Abstract
In the current work, a fully implicit numerical integration scheme is developed for modeling twinning-induced plasticity using a crystal plasticity framework. Firstly, the constitutive formulation of a twin-based micromechanical model is presented to estimate the deformation behavior of steels with low stacking fault energy. Secondly, a numerical integration scheme is developed for discretizing constitutive equations through a fully implicit time integration scheme using the backward Euler method. A time sub-stepping algorithm and the two-norm convergence criterion are used to regulate time step size and stopping criterion. Afterward, a numerical scheme is implemented in finite element software ABAQUS as a user-defined material subroutine. Finally, finite element simulations are executed for observing the validity, performance, and limitations of the numerical scheme. It is observed that the simulation results are in good agreement with the experimental observations with a maximum error of 16% in the case of equivalent stress and strain. It is also found that the developed model is able to estimate well the deformation behavior, magnitude of slip and twin shear strains, and twin volume fraction of three different TWIP steels where the material point is subjected to tension and compression.