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.
(a)
(b)
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.