Combination of modified Lindstedt-Poincare and homotopy perturbation methods

Abstract

Some strongly nonlinear oscillators have been investigated with the help of combined modified Lindstedt–Poincare and homotopy perturbation methods. The solutions are more accurate than those obtained by the modified Lindstedt–Poincare method. Moreover, the determination of solutions is simpler than usual homotopy perturbation methods. Some similarities and dissimilarities between the present method and other existing methods have also been discussed.

Keywords

Nonlinear oscillation Lindstedt-Poincare homotopy perturbation method

Introduction

Mathematical modeling of the physical systems is very important and leads to differential equations many of them are strongly nonlinear such as micro-electromechanical system oscillators, vibration of the trammel pendulum, vibration of a pendulum attached to a rolling wheel, vibration of two-mass spring system etc. Solution of nonlinear differential equations plays an important role in engineering to investigate the physical systems. The methods to solve strongly nonlinear oscillators become inadequate for their analytical treatments. Many researchers have determined the approximate solutions of nonlinear problems by using various methods. Among them the Lindstedt–Poincare (LP) method^1,2 was originally developed to solve

\ddot{x} + ω_{0}^{2} x + ε f (x, \dot{x,} \ddot{x}) = 0, [x (0) = a, \dot{x} (0) = 0], 0 < ε ≪ 1,

(1)

where

ω_{0}

is a constant known as the unperturbed frequency.

Several authors^3–10 extended the method to solve strongly nonlinear problems. Among them Jones³ presented an approximate technique by introducing a new parameter, $α (ε)$ . Such approximate solution is valid for both small and large values of $ε$ . Burton⁴ and Burton et al.⁵ presented a modified version of the LP method. Cheung et al.⁶ further modified the method. However, all approximate solutions obtained by different modified approaches^3–6 are effective for Duffing oscillators with cubical nonlinearity. He^8–10 presented another modification of the LP method. He¹¹ also presented the homotopy perturbation method for solving equation (1) for large amplitudes of oscillation. In a recent article,⁷ Alam et al. provided a generalized form of the modified LP⁶ method and used it to solve various nonlinear oscillators. On the other hand, the homotopy perturbation method has been modified in some recent articles presented by Anjum et al.,¹² He et al.,¹³ Anjum et al.,¹⁴ Ji et al.,¹⁵ and He et al.¹⁶ Moreover, some approximate techniques are available to solve strongly nonlinear oscillators (Qie et al.,¹⁷ Koochi A and Goharimanesh,¹⁸ Anjum et al.,^19,20 Rehman et al.,²¹ Suleman et al.,²² Anjum et al.,²³ Li and Nadeem,²⁴ and Askari et al.²⁵).

The goal of this article is to present a combined technique of the LP and homotopy perturbation methods to solve some strongly nonlinear oscillators. The determination of solutions depends on both parameters $α (ε)$ and $p \to 1$ (well-known as a homotopy parameter¹¹). For some nonlinear problems the combined method helps to obtain significantly improved results. An earlier homotopy perturbation method was used as an alternative to the LP method to solve some truly nonlinear oscillators²⁶ where $ω_{0} = 0$ . A coupling technique of these methods was presented in article²⁷. But the present technique is completely different from article²⁷. Besides these, the present method is more systematic and easier.

The methods

The modified LP method (Cheung et al.⁶)

By a variable transformation $τ = ω t$ , equation (1) becomes

ω^{2} x^{″} + ω_{0}^{2} x + ε f (x, ω x^{'}, ω^{2} x^{″}) = 0, [x (0) = a, x^{'} (0) = 0],

(2)

where

ω

is the frequency of this oscillator and primes denote differentiation with respect to

τ

. Earlier Jones,³ Burton⁴ and Burton, and Hamdan⁵ used the following series for

ω

ω^{2} = ω_{0}^{2} + ε ω_{1} + ε^{2} ω_{2} + \dots,

(3)

where

ω_{1}, ω_{2}, \dots

are constants. In article⁶, equation (3) was replaced by

ω^{2} = (ω_{0}^{2} + ε ω_{1}) (1 + \frac{ε^{2} ω_{2}}{{ω_{0}}^{2} + ε ω_{1}} + \dots),

(4)

and then, it was transformed to

ω^{2} = \frac{ω_{0}^{2}}{1 - α} (1 + δ_{2} α^{2} + \dots), = \frac{ε ω_{1}}{{ω_{0}}^{2} + ε ω_{1}}, ε = \frac{ω_{0}^{2} α}{ω_{1} (1 - α),}

(5)

where

δ_{2}, δ_{3}, \dots

are unknown constants to be determined. Replacing

ε,

equation (3) was expanded in a series of

α

x = x_{0} + α {\tilde{x}}_{1} + α^{2} {\tilde{x}}_{2} + \dots,

(6)

By substituting equation (5) into equation (2), equation (2) then becomes

(1 + δ_{2} α^{2} + \dots) x^{″} + (1 - α) x + \frac{α f (x, ω x^{'}, ω^{2} x^{″})}{ω_{1}} = 0,

(7)

Substituting equation (6) into equation (7) and equating coefficients like powers of $α,$ a set of linear differential equations were obtained

x_{0}^{″} + x_{0} = 0,

(8)

{\tilde{x}}_{1}^{″} + {\tilde{x}}_{1} = x_{0} - \frac{1}{ω_{1}} f (x_{0}, x_{0}^{'}, x_{0}^{″}), \dots,

(9)

The initial conditions were replaced by $x_{0} (0) = a$ , $x_{0}^{'} (0) = 0$ , ${\tilde{x}}_{1} (0) = {\tilde{x}}_{1}^{'} (0) = {\tilde{x}}_{2} (0) = 0, \dots$ , and $x_{0}$ , ${\tilde{x}}_{1}$ and $ω_{1}$ , ${\tilde{x}}_{2}$ and $δ_{2}$ etc. were determined sequentially.

Combination of the modified LP and homotopy perturbation methods

According to the homotopy perturbation technique, equation (7) can be written as

σ x^{″} + σ x + p (- σ x + (1 - α) x + \frac{α f (x, ω x^{'}, ω^{2} x^{″})}{ω_{1}}) = 0,

(10)

where

p

is homotopy parameter,

ω^{2} = ω_{0}^{2} σ

(

σ

is a dimensionless quantity introduced in article⁶) and

σ = 1 + σ_{1} p + σ_{2} p^{2} + \dots, σ_{1} = δ_{2} α^{2},

(11)

If $p \to 1$ , equation (10) returns to equation (7). It is noted that the rest of the coefficients $σ_{2}, σ_{3}, \dots$ are not always equal to ${δ_{3} α}^{3}, δ_{4} α^{4}, \dots$ for every order of approximations. Therefore, $σ$ becomes slightly different from the result obtained by the modified LP method.⁶ The solution is also slightly different. A homotopy perturbation solution is considered as

x = x_{0} + p x_{1} + p^{2} x_{2} + \dots,

(12)

By substitution of equation (11) and equation (12) into equation (10) and equating coefficients of like powers of $p$ , a set of linear differential equations is obtained. To prevent secular terms, it is assumed that $x_{1}, x_{2}, \dots$ exclude first harmonic terms. Under this assumption all unknown coefficients of two series of $x$ and $σ$ are determined.

Examples (applications)

Some examples are presented in this section to illustrate the accuracy of the proposed method. Numerical results are provided and discussed for some special cases.

Example-1: Duffing oscillator with cubic nonlinearity

The structure for the buckling of a column motion of a particle on a rotating parabola shown in Figure 1 has been considered (Presented in Marinca et al.²⁸). The mass $m$ moves in the horizontal direction only. Neglecting the weight of all but the mass, the equation for the motion of $m$ is

m \ddot{x} + (k_{1} - \frac{2 p}{l}) x + (k_{3} - \frac{p}{l^{3}}) x^{3} + \dots = 0,

(13)

where the spring force is given by

F_{s p r i n g} = k_{1} x + k_{3} x^{3} + \dots,

(14)

Figure 1.

Model for the buckling of the column motion of a particle on a rotating parabola.

Now equation (13) can be written in the form

m \ddot{x} + α_{1} x + α_{3} x^{3} + \dots = 0,

(15)

which is considered as Duffing equation

\ddot{x} + x + ε x^{3} = 0,

(16)

Let us consider the initial condition as x $(0) = a, \dot{x} (0) = 0$ .

For equation (16), $f = x^{3}$ and equation (10) becomes

x^{″} + x + p (- x + ((1 - α) x + \frac{α x^{3}}{ω_{1}}) / σ) = 0,

(17)

According to the processes discussed in the previous section, the linear differential equations for $x_{1}, x_{2}, \dots$ have been obtained

x_{0}^{″} + x_{0} = 0,

(18)

x_{1}^{″} + x_{1} = α (- x_{0} + \frac{x_{0}^{3}}{ω_{1}}) \dots,

(19)

x_{2}^{″} + x_{2} = α (- x_{1} + \frac{3 x_{0}^{2} x_{1}}{ω_{1}}) - σ_{1} ((1 - α) x_{0} + \frac{x_{0}^{3}}{ω_{1}}), \dots,

(20)

Solving the above set of equations with the initial conditions described in the previous section, we obtained the following results

σ_{1} = - \frac{α^{2}}{24}, σ_{2} = - \frac{α^{4}}{576}, σ_{3} = \frac{7 α^{4}}{13824} - \frac{α^{6}}{6912}, σ_{4} = \frac{7 α^{6}}{82944} - \frac{5 α^{8}}{331776}, \dots,

(21)

and

x_{0} = a \cos τ, x_{1} = \frac{a α}{24} (\cos 3 τ - \cos τ),

x_{2} = \frac{a α^{3}}{72} (\cos 3 τ - \cos τ) + \frac{a α^{2}}{24} (\cos 5 τ - \cos τ),

(22)

x_{3} = \frac{a (3 - α^{2})}{6912} (\cos 3 τ - \cos τ) - \frac{a α^{4}}{6912} (\cos 5 τ - \cos τ) - \frac{a α^{3}}{13824} (\cos 7 τ - \cos τ), \dots,

Substituting $p = 1$ together with the values of $σ_{1}, σ_{2}, \dots$ , equation (22) becomes

σ = 1 + \sum_{i = 1}^{4} σ_{i} = 1 - \frac{α^{2}}{24} - \frac{17 α^{4}}{13824} - \frac{5 α^{6}}{82944} - \frac{5 α^{8}}{331776} .

(23)

By the modified LP method,⁶ the following results was obtained

δ_{2} = - \frac{α^{2}}{24}, δ_{3} = 0, δ_{4} = - \frac{17 α^{4}}{13824}, δ_{5} = 0, \dots,

σ = 1 + \sum_{i = 1}^{4} δ_{i + 1} = 1 - \frac{α^{2}}{24} - \frac{17 α^{4}}{13824} .

(24)

In both relations of $σ$ , the results are the same until $α^{4}$ and it is the correct result. If the next terms $σ_{5}$ and $δ_{6}$ are added respectively in the series equation (23) and Eq. (24), the results will be the same until $α^{6}$ and later series (presented by equation (24)) do not contain more than $α^{6}$ ordered terms. On the contrary, the former contains more two terms up to $α^{10}$ . Again, all the coefficients are negative. It is observed that all approximate results of $σ$ are smaller than its exact value and the former (obtained by a combined technique) provides better results for this oscillator.

Example-2: nonlinear oscillator with inertia-type cubic nonlinearity

The motion of a ring of mass $m$ sliding freely on a wire described by the parabola $y = q x^{2}$ which rotates with a constant angular velocity $Ω$ about the y-axis as shown in Figure 2 (Nayfeh and Mook²). The equation describing the motion of the ring is (Nayfeh and Mook,² Marinca and Herinasu²⁸).

\ddot{x} + ω^{2} x = - 4 q x (x \ddot{x} + {\dot{x}}^{2}),

(25)

where

ω^{2} = 2 g q - Ω^{2}

and the initial conditions are

x (0) = A, \dot{x} (0) = 0 .

Figure 2.

Motion of a particle on a rotating parabola.

Equation (25) can be written as²⁹

\ddot{x} + x + ε (x^{2} \ddot{x} + {\dot{x}}^{2} x) = 0,

(26)

with its transformed form

ω^{2} x^{″} + x + {ε ({{ω^{2} x}^{2} x^{″} + ω}^{2} x^{'}}^{2} x) = 0,

(27)

or,

x^{″} + ω^{- 2} x + {ε (x^{2} x^{″} + x^{'}}^{2} x) = 0,

(28)

For this oscillator, by choosing the following series (see article⁷ for details)

ω^{- 2} = 1 + ε ω_{1} + ε^{2} ω_{2} + \dots,

(29)

and the corresponding equation of equation (7) becomes

(1 - α) x^{″} + σ x + \frac{α}{ω_{1}} {(x^{2} x^{″} + x^{'}}^{2} x) = 0,

(30)

where

σ

follows equation (11).

Now adding a homotopy parameter it becomes

σ x^{″} + σ x + p (- σ x^{″} + (1 - α) x^{″} + \frac{α}{ω_{1}} {(x^{2} x^{″} + x^{'}}^{2} x)) = 0,

(31)

For equation (31), the following results have been obtained

σ_{1} = - \frac{α^{2}}{8}, σ_{2} = - \frac{α^{4}}{64}, σ_{3} = - \frac{5 α^{4}}{512} - \frac{α^{6}}{256}, σ_{4} = - \frac{5 α^{4}}{512} - \frac{5 α^{8}}{4096}, \dots,

and then,

σ = 1 + \sum_{i = 1}^{4} σ_{i} = 1 - \frac{α^{2}}{8} - \frac{13 α^{4}}{512} - \frac{9 α^{6}}{1024} - \frac{5 α^{8}}{4096},

(32)

On the other hand, the calculated result by the modified LP method⁷ is

σ = 1 + \sum_{i = 1}^{4} δ_{i + 1} = 1 - \frac{α^{2}}{8} - \frac{13 α^{4}}{512},

(33)

All these results are similar to those obtained for the cubic Duffing oscillator (i.e., the coefficients of $α^{3}, α^{5}, \dots$ vanish). For both, oscillators $δ_{2 i + 1} = 0, i = 0, 1, 2, \dots$ and coefficients are much smaller than $1$ .

Example-3: duffing oscillator with quintic nonlinearity

Let us consider a highly nonlinear oscillator, that is, Duffing oscillator with quintic nonlinearity

\ddot{x} + ω_{0}^{2} x + {ε x}^{5} = 0,

(34)

For this oscillator the following results have been obtained

σ_{1} = - \frac{{13 α}^{2}}{120}, σ_{2} = \frac{2 α^{3}}{225} - \frac{169 α^{4}}{14400}, σ_{3} = \frac{1673 α^{4}}{575000} + \frac{13 α^{3}}{4500} - \frac{2197 α^{6}}{864000},

σ_{4} = - \frac{253 α^{5}}{270000} + \frac{11411 α^{6}}{10368000} + \frac{169 α^{7}}{162000} - \frac{28561 α^{8}}{41472000}, \dots,

and

σ = 1 + \sum_{i = 1}^{4} σ_{i} = 1 - \frac{13 α^{2}}{120} + \frac{2 α^{3}}{225} - \frac{5087 α^{4}}{576000} + \frac{527 α^{5}}{270000} - \frac{14953 α^{6}}{10368000} + \frac{169 α^{7}}{162000} - \frac{28561 α^{8}}{41472000},

(35)

On the other hand, the calculated result by the modified LP method⁷ is

σ = 1 + \sum_{i = 1}^{4} δ_{i + 1} = 1 - \frac{13 α^{2}}{120} + \frac{2 α^{3}}{225} - \frac{5087 α^{4}}{576000} + \frac{527 α^{5}}{270000},

(36)

The coefficients of equation (36) are also much smaller than $1$ , but larger than those of equation (23) and equation (24). Moreover, all $δ_{2 i + 1}, i = 0, 1, 2, \dots$ are positive. For this reason, the results oscillate and gradually approach its exact value. However, the errors of results obtained by equation (35) are much smaller than those obtained by equation (36).

Example-4: duffing oscillator with septic nonlinearity

Let us consider another highly nonlinear oscillator, that is, Duffing oscillator with septic nonlinearity

\ddot{x} + ω_{0}^{2} x + {ε x}^{7} = 0,

(37)

For this oscillator the following results have been obtained

σ_{1} = - \frac{{99 α}^{2}}{560}, σ_{2} = \frac{141 α^{3}}{4480} - \frac{9801 α^{4}}{313600}, σ_{3} = \frac{695357 α^{4}}{175616000} + \frac{41877 α^{5}}{2508800} - \frac{970299 α^{6}}{87808000},

σ_{4} = - \frac{2917107 α^{5}}{614656000} + \frac{20131893 α^{6}}{24586240000} + \frac{1381941 α^{7}}{140492800} - \frac{96059601 α^{8}}{19668992000}, \dots,

and

σ = 1 + \sum_{i = 1}^{4} σ_{i} = 1 - \frac{99 α^{2}}{560} + \frac{141 α^{3}}{4480} - \frac{4793203 α^{4}}{175616000} + \frac{3671379 α^{5}}{307328000} - \frac{251551827 α^{6}}{24586240000} + \frac{1381941 α^{7}}{140492800} - \frac{96059601 α^{8}}{19668992000},

(38)

On the other hand, the calculated result by the modified LP method⁷ is

σ = 1 + \sum_{i = 1}^{4} δ_{i + 1} = 1 - \frac{99 α^{2}}{560} + \frac{141 α^{3}}{4480} - \frac{4793203 α^{4}}{175616000},

(39)

The coefficients of equation (39) are also much smaller than $1$ but larger than those of equation (23) and equation (24). Moreover, all $δ_{2 i + 1}, i = 0, 1, 2, \dots$ are positive. For this reason, the results oscillate and gradually approach its exact value. However, the errors of results obtained by equation (38) are much smaller than those obtained by equation (39).

Convergency of the $α$ -series

For a lot of nonlinear oscillators, the $α$ -series of the modified LP method⁷ converge rapidly when $α ≪ 1$ (i.e., in the case of small oscillations). In the case of large amplitude of oscillations (i.e., $α \to 1$ ) the series also converges since the coefficients of the series, $δ_{i}, i = 2, 3, \dots$ gradually decrease. For some of the oscillators, the coefficients $δ_{i}, i = 3, 5, \dots$ vanish. According to the recently presented homotopy perturbation method,¹⁶ the coefficients $ω_{i}, i = 3, 5, \dots$ vanish for oscillator presented by equation (16). In these situations, the frequency of oscillation is same for second and third approximations; and then for fourth and fifth approximations. On the other hand, the coefficients $δ_{i}, i = 3, 5, \dots$ are different in sign of the coefficients $δ_{i}, i = 2, 4, \dots$ for oscillators presented by equation (34) and equation (37). In these situations, approximate results of frequencies oscillate. For these nonlinear oscillators, the calculated frequencies obtained by the homotopy perturbation method¹⁶ also oscillates. Moreover, homotopy perturbation method¹⁶ fails to measure frequency as well as solution of oscillator presented by equation (26) especially in the case of amplitude of oscillation. All these problems have been satisfactorily removed in this article.

Results and discussion

Recently the modified Lindstedt–Poincare method⁶ has been generalized for solving strongly nonlinear oscillators.⁷ The generalized method covers wide varieties of nonlinear problems where $ω_{0} \neq 0$ . (It is noted that the modified technique⁶ is useless for a class of nonlinear oscillators, that is, second example equation (26)). Moreover, the generalized version⁷ has been combined with the homotopy perturbation technique for some truly nonlinear oscillators where $ω_{0} = 0$ . In this article the method has been combined with homotopy perturbation for all values of $ω_{0}$ whether it vanishes or not. Earlier truly nonlinear oscillators were solved by the homotopy method,¹¹ but with a different approach. The new technique is more systematic than all the past methods mentioned above. For large amplitudes of oscillations $α \to 1$ and errors of frequency become the largest. To compare the results with existing results, the ratios of frequency with their exact frequency are presented below.

For the cubic Duffing oscillator equation (16), the ratios respectively become:

By the present method with equation (23): $\underset{.}{\begin{array}{c} \lim \\ ε a^{2} \to \infty \end{array}} \frac{ω}{ω_{E x a c t}} ≅ 1.0000008 .$

and by the modified LP method⁶ with equation (24): $\underset{\begin{array}{c} \lim \\ ε a^{2} \to \infty \end{array}}{.} \frac{ω}{ω_{E x a c t}} ≅ 1.00004 .$

For the nonlinear oscillators equation (26), the ratios respectively become

by the present method with equation (32): $\underset{\begin{array}{c} \lim \\ ε a^{2} \to \infty \end{array}}{.} \frac{ω}{ω_{E x a c t}} ≅ 0.98256,$

and by the modified LP method,⁷ with equation (33): $\underset{\begin{array}{c} \lim \\ ε a^{2} \to \infty \end{array}}{.} \frac{ω}{ω_{E x a c t}} ≅ 0.97676,$

For the highly nonlinear oscillators equation (34), the ratios respectively become

by the present method with equation (35): $\underset{.}{\begin{array}{c} \lim \\ ε a^{4} \to \infty \end{array}} \frac{ω}{ω_{E x a c t}} ≅ 1.000095 \dots,$

and by the modified LP method,⁷ with equation (36) $\underset{\begin{array}{c} \lim \\ ε a^{4} \to \infty \end{array}}{:} \frac{ω}{ω_{E x a c t}} ≅ 1.000704 \dots .$

For the Duffing oscillator with cubic nonlinearity, the second, third, fourth, and fifth approximated frequencies with their percentage error are given in Table 1.

Table 1.

Comparison of the frequencies among the present method, modified Lindstedt–Poincare method,⁷ and the homotopy method¹⁶ for the Duffing oscillator with cubic nonlinearity for $a = 100$ .

Approx.	Present method (% Err.)	MLP method⁷ (% Err.)	Homotopy method¹⁶ (%Err.)
2nd	84.7853 (0.06820)	84.7853 (0.06820)	84.7013 (0.03087)
3rd	84.7085 (0.02243)	84.7853 (0.06820)	84.7013 (0.03087)
4th	84.7245 (0.00354)	84.7309 (0.00401)	84.7281 (0.00073)
5th	84.7275 (0.00008)	84.7309 (0.00401)	84.7281 (0.00073)

Here % Error denotes the absolute percentage error.

Moreover, for the highly nonlinear oscillators equation (34), the second, third, fourth, and fifth approximated frequencies with their percentage error are given in Table 2.

Table 2.

Comparison of the frequencies among the present method, modified Lindstedt–Poincare method,⁷ and the homotopy method¹⁶ for the quintic oscillator $(\ddot{x} + ω_{0}^{2} x + {ε x}^{5} = 0)$ for $a = 100$ .

Approx.	Present method (% Err.)	MLP method⁷ (% Err.)	Homotopy method¹⁶ (% Err.)
2nd	7465.2 (0.0420)	7465.2 (0.0420)	7400.95 (0.90234)
3rd	7465.44 (0.03883)	7502.31 (0.4548)	7455.8541 (0.16718)
4th	7467.57 (0.01031)	7465.44 (0.0388)	7474.9977 (0.089146
5th	7469.05 (0.0095)	7473.60 (0.0704)	7468.15 (0.00254)

Here, % error denotes the absolute percentage error.

The present solution along with numerical, Modified Lindstedt–Poincare⁶ and generalized Lindstedt–Poincare⁷ solutions are presented in Figures 3–7. Also, the frequency analysis of the present method with existing methods is provided in Figures 8 and 9.

Figure 3.

Comparison between the time versus displacement results obtained by the present method, numerical method, Cheung et al.⁶ method, and Alam et al.⁷ method for the Duffing oscillator with cubic nonlinearity when $a = 0.5, ε = 1$ .

Figure 4.

Comparison between the time versus displacement results obtained by the present method, numerical method, and Cheung et al.⁶ method for the Duffing oscillator with cubic nonlinearity when $a = 100, ε = 1$ .

Figure 5.

Comparison between the time versus displacement results obtained by the present method, numerical method, Cheung et al.⁶ method, and Alam et al.⁷ method for equation (26) when $a = 0.7, ε = 1$ .

Figure 6.

Comparison between the time versus displacement results obtained by the present method, numerical method, and Alam et al.⁷ method for equation (26) when $a = 5.0, ε = 1$ .

Figure 7.

Comparison between the time versus displacement results obtained by the present method, numerical method, Cheung et al.⁶ method, and Alam et al.⁷ method for the quintic oscillator $(\ddot{x} + ω_{0}^{2} x + {ε x}^{5} = 0)$ when $a = 0.7, ε = 1$ .

Figure 8.

Comparison of calculated frequencies measured by the present technique and earlier modified technique (Alam et al.⁷), homotopy method¹⁶ together with the numerical results for the quintic oscillator $(\ddot{x} + ω_{0}^{2} x + {ε x}^{5} = 0)$ when $a = 100, ε = 1$ .

Figure 9.

Comparison of calculated frequencies measured by the present technique and earlier modified technique (Alam et al.⁷) together with the numerical results for the septic oscillator $(\ddot{x} + ω_{0}^{2} x + {ε x}^{7} = 0)$ when $a = 10.0, ε = 1$ .

All these figures show that the present results have good agreement with the numerical results; on the other hand, in some cases, the present results are better (Figure 3, Figures 5–9) than other existing results. Also, for a large amplitude, the present results are converging but other results are diverging, that is, the results obtained by Cheung et al.⁶ diverge for equation (26) and the results obtained by Alam et al.⁷ diverge for Duffing cubic, quintic, and septic oscillators.

Conclusion

Earlier modified LP method and homotopy perturbation method were presented for solving strongly nonlinear oscillators especially in case of large amplitude of oscillation. Based on these methods, a combined technique has been presented which provides more satisfactory results. The necessity of present technique has been illustrated by considering four different nonlinear oscillators. Calculated frequencies by the modified LP method of cubical Duffing oscillator (i.e., equation (24)) gradually closes to numerical result with increasing order of approximation. But present frequency (equation (23)) closes faster than that of the modified LP method and a recent homotopy perturbation method.¹⁶ For the quintic oscillator, coefficients $δ_{i}, i = 3, 5, \dots$ are different in sign of $δ_{i}, i = 2, 4, \dots$ . It is noted that $δ_{i} = 0, i = 3, 5, \dots$ for cubical oscillator. So, major problem arises in case of quintic oscillator. In this situation error of the frequency (equation (36), which is obtained by modified LP method) oscillates. Similar problem arises in septic oscillator. Homotopy perturbation method¹⁶ measures also oscillating results for these oscillators. But present method provides more accurate result. Moreover, the present method covers various nonlinear oscillators involving $x$ , $\dot{x}$ and $\ddot{x}$ (e.g., equation (26)). It is noted that the homotopy perturbation method¹⁶ is useless for solving equation (26). Homotopy perturbation method¹⁶ has another limitation. According to this method, n-th ( $n = 2, 3, \dots)$ approximate frequency are calculated from n-th order algebraic equation. On the contrary, modified LP method and present method are straightforward. However, the present method has a limitation for solving the quadratic oscillators. It measures similar results of modified LP method.⁶

Footnotes

Acknowledgements

The authors are grateful to the honorable reviewers for their helpful comments and suggestions in improving the manuscript.

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) received no financial support for the research, authorship, and/or publication of this article.

ORCID iD

Md Helal Uddin Molla

References

Nayfeh

. Perturbation methods. New York, NY: Wiley & Sons, 1973.

Nayfeh

Mook

. Nonlinear oscillations. New York, NY: Wiley & Sons, 1979.

Jones

. Remarks on the perturbation process for certain conservative systems. Int J Non-linear Mech 1978; 13: 125–128.

Burton

. A perturbation method for certain nonlinear oscillators. Int J Non-linear Mech 1984; 19: 397–407.

Burton

Hamdan

. Analysis of nonlinear autonomous conservative oscillators by time transformation method. J Sound Vib 1983; 87: 543–554.

Cheung

Chen

Lau

. A modified Lindstedt-Poincare method for certain strongly nonlinear oscillators. Int J Non-linear Mech 1991; 26: 367–378.

Alam

Yeasmin

Ahamed

. Generalization of the modified Lindstedt-Poincare method for solving some strong nonlinear oscillators, Ain Shams Eng J. 2019; 10: 195–201.

. Modified Lindstedt-Poincare methods for some strongly nonlinear oscillations-I: expansion of a constant. Int J Nonlinear Mech 2002; 37: 309–314.

. Modified Lindstedt-Poincare methods for some strongly non-linear oscillations Part II: a new transformation. Int J Non-linear Mech 2002; 37: 315–320.

10.

. Modified Lindsted-Poincare methods for some strongly nonlinear oscillations Part III: double series expansion. Inter J Nonlinear Sci Numer Simul 2001; 2: 317–320.

11.

. Homotopy perturbation method for bifurcation and nonlinear problems. Int J Non-linear Sci Numer Simul 2005; 6: 207–208.

12.

Anjum

Ain

, et al. Li-He’s modified homotopy perturbation method for doubly-clamped electrically actuated microbeams-based microelectromechanical system. Facta Univ-Ser Mech Eng 2021; 19(4): 601–612.

13.

El-Dib

. The enhanced homotopy perturbation method for axial vibration of strings. Facta Univ-Ser Mech Eng 2021; 19(4): 735–750.

14.

Anjum

. Homotopy perturbation method for N/MEMS oscillators. In: Mathematical methods in the applied sciences. Hoboken, NJ: Wiley, 2020.

15.

Wang

, et al. Li–He’s modified homotopy perturbation method coupled with the energy method for the dropping shock response of a tangent nonlinear packaging system. J Low Freq Noise Vib Act Control 2021; 40(2): 675–682.

16.

J-H

Jiao

M-L

Gepreeld

, et al. Homotopy perturbation method for strongly nonlinear oscillators. J Math Comput Simul 2023; 204: 243–258.

17.

Qie

Houa

. The fastest insight into the large amplitude vibration of a string. Rep Mech Eng 2021; 2(1): 1–5.

18.

Koochi

Goharimanesh

. Nonlinear oscillations of CNT nano-resonator based on nonlocal elasticity: the energy balance method. Rep Mech Eng 2021; 2(1): 41–50.

19.

Anjum

, et al. A brief review on the asymptotic methods for the periodic behaviour of microelectromechanical systems. J Appl Comput Mech 2022; 8(3): 1120–1140.

20.

Anjum

Ul Rahman

, et al. An efficient analytical approach for the periodicity of nano/microelectromechanical systems’ oscillators. Math Probl Eng 2022; 2022: 9712199.

21.

Rehman

Hussain

Rahman

JUl

, et al. Modified laplace based variational iteration method for the mechanical vibrations and its applications. Act Mech et Autom 2022; 16(2): 98–102.

22.

Suleman

Elzaki

, et al. New application of Elzaki projected differential transform method. J Comput Theoret Nanosci 2017; 14(1): 631–639.

23.

Anjum

Suleman

, et al. Numerical iteration for nonlinear oscillators by Elzaki transform. J Low Freq Noise Vib Act Control 2020; 39(4): 879–884.

24.

Nadeem

. He–laplace method for nonlinear vibration in shallow water waves. J Low Freq Noise Vib Act Control 2019; 38(3–4): 1305–1313.

25.

Askari

Younesian

Saadatnia

. Nonlinear oscillations analysis of the elevator cable in a drum drive elevator system. Adv Appl Math Mech 2015; 7(1): 43–57.

26.

Belendez

Pascual

Marquez

, et al. Application of He’s homotopy perturbation method to the Duffing-harmonic oscillator. Int J Non-linear Sci Numer Simul 2007; 8: 79–88.

27.

. A coupling method of homotopy technique and perturbation technique for nonlinear problems. Int J Non-linear Mech 2000; 35: 37–43.

28.

Marinca

Herinasu

. A modified iteration perturbation method for some nonlinear oscillation problems. Acta Mech 2006; 184: 142–231.

29.

Liao

. An approximate solution technique not depending on small parameters: a special example. Int J Non-linear Mech 1995; 30: 371–380.

Combination of modified Lindstedt-Poincare and homotopy perturbation methods

Abstract

Keywords

Introduction

The methods

The modified LP method (Cheung et al. 6 )

Combination of the modified LP and homotopy perturbation methods

Examples (applications)

Example-1: Duffing oscillator with cubic nonlinearity

Example-2: nonlinear oscillator with inertia-type cubic nonlinearity

Example-3: duffing oscillator with quintic nonlinearity

Example-4: duffing oscillator with septic nonlinearity

Convergency of the α -series

Results and discussion

Conclusion

Footnotes

Acknowledgements

Declaration of conflicting interests

Funding

ORCID iD

References

The modified LP method (Cheung et al.⁶)

Convergency of the $α$ -series