Shock and Vibration

Volume 2018 (2018), Article ID 7213606, 12 pages

https://doi.org/10.1155/2018/7213606

## Heteroclinic Bifurcation Behaviors of a Duffing Oscillator with Delayed Feedback

^{1}Department of Traffic and Transportation, Shijiazhuang Tiedao University, Shijiazhuang 050043, China^{2}Department of Mathematics and Physics, Shijiazhuang Tiedao University, Shijiazhuang 050043, China^{3}Department of Engineering Mechanics, Shijiazhuang Tiedao University, Shijiazhuang 050043, China

Correspondence should be addressed to Shao-Fang Wen

Received 17 July 2017; Revised 22 November 2017; Accepted 21 December 2017; Published 21 January 2018

Academic Editor: Gianluca Gatti

Copyright © 2018 Shao-Fang Wen 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.

#### Abstract

The heteroclinic bifurcation and chaos of a Duffing oscillator with forcing excitation under both delayed displacement feedback and delayed velocity feedback are studied by Melnikov method. The Melnikov function is analytically established to detect the necessary conditions for generating chaos. Through the analysis of the analytical necessary conditions, we find that the influences of the delayed displacement feedback and delayed velocity feedback are separable. Then the influences of the displacement and velocity feedback parameters on heteroclinic bifurcation and threshold value of chaotic motion are investigated individually. In order to verify the correctness of the analytical conditions, the Duffing oscillator is also investigated by numerical iterative method. The bifurcation curves and the largest Lyapunov exponents are provided and compared. From the analysis of the numerical simulation results, it could be found that two types of period-doubling bifurcations occur in the Duffing oscillator, so that there are two paths leading to the chaos in this oscillator. The typical dynamical responses, including time histories, phase portraits, and Poincare maps, are all carried out to verify the conclusions. The results reveal some new phenomena, which is useful to design or control this kind of system.

#### 1. Introduction

It is well known that many kinds of nonlinearities exist in engineering systems, such as parametric excitation, nonsmoothness, time delay, discontinuity, and large deformation [1–6]. These nonlinearities will make the system responses more complicated and the system performance deteriorated. In recent years, the complicated dynamics of typical nonlinear systems have been more and more concerned. Duffing oscillator is one of the most common and typical models in nonlinear dynamical systems. Among these nonlinear systems the well-known Duffing oscillator is quite suitable to model the large deformation structure in many physics and engineering fields [7–10]. Moreover, time delay is an unavoidable and very common problem when those systems were controlled. It can influence the system dynamical characteristic and even destroy the system stability [11, 12]. The complicated Duffing nonlinear system with time delay is possible to generate more complex bifurcation and chaotic dynamic phenomena [11–14]. For example, Luo [15] studied the bifurcation trees of time-delay Duffing oscillator by a semianalytical method. Amer et al. [16] investigated the Duffing oscillator with parametric excitation under time-delay feedback based on the multiple scales perturbation method and analyzed the influences of the system parameters. Ji and Leung [17] studied the bifurcation control of a Duffing oscillator with parametrical excitation by a linear time-delayed feedback and discovered that the stable region could be broadened through choosing an appropriate feedback control. Zhang et al. [18] investigated the multipulse global bifurcations and chaotic behaviors of a cantilever beam by extended Melnikov method. Theodossiades and Natsiavas [19] and Van Dooren [20] studied the period-doubling bifurcations for gear-pair systems with periodic stiffness and backlash. Rusinek et al. [21, 22] investigated the dynamics of the time-delay Duffing oscillator. Besides, they also studied the chaos and feedback control of the time-delay Duffing system and found the suitable time delay and feedback gain would destroy the chaotic attractor for Duffing system. Nana Nbendjo et al. [23] studied the effects of the control parameters in a double-well Duffing oscillator with the time delay by Melnikov theory.

The heteroclinic bifurcation and chaos behaviors are two of the most important characteristics in nonlinear dynamical systems. Many scholars paid attention to the heteroclinic bifurcation and chaos phenomena of nonlinear Duffing oscillator [24–26]. For example, Yang and Sun et al. [27, 28] investigated the necessary condition for the generating chaos of a double-well Duffing oscillator with bounded noise excitations and time-delay feedback by Melnikov theory. Siewe et al. [29] investigated the necessary condition for chaotic behavior of Duffing-Rayleigh system based on Melnikov theory. Cao et al. [30] presented a novel construction of homoclinic and heteroclinic orbits for nonlinear systems with a perturbation-incremental method. Chen and Yan [31] studied the heteroclinic bifurcation behavior in the Duffing-VDP oscillator by the hyperbolic Lindstedt-Poincare method and obtained the analytical heteroclinic solution. Chacon [32] and Maki et al. [33] investigated heteroclinic bifurcation phenomenon for different types of nonlinear systems based on Melnikov method, respectively. Lei and Zhang [34] investigated the chaotic motion of the Duffing system and came to the conclusion that the threshold for generating chaos could be changed by choosing the internal parameters of trichotomous noise.

All the above analyses mainly focused on qualitative, numerical, or simplified analytical solutions about the necessary condition for generating chaos. In this paper, the first-order exact analytical solution of the necessary condition for generating chaos in sense of Smale horseshoes in a Duffing oscillator with both delayed displacement feedback and delayed velocity feedback is obtained based on Melnikov theory. Besides, the numerical heteroclinic bifurcation results were presented and some new phenomena were found in the Duffing oscillator with time-delay feedback. The basic structure of this paper is arranged as follows. The Melnikov function is obtained based on Melnikov method, and the analytical necessary condition for generating chaos is also obtained in Section 2. In Section 3, the bifurcation curves and the largest Lyapunov exponents by numerical method are investigated. It can be found that there are two paths leading to chaos via period-doubling bifurcation in this system. Then the time histories, phase portraits, and Poincare maps with the typical system parameters are all presented to verify the new phenomenon. The influences of the delayed displacement and velocity feedback parameters on generating chaos are analyzed. The comparisons of the numerical results with the analytical results are also fulfilled in this section. Finally, the conclusions of this paper are summarized.

#### 2. Analytical Necessary Condition for Chaos in Sense of Smale Horseshoes

Duffing oscillator is one of the most familiar systems in nonlinear dynamics. Under the function of both delayed displacement feedback and delayed velocity feedback, a Duffing oscillator with forcing excitation would be investigated in this section. The dynamic equation is where is the linear stiffness coefficient, is nonlinear stiffness coefficient, is linear damping coefficient, and are excitation amplitude and frequency, respectively, and are displacement and velocity feedback coefficients, respectively, and and are time delays of displacement and velocity feedback, respectively. Here all the system parameters are positive and dimensionless.

Introducing the transformationwhere is a small real parameter, (1) turns intoSupposing , the unperturbed system is There are three equilibrium points, where is a center, and are two saddle points. Generally speaking, if the unperturbed system is conservative system and the number of saddle points is 1 larger than that of the centers, there may be heteroclinic orbit in this system. Here the heteroclinic orbit connecting the two saddle points satisfies the formulaSupposing at , the calculating result isIntegrating (5), one could obtainCalculating (7), it yieldsSo (3) can be rewritten aswhereThen one can establish the heteroclinic orbit as follows:Melnikov theory [35], 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 technique was firstly applied to study a periodically driven Duffing oscillator by Holmes [36]. The Melnikov technique is also valid in time-delay systems, which was proved in [37]. According to Melnikov theory, the perturbed function must be a periodicity of function. In (9a), the perturbed function with time delay is also a periodic function which is satisfied with the basic condition of Melnikov theory. Therefore, many scholars have studied the dynamical characteristics of the time-delay Duffing system by Melnikov theory such as [23, 27, 28, 37]. All the existing analyses were mainly focused on qualitative, numerical, or simplified analytical solutions about the necessary condition for generating chaos. Here, Melnikov method is also applied to yield the Melnikov function as follows:whereMaking use of the odevity of integrand, becomesAnd the other parts in (11a) are Through the analysis of the above results, the necessary condition for generating chaos in sense of Smale horseshoes can be obtained as follows:Replacing the parameters in (12) with the original system parameters, one can get

#### 3. Numerical Simulation and Influences of Delayed Feedback on Bifurcation and Chaotic Behaviors

If the delayed displacement feedback coefficient and velocity feedback coefficient are all small enough, then (13) becomesThrough the analysis of (14), we could find that the analytical necessary condition for generating chaos is influenced by the delayed displacement feedback and delayed velocity feedback, respectively. That is to say, the coupling relationship does not exist between delayed displacement feedback and delayed velocity feedback. According to the above analysis of the relationship between the two kinds of feedback, the influences of the two kinds of delayed feedback on bifurcation and chaos behaviors will be investigated separately in the following sections.

##### 3.1. The Influences of Delayed Displacement Feedback Parameters

If there is only the delayed displacement feedback in (1), (14) becomes

###### 3.1.1. The Influence of Feedback Coefficient

By the analysis of (15), we can findso thatIt could be found that the right-hand parts of (15) would become larger with the increasing of the delayed displacement feedback coefficient . Accordingly, the critical excitation amplitude for generating chaos will become larger with the increasing of . In other words, the threshold value for chaotic motion in (1) will become larger, so that it is harder to generate chaos with larger in (1).

In order to verify the validity of the analytical necessary condition for generating chaos, a set of illustrative system parameters is chosen as , , , , and . The numerical solutions for of (1) are also investigated by numerical iterative method. When , the bifurcation diagram and the corresponding largest Lyapunov exponents are shown in Figures 1(a) and 1(b), respectively. The numerical simulation method for the largest Lyapunov exponents about delayed differential equation can be found in [38]. From the analysis of Figures 1(a) and 1(b), it could be found that the values of generating chaos in two figures are consistent, which implies that the results by numerical iterative method are correct. Through the numerical simulation, the authors also observe some special dynamical phenomenon in this complex system. We find that there exist two paths leading to chaos via different period-doubling bifurcation in (1) which are shown in Figures 2(a) and 3(a), respectively. If the two paths are drawn on the same diagram which is shown in Figure 1, it could be found that the intersection point of the two paths is . When , there is a single periodic solution. If , there will be two periodic-1 solutions, which depend on the initial values. The path shown in Figure 2(a) is simply named as “Type 1” and the other path shown in Figure 3(a) is “Type 2.” In order to verify the conclusion, the time history, phase portrait, and Poincare maps at different typical points are presented and analyzed in the following part.