0
Research Papers

Intrinsic Localized Modes of Harmonic Oscillations in Nonlinear Oscillator Arrays

[+] Author and Article Information
Takashi Ikeda

e-mail: tikeda@hiroshima-u.ac.jp

Keisuke Nishimura

Department of Mechanical Systems Engineering Institute of Engineering,
Hiroshima University
1-4-1, Kagamiyama,
Higashi-Hiroshima, Hiroshima 739-8527 Japan

Contributed by Design Engineering Division of ASME for publication in the Journal of Computational and Nonlinear Dynamics. Manuscript received November 27, 2012; final manuscript received February 18, 2013; published online March 26, 2013. Assoc. Editor: D. Dane Quinn.

J. Comput. Nonlinear Dynam 8(4), 041009 (Mar 26, 2013) (12 pages) Paper No: CND-12-1211; doi: 10.1115/1.4023866 History: Received November 27, 2012; Revised February 18, 2013

Intrinsic localized modes (ILMs) are investigated in an array with N Duffing oscillators that are weakly coupled with each other when each oscillator is subjected to sinusoidal excitation. The purpose of this study is to investigate the behavior of ILMs in nonlinear multi-degree-of-freedom (MDOF) systems. In the theoretical analysis, van der Pol's method is employed to determine the expressions for the frequency response curves for fundamental harmonic oscillations. In the numerical calculations, the frequency response curves are shown for N = 2 and 3 and compared with the results of the numerical simulations. Basins of attraction are shown for a two-oscillator array with hard-type nonlinearities to examine the possibility of appearance of ILMs when an oscillator is disturbed. The influences of the connecting springs for both hard- and soft-type nonlinearities on the appearance of the ILMs are examined. Increasing the values of the connecting spring constants may cause Hopf bifurcation followed by amplitude modulated motion (AMM) including chaotic vibrations. The influence of the imperfection of an oscillator is also investigated. Bifurcation sets are calculated to show the influence of the system parameters on the excitation frequency range of ILMs. Furthermore, time histories are shown for the case of N = 10, and many patterns of ILMs may appear depending on the initial conditions.

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

References

Sievers, A. J., and Takeno, S., 1988, “Intrinsic Localized Mode in Anharmonic Crystal,” Phys. Rev. Lett., 61(8), pp. 970–973. [CrossRef] [PubMed]
Flash, S., and Willis, C. R., 1998, “Discrete Breathers,” Phys. Rep., 295(5), pp. 181–264. [CrossRef]
Campbell, D. K., Flach, S., and Kivshar, Y. S., 2004, “Localizing Energy Through Nonlinearity and Discreteness,” Phys. Today, 57(1), pp. 43–49. [CrossRef]
Dauxois, T., and Peyrard, M., 1993, “Discreteness Effects on the Formation and Propagation of Breathers in Nonlinear Klein-Gordon Equations,” Phys. Rev. E, 48(6), pp. 4768–4778. [CrossRef]
Sato, M., Hubbard, B. E., Sievers, A. J., Ilic, B., Czaplewski, D. A., and Craighead, H. G., 2003,“Observation of Locked Intrinsic Localized Vibrational Modes in a Micromechanical Oscillator Array,” Phys. Rev. Lett., 90(4), p. 044102. [CrossRef] [PubMed]
Sato, M., Hubbard, B. E., English, L. Q., and Sievers, A. J., 2003, “Study of Intrinsic Localized Vibrational Modes in Micromechanical Oscillator Arrays,” Chaos, 13(2), pp. 702–715. [CrossRef] [PubMed]
Sato, M., Hubbard, B. E., and Sievers, A. J., 2006, “Colloquium: Nonlinear Energy Localization and Its Manipulation in Micromechanical Oscillator Arrays,” Rev. Modern Phys., 78(1), pp. 137–157. [CrossRef]
Dick, A. J., Balachandran, B., and Mote, C. D., Jr., 2008, “Intrinsic Localized Modes in Microresonator Arrays and Their Relationship to Nonlinear Vibration Modes,” Nonlinear Dyn., 54(1,2), pp. 13–29. [CrossRef]
Dick, A. J., Balachandran, B., and Mote, Jr., C. D, 2010, “Localization in Microresonator Arrays: Influence of Natural Frequency Tuning,” Trans. ASME J. Comput Nonlinear. Dyn.5(1), p. 11002. [CrossRef]
Chen, Q. F., Huang, L., Lai, Y.-C., and Dietz, D., 2009, “Dynamical Mechanism of Intrinsic Localized Modes in Microelectromechanical Oscillator Arrays,” Chaos, 19(1), p. 013127. [CrossRef] [PubMed]
Vakakis, A. F., and Cetinkaya, C., 1993, “Mode Localization in a Class of Multidegree-of-Freedom Nonlinear Systems With Cyclic Symmetry,” SIAM J. Appl. Math., 53(1), pp. 265–282. [CrossRef]
Vakakis, A., Nayfeh, T., and King, M., 1993, “A Multiple-Scales Analysis of Nonlinear, Localized Modes in a Cyclic Periodic System,” Trans. ASME, Appl. Mech., 60(2), pp. 388–397. [CrossRef]
King, M. E., and Vakakis, A. F., 1995, “A Very Complicated Structure of Resonances in a Nonlinear System With Cyclic Symmetry: Nonlinear Forced Localization,” Nonlinear Dyn., 7(1), pp. 85–104.
Watanabe, Y., Hamada, K., and Sugimoto, N., 2007, “Localized Oscillations of a Spatially Periodic and Articulated Structure,” Wave Motion, 45(1–2), pp. 100–117. [CrossRef]
King, M. E., Aubrecht, J., and Vakakis, A. F., 1995, “Experimental Study of Steady-State Localization in Coupled Beams With Active Nonlinearities,” J. Nonlinear Sci., 5(6), pp. 485–502. [CrossRef]
Emad, J., Vakakis, A. F., and Miller, N., 2000, “Experimental Nonlinear Localization in a Periodically Forced Repetitive System of Coupled Magnetoelastic Beams,” Physica D, 137(1–2), pp. 192–201. [CrossRef]
Kimura, M., and Hikihara, T., 2009, “Coupled Cantilever Array With Tunable On-Site Nonlinearity and Observation of Localized Oscillations,” Phys. Lett. A, 373, pp. 1257–1260. [CrossRef]
Thakur, R. B., English, L. Q., and Sievers, A. J., 2008, “Driven Intrinsic Localized Modes in a Coupled Pendulum Array,” J. Phys. D, 41(1), p. 015503. [CrossRef]
Stoker, J. J., 2005, Nonlinear Vibrations, John Wiley & Sons, New York.
Brent, R. P., 1973, Algorithms for Minimization Without Derivatives, Chapter 4, Prentice-Hall, Englewood Cliffs, NJ.
Wolf, A., Swift, J. B., Swinney, H. L., and Vastano, J. A., 1985, “Determining Lyapunov Exponents From a Time Series,” Physica, 16D, pp. 285–317. [CrossRef]
Doedel, E. J., Champneys, A. R., Fairgrieve, T. F., Kuznetsov, Y. A., Sandstede, B., and Wang, X., 1997, Continuation and Bifurcation Software for Ordinary Differential Equations (with HomCont), AUTO97, Concordia University, Canada.

Figures

Grahic Jump Location
Fig. 2

Frequency response curves of the amplitudes for N = 2 when μi = 1.0, ki = 1.0, ci = 0.02, βi = 0.1, K1 = 0.01, and F = 0.05

Grahic Jump Location
Fig. 3

Frequency response curves of the phase angles corresponding to Fig. 2

Grahic Jump Location
Fig. 4

Stationary time histories at ω = 1.10 in Fig. 2. (a) Pattern I on branches ai; (b) pattern II-1 on branches bi; (c) pattern II-2 on branches ci; and (d) pattern III on branches di.

Grahic Jump Location
Fig. 5

Basins of attraction for Fig. 2 when oscillator 1 is disturbed on branches ai at (a) ω = 1.064, (b) ω = 1.10, and (c) ω = 1.15, and on branches di at (d) ω = 1.15, (e) ω = 1.10, and (f) ω = 1.064. Basin of attraction for ILMs.

Grahic Jump Location
Fig. 6

Same as Fig. 2, but K1 = 0.02

Grahic Jump Location
Fig. 7

Same as Fig. 2, but K1 = 0.04

Grahic Jump Location
Fig. 8

Stationary time histories on branches bi in Fig. 7(a)ω = 1.090; and (b) ω = 1.0871

Grahic Jump Location
Fig. 9

Lyapunov exponents along branches bi in Fig. 7

Grahic Jump Location
Fig. 10

Bifurcation sets in the (ω, K1) plane, including Figs. 2, 6, and 7. — Saddle-node bifurcation set; – – - pitchfork bifurcation set; – - – Hopf bifurcation set.

Grahic Jump Location
Fig. 11

Frequency response curves of the amplitudes for N = 2 when μi = 1.0, ki = 1.0, ci = 0.02, βi = -0.1, K1 = 0.005, and F = 0.03

Grahic Jump Location
Fig. 12

Same as Fig. 11, but K1 = 0.01

Grahic Jump Location
Fig. 13

Same as Fig. 11, but K1 = 0.04

Grahic Jump Location
Fig. 14

Bifurcation sets in the (ω, K1) plane including Figs. 11, 12, and 13. — Saddle-node bifurcation set; – – - pitchfork bifurcation set; – - – Hopf bifurcation set.

Grahic Jump Location
Fig. 15

Same as Fig. 7, but k2 = 0.995

Grahic Jump Location
Fig. 16

Same as Fig. 13, but k2 = 0.995

Grahic Jump Location
Fig. 17

Frequency resonance curves for N = 3 when μi = 1.0, ki = 1.0, ci = 0.02, βi = 0.1, Kj = 0.01, and F = 0.05 (i = 1, 2, 3; j = 1, 2)

Grahic Jump Location
Fig. 18

Enlarged view of Fig. 17

Grahic Jump Location
Fig. 19

Stationary time histories at ω = 1.08 in Fig. 17(a) pattern II-1 on branches bi; and (b) pattern III-1 on branches ei

Grahic Jump Location
Fig. 20

Same as Fig. 17, but Kj = 0.03

Grahic Jump Location
Fig. 21

Same as Fig. 17, but Kj = 0.066

Grahic Jump Location
Fig. 22

Bifurcation sets in the (ω, Kj) plane including Figs. 17, 20, and 21. — Saddle-node bifurcation set; – – - pitchfork bifurcation set; – - – Hopf bifurcation set

Grahic Jump Location
Fig. 23

Stationary time histories for ILMs at ω = 1.08 when μi = 1.0, ki = 1.0, ci = 0.02, βi = 0.1, Kj = 0.01, and F = 0.05. (a) Six oscillators with high amplitudes; (b) three oscillators with high amplitudes; (c) one oscillator with high amplitudes

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