Recent Advances on Methods and Applications of Nonlinear Differential Equations
View this Special IssueResearch Article  Open Access
Equivalent Mathematical Representation of SecondOrder Damped, Driven Nonlinear Oscillators
Abstract
The aim of this paper focuses on applying a nonlinearization method to transform forced, damped nonlinear equations of motion of oscillatory systems into the wellknown forced, damped Duffing equation. The accuracy obtained from the derived equivalent equations of motion is evaluated by studying the amplitudetime, the phase portraits, and the continuous wavelet transform diagrams of the cubicquintic Duffing equation, the generalized pendulum equation, the powerform elastic term oscillator, the Duffing equation with linear and cubic damped terms, and the pendulum equation with a cubic damped term.
1. Introduction
Here, in this paper, we illustrate how the nonlinearization approach can be used to obtain equivalent equations of motion of forced, damped nonlinear oscillators of the form where is the initial oscillation amplitude, and are the system conservative and nonconservative restoring forces, respectively, is a damping parameter, and is a periodic external force with driving frequency . In accordance with the nonlinear transformation approach [1], we first write the conservative force terms as a polynomial expression by using the Chebyshev polynomials of the first kind [2–4]: where and are the Chebyshev polynomials of the first kind. We can see that the usage of (3) could transform (2) into a fifth or higher order polynomial expression. In the case for which a fifthorder Chebyshev polynomial is used, the conservative force in (2) becomes where , and will be defined later on. Therefore, the equivalent representation form of (1) is given as By following the nonlinearization method, we now find the equivalent representation form of (1) as a function of a cubiclike polynomial equation. This procedure leads to where , , and are determined from in which Notice that in our proposed procedure we are assuming that the magnitude of the external force and its driving frequency remain constants during the transformation process. Thus, (1) can be written in equivalent form as We will next explore the applicability of our proposed approach and derive the equivalent representation form of some forced, damped nonlinear systems.
2. The Forced, Damped CubicQuintic Oscillator
The equation of motion that describes the dynamical response of the forced, damped cubicquintic oscillator is given as where denotes the displacement of the system, is the damping coefficient, , , and are system constant parameters, is the magnitude of the external force, and is the driving frequency [6, 7]. We next use and write (11) as where , , , and . By following our proposed nonlinear method and by using (4) and (7)–(9), we obtain the equivalent representation form of (12) as where , , and can be determined from the following equations: Here, and are fitting parameters that satisfy (7)(8). To examine the accuracy of (13), we next compare its solution with the one obtained from (12) by using the fourthorder RungeKutta numerical integration method. Let us consider the parameter values of , , , , , , and . In this case, the parameter values assigned to , , and provide a triplewell potential to the cubicquintic Duffing oscillator that can have up to four resonance frequencies [8]. Figure 1 illustrates the comparison between the amplitudetime response curves of (12) and (13) obtained from their corresponding numerical integration solutions. As one can see form Figure 1, both solutions are almost the same. In fact, the computed rootmeansquare error (RMSE) value does not exceed on with , , , , and . The accuracy of the numerical simulations is surprisingly good if we consider that the potential of a cubic Duffing oscillator cannot have triplewell form.
As a second example, let us consider the parameter values of , , , , , , and . Figure 2 shows the amplitudetime response curves and the corresponding phase portraits, as well as the Morlet continuos wavelet transforms (CWT) obtained from the numerical integration solutions of (12) and (13). Here, the values of and were computed from (15) which provides good agreement between (12) and (13). Notice that the numerical integration solutions of (12) and (13) are almost the same. In this particular problem, the Morlet CWT was used to extract system dynamics effects such as the one shown at the system transient motion in which the transient frequency has strong influence on the system dynamic behavior. In fact on the time interval , the transient frequency dominates the system motion. When , the system oscillates at the driving frequency . Besides, we have computed the RMSE value between both numerical solutions and found that it has the value of . Here, , , and . Of course, we can consider other parameter values, as those shown in Table 1, to describe the dynamic response of (12) by using (13). Therefore, we can conclude that our nonlinear method leads to the derivation of an equivalent equation of motion that follows well the qualitative and quantitative numerical response of the original equation (12).

We next determine the equivalent representation form of the forced, damped general pendulum equation.
3. The Forced, Damped General Pendulum Equation
We now proceed to derive the equivalent representation form of the forced, damped pendulum equation where and represent system constant parameter values, is the damping coefficient, is the driving frequency, and is the external force magnitude [9]. If we introduce the transformation , then (16) can be rewritten as with . By applying our proposed transformation method to (17), we obtain the following expression: where and , , and are given by (14) and (15). Here, and are the first and second order Bessel functions of the first kind. To illustrate the degree of accuracy attained by our derived solution (18), let us consider the system parameter values of , , , , and with or and . One can notice from Figure 3 that the numerical integration solutions of (17) and (18) are almost the same. In this case, ,, , , , and , and the values of , and were fitted by using (15), since these expressions provide the best predictions with a RMSE value of . The same degree of accuracy was found by considering different system parameter values, as those illustrated in Table 1 in which the RMSE values are close to zero.
To further assess the applicability of our nonlinear cubication approach, we next derive the equivalent representation form of a forced, damped oscillator with a powerform elastic term.
4. A Generalized Forced, Damped PowerForm Elastic Term Oscillator
The equation of motion of this oscillator is given as where and are constant parameters and can take any nonnegative real value, such as odd, even, rational, or irrational, that is, [10]. As usual, let us use the following coordinate transformation and write (20) as, We next use Chebyshev polynomial expansion to write the restoring forces as a nonic polynomial expression where Notice that in (23)–(27) the terms represent the Euler gamma function. It is important to point out that in this particular problem we have used five Chebyshev expansion coefficient terms that provide, for the system restoring force, an equivalent representation form that is based on a ninthorder polynomial expression. This example illustrates the applicability of our procedure in using more than three terms in (4). We next follow our solution procedure and find, by using (7) and (8), that where To assess the accuracy of our derived equivalent representation form (28) of (21), we shall consider the following data values: , , , , , and with a driving frequency value of [11]. Figure 4 illustrates the amplitudetime response curves obtained by numerically integrating (21) and (28). As we can see from Figure 4, the numerical integration of (28) follows closely the amplitudetime response curve obtained from (21). In this case, the RMSE value of is obtained by using equations (30)–(32). Here, the red solid and black dashed lines represent, respectively, the numerical integration solution of (21) and (28). The computed parameter values are , , , , , , , , , and . Also, Figure 5 provides a comparison of the numerical solutions of (21) and (28) with respect to the approximate general solution of (28) derived by using Jacobi elliptic functions [5]. One can notice from Figure 5 that all solutions are almost the same. Therefore, we can conclude that our derived equivalent representation form (28) describes well the qualitative and quantitative behavior of (21). The amplitudefrequency response curve of (28) can be obtained by using, for instance, the approximate solutions developed in [5, 11]. As a second case, we now use our equivalent representation form (28) and consider the following parameter values of , , , , , , and in (21) and compute the corresponding amplitudetime response curve. We can see from Figure 6 that the amplitudetime curve obtained from (28) follows well the curve obtained from (21). In this case, the RMSE value is about for which the parameter values are , , , , , , , , , and . For illustrative purposes, we show in Table 2 some values of the exponent with their fitting parameter values of and that can be used to study the dynamical behavior of some nonlinear oscillator with a rational or irrational power restoring forces.

(a)
(b)
We next develop the equivalent representation form of the Duffing equation with linear and cubic damped terms.
5. The Forced Duffing Equation with Linear and Cubic Damped Terms
We now explore the applicability of our method to derive the equivalent representation form of the following equation of motion which has a linear damped term, , and a cubic one, [12, 13]. Let ; then, (33) can be written as where , , and .
Since (34) has a damped nonlinear term of the cubictype, we need to modify our nonlinear method to take into account its effects on the solution response of (34). Therefore, we now assume that (7) and (8) can be rewritten as which yield the equivalent representation form of (34) as where Before we evaluate the accuracy achieved by our derived expression (37), we first recall that Trueba and coworkers in [13], by using Melnikov analysis, found an equivalent equation of motion for (34) given as where is defined as In what follows, we will use (37) and (39) to compare their numerical predictions with those provided by (34). First, let us consider the parameter values of , , , , , and with , and and use our derived expressions to compute the values of , , , and which are given as , , , , and , respectively. Figure 7 shows a comparison of the amplitudetime curves, the phase portrait plots, and the Morlet CWT diagram obtained from the numerical integrations of (34), (37), and (39). Notice from Figure 7, that our equivalent equation of motion (37) closely follows the numerical integration curve of (34). Here, the RMSE value is close to , while the numerical predictions obtained from (39) show some discrepancies in the amplitudetime curve at the time interval of . In this solution, the computed RMSE value is .
(a)
(b)
(c)
As a second case, we now explore the accuracy of our equivalent representation form (37) by assuming that , , , , , and with and . As we can see from Figure 8, the numerical integration solutions of (37) and (39) agree well with the solution of (34). Furthermore, the Morlet CWT of (39) shown in Figure 8, exhibits the subharmonic effects that the original system (34) experiences at the frequencies values of , , and which correspond to the Morlet continuous wavelet transform scale values of , , and , respectively. In this case, the estimated RMSE values are for (37) and for (39). The values of , , , and are found to be , , , , and . Notice that in both examples the values of , , and remain unchanged. In Figures 7 and 8, the black and the red dashed lines describe the amplitudetime and phase portrait curves of (34) and (39), respectively.
(a)
(b)
(c)
As a final example, we now derive the equivalent equation of the forced pendulum equation with a cubic damped term where is the initial oscillation amplitude, , , with . We first use Chebyshev polynomial expansion and then (35) and (36) to get that where and , , and are determined from By using Melnikov analysis, Trueba and coworkers [13] developed the equivalent representation form of (41) which is given as We next consider the parameter values of , , , , , , and , and plot the numerical integration solutions of (41), (42), and (45). The corresponding amplitudetime, phase portrait, and Morlet CWT plots are shown in Figure 9 for which the computed parameter values are , , , , , , , with , and . We can see from Figure 9 that our solution closely follows the numerical simulations of (41). In this case, the computed RMSE values from (42) and (45) are and , respectively. In Figure 9, the black, the purple, and the red solid lines represent, respectively, the numerical integration solutions of (41), (42), and (45).
(a)
(b)
(c)
This confirms the usefulness of our proposed nonlinear method to obtain equivalent equations of motion of nonlinear oscillators.
6. Conclusions
In this paper, we have obtained the equivalent representation form of some driven, damped nonlinear oscillators by using a nonlinearization approach. We have found that in all cases, the numerical predictions obtained from the corresponding equivalent representation form of the cubicquintic, the general pendulum, the powerform elastic term, and the cubic damped nonlinear oscillators describe well the qualitative and quantitative behavior of their original equations of motion.
During the solution processes of the forced Duffing equation with a cubic damped term, we have found that our equivalent solution (37) provides numerical estimates that are similar to those obtained from (39) which was derived from Melnikov analysis in [13]. Besides, the numerical predictions of (37) are closer to the numerical integration values of (34) than those obtained from (39). In this oscillator, we have found that the values of , , and could remain unchanged even at larger values of , , , , and . To further evaluate the accuracy of our proposed method, we have developed the solution of a pendulum equation with a cubic damped term and compared its numerical estimated values with those obtained from the original equation of motion and with respect to those of its equivalent representation form derived by Melnikov analysis [13]. We found that our derived equation (42) describes well the numerical estimated values of (41). Moreover, and based on the numerical simulations performed on the last example of this work, it is clear that the derived equivalent equation of motion, when compared to other solutions such as the one derived by using the wellestablished Melnikov analysis, exhibits good accuracy for a wide range of system parameter values [13]. Of course, one must be careful when using the equivalent transformation forms previously derived, since their degree of accuracy depends not only on the system parameter values but also on the method used to replace the corresponding restoring forces for equivalent ones of the cubic type.
Finally, it is evident that our proposed nonlinearization method can be used to derive equivalent representation forms of other nonlinear oscillators such as the ones examined by the authors in [14, 15] and references cited therein, in which rational or irrational restoring forces, as well as damping terms, are used to model the dynamics behavior of common problems that arise in the physical sciences and engineering fields.
Acknowledgments
This work was funded by the Tecnológico de Monterrey, Campus Monterrey —Campus Monterrey, through the Research Chair in Nanomaterials for Medical Devices and Research Chair in Intelligent Machines. Additional support was provided by the European Union Seventh Framework Programme (FP7PEOPLE2009) under the grant agreement IRSES no. 247476 and from Consejo Nacional de Ciencia y Tecnología (Conacyt), México.
References
 A. ElíasZúñiga, O. MartínezRomero, and R. K. CórdobaDíaz, “Approximate solution for the Duffingharmonic oscillator by the enhanced cubication method,” Mathematical Problems in Engineering, Article ID 618750, 12 pages, 2012. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 A. Beléndez, M. L. Lvarez, E. Fernndez, and I. Pascual, “Cubication of conservative nonlinear oscillators,” European Journal of Physics, vol. 30, no. 5, pp. 973–981, 2009. View at: Publisher Site  Google Scholar
 A. Beléndez, D. I. Méndez, E. Fernández, S. Marini, and I. Pascual, “An explicit approximate solution to the Duffingharmonic oscillator by a cubication method,” Physics Letters, Section A: General, Atomic and Solid State Physics, vol. 373, no. 32, pp. 2805–2809, 2009. View at: Publisher Site  Google Scholar
 A. Beléndez, G. Bernabeu, J. Francés, D. I. Méndez, and S. Marini, “An accurate closedform approximate solution for the quintic Duffing oscillator equation,” Mathematical and Computer Modelling, vol. 52, no. 34, pp. 637–641, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. ElíasZúñiga, “A general solution of the Duffing equation,” Nonlinear Dynamics, vol. 45, no. 34, pp. 227–235, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. S. Siewe and U. H. Hegazy, “Homoclinic bifurcation and chaos control in MEMS resonators,” Applied Mathematical Modelling. Simulation and Computation for Engineering and Environmental Systems, vol. 35, no. 12, pp. 5533–5552, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Taylan, “The effect of nonlinear damping and restoring in ship rolling,” Ocean Engineering, vol. 27, no. 9, pp. 921–932, 2000. View at: Publisher Site  Google Scholar
 S. Jeyakumari, V. Chinnathambi, S. Rajasekar, and M. A. F. Sanjuan, “Analysis of vibrational resonance in a quintic oscillator,” Chaos, vol. 19, no. 4, Article ID 043128, 2009. View at: Publisher Site  Google Scholar
 E. I. Butikov, “Extraordinary oscillations of an ordinary forced pendulum,” European Journal of Physics, vol. 29, no. 2, pp. 215–233, 2008. View at: Publisher Site  Google Scholar
 V. N. Pilipchuk, “Analytical study of vibrating systems with strong nonlinearities by employing sawtooth time transformations,” Journal of Sound and Vibration, vol. 192, no. 1, pp. 43–64, 1996. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 I. Kovacic, “Forced vibrations of oscillators with a purely nonlinear powerform restoring force,” Journal of Sound and Vibration, vol. 330, no. 17, pp. 4313–4327, 2011. View at: Publisher Site  Google Scholar
 A. K. Mallik, V. Kher, M. Puri, and H. Hatwal, “On the modelling of nonlinear elastomeric vibration isolators,” Journal of Sound and Vibration, vol. 219, no. 2, pp. 239–253, 1999. View at: Google Scholar
 J. L. Trueba, J. Rams, and M. A. F. Sanjuán, “Analytical estimates of the effect of nonlinear damping in some nonlinear oscillators,” International Journal of Bifurcation and Chaos in Applied Sciences and Engineering, vol. 10, no. 9, pp. 2257–2267, 2000. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Alex ElíasZúñiga and Oscar MartínezRomero, “Accurate solutions of conservative nonlinear oscillators by the enhanced cubication method,” Mathematical Problems in Engineering, vol. 2013, Article ID 842423, 9 pages, 2013. View at: Publisher Site  Google Scholar
 Alex ElíasZúñiga and Oscar MartínezRomero, “Investigation of the equivalent representation form of strongly damped nonlinear oscillators by a nonlinear transformation approach,” Journal of Applied Mathematics, vol. 2013, Article ID 245092, 7 pages, 2013. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2013 Alex ElíasZúñiga and Oscar MartínezRomero. 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.