Complex dynamics in Duffing–Van der Pol equation

Duffing–Van der Pol equation with fifth nonlinear-restoring force and two external forcing terms is investigated. The threshold values of existence of chaotic motion are obtained under the periodic perturbation. By second-order averaging method and Melnikov method, we prove the criterion of existence of chaos in averaged system under quasi-periodic perturbation for ω₂ = nω₁ + εσ, n = 1, 3, 5, and cannot prove the criterion of existence of chaos in second-order averaged system under quasi-periodic perturbation for ω₂ = nω₁ + εσ, n = 2, 4, 6, 7, 8, 9, 10, where σ is not rational to ω₁, but can show the occurrence of chaos in original system by numerical simulation. Numerical simulations including heteroclinic and homoclinic bifurcation surfaces, bifurcation diagrams, Lyapunov exponent, phase portraits and Poincaré map, not only show the consistence with the theoretical analysis but also exhibit the more new complex dynamical behaviors. We show that cascades of interlocking period-doubling and reverse period-doubling bifurcations from period-2 to -4 and -6 orbits, interleaving occurrence of chaotic behaviors and quasi-periodic orbits, transient chaos with a great abundance of period windows, symmetry-breaking of periodic orbits in chaotic regions, onset of chaos which occurs more than one, chaos suddenly disappearing to period orbits, interior crisis, strange non-chaotic attractor, non-attracting chaotic set and nice chaotic attractors. Our results show many dynamical behaviors and some of them are strictly departure from the behaviors of Duffing–Van der Pol equation with a cubic nonlinear-restoring force and one external forcing.