# Chaos Analysis and Control of Relative Rotation System with Mathieu-Duffing Oscillator

Authors: Yu Zhang; Longsuo Li
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

Chaos analysis and control of relative rotation nonlinear dynamic system with Mathieu-Duffing oscillator are investigated. By using Lagrange equation, the dynamics equation of relative rotation system has been established. Melnikov’s method is applied to predict the chaotic behavior of this system. Moreover, the chaotic dynamical behavior can be controlled by adding the Gaussian white noise to the proposed system for the sake of changing chaos state into stable state. Through numerical calculation, the Poincaré map analysis and phase portraits are carried out to confirm main results.


---

## Body

## 1. Introduction

Chaos, as a kind of physical phenomenon, exists widely in various nonlinear dynamic systems [1–5]. The main ideal of chaos control approach is that one can use given unstable periodic orbits or external excitations to force chaotic system into stable system [6–9]. There are many ways to be used to suppress chaos, such as parametric driven [10–12] and time delayed approach [13, 14], and the Gaussian white noise as random phase control is an interesting one. The control of chaos we discuss in this paper is based on this case.The theory of mechanics of relative rotation system was established by Carmeli in 1986 [15]; bifurcation and chaos of this system got rapid development in variety of areas. Shi et al. [16] studied the chaotic behavior and its control for a class of nonlinear dynamics equation of relative rotation system. Liu et al. [17] found stability and bifurcation for a coupled nonlinear relative rotation system with multitime delay feedbacks. Therefore, this system plays an important role in nonlinear dynamic systems.Stochastic forces or random noise are the most frequently control strategies to be used in suppressing chaos. Ramesh and Narayanan [18] investigated the robustness in presence of uniform noise and found that the system would lose control while noise intensity was raided to a threshold level. Liu et al. [19] explored the effect of bounded noise on chaotic motion of Duffing oscillator under parametric excitation. Wu et al. [20] studied stochastic chaos and its control by the top Lyapunov exponent. Yin et al. [21] proved the peculiar solitary waves are more likely to be chaos by using the Melnikov theory and found that the system can be well controlled when the frequency of the perturbation surpasses the peculiar perturbation frequency with fixed parameters of the unperturbed system. Noise, as random phase, has been used in studying the control of chaos in the paper. So as a kind of effective method of chaos control, whether in theory or in the practical application, the noise has some of the research significance.In this paper, we explore chaos analysis and control in relative rotation system with Mathieu-Duffing oscillator. The paper is organized as follows. The dynamics equation of relative rotation system is established in Section2. In Section 3, the prediction of chaotic motion is given by Melnikov’s method. In Section 4, the Gaussian white noise can increase or decrease the threshold values of chaos, so as to realize the control of chaos.


## 2. Dynamics Equation of Relative Rotation System with Mathieu-Duffing Oscillator

Mathieu-Duffing system is a class of typical vibration system, and researching its nonlinear dynamic characteristics shows to be extremely essential. Considering Mathieu-Duffing oscillator expression form as follows [22, 23](1)ft=N1tx+N2tx3,where (2)N1t=cKt,N2t=bKt,Kt=k02+k1cos⁡⁡ω1tN1(t), N2(t) are linear torsional stiffness of the system. For two quality relative rotation systems of Mathieu-Duffing oscillator, the kinetic energy is(3)E=12J1θ˙12+12J2θ˙22,where J1, J2 are moment of inertia and θi, θ˙i(i=1,2) are rotational angle and speed of rotational angle, respectively. The potential energy of system is(4)U=12N1tθ2+12N2tθ4.Generalized forces are(5)Q1=F11∂θ1∂q1+F22∂θ2∂q1,Q2=F11∂θ1∂q2+F22∂θ2∂q2,where (6)F11=T1-Aθ˙1-θ˙2,F22=T2-Aθ˙2-θ˙1,where T1, T2 are generalized external force, q1, q2 are generalized coordinates, and A is linear damping coefficient. The Lagrange equation is(7)ddt∂E∂q˙j-∂E∂qj+∂U∂qj.Substituting (1)–(5) into Lagrange equation (7), we have(8)J1θ¨1+cKtθ1-θ2+bKtθ1-θ23+Aθ˙1-θ˙2=T1,(9)J2θ¨2+cKtθ2-θ1+bKtθ2-θ13+Aθ˙2-θ˙1=T2,where θ¨1, θ¨2 are acceleration of rotational angle of moment of inertia. Combining (8) and (9), we have (10)θ¨1-θ¨2+J1+J2J1J2cKtθ1-θ2+J1+J2J1J2bKtθ1-θ23+J1+J2J1J2Aθ˙1-θ˙2=J2T1-J1T2J1J2.Let (11)x=θ1-θ2,x˙=θ˙1-θ˙2,x¨=θ¨1-θ˙2,J2T1-J1T2J1J2=Fcos⁡⁡ω2t,a=J1+J2J1J2,λ=J1+J2J1J2A.Hence, the above equation turns to(12)x¨+ak02+k1cos⁡⁡ω1tcx+bx3+λx˙=Fcos⁡⁡ω2t.Equation (12) is dynamics equation of relative rotation system with Mathieu-Duffing oscillator.


## 3. Melnikov’s Method Analysis

Melnikov’s method [24–26] is an effective approach to predict chaotic behavior, its basic idea mainly makes the dynamics system into a Poincaré map system, and we obtain chaotic properties by studying conditions of the mapping system whether there exist homoclinic orbits or heteroclinic orbits.Equation (12) can be transformed into first-order equation(13)x˙1=x2,x˙2=-ack02x1-abk02x13-ak1cos⁡⁡ω1tcx1+bx13-λx2+Fcos⁡⁡ω2t.Small parameter ɛ is added to the nonlinear term of (13) and it may be written as(14)x˙1=x2,x˙2=-ack02x1-abk02x13-ɛak1cos⁡⁡ω1tcx1+bx13+λx2-Fcos⁡⁡ω2t.When ɛ=0, the above equation turns into(15)x˙1=x2,x˙2=-ack02x1-abk02x13.Equation (15) is a unperturbed Hamilton system and its Hamiltonian is(16)Hx1,x2=12x22+12ack02x12+14abk02x14.There exist heteroclinic orbits and they satisfy(17)x˙1=x2,H±-cb,0=12x22+12ack02x12+14abk02x14.Then the parametric equations of two heteroclinic orbits are(18)x1t=±-cbtanh⁡2acak0t,x2t=±ck02-2ab1cosh⁡22acak0t.Melnikov function can be given by (19)Mt0=∫-∞+∞x2t·-ak1cos⁡⁡ω1t+t0cx1t+bx1t3-λx2t+Fcos⁡⁡ω2t+t0dt=±ack1∫-∞+∞x1tx2tcos⁡⁡ω1t+t0dt±abk1∫-∞+∞x1t3x2tcos⁡⁡ω1t+t0dt-λ∫-∞+∞x2t2dt±F∫-∞+∞x2tcos⁡⁡ω2t+t0dt,where (20)∫-∞+∞x1tx2tcos⁡⁡ω1t+t0dt=∫-∞+∞x1tx2tsin⁡ω1tdtsin⁡ω1t0,∫-∞+∞x1t3x2tcos⁡⁡ω1t+t0dt=∫-∞+∞x1t3x2tsin⁡ω1tdtsin⁡ω1t0,∫-∞+∞x2tcos⁡⁡ω2t+t0dt=∫-∞+∞x2tcos⁡⁡ω2tdtcos⁡⁡ω2t0.Let (21)Z1=∫-∞+∞x2t2dt,Z2=∫-∞+∞x1tx2tsin⁡ω1tdt,Z3=∫-∞+∞x1t3x2tsin⁡ω1tdt,Z4=∫-∞+∞x2tcos⁡⁡ω2tdt.Substituting Z1-Z4 into Melnikov function, we have(22)Mt0=-λZ1±ack1Z2+abk1Z3sin⁡ω1t0±FZ4cos⁡⁡ω2t.When ω1=ω2=ack02, (22) can be written as(23)Mt0=-λZ1±ack1Z2+abk1Z32+FZ42sin⁡ω1t0+φ,where (24)cos⁡⁡φ=ack1Z2+abk1Z3ack1Z2+abk1Z32+FZ42,sin⁡φ=FZ4ack1Z2+abk1Z32+FZ42.Because of M(t0)=0, we have (25)sin⁡ω1t0+φ=±λZ1ack1Z2+abk1Z32+FZ42.Due to |sin⁡(ω1t0+φ)|≤1, then (26)λZ1ack1Z2+abk1Z32+FZ42≤1,∂Mt0∂t0=±ω1ack1Z2+abk1Z32+FZ42·cos⁡⁡ω1t0+φ≠0.That is, cos⁡(ω1t0+φ)≠0, so sin(ω1t0+φ)≠0; we can know (27)λZ1ack1Z2+abk1Z32+FZ42<1.There exists a sufficient small ɛ, making M(t0)=0 and ∂M(t0)/∂t0≠0; therefore, the system produces chaotic behavior.For system (13), taking parameters ack02=4, abk02=-4, ack1=1, abk1=-1, λ=0.5, F=1.42, and ω1=ω2=2 with the initial conditions x(0)=1, x˙(0)=0, the Poincaré map and phase portraits were plotted by using the numerical simulation in Figure 1.Figure 1
(a) Poincaré map; (b) phase portraits.
(a)
(b)From Figure1, we can see that the Poincaré map has chaotic attractors, phase portraits show nonoverlapping and are disorganized, and these can prove that relative rotation system with Mathieu-Duffing oscillator is chaos.


## 4. Chaos Control of Relative Rotation System with Mathieu-Duffing Oscillator Using Noise

Adding Gaussian white noise, (13) can be written as(28)x˙1=x2,x˙2=-ack02x1-abk02x13-ak1cos⁡⁡ω1tcx1+bx13-λx2+Fcos⁡⁡ω2t+σξt,where ξ(t) denotes standard Gaussian white noise; it satisfies (29)Eξt=0,Eξtξt+τ=ζτ,where ζ(τ) is Dirac-Delta function and σ is the intensity of noise.Taking the same parameters and initial conditions as Section3. Equation (28) can be written:(30)x˙1=x2,x˙2=-4x1+4x13-cos⁡⁡2tx1+cos⁡⁡2tx13-0.5x2+1.42cos⁡⁡2t+σξt.Making the Poincaré map as(31)Σ⟶Σ,Σxt,x˙t∣t=0,2πω2,4πω2,…⊆R2.Under the given initial conditions, the solution of differential equation (30) is solved by the fourth-order Runge-Kutta method and the solution is plotted for every T=2π/ω2. After deleting the first 100 dots, we plot the Poincaré map by using surplus 200 iteration dots, when the intensity of noise σ=0.4; it is shown in Figure 2(a).Figure 2
(a) Poincaré map withσ=0.4; (b) phase portraits with σ=0.4.
(a)
(b)The Poincaré map appears to be a stable stochastic attractor, which means that the system is stable. To further verify the obtained results, we plotted the phase portraits. Taking the same intensity of noise, the phase portraits are shown in Figure2(b). We can see from Figure 2(b) that the phase portraits turned into a regular annulus. It is proved that chaotic behavior was suppressed, and the system changed chaos state into stable state.


## 5. Concluding Remarks

We investigate the chaos analysis and control of relative rotation system with Mathieu-Duffing oscillator in this paper. By analysis of Melnikov’s method, we confirm chaotic behavior of this system under the given parameters and initial conditions. Adding the Gaussian white noise to phase of system, we plot the Poincaré map and phase portraits under the intensity of noiseσ=0.4. We get that the Poincaré map changed chaotic attractors into a stable attractor, and phase portraits also displayed regular state, so the relative rotation system with Mathieu-Duffing oscillator eventually becomes stable.



---

*Source: 348462-2015-04-30.xml*