Dynamics of a mass-spring-pendulum system with vastly different frequencies

Save this PDF as:

Size: px
Start display at page:

Transcription

1 Dynamics of a mass-spring-pendulum system with vastly different frequencies Hiba Sheheitli, Richard H. Rand, 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.

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

CALCULATION 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

One 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

7 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

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

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

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

Phase 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

Nonlinear 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

Three 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 *,

Lecture 5. Outline: Limit Cycles. Definition and examples How to rule out limit cycles. Poincare-Bendixson theorem Hopf bifurcations Poincare maps

Lecture 5 Outline: Limit Cycles Definition and examples How to rule out limit cycles Gradient systems Liapunov functions Dulacs criterion Poincare-Bendixson theorem Hopf bifurcations Poincare maps Limit

TWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY. P. Yu 1,2 and M. Han 1

COMMUNICATIONS ON Website: http://aimsciences.org PURE AND APPLIED ANALYSIS Volume 3, Number 3, September 2004 pp. 515 526 TWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY P. Yu 1,2

TECHNICAL NOTE: PREDICTION OF THE THRESHOLD OF GLOBAL SURF-RIDING BY AN EXTENDED MELNIKOV METHOD

10 th International Conference 441 TECHNICAL NOTE: PREDICTION OF THE THRESHOLD OF GLOBAL SURF-RIDING BY AN EXTENDED MELNIKOV METHOD Wan Wu, Virginia Polytechnic Institute and State University, wanwu@vt.edu

Dynamical behaviour of a controlled vibro-impact system

Vol 17 No 7, July 2008 c 2008 Chin. Phys. Soc. 1674-1056/2008/17(07)/2446-05 Chinese Physics B and IOP Publishing Ltd Dynamical behaviour of a controlled vibro-impact system Wang Liang( ), Xu Wei( ), and

APPPHYS217 Tuesday 25 May 2010

APPPHYS7 Tuesday 5 May Our aim today is to take a brief tour of some topics in nonlinear dynamics. Some good references include: [Perko] Lawrence Perko Differential Equations and Dynamical Systems (Springer-Verlag

SYNTHESIS OF SINEWAVE OSCILLATOR BASED ON THE MODIFIED VAN

ELECTRONICS 7 19 1 September, Sozopol, BULGARIA SYNTHESIS OF SINEWAVE OSCILLATOR BASED ON THE MODIFIED VAN DER POL EQUATION USING MELNIKOV THEORY Zhivko Dimitrov Georgiev 1, Todor Georgiev Todorov, Emil

The Stick-Slip Vibration and Bifurcation of a Vibro-Impact System with Dry Friction

Send Orders for Reprints to reprints@benthamscience.ae 308 The Open Mechanical Engineering Journal, 2014, 8, 308-313 Open Access The Stick-Slip Vibration and Bifurcation of a Vibro-Impact System with Dry

Improving convergence of incremental harmonic balance method using homotopy analysis method

Acta Mech Sin (2009) 25:707 712 DOI 10.1007/s10409-009-0256-4 RESEARCH PAPER Improving convergence of incremental harmonic balance method using homotopy analysis method Yanmao Chen Jike Liu Received: 10

DIFFERENTIAL EQUATIONS, DYNAMICAL SYSTEMS, AND AN INTRODUCTION TO CHAOS

DIFFERENTIAL EQUATIONS, DYNAMICAL SYSTEMS, AND AN INTRODUCTION TO CHAOS Morris W. Hirsch University of California, Berkeley Stephen Smale University of California, Berkeley Robert L. Devaney Boston University

An 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

OSCILLATORS WITH DELAY COUPLING

Proceedings of DETC'99 1999 ASME Design Engineering Technical Conferences September 1-15, 1999, Las Vegas, Nevada, USA D E T C 9 9 / V IB -8 3 1 8 BIFURCATIONS IN THE DYNAMICS OF TWO COUPLED VAN DER POL

THREE PROBLEMS OF RESONANCE IN COUPLED OR DRIVEN OSCILLATOR SYSTEMS

THREE PROBLEMS OF RESONANCE IN COUPLED OR DRIVEN OSCILLATOR SYSTEMS A Dissertation Presented to the Faculty of the Graduate School of Cornell University in Partial Fulfillment of the Requirements for the

Nonlinear systems, chaos and control in Engineering

Nonlinear systems, chaos and control in Engineering Module 1 block 3 One-dimensional nonlinear systems Cristina Masoller Cristina.masoller@upc.edu http://www.fisica.edu.uy/~cris/ Schedule Flows on the

STRUCTURAL CONTROL USING MODIFIED TUNED LIQUID DAMPERS

STRUCTURAL CONTROL USING MODIFIED TUNED LIQUID DAMPERS A. Samanta 1 and P. Banerji 2 1 Research Scholar, Department of Civil Engineering, Indian Institute of Technology Bombay, Mumbai, India, 2 Professor,

Multiple scale methods

Multiple scale methods G. Pedersen MEK 3100/4100, Spring 2006 March 13, 2006 1 Background Many physical problems involve more than one temporal or spatial scale. One important example is the boundary layer

Figure 5.16 Compound pendulum: (a) At rest in equilibrium, (b) General position with coordinate θ, Freebody

Lecture 27. THE COMPOUND PENDULUM Figure 5.16 Compound pendulum: (a) At rest in equilibrium, (b) General position with coordinate θ, Freebody diagram The term compound is used to distinguish the present

Nonlinear dynamics of structures with propagating cracks

Loughborough University Institutional Repository Nonlinear dynamics of structures with propagating cracks This item was submitted to Loughborough University's Institutional Repository by the/an author.

Columbia University Department of Physics QUALIFYING EXAMINATION

Columbia University Department of Physics QUALIFYING EXAMINATION Monday, January 9, 2017 11:00AM to 1:00PM Classical Physics Section 1. Classical Mechanics Two hours are permitted for the completion of

LIMITATIONS OF THE METHOD OF MULTIPLE-TIME-SCALES

LIMITATIONS OF THE METHOD OF MULTIPLE-TIME-SCALES Peter B. Kahn Department of Physics and Astronomy State University of New York Stony Brook, NY 794 and Yair Zarmi Department for Energy & Environmental

Physics 235 Chapter 4. Chapter 4 Non-Linear Oscillations and Chaos

Chapter 4 Non-Linear Oscillations and Chaos Non-Linear Differential Equations Up to now we have considered differential equations with terms that are proportional to the acceleration, the velocity, and

Phase Model for the relaxed van der Pol oscillator and its application to synchronization analysis

Phase Model for the relaxed van der Pol oscillator and its application to synchronization analysis Mimila Prost O. Collado J. Automatic Control Department, CINVESTAV IPN, A.P. 4 74 Mexico, D.F. ( e mail:

Examination paper for TMA4195 Mathematical Modeling

Department of Mathematical Sciences Examination paper for TMA4195 Mathematical Modeling Academic contact during examination: Elena Celledoni Phone: 48238584, 73593541 Examination date: 11th of December

Imaginary. Axis. Real. Axis

Name ME6 Final. I certify that I upheld the Stanford Honor code during this exam Monday December 2, 25 3:3-6:3 p.m. ffl Print your name and sign the honor code statement ffl You may use your course notes,

Dynamics of multiple pendula without gravity

June 13, 2013 Institute of Physics University of Zielona Góra Chaos 2013 Istambul, 11-15 June Dynamics of multiple pendula without gravity Wojciech Szumiński http://midori-yama.pl For Anna Finding new

Vibrations Qualifying Exam Study Material

Vibrations Qualifying Exam Study Material The candidate is expected to have a thorough understanding of engineering vibrations topics. These topics are listed below for clarification. Not all instructors

1 Lyapunov theory of stability

M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability

C. points X and Y only. D. points O, X and Y only. (Total 1 mark)

Grade 11 Physics -- Homework 16 -- Answers on a separate sheet of paper, please 1. A cart, connected to two identical springs, is oscillating with simple harmonic motion between two points X and Y that

28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod)

28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod) θ + ω 2 sin θ = 0. Indicate the stable equilibrium points as well as the unstable equilibrium points.

11 The Prandtl-Tomlinson Model for Dry Friction

The Prandtl-Tomlinson Model for Dry Friction Photodetector Laser Beam Cantilever. Introduction The development of experimental methods to investigate frictional processes on the atomic scale and numerical

Analytical Approaches for Friction Induced Vibration and Stability Analysis

Approaches for Friction Induced Vibration and Stability Analysis Anooshiravan FarshidianFar a, Mehrdad Noei Aghaee a Mechanical Engineering Department, Ferdowsi University, Azadi Square, P.O. Box: 91775-1111,

A Classical Approach to the Stark-Effect. Mridul Mehta Advisor: Prof. Enrique J. Galvez Physics Dept., Colgate University

A Classical Approach to the Stark-Effect Mridul Mehta Advisor: Prof. Enrique J. Galvez Physics Dept., Colgate University Abstract The state of an atom in the presence of an external electric field is known

Introduction to Vibrational Spectroscopy

Introduction to Vibrational Spectroscopy Harmonic oscillators The classical harmonic oscillator The uantum mechanical harmonic oscillator Harmonic approximations in molecular vibrations Vibrational spectroscopy

Unforced Mechanical Vibrations

Unforced Mechanical Vibrations Today we begin to consider applications of second order ordinary differential equations. 1. Spring-Mass Systems 2. Unforced Systems: Damped Motion 1 Spring-Mass Systems We

arxiv: v1 [nlin.ao] 19 May 2017

Feature-rich bifurcations in a simple electronic circuit Debdipta Goswami 1, and Subhankar Ray 2, 1 Department of Electrical and Computer Engineering, University of Maryland, College Park, MD 20742, USA

Chaotic motion. Phys 750 Lecture 9

Chaotic motion Phys 750 Lecture 9 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t =0to

Separatrix Map Analysis for Fractal Scatterings in Weak Interactions of Solitary Waves

Separatrix Map Analysis for Fractal Scatterings in Weak Interactions of Solitary Waves By Yi Zhu, Richard Haberman, and Jianke Yang Previous studies have shown that fractal scatterings in weak interactions

7 Planar systems of linear ODE

7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution

Lab M5: Hooke s Law and the Simple Harmonic Oscillator

M5.1 Lab M5: Hooke s Law and the Simple Harmonic Oscillator Most springs obey Hooke s Law, which states that the force exerted by the spring is proportional to the extension or compression of the spring

Hamiltonian Systems with Widely Separated. Frequencies. Mathematisch Instituut, Utrecht University

Hamiltonian Systems with Widely Separated Frequencies J.M. Tuwankotta F. Verhulst e-mail: tuwankotta@math.uu.nl, verhulst@math.uu.nl. Mathematisch Instituut, Utrecht University PO Box 80.010, 3508 TA Utrecht

A new look on the stabilization of inverted pendulum with parametric excitation and large random frequencies: analytical and numerical approaches

A new look on the stabilization of inverted pendulum with parametric excitation and large random frequencies: analytical and numerical approaches Roberto da Silva 1, Sandra D. Prado 1, Henrique A. Fernandes

Two-Body Problem. Central Potential. 1D Motion

Two-Body Problem. Central Potential. D Motion The simplest non-trivial dynamical problem is the problem of two particles. The equations of motion read. m r = F 2, () We already know that the center of

NUMERICAL STUDY OF INTRINSIC FEATURES OF ISOLAS IN A 2-DOF NONLINEAR SYSTEM

NUMERICAL STUDY OF INTRINSIC FEATURES OF ISOLAS IN A 2-DOF NONLINEAR SYSTEM T. Detroux 1, J.P. Noël 1, L. Masset 1, G. Kerschen 1, L.N. Virgin 2 1 Space Structures and Systems Lab., Department of Aerospace

Copyright 2008, University of Chicago, Department of Physics. Experiment I. RATIO OF SPECIFIC HEATS OF GASES; γ C p

Experiment I RATIO OF SPECIFIC HEATS OF GASES; γ C p / C v 1. Recommended Reading M. W. Zemansky, Heat and Thermodynamics, Fifth Edition, McGraw Hill, 1968, p. 122-132, 161-2. 2. Introduction You have

Physics 115C Homework 3

Physics 115C Homework 3 Problem 1 In this problem, it will be convenient to introduce the Einstein summation convention. Note that we can write S = i S i i where the sum is over i = x,y,z. In the Einstein

Generalized projective synchronization between two chaotic gyros with nonlinear damping

Generalized projective synchronization between two chaotic gyros with nonlinear damping Min Fu-Hong( ) Department of Electrical and Automation Engineering, Nanjing Normal University, Nanjing 210042, China

Principles of Molecular Spectroscopy

Principles of Molecular Spectroscopy What variables do we need to characterize a molecule? Nuclear and electronic configurations: What is the structure of the molecule? What are the bond lengths? How strong

1. < 0: the eigenvalues are real and have opposite signs; the fixed point is a saddle point

Solving a Linear System τ = trace(a) = a + d = λ 1 + λ 2 λ 1,2 = τ± = det(a) = ad bc = λ 1 λ 2 Classification of Fixed Points τ 2 4 1. < 0: the eigenvalues are real and have opposite signs; the fixed point

Physics 351, Spring 2015, Homework #5. Due at start of class, Friday, February 20, 2015 Course info is at positron.hep.upenn.

Physics 351, Spring 2015, Homework #5. Due at start of class, Friday, February 20, 2015 Course info is at positron.hep.upenn.edu/p351 When you finish this homework, remember to visit the feedback page

BIFURCATIONS OF PERIODIC ORBITS IN THREE-WELL DUFFING SYSTEM WITH A PHASE SHIFT

J Syst Sci Complex (11 4: 519 531 BIFURCATIONS OF PERIODIC ORBITS IN THREE-WELL DUFFING SYSTEM WITH A PHASE SHIFT Jicai HUANG Han ZHANG DOI: 1.17/s1144-1-89-3 Received: 9 May 8 / Revised: 5 December 9

Why are Discrete Maps Sufficient?

Why are Discrete Maps Sufficient? Why do dynamical systems specialists study maps of the form x n+ 1 = f ( xn), (time is discrete) when much of the world around us evolves continuously, and is thus well

Structural Dynamics Lecture 4. Outline of Lecture 4. Multi-Degree-of-Freedom Systems. Formulation of Equations of Motions. Undamped Eigenvibrations.

Outline of Multi-Degree-of-Freedom Systems Formulation of Equations of Motions. Newton s 2 nd Law Applied to Free Masses. D Alembert s Principle. Basic Equations of Motion for Forced Vibrations of Linear

Seminar 8. HAMILTON S EQUATIONS. p = L q = m q q = p m, (2) The Hamiltonian (3) creates Hamilton s equations as follows: = p ṗ = H = kq (5)

Problem 31. Derive Hamilton s equations for a one-dimensional harmonic oscillator. Seminar 8. HAMILTON S EQUATIONS Solution: The Lagrangian L = T V = 1 m q 1 kq (1) yields and hence the Hamiltonian is

Marcus Theory for Electron Transfer a short introduction

Marcus Theory for Electron Transfer a short introduction Minoia Andrea MPIP - Journal Club -Mainz - January 29, 2008 1 Contents 1 Intro 1 2 History and Concepts 2 2.1 Frank-Condon principle applied to

Dynamical analysis and circuit simulation of a new three-dimensional chaotic system

Dynamical analysis and circuit simulation of a new three-dimensional chaotic system Wang Ai-Yuan( 王爱元 ) a)b) and Ling Zhi-Hao( 凌志浩 ) a) a) Department of Automation, East China University of Science and

8.1 Bifurcations of Equilibria

1 81 Bifurcations of Equilibria Bifurcation theory studies qualitative changes in solutions as a parameter varies In general one could study the bifurcation theory of ODEs PDEs integro-differential equations

Final Exam Solution Dynamics :45 12:15. Problem 1 Bateau

Final Exam Solution Dynamics 2 191157140 31-01-2013 8:45 12:15 Problem 1 Bateau Bateau is a trapeze act by Cirque du Soleil in which artists perform aerial maneuvers on a boat shaped structure. The boat

DYNAMICS OF AN UMBILICAL CABLE FOR UNDERGROUND BORE-WELL APPLICATIONS

DYNAMICS OF AN UMBILICAL CABLE FOR UNDERGROUND BORE-WELL APPLICATIONS 1 KANNAN SANKAR & 2 JAYA KUMAR.V.K Engineering and Industrial Services, TATA Consultancy Services, Chennai, India E-mail: sankar.kannan@tcs.com,

Oscillations. PHYS 101 Previous Exam Problems CHAPTER. Simple harmonic motion Mass-spring system Energy in SHM Pendulums

PHYS 101 Previous Exam Problems CHAPTER 15 Oscillations Simple harmonic motion Mass-spring system Energy in SHM Pendulums 1. The displacement of a particle oscillating along the x axis is given as a function

Physics 2001/2051 The Compound Pendulum Experiment 4 and Helical Springs

PY001/051 Compound Pendulum and Helical Springs Experiment 4 Physics 001/051 The Compound Pendulum Experiment 4 and Helical Springs Prelab 1 Read the following background/setup and ensure you are familiar

Dynamics of a double pendulum with distributed mass. Abstract

The following article has been accepted by the American Journal of Physics. After it is published, it will be found at http://scitation.aip.org/ajp. Dynamics of a double pendulum with distributed mass

Chap. 15: Simple Harmonic Motion

Chap. 15: Simple Harmonic Motion Announcements: CAPA is due next Tuesday and next Friday. Web page: http://www.colorado.edu/physics/phys1110/phys1110_sp12/ Examples of periodic motion vibrating guitar

Report E-Project Henriette Laabsch Toni Luhdo Steffen Mitzscherling Jens Paasche Thomas Pache

Potsdam, August 006 Report E-Project Henriette Laabsch 7685 Toni Luhdo 7589 Steffen Mitzscherling 7540 Jens Paasche 7575 Thomas Pache 754 Introduction From 7 th February till 3 rd March, we had our laboratory

Global Analysis of Dynamical Systems Festschrift dedicated to Floris Takens for his 60th birthday

Global Analysis of Dynamical Systems Festschrift dedicated to Floris Takens for his 6th birthday Edited by Henk W. Broer, Bernd Krauskopf and Gert Vegter January 6, 3 Chapter Global Bifurcations of Periodic

AA242B: MECHANICAL VIBRATIONS

AA242B: MECHANICAL VIBRATIONS 1 / 50 AA242B: MECHANICAL VIBRATIONS Undamped Vibrations of n-dof Systems These slides are based on the recommended textbook: M. Géradin and D. Rixen, Mechanical Vibrations:

Penning Traps. Contents. Plasma Physics Penning Traps AJW August 16, Introduction. Clasical picture. Radiation Damping.

Penning Traps Contents Introduction Clasical picture Radiation Damping Number density B and E fields used to increase time that an electron remains within a discharge: Penning, 936. Can now trap a particle

Shallow Water Gravity Waves: A Note on the Particle Orbits

Journal of Oceanography Vol. 5, pp. 353 to 357. 1996 Shallow Water Gravity Waves: A Note on the Particle Orbits KERN E. KENYON 463 North Lane, Del Mar, CA 9014-4134, U.S.A. (Received 4 July 1995; in revised

Vibrations: Second Order Systems with One Degree of Freedom, Free Response

Single Degree of Freedom System 1.003J/1.053J Dynamics and Control I, Spring 007 Professor Thomas Peacock 5//007 Lecture 0 Vibrations: Second Order Systems with One Degree of Freedom, Free Response Single

= w. These evolve with time yielding the

1 Analytical prediction and representation of chaos. Michail Zak a Jet Propulsion Laboratory California Institute of Technology, Pasadena, CA 91109, USA Abstract. 1. Introduction The concept of randomness

Solutions for homework 5

1 Section 4.3 Solutions for homework 5 17. The following equation has repeated, real, characteristic roots. Find the general solution. y 4y + 4y = 0. The characteristic equation is λ 4λ + 4 = 0 which has

Seakeeping Models in the Frequency Domain

Seakeeping Models in the Frequency Domain (Module 6) Dr Tristan Perez Centre for Complex Dynamic Systems and Control (CDSC) Prof. Thor I Fossen Department of Engineering Cybernetics 18/09/2007 One-day

Relative Stability of Mathieu Equation in Third Zone

Relative tability of Mathieu Equation in Third Zone N. Mahmoudian, Ph.D. tudent, Nina.Mahmoudian@ndsu.nodak.edu Tel: 701.231.8303, Fax: 701.231-8913 G. Nakhaie Jazar, Assistant Professor, Reza.N.Jazar@ndsu.nodak.edu

MIPSI: Classic particle pendulum

Chapter 27 MIPSI: Classic particle pendulum (www.motiongenesis.com Textbooks Resources) The pendulum to the right consists of a small heavy metallic ball tied to a light cable which is connected to the

= y(x, t) =A cos (!t + kx)

A harmonic wave propagates horizontally along a taut string of length L = 8.0 m and mass M = 0.23 kg. The vertical displacement of the string along its length is given by y(x, t) = 0. m cos(.5 t + 0.8

Melnikov s Method Applied to a Multi-DOF Ship Model

Proceedings of the th International Ship Staility Workshop Melnikov s Method Applied to a Multi-DOF Ship Model Wan Wu, Leigh S. McCue, Department of Aerospace and Ocean Engineering, Virginia Polytechnic

Investigating Springs (Simple Harmonic Motion)

Investigating Springs (Simple Harmonic Motion) INTRODUCTION The purpose of this lab is to study the well-known force exerted by a spring The force, as given by Hooke s Law, is a function of the amount

Lecture 8. The Second Law of Thermodynamics; Energy Exchange

Lecture 8 The Second Law of Thermodynamics; Energy Exchange The second law of thermodynamics Statistics of energy exchange General definition of temperature Why heat flows from hot to cold Reading for

5.6 Unforced Mechanical Vibrations

5.6 Unforced Mechanical Vibrations 215 5.6 Unforced Mechanical Vibrations The study of vibrating mechanical systems begins here with examples for unforced systems with one degree of freedom. The main example

Oscillations in Damped Driven Pendulum: A Chaotic System

International Journal of Scientific and Innovative Mathematical Research (IJSIMR) Volume 3, Issue 10, October 2015, PP 14-27 ISSN 2347-307X (Print) & ISSN 2347-3142 (Online) www.arcjournals.org Oscillations

Lecture Notes for PHY 405 Classical Mechanics

Lecture Notes for PHY 405 Classical Mechanics From Thorton & Marion s Classical Mechanics Prepared by Dr. Joseph M. Hahn Saint Mary s University Department of Astronomy & Physics September 1, 2005 Chapter

Computer Algebra for Analyzing the Vibrations of Nonlinear Structures

Computer Algebra for Analyzing the Vibrations of Nonlinear Structures Algimantas Čepulkauskas, Regina Kulvietienė, and Genadijus Kulvietis Department of Information Technologies, Vilnius Gediminas Technical

In-Class Problems 30-32: Moment of Inertia, Torque, and Pendulum: Solutions

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Physics Physics 8.01 TEAL Fall Term 004 In-Class Problems 30-3: Moment of Inertia, Torque, and Pendulum: Solutions Problem 30 Moment of Inertia of a

Modeling and Experimentation: Mass-Spring-Damper System Dynamics

Modeling and Experimentation: Mass-Spring-Damper System Dynamics Prof. R.G. Longoria Department of Mechanical Engineering The University of Texas at Austin July 20, 2014 Overview 1 This lab is meant to

Orbital Motion in Schwarzschild Geometry

Physics 4 Lecture 29 Orbital Motion in Schwarzschild Geometry Lecture 29 Physics 4 Classical Mechanics II November 9th, 2007 We have seen, through the study of the weak field solutions of Einstein s equation

Graduate Preliminary Examination, Thursday, January 6, Part I 2

Graduate Preliminary Examination, Thursday, January 6, 2011 - Part I 2 Section A. Mechanics 1. ( Lasso) Picture from: The Lasso: a rational guide... c 1995 Carey D. Bunks A lasso is a rope of linear mass

Application of Viscous Vortex Domains Method for Solving Flow-Structure Problems

Application of Viscous Vortex Domains Method for Solving Flow-Structure Problems Yaroslav Dynnikov 1, Galina Dynnikova 1 1 Institute of Mechanics of Lomonosov Moscow State University, Michurinskiy pr.

Notes for Expansions/Series and Differential Equations

Notes for Expansions/Series and Differential Equations In the last discussion, we considered perturbation methods for constructing solutions/roots of algebraic equations. Three types of problems were illustrated

Chapter 12. Recall that when a spring is stretched a distance x, it will pull back with a force given by: F = -kx

Chapter 1 Lecture Notes Chapter 1 Oscillatory Motion Recall that when a spring is stretched a distance x, it will pull back with a force given by: F = -kx When the mass is released, the spring will pull

NONLINEAR OPTICS. Ch. 1 INTRODUCTION TO NONLINEAR OPTICS

NONLINEAR OPTICS Ch. 1 INTRODUCTION TO NONLINEAR OPTICS Nonlinear regime - Order of magnitude Origin of the nonlinearities - Induced Dipole and Polarization - Description of the classical anharmonic oscillator