Research Article  Open Access
A. Beléndez, M. L. Alvarez, J. Francés, S. Bleda, T. Beléndez, A. Nájera, E. Arribas, "Analytical Approximate Solutions for the CubicQuintic Duffing Oscillator in Terms of Elementary Functions", Journal of Applied Mathematics, vol. 2012, Article ID 286290, 16 pages, 2012. https://doi.org/10.1155/2012/286290
Analytical Approximate Solutions for the CubicQuintic Duffing Oscillator in Terms of Elementary Functions
Abstract
Accurate approximate closedform solutions for the cubicquintic Duffing oscillator are obtained in terms of elementary functions. To do this, we use the previous results obtained using a cubication method in which the restoring force is expanded in Chebyshev polynomials and the original nonlinear differential equation is approximated by a cubic Duffing equation. Explicit approximate solutions are then expressed as a function of the complete elliptic integral of the first kind and the Jacobi elliptic function cn. Then we obtain other approximate expressions for these solutions, which are expressed in terms of elementary functions. To do this, the relationship between the complete elliptic integral of the first kind and the arithmeticgeometric mean is used and the rational harmonic balance method is applied to obtain the periodic solution of the original nonlinear oscillator.
1. Introduction
The cubicquintic Duffing equation is a differential equation with third and/or fifthpower nonlinearity [1]. Due to the presence of this fifthpower nonlinearity added to the third nonlinearity of the common Duffing equation, this oscillator is difficult to handle. As Lai et al. pointed out [1], the cubicquintic Duffing equation can be found in the modelling of different systems and we advise readers to consult this paper to learn more about these applications and a good set of references in which this strongly nonlinear equation can be found. Among these systems modelled by the cubicquintic Duffing equation we would like to mention the nonlinear dynamics of a slender elastica, the compound KdV equation in nonlinear wave systems, or the propagation of a short electromagnetic pulse in a nonlinear medium [1–3].
The purpose of this paper is to present analytical approximations to the periodic solutions to the cubicquintic Duffing oscillator [1–3] in terms of elementary functions. To do this, we use the results previously obtained [4] using the Chebyshev series expansion of the restoring force [5, 6] for a quintic Duffing oscillator. Using this approximation, the original nonlinear differential equation is replaced by a cubic Duffing equation, which can be exactly solved. The replacement of the original nonlinear equation by an “approximate cubic Duffing equation” allows us to obtain an approximate frequencyamplitude relation as a function of the complete elliptic integral of the first kind and the solution in terms of the Jacobi elliptic function cn [4]. Using these first approximate equations, closedform expressions for the approximate frequency and the solution in terms of elementary functions are obtained using the arithmeticgeometric mean combined with the (secondorder) rational harmonic balance method.
In relation to cubication methods, BravoYuste [7] introduced a cubication procedure based on the harmonic balance method and BravoYuste and MartínSánchez [8] also developed a weighted mean square method of cubication. Recently, Beléndez et al. [4, 9, 10] and ElíasZùñiga et al. [11] employed similar procedures based on the use of Chebyshev polynomials. In particular, ElíasZùñiga et al. combined the cubication method with the equivalent nonlinearization method developed by Cai et al. [12].
In the early 1980’s and 1990’s, much effort was put on renormalization techniques and the use of the wellknown arithmeticgeometric mean to analyze nonlinear oscillators and, in particular, the old famous simple pendulum [13]. Recently, and also analyzing the nonlinear pendulum, it has been shown that it is possible to obtain an approximate expression for the complete elliptic integral of the first kind by using a Taylor series expansion [14], which is similar to the second approximation derived using the arithmeticgeometric mean. Finally, the rational harmonic balance method (RHBM) [15] is more difficult to apply than the common harmonic balance method (HBM) for oscillators with complex nonlinearities. However, we will see that the RHBM is easy to apply to the cubic Duffing equation obtained when the Chebyshev polynomials expansion is taken into account. The RHBM presents the advantage versus the usual HBM that the RHBM includes approximations for all the harmonics contributing to the periodic solution [16].
2. Solution Procedure Using a Chebyshev Polynomials Expansion
A cubicquintic Duffing oscillator is a conservative autonomous oscillatory system, which can be described by the following secondorder differential equation with cubicquintic nonlinearities [1]: with initial conditions where is an odd function; and are dimensionless displacement and time variables, respectively; , , and are positive constant parameters [1]. The angular frequency of these oscillations is denoted by , which is a function of the initial amplitude .
The cubicquintic Duffing equation is difficult to handle because of the presence of strong nonlinearity and it has no known closedform solution. However, the quintic Duffing equation (, , and in (2.1)) has been recently solved approximately by means of a cubication procedure based on the expansion of the nonlinear restoring force in terms of the Chebyshev polynomials [4]. To apply this procedure, firstly a reduced variable is introduced in (2.1) and (2.2) and the nonlinear function is expanded in terms of Chebyshev polynomials of the first kind . Retaining only the first two terms in this expansion, we obtain where Then the nonlinear differential equation in (2.1) can be then approximated by the (cubic) Duffing equation where
Equation (2.5) can be exactly solved and its solution can be written in terms of the Jacobian elliptic function cn [15]. Then the approximate frequency and solution for the cubicquintic equation in (2.1) will be the exact frequency and solution for the cubic equation in (2.5). Taking into account the wellknown expressions for the frequency and the solution for the Duffing oscillator in (2.5), we can easily obtain the following analytical approximate expressions for the frequency and solution for the cubicquintic Duffing oscillator as follows: where is the complete elliptic integral of the first kind
3. ClosedForm Expression for the Approximate Frequency in Terms of Elementary Functions Using the ArithmeticGeometric Mean
As we can see from (2.7), the approximate frequency for the cubicquintic Duffing oscillator is expressed in terms of the complete elliptic integral of the first kind, . However, it is also possible to obtain a closedform expression for the approximate frequency in terms of elementary functions. The complete elliptic integral of the first kind (2.9) cannot be expressed in terms of elementary functions, but can be numerically evaluated with high precision by a simple procedure based on the arithmeticgeometric mean. This is due to the fact that the arithmeticgeometric mean is the basis of Gauss’ method for the calculation of elliptic integrals [16–18]. The arithmeticgeometric mean of two numbers and , can be expressed in closedform in terms of the complete elliptic integral of the first kind as [16] which allows us to write in (2.9) as follows [16]:
The Legendre form of the arithmeticgeometric mean is given by [16] where and . From (3.3) it is possible to approximately calculate the arithmeticgeometric mean by choosing a finite number of terms in this equation as follows: From (2.7), (3.2), and (3.4) we obtain the following expression for the approximate frequency of the cubicquintic Duffing oscillator: where only two terms () have been considered in (3.4) where and (2.7). The maximum value of reached for the cubicquintic Duffing oscillator is because for the case in which and we obtain . We can see that (3.5) is a very simple and easily computable expression.
From (3.2) and (3.6) it can be seen that it is possible to take into account the following approximation for the complete elliptic integral of the first kind: In Table 1 we present the percentage errors for comparison of the approximate expression for the complete elliptic integral of the first kind, (3.7), with the exact one, , for the interval of values of corresponding to the cubicquintic Duffing oscillator (). As we can see from this table, is a very good approximation for and the relative error is less than 0.087% for .

To verify if it is sufficient to consider only a twoterm approximation in the Legendre form of the arithmeticgeometric mean, we compare (2.7) and (3.5) for and for and .
For we obtain the following power series expansions for these equations for very large values of the oscillation amplitude
Furthermore, we have the following: The maximum relative error for the frequency using the arithmeticgeometric mean approximation in (3.5) instead of (2.7) is 0.0087%, which allows us to conclude that (3.6) is a suitable approximation to (3.3).
For and we obtain the following power series expansions for these equations for very large values of the oscillation amplitude : where is the Euler gamma function. Furthermore, we have the following equations: Now the maximum relative error for the frequency using the arithmeticgeometric mean approximation in (3.5) instead of (2.7) is as low as 0.0014%, which allows us to conclude that (3.6) is a suitable approximation to (3.3).
4. ClosedForm Expression for the Solution in Terms of Elementary Functions Using the Rational Harmonic Balance Method
In order to obtain an expression for the approximate solution in terms of elementary functions we consider that the approximate frequency is given by (3.5) and we apply the rational harmonic balance method (RHBM). By defining a new independent variable replacing the time variable, , (2.5) can be expressed as follows: where we have written for simplicity. We propose the following approximation to the periodic solution of (4.1) Note that this expression satisfies the prescribed initial conditions independently of parameter , which is to be determined. It is easy to verify that substituting in (4.2), we obtain the rational representation which is the expression usually considered in the bibliography for the secondorder rational harmonic balance method [15]. However, (4.2) has the advantage that it allows us to simplify the procedure to obtain the coefficients of the Fourier series expansion of and . These series can be written as follows: where Substituting (4.4) into (4.1) the following result is obtained: where HOH stands for higherorder harmonics. Setting the coefficient of the lowest harmonic to zero gives which can be solved to obtain as a function of . Before doing this, it is necessary to obtain the values of the first Fourier coefficients and as a function of the unknown parameter . To do this we firstly write the expression for as follows and then Then Fourier coefficients and take the expressions Substitution of (4.10) into (4.7), yields the following expression Solving for gives where , and are given in (2.6) and (3.5). From (4.12) we obtain for and for and In Figure 1 we have plotted parameter as a function of for different values of parameters , , and . Equations (4.13) and (4.14) allow us to know the maximum value of parameter for and (), respectively. In fact, for and we have obtained an accurate approximate solution for the cubic Duffing equation.
5. Comparison with the Exact Solution
The exact frequency for the cubicquintic Duffing oscillator can be obtained by direct integration of (2.1) using the initial conditions in (2.2). The result is (see Appendix B in [1]) where parameters , , and are given as follows: For it is possible to obtain the following expressions for very small values of : For it is possible to obtain the following expression for very large values of : Based on (3.9) and (5.4), the approximation compared with respect to the exact frequency is As can be seen, the relative error of the analytical approximate frequency is 0.38% in the limit as , which implies that the method proposed is suitable for solving (2.1) with much larger than unity. In Figure 2 we have plotted the comparison between the approximate and exact frequencies for different values of parameters , , and .
For this cubicquintic Duffing oscillator, Lai et al. obtained the following limits applying the Newtonharmonic balance method (NHBM) [1]: for the first, second, and thirdorder analytical approximations, respectively. As we can see, the approximate frequency obtained in this paper, combining a cubication procedure (based on Chebyshev polynomials expansion of the restoring force), the arithmetic geometric mean (for approximating the exact frequency of the cubic Duffing oscillator), and the secondorder rational harmonic balance method, is more accurate than the firstorder and secondorder frequencies obtained using the NHBM and slightly less accurate (0.38% versus 0.23%) than the thirdorder frequency obtained using the NHBM.
In Table 2 we have presented the percentage errors for comparison of the approximate frequency with exact frequency for . For these parameters (3.5) can be written as follows: Taking into account (3.9), (4.2), and (4.13) it is easy to verify that the analytical approximate solution for the cubicquintic Duffing oscillator has the following asymptotic expression for : We analyze this example because this situation corresponds to the worst case, since it has the approximate frequency whose relative error is the greatest.

Approximate and exact solutions for a few examples of the cubicquintic Duffing oscillators are presented in Figures 3, 4, 5, and 6 for one complete cycle. The exact periodic solutions have been obtained by numerically integrating (2.1), and the proposed approximate periodic solutions have been computed using (4.2). In these figures parameter is defined as which tends to as . The approximate solution, expressed as a function of , which has been plotted in Figure 6 is Parameters for oscillators analyzed in Figures 3–6 are presented in Table 3 in which the values for the norm errors have been also included in order to produce a global estimation of the accuracy of the approximate solutions [18] The maximum value of norm is 0.0100931 and it is obtained for . In Figure 7 we have plotted the differences for the oscillators presented in Figures 3–6. Figures 3–7 show that (3.9) and (4.2) provide a good approximation to the exact periodic solution to the cubicquintic Duffing oscillator.

6. Conclusions
Using the results obtained by applying a cubication method for the cubicquintic Duffing oscillator, very simple and accurate expression for the frequency and the solution are obtained for the cubicquintic Duffing oscillator. In this procedure the original nonlinear equation is replaced by a related “cubic Duffing equation” that approximates the original “cubicquintic Duffing equation” closely enough to provide useful solutions. The cubic Duffing equation is exactly solved in terms of the complete elliptic integral of the first kind and the Jacobi elliptic function cn and these exact solutions are the approximate solutions to the cubicquintic Duffing equation. Using the relationship between the complete elliptical integral of the first kind and the arithmeticgeometric mean and the secondorder rational harmonic balance method, closedform expressions for the approximate frequency and the solution are obtained in terms of elementary functions with a relative error for the approximate frequency less than 0.38% for all values of the oscillation amplitude. The results of the method presented in this paper have been found to be in good agreement with those of the exact solutions obtained numerically.
Acknowledgments
This work was supported by the “Generalitat Valenciana” of Spain, under Project PROMETEO/2011/021, and by the “Vicerrectorado de Tecnología e Innovación Educativa” of the University of Alicante, Spain, under Project GITE09006UA.
References
 S. K. Lai, C. W. Lim, B. S. Wu, C. Wang, Q. C. Zeng, and X. F. He, “Newtonharmonic balancing approach for accurate solutions to nonlinear cubicquintic Duffing oscillators,” Applied Mathematical Modelling. Simulation and Computation for Engineering and Environmental Systems, vol. 33, no. 2, pp. 852–866, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J. I. Ramos, “On LinstedtPoincaré technique for the quintic Duffing equation,” Applied Mathematics and Computation, vol. 193, no. 2, pp. 303–310, 2007. View at: Publisher Site  Google Scholar
 C. W. Lim, R. Xu, S. K. Lai, Y. M. Yu, and Q. Yang, “Nonlinear free vibration of an elasticallyrestrained beam with a point mass via the Newtonharmonic balancing approach,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 10, no. 5, pp. 661–674, 2009. View at: Publisher Site  Google Scholar
 A. Beléndez, G. Bernabeu, J. Francés, D. I. Méndez, and S. Marini, “An accurate closedform approximate solution for the quintic Duffing oscillator equation,” Mathematical and Computer Modelling, vol. 52, no. 34, pp. 637–641, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J. H. Denman, “An approximate equivalent linearization technique for nonlinear oscillations,” Journal of Applied Mechanics, vol. 36, no. 2, pp. 358–360, 1969. View at: Publisher Site  Google Scholar
 R. E. Jonckheere, “Determination of the period of nonlinear oscillations by means of Chebyshev polynomials,” Zeitschrift fur Angewandte Mathematik und Mechanik, vol. 51, pp. 389–393, 1971. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 S. BravoYuste, “Cubication of nonlinear oscillators using the principle of harmonic balance,” International Journal of NonLinear Mechanics, vol. 27, no. 3, pp. 347–356, 1992. View at: Publisher Site  Google Scholar
 S. BravoYuste and A. MartínSánchez, “A weighted meansquare method of “cubication” for nonlinear oscillators,” Journal of Sound and Vibration, vol. 134, no. 3, pp. 423–433, 1989. View at: Publisher Site  Google Scholar
 A. Beléndez, M. L. Álvarez, E. Fernández, and I. Pascual, “Cubication of conservative nonlinear oscillators,” European Journal of Physics, vol. 30, no. 5, pp. 973–981, 2009. View at: Publisher Site  Google Scholar
 A. Beléndez, D. I. Méndez, E. Fernández, S. Marini, and I. Pascual, “An explicit approximate solution to the Duffingharmonic oscillator by a cubication method,” Physics Letters A, vol. 373, pp. 2805–2809, 2009. View at: Publisher Site  Google Scholar
 A. ElíasZúñiga, O. MartínezRomero, and R. CórdobaDíaz, “Approximate solution for the Duffingharmonic oscillator by the enhanced cubication method,” Mathematical Problems in Engineering, vol. 2012, Article ID 618750, 12 pages, 2012. View at: Publisher Site  Google Scholar
 J. Cai, X. Wu, and Y. P. Li, “An equivalent nonlinearization method for strongly nonlinear oscillations,” Mechanics Research Communications, vol. 32, no. 5, pp. 553–560, 2005. View at: Publisher Site  Google Scholar
 S. Bullett and J. Stark, “Renormalizing the simple pendulum,” SIAM Review, vol. 35, no. 4, pp. 631–640, 1993. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A. Beléndez, E. Arribas, A. Márquez, M. Ortuño, and S. Gallego, “Approximate expressions for the period of a simple pendulum using a Taylor series expansion,” European Journal of Physics, vol. 32, no. 5, pp. 1303–1310, 2011. View at: Publisher Site  Google Scholar
 R. E. Mickens, Oscillations in Planar Dynamic Systems, World Scientific, Singapore, Singapore, 1996.
 E. W. Weisstein, “ArithmeticGeometric Mean,” From MathWorld—A Wolfram Web Resource, http://mathworld.wolfram.com/ArithmeticGeometricMean.html. View at: Google Scholar
 C. G. Carvalhaes and P. Suppes, “Approximation for the period of the simple pendulum based on the arithmeticgeometric mean,” American Journal of Physics, vol. 76, no. 12, pp. 1150–1154, 2008. View at: Publisher Site  Google Scholar
 M. Febbo, “A finite extensibility nonlinear oscillator,” Applied Mathematics and Computation, vol. 217, no. 14, pp. 6464–6475, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH
Copyright
Copyright © 2012 A. Beléndez et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.