This study is aimed at developing an anisotropic elastoplastic constitutive model for geomaterials at finite strain and its stress calculation algorithm based on the fully implicit return-mapping scheme. The Cam-clay plasticity model is adopted as a specific prototype model of geomaterials. As a pertinent representation of deformation-induced anisotropy in geomaterials, nonlinear rotational hardening is incorporated into the model in a theoretically reasonable manner by introducing the dual multiplicative decompositions of the deformation gradient tensor. In addition to the usual decomposition into elastic and plastic parts, the plastic part is decomposed further into a part contributing to the rotational hardening and a remainder part. The former part leads to a back stress ratio tensor related to the rotational hardening via a hyperelastic-type hardening rule. The constitutive theory is thereby formulated on proper intermediate configurations entirely in terms of deformation-like tensorial variables possessing invariance property, without resort to any objective rates of stress or stress-like variables. Combining the Cam-clay plasticity with the concept of subloading surface, a class of unconventional plasticity, enables the model to be capable of reproducing complex hardening/softening accompanied by volumetric contractive/dilative responses.Basic characteristics and predictive capability of the proposed model, as well as the accuracy of the developed numerical scheme, are verified through several numerical examples including monotonic and cyclic loadings.
K E Y W O R D Sanisotropy, cam-clay plasticity, finite strain, return-mapping algorithm, rotational hardening, subloading surface conceptThis is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.