This report proposes the iterative method for the dynamic thermo-plastic-elastic stress analysis by the finite element method (FEM), of which scheme satisfies the condition that the stress point must move on the yield surface during plastic deformation. The common engineering condition of a von Mises, isotropic hardening model and the constitutive equatin expressed as a function of a strain rate-temperature parameter R are used in this FEM program. As a example of calculation using by this FEM, we show the strain and temperature distribution around a hole for a plate loaded in the uni-axial direction.