Jiles–Atherton model

The Jiles–Atherton model of magnetic hysteresis was introduced in 1984 by David Jiles and D. L. Atherton.[1] This is one of the most popular models of magnetic hysteresis. Its main advantage is the fact that this model enables connection with physical parameters of the magnetic material.[2] Jiles–Atherton model enables calculation of minor and major hysteresis loops.[1] The original Jiles–Atherton model is suitable only for isotropic materials.[1] However, an extension of this model presented by Ramesh et al.[3] and corrected by Szewczyk [4] enables the modeling of anisotropic magnetic materials.

Principles

Magnetization ${\displaystyle M}$ of the magnetic material sample in Jiles–Atherton model is calculated in the following steps [1] for each value of the magnetizing field ${\displaystyle H}$:

• effective magnetic field ${\displaystyle H_{\text{e}}}$ is calculated considering interdomain coupling ${\displaystyle \alpha }$ and magnetization ${\displaystyle M}$,
• anhysteretic magnetization ${\displaystyle M_{\text{an}}}$ is calculated for effective magnetic field ${\displaystyle H_{\text{e}}}$,
• magnetization ${\displaystyle M}$ of the sample is calculated by solving ordinary differential equation taking into account sign of derivative of magnetizing field ${\displaystyle H}$ (which is the source of hysteresis).

Parameters

Original Jiles–Atherton model considers following parameters:[1]

Parameter	Units	Description
${\displaystyle \alpha }$		Quantifies interdomain coupling in the magnetic material
${\displaystyle a}$	A/m	Quantifies domain walls density in the magnetic material
${\displaystyle M_{\text{s}}}$	A/m	Saturation magnetization of material
${\displaystyle k}$	A/m	Quantifies average energy required to break pinning site in the magnetic material
${\displaystyle c}$		Magnetization reversibility

Extension considering uniaxial anisotropy introduced by Ramesh et al.[3] and corrected by Szewczyk [4] requires additional parameters:

Parameter	Units	Description
${\displaystyle K_{\text{an}}}$	J/m^3	Average anisotropy energy density
${\displaystyle \psi }$	rad	Angle between direction of magnetizing field ${\displaystyle H}$ and direction of anisotropy easy axis
${\displaystyle t}$		Participation of anisotropic phase in the magnetic material

Modelling the magnetic hysteresis loops

Effective magnetic field

Effective magnetic field ${\displaystyle H_{\text{e}}}$ influencing on magnetic moments within the material may be calculated from following equation:[1]

${\displaystyle H_{\text{e}}=H+\alpha M}$

This effective magnetic field is analogous to the Weiss mean field acting on magnetic moments within a magnetic domain.[1]

Anhysteretic magnetization

Anhysteretic magnetization can be observed experimentally, when magnetic material is demagnetized under the influence of constant magnetic field. However, measurements of anhysteretic magnetization are very sophisticated due to the fact, that the fluxmeter has to keep accuracy of integration during the demagnetization process. As a result, experimental verification of the model of anhysteretic magnetization is possible only for materials with negligible hysteresis loop.[4]
Anhysteretic magnetization of typical magnetic material can be calculated as a weighted sum of isotropic and anisotropic anhysteretic magnetization:[5]

${\displaystyle M_{\text{an}}=(1-t)M_{\text{an}}^{\text{iso}}+tM_{\text{an}}^{\text{aniso}}}$

Isotropic

Isotropic anhysteretic magnetization ${\displaystyle M_{\text{an}}^{\text{iso}}}$ is determined on the base of Boltzmann distribution. In the case of isotropic magnetic materials, Boltzmann distribution can be reduced to Langevin function connecting isotropic anhysteretic magnetization with effective magnetic field ${\displaystyle H_{\text{e}}}$:[1]

${\displaystyle M_{\text{an}}^{\text{iso}}=M_{\text{s}}\left(\coth \left({\frac {H_{\text{e}}}{a}}\right)-{\frac {a}{H_{\text{e}}}}\right)}$

Anisotropic

Anisotropic anhysteretic magnetization ${\displaystyle M_{\text{an}}^{\text{aniso}}}$ is also determined on the base of Boltzmann distribution.[3] However, in such a case, there is no antiderivative for Boltzmann distribution function.[4] For this reason, integration has to be made numerically. In the original publication, anisotropic anhysteretic magnetization ${\displaystyle M_{\text{an}}^{\text{aniso}}}$ is given as:[3]

${\displaystyle M_{\text{an}}^{\text{aniso}}=M_{\text{s}}{\frac {\int _{0}^{\pi }\!e^{E(1)+E(2)}\sin(\theta )\cos(\theta )\,d\theta }{\int _{0}^{\pi }\!e^{E(1)+E(2)}\sin(\theta )\,d\theta }}}$

where

${\displaystyle E(1)={\frac {H_{\text{e}}}{a}}\cos \theta -{\frac {K_{\text{an}}}{M_{\text{s}}\mu _{0}a}}\sin ^{2}(\psi -\theta )}$

${\displaystyle E(2)={\frac {H_{\text{e}}}{a}}\cos \theta -{\frac {K_{\text{an}}}{M_{\text{s}}\mu _{0}a}}\sin ^{2}(\psi +\theta )}$

It should be highlighted, that typing mistake happened in the original Ramesh et al. publication.[4] As a result, for isotropic material (where ${\displaystyle K_{\text{an}}=0)}$), presented form of anisotropic anhysteretic magnetization ${\displaystyle M_{\text{an}}^{\text{aniso}}}$ is not coherent with isotropic anhysteretic magnetization ${\displaystyle M_{\text{an}}^{\text{iso}}}$ given by Langevin equation. Physical analysis leads to conclusion, that equation for anisotropic anhysteretic magnetization ${\displaystyle M_{\text{an}}^{\text{aniso}}}$ has to be corrected to the following form:[4]

${\displaystyle M_{\text{an}}^{\text{aniso}}=M_{\text{s}}{\frac {\int _{0}^{\pi }\!e^{0.5(E(1)+E(2))}\sin(\theta )\cos(\theta )\,d\theta }{\int _{0}^{\pi }\!e^{0.5(E(1)+E(2))}\sin(\theta )\,d\theta }}}$

In the corrected form, model for anisotropic anhysteretic magnetization ${\displaystyle M_{\text{an}}^{\text{aniso}}}$ was confirmed experimentally for anisotropic amorphous alloys.[4]

Magnetization as a function of magnetizing field

In Jiles–Atherton model, M(H) dependence is given in form of following ordinary differential equation:[6]

${\displaystyle {\frac {dM}{dH}}={\frac {1}{1+c}}{\frac {M_{\text{an}}-M}{\delta k-\alpha (M_{\text{an}}-M)}}+{\frac {c}{1+c}}{\frac {dM_{\text{an}}}{dH}}}$

where ${\displaystyle \delta }$ depends on direction of changes of magnetizing field ${\displaystyle H}$ (${\displaystyle \delta =1}$ for increasing field, ${\displaystyle \delta =-1}$ for decreasing field)

Flux density as a function of magnetizing field

Flux density ${\displaystyle B}$ in the material is given as:[1]

${\displaystyle B(H)=\mu _{0}M(H)}$

where ${\displaystyle \mu _{0}}$ is magnetic constant.

Vectorized Jiles–Atherton model

Vectorized Jiles–Atherton model is constructed as the superposition of three scalar models one for each principal axe.[7] This model is especially suitable for finite element method computations.

Numerical implementation

The Jiles-Atherton model is implemented in JAmodel, a MATLAB/OCTAVE toolbox. It uses the Runge-Kutta algorithm for solving ordinary differential equations. JAmodel is open-source is under MIT license.[8]

The two most important computational problems connected with the Jiles–Atherton model were identified:[8]

• numerical integration of the anisotropic anhysteretic magnetization ${\displaystyle M_{\text{an}}^{\text{aniso}}}$
• solving the ordinary differential equation for ${\displaystyle M(H)}$ dependence.

For numerical integration of the anisotropic anhysteretic magnetization ${\displaystyle M_{\text{an}}^{\text{aniso}}}$ the Gauss–Kronrod quadrature formula has to be used. In GNU Octave this quadrature is implemented as quadgk() function.

For solving ordinary differential equation for ${\displaystyle M(H)}$ dependence, the Runge–Kutta methods are recommended. It was observed, that the best performing was 4-th order fixed step method.[8]

Further development

Since its introduction in 1984, Jiles–Atherton model was intensively developed. As a result, this model may be applied for the modeling of:

• frequency dependence of magnetic hysteresis loop in conductive materials [9][10]
• influence of stresses on magnetic hysteresis loops [11][12][13]
• magnetostriction of soft magnetic materials [11][14]

Moreover, different corrections were implemented, especially:

• to avoid unphysical states when reversible permeability is negative [15]
• to consider changes of average energy required to break pinning site [16]

Applications

Jiles–Atherton model may be applied for modeling:

• rotating electric machines [17]
• power transformers [18]
• magnetostrictive actuators [19]
• magnetoelastic sensors [20][21]
• magnetic field sensors (e. g. fluxgates) [22][23]

It is also widely used for electronic circuit simulation, especially for models of inductive components, such as transformers or chokes.[24]

References

1. Jiles, D. C.; Atherton, D.L. (1984). "Theory of ferromagnetic hysteresis". Journal of Applied Physics. 55 (6): 2115. Bibcode:1984JAP....55.2115J. doi:10.1063/1.333582.
2. Liorzou, F.; Phelps, B.; Atherton, D. L. (2000). "Macroscopic models of magnetization". IEEE Transactions on Magnetics. 36 (2): 418. Bibcode:2000ITM....36..418L. doi:10.1109/20.825802.
3. Ramesh, A.; Jiles, D. C.; Roderick, J. M. (1996). "A model of anisotropic anhysteretic magnetization". IEEE Transactions on Magnetics. 32 (5): 4234. Bibcode:1996ITM....32.4234R. doi:10.1109/20.539344.
4. Szewczyk, R. (2014). "Validation of the anhysteretic magnetization model for soft magnetic materials with perpendicular anisotropy". Materials. 7 (7): 5109–5116. Bibcode:2014Mate....7.5109S. doi:10.3390/ma7075109. PMC 5455830. PMID 28788121.
5. Jiles, D.C.; Ramesh, A.; Shi, Y.; Fang, X. (1997). "Application of the anisotropic extension of the theory of hysteresis to the magnetization curves of crystalline and textured magnetic materials". IEEE Transactions on Magnetics. 33 (5): 3961. Bibcode:1997ITM....33.3961J. doi:10.1109/20.619629.
6. Jiles, D. C.; Atherton, D.L. (1986). "A model of ferromagnetic hysteresis". Journal of Magnetism and Magnetic Materials. 61 (1–2): 48. Bibcode:1986JMMM...61...48J. doi:10.1016/0304-8853(86)90066-1.
7. Szymanski, Grzegorz; Waszak, Michal (2004). "Vectorized Jiles–Atherton hysteresis model". Physica B. 343 (1–4): 26–29. Bibcode:2004PhyB..343...26S. doi:10.1016/j.physb.2003.08.048.
8. Szewczyk, R. (2014). Computational problems connected with Jiles–Atherton model of magnetic hysteresis. Advances in Intelligent Systems and Computing. 267. pp. 275–283. doi:10.1007/978-3-319-05353-0_27. ISBN 978-3-319-05352-3.
9. Jiles, D.C. (1994). "Modelling the effects of eddy current losses on frequency dependent hysteresis in electrically conducting media". IEEE Transactions on Magnetics. 30 (6): 4326–4328. Bibcode:1994ITM....30.4326J. doi:10.1109/20.334076.
10. Szewczyk, R.; Frydrych, P. (2010). "Extension of the Jiles–Atherton model for modelling the frequency dependence of magnetic characteristics of amorphous alloy cores for inductive components of electronic devices". Acta Physica Polonica A. 118 (5): 782. doi:10.12693/aphyspola.118.782.
11. Sablik, M.J.; Jiles, D.C. (1993). "Coupled magnetoelastic theory of magnetic and magnetostrictive hysteresis". IEEE Transactions on Magnetics. 29 (4): 2113. Bibcode:1993ITM....29.2113S. doi:10.1109/20.221036.
12. Szewczyk, R.; Bienkowski, A. (2003). "Magnetoelastic Villari effect in high-permeability Mn-Zn ferrites and modeling of this effect". Journal of Magnetism and Magnetic Materials. 254: 284–286. Bibcode:2003JMMM..254..284S. doi:10.1016/S0304-8853(02)00784-9.
13. Jackiewicz, D.; Szewczyk, R.; Salach, J.; Bieńkowski, A. (2014). "Application of extended Jiles–Atherton model for modelling the influence of stresses on magnetic characteristics of the construction steel". Acta Physica Polonica A. 126 (1): 392. doi:10.12693/aphyspola.126.392.
14. Szewczyk, R. (2006). "Modelling of the magnetic and magnetostrictive properties of high permeability Mn-Zn ferrites". Pramana. 67 (6): 1165–1171. Bibcode:2006Prama..67.1165S. doi:10.1007/s12043-006-0031-z.
15. Deane, J.H.B. (1994). "Modeling the dynamics of nonlinear inductor circuits". IEEE Transactions on Magnetics. 30 (5): 2795–2801. Bibcode:1994ITM....30.2795D. doi:10.1109/20.312521.
16. Szewczyk, R. (2007). "Extension of the model of the magnetic characteristics of anisotropic metallic glasses". Journal of Physics D: Applied Physics. 40 (14): 4109–4113. Bibcode:2007JPhD...40.4109S. doi:10.1088/0022-3727/40/14/002.
17. Du, Ruoyang; Robertson, Paul (2015). "Dynamic Jiles–Atherton Model for Determining the Magnetic Power Loss at High Frequency in Permanent Magnet Machines". IEEE Transactions on Magnetics. 51 (6): 7301210. Bibcode:2015ITM....5182594D. doi:10.1109/TMAG.2014.2382594.
18. Huang, Sy-Ruen; Chen, Hong-Tai; Wu, Chueh-Cheng; et al. (2012). "Distinguishing internal winding faults from inrush currents in power transformers using Jiles–Atherton model parameters based on correlation voefficient". IEEE Transactions on Magnetics. 27 (2): 548. doi:10.1109/TPWRD.2011.2181543.
19. Calkins, F.T.; Smith, R.C.; Flatau, A.B. (2008). "Energy-based hysteresis model for magnetostrictive transducers". IEEE Transactions on Magnetics. 36 (2): 429. Bibcode:2000ITM....36..429C. CiteSeerX 10.1.1.44.9747. doi:10.1109/20.825804.
20. Szewczyk, R.; Bienkowski, A. (2004). "Application of the energy-based model for the magnetoelastic properties of amorphous alloys for sensor applications". Journal of Magnetism and Magnetic Materials. 272: 728–730. Bibcode:2004JMMM..272..728S. doi:10.1016/j.jmmm.2003.11.270.
21. Szewczyk, R.; Salach, J.; Bienkowski, A.; et al. (2012). "Application of extended Jiles–Atherton model for modeling the magnetic characteristics of Fe41.5Co41.5Nb3Cu1B13 alloy in as-quenched and nanocrystalline State". IEEE Transactions on Magnetics. 48 (4): 1389. Bibcode:2012ITM....48.1389S. doi:10.1109/TMAG.2011.2173562.
22. Szewczyk, R. (2008). "Extended Jiles–Atherton model for modelling the magnetic characteristics of isotropic materials". Acta Physica Polonica A. 113 (1): 67. Bibcode:2008JMMM..320E1049S. doi:10.12693/APhysPolA.113.67.
23. Moldovanu, B.O.; Moldovanu, C.; Moldovanu, A. (1996). "Computer simulation of the transient behaviour of a fluxgate magnetometric circuit". Journal of Magnetism and Magnetic Materials. 157-158: 565–566. Bibcode:1996JMMM..157..565M. doi:10.1016/0304-8853(95)01101-3.
24. Cundeva, S. (2008). "Computer simulation of the transient behaviour of a fluxgate magnetometric circuit" (PDF). Serbian Journal of Electrical Engineering. 5 (1): 21–30. doi:10.2298/sjee0801021c. Archived from the original (PDF) on 2014-07-24.

External links

• Jiles–Atherton model for Octave/MATLAB - open-source software for implementation of Jiles–Atherton model in GNU Octave and Matlab
• GitHub repository for open-source implementation of Jiles–Atherton model - GitHub repository for open-source implementation of Jiles–Atherton model in GNU Octave and Matlab