Quickly convergent integration methods for plane stress plasticity

Explicit finite-element programs, such as DYNA3D, are frequently used to solve large structural problems, such as automobile impacts, that involve tens of thousands of shell elements. Since the explicit finite-element programs do not use a stiffness matrix, the major cost of the calculation is the evaluation of the force vector. A significant part of that cost is the evaluation of the constitutive models. The standard algorithm for integrating J2 plasticity at the stress points is the radial return method. An iterative variation of radial return is commonly used to impose the zero normal stress condition within the shell elements. In this paper we explore several different non-iterative methods that impose the zero normal stress condition in an approximate manner. We show that the cost of the constitutive evaluation can be reduced by 25 per cent without introducing significant errors.