Analytical solution of strongly nonlinear Duffing oscillators

Abstract

In this paper, a new perturbation technique is employed to solve strongly nonlinear Duffing oscillators, in which a new parameter Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \alpha =\alpha (\epsilon )}

is defined such that the value of α   is always small regardless of the magnitude of the original parameter Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \epsilon }

. Therefore, the strongly nonlinear Duffing oscillators with large parameter ε   are transformed into a small parameter system with respect to Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \alpha } . Approximate solution obtained by the present method is compared with the solution of energy balance method, homotopy perturbation method, global error minimization method and lastly numerical solution. We observe from the results that this method is very simple, easy to apply, and gives a very good accuracy not only for small parameter εbut also for large values of ε.

Keywords

Analytical solution; Perturbation technique; Strongly nonlinear Duffing oscillators

1. Introduction

Duffing oscillators are described by nonlinear differential equations that modeled the behavior of many practical problems that arise in engineering, physics, and in many real world applications [1], [2], [3], [4] and [5]. It is well known that Duffing oscillators can be found in the modeling of free vibrations of a restrained uniform beam with intermediate lumped mass, the nonlinear dynamics of slender elastica, the generalized Pochhammer–Chree (PC) equation, the generalized compound KdV equation in nonlinear wave systems, among others [6].

The study of Duffing oscillators has received considerable attention in recent years due to a variety of engineering applications. Several approaches have been proposed so far dealing with the nonlinear Duffing oscillators. Variational iteration method [7], homotopy perturbation method [8], He’s energy balance method [9], He’s parameter-expanding method [10], He’s max–min approach [11] and global error minimization method [12] are some examples.

In this paper, we will apply He’s modified perturbation technique [13] and [14] to solve nonlinear Duffing oscillators of fifth order in two cases first, without forced term and second with forced term, which hold for all the values of amplitude of the oscillator [15], [16] and [17].

2. Basic idea of the method

To illustrate the basic idea of the present note, we will consider the following nonlinear Duffing equation with nonlinearity of fifth order:

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \ddot{u}+u+\epsilon u^5=0\mbox{,}\quad u(0)=A\mbox{,}\quad \overset{\cdot}{u}(0)= Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): 0\mbox{.}

(1)

Re-write Eq. (1) in the form

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \ddot{u}+{\beta }^2u+\epsilon (\eta u+u^5)=0\mbox{,}\quad u(0)= Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): A\mbox{,}\quad \overset{\cdot}{u}(0)=0\mbox{,}

(2)

where

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\beta }^2+\epsilon \eta =1\mbox{,}\quad \epsilon \eta <1\mbox{.}

(3)

Here β is the angular frequency which is unknown to be further determined and ε is a large parameter.

Now define a new parameter

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \alpha =\frac{\epsilon }{1+\epsilon }\mbox{,}

(4)

such that

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \epsilon =\frac{\alpha }{1-\alpha }\mbox{.}

(5)

Then Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \alpha \in \left[0\mbox{,}1\right)} , when Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \epsilon \in \left[0\mbox{,}\infty \right)} , so ε can be written as follows:

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \epsilon =\alpha (1+\alpha +{\alpha }^2+{\alpha }^3+ Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \cdots )\mbox{.}

(6)

Assume that η and u can be written as follows:

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \eta ={\eta }_0+\alpha {\eta }_1+{\alpha }^2{\eta }_2+ Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \cdots \mbox{.}

(7)

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): u=u_0+\alpha u_1+{\alpha }^2u_2+\cdots \mbox{.}

(8)

Substituting Eqs. (6), (7) and (8) into Eq. (2) and equating coefficients of like powers of α yield the following equations:

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\ddot{u}}_0+{\beta }^2u_0=0\mbox{,}\quad u_0(0)=A\mbox{,}\quad {\overset{\cdot}{u}}_0(0)= Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): 0\mbox{,}

(9)

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\ddot{u}}_1+{\beta }^2u_1+{\eta }_0u_0+u_0^5=0\mbox{,}\quad u_1(0)= Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): 0\mbox{,}\quad {\overset{\cdot}{u}}_1(0)=0\mbox{,}

(10)

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\ddot{u}}_2+{\beta }^2u_2+{\eta }_1u_0+{\eta }_0u_1+ Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\eta }_0u_0+5u_0^4u_1+u_0^5=0\mbox{,}\quad u_2(0)= Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): 0\mbox{,}\quad {\overset{\cdot}{u}}_2(0)=0\mbox{.}

(11)

Solving Eqs. (9), (10) and (11), we obtain,

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): u_0=Acos\beta t\mbox{,}\quad {\eta }_0=-\frac{5A^4}{8}\mbox{,}

(12)

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): u_1=\frac{5A^5}{128{\beta }^2}cos3\beta t+\frac{A^5}{384{\beta }^2}cos5\beta t\mbox{,}\quad {\eta }_1= Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): -\frac{95A^8}{1536{\beta }^2}

(13)

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): u_2=\left(\frac{320A^9}{49152{\beta }^4}+\frac{5A^5}{128{\beta }^2}\right)cos3\beta t+ Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \left(\frac{320A^9}{147456{\beta }^4}+\frac{A^5}{384{\beta }^2}\right)cos5\beta t+ Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \frac{95A^9}{294912{\beta }^4}cos7\beta t+\frac{5A^9}{491520{\beta }^4}cos9\beta t\mbox{.}

(14)

Substituting Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle {\eta }_0}

and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle {\eta }_1}
into Eq. (3), we have

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\beta }^4-\left(1+\frac{5\epsilon A^4}{8}\right){\beta }^2- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \frac{95\epsilon \alpha A^8}{1536}=0\mbox{.}

(15)

Then, the second approximate solution to Eq. (1) becomes,

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): u(t)=Acos\beta t+\alpha \left(\frac{5A^5}{128{\beta }^2}cos3\beta t+\right. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \left. \frac{A^5}{384{\beta }^2}cos5\beta t\right)+ Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\alpha }^2\left\{\left(\frac{320A^9}{49152{\beta }^4}+\right. \right. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \left. \left. \frac{5A^5}{128{\beta }^2}\right)cos3\beta t+\right. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \left. \left(\frac{320A^9}{147456{\beta }^4}+\frac{A^5}{384{\beta }^2}\right)cos5\beta t+\right. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \left. \frac{95A^9}{294912{\beta }^4}cos7\beta t+\frac{5A^9}{491520{\beta }^4}cos9\beta t\right\}\mbox{.}

(16)

To show the remarkable accuracy of the obtained result, the approximate solution with the energy balance method [18], homotopy perturbation method [19], global error minimization method [20] and numerical solution are compared at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \epsilon =1}

in Table 1 (see Fig. 1).

Table 1. Comparison of the present method with the energy balance method, homotopy perturbation method, global error minimization method and numerical solution.

A	EBM [18]	HPM [19]	GEMM [20]	Present	Numerical
0.1	1.00003	1.00003	1.00003	1.00003	1.00003
0.2	1.00047	1.00050	1.00050	1.00050	1.00050
0.3	1.00236	1.00253	1.00254	1.00253	1.00253
0.5	1.01807	1.01934	1.01941	1.01940	1.01940
1	1.25831	1.27475	1.28082	1.28211	1.28079
5	19.1202	19.7895	20.3919	20.5027	19.9719
10	76.3828	79.0633	81.4807	81.9253	80.0801

Figure 1. Comparison of the approximate solution (solid line) with the numerical solution (dashed line).

3. Forced Duffing equation of fifth power nonlinear function

We consider next an oscillator with fifth power nonlinearity in the form

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \ddot{u}+u+\epsilon u^5=pcos\Omega t\mbox{,}\quad u(0)= Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): A\mbox{,}\quad \overset{\cdot}{u}(0)=0\mbox{,}

(17)

which represents a forced Duffing equation. Changing Eq. (17) into another case by using Eq. (3) to obtain

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \ddot{u}+{\beta }^2u+\epsilon (\eta u+u^5)=pcos\Omega t\mbox{.}

(18)

Substituting Eqs. (6) and (8) into Eq. (18), and setting the coefficients of the powers of α equal to zero, resulting in

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\ddot{u}}_0+{\beta }^2u_0=pcos\Omega t\mbox{,}\quad u_0(0)= Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): A\mbox{,}\quad {\overset{\cdot}{u}}_0(0)=0\mbox{,}

(19)

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\ddot{u}}_1+{\beta }^2u_1+\eta u_0+u_0^5=0\mbox{,}\quad u_1(0)= Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): 0\mbox{,}\quad {\overset{\cdot}{u}}_1(0)=0\mbox{.}

(20)

Solving Eq. (19), results in

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): u_0=\frac{p}{\left({\beta }^2-{\Omega }^2\right)}cos\Omega t+ Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \left(A-\frac{p}{{\beta }^2-{\Omega }^2}\right)cos\beta t= Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): mcos\beta t+ncos\Omega t\mbox{.}

(21)

Substituting Eq. (21) into Eq. (20), to obtain

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \eta =-\left(\frac{5}{8}m^4+\frac{30}{8}m^2n^2+\frac{15}{8}n^4\right)\mbox{.}

(22)

and

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): u_1=\frac{\left(8\eta n+5n^5+30m^2n^3+15m^4n\right)}{8\left({\beta }^2-{\Omega }^2\right)}\left(cos\beta t-\right. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \left. cos\Omega t\right)-\frac{\left(5m^5+20m^3n^2\right)}{128{\beta }^2}\left(cos\beta t-\right. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \left. cos3\beta t\right)+\frac{\left(5n^5+20m^2n^3\right)}{16\left({\beta }^2-9{\Omega }^2\right)}\left(cos\beta t-\right. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \left. cos3\Omega t\right)-\frac{m^5}{384{\beta }^2}\left(cos\beta t-\right. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \left. cos5\beta t\right)+\frac{n^5}{16\left({\beta }^2-25{\Omega }^2\right)}\left(cos\beta t-\right. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \left. cos5\Omega t\right)+\frac{5m^4n}{16\left[{\beta }^2-{\left(4\beta +\Omega \right)}^2\right]}\left(cos\beta t-\right. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \left. cos(4\beta +\Omega )t\right)+\frac{20m^4n}{16\left[{\beta }^2-{\left(2\beta +\Omega \right)}^2\right]}\left(cos\beta t-\right. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \left. cos(2\beta +\Omega )t\right)+\frac{20m^4n}{16\left[{\beta }^2-{\left(2\beta -\Omega \right)}^2\right]}\left(cos\beta t-\right. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \left. cos(2\beta -\Omega )t\right)+\frac{5m^4n}{16\left[{\beta }^2-{\left(4\beta -\Omega \right)}^2\right]}\left(cos\beta t-\right. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \left. cos(4\beta -\Omega )t\right)+\frac{5m^3n^2}{8\left[{\beta }^2-{\left(5m3\beta +2\Omega \right)}^2\right]}\left(cos\beta t-\right. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \left. cos(3\beta +2\Omega )t\right)+\frac{15m^3n^2}{8\left[{\beta }^2-{\left(\beta +2\Omega \right)}^2\right]}\left(cos\beta t-\right. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \left. cos(\beta +2\Omega )t\right)+\frac{15m^3n^2}{8\left[{\beta }^2-{\left(\beta -2\Omega \right)}^2\right]}\left(cos\beta t-\right. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \left. cos(\beta -2\Omega )t\right)+\frac{5m^3n^2}{8\left[{\beta }^2-{\left(3\beta -2\Omega \right)}^2\right]}\left(cos\beta t-\right. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \left. cos(3\beta -2\Omega )t\right)+\frac{5m^2n^3}{8\left[{\beta }^2-{\left(2\beta +3\Omega \right)}^2\right]}\left(cos\beta t-\right. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \left. cos(2\beta +3\Omega )t\right)+\frac{15m^2n^3}{8\left[{\beta }^2-{\left(2\beta +\Omega \right)}^2\right]}\left(cos\beta t-\right. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \left. cos(2\beta +\Omega )t\right)+\frac{15m^2n^3}{8\left[{\beta }^2-{\left(2\beta -\Omega \right)}^2\right]}\left(cos\beta t-\right. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \left. cos(2\beta -\Omega )t\right)+\frac{5m^2n^3}{8\left[{\beta }^2-{\left(2\beta -3\Omega \right)}^2\right]}\left(cos\beta t-\right. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \left. cos(2\beta -3\Omega )t\right)+\frac{5m\quad n^4}{16\left[{\beta }^2-{\left(\beta +4\Omega \right)}^2\right]}\left(cos\beta t-\right. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \left. cos(\beta +4\Omega )t\right)+\frac{20m\quad n^4}{16\left[{\beta }^2-{\left(\beta +2\Omega \right)}^2\right]}\left(cos\beta t-\right. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \left. cos(\beta +2\Omega )t\right)+\frac{20m\quad n^4}{16\left[{\beta }^2-{\left(\beta -2\Omega \right)}^2\right]}\left(cos\beta t-\right. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \left. cos(\beta -2\Omega )t\right)+\frac{5m\quad n^4}{16\left[{\beta }^2-{\left(\beta -4\Omega \right)}^2\right]}\left(cos\beta t-\right. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \left. cos(\beta -4\Omega )t\right)\mbox{.}

(23)

Using Eq. (8) the perturbation solution can be expressed as follows:

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): u(t)=u_0+\alpha u_1\mbox{,}

(24)

where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle u_0}

and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle u_1}
are given by Eqs. (21) and (23). As we did before we check numerically the validity of Eq. (17) by calculating the solution of Eq. (24) and also, a good agreement is found for different values of A and ε as shown in  Fig. 2.

Figure 2. Comparison of the approximate solution (solid line) with the numerical solution (dashed line).

4. Conclusion

A new perturbation technique valid for large parameter is employed in the present paper. It is an effective method to determine analytical solution for strongly nonlinear Duffing oscillators. All analytical solutions computed in the present paper are in pretty good agreement with those obtained by the numerical solutions. It is worthy to mention that the present method is an extremely simple and leads to high accuracy of the obtained results.

References

1. [1] A.H. Nayfeh, D.T. Mook; Nonlinear Oscillators; John Wiley, New York (1973)
2. [2] A.L. Maimistov; Propagation of an ultimately short electromagnetic in a nonlinear medium described by the fifth-order Duffing model; Opt. Spectrosc., 94 (2003), pp. 251–257
3. [3] A.E. Zuniga; Exact solution of the cubic-quintic Duffing oscillator; Appl. Math. Model., 37 (2013), pp. 2574–2579
4. [4] T. Öziş, A. Yildirim; Determination of limit cycles by modified straightforward expansion for nonlinear oscillators; Chaos Soliton Fract., 32 (2007), pp. 445–448
5. [5] W.U. Baisheng, L.I. Pengongli; A method for obtaining approximate analytic periods for a class of nonlinear oscillators; Meccanica, 36 (2001), pp. 167–176
6. [6] B.S. Wu, W.P. Sun, C.W. Lim; An analytical approximate technique for a class of strongly non-linear oscillators; Int. J. Non-Linear Mech., 41 (2006), pp. 766–774
7. [7] J.H. He; Variational iteration method – a kind of nonlinear analytical technique: some examples; Int. J. Non-Linear Mech., 34 (1999), pp. 699–704
8. [8] J.H. He; Homotopy perturbation method: a new nonlinear analytical technique; Appl. Math. Comput., 135 (2003), pp. 73–79
9. [9] H. Babazadeh, D.D. Ganji, M. Akbarzade; He’s energy balance method to evaluate the effect of amplitude on the nonlinear frequency in nonlinear vibration systems; Prog. Electromag. Res. M, 4 (2008), pp. 143–154
10. [10] F. Ozen Zengin, M.O. Kaya, S.A. Demirbag; Approximate period calculation for some strongly nonlinear oscillation by He’s parameter-expanding methods; Nonlinear Anal.: Real World Appl., 10 (2009), pp. 2177–2182
11. [11] L.B. Ibsen, A. Barari, A. Kimiaeifar; Analysis of highly nonlinear oscillation systems using He’s max–min method and comparison with homotopy analysis and energy balance methods; Sādhanā, 35 (2010), pp. 433–448
12. [12] M. Akbarzade, J. Langari; Solutions of nonlinear oscillators using global error minimization method; Adv. Stud. Theor. Phys., 5 (2011), pp. 349–356
13. [13] J.H. He; Analytical solution of a nonlinear oscillator by the linearized perturbation technique; Commun. Nonlinear Sci. Numer. Simul., 4 (1999), pp. 109–113
14. [14] J.H. He; Modified straightforward expansion; Meccanica, 34 (1999), pp. 287–289
15. [15] J.H. He; A new perturbation technique which is also valid for large parameters; J. Sound Vib., 229 (2000), pp. 1257–1263
16. [16] J.H. He; Iteration perturbation method for strongly nonlinear oscillations; J. Vib. Control, 7 (2001), pp. 631–642
17. [17] J.H. He; A modified perturbation technique depending upon an artificial parameter; Meccanica, 35 (2000), pp. 299–311
18. [18] H. Babazadeh, D.D. Ganji, M. Akbarzade; He’s energy balance method to evolute the effect of amplitude on the natural frequency in nonlinear vibration systems; Prog. Electromag. Res. M, 4 (2008), pp. 143–154
19. [19] J.H. He; Homotopy perturbation method: a new nonlinear analytical technique; Appl. Math. Comput., 135 (2003), pp. 73–79
20. [20] M. Akbarzade, J. Langari; Solution of nonlinear oscillator using global error minimization method; Adv. Stud. Theor. Phys., 5 (2011), pp. 349–356