A general yield function with differential and anisotropic hardening for strength modelling under various stress states with non-associated flow rule

•A function is proposed to model differential hardening under various states.•Convexity is calculated by geometric-inspired numerical convex analysis method.•Parameters are analytically computed to guarantee accuracy and user-friendliness.•Yield surface evolution is automatically captured in distinct stress states.