We have used a stress integration scheme, based on numerical approximation of the yield function gradients, to implement in the finite element code ABAQUS three elastic isotropic, plastic anisotropic constitutive models with yielding described by Yld2004-18p (Barlat et al., 2005), CPB06ex2 (Plunkett et al., 10 2008) and Yld2011-27p (Aretz and Barlat, 2013) criteria, respectively. We have developed both VUMAT and UMAT subroutines for the three constitutive models.

You can download the codes here: https://edatos.consorciomadrono.es/dataset.xhtml?persistentId=doi:10.21950/IBUQ4Q