We consider the quasi-static magnetic hysteresis model based on a dry-friction like representation of magnetization.The model has a consistent energy interpretation, is intrinsically vectorial, and ensures a direct calculation of the stored and dissipated energies at any moment in time, and hence not only on the completion of a closed hysteresis loop. We discuss the variational formulation of this model and derive an efficient numerical scheme, avoiding the usually employed approximation which can be inaccurate in the vectorial case. The parameters of this model for a nonoriented steel are identified using a set of first order reversal curves. Finally, the model is incorporated as a local constitutive relation into a 2D finite element simulation accounting for both the magnetic hysteresis and the eddy current.Index Terms-vector magnetic hysteresis, energy-based model, variational formulation, finite element simulation.Corresponding author: L. Prigozhin (email: leonid@math.bgu.ac.il).Both the Preisach and Jiles-Atherton models are scalar and, although there exist numerous vector modifications, these also lack a true physical justification. Furthermore, in a general situation, the use of these models to predict the evolving magnetization does not make computing the accompanying energy loss straightforward (see, e.g., [6]). In a seminal work [7], Bergqvist proposed a new quasi-static magnetic hysteresis model, phenomenological but having a consistent and genuine energy interpretation, intrinsically vectorial, and ensuring a direct calculation of the stored magnetic energy and the dissipated energy at any moment in time, and not only after the completion of a closed hysteresis loop as is usually the case. This model differs significantly from the previous