sdx (2*sx_smsld-sy_smsld-sz_smsld)/3 Stress deviator x-component sdy (2*sy_smsld-sx_smsld-sz_smsld)/3 Stress deviator y-component sdz (2*sz_smsld-sx_smsld-sy_smsld)/3 Stress deviator z-component ep_a A_a*exp(-Q_R_a/T)*sinh(xi_a*mises_smsld/s_a)^(1/m_a) Effective inelastic strain rate s_star_a sh_a*(ep_a/A_a*exp(Q_R_a/T))^n_a+eps Saturation value for s_a h_a h0_a*abs(1-s_a/s_star_a)^(a_a-1)*(1-s_a/s_star_a) Hardening function