Abstract

The cube-root truly nonlinear oscillator and the inverse cube-root truly nonlinear oscillator are the most meaningful and classical nonlinear ordinary differential equations on behalf of its various applications in science and engineering. Especially, the oscillators are used widely in the study of elastic force, structural dynamics, and elliptic curve cryptography. In this paper, we have applied modified Mickens extended iteration method to solve the cube-root truly nonlinear oscillator, the inverse cube-root truly nonlinear oscillator, and the equation of pendulum. Comparison is made among iteration method, harmonic balance method, He’s amplitude-frequency formulation, He’s homotopy perturbation method, improved harmonic balance method, and homotopy perturbation method. After comparison, we analyze that modified Mickens extended iteration method is more accurate, effective, easy, and straightforward. Also, the comparison of the obtained analytical solutions with the numerical results represented an extraordinary accuracy. The percentage error for the fourth approximate frequency of cube-root truly nonlinear oscillator is 0.006 and the percentage error for the fourth approximate frequency of inverse cube-root truly nonlinear oscillator is 0.12.

1. Introduction

Nonlinear systems are widespread around us. Nonlinear systems are widely involved in science, engineering, medical science, etc. So, research on the nonlinear systems has enriched science, engineering, medical science, etc. Research on nonlinear systems is complex and sensitive because most of the nonlinear systems suddenly change their characteristics due to some small changes of some valid parameters. The development of the theorems of dynamical systems has been derived by the modeling and formulating of nonlinear oscillators. Thus, nonlinear oscillators are one of the most important parts of nonlinear dynamical systems. Recently, many scientists have made significant improvement in finding a new mathematical tool which would be related to nonlinear dynamical systems whose understanding will rely not on analytic techniques but also on numerical and asymptotic methods. They set up many fruitful and potential methods to operate the nonlinear systems. There are many analytical methods to solve nonlinear oscillators such as perturbation [13]; standard and modified Poincare–Linstedt [4]; harmonic balance [57]; multiple scale [8]; homotopy perturbation [914]; modified He’s homotopy perturbation [15]; He’s frequency-amplitude formulation [16, 17]; cubication method [18, 19]; energy balance method [20]; He’s energy balance method [21]; Mickens iteration method [7, 2225]; Hu’s iteration method [26]; Haque’s iteration method [2729]; Haque’s extended iteration method [3032]; variation iteration method [33]; homotopy analysis method [34]; finite element method and Akbari Ganji method [35]; and Taguchi method [36]. Besides, a few numbers of researchers have done work on the cube-root nonlinear oscillator using different methods such as the methods by Beléndez [9], Beléndez et al. [15], Ganji et al. [21], Mickens [25], Zhang [17], and Lim and Wu [6]. Among them, the proposed method is more easy, simple, and also valid for higher order approximation solutions. Perturbation method is the most widely utilized established method in which the nonlinear terms are considered as small. Mickens has developed a technique named harmonic balance (HB) method and further work has been done by Lim, Hu, Wu, and Alam, and so forth, for handling the oscillators in which nonlinear terms are strong. Mickens has developed another update technique named iteration method. The method of iteration is applicable for small as well as strong nonlinearities. Mickens has also developed the iteration method named extended iteration method; sometimes, this is faster than direct or simple iteration method. Further, the method of iteration and extended iteration has been developed by Lim, Hu, and Haque. Ikramul Haque et al. applied the proposed method to obtain the suitable periodic solution of nonlinear oscillators [3032].

The purpose of this article is to exhibit an extended iteration method which gives us a simple path line for obtaining the approximate angular frequencies and corresponding analytic solutions and a significantly improved results of the “cube-root truly nonlinear oscillator (TNL),” “ the inverse cube-root TNL oscillator,” and “the pendulum equation.”

2. The Method

The fundamental base of the method of extended iteration is to reconvey the prestage’s solution to obtain the present stage’s solution and its sequence of extensions. The most important concern is(i)How to rearrange the preliminary shape of the given nonlinear oscillator so that it will help us to handle step by step iterations with simplification of a term including large number of harmonics.(ii)How to reuse the obtained solutions for the case of harmonic terms so that there is no algebraic complexity.

The outline of the presented procedure is as follows:(i)Step 1: considering the general form of oscillator by the following way:(ii)Step 2: setting initial conditions as(iii) Step 3: making standard form by adding on both sides, we have(iv)Step 4: formulizing the suitable iteration scheme as(v)Step 5: setting initial supposition as(vi)Step 6: setting iterated initial conditions as(vii)Step 7: executing the successive iteration from .

3. Solution Procedure of Cube-Root TNL Oscillator

Consider the cube-root TNL oscillator [25]

Now, introducing the term on equation (7), then we have the standard form:where .

Therefore, .

Using equation (4), the iteration scheme of equation (8) is

First Iteration Step. Putting in equation (9) and using initial supposition (5), we get

Now, expanding the right hand side of equation (10), then it reduces to

To avoid secular term from equation (11), we obtain

Then, equation (11) is transformed into

The required solution of equation (13) is

This is known as the first iterated approximate solution of the oscillator (7).

Second Iteration Step. For the second level of iteration, we have

Now, substituting the value of and from equations (5) and (14) into (15) and expanding the right side of equation (15), it reduces to

To eliminate secular term, we get

Then, equation (16) is transferred into

The required solution of equation (18) is

This is known as the second iterated approximate solution of the oscillator (7).

Third Iteration Step. For the third level of iteration, we have

Now, substituting the value of and from equations (5) and (19) into (20) and expanding the right side of equation (20), it reduces to

To eliminate secular term, we get

Then, equation (21) is transferred into

The required solution of equation (23) is

This is known as the third iterated approximate solution of the oscillator (7).

Fourth Iteration Step. For the fourth level of iteration, we have

Now, substituting the value of and from equations (5) and (24) into (25) and expanding the right side of equation (25), it reduces to

To eliminate secular term, we get

Equations (12), (17), (22), and (27) represent the first, second, third, and fourth approximate frequencies of the cube-root TNL oscillator (7), respectively.

4. Solution of Inverse Cube-Root TNL Oscillator

Consider the inverse cube-root TNL oscillator [25]

Applying the same iteration procedure, we obtain

Equation (29) represents the first, second, third, and fourth approximate frequencies of the inverse cube-root TNL oscillator (28) which are denoted by , , , and , respectively. Also, equation (30) represents the first, second, and third approximate analytic solutions of the inverse cube-root TNL oscillator (28) which are denoted by , and , respectively.

5. Equation of Pendulum

Consider the following equation as a pendulum equation where friction is neglected, as shown in Figure 1 [13].where , is the gravitational acceleration, is the length of the pendulum, and is the angle from the vertical to the pendulum.

To solve the pendulum equation, we consider

, then equation (30) reduces to

The iteration form of equation (32) iswhere

Applying the approximate method (6), we get

For the first iteration, we get

Now, expanding the right side of equation (36) using a Fourier series, we get

To avoid secular term, we obtain

This is the first approximate frequency of the oscillator.

After simplification in equation (37), we have

Thus, the first iterated complete solution of the oscillator (31) iswhere

For the second iteration, we get

Now, expanding the right side of equation (43) using truncation principle, we get

To avoid secular term, we obtain

This is the second approximate frequency of the oscillator.

After simplification in equation (44), we havewhere

Thus, the second iterated complete solution of the oscillator (31) is

6. Results and Discussion

To obtain approximate solutions of cube-root truly nonlinear oscillator, inverse cube-root truly nonlinear oscillator, and pendulum equation, we have applied an extended iteration method. Here, we have calculated the first, second, third, and fourth approximate frequencies of both oscillators and the first and second approximate frequencies of pendulum equation. All the results are given in Tables 1 and 2 and the frequency-amplitude relationships are given in Tables 35. To show the validity of the obtained solutions, we have compared these with the existing results determined by Mickens iteration method [25], Mickens HB method [25], He’s amplitude-frequency formulation [34], modified He’s homotopy perturbation method [15], improved harmonic balance method [6] and homotopy perturbation method [9] in Table 1, He’s energy balanced method [21] in Table 3, Mickens iteration method [25], Mickens HB method [25] in Table 2, and homotopy perturbation method [13] in Table 5. The comparison between the third-order approximate solution of equation (8) for and together with the corresponding exact solution are presented in Figures 2 and 3. The comparison between the third-order approximate solution of equation (28) for together with the corresponding exact solution are presented in Figure 4 and the comparison between the second-order approximate solution of equation (31) for and together with the corresponding exact solution are presented in Figure 5.

Here, the method of Mickens’ iteration [25] is divergent, as error pattern shows an increasing phenomenon for corresponding iteration steps. The process of Mickens’ HB method [13] is too much complicated for higher order approximations for its algebraic complexity, that is why the process was terminated up to second approximation and He could not determine the better solutions. Not only that, there is no flexibility to take any ordered harmonic terms in each iterated step. The improvement of solutions depends on only starting consideration. It is a disadvantage of this method.

For both of the oscillators, we can say, the adopted method is convergent because the solution yields lesser error in each forward iterative step:(1)For cube-root TNL oscillator ,(2)For inverse cube-root TNL oscillator ,(3)For pendulum equation ,(i)where is taken to be of very small positive quantity.(ii)Also, the consistency of the adopted method is acquired because(4)For cube-root TNL oscillator ,(5)For inverse cube-root TNL oscillator ,(6)For pendulum equation ,

7. Conclusion

After reflecting on each point of view about all the methods compared in Tables 13 and 5, we conclude that the adopted method is significantly better than each comparable stages which have been shown by other methods. Finally, we summarized:(i)The proposed method is a powerful technique for analyzing random oscillations. The technique is also easy and useful to obtain approximate frequencies and the corresponding periodic solutions of strongly nonlinear oscillator.(ii)The proposed method gives high validity for both small and large initial amplitudes of oscillations and better results of the approximate frequencies and corresponding periodic solutions than those obtained by other existing results.(iii)The percentage errors for the fourth approximate frequencies of cube-root truly nonlinear oscillator and inverse cube-root truly nonlinear oscillator are 0.006 and 0.12, respectively. Also, the percentage error for the second approximate frequency (for amplitude = 1) of equation of pendulum is 0.0065.(iv)It is investigated that most of the researchers have used the procedure to modify the method to improve the solutions in the iteration method, but we have paid attention to rearranging the leading oscillators with their own merit and choosing suitable harmonic terms from trigonometric expansion. It has been concluded that these two are also imperative matter for obtaining improved solutions.(v)The proposed method is also convergent for the considered equations.

Nomenclature

:Initial oscillation amplitude
:Dimensionless displacement
:Constant parameters
:Time
:General nonlinear functions
:Gravitational acceleration
:Length of the pendulum
:First, second, and third approximate analytic solutions of the cube-root TNL oscillator, respectively
:First, second, and third approximate analytic solutions of the inverse cube-root TNL oscillator, respectively
:First and second approximate analytic solutions of the pendulum equation, respectively
Er (%):Percentage error.
Greek Letters
:Angle from the vertical to the pendulum (°)
:First, second, third, and fourth approximate frequencies of the cube-root TNL oscillator, respectively
:First, second, third, and fourth approximate frequencies of the inverse cube-root TNL oscillator, respectively
:First and second approximate frequencies of the pendulum equation, respectively
:Very small positive quantity
:Exact frequency.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Both authors have equal contributions. They have read and approved the final version of the paper.