# Failure Diagnosis for Demagnetization in Interior Permanent Magnet Synchronous Motors

Authors: Takeo Ishikawa; Naoto Igarashi; Nobuyuki Kurita
Journal: International Journal of Rotating Machinery
 (2017)

---

## Abstract

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 thed-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.


---

## Body

## 1. 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–5]. They analyzed the stator current of a PM motor using several methods: Hilbert-Huang transform, method based on Box-counting 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–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, high-frequency 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.


## 2. Rotor Configuration of PMSM with Permanent Magnet Defect

Figure1(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 1
Rotor of the experimental PMSM.
(a)
Rotor configuration
(b)
Photograph of rotorFigure2(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, healthy means motor with fully magnetized PMs, R-2.5% means motor with PMs demagnetized in the radial direction by 2.5%, and Z-2.5% means motor with PMs demagnetized in the axial direction by 2.5%. PMs were manufactured by 0.1 mm step.Table 1
Dimension of PMs.
PM
L [mm]
W [mm]
H [mm]
V [mm3]
V
/
V
healthy
Healthy
22.6
20.2
2.2
16070
—
R-2.5%
22.6
20.2
2.0
15704
0.977
R-5.0%
22.6
20.2
1.8
15339
0.955
R-7.5%
22.6
20.2
1.6
14974
0.932
Z-2.5%
18.0
20.2
2.2
15661
0.975
Z-5.0%
13.4
20.2
2.2
15252
0.950
Z-7.5%
8.8
20.2
2.2
14843
0.924Figure 2
PMs for radial demagnetization.
(a)
Magnet configuration
(b)
Photograph of magnetFigure 3
PMs for axial demagnetization.
(a)
Magnet configuration
(b)
Photograph of magnet


## 3. 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 Figure4.Figure 4
Finite element mesh for the analyzed motor with demagnetized magnet.Figure5 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.Table 2
Measured flux linkage.
Motor
Φ [Wb]
Φ
/
Φ
healthy
Healthy
0.2157
—
R-2.5%
0.2108
0.978
R-5.0%
0.2055
0.954
R-7.5%
0.2002
0.930
Z-2.5%
0.2104
0.976
Z-5.0%
0.2048
0.950
Z-7.5%
0.1994
0.925Figure 5
Calculated EMF waveform.
(a)
Radial demagnetization
(b)
Axial demagnetizationFigure 6
Measured EMF waveform.
(a)
Radial demagnetization
(b)
Axial demagnetizationFigure7 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 z-direction. Even for axial demagnetization the magnetic flux density is 0.479 T at z = 0 and 0.477 T at z = 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.Figure 7
Calculated magnetic flux density at the center of the air-gap and axial location (z = 0).Figure 8
Calculated magnetic flux density at the center point of air-gap and axial direction for ratio of PM volume.Figure 9
Calculated flux density at the center of the air-gap andθ = 315°.Figure 10
Calculated magnetic flux density at the center point of air-gap andθ = 315° for ratio of PM volume.Figure11 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 R-2.5% is almost same as that for z-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.Figure 11
Calculated magnetic flux density waveform at the stator teeth tip.Figure 12
Fourier analysis of the magnetic flux density at the stator teeth tip.
(a)
Healthy
(b)
Radial demagnetization
(c)
Axial demagnetization


## 4. 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 u and w 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 R and L are calculated. Figures 15(a) and 15(b) show the relative error of L when the rotor position is 65 deg. and 20 deg., respectively. Figure 16 shows the relative error of resistance R. 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 R 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 d-axis is greater than that at q-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 d-axis position is better than that at the q-axis position.Figure 13
Measured impedance betweenu and w terminal at 100 Hz.
(a)
Inductance
(b)
ResistanceFigure 14
Measured impedance betweenu and w terminal at 2000 Hz.
(a)
Inductance
(b)
ResistanceFigure 15
Relative error of measured inductance.
(a)
Atd-axis
(b)
Atq-axisFigure 16
Relative error of measured resistance.
(a)
Atd-axis
(b)
Atq-axis


## 5. Stator Voltage and Current under Vector Control

One of the authors proposed a simulation model using MATLAB/Simulink as shown in Figure17 [12]. The model is a typical block diagram for the vector-controlled PMSM system. In this diagram, two control loops are used; one is an inner loop to regulate the stator currents by converting to the d- and q-axes with the rotor position, and the other is an outer loop to control the motor speed. The PWM voltage is generated by comparing the voltage reference and a truncated triangular carrier waveform as shown in Figure 17(b). In this block, the stator voltage drop is taken into account by considering an input voltage drop ΔVdc and an inverter voltage drop. Figure 18 shows the measured input voltage drop ΔVdc. This voltage drop can be represented by a linear equation and is represented by a blockVdcdrop as shown in Figure 17(c). It is well known that the forward characteristics of the IGBT and diode are nonlinear and vary with temperature. In order to simplify the analysis, this paper assumes that thev-i characteristics of the IGBT and diode are represented by constant forward voltage drops,v_Tr andv_D, respectively, and resistance,r_Tr andr_D, respectively. The inverter situations with positive current and negative current are shown in Figure 19. The letter s denotes a switching function for the transistor. When s = 1, the gate signal of the upper transistor is turned on, and when s = 0, the gate signal of the lower transistor is turned on. When s = 1 and the stator current is positive, the upper transistor is turned on, and then output voltage is shifted down byv_Tr+ i∗r_Tr as shown in case 1. Similarly, when the current is negative, output voltage is shifted up byv_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, d and q inductances Ld and Lq, 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 Ld and Lq 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 17
Block diagram for a brushless DC motor controlled using the vector control strategy.
(a)
Whole system
(b)
Subsystem forPWM inv in (a)
(c)
Subsystem for va,vb, andvc in (b)
(d)
Subsystemvoltage drop in (c)Figure 18
DC voltage drop across the inverter.Figure 19
Analysis of the output voltage of a phase inverter leg.Figure20(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 d- and q-axis currents id and iq. As id is controlled to be 0, the characteristic of iq 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 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.Figure 20
Simulated results using vector control.
(a)
Fundamental component of the stator current
(b)
Stator currentsid, iq
(c)
Fundamental component of the stator voltageExperimental 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 modelPM provided in SimPowerSystem toolbox of MATLAB/Simulink. Therefore the space harmonics are ignored in this mode. Although the actual inductances Ld and Lq 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.Figure 21
Measurement results under vector control.
(a)
Fundamental component of the stator current
(b)
Stator currentsid, iq
(c)
Fundamental component of the stator voltageHere, we will discuss these characteristics. Whenid is controlled to be 0, the voltage and torque equations are expressed by(1)vd=-ωLqiq,(2)vq=R+PLqiq+ωϕ,(3)T=pϕiq,where vd, vq, ω, ϕ, P, and p are d-axis voltage, q-axis voltage, angular velocity, flux linkage, differential operator, and pole pairs, respectively. From (3), as iq is inversely proportional to flux linkage ϕ, iq 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 id = 0. Note that a variation of the stator resistance due to stator temperature variation produces the same effect.From (1) and (2), vq is proportional to ωϕ and vd becomes 0, when the load torque is small; that is, iq is small. In contrast, when the load torque becomes large, that is, iq becomes large, vq is affected by iq 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(4)T=pωΦhealthyΦfaultyω2Lq2+R2,where Φhealthy and Φfaulty are the flux linkages of healthy motor and demagnetized motor, respectively. Here it is assumed that R and Lq do not change by the demagnetization. Therefore, since vq is proportional to ωϕ when iq = 0 at no load, the stator voltage at no load is useful for the diagnosis of demagnetization.


## 6. Conclusions

This paper has investigated the diagnosis of very slight PM demagnetization by both simulations and experiments. The reduction of the flux linkage, that is, back-EMF of the motor with demagnetized PMs, is equal to the reduction of PM volume, and there is no difference in the demagnetization in the radial direction and in the axial direction. Magnetic flux density at the teeth tip of the motor with demagnetized PMs includes even harmonics and subharmonics as well as 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 thed-axis position is better than that at the q-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.



---

*Source: 2716814-2017-10-16.xml*