Journal of Applied Mathematics

Volume 2013 (2013), Article ID 620267, 6 pages

http://dx.doi.org/10.1155/2013/620267

Research Article

Optimal Variational Method for Truly Nonlinear Oscillators

Vasile Marinca1,2 and Nicolae Herişanu1,2

¹Faculty of Mechanical Engineering, “Politehnica” University of Timişoara, Bd. M. Viteazu, 1, 300222 Timişoara, Romania

²Center of Advanced Research in Engineering Sciences, Romanian Academy, Timişoara Branch, Bd. M. Viteazu, 24, 300223, Timişoara, Romania

Received 6 August 2012; Accepted 19 November 2012

Academic Editor: Kale Oyedeji

Copyright © 2013 Vasile Marinca and Nicolae Herişanu. 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.

The Optimal Variational Method (OVM) is introduced and applied for calculating approximate periodic solutions of “truly nonlinear oscillators”. The main advantage of this procedure consists in that it provides a convenient way to control the convergence of approximate solutions in a very rigorous way and allows adjustment of convergence regions where necessary. This approach does not depend upon any small or large parameters. A very good agreement was found between approximate and numerical solution, which proves that OVM is very efficient and accurate.

1. Introduction

It is known that the study of nonlinear differential equations is restricted to a variety of special classes of equations and the method of solution usually involves one or more techniques to achieve analytical approximations to the solutions. Many researchers have recently paid much attention to find and develop approximate solutions. Perturbation methods have been successfully employed to determine approximate solutions to weakly nonlinear problems [1–3]. But the use of perturbation theory in many problems is invalid or it simply breaks down for parameters beyond a certain specified range. Therefore, new analytical techniques should be developed to overcome these shortcomings. Such a new technique should work over a large range of parameters and yield accurate analytical approximate solutions beyond the coverage and ability of the classical perturbation methods. There has been a great need for effective algorithms to avoid the work required by traditional techniques, but it is difficult to obtain convergent results in the cases of strong nonlinearity.

Recently, many new approaches have been proposed for this purpose, such as various modified Lindstedt-Poincare methods [4], some linearization methods [5, 6], the Adomian decomposition method [7], the optimal homotopy asymptotic method [8, 9], the optimal variational iteration method [10], the energy balance method [11], and so on. A variational principle for nonlinear oscillations by constructing the Hamiltonian was studied in [11]. Variational principles have a great importance in physics and engineering since they establish connections between these disciplines and their applications are useful in devising various approximate techniques. Variational methods have been and continue to be popular tools for nonlinear analysis. They combine physical insights into the nature of the solution of the problem and the solutions obtained using possible trial functions are the best. It is known that computing a Lagrangian for dynamical systems with more general Newtonian forces is nowadays applicable only to systems with force derivable from a potential function (basically, conservative systems). Strictly speaking, conservative dynamical systems do not exist in our Newtonian environment. As a result, the Lagrangian representation of conservative Newtonian systems is in general only a crude approximation of physical reality. The problem of the existence of a Lagrangian, Hamiltonian, or Routhian can be studied today with a variety of modern and sophisticated mathematical tools which include the use of functional analysis, prolongation theory, and differential geometry, to cite only a few. This question is called “Inverse Problem in Newtonian Mechanics” [12] and consists in the identification of the methods for the construction of a Lagrangian, Hamiltonian, or Routhian form given equations of motion.

In this paper we construct accurate approximations to periodic solutions and frequencies of the so-called truly nonlinear oscillator (TNO). Following Mickens and Oyedeji [13, 14], the most general form of a TNO is given by the following differential equation:

where the dot denotes the derivative with respect to variable , is a positive arbitrary parameter, and the functions and have the properties:

and does not have for small a dominant term proportional to .

In the present work, we consider , and therefore the truly nonlinear oscillator is modeled by the following nonlinear differential equation:

subject to initial conditions

where , , and are known parameters.

In (3), there exists no small or large parameter.

2. Basic Idea of Optimal Variational Method and Solution

In order to develop an application of the OVM, we consider the following differential equation:

with the initial conditions given by (4) and being an arbitrary odd function.

Introducing a new independent variable and a new unknown as

where is the frequency of the system (5), then this becomes

where prime denotes derivative with respect to the new variable .

The initial conditions (4) become

The variational principle for (7) can be easily established if there exists a function

which admits as extremals the solutions of (7) and (8) where is the Lagrangian of the system (7):

The function is given by the equation

We assume that the approximate periodic solutions of (7) and (8) can be expressed as

where are arbitrary unknowns at this moment and is a positive integer number. Choosing the solution (12) has been made in accordance to the properties (2).

Substituting (12) into (9) results in

Applying the Ritz method [9], we require

From (14) and from the initial condition (8₁) which becomes

we can obtain optimally the parameters , , and the frequency .

We remark that the choice of the approximate solution (12) is not unique. We can alternatively choose another expression of the approximate periodic solution in the form

and so on. With the parameters (called convergence-control parameters) and the frequency known, the approximate periodic solutions is well determined.

The validity of the proposed approach is illustrated on the TNO given by (3). Using the transformations (6), equation (3) can be written in the form

and the Lagrangian becomes

If we consider into (12), then the approximate periodic solutions become

Now, substituting (19) into (18) and (13), we have successively

To calculate the middle term into (18) we use the series expansion [15]:

where . If we denote where is given by (19), then using the relation

with , , we can write successively

Substituting (19) and (24) into (13) we obtain

From (25), (14), and (15) we can obtain the unknowns , , , and .

We remark that with the approximation given by (19), there are four things to be calculated: the optimal values of the convergence-control parameters , , , and the frequency . As it can be seen, the complexity of the equations is such that only numerical solutions can be found for particular values of the parameters , , and , which will be further developed.

3. Numerical Examples

In order to show the validity and accuracy of the OVM, we consider the following cases.

(1) In the first case we consider and and we obtain

 The approximate solution of (3) becomes

 The value of obtained by numerical integration is .(2) In the second case, for and we have

 and therefore the approximate solution becomes

 In this case, the value of obtained by numerical integration is .(3) In the third case, we consider and and thus

 The approximate solution of (3) can be written as

 The result of numerical integration for the frequency is in this case .(4) In the case and we get

 and the approximate solution of (3) will be

 By numerical integration, we get in this case the frequency .(5) For and , we obtain

 and the approximate solution of (3) becomes

 In this case the numerical integration result for the frequency is .(6) In the last case, we consider and and therefore

 and the approximate solution of (3) becomes

 For comparison purposes, the frequency of the system obtained directly by numerical integration is .

Figures 1, 2, 3, 4, 5, and 6 present a comparison between the present solutions (27)–(37), respectively, and the numerical integration results for (3) and (4). Thus, it is easier to emphasize the accuracy of the obtained results, since within these graphical representations, the analytical results obtained through OVM are nearly identical with numerical ones. Also, the approximate frequencies obtained by OVM are in very good agreement with those obtained by numerical integration, which also proves the validity of the approximate results.

Figure 1: Comparison between the approximate solution (27) and numerical solution of (3) in the case , , red solid line: numerical integration results, blue dashed line: approximate solution.

Figure 2: Comparison between the approximate solution (29) and numerical solution of (3) in the case , , red solid line: numerical integration results, blue dashed line: approximate solution.

Figure 3: Comparison between the approximate solution (31) and numerical solution of (3) in the case , , red solid line: numerical integration results, blue dashed line: approximate solution.

Figure 4: Comparison between the approximate solution (33) and numerical solution of (3) in the case , , red solid line: numerical integration results, blue dashed line: approximate solution.

Figure 5: Comparison between the approximate solution (35) and numerical solution of (3) in the case , , red solid line: numerical integration results, blue dashed line: approximate solution.

Figure 6: Comparison between the approximate solution and numerical solution of (37) in the case , , red solid line: numerical integration results, blue dashed line: approximate solution.

4. Conclusions

In this paper we introduce the Optimal Variational Method to propose a new analytic approximate periodic solution to TNO. Our procedure is valid even if the nonlinear equation does not contain small or large parameters. OVM provides us with a simple and rigorous way to control and adjust the convergence of a solution through several convergence-control parameters whose values are optimally determined. This new method is very rapid and effective, and we prove it by comparing the approximate periodic solutions and frequencies obtained through the proposed method with numerical integration results. An excellent agreement has been demonstrated between the analytical and numerical integration results, and also for large values of the initial amplitudes, which validates the effectiveness of the proposed method. The proposed procedure can also be used to find analytical approximate solutions to other classes of conservative oscillators.

The main advantage of the method consists in that it provides us with a great freedom in choosing the approximate periodic solution dependent on an arbitrary number of initially unknown parameters. It is interesting to remark that unlike other known approximate methods applicable for such problems (such as the harmonic balance method, the multiple scales method, and so on), the frequency of the system is not obtained imposing some conditions to avoid secular terms, which is a usual procedure within other methods, but this frequency results directly from the conditions to minimize the residual functional along with the initial conditions.

Acknowledgment

This work was supported by CNCS-UEFISCDI, Project no. PN-II-ID-PCE-2012-4-0358.

References

A. H. Nayfeh and D. T. Mook, Nonlinear Oscillations, John Wiley and Sons, New York, NY, USA, 1979.

P. Hagedorn, Nonlinear Oscillations, vol. 10, Clarendon Press, Oxford, UK, 1988.

R. E. Mickens, Oscillations in Planar Dynamical Systems, vol. 37, World Scientific, Singapore, 1996.

J.-H. He, “Modified Lindstedt-Poincare methods for some strongly nonlinear oscillators. Part.I. expansion of a constant,” International Journal of Non-Linear Mechanics, vol. 37, no. 2, pp. 309–314, 2002.

A. Beléndez, C. Pascual, C. Neipp, T. Beléndez, and A. Hernández, “An equivalent linearization method for conservative nonlinear oscillators,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 9, no. 1, pp. 9–17, 2008.

J. I. Ramos, “Linearized Galerkin and artificial parameter techniques for the determination of periodic solutions of nonlinear oscillators,” Applied Mathematics and Computation, vol. 196, no. 2, pp. 483–493, 2008.

G. Adomian, “A review of the decomposition method in applied mathematics,” Journal of Mathematical Analysis and Applications, vol. 135, no. 2, pp. 501–544, 1988.

V. Marinca and N. Herişanu, “Determination of periodic solutions for the motion of a particle on a rotating parabola by means of the optimal homotopy asymptotic method,” Journal of Sound and Vibration, vol. 329, no. 9, pp. 1450–1459, 2010.

V. Marinca and N. Herişanu, Nonlinear Dynamical Systems in Engineering. Some Approximate Approaches, Springer, Berlin, Germany, 2011.

N. Herişanu and V. Marinca, “A modified variational iteration method for strongly nonlinear problems,” Nonlinear Science Letters A, vol. 1, pp. 183–192, 2010.

J. H. He, “Preliminary report on the energy balance for nonlinear oscillations,” Mechanics Research Communications, vol. 29, pp. 107–111, 2003.

V. Obădeanu and V. Marinca, The Inverse Problem in Analytic Mechanics, University of Timişoara, Timişoara, Romania, 1992.

R. E. Mickens and K. Oyedeji, “Construction of approximate analytical solutions to a new class of nonlinear oscillator equation,” Journal of Sound and Vibration, vol. 162, no. 4, pp. 579–582, 1985.

R. E. Mickens, “A generalized iteration procedure for calculating approximations to periodic solutions of ‘truly nonlinear oscillators’,” Journal of Sound and Vibration, vol. 287, no. 4-5, pp. 1045–1051, 2005.

R. E. Mickens, “Iteration method solutions for conservative and limit-cycle X1/3 force oscillators,” Journal of Sound and Vibration, vol. 292, no. 3–5, pp. 964–968, 2006.