ABSTRACT
A return mapping (fully implicit) algorithm suitable for general isotropic and anisotropic hyperelastoplasticity is developed. The algorithm is cast in the framework of multiplicative decomposition and has no restriction to hyperelastic isotropy, isotropic flow function, or isotropic yield function. As such, it can be applied to both isotropic and anisotropic large-deformation problems. Numerical examples illustrating good performance are shown.