A permanent magnet synchronous motor (PMSM) model with smooth air gap and an exogenous periodic input is introduced and analyzed in this paper. With a simple mathematical transformation, a new nonautonomous Lorenz-like system is derived from this PMSM model, and this new three-dimensional system can display the complicated dynamics such as the chaotic attractor and the multiperiodic orbits by adjusting the frequency and amplitude of the exogenous periodic inputs. Moreover, this new system shows a double-deck chaotic attractor that is completely different from the four-wing chaotic attractors on topological structures, although the phase portrait shapes of the new attractor and the four-wing chaotic attractors are similar. The exotic phenomenon has been well demonstrated and investigated by numerical simulations, bifurcation analysis, and electronic circuit implementation.
1. Introduction
Due to simple structure, low manufacturing cost, and high performance [1], permanent magnet synchronous motors (PMSMs) are widely used in modern direct drive motor systems [2, 3], power conversion [4, 5], computer disk drives [6, 7], and domestic applications [8, 9]. Compared with DC motors and induction motors, PMSMs have the advantages of high speed operational capability, precise torque control, high torque to current ratio, high power to weight ratio, high efficiency, and high robustness to exogenous disturbances [10]. Recently, the studies about the secure and stable operation of PMSMs, which is an essential requirement of industrial automation manufacturing, have attracted more and more attentions. Some research results showed that PMSMs are chaotic when system parameters change within certain ranges [11–13]. The chaotic behavior in PMSMs, which appears mainly as intermittent ripples of torque, low-frequency oscillations of rotational speed of motor, can extremely destroy the stability of the motor and even induce drive system collapse.
Because the mathematical model for a PMSM is a typical nonlinear system and it cannot be analyzed based on the linear theory, the modern nonlinear theory such as bifurcation and chaos, which has been widely used to study the stability of nonlinear systems, can be used to study the permanent magnet machines. Li et al. deeply analyzed the bifurcation and chaos in PMSM with smooth air gap in various cases that stator voltages on rotating frame and external load torque are equal to zero in whole or part [11]. However, their studies about the nonzero external inputs of PMSM are only limited to DC power supply, without considering AC power supply. For the PMSM with nonsmooth air gap, Jing et al. studied the complex dynamics including the stability, the pitchfork, and Hopf bifurcations using the bifurcation theory and the center manifold [14]. In order to avoid chaos in PMSM, Choi gave an adaptive control scheme for the permanent magnets based on the bifurcation analysis [15, 16]. Therefore, the above research results [11, 14–16] show that the permanent magnet machines can experience chaotic oscillations at some situations.
Because the chaotic dynamical behavior has been widely observed in practice, how to apply chaos to the engineering practice and generate strange chaotic attractors with complex topological structure by a simple dynamical system has become the focus in the field of nonlinear dynamics recently. Moreover, sometimes it is necessary to avoid the chaotic dynamics in some systems. Hence, it would be useful to investigate the dynamical behavior of PMSM and find out whether there are strange attractors. In this paper, the focuses of this study are as follows:
investigating the dynamical behavior of a PMSM model with smooth air gap with external periodic excitation to the q-axes stator on rotating frame, which has not been carefully studied;
analyzing the bifurcation and chaos of the PMSM model by phase portrait, Poincare map, Lyapunov exponents, and bifurcation diagram, which are important tools to find strange double-deck butterfly chaotic attractors and multiperiodic orbits and analyze the dynamical behavior of the PMSM model;
analyzing the power spectra of the PMSM model since the power spectrum analysis is an effective tool to explore the frequency components of state variable;
besides, validating the findings by the simulations and experimental circuits.
This paper is organized as follows. In Section 2, the mathematical model for a PMSM with external excitations is proposed. Section 3 focuses on the bifurcation and chaos in the nonautonomous PMSM system by the methods of phase portrait, Poincare map, Lyapunov exponents, bifurcation diagram, and power spectrum. In Section 4, an experimental circuit is implemented to verify the obtained numerical results. Finally, we draw the conclusions.
2. The Mathematical Model
The d-q machine model for a permanent magnet synchronous motor (PMSM) in the rotor rotating reference frame [11, 17] is characterized by
(1)diddt=1Ld(-Rid+Lqϖiq+υd)diqdt=1Lq(-Riq-Ldidϖ-ψrϖ+υq)dϖdt=1J[nψriq+(Ld-Lq)idiq-Bϖ-TL],
where id, iq, and ϖ are the state variables of the nonlinear dynamical system (1). The state variables id, iq, and ϖ and system parameters Ld, Lq, R, J, B, TL, n, and ψr represent the physical meanings, respectively, as follows:
id, iq: d- and q-axes stator currents on rotating frame;
vd, vq: d- and q-axes stator voltages on rotating frame;
Ld, Lq: d- and q-axes stator self-inductances;
ϖ: rotor speed at electrical angle;
R: stator per-phase residence;
J: polar moment of inertia;
B: viscous damping coefficient;
TL: external load torque;
n: number of pole-pairs;
ψr: permanent magnetic flux.
If Ld=Lq=L, the system (1) is the model of smooth air-gap synchronous machines. Applying a time-scaling transformation t=τt~, where parameter τ=L/R is a scaling factor, and further letting
(2)i~d=τnψrBid,i~q=τnψrBiq,ϖ~=τϖ,
the system (1) can be rewritten as
(3)di~ddt~=-i~d+i~qϖ~+υ~ddi~qdt~=-i~q+γϖ~-i~dϖ~+υ~qdϖ~dt~=a(i~q-ϖ~)-T~L.
Here the updated parameters are
(4)γ=-τnψr2BL,a=τBJ,υ~q=τnψrBRυq,υ~d=τnψrBRυd,T~L=τ2TLJ.
Let x=ϖ~, y=i~q, and z=i~d; one obtains that
(5)x˙=a(y-x)-T~Ly˙=γx-y-xz+υ~qz˙=-z+xy+υ~d.
If T~L=υ~q=υ~d=0, the system (5) evolves into the famous Lorenz system with the parameter ρ=1, whose dynamical behavior has been deeply explored [18, 19]. For the system (5) with υ~q≠0, there are two main cases: T~L=υ~d≠0 and T~L=υ~d=0. According to the numerical simulations and experimental investigations, it was found that the dynamics of the system (5) with υ~q≠0 and T~L=υ~d≠0 is similar to that of the system (5) with υ~q≠0 and T~L=υ~d=0. For the ease of analysis, this paper will only focus on the case with υ~q≠0 and T~L=υ~d=0, and a new type of double-deck butterfly chaotic attractor is obtained by adjusting the control signal υ~q, which has not been found in the existing literatures.
3. Bifurcation and Chaos Analysis of the PMSM Model
To ensure that the system (5) can be realized by electrical circuits, the variation of state variables x, y, and z of the system (5) must be limited into a proper range by the method of linear transformation and the linear transformation (x,y,z)→(10x,10y,100z) can be used without changing the properties of the original system. Then, the system (5) can be changed to
(6)x˙=a(y-x)y˙=γx-y-100xz+0.1υ~qz˙=-z+xy,
where υ~q=Asin(ωt). A and ω denote the amplitude and angular frequency of the driving signal υ~q, respectively. To analyze the system (6), there are many methods such as Lyapunov exponent spectrum [20], bifurcation diagram, phase portrait, and power spectrum, and they can be used to analyze the dynamics of nonlinear systems. The fourth-order Runge-Kutta method is used to solve the system (6) in MATLAB.
3.1. Phase Portraits and Poincare Sections
Considering the system (6) with υ~q=Asin(ωt), we set a=10 and γ=100, respectively. Varying parameters A and ω, the dynamics of the system (6) is to be explored.
When the amplitude A and the angular frequency ω of the driving signal υ~q are not equal to zero, the system (6) is a nonautonomous three-dimensional system. By adjusting parameters A and ω, a new chaotic attractor, which was not found in the existing dynamical systems, occurs when A=0.2 and ω=0.6, and meanwhile, uncommon multiperiodic orbits are observed when A=1 and ω=0.6.
Case 1 (A=0.2 and ω=0.6).
When A=0.2 and ω=0.6, the system is chaotic, but the new chaotic attractor is different from the famous attractor of the Lorenz-family systems [21] and is different from any other attractors contained in the known nonlinear dynamical system. Currently, four-wing chaotic attractors have been investigated by many researchers since the four-wing chaotic attractor shows complicated dynamics. From the shape of attractor shown in Figure 2, it can be seen that the new attractor is similar to four-wing chaotic attractor [22, 23], but the topologies are completely different between the new attractor and the known four-wing chaotic attractors since the new attractor contains two double-wing chaotic attractors, and we call this new chaotic attractor “double-deck butterfly chaotic attractor.” To check whether the system (6) is really chaotic when A=0.2 and ω=0.6, we calculated the Lyapunov exponents of the system (6) via the Jacobian method and obtained the Lyapunov exponents (0.042, 0, −0.341, −11.701). Owing to the system (6) being a three-dimensional nonautonomous system, the calculating results show four exponents, among which one is greater than zero, so the system (6) is chaotic.
From Figures 1(a)–1(d), it is very clear that the new chaotic attractor is constructed by double-deck butterfly chaotic attractor, which is different from the common four-wing chaotic attractors. The projections of phase portrait on x-y, x-z, and y-z plane are shown in Figures 1(a)–1(c) and the 3D shape is in Figure 1(d), respectively.
Phase portraits of the system (6) with A=0.2, ω=0.6, a=10, and γ=100: (a) x-y plane; (b) x-z plane; (c) y-z plane; (d) x-y-z space.
Poincare sections: (a) z=0; (b) y=0; (c) x=0.
x-y plane
x-z plane
y-z plane
There is another method to describe the dynamical behavior of the system (6), that is, Poincare map. Via Poincare sections, we can understand the motion of phase trajectories. Generally, the Poincare maps present several points that mean that the dynamics is periodic orbits, and countless points that form a closed curve mean that the dynamics is quasiperiodic orbits, and countless points that are self-similarly distributed at a specific region and cannot form one or more closed curves mean that the dynamics is chaotic.
On Case 1, the corresponding Poincare maps of the system (6) projected on different planes are exhibited in Figures 2(a)–2(c). Since the system (6) is chaotic, the distinct points on the Poincare sections cannot form one or more closed curves.
Case 2 (A=1 and ω=0.6).
When A=1 and ω=0.6, a multiperiodic orbit can be found, as shown in Figures 3(a)–3(c). It can be seen from Figure 3 that the type of multiperiodic orbit has the feature of helix which is rare and different from that of periodic orbits in the known dynamical systems.
Phase portraits of the system (6) with A=1, ω=0.6, a=10, and γ=100: (a) x-y plane; (b) x-z plane; (c) y-z plane.
3.2. Analysis of Lyapunov Spectra and Bifurcation Diagrams
The observations in Section 3.1 are several cases under different parameters. In order to deeply explore the dynamics of the system (6) in a fixed parameter range, it is necessary to investigate the system (6) using Lyapunov spectra and bifurcation diagrams. For a nonlinear dynamical system, the largest Lyapunov exponent can denote dynamical behavior. The largest Lyapunov exponent is expressed by λmax, and then λmax>0, =0 and <0 represent the dynamics of the system to be chaos, periodic orbit, and convergence, respectively.
For υ~q=Asin(ωt), we set that the amplitude A is in [0, 0.6] and the angle frequency ω is in [0.4, 1], respectively. Fixing A=0.2 and changing the angle frequency ω, the dynamical behavior of the system (6) can be shown in Figure 4. Similarly, fixing ω=0.6 and adjusting the amplitude A, the dynamical behavior of the system (6) is exhibited in Figure 5. Here, Figures 4(a) and 5(a) are the bifurcation diagrams of the state variable x for increasing ω and A, respectively. The corresponding largest exponent spectra are shown in Figures 4(b) and 5(b). From Figures 4 and 5, we can investigate the dynamics of the system (6) along with A and ω, respectively.
Dynamical behavior of the system (6) with ω in [0.4, 1] and A=0.2: (a) bifurcation diagram for increasing w; (b) the largest Lyapunov exponent versus ω.
Dynamical behavior of the system (6) with A in [0, 0.6] and ω=0.6: (a) bifurcation diagram for increasing A; (b) the largest Lyapunov exponent versus A.
Further, it is necessary to characterize the dynamical behavior of the system (6) if parameters A and ω simultaneously change. Therefore, we give a two-dimensional dynamics distribution shown in Figure 6 along with the change of A in [0, 0.6] and ω in [0.4, 1], where the dark-cyan region denotes chaotic motion and the magenta region represents periodic motion.
Two-dimensional dynamical behavior distribution on A-ω plane (colored figure online).
3.3. Power Spectra Analysis
The different dynamics subjected to different system parameters are described by the phase portraits and the bifurcation diagrams that can be used to intuitively understand dynamics of the system (6). However, the phase portraits and the bifurcation diagrams cannot provide more information concerning frequency characteristics of a nonlinear system, and the power spectrum analysis is a useful tool to explore frequency characteristics.
As we all know, the power spectrum of periodic motion gives peaks at a primary frequency and its harmonics, and the power spectrum of quasiperiodic motion gives peaks at linear combinations of two or more irrationally related frequencies, and the power spectrum of chaotic motion gives continuous broadband components to the spectrum. When A=1 and ω=0.6, the time series and power spectrum of variable x of the system (6) are shown in Figures 7(a) and 7(b), respectively. From Figure 7(a), it is obvious that the time series of variable x is periodic and the frequency f is approximately equal to 9/95≈0.095 Hz, which is in concordance with ω=0.6; namely, f=ω/2π≈0.095 Hz. Moreover, the primary frequency 0.095 Hz and its higher harmonics 0.287 Hz, 0.475 Hz, 0.67 Hz, 0.86 Hz, 1.05 Hz, 1.24 Hz, and 1.43 Hz are shown in Figure 7(b) which further validate that the dynamics of the system (6) is periodic.
Time series and power spectrum of the system (6) with A=1 and ω=0.6.
Time series of x
Local single-sided amplitude spectrum of x
Chaotic signal, in contrast, has the properties of nonperiodicity, so the power spectrum is a continuous broadband. For the system (6), if A=0.2 and ω=0.6, the time series and power spectrum of variable x are shown in Figures 8(a) and 8(b), respectively. Figure 8(a) shows random-like dynamics of variable x and the corresponding power spectrum shown in Figure 8(b) is continuous spectrum with noise-like background, which is one of the characteristics of chaotic signal. Through further analyzing of the power spectrum shown in Figure 8(b) and comparing it with the periodic power spectrum shown in Figure 7(b), it is easy to find that the primary frequency and its higher harmonics of the periodic signal generated from the system (6) at parameters A=1 and ω=0.6 are still existing and act as the salient frequency components, but some new frequency components between the primary frequency and its higher harmonics are excited because of the change of the parameters of the driving signal υ~q. The excited frequency components are filled in other spectra of frequency, which forms a continuous frequency band. Because of its inadequate information of chaotic time series, richer frequency components are not wholly exhibited in Figure 8(b).
Time series and power spectrum of the system (6) with A=0.2 and ω=0.6.
Time series of x
Local single-sided amplitude spectrum of x
4. Circuit Implementation and Experimental Results4.1. Experimental Circuit
A simpler implementation of a chaotic system can be achieved by an electronic circuit. With the help of the oscillators, the dynamics of many chaotic systems constructed by electronic circuit can be well demonstrated, which has been deeply researched in many chaotic systems. How to realize a chaotic experimental circuit has been investigated fully in [24, 25]. Here, the method is not introduced in details, but the schematic circuit and some necessary explanations are given.
The analog circuit for the nonautonomous system (6) is shown in Figure 9 and this circuit consists of three channels to conduct the integration of the three state variables x, y, and z, respectively. The operational amplifiers LF347 (A1–A6) and associated circuitry perform the basic operations of addition, subtraction, and integration, and the values of resistors and capacitors are labeled in Figure 9. The nonlinear terms of the system (6) are implemented with two analogue multipliers AD633 (U1 and U2).
Circuit implementation of the system (6).
The periodic excitation signal υ~q=Asin(ωt) is directly derived from a function generator YB1638. The input voltage frequency w of the periodic excitation signal is 385 Hz and is fixed. It is necessary to point out that the angular frequency ω=2πf, where f is the frequency, can be directly and freely changed in the experimental environment.
Therefore, it is very convenient for us to change the value of parameters A and ω by adjusting the corresponding knobs in control panel. When we adjusted the amplitude A of the function generator YB1638, rich bifurcation phenomena emerged in the analogue oscilloscope YB4360.
4.2. Experimental Results
The bifurcation processes of chaos and periodic orbits can be obtained by changing the magnitude A of YB1638, and the corresponding strange chaotic attractors and periodic orbits are shown in Figures 10(a)–10(c) and Figures 11(a)–11(c). As compared with Figures 1(a)–1(c) and Figures 3(a)–3(c), a good qualitative agreement between the numerical simulation and the experimental data was obtained. All results were recorded by YB4360.
Chaotic motion shown in the phase portraits of the system (6) with A=0.8 V and f=385 Hz: (a) x-y plane; (b) x-z plane; (c) y-z plane.
Multiperiodic orbit shown in the phase portraits of the system (6) with A=1.4 V and f=385 Hz: (a) x-y plane; (b) x-z plane; (c) y-z plane.
But it should be pointed out that the obtained experimental results given above are not exactly the same as the numerical solutions in MATLAB even if the parameters were the same, since the chaos pattern is very sensitive to tiny changes in initial conditions (or during the evolution) and the modeling errors and environmental noise are unavoidable in hardware experiments.
5. Conclusions
In this paper, we put the emphasis on the investigation of a PMSM model with smooth air gap. By setting υ~q=Asin(ωt), a new nonlinear system having complex dynamical behavior was obtained. The resulting nonautonomous system with given parameters, which can be adjusted shrewdly and accurately by controlling the frequency and amplitude of the external periodic excitation, has shown significant characteristics that there exist a strange double-deck butterfly attractor and multiperiodic orbits. These phenomena have also been validated by the numerical simulations and an electronic circuit, and the experimental circuit showed that the phenomenon really exists in nature. Moreover, because the obtained chaotic attractor has a more complex topological structure that is not found in existing systems and is different from the existing four-wing chaotic attractors, the system (6) may have a good application value in the field of information technology such as secure communication and encryption. In the future we will conduct research on how to control the PMSM to escape from the chaotic behavior to protect the motors in the applications.
Conflict of Interests
The authors declare that they have no conflict of interests.
Acknowledgments
This work was supported in part by the Natural Science Foundation of China under Grant nos. 61174094 and 11202148, the Tianjin Natural Science Foundation under Grant no. 14JCYBJC18700, the Shandong Natural Science Foundation under Grant no. ZR2012FM034, the China/South Africa Research Cooperation Programme under Grant nos. 78673 and CS06-L02, and the South African National Research Foundation under Incentive Grant no. 81705.
WeiD. Q.LuoX. S.WangB. H.FangJ. Q.Robust adaptive dynamic surface control of chaos in permanent magnet synchronous motor20073631-2717710.1016/j.physleta.2006.10.0742-s2.0-33846623353MademlisC.AgelidisV. G.On considering magnetic saturation with maximum torque to current control in interior permanent magnet synchronous motor drives200116324625210.1109/60.9372042-s2.0-0035451139BatzelT. D.LeeK. Y.Electric propulsion with sensorless permanent magnet synchronous motor: implementation and performance200520357558310.1109/TEC.2005.8529562-s2.0-24644508055JollyL.JabbarM. A.QinghuaL.Optimization of the constant power speed range of a saturated permanent-magnet synchronous motor20064241024103010.1109/TIA.2006.8760672-s2.0-33746343020LeeJ. H.JangY. J.HongJ. P.Characteristic analysis of permanent magnet-assisted synchronous reluctance motor for high power application2005971010Q50310.1063/1.18518312-s2.0-20944434851 ZhangX.WangX.DuJ.TangR.Amelioration of coreless permanent-magnet disk synchronous motor based on FEM: motor with wedge airgap200613142743210.1016/S1006-706X(08)60223-X2-s2.0-55249116506QiL.WangX.DuJ.ZHAOF.TANGR.Magnetic field analysis of disk coreless permanent magnet synchronous motor based on thickness variation Halbach array200613143944310.1016/S1006-706X(08)60225-32-s2.0-55249124032ChiS.ZhangZ.XuL.Sliding-mode sensorless control of direct-drive PM synchronous motors for washing machine applications200945258259010.1109/TIA.2009.20135452-s2.0-64349122701ParkN.JangM.LeeJ.HongK.KimJ.Performance improvement of a PMSM sensorless control algorithm using a stator resistance error compensator in the low speed region201010548549010.6113/JPE.2010.10.5.4852-s2.0-77956909589ChakrabortyS.KellerE.RayA.MayerJ.Detection and estimation of demagnetization faults in permanent magnet synchronous motors20139622523610.1016/j.epsr.2012.11.0052-s2.0-84871725765LiZ.ParkJ. B.JooY. H.ZhangB.ChenG.Bifurcations and chaos in a permanent-magnet synchronous motor200249338338710.1109/81.9891762-s2.0-0036506153ZribiM.OteafyA.SmaouiN.Controlling chaos in the permanent magnet synchronous motor20094131266127610.1016/j.chaos.2008.05.0192-s2.0-67249116860XueW.GuoY. L.ChenZ. Q.Analysis of chaos and circuit implementation of a permanent magnet synchronous motor2009581281468151MR26669122-s2.0-75349111859JingZ. J.YuC.ChenG. R.Complex dynamics in a permanent-magnet synchronous motor model200422483184810.1016/j.chaos.2004.02.0542-s2.0-2442603766ChoiH. H.Adaptive control of a chaotic permanent magnet synchronous motor20126931311132210.1007/s11071-012-0349-7MR29433872-s2.0-84863984998WeiD.ZhangB.LuoX.Adaptive synchronization of chaos in permanent magnet synchronous motors based on passivity theory201221303050410.1088/1674-1056/21/3/0305042-s2.0-84863351732AtaeiM.KiyoumarsiA.GhorbaniB.Control of chaos in permanent magnet synchronous motor by using optimal Lyapunov exponents placement2010374414226423010.1016/j.physleta.2010.08.047ZBL1238.341182-s2.0-77956393823MischaikowK.MrozekM.Chaos in the Lorenz equations: a computer assisted proof. Part II: details1998672231023104610.1090/S0025-5718-98-00945-4MR14593922-s2.0-0032393729MischaikowK.MrozekM.SzymczakA.Chaos in the Lorenz equations: a computer assisted proof Part III: classical parameter values20011691175610.1006/jdeq.2000.3894MR18084602-s2.0-0035216342BaskM.GençayR.Testing chaotic dynamics via Lyapunov exponents19981141-21210.1016/S0167-2789(97)00306-0MR16120392-s2.0-0002479384A. AlgabaF. Fernández-SánchezM. MerinoA. J. Rodríguez-LuisCenters on center manifolds in the Lorenz, Chen and Lü systems201419477277510.1016/j.cnsns.2013.08.003MR31192572-s2.0-84883090537CangS.QiG.ChenZ.A four-wing hyper-chaotic attractor and transient chaos generated from a new 4-D quadratic autonomous system201059351552710.1007/s11071-009-9558-0MR25929342-s2.0-77949274584WangZ. H.QiG.SunY.van WykB. J.van WykM. A.A new type of four-wing chaotic attractors in 3-D quadratic autonomous systems201060344345710.1007/s11071-009-9607-8MR26450352-s2.0-77953358076CangS.ChenZ.YuanZ.Analysis and circuit implementation of a new four-dimensional non-autonomous hyper-chaotic system2008573149315012-s2.0-41949087109LiuL.LiuC.ZhangY.Theoretical analysis and circuit implementation of a novel complicated hyperchaotic system201166470771510.1007/s11071-011-9943-3MR28595962-s2.0-82255182193