Mathematical Problems in Engineering

Volume 2017, Article ID 3769870, 8 pages

https://doi.org/10.1155/2017/3769870

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

Department of Construction Engineering, National Kaohsiung First University of Science and Technology, Kaohsiung, Taiwan

Correspondence should be addressed to Tai-Ping Chang; wt.ude.tsufkn.smcc@gnahcpt

Received 6 December 2016; Accepted 16 January 2017; Published 31 January 2017

Academic Editor: Jonathan N. Blakely

Copyright © 2017 Tai-Ping Chang. 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.

#### 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.

#### 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: where denotes the linear stiffness constant of the system, is linear damping coefficient, is the nonlinear damping coefficient, is the nonlinear stiffness constant, and and are the external amplitude and driving frequency, respectively. In the present study, , and are assumed to be small parameters. Hence, a transformation is adopted in order to apply the first-order perturbation scheme of the Melnikov method. Therefore, (1) can be rewritten asWhen , an unperturbed system can be obtained as follows:The system of (3) corresponds to a Hamiltonian system with a potential functionThe unperturbed system of (4) has three equilibrium points: one saddle and two centers .

#### 3. Chaotic Motion of Forced Duffing Oscillator

By integrating system (3), the solution of unperturbed system (3) can be written as follows: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 , if 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:where is the cross-section time of the Poincare map and 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:whereThe simple zeros of (7) give the critical value of the forcing amplitude (), 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), will vanish if a solution can be detected for , and ; therefore, the critical value of the forcing amplitude () can be obtained as follows: where 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:

In addition, the initial conditions in (1) are considered as follows:Based on (9), the variations of critical values of forcing amplitude versus linear damping coefficient , driving frequency , and nonlinear coefficient are depicted in Figures 1(a)–1(c). As can be seen from Figure 1(a), increases linearly with the linear damping coefficient; however, 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 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: , , , , . As can be seen from Figures 1(a)–1(c), the critical values of forcing amplitude 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).