# Dynamic Behaviors and Mechanisms of Permanent Magnet Synchronous Motor with Excitation

Authors: Shaohua Liu; Tingting Shao; Hai-dong Wu; Dong-hui Zhang; Qing-zhen Han
Journal: Mathematical Problems in Engineering
 (2021)

---

## Abstract

The aim of this study is to make a general exploration of the dynamic characteristics of the permanent magnet synchronous motor (PMSM) with parametric or external perturbation. The pitchfork, fold, and Hopf bifurcations are derived by using bifurcation theory. Simulation results not only confirm the theoretical analysis results but also show the Bogdanov–Takens bifurcation of the equilibrium. Dynamic behaviors, such as period three and chaotic motion of PMSM, are analyzed by using bifurcation diagram and phase portraits. The symmetric fold/fold bursting oscillation as well as two kinds of delayed pitchfork bursting oscillations is obtained, and different mechanisms are presented.


---

## Body

## 1. Introduction

Since the permanent magnet synchronous motor (PMSM) has low inertia, low noise, high power density, and high efficiency, it is one of most important development trends in a variety of driving motors [1–3]. PMSM is a multivariable, strongly nonlinear, and strongly coupled electromechanical system which is widely used in industry and may exhibit different dynamic behaviors as the parameters and the external load are not fixed but may fluctuate within a range according to the working environment. So, the dynamic behaviors of the PMSM with parametric perturbation as well as with external perturbation should be discussed deeply to provide theoretical basis.Scholars have done a lot of research on the dynamic characteristics of PMSM. Li et al. [4] derived a dynamic model of PMSM by applying an affine transformation and a time-scaling transformation and studied the complicated dynamic behaviors such as limit cycle and chaos in the smooth air gap PMSM. Jing et al. [5] studied the dynamics of nonsmooth air gap PMSM and derived the bifurcation conditions of pitchfork bifurcation and Hopf bifurcation by using bifurcation theory and center manifold theory. The dynamic behaviors such as the period doubling bifurcation, cyclic fold bifurcation, and chaotic motion are obtained by using numerical simulation. Rasoolzadeh and Tavazoei [6] discussed the dynamic behaviors of PMSM with external excitation and analyzed the effect of excitation angle-frequency on the dynamic behaviors by using bifurcation diagram. Wang et al. [7] analyzed the dynamic response of PMSM under FOC control and obtained different attractors with the change of single parameter. In reference [8], the parameter region is divided into the chaotic area and nonchaotic area by a lot of simulation results, and the maximum torque output of the stable region is estimated based on the nonchaotic area.Bursting oscillation is a kind of fast-slow dynamic behavior when the dynamic system has two or more time scales [9–12]. At present, codimension one bursting oscillations of the dynamic system are studied deeply in the nervous system. Izhikevich [13] summarized the codimension one bursting oscillations that may happened in the nervous system and studied its classification and naming method, and some of them have already been discovered in the nervous system. Bursting oscillation is also studied in other fields. In reference [14], bursting behaviors of a piecewise mechanical system with parameter perturbation in stiffness is studied. Wu et al. [15] discussed the chaotic and periodic bursting phenomena in a memristive Wien-bridge oscillator. The time-delayed bursting oscillations in a two dimensional oscillator is studied in reference [16]. Makouo and Woafo [12] studied the bursting patterns in Van der Pol oscillators by experimental observation. Han [17] discussed the fold/fold bursting oscillation and fold/subHopf bursting oscillation in a hyperchaotic Lorenz system. The high speed spiking and bursting oscillations in a long-delayed broadband optoelectronic oscillator is studied in reference [18]. Bi et al. [19] discussed the mechanism of bursting oscillations with different codimensional bifurcations and nonlinear structures. Analysis of the above literature shows that domestic and foreign scholars have done a lot of research on the fast-slow dynamics, but there are still few reports on PMSM which need to be discussed for the deep exploration of the dynamics of PMSM.This paper gives a general exploration of the dynamic characteristics of PMSM with parametric excitation or load excitation. PMSM parameters are not limited to a special case. The stability conditions and the bifurcation conditions of the equilibrium have been derived. The local dynamic characteristics of the equilibrium points and the global dynamic characteristics of PMSM are discussed by bifurcation set and bifurcation diagram, respectively. The dynamic behaviors at some specific parameters are studied by phase portraits and Poincaré sections, and three kinds of multiple time scale dynamic behaviors are obtained.


## 2. Dynamic Model of PMSM

The dynamic model of PMSM can be described as a three-dimensional nonlinear system [20]:(1)diddt=ud−R1id+ωLqiqLd,diqdt=uq−R1iq−ωLdid+ωψfLq,dωdt=npψfiq+npLd−Lqidiq−TL−βωJ,where id,iq,ω represent the d and q axis stator current and motor angular speed, respectively; ud and uq are the d-q axis stator voltage components, respectively, Ld and Lq are the d-q axis stator inductors, respectively, R1,ψf,β,J,np, and TL are the stator winding resistance, the permanent-magnet flux, the viscous damping coefficient, the polar moment of inertia, the number of pole-pairs, and the external load torque, respectively.By applying an affine transformation and a time-scaling transformation, system (1) can be transformed into the following dimensionless model:(2)di˜ddt˜=−LqLdi˜d+ω˜i˜q+u˜d,di˜qdt˜=−i˜q−ω˜i˜d+γω˜+u˜q,dω˜dt˜=σi˜q−ω˜+εi˜di˜q−T˜L,where γ=npψf2/R1β; σ=Lqβ/R1J; u˜q=npLqψfuq/R12β; u˜d=npLqψfud/R12β; ε=Lqβ2Ld−Lq/LdJnpψf2; and T˜L=Lq2TL/R12J.


## 3. Equilibrium Point and Bifurcations

### 3.1. The Inherent Characteristics of the Smooth-Air-Gap PMSM, Namely,Ld=Lq, u˜d=0, u˜q=0, and T˜L=0

System (2) can be written as(3)di˜ddt˜=−i˜d+ω˜i˜q,di˜qdt˜=−i˜q−ω˜i˜d+γω˜,dω˜dt˜=σi˜q−ω˜.There are only two combined parameters,γ and σ, in (3), and each parameter represents a group of parameters of the actual PMSM. These two parameters are assumed to be totally independent of each other for further discussing the dynamic characteristics, and the analytical results obtained from the dynamics analysis procedures will encompass a general class of machines rather than being limited to a specific machine with given parameters. The actual dynamic behaviors can be obtained by substituting the actual parameters into γ and σ.By the analysis of the equilibrium points,i˜d0,i˜q0,ω˜0, of system (3), the equilibrium points equations are obtained:(4)−i˜d0+ω˜0i˜q0=0,−i˜q0−ω˜0i˜d0+γω˜0=0,σi˜q0−ω˜0=0,⇒i˜d0=ω˜02,i˜q0=ω˜0,ω˜0ω˜02−γ+1=0.Solving the equilibrium points equations, we obtain the following:(1)
Ifγ>1, system (3) would have three equilibrium points, including 0,0,0, γ−1,γ−1,γ−1, and γ−1,−γ−1,−γ−1(2)
Ifγ≤1, system (3) would only have one equilibrium point, 0,0,0The local stability of the equilibrium point is determined by the roots of the characteristic equation, and the equilibrium point is stable if the all the roots of the characteristic equation have negative real parts. The characteristic equation can be expressed as det (λI − J) = 0, where λ, I, and J represent the eigenvalue, identity matrix, and Jacobian matrix, respectively. In fact, the stability of the equilibrium point is usually judged using Routh–Hurwitz stability criterion [21, 22] as solving the characteristic equation is very difficult for high-dimensional dynamic system.The Jacobian matrix at the equilibrium points of system (3) is(5)J=−1ω˜0i˜q0−ω˜0−1−i˜d0+γ0σ−σ.Substituting the Jacobian matrix and equilibrium,ω˜02,ω˜0,ω˜0, into the characteristic equation, detλI−J=0, we obtain(6)λ3+σ+2λ2+1+σω˜02−σγ+ω˜02+2σλ+3σω˜02−σγ+σ=0.According to Routh–Hurwitz stability criterion, the sufficient condition of the local stability of the equilibrium points is(7)σ>−2,3σω˜02−γσ+σ>0,σ+21+σω˜02−σγ+ω˜02+2σ−3ω˜02+γ−1σ>0.By the bifurcation theory [23–25] of the equilibrium points, the equilibrium points may lose stability, and the bifurcation behavior occurs when the parameters pass through the key values. In order to obtain the conditions of Hopf bifurcation, setting λ = nj (n ≠ 0) and substituting it into characteristic equation (6), we obtain(8)−n3j−σ+2n2+1+σω˜02−σγ+ω˜02+2σnj+3σω˜02−σγ+σ=0.Equating real and imaginary parts of the equation,(9)n2=1+σω˜02−σγ+ω˜02+2σ,n2=3σω˜02−σγ+σσ+2.Then, if Hopf bifurcation occurs at the equilibrium point of (3), the following conditions must be satisfied:(10)σ+2>0,1+σω˜02−σγ+ω˜02+2σ>0,−σ2−2ω02+γ−2σ2+γ−4σ−2=0.In addition, ifγ > 1, system (3) would have two stable equilibrium points and one unstable equilibrium point, and if γ < 1, system (3) would only have one stable equilibrium point, the pitchfork bifurcation occurs at γ=1.To verify the theoretical analysis results, the double-parameter bifurcation curves of system (3) in the γ-σ plane are obtained by continuation method (see Figure 1), where “BP” and “subH” represent pitchfork bifurcation curve and subcritical Hopf bifurcation curve, respectively. The pitchfork bifurcation curve together with subcritical Hopf bifurcation curve divide the equilibrium parameter plane into three regions, namely, “A,” “B,” and “C,” where the equilibrium has different characteristics in different regions. There is only one stable equilibrium point in “A” and two stable equilibrium points and one unstable equilibrium point in “B” and three unstable equilibrium points in “C.” To verify the above research, Table 1 lists a set of parameters of PMSM [26], and the multiequilibrium coexistence phenomenon (see Figure 2), is obtained by simulation as the parameters of the PMSM are fall in region “B.”Figure 1
Bifurcation set inγ-σ plane.Table 1
Parameters of the PMSM [26].
ParameterValueUnitLd/Lq1.01mHΨR0.06784Nm/AR10.24ΩJ4.8 × 10−5kg m2pn4β0.01619N/rad/sFigure 2
Multiequilibrium coexistence phenomenon.


### 3.2. When Considering the External Load Namely,Ld=Lq, u˜d=0, u˜q=0, and T˜L≠0

This is a special case that the control inputs of the system are removed after the motor runs for a period of operation. Letu=T˜L, then system (2) becomes(11)di˜ddt˜=−i˜d+ω˜i˜q,di˜qdt˜=−i˜q−ω˜i˜d+γω˜,dω˜dt˜=σi˜q−ω˜−u.By analyzing the equilibrium,i˜d0,i˜q0,ω˜0, of system (11), we obtain(12)−i˜d0+ω˜0i˜q0=0,−i˜q0−ω˜0i˜d0+γω˜0=0,σi˜q0−ω˜0−u=0,⇒i˜d0=ω˜02+uω˜0σ,i˜q0=ω˜0+uσ,ω˜03+uσω02+1−γω˜0−uσ=0.According to the relationship amongi˜d0,i˜q0,ω˜0, the equilibrium point of (11) can be written as ω˜02+uω˜0/σ,ω˜0+u/σ,ω˜0.According to the Routh–Hurwitz criteria, the local stability condition of the equilibrium can be written as(13)σ+2>0,3σω˜02+uω˜0+u−γ+1σ>0,σ2+2σ+σ+2ω˜02+uω˜0−γσ+1+2>0.If fold bifurcation occurs at the equilibrium of (11), the following conditions must be satisfied:(14)ω˜03+uσω02+1−γω˜0+uσ=0,3σω˜02+uω˜0+u−γ+1σ=0.If Hopf bifurcation occurs at the equilibrium of (11), the following conditions must be satisfied:(15)σ+2>0,1+σω˜02+uω˜0+2−γσ+1>0,ω˜03+uσω02+1−γω˜0+uσ=0,σ2+2σ+σ+2ω˜02+uω˜0−γσ+1+2=0.The bifurcation set of equilibrium of (11) (see Figure 3) describes the bifurcation and stability characteristics of the equilibrium, where the parameter γ is fixed as 4.378 obtained by calculating the data from Table 1, and the PMSM itself always has one unstable and two stable equilibrium as shown in Figure 1. “BT” represents the Bogdanov–Takens, one kind of codimension two bifurcation of the equilibrium. The fold bifurcation curve, labelled as “LP,” together with subcritical Hopf bifurcation curve, labelled as “subH,” divided the equilibrium parameter plane into three regions, namely, “A”–“C.” There is only one stable equilibrium in “A,” and the characteristics of this stable equilibrium are not changed as well as two unstable equilibrium appear when the parameters pass fold bifurcation curve from “A” to “B.” There are one stable and two unstable equilibrium points in “B.” When the parameters pass through subcritical Hopf bifurcation line from “B” to “C,” one of the unstable equilibriums become stable, namely, there are two stable and one unstable equilibrium in “C.”Figure 3
Bifurcation set in theu-σ plane.


## 3.1. The Inherent Characteristics of the Smooth-Air-Gap PMSM, Namely,Ld=Lq, u˜d=0, u˜q=0, and T˜L=0

System (2) can be written as(3)di˜ddt˜=−i˜d+ω˜i˜q,di˜qdt˜=−i˜q−ω˜i˜d+γω˜,dω˜dt˜=σi˜q−ω˜.There are only two combined parameters,γ and σ, in (3), and each parameter represents a group of parameters of the actual PMSM. These two parameters are assumed to be totally independent of each other for further discussing the dynamic characteristics, and the analytical results obtained from the dynamics analysis procedures will encompass a general class of machines rather than being limited to a specific machine with given parameters. The actual dynamic behaviors can be obtained by substituting the actual parameters into γ and σ.By the analysis of the equilibrium points,i˜d0,i˜q0,ω˜0, of system (3), the equilibrium points equations are obtained:(4)−i˜d0+ω˜0i˜q0=0,−i˜q0−ω˜0i˜d0+γω˜0=0,σi˜q0−ω˜0=0,⇒i˜d0=ω˜02,i˜q0=ω˜0,ω˜0ω˜02−γ+1=0.Solving the equilibrium points equations, we obtain the following:(1)
Ifγ>1, system (3) would have three equilibrium points, including 0,0,0, γ−1,γ−1,γ−1, and γ−1,−γ−1,−γ−1(2)
Ifγ≤1, system (3) would only have one equilibrium point, 0,0,0The local stability of the equilibrium point is determined by the roots of the characteristic equation, and the equilibrium point is stable if the all the roots of the characteristic equation have negative real parts. The characteristic equation can be expressed as det (λI − J) = 0, where λ, I, and J represent the eigenvalue, identity matrix, and Jacobian matrix, respectively. In fact, the stability of the equilibrium point is usually judged using Routh–Hurwitz stability criterion [21, 22] as solving the characteristic equation is very difficult for high-dimensional dynamic system.The Jacobian matrix at the equilibrium points of system (3) is(5)J=−1ω˜0i˜q0−ω˜0−1−i˜d0+γ0σ−σ.Substituting the Jacobian matrix and equilibrium,ω˜02,ω˜0,ω˜0, into the characteristic equation, detλI−J=0, we obtain(6)λ3+σ+2λ2+1+σω˜02−σγ+ω˜02+2σλ+3σω˜02−σγ+σ=0.According to Routh–Hurwitz stability criterion, the sufficient condition of the local stability of the equilibrium points is(7)σ>−2,3σω˜02−γσ+σ>0,σ+21+σω˜02−σγ+ω˜02+2σ−3ω˜02+γ−1σ>0.By the bifurcation theory [23–25] of the equilibrium points, the equilibrium points may lose stability, and the bifurcation behavior occurs when the parameters pass through the key values. In order to obtain the conditions of Hopf bifurcation, setting λ = nj (n ≠ 0) and substituting it into characteristic equation (6), we obtain(8)−n3j−σ+2n2+1+σω˜02−σγ+ω˜02+2σnj+3σω˜02−σγ+σ=0.Equating real and imaginary parts of the equation,(9)n2=1+σω˜02−σγ+ω˜02+2σ,n2=3σω˜02−σγ+σσ+2.Then, if Hopf bifurcation occurs at the equilibrium point of (3), the following conditions must be satisfied:(10)σ+2>0,1+σω˜02−σγ+ω˜02+2σ>0,−σ2−2ω02+γ−2σ2+γ−4σ−2=0.In addition, ifγ > 1, system (3) would have two stable equilibrium points and one unstable equilibrium point, and if γ < 1, system (3) would only have one stable equilibrium point, the pitchfork bifurcation occurs at γ=1.To verify the theoretical analysis results, the double-parameter bifurcation curves of system (3) in the γ-σ plane are obtained by continuation method (see Figure 1), where “BP” and “subH” represent pitchfork bifurcation curve and subcritical Hopf bifurcation curve, respectively. The pitchfork bifurcation curve together with subcritical Hopf bifurcation curve divide the equilibrium parameter plane into three regions, namely, “A,” “B,” and “C,” where the equilibrium has different characteristics in different regions. There is only one stable equilibrium point in “A” and two stable equilibrium points and one unstable equilibrium point in “B” and three unstable equilibrium points in “C.” To verify the above research, Table 1 lists a set of parameters of PMSM [26], and the multiequilibrium coexistence phenomenon (see Figure 2), is obtained by simulation as the parameters of the PMSM are fall in region “B.”Figure 1
Bifurcation set inγ-σ plane.Table 1
Parameters of the PMSM [26].
ParameterValueUnitLd/Lq1.01mHΨR0.06784Nm/AR10.24ΩJ4.8 × 10−5kg m2pn4β0.01619N/rad/sFigure 2
Multiequilibrium coexistence phenomenon.


## 3.2. When Considering the External Load Namely,Ld=Lq, u˜d=0, u˜q=0, and T˜L≠0

This is a special case that the control inputs of the system are removed after the motor runs for a period of operation. Letu=T˜L, then system (2) becomes(11)di˜ddt˜=−i˜d+ω˜i˜q,di˜qdt˜=−i˜q−ω˜i˜d+γω˜,dω˜dt˜=σi˜q−ω˜−u.By analyzing the equilibrium,i˜d0,i˜q0,ω˜0, of system (11), we obtain(12)−i˜d0+ω˜0i˜q0=0,−i˜q0−ω˜0i˜d0+γω˜0=0,σi˜q0−ω˜0−u=0,⇒i˜d0=ω˜02+uω˜0σ,i˜q0=ω˜0+uσ,ω˜03+uσω02+1−γω˜0−uσ=0.According to the relationship amongi˜d0,i˜q0,ω˜0, the equilibrium point of (11) can be written as ω˜02+uω˜0/σ,ω˜0+u/σ,ω˜0.According to the Routh–Hurwitz criteria, the local stability condition of the equilibrium can be written as(13)σ+2>0,3σω˜02+uω˜0+u−γ+1σ>0,σ2+2σ+σ+2ω˜02+uω˜0−γσ+1+2>0.If fold bifurcation occurs at the equilibrium of (11), the following conditions must be satisfied:(14)ω˜03+uσω02+1−γω˜0+uσ=0,3σω˜02+uω˜0+u−γ+1σ=0.If Hopf bifurcation occurs at the equilibrium of (11), the following conditions must be satisfied:(15)σ+2>0,1+σω˜02+uω˜0+2−γσ+1>0,ω˜03+uσω02+1−γω˜0+uσ=0,σ2+2σ+σ+2ω˜02+uω˜0−γσ+1+2=0.The bifurcation set of equilibrium of (11) (see Figure 3) describes the bifurcation and stability characteristics of the equilibrium, where the parameter γ is fixed as 4.378 obtained by calculating the data from Table 1, and the PMSM itself always has one unstable and two stable equilibrium as shown in Figure 1. “BT” represents the Bogdanov–Takens, one kind of codimension two bifurcation of the equilibrium. The fold bifurcation curve, labelled as “LP,” together with subcritical Hopf bifurcation curve, labelled as “subH,” divided the equilibrium parameter plane into three regions, namely, “A”–“C.” There is only one stable equilibrium in “A,” and the characteristics of this stable equilibrium are not changed as well as two unstable equilibrium appear when the parameters pass fold bifurcation curve from “A” to “B.” There are one stable and two unstable equilibrium points in “B.” When the parameters pass through subcritical Hopf bifurcation line from “B” to “C,” one of the unstable equilibriums become stable, namely, there are two stable and one unstable equilibrium in “C.”Figure 3
Bifurcation set in theu-σ plane.


## 4. Simulation Results

The dynamics of the PMSM model with parametric or external excitation will be discussed in this section. Based on the different order of magnitude of the excitation frequency, the dynamic behaviors are divided into two parts. The bifurcation and chaos are discussed in4.1, and the fast-slow dynamics are discussed in 4.2. For the sake of simplicity, let x=i˜d, y=i˜q, and z=ω˜.


### 4.1. Bifurcation and Chaos

In order to discuss the dynamic characteristics of the PMSM with parametric and external excitation, letγ=γ0+Asinw′t, u=Asinw′t, A = 4, and using the inherent parameters in Table 1. The bifurcation diagram (see Figure 4) is drawn which describes the global nonlinear dynamic behaviors of the PMSM. The period motion and chaotic motion can be easily identified in this figure. For further analysis of the dynamic characteristics, the phase portraits and poincaré sections (see Figures 5–8) of different frequencies are drawn where the blue curve and the red dots represent the orbit and the poincaré section, respectively. It can be seen that, with the change of frequency, the PMSM can exhibit different dynamic characteristics, such as the period one motion (see Figure 5), period two motion (see Figure 6), period four motion (see Figure 7), and period three motion (see Figure 8).Figure 4
The Poincaré bifurcation diagram.Figure 5
Period one motion,w′=1.41.Figure 6
Period two motion,w′=1.425.Figure 7
Period four motion,w′=1.429.Figure 8
Period three motion,w′=2.195.


### 4.2. Fast-Slow Dynamics

An interesting two time scale fast-slow dynamic [14, 16, 17] behavior (see Figure 9(a)) of PMSM with load excitation is obtained by using the inherent parameters in Table 1 and letting u=Asinw′t, and the corresponding time history map is given in Figure 9(b). In order to analyze the mechanism of these phenomena, the equilibrium curve and the overlap of transformed phase portraits (TPP) and equilibrium curve are drawn in Figures 9(c) and 9(d), respectively.Figure 9
Symmetric fold/fold bursting. (a) Phase portrait. (b) Time history. (c) Equilibrium curve. (d) Overlap of TPP and equilibrium curve.
(a)(b)(c)(d)Suppose the trajectory starts from “A” of the up piece of equilibrium curve “E+” in Figure 9(d), where the red curve is the TPP and the blue curve is the equilibrium curve, and moves strictly until the fold bifurcation point “LP1” and then the trajectory falls into the attraction domain of low piece of equilibrium curve “E−” as the stable equilibrium curve “E+” disappears. Similar dynamics take place at the fold bifurcation point “LP2.” A cycle is completed when the trajectory is back to the initial point, and this dynamic behavior is caused by the fold bifurcation of the equilibrium curve and can be named as “symmetric fold/fold bursting.”When considering about the parametric excitation, i.e.,γ = γ0 + Asin (w′t) and u=0, two kinds of delayed pitchfork bursting oscillations are obtained (see Figures 10 and 11). In order to analyze the mechanisms of these two kinds of dynamic behaviors, the overlap of TPP and equilibrium curve is drawn in Figures 10(b) and 11(b), where the red and blue curve represent the TPP and equilibrium curve, respectively. For the delayed pitchfork bursting I (see Figure 10(b)), suppose the orbit starts from point “A” and moves to “B” point as indicated by the arrow, the orbit is not jumped from the up piece of the equilibrium curve to the middle piece of the equilibrium curve at “BP” point because of the delay of pitchfork bifurcation. When the orbit moves from “B” to “C,” for the same reason, the orbit is not jumped from the middle piece of the equilibrium curve to the up piece of the equilibrium curve at point “BP,” and the orbit is back to point “A” finally. As caused by the delay of pitchfork bifurcation, this dynamic behavior can be named as “delayed pitchfork bursting.” For the delayed pitchfork bursting II, the orbit may be attracted by the up piece or the down piece of the equilibrium curve, which is different from the delayed pitchfork bursting I.Figure 10
Delayed pitchfork bursting I;w′=0.03, σ = 10, and γ = 1.5. (a) Phase portrait. (b) Overlap of TPP and equilibrium curve.
(a)(b)Figure 11
Delayed pitchfork bursting II;w′=0.02, σ = 10, and γ0 = 1.5. (a) Phase portrait. (b) Overlap of TPP and equilibrium curve.
(a)(b)


## 4.1. Bifurcation and Chaos

In order to discuss the dynamic characteristics of the PMSM with parametric and external excitation, letγ=γ0+Asinw′t, u=Asinw′t, A = 4, and using the inherent parameters in Table 1. The bifurcation diagram (see Figure 4) is drawn which describes the global nonlinear dynamic behaviors of the PMSM. The period motion and chaotic motion can be easily identified in this figure. For further analysis of the dynamic characteristics, the phase portraits and poincaré sections (see Figures 5–8) of different frequencies are drawn where the blue curve and the red dots represent the orbit and the poincaré section, respectively. It can be seen that, with the change of frequency, the PMSM can exhibit different dynamic characteristics, such as the period one motion (see Figure 5), period two motion (see Figure 6), period four motion (see Figure 7), and period three motion (see Figure 8).Figure 4
The Poincaré bifurcation diagram.Figure 5
Period one motion,w′=1.41.Figure 6
Period two motion,w′=1.425.Figure 7
Period four motion,w′=1.429.Figure 8
Period three motion,w′=2.195.


## 4.2. Fast-Slow Dynamics

An interesting two time scale fast-slow dynamic [14, 16, 17] behavior (see Figure 9(a)) of PMSM with load excitation is obtained by using the inherent parameters in Table 1 and letting u=Asinw′t, and the corresponding time history map is given in Figure 9(b). In order to analyze the mechanism of these phenomena, the equilibrium curve and the overlap of transformed phase portraits (TPP) and equilibrium curve are drawn in Figures 9(c) and 9(d), respectively.Figure 9
Symmetric fold/fold bursting. (a) Phase portrait. (b) Time history. (c) Equilibrium curve. (d) Overlap of TPP and equilibrium curve.
(a)(b)(c)(d)Suppose the trajectory starts from “A” of the up piece of equilibrium curve “E+” in Figure 9(d), where the red curve is the TPP and the blue curve is the equilibrium curve, and moves strictly until the fold bifurcation point “LP1” and then the trajectory falls into the attraction domain of low piece of equilibrium curve “E−” as the stable equilibrium curve “E+” disappears. Similar dynamics take place at the fold bifurcation point “LP2.” A cycle is completed when the trajectory is back to the initial point, and this dynamic behavior is caused by the fold bifurcation of the equilibrium curve and can be named as “symmetric fold/fold bursting.”When considering about the parametric excitation, i.e.,γ = γ0 + Asin (w′t) and u=0, two kinds of delayed pitchfork bursting oscillations are obtained (see Figures 10 and 11). In order to analyze the mechanisms of these two kinds of dynamic behaviors, the overlap of TPP and equilibrium curve is drawn in Figures 10(b) and 11(b), where the red and blue curve represent the TPP and equilibrium curve, respectively. For the delayed pitchfork bursting I (see Figure 10(b)), suppose the orbit starts from point “A” and moves to “B” point as indicated by the arrow, the orbit is not jumped from the up piece of the equilibrium curve to the middle piece of the equilibrium curve at “BP” point because of the delay of pitchfork bifurcation. When the orbit moves from “B” to “C,” for the same reason, the orbit is not jumped from the middle piece of the equilibrium curve to the up piece of the equilibrium curve at point “BP,” and the orbit is back to point “A” finally. As caused by the delay of pitchfork bifurcation, this dynamic behavior can be named as “delayed pitchfork bursting.” For the delayed pitchfork bursting II, the orbit may be attracted by the up piece or the down piece of the equilibrium curve, which is different from the delayed pitchfork bursting I.Figure 10
Delayed pitchfork bursting I;w′=0.03, σ = 10, and γ = 1.5. (a) Phase portrait. (b) Overlap of TPP and equilibrium curve.
(a)(b)Figure 11
Delayed pitchfork bursting II;w′=0.02, σ = 10, and γ0 = 1.5. (a) Phase portrait. (b) Overlap of TPP and equilibrium curve.
(a)(b)


## 5. Conclusion

The nonlinear dynamic behaviors of the PMSM with parametric or load perturbation are generally studied. The primary conclusions of this paper are as follows:(1)
Different kinds of bifurcations, such as pitchfork bifurcation, fold bifurcation, and Hopf bifurcation, are deduced. The bifurcation curves are also simulated.(2)
The nonlinear dynamics of the PMSM with both parametric excitation and load excitation are simulated. Different kinds of motions, such as period one, period two, and period four, are obtained. The period three motion is also obtained which can prove the existence of chaotic motion.(3)
When the frequency of the excitation is superlow, the fast-slow dynamics are obtained. Simulation results show three kinds of fast-slow dynamics: symmetric fold/fold bursting oscillation of the PMSM with external perturbation, two kinds of delayed pitchfork bursting oscillations of the PMSM with parametric perturbation. The mechanisms of those three kinds of fast-slow dynamics are also analyzed by using bifurcation curve of the equilibrium.



---

*Source: 5535856-2021-06-04.xml*