This paper deals with the implementation of an anisotropic plasticity constitutive equations exhibiting non-linear isotropic, kinematic hardening and coupled continuum ductile damage models. A fully implicit integration of the coupling constitutive equations is adopted and leads to two non-linear local scalar equations solved using the Newton method. The consistent local tangent modulus is obtained in a closed form by exact linearization of the algorithm. The numerical treatment of the proposed algorithm is implemented on ABAQUS using user interface material subroutines (UMAT and VUMAT). The performance of the present algorithm is assessed and bared out by numerical examples.
