PurposeMagnetic hysteresis holds significant technical and physical importance in the design of electromagnetic components. Despite extensive research in this area, modeling magnetic hysteresis remains a challenging task that is yet to be fully resolved. The purpose of this paper is to study vector hysteresis play models for anisotropic ferromagnetic materials in a physical, thermodynamical approach.Design/methodology/approachIn this work, hysteresis play models are implemented to interpret magnetic properties, drawing upon classical rate-independent plasticity principles derived from continuum mechanics theory. By conducting qualitative and quantitative verification and validation, various aspects of ferromagnetic vector hysteresis were thoroughly examined. By directly incorporating the hysteresis play models into the primal formulations using fixed point method, the proposed model is validated with measurements in a finite element (FE) environments.FindingsThe proposed vector hysteresis play model is verified with fundamental properties of hysteresis effects. Numerical analysis is performed in an FE environment. Measured data from a rotational single sheet tester (RSST) are validated to the simulated results.Originality/valueThe results of this work demonstrates that the essential properties of the hysteresis effects by electrical steel sheets can be represented by the proposed vector hysteresis play models. By incorporation of hysteresis play models into the weak formulations of the magnetostatic problem in the h-based magnetic scalar potential form, magnetic properties of electrical steel sheets can be locally analyzed and represented.