Existence of Chaotic Phenomena in a Single-Well Duffing Oscillator under Parametric Excitations

In this paper, the existence of chaotic behavior in the single-well Duffing Oscillator was examined under parametric excitations using Melnikov method and Lyapunov exponents. The minimum and maximum values were obtained and the dynamical behaviors showed the intersections of manifold which was illustrated using the MATCAD software. This extends some results in the literature. Simulation results indicate that the single-well oscillator is sensitive to sinusoidal signals in high frequency cases and with high damping factor, the amplitude of the oscillator was reduced.


Introduction
Duffing oscillators have received remarkable attention in the recent decades due to the variety of their Engineering applications, for example magneto-elastic mechanical system [1], large amplitude oscillator of centrifugal governing system [2], nonlinear vibration beans and plates [3] [4] and fluid flow induced vibration [5].It is famous for the existence of chaos behavior in recent decades [6].
In 1979, the chaotic phenomena in Duffing equation had been investigated by Ueda [7].Chaos occurs when the behavior of the dynamical system is extremely sensitive to initial conditions.In mechanical system, it means a motion which trajectories starting from slightly different initial conditions diverge exponen- Thus, according to Smale-Birkhoff theorem in [1], it implies the existence of the chaotic behavior in the Smale-horseshoe sense.The Melnikov theory was firstly used to study chaos in Duffing system by Holmes [2] and generalized Melnikov function was developed by Wiggins [13] [14].This criterion is just the necessary condition for chaos but not sufficient for chaos, therefore, it must be sufficient conditions for the suppression of chaos [15].The Lyapunov exponent is an important indicator in determining the sensitivity of chaotic behavior which characterizes the average rate of the system in phase space between adjacent tracks of convergence and divergence [16].Whether the Lyapunov exponent is greater than zero or not is one of the most straight forward criterions to distinguish the chaotic systems [17].In other to calculate the Lyapunov exponent, some methods of solutions includes [16], nonlinear adaptive filtering method, QR matrix factorization method and its improvement methods.This paper makes use of two methods, the Melnikov method and improved QR matrix factorization.
The objectives of this paper therefore are to investigate the existence of chaotic behavior in a single-well Duffing oscillator forced by parametric excitations.
The rest of the paper is organized as follows: Section 2, explained the preliminaries to the results, Section 3 gives the main results using the Melnikov method and the calculation of the Lyapunov exponent and Section 4 presents the numerical simulations and finally some conclusions are given in Section 5.

Melnikov Method
One of the main tools for determining the existence or non-existence of chaos in a perturbed Hamiltonian system is Melnikov.In his theory, the distance between stable and unstable manifolds of the perturbed system were calculated up to the first order term.Melnikov method is a procedure which gives a bound on the parameters of a system such that chaos is predicted not to occur.The Melnikov method investigate the homoclinic bifurcation in the forced Duffing oscillator system with linear and non-damping.It measures the distance between stable and unstable manifolds in the Poincare section [6] and to preserve the homoclinic loops un- M t , that is the Melnikov function has a simple zero, then a homoclinic bifurcation occurs, implying that the chaotic motion occurs.

Melnikov Method for Predicting Chaos
Melnikov method gives an analytic tool for establishing the existence of transverse homoclinic points of the Poincare map for a periodic orbit of a perturbed dynamical system of the form; ( ) ( ) with n x R ∈ .It can also be used to establish the existence of sub-harmonic pe- riodic orbits of perturbed system of the form in (1).Furthermore, it can be used to show the existence of limit cycles and separatix cycles of perturbed planar system with 2 x R ∈ .For periodically perturbed planar systems, we have the form; ( ) ( ) where x R ∈ and g is periodic with period t in T. We assume that ( ) and we make the assumption; 1) For 0 ε = , the system (2) has a homoclinic orbit; ( ) 2) For 0 ε = , the system has a non-parametric family of periodic orbit.Then the Melnikov function ( ) The Melnikov method can be interpreted as a derivation in energy from the value on the perturbed separatix.Before stating main result established by Melnikov concerning the existence of transverse homoclinic point of the Poincare section, we need the following lemma and theory which establish the existence of a periodic orbit and hence the existence of the Poincare map with sufficient ε .Lemma 2.1 Under assumption 1) and 2), for ε sufficiently small, the system (2) has a unique hyperbolic periodic orbit; Correspondingly, the Poincare map P ε has a unique hyperbolic fixed point of saddle type; Theorem 2.1 Under the assumption 1) and 2), if the Melnikov function ( ) 0 M t has a simple zero in [0,1], then for all sufficiently small 0 ε ≠ , the stable and unstable manifold of the Poincare map P ε intersect transversally, that is, P ε has a transverse homoclinic point.
This theorem was established by Melnikov [1].The idea of the proof is that ( ) M t is a measure of the separation of the stable and unstable manifold of the Poincare map.The theory is an important result because it establishes the existence of transverse homoclinic point for P ε .It implies the existence of strange DOI: 10.4236/wjm.2019.9400570 World Journal of Mechanics invariant set for some iterate of P ε and the same type of chaotic dynamics for system (2) as for the Smale horseshoe map.Generally, the Melnikov method is very useful for detecting the presence of transverse homoclinic orbits and the occurance of homoclinic bifurcations.
Theorem 2.2 (Smale-Birkhoff Homoclinic Theorem) [18] Let f be a diffeomorphism ( C′ ) and suppose p is a hyperbolic fixed point.A homoclinic point is a point q p ≠ which is in the stable and unstable manifolds.If the stable and unstable manifolds intersect transversally at q, then q is called transverse.This implies that there is a homoclinic orbit ( ) n q q γ = such that lim lim . Since the stable and unstable manifolds are invariant, we have; ( ) ( ) for all n ∈  .Moreover, if q is transversal, so are all n q since f is diffeomorphism.

Method of Lyapunov Exponent
The method of Lyapunov exponent serves as a useful tool to qualify chaos.Specifically, Lyapunov exponent measures the rate of convergence or divergence of nearby trajectories [16] [18].Negative Lyapunov exponents indicate convergence while positive Lyapunov exponents demonstrate divergence and chaos.
The magnitude of the Lyapunov exponents is an indicator of the time scale on which chaotic behavior can be predicted or transients for the positive and negative cases respectively [19].
Physically, the Lyapunov exponent measures average exponential divergence or convergence between trajectories that differ only in having an infinitesimally small difference in their initial condition.The system is said to be chaotic if the trajectories remain within a bounded set of the dynamics.If one considers a ball of points in N-dimensional phase space in which each point follows its own trajectory based upon the system equations of motion over time, the ball of points will collapse to a simple point, will stay a ball or will become ellipsoid in shape [20].The measure of the rate at which this infinitesimal ball collapse or expands is the Lyapunov exponent.For a system written in the state-space form ( ) , small derivation from trajectory can be expressed by the equation . The maximal Lyapunov exponent is then defined by this equation.
Other useful quantities are the short time Lyapunov exponent and the local Lyapunov exponent.A short time Lyapunov exponent is simply a Lyapunov exponent defined over some finite time interval.The local Lyapunov exponent is a short time Lyapunov exponent in the limit where the time interval approaches zero.Both are dependent on starting points and the short time Lyapunov exponent is also independent on the magnitude of the time interval.If all points in the neighborhood of a trajectory converge towards the same orbit, the attractor is a fixed point or a limit cycle.However, if the attractor is strange, any two trajectories ( ) ( ) that starts over very close to each other separate exponentially with time.This sensitive initial condi-( ) where λ, the mean rate of separation of trajectories of the system is called the Lyapunov exponent, which can be estimated for long time t as; ( ) Equations ( 7) and ( 8) are for short and local Lyapunov exponent.The exponent can be positive or negative but at least one must be positive for an attractor to be classified as chaotic.In particular, if 0 λ < , the system converges to a stable fixed point or periodic orbits.A negative value of the Lyapunov exponent is characteristic of dissipative or non-conservative systems.If 0 λ = , the system is conservative and converges to a stable cycle limit.If 0 λ > , the system is unstable and chaotic.Hence, if the system is chaotic, it will have at least one positive Lyapunov exponent.Thus, the definition of chaotic system is based on a positive Lyapunov exponent.Finally, If λ = ∞ , the system is random.
Generally, the most used measure of sensitive initial condition is a system characterization by the Lyapunov exponent, which quantifies the rate of separation of infinitesimal close trajectories.For example, consider a one-dimensional system with two trajectories ( ) x t and ( ) x t which at some point 0 t are ar- bitrary close together and their difference in time tracked by the function; ( ) ( ) ( ) x t x t x t δ = − .The sign of the lyapunov exponent characterizes whether or not the system is exhibiting chaotic behavior.If the exponent is negative, the system, at least in that set of initial conditions is said to be stable (like trajectories go to like trajectories).A Lyapunov exponent of zero implies an unstable system which is essentially on the edge stable and chaotic.And of course a positive exponent implies the system is chaotic where trajectories exhibit exponential divergence.

The Single-Well Duffing Oscillator
The single-well Duffing equation under parametrical excitation is shown below; ( ) The System (9) has a unique hyperbolic limit cycle.Using the Melnikov theory, an analysis has been performed of the limit circles in oscillator systems described by single-well Duffing equation under perturbation.
Briefly, we describe Melnikov function and the bifurcations in perturbed Hamiltonian system as; where H the Hamiltonian ( ) , H H x y = is the analytic function.Also the perturbation functions ( ) , p x y and ( ) , q x y are analytic, ε is a small parameter.
Let ( ) ( ) ( ) , , x y x t y t ε ε = be the solution of (3.1).Then the solution of the unperturbed system at ( x t y t = . Further, we note that the unperturbed system at 0 ε = has one equilibrium point i.e. the center surrounded by a closed trajectories.

Melnikov Function for the Perturbed Single-Well Duffing Equation
In this work, the single-well Duffing equation is represented by; ( ) This equation can be rewritten in the following perturbed Hamiltonian system; ( ) where , 0 x t y t ε ε = be the solution of ( 11).The unperturbed system (11) has a Hamiltonian; and one equilibrium point surrounded by a closed trajectories.
The solution of the unperturbed system is expressed as; ( ) where sd is a Jacobian function.
Then, the Melnikov function for the System (3) is given as; now taking ( 13) into ( 14), we get; ( ) Then, using the following properties; The Melnikov function becomes; ( ) ( ) after a long calculation and introducing the notation ( ) We obtain the final expression as; ( ) ( ) where a b λ = .

The Lyapunov Exponent of a Single-Well Duffing Oscillator
Consider the Duffing equation below; ( ) where x  and x  are second-order and first-order derivative, n α ∈  , n β ∈  , δ is the damping, γ is the amplitude of the circle, ω is the angular frequency of the driven circle.In other to Type equation here, determine whether the system is in chaotic state, we need to calculate the Lyapunov exponent using the QR factorization method.
Let y x =  , ( ) which is written in matrix form as; ( ) ( ) According to the variational principle, its variational equations are; where ( ) Y t is a 2 by 2 matrix, I is a 2 by 2 unit matrix, ( ) J t is the Jacobian matrix of the system and its expression is; DOI: 10.4236/wjm.2019.9400574 World Journal of Mechanics Then, QR factorization of ( ) Y t can be written as; ( ) where Q is orthogonal matrix, R is upper triangular matrix.Substituting (21) in (20), we obtain the variational equation; ( ) ( )

Left multiply Equation (22) by T
Q and right multiply by 1 R − , we have; The orthogonal matrix Q is written as a function of angle variables.To the Duffing equation, its orthogonal matrix Q can be expressed by one angle θ .
The upper triangular matrix R can be expressed as; where θ is the angle variable, ( ) λ is the value associated with the Lyapunov exponent.Then, The correspondent matrix elements on both sides of (23) are equal, so we get; ( We add and subtract the first two differential equations and get a new differential equation.Together with the third differential equation, we obtain three new equations; Then from; We obtain; ( ) ( ) ( ) ( ) ( ) The time evolution of the Lyapunov exponent is; Then, the Lyapunov exponent is; ( ) ( )

Numerical Simulation of Single-Well Duffing Oscillator
In this section, we compare the numerical solution of Equation ( 9    represent our chaotic solutions.1.The values were generated using the vector initial function values and the constant.However, the solutions are in .The solid lines represent the numerical solution and the dash lines represent the chaotic behavior., : cos Define an additional argument for the ODE solver: 0 : 0 t = Initial value of independent variable

Conclusions
In the present study, the chaotic behavior in single-well Duffing oscillator is investigated using Melnikov approach and Lyapunov exponent.The distance between the stable and unstable manifold of the nonlinear system is calculated by Melnikov approach.The Lyapunov exponent of the nonlinear system is evaluated by QR factorization to determine whether the chaotic phenomenon of the nonlinear system actually occurs.
As a result, threshold values were obtained and the dynamical behaviors showing the intersections of manifold were illustrated.To detect the chaotic phenomena of the nonlinear system, the Melnikov approach, Lyapunov exponent, the time history, phase portrait of the nonlinear system were presented for various cases.
How to cite this paper: Eze, E.O., Goodluck, O.C., Ngozi, U.R. and Oko, T.O.(2019) Existence of Chaotic Phenomena in a Single-Well Duffing Oscillator under ) using MATCAD simulation.In Figures 1-6, the trajectory versus time response curves are plotted for different sets of parameter values noted in the figure captions.In all figures, the solid lines represent the numerical solution and the dashed lines

Figure 3 .
Figure 3.The phase portrait orbits in the chaotic state at 0.5 k = .

Figure 1 and
Figure 1 and Figure 2 compare solutions by considering a strong nonlinearity value of

Figure 6 .
Figure 6.The phase portrait orbits in the chaotic state at 2 k = .
Numbers of solution values on [t 0 , t 1 ] E. O. Eze et al.

Table 1 .
Solution matrix table for solution functions.