Systems and methods for vector hysteresis analysis

Abstract

Systems and methods are provided for analyzing magnetic hysteresis. A reversible component associated with an applied magnetic field is determined for establishing a series-distributed model. One or more hysteresis loops are generated for analyzing magnetic hysteresis of a magnetic material.

Claims

1. A processor-implemented system for analyzing magnetic hysteresis of a magnetic material, the system comprising: one or more data processors configured to: determine an initial reversible component of a magnetic field applied to a magnetic material based on physical characteristics of the magnetic material that is to be built or modified; determine a first reversible component based at least in part on the initial reversible component using a vector play model by selecting a specific vector play operator that (1) limits a magnitude of an irreversible component of the magnetic field to a maximum value to ensure that a rotational loss property is maintained in saturation and (2) maintains a correct total field for the determined first reversible component; determine a flux density; and generate one or more hysteresis loops based at least in part on the flux density for analyzing magnetic hysteresis of the magnetic material; a non-transitory computer-readable storage medium configured to store data related to the first reversible component, data related to the flux density, and data related to the hysteresis loops for performing a simulation of behavior of the magnetic material, wherein a physical version of the magnetic material is built or modified based on the simulation.

2. The system of claim 1, wherein the one or more data processors are further configured to determine the first reversible component using the specific vector play operator as follows: $h_{\mathrm{rek}}=\{\begin{array}{cc}{h}_{\mathrm{rek}0}& \mathrm{if}.Math.h_k-h_{\mathrm{rek}0}.Math.<r_k\left(h_{\mathrm{rek}}\right)\\ h_k-r_k\left(h_{\mathrm{rek}}\right).Math.\frac{h_k-h_{\mathrm{rek}0}}{.Math.h_k-h_{\mathrm{rek}0}.Math.}& \mathrm{if}.Math.h_k-h_{\mathrm{rek}0}.Math.r_k\left(h_{\mathrm{rek}}\right)\end{array}$ where h.sub.k represents the magnetic field, h.sub.rek0 represents the initial reversible component, h.sub.rek represents the first reversible component, h.sub.s represents the saturation threshold, and r.sub.k represents a radius associated with the irreversible component.

3. The system of claim 2, wherein: $r_k\left(h_{\mathrm{rek}}\right)=\{\begin{array}{cc}\frac{1}{1-k_{}}{}_k& \mathrm{if}.Math.h_{\mathrm{rek}}.Math.\left(h_s-{}_k\right)\\ \frac{1}{1-k_{}}\left(h_s-.Math.h_{\mathrm{rek}}.Math.\right)& \mathrm{if}\left(h_s-{}_k\right)<.Math.h_{\mathrm{rek}}.Math.h_s\\ 0& \mathrm{if}.Math.h_{\mathrm{rek}}.Math.>h_s\end{array};$ where .sub.k represents an intrinsic coercivity for a k-th play hysteron, and k.sub. represents a parameter associated with slopes of one or more recoil lines.

4. The system of claim 3, wherein the one or more hysteresis loops include a major hysteresis loop and one or more minor hysteresis loops, the minor hysteresis loops being formed based at least in part on the major hysteresis loop and the one or more recoil lines.

5. The system of claim 3, wherein the one or more data processors are further configured to determine a slope of a recoil line to be a linear function of a slope of a tangent line of a reversible curve at an intersection with the recoil line, the reversible curve being associated with the first reversible component.

6. The system of claim 3, wherein the slopes of the one or more recoil lines are the same if k.sub. is equal to 0.

7. The system of claim 3, wherein the slopes of the one or more recoil lines are different.

8. The system of claim 2, wherein the one or more data processors are further configured to determine the first reversible component using an iteration algorithm based at least in part on the specific vector play operator.

9. The system of claim 8, wherein the iteration algorithm comprises: 1) assuming h.sub.rek=h.sub.rek0; 2) obtaining h.sub.irk from r.sub.k and a direction of (h.sub.kh.sub.rek0); 3) calculating h=h.sub.k(h.sub.rek+h.sub.irk); 4) letting h.sub.rek=h.sub.rek+h; and 5) repeating operations 2)-4) until |h|/h.sub.s<, wherein h.sub.irk represents the irreversible component.

10. The system of claim 2, wherein the one or more data processors are further configured to determine the magnetic field based at least in part on the flux density using the specific vector play operator.

11. The system of claim 2, wherein the one or more data processors are further configured to determine an irreversible component based at least in part on the first reversible component.

12. The system of claim 11, wherein the one or more data processors are further configured to determine the irreversible component in the direction of h.sub.kh.sub.rek0 with a length of r.sub.k.

13. The system of claim 11, wherein the irreversible component is approximately equal to zero when the magnetic field exceeds a saturation threshold.

14. The system of claim 11, wherein: the magnetic field corresponds to a rotating magnetic field; and a magnetic hysteresis loss of the magnetic field is approximately equal to zero when the magnetic field exceeds the saturation threshold.

15. The system of claim 1, wherein a maximum value of the irreversible component is selected to equal an intrinsic coercivity value when the reversible component is substantially zero, and wherein the maximum value of the irreversible component is selected to decrease to zero as the reversible component approaches saturation.

16. A processor-implemented system for analyzing magnetic hysteresis of a magnetic material, the system comprising: one or more data processors configured to: determine a plurality of flux densities associated with a magnetic material based on physical characteristics of the magnetic material that is to be built or modified; determine a plurality of first magnetic fields based at least in part on the flux densities; determine a total magnetic field based on a weighted sum of the first magnetic fields, wherein the total magnetic field is determined according to a series distributed model that indicates fluctuations in the total magnetic field over time as a function of weighted (w) flux densities (b) according to: $h\left(t\right)=\underset{k=1}{\overset{n}{.Math.}}{w}_k{h}_k\left(b\right);$ generate one or more hysteresis loops based at least in part on the total magnetic field for analyzing magnetic hysteresis of the magnetic material; a non-transitory computer-readable storage medium configured to store data related to the flux densities, data related to the first magnetic fields, and data related to the total magnetic field for performing a simulation of behavior of the magnetic material, wherein a physical version of the magnetic material is built or modified based on the simulation.

17. The system of claim 16, wherein the one or more data processors are further configured to determine the first magnetic fields using an iteration algorithm based on the flux densities.

18. The system of claim 16, wherein the one or more hysteresis loops include a major hysteresis loop and one or more minor hysteresis loops, the minor hysteresis loops being formed by the major hysteresis loop and one or more recoil lines.

19. The system of claim 18, wherein the recoil lines correspond to different slopes.

20. The system of claim 18, wherein the data processors are further configured to identify one or more parameters based at least in part on an input flux density curve and an input hysteresis loop.

21. The system of claim 20, wherein the one or more parameters include: a simulation flux density curve, information associated with slopes of the recoil lines, and a plurality of weighting factors associated with the first magnetic fields.

22. A processor-implemented method for analyzing magnetic hysteresis of a magnetic material, the method comprising: determining, using one or more data processors, an initial reversible component of a magnetic field applied to a magnetic material based on physical characteristics of the magnetic material that is to be built or modified; determining, using the one or more data processors, a first reversible component based at least in part on the initial reversible component using a vector play model by selecting a specific vector play operator that (1) limits a magnitude of an irreversible component of the magnetic field to a maximum value to ensure that a rotational loss property is maintained in saturation and (2) maintains a correct total field for the determined first reversible component; determining, using the one or more data processors, a flux density; and generating, using the one or more data processors, one or more hysteresis loops based at least in part on the flux density for analyzing magnetic hysteresis of the magnetic material for performing a simulation of behavior of the magnetic material, wherein a physical version of the magnetic material is built or modified based on the simulation.

23. A processor-implemented method for analyzing magnetic hysteresis of a magnetic material, the method comprising: determining, using one or more data processors, a plurality of flux densities associated with a magnetic material based on physical characteristics of the magnetic material that is to be built or modified; determining, using the one or more data processors, a plurality of first magnetic fields based at least in part on the flux densities; determining, using the one or more data processors, a total magnetic field based on a weighted sum of the first magnetic fields, wherein the total magnetic field is determined according to a series distributed model that indicates fluctuations in the total magnetic field over time as a function of weighted (w) flux densities (b) according to: $h\left(t\right)=\underset{k=1}{\overset{n}{.Math.}}{w}_k{h}_k\left(b\right);$ and generating, using the one or more data processors, one or more hysteresis loops based at least in part on the total magnetic field for analyzing magnetic hysteresis of the magnetic material for performing a simulation of behavior of the magnetic material, wherein a physical version of the magnetic material is built or modified based on the simulation.

24. A non-transitory computer readable storage medium comprising programming instructions for analyzing magnetic hysteresis of a magnetic material, the programming instructions being configured to cause one or more data processors to execute operations comprising: determining an initial reversible component of a magnetic field applied to a magnetic material based on physical characteristics of the magnetic material that is to be built or modified; determining a first reversible component based at least in part on the initial reversible component using a vector play model by selecting a specific vector play operator that (1) limits a magnitude of an irreversible component of the magnetic field to a maximum value to ensure that a rotational loss property is maintained in saturation and (2) maintains a correct total field for the determined first reversible component; determining a flux density; and generating one or more hysteresis loops based at least in part on the flux density for analyzing magnetic hysteresis of the magnetic material for performing a simulation of behavior of the magnetic material, wherein a physical version of the magnetic material is built or modified based on the simulation.

25. A non-transitory computer readable storage medium comprising programming instructions for analyzing magnetic hysteresis of a magnetic material, the programming instructions being configured to cause one or more data processors to execute operations comprising: determining a plurality of flux densities associated with a magnetic material based on physical characteristics of the magnetic material that is to be built or modified; determining a plurality of first magnetic fields based at least in part on the flux densities; determining a total magnetic field based on a weighted sum of the first magnetic fields, wherein the total magnetic field is determined according to a series distributed model that indicates fluctuations in the total magnetic field over time as a function of weighted (w) flux densities (b) according to: $h\left(t\right)=\underset{k=1}{\overset{n}{.Math.}}{w}_k{h}_k\left(b\right);$ generating one or more hysteresis loops based at least in part on the total magnetic field for analyzing magnetic hysteresis of the magnetic material for performing a simulation of behavior of the magnetic material, wherein a physical version of the magnetic material is built or modified based on the simulation.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1 depicts an example computer-implemented environment wherein users can interact with a vector hysteresis analysis system hosted on one or more servers through a network.

(2) FIG. 2 depicts an example diagram showing a scalar play operator.

(3) FIG. 3 depicts an example diagram showing a b-h curve of a play hysteron.

(4) FIG. 4 depicts an example diagram showing a parallel-connected circuit for scalar play hysterons.

(5) FIG. 5 depicts an example diagram showing an ordinary vector play operator.

(6) FIG. 6 depicts an example diagram showing an irreversible component as a function of a reversible component.

(7) FIG. 7 depicts an example diagram showing a specific vector play operator.

(8) FIG. 8 depicts an example diagram showing recoil lines with variable slopes.

(9) FIG. 9 depicts an example diagram showing irreversible limit modified based on recoil line slopes.

(10) FIG. 10 depicts an example diagram showing a vector diagram for deriving a magnetic field from a flux density.

(11) FIG. 11 depicts an example diagram showing a series-connected circuit for play hysterons.

(12) FIG. 12 depicts an example diagram showing a comparison between simulated results derived from the vector hysteresis model and measured data.

(13) FIG. 13 depicts an example diagram showing a comparison between simulated results of rotational losses derived from the vector hysteresis model and simulated results derived from the ordinary vector play model.

(14) FIG. 14 depicts an example diagram showing a reversible nonlinear m-h curve identified from ascending and descending curves.

(15) FIG. 15 depicts an example diagram showing an input normal b-h curve, an input ascending curve and an input descending curve for parameter identification of the vector hysteresis model.

(16) FIG. 16 depicts an example diagram showing symmetric minor loops and major loops derived from the vector hysteresis model based on identified parameters.

(17) FIG. 17 depicts an example flow chart for vector hysteresis analysis of one hysteron.

(18) FIG. 18 depicts another example flow chart for vector hysteresis analysis.

(19) FIG. 19 depicts an example diagram showing a system for vector hysteresis analysis.

(20) FIG. 20 depicts an example diagram showing a computing system for vector hysteresis analysis.

DETAILED DESCRIPTION

(21) FIG. 1 depicts an example computer-implemented environment wherein users 102 can interact with a vector hysteresis analysis system 104 hosted on one or more servers 106 through a network 108. The vector hysteresis analysis system 104 can assist the users 102 to generate a vector hysteresis model for analyzing magnetic hysteresis of a magnetic material (e.g., to predict the magnetic hysteresis behavior of the magnetic material).

(22) As shown in FIG. 1, the users 102 can interact with the vector hysteresis analysis system 104 through a number of ways, such as over one or more networks 108. One or more servers 106 accessible through the network(s) 108 can host the vector hysteresis analysis system 104. The one or more servers 106 can also contain or have access to one or more data stores 110 for storing data for the vector hysteresis analysis system 104.

(23) The vector hysteresis analysis system 104 determines the vector hysteresis model that satisfy certain properties or basic rules: i) a saturation property (e.g., magnetization being limited to saturation for any magnetizing process); ii) a reduction property (e.g., the vector hysteresis model being reduced to a scalar model for a large field applied in a fixed direction if a magnetic material is originally demagnetized); iii) a rotational symmetry property (e.g., a locus of magnetization vector tip tracing out a circle for any rotating field); iv) a rotational loss property (e.g., the hysteresis loss approaching zero for large rotating fields).

(24) Specifically, the vector hysteresis analysis system 104 determines a specific vector play operator for the vector hysteresis model to satisfy the rotation loss property beyond saturation. The specific vector play operator for the vector hysteresis model determined by the vector hysteresis analysis system 104 is different from the ordinary vector play operator for the ordinary vector play model which is derived from a scalar play model.

(25) A. Scalar Play Model

(26) The output of the scalar play model is given by:

(27) $\begin{array}{cc}m\left(t\right)=\underset{k=1}{\overset{n}{.Math.}}{f}_k\left(h_{\mathrm{rek}}\left(t\right)\right)& \left(1\right)\end{array}$

(28) where .sub.k(.Math.)'s correspond to one or more anhysteretic nonlinear functions. The scalar model implements a scalar play operator h.sub.rek given by:
h.sub.rek(t)=P.sub.k[h.sub.k(t)](2)

(29) The scalar play operator h.sub.rek is shown in FIG. 2. The input of the scalar play model includes an applied field intensity h.sub.k given by:
h.sub.k(t)=h(t) k=1,2, . . . ,n(3)

(30) One or more returned values of the scalar play operator are bounded within two parallel limit lines (e.g., an ascending limit line and a descending limit line), as shown in FIG. 2. A flux density b(t) is obtained as follows:

(31) $\begin{array}{cc}b\left(t\right)={}_0\left(m\left(t\right)+h\left(t\right)\right)=\underset{k=1}{\overset{n}{.Math.}}{b}_k\left(t\right)& \left(4\right)\end{array}$
where
b.sub.k(t)=.sub.0(.sub.k(h.sub.rek(t))+h(t))=b.sub.rek(t)+.sub.0(h.sub.k(t)h.sub.rek(t))(5)
with
b.sub.rek(t)=.sub.0(.sub.k(h.sub.rek(t))+h.sub.rek(t))(6)

(32) The flux density b.sub.k(t) of a play hysteron is shown in FIG. 3. When the applied field h.sub.k(t) increases, b.sub.k(t) tracks the ascending curve. If h.sub.k(t) decreases, b.sub.k(t) tracks a recoil line first, then tracks the descending curve. All recoil lines have a same slope of .sub.0. The output b(t) can be seen as the total contribution of n play hysterons. If each hysteron is represented by a circuit element with h denoting a cross quantity and b a through quantity, Equation (3) may be expressed by a parallel-connected circuit, as shown in FIG. 4. If the applied field h.sub.k is decomposed into two componentsa reversible component and an irreversible component, then the scalar play operator h.sub.rek represents the reversible component, and h.sub.kh.sub.rek stands for the irreversible component h.sub.irk which satisfies:
|h.sub.kh.sub.rek|.sub.k(7)
where .sub.k represents an intrinsic coercivity for a k-th play hysteron.

(33) The parameters of the scalar play model to be identified include one or more anhysteretic curves .sub.k(.Math.), k=1, 2, . . . , n, which are identified by letting the curves derived from the scalar play model best match the measured ones. In order to identify all these anhysteretic curves, a prohibitive experimental effort is often needed, and the identification process may also be quite complicated.

(34) B. Ordinary Vector Play Operator

(35) In some embodiments, the scalar play operator h.sub.rek as shown in FIG. 2 is expressed as:
h.sub.rek=max(min(h.sub.rek0,h.sub.k+.sub.k),h.sub.k.sub.k)(8)
or

(36) $\begin{array}{cc}{h}_{\mathrm{rek}}=h_k-\frac{{}_k\left(h_k-h_{\mathrm{rek}0}\right)}{\max \left({}_k,.Math.h_k-h_{\mathrm{rek}0}.Math.\right)}& \left(9\right)\end{array}$
where h.sub.rek0 represents an initial value of h.sub.rek.

(37) The scalar play operator h.sub.rek expressed in Equation (9) is extended for a vector field as:

(38) $\begin{array}{cc}{h}_{\mathrm{rek}}=h_k-\frac{{}_k\left(h_k-h_{\mathrm{rek}0}\right)}{\max \left({}_k,.Math.h_k-h_{\mathrm{rek}0}.Math.\right)}& \left(10\right)\end{array}$
or

(39) $\begin{array}{cc}{h}_{\mathrm{rek}}=\{\begin{array}{cc}{h}_{\mathrm{rek}0}& \mathrm{if}.Math.h_k-h_{\mathrm{rek}0}.Math.<{}_k\\ h_k-{}_k.Math.\frac{h_k-h_{\mathrm{rek}0}}{.Math.h_k-h_{\mathrm{rek}0}.Math.}& \mathrm{if}.Math.h_k-h_{\mathrm{rek}0}.Math.{}_k\end{array}& \left(11\right)\end{array}$

(40) The ordinary vector play operator h.sub.rek is illustrated by a vector diagram as shown in FIG. 5. A circle can be drawn at the tip of the vector h.sub.rek0 with a radius r.sub.k=.sub.k. If the tip of the applied field h.sub.k falls inside the circle, a reversible component h.sub.rek is kept unchanged (e.g., equal to an initial vector h.sub.rek0), as shown in FIG. 5(A). On the other hand, if the tip of the applied field h.sub.k falls outside the circle, an irreversible component h.sub.irk is set in the direction of h.sub.kh.sub.rek0 with a length of r.sub.k, and then the reversible component h.sub.rek is set as h.sub.rek=h.sub.kh.sub.irk, as shown in FIG. 5(B). Parameters for the ordinary vector play model associated with FIG. 5 are the same as those for the scalar play model described above, and therefore can be identified based on the scalar play model.

(41) As shown in FIG. 5, if the applied field h.sub.k rotates, at a steady state, the irreversible component h.sub.irk is perpendicular to the reversible component h.sub.rek. As an example, in the ordinary vector play operator as shown in FIG. 5, the magnitude of h.sub.irk is constant no matter how large the applied field h.sub.k is, which means m.sub.k, in the same direction of h.sub.rek, lags h.sub.k a certain angle. The ordinary vector play model described above may not satisfy the rotational loss property that the magnetic hysteresis loss of any applied rotating magnetic field approaches zero when the magnitude of the applied rotating magnetic field becomes saturated.

(42) C. Specific Vector Play Operator

(43) In order to satisfy the rotational loss property for any vector play hysteron, the irreversible field component may need to be zero when the applied field is beyond a saturation field h.sub.s. The vector hysteresis analysis system 104 determines a specific vector play operator h.sub.rek as:

(44) $\begin{array}{cc}{h}_{\mathrm{rek}}=\{\begin{array}{cc}{h}_{\mathrm{rek}0}& \mathrm{if}.Math.h_k-h_{\mathrm{rek}0}.Math.<r_k\left(h_{\mathrm{rek}}\right)\\ h_k-r_k\left(h_{\mathrm{rek}}\right).Math.\frac{h_k-h_{\mathrm{rek}0}}{.Math.h_k-h_{\mathrm{rek}0}.Math.}& \mathrm{if}.Math.h_k-h_{\mathrm{rek}0}.Math.r_k\left(h_{\mathrm{rek}}\right)\end{array}& \left(12\right)\end{array}$
where

(45) $\begin{array}{cc}{r}_k\left(h_{\mathrm{rek}}\right)=\{\begin{array}{cc}{}_k& \mathrm{if}.Math.h_{\mathrm{rek}}.Math.\left(h_s-{}_k\right)\\ h_s-.Math.h_{\mathrm{rek}}.Math.& \mathrm{if}\left(h_s-{}_k\right)<.Math.h_{\mathrm{rek}}.Math.h_s\\ 0& \mathrm{if}.Math.h_{\mathrm{rek}}.Math.>h_s\end{array}& \left(13\right)\end{array}$
r.sub.k(h.sub.rek) is shown in FIG. 6. The specific vector play operator h.sub.rek (e.g., considering the saturation behavior) is shown in FIG. 7.
D. Variable Slope for Recoil Lines

(46) Referring back to FIG. 3, for the scalar play operator (or the ordinary vector play operator), all recoil lines inside the hysteresis loop have the same slope of .sub.0, which may cause the area of a minor hysteresis loop bounded within a major hysteresis loop and two parallel recoil lines to be larger than a measured one. In order to more accurately represent individual minor loop behavior, the vector hysteresis analysis system 104 introduces a variable slope algorithm so that the area of a derived minor loop can best match the area of a measured minor loop.

(47) Particularly, the hysteresis analysis system 104 constructs a slope of a recoil line as a linear function of a slope of a tangent line of a reversible curve at an intersection with the recoil line, as shown in FIG. 8. For example, assuming that a recoil line intersects the reversible curve at h.sub.rek, the slope of the tangent line is given by:

(48) $\begin{array}{cc}\stackrel{~}{}=\frac{{\mathrm{df}}_k\left(h\right)}{\mathrm{dh}}{|}_{h=h_{\mathrm{rek}}}={\stackrel{~}{}}_r{}_0& \left(14\right)\end{array}$
The slope of the recoil line is given by:
=k.sub.({tilde over ()}.sub.0)+.sub.0=(k.sub.({tilde over ()}.sub.r1)+1).sub.0(15)
where k.sub. represents a parameter to be identified. When k.sub.=0, the slope of all recoil lines becomes constant, the vector hysteresis model determined by the hysteresis analysis system 104 degrades to the ordinary vector play model associated with FIG. 3. Accordingly, r.sub.k is modified as:

(49) $\begin{array}{cc}{r}_k\left(h_{\mathrm{rek}}\right)=\{\begin{array}{cc}\frac{1}{1-k_{}}{}_k& \mathrm{if}.Math.h_{\mathrm{rek}}.Math.\left(h_s-{}_k\right)\\ \frac{1}{1-k_{}}\left(h_s-.Math.h_{\mathrm{rek}}.Math.\right)& \mathrm{if}\left(h_s-{}_k\right)<.Math.h_{\mathrm{rek}}.Math.h_s\\ 0& \mathrm{if}.Math.h_{\mathrm{rek}}.Math.>h_s\end{array}& \left(16\right)\end{array}$
The modified r.sub.k(h.sub.rek) is shown in FIG. 9.

(50) After h.sub.rek is determined by the specific vector play operator (e.g., as determined by Equation (12)), a reversible flux density b.sub.rek is computed from Equation (6) with the direction of h.sub.rek. The flux density is given by:
b.sub.k(t)=b.sub.rek(t)+(h.sub.rek)(h.sub.k(t)h.sub.rek(t))(17)
E. Local Iteration Algorithm

(51) The vector hysteresis analysis system 104 introduces a local iteration algorithm, which allows either deriving the applied magnetic field h.sub.k from the flux density b.sub.k or deriving the flux density b.sub.k from the applied magnetic field h.sub.k, for numerical stability that is considered as one of the most challenging issues for practical applications. For example, the vector hysteresis analysis system 104 introduces the local iteration algorithm to efficiently locate an operating point on hysteresis loops in addition to a global iteration algorithm (e.g., a Newton-Raphson global nonlinear iteration algorithm). Specifically, Equation (12) is solved to derive h.sub.rek from h.sub.k. Since r.sub.k in Equation (16) depends on h.sub.rek, a local iterating process is performed.

(52) In some embodiments, the vector hysteresis analysis system 104 derives the flux density b.sub.k from the applied magnetic field h.sub.k. When the applied field h.sub.k locates inside the circle as shown in FIG. 5(A), h.sub.rek is kept the same as the initial value h.sub.rek0. On the other hand, when the applied field h.sub.k locates outside the circle as shown in FIG. 5(B), h.sub.rek is derived according to the following iteration process:

(53) 1) Assume h.sub.rek=h.sub.rek0;

(54) 2) Obtain h.sub.irk from r.sub.k by Equation (16) and the direction of (h.sub.kh.sub.rek0);

(55) 3) Calculate h=h.sub.k(h.sub.rek+h.sub.irk);

(56) 4) Let h.sub.rek=h.sub.rek+h;

(57) 5) Repeat operations 2)-4) until |h|/h.sub.s<.

(58) In the above iterating process, represents a relaxation factor which can be optimized based on the historic iterating results, and represents a given tolerance. After h.sub.rek is obtained, b.sub.k is computed according to Equation (17) and Equation (6).

(59) In certain embodiments, the vector hysteresis analysis system 104 derives the applied magnetic field h.sub.k from the flux density b.sub.k. That is, Equation (12) is solved inversely. A vector b.sub.rek0 and are obtained from h.sub.rek0, and a circle is drawn at the tip of the vector b.sub.rek0 with a radius of r.sub.b=r.sub.k(h.sub.rek0), as shown in FIG. 10. When the tip of the vector b.sub.k locates inside the circle, h.sub.rek is kept the same as the initial value h.sub.rek0. On the other hand, when the tip of the vector b.sub.k locates outside the circle as shown in FIG. 10, another iteration process is performed:

(60) 1) Assume b.sub.rek=b.sub.k;

(61) 2) Obtain h.sub.rek from b.sub.rek based on a reversible b-h curve according to Equation (6);

(62) 3) Obtain from h.sub.rek based on Equation (15);

(63) 4) Obtain h.sub.irk from r.sub.k by Equation (16) and the direction of (h.sub.rekh.sub.rek0);

(64) 5) Let b.sub.irk=h.sub.irk;

(65) 6) Calculate b=b.sub.k(b.sub.rek+b.sub.irk);

(66) 7) Let b.sub.rek=b.sub.rek+b;

(67) 8) Repeat operations 2)-7) until |A|/b.sub.s<, where b.sub.s represents the flux density at h.sub.s.

(68) After h.sub.rek is obtained, h.sub.k is computed as follows:
h.sub.k(b.sub.k)=h.sub.rek+(b.sub.kb.sub.rek)/(18)
F. Series-Distributed Model

(69) Furthermore, the vector hysteresis analysis system 104 generates a series-distributed hysteron model for analyzing magnetic materials (e.g., predicting magnetization behavior of the magnetic materials). Particularly, to satisfy the rotational loss property, the b-h hysteresis loop associated with the specific play operator is discontinued at h=h.sub.s for each play hysteron, as shown in FIG. 7. The series-distributed hysteron model based on series-distributed play hysterons accommodates the specific vector play operator and the variable slope algorithm described above so as to reduce any possible discontinuity. For example, the series-distributed play model is used according to the variable slope algorithm for recoil lines and the local iteration algorithm for deriving h.sub.k from b.sub.k for each play hysteron.

(70) In Equation (18), the applied field h.sub.k is expressed as a function of the flux density b.sub.k for the k-th play hysteron. For series connection, b.sub.1=b.sub.2= . . . =b.sub.n=b, and thus, the total applied field for the series-distributed play hysterons can be expressed as follows:

(71) 0$\begin{array}{cc}h\left(t\right)=\underset{k=1}{\overset{n}{.Math.}}{w}_k{h}_k\left(b\right)& \left(19\right)\end{array}$
where the parameters w.sub.k represent weighting factors for all play hysterons. The series-connected circuit is shown in FIG. 11. FIG. 12 depicts an example diagram showing a comparison between simulated results derived from the vector hysteresis model and measured data. FIG. 13 depicts an example diagram showing a comparison between simulated results of rotational losses derived from the vector hysteresis model and simulated results derived from the ordinary vector play model.
G. Parameter Identification

(72) With the introduction of the series-distributed model, all play hysterons can have a same reversible nonlinear b-h curve. As a result, a parameter identification process, as an integral function of the vector hysteresis analysis system 104, may be simplified. In some embodiments, the vector hysteresis analysis system 104 performs parameter identification for the series-distributed hysteron model (e.g., based on available measured data). For example, the measured data are directly available from the manufacture of the magnetic materials. In some embodiments, all play hysterons associated with the series-distributed model can have a same reversible nonlinear b-h curve so that the parameter identification process is greatly simplified.

(73) The parameters of the vector hysteresis model include: i) a reversible nonlinear b-h curve for all play hysterons; ii) a coefficient k.sub. for variable slope of recoil lines; iii) a weighting factor for each play hysteron. These parameters can be identified based on a normal b-h curve and a major hysteresis loop.

(74) The reversible nonlinear b-h curve is obtained from a reversible nonlinear m-h curve which can in turn be derived from a center line of an m-h major hysteresis loop. The major hysteresis loop includes an ascending curve m.sub.asd(h) and a descending curve m.sub.dsc(h). The ascending curve, or the descending curve can be directly obtained from each other based on the odd symmetry condition. Only one curve is needed from input, according to some embodiments.

(75) In certain embodiments, the inverse functions of m.sub.asd(h) and m.sub.dsc(h) are denoted as h.sub.asd(m) and h.sub.dsc(m), respectively, as shown in FIG. 14. The reversible nonlinear m-h curve is determined from:
h.sub.rev(m)=(h.sub.asd(m)+h.sub.dsc(m))/2(20)

(76) If the inverse function of h.sub.rev(m) is expressed as m.sub.rev(h), then the reversible nonlinear b-h curve can be obtained from:
b.sub.rev(h)=.sub.0(m.sub.rev(h)+h)(21)

(77) After the reversible nonlinear b-h curve is derived, the remaining parameters can be derived by making a derived (e.g., simulated) normal b-h curve and a derived (e.g., simulated) ascending curve of the major hysteresis loop best match to the input normal b-h curve and the input ascending curve of the major hysteresis loop respectively.

(78) If h.sub.rek0 starts from 0 and b sweeps from 0 to b.sub.s, the derived normal b-h curve is derived based on the local iteration algorithm for deriving h.sub.k from b.sub.k with b.sub.k=b for all hysterons. If a sweep index is denoted as i, for a given value of k.sub., h.sub.ki is obtained from b.sub.i, and a total field h.sub.i of all play hysterons is determined as follows:

(79) $\begin{array}{cc}{h}_i=\underset{k=1}{\overset{n}{.Math.}}{w}_k{h}_{\mathrm{ki}}& \left(22\right)\end{array}$

(80) If a field on the input normal b-h curve at b.sub.i is denoted as h.sub.i, then a total error between the derived b-h curve and the input normal b-h curve for all sweepings is determined as follows:

(81) $\begin{array}{cc}{e}_{\mathrm{nor}}=\underset{i=1}{\overset{m_1}{.Math.}}{\left(\underset{k=1}{\overset{n}{.Math.}}{w}_k{h}_{\mathrm{ki}}-h_i^{}\right)}^2& \left(23\right)\end{array}$

(82) where m.sub.1 represents the number of sweeps for the normal b-h curve.

(83) Similarly, if h.sub.rek0 starts from h.sub.s and b sweeps from b.sub.s to b.sub.s, the field on the derived ascending curve is obtained from the vector hysteresis model. If an index for the ascending curve sweep is counted from m.sub.1+1 to m=m.sub.1+m.sub.2, the error between the derived ascending curve and the input ascending curve for all sweepings is determined as follows:

(84) $\begin{array}{cc}{e}_{\mathrm{asd}}=\underset{i=m_1+1}{\overset{m}{.Math.}}{\left(\underset{k=1}{\overset{n}{.Math.}}{w}_k{h}_{\mathrm{ki}}-h_i^{}\right)}^2& \left(24\right)\end{array}$
where m.sub.2 represents the number of sweeps for the ascending curve.

(85) To minimize the total error, that is, let

(86) $\begin{array}{cc}\begin{array}{cc}\frac{}{w_j}\left(e_{\mathrm{nor}}+e_{\mathrm{asd}}\right)=0& j=1,2,\mathrm{.Math.},n\end{array}& \left(25\right)\end{array}$

(87) The following equation is obtained:

(88) $\begin{array}{cc}\begin{array}{cc}\underset{i=1}{\overset{m}{.Math.}}\left(\underset{k=1}{\overset{n}{.Math.}}{w}_k{h}_{\mathrm{ki}}-h_i^{}\right)h_{\mathrm{ji}}=0& j=1,2,\mathrm{.Math.},n\end{array}& \left(26\right)\end{array}$

(89) Or the following equation is obtained:

(90) $\begin{array}{cc}\left[\begin{array}{cccc}{a}_{11}& a_{12}& \mathrm{.Math.}& a_{1n}\\ a_{21}& a_{22}& \mathrm{.Math.}& a_{2n}\\ \mathrm{.Math.}& \mathrm{.Math.}& a_{\mathrm{kj}}& \mathrm{.Math.}\\ a_{n1}& n_{n2}& \mathrm{.Math.}& a_{\mathrm{nn}}\end{array}\right].Math.\left[\begin{array}{c}{w}_1\\ w_2\\ \mathrm{.Math.}\\ w_n\end{array}\right]=\left[\begin{array}{c}{y}_1\\ y_2\\ \mathrm{.Math.}\\ y_n\end{array}\right]& \left(27\right)\end{array}$
where

(91) $\begin{array}{cc}\begin{array}{cc}{a}_{\mathrm{kj}}=\underset{i=1}{\overset{m}{.Math.}}\left(h_{\mathrm{ki}}{h}_{\mathrm{ji}}\right)& k=1,2,\mathrm{.Math.},n;j=1,2,\mathrm{.Math.},n\end{array}& \left(28\right)\end{array}$
and

(92) $\begin{array}{cc}\begin{array}{cc}{y}_j=\underset{i=1}{\overset{m}{.Math.}}\left(h_i^{}{h}_{\mathrm{ji}}\right)& j=1,2,\mathrm{.Math.},n\end{array}& \left(29\right)\end{array}$

(93) After w.sub.k is solved from Equation (27), the total error, as a function of k.sub., is obtained from:

(94) $\begin{array}{cc}\mathrm{err}\left(k_{}\right)=\underset{i=1}{\overset{m}{.Math.}}{\left(\underset{k=1}{\overset{n}{.Math.}}{w}_k{h}_{\mathrm{ki}}-h_i^{}\right)}^2& \left(30\right)\end{array}$

(95) For example, an optimal k.sub. can be obtained by minimizing the total error using a one-variable numerical optimal process. In some embodiments, the vector hysteresis analysis system 104 combines a linear regression with a one-dimensional numerical optimization to minimize the total error to simplify the parameter identification of the vector hysteresis model.

(96) FIG. 15 depicts an example diagram showing an input normal b-h curve, an input ascending curve and an input descending curve for parameter identification of the vector hysteresis model. FIG. 16 depicts an example diagram showing symmetric minor loops and major loops derived from the vector hysteresis model based on identified parameters.

(97) FIG. 17 depicts an example flow chart for vector hysteresis analysis of one hysteron. At 1702, a magnetic field (e.g., h.sub.k) applied to a magnetic material is received. The magnetic field includes an initial reversible component (e.g., h.sub.rek0). At 1704, a first reversible component (e.g., h.sub.rek) is determined based at least in part on the initial reversible component (e.g., h.sub.rek0). At 1706, an irreversible component (e.g., h.sub.irk) is determined based at least in part on the first reversible component (e.g., h.sub.rek).

(98) For example, a circle can be drawn at the tip of the initial reversible component (e.g., h.sub.rek0) with a radius r.sub.k (e.g., r.sub.k=.sub.k). If the tip of the magnetic field (e.g., h.sub.k) falls inside the circle, the first reversible component (e.g., h.sub.rek) is kept unchanged (e.g., equal to the initial vector h.sub.rek0). The irreversible component (e.g., h.sub.irk) is determined to be h.sub.kh.sub.rek. On the other hand, if the tip of the magnetic field (e.g., h.sub.k) falls outside the circle, the irreversible component h.sub.irk is set in the direction of h.sub.kh.sub.rek0 with a length of r.sub.k, and then the first reversible component h.sub.rek is set as h.sub.kh.sub.irk.

(99) As an example, the irreversible component is approximately equal to a predetermined value (e.g., zero) when the applied magnetic field exceeds a saturation threshold (e.g., a saturation field h.sub.s). At 1708, a flux density is determined. At 1710, one or more hysteresis loops are generated based at least in part on the flux density for analyzing magnetic hysteresis of the magnetic material.

(100) FIG. 18 depicts another example flow chart for vector hysteresis analysis. At 1802, a plurality of flux densities associated with a magnetic material are determined. At 1804, a plurality of first magnetic fields are determined based at least in part on the flux densities. At 1806, a total magnetic field is determined based on a weighted sum of the first magnetic fields. At 1808, one or more hysteresis loops are based at least in part on the total magnetic field for analyzing magnetic hysteresis of the magnetic material.

(101) FIG. 19 depicts an example diagram showing a system for vector hysteresis analysis. As shown in FIG. 19, the system 10 includes a computing system 12 which contains a processor 14, a storage device 16 and a vector hysteresis analysis module 18. The computing system 12 includes any suitable type of computing device (e.g., a server, a desktop, a laptop, a tablet, a mobile phone, etc.) that includes the processor 14 or provide access to a processor via a network or as part of a cloud based application. The vector hysteresis analysis module 18 includes tasks and is implemented as part of a user interface module (not shown in FIG. 19).

(102) FIG. 20 depicts an example diagram showing a computing system for vector hysteresis analysis. As shown in FIG. 20, the computing system 12 includes a processor 14, memory devices 1902 and 1904, one or more input/output devices 1906, one or more networking components 1908, and a system bus 1910. In some embodiments, the computing system 12 includes the vector hysteresis analysis module 18, and provides access to the vector hysteresis analysis module 18 to a user as a stand-alone computer.

(103) This written description uses examples to disclose the invention, including the best mode, and also to enable a person skilled in the art to make and use the invention. The patentable scope of the invention may include other examples. For example, the systems and methods disclosed herein are configured to improve computational efficiency by using less play hysterons in a vector hysteresis model.

(104) For example, the systems and methods may include data signals conveyed via networks (e.g., local area network, wide area network, internet, combinations thereof, etc.), fiber optic medium, carrier waves, wireless networks, etc. for communication with one or more data processing devices. The data signals can carry any or all of the data disclosed herein that is provided to or from a device.

(105) Additionally, the methods and systems described herein may be implemented on many different types of processing devices by program code comprising program instructions that are executable by the device processing subsystem. The software program instructions may include source code, object code, machine code, or any other stored data that is operable to cause a processing system to perform the methods and operations described herein. Other implementations may also be used, however, such as firmware or even appropriately designed hardware configured to carry out the methods and systems described herein.

(106) The systems' and methods' data (e.g., associations, mappings, data input, data output, intermediate data results, final data results, etc.) may be stored and implemented in one or more different types of non-transitory computer-readable storage medium that is stored at a single location or distributed across multiple locations. The medium can include computer-implemented data stores, such as different types of storage devices and programming constructs (e.g., RAM, ROM, Flash memory, flat files, databases, programming data structures, programming variables, IF-THEN (or similar type) statement constructs, etc.). It is noted that data structures describe formats for use in organizing and storing data in databases, programs, memory, or other computer-readable media for use by a computer program.

(107) The systems and methods may be provided on many different types of computer-readable media including computer storage mechanisms (e.g., CD-ROM, diskette, RAM, flash memory, computer's hard drive, etc.) that contain instructions (e.g., software) for use in execution by a processor to perform the methods' operations and implement the systems described herein.

(108) The computer components, software modules, functions, data stores and data structures described herein may be connected directly or indirectly to each other in order to allow the flow of data needed for their operations. It is also noted that a module or processor includes but is not limited to a unit of code that performs a software operation, and can be implemented for example, as a subroutine unit of code, or as a software function unit of code, or as an object (as in an object-oriented paradigm), or as an applet, or in a computer script language, or as another type of computer code. The software components and/or functionality may be located on a single computer or distributed across multiple computers depending upon the situation at hand.

(109) It should be understood that as used in the description herein and throughout the claims that follow, the meaning of a, an, and the includes plural reference unless the context clearly dictates otherwise. Also, as used in the description herein and throughout the claims that follow, the meaning of in includes in and on unless the context clearly dictates otherwise. Finally, as used in the description herein and throughout the claims that follow, the meanings of and and or include both the conjunctive and disjunctive and may be used interchangeably unless the context expressly dictates otherwise; the phrase exclusive or may be used to indicate situation where only the disjunctive meaning may apply.