Approximate Periodic Solution for the Nonlinear Helmholtz-Duffing Oscillator via Analytical Approaches

Abstract

The conservative Helmholtz-Duffing oscillator is analyzed by means of three analytical techniques. The max-min, second-order of the Hamiltonian, and the global error minimization approaches are applied to achieve natural frequencies. The obtained results are compared with the homotopy perturbation method and numerical solutions. The results show that second-order of the global error minimization method is very accurate, so it can be widely applicable in engineering problems.

1. Introduction

Mathematical modeling and frequency analysis of the nonlinear vibrational systems are an important and interesting field of mechanics. A lot of researchers have worked in this field and have proposed a lot of methods for demonstrating the dynamics responses of these systems [1–4]. They have developed this field of science and have analyzed the responses of the nonlinear vibration problems such as Duffing oscillators [5–10], nonlinear dynamics of a particle on a rotating parabola [11], nonlinear oscillators with discontinuity [12], oscillators with noninteger order nonlinear connection [13], the plasma physics equation [14], and van der Pol oscillator [15, 16]. The Helmholtz-Duffing equation is a nonlinear problem with the quadratic and cubic nonlinear terms. Surveying the literature shows that this equation has wide applications in the engineering problems. For example, due to different vibration behavior of functionally graded materials (FGMs) at positive and negative amplitudes, the governing equations of FGM beams, plates, and shells are conduced to a second-order nonlinear ordinary equation with quadratic and cubic nonlinear terms [17–20]. Moreover, Sharabiani and Yazdi [21] obtained a Helmholtz-Duffing type equation within studying of nonlinear free vibrations of functionally graded nanobeams with surface effects. On the other hand, they revealed application of this equation in FG nanostructures.

In this paper, the frequency-amplitude relationship of the conservative Helmholtz-Duffing oscillator is obtained by means of the max-min [22–26], Hamiltonian [27–31], and global error minimization methods [32–35]. The Hamiltonian approach is a kind of energy method and is proposed by He [27]. It is a simple method and can be used for the conservative nonlinear equations. Recently, it is applied for dynamic analysis of an electromechanical resonator [36]. Moreover, Akbarzade and Khan [37] employed the second-order Hamiltonian approach for nonlinear dynamic analysis of conservative coupled systems of mass-spring. The max-min approach is made on the base of Chengtian’s inequality [38]. It is a valuable method for obtaining the frequency responses of the nonlinear problems, and many researchers are attracted to use this method for studying the nonlinear systems. The global error minimization method is a modified type of variational approach and converts the nonlinear equation to an equivalent minimization problem.

The Helmholtz-Duffing equation is analyzed by many researchers. For instance, Leung and Guo [39] have applied the homotopy perturbation method (HPM) to this equation and have obtained accurate responses. Askari et al. [40] studied approximate periodic solutions of this equation using He’s energy balance method (HEBM) and He’s frequency-amplitude formulation (HFAF). Akbarzade et al. [41] used the first-order of the Hamiltonian approach and coupled homotopy-variational formulation to study the periodic solutions of the Helmholtz-Duffing oscillator; they also discussed the stability of the system for selected constant parameters. Recently, Li et al. [42] determined limit cycles and homoclinic orbits of this oscillator via a generalized harmonic function perturbation method. In the next sections, the max-min, the Hamiltonian, and the global error minimization methods are applied for evaluating the dynamics responses of the Helmholtz-Duffing equation. The results of these methods are compared with the exact ones and HPM solutions.

2. Solution Procedure

Let us consider the Helmholtz-Duffing equation The response oscillates between an asymmetric limit zone , where both and are positive. One can determine in the form of [39] where

2.1. Max-Min Approach (MMA)

We can rewrite (1) in the following form:

By choosing as a trial function that satisfied the initial conditions, the maximum and minimum values of can be calculated with Maple software as and , respectively, so we can write

By using Chengtian interpolation [38], we have where , are weighting factor, . Therefore the frequency can be approximated as Then the solution of (1) can be written:

By using the approximate solution, (4) can be written in the following form:

If by chance (8) is the exact solution of (1), then the right-hand side of (9) vanishes. According to [22], we set where . By substituting (8) into (10), we obtain the following expression for :

Substituting (11) into (7) and after some simplification, we have

2.2. Hamiltonian Approach

2.2.1. The First-Order Hamiltonian Approach (FHA)

The Hamiltonian of (1) is constructed as

Assume the first approximate solution of (1) as Then integrating (13) with respect to time from 0 to , we have Substituting (14) in (15) yields

Set

Finally, The result of first-order Hamiltonian and max-min approaches is alike.

2.2.2. The Second-Order Hamiltonian Approach (SHA)

For the second-order Hamiltonian approach, we consider the following equation as the response of the system: Whereas (19) must satisfy the initial condition we have Substituting (19) into (15) yields

Set Since finding a relation between different parameters from (22) is not easy, the value of other parameters is determined in Table 1 for several values of .

2.3. Global Error Minimization Method

2.3.1. The First-Order Global Error Minimization Method (FGEM)

Based on standard procedure of modified variational approach, the minimization problem is For first-order approximation, consider a trial function as follows: Substituting (24) into (23) yields The solution of (25) could be found through By some simplifications, the following equation is obtained: One may find the first-order approximation by solving (27) as

2.3.2. The Second-Order Global Error Minimization Method (SGEM)

For the second-order approximation, consider a trial function as Whereas (29) must satisfy the initial condition, we have Substituting (29) into (23) yields

The solution of (31) could be found through which yield the following:

Finally, by simultaneously solving (30) and (33) unknown parameters are determined for different values of and as the second-order global error minimization method (SGEM) solution.

3. Results

The relative errors of the max-min, first- and second-orders of the Hamiltonian, and the global error minimization methods are shown in Tables 2, 3, and 4 for , , and , respectively. It is seen that the second-order global error minimization method results are closer to the exact ones than the other mentioned methods for large initial amplitudes.

The differences of obtained responses and velocities with exact ones are plotted in Figures 1 and 2, respectively. It is found that the second-order global error minimization method illustrates a better accuracy than the aforementioned techniques.

Figure 1 

Difference between obtained responses and the numerical solution (,  ).

Figure 2 

Difference between obtained velocities and the numerical solution (,  ).