Failure Diagnosis for Demagnetization in Interior Permanent Magnet Synchronous Motors

Since a high degree of reliability is necessary for permanent magnet synchronous motors, the detection of a precursor for the demagnetization of permanent magnets is very important.This paper investigates the diagnosis of very slight PM demagnetization. The permanent magnet volume is altered so as to mimic the effect of demagnetization. This paper investigates the influence of demagnetization by using several methods: the 3D finite element analysis (FEA) of the motor, the measurement of high-frequency impedance, and the measurement and FEA of the stator voltage and current under vector control. We have obtained the following results.The back-EMF is proportional to permanent magnet volume, and there is no difference in the demagnetization in the radial direction and in the axial direction. Even harmonics and subharmonics of flux density at the teeth tip could be useful for diagnosis if a search coil is installed there. The relatively low frequency resistance at the d-axis position is useful for diagnosis. Under vector control, the stator voltage is useful except in an intermediate torque range, and the intermediate torque is expressed by a simple equation.


Introduction
Permanent magnets (PMs) of permanent magnet synchronous motors (PMSMs) can be demagnetized by temperature, large stator currents, large short-circuit currents produced by the inverter or stator faults, and the aging of magnet itself.Moreover, because the PM material is often magnetized after it is inserted into the rotor assembly, there is a possibility that the PM is not in a state of complete magnetization.Since a high degree of reliability is necessary for PMSMs, the detection of a precursor to the demagnetization of PMs is very important [1].
There have been several research studies that explored methods for detecting the demagnetization of a PM.Rajagopalan et al. created a magnet defect by chipping off a part of the magnet and then measured the stator current [2].Urresty et al. made several contributions [1,[3][4][5].They analyzed the stator current of a PM motor using several methods: Hilbert-Huang transform, method based on Boxcounting fractal dimension, zero-sequence voltage component, and back-EMF based method.Liu and Zhu also showed that the magnetization can be estimated from the back-EMF [6].Moreover, in order to determine the magnet temperature caused by the demagnetization, some approaches have been proposed [7].One possible way to determine the rotor temperature under real-time conditions is a precise flux observer [8].The other method is an invasive method, which estimates the impact of the flux variation on the motor's magnetic circuit from the current response to injected voltage signals [9,10].However, in [1][2][3][4][5] the volume of the PM was decreased by 16.7%.The degree of PM demagnetization in these past studies was relatively large.It would be preferable to be able to detect the onset of demagnetization to provide the user with advance notice that the magnet is beginning to be demagnetized.
Although all poles may be uniformly demagnetized owing to the causes of demagnetization described above, the previous research studies created the magnet defect by chipping off a part of the permanent magnet or by reducing the volume of only one pole.In this study, the volume of only one pole is also reduced to simulate the demagnetization.In order to deal with a less demagnetized situation compared to those in the previous research studies, the volume of one of the four PMs is reduced by 10%, 20%, and 30%; that is, the total volume of PM material used is decreased by 2.5%, 5.0%, and 7.5%, respectively, in separate experiments [11].Two types of demagnetization are taken into account, that is, demagnetization in the radial direction and axial direction.In this paper, physical variables to be available for diagnosis for less demagnetization in interior permanent magnet synchronous motors are clarified by simulations and experiments.First, electromotive force, inductances, highfrequency impedance, and magnetic flux density in the air-gap are investigated.Next, the voltage and current are investigated when the motor is driven under a vector control strategy.

Rotor Configuration of PMSM with Permanent Magnet Defect
Figure 1(a) shows the rotor configuration of an experimental PMSM.This motor's rating is as follows: 1.5 kW, 3000 min −1 , 4.8 N m, 5.6 A, four-pole machine.The PM volume of one of four poles is altered so as to mimic the effect of demagnetization.Figure 1(b) shows a photograph of the rotor of the experimental PMSM.When performing experiments with several motors, slight differences between motors can affect the relative motor performance especially if the amount of demagnetization is small.In order to avoid this problem, we perform experiments on the same stator and the same rotor using the different PMs. Figure 2(a) shows the demagnetization in the radial direction.We define this operating condition as radial demagnetization hereafter.One pole is composed of four PMs, and the thickness of the four PMs which compose one pole is reduced by 10%, 20%, and 30%, in separate experiments, in order to imitate demagnetization or imperfect magnetization as shown in Figure 2(b).Figure 3(a) shows the demagnetization in the axial direction.We define this operating condition as axial demagnetization hereafter.The axial length of two of the four PMs per pole is reduced by 20%, 40%, and 60% as shown in Figure 3(b).Therefore, the PM volume of one pole for axial demagnetization is approximately the same as that for radial demagnetization.In order to remove eccentricity, nonmagnetic materials are inserted into the region where the PM is decreased.The volume of PMs is shown in Table 1.Here,

Electromotive Force of Demagnetized Motor
The flux distribution in the motor at no load is calculated using the 3D finite element method.The analysis model is 1/2 of the model in the axial direction, because there is no symmetric condition in the circumferential direction.The numbers of nodes and elements are 132,813 and 225,758, respectively, as shown in Figure 4. Figure 5 shows the calculated electromotive force, when the motor is rotating at 1,500 min −1 .It is shown that the electromotive force is proportional to the volume of PMs, and the waveform of electromotive force of the motor with radial demagnetization is the same as that of axial demagnetization.Figure 6 shows the measured electromotive force, when the motor is rotating at 1500 min −1 .It is shown that the waveforms of the measured electromotive force are approximately same as the calculated ones.In order to show the difference between the demagnetized and healthy situations, Table 2 shows the flux linkage, which is the integral of the electromotive force.It is well known that the magnetization of PMs can be estimated from the back-EMF [5,6,8] In this paper, it is found by both simulations and experiments that the reduction of the flux linkage of the motor with demagnetized PMs is equal to the reduction of PM volume.There is no difference in the back-EMF of the motor with PMs demagnetized in the radial direction and in the axial direction.Figure 7 shows the calculated magnetic flux density at the center of the air-gap and at the center of the axial direction.The magnetic flux density is reduced from 0.516 T to 0.481 T for the radial demagnetization and from 0.516 T to 0.479 T for the axial demagnetization at  = 315 ∘ .Figure 8 shows the magnetic flux density at the center of air-gap and axial direction versus PM volume.It is found that the magnetic flux density is proportional to PM volume irrespective of the direction of demagnetization.Figure 9 shows the calculated magnetic flux density at the center of the air-gap and at  = 315 ∘ .The magnetic flux density is almost flat in the direction.Even for axial demagnetization the magnetic flux density is 0.479 T at  = 0 and 0.477 T at  = 42 mm.Figure 10 shows the magnetic flux density at the center of air-gap and  = 315 ∘ versus the ratio of PM volume.It is also found that the magnetic flux density is proportional to PM volume irrespective of the direction of demagnetization.

International Journal of Rotating Machinery
Figure 11 shows the calculated waveform of magnetic flux density at the stator teeth tip.It looks like a rectangular shape.Note that the calculated magnetic flux density for -2.5% is almost same as that for -2.5%. Figure 12 shows its Fourier components.Magnetic flux density for healthy motor includes the fundamental and odds harmonics.On the other hand, magnetic flux density for the motor with demagnetized PMs includes even harmonics and subharmonics as well as the fundamental and odds harmonics.Therefore, if a search coil is installed at the teeth tip of the stator, the magnetization in one of four poles could be detected by the even harmonics and subharmonics of the flux density.There is no difference between the Fourier components of the motor with PMs demagnetized in the radial direction and in the axial direction.

Impedance for High-Frequency Injection
A high-frequency signal injection has been widely used for the position sensorless control of PMSM, and it was also used PM temperature estimation [10].This paper investigates high-frequency impedance of PMSM for the diagnosis of very slight PM demagnetization.The impedance was measured when a high-frequency voltage is injected between  and  phases.Here, the motor is not rotating; that is, there is no fundamental voltage and current.Figures 13 and 14 show the measured impedance when the high-frequency voltage of 100 Hz and 2000 Hz is injected.The impedance was measured at each rotor position.It is found that the effect of demagnetization strongly appears at the resistance of 100 Hz.In order to show the difference clearly, relative errors of  and  are calculated.Figures 15(a) and 15(b) show the relative error of  when the rotor position is 65 deg.and 20 deg., respectively.Figure 16 shows the relative error of resistance .It is found from Figure 15 that the relative error is approximately constant but the value is smaller than the reduction ratio of PM volume.Figure 16 shows that the relative error of resistance  is larger than the reduction ratio of PM volume in a low frequency region.Therefore, the resistance of low frequency region is useful for the diagnosis of demagnetization of PM.It is shown from Figure 16 that the relative error of resistance at -axis is greater than that at -axis.Reference [10] showed that the PM temperature of surface PM machine can be estimated from the high-frequency resistance.This paper has clarified that the demagnetization of PMs in interior PM machine can be estimated from the relatively low frequency resistance.Moreover, it has been clarified that the resistance at the -axis position is better than that at the -axis position.

Stator Voltage and Current under Vector Control
One of the authors proposed a simulation model using MATLAB/Simulink as shown in Figure 17  characteristics of the IGBT and diode are nonlinear and vary with temperature.In order to simplify the analysis, this paper assumes that the v-i characteristics of the IGBT and diode are represented by constant forward voltage drops, v Tr and v D, respectively, and resistance, r Tr and r D, respectively.The inverter situations with positive current and negative current are shown in Figure 19.The letter  denotes a switching function for the transistor.When  = 1, the gate signal of the upper transistor is turned on, and when  = 0, the gate signal of the lower transistor is turned on.When  = 1 and the stator current is positive, the upper transistor is turned on, and then output voltage is shifted down by v Tr + i * r Tr as shown in case 1.Similarly, when the current is negative,  output voltage is shifted up by v D + i * r D as shown in case 3. Four cases shown in Figure 19 can be represented by a block shown in Figure 17(d).In the simulation block, the IPMSM is modeled in terms of the stator resistance,  and  inductances   and   , and the flux linkage.In order to evaluate the effect of demagnetization of the PM, the value of the flux linkage is reduced in the simulation.On the contrary, the stator resistance is fixed, and the inductances   and   of the demagnetized motor are assumed to be the same as those of a healthy motor, because the difference of inductances from those of healthy machine is very small as shown in Figures 14 and 15.There is no three-phase imbalance of EMF and inductances caused by one-pole demagnetization, because the stator windings are connected in series.Although the small vibration and acoustic noise may be produced, this paper assumes they are negligible.Figure 20(a) shows the fundamental component of the stator current with radial demagnetization.It is shown that the stator current is approximately proportional to the load torque and that the differences due to demagnetization are very small when the load torque is small, and they become large when the load torque becomes large.Figure 20(b) shows the and -axis currents   and   .As   is controlled to be 0, the characteristic of   is equal to the stator current shown in Figure 20(a).Figure 20(c) shows the fundamental component of the stator voltage.It is found that there is a difference due to demagnetization.When the load torque is small, the stator voltage of the motor demagnetized by 7.5% has the lowest value and that of the healthy motor has the highest value.In contrast, when the load torque is large, the stator International Journal of Rotating Machinery voltage of the motor demagnetized by 7.5% has the highest value.The stator voltage and current of the motor with axial demagnetization are the same as those shown in Figure 20.Experimental conditions are 0%, 10%, 20%, and 30% of rated torque and a speed of 1500 min −1 .Figure 21 shows the fundamental component of the stator current and stator voltage of the IPMSM with radial demagnetization.The fundamental component of the stator current becomes large when the load torque becomes large.When the load torque is small, the stator voltage of the motor demagnetized by 7.5% has the lowest value and that of the healthy motor has the highest value.In contrast, when the load torque  is large, the stator voltage of the motor demagnetized by 7.5% has the highest value.These results are approximately the same as those shown in Figure 20.The measured stator voltage and current of the motor with axial demagnetization were also approximately the same as the simulation results shown in Figure 20.We think that the differences between the simulation results and the measured ones are due to the following assumptions.The simulation block diagram shown in Figure 17 uses permanent magnet synchronous machine model PM provided in SimPowerSystem toolbox of MATLAB/Simulink.Therefore the space harmonics are ignored in this mode.Although the actual inductances   and   change according to the value of stator current, they are assumed to be constant.However, since the simulation results are approximately same as the measured ones, the simulation results have been experimentally verified.
Here, we will discuss these characteristics.When   is controlled to be 0, the voltage and torque equations are expressed by where V  , V  , , , , and  are -axis voltage, -axis voltage, angular velocity, flux linkage, differential operator, and pole pairs, respectively.From (3), as   is inversely proportional to flux linkage ,   of the demagnetized motor is large when the load torque is large, as shown in Figure 20(a).And the flux linkage is proportional to the magnet volume, that is, demagnetization situation.Therefore, the demagnetization of PM can be estimated by the stator current under vector control with   = 0. Note that a variation of the stator resistance due to stator temperature variation produces the same effect.From (1) and ( 2), V  is proportional to  and V  becomes 0, when the load torque is small; that is,   is small.In contrast, when the load torque becomes large, that is,   becomes large, V  is affected by   as shown in (2).As a result, the stator voltage of the demagnetized motor becomes large, as shown in Figure 20(c).The torque at the intersection point is driven as where Φ healthy and Φ faulty are the flux linkages of healthy motor and demagnetized motor, respectively.Here it is assumed that  and   do not change by the demagnetization.Therefore, since V  is proportional to  when   = 0 at no load, the stator voltage at no load is useful for the diagnosis of demagnetization.the fundamental and odds harmonics.The demagnetization of PMs in interior PM machine can be estimated from the relatively low frequency resistance, and the resistance at the -axis position is better than that at the -axis position for diagnosis.Under vector control, the stator voltage is useful except in an intermediate torque range, and the intermediate torque is expressed by a simple equation.

Figure 7 :Figure 8 :
Figure 7: Calculated magnetic flux density at the center of the air-gap and axial location ( = 0).

Figure 9 :Figure 10 :
Figure 9: Calculated flux density at the center of the air-gap and  = 315 ∘ .

Figure 11 :
Figure 11: Calculated magnetic flux density waveform at the stator teeth tip.

Figure 12 :Figure 13 :
Figure 12: Fourier analysis of the magnetic flux density at the stator teeth tip.

Figure 18 :
Figure17: Block diagram for a brushless DC motor controlled using the vector control strategy.

Figure 19 :
Figure 19: Analysis of the output voltage of a phase inverter leg.
Figure 20: Simulated results using vector control.
Figure 21: Measurement results under vector control.