Recent Advances on Methods and Applications of Nonlinear Differential Equations
View this Special IssueResearch Article  Open Access
Q. X. Liu, Y. M. Chen, J. K. Liu, "An Iteration Scheme Suitable for Solving Limit Cycles of Nonsmooth Dynamical Systems", Mathematical Problems in Engineering, vol. 2013, Article ID 582865, 5 pages, 2013. https://doi.org/10.1155/2013/582865
An Iteration Scheme Suitable for Solving Limit Cycles of Nonsmooth Dynamical Systems
Abstract
The Mickens iteration method (MIM) is modified to solve selfexcited systems containing nonsmooth nonlinearities and/or nonlinear damping terms. If the MIM is implemented routinely, the unknown frequency and amplitude of limit cycle (LC) would couple to each other in complicated nonlinear algebraic equations at each iteration. It is cumbersome to solve these algebraic equations, especially for nonsmooth systems. In the modified procedures, the unknown frequency is substituted by the determined value obtained at the previous iteration. By this means, the frequency is decoupled from the nonlinear terms. Numerical examples show that the LCs obtained by the modified MIM agree well with numerical results. The presented method is very suitable for solving selfexcited systems, especially those with nonlinear damping and nonsmooth nonlinearities.
1. Introduction
Recent years have witnessed the wide applications of iteration techniques, such as the Mickens iteration method (MIM) [1–6] and the variational iteration method [7–10]. In order to improve the efficiency, Lim and Wu [11], Marinca and Herisanu [12], and Hu [13] modified the MIM, respectively.
In principle, the approximations can be obtained to any desired accuracy by the MIM as long as the iteration proceeds. In the MIM, algebraic equations are introduced at each iteration to eliminate the socalled secular terms. In the applications in conservative oscillators, the algebraic equations are linear. Nonlinear damping and nonsmooth terms appear widely in many dynamical systems [14, 15]. As for the oscillators with nonlinear damping terms, however, very complicated nonlinear algebraic equations have to be solved at each iteration [16]. Moreover, the algebraic equations cannot be deduced for systems with nonsmooth nonlinearities. It is necessary and worthwhile, therefore, to propose some approaches to simplify the MIM. This paper will present a modified iteration algorithm by decoupling the unknown frequency from nonlinear terms.
2. A Modified MIM
Consider a selfexcited oscillator where the superscript denotes the differentiation with respect to time and is a nonlinear term with damping terms. Assume that system (1) has at least one limit cycle (LC) solution. Since the LC frequency and amplitude are independent of initial conditions, they should be considered as unknowns to be determined at every iteration. Denote the angular frequency as and introduce the transformation as ; thus, we rewrite (1) as subject to the following initial conditions: where is the unknown LC amplitude and the superscript denotes the differentiation with respect to . Note that will be approximated as a series by eliminating the secular terms at each iteration stage. In order to obtain the LC, the MIM [16] can be given as with the initial conditions being rewritten at each iteration as Note that the coefficient of the first harmonic in remains still to be an unknown, that is, . This unknown will couple with the unknown frequency, , which will result in a coupled nonlinear term (i.e., ) in the right side of (4). If higher powers of exist, these terms will lead to very complicated functions in and . Different from conservative systems, can no longer be considered as an independent unknown. In order to simplify the MIM, therefore, a modified scheme is proposed as As , we choose . In the th iteration, is a given constant that is obtained at the ()th iteration. The square of the unknown frequency, that is, , can be treated as an independent parameter, because appears only in . According to the initial conditions, the starting iteration solution can be chosen as It is obvious that as long as the series and are convergent, they must converge to the exact solutions. The righthand side of (6) can be expressed by Fourier series as where the harmonic coefficients and are functions in and . Here, is a positive integer denoting the order of the highest harmonic. Approximations and are determined by eliminating the secular terms, that is, letting These equations can be solved analytically if is considered as an independent unknown. They can also be numerically solved by NewtonRaphson method. The latter is employed in this study.
Different from the existing procedures [16], as increases, is always an independent unknown in the modified MIM. Moreover, unknown does not couple with nonlinear terms. It simplifies the MIM to a large extent, as shown later.
3. Numerical Examples
Example 1 (system with nonlinear damping terms). The van der Pol equation is chosen to illustrate the previous procedures more clearly:
where is a given constant. As known, (10) has a stable LC solution when while an unstable one when .
According to the modified MIM, the corresponding iteration scheme is given as
Introducing a new time variable at each iteration stage, we rewrite (11) as
where the superscript denotes the derivative with respect to . The iteration algorithm begins with an initial solution
Then, we obtain the governing equations in as
Equating the coefficients of and to zeros results into
which yields that and . Substituting them into (16), we have
Considering initial conditions (14), we can obtain
where is to be determined at the next iteration stage.
According to iterative scheme (11), the equation in is deduced as
Equate the coefficients of and to zeros:
By solving (19) numerically, we can determine and . Here, we obtain the secondorder approximation and expand it as
According to [17], the LindstedtPoincare (LP) method provides the second and forthorder approximate frequency and , respectively. The attained approximation agrees well with the 4thorder LP solution. Table 1 indicates that is more accurate than when compared with .

Figure 1 shows the comparison of the phase planes between iteration solutions and numerical result. Rapid convergence of to the numerical result can be observed. Note that all numerical solutions are obtained by the fourthorder RungeKutta (RK) integration method. When , the iteration procedure presented by Chen and Liu [16] does not converge. This is probably the difference between the starting function () and the exact solution is too large. The modified MIM is still effective for . As Figure 2 shows, the LC solution with obtained by the presented method is in excellent agreement with numerical one. It is necessary to point out that the presented method is able to track unstable LCs, whereas the RK method is not.
Also plotted in Figure 2 are the results provided by the LP method [17]. The iteration results are much more precise than the 2nd and 4thorder LP approximations.
It is necessary to point out that the presented method is able to track unstable LCs, whereas the RK method is not. Figure 3 shows an unstable LC of the van der Pol equation with obtained by the presented method. As shown, the RK begins at the LC; however, the solution curve converges to the equilibrium.
In order to further demonstrate the merit of the modified MIM when applied to problems with nonlinear damping terms, we consider the following selfexcited system [18]: The nonlinear term contains high powers of , that is, . If the original MIM is employed, the algebraic equations governing will become very complicated. Therefore, it is necessary to employ the modified approach. Figure 4 indicates that the approximations obtained by the presented method converge rapidly to the numerical solution as increased.
Example 2 (system with nonsmooth nonlinearity). The modified MVIM is further applied to nonsmooth dynamical system expressed as
Here, is a nonlinear damping term, and is a nonsmooth function. If substituting the into , on account of contained unknown quantities , so can not be expanded as Fourier progression by numerical integration. To this end, (9) cannot be deduced by eliminating the secular terms. Likewise, we present the following iteration scheme:
In this scheme, can be expanded as a Fourier series since has been determined at the previous iteration.
Let us consider a van der Pol type oscillator with a nonsmooth function as
with
Figure 5 shows the LC of system (14) with and . The 5thorder approximations obtained by the presented method agree well with the numerical solution when the nonsmooth term is given as
The LC can also be obtained very accurate, as Figure 6 shows.
4. Conclusions
The Mickens iteration method (MIM) has been modified, so that it is suitable for solving LC solutions of selfexcited systems with nonsmooth and/or damping nonlinearities. Different from the routinelyused MIM, the modified method decouples the unknown frequency from nonlinear terms. This modification simplifies the MIM significantly. Numerical examples show the feasibility and validity of the presented method, which implies that it could be applicable to more nonlinear dynamical systems, especially those with nonlinear damping terms and nonsmooth nonlinearities.
Acknowledgment
This work is supported by the National Natural Science Foundation of China (11002088, 11272361, 11172333), Doctoral Program Foundation of Ministry of Education of China (20130171110039), Guangdong Province Natural Science Foundation (S2012040007920, S2013010013802), Fundamental Research Funds for the Central Universities (13lgzd06), and the Guangdong Province Science and Technology Program (2012A030200011).
References
 R. E. Mickens, “Iteration procedure for determining approximate solutions to nonlinear oscillator equations,” Journal of Sound and Vibration, vol. 116, no. 1, pp. 185–187, 1987. View at: Publisher Site  Google Scholar  MathSciNet
 R. E. Mickens, “Harmonic balance and iteration calculations of periodic solutions to ${y}^{\u2033}+{y}^{1}=0$,” Journal of Sound and Vibration, vol. 306, no. 3–5, pp. 968–972, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 R. E. Mickens, “A generalized iteration procedure for calculating approximations to periodic solutions of “truly nonlinear oscillators”,” Journal of Sound and Vibration, vol. 287, no. 45, pp. 1045–1051, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. Bhattacharjee and J. K. Bhattacharjee, “Lindstedt Poincare technique applied to molecular potentials,” Journal of Mathematical Chemistry, vol. 50, no. 6, pp. 1398–1410, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 H. Hu and J. H. Tang, “A classical iteration procedure valid for certain strongly nonlinear oscillators,” Journal of Sound and Vibration, vol. 299, no. 12, pp. 397–402, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. I. Ramos, “On LinstedtPoincaré technique for the quintic Duffing equation,” Applied Mathematics and Computation, vol. 193, no. 2, pp. 303–310, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 F. K. Yin, J. Q. Song, and X. Q. Cao, “Couple of the variational iteration method and Legendre wavelets for nonlinear partial differential equations,” Journal of Applied Mathematics, vol. 2013, Article ID 157956, 11 pages, 2013. View at: Publisher Site  Google Scholar
 M. T. Atay and O. Kilic, “The semianalytical solutions for stiff systems of ordinary differential equations by using variational iteration method and modified variational iteration method with comparison to exact solutions,” Mathematical Problems in Engineering, vol. 2013, Article ID 143915, 11 pages, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 A.J. Chen, “Resonance analysis for tilted support spring coupled nonlinear packaging system applying variational iteration method,” Mathematical Problems in Engineering, vol. 2013, Article ID 384251, 4 pages, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 V. Marinca, N. Herişanu, and C. Bota, “Application of the variational iteration method to some nonlinear onedimensional oscillations,” Meccanica, vol. 43, no. 1, pp. 75–79, 2008. View at: Publisher Site  Google Scholar
 C. W. Lim and B. S. Wu, “A modified Mickens procedure for certain nonlinear oscillators,” Journal of Sound and Vibration, vol. 257, no. 1, pp. 202–206, 2002. View at: Publisher Site  Google Scholar  MathSciNet
 V. Marinca and N. Herisanu, “A modified iteration perturbation method for some nonlinear oscillation problems,” Acta Mechanica, vol. 184, no. 1–4, pp. 231–242, 2006. View at: Publisher Site  Google Scholar
 H. Hu, “Solutions of the Duffingharmonic oscillator by an iteration procedure,” Journal of Sound and Vibration, vol. 298, no. 12, pp. 446–452, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Y. M. Chen, G. Meng, and J. K. Liu, “A new method for Fourier series expansions: applications in rotorseal systems,” Mechanics Research Communications, vol. 38, no. 5, pp. 399–403, 2011. View at: Publisher Site  Google Scholar
 T. Pirbodaghi, M. T. Ahmadian, and M. Fesanghary, “On the homotopy analysis method for nonlinear vibration of beams,” Mechanics Research Communications, vol. 36, no. 2, pp. 143–148, 2009. View at: Publisher Site  Google Scholar
 Y. M. Chen and J. K. Liu, “A modified Mickens iteration procedure for nonlinear oscillators,” Journal of Sound and Vibration, vol. 314, no. 3–5, pp. 465–473, 2008. View at: Publisher Site  Google Scholar
 A. H. Nayfeh, Introduction to Perturbation Techniques, WileyInterscience, New York, NY, USA, 1993. View at: MathSciNet
 W. J. Ding, SelfExcited Vibration, Tsinghua University Press, Beijing, China, 2009.
Copyright
Copyright © 2013 Q. X. Liu 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.