This page shows my version of a frequency dependent hysteresis model showing a typical behavior like it can be observed with special magnetic materials. Unlike the current pSPICE model, it has full frequency dependency and even contains Barkhausen's effect.

static low frequency behavior with moderate saturation

The coarseness of the Barkhausen steps, the steepness of the saturation region and the saturation knee are entirely parametric and can be perfectly adjusted to the four main field parameters – coercive field strength, saturation magnetisation and remanence. The Barkhausen steps are deterministically coupled to the settings and the strength of the magnetic field and are also generated by stochastic processes using a noise generator. The jumping of the domain walls is implemented using resonator models. Frequency dependence is firmly integrated into the model and exhibits natural behaviour across a wide frequency range without the need for parameter adjustments. However, it can be adjusted to match experimental measurements though.

dynamic frequency behavior with stronger saturation  

Typical elliptic behavior of the B-H-curve when it comes to higher frequencies is e.g. because of delayed reaction of the Weiss areas and the magnetic domain wall pinning. The increasingly elliptical to round B-H loop at higher frequencies takes into account, the various causes, including dynamic magnetisation processes, turbulence currents, and stray losses.

Read more about my magnetic research here Optimized Parameters for the Jiles Atherton Model