Bifurcations and Chaos in the Duffing Equation with One Degenerate Saddle Point and Single External Forcing ()
1. Introduction
Since in 1918, the German electrical engineer Georg Duffing introduced the Duffing equation, many scientists have been widely studied the equation in physics, economics, engineering, and found many other physical phenomena. The Duffing oscillator, is normally written as
(1)
Depending on the parameters chosen, the equation can take a number of special forms. For example, Bender and Orszag [1] and Zwillinger [2] took the parameters
and studied the Duffing Equation (1) with no damping and no forcing,
(2)
Wiggins [3] took
and studied the following Duffing equation
(3)
Ravichandran et al. [4] replaced the external forcing
as various periodic external forcing. Equation (3) with one external forcing and two potential wells has many different types of oscillations such as chaos and limit cycles [5]. And Huang and Jing [6] studied the three well Duffing equation with one external forcing
(4)
and obtained the conditions of existence and bifurcations for harmonics, subharmonics and superharmonics under small perturbations and the threshold values of chaotic motion under periodic perturbation. And Jing et al. [7] [8] obtained complex dynamics of the three well Duffing equation with two external forcings,
(5)
Wang [9] presented analytical and numerical results concerning the inhibition of chaos in the Duffing equation with two weak forcing excitations. Jiang et al. [10], studied bifurcation and chaos of the three well Duffing equation with parametric excitation and one external forcing
(6)
But less attention was focused on the two well Duffing equation with one degenerated saddle. In this paper we studied the following Duffing equation
(7)
where
are real parameters. Physically,
can be regarded as dissipation or damping;
and
is the amplitude and frequency of the external force.
The structure of the paper is as follows. In Section 2, the fixed points and phase portraits are obtained for the unperturbed system of (7). In Section 3, the conditions of existence of chaos under periodic perturbation resulting from the homoclinic bifurcations are performed by Melnikov method. Finally, we make some numerical computations which give support to the theoretical analysis and some complex dynamics in Section 4.
2. Fixed Points and Phase Portrait of Unperturbed System of (7)
In this section, we obtain the stability of fixed points and phase portrait of unperturbed system of (7).
For we take
and obtain the unperturbed system of (7) as follows
(8)
The unperturbed system (7) can be easily obtained three fixed points: a degenerate saddle
and two centers
and
. The phase portrait of the unperturbed system (7) is plotted in Figure 1. The degenerate saddle is connected by two homoclinic orbits
and
, respectively.
In essence we use perturbation methods to study the system (7), we therefore study how the dynamics of unperturbed system (8) are changed under the periodic perturbation in the following parts.
3. Chaos for Periodic Perturbations
In this section, we consider the chaotic behaviors of system (7) in which
are assumed to be small parameters with order
. The Duffing system can be written an as follows:
(9)
where
.
The closed homoclinic orbits break when the perturbation is added, and system (7) may have transverse homoclinic orbits. By the Smale-Birkhoff Theorem [3] [11] [12], the existence of such orbits may results in chaotic dynamics. Therefore, we apply the Melnikov method to system (7) for finding the criteria of the existence of homoclinic bifurcation and chaos.
For the homoclinic orbit
, we have the Melnikov function,
(10)
where
is Bessel functions of the second. If we define
(11)
then it follows from Theorem 4.5.3 in [11] that if
, the stable manifold of the fixed point (0,0) intersects the unstable manifold for
sufficiently small, and if
, the stable manifold doesn’t intersect the unstable manifold. Moreover, since
has quadratic zeros when
, there is a bifurcation curve of system (9) in the
plane for each fixed
, tangent to
at
. This implies that if
is sufficiently small, the transverse heteroclinic orbits exist and system (9) may be chaotic. In Figure 2, we give the diagram of (11) in
plane for
.
For the homoclinic orbit
, the computation is identical and the similar result is obtained.
4. Numerical Simulations
In this section we give numerical simulations to look for other new dynamics. In the process of numerical simulation, we vary one parameter and fix the other parameters of system (7) as follows:
1) Varying
in the range
and fixing
and for rational and irrational values of
.
2) Varying
in the range
and fixing
and for rational and irrational values of
.
For case 1). The bifurcation diagram of system (7) in
plane and the corresponding Lyapunov exponents for
are given in Figure 3(a) and Figure 3(b), respectively. We observe that a period-1 window for
becomes chaos at
and chaos at 1.391 becomes a period-1 window for
. From Figure 3(c), the local amplification of Figure 3(a) for
, there exist three period-doubling bifurcation to chaos for
.
The bifurcation diagram of system (7) in
plane and the corresponding Lyapunov exponents for
are given in Figure 4(a) and Figure 4(b), respectively. A period-1 window suddenly becomes chaos at
and chaoticmotion suddenly becomes a period-doubling bifurcation at
. And for
period-doubling bifurcation to chaos and chaos to period-doubling bifurcations alternatively appear. Poincaré maps of chaotic attractors for
and
are shown in Figures 4(c)-(e), respectively.
Figure 2. The diagram of (11) in
plane.
(a) (b) (c)
Figure 3. (a) Bifurcation diagram of system (7) for
; (b) Lyapunov exponents corresponding to (a); (c) Local amplification of (a) for
.
For case 2). The bifurcation diagram of system (7) in
plane and the corresponding Lyapunov exponents for
are given in Figure 5(a) and (b). There exist onset of chaos at
from a period-1 window for
. And chaos suddenly becomes three inverse period-doubling bifurcation at
. We observe that chaos to inverse period-doubling bifurcations and inverse period-doubling bifurcations to chaos alternatively appear. And the size of chaotic attractors becomes smaller at
that an
interior crisis occurs. Poincaré maps of chaotic attractors for
and
are shown in Figure 5(c) and Figure 5(d), respectively.
The bifurcation diagram of system (7) in
plane and the corresponding Lyapunov exponents for
are given in Figure 6(a) and Figure 6(b). Period-2 orbit becomes period-1 orbit at
. From the local amplification Figure 6(c) of Figure 6(a), a period-1 window disappears and chaos appear at
. There exists an interior crisis of chaos at
and chaos regions becomes inverse period-doubling bifurcations to period-doubling
bifurcations for
. At
an intermittence of chaos occurs and chaotic motion becomes inverse period-doubling bifurcation. There is a bubble for
in the local amplification Figure 6(d) of Figure 6(a). Poincaré map for
is shown in Figure 6(e). Poincaré map of chaotic attractor for
is shown in Figure 6(f).
Acknowledgements
This work was supported by the National Science Foundations of China (10671063, 10801135 and 61571052).