Most of the oscillation problems in engineering sciences are nonlinear, and it is difficult to solve such equations analytically. Recently, nonlinear oscillator models have been widely considered in physical and chemical sciences. Due to the limitation of existing exact solutions, many approximate analytical and numerical approaches have been investigated. Many real-life problems that arise in several branches of pure and applied science can be expressed using the nonlinear differential equations. Therefore, these nonlinear equations must be solved using analytical/numerical methods. Many researchers have been working on various analytical methods for solving nonlinear oscillation systems in the last decades. Nowadays, the computational experience is significant, and several numerical methods have been suggested and analyzed under certain conditions. These numerical methods have been developed using different techniques such as Taylor series, homotopy perturbation method, quadrature formula, variation iteration method and decomposition method [1, 2, 3, 4, 5, 6, 7]. Noor et al. [8] have applied a sixth-order predictor–corrector iterative method for solving the nonlinear equations. Chun [9] proposed a class of fifth-order and sixth-order iterative methods. In this paper, we apply the wavelet transform method and homotopy perturbation method for solving the nonlinear oscillation equations. The homotopy perturbation method was developed by He [10, 11] and has been applied to a wide class of nonlinear and linear problems arising in various branches.
Recently, the homotopy perturbation method (HPM) [12, 13, 14, 15, 16] with expanding parameter is applied to solve some of the nonlinear equations. Ren and Hu [17], Yu et al. [18] and Wu and He [19] developed the nonlinear oscillation problem using HPM. Abbasbandy et al. [20] developed the homotopy perturbation method and the decomposition method [21] for solving nonlinear equations. Ganji et al. [22] applied the oscillation systems with nonlinearity terms such as the motion of a rigid rod rocking. Khudayarov and Turaev [23] developed the mathematical model of the problem of nonlinear oscillations of a viscoelastic pipeline conveying fluid. Nasab et al. [24] solved the nonlinear singular boundary value problems.
To our knowledge, no exact analytical expressions for oscillation problems are reported. However, it is difficult for us to obtain the exact solution for these problems. The purpose of this paper was to derive the approximate analytical expressions for some oscillation problems in engineering sciences.
The general form of the differential equation describing the oscillations of single-degree-of-freedom systems can be written as [25]:
with initial conditions:
where
In this section, we discuss some relevant properties of the function
These polynomials are orthogonal on [–1, 1]
The following properties of second kind Chebyshev polynomials [26] are of fundamental importance in the sequel. They are eigenfunctions of the following singular Sturm–Liouville equation:
where
starting from
The shifted second kind Chebyshev polynomials are defined on the interval [0, 1] by
The first derivative
The first derivative of the shifted second kind Chebyshev polynomial is given by
Wavelets constitute a family of functions constructed from dilation and translation of a single function called the mother wavelet. When the dilation parameter
Shifted second kind Chebyshev wavelets
A function
where
and
where C and
In this section, we give some numerical results obtained using the algorithms presented in the previous sections. We consider the following examples.
We consider the nonlinear differential equation [27]
with initial conditions:
Eq. (17) is solved using the wavelet scheme,
Moreover, the initial conditions are
Using Eq. (19) and Eq. (20), the following system of algebraic equations can be obtained
The suggested method of solution
Solving the nonlinear Eq. (17) using a new approach of homotopy perturbation method (Appendix A)
The velocity becomes
We consider equation [28]
with the initial conditions
The wavelet scheme of Eq. (27) is
The initial conditions are
Using Eq. (29) and Eq. (30), the following system of algebraic equations can be obtained
The solution of
Solving the nonlinear Eq. (27) using a new approach of homotopy perturbation method, we get
The velocity becomes
We consider the differential equation [29]
with initial conditions:
Eq. (37) is solved using the wavelet scheme,
Moreover, the initial conditions are
Using Eq. (39) and Eq. (40), the following system of algebraic equations can be obtained
The solution of
Solving the nonlinear Eq. (37) using a new approach of homotopy perturbation method, we get
The velocity becomes
To illustrate the applicability, accuracy and effectiveness of the proposed method, we have compared the approximate analytical solution of the nonlinear differential equations with numerical data. The function ode 45 (Runge–Kutta method) in MATLAB software, which is a function of solving the initial value problems, is used. In Figure 1(a–c), we have plotted the numerical solution and the approximate solution derived by our proposed method using HPM and CWM. The figure shows the behaviour of the solution for various values of the parameter. From the figure, it is observed that the variation in the approximate solution is small, when
(a–c) Comparison of CWM (Eq.(24), HPM (Eq.(25) and numerical method (MATLAB result) for various parameter values. Fig.1(a) Fig. 1
Plot of displacement and velocity for oscillator Eq. (26) with weak nonlinearity and small amplitude oscillations Fig. 2
Figure 3(a–c) denotes the displacement versus time for various values of the parameters
Comparison of CWM (Eq. (34), HPM (Eq. (35) and numerical method (MATLAB result) for various parameter values. Fig. 3(a) Fig. 3
Plot of displacement and velocity for oscillator Eq. (36) with weak nonlinearity and small amplitude oscillations Fig. 4
Comparison of CWM (Eq. (44), HPM (Eq. (45) and numerical method (MATLAB result) for fixed parameter values Fig. 5
Plot of displacement and velocity for oscillator Eq. (46) with weak nonlinearity and small amplitude oscillations Fig. 6
In this paper, a wavelet technique has been employed for the approximate solution successfully to solve the well-known nonlinear oscillator differential equations such as Duffing equation with different parameters. There is no need for iterations for achieving sufficient accuracy in numerical results. The results are also obtained via CWM, HPM and numerical solution. Moreover, the proposed method is used to compare CWM, HPM, and NM iteration with the nonlinear part. The effects of constant parameters on responses of the system for approximate solution are also shown in figures.