Research Article  Open Access
Threshold for Chaos of a Duffing Oscillator with FractionalOrder Derivative
Abstract
In this paper, the necessary condition for the chaotic motion of a Duffing oscillator with the fractionalorder derivative under harmonic excitation is investigated. The necessary condition for the chaos in the sense of Smale horseshoes is established based on the Melnikov method, and then the chaotic threshold curve is obtained. The largest Lyapunov exponents are provided, and some other typical numerical simulation results, including the time histories, frequency spectrograms, phase portraits, and Poincare maps, are presented and compared. From the analysis of the numerical simulation results, it could be found that, near the chaotic threshold curve, the system generates chaos via the perioddoubling bifurcation, from single periodic motion to period2 motion and period4 motion to chaotic motion. The effects of fractionalorder parameters, the stiffness coefficient, and the damping coefficient on the threshold value of the chaotic motion are analytically discussed. The results show that the coefficient of the fractionalorder derivative has great effect on the threshold value of the chaotic motion, while the order of the fractionalorder derivative has less. The analysis results reveal some new phenomena, and it could be useful for designing or controlling dynamic systems with the fractionalorder derivative.
1. Introduction
The development of fractionalorder calculus has a history of more than 300 years, but until recently, the research in this field still focused on pure mathematical theory. In recent years, the research and application of fractional calculus in the mechanical engineering field has gradually increased [1], especially to accurately characterize the true constitutive relationship of viscoelastic materials [2, 3] and controlling engineering [4, 5].
In many nonlinear dynamical systems, such as van der Pol, Mathieu, Hill’s oscillators, Duffing oscillators, and other systems, large number of results can be found in the literature [6–10]. The complex dynamic phenomena of typical nonlinear systems, such as periodic solutions [11, 12], bifurcation, and chaos [13], have been more and more concerned. Among these nonlinear dynamical phenomena, chaos is a unique dynamic behavior that involves nonperiodic motion originating from deterministic systems and depending on initial conditions. Its effect on mechanical systems cannot be ignored. Therefore, the study of chaotic threshold is of great significance to the understanding of the mechanism of chaos generation, which can lead to the control and anticontrol of chaos [14, 15]. A chaotic threshold is the boundary between the periodic motion and chaotic motion of the nonlinear dynamical system. A Duffing oscillator with the fractionalorder derivative is a typical nonlinear dynamic system with complex dynamic behaviors. At present, the research on fractionalorder Duffing oscillators mainly focuses on the analysis of dynamic response characteristics. Shen et al. [16] and Mobin and Saeed [17] studied Duffing oscillators with fractional order and obtained the amplitudefrequency response curve of the systems by using the average method. They also analyzed the effect of fractionalorder parameters on the primary resonance amplitude value and frequency. Wen et al. [18, 19] derived the approximate analytical solution of the fractionalorder Mathieu–Duffing equation based on the incremental harmonic balance method. Xu et al. [20] considered a complex Duffing system subjected to nonstationary random excitation. Chen et al. [21–23] studied the random jumping and bifurcation behaviors of fractionalorder Duffing oscillators under combined harmonic and white noise excitation, and they designed a controller based on the fractionalorder PID.
There has been some research on chaotic motions and synchronization control of chaos to the fractionalorder Duffing system. Liu et al. [24, 25] proposed a numerical simulation method called extended generalized unit mapping for a nonautonomous fractionalorder Duffing system, and they studied the chaotic turbulence inside the system. Zufiria and Jiménez [26] described the chaotic behavior of fractionalorder Duffing oscillators through the Lyapunov exponents and strange attractors in the effective phase space of the systems. Mohammad et al. [27] suggested a technique for finding unstable periodic orbits in chaotic fractionalorder systems, using a fractionalorder Duffing system as an example. Some scholars [28–30] analyzed the chaotic synchronization behavior of a modified van der PolDuffing system and revealed the effect of fractional order on chaotic synchronization controls. These examples show that applications of the fractionalorder derivative in engineering fields are becoming increasingly important.
The Melnikov method has certain advantages when it is used to analyze the chaotic motion threshold value of multiparameter nonlinear systems. For getting threshold values of the chaotic motion, many scholars have used the analytic Melnikov method to study the necessary condition for chaos of Duffing oscillators without a fractionalorder derivative [31–35]. However, due to the complexity of the fractionalorder derivative, it is rare to discuss its effect on the chaotic threshold of nonlinear fractionalorder systems. Anague Tabejieu et al. [36] studied the Rayleigh beams with fractional order, but there was no relevant literature on the threshold value of the chaotic motion of Duffing oscillators with a fractionalorder derivative.
In this paper, the necessary condition for generating chaos in sense of Smale horseshoes in a typical Duffing oscillator with the fractionalorder derivative under harmonic excitation is investigated based on the Melnikov method. The paper is organized as follows. In Section 2, the Melnikov function is established based on Melnikov theory, and the necessary conditions for chaos are obtained. The fractionalorder derivative is calculated by the numerical method, and the analytical expressions are obtained in the rest part. In Section 3, the necessary condition is verified by the maximum Lyapunov exponents, and the typical dynamical responses, including the time histories, frequency spectrograms, phase portraits, and Poincare maps of the system, are also investigated. The results show that this system has chaotic motions. The result obtained based on the numerical solution satisfies the chaos threshold curve, which fully verifies the correctness of the necessary conditions. In Section 4, the effects of fractionalorder parameters, stiffness coefficient, and the damping coefficient of the system are analyzed. Finally, the detailed conclusions are summarized and analyzed.
2. Necessary Condition for Chaos of FractionalOrder Duffing Oscillators
Duffing oscillators are one of the most typical and important objects in nonlinear dynamics. In this paper, the ordinary Duffing oscillator with a fractionalorder derivative under harmonic excitation is described withwhere m, k_{1}, k_{3}, c_{1}, F, and ω are the system mass, linear stiffness coefficient, nonlinear stiffness coefficient, linear viscous damping coefficient, excitation amplitude, and excitation frequency, respectively. The expression is the p^{th}order derivative from x(t) to t for the fractional coefficient h(h > 0), and the fractional order is restricted to 0 < p < 1. All the system parameters are positive and dimensionless.
Using the transformationswhere ε is a small real parameter, , equation (1) becomes
2.1. Homoclinic Orbits of FractionalOrder Duffing Oscillators
Equation (1) can be transformed into the form of a state equation:
If ε = 0, meaning that neither system damping nor external excitation exists, then equation (3) can be transformed into
Now the system is a Hamilton system whose Hamilton function satisfies f_{1} = ∂H/∂y and f_{2} = −∂H/∂x.
In our case,
Letting and leads to three fixed points: , , and .
According to equation (5), the characteristic equation of the system is
Using equation (7), one can obtain
When the fixed point is (0, 0), the eigenvalue of the system is . When the fixed point is or , the eigenvalue of the system is . According to the properties of a fixed point [37], fixed point (0, 0) is the saddle point of the system and fixed points and are the center of the system. Using equation (6), one can obtain the homoclinic orbit, and the result is shown in Figure 1. When the Hamilton function H(x, y) = 0, there are two homoclinic orbits connecting saddle points. The positions of a saddle point and two centers are indicated in Figure 1. The saddle point is unstable, and the centers are stable.
The homoclinic orbits of the system described in equation (1) can be obtained by solving the following expression:
Equation (8) leads to
The variables in equation (9) can be separated and the integral done, and this leads to the parametric equations of the two homoclinic orbits:
An analysis of equation (10) leads to the conclusion that a plus sign of x(t) represents the positive axis of a homoclinic orbit and a minus sign represents the negative part. The plus sign of y(t) represents the upper half of a homoclinic orbit, and minus sign represents the lower half.
2.2. Melnikov Function Determines the Chaotic System’s Necessary Condition
If , where and , equation (3) can be transformed into
Melnikov theory [38], proposed by Melnikov, is a perturbation method originally for proving the existence of transverse homoclinic or heteroclinic orbits. Melnikov theory has been applied in many researches. The Melnikov theory was firstly applied to study a periodically driven Duffing oscillator by Holmes [39]. According to the Melnikov process, one can determine the Melnikov function with y(t) as the odd function:
Substituting , , and into M(t_{0}) yields
According to Melnikov theory [37], this leads to
Since , equation (15) can be used to the get necessary condition:
The integral in equation (11) can be divided into three parts: I_{1}, I_{2}, and I_{3}. These three parts were calculated in the following.
2.2.1. Calculating I_{1} and I_{3}
We substituted the parametric equation (11) of the two homoclinic orbits into equation (14). Since I_{1} and I_{3} do not contain fractionalorder derivatives, the residue theorem and contour integration can be used to calculate I_{1} and I_{3} [38, 40]:
2.2.2. Calculating I_{2}
Since I_{2} in equation (14) contains a fractionalorder derivative, it is impossible to integrate directly, so it can only be calculated using the definition of a fractionalorder derivative. The Caputo definition and the Riemann–Liouville definition are as follows:where Γ(x) is the gamma function satisfying Γ(x + 1) = xΓ(x).
If the initial value of x(t) is nonzero, it is relatively easy to use the Caputo definition [41] to calculate the fractionalorder derivative. Thus, we used Caputo’s definition to calculate I_{2}. Comparing the two definitions of the fractionalorder derivative, the relationship between them can be obtained as follows:where m = ΓpΊ and f(x) is a penalty function.
The specific process of calculating the highaccuracy Caputo fractionalorder derivative is below. First, the improved direct recursive algorithm [42] is used to calculate the Riemann–Liouville fractionalorder derivative with high accuracy. Then, the penalty function f(t) is calculated, and finally, the penalty function is combined with the Riemann–Liouville fractionalorder derivative, leading to the highaccuracy Caputo fractionalorder derivative. In this way, the integral with the fractionalorder derivative can be calculated, and the value of I_{2} can be obtained:where
2.2.3. Discussion of the Parameter t_{0}
Equation (20) contains the parameter t_{0}, which is the initial moment when the stable manifold intersects the unstable manifold, that is to say, when the system generates chaos in the sense of the Smale horseshoe. The variable t_{0} can be regarded as a variable about time, , where T is the period of the system motion. We analyzed the effect of t_{0} on the chaos threshold value (the system parameters were m = 1, k_{1} = –1, c_{1} = 0.5, k_{3} = 1, and h = 0.5). Figure 2 shows the results of this analysis.
Figure 2 shows that an increase of t_{0} would cause the threshold extremum of the chaotic motion to decrease, reducing the chaotic area and making it more difficult to generate the chaos, but the parameter t_{0} has little effect on the threshold extremum of chaotic motion (the chaos thresholds of point A and point B are 2.807 and 2.775, respectively). Therefore, we selected t_{0} = 1 when analyzing the chaotic threshold of system.
2.2.4. Necessary Condition and Threshold Curves of a Chaotic System
Substituting the result of equations (17), (18), and (21) into equation (16), the necessary condition for generating chaos in the sense of Smale horseshoes can be obtained as follows:
Based on the relationship between the excitation amplitude F and the excitation frequency ω obtained from equation (23), the necessary condition of the system for generating chaos can be determined. The chaotic threshold curve is shown in Figure 3 (the system parameters are m = 1, k_{1} = 1, c_{1} = 0.5, k_{3} = 1, h = 0.5, and ).
Figure 3 shows that when the value of the parameter 1/F is located in the region below the curve (gray region), the system might generate chaos or not. Equation (23) shows a condition for chaos to occur, but it is a sufficient but not a necessary condition. When the value of the parameter 1/F is located in the region above the curve (white region), the system will have a stable periodic motion, which may be a singlecycle motion or a multicycle motion.
From the analysis of Figure 3, it could be concluded that whether chaotic motion occurs is affected by the force and frequency of external excitation. The following is an analysis of these two aspects:(1)When the external excitation force of the system is small, the system will not generate chaotic motions, and the system will have regular periodic motions. When the external excitation force increases and exceeds the critical value for chaos, the chaotic motion may occur in the system, and the larger the excitation force, the higher the possibility of the chaotic motion.(2)The chaotic threshold value first increases with the increasing external excitation frequency, then decreases, and finally approaches zero. At the resonance frequency, the system chaotic threshold value reaches a maximum. In the resonance region of the system, the vibration amplitude is larger, and chaos is more likely to occur.
Therefore, in the practical engineering of the systems, we should try to avoid the resonance region of the system in order to avoid the irregular chaos phenomenon.
3. Numerical Simulation Analysis of a FractionalOrder Duffing System
In order to verify the validity of the necessary condition for generating chaos, a set of illustrative system parameters is chosen as m = 1, k_{1} = 1, k_{3} = 1, c_{1} = 0.5, h = 0.5, and . The chaotic threshold curve of the external force F with respect to the excitation frequency ω is obtained. Above the curve is the region of the periodic motion, and below the curve is the region where chaos may occur. Different points in the two regions are selected, and the location of each point is shown in Figure 4. When F = 1, F = 1.5, and F = 1.75, the largest Lyapunov exponents are shown in Figures 5–7, respectively. The numerical simulation method for the largest Lyapunov exponents about the differential equation can be found in [43]. From the analysis of Figures 4–7, it can be found that the maximum Lyapunov exponents of points in the chaotic region is greater than zero, while the maximum Lyapunov exponents of points in the periodic motion region is less than zero. The results imply the chaotic threshold curve is correct.
Some typical points are shown in Figure 4 with solid dots. The points from A to E are in the region of the periodic motion, and points F, G, and H are in the region of the chaotic motion. Among them, points A and G would occur when the excitation amplitude F = 1, points H, E, D, and B would occur when the excitation amplitude F = 1.5, while points F and C would occur when F = 1.75. The horizontal abscissa of the intersection points of the three straight lines and chaotic edge curves are 1.45, 1.86, and 2, respectively, and the Lyapunov exponents are equal to zero at these frequencies. In order to further explain Figure 4, the largest Lyapunov exponents and the system responses at each point are analyzed, respectively, in the following text.
3.1. The Largest Lyapunov Exponents
The Lyapunov exponents (LEs) measure the exponential growth, or decay, of infinitesimal phase space perturbations of a chaotic dynamical system. We use the Benettin–Wolf algorithm to determine the Lyapunov exponent. The main steps to determine numerically the Lyapunov exponent are as follows: using the numerical method of power series expansion (PSE) to numerically integrate the fractionalorder differential equations to obtain the time series of x(t) [41], the Gram–Schmidt procedure and picking up the exponents during the renormalization procedure, and the Lyapunov exponents being determined as the average of the logarithm of the stretching factor of each perturbation. The numerical simulation method for the largest Lyapunov exponents about the Duffing oscillator with the fractionalorder derivative can be found in [44]. The largest Lyapunov exponents at F = 1 N, F = 1.5 N, and F = 1.75 N are shown in Figures 5–7.
Figure 5 shows that the maximum Lyapunov exponents of point G is greater than zero, the maximum Lyapunov exponents of point A is less than zero, and the critical excitation frequency for chaos is 1.51. Accordingly, the system performs periodic motions at point A and chaotic motions at point G.
Figure 6 shows that the maximum Lyapunov exponents of point H is greater than zero, the maximum Lyapunov exponents of points B, D, and E are less than zero, and the critical excitation frequency for chaos is 1.92. Accordingly, the system performs periodic motion at points B, D, and E and chaotic motion at point H.
Figure 7 shows that the maximum Lyapunov exponents of point F is greater than zero, the maximum Lyapunov exponents of point C is less than zero, and the critical excitation frequency for chaos is 2.06. Accordingly, the system performs periodic motions at point C and chaotic motions at point F.
Compared with Figures 4–7, it could be found that the critical excitation frequency based on the Melnikov method is little less than the numerical simulation method for the largest Lyapunov exponents (the difference value of critical excitation frequencies between the two methods is 0.06, 0.06, and 0.06, respectively). Although there is difference between the two results quantitatively, this difference is reasonable because the Melnikov method is a firstorder approximate way to find chaotic motions [32]. This could also be useful to design or revise the similar dynamical system. From the observation of Figures 5–7, it could be found that increase of the critical excitation frequency would be against the chaos generation. In addition, it could be found that the maximum Lyapunov exponents in the chaotic region may appear as the negative value (when the system is in the state of periodic motion) from the region surrounded by the black circles 1, 2, and 3 in Figures 5–7, which show that the chaotic threshold curve determined based on Melnikov theory is the necessary condition for the chaotic motion of the system.
3.2. Analysis of the System Responses at Each Point
In order to further verify the conclusion, the typical dynamical responses, including the time history, frequency spectrograms, phase portraits, and Poincare maps, at different typical points are presented and analyzed in the following text. During the numerical simulation process [19, 41], the adopted numerical scheme iswhere t_{l} = lh is the time sample points, h is the sample step, and C_{j}^{p} is the fractional binomial coefficient with the iterative relationship as
Here, we select h = 0.01, the total computation time is 500 excitation periods, and 500 points were sampled in each period. The frontal transient response is omitted, and the posterior 100 periods are selected for analysis. In the following section, the points (in Figure 4) in the region of the periodic motion and chaotic motion are simulated and analyzed. Some typical samples in different ranges are taken to analyze the dynamical phenomena. Through the analysis of these dynamical characteristics at different typical points, the conclusion is that the curve obtained in Figure 4 is the boundary between chaos and periodic motion.
3.2.1. Region of Periodic Motions
The graphs in Figure 8 illustrate the periodic motion when ω = 2 rad/s and F = 1 N (point A).
(a)
(b)
(c)
(d)
The graphs in Figure 9 illustrate the periodic motion when ω = 3 rad/s and F = 1.5 N (point B).
(a)
(b)
(c)
(d)
The graphs in Figure 10 illustrate the periodic motion when ω = 2.5 rad/s and F = 1.75 N (point C).
(a)
(b)
(c)
(d)
Figures 8–10 show that, for different external excitation frequencies and excitation amplitudes, the phase portraits have a stable phase trajectory, and the Poincare maps have only one fixed point. The dominant frequency on the frequency spectrograms indicates that the system has forced vibration when the frequency is equal to the excitation frequency. Therefore, the system performs the steadystate periodic motion.
The graphs in Figure 11 illustrate the period2 motion when ω = 2.2 rad/s and F = 1.5 N (point D).
(a)
(b)
(c)
(d)
Figure 11 shows that the Poincare maps have two fixed points, and there are dominant frequencies, a half frequency division, and weak frequency multiplication in the frequency spectrogram, showing that there was period2 motion in this case.
The graphs in Figure 12 illustrate the period4 motion when ω = 1.95 rad/s and F = 1.5 N (point E).
(a)
(b)
(c)
(d)
Figure 12 shows that the Poincare maps have four fixed points, and there are dominant frequencies, quarter frequency division, half frequency division, and weak frequency multiplication in the frequency spectrogram. The frequency multiplication will reduce the smoothness of the timedomain curve of the system, showing that there is period4 motion in this case.
3.2.2. Region of Chaotic Motions
The graphs in Figure 13 illustrate the chaotic motion when ω = 1.85 rad/s and F = 1.75 N (point F).
(a)
(b)
(c)
(d)
The graphs in Figure 14 illustrate the chaotic motion when ω = 1 rad/s and F = 1 N (point G).
(a)
(b)
(c)
(d)
The graphs in Figure 15 illustrate the chaotic motion when ω = 1.5 rad/s and F = 1.5 N (point H).
(a)
(b)
(c)
(d)
In Figures 13 and 14, the time histories show that the motion of the system is irregular. In the frequency spectrograms, except for a dominant frequency that is equal to the excitation frequency, there are continuous frequency components. Phase trajectories show disorder. The Poincare maps show a collection of discrete points. These all show that the system has the chaotic motion in these cases. In Figure 15, the timehistory response curve of the system shows the paroxysmal motion state with a random alternation between the regular and irregular motion. This state is one of the paths through which the system moves from the steadystate periodic motion to chaotic motion. With the change of system parameters, the time of the irregular motion becomes longer and longer. Finally, the system completely enters the chaotic motion state.
The system responses for different typical points are shown in Tables 1–3 when F = 1 N, F = 1.5 N, and F = 1.75 N.



Based on the above analytical results, we concluded that when the parameters of the system are above the chaotic threshold curve, the system shows the stable periodic motion, and when they are under the chaotic threshold curve, the system may have chaotic motions. From the dynamic response diagram of the system at the points B, D, E, and H, it can be seen that as the points gradually approach the chaotic threshold curve, the system changes from single periodic motion to period2 motion and period4 motion to chaotic motion. When it enters the chaotic motion region, it has reciprocating aperiodic motion, so there are abundant nonlinear dynamic phenomena in the region near the chaotic threshold curve. Accordingly, the numerical simulation results verify the correctness of the chaotic threshold curve determined by the previous calculations, and the results further illustrate the validity of the necessary condition obtained by the Melnikov method.
4. Effect of System Parameters on the Chaotic Threshold Curve
The main purpose of this section is to describe the effects of four parameters k_{1}, c_{1}, h, and p of the fractionalorder Duffing system on the chaotic threshold curve. Based on the necessary condition determined by equation (21), the chaotic threshold curve can be obtained. By changing the value of each parameter, the changes for the chaotic threshold curves can be analyzed. The results are presented individually in Figures 16–19. The natural frequency of the system discussed below is .
Figure 16 shows that increasing linear stiffness coefficient k_{1} leads to smaller threshold values for generating chaos, so the generation of chaos becomes more difficult. In addition, when the external excitation frequency is greater than a certain value, the attenuation speed of the chaotic threshold curve decreases with the increasing system stiffness, and it eventually approaches to zero. When the external excitation frequency is relatively large, the chaotic motion of the system is not likely to occur.
Figure 18 shows that when the linear damping coefficient c_{1} increases, the peak value of the chaotic threshold curve decreases, the area where chaos occurs decreases, and the possibility of chaos in the system decreases.
Figures 16 and 18 show that, with increasing stiffness and damping of the system, the peak value of the chaotic threshold curve decreases, the area where chaos may occur decreases, and the possibility of chaos in the system decreases. Compared with the damping of the system, the stiffness of the system has a greater impact on the chaotic region of the system. Therefore, in practical work, the chaotic phenomena can be avoided by adjusting the stiffness of the system.
Figure 19 shows that increasing coefficient h of the fractionalorder derivative leads to decreasing the peak value of the chaotic threshold curve. The area where chaos occurs decreases, and the possibility of chaos in the system decreases.
Figure 17 shows that increasing the order p of the fractionalorder derivative leads to increasing threshold values for generating chaos, increasing areas where chaos occurs, and increasing likelihood of chaos generation.
An analysis of Figures 16 and 17 shows that the parameters of the fractionalorder derivative have effects on the threshold value of chaos, especially the coefficient of the fractionalorder derivative term. We also found that the system would be less chaotic with the smaller order and larger coefficient of the fractionalorder derivative.
The fractionalorder derivative was introduced to more accurately describe the nonlinear dynamic behaviors of physical objects with fractionalorder characteristics, such as viscoelastic materials, damping, dynamic characteristics of gear clearance, friction, and impact in mechanical systems [45–48]. The analysis shows that the parameters of the fractionalorder derivative have an important effect on the chaotic threshold of the system.
The above analysis shows that if the system is under stationary excitation amplitude with a small amplitude, chaotic motion will not occur, but under large excitation amplitude, such as impact and chaotic motion, may occur. Since chaotic motion may increase the uncertainty of the system, appropriate parameters should be selected to avoid the occurrence of chaotic phenomena.
5. Conclusions
In this paper, threshold for chaos of a Duffing oscillator with the fractionalorder derivative is investigated. Furthermore, the effects of various parameters on the threshold value of chaotic motions are analyzed.(1)The necessary condition for generating chaos in the sense of Smale horseshoes is obtained based on the Melnikov method. The relationship between each parameter is established, and the effect of each parameter on the chaotic threshold is studied, respectively.(2)The typical dynamical responses, including the largest Lyapunov exponents, time histories, frequency spectrograms, phase portraits, and Poincare maps, are obtained by a numerical simulation method to analyze the dynamic response of the system. It could be found that the system has chaotic motions under certain conditions, and the transition from the steady periodic motion to chaotic motion is a gradual process, which verified the correctness of the necessary condition indirectly.(3)The effects of the system parameters on the chaotic threshold curve are analyzed in detail. It concluded that the stiffness coefficient and damping coefficient of the system are not only effects on the chaotic threshold curve, and the fractionalorder parameters also have a significant effect on the chaotic threshold curve.
The corresponding results could provide a reference for analyzing dynamic systems with a fractionalorder derivative, which is important to design or control this kind of system.
Data Availability
The authors declare that the underlying data related to our submission are available.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors are grateful to the support by the National Natural Science Foundation of China (no. 11172183).
References
 J. T. Katsikadelis, “Generalized fractional derivatives and their applications to mechanical systems,” Archive of Applied Mechanics, vol. 85, no. 910, pp. 1307–1320, 2015. View at: Publisher Site  Google Scholar
 H. Sun, Y. Zhang, D. Baleanu, W. Chen, and Y. Chen, “A new collection of real world applications of fractional calculus in science and engineering,” Communications in Nonlinear Science and Numerical Simulation, vol. 64, pp. 213–231, 2018. View at: Publisher Site  Google Scholar
 M. Caputo and F. Mainardi, “A new dissipation model based on memory mechanism,” Pure and Applied Geophysics PAGEOPH, vol. 91, no. 1, pp. 134–147, 1971. View at: Publisher Site  Google Scholar
 E. Emmanuel and K. Reza, “On fractionalorder PID controllers,” IFAC PapersOnLine, vol. 51, no. 4, pp. 739–744, 2018. View at: Google Scholar
 T. Sathiyaraj and P. Balasubramaniam, “Controllability of fractional higher order stochastic integrodifferential systems with fractional brownian motion,” ISA Transactions, vol. 82, pp. 107–119, 2017. View at: Publisher Site  Google Scholar
 G. M. Mahmoud, E. E. Mahmoud, and M. E. Ahmed, “A hyperchaotic complex chen system and its dynamics,” International Journal of Applied Mathematics and Statistics, vol. 12, no. 7, pp. 90–100, 2007. View at: Google Scholar
 Y. Xu, W. Xu, and G. M. Mahmoud, “On a complex beambeam interaction model with random forcing,” Physica A: Statistical Mechanics and its Applications, vol. 336, no. 34, pp. 347–360, 2004. View at: Publisher Site  Google Scholar
 G. M. Mahmoud and S. A. Aly, “Periodic attractors of complex damped nonlinear systems,” International Journal of NonLinear Mechanics, vol. 35, no. 2, pp. 309–323, 2000. View at: Publisher Site  Google Scholar
 G. M. Mahmoud, “Stability regions for coupled Hill’s equations,” Physica A: Statistical Mechanics and its Applications, vol. 242, no. 12, pp. 239–249, 1997. View at: Publisher Site  Google Scholar
 K. W. Chung, C. L. Chan, Z. Xu, and G. M. Mahmoud, “A perturbationincremental method for strongly nonlinear autonomous oscillators with many degrees of freedom,” Nonlinear Dynamics, vol. 28, no. 34, pp. 243–259, 2002. View at: Publisher Site  Google Scholar
 G. M. Mahmoud and S. A. H. Aly, “On periodic solutions of parametrically excited complex nonlinear dynamical systems,” Physica A: Statistical Mechanics and its Applications, vol. 278, no. 34, pp. 390–404, 2000. View at: Publisher Site  Google Scholar
 G. M. Mahmoud, “Periodic solutions of strongly nonlinear Mathieu oscillators,” International Journal of NonLinear Mechanics, vol. 32, no. 6, pp. 1177–1185, 1997. View at: Publisher Site  Google Scholar
 G. M. Mahmoud, E. E. Mahmoud, and A. A. Arafa, “Passive control of ndimensional chaotic complex nonlinear systems,” Journal of Vibration and Control, vol. 19, no. 7, pp. 1061–1071, 2013. View at: Publisher Site  Google Scholar
 G. M. Mahmoud, A. A. Arafa, T. M. Abed–Elhameed, and E. E. Mahmoud, “Chaos control of integer and fractional orders of chaotic burkeshaw system using time delayed feedback control,” Chaos, Solitons and Fractals, vol. 104, pp. 680–692, 2017. View at: Publisher Site  Google Scholar
 M. M. AlSawalha and A. AlSawalha, “Antisynchronization of fractional order chaotic and hyperchaotic systems with fully unknown parameters using modified adaptive control,” Open Physics, vol. 14, no. 1, pp. 304–313, 2016. View at: Publisher Site  Google Scholar
 Y. Shen, S. Yang, H. Xing, and G. Gao, “Primary resonance of Duffing oscillator with fractionalorder derivativefing oscillator with fractionalorder derivative,” Communications in Nonlinear Science and Numerical Simulation, vol. 17, no. 7, pp. 3092–3100, 2012. View at: Publisher Site  Google Scholar
 K. Mobin and S. Saeed, “Dynamic behaviors of a fractional order nonlinear oscillator,” Journal of King Saud University: Science, vol. 31, no. 1, pp. 14–20, 2019. View at: Google Scholar
 S. F. Wen, Y. J. Shen, X. N. Wang, S. P. Yang, and H. J. Xing, “Dynamical analysis of strongly nonlinear fractionalorder MathieuDuffing equation,” Chaos: An Interdisciplinary Journal of Nonlinear Science, vol. 26, no. 8, Article ID 084309, 2016. View at: Publisher Site  Google Scholar
 Y.J. Shen, S.F. Wen, X.H. Li, S.P. Yang, and H.J. Xing, “Dynamical analysis of fractionalorder nonlinear oscillator by incremental harmonic balance method,” Nonlinear Dynamics, vol. 85, no. 3, pp. 1457–1467, 2016. View at: Publisher Site  Google Scholar
 Y. Xu, W. Xu, and G. M. Mahmoud, “On a complex duffing system with random excitation,” Chaos, Solitons & Fractals, vol. 35, no. 1, pp. 126–132, 2008. View at: Publisher Site  Google Scholar
 L. Chen and W. Zhu, “Stochastic jump and bifurcation of duffing oscillator with fractional derivative damping under combined harmonic and white noise excitations,” International Journal of NonLinear Mechanics, vol. 46, no. 10, pp. 1324–1329, 2011. View at: Publisher Site  Google Scholar
 L. C. Chen, T. L. Zhao, W. Li, and J. Zhao, “Bifurcation control of bounded noise excited duffing oscillator by a weakly fractionalorder PIλDμ feedback controller,” Nonlinear Dynamics, vol. 83, no. 12, pp. 529–539, 2016. View at: Publisher Site  Google Scholar
 L. C. Chen, F. Hu, and W. Q. Zhu, “Stochastic dynamics and fractional optimal control of quasi integrable Hamiltonian systems with fractional derivative damping,” Fractional Calculus and Applied Analysis, vol. 16, no. 1, pp. 189–225, 2013. View at: Publisher Site  Google Scholar
 X. J. Liu, L. Hong, and J. Jiang, “Crises in a nonautonomous fractionalorder duffing system,” Acta Physica Sinica, vol. 65, no. 18, Article ID 180502, 2016. View at: Google Scholar
 X. J. Liu, L. Hong, and J. Jiang, “Global bifurcations in fractionalorder chaotic systems with an extended generalized cell mapping method,” Chaos: An Interdisciplinary Journal of Nonlinear Science, vol. 26, no. 8, Article ID 084304, 2016. View at: Publisher Site  Google Scholar
 P. J. Zufiria and S. Jiménez, “Characterizing chaos in a type of fractional duffing’s equation,” Discrete & Continuous Dynamical Systems, vol. 35, pp. 660–669, 2015. View at: Google Scholar
 A. R. Mohammad, S. Hassan, and A. Aria, “Stabilizing periodic orbits of fractional order chaotic systems via linear feedback theory,” Applied Mathematical Modelling, vol. 36, no. 3, pp. 863–877, 2011. View at: Google Scholar
 M. Buslowicz and A. Makarewicz, “Chaos synchronization of the modified van der polduffing oscillator of fractional order,” Recent Advances in Automation, Robotics and Measuring Techniques, vol. 267, pp. 33–43, 2014. View at: Google Scholar
 T. A. Giresse and K. T. Crépin, “Chaos generalized synchronization of coupled mathieuvan der pol and coupled duffingvan der pol systems using fractional orderderivative,” Chaos, Solitons & Fractals, vol. 98, no. 33, pp. 88–100, 2017. View at: Publisher Site  Google Scholar
 A. E. Matouk, “Chaos, feedback control and synchronization of a fractionalorder modified autonomous van der polduffing circuit,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 2, pp. 975–986, 2011. View at: Publisher Site  Google Scholar
 M.X. Cai, J.P. Yang, and J. Deng, “Bifurcations and chaos in Duffing equation with damping and external excitations,” Acta Mathematicae Applicatae Sinica, English Series, vol. 30, no. 2, pp. 483–504, 2014. View at: Publisher Site  Google Scholar
 Y.J. Shen, S.F. Wen, S.P. Yang, S.Q. Guo, and L.R. Li, “Analytical threshold for chaos in a duffing oscillator with delayed feedbacks,” International Journal of NonLinear Mechanics, vol. 98, pp. 173–179, 2018. View at: Publisher Site  Google Scholar
 S. F. Wen, J. F. Chen, and S. Q. Guo, “Heteroclinic bifurcation behaviors of a duffing oscillator with delayed feedback,” Shock and Vibration, vol. 2018, Article ID 7213606, 12 pages, 2018. View at: Publisher Site  Google Scholar
 Y.G. Wang, H.F. Song, D. Li, and J. Wang, “Bifurcations and chaos in a periodic timevarying temperatureexcited bimetallic shallow shell of revolution,” Archive of Applied Mechanics, vol. 80, no. 7, pp. 815–828, 2010. View at: Publisher Site  Google Scholar
 L. Du, Y. Zhao, Y. Lei, J. Hu, and X. Yue, “Suppression of chaos in a generalized Duffing oscillator with fractionalorder deflection,” Nonlinear dynamics, vol. 92, no. 4, pp. 1921–1933, 2018. View at: Publisher Site  Google Scholar
 L. M. Anague Tabejieu, B. R. Nana Nbendjo, and P. Woafo, “On the dynamics of Rayleigh beams resting on fractionalorder viscoelastic Pasternak foundations subjected to moving loads,” Chaos, Solitons & Fractals, vol. 93, pp. 39–47, 2016. View at: Publisher Site  Google Scholar
 H. Hu, Applied Nonlinear Dynamics, Aviation Industry Press, Beijing, China, 2000.
 J. Guckenheimer and P. Holmes, Nonlinear Oscillations, Dynamical System and Bifurcations of Vector Fields, SpringerVerlag, New York, NY, USA, 1983.
 P. Holmes, “A nonlinear oscillator with a strange attractor,” Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, vol. 292, no. 1394, pp. 419–448, 1979. View at: Publisher Site  Google Scholar
 P. P. Bruce, An Introduction to Complex Function Theory, SpringerVerlag, New York, NY, USA, 1991.
 I. Petras, FractionalOrder Nonlinear Systems: Modeling, Analysis and Simulation, Higher Education Press, Beijing, China, 2011.
 D. Y. Xue, Fractional Calculus and FractionalOrder Control, Science Press, Beijing, China, 2018.
 J. C. Sprott, “A simple chaotic delay differential equation,” Physics Letters A, vol. 366, no. 45, pp. 397–402, 2007. View at: Publisher Site  Google Scholar
 M.F. Danca and N. Kuznetsov, “Matlab code for Lyapunov exponents of fractionalorder systems,” International Journal of Bifurcation and Chaos, vol. 28, no. 5, Article ID 1850067, 2018. View at: Publisher Site  Google Scholar
 J.H. Jia, X.Y. Shen, and H.X. Hua, “Viscoelastic behavior analysis and application of the fractional derivative maxwell model,” Journal of Vibration and Control, vol. 13, no. 4, pp. 385–401, 2007. View at: Publisher Site  Google Scholar
 Y. Yang and W. Zhao, “Curvelet transformbased identification of void diseases in ballastless track by groundpenetrating radar,” Structural Control and Health Monitoring, vol. 29, no. 4, p. e2322, 2019. View at: Google Scholar
 A. Shokooh and L. Suárez, “A comparison of numerical methods applied to a fractional model of damping materials,” Journal of Vibration and Control, vol. 5, no. 3, pp. 331–354, 1999. View at: Publisher Site  Google Scholar
 R. S. Barbosa and J. A. T. Machado, “Describing function analysis of systems with impacts and backlash,” Nonlinear Dynamics, vol. 29, no. 1–4, pp. 235–250, 2002. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2019 Wuce Xing et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.