Homotopy Analysis Method for a Conservative Nonlinear Oscillator with Fractional Power

In this paper, the homotopy analysis method is applied to deduce the periodic solutions of a conservative nonlinear oscillator for which the elastic force term is proportional to u 1/3 . By introducing the auxiliary linear operator and the initial guess of solution, the homotopy analysis solving is set up. By choosing the suitable convergence-control parameter, the accurate high-order approximations of solution and frequency for the whole range of initial amplitudes can easily be obtained. Comparison of the results obtained using this method with those obtained by different methods reveals that the former is more accurate, effective and convenient for these types of nonlinear oscillators.

Among these strongly nonlinear methods, the HAM (the abbreviation for homotopy analysis method) proposed by Liao in 1992 [9] is a general analytical method. The advantage of the HAM as compared to other perturbation methods is that the HAM is independent of the small/large parameters. More importantly, by introducing the non-zero convergence-control parameter h and the nonzero auxiliary function H(t), this method provides a simple way to control and to ensure the convergence of approximation series. Another advantage of the HAM is that one can construct a continuous mapping of the initial guess approximation to the exact solution of the given problem through an auxiliary linear operator. Hence, the HAM attracts a great deal of interest of more and more researchers and has found applications in various nonlinear problems. On the basic ideas and the applications of HAM, one can refer to [9] [10] [11].
In this paper, Liao's HAM is applied to obtain the high-order analytical periodic solutions of a conservative nonlinear oscillator for which the elastic force term is proportional to u 1/3 . This nonlinear oscillator has been recently studied by Beléndez [12] [13] using the HPM (the abbreviation for homotopy perturbation method), and by Lim [14] and Mickens [15] using the HBM (the abbreviation for harmonic balance method). By introducing the auxiliary linear operator and the initial guess of solution, the homotopy analysis solving for this nonlinear oscillator is set up, the high-order approximations of frequency and solution can easily be obtained. Excellent agreement of the approximate frequency with the exact one has been demonstrated and discussed, and we show that the secondorder approximation for frequency in this paper is better than the results by the preceding methods, the maximal relative error of the sixth-order approximate frequency by HAM is less than 0.0045%. As can be seen, the results presented in this paper reveal that the HAM is very effective and convenient for conservative nonlinear oscillators with non-polynomial elastic terms.

Solution Procedure
Considering the following nonlinear oscillator which was introduced as a model "truly nonlinear oscillator" by Mickens [15] [16] [17] where the over-dot denotes differentiation with respect to time t and A is the amplitude of the oscillation. System (1) is a conservative nonlinear oscillator with a fractional power restoring force, it is not amenable to exact treatment and, therefore, the HAM can be applied to solve it. In order to apply the HAM effectively, we rewrite (1) in the following form ( ) with initial conditions where 6 ω Ω = , prime denotes the derivative with respect τ , and ω represents the angular frequency.
According to [12], the periodic solution of (5) does exist, the solution ( ) y τ can be expressed by such a set of base functions where m α are coefficients that should be determined. Considering the initial conditions (6) and the rule of solution expression described by (8), a good initial guess of ( ) y τ can be determined as According to Liao [9], the auxiliary linear operator can be given easily Now, combining (5) with (10), we can construct the zeroth-order deformation equation with the initial conditions where the non-zero constant h is called the convergence-control parameter, is called the homotopy parameter,

( )
H τ is called the non-zero auxiliary function, and is the nonlinear operator. When 0 When 1 p = , the system (11) becomes (5). Thus, by ranging p from 0 p = to 1 p = , it turns out that the initial guess of solution (9) deforms continuously to the exact solution of system (5) provided the following series converge when 1 p = . The series in (14), (15) are called the homotopy series, in which, ( ) 0 y τ is defined in (9). By selecting the convergence-control parameter h suitably, in general, the convergence of the series (14), (15) can be guaranteed.
If the series (14), (15) (14), (15), according to Liao [9], we get the higher-order deformation equation With the perturbation procedure described above going on, the higher-order approximations for

Results and Discussion
In this section, the HAM is applied to obtain the frequency and periodic solution of a conservative nonlinear oscillator for which the elastic force term is proportional to 1 3 u and its accuracy and efficiency are illustrated by comparing the approximate frequency obtained by HAM with the exact one ex ω and other results in the literature.
The exact frequency for system (1) is given by the following expression [12] ( ) ( ) ( ) A Ω becomes independent of h and remains fixed. From Figure 1, according to Liao [9], the derivative of 2 A Ω with respect to h is 0. So, we can choose . Under this condition, the second-order approximations of frequency and solution of system (1) therefore, it can be easily obtained that the relative error of frequency is 0.039%, that is, the second-order approximations give the frequency with the highest accuracy. This indicates that the convergence-control parameter h plays an important role in the HAM.     the HAM provides excellent approximations to the exact periodic solution for the wide range of initial amplitudes for this case study.
In Table 1, we present the comparison between the approximate and exact frequency for the second-order approximation using different methods. It is clear that for the second-order approximation, the result obtained by HAM is better than those obtained previously by other authors.
In order to make the solutions more accurate, we get the higher-order approximations easily by the perturbation procedure described in section 2. Here, the sixth-order approximations of frequency and solution can be expressed in the following