ABSTRACT

We applied multiple parameters method (MPM) to obtain natural frequency of the nonlinear oscillator with rational restoring force. A frequency analysis is carried out and the relationship between the angular frequency and the initial amplitude is obtained in analytical/numerical form. This equation is analyzed in three cases: the relativistic harmonic oscillator, a mass attached of a stretched elastic wire and oscillations of a punctual charge in the electric field of charged ring. The three and four parameters solutions are obtained. The results obtained are compared with the numerical solution, showing good agreement.

1. Introduction

Nonlinear oscillators play a pivotal role in physics and engineering. Recently, considerable attention has been directed towards analytical/numerical solutions of nonlinear equations. Many new techniques have been appeared in writing, for example, max-min approach [1-3], frequency amplitude formulation [4], homotopy methods [5-7], harmonic balance method [8], parameter-expanding method [9], variational approach [10-11], Hamiltonian methods [12-14] and Lindstedt-Poincaré methods with modification etc. [15-16].

In this paper, we consider a generalized nonlinear oscillator

(1)

with initial conditions

and             (2)

This equation occurs in certain phenomenon in relativistic physics, vibration of a stretched elastic wire due to mass attached to the centre and oscillation of a punctual charge in the electric field of charged ring. This equation has been investigated by various authors [17,18] for special cases.

It is interesting to note that reduce to the oscillations of a charge in the electric field of a charged ring equation. This connection is given as follows: We consider a ring of radius r with a charge spread uniformly around the ring. The electric field E on the x-axis of the ring is given by

(3)

where x is the distance along the axis. If a negative punctual charge is placed at a point on the ring axis, the charge will experience a force

(4)

The equation of motion of the punctual charge is given by the following nonlinear differential equation

(5)

Equation (5) can be written as

(6)

where

Now, dimensionlize the Equation (6) by taking

and               (7)

Substituting these dimensionless variables into Equation (6) gives

(8)

This is a special case of generalized oscillator Equaiton (1).

We assume the solution in the following form:

(9)

where are the angular frequency of motion and Fourier coefficients, respectively.

2. Governing Equation

Consider the following nonlinear equation:

(10)

where a rational function and u is the displacement. The imposed initial conditions take the forms

(11)

After multiplying both sides of Equation (10) by, integrating and using the initial condition (11), we obtain

(12)

(13)

The exact frequency of the motion can be expressed by the relation (see Equation (14)).

A general scheme of the procedure of three and four parameters is depicted in Figures 1 and 2.

3. Applications of MPM to Oscillators with Rational Restoring Force

3.1. Example for

In this case, the nonlinear equation reduces to equation of relativistic oscillator

(15)

with initial conditions

and            (16)

Multiplying both sides of Equation (15) by and integrating, with initial conditions, we get

(17)

In MPM the solution of the problem is assumed to be

(18)

Differentiating Equation (18) leads to the results

(19)

(20)

From the initial condition Equations (16) and (18), we have

(21)

Substituting Equations (18) and (19) into Equaiton

(17) at, we will find the following equation:

(22)

Considering the acceleration at the time, from Equations (15), (16) and (20), we get the following equation

(23)

From the three Equations (21)-(23), three unknowns can be solved analytically/ numerically.

(24)

The computed results of three parameters methodsfor

(25)

(14)

The frequency of three parameters method is not highly perfect. In order to find more accuracy of the solution the four parameters technique is introduce

(26)

According to the initial conditions:

(27)

Substituting Equation (26) and its derivative at and, the following equations are obtained

(28)

(29)

The acceleration at, from Equaiton (15), we will find the following equation:

(30)

After some mathematical simplification withwe achieve the numerical values

(31)

From Equation (14), we obtain the exact frequency for is

(32)

The frequency-amplitude relationship of relativistic oscillator Equation (15) obtained by Zhao [18] and Beléndez et al. [19] by frequency-amplitude formulation (FAF) and homotopy perturbation method is given by

(33)

The frequency of Equation (15) found by Shen and Mo [17] using max-min approach is given as

(34)

where

where are the complete elliptic integral and elliptic integral of the first kind respectively. The accuracy of FAF/HPM reaches 0.332517%, the accuracy of max-min approach is 3.09685% and the accuracy of four parameter lower than 0.004816% for. After comparison between the exact frequency with these methods, we conclude that the four parameter approach is better than the FAF, max-min approach and HPM.

3.2. Example for

In this case, Equation (1) is reduces to the equation of motion of a mass attached to the center of a stretched elastic wire

(35)

with initial conditions

and            (36)

Multiplying both sides of Equation (35) by and integrating, taking into account the initial conditions, we get

(37)

In MPM the solution of the problem is assumed to be

(38)

Differentiating Equaiton (38) leads to the results

(39)

(40)

From the initial condition Equations (36) and (38), we have

(41)

Substituting Equations (38) and (39) into Equation (37)

at, we will find the following equation:

(42)

Considering the acceleration at the time, from Equations (35) and (40), we get the following equation

(43)

From, Equations (41)-(43), three unknowns can be solved analytically/numerically (see Equation (44)).

In order to find more accuracy of the solution the four parameters technique is introduce

(45)

According to the initial conditions:

(46)

Substituting Equation (45) and its derivative at

and                the following equations are obtained (see Equation (47)).

(44)

(47)

(48)

From Equation (45) and (36), we will find the following equation:

(49)

After some mathematical simplification with, we achieved the numerical values

(50)

which is very close to the exact solution.

3.3. Example for

In this case, Equation (1) is reduces to the equation of a punctual charged ring, which is a free charge [21,22]

(51)

with initial conditions

and            (52)

After utilizing the same procedure in the previous examples, we obtained the approximate solutions of the problem are portrayed in Figure 5.

4. Conclusion

The multiple parameter approach has proved to be a powerful mathematical tool to find an approximate analytical solution for relativistic harmonic oscillator, a mass attached of a stretched elastic wire and oscillation of a charged ring. Comparison of the results obtained with previous methods shows that the approximate solutions are accurate and valid for the whole solution domain, and very convenient and effective. It is also found the method gives the better results if the number of parameters will be increased. It is also observed that multiple parameter approach converts the original differential equation into system of nonlinear algebraic equations. It is found that an iterative procedure for solving the corresponding system of algebraic equations creates an extremely effective method for constructing periodic solutions for nonlinear oscillators. Figures 3-5 are plotted for three, four para-

meters and exact trajectories and Table 1 shows the numerical comparisons for second examples.

5. Acknowledgements

The corresponding author is highly thankful for reviewers and the Dean faculty of Sciences for supporting this research work.

REFERENCES