Novel Hyperbolic Homoclinic Solutions of the Helmholtz-Duffing Oscillators

The exact and explicit homoclinic solution of the undamped Helmholtz-Duffing oscillator is derived by a presented hyperbolic function balance procedure. The homoclinic solution of the self-excited Helmholtz-Duffing oscillator can also be obtained by an extended hyperbolic perturbation method. The application of the present homoclinic solutions to the chaos prediction of the nonautonomous Helmholtz-Duffing oscillator is performed. Effectiveness and advantage of the present solutions are shown by comparisons.


Introduction
It has been widely accepted that homoclinic solutions play a fundamental role in global bifurcations and chaos predictions of dynamical systems [1,2].For instance, the experimental study of certain magnetic pendulum verified the homoclinic solutions as the precursors to chaotic vibration [3].Some occurrences of homoclinic solutions can be regarded as the criterion from single well chaos to cross well chaos motion of oscillators [4], or as the onsets of chaotic vibrations of asymmetric nonconservative oscillators [5].Homoclinic solutions were also adopted in bifurcation and chaotic vibration controls for beam structures [6,7].Another typical application of homoclinic solutions aims at solitary wave studies.For instance, a proper homoclinic solution can govern the solitary roll waves down an open inclined channel [8], or optical solitary waves propagating in fibers [9,10].The association between the singular solitary waves and homoclinic solutions can be interpreted based on phase plane analysis [11].
Because of their importance in nonlinear systems, many homoclinic solutions have been derived in the past few decades.Such works include but are not limited to the following: Xu et al. [12] proposed the perturbation-incremental method for homoclinic solutions; Chan et al. [13] applied the perturbation-incremental method to study the stability and the homoclinic bifurcations of limit cycles; Belhaq et al. [14] analytically developed criterions for predicting homoclinic connection of limit cycle.Mikhlin and Manucharyan [15] and Manucharyan and Mikhlin [16] applied the Padé and quasi-Padé approximants for homo-and heteroclinic solutions.Y. Y. Chen and S. H. Chen [17] and Chen et al. [18] developed perturbation techniques by hyperbolic functions for homoclinic solutions of strongly nonlinear oscillators.Cao et al. [19] improved the perturbation-incremental homoclinic solutions for strongly nonlinear oscillators.Recently, Li et al. [20] improved the perturbation method based on harmonic functions to derive homoclinic solutions of Helmholtz-Duffing oscillators.
Nevertheless, to the best of our knowledge, the completely analytical, exact, and explicit homoclinic solution of the strongly nonlinear Helmholtz-Duffing oscillators has not yet been derived, in spite of the wide application of its equation for many engineering problems such as ship dynamics, oscillation of the human ear drum, oscillations of onedimensional structural system with an initial curvature, some electrical circuits, microperforated panel absorber, and heavy symmetric gyroscope [21][22][23][24][25][26][27][28][29].It should be pointed out that the previous typical solutions [ [14][15][16][17][18] become invalid for such mix-parity systems.Even for the conservative Helmholtz-Duffing oscillator, solutions by the perturbation methods [12,13,19,20] based on generalized harmonic functions can only be obtained implicitly, in which the infinite time domain of a homoclinic motion has to be transformed into a finite period of the harmonic.Moreover, for strongly nonlinear oscillators, as the perturbation-incremental method [12,13,19] consists of perturbation procedure with the incremental harmonic balance method, their solutions are always expressed by harmonic functions with numerical coefficients.That means such implicit solutions are semianalytical and seminumerical and cumbersome for practical application.
This paper aims to present new homoclinic solutions of the Helmholtz-Duffing oscillators.The completely analytical, exact, and explicit homoclinic solution of the conservative Helmholtz-Duffing oscillator will be derived by a hyperbolic function balance procedure.Then the homoclinic solution of the self-excited Helmholtz-Duffing oscillator will also be obtained by an extended hyperbolic perturbation method.The application of the present solutions to the chaos prediction of the nonautonomous Helmholtz-Duffing oscillator is performed.The preference of the present solution will be illustrated by comparison.

The Explicit and Exact Homoclinic Solution of the Undamped Helmholtz-Duffing Oscillator
Consider the homoclinic solution of the undamped Helmholtz-Duffing equation If  2 = 0, (1) becomes the classical Duffing equation, which possesses a homoclinic solution with  1 < 0 and  3 > 0. Such homoclinic solution of classical Duffing equation has been discussed in detail in [17], in which the solution can be written as If  3 = 0, (1) becomes the classical Helmholtz equation, which possesses a homoclinic solution.Such homoclinic solution of classical Helmholtz equation has been discussed in detail in [18], in which the solution can be written as Noting the relationship as below, sech Here, to find a proper trial solution form for (1), we can observe the two special cases above.It can be seen that the two expressions above are similar, because they have the common form expressed as in which, when  3 = 0, the constants of ( 6) are While when  2 = 0,  1 < 0, and  3 > 0, the constants of ( 6) are  0 = ±√−2 1 / 3 ,  = 0,  0 = √ 1 , and  = 0. Thus, the time derivative of ( 6) is Note that (, 0) is the homoclinic point.For  2 ̸ = 0 and  3 ̸ = 0, we adopt (6) as a trial solution for the homoclinic solution of (1) and try to determine all its constants by substituting ( 6) into (1); that is, In order to balance (8) for all time , we equate coefficients of like powers of the hyperbolic function term ( + cosh  0 ) and get the following nonlinear algebraic equations: From ( 12), or The left hand side of ( 12) can be regarded as the restoring force of the oscillator with  = .In other words, (12) means that the displacement derivative of the potential energy curve at  =  is zero.Furthermore, we have to make sure that the potential energy curve at  =  is not concave.Thus, the displacement derivative of the restoring force at  =  will not be positive; that is, Therefore,  can be determined by ( 13)-( 15) and then, ( 9)-( 11) can be discussed, respectively, in the two cases as follows.
Example 1.Here we apply the method for equation which is a case of (1) with  1 = −1,  2 = −3, and  3 = 1.From ( 14), (21), and ( 22), we can determine all the constants and get the homoclinic solution as The time histories and the phase portraits of the solutions by different methods are shown in Figures 1 and 2, respectively.It can be seen from the figures that the present method yields accurate and explicit solutions in both the figures, while the generalized harmonic function perturbation method can only provide valid solution in Figure 2. The reason is that based on harmonic functions [12,13,19,20] the homoclinic solutions can only be expressed implicitly by the nonlinear time scale they adopted and be investigated only in phase planes.Such implicit solutions are too abstract or cumbersome to use in some practical problems.Therefore, the present explicit solutions in respect to time  are more applicable.

Perturbation Homoclinic Solution of the Self-Excited Helmholtz-Duffing Oscillator
Consider the homoclinic solution of the self-excited Helmholtz equation Shock and Vibration  where  denotes constants.We assume the homoclinic solution of (26) can still be expressed in the similar form of (6); that is, however, the amplitude  and the nonlinear time scale  will depend upon the perturbation parameter .Thus  and , respectively, can be expanded in the powers of ; that is, Then ( 27) can be rewritten as where Then After substituting ( 29) and ( 30) into ( 26), equating coefficients of like powers of  yields the following equations: Then solutions  0 ,  1 , . . .can be determined by solving linear equations ( 33), (34), . . .one by one.It can be seen that ( 33) is obtained from (1) via the transformation in (29).Therefore, the homoclinic solution of (33) can be given by (6).Multiplying (34) by   0 and integrating it from  0 to , we obtain where Noting the properties of hyperbolic functions, we have 0 (±∞) = , (37b) Thus letting  0 = −∞ and  = +∞ in (35), we derive Equation ( 38), which can also be derived by the classical Melnikov method, represents the critical condition under which the homoclinic bifurcation occurs.In other words, there exists a homoclinic solution once all the constants in ( 26) satisfy (38).Letting  0 = 0 and  = +∞ in (35) gives 1 ( 0 / ( + 1) + ) +  2 ( 0 / ( + 1) + ) 2 +  3 ( 0 / ( + 1) + ) 3 .
(39) Furthermore, substituting  0 = 0 into (35) yields The three equations above allow   ,  1 , and  1 to be determined one by one.As an illustration, here we consider in which  0 ,  1 , and  2 are constants.Noting (41), and substituting ( 6) and ( 7) into (36), the latter becomes where the expressions of the constants , ,   ,   , and   are listed in Appendix.Noting that, in (41),  0 is an even function and   0 is an odd function with respect to , (38) can be rewritten as Thus substituting (42) into (43), one derives Equation ( 44) is the condition under which homoclinic bifurcation occurs.Substituting (43) into (39) gives Then (40) can be rewritten as Finally, the expression for homoclinic solution can be expressed by It can be seen that once  2 or  3 becomes zero, the present procedures can be reduced to the methods and solutions presented in [17] or in [18].
where  cos Ω is the external harmonic excitation with the amplitude  and the frequency Ω.According to the Melnikov method [1,2], the Melnikov function of (52) can be written as where  0 is the solution presented in (6) and ∧ denotes the vector cross product.Then the chaotic response of the system may occur if there exists a simple zero point of ().Thus, we substitute ( 7) into (53), and let () = 0; the latter yields be pointed out that the Melnikov method is only regarded as one of the conditions for chaotic prediction.At present, a chaotic motion should still be evaluated more thoroughly with qualitative theory and numerical method.Below are two examples which satisfy the Melnikov condition but with different characters of chaotic motions.
Example 3. Consider the equation which is the case of (52) with  1 = −1,  2 = −1,  3 = 2,  = 0.1,  0 = −1, and Ω = 2. From ( 17), (55), and (56) we can derive that   = 1.43 when  0 = 0.94868 and   = 1.11 when  0 = −0.94868.Therefore it can be estimated that chaotic motion may happen if  > 1.43.The numerical results of homoclinic bifurcation by AUTO numerical method [30,31] show that, with  = 3, the Lyapunov exponent value shows chaos behavior from  = 0 to about  = 145 and then converts to 0 gradually.That means the chaotic motion possesses dissipative chaos property.The Lyapunov exponent diagram is shown in Figure 5.The Lyapunov exponent value converts to less than 0.01 after  = 1000.The phase portrait of the system after  = 1000 is shown in Figure 6, which shows that the motion converts to a limit cycle.17), (55), and (56) we can derive that   = 1.94 when  0 = 1.4138 and   = 1.82 when  0 = −1.4138.Therefore it can be estimated that chaotic motion may happen if  > 1.94.The numerical results of homoclinic bifurcation by AUTO numerical method [30,31] show that the Lyapunov exponent value stays more than 0 when  = 2.90.The Lyapunov exponent diagram is shown in Figure 7.The Poincaré projection of the system from  = 500 to  = 5000 is shown in Figure 8.In the figure, the fractal character of a strange attractor can be observed, which supports the prediction of chaotic motion.

Conclusions
The present procedures are efficient for constructing homoclinic solutions of the Helmholtz-Duffing oscillator.The exact and explicit homoclinic solution of the undamped Helmholtz-Duffing oscillator is derived by a hyperbolic function balance procedure.The homoclinic solution of the self-excited system is then obtained by the extension of the hyperbolic perturbation procedure.The application to the chaos prediction of the nonautonomous Helmholtz-Duffing oscillator can also be conducted.