Journal of Applied Mathematics

Volume 2012 (2012), Article ID 518684, 7 pages

http://dx.doi.org/10.1155/2012/518684

## High-Order Energy Balance Method to Nonlinear Oscillators

Faculty of Aeronautics and Astronautics, Istanbul Technical University, Maslak, 34469 Istanbul, Turkey

Received 1 April 2012; Revised 25 June 2012; Accepted 9 July 2012

Academic Editor: Livija Cveticanin

Copyright © 2012 Seher Durmaz and Metin Orhan Kaya. 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.

#### Abstract

Energy balance method (EBM) is extended for high-order nonlinear oscillators. To illustrate the effectiveness of the method, a cubic-quintic Duffing oscillator was chosen. The maximum relative errors of the frequencies of the oscillator read 1.25% and 0.6% for the first- and second-order approximation, respectively. The third-order approximation has an accuracy as high as 0.008%. Excellent agreement of the approximated frequencies and periodic solutions with the exact ones is demonstrated for several values of parameters of the oscillator.

#### 1. Introduction

A great deal of work has been devoted to the nonlinear problems encountered in the fields of applied mathematics, physics, and engineering sciences. In general, the analytical approximation to the solution of a given nonlinear problem is difficult, sometimes impossible; plenty of techniques based on numerical methods have been implemented. Among them are the variational iteration method [1, 2], the harmonic balance method [3, 4], and energy balance method [5–8] developed to solve nonlinear differential equations. In this study, we have investigated the application of high-order energy balance method to cubic-quintic Duffing oscillator. The nonlinear frequencies are calculated for the first-, second-, and third-order EBM and compared with the results of different techniques.

#### 2. The Basic Idea of Energy Balance Method

This section briefly introduces energy balance method proposed by He [5]. In this method, a variational principle for the oscillation is established, then the corresponding Hamiltonian is considered from which the angular frequency can be easily obtained by Galerkin method.

Let us consider the motion of a general oscillator with the initial conditions in the form where is the initial amplitude.

Its variational can be written as Here is the period of the nonlinear oscillation and .

The Hamiltonian of (2.1) can be written in the form: In (2.2) the kinetic energy and potential energy can be, respectively, expressed as , .

Throughout the oscillation since the system is conservative, the total energy remains unchanged during the motion; the Hamiltonian of the oscillator becomes a constant value, For the first-order approximation, the following trial function can be assumed: Substituting (2.5) into (2.3) yields the following residual: The residual is forced to zero, in an average sense, by setting weighted integrals of the residual to zero where is a set of weighting function (or test).

There are lots of weighting functions, that is, Galerkin, least squares, collocation and so forth. In this study, we used Galerkin method as a weighting function.

#### 3. High-Order Energy Balance Method

In order to extend He’s energy balance method, let us assume that the solution of (2.1) can be expressed as From the initial conditions, the coefficients should satisfy the following constrain: One of these parameters can be chosen as a dependent parameter. Hence, By inserting (3.1) into (2.7), the following systems can be obtained:

#### 4. Example

A cubic-quintic Duffing oscillator is considered. In the following sections, the nonlinear frequencies will be compared with the results of different techniques to illustrate the efficiency and accuracy of energy balance method.

The governing differential equation of this oscillator is in the form of with the initial conditions The Hamiltonian of (4.1a) is given as follows: For the first-order approximation, assume that is in the following form: Substituting the first approximation into (4.2) yields First-order approximation can be obtained by setting The amplitude-frequency relationship for the first approximation is obtained as To obtain a more accurate result, let us define as follows: Substituting (4.7) into (4.2) results in the following residual: We set

By solving (4.9a)-(4.9b) simultaneously, one can obtain the second-order approximate amplitude-frequency relation. For different values, the approximate frequencies are given in Table 1.

Moreover, the accuracy of results will be further improved by defining in the following form: Substituting (4.10) into (4.2), we get the following residual for the third-order approximation: Inserting (4.11) into (3.3) for and using the same procedure explained above, we get three weighted integrals. Solving these three equations simultaneously, the amplitude-frequency relation for the third-order approximation is obtained. For higher-order approximations, the similar procedures can be applied, however, the accuracy of the third-order approximation is appropriate for several values of parameters , , and .

In the following, the nonlinear frequencies for the cubic-quintic oscillator are calculated for two different cases: (i) and (ii) .

The first case considered here corresponds to the cubic Duffing oscillator. For , the results of the nonlinear frequencies are obtained by the first-, second-, and third-order energy balance method and compared with the results of [4, 9]. In the second case, for , the nonlinear frequencies are given in Table 2. Additionally, the numerical solution for all cases is acquired by standard Runge-Kutta method (R-K).

In Table 1, the relative errors for the first-order approximation read 1.25%, while this error reduces to 0.59% in the second approximation. We observe that the differences between the third-order and the exact frequencies are sufficiently small.

As seen in Table 2, the results of third order approximation are in very good agreement between the numerical results. The comparison of approximate and numerical solutions can also be found in Figure 1. It can be seen that the first-, and the second-order results have slight differences compared to the numerical solution. However, the third-order approximation is overlapping with the numerical solution.

#### 5. Conclusion

In this paper, energy balance method is extended for high-order solutions. The first-order approximate frequency for Duffing oscillator gives 1.25% relative error, while the second-, and the third-order approximated frequencies reach 0.59% and 0.008% relative errors, respectively. Moreover, relative errors in high-order energy balance reduce to smaller values than global error minimization and harmonic balance methods. Consequently, we can state that extended method is very effective and convenient for the cubic-quintic Duffing oscillator.

#### References

- J.-H. He, “Variational iteration method—some recent results and new interpretations,”
*Journal of Computational and Applied Mathematics*, vol. 207, no. 1, pp. 3–17, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - V. Marinca and N. Herisanu, “Periodic solutions for some strongly nonlinear oscillations by He's variational iteration method,”
*Computers & Mathematics with Applications*, vol. 54, no. 7-8, pp. 1188–1196, 2007. View at Publisher · View at Google Scholar - A. Beléndez, E. Gimeno, T. Beléndez, and A. Hernández, “Rational harmonic balance based method for conservative nonlinear oscillators: application to the Duffing equation,”
*Mechanics Research Communications*, vol. 36, no. 6, pp. 728–734, 2009. View at Publisher · View at Google Scholar - B. S. Wu, W. P. Sun, and C. W. Lim, “An analytical approximate technique for a class of strongly non-linear oscillators,”
*International Journal of Non-Linear Mechanics*, vol. 41, no. 6-7, pp. 766–774, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - J. H. He, “Preliminary report on the energy balance for nonlinear oscillations,”
*Mechanics Research Communications*, vol. 29, no. 2-3, pp. 107–111, 2002. View at Publisher · View at Google Scholar · View at Scopus - H.-L. Zhang, “Periodic solutions for some strongly nonlinear oscillations by He's energy balance method,”
*Computers & Mathematics with Applications*, vol. 58, no. 11-12, pp. 2480–2485, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - M. Momeni, N. Jamshidi, A. Barari, and D. D. Ganji, “Application of He's energy balance method to Duffing-harmonic oscillators,”
*International Journal of Computer Mathematics*, vol. 88, no. 1, pp. 135–144, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - A. Yildirim, H. Askari, Z. Saadatnia, M. KalamiYazdi, and Y. Khan, “Analysis of nonlinear oscillations of a punctual charge in the electric field of a charged ring via a Hamiltonian approach and the energy balance method,”
*Computers & Mathematics with Applications*, vol. 62, no. 1, pp. 486–490, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - Y. Farzaneh and A. Akbarzadeh Tootoonchi, “Global error minimization method for solving strongly nonlinear oscillator differential equations,”
*Computers & Mathematics with Applications*, vol. 59, no. 8, pp. 2887–2895, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH