An algorithm for finite element analysis of problems in elastoplastici
ty with continuous stress and strain approximation is presented. By a
global iteration procedure, equilibrium is preserved at the nodes in a
weak sense, and the local constitutive relation between stresses and
strains is satisfied. A high order numerical integration is used to ac
hieve a good quality stiffness matrix and to evaluate the boundary bet
ween elastic and plastic regions in the case of partly plastic element
s.