Energy Method to Obtain Approximate Solutions of Strongly Nonlinear Oscillators

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.


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 cubicor 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][4][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 cubicquintic 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%.

A Nonlinearization Energy Procedure
Conservative nonlinear oscillators are commonly modelled by equations of the form where () is the conservative system restoring force and  10 is the initial oscillation amplitude.Integration of (1) with Mathematical Problems in Engineering respect to the displacement  leads to its energy equation of motion Here, ẋ and () represent the velocity and the system potential energy terms, respectively, and  0 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][4][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 ,   ,  = 1, 2, . . ., , depend on the Chebyshev coefficient terms.Therefore, the energy equation ( 2) can be equivalently represented by the following expression: where  0 ≡  0 + .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  0 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).

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  = / 10 , then (16) can be written as Here  = ,  =  2 10 , and  =  4 10 .To apply our nonlinear energy approach, we first integrate (17) with respect to .This step provides the following energy equation: where the integration constant  0 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  = √ −1.The integration of (23) provides the following time expression: Here  represents the incomplete elliptic integral of the first kind with modulus ( 3 − 2 )/( 1 − 2 ).Thus, the approximate angular frequency value of the cubic-quintic oscillator is given as where   = 4  on  ∈ [0, 1].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  1 ,  1 , and Δ 1 can be determined by using ( 14) and (15) which yields where the value of  1 must satisfy (14).Now the solution of yields where Thus, the approximate angular frequency expression is given as   = 2/  , where   = 4  .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  = 1,  = 10, and  = 1 with  = 1.0854.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  1 = 1.13731, 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 fourthorder 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).

A Finite Extensibility Nonlinear Oscillator
This oscillator is used to describe chain dynamics of polymerlike materials subjected to high extensibility values [12].Its equation of motion has the form where  is the initial oscillation amplitude with  ∈ (0, 1).If we introduce the transformation  = /, then (33) can be written as It is well known that when  → 1, 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][8][9][10][11][12][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 6  , 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  1 =  1 ,  1 = 2 1 , and Δ 1 =  1 /2. 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,  = 0.9999999, 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][8][9][10][11][12][13]. Furthermore, the maximum error values of both curves do not exceed 3.2%.Here the fitting parameter  has the value of 1.01195.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 0 ≤  ≤ |1|. 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,  = 0.999.

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  = 0.9999999, 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.

Figure 3 :
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 0 <  ≤ 0.9999999.

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

Figure 5 :
Figure5: 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  = 0.999.