# Chaotic Motions of the Duffing-Van der Pol Oscillator with External and Parametric Excitations

Authors: Liangqiang Zhou; Fangqi Chen
Journal: Shock and Vibration
 (2014)

---

## Abstract

The chaotic motions of the Duffing-Van der Pol oscillator with external and parametric excitations are investigated both analytically and numerically in this paper. The critical curves separating the chaotic and nonchaotic regions are obtained. The chaotic feature on the system parameters is discussed in detail. Some new dynamical phenomena including the controllable frequency are presented for this system. Numerical results are given, which verify the analytical ones.


---

## Body

## 1. Introduction

Dufffing-Van der Pol oscillator has a wide usage in many fields. The dynamics of Duffing-Van der Pol oscillator has been investigated widely in these years. Using the Melnikov method, Ravisankar et al. [1] studied horseshoe chaos in Duffing-Van der Pol oscillator driven by different periodic forces. Melnikov threshold curve was drawn in a parameter space. With the second-order averaging method and Melnikov method, Jing et al. [2] investigated chaotic motions in Duffing-Van der Pol equation with fifth nonlinear-restoring force and two external forcing terms. Numerical simulations were given to show the consistence with the theoretical analysis and exhibit the more new complex dynamical behaviors. With the singularity analysis, bifurcation properties of Duffing-Van der Pol system with two parameters under multifrequency excitations were studied by Qin and Chen [3]. Using the residue harmonic method, Leung et al. [4] investigated periodic bifurcation of Duffing-Van der Pol oscillators having fractional derivatives and time delay. It was shown that jumps and hysteresis phenomena can be delayed or removed. By using a simple transformation, the first integrals and the solutions of the Duffing-Van der Pol type equation under certain conditions were obtained by Udwadia and Cho [5]. Using the residue harmonic homotopy, a generalized Duffing-Van der Pol oscillator with nonlinear fractional order damping was investigated by Leung et al. [6]. Nonlinear dynamic behaviors of the harmonically forced oscillator were further explored by the harmonic balance method along with the polynomial homotopy continuation technique. By using the GYC partial region stability theory, Ge and Li [7] studied the synchronization of new Mathieu-Van der Pol systems with new Duffing-Van der Pol systems. With the second-order averaging method, the conditions for the existence and the bifurcations of harmonics for the damped and driven Duffing-Van der Pol system were obtained by Zhang et al. [8]. By the averaging method together with truncation of Taylor expansions, Li et al. [9] investigated the dynamics of Duffing-Van der Pol oscillators under linear-plus-nonlinear position feedback control with two time delays. By applying phase diagrams, potential diagram, Poincaré maps, bifurcation diagrams, and maximal Lyapunov exponent diagrams, the nonlinear behavior and the complex state of the Duffing-Van der Pol equation with fifth nonlinear-restoring force and two external periodic excitations were investigated by Shi et al. [10]. By using Melnikov analysis and numerical simulations, the dynamical behaviors including chaos, period-doubling cascades, and strange attractors of the extended Duffing-Van der Pol system were investigated by Yu et al. [11, 12]. With numerical methods, Yu et al. [13] also investigated the dynamical behavior of the extended Duffing-Van der Pol oscillator with ϕ
6 potential. Different routes to chaos and rich dynamical phenomena were observed. Patel and Sharma [14] revisited the stochastic Duffing-Van der Pol “filtering” in the Fokker-Planck setting in lieu of the filtering in the Kushner setting. With the help of the modified quasiconservative averaging, Li et al. [15] investigated the stochastic responses of Duffing-Van der Pol vibroimpact system under additive colored noise excitation.In this paper, the chaotic motions of the Duffing-Van der Pol oscillator with external and parametric excitations are studied analytically with the Melnikov method. The critical curves separating the chaotic and nonchaotic regions are obtained. The chaotic feature on the system parameters is discussed in detail and some new dynamical phenomena are presented. Numerical simulations verify the analytical results.


## 2. Formulation of the Problem

Consider the Duffing-Van der pol oscillator with external and parametric excitations(1)
x
¨
+
p
x
˙
(
1
-
x
2
)
-
α
x
+
β
x
3
=
F
(
t
)
,
where p is a damping parameter, α
>
0, β
>
0 are real parameters, and F
(
t
)
=
f
(
1
+
δ
x
)
cos
⁡
⁡
ω
t is the external and parametric forces.Assume the damping and excitation termsp, f are small, setting p
=
ε
p
-, f
=
ε
f
-, where ε is a small parameter; then (1) can be written as
(2)
x
˙
=
y
,
y
˙
=
α
x
-
β
x
3
+
ε
[
-
p
-
y
(
1
-
x
2
)
+
f
-
(
1
+
δ
x
)
cos
⁡
⁡
ω
t
]
.Using the transformations(3)
x
=
u
α
β
,
t
=
1
α
τ
,
then (2) can be written as
(4)
u
′
=
v
,
v
′
=
u
-
u
3
-
ε
p
~
v
(
1
-
α
β
u
2
)
+
ε
f
~
(
1
+
δ
u
α
β
)
cos
⁡
⁡
(
ω
~
τ
)
,
where p
~
=
p
-
/
α, f
~
=
f
-
/
(
α
α
/
β
), ω
~
=
ω
/
α, and ′ represents d
/
d
τ.Whenε
=
0, the unperturbed system of (4) is
(5)
u
′
=
v
,
v
′
=
u
-
u
3
which is a planar Hamiltonian system with the Hamiltonian
(6)
H
(
u
,
v
)
=
v
2
2
-
u
2
2
+
u
4
4
.System (6) has three equilibrium points, where (
0,0
) is a saddle point, and (
±
1,0
) are all centers.There exist homoclinic orbits connecting(
0,0
) to itself with the expressions [16]
(7)
u
hom
(
τ
)
=
±
2
sech
⁡
(
τ
)
,
v
hom
(
τ
)
=
∓
2
sech
⁡
(
τ
)
tanh
⁡
(
τ
)
and closed periodic orbits around (
±
1,0
) with the expressions [16]
(8)
u
k
(
τ
)
=
k
2
k
2
-
1
c
n
(
1
2
k
2
-
1
τ
,
k
)
,
v
k
(
τ
)
=
-
k
2
k
2
-
1
s
n
(
1
2
k
2
-
1
τ
,
k
)
d
n
(
1
2
k
2
-
1
τ
,
k
)
;
see Figure 1, where s
n
,
c
n
,
d
n are Jacobi elliptic functions, and 1
/
2
<
k
<
1 is the modulus of the Jacobi elliptic functions. The period of the closed orbit is T
k
=
4
2
k
2
-
1
K
(
k
), where K
(
k
) is the complete elliptic integral of the first kind.Figure 1
The phase portrait of system (5).


## 3. Chaotic Motions of the System

Melnikov method [17] is an analytical tool to study chaotic systems. Recently, chaotic motions of many systems, for example, ϕ
6-Rayleigh oscillator [18], Duffing oscillator [19], Gylden’s problem [20], and nonsmooth systems [21], have been investigated by the Melnikov method. In this section, we use the Melnikov method to investigate the chaotic motions of system (4). We compute the Melnikov functions of system (4) along the homoclinic orbit (7) as follows:
(9)
M
(
τ
0
)
=
∫
-
∞
+
∞
‍
-
p
~
v
hom
2
(
τ
)
(
1
-
α
β
u
hom
2
)
(
τ
)
d
τ
+
∫
-
∞
+
∞
‍
f
~
cos
⁡
⁡
ω
~
(
τ
+
τ
0
)
v
hom
(
τ
)
d
τ
+
δ
∫
-
∞
+
∞
‍
f
~
α
β
u
hom
(
τ
)
v
hom
(
τ
)
cos
⁡
⁡
ω
~
(
τ
+
τ
0
)
d
τ
≡
-
μ
~
I
0
+
g
~
(
I
1
+
δ
I
2
)
sin
⁡
ω
~
τ
0
,
where
(10)
I
0
=
∫
-
∞
+
∞
‍
v
hom
2
(
τ
)
(
1
-
α
β
u
hom
2
(
τ
)
)
d
τ
=
4
3
-
16
α
15
β
,
I
1
=
∫
-
∞
+
∞
‍
v
hom
(
τ
)
sin
⁡
ω
~
τ
d
τ
=
2
π
ω
~
sech
⁡
(
π
ω
~
2
)
=
2
π
ω
α
sech
⁡
(
π
ω
2
α
)
,
I
2
=
∫
-
∞
+
∞
‍
α
β
u
hom
(
τ
)
v
hom
(
τ
)
sin
⁡
ω
~
τ
d
τ
=
-
2
α
β
∫
-
∞
+
∞
‍
sech
2
⁡
(
τ
)
tanh
2
(
τ
)
sin
⁡
ω
~
τ
d
τ
=
2
(
α
/
β
)
e
-
π
ω
/
2
α
(
ω
2
-
α
)
π
(
e
-
π
ω
/
α
-
1
)
α
.
By Melnikov analysis, we estimate the condition for transverse intersection and chaotic separatrix motion as follows:
(11)
p
~
f
~
<
|
I
1
+
δ
I
2
|
|
I
0
|
;
that is,
(12)
p
-
f
-
<
|
I
1
+
δ
I
2
|
|
I
0
|
·
β
α
.First, takingδ
=
0, which is the case of periodic external excitation, letting α
=
1, for different values of β, we get the critical curves separating the chaotic regions (below) and nonchaotic regions (above) as in Figure 2(a). Next, Letting β
=
1, the critical curves for different values of α are shown in Figure 2(b).The critical curves for chaotic motions of system (2) in the case of δ
=
0.
(a)
(b)Secondly, takingδ
=
1, which is the case of both periodic external and parametric excitations, letting α
=
1, for different values of β, we get the critical curves as in Figure 3(a). Next, Letting β
=
1, the critical curves for different values of α are shown in Figure 3(b).The critical curves for chaotic motions of system (2) in the case of δ
=
1.
(a)
(b)From Figures2-3 we can obtain the following conclusions.(1)
For the case of periodic external excitation, the critical curves have the classical bell shape; this means that, with the excitations possessing sufficiently small or very large periods, the systems are not chaotically excited. Whenα is fixed, for each β
∈
(
0,0.8
), the larger the values of β, the larger the critical values for chaotic motions, while for β
>
0.8, the larger the values of β, the smaller the critical values for chaotic motions. On the other hand, when β is fixed, if α
<
0.9, for the case of small values of ω, that is, the period of the excitation is large, the critical value for chaotic motions decreases as α increases; when ω crosses a critical value, the case is opposite, so for the case of large values of ω, that is, the period of the excitation is small, the critical value for chaotic motions increases as α increases; if α
≥
0.9, for each excitation frequency ω, the critical value increases as α increases.
(2)
For the case of parametric excitations, the critical curve first decreases quickly to zero and then increases; at last it decreases to zero asω increases from zero. There exists a controllable frequency ω excited at which chaotic motions do not take place no matter how large the excitation amplitude is. When α (β) is fixed, the controllable frequency increases as β (α) increases.


## 4. Numerical Simulations

First, choosing the system parametersε
=
0.01, α
=
1, δ
=
0, β
=
0.5, and ω
∈
(
0,3
), which is the case of external excitations, the critical values p
-
/
f
- for chaotic motions are shown in Figure 4(a), where the red curve is the theoretical predictions, and the black points are the numerical predictions.The theoretical and numerical critical values for chaos of system (2) in the case of (a) δ
=
0 and (b) δ
=
1.
(a)
(b)Next, choosing the system parametersε
=
0.01, α
=
1, δ
=
1, β
=
0.5, and ω
∈
(
0,3
), which is the case of both external and parametric excitations, the critical values p
-
/
f
- for chaotic motions are shown in Figure 4(b), where the red curve is the theoretical predictions and the black points are the numerical predictions. From Figure 4, we can see that the difference of the critical values for chaotic motions between the theoretical and numerical predictions is very small, so numerical simulations agree with the analytical results.


## 5. Conclusions

Using the Melnikov and numerical methods, the chaotic motions for the Duffing-Van der Pol oscillator with external and parametrical excitations are investigated in this paper. The critical curves separating the chaotic and nonchaotic regions are obtained. It is shown that there exists a controllable frequencyω for the system with parametric excitations. When the system parameter α (β) is fixed, the controllable frequency increases as β (α) increases. Some new dynamical behaviors are presented.



---

*Source: 131637-2014-07-15.xml*