Computational modelling of plasticity induced by martensitic phase transformations

A time integration scheme is presented for the martensitic phase transformation model developed in recent theoretical work of Turteltaub and Suiker (A multi-scale thermomechanical model for cubic to tetragonal martensitic phase transformations 2005; Transformation-induced plasticity in ferrous alloys 2005). The phase transformation model can be used for analysing transformation-induced plasticity (TRIP) phenomena in ferrous alloys. The microstructural information for the phase transformation model is provided by the crystallographic theory of martensitic transformations. The transformation characteristics depend on the specific transformation systems activated during a loading process. The time integration scheme is formulated within a framework of finite deformations, where the stress-update algorithm is based on a fully implicit Euler backward discretization. A robust search algorithm is used for detecting the transformation systems activated during loading. The completion of the transformation process is prescribed by a constraint on the total martensitic volume fraction, which is accurately satisfied in the converged state using a sub-stepping algorithm. The computation of the consistent tangent operator is performed through a numerical differentiation method, which avoids the determination of extensive analytical derivatives and allows the model to be easily adapted if necessary. The ability of the algorithm to solve complex transformation-induced plasticity problems is illustrated with the aid of three-dimensional analyses, in which an aggregate of single-crystal grains of retained austenite embedded in a ferrite-based matrix is subjected to uniaxial tension. The response characteristics are in agreement with experimental findings of Jacques et al.