Abstract
This paper is a step forward to generalize the fundamentals of the relationship between fractional-order system modeling and integer order system modeling. First, we establish a fractional-order mathematical modeling of permanent magnet synchronous motor from the view of engineering applications. Obviously, the effect of fractional order is the key factor for extra freedom, more flexibility and novelty. As a pioneering work, we analyze the effect of fractional order on the bifurcation point, and give its basic law, which is a fundamental work to optimize the use of the previous integer-order research results, and a bridge between fractional-order system and integer order system. Furthermore, we discuss the beginning point and ending point of chaos and the size of the attractors with varying fractional-order, respectively. In addition, nonlinear dynamics analyses of the presented fractional order system are also given. Finally, numerical simulations are given which match the analytical discussion.
1. Introduction
Permanent magnet synchronous motor (PMSM) is widely used in high performance applications due to high efficiency, high power density, small size, and simple structure (Chang et al., 2011; Wei and Wang, 2013; Liu and Zhu, 2014; Ananthamoorthy and Baskaran, 2013). However, for those cases where low vibration, low acoustic noise, and high accurateness are highly demanded, magnetically induced vibration is a major concern (Shanmugasundram et al., 2013). Consequently, most of the prior contributions are major in new control methods with a good ability of robustness, and not try to solve its essential problem-basic mathematical model and vibration analysis (Choi and Jung, 2012; Yu et al., 2013; Arumugam et al., 2014). Because of this concern, it is essential to build new models approaching the engineering application and understand the excitation and vibration behaviors.
Fractional calculus dates back to more than four centuries, however, great improvements in the study of fractional calculus were made in the most recent twenty years (Li et al., 2013), because fractional calculus depends on the history of the function, which is more realistic and suitable for modeling, analyzing, and synthesizing and for solving many problems in basic natural sciences. Now, the fractional-order generalization has been applied to nonlinear dynamics, mathematics, physics, mechanics, control theory and engineering, agriculture, electromagnetic (Bhalekar et al., 2012; Chen et al., 2012, 2013, 2014; Jia et al., 2013; Golmankhaneh et al., 2013; Ho et al., 2013; Kayedi-Bardeh et al., 2012; Irandoust-pakchin et al., 2013; Grahovac and Spasic, 2013; Lopes and Machado, 2013; Romanovas et al., 2013; Jalloul et al., 2013; Wu and Baleanu, 2014; Liu et al., 2014; Zhang et al., 2014; Teng et al., 2014; Faieghi et al., 2014a, 2014b; Wu et al., 2014). The challenges in realizing the fractional element prevail through research and practice. Many realizations and patents of the fractional-order element using different technologies have been investigated. Therefore, we try to establish a fractional-order mathematical modeling of PMSM.
Motivated by the above discussions, this paper brings attention to establishing a fractional-order mathematical modeling of PMSM, and discussing nonlinear dynamic behavior. This first point is to build a novel fractional-order mathematical modeling of PMSM from the view of engineering applications. Second, we study the nonlinear dynamics of the above fractional-order modeling of PMSM in detail, especially the effect of fractional order on bifurcation point, which is a prior work. Furthermore, we discuss the beginning point and ending point of chaos and the size of the attractors with varying fractional-order, respectively.
This paper is organized as follows: section 2 build a fractional-order nonlinear mathematical modeling of PMSM, in section 3, nonlinear dynamics of the above PMSM are studied in detail and section 4 concludes this paper.
2. System modelling
Here, we suppose that PMSM have uniform air gap. Also, the other assumptions are:
The stator winding has wye connection, and balanced three-phase sine currents flow in the three-phase symmetric winding. Stator-flux-orientation has sine distribution, for that we neglect the effect of harmonic distortion and saturation flux. We also neglect the effect of magnetic hysteretic losses and eddy current loss. We set the damping of the rotor winding is zero. A diagram of the space vector of PMSM.
Usually, we use d-q rotor coordinate system, which will synchronous rotate with rotor. Therefore, we can easily get the space vector of PMSM as shown in Figure 1. Thus, we can get the voltage equation, ferromagnetic chain equation and torque equation as follows

The PMSM parameters.
From equation (1) to equation (6), we can get the equation of state as
We set
Using affine invariance and time scaling property, one gets
By simplifying equation (9), we get
3. Nonlinear dynamic analysis
In the following section we will study nonlinear dynamics of the above fractional-order PMSM in detail. Here, the numerical simulation results were carried out by using algorithm Adams-Bashforth-Moulton with a fixed step size of 0.01.
From Figures 2 and 3, the system goes directly through Hopf bifurcation to chaos. Moreover, the time domain waveforms and the diagram of phase trajectory agree with each other. When q1 = q2 = 1.00, there are period three between chaos from Figure 4. With the increasing of fractional order q1 = q2, a wide variety of bifurcation phenomena occurs. For example, from Figure 8, there is paroxysmal flashes of multi-periods near period two; from Figure 9, the system goes through period two to chaos omission of period one; from Figures 12 and 13, special bifurcation occurs.
Bifurcation diagrams, time domain waveform and phase portrait of fractional-order PMSM with fractional order q1 = q2 = 0.95, q3 = 1.00, parameter a = 50, initial condition x0 = 35, y0 = 0.02, z0 = 0.01. (a) Bifurcation diagram versus parameter b for fractional-order PMSM with fractional-orders q1 = q2 = 0.95, q3 = 1; (b1) time domain waveform with b = 9.0; (b2) phase portrait with b = 9.0, quasi periodic orbit; (c1) time domain waveform with b = 9.6; (c2) phase portrait with b = 9.6, double-scroll attractor. Bifurcation diagrams, time domain waveform and phase portrait of fractional-order PMSM with q1 = q2 = 0.98, q3 = 1.00, parameter a = 50, initial condition x0 = 35, y0 = 0.02, z0 = 0.01. (a) Bifurcation diagrams versus parameter b for fractional-order PMSM with fractional-orders q1 = q2 = 0.98 and q3 = 1.00; (b1) time domain waveform with b = 2.5; (b2) phase portrait with b = 2.5, quasi periodic orbit; (c1) time domain waveform with b = 3.0; (c2) phase portrait with b = 3.0, double-scroll attractor. Bifurcation diagrams, time domain waveform and phase portrait of fractional-order PMSM for fractional-orders q1 = q2 = 1.00, q3 = 1.00 with parameter a = 50, initial condition x0 = 35, y0 = 0.02, z0 = 0.01. (a) Bifurcation diagrams versus parameter b for fractional-order PMSM with fractional-orders q1 = q2 = 1.00, q3 = 1.00; (b1) time domain waveform with b = 1.5; (b2) phase portrait with b = 1.5, quasi periodic orbit; (c1) time domain waveform with b = 2.0; (c2) phase portrait with b = 2.0, double-scroll attractor. Bifurcation diagrams, time domain waveform and phase portrait of fractional-order PMSM for fractional-orders q1 = q2 = 1.02, q3 = 1.00 with parameter a = 50, initial condition x0 = 35, y0 = 0.02, z0 = 0.01. (a) bifurcation diagram versus parameter b for fractional-order PMSM with fractional-orders q1 = q2 = 1.02, q3 = 1; (b1) time domain waveform with b = 2.0; (b2) phase portrait with b = 2.0, period-2 orbit; (c1) time domain waveform with b = 2.6; (c2) phase portrait with b = 2.6, double-scroll attractor; (d1) time domain waveform with b = 2.8; (d2) phase portrait with b = 2.8, double-scroll attractor; (e1) time domain waveform with b = 2.915; (e2) phase portrait with b = 2.915, period-4 orbit; (f1) time domain waveform with b = 4.17; (f2) phase portrait with b = 4.17, period-3 orbit. Bifurcation diagrams, time domain waveform and phase portrait of fractional-order PMSM for fractional-orders q1 = q2 = 1.05, q3 = 1.00 with parameter a = 50, initial condition x0 = 35, y0 = 0.02, z0 = 0.01. (a) bifurcation diagram for fractional-order PMSM with fractional-orders q1 = q2 = 1.05, q3 = 1; (b1) time domain waveform with b = 1.7; (b2) phase portrait with b = 1.7, period-1 orbit; (c1) time domain waveform with b = 2.5; (c2) phase portrait with b = 2.5, double-scroll attractor; (d1) time domain waveform with b = 2.8; (d2) phase portrait with b = 2.8, double-scroll attractor. Bifurcation diagrams, time domain waveform and phase portrait of fractional-order PMSM for fractional-orders q1 = q2 = 1.10, q3 = 1.00 with parameter a = 50, initial condition x0 = 35, y0 = 0.02, z0 = 0.01. (a) bifurcation diagram for fractional-order PMSM with fractional-orders q1 = q2 = 1.10, q3 = 1; (b1) time domain waveform with b = 2.0; (b2) phase portrait with b = 2.0, period-1 orbit; (c1) time domain waveform with b = 2.7; (c2) phase portrait with b = 2.7, double-scroll attractor; (d1) time domain waveform with b = 2.9; (d2) phase portrait with b = 2.9, double-scroll attractor. Bifurcation diagrams, time domain waveform and phase portrait of fractional-order PMSM for fractional-orders q1 = q2 = 1.14, q3 = 1.00 with parameter a = 50, initial condition x0 = 35, y0 = 0.02, z0 = 0.01. (a) bifurcation diagram for fractional-order PMSM with fractional-orders q1 = q2 = 1.14, q3 = 1; (b1) time domain waveform with b = 2.2; (b2) phase portrait with b = 2.2, period-1 orbit; (c1) time domain waveform with b = 2.556; (c2) phase portrait with b = 2.556, period-3 orbit; (d1) time domain waveform with b = 2.57; (d2) phase portrait with b = 2.57, period-3 orbit; (e1) time domain waveform with b = 2.8; (e2) phase portrait with b = 2.8, double-scroll attractor; (f1) time domain waveform with b = 3.0; (f2) phase portrait with b = 3.0, double-scroll attractor. Bifurcation diagrams, time domain waveform and phase portrait of fractional-order PMSM for fractional-orders q1 = q2 = 1.15, q3 = 1.00 with parameter a = 50, initial condition x0 = 35, y0 = 0.02, z0 = 0.01. (a) bifurcation diagram for fractional-order PMSM with fractional-orders q1 = q2 = 1.15, q3 = 1; (b1) time domain waveform with b = 2.606; (b2) phase portrait with b = 2.606, period-3 orbit; (c1) time domain waveform with b = 2.87; (c2) phase portrait with b = 2.87, double-scroll attractor; (d1) time domain waveform with b = 2.946; (d2) phase portrait with b = 2.946, double-scroll attractor. Bifurcation diagrams, time domain waveform and phase portrait of fractional-order PMSM for fractional-orders q1 = q2 = 1.20, q3 = 1.00 with parameter a = 50, initial condition x0 = 35, y0 = 0.02, z0 = 0.01. (a) bifurcation diagram for fractional-order PMSM with fractional-orders q1 = q2 = 1.20, q3 = 1; (b1) phase portrait with b = 7.0; (b2) attractor projection on y-z plane with b = 7.0; (b3) attractor projection on x-z plane with b = 7.0; (b4) attractor projection on x-y plane with b = 7.0. Bifurcation diagrams, time domain waveform and phase portrait of fractional-order PMSM for fractional-orders q1 = q2 = 1.25, q3 = 1.00 with parameter a = 50, initial condition x0 = 35, y0 = 0.02, z0 = 0.01. (a) bifurcation diagram for fractional-order PMSM with fractional-orders q1 = q2 = 1.25, q3 = 1; (b1) phase portrait with b = 12.5; (b2) attractor projection on y-z plane with b = 12.5; (b3) attractor projection on x-z plane with b = 12.5; (b4) attractor projection on x-y plane with b = 12.5. Bifurcation diagrams, time domain waveform and phase portrait of fractional-order PMSM for fractional-orders q1 = q2 = 1.30, q3 = 1.00 with parameter a = 50, initial condition x0 = 35, y0 = 0.02, z0 = 0.01. (a) bifurcation diagram for fractional-order PMSM with fractional-orders q1 = q2 = 1.30, q3 = 1; (b1) phase portrait with b = 26.4; (b2) attractor projection on y-z plane with b = 26.4; (b3) attractor projection on x-z plane with b = 26.4; (b4) attractor projection on x-y plane with b = 26.4. Bifurcation diagrams, time domain waveform and phase portrait of fractional-order PMSM for fractional-orders q1 = q2 = 1.35, q3 = 1.00 with parameter a = 50, initial condition x0 = 35, y0 = 0.02, z0 = 0.01. (a) bifurcation diagram for fractional-order PMSM with fractional-orders q1 = q2 = 1.35, q3 = 1; (b1) phase portrait with b = 48.07; (b2) attractor projection with b = 48.07; (b3) phase orbits on x-z plane with b = 48.07; (b4) phase orbits on x-y plane with b = 48.07; (c1) time domain waveform with b = 48.2; (c2) phase portrait with b = 48.2; (d1) time domain waveform with b = 50.1; (d2) phase portrait with b = 50.1.











On the other hand, for the shape of the attractors, it is similar to Chen's attractor at first with q1 = q2 = 1.14. However, from q1 = q2 = 1.14 to q1 = q2 = 1.35, the attractors are not similar to Lorenz’s attractors or Chen's attractors.
Furthermore, Bifurcation diagrams of x versus fractional orders q1 = q2 are also presented, from Figure 14. Interestingly, with the decreasing of fractional order q1 = q2, the system routs period-doubling bifurcation to chaos.
Bifurcation diagrams of x versus fractional orders q1 = q2 for fractional-order PMSM with parameter a = 50, initial condition x0 = 35, y0 = 0.02, z0 = 0.01; (a) b = 2.5; (b) b = 2.7.
From Figures 8 and 9, we also can observe a typical transient phenomena, which are that the maximum amplitude abruptly become larger near b = 3.0.
In the following contents, we will study nonlinear dynamical analysis in detail.
From Figure 2, the system with q1 = q2 = 0.95, q3 = 1 goes into chaotic state at b = 9.2, and there is not any period window. Interestingly, the amplitude of x abruptly increases near b = 9.5, which can be observed obviously comparing Figure 2(b1) and (c 1 ). From the viewpoint of an electrical engineer, the current appears ripple current at b = 9.0 while it becomes chaotic vibration at b = 9.6 as shown in Figure 2(b1) and (c1), respectively.
From Figure 3, the system with q1 = q2 = 0.98, q3 = 1.00 become chaos from b = 2.75. For b = 2.5, the system goes through quasi-periodicity to a stable point, as shown in Figure 3(a). For b = 3.0, the chaotic attractor of the system is similar to the Lorentz’s attractor. From the viewpoint of an electrical engineer, the current stays a current value at b = 2.5 while it becomes chaotic vibration with high value at b = 3.0 as shown in Figures 3(b1) and (c1), respectively.
From Figure 4, the system with q1 = q2 = 1.00, q3 = 1.00 abruptly goes into chaos without any transient process at b = 1.85, although there are rich period windows from b = 1.85 to b = 1.85. Moreover, the attractor similar to Lorentz system as shown in Figure 4(c2) is bigger than that is shown in Figure 3(c2).
From Figure 5, there is a wide period window, which can be also called periodic zone, from b = 1.29 to b = 2.72. In this regions, the system goes through inversed period doubling bifurcation to period two from chaos. Furthermore, period four and period three occur at b = 2.915 and b = 4.17 as shown in Figure 5(e) and (f), respectively. In other words, the current appears periodic vibration.
The bifurcation changes considerably, although the fractional order change little, from Figure 6. The system with q1 = q2 = 1.05, q3 = 1 goes through inverse period doubling bifurcation to period one from chaos at b = 0.77. Interestingly, the chaotic attractor of the system is not the same as before at b = 2.5 as shown in Figure 6(c). From b = 3.2 to b = 3.4, there is period two window. And then, the system goes through period doubling bifurcation to chaos. From an engineer's viewpoint, the current alternate periodic and chaotic vibration.
From Figure 7, the initial chaotic region has now narrowed. As the increasing of fractional order, there is not period windows. Comparing Figure 7(c) and (d), the magnitude of x obviously changes. From Figure 7(d), the divergence of the attractor occurs, which shows that the shape of the attractor is changing.
From Figure 8, the prior chaotic region disappears for the system with q1 = q2 = 1.14, q3 = 1 the system goes through period doubling bifurcation to chaos. Interestingly, there are flashes of the branch of multi-cycle on the branch of period two, which is extremely rare. For b = 2.556 and b = 2.57, the flashes of period three occurs as shown in Figure 8(c) and (d). Comparing Figure 8(e) and (f), the magnitude of x abruptly becomes bigger, which is a typical transient dynamical phenomena. From an engineer's viewpoint, the current suddenly becomes bigger. Now, the shape of the attractor is similar to Chen’s system as shown in Figure 8(f2).
From Figure 9, the period one disappears for the system with q1 = q2 = 1.15, q3 = 1. So, the system goes though period doubling bifurcation to chaos from period two. Interestingly, the flashes of multi-cycle occurs near period two, which means that the current alternates between different states.
From Figure 10, the bifurcation changes significantly for the system with q1 = q2 = 1.20, q3 = 1, which means that the system goes into chaos directly. Moreover, the shape of the chaotic attractor becomes surprising, which is not similar to typical attractors. From Figure 10(b2), the attractor on y-z plane is centrosymmetric. While, the attractor is axisymmetric from Figure 10(b3) and (b4).
From Figure 11, the system with q1 = q2 = 1.25, q3 = 1 goes into chaos directly, and there is obviously period window. Moreover, the shape of the attractor also changes as shown in Figure 11(b2). From Figure 11(b2), the attractor on y-z plane is centrosymmetric. While, the attractor is axisymmetric from Figure 11(b3) and (b4).
Form Figure 12, the bifurcation phenomena changes momentously for the system with q1 = q2 = 1.30, q3 = 1. More specially, the system goes through turbulence to period one from a quick chaos. Moreover, the shape of the attractor also changes as shown in Figure 12(b). From Figure 12(b2), the attractor on y-z plane is centrosymmetric. While, the attractor is axisymmetric from Figure 12(b3) and (b4).
From Figure 13, the presented system with q1 = q2 = 1.35, q3 = 1 goes from chaos through quasi-periodicity to period one in a strange way. From Figure 13(c1), the current abruptly tends to infinity from period vibration at b = 48.2. It, however, holds periodic vibration at b = 50.1.
For both of b = 2.5 and b = 2.7, the general trend of the bifurcation are the same, which is that the system abruptly becomes chaos, and then goes through inverse period doubling bifurcation to period two. The same process will be occurs again during the next region. The difference is that the flashes of the multi-cycle bifurcation occurs for b = 2.5. Moreover, the system ends with period two for b = 2.5, while it ends with period four for b = 2.7.
We will discuss how the changes of fractional order affect the bifurcation. From Figures 2(a) to 13(a), obvious shifts can be observed, which is that the bifurcation point moves to the left as the fractional-orders of the system increase. In other words, the value of bifurcation parameter b decreases where bifurcation occurs with the increasing of fractional orders. More specially, when the system order satisfies q1 = q2 = 0.95, q3 = 1, the bifurcation point is b = 9.30 (Figure 2(a)), while q1 = q2 = 0.98, q3 = 1, the bifurcation point dramatically shifts to b = 2.75 (Figure 3(a)). When q1 = q2 = q3 = 1, bifurcation happens at b = 1.85 (Figure 4(a)); When q1 = q2 = 1.02, q3 = 1, bifurcation happens at b = 1.3 (Figure 5(a)); When q1 = q2 = 1.05, q3 = 1, bifurcation happens at b = 0.77 (Figure 6(a)); When q1 = q2 = 1.10, q3 = 1, bifurcation happens at b = 0.04, bur soon get rid of chaos and enter period 1 orbit (Figure 7(a)). These sorts of shifts can be observed for a relatively long time before the Hopf bifurcation disappears.
This interesting phenomena can be explained by the stability theory of fractional-order nonlinear system. As we all know that if and only if
Furthermore, we discuss the beginning point and ending point of chaos with varying fractional-order. From Figure 15(a), the beginning point and ending point of main chaotic region both decrease first and then increase with the increasing of fractional-order q1 = q2, whose minimum is at about q1 = q2 = 1. From Figure 15(b), the beginning point and ending point of the secondary chaotic region decrease until the chaotic region disappear increasing the fractional-order q1 = q2.
The values of parameter b for the beginning and end of chaos versus different fractional-order q1 = q2 with parameter a = 50, initial conditions x0 = 35, y0 = 0.02, z0 = 0.01. (a) regular figure; (b) detail with enlarged scale.
From Figure 16, we also find another interesting phenomena, which is that the size of the attractor becomes bigger with the increasing of fractional-order q1 = q2.
The attractors of the above PMSM system for different fractional-orders q1 = q2, with parameter a = 50, initial conditions x0 = 35, y0 = 0.02, z0 = 0.01. (a) q1 = q2 = 0.95; (b) q1 = q2 = 1.00; (c) q1 = q2 = 1.05; (d) q1 = q2 = 1.10; (e) q1 = q2 = 1.15; (f) q1 = q2 = 1.20.
4. Conclusion and discussions
This paper demonstrates the generalized concepts from the very narrow integer order scope system to fractional-order mathematical modeling. More specially, we established a fractional-order mathematical modeling of PMSM from the view of engineering applications. Moreover, we studied the nonlinear dynamics of the above fractional-order nonlinear system including bifurcation diagrams, waveforms in time domain and phase orbits. As a pioneering work, we analyzed the effect of fractional order on the bifurcation point, and gave its basic law, which is a fundamental work to optimize the use of the previous integer-order research results, and a bridge between fractional-order system and integer order system. Moreover, we discuss the beginning point and ending point of chaos and the size of the attractors with varying fractional-order, respectively. Finally, numerical simulations are given which match the analytical discussion.
In the future, we will try to propose new nonlinear dynamical indicators to analyze the effect of the fractional-order on the stability of nonlinear systems, and describe dynamic characters from a different viewpoint. Moreover, we will test the presented fractional-order mathematical modeling by conducting scientific experiments.
Footnotes
Funding
This work was supported by the Scientific Research Foundation of the National Natural Science Foundation (grant numbers 51109180 and 51279167), the National Science and Technology Supporting Plan from the Ministry of Science and Technology of China (grant numbers 2011BAD29B08 and 2012BAD10B02), the Fundamental Research Funds for the Central Universities (grant number 201304030577) and Scientific Research Funds of Northwest A&F University (grant number 2013BSJJ095).
