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 amplitude-time 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 double-well 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 closed-form 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, [710] 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 [1115]: 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 cubic-like 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 cubic-quintic Duffing equation, the damped general pendulum equation of motion, the damped rational-form elastic term oscillator, and the nonlinear damped cubic term oscillator.

2. The Damped Cubic-Quintic 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 hinged-hinged 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 radio-frequency 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 cubic-like 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 root-mean-square 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 re-write (18) as follows: We next use our proposed approach and write (19) in its equivalent cubic-like form where Here, and are the first and second-order 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 cubic-like 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 two-degree-of-freedom 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 single-degree-of-freedom 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 cubic-like 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 amplitude-time 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 cubic-quintic 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.