Implicit numerical integration of highly nonlinear plasticity models

When using an implicit numerical integration algorithm to implement a highly nonlinear elasto-plastic constitutive models in FEM simulations, some difficulties would encounter including Jacobian matrix singularity and non-convergence. Based on the composition of the elastoplastic constitutive model and the characteristics of the traditional implicit algorithms, this paper proposes three improved implicit algorithms, namely the Homotopy-Newton-CPPM algorithm, the CG-Newton-CPPM algorithm and the two-stage algorithm, which would contribute to improve the convergence and deal with singularity of the Jacobian matrix. Besides, the convergence and the computational efficiency as well as the accuracy of five algorithms, including a conventional explicit algorithm, a conventional implicit algorithm and three improved implicit algorithms, are compared with reference to numerical simulations of single element tests and multi-element analysis.