This work is licensed under the Creative Commons Attribution 4.0 Public License.
Introduction
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.
Mathematical formulation of the problems
The general form of the differential equation describing the oscillations of single-degree-of-freedom systems can be written as [25]:
$$\begin{array}{}
\displaystyle
\ddot x + F(x) = 0
\end{array}$$
with initial conditions:
$$\begin{array}{}
\displaystyle
x(0) = l \mathrm{\;\; and \;\;}\dot x(0) = 0
\end{array}$$
where F(x) represents the linear and nonlinear terms. A typical nonlinear conservative system that has been the subject of many investigations is Duffing-type oscillator. This nonlinear oscillator represents the dynamic behaviour of many engineering problems. Exact analytical solutions to the oscillatory problem in the form of Eq. (1) are generally impossible, and therefore, some numerical solution is to be reported. The main goal is to apply the discretization and linearization concepts to develop the HPM and CWM which can provide periodic solutions for oscillatory problems.
Some properties of shifted second kind Chebyshev polynomials [26]
In this section, we discuss some relevant properties of the function
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:
The shifted second kind Chebyshev polynomials are defined on the interval [0, 1] by $\begin{array}{}
\displaystyle
U_n^*
\end{array}$(x) = Un(2x – 1). All the results of second kind Chebyshev polynomials can be easily transformed to give the corresponding results for their shifted forms. The orthogonality relation with respect to the weight function $\begin{array}{}
\displaystyle
\sqrt {x - {x^2}}
\end{array}$ is given by
Shifted second kind Chebyshev operational matrix of derivatives (S2KCOM)
Wavelets constitute a family of functions constructed from dilation and translation of a single function called the mother wavelet. When the dilation parameter a and the translation parameter b vary continuously, we have the following family of continuous wavelet:
Shifted second kind Chebyshev wavelets ψnm(t) = ψ(k, n, m, t) have four arguments: k, n can assume any positive integer, m is the order of second kind Chebyshev polynomials and tis the normalized time. They are defined on the interval [0, 1] by
Solution pertaining to the nonlinear differential equations using shifted second kind Chebyshev wavelet method (S2KCWM) and homotopy perturbation method
In this section, we give some numerical results obtained using the algorithms presented in the previous sections. We consider the following examples.
Example 1
We consider the nonlinear differential equation [27]
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 μ ≤ 1. From Figure 1(a–c), it is also observed that the amplitude depends upon the initial conditions. Figure 2 represents the displacement and velocity versus time t. From the figure, it is noted that the amplitude of displacement and velocity are equal.
Figure 3(a–c) denotes the displacement versus time for various values of the parameters α and β. The numerical solution is compared with our analytical results in Figure 3(a–c) and found to be satisfactory. Displacement and velocity are shown in Figure 4. Here also, the amplitude of displacement and velocity are equal. Figure 5 shows the displacement versus time for the oscillation problem (Eq. (45)). Displacement and velocity versus time are shown in Figure 6. All the result also confirmed for the problem Eq. (37).
Conclusion
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.