Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 842423, 9 pages

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

Research Article

Accurate Solutions of Conservative Nonlinear Oscillators by the Enhanced Cubication Method

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

Departamento de Ingeniería Mecánica, Tecnológico de Monterrey, Campus Monterrey, E. Garza Sada 2501 Sur, 64849 Monterrey, NL, Mexico

Received 2 December 2012; Accepted 8 February 2013

Academic Editor: Erwei Bai

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.

The enhanced cubication method is applied to develop approximate solutions for the most common nonlinear oscillators found in the literature. It is shown that this procedure leads to amplitude-time response curves and angular frequency values with maximum relative errors lower than those found by previously developed approximate solutions.

1. Introduction

In this paper, we use the enhanced-cubication method to develop the approximate solution of some conservative nonlinear oscillators of the form

where is the system restoring force which could have rational or irrational elastic terms, and is the oscillation amplitude. We recall that the enhanced-cubication method uses the Chebyshev polynomial expansion to replace the system restoring force by an equivalent one. Then, the nonlinearization method is used to find and equivalent Duffing equation for (1) that is valid for the complete range of oscillation amplitudes [1]. Elías-Zúñiga and coworkers used this technique to derive the approximate solution of the Duffing-harmonic oscillator of the form

where , , and are constant system parameters [1]. They found by using this method that the approximate angular frequency values agree well with the exact ones for the whole range of oscillation amplitude values with a maximum relative error value that is less than %.

Here in this paper, we examine the application of the enhanced-cubication method to obtain approximate solutions of the general pendulum equation, the generalized power-form elastic term oscillator, the general restoring force term oscillators, and of the singular oscillator that arises in plasma physics. To assess the accuracy of the enhanced-cubication method, we compare our derived approximate angular frequency values with those obtained by numerical integration. It is shown that in most cases, our derived solutions have maximum relative errors that are lower than those obtained by using other solution techniques.

We will next start with the application of our enhanced cubication approach by deriving the approximate solution of the general pendulum equation.

2. A General Pendulum Equation

We first consider the general pendulum equation of the form [2]

with initial conditions and , and system parameters and . Next, we use the following transformation and rewrite (3) as

Then, we write in equivalent form by using the Chebyshev polynomials of the first kind [3–6]:

where

Thus, the equivalent restoring force can be written as

where the coefficients , , and depend on , the amplitude of oscillation, , and the Chebyshev coefficient terms. Therefore, the nonlinear differential equation (4) can be replaced by an equivalent cubic-quintic Duffing equation of the form

Since the cubication procedure requires a cubic polynomial representation of the restoring force, we shall next transform the restoring force term of (9) into an equivalent cubic polynomial

where and can be found by using the equivalent nonlinearization method [2]. In this method, we replace the terms by the cubic polynomial that satisfies

By following this procedure, the coefficients and can be determined by using the expressions

where the value of is fitted to satisfy (11). Thus, the equivalent representation form of (4) is given by

whose exact angular frequency-amplitude relationship is given by

where

Here, the coefficients , , and are obtained by using (6)–(8). It is easy to show that these coefficients are given as

where and are the first and second order Bessel functions of the first kind, and and are the complete elliptic integral of the first and second kinds for the modulus , respectively. Figure 1 illustrates the relative errors achieved by comparing the approximate angular frequency values computed from (14) with the exact ones, obtained by using Radhakrishnan et al. procedure introduced in [7]. As we may see from Figure 1 and for , the relative error values are not bigger than 0.2%.

Figure 1: Prediction of the error attained by computing the general pendulum angular frequency oscillator values by the enhanced cubication method. Notice that the maximum relative error attained between the approximate and the exact angular frequency values does not exceed of 0.2% for .

To further assess the applicability of the enhanced cubication approach, we next derive the solution of an oscillator with fractional restoring form.

3. A Generalized Power-Form Elastic Term Oscillator

In this case, the restoring force term in (3) has the form

The restoring force given by (18) is similar in form to the one proposed by Pilipchuk [9], and due to its oddness with respect to time, the exponent can continuously take any nonnegative real value, such as odd, even, rational, or irrational that is, . Notice that is a constant parameter whose value depends on the physical system under consideration. We next use the transformation and write the corresponding equation of motion in the form

Then, we use Chebyshev polynomial expansion and write the restoring force as

where

Here, is the Euler gamma function. Table 1 shows the corresponding values of the fitting parameter that satisfies (11), for selected values of exponent. In all cases, we found that the maximum relative errors attained by comparing the approximate angular frequency values computed from (14) with the exact ones never exceed for the interval . The amount of error attained by our enhanced cubication approach is remarkably lower than those previously found in the literature. See, for instance, [10–12] and references cited therein. This confirms the usefulness of our enhanced cubication approach to obtain approximate solutions of nonlinear oscillators with generalized power-form elastic terms. For illustrative purposes, Figure 2 shows the corresponding amplitude-time curve by assuming that , , , and . There, the solid line represents the numerical integration solution of (19) obtained by using the fourth-order Runge-Kutta algorithm provided by the Mathematica symbolic package, while the red dashed lines represent the enhanced cubication approximate solution. As we can see from Figure 2, both solutions agree well. In fact, the relative error between both solutions is less than .

Table 1: Fitting parameter for selected values of the exponent .

Exponent value
Fitting parameter value

α
σ

1/3
1.0481977

2/3
1.045645

4/3
1.041253

5/3
1.039345

7/5
1.040857

0.99608975

Figure 2: Amplitude-time response curve by assuming that , , , and . Here, the solid lines represents the numerical integration solution of (19), while the dashed lines represent the enhanced cubication approximate solution.

We next develop the approximate solution of a most general restoring force-type that involves terms that describing some of the most common conservative nonlinear oscillators found in the literature.

4. A General Restoring Force Term

Here, we use the enhanced cubication approach to develop the approximate solution of a general oscillator described by (1) with a restoring force given by

where , , , , and are system parameter values. If we introduce the transformation , then the general oscillator differential equation becomes

with initial conditions and . We next use Chebyshev polynomials to find the equivalent representation of as

where

Here,

is the Gauss hypergeometric function. Thus, the approximate solution to (23) is obtained by using (13), (14), (15), and (16).

To assess the accuracy of our enhanced cubication solution, let us consider first the case of a nonlinear dynamical oscillator with an irrational restoring force for which , , , , and [13]. Figure 3 shows that the maximum error attained by comparing our approximate angular frequency values with the exact ones is less than 0.05%. This error is bigger than the error value of 0.04% computed by Lai and Xiang in [13] by using the first two order approximations of the generalized Senator Bapat perturbation technique. Both solutions as well as the numerical integration of (23) are plotted in Figure 4 by considering that . As we can see from Figure 4, all solutions have the same qualitative and quantitative behaviors. Here, the solid black line represents the numerical integration solution, the red dashed line our proposed solution, while the black dots describe the solution derived by Lai and Xiang [13]. For convenience, we have plotted in Figure 5 the frequency ratio versus the amplitude . We may notice that both approximate angular frequency values are close to the exact one. However, we found that Senator Bapat solution fail for small values of , that is, for . Of course, we can improve the percentage error value computed from our enhanced cubication method by considering additional terms in the Chebyshev polynomial expansion of the restoring force term. In the present paper, we shall not consider these additional Chebyshev polynomial expansion terms.

Figure 3: Prediction of the error attained by computing the elastic wire angular frequency oscillator values by the enhanced cubication method. Notice that the maximum relative error between the approximate and the exact angular frequency values does not exceed of 0.05%.

Figure 4: Amplitude-time response curves of a nonlinear oscillator with an irrational restoring force for which , , , , and with a fitting parameter value of . Here, the red dashed line, the black dots, and the green solid line represent the enhance cubication solution, Senator Bapat perturbation solution, and the numerical integration solution of (23), respectively.

Figure 5: Comparison of the frequency ratios versus the amplitude . Here, .

We next study the applicability of the enhanced cubication method by deriving the approximate solution of the finite extensibility nonlinear oscillator for which , , , , and [14, 15]. As we can see from Figure 6, the maximum relative error between the approximate angular frequency value and the exact one is not bigger than % at with . This value is lower than the error value of for reported in [14, 15]. Figure 7 shows the amplitude-time response curves with an initial oscillation amplitude value of . Here, we compare our cubication solution (red dashed line) with respect to the approximate Beléndez et al. solution (black dots) and with the numerical integration solution of (23) (black solid line). As one can see from Figure 7, our derived solution is closer to the numerical integration one. Furthermore, we have computed the approximate angular frequency values of our derived solution for several values of and compared these with respect to the exact values obtained from numerical integrations and with those values computed from the approximate solution derived by Beléndez and coworkers in [15]. Table 2 shows that in general, our enhanced angular frequency value is closer to the exact numerical one even for larger values of . Finally, we have plotted the angular frequency ratio versus the amplitude . It is clear from Figure 8 that both approximate solutions follow well the exact numerical value, however, Beléndez et al. solution shows higher deviation on . This confirms the applicability of our enhanced-cubication approach to obtain accurate solutions of strongly nonlinear oscillators.

Table 2: Comparison of the exact and the approximate angular frequency values for a finite extensibility nonlinear oscillator.

0.5
1.1113
0.9988
1.0000

0.75
1.3369
0.9922
1.0004

0.85
1.5412
0.9848
0.9991

0.95
2.0355
0.9723
0.9707

0.98
2.4600
0.9776
0.9098

0.985
2.5901
0.9836
0.8861

0.9953
3.0856
1.0345
0.7875

Figure 6: Finite extensibility nonlinear oscillator for which , , , , and . The maximum relative error value of % between the approximate angular frequency value and the exact one occurs at the oscillation amplitude of with a fitting parameter value of .

Figure 7: Amplitude-time response curves of the finite extensibility nonlinear oscillator for which , , , , , and with a fitting parameter value of . Here, the red dashed line, the black dots, and the black solid line represent the enhance cubication solution, Beléndez et al. approximate solution, and the numerical integration solution of (23), respectively.

Figure 8: Comparison of the frequency ratios versus the amplitude .

Of course, other equivalent restoring force representations can be studied by assuming different parameters values of , , , , and in (22).

5. A Singular Oscillator

As a final example, we now focus our attention on deriving the approximate solution of the singular oscillator that describes the motion of injected electrons in plasma physics. In this case, the corresponding equation of motion is given as

where is a positive entire number that usually takes the value of one [8, 16]. By using (6)–(8), we can write the approximate Chebyshev equivalent representation form of the restoring force as

It was found in [8, 16, 17] that the exact angular frequency value of (27) when is given as

To find the approximate angular frequency value by using the enhanced cubication approach, we now substitute (28) into (14) and (16), to get that

where is a fitting parameter that satisfies (11). In this particular case, we note that can be found by equating (29) with (30) since this provides an equation that depends only on . By numerically solving the resulting equation, we found two real roots for the fitting parameter . These real roots have the values of and . In this way, we have found two possible approximate solutions to (27) since is equal to (29), that is, . Figure 9 shows the amplitude-time plots obtained by using the values of , , and the solution of (13)

with initial conditions , and . Notice from Figure 9(a) that when we use the root value of , the estimated amplitude-time curve follows surprisingly well the exact numerical solution curve. In this case, the modulus of the elliptic function, , has the value of , where . Also, we have plotted in Figure 9(a) the approximate solution of (27) derived in [8] for which its predicted numerical solution has a dispersion value of when compared to the numerical integration one. As we can see from Figure 9(a), these solutions are almost the same. In Figure 9(a), the solid line represents the derived enhanced-cubication solution, the dashed line represents the numerical integration solution, while the black dots are the estimated prediction values obtained by using the approximate solution derived in [8]. We next use the parameter value of and plot the amplitude-time curves by using (32). In this case, although our enhanced-cubication solution provides the same angular frequency value , the amplitude-time response curve showed in Figure 9(b) differs from the curve obtained by numerical integration. Therefore, we must conclude that the value of is not a valid solution to (27). Furthermore, one can numerically prove that this solution does not satisfy (27) at all times . Thus, we can conclude that only the fitting parameter value of provides a solution that predicts well the trajectory of the singular oscillator. Finally, we must mention that during the derivation of the equivalent restoring force representation of (27), Chebyshev polynomial expansion provides a poor fit to this force term. Due to this finding, we ran several computer simulations by considering different initial condition values. Then, we plotted the amplitude-time response curve of our derived solution and compared this with the exact numerical one. We found in all cases that for , our enhanced solution provided the same qualitative and quantitative behaviors than those of the numerical integration solutions of (27).

(a)

(b)

Figure 9: Amplitude-time response curves of a singular oscillator with initial conditions ,   (a) , and (b) . Here, the solid line represent the derived enhanced-cubication solution, the dashed lines represents the numerical integration solution, while the black dots represent the approximate solution derived in [8] by using Jacobi elliptic functions.

6. Conclusions

In this paper, we have used the enhanced cubication procedure introduced by Elías-Zúñiga et al. in [1] to develop approximate solutions of some strongly nonlinear oscillators that described the dynamical behavior response of several physical and engineering systems. During the solution processes, we used the enhanced cubication approach to replace the original equation of motion by the homogeneous Duffing equation. This transformation is achieved by computing the Chebyshev polynomial representation of the restoring force terms. Then, we used the homogenous Duffing equation to find the angular frequency expression since it is well known that the Duffing equation has an exact solution that depends on Jacobi elliptic functions. Based on these expressions, we have assessed the accuracy of the angular frequency equations for the complete range of oscillations amplitudes. For instance, we have examined the general pendulum response, given by (4), at the interval oscillation values of . In this case, we have found that the maximum relative error does not exceed 0.2%. Then, we used our enhanced cubication method to derive the approximate solution of a generalized power-form elastic term oscillator and found, in all cases, that the maximum error attained is less than which is remarkably lower when compared to other approximate solutions [10–12]. However, when we examined the nonlinear oscillator described by (23) for which , , , , and , the percentage error value found by the enhanced cubication method is bigger than that obtained by using the generalized Senator Bapat perturbation technique [13]. Of course, further improvement could be achieved in our enhanced perturbation method if we consider more terms on the Chebyshev polynomial expansion (6).

We next used the enhanced cubication approach and derived the approximate solution of the finite extensibility nonlinear oscillator and found a maximum error of about % for which represents an improvement of the error values estimated by using other perturbation techniques [14, 15]. Finally, we examined the dynamical response of a singular oscillator (27) and found that our enhanced cubication technique provides two approximate solutions by using the same initial conditions. However, the solution related to the fitting parameter value of provides an amplitude-time curve whose trajectory differs from the numerical one, and then it was discarded as a possible solution to (27).

Based on the previous results, we can conclude that our proposed enhanced-cubication method provides good approximate solutions to strongly nonlinear oscillators even for irrational restoring forces.

In a forthcoming paper, we shall introduce some modifications to our enhanced cubication method to account for linear or nonlinear damping effects of strongly nonlinear oscillators.

Acknowledgment

This work was partially funded by the Tecnológico de Monterrey, Campus Monterrey through the Research Chairs in Nanotechnology and Intelligent Machines.

References

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.

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.

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.

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, M. L. Alvarez, J. Francés et al., “Analytical approximate solutions for the cubic-quintic Duffing oscillator in terms of elementary functions,” Journal of Applied Mathematics, vol. 2012, Article ID 286290, 16 pages, 2012.

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.

G. Radhakrishnan, B. N. Rao, and M. S. Sarma, “On the uniqueness of angular frequency using harmonic balance from the equation of motion and the energy relation,” Journal of Sound and Vibration, vol. 200, no. 3, pp. 367–370, 1997.

A. Elías-Zúñiga, C. A. Rodríguez, and O. Martínez-Romero, “Application of the elliptic balance method to a nonlinear singular oscillator,” Applied Mathematics and Computation, vol. 218, no. 22, pp. 11112–11117, 2012.

V. N. Pilipchuk, “Analytical study of vibrating systems with strong non-linearities by employing saw-tooth time transformations,” Journal of Sound and Vibration, vol. 192, no. 1, pp. 43–64, 1996.

A. Yildirim, Z. Saadatnia, and H. Askari, “Application of the Hamiltonian approach to nonlinear oscillators with rational and irrational elastic terms,” Mathematical and Computer Modelling, vol. 54, no. 1-2, pp. 697–703, 2011.

L. Cveticanin, “Oscillator with fraction order restoring force,” Journal of Sound and Vibration, vol. 320, pp. 1064–1077, 2009.

A. Beléndez, C. Pascual, S. Gallego, M. Ortuño, and C. Neipp, “Application of a modified He's homotopy perturbation method to obtain higher-order approximations of an force nonlinear oscillator,” Physics Letters A, vol. 371, no. 5-6, pp. 421–426, 2007.

S. K. Lai and Y. Xiang, “Application of a generalized Senator-Bapat perturbation technique to nonlinear dynamical systems with an irrational restoring force,” Computers & Mathematics with Applications, vol. 60, no. 7, pp. 2078–2086, 2010.

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.

J. I. Ramos, “Generalized decomposition methods for singular oscillators,” Chaos, Solitons & Fractals, vol. 42, no. 2, pp. 1149–1155, 2009.

A. Beléndez, D. I. Méndez, T. Beléndez, A. Hernández, and M. L. Álvarez, “Harmonic balance approaches to the nonlinear oscillators in which the restoring force is inversely proportional to the dependent variable,” Journal of Sound and Vibration, vol. 314, pp. 775–782, 2008.