A new finite strain elatoplastic J2-flow model with coupling effects of both isotropic and anisotropic hardening is proposed with the co-rotational logarithmic rate.In terms of certain single-variable shape functions representing uniaxial loading and unloading curves,explicit multi-axial expressions for the three hardening quantities incorporated in the new model proposed are derived in unified forms for the purpose of automatically and accurately simulating complex pseudoelastic-to-plastic transition effects of shape memory alloys (SMAs) under multiple loading-unloading cycles.Numerical examples show that with only a single parameter of direct physical meaning for each cycle,accurate and explicit simulations may be achieved for extensive data from multiple cycle tests.