Duffing-Type Oscillators with Amplitude-Independent Period
|
|
- Arthur Atkins
- 6 years ago
- Views:
Transcription
1 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 classical conservative oscillators are known to have an amplitude-dependent period. In this work, we design oscillators with the Duffing-type restoring force but an amplitude-independent period. We present their Lagrangians, equations of motion, conservation laws, as well as solutions for motion. 1 Introduction Classical Duffing oscillators are governed by Rx C x x 3 D : (1) Their restoring force F D x x 3 includes a linear geometric term as well as a cubic geometric term: a positive sign in front of the cubic term corresponds to a hardening Duffing oscillator (HDO) and the negative one to a softening Duffing oscillator (SDO) [4, 6]. Unlike the majority of nonlinear oscillators, both of them have a closed-form exact solution, which is expressed in terms of Jacobi cn or sn elliptic functions. These solutions corresponding to the following initial conditions x./ D A; Px./ D ; () I.N. Kovacic ( ) Department of Mechanics, Faculty of Technical Sciences, University of Novi Sad, 115 Novi Sad, Serbia ivanakov@uns.ac.rs R.H. Rand Department of Mathematics, Cornell University, Ithaca, NY 14853, USA Department of Mechanical and Aerospace Engineering, Cornell University, Ithaca, NY 14853, USA rhr@cornell.edu Springer International Publishing Switzerland 14 J. Awrejcewicz (ed.), Applied Non-Linear Dynamical Systems, Springer Proceedings in Mathematics & Statistics 93, DOI 1.17/
2 I.N. Kovacic and R.H. Rand Table 1 Solutions for motion of the HDO and SDO ( in case of the SDO one requires jaj <1 for the closed orbits around the stable origin) Solution for motion Frequency Elliptic modulus HDO: x HDO D Acn.! HDO t;k HDO /! HDO D 1 C A khdo D A D! HDO SDO : x SDO D Asn.! SDO t;k SDO /! SDO A D 1 k SDO D A D! SDO A.1CA / A 1 A 1 T 8 SDO p HDO A Fig. 1 Period of the HDO (blue solid line)andsdo(red dashed line) as a function of amplitude A are given in Table 1 together with their frequencies! and elliptic moduli k (see, for example, [4, 6] for more detail). Given the fact that the period of cn and sn is T D 4K.k/ =!, where K.k/ stands for the complete elliptic integral of the first kind, the following expressions are obtained for the period of the HDO and SDO: T HDO D r A 4K p 1 C A.1CA / ; T SDO D 4K r A 1 A q 1 A! : (3) As seen from these expressions and Fig. 1, the period of the HDO and SDO depends on the amplitude A. This leads us to the question of designing an oscillator having the Duffing-type restoring force F D x x 3, but a constant, amplitudeindependent period, corresponding to the so-called isochronous oscillators [1, ]. In what follows we present such Duffing-type oscillators, both hardening and softening, modelled by Rx C G.x; Px/ C x x 3 D ; (4) finding also the corresponding Lagrangians, conservation laws, as well as solutions for motion.
3 Duffing-Type Oscillators with Amplitude-Independent Period 3 Derivation The derivation of the mathematical model (4) is based on the transformation approach [5], in which the kinetic energy E k and potential energy E p of nonlinear oscillators are made equal to the one of the simple harmonic oscillator (SHO) E ksho D 1 P X ; E psho D 1 X ; (5) which is known to have a constant, amplitude-independent period (note that its generalized coordinate is labelled here by X)..1 Case I To find the first mathematical model of the form (4), we assume that E p D E p.x/ and let E p E psho D X =, obtaining X D p E p : (6) Then, we also make the kinetic energy E k of nonlinear oscillators equal to the one of the SHO and use Eq. (6) to derive E k D X P E D p Px ; (7) 4E p where Ep D de p=dx. The differential equation of motion stemming from the Lagrangian L D E k E p has a general form Rx C! E p E p Px C E p Ep E p Ep D : (8) The last term on the left-hand side E p =Ep is required to correspond to the Duffing restoring force F =x x 3, which gives the potential energy E p D Equation (7) now yields the kinetic energy E k D x.1 x / : (9) 1.1 x / 3 Px ; (1)
4 4 I.N. Kovacic and R.H. Rand so that Eq. (8) becomes Rx 3x 1 x Px C x x 3 D : (11) Based on X = C P X = D const:, the system exhibits the first integral 1.1 x / 3 Px C x 1 x D A 1 A : (1) Taking the solution for motion of the SHO in the form X D a cos.t C /, where a and are constants, and using Eqs. (6), (9) and (), the solution for motion of Eq. (11) is obtained A p cos t D x p : (13) 1 A 1 x Numerical verifications of the analytical results for the motion (13) and phase trajectories (1) are shown in Figs. and 3 for the HDO and SDO, respectively. These figures confirm that the analytical results coincide with the solutions obtained by solving the equation of motion (11) numerically. In addition, they illustrate the fact that the period stays constant despite the fact that the amplitude A changes, i.e., that the systems perform isochronous oscillations. a b ẋ.5 x. - p 4p 6p t.5 x Fig. Isochronous oscillations of the HDO, Eq. (11) for A D :5;.5;.75:(a) time histories obtained numerically from Eq. (11) (black dots) and from Eq. (13) (blue solid line); (b) phase trajectories obtained numerically from Eq. (11) (black dots) and from Eq. (1) (blue solid line) (upper signs are used in all these equations)
5 Duffing-Type Oscillators with Amplitude-Independent Period 5 a.5 b.5 ẋ x. x p 4p 6p t - Fig. 3 Isochronous oscillations of the SDO, Eq. (11) for A D :5;.5;.75: (a) time histories obtained numerically from Eq. (11) (black dots) and from Eq. (13) (red dashed line); (b) phase trajectories obtained numerically from Eq. (11) (black dots) and from Eq. (1) (red dashed line) (lower signs are used in all these equations). Case II In this case we consider the system whose potential and kinetic energies are E p D 1 X D 1.xf/ ; E k D 1 X P D 1 Pxf C x f ; (14) where f f.i /, I D R t x.t/ dt, and f D df =d I. The corresponding Lagrange s equation is Rx C 3x Px f f C x C f f x3 D : (15) This system has two independent first integrals. The first one is the energy conservation law stemming from XP = C X = D const: Pxf C x f C.x f/ D h 1 ; h 1 D const: (16) The other first integral is related to the principle of conservation of momentum for the SHO XP C R t Xdt D X P./ D const.byusingx and XP from (14) and knowing that di=dt D x, we obtain Z Px f C x f C f.i / di D h ; h D const: (17)
6 6 I.N. Kovacic and R.H. Rand In addition, as the solution for motion for the SHO can be written down as X D a sin.t C /, the following should be satisfied: a sin.t C / D xf; acos.t C / DPxf C x f : (18) For the hardening-type nonlinearity in Eq. (15), one requires f =f D 1, which leads to Z t f HDO D exp x.t/ dt ; (19) and the equation of motion takes the form Rx C 3x Px C x C x 3 D : () Two first integrals (16) and (17) are h x3 exp.x 1 / C x i C x D h 1 ; (1) and where x 1 D Z t with initial conditions being [see Eq. ()] exp.x 1 / x 3 C x C 1 D h ; () x.t/ dt; x DPx 1 D x; x 3 DPx DPx; (3) x 1./ D ; x./ DPx 1./ D A; x 3./ DRx 1./ D : (4) Equation (1) is plotted in Fig. 4aforh 1 D 1. To analyze phase trajectories in more detail, Eq. () is squared and divided by Eq. (1) to obtain x3 C x C 1 x3 C x D B; B D const: (5) C x This expression agrees with the first integral obtained and studied in [3] and is plotted in Fig. 4b, where periodic solutions correspond to the case B > 1. Note that for the initial conditions (4), one has B D 1 C 1=A, which implies that B is always higher than unity.
7 Duffing-Type Oscillators with Amplitude-Independent Period 7 a x 3 x 1 x b 1 5 B=1.5 B=1.1 x 3 B= B=5-5 B=.1 B= x Fig. 4 (a) 3D plot of Eq.(1); (b) phase trajectories obtained from Eq. (5) for different values of B By using (18) and (19) one can derive Px sin.t C / C x sin.t C / x cos.t C / D : (6) Its solution satisfying Eq. () is sin t C arc tan 1 A x D q 1 C 1 cos : (7) t C arc tan 1 A A
8 8 I.N. Kovacic and R.H. Rand.5 x. - p 4p 6p t Fig. 5 Time response of the HDO, Eq. () for A D :5;.5;.75: numerically obtained solution from Eq. (8) (black dots) and from Eq. (7) (blue solid line) To compare the analytical solution for motion (7) with a numerically obtained solution of the equation of motion (), the latter is written down in the form «x 1 C 3 Px 1 Rx 1 CPx 1 CPx 1 3 D ; (8) and numerically integrated by using the initial conditions (4). This comparison, plotted in Fig. 5, shows that these two types of solution are in full agreement as well as that the period is amplitude-independent. The equation of motion (15) corresponds to the SDO if f =f D 1, which is satisfied for Z t f SDO D cos x.t/ dt : (9) This equation of motion is now given by Z t Rx 3x Px tan x.t/ dt C x x 3 D : (3) By using the notation given in Eq. (3), the equation of motion (3) transforms to «x 1 3 Px 1 Rx 1 tan x 1 CPx 1 Px 1 3 D ; (31) with the initial conditions given in Eq. (4). Two first integrals (16) and (17) become x3 cos x 1 x sin x 1 C x cos x 1 D h 1 ; (3)
9 Duffing-Type Oscillators with Amplitude-Independent Period 9 and x 3 cos x 1 x sin x 1 C sin x 1 D h : (33) These two integrals can be manipulated to exclude x 1 and to derive.x 1/.h 1 C. 1 h 1 C h /x C x4 / C.1 h 1 h C. 1 C h 1 h /x /.h 1.x 1/ x.h C x 1//x 3 C.. 1 C h 1 /.1 C h 1 /h C h4 /x4 3 D : (34) For the initial conditions (4) one has h 1 D A and h D. Introducing these values into Eq. (34) and solving it with respect to x 3, the following explicit solution for phase trajectories is obtained: s x 3 D.x 1/ A x 1 A : (35) Combining equations in (18) and using a D A and D =, we derive Its solution satisfying Eq. () is Px Acos t x p x A cos t C xasin t D : (36) x D A cos t p 1 A sin t : (37) This solution is plotted in Fig. 6 together with the numerical solution of Eq. (31) with Eq. (4) for different values of A. These solutions coincide and confirm isochronicity..5 x. - p 4p 6p t Fig. 6 Time response of the SDO, Eq. (3) with A D :5;.5;.75: numerically obtained solution from Eq. (31) (black dots) and from Eq. (37) (red dashed line)
10 1 I.N. Kovacic and R.H. Rand 3 Conclusions In this work we have considered nonlinear oscillators with a hardening and softening Duffing restoring force. Unlike classical conservative Duffing oscillators, which have an amplitude-dependent period, the designed Duffing-type oscillators have the period that does not change with their amplitude and are, thus, isochronous. Two separate cases are considered with respect to the form of their potential and kinetic energy, which are made equal to the corresponding energies of the SHO, which is known to be isochronous. Corresponding equations of motions are derived, as well as their solutions for motion. Numerical verifications of these isochronous solutions are provided. In addition, two independent first integrals are presented: the energyconservation law and the principle of conservation of momentum. Acknowledgements Ivana Kovacic acknowledges support received from the Provincial Secretariat for Science and Technological Development, Autonomous Province of Vojvodina (Project No ). References 1. Calogero, F.: Isochronous Systems. Oxford University Press, Oxford (8). Calogero, F.: Isochronous dynamical systems. Phil. Trans. R. Soc. A 369, (11) 3. Iacono,R., Russo, F.: Class of solvable nonlinear oscillators with isochronous orbits. Phys. Rev. E 83, 761 (4 pages) (11) 4. Kovacic, I., Brennan, M.J.: The Duffing Equation: Nonlinear Oscillators and their Behaviour. Wiley, Chichester (11) 5. Kovacic, I., Rand, R.: About a class of nonlinear oscillators with amplitude-independent frequency. Nonlinear Dyn. 74, (13) 6. Rand, R.H.: Lecture notes on nonlinear vibrations (version 53). handle/1813/8989
Dynamics of a mass-spring-pendulum system with vastly different frequencies
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
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 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 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 informationNUMERICAL METHODS. x n+1 = 2x n x 2 n. In particular: which of them gives faster convergence, and why? [Work to four decimal places.
NUMERICAL METHODS 1. Rearranging the equation x 3 =.5 gives the iterative formula x n+1 = g(x n ), where g(x) = (2x 2 ) 1. (a) Starting with x = 1, compute the x n up to n = 6, and describe what is happening.
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 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 informationLecture 2: Controllability of nonlinear systems
DISC Systems and Control Theory of Nonlinear Systems 1 Lecture 2: Controllability of nonlinear systems Nonlinear Dynamical Control Systems, Chapter 3 See www.math.rug.nl/ arjan (under teaching) for info
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 informationCanonical transformations (Lecture 4)
Canonical transformations (Lecture 4) January 26, 2016 61/441 Lecture outline We will introduce and discuss canonical transformations that conserve the Hamiltonian structure of equations of motion. Poisson
More informationSimple Harmonic Motion Concept Questions
Simple Harmonic Motion Concept Questions Question 1 Which of the following functions x(t) has a second derivative which is proportional to the negative of the function d x! " x? dt 1 1. x( t ) = at. x(
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 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 informationProblem 1: Lagrangians and Conserved Quantities. Consider the following action for a particle of mass m moving in one dimension
105A Practice Final Solutions March 13, 01 William Kelly Problem 1: Lagrangians and Conserved Quantities Consider the following action for a particle of mass m moving in one dimension S = dtl = mc dt 1
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 informationPhysics 115/242 Comparison of methods for integrating the simple harmonic oscillator.
Physics 115/4 Comparison of methods for integrating the simple harmonic oscillator. Peter Young I. THE SIMPLE HARMONIC OSCILLATOR The energy (sometimes called the Hamiltonian ) of the simple harmonic oscillator
More information18.12 FORCED-DAMPED VIBRATIONS
8. ORCED-DAMPED VIBRATIONS Vibrations A mass m is attached to a helical spring and is suspended from a fixed support as before. Damping is also provided in the system ith a dashpot (ig. 8.). Before the
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 informationUSING CARLEMAN EMBEDDING TO DISCOVER A SYSTEM S MOTION CONSTANTS
(Preprint) AAS 12-629 USING CARLEMAN EMBEDDING TO DISCOVER A SYSTEM S MOTION CONSTANTS John E. Hurtado and Andrew J. Sinclair INTRODUCTION Although the solutions with respect to time are commonly sought
More informationFree Vibration of Single-Degree-of-Freedom (SDOF) Systems
Free Vibration of Single-Degree-of-Freedom (SDOF) Systems Procedure in solving structural dynamics problems 1. Abstraction/modeling Idealize the actual structure to a simplified version, depending on the
More informationMassachusetts Institute of Technology - Physics Department
Massachusetts Institute of Technology - Physics Department Physics - 8.01 Assignment #4 October 6, 1999. It is strongly recommended that you read about a subject before it is covered in lectures. Lecture
More informationHamilton-Jacobi theory
Hamilton-Jacobi theory November 9, 04 We conclude with the crowning theorem of Hamiltonian dynamics: a proof that for any Hamiltonian dynamical system there exists a canonical transformation to a set of
More informationApplication of He s Amplitude - Frequency. Formulation for Periodic Solution. of Nonlinear Oscillators
Adv. heor. Appl. Mech., Vol.,, no. 7, - 8 Application of He s Amplitude - Frequency Formulation for Periodic Solution of Nonlinear Oscillators Jafar Langari* Islamic Azad University, Quchan Branch, Sama
More informationModeling the Duffing Equation with an Analog Computer
Modeling the Duffing Equation with an Analog Computer Matt Schmitthenner Physics Department, The College of Wooster, Wooster, Ohio 44691, USA (Dated: December 13, 2011) The goal was to model the Duffing
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 informationThe Big Picture. Python environment? Discuss Examples of unpredictability. Chaos, Scientific American (1986)
The Big Picture Python environment? Discuss Examples of unpredictability Email homework to me: chaos@cse.ucdavis.edu Chaos, Scientific American (1986) Odds, Stanislaw Lem, The New Yorker (1974) 1 Nonlinear
More informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 5. Global Bifurcations, Homoclinic chaos, Melnikov s method Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Motivation 1.1 The problem 1.2 A
More informationIntroduction to Vibration. Professor Mike Brennan
Introduction to Vibration Professor Mie Brennan Introduction to Vibration Nature of vibration of mechanical systems Free and forced vibrations Frequency response functions Fundamentals For free vibration
More informationFaculty of Computers and Information. Basic Science Department
18--018 FCI 1 Faculty of Computers and Information Basic Science Department 017-018 Prof. Nabila.M.Hassan 18--018 FCI Aims of Course: The graduates have to know the nature of vibration wave motions with
More informationOscillations. Phys101 Lectures 28, 29. Key points: Simple Harmonic Motion (SHM) SHM Related to Uniform Circular Motion The Simple Pendulum
Phys101 Lectures 8, 9 Oscillations Key points: Simple Harmonic Motion (SHM) SHM Related to Uniform Circular Motion The Simple Pendulum Ref: 11-1,,3,4. Page 1 Oscillations of a Spring If an object oscillates
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 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 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 information13 Path Planning Cubic Path P 2 P 1. θ 2
13 Path Planning Path planning includes three tasks: 1 Defining a geometric curve for the end-effector between two points. 2 Defining a rotational motion between two orientations. 3 Defining a time function
More informationFrom quantum to classical statistical mechanics. Polyatomic ideal gas.
From quantum to classical statistical mechanics. Polyatomic ideal gas. Peter Košovan peter.kosovan@natur.cuni.cz Dept. of Physical and Macromolecular Chemistry Lecture 5, Statistical Thermodynamics, MC260P105,
More informationResearch Article Analytical Approximate Solutions for the Cubic-Quintic Duffing Oscillator in Terms of Elementary Functions
Hindawi Publishing Corporation Journal of Applied Mathematics Volume 22, Article ID 28629, 6 pages doi:.55/22/28629 Research Article Analytical Approximate Solutions for the Cubic-Quintic Duffing Oscillator
More informationS13 PHY321: Final May 1, NOTE: Show all your work. No credit for unsupported answers. Turn the front page only when advised by the instructor!
Name: Student ID: S13 PHY321: Final May 1, 2013 NOTE: Show all your work. No credit for unsupported answers. Turn the front page only when advised by the instructor! The exam consists of 6 problems (60
More informationAn analytic approach to Van der Pol limit cycle
n analytic approach to Van der Pol limit cycle lberto Strumia Istituto Nazionale di lta Matematica F. Severi (i.n.d.a.m.) - Italy c CRS Registration Number: 6193976149 February 3, 2018 bstract We propose
More informationPHY321 Homework Set 10
PHY321 Homework Set 10 1. [5 pts] A small block of mass m slides without friction down a wedge-shaped block of mass M and of opening angle α. Thetriangular block itself slides along a horizontal floor,
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 informationLecture 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
More informationSeminar 6: COUPLED HARMONIC OSCILLATORS
Seminar 6: COUPLED HARMONIC OSCILLATORS 1. Lagrangian Equations of Motion Let consider a system consisting of two harmonic oscillators that are coupled together. As a model, we will use two particles attached
More informationMS 3011 Exercises. December 11, 2013
MS 3011 Exercises December 11, 2013 The exercises are divided into (A) easy (B) medium and (C) hard. If you are particularly interested I also have some projects at the end which will deepen your understanding
More informationCoordinates, phase space, constraints
Coordinates, phase space, constraints Sourendu Gupta TIFR, Mumbai, India Classical Mechanics 2011 August 1, 2011 Mechanics of a single particle For the motion of a particle (of constant mass m and position
More informationAppendix to Lecture 2
PHYS 652: Astrophysics 1 Appendix to Lecture 2 An Alternative Lagrangian In class we used an alternative Lagrangian L = g γδ ẋ γ ẋ δ, instead of the traditional L = g γδ ẋ γ ẋ δ. Here is the justification
More informationClassical mechanics of particles and fields
Classical mechanics of particles and fields D.V. Skryabin Department of Physics, University of Bath PACS numbers: The concise and transparent exposition of many topics covered in this unit can be found
More informationThe He s Amplitude-Frequency Formulation. for Solving Strongly Nonlinear Oscillators. Differential Equations
Adv. Studies heor. Phys., Vol. 5,, no. 8, 33-38 he He s Amplitude-Frequency Formulation for Solving Strongly Nonlinear Oscillators Differential Equations Jafar Langari Department of Mechanical Engineering
More informationLogistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations
Logistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations S. Y. Ha and J. Park Department of Mathematical Sciences Seoul National University Sep 23, 2013 Contents 1 Logistic Map 2 Euler and
More informationLecture 3 Dynamics 29
Lecture 3 Dynamics 29 30 LECTURE 3. DYNAMICS 3.1 Introduction Having described the states and the observables of a quantum system, we shall now introduce the rules that determine their time evolution.
More informationDIFFERENTIAL EQUATIONS AND FAMILIES OF CURVES
DIFFERENTIAL EQUATIONS AND FAMILIES OF CURVES General solution of a first-order ODE y = f(x,y) contains an arbitrary constant: y = (x,c) one curve in x,y plane for each value of c general solution can
More informationTHEODORE VORONOV DIFFERENTIABLE MANIFOLDS. Fall Last updated: November 26, (Under construction.)
4 Vector fields Last updated: November 26, 2009. (Under construction.) 4.1 Tangent vectors as derivations After we have introduced topological notions, we can come back to analysis on manifolds. Let M
More informationChapter 13 Lecture. Essential University Physics Richard Wolfson 2 nd Edition. Oscillatory Motion Pearson Education, Inc.
Chapter 13 Lecture Essential University Physics Richard Wolfson nd Edition Oscillatory Motion Slide 13-1 In this lecture you ll learn To describe the conditions under which oscillatory motion occurs To
More informationClosed form integration of a hyperelliptic, odd powers, undamped oscillator
Closed form integration of a hyperelliptic, odd powers, undamped oscillator Giovanni Mingari Scarpello and Daniele Ritelli Abstract A known one-dimensional, undamped, anharmonic oscillator whose restoring
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 informationElectric and Magnetic Forces in Lagrangian and Hamiltonian Formalism
Electric and Magnetic Forces in Lagrangian and Hamiltonian Formalism Benjamin Hornberger 1/26/1 Phy 55, Classical Electrodynamics, Prof. Goldhaber Lecture notes from Oct. 26, 21 Lecture held by Prof. Weisberger
More informationPHYS 705: Classical Mechanics. Euler s Equations
1 PHYS 705: Classical Mechanics Euler s Equations 2 Euler s Equations (set up) We have seen how to describe the kinematic properties of a rigid body. Now, we would like to get equations of motion for it.
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 informationPHYS 314 FIRST HOUR EXAM
PHYS 314 FIRST HOUR EXAM Spring 2017 This is a closed book, closed note exam. You will not need nor be allowed to use calculators or other electronic devices on this test. At this time, store all electronic
More informationPeriodic Solutions of the Duffing Harmonic Oscillator by He's Energy Balance Method
Journal of pplied and Computational Mechanics, Vol, No 1, (016), 5-1 DOI: 10055/jacm016169 Periodic Solutions of the Duffing Harmonic Oscillator by He's Energy Balance Method M El-Naggar 1, G M Ismail
More informationChapter 2 Optimal Control Problem
Chapter 2 Optimal Control Problem Optimal control of any process can be achieved either in open or closed loop. In the following two chapters we concentrate mainly on the first class. The first chapter
More informationElliptical Integrals and Functions
Appendix A Elliptical Integrals and Functions Elliptic integrals and functions are mathematical objects, which nowadays are often omitted in the mathematical curricula of universities. One quite trivial
More informationMCE693/793: Analysis and Control of Nonlinear Systems
MCE693/793: Analysis and Control of Nonlinear Systems Lyapunov Stability - I Hanz Richter Mechanical Engineering Department Cleveland State University Definition of Stability - Lyapunov Sense Lyapunov
More informationLecture 13 REVIEW. Physics 106 Spring What should we know? What should we know? Newton s Laws
Lecture 13 REVIEW Physics 106 Spring 2006 http://web.njit.edu/~sirenko/ What should we know? Vectors addition, subtraction, scalar and vector multiplication Trigonometric functions sinθ, cos θ, tan θ,
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 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 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 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 informationVibrational Motion. Chapter 5. P. J. Grandinetti. Sep. 13, Chem P. J. Grandinetti (Chem. 4300) Vibrational Motion Sep.
Vibrational Motion Chapter 5 P. J. Grandinetti Chem. 4300 Sep. 13, 2017 P. J. Grandinetti (Chem. 4300) Vibrational Motion Sep. 13, 2017 1 / 20 Simple Harmonic Oscillator Simplest model for harmonic oscillator
More informationPhys101 Lectures 28, 29. Wave Motion
Phys101 Lectures 8, 9 Wave Motion Key points: Types of Waves: Transverse and Longitudinal Mathematical Representation of a Traveling Wave The Principle of Superposition Standing Waves; Resonance Ref: 11-7,8,9,10,11,16,1,13,16.
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 informationControl of Chaos in Strongly Nonlinear Dynamic Systems
Control of Chaos in Strongly Nonlinear Dynamic Systems Lev F. Petrov Plekhanov Russian University of Economics Stremianniy per., 36, 115998, Moscow, Russia lfp@mail.ru Abstract We consider the dynamic
More informationPhysics 5153 Classical Mechanics. Solution by Quadrature-1
October 14, 003 11:47:49 1 Introduction Physics 5153 Classical Mechanics Solution by Quadrature In the previous lectures, we have reduced the number o eective degrees o reedom that are needed to solve
More informationPhysics 106a, Caltech 13 November, Lecture 13: Action, Hamilton-Jacobi Theory. Action-Angle Variables
Physics 06a, Caltech 3 November, 08 Lecture 3: Action, Hamilton-Jacobi Theory Starred sections are advanced topics for interest and future reference. The unstarred material will not be tested on the final
More informationIntroduction to Seismology Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 12.510 Introduction to Seismology Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 12.510 Introduction
More informationMath 141: Lecture 11
Math 141: Lecture 11 The Fundamental Theorem of Calculus and integration methods Bob Hough October 12, 2016 Bob Hough Math 141: Lecture 11 October 12, 2016 1 / 36 First Fundamental Theorem of Calculus
More informationPoincaré Map, Floquet Theory, and Stability of Periodic Orbits
Poincaré Map, Floquet Theory, and Stability of Periodic Orbits CDS140A Lecturer: W.S. Koon Fall, 2006 1 Poincaré Maps Definition (Poincaré Map): Consider ẋ = f(x) with periodic solution x(t). Construct
More informationLYAPUNOV EXPONENTS AND STABILITY FOR THE STOCHASTIC DUFFING-VAN DER POL OSCILLATOR
LYAPUNOV EXPONENTS AND STABILITY FOR THE STOCHASTIC DUFFING-VAN DER POL OSCILLATOR Peter H. Baxendale Department of Mathematics University of Southern California Los Angeles, CA 90089-3 USA baxendal@math.usc.edu
More informationA FIELD METHOD FOR SOLVING THE EQUATIONS OF MOTION OF EXCITED SYSTEMS UDC : (045) Ivana Kovačić
FACTA UNIVERSITATIS Series: Mechanics, Automatic Control and Robotics Vol.3, N o,, pp. 53-58 A FIELD METHOD FOR SOLVING THE EQUATIONS OF MOTION OF EXCITED SYSTEMS UDC 534.:57.98(45) Ivana Kovačić Faculty
More informationIllustrating Dynamical Symmetries in Classical Mechanics: The Laplace-Runge-Lenz Vector Revisited
Illustrating Dynamical Symmetries in Classical Mechanics: The Laplace-Runge-Lenz Vector Revisited Ross C. O Connell and Kannan Jagannathan Physics Department, Amherst College Amherst, MA 01002-5000 Abstract
More informationThe Accelerator Hamiltonian in a Straight Coordinate System
Hamiltoninan Dynamics for Particle Accelerators, Lecture 2 The Accelerator Hamiltonian in a Straight Coordinate System Andy Wolski University of Liverpool, and the Cockcroft Institute, Daresbury, UK. Given
More informationIsometric elastic deformations
Isometric elastic deformations Fares Al-Azemi and Ovidiu Calin Abstract. This paper deals with the problem of finding a class of isometric deformations of simple and closed curves, which decrease the total
More informationTheoretical physics. Deterministic chaos in classical physics. Martin Scholtz
Theoretical physics Deterministic chaos in classical physics Martin Scholtz scholtzzz@gmail.com Fundamental physical theories and role of classical mechanics. Intuitive characteristics of chaos. Newton
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 informationAdvanced Higher Mathematics of Mechanics Course Assessment Specification (C702 77)
Advanced Higher Mathematics of Mechanics Course Assessment Specification (C70 77) Valid from August 015 This edition: April 016, version.4 This specification may be reproduced in whole or in part for educational
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 informationFeedback Control and Stability of the Van der Pol Equation Subjected to External and Parametric Excitation Forces
International Journal of Applied Engineering Research ISSN 973-456 Volume 3, Number 6 (8) pp. 377-3783 Feedback Control and Stability of the Van der Pol Equation Subjected to External and Parametric Excitation
More informationTwo-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
More informationReport 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
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 informationPhase Synchronization
Phase Synchronization Lecture by: Zhibin Guo Notes by: Xiang Fan May 10, 2016 1 Introduction For any mode or fluctuation, we always have where S(x, t) is phase. If a mode amplitude satisfies ϕ k = ϕ k
More informationLecture 6. Angular Momentum and Roataion Invariance of L. x = xcosφ + ysinφ ' y = ycosφ xsinφ for φ << 1, sinφ φ, cosφ 1 x ' x + yφ y ' = y xφ,
Lecture 6 Conservation of Angular Momentum (talk about LL p. 9.3) In the previous lecture we saw that time invariance of the L leads to Energy Conservation, and translation invariance of L leads to Momentum
More informationVibrations and Waves Physics Year 1. Handout 1: Course Details
Vibrations and Waves Jan-Feb 2011 Handout 1: Course Details Office Hours Vibrations and Waves Physics Year 1 Handout 1: Course Details Dr Carl Paterson (Blackett 621, carl.paterson@imperial.ac.uk Office
More informationLocal time path integrals and their application to Lévy random walks
Local time path integrals and their application to Lévy random walks Václav Zatloukal (www.zatlovac.eu) Faculty of Nuclear Sciences and Physical Engineering Czech Technical University in Prague talk given
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 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 informationImproving 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
More informationSolution of an anti-symmetric quadratic nonlinear oscillator by a modified He s homotopy perturbation method
Nonlinear Analysis: Real World Applications, Vol. 10, Nº 1, 416-427 (2009) Solution of an anti-symmetric quadratic nonlinear oscillator by a modified He s homotopy perturbation method A. Beléndez, C. Pascual,
More information3. Riley 12.9: The equation sin x dy +2ycos x 1 dx can be reduced to a quadrature by the standard integrating factor,» Z x f(x) exp 2 dt cos t exp (2
PHYS 725 HW #4. Due 15 November 21 1. Riley 12.3: R dq dt + q C V (t); The solution is obtained with the integrating factor exp (t/rc), giving q(t) e t/rc 1 R dsv (s) e s/rc + q() With q() and V (t) V
More informationRelative 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
More informationEN Nonlinear Control and Planning in Robotics Lecture 3: Stability February 4, 2015
EN530.678 Nonlinear Control and Planning in Robotics Lecture 3: Stability February 4, 2015 Prof: Marin Kobilarov 0.1 Model prerequisites Consider ẋ = f(t, x). We will make the following basic assumptions
More information