Transcendental cases in stability problem. Hamiltonian systems
|
|
- Abigail Shields
- 5 years ago
- Views:
Transcription
1 of Hamiltonian systems Boris S. Bardin Moscow Aviation Institute (Technical University) Faculty of Applied Mathematics and Physics Department of Theoretical Mechanics Hamiltonian Dynamics and Celestial Mechanics To honour Professor Kenneth Meyer in his 75th year 30. May 3. June, 2011
2 Contents 1 Linear problem of stability 2 3 4
3 Problem statement Let us consider Hamiltonian system dq dt = H p, dp dt = H q, (1) where q = (q 1,..., q n ) T and p = (p 1,..., p n ) T Suppose that Hamiltonian H is analytic in a small neighborhood of the origin q = p = 0, which is an equilibrium H = H 2 + H H m +... (2) Hamiltonian H is autonomous or periodically depends on t
4 Linear system Linear Hamiltonian system dx dt = Ax (3) where A = JH 2, x = (q 1,..., q n, p 1,..., p n ) T
5 Linear system Linear Hamiltonian system dx dt = Ax (3) where A = JH 2, x = (q 1,..., q n, p 1,..., p n ) T Autonomous case: A is a constant matrix Periodic case: A(t + 2π) = A(t) Lyapunov-Floquet Theorem.There exists a linear 2π periodic change of variables that transforms the 2π periodic linear system to an autonomous linear system.
6 Let us perform linear transformation x = X(t)e Bt y where X(t) is the fundamental matrix of solution of the system linear system, X(0) = E. In the new variables we have dy dt = By. (4) The eigenvalues λ i of matrix B are called a characteristic exponent. The stability problems for systems are equivalent for both x and y variables.
7 Stability in the First Approximation Theorem 1. If at least one of the characteristic exponent of the linear system has a non-zero real part, then irrespective of terms of order higher than one, the equilibrium is unstable. Theorem 2. If the Hamiltonian system autonomous and the quadratic part H 2 of the Hamiltonian is positive (or negative) definite, then the equilibrium is stable in the sense of Liapunov.
8 When nonlinear stability analysis is necessary? In what follows We suppose that real parts of all characteristic exponents are equal to zero ±iω 1, ±iω 2,..., ±iω n, In autonomous case the quadratic part H 2 of the Hamiltonian is not definite We consider autonomous Hamiltonian systems with two degrees of freedom or periodic systems with one degree of freedom
9 Nonresonant case Arnold Moser Theorem. If the Hamiltonian satisfied the following conditions: 1 There is no resonance up to fourth order, i.e. m 1 ω 1 + m 2 ω 2 0, 0 < m 1 + m 2 4 where m 1, m 2 are integer numbers 2 The coefficients of the Hamiltonian normal form H = ω 1 r 1 ω 2 r 2 + a 11 r a 12 r 1 r 2 + a 22 r O((r 1 + r 2 ) 5 2 ) q i = 2r i sin ϕ i, p i = 2r i cos ϕ i (i = 1, 2) satiesfy the following inequality a 11 ω a 12 ω 1 ω 2 + a 22 ω 2 2 0, then the Hamiltonian system is stable in the sense of Lyapunov. If a 11 ω a 12 ω 1 ω 2 + a 22 ω 2 2 = 0, then nonlinear analysis up to terms of order higher then four (with respect to q i, p i ) is necessary.
10 Third order resonance ω 1 = 2ω 2. Let us suppose that the third order resonance takes place, i.e. ω 1 = 2ω 2 then the Hamiltonian normal form reads H = 2ω 2 r 1 ω 2 r 2 + a 11 r a 12 r 1 r 2 + a 22 r ar 2 r1 cos(ϕ 1 + 2ϕ 2 ) + H (ϕ 1, r 1, ϕ 2, r 2 ), (5) Markeev s Theorem. If a 0 of Hamiltonian (5) is not equal to zero, then equilibrium is unstable in the sense of Liapunov. If a = 0 and the second condition of the Arnold - Moser Theorem is fulfilled a 11 ω a 12 ω 1 ω 2 + a 22 ω 2 2 0, then equilibrium is stable. If a = 0 and the second condition of the Arnold - Moser Theorem is not fulfilled, then additional study of stability is necessary.
11 Fourth order resonance ω 1 = 3ω 2. Suppose that the fourth order resonance takes place, i.e. ω 1 = 3ω 2 Hamiltonian takes the form H = 3ω 2 r 1 ω 2 r 2 + a 11 r a 12 r 1 r 2 + a 22 r dr 2 r1 r 2 cos(ϕ 1 + 3ϕ 2 ) + H (ϕ 1, r 1, ϕ 2, r 2 ), (6) Markeev s Theorem. If the coefficients of the Hamiltonian (6) satisfy the following inequality then equilibrium is stable; a a a 22 > 3 3d, If the coefficients of the Hamiltonian satisfy the following inequality then equilibrium is unstable. a a a 22 < 3 3d, In the case a a a 22 = 3 3d an additional study of stability is necessary.
12 How to go on? In nonresonant case the Arnold - Moser Theorem is still valid; In the case of higher order resonance (the order higher then four) Markeev s theorem 1 or Cabral Meyer theorem 2 can be applied; For strong degeneration in the case of fourth order resonance Markeev formulated stability condition. For other resonances the stability question is in the case of strong degeneration still open. 1 A. P. Markeev. Libration Points in Celestial Mechanics and Space Dynamics. Nauka, Moscow, (in Russian). 2 H. E. Cabral and K. R. Meyer. Stability of equilibria and fixed points of conservative systems. Nonlinearity, 12(5): , 1999.
13 Low order resonances ω 1 = ω 2, ω 1 = 0, or ω 2 = 0 In general position the Jordan canonical form of the matrix of the linear system is not diagonal. Let us consider resonance ω 1 = 0. The Hamiltonian takes the form H = δ 2 η ω(ξ2 2 + η 2 2) + [ M m 2 ] a m 2l,2l (ξ2 2 + η2) 2 l ξ m 2L 1 + H, (7) m=3 l=0 Sokolskii s Theorem 3. Let us suppose that a k,0 = 0 for k = 1,..., M 1 and a M,0 0. If M is odd, then the equilibrium is unstable. If M is even, then the equilibrium is unstable for δa M,0 < 0 and stable for δa M,0 > 0. 3 Sokolsky, A. G. On stability of an autonomous Hamiltonian system with two degrees of freedom under first-order resonance J. Appl. Math. Mech., 1977, 41, 20 28
14 Transcendental cases M = 3: M = 4: H = δ 2 η ω(ξ2 2 + η 2 2) + a 3,0 ξ a 1,2 ξ 1 (ξ η 2 2) + H, (8) H = δ 2 η ω(ξ2 2 + η 2 2) + a 4,0 ξ a 2,2 ξ 2 1(ξ η 2 2) + a 0,4 (ξ η 2 2) 2 + H, (9) Definition. If for any integer m the coefficients of the normal form a m,0 = 0 for any integer m, then we shall say that a transcendental case takes place.
15 Transcendental cases When transcendental cases can appear? Is equilibrium position of system stable or unstable in transcendental case?
16 When transcendental cases can appear? A. J. Maciejewski 4 has shown that transcendental case takes place if the canonical system with Hamiltonian (1) has invariant manifold M = {(ξ 1, ξ 2, η 1, η 2 ) : ξ 1 = constant, ξ 2 = η 1 = η 2 = 0} of equilibrium positions. The above sufficient condition for the existence of the transcendental case can be generalized. In fact the following lemma takes place. 4 A.J. Maciejewski. Ruch obrotowy bryly sztywnej umieszczonej w punkcie libracyjnym ograniczonego zagadnienia trzech cial. PhD thesis, Torun, 1987
17 When transcendental cases can appear? Lemma. If the equilibrium position ξ 1 = η 1 = ξ 2 = η 2 = 0 is not isolated then the transcendental case takes place. We shall prove the Lemma by contradiction. Let us suppose that the equilibrium ξ 1 = η 1 = ξ 2 = η 2 = 0 is not isolated, but the exists an integer number M such that the coefficients a k,0 (k = 1,..., M) of the normal form H = δ 2 η ω(ξ2 2 + η 2 2) + [ M m 2 ] a m 2l,2l (ξ2 2 + η2) 2 l ξ m 2L 1 m=3 l=0 satisfy the conditions a k,0 = 0 for k = 1,..., M 1 and a M,0 0 First we calculate the Jacobian ( H, H, H ) η 1 ξ 2 η 2 δ 0 0 = (η 1, ξ 2, η 2 ) = 0 ω ω = δω2 where H = H + H η1 =ξ 2 =η 2 =0
18 Hence 0 then by using implicate function theorem for for small values of ξ 1 the system of equations H = 0 η 1 H = 0 ξ 2 H = 0 η 2 can be solved with respect to η 1, ξ 2, η 2. For ξ 1 small enough the solution is analytic with respect to ξ 1 and the following asymptotic expressions take place η 1 = η 1 (ξ 1 ) = Aξ α 1 + O(ξ α+1 1 ) ξ 2 = ξ 2 (ξ 1 ) = Bξ β 1 + O(ξβ+1 1 ) η 2 = η 2 (ξ 1 ) = Cξ γ 1 + O(ξγ+1 1 ) where A, B, C are constant, α, β, γ N. (10)
19 Moreover, the truncated system of equations H = 0 η 1 H = 0 ξ 2 H = 0 η 2 where H is the normalized part of Hamiltonian (9), has the solution η 1 = ξ 2 = η 2 = 0. It follows from the above fact that α M, β M, γ M.
20 The equilibrium is not isolated, i.e. there exists a nontrivial equilibrium position in a small neighbourhood of η 1 = ξ 2 = η 2 = 0. Let us look for such an equilibrium. To this end we substitute in the equation Then we get η 1 = η 1 (ξ 1 ), ξ 2 = ξ 2 (ξ 1 ), η 2 = η 2 (ξ 1 ). H ξ 1 = 0 H ξ 1 (ξ 1, η 1 (ξ 1 ), ξ 2 (ξ 1 ), η 2 (ξ 1 )) = ξ M 1 1 (a M,0 + O(ξ 1 )) = 0 Hence a M,0 0, then for ξ 1 small enough the last equation has a unique root ξ 1 = 0. Thus, at a M,0 0 system has isolated equilibrium ξ 1 = η 1 = ξ 2 = η 2 = 0. This contradicts our assumption.
21 Suppose that the transcendental case takes place and the Hamiltonian normal form reads H = 1 2 η2 1 + ωr 2 + M [m/2] m=3 l=0 where ξ 2 = 2r 2 sin ϕ 2, η 2 = 2r 2 cos ϕ 2. a m 2l,2l ξ m 2l 1 r l 2 + H (M+1) (ξ 1, η 1, ϕ, r) Isoenergetic reduction. Consider the motion on fixed energy level H = 0. For small values of ξ 1, η 1 r 2 = K (ξ 1, η 1, ϕ 2 ) = 1 2ω η2 1 + K (3) (ξ 1, η 1, ϕ 2 ) (11) K (3) is periodic function of variable ϕ. Whittaker equations dξ 1 = K, dϕ 2 η 1 dη 1 = K (12) dϕ 2 ξ 1
22 By a proper choice of variables ξ 1, η 1, ξ 2, η 2 we can represent the invariant manifold of equilibrium positions in the form M = {(ξ 1, ξ 2, η 1, η 2 ) : ξ 1 = constant, ξ 2 = η 1 = η 2 = 0} This manifold lies on the energy level H = 0. That is why equations (12) have solution ξ 1 = const, η 1 = 0. It means that ( K (ξ 1, η 1, ϕ 2 ) = η ω a ) 1,2 2ω ξ 1η = η1 2 G(ξ 1, η 1, ϕ 2 ), (13) where G = 1 2ω a 1,2 2ω ξ 1η
23 Let us consider Chetaev function V = ξ 1 η 1 V is positive in the domain ξ 1 > 0, η 1 > 0 and is equal to zero on its border. Let us calculate the derivative of V with respect to equations (12). dv = dξ ( ) ( ) 1 dη 1 η 1 + ξ 1 = η 1 2η 1 G + η 2 G 1 + ξ 1 η 2 G 1 = dϕ 2 dϕ 2 dϕ 2 η 1 ξ 1 ( ) = η1 2 G G 2G + η 1 ξ 1 = η1 2 F (ξ 1, η 1, ϕ 2 ) η 1 ξ 1 The above derivative is positive in the Domain V > 0 because F = 1 ω + > 0 where by dots we denote terms of order one and higher with respect of canonical variables ξ, η.
24 Accourding to Chetaev theorem the equilibrium ξ 1 = η 1 = 0 of the reduced system is unstable. It follows that the equilibrium ξ 1 = η 1 = ξ 2 = η 2 = 0 of the original system is also unstable. By summarizing the results we can formulate Theorem. If the equilibrium position ξ 1 = η 1 = ξ 2 = η 2 = 0 of the canonical system with Hamiltonian H = 1 2 η2 1 + ω(ξ η 2 2) + H (3) (ξ 1, η 1, ϕ, r) is not isolated, then it is unstable.
25 Motion of a Satellite with respect to its Center of Mass Main assumptions: Satellite is a rigid body Center of mass moves in circular or elliptic orbit Motion of mass center does not depend on motion with respect to center of mass Reference frames orbital frame (axes OX, OY, OZ directed along the radius vector of the center of mass O, the transversal and normal of the orbit, respectively) principal axes frame Oxyz We use the Euler angles ψ, θ, ϕ to describe the position of the the satellite in orbital system
26 Euler angles
27 Hamiltonian of the problem The Hamiltonian of the problem reads H = p2 ψ 2 sin 2 θ + p2 θ 2 p cos θ ψ cot θ cos ψ αβp ψ sin 2 θ p θ sin ψ + αβ cos ψ sin θ + + α2 β 2 2 sin 2 θ (α 1) cos2 θ. (14) Parameters: α = C/A and A, B, C(A = B) are the principal moments of inertia (0 < α 2). β = Ω 0 /ω 0 ; ω 0 and Ω 0 are the orbital angular velocity and the projection of the angular velocity of the satellite onto its symmetry axis respectively. The parameter β can take any real value. The angle ϕ is a cyclic coordinate so in the Hamiltonian we have put p ϕ = αβ.
28 The equilibria solution (4 3α) cos θ 0 = αβ, ψ 0 = π 2, p θ 0 = 0, p ψ0 = 1 3(α 1) cos 2 θ 0. (15) Solution (15) describes the so-called conical precession of the satellite representing the steady-state rotation in which the satellite symmetry axis is perpendicular to the velocity vector of the center of mass and forms the angle θ 0 with the normal to the orbit plane.
29 The stability problem of the conical precession has been studied in detail. The stability question of the conical precession is still open 5 only for the special case α = 4/3, β = 0, when the angle θ 0 may take any constant value. In this case relations (15) describe not a unique equilibria solution, but a one parametric family of the following equilibria solutions 0 < θ 0 π 2, ψ 0 = π 2, p θ 0 = 0, p ψ0 = sin 2 θ 0. (16) 5 A. P. Markeev, A. G. Skolskii, and T. N. Chekhovskaia. Stability of the conical precession of a dynamically symmetric rigid body. Soviet Astronomy Letters, 3: , August 1977.
30 The local canonical variables q i, p i (i = 1, 2) ψ = π 2 + q 1, p ψ = sin 2 θ 0 + p 1, θ = θ 0 + q 2, p θ0 = p 2. (17) The Hamiltonian quadratic part H 2 = p sin 2 θ 0 + p sin2 θ 0 q cot θ 0 p 1 q 2 (18) The CE of the linear system with Hamiltonian H 2 λ 4 + (4 cos 2 θ 0 1)λ 2 = 0 (19) If π/3 < θ 0 π/2, then the CE has two real roots Instability. If 0 < θ 0 < π/3 then roots of the CE are pure imaginary. Matrix of the linear system is not diagonal Nonlinear stability analysis is necessary. The case θ 0 = π/3, when all roots of the CE are zero we do not consider here.
31 The linear canonical change of variables q 1 = 2 cot θ 0 x 1 + ω 1 ω sin ω y 2, q 2 = 2 cos θ 0 ω x ω y 1, p 1 = sin θ 0 x 2 + sin 2θ 0 y 1, p 2 = 1 ω ω ω x cos θ (20) 0 y 2 ω The normal form of the Hamiltonian quadratic part H 2 = 1 2 x ω(x y 2 2 ) (21) Conclusion: The conical precession is unstable for the case 0 < θ 0 < π/3
32 Other problems Three body problem. Stability of L 4 at µ = 0. Stability of regular precession of a satellite at α = 0 (spherically symmetric satellite). Stability of planar periodic pendulum-like motions of a rigid body with a fixed point in the case of Lagrange and in the case of Gorychev-Chaplygin 6. 6 Bardin, B. S. On stability problem of rigid body penduliform motions in Goraychev-Chaplygin s case Akademiia Nauk SSSR, Izvestiia, Mekhanika Tverdogo Tela, 2007, 14-21
33 Open questions What are conditions for transcendental situation nonresonant case and in cases of other resonances? Is equilibrium can be stable in the transcendental case? What are relation between transcendental situation and integrability of the system? What is dynamical behavior of the system in the neighborhood of equilibria in transcendental case?
34 Thank you for attention and Happy Birthday, Ken!!!
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
More informationLinearization at equilibrium points
TMA4165 2007 Linearization at equilibrium points Harald Hanche-Olsen hanche@math.ntnu.no This note is about the behaviour of a nonlinear autonomous system ẋ = f (x) (where x(t) R n ) near an equilibrium
More informationLectures on Dynamical Systems. Anatoly Neishtadt
Lectures on Dynamical Systems Anatoly Neishtadt Lectures for Mathematics Access Grid Instruction and Collaboration (MAGIC) consortium, Loughborough University, 2007 Part 3 LECTURE 14 NORMAL FORMS Resonances
More informationEntrance Exam, Differential Equations April, (Solve exactly 6 out of the 8 problems) y + 2y + y cos(x 2 y) = 0, y(0) = 2, y (0) = 4.
Entrance Exam, Differential Equations April, 7 (Solve exactly 6 out of the 8 problems). Consider the following initial value problem: { y + y + y cos(x y) =, y() = y. Find all the values y such that the
More informationPoint vortices in a circular domain: stability, resonances, and instability of stationary rotation of a regular vortex polygon
Point vortices in a circular domain: stability, resonances, and instability of stationary rotation of a regular vortex polygon Leonid Kurakin South Federal University Department of Mathematics, Mechanics
More informationDynamical Systems & Lyapunov Stability
Dynamical Systems & Lyapunov Stability Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University Outline Ordinary Differential Equations Existence & uniqueness Continuous dependence
More informationLecture 38: Equations of Rigid-Body Motion
Lecture 38: Equations of Rigid-Body Motion It s going to be easiest to find the equations of motion for the object in the body frame i.e., the frame where the axes are principal axes In general, we can
More informationIntroduction to Applied Nonlinear Dynamical Systems and Chaos
Stephen Wiggins Introduction to Applied Nonlinear Dynamical Systems and Chaos Second Edition With 250 Figures 4jj Springer I Series Preface v L I Preface to the Second Edition vii Introduction 1 1 Equilibrium
More informationThe 3D restricted three-body problem under angular velocity variation. K. E. Papadakis
A&A 425, 11 1142 (2004) DOI: 10.1051/0004-661:20041216 c ESO 2004 Astronomy & Astrophysics The D restricted three-body problem under angular velocity variation K. E. Papadakis Department of Engineering
More informationLecture 38: Equations of Rigid-Body Motion
Lecture 38: Equations of Rigid-Body Motion It s going to be easiest to find the equations of motion for the object in the body frame i.e., the frame where the axes are principal axes In general, we can
More informationAVERAGING AND RECONSTRUCTION IN HAMILTONIAN SYSTEMS
AVERAGING AND RECONSTRUCTION IN HAMILTONIAN SYSTEMS Kenneth R. Meyer 1 Jesús F. Palacián 2 Patricia Yanguas 2 1 Department of Mathematical Sciences University of Cincinnati, Cincinnati, Ohio (USA) 2 Departamento
More informationExistence and stability of collinear equilibrium points in elliptic restricted three body problem with radiating primary and triaxial secondary
Modelling, Measurement and Control A Vol. 9, No., March, 08, pp. -8 Journal homepage: http://iieta.org/journals/mmc/mmc_a Existence and stability of collinear equilibrium points in elliptic restricted
More informationPhys 7221 Homework # 8
Phys 71 Homework # 8 Gabriela González November 15, 6 Derivation 5-6: Torque free symmetric top In a torque free, symmetric top, with I x = I y = I, the angular velocity vector ω in body coordinates with
More informationTRANSFER TO THE COLLINEAR LIBRATION POINT L 3 IN THE SUN-EARTH+MOON SYSTEM
TRANSFER TO THE COLLINEAR LIBRATION POINT L 3 IN THE SUN-EARTH+MOON SYSTEM HOU Xi-yun,2 TANG Jing-shi,2 LIU Lin,2. Astronomy Department, Nanjing University, Nanjing 20093, China 2. Institute of Space Environment
More informationStability Implications of Bendixson s Criterion
Wilfrid Laurier University Scholars Commons @ Laurier Mathematics Faculty Publications Mathematics 1998 Stability Implications of Bendixson s Criterion C. Connell McCluskey Wilfrid Laurier University,
More informationCP1 REVISION LECTURE 3 INTRODUCTION TO CLASSICAL MECHANICS. Prof. N. Harnew University of Oxford TT 2017
CP1 REVISION LECTURE 3 INTRODUCTION TO CLASSICAL MECHANICS Prof. N. Harnew University of Oxford TT 2017 1 OUTLINE : CP1 REVISION LECTURE 3 : INTRODUCTION TO CLASSICAL MECHANICS 1. Angular velocity and
More informationStudy of the Restricted Three Body Problem When One Primary Is a Uniform Circular Disk
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 93-9466 Applications and Applied Mathematics: An International Journal (AAM) Vol. 3, Issue (June 08), pp. 60 7 Study of the Restricted Three Body
More informationThere is a more global concept that is related to this circle of ideas that we discuss somewhat informally. Namely, a region R R n with a (smooth)
82 Introduction Liapunov Functions Besides the Liapunov spectral theorem, there is another basic method of proving stability that is a generalization of the energy method we have seen in the introductory
More informationINVARIANT TORI IN THE LUNAR PROBLEM. Kenneth R. Meyer, Jesús F. Palacián, and Patricia Yanguas. Dedicated to Jaume Llibre on his 60th birthday
Publ. Mat. (2014), 353 394 Proceedings of New Trends in Dynamical Systems. Salou, 2012. DOI: 10.5565/PUBLMAT Extra14 19 INVARIANT TORI IN THE LUNAR PROBLEM Kenneth R. Meyer, Jesús F. Palacián, and Patricia
More informationPart II. Classical Dynamics. Year
Part II Year 28 27 26 25 24 23 22 21 20 2009 2008 2007 2006 2005 28 Paper 1, Section I 8B Derive Hamilton s equations from an action principle. 22 Consider a two-dimensional phase space with the Hamiltonian
More informationAngular momentum. Quantum mechanics. Orbital angular momentum
Angular momentum 1 Orbital angular momentum Consider a particle described by the Cartesian coordinates (x, y, z r and their conjugate momenta (p x, p y, p z p. The classical definition of the orbital angular
More information1.7. Stability and attractors. Consider the autonomous differential equation. (7.1) ẋ = f(x),
1.7. Stability and attractors. Consider the autonomous differential equation (7.1) ẋ = f(x), where f C r (lr d, lr d ), r 1. For notation, for any x lr d, c lr, we let B(x, c) = { ξ lr d : ξ x < c }. Suppose
More information1/30. Rigid Body Rotations. Dave Frank
. 1/3 Rigid Body Rotations Dave Frank A Point Particle and Fundamental Quantities z 2/3 m v ω r y x Angular Velocity v = dr dt = ω r Kinetic Energy K = 1 2 mv2 Momentum p = mv Rigid Bodies We treat a rigid
More informationIn Arnold s Mathematical Methods of Classical Mechanics (1), it
Near strongly resonant periodic orbits in a Hamiltonian system Vassili Gelfreich* Mathematics Institute, University of Warwick, Coventry CV4 7AL, United Kingdom Communicated by John N. Mather, Princeton
More informationEulerian equilibria of a rigid body in the three body problem
Proceedings of the nd WSEAS Int. Conference on Applied and Theoretical Mechanics, Venice, Italy, November 0-, 006 68 Eulerian equilibria of a rigid body in the three body problem J.A. VERA Universidad
More informationProf. Krstic Nonlinear Systems MAE281A Homework set 1 Linearization & phase portrait
Prof. Krstic Nonlinear Systems MAE28A Homework set Linearization & phase portrait. For each of the following systems, find all equilibrium points and determine the type of each isolated equilibrium. Use
More informationNonlinear Control of Electrodynamic Tether in Equatorial or Somewhat Inclined Orbits
Proceedings of the 15th Mediterranean Conference on Control & Automation, July 7-9, 7, Athens - Greece T-5 Nonlinear Control of Electrodynamic Tether in Equatorial or Somewhat Inclined Martin B. Larsen
More informationMoment of inertia. Contents. 1 Introduction and simple cases. January 15, Introduction. 1.2 Examples
Moment of inertia January 15, 016 A systematic account is given of the concept and the properties of the moment of inertia. Contents 1 Introduction and simple cases 1 1.1 Introduction.............. 1 1.
More informationHandout 6: Rotational motion and moment of inertia. Angular velocity and angular acceleration
1 Handout 6: Rotational motion and moment of inertia Angular velocity and angular acceleration In Figure 1, a particle b is rotating about an axis along a circular path with radius r. The radius sweeps
More informationDYNAMICAL SYSTEMS WITH A CODIMENSION-ONE INVARIANT MANIFOLD: THE UNFOLDINGS AND ITS BIFURCATIONS
International Journal of Bifurcation and Chaos c World Scientific Publishing Company DYNAMICAL SYSTEMS WITH A CODIMENSION-ONE INVARIANT MANIFOLD: THE UNFOLDINGS AND ITS BIFURCATIONS KIE VAN IVANKY SAPUTRA
More informationSOLUTIONS TO THE GINZBURG LANDAU EQUATIONS FOR PLANAR TEXTURES IN SUPERFLUID 3 He
SOLUTIONS TO THE GINZBURG LANDAU EQUATIONS FOR PLANAR TEXTURES IN SUPERFLUID 3 He V. L. GOLO, M. I. MONASTYRSKY, AND S. P. NOVIKOV Abstract. The Ginzburg Landau equations for planar textures of superfluid
More informationHamiltonian Dynamics from Lie Poisson Brackets
1 Hamiltonian Dynamics from Lie Poisson Brackets Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 12 February 2002 2
More informationDifferential Kinematics
Differential Kinematics Relations between motion (velocity) in joint space and motion (linear/angular velocity) in task space (e.g., Cartesian space) Instantaneous velocity mappings can be obtained through
More informationANALYTICAL MECHANICS. LOUIS N. HAND and JANET D. FINCH CAMBRIDGE UNIVERSITY PRESS
ANALYTICAL MECHANICS LOUIS N. HAND and JANET D. FINCH CAMBRIDGE UNIVERSITY PRESS Preface xi 1 LAGRANGIAN MECHANICS l 1.1 Example and Review of Newton's Mechanics: A Block Sliding on an Inclined Plane 1
More informationThe phenomenon of reversal in the Euler Poincaré Suslov nonholonomic systems
The phenomenon of reversal in the Euler Poincaré Suslov nonholonomic systems Valery V.Kozlov arxiv:1509.01089v1 [nlin.si] 3 Sep 2015 Steklov Mathematical Institute, Russian Academy of Sciences Gubkina
More informationPhysical Dynamics (SPA5304) Lecture Plan 2018
Physical Dynamics (SPA5304) Lecture Plan 2018 The numbers on the left margin are approximate lecture numbers. Items in gray are not covered this year 1 Advanced Review of Newtonian Mechanics 1.1 One Particle
More informationStabilization and Passivity-Based Control
DISC Systems and Control Theory of Nonlinear Systems, 2010 1 Stabilization and Passivity-Based Control Lecture 8 Nonlinear Dynamical Control Systems, Chapter 10, plus handout from R. Sepulchre, Constructive
More informationClassical Mechanics III (8.09) Fall 2014 Assignment 3
Classical Mechanics III (8.09) Fall 2014 Assignment 3 Massachusetts Institute of Technology Physics Department Due September 29, 2014 September 22, 2014 6:00pm Announcements This week we continue our discussion
More information2 Lyapunov Stability. x(0) x 0 < δ x(t) x 0 < ɛ
1 2 Lyapunov Stability Whereas I/O stability is concerned with the effect of inputs on outputs, Lyapunov stability deals with unforced systems: ẋ = f(x, t) (1) where x R n, t R +, and f : R n R + R n.
More informationTHE CELT RATTLEBACK DYNAMICS WITH THE FRICTION INFLUENCE Institute of Mechanical Engineering Problems RAS
THE CELT RATTLEBACK DYNAMICS WITH THE FRICTION INFLUENCE Institute of Mechanical Engineering Problems RAS e-mail . Introduction. The Celt stone dynamics investigation is the classical
More informationPhysical Dynamics (PHY-304)
Physical Dynamics (PHY-304) Gabriele Travaglini March 31, 2012 1 Review of Newtonian Mechanics 1.1 One particle Lectures 1-2. Frame, velocity, acceleration, number of degrees of freedom, generalised coordinates.
More informationKinematics (special case) Dynamics gravity, tension, elastic, normal, friction. Energy: kinetic, potential gravity, spring + work (friction)
Kinematics (special case) a = constant 1D motion 2D projectile Uniform circular Dynamics gravity, tension, elastic, normal, friction Motion with a = constant Newton s Laws F = m a F 12 = F 21 Time & Position
More informationLecture 41: Highlights
Lecture 41: Highlights The goal of this lecture is to remind you of some of the key points that we ve covered this semester Note that this is not the complete set of topics that may appear on the final
More informationDynamics and Stability application to submerged bodies, vortex streets and vortex-body systems
Dynamics and Stability application to submerged bodies, vortex streets and vortex-body systems Eva Kanso University of Southern California CDS 140B Introduction to Dynamics February 5 and 7, 2008 Fish
More information27. Euler s Equations
27 Euler s Equations Michael Fowler Introduction We ve just seen that by specifying the rotational direction and the angular phase of a rotating body using Euler s angles, we can write the Lagrangian in
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 information154 Chapter 9 Hints, Answers, and Solutions The particular trajectories are highlighted in the phase portraits below.
54 Chapter 9 Hints, Answers, and Solutions 9. The Phase Plane 9.. 4. The particular trajectories are highlighted in the phase portraits below... 3. 4. 9..5. Shown below is one possibility with x(t) and
More informationA qualitative analysis of bifurcations to halo orbits
1/28 spazio A qualitative analysis of bifurcations to halo orbits Dr. Ceccaroni Marta ceccaron@mat.uniroma2.it University of Roma Tor Vergata Work in collaboration with S. Bucciarelli, A. Celletti, G.
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 informationLinear ODEs. Existence of solutions to linear IVPs. Resolvent matrix. Autonomous linear systems
Linear ODEs p. 1 Linear ODEs Existence of solutions to linear IVPs Resolvent matrix Autonomous linear systems Linear ODEs Definition (Linear ODE) A linear ODE is a differential equation taking the form
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 informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 4. Bifurcations Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Local bifurcations for vector fields 1.1 The problem 1.2 The extended centre
More informationRigid bodies - general theory
Rigid bodies - general theory Kinetic Energy: based on FW-26 Consider a system on N particles with all their relative separations fixed: it has 3 translational and 3 rotational degrees of freedom. Motion
More informationNonlinear systems. Lyapunov stability theory. G. Ferrari Trecate
Nonlinear systems Lyapunov stability theory G. Ferrari Trecate Dipartimento di Ingegneria Industriale e dell Informazione Università degli Studi di Pavia Advanced automation and control Ferrari Trecate
More informationLyapunov Stability Theory
Lyapunov Stability Theory Peter Al Hokayem and Eduardo Gallestey March 16, 2015 1 Introduction In this lecture we consider the stability of equilibrium points of autonomous nonlinear systems, both in continuous
More informationChaotic motion. Phys 420/580 Lecture 10
Chaotic motion Phys 420/580 Lecture 10 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
More informationNumeric Deduction in Symbolic Computation. Application to Normalizing Transformations
arxiv:1607.02016v1 [cs.sc] 27 May 2016 Numeric Deduction in Symbolic Computation. Application to Normalizing Transformations Ivan I. Shevchenko Institute of Theoretical Astronomy, Russian Academy of Sciences,
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 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 information415. On evolution of libration points similar to Eulerian in the model problem of the binary-asteroids dynamics
45 On evolution of libration points similar to Eulerian in the model problem of the binary-asteroids dynamics V V Beletsky and A V Rodnikov Keldysh Institute of Applied Mathematics of Russian Academy of
More informationCDS 101 Precourse Phase Plane Analysis and Stability
CDS 101 Precourse Phase Plane Analysis and Stability Melvin Leok Control and Dynamical Systems California Institute of Technology Pasadena, CA, 26 September, 2002. mleok@cds.caltech.edu http://www.cds.caltech.edu/
More informationCeres Rotation Solution under the Gravitational Torque of the Sun
Ceres Rotation Solution under the Gravitational Torque of the Sun Martin Lara, Toshio Fukushima, Sebastián Ferrer (*) Real Observatorio de la Armada, San Fernando, Spain ( ) National Astronomical Observatory,
More informationON COLISSIONS IN NONHOLONOMIC SYSTEMS
ON COLISSIONS IN NONHOLONOMIC SYSTEMS DMITRY TRESCHEV AND OLEG ZUBELEVICH DEPT. OF THEORETICAL MECHANICS, MECHANICS AND MATHEMATICS FACULTY, M. V. LOMONOSOV MOSCOW STATE UNIVERSITY RUSSIA, 119899, MOSCOW,
More information4 Second-Order Systems
4 Second-Order Systems Second-order autonomous systems occupy an important place in the study of nonlinear systems because solution trajectories can be represented in the plane. This allows for easy visualization
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 informationSymplectic maps. James D. Meiss. March 4, 2008
Symplectic maps James D. Meiss March 4, 2008 First used mathematically by Hermann Weyl, the term symplectic arises from a Greek word that means twining or plaiting together. This is apt, as symplectic
More informationB. Differential Equations A differential equation is an equation of the form
B Differential Equations A differential equation is an equation of the form ( n) F[ t; x(, xʹ (, x ʹ ʹ (, x ( ; α] = 0 dx d x ( n) d x where x ʹ ( =, x ʹ ʹ ( =,, x ( = n A differential equation describes
More informationSPACECRAFT DYNAMICS NEAR A BINARY ASTEROID. F. Gabern, W.S. Koon and J.E. Marsden
PROCEEDINGS OF THE FIFTH INTERNATIONAL CONFERENCE ON DYNAMICAL SYSTEMS AND DIFFERENTIAL EQUATIONS June 16 19, 2004, Pomona, CA, USA pp. 1 10 SPACECRAFT DYNAMICS NEAR A BINARY ASTEROID F. Gabern, W.S. Koon
More informationPHY 5246: Theoretical Dynamics, Fall Assignment # 10, Solutions. (1.a) N = a. we see that a m ar a = 0 and so N = 0. ω 3 ω 2 = 0 ω 2 + I 1 I 3
PHY 54: Theoretical Dynamics, Fall 015 Assignment # 10, Solutions 1 Graded Problems Problem 1 x 3 a ω First we calculate the moments of inertia: ( ) a I 1 = I = m 4 + b, 1 (1.a) I 3 = ma. b/ α The torque
More informationPERIODIC ORBITS OF A COLLINEAR RESTRICTED THREE BODY PROBLEM
PERIODIC ORBITS OF A COLLINEAR RESTRICTED THREE BODY PROBLEM MONTSERRAT CORBERA Departament d Informàtica i Matemàtica, Escola Politècnica Superior, Universitat de Vic, C/ Laura 13, 85 Vic, Barcelona,
More informationOscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. Oscillatory solution
Oscillatory Motion Simple pendulum: linear Hooke s Law restoring force for small angular deviations d 2 θ dt 2 = g l θ θ l Oscillatory solution θ(t) =θ 0 sin(ωt + φ) F with characteristic angular frequency
More informationOscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. small angle approximation. Oscillatory solution
Oscillatory Motion Simple pendulum: linear Hooke s Law restoring force for small angular deviations d 2 θ dt 2 = g l θ small angle approximation θ l Oscillatory solution θ(t) =θ 0 sin(ωt + φ) F with characteristic
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 informationHomework 1. Due Tuesday, January 29.
Homework 1. Due Tuesday, January 29. Problem 1. An ideal rope (no friction) lying on a table slides from its edge down to a scales lying on the floor. The table s height is h. Find a stationary velocity
More information[#1] R 3 bracket for the spherical pendulum
.. Holm Tuesday 11 January 2011 Solutions to MSc Enhanced Coursework for MA16 1 M3/4A16 MSc Enhanced Coursework arryl Holm Solutions Tuesday 11 January 2011 [#1] R 3 bracket for the spherical pendulum
More informationNumerical Continuation and Normal Form Analysis of Limit Cycle Bifurcations without Computing Poincaré Maps
Numerical Continuation and Normal Form Analysis of Limit Cycle Bifurcations without Computing Poincaré Maps Yuri A. Kuznetsov joint work with W. Govaerts, A. Dhooge(Gent), and E. Doedel (Montreal) LCBIF
More informationChemistry 532 Practice Final Exam Fall 2012 Solutions
Chemistry 53 Practice Final Exam Fall Solutions x e ax dx π a 3/ ; π sin 3 xdx 4 3 π cos nx dx π; sin θ cos θ + K x n e ax dx n! a n+ ; r r r r ˆL h r ˆL z h i φ ˆL x i hsin φ + cot θ cos φ θ φ ) ˆLy i
More informationConnecting orbits and invariant manifolds in the spatial three-body problem
C C Dynamical A L T E C S H Connecting orbits and invariant manifolds in the spatial three-body problem Shane D. Ross Control and Dynamical Systems, Caltech Work with G. Gómez, W. Koon, M. Lo, J. Marsden,
More informationIm + α α. β + I 1 I 1< 0 I 1= 0 I 1 > 0
ON THE HAMILTONIAN ANDRONOV-HOPF BIFURCATION M. Olle, J. Villanueva 2 and J. R. Pacha 3 2 3 Departament de Matematica Aplicada I (UPC), Barcelona, Spain In this contribution, we consider an specic type
More informationISOLATING BLOCKS NEAR THE COLLINEAR RELATIVE EQUILIBRIA OF THE THREE-BODY PROBLEM
ISOLATING BLOCKS NEAR THE COLLINEAR RELATIVE EQUILIBRIA OF THE THREE-BODY PROBLEM RICHARD MOECKEL Abstract. The collinear relative equilibrium solutions are among the few explicitly known periodic solutions
More informationNBA Lecture 1. Simplest bifurcations in n-dimensional ODEs. Yu.A. Kuznetsov (Utrecht University, NL) March 14, 2011
NBA Lecture 1 Simplest bifurcations in n-dimensional ODEs Yu.A. Kuznetsov (Utrecht University, NL) March 14, 2011 Contents 1. Solutions and orbits: equilibria cycles connecting orbits other invariant sets
More informationDYNAMICS OF GENERALIZED EULER TOPS WITH CONSTRAINTS. Dmitry V. Zenkov, Anthony M. Bloch
PROCEEDINGS OF THE INTERNATIONAL CONFERENCE ON DYNAMICAL SYSTEMS AND DIFFERENTIAL EQUATIONS May 18 21, 2000, Atlanta, USA pp. 398 405 DYNAMICS OF GENERALIZED EULER TOPS WITH CONSTRAINTS Dmitry V. Zenkov,
More informationNow I switch to nonlinear systems. In this chapter the main object of study will be
Chapter 4 Stability 4.1 Autonomous systems Now I switch to nonlinear systems. In this chapter the main object of study will be ẋ = f(x), x(t) X R k, f : X R k, (4.1) where f is supposed to be locally Lipschitz
More informationPart II. Dynamical Systems. Year
Part II Year 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2017 34 Paper 1, Section II 30A Consider the dynamical system where β > 1 is a constant. ẋ = x + x 3 + βxy 2, ẏ = y + βx 2
More informationMathematical Foundations of Neuroscience - Lecture 7. Bifurcations II.
Mathematical Foundations of Neuroscience - Lecture 7. Bifurcations II. Filip Piękniewski Faculty of Mathematics and Computer Science, Nicolaus Copernicus University, Toruń, Poland Winter 2009/2010 Filip
More information7 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
More informationInstability of Periodic Orbits of Some Rhombus and Parallelogram Four Body Problems
Instability of Periodic Orbits of Some Rhombus and Parallelogram Four Body Problems by Abdalla Mansur A thesis submitted to the Department of Mathematics & Statistics in conformity with the requirements
More informationDynamics of a rigid rotor in the elastic bearings
Theoret. Appl. Mech., Vol.31, No.1, pp.73-83, Belgrade 2004 Dynamics of a rigid rotor in the elastic bearings Inga M. Arkhipova Abstract As a rule in the studies of a rigid rotor in the elastic bearings
More informationSolution of Additional Exercises for Chapter 4
1 1. (1) Try V (x) = 1 (x 1 + x ). Solution of Additional Exercises for Chapter 4 V (x) = x 1 ( x 1 + x ) x = x 1 x + x 1 x In the neighborhood of the origin, the term (x 1 + x ) dominates. Hence, the
More informationBIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs
BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits. 2. Equilibria. 3. Periodic orbits and limit cycles. 4. Homoclinic orbits.
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 37: Principal Axes, Translations, and Eulerian Angles
Lecture 37: Principal Axes, Translations, and Eulerian Angles When Can We Find Principal Axes? We can always write down the cubic equation that one must solve to determine the principal moments But if
More informationSteering the Chaplygin Sleigh by a Moving Mass
Steering the Chaplygin Sleigh by a Moving Mass Jason M. Osborne Department of Mathematics North Carolina State University Raleigh, NC 27695 jmosborn@unity.ncsu.edu Dmitry V. Zenkov Department of Mathematics
More informationWIDELY SEPARATED FREQUENCIES IN COUPLED OSCILLATORS WITH ENERGY-PRESERVING QUADRATIC NONLINEARITY
WIDELY SEPARATED FREQUENCIES IN COUPLED OSCILLATORS WITH ENERGY-PRESERVING QUADRATIC NONLINEARITY J.M. TUWANKOTTA Abstract. In this paper we present an analysis of a system of coupled oscillators suggested
More informationFINAL EXAM MATH303 Theory of Ordinary Differential Equations. Spring dx dt = x + 3y dy dt = x y.
FINAL EXAM MATH0 Theory of Ordinary Differential Equations There are 5 problems on 2 pages. Spring 2009. 25 points Consider the linear plane autonomous system x + y x y. Find a fundamental matrix of the
More informationQualifying Exam. Aug Part II. Please use blank paper for your work do not write on problems sheets!
Qualifying Exam Aug. 2015 Part II Please use blank paper for your work do not write on problems sheets! Solve only one problem from each of the four sections Mechanics, Quantum Mechanics, Statistical Physics
More informationHamiltonian Dynamics In The Theory of Abstraction
Hamiltonian Dynamics In The Theory of Abstraction Subhajit Ganguly. email: gangulysubhajit@indiatimes.com Abstract: This paper deals with fluid flow dynamics which may be Hamiltonian in nature and yet
More information13. Rigid Body Dynamics II
University of Rhode Island DigitalCommons@URI Classical Dynamics Physics Course Materials 2015 13. Rigid Body Dynamics II Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative Commons License
More informationOn the influence of sliding on the Celt rattleback motion
432 Proceedings of XXXV International Summer School-Conference APM 2007 On the influence of sliding on the Celt rattleback motion Tatiana P. Tovstik tptovstik@mail.ru Abstract The problem of the Celt stone
More informationChaotic 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
More information