The paper presents a numerical method for transient magneto-quasistatic field analysis including hysteretic ferromagnetic materials. The procedure is based on a Finite Integration Implicit Time Domain(FI
^2
TD) formulation using the classical Preisach model. An hybrid scheme using the magnetic polarization update scheme for hysteretic materials and the successive approximation scheme for non-hysteretic nonlinear materials is given. Relaxation parameters are introduced to avoid convergence problems of the iterative scheme around the critical zones of the hysteresis curve. In numerical experiments with a test configuration the simulated hysteresis loops are compared with those achieved by measurements and the corresponding hysteretic losses are presented.
