A crystal plasticity model exhibiting deformation twinning is introduced based on the classical crystal plasticity theory,and its numerical implementation is conducted in which a fully implicit integration procedure is employed. In this constitutive model,the saturation-type hardening law is adopted and the direct proportion relationship between twin resistance and slip hardening is used to describe the effect of twinning on slip hardening and twin hardening. By considering the physical meaning of all the 1...