# Parameters Approach Applied on Nonlinear Oscillators

Authors: Najeeb Alam Khan; Fatima Riaz; Nadeem Alam Khan
Journal: Shock and Vibration
 (2014)

---

## Abstract

We applied an approach to obtain the natural frequency of the generalized Duffing oscillatoru
¨
+
u
+
α
3
u
3
+
α
5
u
5
+
α
7
u
7
+
⋯
+
α
n
u
n
=
0 and a nonlinear oscillator with a restoring force which is the function of a noninteger power exponent of deflection u
¨
+
α
u
|
u
|
n
−
1
=
0. This approach is based on involved parameters, initial conditions, and collocation points. For any arbitrary power of n, the approximate frequency analysis is carried out between the natural frequency and amplitude. The solution procedure is simple, and the results obtained are valid for the whole solution domain.


---

## Body

## 1. Introduction

Although a large amount of the efforts on dynamical systems are related to second-order differential equations, some dynamical systems can be described by nonlinear (second-order) differential equations. Attention in nonlinear oscillator equations involving the second temporal derivative of displacement has recently been focused on the existence of periodic solutions. The study of nonlinear periodic oscillator is of interest to many researchers and various methods of solution have been suggested. Several approaches have been proposed to deal with different kinds of oscillator equations, for example, [1–7]. He in [8] used Hamiltonian method to calculate the analytical approximate periodic solutions of nonlinear oscillator equations. The approximations to the periodic solution and the angular frequency obtained by He were not accurate enough.Yildirim et al. [9] and Khan et al. [10], respectively, applied a higher order Hamiltonian formulation combined with parameters for nonlinear oscillators. Our concern in this work is the derivation of amplitude-frequency relationship for the nonlinear oscillator equations u
¨
+
u
+
α
3
u
3
+
α
5
u
5
+
α
7
u
7
+
⋯
+
α
n
u
n
=
0 and u
¨
+
α
u
|
u
|
n
-
1
=
0. The attention here has been restricted primarily to odd positive integer power for the first equation and rational powers greater than unity for the second equation. There are examples of systems, however, for which these exponents can be of noninteger order, for instance, the flexible elements of vibration isolators made of wire-mesh and felt materials, cable isolators, and radially loaded rubber cylinder.In the present work, the mentioned parameters are the undetermined values in the assumed solution. In the parameters technique, the motion has been assumed asu
=
∑
k
=
0
N
A
2
k
+
1
cos
⁡
(
2
k
+
1
)
ω
t where ω, A
k, and k
=
0,1
,
2
,
… are the angular frequency of motion and Fourier coefficients, respectively. The method in this approach to obtain the parameters is quite different from the method in He’s Hamiltonian technique. Hence, the present technique is not similar to He’s Hamiltonian technique. Finally, the paper provides some accurate results for the angular frequency ω of the motion.


## 2. Analysis of the Method

### 2.1. Generalized Duffing Oscillator

First, we consider a general form of nonlinear oscillator(1)
u
¨
+
u
+
α
3
u
3
+
α
5
u
5
+
α
7
u
7
+
⋯
+
α
n
u
n
=
0
with initial conditions
(2)
u
(
0
)
=
A
,
u
˙
(
0
)
=
0
,
where A, α
2
k
+
1, and k
=
1,2
,
3
,
… are constants. Multiplying both sides of (1) by 2
u
˙ and integrating, with initial conditions, we get
(3)
u
˙
2
+
u
2
+
α
3
u
4
2
+
α
5
u
6
3
+
⋯
+
2
α
n
u
n
+
1
n
+
1
=
A
2
+
α
3
A
4
2
+
α
5
A
6
3
+
⋯
+
2
α
n
A
n
+
1
n
+
1
.
In this approach the solution of the problem is assumed to be
(4)
u
=
∑
k
=
0
N
A
2
k
+
1
cos
⁡
(
2
k
+
1
)
ω
t
.
Differentiating (4) leads to the results
(5)
u
˙
=
-
ω
∑
k
=
0
N
A
2
k
+
1
(
2
k
+
1
)
sin
(
2
k
+
1
)
ω
t
(6)
u
¨
=
-
ω
2
∑
k
=
0
N
A
2
k
+
1
(
2
k
+
1
)
2
cos
⁡
(
2
k
+
1
)
ω
t
.
From the initial condition equations (2) and (4), we have
(7)
A
=
∑
i
=
0
N
A
2
k
+
1
.
Substituting (4) and (5) into (3) at different values of ω
t is mentioned in Table 1.Table 1
Collocation points for differentω
t.
N
=
1
π
/
2
—
—
—
N
=
2
π
/
4
π
/
2
—
—ForN
=
1 the following equation is obtained:
(8)
ω
2
(
∑
k
=
0
N
A
2
k
+
1
(
2
k
+
1
)
)
2
=
A
2
+
α
3
A
4
2
+
α
5
A
6
3
+
⋯
+
2
α
n
A
n
+
1
n
+
1
.
For N
=
2, corresponding number of equations will be obtained, respectively. Considering the acceleration at ω
t
=
0, from (1), (2), and (6), we get the following equation:
(9)
ω
2
∑
k
=
0
N
A
2
k
+
1
(
2
k
+
1
)
2
=
(
A
+
α
3
A
3
+
α
5
A
5
+
α
7
A
7
+
⋯
+
α
n
A
n
)
.
The numerical solution of these algebraic equations can be obtained by symbolic software or an iterative scheme. From these system of equations, the approximate frequency-amplitude relationship of a nonlinear oscillator up to higher order will be attained. From (3), the period of the motion is given by(10)
T
=
∫
-
A
A
2
d
u
A
2
+
α
3
(
A
4
/
2
)
+
α
5
(
A
6
/
3
)
+
⋯
+
2
α
n
(
A
n
+
1
/
(
n
+
1
)
)
-
(
u
2
+
α
3
(
u
4
/
2
)
+
α
5
(
u
6
/
3
)
+
⋯
+
2
α
n
(
u
n
+
1
/
(
n
+
1
)
)
)
.The angular frequency of the motion can be expressed by the relation
(11)
ω
=
2
π
T
.


### 2.2. Noninteger Order Force-Deflection Oscillator

Now, consider a nonlinear oscillator of the form(12)
u
¨
+
α
u
|
u
|
n
-
1
=
0
,
n
>
1
with initial conditions
(13)
u
(
0
)
=
A
,
u
˙
(
0
)
=
0
.
Multiply both sides of (12) by 2
u
˙ and integrate, with initial conditions
(14)
u
˙
2
+
2
α
u
2
|
u
|
n
-
1
n
+
1
=
2
α
A
2
|
A
|
n
-
1
n
+
1
.
Proceeding in the same way as in the previous section, the obtained algebraic equations can be solved by symbolic software. From these system of equations, one can obtain approximate frequency-amplitude relationship of a nonlinear oscillator up to higher order. From (14), the period of the motion is obtained as
(15)
T
=
∫
-
A
A
2
d
u
2
α
(
A
2
|
A
|
n
-
1
/
(
n
+
1
)
)
-
2
α
(
u
2
|
u
|
n
-
1
/
(
n
+
1
)
)
.
The angular frequency of the motion can be expressed by the relation
(16)
ω
=
2
π
(
∫
-
A
A
2
d
u
×
(
2
α
A
2
|
A
|
n
-
1
n
+
1
-
2
α
u
2
|
u
|
n
-
1
n
+
1
)
-
1
/
2
)
-
1
.
It should be pointed out that the square root under the integral sign in (16) should be positive.


## 2.1. Generalized Duffing Oscillator

First, we consider a general form of nonlinear oscillator(1)
u
¨
+
u
+
α
3
u
3
+
α
5
u
5
+
α
7
u
7
+
⋯
+
α
n
u
n
=
0
with initial conditions
(2)
u
(
0
)
=
A
,
u
˙
(
0
)
=
0
,
where A, α
2
k
+
1, and k
=
1,2
,
3
,
… are constants. Multiplying both sides of (1) by 2
u
˙ and integrating, with initial conditions, we get
(3)
u
˙
2
+
u
2
+
α
3
u
4
2
+
α
5
u
6
3
+
⋯
+
2
α
n
u
n
+
1
n
+
1
=
A
2
+
α
3
A
4
2
+
α
5
A
6
3
+
⋯
+
2
α
n
A
n
+
1
n
+
1
.
In this approach the solution of the problem is assumed to be
(4)
u
=
∑
k
=
0
N
A
2
k
+
1
cos
⁡
(
2
k
+
1
)
ω
t
.
Differentiating (4) leads to the results
(5)
u
˙
=
-
ω
∑
k
=
0
N
A
2
k
+
1
(
2
k
+
1
)
sin
(
2
k
+
1
)
ω
t
(6)
u
¨
=
-
ω
2
∑
k
=
0
N
A
2
k
+
1
(
2
k
+
1
)
2
cos
⁡
(
2
k
+
1
)
ω
t
.
From the initial condition equations (2) and (4), we have
(7)
A
=
∑
i
=
0
N
A
2
k
+
1
.
Substituting (4) and (5) into (3) at different values of ω
t is mentioned in Table 1.Table 1
Collocation points for differentω
t.
N
=
1
π
/
2
—
—
—
N
=
2
π
/
4
π
/
2
—
—ForN
=
1 the following equation is obtained:
(8)
ω
2
(
∑
k
=
0
N
A
2
k
+
1
(
2
k
+
1
)
)
2
=
A
2
+
α
3
A
4
2
+
α
5
A
6
3
+
⋯
+
2
α
n
A
n
+
1
n
+
1
.
For N
=
2, corresponding number of equations will be obtained, respectively. Considering the acceleration at ω
t
=
0, from (1), (2), and (6), we get the following equation:
(9)
ω
2
∑
k
=
0
N
A
2
k
+
1
(
2
k
+
1
)
2
=
(
A
+
α
3
A
3
+
α
5
A
5
+
α
7
A
7
+
⋯
+
α
n
A
n
)
.
The numerical solution of these algebraic equations can be obtained by symbolic software or an iterative scheme. From these system of equations, the approximate frequency-amplitude relationship of a nonlinear oscillator up to higher order will be attained. From (3), the period of the motion is given by(10)
T
=
∫
-
A
A
2
d
u
A
2
+
α
3
(
A
4
/
2
)
+
α
5
(
A
6
/
3
)
+
⋯
+
2
α
n
(
A
n
+
1
/
(
n
+
1
)
)
-
(
u
2
+
α
3
(
u
4
/
2
)
+
α
5
(
u
6
/
3
)
+
⋯
+
2
α
n
(
u
n
+
1
/
(
n
+
1
)
)
)
.The angular frequency of the motion can be expressed by the relation
(11)
ω
=
2
π
T
.


## 2.2. Noninteger Order Force-Deflection Oscillator

Now, consider a nonlinear oscillator of the form(12)
u
¨
+
α
u
|
u
|
n
-
1
=
0
,
n
>
1
with initial conditions
(13)
u
(
0
)
=
A
,
u
˙
(
0
)
=
0
.
Multiply both sides of (12) by 2
u
˙ and integrate, with initial conditions
(14)
u
˙
2
+
2
α
u
2
|
u
|
n
-
1
n
+
1
=
2
α
A
2
|
A
|
n
-
1
n
+
1
.
Proceeding in the same way as in the previous section, the obtained algebraic equations can be solved by symbolic software. From these system of equations, one can obtain approximate frequency-amplitude relationship of a nonlinear oscillator up to higher order. From (14), the period of the motion is obtained as
(15)
T
=
∫
-
A
A
2
d
u
2
α
(
A
2
|
A
|
n
-
1
/
(
n
+
1
)
)
-
2
α
(
u
2
|
u
|
n
-
1
/
(
n
+
1
)
)
.
The angular frequency of the motion can be expressed by the relation
(16)
ω
=
2
π
(
∫
-
A
A
2
d
u
×
(
2
α
A
2
|
A
|
n
-
1
n
+
1
-
2
α
u
2
|
u
|
n
-
1
n
+
1
)
-
1
/
2
)
-
1
.
It should be pointed out that the square root under the integral sign in (16) should be positive.


## 3. Numerical Experiments

### 3.1. For Duffing Oscillatorn
=
3

Consider the following nonlinear Duffing oscillator:(17)
u
¨
+
u
+
α
3
u
3
=
0
with initial conditions u
(
0
)
=
A and u
˙
(
0
)
=
0. Assume that the solution can be expressed as
(18)
u
=
A
1
cos
⁡
ω
t
+
A
3
cos
⁡
3
ω
t
.
According to the initial conditions,
(19)
A
1
+
A
3
=
A
.
Substituting (18) and its derivative into (3) at ω
t
=
π
/
2, the following equation is obtained:
(20)
ω
2
(
A
1
-
3
A
3
)
2
=
A
2
+
α
3
A
4
2
.
From (17) the acceleration at ω
t
=
0 has been obtained as follows:
(21)
ω
2
(
A
1
+
9
A
3
)
=
A
+
α
3
A
3
.
UsingMATHEMATICA the simplified values will be in the following form:
(22)
A
1
=
5
A
8
-
A
8
(
1
+
α
3
A
2
)
+
A
16
+
22
α
3
A
2
+
7
α
3
2
A
4
8
(
1
+
α
3
A
2
)
,
A
3
=
3
A
8
+
A
8
(
1
+
α
3
A
2
)
-
A
16
+
22
α
3
A
2
+
7
α
3
2
A
4
8
(
1
+
α
3
A
2
)
,
ω
(
A
)
=
1
3
(
1
+
α
3
2
A
2
+
4
+
7
2
α
3
A
2
)
.
The frequency-amplitude relationship of nonlinear oscillator for α
3
=
1, A
=
1, has been obtained as
(23)
ω
2
(
app
)
=
1.32111921
.From (3), the frequency of the motion for A
=
1 is
(24)
ω
(
Exact
)
=
(
3
/
2
)
π
2
Elliptic
K
[
-
1
/
3
]
≈
1.31777606
,
where Elliptic
K
[
m
] is the complete elliptic integral of the first kind. The three-parameter approach provides good approximations to the exact frequency and the relative error lower than 0.2536%.Since the accuracy of the obtained results in three-parameter technique is not so high, the four-parameter technique has been introduced as follows:(25)
u
=
A
1
cos
⁡
ω
t
+
A
3
cos
⁡
3
ω
t
+
A
5
cos
⁡
5
ω
t
,
where A
1, A
3, A
5, and ω are four undetermined parameters. Four equations can be formulated for the solution of four parameters. According to the initial conditions,
(26)
A
=
A
1
+
A
3
+
A
5
.
Substituting (25) and its derivative into (3) at the times ω
t
=
π
/
2 and ω
t
=
π
/
4 as mentioned in Table 1, the following equations have been obtained:
(27)
ω
2
2
(
A
1
+
3
A
3
-
5
A
5
)
2
+
u
1
2
+
α
3
u
1
4
2
=
A
2
+
α
3
A
4
2
ω
2
(
A
1
-
3
A
3
+
5
A
5
)
2
=
A
2
+
α
3
A
4
2
,
where u
1
=
(
1
/
2
)
(
A
1
-
A
3
-
A
5
).The acceleration atω
t
=
0, from (17), resulted with the following equation:
(28)
ω
2
(
A
1
+
9
A
3
+
25
A
5
)
2
=
A
+
α
3
A
3
.
It is difficult to obtain the analytical expression for the unknown values; we derived the numerical results for the unknown parameters ω, A
1, A
3, and A
5 by using (26)–(28). After some mathematical simplification by taking A
=
1
,
α
3
=
1, the following values have been obtained:
(29)
A
1
=
0.981716
,
A
3
=
0.0179609
,
A
5
=
0.000322698
,
ω
=
1.31777939
which is very close to the exact solution and highly accurate in comparison with [11, 12, 17, 18].


### 3.2. For Duffing Oscillatorn
=
5

The following nonlinear Duffing oscillator is considered:(30)
u
¨
+
u
+
α
3
u
3
+
α
5
u
5
=
0
with initial conditions u
(
0
)
=
A and u
˙
(
0
)
=
0. Let the solution be of the form
(31)
u
=
A
1
cos
⁡
ω
t
+
A
3
cos
⁡
3
ω
t
.
According to the initial conditions
(32)
A
1
+
A
3
=
A
.
Substituting (31) and its derivatives into (3) at ω
t
=
π
/
2, the following equation is obtained:
(33)
ω
2
(
A
1
-
3
A
3
)
2
=
A
2
+
α
3
A
4
2
+
α
5
A
6
3
.
The acceleration will be procured from (30) at ω
t
=
0 as follows:
(34)
ω
2
(
A
1
+
9
A
3
)
=
A
+
α
3
A
3
+
α
5
A
5
.
After some mathematical simplification usingMATHEMATICA, the following values were achieved:(35)
A
1
=
5
A
8
+
-
3
A
+
α
5
A
5
24
(
1
+
α
3
A
2
+
α
5
A
4
)
+
A
144
+
198
α
3
A
2
+
63
α
3
2
A
4
+
168
α
5
A
4
+
102
α
3
α
5
A
6
+
40
α
5
2
A
8
24
(
1
+
α
3
A
2
+
α
5
A
4
)
,
A
3
=
3
A
8
+
3
A
-
α
5
A
5
24
(
1
+
α
3
A
2
+
α
5
A
4
)
-
A
144
+
198
α
3
A
2
+
63
α
3
2
A
4
+
168
α
5
A
4
+
102
α
3
α
5
A
6
+
40
α
5
2
A
8
24
(
1
+
α
3
A
2
+
α
5
A
4
)
,
ω
=
5
9
+
4
A
2
9
α
3
+
11
A
4
27
α
5
+
1
27
144
+
198
α
3
A
2
+
63
α
3
2
A
4
+
168
α
5
A
4
+
102
α
3
α
5
A
6
+
40
α
5
2
A
8
.From (3), the exact period of the motion for A
=
1 has been obtained as
(36)
ω
(
Exact
)
=
1.52358602
.
The three-parameter approach provides good approximations to the exact frequency and the relative error lower than 1.6332%.Since the accuracy of the obtained results in three-parameter technique is not so high, the four-parameter technique has been introduced as follows:(37)
u
=
A
1
cos
⁡
ω
t
+
A
3
cos
⁡
3
ω
t
+
A
5
cos
⁡
5
ω
t
,
where A
1, A
3, A
5, and ω are four undetermined parameters. Four equations can be formulated for the solution of four parameters. According to the initial conditions,
(38)
A
=
A
1
+
A
3
+
A
5
.Substituting (37) and its derivative into (3) at the times ω
t
=
π
/
2 and ω
t
=
π
/
4 as mentioned in Table 1, the following equations have been obtained:
(39)
ω
2
(
A
1
+
3
A
3
-
5
A
5
)
2
+
u
1
2
+
α
3
u
1
4
2
+
α
5
u
1
6
3
=
A
2
+
α
3
A
4
2
+
α
5
A
6
3
ω
2
(
A
1
-
3
A
3
+
5
A
5
)
2
=
A
2
+
α
3
A
4
2
+
α
5
A
6
3
,
where u
1
=
(
1
/
2
)
(
A
1
-
A
3
-
A
5
).By the acceleration atω
t
=
0, from (30), the following equation will get
(40)
ω
2
(
A
1
+
9
A
3
+
25
A
5
)
2
=
A
+
α
3
A
3
+
α
5
A
5
.
Similarly, from (38)–(40), four unknowns ω, A
1, A
3, and A
5 can be solved numerically by using MATHEMATICA taking A
=
1 to get
(41)
A
1
=
0.968372
,
A
3
=
0.0296105
,
A
5
=
0.00201743
,
ω
=
1.52376242
which is very close to the exact solutions. The four-parameter approach provides good approximations to the exact frequency and the relative error lower than 0.0115%.


### 3.3. Force-Deflection Oscillator (for Nonintegern)

Consider the following nonlinearnon-integer force-deflection oscillator:
(42)
u
¨
+
u
|
u
|
n
-
1
=
0
,
n
>
1
with initial conditions
(43)
u
(
0
)
=
A
,
u
˙
(
0
)
=
0
.
The differential (42) with (43) has an exact analytical solution in the form of the Ateb cam function [13, 14]:
(44)
u
=
Acam
(
t
|
A
|
n
-
1
(
n
+
1
)
2
)
which is the inverse incomplete Euler Beta function. The exact period of the oscillation is
(45)
T
=
4
2
|
A
|
n
-
1
(
n
+
1
)
B
(
1
n
+
1
,
1
2
)
,
where B
(
m
,
n
) is the complete Beta function. Consider that the solution can be written as
(46)
u
=
A
1
cos
⁡
ω
t
+
A
3
cos
⁡
3
ω
t
.
According to the initial conditions,
(47)
A
1
+
A
3
=
A
.
Substituting (46) into (14) at ω
t
=
π
/
2, the following equation has been obtained:
(48)
ω
2
(
A
1
-
3
A
3
)
2
=
A
2
|
A
|
n
-
1
n
+
1
.
By the acceleration at ω
t
=
0, from (46) and (42), the following equation will be achieved:
(49)
ω
2
(
A
1
+
9
A
3
)
2
=
A
n
.
An analytical frequency-amplitude relation has been heeded from (47)–(49), byMATHEMATICA 8 built-in utilities, as
(50)
A
1
=
3
A
8
-
A
n
+
1
-
(
7
+
3
n
)
A
8
(
n
+
1
)
-
2
2
A
|
A
|
1
-
n
(
5
+
3
n
)
|
A
|
2
n
+
2
(
n
+
1
)
|
A
|
2
,
A
3
=
5
A
8
+
A
n
+
1
-
(
7
+
3
n
)
A
8
(
n
+
1
)
-
A
|
A
|
1
-
n
(
5
+
3
n
)
|
A
|
2
n
+
2
2
2
(
n
+
1
)
|
A
|
2
,
ω
=
1
3
|
A
|
n
-
1
(
n
+
1
)
(
1
+
4
+
3
n
)
.
The four-parameter approach provides good approximations to the exact frequency. The computed results and its comparison with exact frequency, Lindstedt-Poincaré (LP) method “ω
LP”, and modified Lindstedt-Poincaré (MLP) method “ω
MLP” of second order have been tabulated in Table 2.Table 2
The frequenciesω
LP and ω
MLP given by [15], ω
Present
* and ω
Present
*
* given by three-parameter and four-parameter approaches, ω
Exact given by (16) as a function of various values of the parameters n and A
=
1.
n
ω
LP
ω
MLP
ω
Present
*
ω
Present
*
*
ω
Exact
1
1
1
1
1
1
4/3
0.97013
0.95452
0.83543
0.96326
0.96916
3/2
0.95671
0.93298
0.82546
0.94793
0.95469
5/3
0.94418
0.91212
0.81649
0.93417
0.94081
2
0.92130
0.87214
0.80103
0.91045
0.91468


## 3.1. For Duffing Oscillatorn
=
3

Consider the following nonlinear Duffing oscillator:(17)
u
¨
+
u
+
α
3
u
3
=
0
with initial conditions u
(
0
)
=
A and u
˙
(
0
)
=
0. Assume that the solution can be expressed as
(18)
u
=
A
1
cos
⁡
ω
t
+
A
3
cos
⁡
3
ω
t
.
According to the initial conditions,
(19)
A
1
+
A
3
=
A
.
Substituting (18) and its derivative into (3) at ω
t
=
π
/
2, the following equation is obtained:
(20)
ω
2
(
A
1
-
3
A
3
)
2
=
A
2
+
α
3
A
4
2
.
From (17) the acceleration at ω
t
=
0 has been obtained as follows:
(21)
ω
2
(
A
1
+
9
A
3
)
=
A
+
α
3
A
3
.
UsingMATHEMATICA the simplified values will be in the following form:
(22)
A
1
=
5
A
8
-
A
8
(
1
+
α
3
A
2
)
+
A
16
+
22
α
3
A
2
+
7
α
3
2
A
4
8
(
1
+
α
3
A
2
)
,
A
3
=
3
A
8
+
A
8
(
1
+
α
3
A
2
)
-
A
16
+
22
α
3
A
2
+
7
α
3
2
A
4
8
(
1
+
α
3
A
2
)
,
ω
(
A
)
=
1
3
(
1
+
α
3
2
A
2
+
4
+
7
2
α
3
A
2
)
.
The frequency-amplitude relationship of nonlinear oscillator for α
3
=
1, A
=
1, has been obtained as
(23)
ω
2
(
app
)
=
1.32111921
.From (3), the frequency of the motion for A
=
1 is
(24)
ω
(
Exact
)
=
(
3
/
2
)
π
2
Elliptic
K
[
-
1
/
3
]
≈
1.31777606
,
where Elliptic
K
[
m
] is the complete elliptic integral of the first kind. The three-parameter approach provides good approximations to the exact frequency and the relative error lower than 0.2536%.Since the accuracy of the obtained results in three-parameter technique is not so high, the four-parameter technique has been introduced as follows:(25)
u
=
A
1
cos
⁡
ω
t
+
A
3
cos
⁡
3
ω
t
+
A
5
cos
⁡
5
ω
t
,
where A
1, A
3, A
5, and ω are four undetermined parameters. Four equations can be formulated for the solution of four parameters. According to the initial conditions,
(26)
A
=
A
1
+
A
3
+
A
5
.
Substituting (25) and its derivative into (3) at the times ω
t
=
π
/
2 and ω
t
=
π
/
4 as mentioned in Table 1, the following equations have been obtained:
(27)
ω
2
2
(
A
1
+
3
A
3
-
5
A
5
)
2
+
u
1
2
+
α
3
u
1
4
2
=
A
2
+
α
3
A
4
2
ω
2
(
A
1
-
3
A
3
+
5
A
5
)
2
=
A
2
+
α
3
A
4
2
,
where u
1
=
(
1
/
2
)
(
A
1
-
A
3
-
A
5
).The acceleration atω
t
=
0, from (17), resulted with the following equation:
(28)
ω
2
(
A
1
+
9
A
3
+
25
A
5
)
2
=
A
+
α
3
A
3
.
It is difficult to obtain the analytical expression for the unknown values; we derived the numerical results for the unknown parameters ω, A
1, A
3, and A
5 by using (26)–(28). After some mathematical simplification by taking A
=
1
,
α
3
=
1, the following values have been obtained:
(29)
A
1
=
0.981716
,
A
3
=
0.0179609
,
A
5
=
0.000322698
,
ω
=
1.31777939
which is very close to the exact solution and highly accurate in comparison with [11, 12, 17, 18].


## 3.2. For Duffing Oscillatorn
=
5

The following nonlinear Duffing oscillator is considered:(30)
u
¨
+
u
+
α
3
u
3
+
α
5
u
5
=
0
with initial conditions u
(
0
)
=
A and u
˙
(
0
)
=
0. Let the solution be of the form
(31)
u
=
A
1
cos
⁡
ω
t
+
A
3
cos
⁡
3
ω
t
.
According to the initial conditions
(32)
A
1
+
A
3
=
A
.
Substituting (31) and its derivatives into (3) at ω
t
=
π
/
2, the following equation is obtained:
(33)
ω
2
(
A
1
-
3
A
3
)
2
=
A
2
+
α
3
A
4
2
+
α
5
A
6
3
.
The acceleration will be procured from (30) at ω
t
=
0 as follows:
(34)
ω
2
(
A
1
+
9
A
3
)
=
A
+
α
3
A
3
+
α
5
A
5
.
After some mathematical simplification usingMATHEMATICA, the following values were achieved:(35)
A
1
=
5
A
8
+
-
3
A
+
α
5
A
5
24
(
1
+
α
3
A
2
+
α
5
A
4
)
+
A
144
+
198
α
3
A
2
+
63
α
3
2
A
4
+
168
α
5
A
4
+
102
α
3
α
5
A
6
+
40
α
5
2
A
8
24
(
1
+
α
3
A
2
+
α
5
A
4
)
,
A
3
=
3
A
8
+
3
A
-
α
5
A
5
24
(
1
+
α
3
A
2
+
α
5
A
4
)
-
A
144
+
198
α
3
A
2
+
63
α
3
2
A
4
+
168
α
5
A
4
+
102
α
3
α
5
A
6
+
40
α
5
2
A
8
24
(
1
+
α
3
A
2
+
α
5
A
4
)
,
ω
=
5
9
+
4
A
2
9
α
3
+
11
A
4
27
α
5
+
1
27
144
+
198
α
3
A
2
+
63
α
3
2
A
4
+
168
α
5
A
4
+
102
α
3
α
5
A
6
+
40
α
5
2
A
8
.From (3), the exact period of the motion for A
=
1 has been obtained as
(36)
ω
(
Exact
)
=
1.52358602
.
The three-parameter approach provides good approximations to the exact frequency and the relative error lower than 1.6332%.Since the accuracy of the obtained results in three-parameter technique is not so high, the four-parameter technique has been introduced as follows:(37)
u
=
A
1
cos
⁡
ω
t
+
A
3
cos
⁡
3
ω
t
+
A
5
cos
⁡
5
ω
t
,
where A
1, A
3, A
5, and ω are four undetermined parameters. Four equations can be formulated for the solution of four parameters. According to the initial conditions,
(38)
A
=
A
1
+
A
3
+
A
5
.Substituting (37) and its derivative into (3) at the times ω
t
=
π
/
2 and ω
t
=
π
/
4 as mentioned in Table 1, the following equations have been obtained:
(39)
ω
2
(
A
1
+
3
A
3
-
5
A
5
)
2
+
u
1
2
+
α
3
u
1
4
2
+
α
5
u
1
6
3
=
A
2
+
α
3
A
4
2
+
α
5
A
6
3
ω
2
(
A
1
-
3
A
3
+
5
A
5
)
2
=
A
2
+
α
3
A
4
2
+
α
5
A
6
3
,
where u
1
=
(
1
/
2
)
(
A
1
-
A
3
-
A
5
).By the acceleration atω
t
=
0, from (30), the following equation will get
(40)
ω
2
(
A
1
+
9
A
3
+
25
A
5
)
2
=
A
+
α
3
A
3
+
α
5
A
5
.
Similarly, from (38)–(40), four unknowns ω, A
1, A
3, and A
5 can be solved numerically by using MATHEMATICA taking A
=
1 to get
(41)
A
1
=
0.968372
,
A
3
=
0.0296105
,
A
5
=
0.00201743
,
ω
=
1.52376242
which is very close to the exact solutions. The four-parameter approach provides good approximations to the exact frequency and the relative error lower than 0.0115%.


## 3.3. Force-Deflection Oscillator (for Nonintegern)

Consider the following nonlinearnon-integer force-deflection oscillator:
(42)
u
¨
+
u
|
u
|
n
-
1
=
0
,
n
>
1
with initial conditions
(43)
u
(
0
)
=
A
,
u
˙
(
0
)
=
0
.
The differential (42) with (43) has an exact analytical solution in the form of the Ateb cam function [13, 14]:
(44)
u
=
Acam
(
t
|
A
|
n
-
1
(
n
+
1
)
2
)
which is the inverse incomplete Euler Beta function. The exact period of the oscillation is
(45)
T
=
4
2
|
A
|
n
-
1
(
n
+
1
)
B
(
1
n
+
1
,
1
2
)
,
where B
(
m
,
n
) is the complete Beta function. Consider that the solution can be written as
(46)
u
=
A
1
cos
⁡
ω
t
+
A
3
cos
⁡
3
ω
t
.
According to the initial conditions,
(47)
A
1
+
A
3
=
A
.
Substituting (46) into (14) at ω
t
=
π
/
2, the following equation has been obtained:
(48)
ω
2
(
A
1
-
3
A
3
)
2
=
A
2
|
A
|
n
-
1
n
+
1
.
By the acceleration at ω
t
=
0, from (46) and (42), the following equation will be achieved:
(49)
ω
2
(
A
1
+
9
A
3
)
2
=
A
n
.
An analytical frequency-amplitude relation has been heeded from (47)–(49), byMATHEMATICA 8 built-in utilities, as
(50)
A
1
=
3
A
8
-
A
n
+
1
-
(
7
+
3
n
)
A
8
(
n
+
1
)
-
2
2
A
|
A
|
1
-
n
(
5
+
3
n
)
|
A
|
2
n
+
2
(
n
+
1
)
|
A
|
2
,
A
3
=
5
A
8
+
A
n
+
1
-
(
7
+
3
n
)
A
8
(
n
+
1
)
-
A
|
A
|
1
-
n
(
5
+
3
n
)
|
A
|
2
n
+
2
2
2
(
n
+
1
)
|
A
|
2
,
ω
=
1
3
|
A
|
n
-
1
(
n
+
1
)
(
1
+
4
+
3
n
)
.
The four-parameter approach provides good approximations to the exact frequency. The computed results and its comparison with exact frequency, Lindstedt-Poincaré (LP) method “ω
LP”, and modified Lindstedt-Poincaré (MLP) method “ω
MLP” of second order have been tabulated in Table 2.Table 2
The frequenciesω
LP and ω
MLP given by [15], ω
Present
* and ω
Present
*
* given by three-parameter and four-parameter approaches, ω
Exact given by (16) as a function of various values of the parameters n and A
=
1.
n
ω
LP
ω
MLP
ω
Present
*
ω
Present
*
*
ω
Exact
1
1
1
1
1
1
4/3
0.97013
0.95452
0.83543
0.96326
0.96916
3/2
0.95671
0.93298
0.82546
0.94793
0.95469
5/3
0.94418
0.91212
0.81649
0.93417
0.94081
2
0.92130
0.87214
0.80103
0.91045
0.91468


## 4. Conclusions

The parameter method [16] gives the approximate solution for the generalized Duffing and non-integer order oscillator equations. Accuracy and validity of the obtained results have been examined by comparing it with the exact ones in time histories and table. Figures 1 and 2 are depicted for cubic and quanta oscillator equations and Figure 3 for non-integer oscillator equation. The frequency of vibration depends on the initial amplitude, the coefficient of nonlinearity, and the value of the fractional power. The main results of the paper obtained by this method can be summarized as follows.(i)
It has been observed that if numbers of parameters are increased, then the method will give the better results.
(ii)
The nonlinear oscillator equation has converted to nonlinear algebraic equation by this approach which can be solved by numerical methods that leads to a better result.
(iii)
It has also been observed that presented results for oscillator equations (as in numerical experiments Sections3.1 and 3.2) are accurate in comparison with [11, 12, 17, 18].
(iv)
The parameters method may also provide for large values of amplitude of the motion. The method is also valid for any arbitrary values ofα
3, α
5, α, and n. It will provide liberty to solve any kind of nonlinear oscillator problem which is not suitable for established methods.
(v)
Figures4–6 have been plotted for comparison of the numerical results with presented and existing methods in the literature. The LP and MLP methods give the same results for motion in the first approximation.
(vi)
The higher order MLP gives minor corrections in LP method. It has been observed that four-parameter approximations give the better results as shown in Figure5.(a) Comparison for theu versus u
˙ trajectory for the case of A
=
1, (b) t versus u (displacement A
=
1), (c) t versus u (displacement A
=
100), (d) t versus u (displacement A
=
1000).
(a)
(b)
(c)
(d)(a) Comparison for theu versus u
˙ trajectory for the case of A
=
1 and (b) t versus u (displacement A
=
1).
(a)
(b)(a) Comparison for theu versus u
˙ trajectory for the case of A
=
1 and (b) t versus u (displacement A
=
1).
(a)
(b)Comparison of the present method with numerical values and [18] for ω versus A of Duffing Oscillator (a) with n
=
3, α
3
=
1 (b) with n
=
3, α
3
=
1, α
5
=
1.
(a)
(b)Figure 5
Comparison of various methods forω versus A of force-deflection oscillator.Figure 6
Comparison of LP method and numerical values forω versus A of force-deflection oscillator.



---

*Source: 624147-2014-03-13.xml*