Review Article

Mathieu's Equation and Its Generalizations: Overview of Stability Charts and Their Features

[+] Author and Article Information
Ivana Kovacic

Centre of Excellence for Vibro-Acoustic
Systems and Signal Processing,
Faculty of Technical Sciences,
University of Novi Sad,
Novi Sad 21215, Serbia

Richard Rand

Department of Mathematics,
Cornell University,
Ithaca, NY 14853;
Department of Mechanical and
Aerospace Engineering,
Cornell University,
Ithaca, NY 14853

Si Mohamed Sah

Technical University of Denmark,
Department of Mechanical Engineering,
Section of Solid Mechanics,
Lyngby 2800, Denmark

Manuscript received August 7, 2017; final manuscript received January 18, 2018; published online February 14, 2018. Editor: Harry Dankowicz.

Appl. Mech. Rev 70(2), 020802 (Feb 14, 2018) (22 pages) Paper No: AMR-17-1054; doi: 10.1115/1.4039144 History: Received August 07, 2017; Revised January 18, 2018

This work is concerned with Mathieu's equation—a classical differential equation, which has the form of a linear second-order ordinary differential equation (ODE) with Cosine-type periodic forcing of the stiffness coefficient, and its different generalizations/extensions. These extensions include: the effects of linear viscous damping, geometric nonlinearity, damping nonlinearity, fractional derivative terms, delay terms, quasiperiodic excitation, or elliptic-type excitation. The aim is to provide a systematic overview of the methods to determine the corresponding stability chart, its structure and features, and how it differs from that of the classical Mathieu's equation.

Copyright © 2018 by ASME
Your Session has timed out. Please sign back in to continue.


Rand, R. H. , 2017, “ Lecture Notes on Nonlinear Vibrations (Version 53),” Cornell University, Ithaca, NY, accessed July 23, 2017, http://dspace.library.cornell.edu/handle/1813/28989
Magnus, W. , and Winkler, S. , 1961, “ Hill's Equation—Part II: Transformations, Approximation, Examples,” New York University, New York, Report No. BR-38. https://archive.org/details/hillsequationiit00magn
McLachlan, N. W. , 1947, Theory and Applications of Mathieu Functions, Clarendon Press, Oxford, UK.
Erdelyi, A. , 1955, Higher Transcendental Functions, Vol. III, McGraw-Hill Book Company, New York.
Stoker, J. J. , 1950, Nonlinear Vibrations in Mechanical and Electrical Systems, Interscience Publishers, New York.
Cartmell, M. , 1990, Introduction to Linear, Parametric and Nonlinear Vibrations, Chapman and Hall, London.
Ruby, L. , 1996, “ Applications of the Mathieu Equation,” Am. J. Phys., 64(1), pp. 39–44. [CrossRef]
Mathieu, E. , 1868, “ Mémoire sur Le Mouvement Vibratoire d'une Membrane de forme Elliptique,” J. Math. Pures Appl., 13, pp. 137–203.
Heine, E. , 1878, Hanbuch Der Kugelfunktionen, Vol. 2, Georg Reimer, Berlin, p. 81.
Floquet, G. , 1883, “ Sur les equations differetielles lineaires,” Ann. de L' Cole Normale Super., 12, pp. 47–88. [CrossRef]
Hill, G. W. , 1886, “ Mean Motion of the Lunar Perigee,” Acta Math., 8, pp. 1–36. [CrossRef]
Lord Rayleigh , 1887, “ Maintenance of Vibrations by Forces of Double Frequency and Propagation of Waves Through a Medium With a Periodic Structure,” Philo Magaz., 24(147), pp. 145–159. http://optoelectronics.eecs.berkeley.edu/LordRayleighA.pdf
Sieger, B. , 1908, “ Die Beugung einer ebened elektrischen Weilen an einem Schirm von elliptischen Querschnitt,” A. der P., 27(13), p. 626.
Whittaker, E. T. , 1912, “ Elliptic Cylinder Functions in Harmonic Analysis,” Paediatric Intensive Care Medicine, Vol. 1, Springer, London, p. 366.
Ince, E. , 1927, “ Research Into the Characteristic Numbers of Mathieu Equation,” Proc. R. Soc. Edinburgh, 46, pp. 20–29. [CrossRef]
Strutt, M. J. O. , 1928, “ Zur Wellenmechanik des Atomgitters,” Ann. Phys., 391(10), pp. 319–324. [CrossRef]
Stephenson, A. , 1908, “ On a New Type of Dynamic Stability,” Memoirs and Proceedings of the Manchester Literary and Philosophical Society, Vol. 52, Manchester Literary and Philosophical Society, Manchester, UK, pp. 1–10.
Stephenson, A. , 1908, “ On Induced Stability,” Philos. Mag., 15(86), pp. 233–236. [CrossRef]
Thomsen, J. J. , 2003, Vibrations and Stability: Advanced Theory, Analysis, and Tools, Springer-Verlag, Berlin.
Seyranian, A. P. , and Mailybaev, A. A. , 2003, Multiparameter Stability Theory with Mechanical Applications, Vol. 13, World Scientific, Singapore. [CrossRef]
Verhulst, F. , 2009, “ Perturbation Analysis of Parametric Resonance,” Encyclopedia of Complexity and Systems Science, R. Meyers, ed., Springer, New York.
Vlajic, N. , Liu, X. , Karki, H. , and Balachandran, B. , 2014, “ Torsional Oscillations of a Rotor With Continuous Stator Contact,” Int. J. Mech. Sci., 83, pp. 65–75. [CrossRef]
Yang, T. L. , and Rosenberg, R. M. , 1967, “ On the Vibrations of a Particle in the Plane,” Int. J. Non-Linear Mech., 2(1), pp. 1–25. [CrossRef]
Yang, T. L. , and Rosenberg, R. M. , 1968, “ On the Forced Vibrations of a Particle in the Plane,” Int. J. Non-Linear Mech., 3(1), pp. 47–63. [CrossRef]
Cole, J. D. , 1968, Perturbation Methods in Applied Mathematics, Blaisdell, Waltham, MA.
Nayfeh, A. , 1973, Perturbation Methods, Wiley, New York.
Insperger, T. , and Stépán, G. , 2011, Semi Discretization for Time Delay Systems: Stability and Engineering Applications (Applied Mathematical Science, Vol. 178), Springer Science+Business Media, New York. [CrossRef]
Kovacic, I. I. , and Rand, R. , 2014, “ Duffing-Type Oscillators With Amplitude-Independent Period,” Applied Nonlinear Dynamical Systems, Vol. 93, J. Awrejcewicz , ed., Springer, Berlin, pp. 1–10. [CrossRef]
Kovacic, I. , and Brennan, M. J. , 2011, The Duffing Equation: Nonlinear Oscillators and Their Behaviour, Wiley, Chichester, UK. [CrossRef]
Rand, R. H. , Ramani, D. V. , Keith, W. L. , and Cipolla, K. M. , 2000, “ The Quadratically Damped Mathieu Equation and Its Application to Submarine Dynamics,” Control of Noise and Vibration: New Millenium, AD-Vol. 61, ASME, New York, pp. 39–50.
Ramani, D. V. , Keith, W. L. , and Rand, R. H. , 2004, “ Perturbation Solution for Secondary Bifurcation in the Quadratically-Damped Mathieu Equation,” Int. J. Non-Linear Mech., 39(3), pp. 491–502. [CrossRef]
Morrison, T. M. , and Rand, R. H. , 2007, “ 2:1 Resonance in the Delayed Nonlinear Mathieu Equation,” Nonlinear Dyn., 50(1–2), pp. 341–352. [CrossRef]
Insperger, T. , and Stepan, G. , 2002, “Stability Chart for the Delayed Mathieu Equation,” Proc. R. Soc. A, 458(2024), pp. 1989–1998.
Butcher, E. A. , and Mann, B. P. , 2009, “ Stability Analysis and Control of Linear Periodic Delayed Systems Using Chebyshev and Temporal Finite Element Methods,” Delay Differential Equations: Recent Advances and New Directions, B. Balachandran , D. Gilsinn , and T. Kalmar-Nagy , eds., Springer, New York. [CrossRef]
Atay, F. M. , 1998, “ Van Der Pol's Oscillator Under Delayed Feedback,” J. Sound Vib., 218(2), pp. 333–339. [CrossRef]
Wirkus, S. , and Rand, R. H. , 2002, “ The Dynamics of Two Coupled Van Der Pol Oscillators With Delay Coupling,” Nonlinear Dyn., 30(3), pp. 205–221. [CrossRef]
Sah, S. M. , and Rand, R. , 2016, “ Delay Terms in the Slow Flow,” J. Appl. Nonlinear Dyn., 5(4), pp. 471–484. [CrossRef]
Bernstein, A. , and Rand, R. , 2016, “ Delay-Coupled Mathieu Equations in Synchrotron Dynamics,” J. Appl. Nonlinear Dyn., 5(3), pp. 337–348. [CrossRef]
Rand, R. H. , Sah, S. M. , and Suchorsky, M. K. , 2010, “ Fractional Mathieu Equation,” Commun. Nonlinear Sci. Numer. Simul., 15(11), pp. 3254–3262. [CrossRef]
Ross, B. , 1975, “ A Brief History and Exposition of the Fundamental Theory of Fractional Calculus,” Fractional Calculus and Its Applications (Springer Lecture Notes in Mathematics, Vol. 57), Springer, New York, pp. 1–36. [CrossRef]
Mesbahi, A. , Haeri, M. , Nazari, M. , and Butcher, E. A. , 2015, “ Fractional Delayed Damped Mathieu Equation,” Int. J. Control, 88(3), pp. 622–630. [CrossRef]
Rand, R. , Zounes, R. , and Hastings, R. , 1997, “ Dynamics of a Quasiperiodically Forced Mathieu Oscillator,” Nonlinear Dynamics: The Richard Rand 50th Anniversary Volume, A. Guran , ed., World Scientific, Singapore, pp. 203–221. [CrossRef]
Zounes, R. S. , and Rand, R. H. , 1998, “ Transition Curves in the Quasiperiodic Mathieu Equation,” SIAM J. Appl. Math., 58(4), pp. 1094–1115. [CrossRef]
Abouhazim, N. , Rand, R. H. , and Belhaq, M. , 2006, “ The Damped Nonlinear Quasiperiodic Mathieu Equation Near 2:2:1 Resonance,” Nonlinear Dyn., 45(3–4), pp. 237–247. [CrossRef]
Rand, R. , Guennoun, K. , and Belhaq, M. , 2003, “ 2:2:1 Resonance in the Quasiperiodic Mathieu Equation,” Nonlinear Dyn., 31(4), pp. 367–374. [CrossRef]
Rand, R. , and Morrison, T. , 2005, “ 2:1:1 Resonance in the Quasi-Periodic Mathieu Equation,” Nonlinear Dyn., 40(2), pp. 195–203. [CrossRef]
Sharma, A. , and Sinha, S. C. , 2017, “ An Approximate Analysis of Quasi-Periodic Systems Via Floquet Theory,” ASME J. Comput. Nonlinear Dyn., 13(2), p. 021008. [CrossRef]
Zounes, R. S. , and Rand, R. H. , 2002, “ Global Behavior of a Nonlinear Quasiperiodic Mathieu Equation,” Nonlinear Dyn., 27(1), pp. 87–105. [CrossRef]
Abramowitz, M. , and Stegun, I. , 1965, Handbook of Mathematical Functions, Dover Publications, Mineola, NY.
Kovacic, I. , Cveticanin, L. , Zukovic, M. , and Rakaric, Z. , 2016, “ Jacobi Elliptic Functions: A Review of Nonlinear Oscillatory Application Problems,” J. Sound Vib., 380, pp. 1–36. [CrossRef]
Byrd, P. , and Friedman, M. , 1954, Handbook of Elliptic Integrals for Engineers and Scientists, Springer, Berlin. [CrossRef]
Gradshteyn, I. S. , and Ryzhik, I. M. , 2000, Tables of Integrals, Series and Products, Academic Press, New York.
Kovacic, I. , and Zukovic, M. , 2014, “ A Pendulum With an Elliptic-Type Parametric Excitation: Stability Charts for a Damped and Undamped System,” Commun. Nonlinear Sci. Numer. Simul., 19(4), pp. 1185–1202. [CrossRef]
Sah, S. M. , and Mann, B. , 2012, “ Transition Curves in a Parametrically Excited Pendulum With a Force of Elliptic Type,” Proc. R. Soc. A, 468(2148), pp. 3995–4007. [CrossRef]
Dingle, R. B. , and Müller-Kirsten, H. J. W. , 1962, “ Asymptotic Expansions of Mathieu Functions and Their Characteristic Numbers,” J. Die Reine Angew. Math., 1962(211), pp. 11–32.


Grahic Jump Location
Fig. 1

(a) Mathematical pendulum whose support moves periodically in a vertical direction and (b) “the particle in the plane” problem

Grahic Jump Location
Fig. 2

(a) Two transition curves of Mathieu's equation for the first region of instability called a tongue for an undamped case (solid line), Eq. (25), (b) example of motion of point P1 located inside the tongue, and (c) example of motion of point P2 located outside the tongue

Grahic Jump Location
Fig. 3

A stability chart of classical Mathieu's equation (1) obtained by using numerical integration in conjunction with Floquet theory: gray region = unstable (U), white region = stable (S). Thick solid line-transition curves obtained by harmonic balancing, Eqs. (C1)(C9).

Grahic Jump Location
Fig. 4

(a) Analytically obtained transition curves of Mathieu's equation for the first tongue for an undamped case (solid line), Eq. (25) and the case with linear viscous damping (dashed line), Eq. (52) and (b) numerically obtained stability chart of damped of Mathieu's equation, Eq. (45)

Grahic Jump Location
Fig. 5

The δ-ϵ plane of Mathieu's equation with cubic geometric nonlinearity, Eq. (56): (a) existence of different equilibria and (b) bifurcations

Grahic Jump Location
Fig. 6

Phase portraits of the Poincaré map in the different regions of the parameter plane in a quadratically damped Mathieu equation, Eq. (66)

Grahic Jump Location
Fig. 7

(a) First instability tongue in the parameter space (δ,τ,ϵ), Eq. (99), for β=3/5, α=1, μ=0 and (b) three-dimensional stability chart in the parameter space (δ,b,ϵ) of the delayed undamped Mathieu equation with τ=2π, redrawn based on Ref. [33] (adapted from Ref. [32] with permission from The Royal Society)

Grahic Jump Location
Fig. 8

Hopf bifurcation curve, Eq. (106), for ϵ=0.05, β=3/5,α=1. LC = limit cycle, No LC = no limit cycle.

Grahic Jump Location
Fig. 9

First transition curves, Eq. (119), in the fractional Mathieu equation, Eq. (109) for c = 0.1 and different values of α

Grahic Jump Location
Fig. 10

Plot of the minimum quantity of forcing amplitude ϵmin/c, Eq. (120) necessary to produce instability as a function of fractional derivative order α

Grahic Jump Location
Fig. 11

First transition curves, Eq. (123), for the Mathieu equation with two fractional terms and delay, Eq. (122) for b = c = 0.1, different values of τ, Λ∗ satisfying Eq. (124) and: (a) α = 0.25 and (b) α = 0.75

Grahic Jump Location
Fig. 12

As a first approximation, we may think of the instabilities occurring in Eq. (125) as consisting of the union of the instabilities in Eqs. (126) and (127)

Grahic Jump Location
Fig. 13

Same information as in Fig. 12 plotted for fixed ϵ with varying ω and δ

Grahic Jump Location
Fig. 14

(a) Stability of Eq. (125) as determined directly from numerical integration. ϵ=0.1, (b) enlargement of (a) around δ=0.25,ω=1. ϵ=0.1, and (c) enlargement of (a) around δ=0.25,ω=0.1. ϵ=0.1. Black = stable and white = unstable.

Grahic Jump Location
Fig. 15

Stability of Eq. (125) as determined directly from numerical integration: (a) ϵ=0.01, (b) ϵ=0.05, (c) ϵ=0.1, (d) ϵ=0.5 (Note that Figs. 15(a)15(c) lie in the region bounded by the thick straight lines in this part), and (e) ϵ = 1. Black = stable and white = unstable.

Grahic Jump Location
Fig. 16

Period of the elliptic-type cn excitation as a function of the elliptic parameter m for cases I–III

Grahic Jump Location
Fig. 17

Stability charts of Eq. (129), where gray regions indicate instability: (a) m=−0.5, (b) m = 0.5, and (c) m = 1.5. Results obtained by means of the approximated equation of motion are labeled by thick solid lines (in all cases Ñ=10).

Grahic Jump Location
Fig. 18

Stability chart redrawn based on a figure from page 28 of Ince's paper [15]




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 eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In