The Estimation Method Comparison of Iron Loss Coefficients through the Iron Loss Calculation

Journal of Electrical Engineering and Technology.
2013.
Nov,
8(6):
1409-1414

- Received : January 25, 2013
- Accepted : June 05, 2013
- Published : November 01, 2013

Download

PDF

e-PUB

PubReader

PPT

Export by style

Article

Metrics

Cited by

TagCloud

A new calculation method for iron loss coefficients is proposed by using the Steinmetz equation from Epstein data. The hysteresis loss must have linear characteristic according to the frequency. However, the existing iron loss coefficients are defined by formula of frequency. In this case, the hysteresis loss has non-linear characteristics by frequency. So, in this paper, the iron loss coefficients were defined by a function of the magnetic flux density, and the iron loss calculation is applied for Interior Permanent Magnet Synchronous Motor(IPMSM) of 600(W) and 200(W). The iron loss calculation results and the experimental results are compared according to the various materials.
The shape of the applied model
Specifications of the applied model
where,
f
is the frequency of the external magnetic field,
B_{m}
is the maximum value of magnetic flux density,
k_{h}
is the hysteresis loss coefficient,
k_{e}
is the eddy current loss coefficient,
k_{a}
is the abnormal eddy current loss coefficient, and
n
is the Steinmetz constant
[2]
.
Eq. (1) responds rapidly to the eddy current loss term, but it does not correspond to hysteresis loss term by in the increased magnetic flux density in the core material and PM. Therefore, the Steinmetz constant is fixed to 2. Because, the magnetic flux density was increased by improvement of the PM and iron core material in electrical machine. Eventually, the other coefficients were estimated. In conclusion, the numerical formula of iron loss can be obtained by (2).
The provided Epstein data from the manufacturer are shown in
Fig. 2
. It did not contain the required information by magnetic flux density and frequency. So, the iron loss must be recalculated from the desired frequency band by using the provided data. The Curve Fitting Method(CFM) and iron loss calculation of IPMSM were performed on the basis of the provided Epstein data. The estimated coefficients were satisfied at 50(Hz), 60(Hz), 100(Hz) and 200(Hz). Therefore, the iron loss coefficients were estimated up to the 20th harmonic of operating frequency. The calculation result was obtained by summation of harmonic components in iron loss. The
Table 2
shows the iron loss coefficient fitting results according to frequency.
Epstein data by the manufacturer
Calculation results of iron loss coefficients according to the frequency
For the calculation of the iron loss, the flux density was calculated through the FEM by dividing one stator teeth into 55 areas, as shown in
Fig. 3
. In this paper, the driving frequency is 200 [Hz]; the magnetic flux density in the yoke and the teeth was calculated by using the finite element method from the iron loss data of 50PN1300
[3]
.
Divided stator core shape
The estimated coefficients of
Table 2
were inserted in (2). Thus, the calculated iron loss results and loss ratio were shown in
Tables 3
~
6
.
Iron loss calculation results in 600(W); CFM
Iron loss ratio calculation results in 600(W); CFM
Iron loss calculation results in 200(W); CFM
Iron loss ratio calculation results in 200(W); CFM
In operating frequency, hysteresis loss is 20%~30%, eddy current loss is 60%~70%, and abnormal eddy current losses is 10% in the total iron loss
[4
,
5]
. As shown in (1), Steinmetz constant should be estimated and set as a variable. However, the proposed method of this paper, Steinmetz constant value was fixed, and other coefficients were estimated more accurately. Steinmetz constant was fixed, and the three coefficients were calculated. This result has a strong resemblance to an iron loss aspect of rotating machinery.
Table 4
,
6
shows the calculation results of hysteresis loss and eddy current loss ratio.
The calculation values were analyzed. As a result, the abnormal eddy current losses according to increasing frequency were reduced at 600(W) class motor, on the other hand, that was increased at 200(W) class motor. The difference between the calculation and experimental results is similar in
Table 3
,
5
. But, the results are different with an aspect of abnormal eddy current loss of (2).
Fig. 4
shows results of iron loss divided by frequency
[6]
. In this case, the hysteresis loss should have a constant value according to the frequency. But, the hysteresis loss was reduced, because iron loss coefficients were defined as a function of frequency. To compensate the reduction of the loss, the iron loss coefficients were expressed by function of magnetic flux density in this paper.
The aspect of the iron loss divided by the frequency (CFM)
The result of iron loss according to frequency at each magnetic flux density
The iron loss coefficients were estimated for each magnetic flux density. That result was expressed for iron loss coefficients versus magnetic flux density. By using the calculated results of existing method of iron loss coefficients calculation, Each loss divided by frequency. As a result, the hysteresis loss was reduced, as shown in
Fig. 4
, because iron loss coefficients were defined as a function of frequency. To solve the problem, the iron loss coefficients were expressed by function of magnetic flux density as shown in (3)~(5). Then, the results can be represented as
Fig. 6
[1]
. The same process was applied to IPMSM of 200(W) class motor.
Iron loss coefficients calculation according toflux density(50PN1300)
Process for an iron loss experiment
Iron loss experiment equipment
Table 7
,
9
shows the results of calculation value and experimental value in 600(W) and 200(W). Experimental value and calculation value was compared. The error between experimental value and calculation value was 0.45(W) at operating frequency of 600(W) class motor. And the error of 200(W) class motor was 0.5(W) at the operating frequency.
W_{i} (W)cal : calculation data, W_{i} (W)exp : experimental data
W_{i} (W)cal : calculation data, W_{i} (W)exp : experimental data
The abnormal eddy current loss ratio was increased. However, it shows decreasing aspect at the specific frequency. Because, the iron loss was divided by frequency and the abnormal eddy current was a root function about frequency. Those results are shown in the
Table 8
,
10
.
Iron loss ratio calculation results in 600(W); cubic function
Iron loss ratio calculation results in 200(W); cubic function
The iron loss of (2) was divided by frequency, and the hysteresis loss has a constant value about the frequency. Consequentially, the iron loss coefficients were defined by function of the magnetic flux density. As a result, iron loss analysis results were different from those of
Fig. 4
. The results are shown in
Fig. 9
.
An aspect of the iron loss by the frequency (cubic function)
_{h}
divided by frequency f and the results W
_{h}
/
f
was reduced about increasing frequency, as shown in
Fig. 4
. This is because the iron loss coefficients
k
was defined as a function of frequency. To solve the problem, the iron loss coefficients
k
were proposed by function of magnetic flux density. And the calculated results of the iron loss, the hysteresis loss W
_{h}
was divided by frequency
f
. As a result, W
_{h}
/
f
has a constant value about the frequency, as shown in
Fig. 9
.
And, the iron loss by using proposed method, 600(W) class motor has hysteresis loss components of 37.43%, eddy current loss components of 50.24% and abnormal eddy current loss components of 12.33%. 200(W) class motor has hysteresis loss components of 43.25%, eddy current loss components of 43.48%, abnormal eddy current loss components of 13.26%. In conclusion, computation results well agreed with classically loss components ratio as widely known in existing iron loss calculation method.
Yong-Tae Kim He received B.S degree in electrical engineering from Changwon National university, Korea. He is Currently Pursuing his M.S degree at Changwon National University. His research interests are Electrical machine and FEM analysis.
Gyu-Won Cho He received B.S and M.S degree in electrical engineering from Changwon National University, Korea. He is currently pursuing his Ph.D degree at Changwon National University. His research interests are electrical machine and FEM analysis.
Gyu-Tak Kim He received B.S, M.S and Ph.D degree in electrical engineeering from Hanyang University, Korea. He is presently a Professor of Changwon National University. His research interests are Electrical machine and FEM analysis.

IPMSM
;
Steinmetz equation
;
New calculation method of iron loss coefficients
;
Curve fitting method

1. Introduction

The IPMSM can use a reluctance torque. As a result, the torque per unit volume is relatively large and IPMSM has high magnetic flux density. Therefore the iron loss is greater than the other motors. The formulation about iron loss is defined by Steinmetz. But, the research on iron loss is still being done
[1]
. The residual magnetic flux density of the Permanent Magnet(PM) has been increased and the iron core material has been improved too. Therefore, the previous iron loss coefficients are hard to be applied for recent ferromagnetic material. The existing calculation method of iron loss coefficients assumes that the Steinmetz constant n is unknown quantity. In this paper, it was fixed at 2 by considering enhanced magnetic flux density and iron core material. And, the other coefficients were calculated. So, the new calculation method of the iron loss coefficients is proposed by using the function of magnetic flux density. The coefficients calculation method was verified through experimental value of 600(W) and 200(W) class motor.
2. Application Model

As shown in
Fig. 1
, IPMSM of 600(W) and 200(W) class was used to calculate the iron loss. The specifications of the two models are described in
Table 1
. The 600(W) IPMSM used electrical steel 50PN1300 at iron core. In order to reduce the iron loss, 200(W) IPMSM used 50PN800 at iron core.
PPT Slide

Lager Image

Specifications of the applied model

PPT Slide

Lager Image

3. The Existing Method of Iron Loss Coefficients Calculation

Generally, the Steinmetz equation of iron loss is given as in (1).
PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

Calculation results of iron loss coefficients according to the frequency

PPT Slide

Lager Image

PPT Slide

Lager Image

Iron loss calculation results in 600(W); CFM

PPT Slide

Lager Image

Iron loss ratio calculation results in 600(W); CFM

PPT Slide

Lager Image

Iron loss calculation results in 200(W); CFM

PPT Slide

Lager Image

Iron loss ratio calculation results in 200(W); CFM

PPT Slide

Lager Image

PPT Slide

Lager Image

4. The Proposed Method of Iron Loss Coefficients Calculation

Fig. 5
shows Epstein data divided by frequency
[1]
.
PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

5. The Iron Loss Comparison of Experimented and Calculated Results

The iron loss experiment was performed through the process outlined in
Fig. 7
. And
Fig. 8
shows the iron loss experiment system. So, power analyzer is 99.5% accurate. Variation of copper loss was considered by measuring the input power and current from the inverter to the operating motor. A non-magnetized model was manufactured to measure the mechanical loss.
PPT Slide

Lager Image

PPT Slide

Lager Image

Iron loss calculation results in 600(W); cubic function

PPT Slide

Lager Image

Iron loss calculation results in 200(W); cubic function

PPT Slide

Lager Image

Iron loss ratio calculation results in 600(W); cubic function

PPT Slide

Lager Image

Iron loss ratio calculation results in 200(W); cubic function

PPT Slide

Lager Image

PPT Slide

Lager Image

6. Conclusion

In this paper, the new estimation method of iron loss coefficients were proposed and iron loss was calculated by using proposed method. So, calculation results of iron loss were verified from comparison of experimental and calculation value.
At the existing calculation method of iron loss coefficients, iron loss coefficients were defined as a function of frequency. As a result, the hysteresis loss W
Acknowledgements

This research was financially supported by the Ministry of Education, Science Technology(MEST) and National Research Foundation of Korea(NRF) through the Human Resource Training Project for Regional Innovation.

BIO

Ionel D. M.
,
Popescu M.
,
Dellinger S. J.
,
Miller T. J. E.
,
Heideman R. J.
,
McGilp M. I.
2006
“On the variation with flux and frequency of the coreloss coefficients in electrical machines”
IEEE Transactions on.Industry Applications
42
(3)
658 -
668
** DOI : 10.1109/TIA.2006.872941**

Bertotti G.
1988
“General Properties of Power Losses in Soft Ferromagnetic Materials”
IEEE Transactions on. Magnetics
24
(1)
621 -
630
** DOI : 10.1109/20.43994**

Lee Jeong-Jong
,
Kim Young-Kyoun
,
Nam Hyuk
,
Ha Kyung-Ho
,
Hong Jung-Pyo
2004
“Loss Distribution of Three Phase Induction Motor Fed By Pulse width Modulated Inverter”
IEEE Transactions on. Magnetics
40
(2)
762 -
765
** DOI : 10.1109/TMAG.2004.825445**

Chen Y.
,
Pillay P.
2002
“An Improved formula for lamination core loss calculations in machines operating with high frequency and high flux density excitation”
Industry Application Conference
2
13 -
18

Lee Ping-Kun
,
Kuo Kai-Chen
,
Wu Cheng-Ju
,
Wong Zuo-Tin
,
Yen Jia-Yush
2011
“Prediction of iron losses using the modified Steinmetz equation under the sinusoidal waveform”
Control Conference (ASCC)
579 -
584

Domeki H.
,
et al.
2004
“Investigation of benchmark model for estimating iron loss in rotating machine”
IEEE Transactions on. Magnetics
et al.
40
(2)
794 -
797
** DOI : 10.1109/TMAG.2004.825442**

Citing 'The Estimation Method Comparison of Iron Loss Coefficients through the Iron Loss Calculation
'

@article{ E1EEFQ_2013_v8n6_1409}
,title={The Estimation Method Comparison of Iron Loss Coefficients through the Iron Loss Calculation}
,volume={6}
, url={http://dx.doi.org/10.5370/JEET.2013.8.6.1409}, DOI={10.5370/JEET.2013.8.6.1409}
, number= {6}
, journal={Journal of Electrical Engineering and Technology}
, publisher={The Korean Institute of Electrical Engineers}
, author={Kim, Yong-Tae
and
Cho, Gyu-Won
and
Kim, Gyu-Tak}
, year={2013}
, month={Nov}