Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 620591, 7 pages

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

Research Article

Energy Method to Obtain Approximate Solutions of Strongly Nonlinear Oscillators

Alex Elías-Zúñiga and Oscar Martínez-Romero

Centro de Innovación en Diseño y Tecnología, Tecnológico de Monterrey, Campus Monterrey, E. Garza Sada 2501 Sur, 64849 Monterrey, NL, Mexico

Received 13 April 2013; Revised 15 November 2013; Accepted 18 November 2013

Academic Editor: Xinzhu Meng

Copyright © 2013 Alex Elías-Zúñiga and Oscar Martínez-Romero. 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.

We introduce a nonlinearization procedure that replaces the system potential energy by an equivalent representation form that is used to derive analytical solutions of strongly nonlinear conservative oscillators. We illustrate the applicability of this method by finding the approximate solutions of two strongly nonlinear oscillators and show that this procedure provides solutions that follow well the numerical integration solutions of the corresponding equations of motion.

1. Introduction

This paper focuses on introducing a nonlinearization procedure to derive approximate solutions of strongly conservative nonlinear oscillators by replacing the system potential energy by equivalent representation forms that are based on cubic- or quartic-like polynomial expressions. The main motivation for finding an equivalent representation form of the potential energy terms is based on the different approaches that use equivalent restoring forces to derive approximate solutions of nonlinear oscillators. For instance, Yuste and Sánchez in [1, 2] used a cubication approach to replace the restoring forces by equivalent cubic polynomial expressions by applying the weighted mean-square and the harmonic balance methods. By using Chebyshev polynomial expansions, Beléndez and coworkers replaced the restoring forces by an equivalent cubic-like representation form that transformed the original equation of motion into a Duffing equation [3–5]. Then, Elas-Zúñiga and coworkers used this idea and considered additional Chebyshev polynomial expansion terms to transform the restoring forces into cubic-quintic polynomial representation forms that improved the accuracy of the approximate solutions of different types of strongly nonlinear oscillators [6, 7]. Here in this paper, instead of finding equivalent representation forms of the system restoring forces, we focus our attention on replacing the system potential energy by an equivalent energy expression that will be used to derive the corresponding approximate solutions. The accuracy of our proposed procedure will be illustrated by finding the approximate solutions of two nonlinear oscillators: (a) the cubic-quintic Duffing and (b) the finite extensibility nonlinear oscillator (FENO). We will show that our proposed approach provides solutions that follow well the numerical integration solution of the corresponding equations of motion. In fact, we have found that, for the case of the cubic-quintic oscillator, the nonlinear transformation form that is based on a quartic polynomial provides approximate angular frequency values whose maximum errors, when compared to the exact ones, do not exceed 0.004%. Similarly, we found for the FENO oscillator that when the initial oscillation amplitude approaches to one, the errors achieved do not exceed the percentage value of 0.3787%.

2. A Nonlinearization Energy Procedure

Conservative nonlinear oscillators are commonly modelled by equations of the form

where is the conservative system restoring force and is the initial oscillation amplitude. Integration of (1) with respect to the displacement leads to its energy equation of motion

Here, and represent the velocity and the system potential energy terms, respectively, and is an integration constant whose value is determined from the given initial conditions. Now we assume that the potential energy term can be replaced by an equivalent representation form which can be obtained by applying available techniques such as the nonlinearization approach [8], the Chebyshev polynomials of the first kind [3–5], the weighted mean-square method [2, 9], or other techniques available in the literature. In this paper, we adopt Chebyshev polynomial and used its general definition to write in equivalent form the potential energy terms [10]

where

and the first four polynomials of the first kind are given by

Thus, the equivalent representation form of the system potential energy can be written as

where the forms of the parameters , , , depend on the Chebyshev coefficient terms. Therefore, the energy equation (2) can be equivalently represented by the following expression:

where . However, this representation form (7) needs to be further transformed to a more tractable expression. To achieve this transformation, we can use the nonlinearization processes discussed in [6, 7]. Before we perform this transformation, we need to identify which polynomial order has to be used during the application of the nonlinearization processes, that is, a cubic-like or a quartic-like form. To better understand the implication of the choice of the polynomial order, we recall that the Duffing equation

has an energy equation representation given as

where , , and are known parameter values. Therefore, it could be obvious that a quartic-like polynomial form has to be chosen to replace the potential energy terms of (7). However, we explore here two different options to replace the potential energy terms: (a) a cubic-like polynomial and (b) a quartic-like polynomial expression. If we use option (a), the nonlinearization processes will lead us to the following expression:

where the parameters , , and can be found from

Thus, the equivalent representation of (2) has the form

However, if we adopt the quartic-like polynomial form, (2) can be transformed to

where

with

Of course, other possible polynomial representation forms could be assumed, however, we focus here on the usage of the cubic-like and the quartic-like polynomial forms to determine the angular frequency expressions of nonlinear oscillators under consideration. Notice that in (10) and (14), the value of must be appropriately chosen to ensure that the equivalent potential energy terms are qualitatively and quantitatively similar to those of the original equations of motion. In this sense, we have followed the approach developed by Cai and coworkers in [8].

We will next study the feasibility of our proposed nonlinearization energy approach to derive the equivalent energy equations of the cubic-quintic Duffing equation and then we will explore the solution of the finite extensibility nonlinear oscillator (FENO).

3. The Cubic-Quintic Duffing Oscillator

Although the exact solution of this nonlinear oscillator is well known [11], we choose its equation of motion to illustrate the applicability of our proposed mathematical technique to learn more about some possible issues that could arise during its implementation to solve nonlinear systems. This oscillator has an equation of motion given as

where denotes the displacement of the system and , , and are system constant parameters. If we set , then (16) can be written as

Here , , and . To apply our nonlinear energy approach, we first integrate (17) with respect to . This step provides the following energy equation:

where the integration constant has the value of

First, we follow option (a) and transform the potential energy terms of (18) into a cubic polynomial by following our proposed energy nonlinearization method. Then, we use (10) and (11) to write the potential energy in equivalent form as

where

where is a parameter whose value satisfies (10). Thus, the equivalent energy transformation form of (18) is given as:

To find the angular frequency of (22), we next proceed with its integration and find that

where

with . The integration of (23) provides the following time expression:

Here represents the incomplete elliptic integral of the first kind with modulus . Thus, the approximate angular frequency value of the cubic-quintic oscillator is given as

where on . Before we assess the accuracy of our derived angular frequency expression (26), we next follow option (b) and write the potential energy as a function of a quartic polynomial expression

Notice that the parameters , , and can be determined by using (14) and (15) which yields

where the value of must satisfy (14). Now the solution of

yields

where

Thus, the approximate angular frequency expression is given as , where . Also notice that (25) and (30) can be rewritten in different forms by putting as a function of time and by using the elliptic function identities. However, we will not elaborate on these solutions here.

We next compare our derived solutions with respect to the exact one of (17) found in [11]. Figure 1 illustrates the percentage error attained by comparing the exact angular frequency values of (17) with those obtained by using (26) plotted against the initial oscillation amplitude values. In this case, we use the parameter values of , , and with . We can see from Figure 1 that the maximum error values do not exceed 1.194%. We next use the quartic nonlinearization approximate angular frequency value given by (30) and compute the corresponding error values. One can notice from Figure 2, for the same system parameter values of , , and , with , that now the maximum error values do not exceed 0.004%. This represents a great improvement when compared to the estimated values of . We believe that the improvement achieved by using (30) could be due to the usage of a fourth-order polynomial to approximate the potential energy expression which is closer to the order of the sextic polynomial that is used to derive, through the energy equation, the exact solution of (16). To further assess the precision of our proposed approach, we shall next explore the solution of a finite extensibility nonlinear oscillator (FENO).

Figure 1: Percentage relative angular frequency error values curve versus the initial oscillation amplitude . The parameter values are , , , and .

Figure 2: Percentage relative angular frequency error values curve versus the initial oscillation amplitude . The parameter values are , , and .

4. A Finite Extensibility Nonlinear Oscillator

This oscillator is used to describe chain dynamics of polymer-like materials subjected to high extensibility values [12]. Its equation of motion has the form

where is the initial oscillation amplitude with . If we introduce the transformation , then (33) can be written as

It is well known that when , this oscillator tends to behave as an anharmonic one and thus, depending on the methods used, its approximate solutions to predict its dynamical response could diverge from the numerical one. See, for instance, [7–13]. The aim of this section focuses on exploring the applicability of our proposed nonlinearization energy approach to develop an approximate solution to (34) valid for the whole interval value of .

First, we derive the energy expression of (34) by performing its integration with respect to . This step leads to the following equation:

where

Since we need to write (34) in a polynomial-like form, then we follow our proposed approach and use the Chebyshev polynomial expansion (3) to get that

where

We now adopt options (a) and (b) and determine the equivalent system potential energy. By using the cubic transformation form, we obtain that

where

If we now follow option (b), the equivalent system potential energy is given as

for which , , and .

Figure 3 illustrates the percentage error attained by comparing the angular frequency values of and with respect to the numerical integration ones plotted versus the initial oscillation amplitude . It is interesting to note that when approaches to one, that is, , the percentage error curves tend to the values of 0.3787% and 0.2040%, respectively. Figure 4 shows a comparison among our derived solutions and the angular frequency value obtained by Beléndez et al. [13] in which

It is clear from Figure 4 that when approaches to one, the highest percentage error is obtained from (42). Therefore, we can conclude that our derived angular frequency solutions provide error values that are remarkably lower than those previously reported in [7–13]. Furthermore, the maximum error values of both curves do not exceed 3.2%. Here the fitting parameter has the value of 1.01195.

Figure 3: Percentage relative angular frequency error value curves versus the initial oscillation amplitude . Here the fitting parameter has the value of 1.01195 and .

Figure 4: Comparison of the percentage relative angular frequency error value curves versus the initial oscillation amplitude .

We next focus on the qualitative behavior of our derived solution by plotting its amplitude-time response curve and compare this with its numerical integration solution. In this case, we have used (25) and (30) and considered half periods to ensure that our FENO system satisfies the initial conditions. In fact, we found that, for given values of bigger than one, the estimated values of and computed from (25) and (30), respectively, become imaginary. Therefore and in order to plot the amplitude-time response curves, the values of must satisfy the condition . Figure 5 shows the amplitude-time response curves obtained from (25) and the numerical integration solution of (35). Notice that both solutions have the same qualitative behavior even though the assumed initial system oscillation amplitude has a value close to one; that is, .

Figure 5: Computed amplitude-time response curves of (34) and (25). Here the solid line represents the numerical integration solution of (34) while the black dashed line represents our derived energy approximate solution (25), by considering an initial oscillation amplitude value of .

5. Conclusions

We have introduced a nonlinearization energy method to predict the qualitative and quantitative behavior of strongly nonlinear oscillators that is based on the nonlinear transformation of the system energy equation. We have assumed that the original system potential energy can be equivalently written as a cubic or quartic polynomial representation by using Chebyshev polynomials and the nonlinearization method. We have shown that this procedure leads to equations that can be exactly integrated by using elliptic functions. Then, we have used these equations to study the dynamical response of two nonlinear oscillators: one related to the cubic-quintic Duffing equation and the other related to a finite extensibility nonlinear oscillator (FENO). We have found that our proposed approach follows well the numerical integration solution of the corresponding equations of motion. In fact and for the case of the cubic-quintic oscillator, the nonlinear transformation approach based on a quartic polynomial provides approximate angular frequency values whose maximum errors do not exceed 0.004%. We believe that the improvement achieved by replacing the system potential energy of (18) by an equivalent potential energy form could be due to the usage of a fourth-order polynomial expression which is closer to the order of the sextic polynomial that is needed to derive, through the energy equation, the exact solution of (17).

Similarly, we found, for the FENO oscillator, that when the initial oscillation amplitude approaches to one, that is, when , the errors achieved do not exceed the percentage value of 0.3787%. Furthermore, we have proved that the qualitative behavior of our derived solutions follows well that of the corresponding numerical integration solutions. Therefore, we can conclude that our nonlinearization energy approach provides results that describe well the dynamical response of these strongly nonlinear oscillators.

It is evident that our proposed approach could be used to study the dynamical response of other dynamical oscillators and that it can be expanded to obtain approximate solutions of dynamical system with two or more degrees of freedom.

Conflict of Interests

The authors declare that they have no conflict of interests with any mentioned entities in the paper.

Acknowledgments

This work was funded by Tecnológico de Monterrey, Campus Monterrey, through the Research Chair in Nanomaterials for Medical Devices and the Research Chair in Intelligent Machines. Additional support was provided by the European Union Seventh Framework Programme (FP7 PEOPLE-2009) under the Grant Agreement IRSES no. 247476 and from Consejo Nacional de Ciencia y Tecnología (Conacyt), Mexico.

References

S. B. Yuste and Á. M. Sánchez, “A weighted mean-square method of “cubication” for nonlinear oscillators,” Journal of Sound and Vibration, vol. 134, no. 3, pp. 423–433, 1989.

S. B. Yuste, “‘Cubication’ of nonlinear oscillators using the principle of harmonic balance,” International Journal of Non-Linear Mechanics, vol. 27, no. 3, pp. 347–356, 1992.

A. Belendez, M. L. Alvarez, E. Fernandez, and I. Pascual, “Cubication of conservative nonlinear oscillators,” European Journal of Physics, vol. 30, pp. 973–981, 2009.

A. Beléndez, D. I. Méndez, E. Fernández, S. Marini, and I. Pascual, “An explicit approximate solution to the Duffing-harmonic oscillator by a cubication method,” Physics Letters A, vol. 373, no. 32, pp. 2805–2809, 2009.

A. Beléndez, G. Bernabeu, J. Francés, D. I. Méndez, and S. Marini, “An accurate closed-form approximate solution for the quintic Duffing oscillator equation,” Mathematical and Computer Modelling, vol. 52, no. 3-4, pp. 637–641, 2010.

A. Elías-Zúñiga, O. Martínez-Romero, and R. K. Córdoba-Díaz, “Approximate solution for the Duffing-harmonic oscillator by the enhanced cubication method,” Mathematical Problems in Engineering, vol. 2012, Article ID 618750, 12 pages, 2012.

A. Elías-Zúñiga and O. Martínez-Romero, “Accurate solutions of conservative nonlinear oscillators by the enhanced cubication method,” Mathematical Problems in Engineering, vol. 2013, Article ID 842423, 9 pages, 2013.

J. Cai, X. Wu, and Y. P. Li, “An equivalent nonlinearization method for strongly nonlinear oscillations,” Mechanics Research Communications, vol. 32, pp. 553–560, 2005.

S. C. Sinha and P. Srinivasan, “A weighted mean square method of linearization in non-linear oscillations,” Journal of Sound and Vibration, vol. 16, pp. 139–148, 1971.

A. Gil, J. Segura, and N. M. Temme, Numerical Methods for Special Functions, SIAM Society for industrial and applied mathematics, 2007.

A. Elías-Zúñiga, “Exact solution of the cubic-quintic Duffing oscillator,” Applied Mathematical Modelling, vol. 37, no. 4, pp. 2574–2579, 2013.

M. Febbo, “A finite extensibility nonlinear oscillator,” Applied Mathematics and Computation, vol. 217, no. 14, pp. 6464–6475, 2011.

A. Beléndez, E. Arribas, J. Francés, and I. Pascual, “Comments on a finite extensibility nonlinear oscillator,” Applied Mathematics and Computation, vol. 218, no. 10, pp. 6168–6175, 2012.