EXPONENTIAL MAP ALGORITHM FOR STRESS UPDATES IN ANISOTROPIC MULTIPLICATIVE ELASTOPLASTICITY FOR SINGLE CRYSTALS

This paper presents a new stress update algorithm for large-strain rate-independent single-crystal plasticity. The theoretical frame is the well-established continuum slip theory based on the multiplicative decomposition of the deformation gradient into elastic and plastic parts. A distinct feature of the present formulation is the introduction and computational exploitation of a particularly simple hyperelastic stress response function based on a further multiplicative decomposition of the elastic deformation gradient into spherical and unimodular parts, resulting in a very convenient representation of the Schmid resolved shear stresses on the crystallographic slip systems in terms of a simple inner product of Eulerian vectors. The key contribution of this paper is an algorithmic formulation of the exponential map exp : sI(3) -+ SL(3) for updating the special linear group SL(3) of unimodular plastic deformation maps. This update preserves exactly the plastic incompressibility condition of the anisotropic plasticity model under consideration. The resulting fully implicit stress update algorithm treats the possibly redundant constraints of singlecrystal plasticity by means of an active set search. It exploits intrinsically the simple representation of the Schmid stresses by formulating the retum algorithm and the associated consistent elastoplastic moduli in terms of Eulerian vectors updates. The performance of the proposed algorithm is demonstrated by means of a representative numerical example.