An Easily Used Phenomenological Magnetization Model and Its Empirical Expressions Based on Jiles–Atherton Parameters

In this paper, a simple magnetization model convenient for engineering applications is presented based on the expressions of the first-order LTI system model. Considering the trade-off between the nonlinearity of anhysteretic magnetization and the hysteresis width, the proposed model employs two different equations with different magnetic field amplitudes. Furthermore, the proposed model utilizes the first-order LTI system model with a low magnetic field amplitude and a simple nonlinear function, based on the amplitude–frequency function, with a high magnetic field amplitude. Two important characteristic parameters for engineering applications, namely, amplitude and the equivalent phase lag, were exacted and analyzed to validate the computation precision of the proposed model. Then, the model was verified through comparisons to the validated Jiles–Atherton model. For easy use, similar to a physics-based model instead of a fitting method, empirical expressions for the model parameters were given, and applicable ranges of these equations were determined using the parameters of the Jiles–Atherton model. Finally, an example of the magnetization model applied to an on/off type device was computed to further verify the effectiveness of the proposed model with quite a simple expression.


Introduction
It seems that more complex magnetization models have become more popular in recent years [1][2][3][4][5][6][7][8][9][10][11][12], but this is not good news for the engineering applications of the models.These models directly used the classical model as a sub-model, or performed some adaptive improvements to the classical model to improve the applicability.With quite a strong nonlinearity, especially for the hysteresis, magnetization models are difficult to describe and solve.Moreover, it will bring disaster to the implementation of online control.
Current magnetization models are divided into two types: physics-based models and phenomenological models.Physics-based models consider the magnetizing mechanism, and most parameters of this type of model have certain physical means.A commonly used physics-based magnetization model was proposed by Jiles and Atherton [13][14][15][16], named the Jiles-Atherton (J-A) model.The quasi-static J-A model utilizes five equations, including the differential element and Langevin function, to describe magnetization when complicated expressions provide no possibilities for any analytical solution.Even a high-precision numerical solution was not easily achieved in previous papers when the magnetic field intensity was not high enough [17][18][19].The Smith free-energy model [20,21] employs the mean value of the local magnetization derived from the quadratic Helmholtz model and executes a double integral on the core function.Although it has a simpler expression than the J-A model, the Smith free-energy model is hardly solved, and its accuracy relies on the discontinuous core function [22,23].Without analyzing the magnetizing process, a phenomenological model generally has more concise expressions than a physics-based model.The Preisach model [24], the Duhem model [25,26] and the neural network model [6,[27][28][29][30][31] are the three most commonly used phenomenological models.The Preisach model utilizes a surface integral to describe magnetization, while the Duhem model uses a first-order differential equation and is more easily imposed with identification algorithms and compensation control.The neural network provides many choices for describing magnetization while considering the trade-off between high precision and high efficiency.Besides the above models, physics-based models, such as the Armstrong model [32,33] and the Bouc-Wen model [34,35], and phenomenological models, such as the polynomial model [36], the Krasnoselskii-Pokrovskii (KP) model [37] and the Prandtl-Ishlinskii (PI) model [38,39], have also been employed to describe magnetization [40].
The abovementioned magnetization models have been commonly used in many areas; however, two key problems in engineering applications have not been solved.Firstly, these models have not reached the simplification degree required by engineering applications.Among the abovementioned models, the Duhem model, the Bouc-Wen model (a nonlinear second-order differential equation) and the polynomial model have the most concise forms, while the sign function or absolute value function always remains to distinguish the rising and falling processes that the first derivative at the turning point indicate are not continuous.It seems that these models are "simplified while not completely simplified" and that they are unhelpful for future control.This problem is difficult to solve, resulting in many studies directly covering up the magnetization-magnetic field sub-models and constructing a more macroscopic output-input relationship.Besides the engineering expectation of simplification described below, the verification index for the magnetization model is highly confusing.It seems that either the model supporter or the user does not know how precise the model is required to be or what characteristic parameter represents "precision".In many studies, the models have only been verified by putting together the model and the test results, which indicated that the model is in good agreement with the experiment, as the results looked similar from the perspective of the reader.The unclear requirements of magnetization models have led to the blindness of model applications.
This paper proposed a new model for static magnetization, taking inspiration from the first-order LTI system model.The proposed model fully utilizes the different main features of magnetization under various magnetic field amplitudes on the premise of meeting the strict requirements requested by engineering applications for the computation of amplitude and equivalent phase lag.The model is easily controlled, as it employs the first-order LTI system model under a low magnetic field and a nonlinear function based on the amplitudefrequency function of the first-order LTI system under a high magnetic field.The research idea and framework of the article are shown in Figure 1.In Sections 2 and 3, some necessary information on the J-A model and 1st LTI system model is given and analyzed.Then, the simplification idea and format are supplied in Section 4. The parameters of the simplified model were determined in two methods, respectively, from extra knowledge of the maximum magnetization and just from J-A parameters.Additionally, optimal values and empirical equations for parameter determination are given in Sections 5-7 for best computational effects or convenience.At last, a commonly faced condition of an on-off type device is discussed in Section 8 to give a necessary supplement.Comparisons of the calculated results between the proposed model and verified J-A model were carried out to validate the computation precision of the proposed model.From comparisons, the model reached acceptable computation precisions on the amplitude and phase lag (equivalent) of the magnetization with quite simple expression.

Description and Solution
Xue et al. summarized the various expressions of the Jiles-Atherton hysteresis model cited in different references, and they proposed a reasonable expression based on the mathematical properties of the Langevin function and the actual magnetization process [15,35].The sorted formula is expressed as where Ms is the saturation magnetization; Man is the anhysteretic magnetization; α represents the quantified domain interactions; a is the shape parameter for Man; c is the reversibility coefficient; k is the average energy required to break pinning sites; and δ is the sign function of the derivative of the magnetic field strength with respect to time, expressed by sign(dH/dt) and equal to 1 when dH/dt > 0 and equal to −1 when dH/dt < 0. δM is employed to guarantee positive incremental susceptibilities; it is equal to 0 under the condition of sign[(dH/dt)(Man − M)] < 0, and it is equal to 1 under any other condition.
The analytical solution of Equation ( 1) cannot be obtained, and a high−precision solution is not easy.Ref. [15] supplied a fast, high−precision solving method, where Man is solved by the fixed−point iteration, and then M is solved using fourth−order Runge-Kutta equations.The solving method is expressed by where Man (i) is the value of Man after i times of iterations.The iterative initial value can be imposed as Man (0) = 50 H, and the number of iterations imax ≥ 3, or Man(0) is set as any con-

Description and Solution
Xue et al. summarized the various expressions of the Jiles-Atherton hysteresis model cited in different references, and they proposed a reasonable expression based on the mathematical properties of the Langevin function and the actual magnetization process [15,35].The sorted formula is expressed as where M s is the saturation magnetization; M an is the anhysteretic magnetization; α represents the quantified domain interactions; a is the shape parameter for M an ; c is the reversibility coefficient; k is the average energy required to break pinning sites; and δ is the sign function of the derivative of the magnetic field strength with respect to time, expressed by sign(dH/dt) and equal to 1 when dH/dt > 0 and equal to −1 when dH/dt < 0. δ M is employed to guarantee positive incremental susceptibilities; it is equal to 0 under the condition of sign[(dH/dt)(M an − M)] < 0, and it is equal to 1 under any other condition.
The analytical solution of Equation ( 1) cannot be obtained, and a high-precision solution is not easy.Ref. [15] supplied a fast, high-precision solving method, where M an is solved by the fixed-point iteration, and then M is solved using fourth-order Runge-Kutta equations.The solving method is expressed by where M an (i) is the value of M an after i times of iterations.The iterative initial value can be imposed as M an (0) = 50 H, and the number of iterations i max ≥ 3, or M an (0) is set as any constant, and i max ≥ 5.In Equation ( 2), the differential function value ψ(H,M) Figure 2 shows the description effect of the Jiles-Atherton model on magnetic hysteresis with different parameters.Roughly speaking, the values of k and c mainly influence the width of the hysteresis loop, and an increased k value or a decreased c value is helpful for a narrower hysteresis curve.The values of a and α mainly affect the slope of anhysteretic magnetization, and one obtains steeper slopes with higher values of α or lower values of a.As the scale parameter of M an , the saturation magnetization M s mainly influences the magnification degree of the magnetization value.
stant, and imax ≥ 5.In Equation ( 2), the differential function value Figure 2 shows the description effect of the Jiles-Atherton model on magnetic hysteresis with different parameters.Roughly speaking, the values of k and c mainly influence the width of the hysteresis loop, and an increased k value or a decreased c value is helpful for a narrower hysteresis curve.The values of a and α mainly affect the slope of anhysteretic magnetization, and one obtains steeper slopes with higher values of α or lower values of a.As the scale parameter of Man, the saturation magnetization Ms mainly influences the magnification degree of the magnetization value.

Key Characteristic Parameters in Time−Domain Performance
The nonlinear variation and the width of the hysteresis loop are the key characteristics of magnetization in engineering applications, while the performance of these two characteristics has a trade−off relationship under fixed parameters.That is, a strong nonlinear variation in anhysteretic magnetization is always combined with a relatively narrow hysteresis curve when the magnetization amplitude is high.However, when the magnetization amplitude is low, anhysteretic magnetization shows a weak nonlinear variation, while the hysteresis loop is quite wide.
The amplitude and width of the M−H loop approximately correspond to the maximum value and the equivalent phase lag of time−domain magnetization.To demonstrate the trade−off relationship between the nonlinear variation and the width of the hysteresis loop, time−domain curves of the magnetic field and magnetization are shown in Figure 3, and the corresponding M−H curves are shown in Figure 4.The "phase lag" of the magnetization was not directly defined in the hysteresis while could be reflected by the width of hysteresis loop, coercivity and remanence.Just like the harmonic response, equivalent phase lag in the M−H model can be determined by the coordinate difference of M and H intersections with horizontal axis, as shown in Figure 3.For the M−H curve, the equivalent phase lag can be reflected by the proportion of the H value at M = 0 to Hmax.
When employing a low−amplitude magnetic field, shown by the solid black lines in the figures, the loop was wide, indicating that the time−domain curve of magnetization M had an obvious equivalent phase lag compared with the time−domain curve of the magnetic field H.Meanwhile, the variation in anhysteretic magnetization was approxi-

Key Characteristic Parameters in Time-Domain Performance
The nonlinear variation and the width of the hysteresis loop are the key characteristics of magnetization in engineering applications, while the performance of these two characteristics has a trade-off relationship under fixed parameters.That is, a strong nonlinear variation in anhysteretic magnetization is always combined with a relatively narrow hysteresis curve when the magnetization amplitude is high.However, when the magnetization amplitude is low, anhysteretic magnetization shows a weak nonlinear variation, while the hysteresis loop is quite wide.
The amplitude and width of the M-H loop approximately correspond to the maximum value and the equivalent phase lag of time-domain magnetization.To demonstrate the trade-off relationship between the nonlinear variation and the width of the hysteresis loop, time-domain curves of the magnetic field and magnetization are shown in Figure 3, and the corresponding M-H curves are shown in Figure 4.The "phase lag" of the magnetization was not directly defined in the hysteresis while could be reflected by the width of hysteresis loop, coercivity and remanence.Just like the harmonic response, equivalent phase lag in the M-H model can be determined by the coordinate difference of M and H intersections with horizontal axis, as shown in Figure 3.For the M-H curve, the equivalent phase lag can be reflected by the proportion of the H value at M = 0 to H max .
When employing a low-amplitude magnetic field, shown by the solid black lines in the figures, the loop was wide, indicating that the time-domain curve of magnetization M had an obvious equivalent phase lag compared with the time-domain curve of the magnetic field H.Meanwhile, the variation in anhysteretic magnetization was approximately linear, indicating that the time-domain curve of magnetization M was quite similar to a standard sinusoidal curve, and the M-H curve was not obviously curved.On the contrary, when the H value was quite high, as shown by the blue dotted lines in Figures 3 and 4, the loop of the M-H curve seemed narrow, indicating that M had a few equivalent phase lags compared to H in the time-domain; the magnetization variation was quite nonlinear, indicating that the time-domain curve of M was completely different from a sinusoidal curve; and the M-H curve was obviously curved.
Generally speaking, when the magnetic field strength is high, the variation in the nonlinearity of anhysteretic magnetization is strong, while the equivalent phase lag of magnetization to the magnetic field is little.When the magnetic field strength is low, the equivalent phase lag is high, while the variation in anhysteretic magnetization seems to be linear.This provides the new idea that only the main characteristic of hysteresis is considered to simplify the hysteresis model in engineering applications.Moreover, the trade-off between the nonlinear variation in M an and the width of the hysteresis loop provides the possibility for this type of simplification.
mately linear, indicating that the time−domain curve of magnetization M was quite similar to a standard sinusoidal curve, and the M−H curve was not obviously curved.On the contrary, when the H value was quite high, as shown by the blue dotted lines in Figures 3 and 4, the loop of the M−H curve seemed narrow, indicating that M had a few equivalent phase lags compared to H in the time−domain; the magnetization variation was quite nonlinear, indicating that the time−domain curve of M was completely different from a sinusoidal curve; and the M−H curve was obviously curved.
Generally speaking, when the magnetic field strength is high, the variation in the nonlinearity of anhysteretic magnetization is strong, while the equivalent phase lag of magnetization to the magnetic field is little.When the magnetic field strength is low, the equivalent phase lag is high, while the variation in anhysteretic magnetization seems to be linear.This provides the new idea that only the main characteristic of hysteresis is considered to simplify the hysteresis model in engineering applications.Moreover, the trade−off between the nonlinear variation in Man and the width of the hysteresis loop provides the possibility for this type of simplification.mately linear, indicating that the time−domain curve of magnetization M was quite similar to a standard sinusoidal curve, and the M−H curve was not obviously curved.On the contrary, when the H value was quite high, as shown by the blue dotted lines in Figures 3 and 4, the loop of the M−H curve seemed narrow, indicating that M had a few equivalent phase lags compared to H in the time−domain; the magnetization variation was quite nonlinear, indicating that the time−domain curve of M was completely different from a sinusoidal curve; and the M−H curve was obviously curved.
Generally speaking, when the magnetic field strength is high, the variation in the nonlinearity of anhysteretic magnetization is strong, while the equivalent phase lag of magnetization to the magnetic field is little.When the magnetic field strength is low, the equivalent phase lag is high, while the variation in anhysteretic magnetization seems to be linear.This provides the new idea that only the main characteristic of hysteresis is considered to simplify the hysteresis model in engineering applications.Moreover, the trade−off between the nonlinear variation in Man and the width of the hysteresis loop provides the possibility for this type of simplification.

Description and Solution
The first-order LTI system is the simplest LTI system.Inheriting the parameters of an RL series circuit model, the first-order LTI system, describing the relationship of magnetization M(t) and the magnetic field H(t), can be written as where L and R are two parameters independent of time t.The maximum value of M(t) is constrained by R, L, and saturation magnetization M s .
According to the theory of the first-order LTI system, the steady-state response under where ϕ is the equivalent phase lag of M(t) compared to H(t).ϕ satisfies the conditions that tanϕ = ωL/R and 0 ≤ ϕ < π.
From the solution expressed by Equation (4), M s H max / (Lω) 2 + R 2 determines the amplitude, and ϕ determines the lagging phase of steady-state magnetization.
From Appendix A, the first-order LTI system is capable of describing the hysteresis characteristics of the independent and dependent variables, and it can predict various hysteresis characteristics by adjusting R and L.

Description of Magnetic Hysteresis by the First-Order LTI System Model
According to the results shown in the above sections, it seems feasible to predict the magnetic hysteresis using the first-order LTI system model instead of the Jiles-Atherton model.Furthermore, various hysteresis loops with different widths or amplitudes can also be achieved by adjusting the values of R or L, similar to the influences of the parameters of the Jiles-Atherton model on the hysteresis.
Figure 5 shows the description effects of the first-order LTI system model, including the time-domain and M-H curves, on the hysteresis loop when the amplitude of the magnetic field intensity is low.From the calculation results, it is observed that the first-order LTI system can effectively predict the amplitude and width of the hysteresis loop as long as the appropriate parameters are provided.For the corresponding time-domain curve of magnetization, the curve computed by the Jiles-Atherton model is quite similar to standard sine curves, and the first-LTI system model with the appropriate parameters provides a good predicting effect.The main deviations occur during the variation process from 0 to maximum (or minimum) magnetization.These deviations are not so apparent and are further reduced when calculating the final independent response with other linear or nonlinear links being superimposed.However, when the maximum value of the magnetic field intensity Hmax is high, as shown in Figure 6, the deviations between the results calculated from the first-order LTI system model and those calculated from the Jiles-Atherton model are quite high.With a high magnetic field, the time-domain curve of magnetization deviates far from standard sinusoidal curves.Moreover, under this condition, anhysteretic magnetization M an has a strong nonlinearity, indicating that the hysteresis loop has a large bending, and the first-order LTI system cannot describe a curved hysteresis loop, as analyzed in Section 2.2.

Introduction of Amplitude-Frequency Function
Although the first−LTI system model fails to describe the nonlinear characteristic of Man, mainly the saturation nonlinearity, the amplitude-frequency characteristic of the specific first−order LTI system model shows saturation characteristics.Therefore, the amplitude-frequency function can be employed to predict the nonlinear relationship between M and H when the amplitude of the magnetic field is high.
Consider the function of 22 () . Take this function to describe magnetization as where f[H(t)] is an undetermined function meeting lim ( ) In essence, Equation ( 5) is a simplified fitting function ignoring the hysteresis loop to approximate M, which is similar to using Man to approximate M. The difference here is that Man is expressed implicitly by the Langevin function, while Equation ( 5) supplies quite a simple explicit equation.

Simple Hysteresis Model
A simplification idea can be realized from the discussions in Sections 1.2, 2.2, and 2.3.When the magnetic field intensity is low, the hysteresis loop of magnetization is wide, the variation in the nonlinearity of the anhysteretic magnetization is weak, and magnetization can be illustrated by the first−order LTI system model expressed as Equation (3).In contrast, when the magnetic field intensity is high, the variation in the non-

Introduction of Amplitude-Frequency Function
Although the first-LTI system model fails to describe the nonlinear characteristic of M an , mainly the saturation nonlinearity, the amplitude-frequency characteristic of the specific first-order LTI system model shows saturation characteristics.Therefore, the amplitude-frequency function can be employed to predict the nonlinear relationship between M and H when the amplitude of the magnetic field is high.
Consider the function of Take this function to describe magnetization as where f [H(t)] is an undetermined function meeting lim x→∞ f (x)/x = 0.In essence, Equation ( 5) is a simplified fitting function ignoring the hysteresis loop to approximate M, which is similar to using M an to approximate M. The difference here is that M an is expressed implicitly by the Langevin function, while Equation ( 5) supplies quite a simple explicit equation.

Simple Hysteresis Model
A simplification idea can be realized from the discussions in Sections 3.1-3.3.When the magnetic field intensity is low, the hysteresis loop of magnetization is wide, the variation in the nonlinearity of the anhysteretic magnetization is weak, and magnetization can be illustrated by the first-order LTI system model expressed as Equation (3).In contrast, when the magnetic field intensity is high, the variation in the nonlinearity is strong, the loop is narrow, and magnetization can be described using the simple nonlinear function shown in Equation ( 5).Then, one obtains the magnetization model expressed by where K an is a parameter only related to the material properties (the simplest form of Equation ( 5)), and H r is the split magnetic field used to distinguish the application scopes of the two sub-models.The value of r in H r is a real number belonging to [0, 1], and H r represents the required magnetic field intensity, achieving magnetization with a value of r•M s .From the definition of H r , one obtains H 0 = 0, H 1 = +∞.The hysteresis model expressed by Equation ( 6) considers the width while neglecting the bending state of the hysteresis loop under a low magnetic field intensity.From the point of view of the time-domain magnetization curve, the model considers the phase lag of magnetization while neglecting the variation process from 0 to the maximum or minimum value.On the contrary, when the magnetic field intensity is high, the model considers the bending state while neglecting the width of the loop, which means that the model considers the variation process while neglecting the equivalent phase lag.The proposed model was verified from the comparisons of the calculated results with any verified magnetizations from verified models or tests.Considering the verified model is more convenient than tests to acquire magnetization results of various materials, the physics-based Jiles-Atherton model with abundant parameters was chosen to generate accurate magnetization values for comparisons.

Optimal Parameter for Amplitude
The determination of the values of R and L is a complex parameter optimization problem.It was expected to achieve the objective of the amplitude and equivalent phase lag of the computed magnetization being close to those calculated using the Jiles-Atherton model.Furthermore, considering the importance of amplitude in engineering applications, the calculation accuracy of the amplitude should be guaranteed first when there is a trade-off between the calculation accuracies of the amplitude and equivalent phase lag.
Moreover, the equivalent phase lag could be determined by using Equation ( 4) to directly fit the calculated results of the Jiles-Atherton model.Of course, the direct fitting method, using the least square fitting method, takes into account both the amplitude and phase lag, and it indicated that the amplitude of magnetization was not determined so accurately.
To acquire a more accurate amplitude, the values of maximum magnetization M max , the maximum magnetic field intensity H max , and the equivalent phase lag ϕ were first extracted.Based on the amplitude-frequency and phase-frequency relationships shown in Equation ( 4), the optimal parameter value can be directly computed from

√
(Lω) 2 +R 2 H max and tanϕ = ωL/R, where one obtains In short, from the three steps of (1) direct fitting to extract the equivalent phase lag, (2) extracting the maximum value, and (3) computing, the optimal values of R and L for the first-order LTI system model were determined.
Figure 7 shows the computation effects of the different parameters from direct fitting and Equation (7), labeled as "Revised fitting".From the results, it can be observed that the first-order LTI system model with the appropriate parameters effectively described the amplitude and equivalent phase lag of magnetization under a low H max .Furthermore, compared to direct fitting, revised fitting with the parameters computed from Equation (7) preserved the amplitude and equivalent phase lag more precisely.This could indicate that Equation ( 7) achieved the perfect amplitude and equivalent phase lag by sacrificing the calculation accuracy in the changing process of magnetization from 0 to maximum or minimum, which is not so important in engineering applications.
that the first−order LTI system model with the appropriate parameters effectively described the amplitude and equivalent phase lag of magnetization under a low Hmax.Furthermore, compared to direct fitting, revised fitting with the parameters computed from Equation ( 7) preserved the amplitude and equivalent phase lag more precisely.This could indicate that Equation ( 7) achieved the perfect amplitude and equivalent phase lag by sacrificing the calculation accuracy in the changing process of magnetization from 0 to maximum or minimum, which is not so important in engineering applications.

Empirical Equations for the Parameters
Equation ( 7) can predict the exact values of the amplitude and the equivalent phase lag, and the equation relies on a fitting method to predetermine the equivalent phase lag.Here, some empirical expressions are provided for easy use.

Function Format
The values of R and L computed from the equivalent phase lag and Mmax are time−independent and related to the material properties.The main concern is that they are not independent of the maximum magnetic field intensity Mmax. Figure 8 provides the optimal values of R and L at different values of Hmax.

Empirical Equations for the Parameters
Equation ( 7) can predict the exact values of the amplitude and the equivalent phase lag, and the equation relies on a fitting method to predetermine the equivalent phase lag.Here, some empirical expressions are provided for easy use.

Function Format
The values of R and L computed from the equivalent phase lag and M max are timeindependent and related to the material properties.The main concern is that they are not independent of the maximum magnetic field intensity M max .Figure 8 provides the optimal values of R and L at different values of H max .To express the variables of R and L in simple and accurate formulae, the following can be observed: (1) R or L is better expressed by a function of Hmax while not being an Hmax−independent variable.That is, it is better to predetermine the value of Hmax.To express the variables of R and L in simple and accurate formulae, the following can be observed: (1) R or L is better expressed by a function of H max while not being an H max -independent variable.That is, it is better to predetermine the value of H max .(2) To meet lim 2 for a similar amplitude-frequency function indicating that lim The equivalent phase lag ϕ decreases with an increase in H max until 0, and then, one obtains tanϕ = ωL/R.Based on (2), the expression of R must include the linear term of H max and meet lim Lω/H max = 0. (4) The power function and constant are used, while the other function forms are not introduced as far as possible for simple expressions.
Based on the curve shape and the above principles, the functions of the variables R and L were determined as follows: where A L and B L are the parameters of L, and C R , D R , and E R are the parameters of R. All parameters are independent of time, and they are only related to the material properties, such as M S and a.To meet = 0 and lim 1).
ω L was introduced to eliminate the influence of the waveform and the frequency of the magnetic field intensity H.In essence, ω L was used to remove the influence of dH/dt on the calculated magnetization.The calculation results of the static Jiles-Atherton model are only related to the amplitude of H and not to the increase in H, while the first-order LTI system model must be independent of the increase in H.To remove the influence of the frequency or waveform, the expression of 1/ω L can be added to L, and then L•ω can be formed as a parameter independent of the increase in H.The value of ω L can be easily determined from several empirical fittings.For a sinusoidal magnetic field, ω L = ω.For aperiodic signals or rising and falling signals commonly used in switching devices, 1/ω L can be determined without considering the frequency and only considering the influence of the waveform.So, one advantage of employing the first-order LTI system model to describe the magnetic hysteresis is that the amplitude and the phase of the linear steady system response can naturally reflect the influence of the frequency loss and then easily describe the magnetic loss at a high frequency.Instead, the effect of the frequency should be eliminated in the case of static hysteresis, or the hysteresis caused by the frequency should be ignored.

Parameter Determination
The five parameters A L , B L , C L , D R , and E R were only determined by the properties of the material, referring to α, a, c, and k in the Jiles-Atherton model (M S was previously considered).For a specific material, the five parameters are invariant parameters and can be easily determined by any fitting method.
Here, an easy-to-use empirical equation suitable for various materials was presented.It is difficult to determine these multivariate functions, especially when the function form is unknown and there is no clear fitting direction.The univariate fitting of a single independent variable was executed to determine several univariate functions, and then the multivariate function was constructed by the linear combination of these univariate functions.
(1) Select one parameter as the independent variable and fix the other variables as any value.Calculate the magnetization values using the Jiles-Atherton model or any other verified hysteresis model based on the different values of the specified independent variable.(2) Extract the maximum value of magnetization M max and determine the equivalent phase lag ϕ via direct fitting.(3) Calculate the optimal values of R and L under variable parameters using Equation (7) based on the obtained M max and ϕ under different values of H max .(4) Determine the optimal values of A L , B L , C L , D R , and E R using the fitting method based on Equation ( 8).
(5) Fit the functional relationship between the optimal values of A L , B L , C L , D R , and E R and the specified independent variable.The computational precision of amplitude is guaranteed with priority.(6) Change the selected variable, return to (1), and then repeat (1)-( 6) to determine the functional relationship (univariate) between the optimal values of A L , B L , C L , D R , E R , and the other parameters.(7) Employ the linear combination method to transform multiple univariate functions into a multivariate function.
An example is given in Appendix B. Other parameters are determined successively according to the above method.Finally, one obtains  By substituting Equations ( 8) and ( 9) into Equation ( 3), one obtains the first-order LTI system sub-model, which is effective under the magnetic field intensity with a low amplitude.

Verification under Fixed Values
Figure 9 shows the calculation effect of adopting the first-order LTI system model using empirical equations under the parameters of k = 2 kA/m, a = 12 kA/m, c = 0.2, α = −0.001,and M S = 80 kA/m.The maximum values of the magnetic field intensity were taken as H 0.083 = 5 kA/m, H 0.209 = 10 kA/m, H 0.327 = 15 kA/m, and H 0.429 = 20 kA/m, where the corresponding maximum magnetization values at these magnetic field intensities reached 0.083 M S , 0.209 M S , 0.327 M S , and 0.429 M S , respectively.The parameters were calculated using Equation (9).Compared to the Jiles-Atherton model, the first-order LTI system model could effectively predict the amplitude and equivalent phase lag of magnetization, which is consistent with prior designs.The prediction effect deteriorated with an increase in H max .Moreover, the deviation was mainly embodied in the change process between 0 and maximum or minimum and 0. We extracted the relative errors of the first−order LTI system model when computing the amplitude and equivalent phase lag, and they are shown in Figure 10, in which the abscissa variable is used as the r value of Hr to represent the magnetic field intensity with an amplitude from 5 kA/m to 25 kA/m.From the results of the relative errors, it was observed that the first−order LTI system model could describe the amplitude of magnetization with high precision, as the relative error was less than 2%.In contrast, the proposed first−order LTI system model failed to precisely predict the equivalent phase lag, as the relative error exceeded 15% in some cases.Combined with the absolute value We extracted the relative errors of the first-order LTI system model when computing the amplitude and equivalent phase lag, and they are shown in Figure 10, in which the abscissa variable is used as the r value of H r to represent the magnetic field intensity with an amplitude from 5 kA/m to 25 kA/m.From the results of the relative errors, it was observed that the first-order LTI system model could describe the amplitude of magnetization with high precision, as the relative error was less than 2%.In contrast, the proposed first-order LTI system model failed to precisely predict the equivalent phase lag, as the relative error exceeded 15% in some cases.Combined with the absolute value of the equivalent phase lag shown in Figure 10 and the time domain curves shown in Figure 9, this level of calculation accuracy of the equivalent phase lag is acceptable in engineering applications.In fact, the calculation deviation will be weakened when more links are added to the final output.We extracted the relative errors of the first−order LTI system model when computing the amplitude and equivalent phase lag, and they are shown in Figure 10, in which the abscissa variable is used as the r value of Hr to represent the magnetic field intensity with an amplitude from 5 kA/m to 25 kA/m.From the results of the relative errors, it was observed that the first−order LTI system model could describe the amplitude of magnetization with high precision, as the relative error was less than 2%.In contrast, the proposed first−order LTI system model failed to precisely predict the equivalent phase lag, as the relative error exceeded 15% in some cases.Combined with the absolute value of the equivalent phase lag shown in Figure 10 and the time domain curves shown in Figure 9, this level of calculation accuracy of the equivalent phase lag is acceptable in engineering applications.In fact, the calculation deviation will be weakened when more links are added to the final output.Area of hysteretic loop is an important parameter to indicate the energy loss per period.As the amplitude and equivalent phase lag can be acquired with low deviations by the proposed model, it can be concluded approximately that the loop area can be calculated effectively.Table 1 gives the calculated areas of hysteretic loop from the two models.From the calculated results, the first-order LTI system model had relative error not higher than 4.0% compared to the proposed model in calculating the loop area accurately.

Parameter Applicability
Figure 10 shows the relative errors under fixed parameters.We extracted the maximum relative errors under more parameters (univariate change), as shown in Figure 11.According to the error analysis results, the prediction effect of the magnetization amplitude is good, as the calculation errors of the amplitudes are lower than 5%.Moreover, the precision of computing the phase lag is acceptable within the given parameter range.With an acceptable upper error limit of about 5%, the empirical equations proposed in this paper for the first-order LTI system model are approximately suitable for the conditions of k ∈ tude is good, as the calculation errors of the amplitudes are lower than 5%.Mo the precision of computing the phase lag is acceptable within the given parameter With an acceptable upper error limit of about 5%, the empirical equations propo this paper for the first− The relative errors in Figure 11 were obtained by changing a single parameter, and the other parameters were fixed.To implement an effective supplement for the multivariate function expression, Figure 12 shows the maximum relative errors under three boundary conditions when computing the magnetization amplitude.From the calculation results, besides the effective parameter range given above, it was determined that the value of c should not be too low, preferably not lower than 0.12, when the value of k/a is quite low.Moreover, the value of c should not be too high, preferably not higher than 0.24, when the value of k/a is quite high.Meanwhile, the value of α had a constraint effect on the relative error, indicating that a higher value of α is conducive to roughly averaging the errors to different values of H max .The relative errors in Figure 11 were obtained by changing a single parameter, and the other parameters were fixed.To implement an effective supplement for the multivariate function expression, Figure 12 shows the maximum relative errors under three boundary conditions when computing the magnetization amplitude.From the calculation results, besides the effective parameter range given above, it was determined that the value of c should not be too low, preferably not lower than 0.12, when the value of k/a is quite low.Moreover, the value of c should not be too high, preferably not higher than 0.24, when the value of k/a is quite high.Meanwhile, the value of α had a constraint effect on the relative error, indicating that a higher value of α is conducive to roughly averaging the errors to different values of Hmax.

Nonlinear Simple Function Sub-Model under High Magnetic Field Intensity 6.1. Optimal Parameter for Amplitude
The nonlinear simple function sub-model, similar to anhysteretic magnetization M an , is suitable when the hysteresis loop can be neglected.In the expression of K an is the only parameter independent of t or H while being related to the material properties.For the exact amplitude, K an under certain materials can be easily determined from = M s K an ≈ 1 3a/M s −α .By introducing a new parameter β to reduce the deviations of the amplitudes of M and M an , one obtains an empirical expression of K an :

Parameter Determination
By setting Equation (11) to be equal to Equation (10), one obtains the optimal estimate K an is a parameter dependent on the material property, which is represented as the five parameters of k, a, c, α, and M S in the Jiles-Atherton model.As a, α, and M S are considered in the sub-part of K an , as shown in Equation (10), it is reasonable to construct β using the function mainly dependent on k and c.The estimated values of β under different parameters are extracted, and they are plotted in Figure 13.By constructing the multivariate function expression using a basic polynomial fitting of binary function, or by combining the univariate functions as described in Section 5.2.2 and then taking a simple function of α to execute a minor correction, one finally obtains the β expression: By substituting Equation (12) into Equation ( 11), one obtains Kan is a time−independent and Hmax−independent parameter.By substituting Equation (13) into  By constructing the multivariate function expression using a basic polynomial fitting of binary function, or by combining the univariate functions as described in Section 5.2.2 and then taking a simple function of α to execute a minor correction, one finally obtains the β expression: By substituting Equation (12) into Equation (11), one obtains Kan = (3a K an is a time-independent and H max -independent parameter.By substituting Equation ( 13) into M(t) = M s H(t)/ H(t) 2 + K 2 an , a simple nonlinear expression describing the M-H relationship can be achieved and is suitable for the condition of H max being high.Compared to the first-order LTI system model applied to the low-H max condition, the simple nonlinear model employs the parameter K an , independent of the magnetic field intensity, and then predicting the value of H max is no longer required.Furthermore, the influence of the frequency is not naturally considered in the nonlinear model, indicating that some frequency factors should be added to K an (or other terms) under quite a high frequency or that other conditions where the frequency has an influence should be considered.Compared to the first−order LTI system model applied to the low−Hmax condition, the simple nonlinear model employs the parameter Kan, independent of the magnetic field intensity, and then predicting the value of Hmax is no longer required.Furthermore, the influence of the frequency is not naturally considered in the nonlinear model, indicating that some frequency factors should be added to Kan (or other terms) under quite a high frequency or that other conditions where the frequency has an influence should be considered.It can be observed from the calculation results that the model can effectively predict the amplitude and shape of magnetization when the maximum value of the magnetic field intensity belongs to H0.5~H0.9, as the maximum relative errors were only 0.98%, 0.19%, 2.99%, and 4.38%.In contrast, the maximum relative errors were 9.47%, 2.58%, 0.99%, and 0.50% when Man was used for the calculation.Moreover, with the increase in Hmax, the influence of the hysteresis width becomes smaller, meaning that the computation effect of Man will be better.Overall, the nonlinear function in the form of

Parameter Applicability
Figure 15 provides the extracted maximum relative errors under more parameters (univariate change).It can be observed from the calculation results that the proposed empirical equation exhibited good parameter applicability, as the relative errors under various parameters were low.Taking 5% as the acceptable upper limit of the relative error, the simple nonlinear model is roughly applicable to the parameter range of k ∈ [1000, 4000], a ∈ [5000, 20,000], α ∈ [−0.01, 0.002], and c ∈ (0, 1).It can be observed from the calculation results that the model can effectively predict the amplitude and shape of magnetization when the maximum value of the magnetic field intensity belongs to H 0.5 ~H0.9 , as the maximum relative errors were only 0.98%, 0.19%, 2.99%, and 4.38%.In contrast, the maximum relative errors were 9.47%, 2.58%, 0.99%, and 0.50% when M an was used for the calculation.Moreover, with the increase in H max , the influence of the hysteresis width becomes smaller, meaning that the computation effect of M an will be better.Overall, the nonlinear function in the form of M(t) = M s H(t)/ H(t) 2 + K 2 an can accurately describe magnetization under a large range with quite a simple format.

Parameter Applicability
Figure 15 provides the extracted maximum relative errors under more parameters (univariate change).It can be observed from the calculation results that the proposed empirical equation exhibited good parameter applicability, as the relative errors under various parameters were low.Taking 5% as the acceptable upper limit of the relative error, the simple nonlinear model is roughly applicable to the parameter range of k ∈ [1000, 4000], a ∈ [5000, 20,000], α ∈ [−0.01, 0.002], and c ∈ (0, 1). Figure 16 shows the maximum relative errors under three boundary conditi supply the error analysis.From the calculation results, besides the effective para range given above, the proposed nonlinear model should meet k ∈ [a/7.2, a/3.5] w the parameters take extreme values (such cases are rare) and when the magneti intensity is not too large.Figure 16 shows the maximum relative errors under three boundary conditions to supply the error analysis.From the calculation results, besides the effective parameter range given above, the proposed nonlinear model should meet k ∈ [a/7.2, a/3.5] when all the parameters take extreme values (such cases are rare) and when the magnetic field intensity is not too large.Figure 16 shows the maximum relative errors under three boundary conditions to supply the error analysis.From the calculation results, besides the effective parameter range given above, the proposed nonlinear model should meet k ∈ [a/7.2, a/3.5] when all the parameters take extreme values (such cases are rare) and when the magnetic field intensity is not too large.

A Simple Summary
To describe M, Equation (6) gives the method of using the first-order LTI system under a low value of H max and a simple nonlinear function under a high value of H max .
where L and R are dependent on the material property parameters (except for M S ) and H max , and L is dependent on an increase in H at the same time.When H max ≈ H r , both sub-models are acceptable.
All the parameters can be easily determined by using the fitting method for certain materials.Another easily used determination method is given by the following empirical expressions: where ω L is employed to remove the influence of the waveform and the frequency of the magnetic field intensity.Other parameters are based on the J-A model.
It should be repeated that the empirical equations shown in Equation (15) are not necessary for parameter determination.For specific materials, these parameters are fixed as α, k, etc., in the Jiles-Atherton model.Moreover, all the parameters can be determined using any fitting method or computation after extracting some important characteristic values, as described in Sections 5.1 and 6.1.In fact, the parameters determined by the fitting method or computation after extraction can achieve better effects than the empirical equations shown in Equation ( 15), as they do not need to be responsible for other situations.

On/Off Type Device-A Special Example
For on/off type devices, the response is mainly determined by the amplitude and response time, while the change process from minimum to maximum has little influence, as its duration is quite short.Furthermore, the time lag (or phase lag) of magnetization compared to the magnetic field occupies a small proportion of the whole time lag, which is mainly caused by the electrical inductance, mechanical clearance, and damper.Under this condition, the first-order LTI system or the nonlinear function model is capable of describing magnetization more accurately.In fact, any model that can accurately describe the M max -H max relationship and that can supply a low phase lag for M will be able to effectively predict the response of this type of device.
We took the giant magnetostrictive material used on a high-speed on/off device as an example to analyze the calculation differences between the above models.Table 2 shows the Jiles-Atherton model parameters of the giant magnetostrictive material, the first-order LTI system model parameters from the empirical equation, the first-order LTI system model parameters from direct fitting, the nonlinear function model parameters from the empirical equation, and the nonlinear function model parameters from direct fitting.To guarantee a high-enough magnetic field intensity, a coil with a high ratio of inductance to resistance is often used.The relationship between the magnetic field intensity H(t), coil current I(t), and inputted voltage U(t) were given by a linear proportional equation H(t) = CNI(t)/L and a first-order linear ordinary differential equation L elec dI(t)/dt + R elec I(t) = U(t), where C, N, and L are proportional parameters of reluctance, the number of the winding turns, and the material length, respectively, and L elec and R elec are the inductance and resistance of the circuit with values of 50 mH and 5 Ω, respectively.
The coil was stimulated by a DC square wave with a high-voltage duration of π/5 s, and the response calculation results under different models are shown in Figure 17.From the calculated results, it was observed that the computation deviations in the change process were covered up, as the rise and fall times were quite short compared with the duration of the high magnetic field, similar to the above analysis.Moreover, the results verified the effectiveness of the proposed models, as the prediction effects of both models on the variation trend and amplitude of magnetization were acceptable.
ters from the empirical equation, and the nonlinear function model parameters from direct fitting.To guarantee a high−enough magnetic field intensity, a coil with a high ratio of inductance to resistance is often used.The relationship between the magnetic field intensity H(t), coil current I(t), and inputted voltage U(t) were given by a linear proportional equation H(t) = CNI(t)/L and a first−order linear ordinary differential equation LelecdI(t)/dt + RelecI(t) = U(t), where C, N, and L are proportional parameters of reluctance, the number of the winding turns, and the material length, respectively, and Lelec and Relec are the inductance and resistance of the circuit with values of 50 mH and 5 Ω, respectively.
The coil was stimulated by a DC square wave with a high−voltage duration of π/5 s, and the response calculation results under different models are shown in Figure 17.From the calculated results, it was observed that the computation deviations in the change process were covered up, as the rise and fall times were quite short compared with the duration of the high magnetic field, similar to the above analysis.Moreover, the results verified the effectiveness of the proposed models, as the prediction effects of both models on the variation trend and amplitude of magnetization were acceptable.To further analyze the response time, Figure 18 shows the proportion of the magnetization response time to the high−voltage duration.The rise time was defined by the time required for magnetization to increase from 0 (or minimum) to 90% of the maximum value.From the calculations, it was observed that the response times of all the models accounted for very low proportions, less than 6%, which had small effects on the performance of the on/off type device.In addition, the response time calculated using the first−order LTI system sub−model was always slightly higher than the one calculated using the Jiles-Atherton model for the introduction of new differential elements.In con- To further analyze the response time, Figure 18 shows the proportion of the magnetization response time to the high-voltage duration.The rise time was defined by the time required for magnetization to increase from 0 (or minimum) to 90% of the maximum value.From the calculations, it was observed that the response times of all the models accounted for very low proportions, less than 6%, which had small effects on the performance of the on/off type device.In addition, the response time calculated using the first-order LTI system sub-model was always slightly higher than the one calculated using the Jiles-Atherton model for the introduction of new differential elements.In contrast, the nonlinear function model did not introduce a new delay element, indicating that it was better in describing the response time.Figure 19 shows the relative errors of the proposed models in computing the amplitude.According to the calculation results, as long as Hmax was not less than H0.25 (this condition can always be satisfied in engineering applications), the relative errors of the two models were not more than 5%, which shows a good calculation effect of the two models.Furthermore, the first−order LTI system sub−model with direct fitting parameters had the highest calculation accuracy, and the relative errors of this method were less than 1.5%.When describing the magnetization stimulated by the magnetic field in a DC square waveform, both the first−order LTI system sub−model and the nonlinear function model "fail" to describe the residual magnetization, which dramatically reflects the advantages of the two models used in engineering modeling.Residual magnetization is not critical, as its effect is always removed (except for materials working based on their residuals); strictly speaking, it should be removed.Residual magnetization with a low amplitude fails to overcome the preload or interstice of the devices, and it cannot produce any mechanical response, which is the same as the effect of zero magnetization.Moreover, when the residual magnetization is large enough to trigger an effective response, negative magnetic field compensation is generally required to eliminate the influence of the Figure 19 shows the relative errors of the proposed models in computing the ampli-According to the calculation results, as long as H max was not less than H 0.25 (this condition can always be satisfied in engineering applications), the relative errors of the two models were not more than 5%, which shows a good calculation effect of the two models.Furthermore, the first-order LTI system sub-model with direct fitting parameters had the highest calculation accuracy, and the relative errors of this method were less than 1.5%.
trast, the nonlinear function model did not introduce a new delay element, indicating that it was better in describing the response time.Figure 19 shows the relative errors of the proposed models in computing the amplitude.According to the calculation results, as long as Hmax was not less than H0.25 (this condition can always be satisfied in engineering applications), the relative errors of the two models were not more than 5%, which shows a good calculation effect of the two models.Furthermore, the first−order LTI system sub−model with direct fitting parameters had the highest calculation accuracy, and the relative errors of this method were less than 1.5%.When describing the magnetization stimulated by the magnetic field in a DC square waveform, both the first−order LTI system sub−model and the nonlinear function model "fail" to describe the residual magnetization, which dramatically reflects the advantages of the two models used in engineering modeling.Residual magnetization is not critical, as its effect is always removed (except for materials working based on their residuals); strictly speaking, it should be removed.Residual magnetization with a low amplitude fails to overcome the preload or interstice of the devices, and it cannot produce any mechanical response, which is the same as the effect of zero magnetization.Moreover, when the residual magnetization is large enough to trigger an effective response, negative magnetic field compensation is generally required to eliminate the influence of the When describing the magnetization stimulated by the magnetic field in a DC square waveform, both the first-order LTI system sub-model and the nonlinear function model "fail" to describe the residual magnetization, which dramatically reflects the advantages of the two models used in engineering modeling.Residual magnetization is not critical, as its effect is always removed (except for materials working based on their residuals); strictly speaking, it should be removed.Residual magnetization with a low amplitude fails to overcome the preload or interstice of the devices, and it cannot produce any mechanical response, which is the same as the effect of zero magnetization.Moreover, when the residual magnetization is large enough to trigger an effective response, negative magnetic field compensation is generally required to eliminate the influence of the residual magnetization.That is, the final response of the device should be zero with any residual magnetization.To model this practical condition, the Jiles-Atherton model or other hysteresis models with a residual must be used, at least with the cooperation of a dead-zone function.On the contrary, the first-order LTI system sub-model or the nonlinear function model shows their conciseness and efficiency, as both can describe this condition without the introduction of any new sub-models.

Conclusions
The first-order LTI system model and a nonlinear function model were proposed to describe the hysteresis characteristics of a ferromagnetic material with a simple format.Moreover, the parameter determination methods and results were given based on the Jiles-Atherton model parameters.
(1) Neglecting the nonlinearity of anhysteretic magnetization, the first-order LTI system model was suitable for conditions with a magnetic field amplitude not higher than H 0.5 .Moreover, empirical expressions for parameters applicable to various materials were given based on univariate fitting and the linear combination of univariate functions.The proposed equations required the foreknowledge of the maximum and the increase in the magnetic field intensity, and they were approximately suitable for the parameter ranges k ∈ [1000, 4000], a ∈ [8000, 15,500], α ∈ [−0.01, 0.005], and c ∈ (0.1, 0.3).An error analysis showed good performance of the first-order LTI system model, as the calculation errors of the amplitudes were lower than 5%, and the precision of computing the phase lag was acceptable.(2) Neglecting the phase lag, the nonlinear function model was suitable under a magnetic field amplitude not lower than H 0.5 .The nonlinear function model had only one parameter and did not require any foreknowledge.Empirical equations for the parameters were given and suitable for the parameter ranges k ∈ [1000, 4000], a ∈ [5000, 20,000], α ∈ [−0.01, 0.002], and c ∈ (0, 1).The error analysis showed good performance of the nonlinear function model, as the calculation errors were lower than 5%.(3) Both the first-order LTI system model and the nonlinear function model can effectively predict the magnetization of the material employed in an on/off type device.Taking giant magnetostrictive material as an example, the two models with the parameters from the direct fitting or empirical equations showed acceptable calculation effects, as the relative errors when computing the amplitude were not higher than 5%, and the proportions of the responding time under all conditions were less than 6%.Furthermore, the first-order LTI system model with direct fitting parameters performed the best when computing the amplitude, with a relative error less than 1.5%, while the nonlinear function model was better in describing the response time when not introducing a time lag.
values of ϕ under the condition that A ms = 5 × 10 −6 , and Figure A1c shows various M-H curves with different values of H max or M S .The calculation results show that the first-order LTI system model can effectively predict the hysteresis loop with various amplitudes, widths, and slopes (the slope of the line connecting the farthest points and the center of the hysteresis loop).However, the first-order LTI system fails to predict the hysteresis loop with a changing slope, as the mean points of the upper and lower branches of the hysteresis loop are always summarized in a straight line.

Figure 1 .
Figure 1.Structure of this article.

Figure 3 .
Figure 3. Time−domain curves of H and M with different amplitudes and different lagging phases.

Figure 3 .
Figure 3. Time-domain curves of H and M with different amplitudes and different lagging phases.

Figure 3 .
Figure 3. Time−domain curves of H and M with different amplitudes and different lagging phases.

Figure 4 .
Figure 4. M−H curves with different amplitudes.Figure 4. M-H curves with different amplitudes.

Figure 4 .
Figure 4. M−H curves with different amplitudes.Figure 4. M-H curves with different amplitudes.

Figure 5 .Figure 6 .
Figure 5. Magnetization curves under low magnetic field intensity calculated from the Jiles-Atherton model and the 1st−order LTI system model: (a) time−domain curves; (b) M−H curves.

Figure 5 .
Figure 5. Magnetization curves under low magnetic field intensity calculated from the Jiles-Atherton model and the 1st-order LTI system model: (a) time-domain curves; (b) M-H curves.

Figure 5 .Figure 6 .
Figure 5. Magnetization curves under low magnetic field intensity calculated from the Jiles-Atherton model and the 1st−order LTI system model: (a) time−domain curves; (b) M−H curves.

Figure 6 .
Figure 6.Magnetization curves under high magnetic field intensity calculated from the Jiles-Atherton model and the 1st-order LTI system model: (a) time-domain curves; (b) M-H curves.

Figure 7 .
Figure 7. Comparisons of calculated magnetization curves from Jiles-Atherton model, the 1st LTI system model with parameters from directly fitting and the 1st LTI system model with parameters from fitting and extracting (under fixed parameters of k = 2 kA/m, a = 12 kA/m, c = 0.2, α = −0.001and Ms = 80 kA/m during different amplitudes of magnetic field intensity): (a) time−domain curves; (b) M−H curves.

Figure 7 .
Figure 7. Comparisons of calculated magnetization curves from Jiles-Atherton model, the 1st LTI system model with parameters from directly fitting and the 1st LTI system model with parameters from fitting and extracting (under fixed parameters of k = 2 kA/m, a = 12 kA/m, c = 0.2, α = −0.001and M s = 80 kA/m during different amplitudes of magnetic field intensity): (a) time-domain curves; (b) M-H curves.

Figure 8 .
Figure 8. Optimal values of R and L under fixed parameters of k = 2 kA/m, a = 12 kA/m, c = 0.2, α = −0.001and Ms = 80 kA/m during different amplitudes of magnetic field intensity (the Abscissa axis gave the values of amplitude and r value of Hr simultaneously).

4 )
constructed functions of L and R can follow (Lω) 2 + R 2 = f 2 (Hmax) + Hmax 2 for a similar amplitude-frequency function indicating that f x = Kx. (3)The equivalent phase lag φ decreases with an increase in Hmax until 0, and then, one obtains tanφ = ωL/R.Based on (2), the expression of R must include the linear term of Hmax and meet The power function and constant are used, while the other function forms are not introduced as far as possible for simple expressions.

Figure 8 .
Figure 8. Optimal values of R and L under fixed parameters of k = 2 kA/m, a = 12 kA/m, c = 0.2, α = −0.001and M s = 80 kA/m during different amplitudes of magnetic field intensity (the Abscissa axis gave the values of amplitude and r value of H r simultaneously).
80 kA/m) during various magnetic field intensities with low amplitudes: (a) time−domain curves; (b) M−H curves.

Figure 10 .
Figure 10.Relative errors of the 1st LTI system model in calculating the amplitude and equivalent phase lag under fixed parameters (k = 2 kA/m, a = 12 kA/m, c = 0.2, α = −0.001and M S = 80 kA/m): (a) relative errors; (b) absolute value of equivalent phase lag.
Figure10shows the relative errors under fixed parameters.We extracted the maximum relative errors under more parameters (univariate change), as shown in Figure11.According to the error analysis results, the prediction effect of the magnetization amplitude is good, as the calculation errors of the amplitudes are lower than 5%.Moreover, the precision of computing the phase lag is acceptable within the given parameter range.With an acceptable upper error limit of about 5%, the empirical equations proposed in this paper for the first-order LTI system model are approximately suitable for the conditions of k ∈ [1000 kA/m, 4000 kA/m], a ∈ [8000 kA/m, 15,500 kA/m], α ∈ [−0.01, 0.005], and c ∈ [0.1, 0.3].

Figure 12 .Figure 12 .
Figure 12.Relative errors of the 1st LTI system model in calculating the magnetization amplitude under boundary conditions with highest errors.6. Nonlinear Simple Function Sub−Model under High Magnetic Field Intensity 6.1.Optimal Parameter for AmplitudeThe nonlinear simple function sub−model, similar to anhysteretic magnetization Man, is suitable when the hysteresis loop can be neglected.In the expression of

Materials 2022 ,.
15,  x FOR PEER REVIEW 15 of 24 can be used to make a comparison.The derivative function of the nonlinear simple funcnew parameter β to reduce the deviations of the amplitudes of M and Man, one obtains an empirical expression of Kan:(11) to be equal to Equation(10), one obtains the optimal estimate of β as  Kan is a parameter dependent on the material property, which is represented as the five parameters of k, a, c, α, and MS in the Jiles-Atherton model.As a, α, and MS are considered in the sub−part of Kan, as shown in Equation(10), it is reasonable to construct β using the function mainly dependent on k and c.The estimated values of β under different parameters are extracted, and they are plotted in Figure13.
simple nonlinear expression describing the M−H relationship can be achieved and is suitable for the condition of Hmax being high.

Figure 14
Figure 14 shows the prediction effect of the proposed nonlinear model, including timedomain curves and M-H curves under various magnetic field intensities.At the same time, the calculation effect of using M an to approximate M is also given.The maximum values of the magnetic field intensity are H 0.427 = 20 kA/m, H 0.681 = 40 kA/m, H 0.790 = 60 kA/m, and H 0.845 = 80 kA/m.

Figure 14
Figure 14 shows the prediction effect of the proposed nonlinear model, including time−domain curves and M−H curves under various magnetic field intensities.At the same time, the calculation effect of using Man to approximate M is also given.The maximum values of the magnetic field intensity are H0.427 = 20 kA/m, H0.681 = 40 kA/m, H0.790 = 60 kA/m, and H0.845 = 80 kA/m.It can be observed from the calculation results that the model can effectively predict the amplitude and shape of magnetization when the maximum value of the magnetic field intensity belongs to H0.5~H0.9, as the maximum relative errors were only 0.98%, 0.19%, 2.99%, and 4.38%.In contrast, the maximum relative errors were 9.47%, 2.58%, 0.99%, and 0.50% when Man was used for the calculation.Moreover, with the increase in Hmax, the influence of the hysteresis width becomes smaller, meaning that the computation effect of Man will be better.Overall, the nonlinear function in the form of

Figure 16 .
Figure 16.Relative errors of the nonlinear model under boundary conditions with highest e

Figure 16 .
Figure 16.Relative errors of the nonlinear model under boundary conditions with highest errors.

Figure 16 .
Figure 16.Relative errors of the nonlinear model under boundary conditions with highest errors.

Figure 17 .
Figure 17.Time−domain responses calculated from different models (a) low magnetic field intensity, (b) medium magnetic field intensity, (c) high magnetic field intensity.

Figure 17 .
Figure 17.Time-domain responses calculated from different models (a) low magnetic field intensity, (b) medium magnetic field intensity, (c) high magnetic field intensity.

Materials 2022 ,
15,  x FOR PEER REVIEW 20 of 24 trast, the nonlinear function model did not introduce a new delay element, indicating that it was better in describing the response time.

Figure 18 .
Figure 18.Proportions of the response time calculated from different models.

Figure 19 .
Figure 19.Relative errors on computing the magnetization amplitude.

Figure 18 .
Figure 18.Proportions of the response time calculated from different models.

Figure 18 .
Figure 18.Proportions of the response time calculated from different models.

Figure 19 .
Figure 19.Relative errors on computing the magnetization amplitude.

Figure 19 .
Figure 19.Relative errors on computing the magnetization amplitude.

Figure A2 .
Figure A2.Optimal and fitted values of the parameters under univariate analysis.(a) take k as independent variable; (b) take a as independent variable; (c) take c as independent variable; (d) take α as independent variable.

Figure A2 .
Figure A2.Optimal and fitted values of the parameters under univariate analysis.(a) take k as independent variable; (b) take a as independent variable; (c) take c as independent variable; (d) take α as independent variable.

Table 1 .
Calculated areas covered by the hysteresis curves from J-A model and 1st-order LTI system model.

Table 2 .
Parameters of the giant magnetostrictive material from various models or accessing ways.

Table 2 .
Parameters of the giant magnetostrictive material from various models or accessing ways.