The brushless doubly-fed machine (BDFM) has two stator windings with different pole numbers, supplied with different frequencies. Therefore, the distribution of magnetic fields in stator and rotor iron is complex. In addition, the stator flux density distribution is non-sinusoidal and has a DC offset at the natural speed. This makes the use of conventional hysteresis models utilised for sinusoidal fields impractical for the BDFM. In this paper a new hysteresis model is proposed for the BDFM stator iron based on the scalar Preisach model. The rotational characteristics of the magnetic fields in the BDFM are also considered and their effects in generating iron losses are assessed. 2-D time-stepping finite element (FE) models are developed for a prototype D160 BDFM to estimate iron losses and are validated by experiments.