An Accurate Iron Loss Analysis Method based on Finite Element Analysis considering Dynamic Anomalous Loss Katsuyuki Narita JMAG Buisiness Company JSOL Corporation Tokyo, Japan narita.katuyuki@jsol.co.jp Ryosuke Akaki Electric Component Development Dept. SUZUKI Motor Corporation Shizuoka, Japan ryosukeakaki@hhq.suzuki.co.jp Abstract To estimate the iron loss with higher harmonics accurately, the method based on the finite element analysis is proposed. The accuracy of the conventional methods based on the empirical approach deteriorates at the high frequencies and in the case with the minor loops, because they rely on the material measurement that differs from the actual operating conditions. The proposed method can calculate the hysteresis losses with the minor loops accurately using the Play-Hysteron model. It can estimate the classical eddy current losses at higher frequencies using the 1D method. Regarding the anomalous eddy current loss, dynamic anomalous loss factor is introduced in this research. The iron losses due to the minor loops of the lamination steel sheets are measured and compared with the simulation result using the Play-Hysteron model and the 1D method. The iron losses are also measured under the conditions considering the effect of the slot harmonics and PWM carrier harmonics. The proposed method reproduces well the measurement. Keywords Hysteresis, Anomalous loss, finite element analysis (FEA) I. INTRODUCTION In recent years, as demands for the miniaturization and the high speed rotation have been increasing with automobile driving motors etc., the magnetic flux waveform inside a motor tends to be distorted. In addition, the optimization of overall efficiency of a motor drive system including an inverter is required. From such backgrounds, the accuracy of the iron loss calculation under the magnetic flux waveform in motors including the spatial harmonics and time harmonics is required. It is well known that iron losses in a motor greatly increase when the spatial harmonics in magnetomotive force and the time harmonics caused by inverter carriers overlap with the fundamental magnetic flux wave. The loss evaluation methods based on the empirical approach proposed by Steinmetz [1] have been widely used for electrical machine analyses. It is pointed out that such conventional method underestimates iron losses of actual machines with the waveform including harmonics because it utilizes the iron loss properties under sinusoidal magnetic flux excitation condition. On the other hand, the Play-Hysteron model [2] and the 1D method [3][4] have been proposed to calculate iron losses by arbitrary magnetic flux waveforms. We applied the Play-Hysteron Hiroyuki Sano JMAG Buisiness Company JSOL Corporation Tokyo, Japan sano.hiroyuki @jsol.co.jp Masahiro Aoyama Department of Electrical and Electronics Enginnering Shizuoka Uvinversity Shizuoka, Japan aoyama.masahiro@shizuoka.ac.jp 978-1-4799-7312-5/18/$3 218 IEEE 439 Takashi Yamada JMAG Buisiness Company JSOL Corporation Tokyo, Japan yamada.takashi@jsol.co.jp model and the 1D method to the iron loss analysis of the switched reluctance motor and showed that the iron loss calculation accuracy greatly improved [5]. However, in order to put the Play-Hysteron model and the 1D method into practical use, there remains big issue that an anomalous eddy current loss are not taken into consideration in the 1D method. The influence by anomalous eddy current losses are especially remarkable in motors for automobile driving because high grade steel sheets which have large anomalous eddy current losses are often used. First, we state the definition of the anomalous eddy current losses in this research. We define the anomalous eddy current loss as the eddy current loss except classical eddy current loss, a loss due to eddy currents generated in a uniform conductor. One of the cause of the anomalous eddy current loss is magnetic domain wall movement. A magnetic body such as an electromagnetic steel sheet has magnetic domains, and their magnetic domain walls move or bow due to an alternating magnetic field. The anomalous eddy current loss refers to the eddy current loss caused by the movement of the magnetic domain walls. At present, it seems that the practical physical modeling has not been established for the anomalous eddy current loss calculation. The empirical approach have been studied to consider anomalous eddy current losses as the correction factor (anomalous loss factor : κ) to classical eddy current losses given by the 1D method. The conventional method to decide the anomalous loss factor κ is shown in Fig.1. In this method, the anomalous loss factor κ is often determined by the ratio of the eddy current loss including the anomalous eddy current loss separated from the measured iron loss by the twofrequency method and the classical eddy current loss W e by the theoretical formula (1) that does not consider the skin effect [6]. W e = 2 h 2 B 2 m f 2 (1) 6 σ:electric conductivity(1/(ohm m)), h:thickness(m) Bm: magnetic flux density amplitude(t),f :frequency(hz) This separation method needs to be performed at the low frequencies at which the skin effect can be ignored. When the
anomalous loss factor κ at the low frequencies is applied to the iron loss calculation at the high frequencies, the accuracy may be deteriorated. In the proposed method, the anomalous loss factors κ for arbitrary magnetic flux density amplitudes B m and frequencies f are determined by separating the anomalous eddy current loss at each conditions. We report that the accuracy of the iron loss calculation caused by harmonics is greatly improved by using those κ (B m,f). Fig. 1. Conventional anomalous loss separation method To improve above issue, the anomalous eddy current losses are separated by the method shown in Fig. 3, and the anomalous loss factors (B m,f) are obtained. Measure the iron losses for the sinusoidal magnetic flux with each amplitudes B m and frequencies f. Utilize the 1D method. In the 1D method calculation, the non-linear permeability is given by the Play- Hysteron model. Calculate the hysteresis loss and the classical eddy current loss from the result of the 1D method. Calculate the anomalous eddy current loss by subtracting the hysteresis loss and the classical eddy current loss from the measured iron loss. Decide κ (B m,f) as κ (B m,f)=(classical eddy current loss + anomalous eddy current loss) / classical eddy current loss. Fig. 4 shows the results of the calculated κ (B m, f) by applying above method to the electromagnetic steel sheet 35A36. It tends to decrease as the frequency increases. II. PROPOSED METHOD A. Separate anomalous loss from iron loss To separate anomalous eddy current losses for arbitrary magnetic flux density amplitudes and frequencies, it is necessary to calculate accurately hysteresis losses and classical eddy current losses in consideration of the skin effect due to nonlinear permeability. For that purpose, we first calculated the classical eddy current loss by the 1D method with the initial magnetizing curve as the nonlinear permeability. However, that method often gave κ (B m,f) less than, in other words, the anomalous eddy current loss is negative. For example, Fig. 2 shows κ (B m,f) of 5A47 by the initial magnetizing curve. We think that phenomenon is caused by incorrect permeability of the initial magnetizing curve. Because the permeability of the initial magnetizing curve is usually smaller than the permeability of DC hysteresis curve, the skin effect becomes weak and the classical eddy current loss by the 1D method is overestimated. Anomalous loss factor, 1.4 1.3 1.2 1.1.9.8.7.6.5 5Hz 2Hz 1Hz 4Hz.3.6.9 1.2 Flux density(t) Fig. 2. Anomalous factor of 5A47 given by initial magnetizing curve Fig. 3. Proposed anomalous loss separation method B. Decide dynamic anomalous loss factor In order to apply the anomalous loss factor to arbitrary waveforms by the transient response magnetic field analysis, (B m, f) need to be converted to the dependency on the equivalent magnetic flux density variation db/dt. In the case of the sinusoidal magnetic flux which has the amplitude B m and the frequency f, the effective value of db/dt is given by (2). db = 2 B mf dt 2 Those (db/dt) after conversion are called as dynamic anomalous loss factor[7]. Fig.5 shows (db/dt) of 35A36 converted from (B m, f) shown in Fig.4. Fig.5 shows (db/dt) tend to decrease as db/dt increases, and it can be fitted approximately by the logarithmic function of db/dt as shown by (3) (solid line in Fig. 5). When referring to the frequency and the magnetic flux density outside the measurement range, (db/dt) is decided by extrapolating by the function. We attempted the same method for several other grades of steel sheet and those (db/dt) were possible to fit in the same way. (2) 431
db = db (3) dt dt C. Apply anomalous factor to arbitrary magnetic flux waveform For arbitrary magnetic flux waveforms including harmonics, the anomalous loss factors κ ' (db/dt) are determined by (3) from db/dt between the time steps of the FEA, and the anomalous eddy current loss is calculated by multiplying the classical eddy current loss by the average value of κ ' (db/dt) for one electric cycle. Anomalous loss factor, 1.8 1.7 1.6 1.4 1.3 1.2 1.1 5Hz 2Hz 1kHz 1Hz 4Hz.3.6.9 1.2 Flux density(t) harmonics, the calculation accuracy of minor loop losses by the Play-Hysteron model was confirmed by comparing with the measurement. A. Measurement and analysis condition The minor loop losses were measured by the DC B-H analyzer to remove the effect by the eddy currents. The ring shape test piece which has 127mm outer diameter and 12mm inner diameter was made from 35A36. The amplitude of minor loops B m were.5,.1,.4 T and the maximum magnetic flux density of the minor loops B max were.2 1.8 T (see Fig.6). Regarding the set of symmetric loops used for the Play- Hysteron model, the amplitude of the symmetric loops were from.1 to 1.8 T and number of loops was 18(the increments of amplitude is.1 T). Fig. 4. Anomalous loss factor(35a36) Anomalous loss factor, Fig. 5. 1.8 1.7 1.6 1.4 1.3 1.2 1.1 2 4 6 8 db/dt Dynamic anomalous loss factor(35a36) III. VALIDATION OF PLAY-HYSTERON MODEL FOR ESTIMATING MINOR LOOP LOSS The iron losses by harmonics consist of the hysteresis losses, the classical eddy current losses, and the anomalous eddy current losses. Among those components, the classical eddy current losses are calculated by inputting the magnetization characteristics and the electric conductivity to the 1D method, so it is considered that the error factors are relatively small. On the other hand, since the shape functions used in the Play-Hysteron model are identified from the set of measured major loops (symmetric loops), the Play-Hysteron model can not necessarily reproduce the minor loop losses. Regarding the anomalous eddy current loss, it is separated from the iron loss as described above, because no practical physical model has been established yet. Therefore, there are two error components, which makes the error evaluation difficult. Prior to the verification of the iron losses by 4311 Fig. 6. Measurement of minor loop B. Comparison between measurement and analysis Fig.7 shows the result of the minor loop losses with B m=.1 T. There are two lines given by the Play-Hysteron model. One is the result by using the measured symmetric loop set as it is and the other is the result by interpolating the measured symmetric loop set into.5 T increments by using the appropriate interpolation technique [8]. From the comparison with the measurement, there is the deviation from the measurement when using the symmetric loop set with.1 T increments. In the principle of the Play-Hysteron model, the symmetric loop set with.1 T increments can certainly calculate minor loops with.1 T amplitude. However, in this case, the minor loops with.1 T amplitude are reproduced with only four polygonal lines, and the error occurs because resolution is insufficient to sufficiently reproduce the actually measured minor loop shape. On the other hand, the symmetric loop set with.5 T increments created by the interpolation technique can reproduce the measured minor loops. As a result, it was confirmed that by using the symmetric loop set that sufficiently reproduces the amplitude of the minor loops to be calculated, the Play-Hysteron model can calculate minor loop losses at least for this type of steel sheet and excitation conditions.
Hysteresis loss(w/kg).16.14.12.1.8.6.4.2 B=.1T B=.5T(interpolated) Measured.5 1 Bmax(T) Fig. 7. Hysteresis loss of minor loop(b m:.1t) Fig. 8. Measurement system IV. VALIDATION OF ANOMALOUS LOSS CALCULATION A. Measurement and analysis condition The measurement and analysis were done for same test piece of 35A36 for minor loop measurement. Fig. 8 shows the measurement system used to measure the iron losses for the anomalous eddy current loss separation and the iron losses due to harmonics. The magnetic flux is measured from the induced voltage of the secondary coil(12 turn). The magnetic field strength decided from the current of the primary coil(117 turn). The feedback control is performed so that the magnetic flux waveform matches the command waveform. The iron losses are calculated from the area of B-H loops. The iron losses due to the harmonics were measured under the following two conditions. Assuming the spatial harmonics Fundamental waveforms + 2%, 7th order harmonics Assuming the time harmonics due to a PWM invertor Fundamental waveforms + 5%, 6 khz harmonics Fig.9 shows the example of the measured waveforms. In the magnetic field analysis, 2D FEA model for the part of the ring was used for the loss calculation. The 2D FEA model is shown in Fig.1. The magnetic flux density waveforms obtained by the measurement were input as the vector potential boundary. The iron losses were calculated by the Play-Hysteron model and the 1D method, and the anomalous eddy current losses were calculated with the dynamic anomalous loss factor (db/dt) in Fig. 5. In the 1D method, although the influence to the eddy current by the edge of the steel sheet is ignored, we confirmed that the influence is small since the width of the ring is enough. Regarding the set of symmetric loops used for the Play-Hysteron model, the amplitude of the symmetric loops were from.5 to 1.8 T and number of loops was 36(the increments of amplitude is.5 T) and the increments were increased to.25 T by above interpolation technique. (a)magnetic flux density waveform (b)b-h loops Fig. 9. Example of measured waveform(t,2hz+6khz harmonics) Fig. 1. 2D FEA model (2.5deg of the ring) B. Comparison between measurement and analysis Fig. 11 and Fig.12 shows the measured and analyzed iron losses. "Conv." in Fig.11, Fig.12 show the results by the conventional iron loss calculation method using the iron loss values under the sinusoidal magnetic flux excitation. The eddy current losses of "Conv." include the anomalous eddy current 4312
losses. "Prop." shows the results by the proposed method. From the results, when superimposing the 7th harmonics on the 1T fundamental wave, the errors from the measurement are less than 1% even in the conventional method, but the errors increase in the case of 1.4T fundamental wave. As can be seen from Fig. 7, the minor loop losses increase when the fundamental wave amplitude becomes high and the minor loops occurs at the high position in the major loop, but that phenomenon can t be taken into account in the conventional method. The errors in the proposed method are less than several percent regardless of the fundamental wave amplitudes and frequencies. Looking at the breakdown of the iron losses, the contribution of the Play-Hysteron model is large at the low frequencies because the ratio of the hysteresis losses is high, and the contribution of 1D method is large at the high frequencies. "Prop. (Const)" shows the results where the anomalous loss factor κ is fixed at the value at1t/5hz (κ=8) by the conventional separation method in Fig.1. It is greatly overestimated at the high frequencies, resulting in deterioration in the accuracy over the conventional iron loss calculation method. Finally, the influence of the extrapolation by the logarithmic function shown in (3) was investigated. Fig.13 shows the errors to the measurement iron losses(fundamental + 7th harmonics) when narrowing the range of measurement data used for deciding the anomalous loss factors from 5-1kHz to 5-4Hz and 5-1 Hz. Fig.13 shows only low frequencies data such as 5 and 1 Hz for deciding the anomalous loss factor is not enough for calculating the iron loss at 1kHz fundamental wave. It causes almost 1% error. If the data up to 4 Hz exist, the iron losses at the 1kHz fundamental wave can be reproduced by the extrapolation by the logarithmic function (3). V. CONCLUSIONS AND FUTURE WORK By applying the dynamic anomalous loss factors identified from arbitrary magnetic flux density amplitude and frequency conditions by the Play-Hysteron model and the 1D method. The obtained dynamic anomalous loss factors were applied to the calculation of the iron losses due to harmonics. The results match well to the measured losses. In the future, we will apply this method to the iron loss calculation of motors and verify the accuracy. Future tasks include studying the physical background of fitting the dynamic anomalous loss factors with the logarithmic function and improvement of the fact that the factor κ becomes below at the high frequencies. In addition, it is necessary to evaluate the errors by applying the anomalous loss factors κ identified by the sinusoidal magnetic flux excitation conditions to the minor loops. 2.5 2..5. 6 5 4 3 2 1 5 4 3 2 1 12 1 8 6 4 2 (a) T,5Hz+7th order (b) T,4Hz+7th order (c)1.4,1hz+7th order (d) 1.4T,4Hz+7th order Fig. 11. Iron losses due to 7th orderharmonics 4313
12 1 8 6 4 2 (a) T,1Hz+6kHz 3 25 2 15 1 5 25 2 15 1 5 12 1 8 6 4 2 (b) T,4Hz+6kHz (c) 1.4T,1Hz+6kHz (d) 1.4T,4Hz+6kHz Solid circle : data exist, Hollow circle : data don t exist(extrapolated) Fig. 13. Measurement data range to use identification and error REFERENCES [1] G. Bertotti, General properties of power losses in soft ferromagnetic materials, IEEE Transactions on Magnetics, vol. 24, no. 1, pp.621-63, 1988 [2] T. Matsuo, D. Shimode, Y. Terada, M. Shimasaki, Application of stop and play models to the representation of magnetic characteristics of silicon steel sheet", IEEE Transactions on Magnetics, vol.39, no.3, pp.1361-1364, 23 [3] O. Bottauscio, Advanced Model of Laminated Magnetic Cores for Two-Dimensional Field Analysis, IEEE Transactions on Magnetics, vol.36, no.3, pp561-573, 2 [4] K. Yamazaki, N. Fukushima, Experimental Validation of Iron Loss Model for Rotating Machines Based on Direct Eddy Current Analysis of Electrical Steel Sheets, International Electric Machines and Drives Conference (IEMDC), pp.851-857, 29 [5] K. Narita, T. Asanuma, K. Semba, H. Sano, T. Yamada, K. Aiso, K. Akatsu, An accurate iron loss evaluation method based on finite element analysis for switched reluctance motors, Energy Conversion Congress and Exposition (ECCE), 215 [6] K. Yamazaki, N. Fukushima, "Iron-Loss Modeling for Rotating Machines: Comparison Between Bertotti's Three-Term Expression and 3-D Eddy-Current Analysis", IEEE Transactions on Magnetics, vol.46, no.8, pp. 3121-3124, 21 [7] R. Mitsuoka, T. Mifune, T. Matsuo, Dynamic hysteresis modeling of silicon steel having nonuniform magnetic property, J. Appl. Phys. 111, 7D11 (212) [8] K. Narita, H. Sano, T. Yamada, Study of processing and interpolation method of symmetric loops for play model, The Papers of Joint Technical Meeting on Static Apparatus and Rotating Machinery, IEE Japan, SA-17-58, RM-17-89 (217) (in Japanese) Fig. 12. Iron losses due to 6kHz harmonics 4314