Recent Development in Fixed-Point Theory, Optimization, and their ApplicationsView this Special Issue
Research Article | Open Access
Application of Extended Homotopy Analysis Method to the Two-Degree-of-Freedom Coupled van der Pol-Duffing Oscillator
The extended homotopy analysis method (EHAM) is presented to establish the analytical approximate solutions for two-degree-of-freedom (2-DOF) coupled van der Pol-Duffing oscillator. Meanwhile, the comparisons between the results of the EHAM and standard Runge-Kutta numerical method are also presented. The results demonstrate that the analytical approximate solutions of the EHAM agree well with the numerical integration solutions. For EHAM as an analytical approximation method, we are not sure whether it can apply to all of the nonlinear systems; we can only verify its effectiveness through specific cases. As a result of the existence of nonlinear terms, we must study different types of systems, no matter from the complication of calculation and physical significance.
Mathematical methods for the natural and engineering sciences problems have drawn considerable attention in recent years. In normal circumstances, most of the nonlinear dynamical models can be governed by a set of differential equations and auxiliary conditions from modeling processes . Numerous analytical methods have been developed to deal with the nonlinear differential equations, such as the modified perturbation methods [2–5], improved harmonic balance methods [6, 7], energy balance method [8, 9], and the frequency-amplitude formulation [10, 11]. Enlightening from the basic concepts of the homotopy in topology [12, 13], Liao developed the homotopy analysis method (HAM) [14–16] which does not require small parameters as one of the efficient analytical techniques in solving a variety of nonlinear vibration problems.
Recently, great attention is paid to the discussion of coupled oscillators of nonlinear dynamical systems because most of practical engineering problems can be governed by such coupled systems [17–19]. The extended homotopy analysis method (EHAM) is one method based on the HAM envisioned first by Liao . More recently, Qian et al.  extended the HAM to deal with strongly nonlinear coupled van der Pol oscillators. For EHAM as an analytical approximation method, we are not sure whether it can apply to all of the nonlinear systems; we can only verify its effectiveness through specific cases. As a result of the existence of nonlinear terms, we must study different types of systems, no matter from the complication of calculation and physical significance. By solving such example, it is illustrated that the present techniques are not an adhoc approach; it can be generalized to investigate more complicated nonlinear multi-degree-of-freedom (MDOF) dynamical systems.
In the present work, the exact analytical series solutions of the two-degree-of-freedom (2-DOF) coupled van der Pol-Duffing system are obtained by using the EHAM, and we also establish the comparisons between the results of the EHAM and standard Runge-Kutta numerical method. It is shown that the periodic solutions of the EHAM are in excellent agreement with the numerical integration ones, even if time progresses to a certain large domain. In what follows, Section 2 presents the EHAM of the MDOF dynamical system. Moreover, the EHAM is presented to establish the analytical approximate solutions for 2-DOF coupled van der Pol-Duffing oscillator in the next section. In Section 4, numerical comparisons are carried out to authenticate the correctness and accuracy of the present method. Finally, the paper ends with concluding remarks in Section 5.
2. The Extended Homotopy Analysis Method
The MDOF dynamical system is considered by the following equation: where is an -dimensional unknown vector, a dot denotes the derivative with respect to time , , , and are, respectively, mass, damping, and stiffness matrixes, and is the vector function of , , and . Let ; then (1) is an autonomous dynamical system.
From (1), we define a nonlinear operator as where is an unknown vector value function and and are spatial and temporal variables, respectively.
In (2), the unknown vector functions of , , and are, respectively,
According to the fundamental concepts and working procedures of the HAM [1, 2], the zeroth-order deformation equation can be constructed as follows: where is an embedding parameter, is the solution of initial guess, is an auxiliary linear operator, and and are the auxiliary parameters and the functions, respectively.
The operator has the following property: When and , the zeroth-order deformation equation (4) is and , respectively. Hence, as increases from 0 to 1, the solution varies from the initial guess solution to the exact solution . In this paper, where with the initial conditions
Setting and expanding into the Taylor series expansion with respect to in accordance with the theorem of vector-valued function, we obtain
If the auxiliary linear operator, initial guess solution, auxiliary parameters , and auxiliary functions are properly chosen, the series equation (10) converges at , and we arrive at For brevity, the vector of is defined as
Differentiating the zeroth-order deformation equation (4) times with respect to then dividing the equation by and setting yield where
The th-order deformation equation (13) is a linear equation, which can be readily solved by the symbolic software such as Mathematica.
3. Application of the EHAM
In this section, we apply the EHAM for analysis of the two coupled van der Pol-Duffing oscillators:where the superscript denotes the differentiation with respect to time , and are the unknown real functions, and , , , , , , and are parameters.
We introduce a new variable and substitute , , and into (15a) and (15b). Therefore, we havesubject to the initial conditions where a prime denotes the derivative with respect to variable . Provided that the periodic solutions in (16a) and (16b) can be expressed by a set of base functions one obtains
For the initial approximation, and are assumed as and the linear operator is defined as
We can define a nonlinear operator as the following by EHAM:and the nonlinear operator is
In terms of the principle of solution expression, we select the auxiliary functions as and ; thus the zeroth-order deformation equation is given by where with the initial conditions
Obviously, as the embedding parameter varies from 0 to 1, changes from the initial guess to the exact solutions . In addition, changes from the initial guess frequency to the nonlinear physical frequency .
With the help of the Taylor series expansion and (13), we obtainwhere
For simplicity, the following vectors are defined as
By differentiating the zeroth-order deformation equation (24) times with respect to , then dividing the equation by , and setting , the th-order deformation equation is formulated as follows: with the initial conditions in which
Because of the principle of solution expression and the linear operator , the right side of (33) should not contain the terms of and or the secular terms and . The coefficients are set to be zero to yield
The solutions of , , , and () from (33) and (36a), (36b), (36c), and (36d) can be computed successively. To achieve more accurate results, we modify the solution of as follows: where is a small parameter.
4. Numerical Simulation and Discussion
In this section, numerical experiment is conducted to verify the accuracy of the present approach.
Taking , , , , , , , and the initial approximations of , , , and are , , , and , respectively.
For simplicity and accuracy, we set , , and ; then the comparison of the phase curves of the fifth-order approximation with the numerical integration solution is shown in Figure 1.
The initial conditions of the numerical integration method are , , , and .
Moreover, the fifth-order analytical solutions (, , and ) are given as in which
The above results demonstrate that the system has an in-phase solution. While giving the initial approximations of , , , and , we can get an out-of-phase solution to the system.
In this case, the comparison between the phase curves of the fifth-order approximation and the numerical integration solution is portrayed in Figure 2.
The initial conditions of the numerical integration method are , , , and , and the fifth-order analytical solutions are written as in which
In the present paper, the EHAM approach is applied to get asymptotic analytical series solutions of 2-DOF van der Pol-Duffing oscillators with a nonlinear coupling. The basic idea described in this paper is expected to be more employed in solving other dynamical systems in engineering and physical sciences. Comparisons with the numerical results are presented to demonstrate the validity of this method. In summary, compared with some other methods, the EHAM has the following advantages.(1)The EHAM provides an ingenious avenue for controlling the convergences of approximation series. Numerical comparisons demonstrate that the EHAM is an effective and robust analytical method of 2-DOF van der Pol-Duffing oscillators.(2)Because of its flexibility, the present techniques can also be further generalized to analyze more complicated nonlinear MDOF dynamical systems that can only be analyzed by numerical methods.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
All the authors contributed equally and significantly to the writing of this paper. All the authors read and approved the final paper.
The author Y. H. Qian gratefully acknowledges the support of the National Natural Science Foundations of China (NNSFC) through Grants nos. 11202189 and 11304286 and the Natural Science Foundation of Zhejiang Province of China through Grant no. LY12A02002. The author S. M. Chen gratefully appreciates the financial support from the NNSFC through Grants no. 11371326. The author L. Shen gratefully acknowledges the support from open experiment project of Zhejiang Normal University. The authors are also grateful to the anonymous reviewers for their constructive comments and suggestions.
- R. E. Mickens, Mathematical Methods for the Natural and Engineering Sciences, World Scientific, Singapore, 2004.
- Y. K. Cheung, S. H. Chen, and S. L. Lau, “A modified Lindstedt-Poincaré method for certain strongly nonlinear oscillators,” International Journal of Non-Linear Mechanics, vol. 26, no. 3-4, pp. 367–378, 1991.
- M. Senator and C. N. Bapat, “A perturbation technique that works even when the nonlinearity is not small,” Journal of Sound and Vibration, vol. 164, no. 1, pp. 1–27, 1993.
- P. Amore and A. Aranda, “Improved Lindstedt-Poincaré method for the solution of nonlinear problems,” Journal of Sound and Vibration, vol. 283, no. 3–5, pp. 1115–1136, 2005.
- R. R. Pušenjak, “Extended Lindstedt-Poincare method for non-stationary resonances of dynamical systems with cubic nonlinearities,” Journal of Sound and Vibration, vol. 314, no. 1-2, pp. 194–216, 2008.
- J. L. Summers and M. D. Savage, “Two timescale harmonic balance. I. Application to autonomous one-dimensional nonlinear oscillators,” Proceedings of the Royal Society of London A. Mathematical, Physical and Engineering Sciences, vol. 340, no. 1659, pp. 473–501, 1992.
- B. Wu and P. Li, “A method for obtaining approximate analytic periods for a class of nonlinear oscillators,” Meccanica, vol. 36, no. 2, pp. 167–176, 2001.
- S. S. Ganji, D. D. Ganji, Z. Z. Ganji, and S. Karimpour, “Periodic solution for strongly nonlinear vibration systems by He's energy balance method,” Acta Applicandae Mathematicae, vol. 106, no. 1, pp. 79–92, 2009.
- I. Mehdipour, D. D. Ganji, and M. Mozaffari, “Application of the energy balance method to nonlinear vibrating equations,” Current Applied Physics, vol. 10, no. 1, pp. 104–112, 2010.
- L. Zhao, “He's frequency-amplitude formulation for nonlinear oscillators with an irrational force,” Computers & Mathematics with Applications, vol. 58, no. 11-12, pp. 2477–2479, 2009.
- J. Fan, “He's frequency-amplitude formulation for the Duffing harmonic oscillator,” Computers & Mathematics with Applications, vol. 58, no. 11-12, pp. 2473–2476, 2009.
- P. J. Hilton, An Introduction to Homotopy Theory, Cambridge University Press, Cambridge, UK, 1953.
- C. Nash and S. Sen, Topology and Geometry for Physicists, Academic Press, London, UK, 1983.
- S. J. Liao, The proposed homotopy analysis techniques for the solution of nonlinear problems [Ph.D. thesis], Shanghai Jiao Tong University, Shanghai, China, 1992.
- S. J. Liao, “A kind of approximate solution technique which does not depend upon small parameters. II. An application in fluid mechanics,” International Journal of Non-Linear Mechanics, vol. 32, no. 5, pp. 815–822, 1997.
- S. Liao, Beyond Perturbation: Introduction to the Homotopy Analysis Method, Chapman & Hall, Boca Raton, Fla, USA, 2004.
- S. Rajasekar and K. Murali, “Resonance behaviour and jump phenomenon in a two coupled Duffing-van der Pol oscillators,” Chaos, Solitons and Fractals, vol. 19, no. 4, pp. 925–934, 2004.
- B. T. Nohara and A. Arimoto, “Non-existence theorem except the out-of-phase and in-phase solutions in the coupled van der Pol equation system,” Ukrainian Mathematical Journal, vol. 61, no. 8, pp. 1311–1337, 2009.
- Y. J. Li, B. T. Nohara, and S. J. Liao, “Series solutions of coupled van der Pol equation by means of homotopy analysis method,” Journal of Mathematical Physics, vol. 51, no. 6, pp. 1–12, 2010.
- Y. H. Qian, C. M. Duan, S. M. Chen, and S. P. Chen, “Asymptotic analytical solutions of the two-degree-of-freedom strongly nonlinear van der Pol oscillators with cubic couple terms using extended homotopy analysis method,” Acta Mechanica, vol. 223, no. 2, pp. 237–255, 2012.
Copyright © 2014 Y. H. Qian et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.