0
Research Papers

An Approximate Analysis of Quasi-Periodic Systems Via Floquét Theory

[+] Author and Article Information
Ashu Sharma

Department of Mechanical Engineering,
Auburn University,
Auburn, AL 36849
e-mail: azs0111@auburn.edu

S. C. Sinha

Life Fellow ASME
Alumni Professor Emeritus,
Department of Mechanical Engineering,
Auburn University,
Auburn, AL 36849
e-mail: ssinha@eng.auburn.edu

Contributed by Design Engineering Division of ASME for publication in the JOURNAL OF COMPUTATIONAL AND NONLINEAR DYNAMICS. Manuscript received March 6, 2017; final manuscript received August 17, 2017; published online November 9, 2017. Assoc. Editor: Bogdan I. Epureanu.

J. Comput. Nonlinear Dynam 13(2), 021008 (Nov 09, 2017) (18 pages) Paper No: CND-17-1104; doi: 10.1115/1.4037797 History: Received March 06, 2017; Revised August 17, 2017

Parametrically excited linear systems with oscillatory coefficients have been generally modeled by Mathieu or Hill equations (periodic coefficients) because their stability and response can be determined by Floquét theory. However, in many cases, the parametric excitation is not periodic but consists of frequencies that are incommensurate, making them quasi-periodic. Unfortunately, there is no complete theory for linear dynamic systems with quasi-periodic coefficients. Motivated by this fact, in this work, an approximate approach has been proposed to determine the stability and response of quasi-periodic systems. It is suggested here that a quasi-periodic system may be replaced by a periodic system with an appropriate large principal period and thus making it suitable for an application of the Floquét theory. Based on this premise, a systematic approach has been developed and applied to three typical quasi-periodic systems. The approximate boundaries in stability charts obtained from the proposed method are very close to the exact boundaries of original quasi-periodic equations computed numerically using maximal Lyapunov exponents. Further, the frequency spectra of solutions generated near approximate and exact boundaries are found to be almost identical ensuring a high degree of accuracy. In addition, state transition matrices (STMs) are also computed symbolically in terms of system parameters using Chebyshev polynomials and Picard iteration method. Stability diagrams based on this approach are found to be in excellent agreement with those obtained from numerical methods. The coefficients of parametric excitation terms are not necessarily small in all cases.

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

References

Bolotin, V. V. , 1964, The Dynamic Stability of Elastic Systems, Holden-Day, San Francisco, CA.
Johnson, W. , 1980, Helicopter Theory, Princeton University Press, Princeton, NJ.
Roseau, M. , 1987, Vibration in Mechanical Systems: Analytical Methods and Applications, Springer-Verlag, New York. [CrossRef]
Streit, D. A. , Krousgrill, C. M. , and Bajaj, A. K. , 1989, “Nonlinear Response of Flexible Robotic Manipulators Performing Repetitive Tasks,” ASME J. Dyn. Syst., 111(3), pp. 470–479. [CrossRef]
Sanchez, N. E. , and Nayfeh, A. H. , 1990, “Nonlinear Rolling Motion of Ships in Longitudinal Waves,” Int. Shipbuild. Prog., 37(411), pp. 247–272. https://trid.trb.org/view.aspx?id=407297
Mingori, D. L. , 1969, “Effects of Energy Dissipation on the Attitude Stability of Dual Spin Satellite,” AIAA J., 7(1), pp. 20–27. [CrossRef]
Glass, L. , 1991, “Cardiac Arrhythmias and Circle Maps—A Classical Problem,” Chaos, 1(1), pp. 13–19. [CrossRef] [PubMed]
Powell, J. L. , and Crasemann, B. , 1961, Quantum Mechanics, Addison-Wesley, Boston, MA.
Floquét, M. G. , 1883, “Sur Les Équations Différentielles Linéaires à Coefficients Périodiques,” Ann. Sci. Ec. Norm. Superieue, 12(2), pp. 47–88. [CrossRef]
Hsu, C. S. , 1974, “On Approximating a General Linear Periodic System,” J. Math. Anal. Appl., 45(1), pp. 234–251. [CrossRef]
Sinha, S. C. , Chou, C. C. , and Denman, H. H. , 1979, “Stability Analysis of Systems With Periodic Coefficients: An Approximate Approach,” J. Sound Vib., 64(4), pp. 515–527. [CrossRef]
Hammond, C. E. , 1974, “An Application of Floquet Theory to Prediction of Mechanical Instability,” J. Am. Helicopter Soc., 19(4), pp. 14–23. [CrossRef]
Friedmann, P. , Hammond, C. E. , and Woo, T. H. , 1977, “Efficient Numerical Treatment of Periodic Systems With Application to Stability Problems,” Int. J. Numer. Methods Eng., 11(7), pp. 1117–1136. [CrossRef]
Nayfeh, A. H. , and Mook, D. T. , 1979, Nonlinear Oscillations, Wiley, New York.
Sanders, J. A. , and Verhulst, F. , 1985, Averaging Methods in Nonlinear Dynamical Systems, Springer, New York. [CrossRef]
Sinha, S. C. , and Wu, D. H. , 1991, “An Efficient Computational Scheme for the Analysis of Periodic Systems,” J. Sound Vib., 151(1), pp. 91–117. [CrossRef]
Sinha, S. C. , and Butcher, E. A. , 1997, “Symbolic Computation of Fundamental Solution Matrices for Time Periodic Dynamical Systems,” J. Sound Vib., 206(1), pp. 61–85. [CrossRef]
Kirkland, G. W. , and Sinha, S. C. , 2016, “Symbolic Computation of Quantities Associated With Time-Periodic Dynamical Systems,” ASME J. Comput. Nonlinear Dyn., 11(4), p. 041022. [CrossRef]
Sinha, S. C. , Ravindra, B. , and Henrichs, J. T. , 2000, “A General Approach in the Design of Active Controllers for Nonlinear Systems Exhibiting Chaos,” Int. J. Bifurcation Chaos, 10(1), pp. 165–178. [CrossRef]
Dávid, A. , and Sinha, S. C. , 2003, “Bifurcation Control of Nonlinear Systems With Time-Periodic Coefficients,” ASME J. Dyn. Syst., 125(4), pp. 541–548. [CrossRef]
Ko, W. , and Stockie, J. M. , 2016, “Parametric Resonance in Spherical Immersed Elastic Shells,” SIAM J. Appl. Math., 76(1), pp. 58–86. [CrossRef]
Ko, W. , and Stockie, J. M. , 2015, “An Immersed Boundary Model of the Cochlea With Parametric Forcing,” SIAM J. Appl. Math., 75(3), pp. 1065–1089. [CrossRef]
Karavitaki, K. D. , and Mountain, D. C. , 2007, “Evidence for Outer Hair Cell Driven Oscillatory Fluid Flow in the Tunnel of Corti,” Biophys. J., 92(9), pp. 3284–3293. [CrossRef] [PubMed]
Robles, L. , Ruggero, M. A. , and Rich, N. C. , 1997, “ Two-Tone Distortion on the Basilar Membrane of the Chinchilla Cochlea,” J. Neurophysiol., 77(5), pp. 2385–2399. http://jn.physiology.org/content/77/5/2385.long [PubMed]
Ruggero, M. A. , Robles, L. , Rich, N. C. , and Recio, A. , 1992, “Basilar Membrane Responses to Two-Tone and Broadband Stimuli,” Philos. Trans. R. Soc. London, Ser. B, 336(1278), pp. 307–315. [CrossRef]
Ben-Haim, S. , Darvish, N. , Mika, Y. , Rousso, B. , Felzen, B. , and Malonek, D. , 2002, “Local Cardiac Motion Control Using Applied Electrical Signals,” U.S. Patent No. 6,442,424. https://www.google.ch/patents/US6442424
Matheny, R. G. , and Taylor, C. S. , 1997, “Method of Using Vagal Nerve Stimulation in Surgery,” Cardiothoracic Systems, Inc., Cupertino, CA, U.S. Patent No. 5651378A. https://www.google.tl/patents/US5651378
Davis, S. H. , and Rosenblat, S. , 1980, “A Quasiperiodic Mathieu-Hill Equation,” SIAM J. Appl. Math., 38(1), pp. 139–155. [CrossRef]
Johnson, R. A. , and Moser, J. , 1982, “The Rotation Number for Almost Periodic Potentials,” Commun. Math. Phys., 84(3), pp. 403–438. [CrossRef]
Zounes, R. S. , and Rand, R. H. , 1998, “Transition Curves for the Quasi-Periodic Mathieu Equation,” SIAM J. Appl. Math., 58(4), pp. 1094–1115. [CrossRef]
Broer, H. , and Simó, C. , 1998, “Hill’s Equation With Quasi-Periodic Forcing: Resonance Tongues, Instability Pockets and Global Phenomena,” Bull. Braz. Math. Soc., 29(2), pp. 253–293. [CrossRef]
Puig, J. , and Simó, C. , 2011, “Resonance Tongues in Quasi-Periodic Hill-Schrödinger Equation With Three Frequencies,” Regular Chaotic Dyn., 16(1), pp. 61–78. [CrossRef]
Waters, T. J. , 2010, “Stability of a 2-Dimensional Mathieu Type System With Quasiperiodic Coefficients,” Nonlinear Dyn., 60(3), pp. 341–356. [CrossRef]
Grimshaw, R. , 1990, Nonlinear Ordinary Differential Equations, Blackwell Scientific Publications, Boston, MA.
Yakubovich, V. A. , and Starzhinski, V. M. , 1975, Linear Differential Equations With Periodic Coefficients—Part I , Wiley, New York.
Butcher, E. , and Sinha, S. C. , 1998, “Symbolic Computation of Local Stability and Bifurcation Surfaces of Time-Periodic Nonlinear Systems,” Nonlinear Dyn., 17(1), pp. 1–21. [CrossRef]
Bogoljubov, N. N. , Mitropoliskii, J. A. , and Samoilenko, A. M. , 1976, Methods of Accelerated Convergence in Nonlinear Mechanics, Hindustan Publishing, Delhi, India. [CrossRef]
Murdock, J. A. , 1978, “On the Floquet Problem for Quasiperiodic Systems,” Proc. Am. Math. Soc., 68(2), pp. 179–184. [CrossRef]
Broer, H. , and Levi, M. , 1995, “Geometrical Aspects of Stability Theory of Hill’s Equations,” Arch. Ration. Mech. Anal., 131(3), pp. 225–240. [CrossRef]
Broer, H. , and Simó, C. , 2000, “Resonance Tongues in Hill’s Equations: A Geometric Approach,” J. Differ. Equations, 166(2), pp. 290–327. [CrossRef]
Sato, C. , 1966, “Correction of Stability Curves in Hill-Meissner’s Equation,” Math. Comput., 20(93), pp. 98–106.
Brogan, W. L. , 1991, Modern Control Theory, Prentice Hall, Upper Saddle River, NJ.

Figures

Grahic Jump Location
Fig. 1

Stability diagram of the QP Hill equation with  ω1=π and ω2=7.0.  ω min=0.274334. Solid: 2Ta periodic and dashed: Ta periodic.

Grahic Jump Location
Fig. 2

Stability diagram of the QP Hill equation with  ω1=π and ω2=7.0.  ω min=0.0619600. Solid: 2Ta periodic and dashed: Ta periodic.

Grahic Jump Location
Fig. 3

Stability diagram of the QP Hill equation with  ω1=π and ω2=7.0.  ω min=0.0442270. Red/solid: 2Ta periodic and blue/dashed: Ta periodic.

Grahic Jump Location
Fig. 4

Convergence diagram of bifurcation points of the main instability regions of QP Hill equation with  ω1=π and ω2=7.0

Grahic Jump Location
Fig. 5

Poincaré maps of approximate and exact solutions constructed at a few typical points in the stability diagram of QP Hill equation with  ω1=π and ω2=7.0. (ω min=0.0442270; see Fig. 3). Left: approximate system (Eq. (19)), right: original QP system, Eq. (18). (i) Point A: aa=0.268080  b=3.5  and  d=0.0, (ii) point B: aa=0.356044  b=3.6  and  d=0.0, (iii) point C: aa=3.976490   b=2.5  and  d=0.0 , (iv) point D: aa=5.556   b=4.5  and  d=0.0, (v) point F: aa=8.730260   b=5.0  and  d=0.0, (vi) point G: aa=10.52374  b=4.0  and  d=0.0, and (vii) point H: aa=11.1  b=1.0  and  d=0.0.

Grahic Jump Location
Fig. 6

Poincaré map of an exact solution computed near the exact stability boundary (ae=0.358719  b=3.6  and  d=0.0) of QP system with ω1=π andω2=7.0

Grahic Jump Location
Fig. 7

Frequency spectrum of solutions: (i) approximate solution generated at point B (aa=0.356044  b=3.6and  d=0.0 ; see Fig. 3) using Floquét theory with Ta=142.067 and (ii) exact solution generated near the exact boundary (ae=0.358719  b=3.6  and  d=0.0)

Grahic Jump Location
Fig. 8

(i) and (ii) Poincaré maps of approximate and exact solutions, respectively, constructed at point C1 in Fig. 2 with ω1=π and ω2=7.0. ω min=0.0619600. (iii) Frequency spectrum of the approximate solution generated at point C1.

Grahic Jump Location
Fig. 9

Frequency spectrum of solutions: (i) approximate solution generated at point D (aa=5.556b=4.5and  d=0.0 ; see Fig. 3) using Floquét theory with Ta=142.067 and (ii) exact solution generated near the exact boundary (ae=5.593  b=4.5  and  d=0.0)

Grahic Jump Location
Fig. 10

Comparison of instability regions of the QP system (ω1=π and ω2=7.0, ω min=0.0442270) with periodic systems: (i) QP system versus periodic system with ω1=π and (ii) QP system versus periodic system ω2=7.0

Grahic Jump Location
Fig. 11

Comparison between the stability boundaries of the main instability regions computed using Floquét theory (Ta=142.067) and Hill’s type approach

Grahic Jump Location
Fig. 12

Stability diagram of the damped QP Hill equation with  ω1=π and ω2=7.0.  ω min=0.0442270, and d=0.1.

Grahic Jump Location
Fig. 13

Poincaré maps of the approximate and the exact solutions of the damped QP Hill equation with  ω1=π and ω2=7.0. ωmin=0.0442270 (see Fig. 12). (i) Point  : a=2.467, b=0.315, and d=0.1, (ii) point B̂ : a=2.467, b=0.313, and d=0.1, (iii) point Ĉ : a=4.120, b=3.2, and d=0.1, (iv) point D̂ : a=4.130, b=3.2, and d=0.1, (v) point Ê : a=4.140, b=3.2, and d=0.1, (vi) point F̂ : a=11.230, b=2.2, and d=0.1, (vii) point Ĝ : a=11.250, b=2.2, and d=0.1, and (viii) point Ĥ : a=11.280, b=2.2, and d=0.1.

Grahic Jump Location
Fig. 14

Comparison between the stability boundaries of the main instability regions of the damped QP Hill equation ( ω1=π and ω2=7.0) computed using three approaches: maximal Lyapunov exponent, proposed method using Floquét theory (Ta=142.067) and Hill’s type approach

Grahic Jump Location
Fig. 15

Effect of b2 on the main instability region corresponding to ω1=π, i.e., R1 in the damped QP Hill equation with  ω1=π,  ω2=7.0, and ω min=0.0442270

Grahic Jump Location
Fig. 16

Convergence diagram of bifurcation points of the main instability regions of QP Hill equation with  ω1=1 and ω2=(1+5)/2

Grahic Jump Location
Fig. 17

Stability diagram of the QP Hill equation with  ω1=1,  ω2=(1+5)/2,  ω min=0.0557281. Red: 2Ta periodic and blue: Ta periodic.

Grahic Jump Location
Fig. 18

Comparison between the stability boundaries computed using Floquét theory ( ω min=0.0557281) and those computed by Broer and Simó [31]

Grahic Jump Location
Fig. 19

Stability diagram of the damped QP Hill equation with  ω1=1 and ω2=(1+5)/2. ( ω min=0.0557281 and d=0.1).

Grahic Jump Location
Fig. 20

Convergence diagram of bifurcation points of the main instability regions of QP Hill equation with  ω1=1,  ω2=3,  and ω3=11 (see Fig. 21)

Grahic Jump Location
Fig. 21

Stability diagram of the QP Hill equation with  ω1=1,  ω2=3,  ω3=11, ω min=0.0216711. Red: 2Ta periodic and blue: Ta periodic.

Grahic Jump Location
Fig. 22

Stability diagram of the QP Hill equation with ω1=1  and  ω2=(1+5)/2 using symbolic FTM (ω min=0.0557281, mc=15, and p=20). Solid: 2Ta periodic and dashed: Ta periodic.

Grahic Jump Location
Fig. 23

Stability diagram of the QP Hill equation with  ω1=1 and ω2=(1+5)/2 using symbolic FTM (ω min=0.0557281, mc=20, and p=30). Solid: 2Ta periodic and dashed: Ta periodic.

Grahic Jump Location
Fig. 24

Stability diagram of the QP Hill equation with ω1=1.0, ω2=3, and ω3=11 using symbolic FTM (ω min=0.0216711, si=40, mc=25, and p=35). Red: 2Ta periodic and blue: Ta periodic.

Tables

Errata

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 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