Implementation of subloading surface model for hyperelastoplasticity with nonlinear kinematic/isotropic hardening based on reference and intermediate configurations

In this paper, the elastoplastic behavior of isotropic materials under finite strains is modelled and simulated under consideration of the subloading surface concept. Using the subloading surface model has many advantages, such as the ability of modeling more smoothly transition from elastic to plastic state and also more realistic modeling of the Bauschinger effect. For this purpose, three time integration algorithms are utilized in the reference and intermediate configurations in the present study. The first algorithm is based on a time integration scheme, available in the literature, and the second time integration is proposed by tensor exponential-based time integration and the third one is represented by a standard backward Euler time integration scheme. In the second proposed time integration, which is defined in the intermediate configuration, there is no need to use corotational objective time rates. However, in the third one, although it is based on the intermediate configuration, the Zaremba-Jaumann co-rotational rate is employed. Besides, in the second two time integrations, the constitutive equations are derived based on the unknown tensors defined in the intermediate configuration and as a consequence there is no restriction of using different hyperelastic models. For case studies, diverse deformation gradient tensors with and without volume changes are considered. These studied cases are accomplished for SUS 304 stainless steel and a glassy polymer (oriented PET). The obtained results of some cases are compared with available results in the literature. The agreement among the predicted results in this paper and the results from literature is good. & COPY; 2023 Elsevier Inc. All rights reserved.