Research Article  Open Access
Alex ElíasZúñiga, 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. https://doi.org/10.1155/2013/245092
Investigation of the Equivalent Representation Form of Strongly Damped Nonlinear Oscillators by a Nonlinear Transformation Approach
Abstract
We use a nonlinear transformation method to develop equivalent equations of motion of nonlinear homogeneous oscillatory systems with linear and nonlinear odd damping terms. We illustrate the applicability of our approach by using the equations of motion that arise in many engineering problems and compare their amplitudetime curves with those obtained by the numerical integration solutions of the original equations of motion.
1. Introduction
The dynamics response of some systems can be more precisely described when nonlinear damping terms are used to model their dynamics behaviors. For instance, the elastomeric vibration isolators [1], the motion of a rolling ship subjected to the synchronous beam waves [2], the backlash and impact phenomena [3], and the micromechanical oscillators [4], to say a few, have been modeled by considering nonlinear damping terms. In fact, the dynamic behavior of doublewell oscillators in which a nonlinear damping term with a fractional exponent covers the gaps between viscous, dry friction, and turbulent damping phenomena has been used by Litak et al. in [5] to study, by using the Melnikov criterion, the system global homoclinic bifurcation and its transition to chaos. It is evident from the previously mentioned works and references cited therein that the global system dynamics behavior can be accurately described if one is able to identify the order of the nonlinear stiffness and the damping effects that agree with the experimental observations [6]. Of course, the influence of the nonlinear damping terms on the resulting equations of motion increases the difficulty of finding their closedform solutions.
In this paper, a nonlinear transformation of the damped nonlinear equation is proposed to obtain its equivalent damped Duffing’s equation of motion. Here, we assume that is the nonconservative system restoring force which could have rational or irrational conservative force terms as well as linear or nonlinear damping terms, and we assume that is the initial amplitude. The main motivation to find a nonlinear transformation form of (1) is based on the fact that exact or approximate solutions of damped nonlinear oscillators of the Duffing type can be found in the literature. See, for instance, [7–10] and references cited therein. Therefore, if we can transform (1) into the damped Duffing equation, one could find its dynamical response in an easier way.
In order to achieve such transformation, we first assume that the conservative terms of can be written in equivalent forms by using, for instance, the Chebyshev polynomials of the first kind [11–15]: where Since the first three Chebyshev polynomials of the first kind are given by thus, the equivalent conservative restoring force can be written as One must notice that the coefficients , , and depend on the amplitude of oscillation, , and the Chebyshev coefficient terms. Therefore, we assume that the nonlinear differential equation (1) can be replaced by an equivalent equation of the form where represents the system damping terms and is the damping coefficient. We next use a cubication transformation to write the restoring force term of (6) in the form Here, , , and can be found by replacing the terms by the cubic polynomial that satisfies
Thus, the cubiclike equivalent representation form of (1) is given as
We next examine the application of our proposed nonlinear transformation approach to obtain approximate solutions of the damped oscillatory systems such as the damped cubicquintic Duffing equation, the damped general pendulum equation of motion, the damped rationalform elastic term oscillator, and the nonlinear damped cubic term oscillator.
2. The Damped CubicQuintic Duffing Equation
It is well known that this equation is used to describe the dynamical behavior of beams resting on an elastic substrate [16], the nonlinear transverse vibration of a hingedhinged flexible beam subjected to constant excitation at its free end [17], the biological rhythmic movements [18], the propagation of electromagnetic pulses in media with saturable nonlinearity [19], the intermodulation distortion in radiofrequency microelectromechanical systems (MEMS) capacitors [20], the motion of a rolling ship subjected to synchronous beam waves [2], and so forth. This equation has the form where denotes the displacement of the system, is the damping coefficient, and , , and are the system constant parameters. We next introduce the following change of the variable which transforms (13) into an equation of the form Here, , , and . In accordance with our proposed nonlinear transformation approach, we first replace the restoring force by an equivalent cubiclike polynomial expression by using (8), (9), and (10). This provides the following restoring force expression: Thus, (13) can be written in an equivalent form as follows: where Notice that and are the parameter values that must satisfy (9) and (10). Figure 1 illustrates the numerical integration solutions of (12) and (15) by considering the parameter values , , , and with the initial conditions and . To obtain the smallest rootmeansquare error (RMSE) between both numerical solutions, we have fitted the values and in (17). Thus, the parameter values of the damped cubic Duffing equation (15) become , , and with an RMSE value of for the time interval shown in Figure 1. Also, we have computed the RMSE values by considering different system initial conditions as showed in Table 1. As we may see from Table 1, the RMSE values are not bigger than . As a second example, we now assume the following system parameter values , , , and . The RMSE values obtained from the comparison of the numerical integration solutions of (14) and (15) are shown in Table 2. Notice that the maximum RMSE value is now when .


We will next examine the applicability of our proposed approach by deriving the equivalent expression of the damped general pendulum equation.
3. A Damped Pendulum Equation
We now use our nonlinear transformation approach to find the equivalent equation of motion of the damped pendulum equation [21] with the initial conditions and . Here, and represent the system constant parameter values. First, let us introduce the transformation and rewrite (18) as follows: We next use our proposed approach and write (19) in its equivalent cubiclike form where Here, and are the first and secondorder Bessel functions of the first kind, respectively. Once again, and are the fitting parameters that must satisfy (8) and (9). We next illustrate in Figure 2 the numerical integration solutions of (19) and (20) by considering the system parameter values , , and . In this particular case, the values of , , and are , , and . The initial angular displacement amplitude is assumed to have the value of . As we may see from Figure 2, both solutions are almost the same. In fact, the maximum RMSE value attained is with the computed parameter values , , , , and .
To further assess the applicability of our proposed approach, we next derive the equivalent cubiclike representation form of a damped oscillator with a mass attached to two stretched elastic springs.
4. The Damped Nonlinear Oscillator with an Irrational Restoring Force
Oscillators with irrational restoring forces are used to model the oscillations of a mass attached to a stretched wire [22] and are also used to study the dynamical response of vibration isolators or vibration absorbers in twodegreeoffreedom systems [23]. In this case, we assume that, besides the nonlinear irrational restoring force due to the stretched wires, we have a damper attached to the mass system, and hence, the nonlinear differential equation of motion for a singledegreeoffreedom system can be written as where , , , and are the system parameter values. If we introduce the transformation , then, (22) becomes with the initial conditions and . We next apply our proposed nonlinear transformation approach to get the equivalent representation form of (23) as follows: where Here, and represent the complete elliptic integrals of the first and second kinds, respectively, with modulus equal to . To assess the accuracy of our equivalent damped equation of motion (24), let us consider the case for which , , , and . Figure 3 shows the numerical integration solutions of (22) and (24) with . As we can see from Figure 3, both solutions are almost the same for most of the time interval shown. The computed parameter values are , , , , , , and . In this case, the RMSE value is equal to . It is easy to show that the RMSE values do not exceed on . It is evident from our numerical results that the usage of irrational restoring forces with damping effects can be equivalently described by the damped Duffing oscillator. This could help us in better understanding the influence of irrational forces on the dynamical responses of the vibrational systems with two or more degrees of freedom.
We next derive the equivalent equation of motion of a Duffing oscillator with linear and cubic damped terms.
5. The Duffing Equation with Linear and Cubic Damped Terms
The equation of motion with linear and cubic damped terms is used to model the dynamical responses of several engineering applications such as the nanomechanical dynamical response of a doubly clamped beam [24], the nonlinear rolling motion of a ship in random beam seas [25], the modeling of nonlinear elastomeric vibration isolators [1], and the generalized damped general pendulum equation [26], among others. In (26), and represent the magnitude of the linear and nonlinear cubic damped terms, respectively. Of course, the exact solution of (26) is unknown, and thus, the numerical methods or perturbation techniques must be used to obtain its approximate solution.
Trueba and coworkers in [26] used the Melnikov analysis to write the forced version of (26) as the damped Duffing equation in which is the forcing term, , and Our aim in this section is to use our proposed approach to develop an equivalent representation form of (26) and compare its theoretical predictions with respect to those of (26) and (27) with .
First, let us assume that the restoring force in (26) can be written in an equivalent form by using the following equations: where the coefficients , , and can be determined by using the expressions (10). Notice that we have introduced a slight modification in the damping coefficient terms of (29) to take into account the influence of the cubic nonlinear damping term, , in the dynamical system response. We then follow our nonlinear transformation method to obtain the equivalent representation form of (26) which is given as To numerically evaluate the accuracy of our proposed equivalent representation form (30), let us consider the following system parameter values , , , and with . As we may see from Figure 4, the numerical integration solution of (30) with the estimated parameter values , , , , and agrees well with the numerical integration solution of (26). In this case, the computed RMSE value does not exceed 0.0426, while the corresponding RMSE value found by using (27) is about . In Figure 4, the solid line and the red and purple dashed lines represent, respectively, the numerical integration solutions of (26), (27), and (30). It is evident from Figure 4 that our developed cubiclike solution exhibits good accuracy when compared with the numerical predictions of (27). To further assess the precision of our equivalent representation form (30), we next consider the system parameter values , , , and with . In this case, we found that , , , , and , with the computed RMSE values and obtained by using (27) and (30), respectively. Figure 5 illustrates the numerical integration amplitudetime response curves of (26), (27), and (30). As expected, our equivalent equation of motion (30) follows well the numerical solution of (26). This confirms the usefulness of our proposed approach to solve nonlinear oscillators with nonlinear damping terms.
6. Conclusions
In this paper, we have used a nonlinear transformation procedure to obtain equivalent equations of motion of nonlinear oscillators with conservative and dissipative restoring forces. Our solution procedure provides equivalent equations whose numerical predictions follow well the numerical integration solutions of the original equations of motion for small or larger damping coefficient values. To show the feasibility of our proposed approach, we have found the equivalent representation forms of four damped nonlinear oscillators. In fact, we have shown that (15) describes well the damped cubicquintic Duffing equation since the computed RMSE values shown in Tables 1 and 2 are not bigger than on even for larger nonlinear parameter values. Furthermore, we found that the dynamical behavior of a damped pendulum is accurately described by its equivalent damped Duffing equation (20). In this case, the RMSE value does not exceed for the system parameter values , , and with . A similar conclusion can be drawn by using the equivalent equation of motion (24) that models a damped nonlinear oscillator with an irrational restoring force term.
Finally, we have shown that the effect of a nonlinear damping dissipative force is equivalent to a linearly damped nonlinear Duffing oscillator with a modified damping coefficient. This is illustrated in the case of the Duffing equation with linear and nonlinear damping terms in which numerical predictions follow well the responses of the original equations of motion. Here, our numerical predictions coincide with those obtained from the Melnikov analysis even at larger nonlinear damping values. It is clear that our proposed approach could be extended not only to include driving forces or dissipative effects with even or fractional nonlinear damping terms but also to obtain equivalent equations of physical system with two or more degrees of freedom.
Acknowledgment
This work was funded by the Tecnológico de Monterrey, Monterrey Campus, through the Research Chair in Nanomaterials for Medical Devices and the Research Chair in Intelligent Machines.
References
 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
 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
 R. S. Barbosa and J. A. T. Machado, “Describing function analysis of systems with impacts and backlash,” Nonlinear Dynamics, vol. 29, no. 1–4, pp. 235–250, 2002. View at: Publisher Site  Google Scholar
 S. Zaitsev, O. Shtempluck, E. Buks, and O. Gottlieb, “Nonlinear damping in a micromechanical oscillator,” Nonlinear Dynamics, vol. 67, no. 1, pp. 859–883, 2012. View at: Publisher Site  Google Scholar
 G. Litak, M. Borowiec, and A. Syta, “Vibration of generalized double well oscillators,” Zeitschrift für Angewandte Mathematik und Mechanik, vol. 87, no. 89, pp. 590–602, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 C. Ho, Z. Lang, and S. A. Billings, “The benefits of nonlinear cubic viscous damping on the force transmissibility of a Duffingtype vibration isolator,” in Proceedings of the UKACC International Conference on Control, Cardiff, UK, September 2012. View at: Google Scholar
 V. K. Chandrasekar, M. Senthilvelan, and M. Lakshmanan, “New aspects of integrability of forcefree Duffingvan der Pol oscillator and related nonlinear systems,” Journal of Physics A, vol. 37, no. 16, pp. 4527–4534, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Z. Feng and Q. Meng, “First integrals for the damped Helmholtz oscillator,” International Journal of Computer Mathematics, vol. 87, no. 12, pp. 2798–2810, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. K. Lai and K. W. Chow, “Exact solutions for oscillators with quadratic damping and mixedparity nonlinearity,” Physica Scripta, vol. 85, no. 4, Article ID 045006, 6 pages, 2012. View at: Publisher Site  Google Scholar
 B. S. Wu and W. P. Sun, “Construction of approximate analytical solutions to strongly nonlinear damped oscillators,” Archive of Applied Mechanics, vol. 81, no. 8, pp. 1017–1030, 2011. View at: Publisher Site  Google Scholar
 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, vol. 2012, Article ID 618750, 12 pages, 2012. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 A. ElíasZúñiga and O. 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
 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 A, 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
 S. Lenci, G. Menditto, and A. M. Tarantino, “Homoclinic and heteroclinic bifurcations in the nonlinear dynamics of a beam resting on an elastic substrate,” International Journal of NonLinear Mechanics, vol. 34, no. 4, pp. 615–632, 1999. View at: Publisher Site  Google Scholar  MathSciNet
 H. M. Sedighi, K. H. Shirazi, and J. Zare, “An analytic solution of transversal oscillation of quintic nonlinear beam with homotopy analysis method,” International Journal of NonLinear Mechanics, vol. 47, pp. 777–784, 2012. View at: Google Scholar
 D. Delignières, D. Nourrit, T. Deschamps, B. Lauriot, and N. Caillou, “Effects of practice and task constraints on stiffness and friction functions in biological movements,” Human Movement Science, vol. 18, no. 6, pp. 769–793, 1999. View at: Publisher Site  Google Scholar
 A. I. Maǐmistov, “Propagation of an ultimately short electromagnetic pulse in a nonlinear medium described by the fifthorder duffing model,” Optics and Spectroscopy, vol. 94, no. 2, pp. 251–257, 2003. View at: Publisher Site  Google Scholar
 M. S. Siewe and U. H. Hegazy, “Homoclinic bifurcation and chaos control in MEMS resonators,” Applied Mathematical Modelling, vol. 35, no. 12, pp. 5533–5552, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 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. View at: Publisher Site  Google Scholar
 S. K. Lai and Y. Xiang, “Application of a generalized SenatorBapat perturbation technique to nonlinear dynamical systems with an irrational restoring force,” Computers & Mathematics with Applications, vol. 60, no. 7, pp. 2078–2086, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Q. Cao and Y. Xiong, “Dynamic analysis of a two degree of freedom system with irrational nonlinearity,” in Proceedings of the 7th European Nonlinear Dynamics Conference, D. Bernardini, G. Rega, and F. Romeo, Eds., Universita di Roma, Rome, Italy, July 2011. View at: Google Scholar
 S. Zaitsev, R. Almong, O. Shtempluck, and E. Buks, “Nonlinear dynamics in nanomechanical oscillators,” in Proceedings of the International Conference on MEMS, NANO and Smart Systems (ICMENS ’05), pp. 387–391, 2005. View at: Google Scholar
 J.Y. Gu, “Nonlinear rolling motion of ship in random beam seas,” Journal of Marine Science and Technology, vol. 12, no. 4, pp. 273–279, 2004. 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
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.