A penalty finite element method has been developed for accurately predicting stress distributions at the tool-workpiece interfaces. The basic formulation is described, with the emphasis on the algorithm to deal with the normal stress and the frictional stress at the interface. Comparison with the experimental data and the theoretical solutions found in the literature is made for the forming processes selected.