# Dynamic Modeling and Analysis of Demagnetizing Rotor of Permanent Magnet Synchronous Motor

Authors: Bifeng Zhou; Guoning Tang; Yiping Luo
Journal: Shock and Vibration
 (2021)

---

## Abstract

This study performs a dynamic modeling and analysis of the axial uniform demagnetization rotor of a motor to investigate the demagnetization peculiar for permanent magnet synchronous motor (PMSM). First, the air-gap change in the motor is analyzed by constructing a dual coordinate system of stator and rotor, and the unbalanced magnetic pull (UMP) model under uniform demagnetization and eccentricity is constructed. Second, combined with the UMP and the self-gravity, a dynamic model of the rotor under a uniform demagnetization and eccentricity is established. Lastly, the accuracy of the mathematical model of UMP under uniform demagnetization and eccentricity is confirmed by the mutual corroboration of Maxwell simulation and MATLAB calculation; and based on the dynamics model, the dynamic characteristics of the rotor system under different degrees of uniform demagnetization are studied. This study provides a theoretical basis for the accurate demagnetization fault diagnosis and vibration control of PMSMs in the future.


---

## Body

## 1. Introduction

The stability of permanent magnet materials for permanent magnet synchronous motors (PMSMs) is affected by temperature, magnetic field, chemical factors, and time. Considering the actual working condition environment of PMSMs and according to research on the demagnetization mechanism of PMSMs [1–3], the demagnetization of PMSMs is inevitable. Substantial research on the demagnetization of PMSMs has focused on the control of the motor and the performance of the PMSM after demagnetization. Such research has covered the permanent magnet flux linkage of PMSMs [4–7], the fault-tolerant control of PMSM demagnetization [8–11], and the finite element modeling of PMSM demagnetization [12–14].However, PMSMs often exist as a key component in practical industrial systems, such as permanent magnet synchronous wind generators. Not only does the demagnetization of PMSMs affect the performance of the motor itself, but the vibration of the motor due to demagnetization also considerably affects the system where the motor is located. Thus, studying the demagnetization vibration of PMSMs is crucial.Existing studies on the demagnetization vibration of PMSMs are relatively rare. Xiang [15] analyzed the stiffness characteristics of the rotor system of the PMSM and investigated the nonlinear dynamic behavior influenced by the unbalanced magnetic pull (UMP) based on the Jeffcott rotor model. Liu [16] investigated the nonlinear oscillations of a PMSM based on a Jeffcott rotor-bearing system and considered the effects of the UMP, the nonlinear restoration forces due to the Hertz contact force and bearing clearance, the rotor weight, and the rotor mass eccentricity.Some scholars have studied the demagnetization vibration of PMSMs and achieved promising results. However, the existing research results only consider the UMP generated by PMSM demagnetization and do not consider the specific demagnetization form of the permanent magnet of the PMSM. That is, only the demagnetization phenomenon is considered, and the demagnetization state is not quantitatively described. Therefore, the resulting dynamic model does not consider the demagnetization quantitatively.Motivated by the above analysis, this study focuses on the dynamic modeling and analysis of the demagnetizing rotor of PMSMs. First, the uniform demagnetization state of permanent magnets is quantitatively described by the demagnetization angle and the decrease of air-gap magnet motive force (MMF) at the demagnetization angle (air-gap MMF drop percentage); on the basis of this description, the air-gap change in the motor is analyzed by constructing a dual coordinate system of stator and rotor, and the unbalanced magnetic pull (UMP) model under uniform demagnetization and eccentricity is constructed. Second, combined with the UMP and the self-gravity, a dynamic model of the rotor under a uniform demagnetization and eccentricity is established. Lastly, the accuracy of the mathematical model of UMP under uniform demagnetization and eccentricity is confirmed by the mutual corroboration of Maxwell simulation and MATLAB calculation, and based on the dynamics model, the dynamic characteristics of the rotor system under different degrees of uniform demagnetization are studied. This study provides a theoretical basis for the accurate demagnetization fault diagnosis and vibration control of PMSMs in the future.


## 2. Nonlinear Dynamic Model of a Rotor Under Demagnetization

### 2.1. Dynamical Model of a Rotor System

The demagnetization of the permanent magnet in the rotor of the permanent magnet synchronous motor (PMSM) results in an uneven air-gap flux density in the motor, resulting in unbalanced magnetic pull (UMP) acting on the rotor, which makes the rotor eccentricity, and the eccentricity of the rotor further changes the air-gap flux density of the motor, thereby affecting the dynamic characteristics of rotor operation.Figure1 shows a schematic of rotor eccentricity under demagnetization. O and O′ are the geometric centers of the stator and rotor, respectively; εt is the eccentricity between the geometric center of the rotor and the stator; R is the radius of the rotor; δ0 is the average of air-gap length; and δψ is the length of the air-gap corresponding to angle ψ in the rotating coordinate system x′O′y′.Figure 1
The schematic diagram of rotor eccentricity under demagnetization.To establish the rotor dynamic model of a PMSM under demagnetization, this study assumes the following:(1)
All damping is considered as linear viscous damping in the calculation process, and disregard the effects of nonlinearity; the concentrated mass method is adopted without considering the axial motion of each component [17].(2)
The case where the rotor’s permanent magnet is axially uniformly demagnetized is considered. The combined amount of changes in the corresponding air-gap magnet motive force (MMF) caused by demagnetization in any axially uniform state of the permanent magnet can be described by the following two parameters: demagnetization angle and the decrease of air-gap MMF at the demagnetization angle (air-gap MMF drop percentage).On the basis of the above assumptions, the dynamic model of the magnetic field-rotor coupling system can be obtained.(1)mx¨t+cxxx˙t+kxxxt=FMxt,my¨t+cyyy˙t+kyyyt=FMyt−mg,where x,y,x˙,y˙,x¨,y¨ are the displacement, velocity, and acceleration of the rotor in the x and y axes, respectively; m is the structural quality of the rotor; cxx,cyy are the operational damping in the x and y directions; kxx,kyy are the stiffness of the rotor in the positive directions of the x and y axes; and FMxt and FMyt are the UMP in the x and y directions, respectively.


### 2.2. Modeling of Demagnetization

This section analyzesFMxt and FMyt mentioned in the previous section.


#### 2.2.1. Air-Gap Flux Density

The demagnetization of permanent magnet in the rotor of permanent magnet synchronous motor (PMSM) results in unbalanced magnetic pull (UMP), which makes the rotor eccentric. Figure2 shows a cross section diagram of rotor demagnetization. With the center O′ of the rotor taken as the origin, the rotor rectangular coordinate system x′O′y′ is established. In such system, the x′-axis direction is opposite to the center line direction of the demagnetization angle, and the coordinate system rotates with the rotor. With the center of the stator O taken as the origin, the stator stationary coordinate system is constructed with the initial direction of the x′ and y′ coordinate axes in the coordinate system x′O′y′ of the rotor. Thus, the initial stator and rotor coordinate systems are parallel.Figure 2
The cross-section diagram of rotor demagnetization.The displacement variation of eccentricity between the stator and rotor of a PMSM is shown in Figure2. β is half of demagnetization angle, ψ is the space angle associated with the x′ axis in the rotating coordinate system x′O′y′, δ0 is the average air-gap length of the PMSM, and εt is the eccentricity scale of the rotor, and θt is the angle between the connecting line of the stator and rotor center O, O′ and the x axis. When the rotor runs, it satisfies(2)θt=ωt+θ0,where θ0 is the angle between the resultant force of the rotor and the x′ axis in the rotating coordinate system x′O′y′ and ω is the rotational angular velocity of the rotor.ξt=εt/δ0 is set as the eccentricity ratio of the rotor. In the rotating coordinate system x′O′y′, the rotor eccentricity caused by the demagnetization of the permanent magnet is considered, and the air-gap distribution δψ,t can be approximately expressed as follows [18]:(3)δψ,t≈δ0−εtcosψ−θ0=δ01−ξtcosψ−θ0=δ0−x′tcosψ−y′tsinψ,where x′t=εtcosθ0 and y′t=εtsinθ0 are the projection displacements of the eccentric vector in the x′O′y′ coordinate system. The permeability of the air gap is as follows:(4)Λt=μ0δψ,t≈μ0δ01+ξtcosψ−θ0=μ0δ02δ0+x′tcosψ+y′tsinψ,where μ0 is the permeability of vacuum.Considering the demagnetization of the permanent magnets of a PMSM, the fundamental wave of the air-gap MMF established by the current in the torque winding and the permanent magnet field in the rotor of the PMSM is as follows [18]:(5)Fψ,t=Fδcoswt−pψ−u,where p is the number of pole pairs of the torque winding, Fδ is the amplitude of the fundamental wave of the air-gap MMF, u is the spatial initial phase angle of the fundamental wave of the air-gap MMF, and w is the angular frequency of the winding current. The rotation angular velocity ω of the rotor and the angular frequency w of the winding current satisfy ω=60∗2∗π∗w/p.In the rotating coordinate systemx′O′y′ and in consideration of the axial direction uniform demagnetization of the rotor of a PMSM, the change of the amplitude of the air-gap MMF caused by the demagnetization of the permanent magnet is equivalently mapped to demagnetization angle β on both sides of the coordinate −x′ axis, and the decrease of air-gap MMF drops to 1−ζ∗F^ (at the demagnetization angle). At this time, the amplitude of the fundamental wave of the air-gap MMF Fδ is(6)Fδ=F^,ψ∈−π+β,π−β,1−ζF^,ψ∈−π,−π+β∪π−β,π,where β∈0,π and ζ⩽1 is the decrease of air-gap MMF at the demagnetization angle (air-gap MMF drop percentage).In the rotating coordinate systemx′O′y′, in consideration of the local demagnetization of the permanent magnet of a PMSM, the air-gap flux density in the motor is as follows:(7)Bψ,t=Λψ,tFψ,t=μ0δ02Fδδ0+x′tcosψ+y′tsinψcoswt−pψ−u.Assuming that the permanent magnet is in normal condition and the air gap is uniform, the amplitude of flux density of air-gap for pair polep is [18](8)B^=μ0δ0F^=μ0δ0324πNI2p,where N is the number of turns of each phase of the torque winding and I is the amplitude of the excitation current of the torque winding.


#### 2.2.2. Analysis of UMP

The magnetic tension formed on the boundary of a magnetic material with varying permeability (such as air and core) in the magnetic circuit of a PMSM is called Maxwell force, also known as detent force; the direction of the force is perpendicular to the boundary surface of the magnetic material.When analyzing the magnetic field of the air gap of a PMSM, we assume that the magnetic density of the air gap is symmetrically distributed. At this time, the resultant force of the Maxwell force is zero. However, considering the local demagnetization of the permanent magnet in PMSM, the magnetic density distribution in the air gap of the motor is not uniform. Thus, the resultant Maxwell force is not zero.The rotor of a PMSM is subjected to the Maxwell force of the nonzero resultant force; as such, the rotor produces eccentricity, and the magnetic density in the air gap of the PMSM changes, further altering the resultant Maxwell force. The larger the eccentricity of the rotor, the greater the Maxwell force.In accordance with the Maxwell tensor method, the magnetic density in the air gap of a PMSM is set asB, and the Maxwell force acting on the dS area of the rotor surface is as follows:(9)dFM=B2dS2μ0.Then the component of Maxwell force along thex, y direction is as follows:(10)dFMxψ=B2ψ,tl0rdψ2μ0cosψ,(11)dFMyψ=B2ψ,tl0rdψ2μ0sinψ,where Bψ,t is the distribution of magnetic density in the air gap; l0 is the effective core length of the rotor; r is the external diameter of the rotor; and ψ is the space angle.By substituting equation (7) into equations (10) and (11), the integral operation of equations (10) and (11) is completed, and the Maxwell force in the x′,y′ direction can be obtained under the rotor eccentricity caused by the local demagnetization of the permanent magnet of the PMSM:(12)FMx′t=∫02πηκl0r2μ0μ0σ02Fσ2σ0+γx′tcosψ+y′tsinψ2cos2wt−pψ−ucosψdψ,(13)FMy′t=∫02πηκl0r2μ0μ0σ02Fσ2σ0+γx′tcosψ+y′tsinψ2cos2wt−pψ−usinψdψ,where FMx′t and FMy′t are UMPs in the direction of the x′ and y′ coordinate axes in the rotating coordinate system x′O′y′. Initially, the angle between the coordinate axis of the rotating coordinate system x′O′y′ and that of stationary coordinate system xOy is zero. During the operation of the rotating coordinate system x′O′y′, the angle between the two coordinate systems is related to the rotational speed of the rotor, satisfying x=x′cosωt, y=x′sinωt, y=y′cosωt ,and x=−y′sinωt. Under stationary coordinate system xOy, the UMP is as follows:(14)FMx=FMx′tcosωt−FMy′tsinωt,(15)FMy=FMx′tsinωt+FMy′tcosωt.Note 1.
Considering the actual operation of the motor and air-gap flux density loss and other factors, set three parameters ofη, κ, γ in the UMP model (12), (13). For different types of PMSM, the three parameters η, κ, γ are different. Among them, in the UMP model (12), (13), η is the utilization coefficient of air-gap flux density; κ is the influence coefficient of demagnetization on air-gap flux density of motor, when the PMSM is normal, κ=1. The γ coefficient is related to the change of air-gap flux density caused by rotor eccentricity.Note 2.
Equations (12) and (13) are integral formulas, which shows variable Fδ. Fδ is a piecewise function in equation (6) which is related to ψ, and so equations (12) and (13) are also an integral function related to ψ. In the process of integrating equations (12) and (13), the piecewise change of the demagnetization function Fδ should be considered.Note 3.
Based on equations (12) and (13), consider the rotor speed, and perform coordinate transformation on FMx′t and FMy′t. The equations (14) and (15) are the UMP of the rotor under the directions Ox and Oy in the rectangular coordinate system xOy of the stator. Therefore, in equation (1), FMxt and FMyt are functions related to xt, yt and Fδ. In our research, the demagnetization of the permanent magnet of the PMSM results in an uneven air-gap flux density, results in unbalanced magnetic pull (UMP), which makes the rotor eccentric, and the eccentricity of the rotor further changes the air-gap flux density of the motor, resulting in UMP changing again. This further illustrates that FMxt and FMyt are functions related to xt, yt and Fδ. So, equation (1) is the relation between the demagnetization parameters Fδ and the eccentricity parameters xt,yt.
Through the combination of equation (1) with equations (12)–(15), the dynamic equation of the magnetic field-rotor coupling system under demagnetization can be obtained:(16)Mrx¨y¨+Crx˙y˙+Krxy=FMxFMy,where Mr, Cr, and Kr are the quality, damping, and stiffness matrix of the rotor system under the demagnetization. Mr=m00m, Cr=cxx00cyy, and Kr=kxx00kyy.


## 2.1. Dynamical Model of a Rotor System

The demagnetization of the permanent magnet in the rotor of the permanent magnet synchronous motor (PMSM) results in an uneven air-gap flux density in the motor, resulting in unbalanced magnetic pull (UMP) acting on the rotor, which makes the rotor eccentricity, and the eccentricity of the rotor further changes the air-gap flux density of the motor, thereby affecting the dynamic characteristics of rotor operation.Figure1 shows a schematic of rotor eccentricity under demagnetization. O and O′ are the geometric centers of the stator and rotor, respectively; εt is the eccentricity between the geometric center of the rotor and the stator; R is the radius of the rotor; δ0 is the average of air-gap length; and δψ is the length of the air-gap corresponding to angle ψ in the rotating coordinate system x′O′y′.Figure 1
The schematic diagram of rotor eccentricity under demagnetization.To establish the rotor dynamic model of a PMSM under demagnetization, this study assumes the following:(1)
All damping is considered as linear viscous damping in the calculation process, and disregard the effects of nonlinearity; the concentrated mass method is adopted without considering the axial motion of each component [17].(2)
The case where the rotor’s permanent magnet is axially uniformly demagnetized is considered. The combined amount of changes in the corresponding air-gap magnet motive force (MMF) caused by demagnetization in any axially uniform state of the permanent magnet can be described by the following two parameters: demagnetization angle and the decrease of air-gap MMF at the demagnetization angle (air-gap MMF drop percentage).On the basis of the above assumptions, the dynamic model of the magnetic field-rotor coupling system can be obtained.(1)mx¨t+cxxx˙t+kxxxt=FMxt,my¨t+cyyy˙t+kyyyt=FMyt−mg,where x,y,x˙,y˙,x¨,y¨ are the displacement, velocity, and acceleration of the rotor in the x and y axes, respectively; m is the structural quality of the rotor; cxx,cyy are the operational damping in the x and y directions; kxx,kyy are the stiffness of the rotor in the positive directions of the x and y axes; and FMxt and FMyt are the UMP in the x and y directions, respectively.


## 2.2. Modeling of Demagnetization

This section analyzesFMxt and FMyt mentioned in the previous section.


### 2.2.1. Air-Gap Flux Density

The demagnetization of permanent magnet in the rotor of permanent magnet synchronous motor (PMSM) results in unbalanced magnetic pull (UMP), which makes the rotor eccentric. Figure2 shows a cross section diagram of rotor demagnetization. With the center O′ of the rotor taken as the origin, the rotor rectangular coordinate system x′O′y′ is established. In such system, the x′-axis direction is opposite to the center line direction of the demagnetization angle, and the coordinate system rotates with the rotor. With the center of the stator O taken as the origin, the stator stationary coordinate system is constructed with the initial direction of the x′ and y′ coordinate axes in the coordinate system x′O′y′ of the rotor. Thus, the initial stator and rotor coordinate systems are parallel.Figure 2
The cross-section diagram of rotor demagnetization.The displacement variation of eccentricity between the stator and rotor of a PMSM is shown in Figure2. β is half of demagnetization angle, ψ is the space angle associated with the x′ axis in the rotating coordinate system x′O′y′, δ0 is the average air-gap length of the PMSM, and εt is the eccentricity scale of the rotor, and θt is the angle between the connecting line of the stator and rotor center O, O′ and the x axis. When the rotor runs, it satisfies(2)θt=ωt+θ0,where θ0 is the angle between the resultant force of the rotor and the x′ axis in the rotating coordinate system x′O′y′ and ω is the rotational angular velocity of the rotor.ξt=εt/δ0 is set as the eccentricity ratio of the rotor. In the rotating coordinate system x′O′y′, the rotor eccentricity caused by the demagnetization of the permanent magnet is considered, and the air-gap distribution δψ,t can be approximately expressed as follows [18]:(3)δψ,t≈δ0−εtcosψ−θ0=δ01−ξtcosψ−θ0=δ0−x′tcosψ−y′tsinψ,where x′t=εtcosθ0 and y′t=εtsinθ0 are the projection displacements of the eccentric vector in the x′O′y′ coordinate system. The permeability of the air gap is as follows:(4)Λt=μ0δψ,t≈μ0δ01+ξtcosψ−θ0=μ0δ02δ0+x′tcosψ+y′tsinψ,where μ0 is the permeability of vacuum.Considering the demagnetization of the permanent magnets of a PMSM, the fundamental wave of the air-gap MMF established by the current in the torque winding and the permanent magnet field in the rotor of the PMSM is as follows [18]:(5)Fψ,t=Fδcoswt−pψ−u,where p is the number of pole pairs of the torque winding, Fδ is the amplitude of the fundamental wave of the air-gap MMF, u is the spatial initial phase angle of the fundamental wave of the air-gap MMF, and w is the angular frequency of the winding current. The rotation angular velocity ω of the rotor and the angular frequency w of the winding current satisfy ω=60∗2∗π∗w/p.In the rotating coordinate systemx′O′y′ and in consideration of the axial direction uniform demagnetization of the rotor of a PMSM, the change of the amplitude of the air-gap MMF caused by the demagnetization of the permanent magnet is equivalently mapped to demagnetization angle β on both sides of the coordinate −x′ axis, and the decrease of air-gap MMF drops to 1−ζ∗F^ (at the demagnetization angle). At this time, the amplitude of the fundamental wave of the air-gap MMF Fδ is(6)Fδ=F^,ψ∈−π+β,π−β,1−ζF^,ψ∈−π,−π+β∪π−β,π,where β∈0,π and ζ⩽1 is the decrease of air-gap MMF at the demagnetization angle (air-gap MMF drop percentage).In the rotating coordinate systemx′O′y′, in consideration of the local demagnetization of the permanent magnet of a PMSM, the air-gap flux density in the motor is as follows:(7)Bψ,t=Λψ,tFψ,t=μ0δ02Fδδ0+x′tcosψ+y′tsinψcoswt−pψ−u.Assuming that the permanent magnet is in normal condition and the air gap is uniform, the amplitude of flux density of air-gap for pair polep is [18](8)B^=μ0δ0F^=μ0δ0324πNI2p,where N is the number of turns of each phase of the torque winding and I is the amplitude of the excitation current of the torque winding.


### 2.2.2. Analysis of UMP

The magnetic tension formed on the boundary of a magnetic material with varying permeability (such as air and core) in the magnetic circuit of a PMSM is called Maxwell force, also known as detent force; the direction of the force is perpendicular to the boundary surface of the magnetic material.When analyzing the magnetic field of the air gap of a PMSM, we assume that the magnetic density of the air gap is symmetrically distributed. At this time, the resultant force of the Maxwell force is zero. However, considering the local demagnetization of the permanent magnet in PMSM, the magnetic density distribution in the air gap of the motor is not uniform. Thus, the resultant Maxwell force is not zero.The rotor of a PMSM is subjected to the Maxwell force of the nonzero resultant force; as such, the rotor produces eccentricity, and the magnetic density in the air gap of the PMSM changes, further altering the resultant Maxwell force. The larger the eccentricity of the rotor, the greater the Maxwell force.In accordance with the Maxwell tensor method, the magnetic density in the air gap of a PMSM is set asB, and the Maxwell force acting on the dS area of the rotor surface is as follows:(9)dFM=B2dS2μ0.Then the component of Maxwell force along thex, y direction is as follows:(10)dFMxψ=B2ψ,tl0rdψ2μ0cosψ,(11)dFMyψ=B2ψ,tl0rdψ2μ0sinψ,where Bψ,t is the distribution of magnetic density in the air gap; l0 is the effective core length of the rotor; r is the external diameter of the rotor; and ψ is the space angle.By substituting equation (7) into equations (10) and (11), the integral operation of equations (10) and (11) is completed, and the Maxwell force in the x′,y′ direction can be obtained under the rotor eccentricity caused by the local demagnetization of the permanent magnet of the PMSM:(12)FMx′t=∫02πηκl0r2μ0μ0σ02Fσ2σ0+γx′tcosψ+y′tsinψ2cos2wt−pψ−ucosψdψ,(13)FMy′t=∫02πηκl0r2μ0μ0σ02Fσ2σ0+γx′tcosψ+y′tsinψ2cos2wt−pψ−usinψdψ,where FMx′t and FMy′t are UMPs in the direction of the x′ and y′ coordinate axes in the rotating coordinate system x′O′y′. Initially, the angle between the coordinate axis of the rotating coordinate system x′O′y′ and that of stationary coordinate system xOy is zero. During the operation of the rotating coordinate system x′O′y′, the angle between the two coordinate systems is related to the rotational speed of the rotor, satisfying x=x′cosωt, y=x′sinωt, y=y′cosωt ,and x=−y′sinωt. Under stationary coordinate system xOy, the UMP is as follows:(14)FMx=FMx′tcosωt−FMy′tsinωt,(15)FMy=FMx′tsinωt+FMy′tcosωt.Note 1.
Considering the actual operation of the motor and air-gap flux density loss and other factors, set three parameters ofη, κ, γ in the UMP model (12), (13). For different types of PMSM, the three parameters η, κ, γ are different. Among them, in the UMP model (12), (13), η is the utilization coefficient of air-gap flux density; κ is the influence coefficient of demagnetization on air-gap flux density of motor, when the PMSM is normal, κ=1. The γ coefficient is related to the change of air-gap flux density caused by rotor eccentricity.Note 2.
Equations (12) and (13) are integral formulas, which shows variable Fδ. Fδ is a piecewise function in equation (6) which is related to ψ, and so equations (12) and (13) are also an integral function related to ψ. In the process of integrating equations (12) and (13), the piecewise change of the demagnetization function Fδ should be considered.Note 3.
Based on equations (12) and (13), consider the rotor speed, and perform coordinate transformation on FMx′t and FMy′t. The equations (14) and (15) are the UMP of the rotor under the directions Ox and Oy in the rectangular coordinate system xOy of the stator. Therefore, in equation (1), FMxt and FMyt are functions related to xt, yt and Fδ. In our research, the demagnetization of the permanent magnet of the PMSM results in an uneven air-gap flux density, results in unbalanced magnetic pull (UMP), which makes the rotor eccentric, and the eccentricity of the rotor further changes the air-gap flux density of the motor, resulting in UMP changing again. This further illustrates that FMxt and FMyt are functions related to xt, yt and Fδ. So, equation (1) is the relation between the demagnetization parameters Fδ and the eccentricity parameters xt,yt.
Through the combination of equation (1) with equations (12)–(15), the dynamic equation of the magnetic field-rotor coupling system under demagnetization can be obtained:(16)Mrx¨y¨+Crx˙y˙+Krxy=FMxFMy,where Mr, Cr, and Kr are the quality, damping, and stiffness matrix of the rotor system under the demagnetization. Mr=m00m, Cr=cxx00cyy, and Kr=kxx00kyy.


## 2.2.1. Air-Gap Flux Density

The demagnetization of permanent magnet in the rotor of permanent magnet synchronous motor (PMSM) results in unbalanced magnetic pull (UMP), which makes the rotor eccentric. Figure2 shows a cross section diagram of rotor demagnetization. With the center O′ of the rotor taken as the origin, the rotor rectangular coordinate system x′O′y′ is established. In such system, the x′-axis direction is opposite to the center line direction of the demagnetization angle, and the coordinate system rotates with the rotor. With the center of the stator O taken as the origin, the stator stationary coordinate system is constructed with the initial direction of the x′ and y′ coordinate axes in the coordinate system x′O′y′ of the rotor. Thus, the initial stator and rotor coordinate systems are parallel.Figure 2
The cross-section diagram of rotor demagnetization.The displacement variation of eccentricity between the stator and rotor of a PMSM is shown in Figure2. β is half of demagnetization angle, ψ is the space angle associated with the x′ axis in the rotating coordinate system x′O′y′, δ0 is the average air-gap length of the PMSM, and εt is the eccentricity scale of the rotor, and θt is the angle between the connecting line of the stator and rotor center O, O′ and the x axis. When the rotor runs, it satisfies(2)θt=ωt+θ0,where θ0 is the angle between the resultant force of the rotor and the x′ axis in the rotating coordinate system x′O′y′ and ω is the rotational angular velocity of the rotor.ξt=εt/δ0 is set as the eccentricity ratio of the rotor. In the rotating coordinate system x′O′y′, the rotor eccentricity caused by the demagnetization of the permanent magnet is considered, and the air-gap distribution δψ,t can be approximately expressed as follows [18]:(3)δψ,t≈δ0−εtcosψ−θ0=δ01−ξtcosψ−θ0=δ0−x′tcosψ−y′tsinψ,where x′t=εtcosθ0 and y′t=εtsinθ0 are the projection displacements of the eccentric vector in the x′O′y′ coordinate system. The permeability of the air gap is as follows:(4)Λt=μ0δψ,t≈μ0δ01+ξtcosψ−θ0=μ0δ02δ0+x′tcosψ+y′tsinψ,where μ0 is the permeability of vacuum.Considering the demagnetization of the permanent magnets of a PMSM, the fundamental wave of the air-gap MMF established by the current in the torque winding and the permanent magnet field in the rotor of the PMSM is as follows [18]:(5)Fψ,t=Fδcoswt−pψ−u,where p is the number of pole pairs of the torque winding, Fδ is the amplitude of the fundamental wave of the air-gap MMF, u is the spatial initial phase angle of the fundamental wave of the air-gap MMF, and w is the angular frequency of the winding current. The rotation angular velocity ω of the rotor and the angular frequency w of the winding current satisfy ω=60∗2∗π∗w/p.In the rotating coordinate systemx′O′y′ and in consideration of the axial direction uniform demagnetization of the rotor of a PMSM, the change of the amplitude of the air-gap MMF caused by the demagnetization of the permanent magnet is equivalently mapped to demagnetization angle β on both sides of the coordinate −x′ axis, and the decrease of air-gap MMF drops to 1−ζ∗F^ (at the demagnetization angle). At this time, the amplitude of the fundamental wave of the air-gap MMF Fδ is(6)Fδ=F^,ψ∈−π+β,π−β,1−ζF^,ψ∈−π,−π+β∪π−β,π,where β∈0,π and ζ⩽1 is the decrease of air-gap MMF at the demagnetization angle (air-gap MMF drop percentage).In the rotating coordinate systemx′O′y′, in consideration of the local demagnetization of the permanent magnet of a PMSM, the air-gap flux density in the motor is as follows:(7)Bψ,t=Λψ,tFψ,t=μ0δ02Fδδ0+x′tcosψ+y′tsinψcoswt−pψ−u.Assuming that the permanent magnet is in normal condition and the air gap is uniform, the amplitude of flux density of air-gap for pair polep is [18](8)B^=μ0δ0F^=μ0δ0324πNI2p,where N is the number of turns of each phase of the torque winding and I is the amplitude of the excitation current of the torque winding.


## 2.2.2. Analysis of UMP

The magnetic tension formed on the boundary of a magnetic material with varying permeability (such as air and core) in the magnetic circuit of a PMSM is called Maxwell force, also known as detent force; the direction of the force is perpendicular to the boundary surface of the magnetic material.When analyzing the magnetic field of the air gap of a PMSM, we assume that the magnetic density of the air gap is symmetrically distributed. At this time, the resultant force of the Maxwell force is zero. However, considering the local demagnetization of the permanent magnet in PMSM, the magnetic density distribution in the air gap of the motor is not uniform. Thus, the resultant Maxwell force is not zero.The rotor of a PMSM is subjected to the Maxwell force of the nonzero resultant force; as such, the rotor produces eccentricity, and the magnetic density in the air gap of the PMSM changes, further altering the resultant Maxwell force. The larger the eccentricity of the rotor, the greater the Maxwell force.In accordance with the Maxwell tensor method, the magnetic density in the air gap of a PMSM is set asB, and the Maxwell force acting on the dS area of the rotor surface is as follows:(9)dFM=B2dS2μ0.Then the component of Maxwell force along thex, y direction is as follows:(10)dFMxψ=B2ψ,tl0rdψ2μ0cosψ,(11)dFMyψ=B2ψ,tl0rdψ2μ0sinψ,where Bψ,t is the distribution of magnetic density in the air gap; l0 is the effective core length of the rotor; r is the external diameter of the rotor; and ψ is the space angle.By substituting equation (7) into equations (10) and (11), the integral operation of equations (10) and (11) is completed, and the Maxwell force in the x′,y′ direction can be obtained under the rotor eccentricity caused by the local demagnetization of the permanent magnet of the PMSM:(12)FMx′t=∫02πηκl0r2μ0μ0σ02Fσ2σ0+γx′tcosψ+y′tsinψ2cos2wt−pψ−ucosψdψ,(13)FMy′t=∫02πηκl0r2μ0μ0σ02Fσ2σ0+γx′tcosψ+y′tsinψ2cos2wt−pψ−usinψdψ,where FMx′t and FMy′t are UMPs in the direction of the x′ and y′ coordinate axes in the rotating coordinate system x′O′y′. Initially, the angle between the coordinate axis of the rotating coordinate system x′O′y′ and that of stationary coordinate system xOy is zero. During the operation of the rotating coordinate system x′O′y′, the angle between the two coordinate systems is related to the rotational speed of the rotor, satisfying x=x′cosωt, y=x′sinωt, y=y′cosωt ,and x=−y′sinωt. Under stationary coordinate system xOy, the UMP is as follows:(14)FMx=FMx′tcosωt−FMy′tsinωt,(15)FMy=FMx′tsinωt+FMy′tcosωt.Note 1.
Considering the actual operation of the motor and air-gap flux density loss and other factors, set three parameters ofη, κ, γ in the UMP model (12), (13). For different types of PMSM, the three parameters η, κ, γ are different. Among them, in the UMP model (12), (13), η is the utilization coefficient of air-gap flux density; κ is the influence coefficient of demagnetization on air-gap flux density of motor, when the PMSM is normal, κ=1. The γ coefficient is related to the change of air-gap flux density caused by rotor eccentricity.Note 2.
Equations (12) and (13) are integral formulas, which shows variable Fδ. Fδ is a piecewise function in equation (6) which is related to ψ, and so equations (12) and (13) are also an integral function related to ψ. In the process of integrating equations (12) and (13), the piecewise change of the demagnetization function Fδ should be considered.Note 3.
Based on equations (12) and (13), consider the rotor speed, and perform coordinate transformation on FMx′t and FMy′t. The equations (14) and (15) are the UMP of the rotor under the directions Ox and Oy in the rectangular coordinate system xOy of the stator. Therefore, in equation (1), FMxt and FMyt are functions related to xt, yt and Fδ. In our research, the demagnetization of the permanent magnet of the PMSM results in an uneven air-gap flux density, results in unbalanced magnetic pull (UMP), which makes the rotor eccentric, and the eccentricity of the rotor further changes the air-gap flux density of the motor, resulting in UMP changing again. This further illustrates that FMxt and FMyt are functions related to xt, yt and Fδ. So, equation (1) is the relation between the demagnetization parameters Fδ and the eccentricity parameters xt,yt.
Through the combination of equation (1) with equations (12)–(15), the dynamic equation of the magnetic field-rotor coupling system under demagnetization can be obtained:(16)Mrx¨y¨+Crx˙y˙+Krxy=FMxFMy,where Mr, Cr, and Kr are the quality, damping, and stiffness matrix of the rotor system under the demagnetization. Mr=m00m, Cr=cxx00cyy, and Kr=kxx00kyy.


## 3. Simulation Calculation and Dynamic Analysis

In this section, first, we confirm the accuracy of the mathematical model for unbalanced magnetic pull (UMP) under rotor uniform demagnetization and eccentricity by using the mutual verification of Maxwell simulation and MATLAB calculation. Secondly, according to the established dynamic model of PMSMs rotor-magnetic field, the characteristics of the PMSMs rotor system under demagnetization are analyzed. Considering the actual model parameters of PMSM, the parameters are as follows (Table1).Table 1
The main parameters of PMSM.
VariableUnitValueAverage air-gap length,δ0mm1Outer radius of the rotor,rmm244Shaft diameter of the rotor,dmm52Effective core length,l0mm245Shaft length,lmm410Number of pole pairs,pPair2Phase number of motor,MpPhase3Coil turn,NTurn10Rotor quality,mkg76Stator current,IA130Shaft material40 Cr


### 3.1. Simulation and Verification of UMP Model

For the uniform demagnetization of PMSM, this section focuses on the unbalanced magnetic pull (UMP) model (12) (13) under the condition of uniform demagnetization and eccentricity. The uniform demagnetization and eccentricity of the PMSM are set by Maxwell software, and the results of simulation analysis are verified with MATLAB calculation. The verification mainly considers two aspects:(i)
The accuracy of model (12) and (13)(ii)
Considering the causes of rotor eccentricity (mainly considering whether the rotor eccentricity is completely caused by demagnetization), analyzing the changes of UMP under different eccentric states under the demagnetizationFigure3 shows the demagnetization model of PMSM by Maxwell. Demagnetization is set in the model with the negative half axis of rotor rectangular coordinate system as the center line. And Figures 4 to 7 can be obtained through corresponding parameters setting and simulation.Figure 3
The diagram of UMP trend under the change of air-gap MMF.Figure 4
The diagram of UMP trend under the change of air-gap MMF.Figure 5
The diagram of UMP trend under the change of demagnetization angle.Figure 6
The diagram of curve of UMP in directionsx,y.Figure 7
The diagram of UMP trend under the eccentricity change in directionx.In Figure4, set demagnetization angle 2∗β=π/2, the initial eccentric coordinate of the rotor is (0.2, 0.1), and set the decrease of air-gap MMF at the demagnetization angle (air-gap MMF drop percentage) ζ as a variable. The change trend of the resultant force in directions x,y of the UMP under the change of air-gap MMF can be seen that, where the red curve is calculated by MATLAB, and the blue curve is Maxwell data. By comparing the running trend of the two curves, it can be concluded that the two curves are basically consistent.In Figure5, set the decrease of air-gap MMF at the demagnetization angle (air-gap MMF drop percentage) ζ=50%, the initial eccentric coordinate of the rotor is (0.2, 0.1), and let demagnetization angle β as a variable. We can get from Figure 5 the diagram of resultant force in directions x,y of UMP trend under the change of demagnetization angle, where the red curve is calculated by MATLAB according to the mathematical model, and the blue curve is Maxwell data. By comparing the running trend of the two curves, it can be concluded that the two curves are basically consistent.In Figure6, set the decrease of air-gap MMF at the demagnetization angle (air-gap MMF drop percentage) ζ=50%, the demagnetization angle 2∗β=π/2, and the initial eccentric coordinate of the rotor is (0.2, 0.1). We can get from Figure 6 the diagram of curve of UMP in directions x,y, where the red curve is the UMP in direction x, and the blue curve is the UMP in direction y. In the figure, the scatter curve is Maxwell simulation data, and continuous curve is MATLAB calculation data. By comparing the data calculated by MATLAB based on mathematical model and Maxwell simulation data, it can be found that the two data are basically consistent.In Table2, set the decrease of air-gap MMF at the demagnetization angle (air-gap MMF drop percentage) ζ=50%, the demagnetization angle 2∗β=π/2. We can get from Table 2 the data table of UMP variation under eccentric direction change.Table 2
Data table of UMP variation under eccentric direction change.
Initial eccentric coordinateMaxwellMATLABmm (x, y)Data (N)Data (N)(−0.2, 0.1)33553202(0.2, −0.1)39803900(−0.2, −0.1)33893223(0.2, 0.1)40503950(0.225, 0)41364120(−0.225, 0)32803050(0, 0.225)36803628(0, −0.225)37103625In Figure7, set the decrease of air-gap MMF at the demagnetization angle (air-gap MMF drop percentage) ζ=50%, the demagnetization angle 2∗β=π/2. We can get from Figure 7 the diagram of resultant force in directions x,y of UMP trend under the eccentricity change in direction x.By analyzing Table2 and Figure 7, it can be concluded that when demagnetization is constant, the UMP on the rotor of PMSM will also change with the change of eccentricity direction and eccentricity.Therefore, by comparing Figures4–7 and the data in Table 2, it can be concluded that the model (12) (13) of UMP under uniform demagnetization and eccentricity can perfectly describe the change of unbalanced magnetic pull (UMP) under different demagnetization states.


### 3.2. Vibration Characteristics Analysis of Rotor-Magnetic Field System

This section mainly analyzes the influence of UMP on the dynamic characteristics of PMSM caused by axial uniform demagnetization of permanent magnet. Considering the rotor eccentricity caused by uniform demagnetization of PMSM and the unbalanced magnetic pull on the rotor, the changes of critical speed and axis orbit of the rotor system are analyzed.Considering the uniform demagnetization of rotor, the influence of demagnetization angle and decrease of air-gap MMF at the demagnetization angle on rotor dynamic characteristics is studied. In order to further study the dynamic characteristics of rotor under uniform demagnetization, the change of demagnetization angle and decrease of air-gap MMF on rotor dynamic characteristics is studied in this section. Assuming that the initial state is the ideal working state of the motor, the rotor eccentricity caused by the motor assembly is ignored; that is, the rotor coordinate systemO′x,y is consistent with the stator coordinate system Ox,y in the initial state of PMSM.In this section, the dynamic system is solved by the ode45 solver of MATLAB. Due to the UMP, the system (16) is time-varying differential equations. Set 1/300 of the rotor rotation period as a single step; in a single step, the time variables in the time-varying system (16) are linearized by Lipschitz condition, and the ordinary differential equations corresponding to the system (16) under each step are solved in turn, and finally the time-varying differential equation system is solved. Calculate 100 rotation cycles, discard the unsteady state data, taking the solution of steady-state for the next 30 cycles, and discuss the dynamic characteristics of the rotor in the speed ω=3.14−628rad/s.Figure8 shows bifurcation diagram of PMSM’s rotor system with demagnetization angle variation, setting the decrease of air-gap MMF at the demagnetization angle (air-gap MMF drop percentage) is ζ=30%, and the demagnetization angle β takes 0π/6π/3π/2.Figure 8
Bifurcation diagram of rotor system with demagnetization angle variation. (a)β=0, (b) β=π/6, (c) β=π/3, (d) β=π/2.
(a)(b)(c)(d)It can be seen from Figure8 that the PMSM is normal. When demagnetization angle is at β=0, at this time, the resultant force of UMP caused by demagnetization is zero, so the amplitude of the PMSM’s rotor is small. Compared with Figure 8, When the decrease of air-gap MMF is fixed, and demagnetization angle β∈0,π/2, with the increase of β, the nonlinear dynamic characteristics of the rotor system are enhanced, the resonance peak value of the system increases, and the critical speed decreases.Figure9 shows bifurcation diagram of rotor system with air-gap MMF drop percentage variation at demagnetization angle, setting the demagnetization angle 2∗β=2∗π/15, and the decrease of air-gap MMF at the demagnetization angle (air-gap MMF drop percentage) ζ takes 020%40%60%.Figure 9
Bifurcation diagram of rotor system with air-gap MMF drop percentage variation at demagnetization angle. (a)ζ=0, (b) ζ=20%, (c) ζ=40%, (d) ζ=60%.
(a)(b)(c)(d)It can be seen from Figure9 that the PMSM is normal. When air-gap MMF drop percentage ζ=0, at this time, the resultant force of UMP caused by demagnetization is zero, so the amplitude of the PMSM’s rotor is small. Compared with Figure 9, when the demagnetization angle is fixed, with the increase of the decrease of air-gap MMF at the demagnetization angle (air-gap MMF drop percentage) ζ, the nonlinear dynamic characteristics of the rotor system are enhanced, the resonance peak value of the system increases, and the critical speed decreases.Figure10 shows the operating state diagram of the rotor under different demagnetization states at ω=204.2033rad/s. Figure 10(a) shows the axis orbit, Poincare, phase plane, and UMP when demagnetization states ζ=60% and β=π/15 are set. Figure 10(b) shows the axis orbit, Poincare, phase plane, and UMP when demagnetization states ζ=30% and β=π/3 are set. The comparison and analysis of the above figures reveal that in the two kinds of demagnetization states the operating state of the rotor is quasiperiodic; however, the axis orbit, Poincare, phase plane, and UMP are different under two demagnetization states.Figure 10
The rotor operating state diagram under different demagnetization states atω=204.2033rad/s. (a) Demagnetization states ζ=60% and β=π/15. (b) Demagnetization states ζ=30% and β=π/3.
(a)(b)Figure11 shows the operating state diagram of the rotor under different demagnetization states at ω=471.2385rad/s. Figure 11(a) shows the axis orbit, Poincare, phase plane, and UMP when demagnetization states ζ=60% and β=π/15 are set. Figure 11(b) shows the axis orbit, Poincare, phase plane, and UMP when demagnetization states ζ=30% and β=π/3 are set. The comparison and analysis of the above figures reveal that in the two kinds of demagnetization states the operating state of the rotor is quasiperiodic; however, the axis orbit, Poincare, phase plane, and UMP are different under two demagnetization states.Figure 11
The rotor operating state diagram under different demagnetization states atω=471.2385rad/s. (a) Demagnetization states ζ=60% and β=π/15. (b) Demagnetization states ζ=30% and β=π/3.
(a)(b)By comparing Figures8–11, it can be found that the influence of different demagnetization states and speeds on the dynamic characteristics of PMSM’s rotor is different.


## 3.1. Simulation and Verification of UMP Model

For the uniform demagnetization of PMSM, this section focuses on the unbalanced magnetic pull (UMP) model (12) (13) under the condition of uniform demagnetization and eccentricity. The uniform demagnetization and eccentricity of the PMSM are set by Maxwell software, and the results of simulation analysis are verified with MATLAB calculation. The verification mainly considers two aspects:(i)
The accuracy of model (12) and (13)(ii)
Considering the causes of rotor eccentricity (mainly considering whether the rotor eccentricity is completely caused by demagnetization), analyzing the changes of UMP under different eccentric states under the demagnetizationFigure3 shows the demagnetization model of PMSM by Maxwell. Demagnetization is set in the model with the negative half axis of rotor rectangular coordinate system as the center line. And Figures 4 to 7 can be obtained through corresponding parameters setting and simulation.Figure 3
The diagram of UMP trend under the change of air-gap MMF.Figure 4
The diagram of UMP trend under the change of air-gap MMF.Figure 5
The diagram of UMP trend under the change of demagnetization angle.Figure 6
The diagram of curve of UMP in directionsx,y.Figure 7
The diagram of UMP trend under the eccentricity change in directionx.In Figure4, set demagnetization angle 2∗β=π/2, the initial eccentric coordinate of the rotor is (0.2, 0.1), and set the decrease of air-gap MMF at the demagnetization angle (air-gap MMF drop percentage) ζ as a variable. The change trend of the resultant force in directions x,y of the UMP under the change of air-gap MMF can be seen that, where the red curve is calculated by MATLAB, and the blue curve is Maxwell data. By comparing the running trend of the two curves, it can be concluded that the two curves are basically consistent.In Figure5, set the decrease of air-gap MMF at the demagnetization angle (air-gap MMF drop percentage) ζ=50%, the initial eccentric coordinate of the rotor is (0.2, 0.1), and let demagnetization angle β as a variable. We can get from Figure 5 the diagram of resultant force in directions x,y of UMP trend under the change of demagnetization angle, where the red curve is calculated by MATLAB according to the mathematical model, and the blue curve is Maxwell data. By comparing the running trend of the two curves, it can be concluded that the two curves are basically consistent.In Figure6, set the decrease of air-gap MMF at the demagnetization angle (air-gap MMF drop percentage) ζ=50%, the demagnetization angle 2∗β=π/2, and the initial eccentric coordinate of the rotor is (0.2, 0.1). We can get from Figure 6 the diagram of curve of UMP in directions x,y, where the red curve is the UMP in direction x, and the blue curve is the UMP in direction y. In the figure, the scatter curve is Maxwell simulation data, and continuous curve is MATLAB calculation data. By comparing the data calculated by MATLAB based on mathematical model and Maxwell simulation data, it can be found that the two data are basically consistent.In Table2, set the decrease of air-gap MMF at the demagnetization angle (air-gap MMF drop percentage) ζ=50%, the demagnetization angle 2∗β=π/2. We can get from Table 2 the data table of UMP variation under eccentric direction change.Table 2
Data table of UMP variation under eccentric direction change.
Initial eccentric coordinateMaxwellMATLABmm (x, y)Data (N)Data (N)(−0.2, 0.1)33553202(0.2, −0.1)39803900(−0.2, −0.1)33893223(0.2, 0.1)40503950(0.225, 0)41364120(−0.225, 0)32803050(0, 0.225)36803628(0, −0.225)37103625In Figure7, set the decrease of air-gap MMF at the demagnetization angle (air-gap MMF drop percentage) ζ=50%, the demagnetization angle 2∗β=π/2. We can get from Figure 7 the diagram of resultant force in directions x,y of UMP trend under the eccentricity change in direction x.By analyzing Table2 and Figure 7, it can be concluded that when demagnetization is constant, the UMP on the rotor of PMSM will also change with the change of eccentricity direction and eccentricity.Therefore, by comparing Figures4–7 and the data in Table 2, it can be concluded that the model (12) (13) of UMP under uniform demagnetization and eccentricity can perfectly describe the change of unbalanced magnetic pull (UMP) under different demagnetization states.


## 3.2. Vibration Characteristics Analysis of Rotor-Magnetic Field System

This section mainly analyzes the influence of UMP on the dynamic characteristics of PMSM caused by axial uniform demagnetization of permanent magnet. Considering the rotor eccentricity caused by uniform demagnetization of PMSM and the unbalanced magnetic pull on the rotor, the changes of critical speed and axis orbit of the rotor system are analyzed.Considering the uniform demagnetization of rotor, the influence of demagnetization angle and decrease of air-gap MMF at the demagnetization angle on rotor dynamic characteristics is studied. In order to further study the dynamic characteristics of rotor under uniform demagnetization, the change of demagnetization angle and decrease of air-gap MMF on rotor dynamic characteristics is studied in this section. Assuming that the initial state is the ideal working state of the motor, the rotor eccentricity caused by the motor assembly is ignored; that is, the rotor coordinate systemO′x,y is consistent with the stator coordinate system Ox,y in the initial state of PMSM.In this section, the dynamic system is solved by the ode45 solver of MATLAB. Due to the UMP, the system (16) is time-varying differential equations. Set 1/300 of the rotor rotation period as a single step; in a single step, the time variables in the time-varying system (16) are linearized by Lipschitz condition, and the ordinary differential equations corresponding to the system (16) under each step are solved in turn, and finally the time-varying differential equation system is solved. Calculate 100 rotation cycles, discard the unsteady state data, taking the solution of steady-state for the next 30 cycles, and discuss the dynamic characteristics of the rotor in the speed ω=3.14−628rad/s.Figure8 shows bifurcation diagram of PMSM’s rotor system with demagnetization angle variation, setting the decrease of air-gap MMF at the demagnetization angle (air-gap MMF drop percentage) is ζ=30%, and the demagnetization angle β takes 0π/6π/3π/2.Figure 8
Bifurcation diagram of rotor system with demagnetization angle variation. (a)β=0, (b) β=π/6, (c) β=π/3, (d) β=π/2.
(a)(b)(c)(d)It can be seen from Figure8 that the PMSM is normal. When demagnetization angle is at β=0, at this time, the resultant force of UMP caused by demagnetization is zero, so the amplitude of the PMSM’s rotor is small. Compared with Figure 8, When the decrease of air-gap MMF is fixed, and demagnetization angle β∈0,π/2, with the increase of β, the nonlinear dynamic characteristics of the rotor system are enhanced, the resonance peak value of the system increases, and the critical speed decreases.Figure9 shows bifurcation diagram of rotor system with air-gap MMF drop percentage variation at demagnetization angle, setting the demagnetization angle 2∗β=2∗π/15, and the decrease of air-gap MMF at the demagnetization angle (air-gap MMF drop percentage) ζ takes 020%40%60%.Figure 9
Bifurcation diagram of rotor system with air-gap MMF drop percentage variation at demagnetization angle. (a)ζ=0, (b) ζ=20%, (c) ζ=40%, (d) ζ=60%.
(a)(b)(c)(d)It can be seen from Figure9 that the PMSM is normal. When air-gap MMF drop percentage ζ=0, at this time, the resultant force of UMP caused by demagnetization is zero, so the amplitude of the PMSM’s rotor is small. Compared with Figure 9, when the demagnetization angle is fixed, with the increase of the decrease of air-gap MMF at the demagnetization angle (air-gap MMF drop percentage) ζ, the nonlinear dynamic characteristics of the rotor system are enhanced, the resonance peak value of the system increases, and the critical speed decreases.Figure10 shows the operating state diagram of the rotor under different demagnetization states at ω=204.2033rad/s. Figure 10(a) shows the axis orbit, Poincare, phase plane, and UMP when demagnetization states ζ=60% and β=π/15 are set. Figure 10(b) shows the axis orbit, Poincare, phase plane, and UMP when demagnetization states ζ=30% and β=π/3 are set. The comparison and analysis of the above figures reveal that in the two kinds of demagnetization states the operating state of the rotor is quasiperiodic; however, the axis orbit, Poincare, phase plane, and UMP are different under two demagnetization states.Figure 10
The rotor operating state diagram under different demagnetization states atω=204.2033rad/s. (a) Demagnetization states ζ=60% and β=π/15. (b) Demagnetization states ζ=30% and β=π/3.
(a)(b)Figure11 shows the operating state diagram of the rotor under different demagnetization states at ω=471.2385rad/s. Figure 11(a) shows the axis orbit, Poincare, phase plane, and UMP when demagnetization states ζ=60% and β=π/15 are set. Figure 11(b) shows the axis orbit, Poincare, phase plane, and UMP when demagnetization states ζ=30% and β=π/3 are set. The comparison and analysis of the above figures reveal that in the two kinds of demagnetization states the operating state of the rotor is quasiperiodic; however, the axis orbit, Poincare, phase plane, and UMP are different under two demagnetization states.Figure 11
The rotor operating state diagram under different demagnetization states atω=471.2385rad/s. (a) Demagnetization states ζ=60% and β=π/15. (b) Demagnetization states ζ=30% and β=π/3.
(a)(b)By comparing Figures8–11, it can be found that the influence of different demagnetization states and speeds on the dynamic characteristics of PMSM’s rotor is different.


## 4. Conclusion

In this paper, the rotor-magnetic field coupling dynamics Jeffcott model of permanent magnet synchronous motor (PMSM) under axial uniform demagnetization is considered. On this basis, the demagnetization status of permanent magnet in PMSM is described from two aspects of demagnetization angle and decrease of air-gap MMF at the demagnetization angle (air-gap MMF drop percentage), and the influence of demagnetization on the operation of rotor system is analyzed. This study provides a theoretical basis for the accurate fault diagnosis and vibration control of PMSMs in the future.(i)
In view of the UMP produced by the uniform demagnetization of PMSM causing the rotor eccentricity to change the air-gap flux density of the PMSM again, in this paper, the PMSM is modeled and simulated by Maxwell software. By setting the demagnetization state and eccentricity state in the model, the UMP data under different uniform demagnetization and eccentricity states are obtained by simulation. The UMP data generated by Maxwell is compared with the MATLAB data calculated based on the mathematical model under the same uniform demagnetization and eccentricity conditions in this paper, so as to determine the validity of the UMP model under uniform demagnetization and eccentricity.Meanwhile, this paper simulates by setting different eccentric states (eccentric magnitude, eccentric angle) of the rotor under the same demagnetization state. By comparing the UMP data, it can be seen that the UMP generated by the same demagnetization state and different rotor eccentricity state has great difference. Therefore, combined with the model (16) established in this paper, it can provide a theoretical basis for the judgment of demagnetization and other mixed faults of PMSM.(ii)
Based on the rotor-magnetic field coupling dynamics model of PMSM under uniform demagnetization, analyzing the influence of demagnetization of the rotor permanent magnet on the operation of the rotor system, the following conclusions can be drawn:When the decrease of air-gap MMF is fixed, and demagnetization angleβ∈0,π/2, with the increase of β, the nonlinear dynamic characteristics of the rotor system are enhanced, the resonance peak value increases, and the critical speed decreases.When the demagnetization angle is fixed, with the increase of decrease of air-gap MMF at the demagnetization angle (air-gap MMF drop percentage)ζ, the nonlinear dynamic characteristics of the rotor system are enhanced, the resonance peak value increases, and the critical speed decreases.In different demagnetization states (demagnetization angle, decrease of air-gap MMF), the axis orbit, Poincare, phase plane, and UMP have unique characteristics. So, the demagnetization state can be reflected from the axis orbit, Poincare, phase plane, and UMP, which provides a theoretical basis for the accurate fault diagnosis and vibration control of PMSMs in the future.



---

*Source: 6622851-2021-05-13.xml*