An improved return-map stress update algorithm for finite deformation analysis of general isotropic elastoplastic geomaterials

This paper develops a novel return mapping algorithm for the numerical integration of general isotropic finite strain elastoplastic constitutive models for geomaterials. The constitutive formulation is founded on multiplicative decomposition of the deformation gradient. The logarithmic strain measure as well as the exponential approximation of the plastic flow rule is utilized to restore the standard infinitesimal format return mapping algorithm. Central to the algorithm is the exploitation of a set of three mutually orthogonal unit base tensors for the representation of constitutive relations and the corresponding integration of the rate form of the constitutive equations. The base tensors constitute a local cylindrical coordinate system in the principal space, which allows to formulate the return mapping algorithm in the three-dimensional space and reduce the dimension of the problem to be analyzed from six down to three. With the proposed approach, direct determination of the principal axes and the transformation procedure between the general space and the principal space, as required in traditional spectral decomposition, are avoided. Furthermore, the matrices that are involved in the inversion evaluation take simple forms, leading to extremely easy inverse computation. As a result, the consistent tangent operator can be streamlined into a form simpler and more compact than those by conventional integration methods. Following the formulation of the integration procedure, a numerical experiment is performed to assess the accuracy and efficiency of the proposed algorithm. Copyright (c) 2014 John Wiley & Sons, Ltd.