An integration algorithm and the corresponding consistent tangent operator for fully coupled elastoplastic and damage equations

This paper deals with the implementation in finite element calculations of complex elastoplastic constitutive equations exhibiting non-linear kinematic and isotropic hardening, and fully coupled to a continuous ductile damage evolution model. The Newton method is used to solve the non-linear global equilibrium equations as well as the non-linear local equations obtained by fully implicit integration of the constitutive equations. The consistent local tangent modulus is obtained by exact linearization of the algorithm. The procedure described has been implemented in the general purpose code ABAQUS.