A Physical Core-Loss Model for Laminated Magnetic Sheet Steels ()
1. Introduction
The models for the power loss of magnetic materials in the changing magnetic field have been developed for decades from the earliest Steinmetz equation [1].
(1)
denotes the average of core loss,
are coefficients. Steinmetz equation only locates on frequency domain to present limited accuracy for most of magnetic materials unless ferrites. Based on frequency-domain Steinmetz equation, a series of improved time-domain equations [2] [3] [4] [5] are developed but still have not good accuracy for electrical steel sheets.
(2)
Time-domain equations can simulate the complex power loss in non-sinusoidal or non-periodic alternating magnetic field, including high-frequency trapezoidal transformers, PWM drivers, variable-speed motors and so on.
[8] [9] [10] give some accuracy solution for limited ranges of applications, and the necessary integral terms increase the calculation process and deviate from the instantaneous accuracy. [6] [7] use multiple of indexes and coefficients to improve the weakness and provide high accuracy for very wide range of frequencies and magnetic densities. Following the researches, a new equation is deduced from the physical magnetic induction equation to avoid the integral calculation and eliminate complex coefficients. This physical core-loss equation even demonstrates each parts of eddy loss, hysteresis loss and excess loss. Considering a simple electrical-magnetic induction model as the following core and coils:
Figure 1(a) shows a simple magnetic inductance system with a core and coils, the magnetic flux locates inside of the core and electrical current locates on the coil. Figure 1(b) of the section view shows the eddy currents of the laminated structure of the core. If the laminated layers are uniform, then a induction equation can be expressed as
(a)(b)
Figure 1. (a) The magnetic induction system of a laminated core and turning coils; (b) The section view with induction eddy currents.
(3)
is a constant coefficient.
The relationship between eddy current and magnetic flux can be expressed as an approximate equation of [11],
(4)
is a constant coefficient. The independent relation between the coil current and the magnetic flux can be obtained by substituting Equation (4) into Equation (3)
(5)
Therefore, the input power per unit volume is
(6)
where
(7)
2. Physical Core-Loss Equations
The primary part of Equation (7) denotes the inner power of magnetic flux and the secondary part is usually ignored. Actually the permeability should be dependent by time in real cases, the increased permeability presents a loss power in the equation by an irreversible transformation of heat. This phenomenon is similar to the resistivity of the core which produces the eddy-current loss, here the called permeability loss is one kind of core loss and consists of hysteresis loss and excess loss as usual definitions. To integrate the permeability loss and eddy loss, the core loss shows as Equation (8).
(8)
Under a periodic magnetic field, It can be reasonably predicted that the variation of
should be periodic and synchronized with the magnetic flux density. The core losses for most of magnetic materials are positive, therefore the increasing of
will happens on the ranges of large magnetic flux density, and the decreasing of
will happens on the ranges of magnetic flux density close to zero. By the above hypothesis, a simplified permeability loss equation is suggested as following:
(9)
is a constant coefficient of the simplified permeability loss equation,
,
are constant index coefficients.
However, the hypothesis simplification of Equation (9) will not approach accurately enough to the real measuring data by curve fitting, especially on the range of lower frequencies and lower magnetic densities. The phenomenon is similar to hysteresis loss, which exists a saturation limitation and can be approached accurately by the following frequency-domain equation:
(10)
is the constant coefficients.
presents a magnetic point for the saturation of hysteresis loss. The frequency-domain loss energy can be transferred into the time-domain power by differential:
(11)
is the constant coefficients of hysteresis loss equation.
Finally, the tolerance of the theoretic eddy-loss equation can be improved as adding extra one parameter of an exponential magnetic density term to compensate the nonlinear effects:
(12)
is the constant coefficients of eddy loss equation.
are constant index coefficients.
The core loss equation becomes to three parts including eddy loss, permeability loss and hysteresis loss,
(13)
The simplified permeability loss
replaces the excess loss term of the traditional core loss equation with the physical and mathematic interpretations.
3. Simulation Results
65CS400 of China Steel Corporation is a high-grade magnetic steel sheet, composed of Si-Fe-Al alloy and 0.65 mm thickness with good efficiency and permeability properties to be applied at high efficient motors and transformers. Compared to normal grade productions, 65CS400 provides more complete loss data measured in laminated structure. Table 1 shows the measured core loss by Epstein test in the range of frequency 50 - 5000 Hz and magnetic flux density of 0.1 - 1.8 T.
Table 1. The core-loss data of 65CS400 measured by Epstein test.
There are 8 coefficient variables needed to be converged in Equation (13). The convergence error
is calculated by using the root mean square error method,
(14)
is the calculated value of core-loss equation,
is the measured core- loss data of 65CS400. The converging process is to adjust the 8 coefficient variables step by step to make error
minimum. Here sheet thickness d is
0.65 mm
, the resistivity
of 65CS400 is assumed as 55e−8 Ω∙m, the density is
7800 kg
/m3. The data of 50, 100, 200, 400, 1000, 2500 and 5000 Hz is applied in the converging process, the final optimum coefficients are converged to ceddy = 0.84073, r1 = 1.62, r2 = −0.029, cperm = 0.198, n1 = 0.93, n2 = 2.15, chyst = 55.5, B0 = 0.58, m = 2.11, and the average error
is 3.38%.
Figure 2 shows the simulation results of physical core-loss equation in frequency domain, the line curves preset the periodic average core losses of 50, 60, 100, 200, 400, 1000, 2500 and 5000 Hz from left to right, and the circles present the measured core losses.
Figure 3 shows the each parts of total core loss, hysteresis loss and eddy loss in frequency 50, 400 and 5000 Hz from left to right, the dash lines present the hysteresis losses for each frequency, the dot lines present the eddy losses for each frequency.
Figure 4 presents the time-series core loss of 1000Hz-1.0T for one period, the solid line shows the total core loss, the dot line shows the eddy loss, the dash line shows the hysteresis loss and the fine line shows the permeability loss, the maximum loss exists at the peak area of the magnetic density.
Figure 2. Simulation of physical core-loss equation in frequency domain, the line curves preset the periodic average core losses of 50, 60, 100, 200, 400, 1000, 2500 and 5000 Hz from left to right, and the circles present the measured core losses.
Figure 3. The total core loss, hysteresis loss and eddy loss in frequency 50, 400 and 5000 Hz from left to right, the dash lines present the hysteresis losses for each frequency, the dot lines present the eddy losses for each frequency.
Figure 4. Time-series core loss of 1000 Hz - 1.0 T for one period, the solid line shows the total core loss, the dot line shows the eddy loss, the dash line shows the hysteresis loss and the fine line shows the permeability loss.
4. Conclusions
A interesting phenomenon is found when we apply the physical core-loss equation to several electrical steel sheets of China Steel Corporation, the optimum value of coefficient m of the hysteresis-loss equation is equal to 2 precisely, and the others approach a very small difference from the optimum when m is fixed at 2. For example of 65CS400, the coefficients of fixed m = 2 are converged to ceddy = 0.82953, r1 = 1.63, r2 = −0.007, cperm = 0.192, n1 = 0.93, n2 = 2.15, chyst = 53, B0 = 0.62, and the average error
of 3.41% is very closed to the optimum 3.38%. m fixing at 2 which means the hysteresis loss is only related to magnetic energy in physical explanation. To reduce one variable makes the converging process easier and faster.
The physical core-loss equation can serve many kinds of machines from the low-speed motors to high-frequency pulse transformers. It performs so good predictability for wide range of frequencies and magnetic densities, therefore we can consider that each coefficient of the equation might be associated with the physical properties, such as alloy composition, grain size or sheet thickness and so on. To study the influence of sheet thickness on the 8 coefficients will be our next subject.
Acknowledgements
I would like to acknowledge Metal Industries Research & Development Centre which collectively funded this project.
Symbol Definition
: The magnetic flux inside the magnetic core
: The magnetic flux density
: The peak value of magnetic flux density for frequency domain
: The frequency of periodic magnetic field
: The magnetic permeability of the magnetic core
: The electrical resistivity of the magnetic core
: The section area of the magnetic core
: The loop length of the magnetic core
: The sheet thickness of the magnetic core
: The turn number of the coil
: The coil voltage
: The coil current
: The eddy current of the magnetic core
: The input power of unit volume
: The magnetic inner power of unit volume
: The eddy loss of unit volume
: The simplified magnetic permeability loss of unit volume
: The magnetic hysteresis loss of unit volume
: The dynamic core loss of unit volume
: The hysteresis loss energy