0
Research Papers

# An Approximate Solution for Period-1 Motions in a Periodically Forced Van Der Pol Oscillator

[+] Author and Article Information
Albert C. J. Luo

Department of Mechanical and
Industrial Engineering,
Southern Illinois University Edwardsville,
Edwardsville, IL 62026-1805
e-mail: aluo@siue.edu

Arash Baghaei Lakeh

Department of Mechanical and
Industrial Engineering,
Southern Illinois University Edwardsville,
Edwardsville, IL 62026-1805

1Corresponding author.

Contributed by the Design Engineering Division of ASME for publication in the JOURNAL OF COMPUTATIONAL AND NONLINEAR DYNAMICS. Manuscript received January 27, 2013; final manuscript received December 31, 2013; published online February 13, 2014. Assoc. Editor: Claude-Henri Lamarque.

J. Comput. Nonlinear Dynam 9(3), 031001 (Feb 13, 2014) (7 pages) Paper No: CND-13-1021; doi: 10.1115/1.4026425 History: Received January 27, 2013; Revised December 31, 2013

## Abstract

In this paper the approximate analytical solutions of period-1 motion in the periodically forced van der Pol oscillator are obtained by the generalized harmonic balance (HB) method. Such an approximate solution of periodic motion is given by the Fourier series expression, and the convergence of such an expression is guaranteed by the Fourier series theory of periodic functions. The approximate solution is different from traditional, approximate solution because the number of total harmonic terms ($N$) is determined by the precision of harmonic amplitude quantity level, set by the investigator (e.g., $AN≤ɛ$ and $ɛ=10-8$). The stability and bifurcation analysis of the period-1 solutions is completed through the eigenvalue analysis of the coefficient dynamical systems of the Fourier series expressions of periodic solutions, and numerical illustrations of period-1 motions are compared to verify the analytical solutions of periodic motions. The trajectories and analytical harmonic amplitude spectrum for stable and unstable periodic motions are presented. The harmonic amplitude spectrum shows the harmonic term effects on periodic motions, and one can directly know which harmonic terms contribute on periodic motions and the convergence of the Fourier series expression is clearly illustrated.

<>

## References

van der Pol, B., 1920, “A Theory of the Amplitude of Free and Forced Triode Vibrations,” Radio Rev., 1, pp. 701–710.
van der Pol, B., and van der Mark, J., 1927, “Frequency Demultiplication,” Nature, 120 pp. 363–364.
Lagrange, J. L., 1788, Mecanique Analytique, Albert Balnchard, Paris, Vol. 2.
Fatou, P., 1928, “Sure′ le Mouvement d’un Systeme Soumis a′ des Forces a Courte Periode,” Bull. Soc. Math. France, 56, pp. 98–139.
Krylov, N. M., and Bogolyubov, N. N., 1935, “Methodes Approchoes de la Mecanique Non Lineaire Dans Leur Application a l’6tude de la Perturbation des Mouvements Periodiques de Divers Phenomenes de Resonance s’y Rapportant,” Academie des Sciences d’Ukraine, Kiev.
Bogoliubov, N., and Mitropolsky, Y., 1961, Asymptotic Methods in the Theory of Nonlinear Oscillations, Gordon and Breach, New York.
Hayashi, C.1964, Nonlinear Oscillations in Physical Systems, McGraw-Hill, New York, Vol. 33.
Cap, F., 1974, “Averaging Method for the Solution of Nonlinear Differential Equations With Periodic Non-Harmonic Solutions,” Int. J. Non-Linear Mech., 9, pp. 441–450.
Coppola, V., and Rand, R., 1990, “Averaging Using Elliptic Functions: Approximation of Limit Cycles,” Acta Mech., 81, pp. 125–142.
Rand, R., and Armbruster, D., 1987, Perturbation Methods, Bifurcation Theory, and Computer Algebra, Springer, New York.
Garcia-Margallo, J., and Bejarano, J. D., 1987, “A Generalization of the Method of Harmonic Balance,” J. Sound Vib., 116, pp. 591–595.
Xu, Z., and Cheung, Y., 1994, “Averaging Method Using Generalized Harmonic Functions for Strongly Nonlinear Oscillators,” J. Sound Vib., 174, pp. 563–576.
Buonomo, A., 1998, “On the Periodic Solution of the van der Pol Equation for the Small Damping Parameter,” Int. J. Circuit Theory Appl., 26, pp. 39–52.
Buonomo, A., 1998, “The Periodic Solution of van der Pol’s Equation,” SIAM J. Appl. Math., 59, pp. 156–171.
Waluya, S. B., and Van Horssen, W. T., 2003, “On Approximation of First Integrals for a Strongly Nonlinear Forced Oscillator,” Nonlinear Dyn., 33, pp. 225–252.
Luo, A. C. J., and Huang, J., 2011, “Approximate Solutions of Periodic Motions in Nonlinear Systems via a Generalized Harmonic Balance,” J. Vib. Control, 18, pp. 1661–1671.
Luo, A. C. J., and Huang, J.2012, “Analytical Dynamics of Period-m Flows and Chaos in Nonlinear Systems,” Int. J. Bifurcation Chaos, 22(4), p. 1250093.
Luo, A. C. J., 2012, Continuous Dynamical Systems, HEP-L&H Scientific, Beijing/Glen Carbon.

## Figures

Fig. 1

Frequency–amplitude responses of period-1 motion for van der Pol oscillator from the 22 harmonic terms (HB22): (a)–(d) A1-A11 (α1 = 0.5,α2 = 5,α3 = 10,Q0 = 0.5,1.0,2.0,4.0,8.0). Acronym “HB” represents Hopf bifurcation.

Fig. 3

Analytical and numerical solutions of limit cycle with 22 harmonic terms (HB22) with Ω≈3.1574 and Q0 = 0: (a) phase plane and (b) analytical harmonic amplitude spectrum. Solid and symbol curves are for numerical and analytical solutions. The initial condition is x0 = 0.0844137, and y0 = 2.02906 (α1 = 0.5,α2 = 5,α3 = 10).

Fig. 4

Analytical and numerical solutions of stable period-1 motion with 22 harmonic terms (HB22) at Ω = 1.048 and Q0 = 4.0: (a) phase plane and (b) analytical harmonic amplitude spectrum. Solid and symbol curves are for numerical and analytical solutions. The initial condition is x0 = 0.391628, and y0 =  -0.556119 (α1 = 0.5,α2 = 5,α3 = 10).

Fig. 5

Analytical and numerical solutions of unstable period-1 motion with 22 harmonic terms (HB22) at Ω = 0.5 and Q0 = 4: (a) phase plane, (b) displacement, (c) velocity, and (d) analytical harmonic amplitude spectrum. Solid and symbol curves are for numerical and analytical solutions. The initial condition is x0 = 0.409654, and y0 = 6.472250×10-3 (α1 = 0.5,α2 = 5,α3 = 10).

Fig. 2

Parameter map for the period-1 motion of the periodically forced van der Pol oscillator. The gray and white areas are for stable and unstable period-1 motions, respectively (α1 = 0.5,α2 = 5,α3 = 10).

Fig. 6

Analytical and numerical solutions of unstable period-1 motion with 22 harmonic terms (HB22) at Ω = 8 and Q0 = 4: (a) phase plane and (b) analytical harmonic amplitude spectrum. Solid and symbol curves are for numerical and analytical solutions. The initial condition is x0 = -0.073682, and y0 =  -0.042901 (α1 = 0.5,α2 = 5,α3 = 10).

## Discussions

Some tools below are only available to our subscribers or users with an online account.

### Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related Proceedings Articles
Related eBook Content
Topic Collections