Finite element stress formulation for dynamic elastic-plastic analysis

The solution to wave propagation problems in solids with elastic-plastic material properties is obtained by using the finite element method directly in terms of the stresses. A variational principle due to Gurtin is modified by including a plastic strain tensor in the constitutive relationship. The resulting finite element equations, which represent the strain-displacement equations written in terms of the stresses, are simultaneous integral equations in time. With a transformation of variables, a set of simultaneous differential equations is obtained of the form$\text{H}\text{s}\text{?}\text{+ Qs+ Ve}{}_{\mathrm{p}}\text{= q(t)}$, where $\text{H}$ is a symmetric positive-semidefinite matrix, and $\text{Q}$ is a symmetric positive-definite matrix. The stresses and the plastic strains are represented by $\text{s}\text{?}$ and $\text{e}{}_{\mathrm{p}}$, respectively.Finite element equations are developed for an axisymmetric ring element with an arbitrary quadrilateral cross section in which the stresses and the plastic strains vary linearly along the sides of the elements. The equations are numerically integrated with respect to time by Newmark's generalized acceleration method.An iterative procedure is presented, which uses the finite element strain-displacement equations and the plasticity relationships, to determine the state of stress at the end of the time step. Several examples are used to demonstrate the solution technique for elastic and elastic-plastic problems.