# Chaotic Motion in Forced Duffing System Subject to Linear and Nonlinear Damping

Authors: Tai-Ping Chang
Journal: Mathematical Problems in Engineering
 (2017)

---

## Abstract

This paper investigates the chaotic motion in forced Duffing oscillator due to linear and nonlinear damping by using Melnikov technique. In particular, the critical value of the forcing amplitude of the nonlinear system is calculated by Melnikov technique. Further, the top Lyapunov exponent of the nonlinear system is evaluated by Wolf’s algorithm to determine whether the chaotic phenomenon of the nonlinear system actually occurs. It is concluded that the chaotic motion of the nonlinear system occurs when the forcing amplitude exceeds the critical value, and the linear and nonlinear damping can generate pronounced effects on the chaotic behavior of the forced Duffing oscillator.


---

## Body

## 1. Introduction

The study of chaotic motion in nonlinear systems has been a popular area of research during the last few decades. Many investigations have been performed on different nonlinear chaotic systems to understand the complex behavior of these systems. Three of the fundamental forced oscillators, Duffing, Van der Pol, and Rayleigh oscillators, have been extensively examined since lots of dynamic characteristics embedded in the physical systems can be realized from these three systems [1–3]. Among them, forced Duffing oscillator is the most useful nonlinear dynamical systems, which is considered as a prototype model for various physical and engineering problems such as dynamics of a buckled elastic beam, particle in a forced double well, particle in a plasma, and a defect in solids [4]. More recently, Wang et al. [5–7] show that many nonlinear dynamical problems can be reduced to Duffing system for structures subject to both free and forced excitation. Therefore, the Duffing oscillator is still an interesting model to study for discovering the characteristics of chaos in nonlinear physical systems despite lots of investigations have been performed on this model for so many years. The Duffing oscillator can be interpreted as a damped oscillator with a complicated potential. The damping or dissipation here is very important as it decides the border of stability and instability. In the past, most of the studies on Duffing system have been performed by considering the linear viscous damping. Several researches have been reported on the chaotic motion of forced Duffing equations [8–10]. However, the consideration of nonlinear damping is quite necessary in various engineering applications such as drag forces in flow induced vibrations [4] and vibration isolators [11]. Among others, some researches have made the contributions to the chaotic behavior of Duffing oscillator due to nonlinear damping [12–21].Generally speaking, the Melnikov method is very useful for detecting the presence of transverse homoclinic orbits and the occurrence of homoclinic bifurcations. However, the traditional Melnikov methods strongly depend on the small perturbation parameters so that these methods are limited in coping with the systems with strong nonlinearities. Liu et al. [22] presents a procedure to investigate the chaos and subharmonic resonance of strongly nonlinear practical systems by using a homotopy method that is adopted to extend the Melnikov functions to the strongly nonlinear systems.In the present study, the chaotic behavior of a forced Duffing oscillator subjected to linear and nonlinear damping is investigated by using Melnikov technique [23, 24]. In particular, the critical values of the forcing amplitude of the nonlinear system is calculated by Melnikov technique. The Lyapunov exponents of the nonlinear system are evaluated by Wolf’s algorithm [25] to determine whether the chaotic phenomenon of the nonlinear system actually occurs.Although the Melnikov method is used here only for Duffing system without the control action, however, it is known that control action is an important direction of the nonlinear dynamics field. Among others, recently Wang and Li [26] investigated the nonlinear dynamical characteristics of the Duffing–Van der Pol oscillator subject to both external and parametric excitations with time delayed feedback control by using the multiple scale method. Actually, the Melnikov method can also be used for Duffing system with the feedback control. Yagasaki [27] considered a pendulum subjected to linear feedback control with periodic desired motions. He studied local bifurcations of harmonics and subharmonics using the second-order averaging method and Melnikov’s method.


## 2. Governing Equation of Nonlinear System

In many engineering and physics fields, there are lots of models can be converted into the following forced Duffing equation:(1)x¨-ω02x+Dx˙+γx˙3+βx3=Fcos⁡ωt, where -ω02 denotes the linear stiffness constant of the system, D is linear damping coefficient, γ is the nonlinear damping coefficient, β is the nonlinear stiffness constant, and F and ω are the external amplitude and driving frequency, respectively. In the present study, D,γ, and F are assumed to be small parameters. Hence, a transformation D→εD,γ→εγ,F→εF is adopted in order to apply the first-order perturbation scheme of the Melnikov method. Therefore, (1) can be rewritten as(2)x˙=y,y˙=ω02x-βx3+ε-Dy-γy3+Fcos⁡ωt.When ε=0, an unperturbed system can be obtained as follows:(3)x˙=y,y˙=ω02x-βx3.The system of (3) corresponds to a Hamiltonian system with a potential function(4)Vx=-ω022x2+β4x4.The unperturbed system of (4) has three equilibrium points: one saddle 0,0 and two centers (±ω0/β,0).


## 3. Chaotic Motion of Forced Duffing Oscillator

By integrating system (3), the solution of unperturbed system (3) can be written as follows:(5)x0t=2βω0sech⁡ω0t,x˙0t=y0t=2βω02sech⁡ω0ttanh⁡ω0t.Now we can use Melnikov’s method [23, 24] to investigate the homoclinic bifurcation in the forced Duffing oscillator system with linear and nonlinear damping as given in (1). The Melnikov function measures the distance between the stable and unstable manifolds in the Poincare section, and to preserve the homoclinic loops under a perturbation requires that, at t0, if M(t0) has a simple zero, then a homoclinic bifurcation occurs, implying that the chaotic motion occurs. The Melnikov function for the forced Duffing oscillator system shown in (1) can be obtained as follows:(6)Mt0=-D∫-∞∞y02dt-γ∫-∞∞y04dt+F∫-∞∞y0cos⁡ωt-t0dt,where t0 is the cross-section time of the Poincare map and t0 can be interpreted as the initial time of the forcing term. By substituting (5) into (6) and calculating the integral, we can get the Melnikov function as follows:(7)Mt0=-DI1-γI2+FI3sin⁡ωt0,where(8)I1=2ω04β∫-∞∞sech2⁡ω0ttanh2⁡ω0tdt=4ω033β,I2=4ω08β2∫-∞∞sech4⁡ω0ttanh4⁡ω0tdt=16ω0735β2,I3=-2βω0∫-∞∞sech⁡ω0ttanh⁡ω0tsin⁡ωtdt=-πω2ββsech⁡πω2ω0.The simple zeros of (7) give the critical value of the forcing amplitude (Fc), and the first homoclinic bifurcation occurs when we cross the critical value by increasing the forcing amplitude, implying that the chaotic motion occurs in the system.In (7), M(t0) will vanish if a solution can be detected for t0, and sin(ω0t)≤1; therefore, the critical value of the forcing amplitude (Fc) can be obtained as follows: (9)F=Fc=DI1+γI2I3,where I1,I2,I3 are given in (8).The main purpose of the present study is to investigate the chaotic motion of the nonlinear differential equation shown in (1) by evaluating the critical value of the forcing amplitude given in (9). The top Lyapunov exponent of the nonlinear system is evaluated by Wolf’s algorithm [25] to check whether the chaotic phenomenon of the nonlinear system occurs. It is noticed that when the top Lyapunov exponent changes from negative to positive, the chaotic motion of the nonlinear system happens.


## 4. Numerical Examples and Discussions

First of all, the numerical values of the parameters in (1) are fixed as follows:(10)ω0=1,β=1,γ=0,ω=1.In addition, the initial conditions in (1) are considered as follows:(11)x0=1.0,x˙0=0.0.Based on (9), the variations of critical values of forcing amplitude (Fc) versus linear damping coefficient D, driving frequency ω, and nonlinear coefficient β are depicted in Figures 1(a)–1(c). As can be seen from Figure 1(a), Fc increases linearly with the linear damping coefficient; however, Fc varies nonlinearly with respect to the driving frequency as detected from Figure 1(b); first it decreases with the increase in ω and after reaching a minimum, it begins increasing. In Figure 1(c), the critical values of forcing amplitude (Fc) decreases nonlinearly with respect to the nonlinear coefficient β. Now only the linear damping effect of the system is considered, and the parameters are fixed as follows: ω0=1, β=1, D=0.4, γ=0, ω=1. As can be seen from Figures 1(a)–1(c), the critical values of forcing amplitude (Fc) are detected as 0.3012; hence, 0.305 is used for the numerical computations shown in Figures 3(a)–3(d) and 4(a)–4(d).Figure 1
(a) Critical value of forcing amplitudeFc versus linear damping coefficient (ω0=1,β=1,γ=0,ω=1). (b) Critical value of forcing amplitude Fc versus driving frequency (ω0=1,β=1,D=0.4,γ=0). (c) Critical value of forcing amplitude Fc versus nonlinear coefficient ω0=1,D=0.4,γ=0,ω=1.
(a) (b) (c)Now we include the nonlinear damping effect, that is,γ=0.1; while the other parameters remain the same as before, as can be seen in Figures 2(a)–2(c), the variations of critical values of forcing amplitude (Fc) versus linear damping coefficient D, driving frequency ω, and nonlinear coefficient β are presented. It is noticed that the similar phenomena can be found by including the nonlinear damping effect except that the critical values of forcing amplitude (Fc) are larger, which is quite reasonable.Figure 2
(a) Critical value of forcing amplitudeFc versus linear damping coefficient D(ω0=1, β=1, γ=0.1, ω=1). (b) Critical value of forcing amplitude Fc versus driving frequency (ω0=1, β=1, D=0.4, γ=0.1). (c) Critical value of forcing amplitude Fc versus nonlinear coefficient β(ω0=1, D=0.4, γ=0.1, ω=1).
(a) (b) (c)Figure 3
(a) Top Lyapunov exponents versus time(ω0=1, β=1, D=0.4, γ=0, ω=1, F=0.305). (b) Time history of nonlinear system (ω0=1, β=1, D=0.4, γ=0, ω=1, F=0.305). (c) Phase portrait of nonlinear system (ω0=1, β=1, D=0.4, γ=0, ω=1, F=0.305). (d) Poincare map of nonlinear system (ω0=1, β=1, D=0.4, γ=0, ω=1, F=0.305).
(a) (b) (c) (d)Figure 4
(a) Top Lyapunov exponents versus time(ω0=1, β=1, D=0.4, γ=0.1, ω=1, F=0.305). (b) Time history of nonlinear system (ω0=1, β=1, D=0.4, γ=0.1, ω=1, F=0.305). (c) Phase portrait of nonlinear system (ω0=1, β=1, D=0.4, γ=0.1, ω=1, F=0.305). (d) Poincare map of nonlinear system (ω0=1, β=1, D=0.4, γ=0.1, ω=1, F=0.305).
(a) (b) (c) (d)In Figures3(a)–3(d), we plot the top Lyapunov exponents, the time history, the phase portrait, and Poincare map of nonlinear system. As can be detected from Figure 3(a), the top Lyapunov exponent is computed as λ1=0.0088301 which is positive, and the chaotic motion of the system will occur as shown in Figure 3(d). However, if the nonlinear damping effect is considered, namely, γ=0.1, the other parameters remain the same; the top Lyapunov exponents, the time history, the phase portrait, and Poincare map of nonlinear system are presented in Figures 4(a)–4(d). As can be seen from Figure 4(a), the top Lyapunov exponent turns out to be λ1=-0.052549, which is negative, and the periodic motion of the system can be detected as presented in Figure 4(d). Now let us use the numerical results shown in Figures 2(a)–2(c); the critical value of forcing amplitude (Fc) of the nonlinear system with both linear and nonlinear damping is found as 0.3270; hence, 0.400 is adopted for the numerical computations shown in Figures 5(a)–5(d).Figure 5
(a) Top Lyapunov exponents versus time(ω0=1, β=1, D=0.4, γ=0.1, ω=1, F=0.4). (b) Time history of nonlinear system (ω0=1, β=1, D=0.4, γ=0.1, ω=1, F=0.4). (c) Phase portrait of nonlinear system (ω0=1, β=1, D=0.4, γ=0.1, ω=1, F=0.4). (d) Poincare map of nonlinear system (ω0=1, β=1, D=0.4, γ=0.1, ω=1, F=0.4).
(a) (b) (c) (d)In Figures5(a)–5(d), the top Lyapunov exponents, the time history, the phase portrait, and Poincare map of nonlinear system are depicted under forcing amplitude F=0.400. The top Lyapunov exponent now changes from negative to positive as detected from Figure 5(a); as a result, the chaotic motion of the nonlinear system occurs again as shown in Figure 5(d), which is not surprising.Now we switch the driving frequency of the system fromω=1.0 to ω=2.0 and consider first the linear damping effect only; in this case the following parameters are adopted: ω0=1, β=1, D=0.4, γ=0, ω=2. As can be seen from Figure 1(b), the critical values of forcing amplitude (Fc) is found as 0.6599; hence, 0.800 is used for the numerical computations shown in Figures 6(a)–6(d). In Figures 6(a)–6(d), we plot the top Lyapunov exponents, the time history, the phase portrait, and Poincare map of nonlinear system. As can be seen from Figure 6(a), the top Lyapunov exponent is computed as λ1=0.044337 which is positive; the chaotic motion of the system will occur as shown in Figure 6(d). Now if the nonlinear damping effect is considered, namely, γ=0.1, the other parameters remain the same; the top Lyapunov exponents, the time history, the phase portrait, and Poincare map of nonlinear system are presented in Figures 7(a)–7(d). As can be seen from Figure 7(a), the top Lyapunov exponent turns out to be λ1=-0.01584, which is negative, and the periodic motion of the system can be seen as depicted in Figure 7(d).Figure 6
(a) Top Lyapunov exponents versus time(ω0=1, β=1, D=0.4, γ=0.0, ω=2, F=0.8). (b) Time history of nonlinear system (ω0=1, β=1, D=0.4, γ=0.0, ω=2, F=0.8). (c) Phase portrait of nonlinear system (ω0=1, β=1, D=0.4, γ=0.0, ω=2, F=0.8). (d) Poincare map of nonlinear system (ω0=1, β=1, D=0.4, γ=0.0, ω=2, F=0.8).
(a) (b) (c) (d)Figure 7
(a) Top Lyapunov exponents versus time(ω0=1, β=1, D=0.4, γ=0.1, ω=2, F=0.8). (b) Time history of nonlinear system (ω0=1, β=1, D=0.4, γ=0.1, ω=2, F=0.8). (c) Phase portrait of nonlinear system (ω0=1, β=1, D=0.4, γ=0.1, ω=2, F=0.8). (d) Poincare map of nonlinear system (ω0=1, β=1, D=0.4, γ=0.1, ω=2, F=0.8).
(a) (b) (c) (d)Finally we choose the forcing amplitude asF=1.0 which is larger than the critical values of forcing amplitude (Fc) based on the numerical results shown in Figure 2(b). Then the top Lyapunov exponents, the time history, the phase portrait, and Poincare map of nonlinear system are depicted in Figures 8(a)–8(d). The top Lyapunov exponent is estimated as λ1=0.19182 shown in Figure 8(a), which is positive, and the chaotic motion of the system will occur as shown in Figure 8(d).Figure 8
(a) Top Lyapunov exponents versus time(ω0=1, β=1, D=0.4, γ=0.1, ω=2, F=1.0). (b) Time history of nonlinear system (ω0=1, β=1, D=0.4, γ=0.1, ω=2, F=1.0). (c) Phase portrait of nonlinear system (ω0=1, β=1, D=0.4, γ=0.1, ω=2, F=1.0). (d) Poincare map of nonlinear system (ω0=1, β=1, D=0.4, γ=0.1, ω=2, F=1.0).
(a) (b) (c) (d)Based on all the numerical results mentioned above, it is understood that the linear and nonlinear damping can have significant effects on the chaotic behavior of the nonlinear system. That is, the linear and nonlinear damping can influence the nonlinear system into or out of a chaotic regime.


## 5. Conclusions

In the present study, the chaotic behavior of a forced Duffing oscillator subjected to linear and nonlinear damping is investigated by using Melnikov technique. In particular, the critical value of the forcing amplitude of the nonlinear system is calculated by Melnikov technique. The top Lyapunov exponent of the nonlinear system is evaluated by Wolf’s algorithm to determine whether the chaotic phenomenon of the nonlinear system actually occurs. As a result, the critical value of forcing amplitude increases linearly with the linear damping coefficient; however, it varies nonlinearly with respect to the driving frequency. Further, the critical value of forcing amplitude decreases nonlinearly with respect to the nonlinear coefficient. To detect the chaotic phenomena of the nonlinear system, the top Lyapunov exponent, the time history, the phase portrait, and Poincare map of the nonlinear system are presented for various cases. It is concluded that the chaotic motion of the nonlinear system with the linear and nonlinear damping occurs when the forcing amplitude exceeds the critical value, and the linear and nonlinear damping have significant effects on the chaotic behavior of the forced Duffing oscillator.



---

*Source: 3769870-2017-01-31.xml*