Dynamics of a mass-spring-pendulum system with vastly different frequencies
|
|
- Alan Rich
- 6 years ago
- Views:
Transcription
1 Dynamics of a mass-spring-pendulum system with vastly different frequencies Hiba Sheheitli, hs497@cornell.edu Richard H. Rand, rhr2@cornell.edu Cornell University, Ithaca, NY, USA Abstract. We investigate the dynamics of a simple pendulum coupled to a horizontal mass-spring system. The spring is assumed to have a very large stiffness value such that the natural frequency of the mass-spring oscillator, when uncoupled from the pendulum, is an order of magnitude larger than that of the oscillations of the pendulum. The leading order dynamics of the autonomous coupled system is studied using the method of Direct Partition of Motion(DPM), in conjunction with a rescaling of fast time in a manner that is inspired by the WKB method. We particularly study the motions in which the amplitude of the motion of the harmonic oscillator is an order of magnitude smaller than that of the pendulum. In this regime, a pitchfork bifurcation of periodic orbits is found to occur for energy values larger that a critical value. The bifurcation gives rise to non-local periodic and quasi-periodic orbits in which the pendulum oscillates about an angle between zero and π/2 from the down right position. The bifurcating periodic orbits are nonlinear normal modes of the coupled system and correspond to fixed points of a Poincare map. An approximate expression for the value of the new fixed points of the map is obtained. These formal analytic results are confirmed by comparison with numerical integration. Keywords: Coupled oscillators, DPM, Method of direct partition of motion, WKB method, Bifurcations 1. INTRODUCTION The method of direct partition of motion (DPM), formalized by Blekhman (Blekhman, 2), serves to facilitate the study of the non-trivial effects of fast excitation that has been elaborately studied in recent years (Blekhman, 2), (Jensen, 1999), (Thomsen 23, 25). In most problems addressed in the literature, the fast excitation is due to an external source, that is, the system considered is non-autonomous. However, similar non-trivial effects could occur even if the fast excitation is internal to the system, instead of coming from an external source. An example of such a case would be a nonlinear oscillator coupled to a much faster oscillator (Tuwankotta and Verhulst, 23), (Nayfeh and Chin, 1995), (Sheheitli and Rand, 211). In these latter autonomous systems with widely fast frequencies, the leading order dynamics of the fast oscillator is unaffected by the slow oscillator. This is not the case for the system we study in this paper, as the amplitude and frequency of the fast oscillation, to leading order, are found to be a function of the amplitude of the slow oscillation. This is established by observing that the equation of the fast degree of freedom can be treated as a fast oscillator with a slowly varying frequency, for which the WKB method is particularly suited, and thus using a transformation of fast time analogous to that proposed by the WKB method (Wilcox, 1995). 2. THE MASS-SPRING-PENDULUM SYSTEM Figure 1: Schematic for the mass-spring-pendulum system We consider a simple pendulum whose point of suspension is connected to a mass on a spring that is restricted to move
2 horizontally, as shown in Fig.1. Ignoring dissipation, the system is governed by the following equations of motion: ml 2 + ml x cos + mgl sin = (M + m) x + ml cos ml 2 sin + k x = where primes denotes differentiation with respect to time τ. We introduce the following change of variables: / x = x l, t = τ g/l The nondimensionalized equations, which we refer to as the full system, become: + sin = ẍcos ẍ + Ω ) 2 x = µ ( cos 2 sin where µ = where dots represent differentiation with respect to time t. m M + m, Ω2 = kl g (M + m) (1) 2.1 Assumptions We are interested in the case where the linear oscillator has a natural frequency that is an order of magnitude larger than the linearized frequency of the pendulum, and its motion has an amplitude that is an order of magnitude smaller than that of the pendulum. This is implemented through the following rescaling: x = εχ, Ω 2 = Ω 2/ ε 2. Here, Ω and χ are O (1) quantities while ε << 1. Without loss of generality, we take Ω = 1. The rescaled equations become: + sin = ε χ cos χ + 1 ( ε χ = µ cos 2 sin ) 2 ε (2) This system has a conserved quantity that can be expressed as: h = 1 2 ε2 χ µ 2 + µε χ cos χ2 + µ (1 cos) 2.2 The bifurcation of periodic orbits We numerically integrate the full system in Eqs.(1), for typical parameter values and initial conditions (IC s), in order to illustrate the bifurcation. As an example, we take µ =.4 and Ω = 5 (ε =.2). Since the system is conservative, we will look at how the dynamics change as the energy is increased. Fig.2a shows the Poincare map for the energy level h =.5. The shown fixed point (center) of the map corresponds to a periodic orbit in which x. This periodic orbit is a nonlinear normal mode of the coupled system that appears as a nearly straight line through the origin if viewed in the configuration plane vs. x. Fig.2b shows the Poincare map for h =.7. We can see from the Poincare map that the fixed point corresponding to the nonlinear normal mode with x has lost stability and is now a saddle point of the map. Consequently, we can predict that two new fixed points (centers) were born in the process, and closed orbits of the map about the new centers would correspond to oscillations of the pendulum about a non-zero angle. The aim of this paper is to shed light on these latter non-trivial solutions, in which the pendulum oscillates about a non-zero angle, and describe their dependence on initial conditions and the parameter µ. 3. THE APPROXIMATE SOLUTION In Eqs.(2), each of the equations contains the second derivative of both χ and. We can rewrite the system of equations so that each second derivative appears in only one of the equations, as follows: 1 + (sin + µ 1 µ cos 2 2 cos sin 1ε ) χ cos = ( 1 1 χ + 1 µ cos 2 ε χ µ ( ) ) 2 sin + cos sin = 2 ε
3 Proceedings of PACAM XII poincare section x=, xdot >, Ω = 5, µ =.4, h =.5 poincare section x=, xdot >, Ω = 5, µ =.4, h = dot 1.3 dot (a) (b) Figure 2: Poincare map x =, x > for: (a) h =.5, (b) h =.7 In this latter form, the equation appears as that of a nonlinear oscillator parametrically forced by χ, which we expect to be a fast oscillation. Hence this suggests the partitioning of the motion into a slow component overlaid by a fast component, as in the DPM ansatz (Blekhman, 2). Also, we can see that the χ equation appears as that of a fast oscillator with a frequency whose magnitude is modulated by which we expect to be a slow oscillation, that is, it appears as an equation of a fast oscillator with a slowly changing frequency, similar to that which the WKB method (Wilcox, 1995) is well suited for. This suggests rescaling fast time in the following manner: ω (t) dt = dt ε or T = Zt ω (t ) dt ε and the assumed solution becomes: χ = χ (t, T ) (t, T ) = (t) + ε1 (t, T ) Here, ω (t) is to be chosen such that the fast oscillation is a perturbation of a harmonic oscillation on the new timescale T. That is, we will choose ω (t) so that the χ equation has the form: 2χ + χ + O (ε) = T 2 where ω (t) = p 1 1 µ cos2 then, an approximate expression for χ can be found using regular perturbations. After applying the standard DPM procedure (Blekhman, 2), we find that, to leading order, is governed by the following equation (see (Sheheitli and Rand, to appear) for the details): d2 1 cos sin p + sin C 2 = 2 dt 2 (1 µ cos2 ) 1 µ cos2 (3) p 1 is found to be 1 = χ cos and, to leading order, χ is given by χ C ω (t) cos T (see (Shehetli and Rand, to appear) for the details), where C is an arbitrary constant that depends on initial conditions. Then the motion of the pendulum, in the rescaled system described by Eqs.(2), can be expressed as: εχ cos Recall that Eqs.(2) are a rescaled version of the original system of interest given by Eqs.(1), where χ is related to the motion of the mass-spring oscillator as follows: x = εχ. Hence, the solution to Eqs.(1), for the assumed regime of motion, can be expressed in terms of the variables of Eqs.(1) as follows: p x cos, x εc ω (t) cos T (4)
4 4. THE SLOW DYNAMICS At the end of the DPM procedure that is described in (Sheheitli and Rand, to appear), the solution to the two degree of freedom mass-spring-pendulum system is expressed in Eqs.(4) in terms of, the leading order slow motion of the pendulum, which is governed by Eq.(3). The arbitrary constant C that appears in the equation can be expressed in terms of the initial conditions. For initial zero velocities, the initial conditions take the form: { () = () = A, { ẋ () = x () = εb () =, () A x () = εb εc ω () Now, we rewrite the equation governing as a system of two first order equations: C = B ( 1 µ cos 2 A ) 1 4 = φ, φ = sin C2 sin cos (1 µ cos 2 ) 1 µ cos 2 (5) For small enough values of C, the above system has a neutrally stable equilibrium point (center) at the origin(φ =, = ) and two saddle points at (φ =, = π, π), so that the phase portrait resembles that of the simple pendulum. As C increases in value, a pitchfork bifurcation takes place, in which the origin becomes a saddle point and two new centers are born. The critical value of C is related to the parameter µ as follows (see (Sheheitli and Rand, to appear) for the details): C 2 cr = 2 (1 µ) 3 2 In (Sheheitli and Rand, to appear), it is explained how this condition on C translates into the following condition on the energy value h: h cr = 1 µ 4.1 The predicted nonlinear normal modes Note that each value of C leads to a phase portrait filled with closed orbits, however, out of those orbits, the only one which corresponds to a solution of the full system (1) is that associated with the specific IC s that led to that value of C. An interesting case occurs when the choice of IC s results in a phase portrait that has a non-trivial equilibrium point which coincides in value with the initial amplitude, A. That is, we start with IC s of the form: { () = () = () = A, { ẋ() = x () = εb and the corresponding value of C results in non-trivial equilibrium points (centers) for the equation at: = ±E, φ = Then, if E = A, will remain equal to E for all time. It would mean that we are starting at a neutrally stable equilibrium point of the equation, so the solution will remain at that point for all time. In (Sheheitli and Rand, to appear), it is shown that these special values of initial amplitude can be expressed in terms of h and µ as: () = ±A = ±cos 1 µ h ± (h µ) 2 + 8µ (7) 4µ The corresponding value of x is expressed as: x () = εb = ε 2 (h µ (1 cosa)) (8) (6) Hence, we predict that these special initial amplitudes, with zero initial velocities, will lead to a solution in which: A x cosa Such a solution would be a nonlinear normal mode of the coupled mass-spring-pendulum system and corresponds to a non-trivial fixed points of the Poincare map. (9)
5 4.2 Relation of to the Poincare map For a given energy level, the phase portrait of the equation is filled with closed orbits and the picture is topologically similar to that of the Poincare map. That is, for a given initial condition, the resulting orbit in the phase plane corresponds to a closed orbit in the Poincare map, however, the orbits are not identical. This is due to the fact that, while, differs from by an O (1) quantity; as shown in (Sheheitli and Rand, to appear), for the points of the Poincare map, can be expressed in terms of as follows: Pm C ( 1 µ cos 2 ) 3 4 cos (1) So for given IC s, we can obtain the corresponding orbit in the Poincare map by first numerically integrating the equation to obtain and and then generating the orbit in the Poincare map by plotting the corresponding values of Pm vs.. This means that we can generate an approximate picture of the Poincare map of the full system (1) by numerically integrating the slow dynamics equation governing instead of integrating the full system (1) which contains the fast dynamics and thus requires a much smaller step size of integration. Also, by comparing this procedure with the results of numerical integration, we can obtain a check on the accuracy of the various approximations made in this work. 5. COMPARISON TO NUMERICS We compare the solution resulting from the numerical integration of the original equations with that from the integration of the equation. We have set µ =.4 and Ω = 5 (ε =.2). In the plots of Fig.3, the thick line correspond t t (a) (b) Figure 3: Comparison plots of vs. time for IC s with (a) A=π/9, B=.9756 (h=.5) (b) A=π/9, B= (h=.7) to vs. time for the solution of the numerical integration of the full system (1), and the apparent thickness is due to the fast component present in the oscillation of the pendulum; the thin line corresponds to the approximate solution, that is from the numerical integration of the equation, and captures only the leading order slow component of the pendulum oscillation. Fig.4 displays the Poincare map orbits for h = 1. Near each of the orbits, a small arrow points to the orbit which is predicted from the equation for corresponding initial conditions. We can see that that the approximate solution compares well with that from numerical integration of the full system (1). 6. CONCLUSION We have used the method of direct partition of motion to study the dynamics of a mass-spring-pendulum system, in which the harmonic oscillator is restricted to move horizontally. We have considered the case where the stiffness of the spring is very large, so that the frequency of the oscillation of the uncoupled harmonic oscillator is an order of magnitude larger than that of the uncoupled pendulum. We have also limited our attention to the regime of motion where the amplitude of motion of the harmonic oscillator is an order of magnitude smaller than that of the pendulum. Under these assumptions, an approximate expression for the solution of the two degree of freedom system is found in terms of, the leading order slow oscillation of the pendulum. An equation governing is presented and found to undergo a pitchfork bifurcation for a critical value of C which is a parameter related to the initial amplitudes of and x. It is shown.1
6 µ =.4, h = dot Figure 4: Comparison of the predicted Poincare map orbits (arrows) with those from the integration of the full system (1) that the pitchfork bifurcation in the slow dynamics equation corresponds to a pitchfork bifurcation of periodic orbits of the full system (1) that occurs as the energy is increased past a critical value which is expressed in terms of the parameter µ. This bifurcation can be seen to occur in the Poincare map of the full system (1), where the fixed point corresponding to the nonlinear normal mode x loses stability and two new centers are born in the map. The new centers correspond to new periodic motions, which are nonlinear normal modes with A x cosa, where the expression for A is found in terms of µ and h. For these modes, the motion of the pendulum is predicted to be a small fast oscillation about the non-zero value = A. Along with these special motions, quasi-periodic motions exist in which the pendulum undergoes slow oscillation about a non-zero angle, with overlaid fast oscillation. These latter orbits correspond to closed orbits about the new centers in the Poincare map. A relation between and is given for points of the Poincare map, such that the orbits of the map can be generated approximately by numerically integrating the slow dynamics equation. Finally, the approximate solution, as well as the predications made based on it, are checked against numerical integration of the full system (1) and found to agree well. 7. REFERENCES Blekhman, I.I., 2, "Vibrational mechanics-nonlinear dynamic effects, general approach, application", Singapore: World Scientific. Jensen, J.S., 1999, "Non-Trivial Effects of Fast Harmonic Excitation", PhD dissertation, DCAMM Report, S83, Dept. Solid Mechanics, Technical University of Denmark. Nayfeh, A.H., Chin, C.M., 1995, "Nonlinear Interactions in a parametrically excited system with widely spaced frequencies", Nonlinear Dynamics 7: Sheheitli, H., Rand, R.H., 211, "Dynamics of three coupled limit cycle oscillators with vastly different frequencies", Nonlinear Dynamics 64: Sheheitli, H., Rand, R.H., (to appear) "Dynamics of a mass-spring-pendulum system with vastly different frequencies". Thomsen, J.J., 23, "Vibrations and Stability, Advanced Theory, Analysis and Tools", Germany: Springer-Verlag Berlin Heidelberg. Thomsen, J.J., 25, "Slow high-frequency effects in mechanics: problems, solutions, potentials", Int. J. Bifurcation Chaos 15, pp Tuwankotta, J.M., Verhulst, F., 23, "Hamiltonian systems with widely separated frequencies", Nonlinearity 16: Wilcox, D.C., 1995, "Perturbation methods in the computer age", DCW Industries Inc., California.
ON THE DYNAMICS OF SYSTEMS WITH VASTLY DIFFERENT FREQUENCIES
ON THE DYNAMICS OF SYSTEMS WITH VASTLY DIFFERENT FREQUENCIES A Dissertation Presented to the Faculty of the Graduate School of Cornell University in Partial Fulfillment of the Requirements for the Degree
More informationLecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3.12 in Boas)
Lecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3 in Boas) As suggested in Lecture 8 the formalism of eigenvalues/eigenvectors has many applications in physics, especially in
More informationTwo models for the parametric forcing of a nonlinear oscillator
Nonlinear Dyn (007) 50:147 160 DOI 10.1007/s11071-006-9148-3 ORIGINAL ARTICLE Two models for the parametric forcing of a nonlinear oscillator Nazha Abouhazim Mohamed Belhaq Richard H. Rand Received: 3
More informationA plane autonomous system is a pair of simultaneous first-order differential equations,
Chapter 11 Phase-Plane Techniques 11.1 Plane Autonomous Systems A plane autonomous system is a pair of simultaneous first-order differential equations, ẋ = f(x, y), ẏ = g(x, y). This system has an equilibrium
More information2:1 Resonance in the delayed nonlinear Mathieu equation
Nonlinear Dyn 007) 50:31 35 DOI 10.1007/s11071-006-916-5 ORIGINAL ARTICLE :1 Resonance in the delayed nonlinear Mathieu equation Tina M. Morrison Richard H. Rand Received: 9 March 006 / Accepted: 19 September
More informationNONLINEAR NORMAL MODES OF COUPLED SELF-EXCITED OSCILLATORS
NONLINEAR NORMAL MODES OF COUPLED SELF-EXCITED OSCILLATORS Jerzy Warminski 1 1 Department of Applied Mechanics, Lublin University of Technology, Lublin, Poland, j.warminski@pollub.pl Abstract: The main
More informationMethod of Averaging for Differential Equations on an Infinite Interval
Method of Averaging for Differential Equations on an Infinite Interval Theory and Applications SUB Gottingen 7 222 045 71X ;, ' Vladimir Burd Yaroslavl State University Yaroslavl, Russia 2 ' 08A14338 Contents
More information2:2:1 Resonance in the Quasiperiodic Mathieu Equation
Nonlinear Dynamics 31: 367 374, 003. 003 Kluwer Academic Publishers. Printed in the Netherlands. ::1 Resonance in the Quasiperiodic Mathieu Equation RICHARD RAND Department of Theoretical and Applied Mechanics,
More informationNON-STATIONARY RESONANCE DYNAMICS OF THE HARMONICALLY FORCED PENDULUM
CYBERNETICS AND PHYSICS, VOL. 5, NO. 3, 016, 91 95 NON-STATIONARY RESONANCE DYNAMICS OF THE HARMONICALLY FORCED PENDULUM Leonid I. Manevitch Polymer and Composite Materials Department N. N. Semenov Institute
More informationAdditive resonances of a controlled van der Pol-Duffing oscillator
Additive resonances of a controlled van der Pol-Duffing oscillator This paper has been published in Journal of Sound and Vibration vol. 5 issue - 8 pp.-. J.C. Ji N. Zhang Faculty of Engineering University
More information2:1:1 Resonance in the Quasi-Periodic Mathieu Equation
Nonlinear Dynamics 005) 40: 95 03 c pringer 005 :: Resonance in the Quasi-Periodic Mathieu Equation RICHARD RAND and TINA MORRION Department of Theoretical & Applied Mechanics, Cornell University, Ithaca,
More informationA Model of Evolutionary Dynamics with Quasiperiodic Forcing
paper presented at Society for Experimental Mechanics (SEM) IMAC XXXIII Conference on Structural Dynamics February 2-5, 205, Orlando FL A Model of Evolutionary Dynamics with Quasiperiodic Forcing Elizabeth
More informationCALCULATION OF NONLINEAR VIBRATIONS OF PIECEWISE-LINEAR SYSTEMS USING THE SHOOTING METHOD
Vietnam Journal of Mechanics, VAST, Vol. 34, No. 3 (2012), pp. 157 167 CALCULATION OF NONLINEAR VIBRATIONS OF PIECEWISE-LINEAR SYSTEMS USING THE SHOOTING METHOD Nguyen Van Khang, Hoang Manh Cuong, Nguyen
More informationChapter 2 PARAMETRIC OSCILLATOR
CHAPTER- Chapter PARAMETRIC OSCILLATOR.1 Introduction A simple pendulum consists of a mass m suspended from a string of length L which is fixed at a pivot P. When simple pendulum is displaced to an initial
More informationDifference Resonances in a controlled van der Pol-Duffing oscillator involving time. delay
Difference Resonances in a controlled van der Pol-Duffing oscillator involving time delay This paper was published in the journal Chaos, Solitions & Fractals, vol.4, no., pp.975-98, Oct 9 J.C. Ji, N. Zhang,
More informationA New Dynamic Phenomenon in Nonlinear Circuits: State-Space Analysis of Chaotic Beats
A New Dynamic Phenomenon in Nonlinear Circuits: State-Space Analysis of Chaotic Beats DONATO CAFAGNA, GIUSEPPE GRASSI Diparnto Ingegneria Innovazione Università di Lecce via Monteroni, 73 Lecce ITALY giuseppe.grassi}@unile.it
More informationL = 1 2 a(q) q2 V (q).
Physics 3550, Fall 2011 Motion near equilibrium - Small Oscillations Relevant Sections in Text: 5.1 5.6 Motion near equilibrium 1 degree of freedom One of the most important situations in physics is motion
More informationDETC NONLINEAR DYNAMICS OF A SYSTEM OF COUPLED OSCILLATORS WITH ESSENTIAL STIFFNESS NONLINEARITIES
Proceedings of DETC 2003 2003 ASME Design Engineering Technical Conferences September 2-6, 2003, Chicago, Illinois, USA DETC2003-4446 NONLINEAR DYNAMICS OF A SYSTEM OF COUPLED OSCILLATORS WITH ESSENTIAL
More informationA New Circuit for Generating Chaos and Complexity: Analysis of the Beats Phenomenon
A New Circuit for Generating Chaos and Complexity: Analysis of the Beats Phenomenon DONATO CAFAGNA, GIUSEPPE GRASSI Diparnto Ingegneria Innovazione Università di Lecce via Monteroni, 73 Lecce ITALY Abstract:
More informationAutoparametric Resonance in Mechanical Systems. Thijs Ruijgrok Ferdinand Verhulst Radoslav Nabergoj
Autoparametric Resonance in Mechanical Systems Aleš Tondl Thijs Ruijgrok Ferdinand Verhulst Radoslav Nabergoj PUBLISHED BY THE PRESS SYNDICATE OF THE UNIVERSITY OF CAMBRIDGE The Pitt Building, Trumpington
More informationStudies on Rayleigh Equation
Johan Matheus Tuwankotta Studies on Rayleigh Equation Intisari Dalam makalah ini dipelajari persamaan oscilator tak linear yaitu persamaan Rayleigh : x x = µ ( ( x ) ) x. Masalah kestabilan lokal di sekitar
More informationDynamics of Two Coupled van der Pol Oscillators with Delay Coupling Revisited
Dynamics of Two Coupled van der Pol Oscillators with Delay Coupling Revisited arxiv:1705.03100v1 [math.ds] 8 May 017 Mark Gluzman Center for Applied Mathematics Cornell University and Richard Rand Dept.
More informationClearly the passage of an eigenvalue through to the positive real half plane leads to a qualitative change in the phase portrait, i.e.
Bifurcations We have already seen how the loss of stiffness in a linear oscillator leads to instability. In a practical situation the stiffness may not degrade in a linear fashion, and instability may
More informationDuffing-Type Oscillators with Amplitude-Independent Period
Duffing-Type Oscillators with Amplitude-Independent Period Ivana N. Kovacic and Richard H. Rand Abstract Nonlinear oscillators with hardening and softening cubic Duffing nonlinearity are considered. Such
More informationKeble College - Hilary 2014 CP3&4: Mathematical methods I&II Tutorial 5 - Waves and normal modes II
Tomi Johnson 1 Keble College - Hilary 2014 CP3&4: Mathematical methods I&II Tutorial 5 - Waves and normal modes II Prepare full solutions to the problems with a self assessment of your progress on a cover
More informationAvailable online at ScienceDirect. Procedia IUTAM 19 (2016 ) IUTAM Symposium Analytical Methods in Nonlinear Dynamics
Available online at www.sciencedirect.com ScienceDirect Procedia IUTAM 19 (2016 ) 11 18 IUTAM Symposium Analytical Methods in Nonlinear Dynamics A model of evolutionary dynamics with quasiperiodic forcing
More informationBifurcations of phase portraits of pendulum with vibrating suspension point
Bifurcations of phase portraits of pendulum with vibrating suspension point arxiv:1605.09448v [math.ds] 9 Sep 016 A.I. Neishtadt 1,,, K. Sheng 1 1 Loughborough University, Loughborough, LE11 3TU, UK Space
More informationSymmetries 2 - Rotations in Space
Symmetries 2 - Rotations in Space This symmetry is about the isotropy of space, i.e. space is the same in all orientations. Thus, if we continuously rotated an entire system in space, we expect the system
More informationSuppression of the primary resonance vibrations of a forced nonlinear system using a dynamic vibration absorber
Suppression of the primary resonance vibrations of a forced nonlinear system using a dynamic vibration absorber J.C. Ji, N. Zhang Faculty of Engineering, University of Technology, Sydney PO Box, Broadway,
More informationChapter 14 Three Ways of Treating a Linear Delay Differential Equation
Chapter 14 Three Ways of Treating a Linear Delay Differential Equation Si Mohamed Sah and Richard H. Rand Abstract This work concerns the occurrence of Hopf bifurcations in delay differential equations
More informationNON-LINEAR VIBRATION. DR. Rabinarayan Sethi,
DEPT. OF MECHANICAL ENGG., IGIT Sarang, Odisha:2012 Course Material: NON-LINEAR VIBRATION PREPARED BY DR. Rabinarayan Sethi, Assistance PROFESSOR, DEPT. OF MECHANICAL ENGG., IGIT SARANG M.Tech, B.Tech
More informationThe Effect of a Rotational Spring on the Global Stability Aspects of the Classical von Mises Model under Step Loading
Copyright cfl 2001 Tech Science Press CMES, vol.2, no.1, pp.15-26, 2001 The Effect of a Rotational Spring on the Global Stability Aspects of the Classical von Mises Model under Step Loading D. S. Sophianopoulos
More informationStability of Nonlinear Systems An Introduction
Stability of Nonlinear Systems An Introduction Michael Baldea Department of Chemical Engineering The University of Texas at Austin April 3, 2012 The Concept of Stability Consider the generic nonlinear
More information4. Complex Oscillations
4. Complex Oscillations The most common use of complex numbers in physics is for analyzing oscillations and waves. We will illustrate this with a simple but crucially important model, the damped harmonic
More informationMath 4200, Problem set 3
Math, Problem set 3 Solutions September, 13 Problem 1. ẍ = ω x. Solution. Following the general theory of conservative systems with one degree of freedom let us define the kinetic energy T and potential
More information7 Pendulum. Part II: More complicated situations
MATH 35, by T. Lakoba, University of Vermont 60 7 Pendulum. Part II: More complicated situations In this Lecture, we will pursue two main goals. First, we will take a glimpse at a method of Classical Mechanics
More informationNonlinear Oscillators: Free Response
20 Nonlinear Oscillators: Free Response Tools Used in Lab 20 Pendulums To the Instructor: This lab is just an introduction to the nonlinear phase portraits, but the connection between phase portraits and
More informationCoexistence phenomenon in autoparametric excitation of two
International Journal of Non-Linear Mechanics 40 (005) 60 70 wwwelseviercom/locate/nlm Coexistence phenomenon in autoparametric excitation of two degree of freedom systems Geoffrey Recktenwald, Richard
More informationFrequency locking in a forced Mathieu van der Pol Duffing system
Nonlinear Dyn (008) 54:3 1 DOI 10.1007/s11071-007-938-x ORIGINAL ARTICLE requency locking in a forced Mathieu van der Pol Duffing system Manoj Pandey Richard H. Rand Alan T. Zehnder Received: 4 August
More informationfor changing independent variables. Most simply for a function f(x) the Legendre transformation f(x) B(s) takes the form B(s) = xs f(x) with s = df
Physics 106a, Caltech 1 November, 2018 Lecture 10: Hamiltonian Mechanics I The Hamiltonian In the Hamiltonian formulation of dynamics each second order ODE given by the Euler- Lagrange equation in terms
More informationIn the presence of viscous damping, a more generalized form of the Lagrange s equation of motion can be written as
2 MODELING Once the control target is identified, which includes the state variable to be controlled (ex. speed, position, temperature, flow rate, etc), and once the system drives are identified (ex. force,
More informationnonlinear oscillators. method of averaging
Physics 4 Spring 7 nonlinear oscillators. method of averaging lecture notes, spring semester 7 http://www.phys.uconn.edu/ rozman/courses/p4_7s/ Last modified: April, 7 Oscillator with nonlinear friction
More information7 Two-dimensional bifurcations
7 Two-dimensional bifurcations As in one-dimensional systems: fixed points may be created, destroyed, or change stability as parameters are varied (change of topological equivalence ). In addition closed
More informationTWO VAN DER POL-DUFFING OSCILLATORS WITH HUYGENS COUPLING
ENOC-8, Saint Petersburg, Russia, June, 3 July, 8 TWO VAN DER POL-DUFFING OSCILLATORS WITH HUYGENS COUPLING V.N. Belykh, E.V. Pankratova Mathematics Department Volga State Academy Nizhny Novgorod, 3, Russia.
More informationQuasiperiodic phenomena in the Van der Pol - Mathieu equation
Quasiperiodic penomena in te Van der Pol - Matieu equation F. Veerman and F. Verulst Department of Matematics, Utrect University P.O. Box 80.010, 3508 TA Utrect Te Neterlands April 8, 009 Abstract Te Van
More informationCHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD
CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: August 22, 2018, at 08 30 12 30 Johanneberg Jan Meibohm,
More informationHigher Order Averaging : periodic solutions, linear systems and an application
Higher Order Averaging : periodic solutions, linear systems and an application Hartono and A.H.P. van der Burgh Faculty of Information Technology and Systems, Department of Applied Mathematical Analysis,
More informationPeriodic Skeletons of Nonlinear Dynamical Systems in the Problems of Global Bifurcation Analysis
Periodic Skeletons of Nonlinear Dynamical Systems in the Problems of Global Bifurcation Analysis M Zakrzhevsky, I Schukin, A Klokov and E Shilvan PERIODIC SKELETONS OF NONLINEAR DYNAMICAL SYSTEMS IN THE
More informationJournal of Applied Nonlinear Dynamics
Journal of Applied Nonlinear Dynamics 4(2) (2015) 131 140 Journal of Applied Nonlinear Dynamics https://lhscientificpublishing.com/journals/jand-default.aspx A Model of Evolutionary Dynamics with Quasiperiodic
More informationStrange dynamics of bilinear oscillator close to grazing
Strange dynamics of bilinear oscillator close to grazing Ekaterina Pavlovskaia, James Ing, Soumitro Banerjee and Marian Wiercigroch Centre for Applied Dynamics Research, School of Engineering, King s College,
More informationStabilization of Hyperbolic Chaos by the Pyragas Method
Journal of Mathematics and System Science 4 (014) 755-76 D DAVID PUBLISHING Stabilization of Hyperbolic Chaos by the Pyragas Method Sergey Belyakin, Arsen Dzanoev, Sergey Kuznetsov Physics Faculty, Moscow
More informationAPPROXIMATE ANALYTICAL SOLUTIONS TO NONLINEAR OSCILLATIONS OF NON-NATURAL SYSTEMS USING HE S ENERGY BALANCE METHOD
Progress In Electromagnetics Research M, Vol. 5, 43 54, 008 APPROXIMATE ANALYTICAL SOLUTIONS TO NONLINEAR OSCILLATIONS OF NON-NATURAL SYSTEMS USING HE S ENERGY BALANCE METHOD D. D. Ganji, S. Karimpour,
More informationDynamicsofTwoCoupledVanderPolOscillatorswithDelayCouplingRevisited
Global Journal of Science Frontier Research: F Mathematics and Decision Sciences Volume 7 Issue 5 Version.0 Year 07 Type : Double Blind Peer Reviewed International Research Journal Publisher: Global Journals
More informationDETC DELAY DIFFERENTIAL EQUATIONS IN THE DYNAMICS OF GENE COPYING
Proceedings of the ASME 2007 International Design Engineering Technical Conferences & Computers and Information in Engineering Conference IDETC/CIE 2007 September 4-7, 2007, Las Vegas, Nevada, USA DETC2007-3424
More informationNon-linear dynamics of a system of coupled oscillators with essential stiness non-linearities
Available online at www.sciencedirect.com International Journal of Non-Linear Mechanics 39 (2004) 1079 1091 Non-linear dynamics of a system of coupled oscillators with essential stiness non-linearities
More informationDynamics of four coupled phase-only oscillators
Communications in Nonlinear Science and Numerical Simulation 13 (2008) 501 507 www.elsevier.com/locate/cnsns Short communication Dynamics of four coupled phase-only oscillators Richard Rand *, Jeffrey
More informationNonlinear Dynamic Systems Homework 1
Nonlinear Dynamic Systems Homework 1 1. A particle of mass m is constrained to travel along the path shown in Figure 1, which is described by the following function yx = 5x + 1x 4, 1 where x is defined
More informationPendulum with a vibrating base
Pendulum with a vibrating base Chennai, March 14, 2006 1 1 Fast perturbations- rapidly oscillating perturbations We consider perturbations of a Hamiltonian perturbed by rapid oscillations. Later we apply
More informationA Truncated Model for Finite Amplitude Baroclinic Waves in a Channel
A Truncated Model for Finite Amplitude Baroclinic Waves in a Channel Zhiming Kuang 1 Introduction To date, studies of finite amplitude baroclinic waves have been mostly numerical. The numerical models,
More informationM2A2 Problem Sheet 3 - Hamiltonian Mechanics
MA Problem Sheet 3 - Hamiltonian Mechanics. The particle in a cone. A particle slides under gravity, inside a smooth circular cone with a vertical axis, z = k x + y. Write down its Lagrangian in a) Cartesian,
More informationOne Dimensional Dynamical Systems
16 CHAPTER 2 One Dimensional Dynamical Systems We begin by analyzing some dynamical systems with one-dimensional phase spaces, and in particular their bifurcations. All equations in this Chapter are scalar
More informationLinear and Nonlinear Oscillators (Lecture 2)
Linear and Nonlinear Oscillators (Lecture 2) January 25, 2016 7/441 Lecture outline A simple model of a linear oscillator lies in the foundation of many physical phenomena in accelerator dynamics. A typical
More informationIn-class exercises Day 1
Physics 4488/6562: Statistical Mechanics http://www.physics.cornell.edu/sethna/teaching/562/ Material for Week 11 Exercises due Mon Apr 16 Last correction at April 16, 2018, 11:19 am c 2018, James Sethna,
More informationSTABILITY. Phase portraits and local stability
MAS271 Methods for differential equations Dr. R. Jain STABILITY Phase portraits and local stability We are interested in system of ordinary differential equations of the form ẋ = f(x, y), ẏ = g(x, y),
More informationPhase Synchronization of Coupled Rossler Oscillators: Amplitude Effect
Commun. Theor. Phys. (Beijing, China) 47 (2007) pp. 265 269 c International Academic Publishers Vol. 47, No. 2, February 15, 2007 Phase Synchronization of Coupled Rossler Oscillators: Amplitude Effect
More informationTWO DIMENSIONAL FLOWS. Lecture 5: Limit Cycles and Bifurcations
TWO DIMENSIONAL FLOWS Lecture 5: Limit Cycles and Bifurcations 5. Limit cycles A limit cycle is an isolated closed trajectory [ isolated means that neighbouring trajectories are not closed] Fig. 5.1.1
More informationNonlinear Control Lecture 1: Introduction
Nonlinear Control Lecture 1: Introduction Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2011 Farzaneh Abdollahi Nonlinear Control Lecture 1 1/15 Motivation
More informationEdward Lorenz. Professor of Meteorology at the Massachusetts Institute of Technology
The Lorenz system Edward Lorenz Professor of Meteorology at the Massachusetts Institute of Technology In 1963 derived a three dimensional system in efforts to model long range predictions for the weather
More informationLecture 3 : Bifurcation Analysis
Lecture 3 : Bifurcation Analysis D. Sumpter & S.C. Nicolis October - December 2008 D. Sumpter & S.C. Nicolis General settings 4 basic bifurcations (as long as there is only one unstable mode!) steady state
More informationDaba Meshesha Gusu and O.Chandra Sekhara Reddy 1
International Journal of Basic and Applied Sciences Vol. 4. No. 1 2015. Pp.22-27 Copyright by CRDEEP. All Rights Reserved. Full Length Research Paper Solutions of Non Linear Ordinary Differential Equations
More informationThree ways of treating a linear delay differential equation
Proceedings of the 5th International Conference on Nonlinear Dynamics ND-KhPI2016 September 27-30, 2016, Kharkov, Ukraine Three ways of treating a linear delay differential equation Si Mohamed Sah 1 *,
More informationHopf Bifurcation Analysis and Approximation of Limit Cycle in Coupled Van Der Pol and Duffing Oscillators
The Open Acoustics Journal 8 9-3 9 Open Access Hopf ifurcation Analysis and Approximation of Limit Cycle in Coupled Van Der Pol and Duffing Oscillators Jianping Cai *a and Jianhe Shen b a Department of
More informationCOMPLEX DYNAMICS AND CHAOS CONTROL IN DUFFING-VAN DER POL EQUATION WITH TWO EXTERNAL PERIODIC FORCING TERMS
International J. of Math. Sci. & Engg. Appls. (IJMSEA) ISSN 0973-9424, Vol. 9 No. III (September, 2015), pp. 197-210 COMPLEX DYNAMICS AND CHAOS CONTROL IN DUFFING-VAN DER POL EQUATION WITH TWO EXTERNAL
More informationPhysics 41 HW Set 1 Chapter 15 Serway 8 th ( 7 th )
Conceptual Q: 4 (7), 7 (), 8 (6) Physics 4 HW Set Chapter 5 Serway 8 th ( 7 th ) Q4(7) Answer (c). The equilibrium position is 5 cm below the starting point. The motion is symmetric about the equilibrium
More information!T = 2# T = 2! " The velocity and acceleration of the object are found by taking the first and second derivative of the position:
A pendulum swinging back and forth or a mass oscillating on a spring are two examples of (SHM.) SHM occurs any time the position of an object as a function of time can be represented by a sine wave. We
More informationα Cubic nonlinearity coefficient. ISSN: x DOI: : /JOEMS
Journal of the Egyptian Mathematical Society Volume (6) - Issue (1) - 018 ISSN: 1110-65x DOI: : 10.1608/JOEMS.018.9468 ENHANCING PD-CONTROLLER EFFICIENCY VIA TIME- DELAYS TO SUPPRESS NONLINEAR SYSTEM OSCILLATIONS
More informationNon-Linear Response of Test Mass to External Forces and Arbitrary Motion of Suspension Point
LASER INTERFEROMETER GRAVITATIONAL WAVE OBSERVATORY -LIGO- CALIFORNIA INSTITUTE OF TECHNOLOGY MASSACHUSETTS INSTITUTE OF TECHNOLOGY Technical Note LIGO-T980005-01- D 10/28/97 Non-Linear Response of Test
More informationDYNAMICS OF THREE COUPLED VAN DER POL OSCILLATORS WITH APPLICATION TO CIRCADIAN RHYTHMS
Proceedings of IDETC/CIE 2005 ASME 2005 International Design Engineering Technical Conferences & Computers and Information in Engineering Conference September 24-28, 2005, Long Beach, California USA DETC2005-84017
More informationAutoparametric Resonance of Relaxation Oscillations
CHAPTER 4 Autoparametric Resonance of Relaation Oscillations A joint work with Ferdinand Verhulst. Has been submitted to journal. 4.. Introduction Autoparametric resonance plays an important part in nonlinear
More informationHardening nonlinearity effects on forced vibration of viscoelastic dynamical systems to external step perturbation field
Int. J. Adv. Appl. Math. and Mech. 3(2) (215) 16 32 (ISSN: 2347-2529) IJAAMM Journal homepage: www.ijaamm.com International Journal of Advances in Applied Mathematics and Mechanics Hardening nonlinearity
More information8 Example 1: The van der Pol oscillator (Strogatz Chapter 7)
8 Example 1: The van der Pol oscillator (Strogatz Chapter 7) So far we have seen some different possibilities of what can happen in two-dimensional systems (local and global attractors and bifurcations)
More informationFIRST YEAR MATHS FOR PHYSICS STUDENTS NORMAL MODES AND WAVES. Hilary Term Prof. G.G.Ross. Question Sheet 1: Normal Modes
FIRST YEAR MATHS FOR PHYSICS STUDENTS NORMAL MODES AND WAVES Hilary Term 008. Prof. G.G.Ross Question Sheet : Normal Modes [Questions marked with an asterisk (*) cover topics also covered by the unstarred
More informationIM4. Modul Mechanics. Coupled Pendulum
IM4 Modul Mechanics Coupled Pendulum Two pendulums that can exchange energy are called coupled pendulums. The gravitational force acting on the pendulums creates rotational stiffness that drives each pendulum
More informationNORMAL MODES, WAVE MOTION AND THE WAVE EQUATION. Professor G.G.Ross. Oxford University Hilary Term 2009
NORMAL MODES, WAVE MOTION AND THE WAVE EQUATION Professor G.G.Ross Oxford University Hilary Term 009 This course of twelve lectures covers material for the paper CP4: Differential Equations, Waves and
More informationSolution of a Quadratic Non-Linear Oscillator by Elliptic Homotopy Averaging Method
Math. Sci. Lett. 4, No. 3, 313-317 (215) 313 Mathematical Sciences Letters An International Journal http://dx.doi.org/1.12785/msl/4315 Solution of a Quadratic Non-Linear Oscillator by Elliptic Homotopy
More informationProblem Set Number 01, MIT (Winter-Spring 2018)
Problem Set Number 01, 18.377 MIT (Winter-Spring 2018) Rodolfo R. Rosales (MIT, Math. Dept., room 2-337, Cambridge, MA 02139) February 28, 2018 Due Thursday, March 8, 2018. Turn it in (by 3PM) at the Math.
More informationInternal and external synchronization of self-oscillating oscillators with non-identical control parameters
Internal and external synchronization of self-oscillating oscillators with non-identical control parameters Emelianova Yu.P., Kuznetsov A.P. Department of Nonlinear Processes, Saratov State University,
More information9.6 - Energy and the Simple Harmonic Oscillator *
OpenStax-CNX module: m6040 9.6 - Energy and the Simple Harmonic Oscillator * Albert Hall Based on Energy and the Simple Harmonic Oscillator by OpenStax his wor is produced by OpenStax-CNX and licensed
More informationTORSION PENDULUM: THE MECHANICAL NONLINEAR OSCILLATOR
TORSION PENDULUM: THE MECHANICAL NONLINEAR OSCILLATOR Samo Lasič, Gorazd Planinšič,, Faculty of Mathematics and Physics University of Ljubljana, Slovenija Giacomo Torzo, Department of Physics, University
More informationAn introduction to Birkhoff normal form
An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an
More informationThe SD oscillator and its attractors
The SD oscillator and its attractors Qingjie Cao, Marian Wiercigroch, Ekaterina Pavlovskaia Celso Grebogi and J Michael T Thompson Centre for Applied Dynamics Research, Department of Engineering, University
More informationCoupled Parametrically Driven Modes in Synchrotron Dynamics
paper presented at Society for Experimental Mechanics SEM) IMAC XXXIII Conference on Structural Dynamics February -5, 015, Orlando FL Coupled Parametrically Driven Modes in Synchrotron Dynamics Alexander
More informationDYNAMICS OF ASYMMETRIC NONLINEAR VIBRATION ABSORBER
8 Journal of Marine Science and Technology, Vol. 11, No. 1, pp. 8-19 (2003) DYNAMICS OF ASYMMETRIC NONLINEAR VIBRATION ABSORBER Chiou-Fong Chung and Chiang-Nan Chang* Key words: bifurcation, absorber,
More informationChapter 24 BIFURCATIONS
Chapter 24 BIFURCATIONS Abstract Keywords: Phase Portrait Fixed Point Saddle-Node Bifurcation Diagram Codimension-1 Hysteresis Hopf Bifurcation SNIC Page 1 24.1 Introduction In linear systems, responses
More informationComplex Behavior in Coupled Nonlinear Waveguides. Roy Goodman, New Jersey Institute of Technology
Complex Behavior in Coupled Nonlinear Waveguides Roy Goodman, New Jersey Institute of Technology Nonlinear Schrödinger/Gross-Pitaevskii Equation i t = r + V (r) ± Two contexts for today: Propagation of
More informationAn improved brake squeal source model in the presence of kinematic and friction nonlinearities
An improved brake squeal source model in the presence of kinematic and friction nonlinearities Osman Taha Sen, Jason T. Dreyer, and Rajendra Singh 3 Department of Mechanical Engineering, Istanbul Technical
More informationSimple Harmonic Motion Practice Problems PSI AP Physics B
Simple Harmonic Motion Practice Problems PSI AP Physics B Name Multiple Choice 1. A block with a mass M is attached to a spring with a spring constant k. The block undergoes SHM. Where is the block located
More information! 4 4! o! +! h 4 o=0! ±= ± p i And back-substituting into the linear equations gave us the ratios of the amplitudes of oscillation:.»» = A p e i! +t»»
Topic 6: Coupled Oscillators and Normal Modes Reading assignment: Hand and Finch Chapter 9 We are going to be considering the general case of a system with N degrees of freedome close to one of its stable
More informationClassical Mechanics Comprehensive Exam Solution
Classical Mechanics Comprehensive Exam Solution January 31, 011, 1:00 pm 5:pm Solve the following six problems. In the following problems, e x, e y, and e z are unit vectors in the x, y, and z directions,
More informationAutonomous Mechanical systems of 1 degree of freedom
Autonomous Mechanical systems of degree of freedom = f (, ), = y = y = f (, y) = f (, ) : position y : velocity f : force function Solution : = ( t; c, c), y = y( t; c, c ), ci = ci(, y) Integral of motion
More information