Lagrange-d Alembert integrators for constrained systems in mechanics
|
|
- Arabella Lawson
- 5 years ago
- Views:
Transcription
1 Lagrange-d Alembert integrators for constrained systems in mechanics Conference on Scientific Computing in honor of Ernst Hairer s 6th birthday, Geneva, Switzerland Laurent O. Jay Dept. of Mathematics, The University of Iowa, USA June 7-2, 29
2 Lagrangian systems with nonholonomic constraints Coordinates and velocities: (q, v) Lagrangian: L(q, v) Nonholonomic constraints: k(q, v) = Lagrange equations of the second kind: d dt q = v d dt vl(q,v) = q L(q,v) K(q,v) T ψ = k(q,v) They form a system of differential-algebraic equations (DAEs) Index 2 assumption: ( 2 J(q,v) := vv L(q,v) K(q,v) T ) nonsingular k v (q,v) O Ideal nonholonomic constraints when K(q,v) = k v (q,v)
3 Applications with nonholonomic constraints Rolling balls, wheels, and disks Skating and sleighs Trajectory planning for vehicles, e.g., car trailer Robotic mechanisms, e.g., wheeled robots Molecular dynamics simulations with constant kinetic temperature Missile pursuit etc.
4 The underlying forced Lagrangian system where d dt q = v d dt vl(q,v) = q L(q,v) + f L (q,v) f L (q,v) := K(q,v) T Ψ(q,v) and Ψ(q,v) := ψ is the solution with a := d dt v of the system of linear equations ( ) ( ) a 2 J(q,v) = qv L(q,v)v + q L(q,v) ψ k q (q,v)v The nonholonomic constraints = k(q,v) are strong invariants of the underlying forced Lagrangian system.
5 The Lagrange-d Alembert principle The index 2 DAEs can be derived from the Lagrange-d Alembert principle ( d ) dt L v(q, q) + L q (q, q) δ q = and the nonholonomic constraints = k(q, q) where δ q are virtual displacements satisfying the Maurer-Appell-Chetaev- Johnsen-Hamel rule k v (q, q)δ q =
6 Equivalence to a skew critical problem when k(q, v) = K(q)v Define the action A(q) := tn For nonholonomic constraints linear in v t L(q(t), q(t))dt = K(q)v the Lagrange-d Alembert principle is equivalent to a skew critical problem This motivates... δa(q)(δq) = δq C K(q)δq
7 Cortés discrete Lagrange-d Alembert principle (2) Assume that where Consider a discrete Lagrangian = c(q k,q k+ ) for k =,...,N c(q k, q k+ ) := K(q k+ )v k+ (q k,q k+ ) {z } v(t k+ ) L d (q k, q k+ ) Z tk+ t k L(q(t), q(t))dt where q(t k ) = q k, q(t k+ ) = q k+. Cortés discrete Lagrange-d Alembert principle is defined as follows: the discrete action sum N X L d (q k, q k+ ) k= is stationary with respect to variations {δq k } N k= δq k Ker(K(q k )) R n (K(q k )δq k = ) and δq = = δq N. This leads to Cortés discrete Euler-Lagrange equations 2 L d (q k,q k ) + L d (q k,q k+ ) = K(q k ) T ψ k
8 Example of a LDA method For L(q,v) := 2 vt Mv U(q), k(q,v) := K(q)v + b(q) M sym. pos. def., McLachlan-Perlmutter s 2-stage (Lobatto IIIB-A) LDA method (26): Q = q + h 2 v Mv = Mv h( U(Q ) + K(Q ) T Ψ ) q = q + h 2 (v + v ) = Q + h 2 v = K(q )v + b(q )
9 Quote from McLachlan and Perlmutter (26)...much work remains to be done to clarify the nature of discrete nonholonomic mechanics and to pinpoint the correct discrete analog of the Lagrange-d Alembert principle.
10 Forced Lagrangian systems d dt q = v d dt vl(q,v) = q L(q,v) + f L (q,v) the integral Lagrange-d Alembert principle is tn δa(q)(δq) + f L (q(t), q(t)) T δq(t)dt = δq C t Idea: For nonholonomic Lagrangian systems take f L (q,v) := K(q,v) T Ψ(q,v) and consider the nonholonomic constraints = k(q,v) as strong invariants of the underlying forced Lagrangian system.
11 Notation Let q(t) satisfy q(t k ) = q k, q(t k+ ) = q k+. Consider a discrete Lagrangian L d (q k,q k+ ) approximating the exact discrete Lagrangian L E d (q k,q k+ ) := tk+ t k L(q(t), q(t))dt Consider discrete forces f d (q k,q k+ ),f + d (q k,q k+ ) approximating the exact discrete forcing terms f E d (q k,q k+ ) T := f E+ d (q k,q k+ ) T := tk+ t k tk+ t k f L (q(t), q(t)) T qk q(t)dt f L (q(t), q(t)) T qk+ q(t)dt
12 The discrete Lagrange-d Alembert principle of C. Kane, Marsden, Ortiz, and West (2) for forced Lagrangian systems N L d (q k,q k+ )+ N δ k= k= ( f d (q k,q k+ ) T δq k + f + d (q k,q k+ ) T δq k+ ) = {δq k } N k= δq k R n and δq = = δq N. This leads to the forced discrete Euler-Lagrange equations 2 L d (q k,q k ) + L d (q k,q k+ ) + f + d (q k,q k ) + f d (q k,q k+ ) = for k =,...,N and which are satisfied by the exact discrete Lagrangian and the exact discrete forcing terms.
13 A discrete Lagrange-d Alembert principle for Lagrangian systems with nonholonomic constraints Combine the forced discrete principle Lagrange-d Alembert principle with f L (q,v) := K(q,v) T Ψ(q,v) and the constraints = c(q k,q k+ ) = k(q k+,v k+ (q k,q k+ )) for k =,...,N }{{} v(t k+ ) This leads to the nonholonomically constrained discrete Euler-Lagrange equations 2 L d (q k,q k ) + L d (q k,q k+ ) + f + d (q k,q k ) + f d (q k,q k+ ) = c(q k,q k ) = c(q k,q k+ ) for k =,...,N where c(q,q ) = or c(q N,q N ) =.
14 The exact discrete forcing terms One can show that L E d (q k,q k+ ) = p k f E d (q k,q k+ ) 2 L E d (q k,q k+ ) = p k+ f E+ d (q k,q k+ ) where p k := v L(q k,v k (q k,q k+ )), p k+ := v L(q k+,v k+ (q k,q k+ )) and t f E d (q k,q k+ ) T := Ψ(q(t),v(t)) T K(q(t),v(t)) qk q(t)dt f E+ d (q k,q k+ ) T := t t The exact terms L E d,f E d,f E+ d satisfy the new principle do not satisfy Cortés principle t Ψ(q(t),v(t)) T K(q(t),v(t)) qk+ q(t)dt
15 Specialized Partitioned Additive Runge-Kutta (SPARK) methods for index 2 DAEs (J., BIT, 23; Math. Comput., 26) Q i = q + h sx a ij V j i =,...,s j= sx q = q + h b j V j j= sx sx vl(q i,v i ) = vl(q, v ) + h ba ij ql(q j,v j ) h ea ij K(Q j, V j ) T Ψ j i =,...,s j= j= sx sx vl(q, v ) = vl(q, v ) + h b bj ql(q j,v j ) h e bj K(Q j, V j ) T Ψ j j= j= sx = b j c i j k(q j,v j ) i =,..., s j= = k(q,v )
16 Examples of SPARK methods For L(q, v) := 2 vt Mv U(q), M sym. pos. def., k(q, v) := K(q)v + b(q) The -stage Gauss SPARK (midpoint) method: q = q + h 2 (v + v ) Mv = ««q + q q + q T Mv h U hk Ψ 2 2 = K(q )v + b(q ) The 2-stage Lobatto IIIA-B SPARK (Verlet) method: MV = Mv h 2 ( U(q ) + K(q ) T Ψ ) q = q + hv or modified = K = 2 (K(q )V + b(q )) + 2 (K(q )V + b(q )) q + h «2 V V + b q + h «2 V Mv + K(q ) T h 2 Ψ 2 «= MV h 2 U(q ) K(q )v = b(q ) but no more SPARK
17 Main Theorem Let t,q,t,q be given. If b i = b i, bi a ij + b j â ji b i b j = for i,j =,...,s, then we have a nonholonomically constrained discrete Lagrange-d Alembert integrator in the new sense with L d (q,q ) = h s b i L(Q i,v i ) i= f d (q,q ) T = h f + d (q,q ) T = h s i= s i= b i b i h h s ã ij Ψ T j K(Q j,v j ) q V i j= s ( b j ã ij )Ψ T j K(Q j,v j ) q V i j=
18 Main Theorem (cont.) Moreover, if c s+ i + c i = a s+ i,s+ j + a ij â s+ i,s+ j + â ij ã s+ i,s+ j + ã ij = b s+ j = b j = b s+ j = b j = b s+ j = b j for i,j =,...,s then the SPARK method is symmetric and L d (q,q ) = L d (q,q ) f d (q,q ) = f + d (q,q ) f + d (q,q ) = f d (q,q ) as for the exact solution q(t ),q(t ) and exact terms L E E d,fd,f E+ d.
19 The exact discrete forcing terms From q(t) = q R t v(s)ds we have t q q(t) = R t t q v(s)ds. Hence, Similarly f E d (q, q ) T = = = = f E+ d (q, q ) T = Z t Z t «Ψ(q(t), v(t)) T K(q(t), v(t)) q v(s)ds dt t t Z t Z t «Ψ(q(t), v(t)) T K(q(t), v(t)) q v(s)ds dt t t Z t Z s «Ψ(q(t), v(t)) T K(q(t), v(t)) q v(s)dt ds t t Z t Z s «Ψ(q(t), v(t)) T K(q(t), v(t))dt t t q v(s)ds Z t Z t t s «Ψ(q(t), v(t)) T K(q(t), v(t))dt q v(s)ds
20 Remark For forced Lagrangians Kane, Marsden, Ortiz, and West (2) originally proposed to take as a reasonable choice f d (q,q ) T = h corresponding here to s b i f L (Q i,v i ) T q Q i i= f d (q,q ) T = h = h s b i Ψ T i K(Q i,v i ) q Q i i= s b i Ψ T i K(Q i,v i ) I + h i= j= s a ij q V j which is not equivalent
21 Example: The nonholonomic particle L(q,v) = 2 (v2 + v2 2 + v2 3 ) (q2 + q2 2 ) = v 3 q 2 v = ( ) T q v = ( ) T t = h =.2
22 Example: The nonholonomic particle (cont.) 8 x 3 stage Gauss SPARK 6 E k E 4 2 E k E t x 2 stage Gauss SPARK t 2 x 3 stage Gauss SPARK 7.5 E k E t
23 Example: The skate on an inclined plane L(q,v) = m 2 (v2 + v 2 2) + J 2 v2 3 + mg sin(α)q = cos(q 3 )v 2 sin(q 3 )v m =, J =, g sin(α) = = ( ) T q v = ( ) T t = h =.
24 Example: The skate on an inclined plane (cont.) x 3 2 stage Lobatto IIIA B SPARK E k E t 8 x 3 stage Lobatto IIIA B SPARK 7 6 E k E t x 4 stage Lobatto IIIA B SPARK E k E t
25 Example: A mobile robot L(q,v) = m 2 (v2 + v2 2 ) + J 2 v J ωv 2 4 sin(q 4) = v lcos(q 3 )v 4 = v 2 lsin(q 3 )v 4 m =, J =, J ω =, l = = ( ) T q v = ( ) T t = h =.2
26 Example: A mobile robot (cont.). 2 stage Lobatto IIIA B D SPARK E k E t 6 x 3 stage Lobatto IIIA B D SPARK 4 4 E k E 2 E k E t x 4 stage Lobatto IIIA B D SPARK t
27 Example: McLachlan and Perlmutter s particles T(v) = 2 v 2 2, n = 2m + (m 2) U(q) = m! q q2 m+2 q2 m+3 + X q+i 2 q2 m++i i= nx = v + q i v i t = i=m+2 q j = ` cos(jπ/(2j)) T v j = ` sin(jπ/(2j)) T E = 3.6 h =.5 j =,..., J j =,..., J Let µ((e k E ) 2 ) := (J+) P J j= (E jk E ) 2 and J + =
28 Lagrange-d Alembert integrators for constrained systems in mechanics.k 2 µ((e E ) )/h 4 Example: McLachlan and Perlmutter s particles (cont.) 2 stage Lobatto IIIA B SPARK k 2 µ((e E ) )/h 4 x modified 2 stage Lobatto IIIA B SPARK k 2 µ((e E ) )/h 4 x McLachlan Perlmutter 2 stage Lobatto IIIB A 5 We observe.5 h t 3 µ((e k E )2 ) O(tk ) or equivalently like for a random walk x p µ((e k E )2 ) O(h2 tk )
29 Constrained systems in mechanics is a system of DAEs where d dt y = v(y,z) d p(y, z) dt = f (y, z, ψ) + r(y, λ) = g(y) y : coordinates v : velocities p : momenta f,r : forces = g y(y)v(y, z) = k(y, z) g : holonomic constraints λ : Lagrange multipliers associated to g k : nonholonomic constraints ψ : Lagrange multipliers associated to k Lagrange s mechanics: v(y, z) z Hamilton s mechanics: p(y, z) z
30 Constrained systems in mechanics (cont.) This system of DAEs is: overdetermined implicit partitioned additive (split) Manifold of constraints: M := {(y,z) R ny R nz = g(y), = g y(y)v(y, z), = k(y, z)} Assumptions: (q,v ) M, v,p, f,... sufficiently differentiable, and invertibility of pz r λ g yv z O «, = existence and uniqueness of a pz r λ f ψ g yv z O O k z O O A ()
31 SPARK methods Y i = y + h sx α ij v(y j,z j ) i =,...,s + j= y = Y s+ sx esx p(y i,z i ) = p(y, z ) + h bα ij f (Y j,z j, Ψ j ) + h eα ij r( Y e j, Λ e j ) i =,..., s + z = Z s+ ey i = y + h j= sx a ij v(y j,z j ) j= i =,,..., es = g( e Y i ) i =,,..., es ( e Y = y ) = g(y ) (y = e Y es ) = g y(y )v(y, z ) = sx b j c i j k(y j,z j ) i =,..., s j= = k(y, v ) j=
32 Some definitions α := A b T «, bα := ba b b T eq := ` I s s + AM ½ s M := b T b T b T A. b T (s )b T C s 2 A!, eα := b T A C A D(s ) = ea e b T!, ½ s := b T b T C. b T C s C A C A, C := diag(c,..., c s), with simplifying assumption D(s ) : kb T C k A = b T b T C k k =,..., s Q := ˇQ := Q eq γ T!, γ T := γ s+ b T A, Ǎ! a a s s T, Ǎ := B.. C. A a s a ss where we assume es = s
33 Equivalent reformulation of SPARK methods p(y,z ) = (Q I nz ) BB. p(y A ½ s+ p(y,z ) s,z s) p(y,z ) f (Y,Z,Ψ ) r( Y e, e Λ ) B C h(bα I nz A h(eα I nz ) CC. AA f (Y s,z s,ψ s) r( Y e s, Λ e s) h g y + h P s j= a jv(y j,z j ) = (ˇQ I nλ ). h g y + h P s j= a C sjv(y j,z j ) A g y(y )v(y,z ) k(y,z ) = ( Q e I nψ ) B. k(y A s,z s) k(y, z )
34 Modified Newton iterations Modified Jacobian of the reformulated SPARK methods I s+ I ny O O O O Q p z hq eα r λ hq bα f ψ O Q g yv z O O O Q e kz O O For the modified Newton iterations to obtain a block diagonal linear system with the 2 matrices () we define the intermediate quantities C A z.. z s z s+ ψ.. ψ s Z C A := (Q Inz ) Z s z C A := h( Q e bα I nψ ) C A, Ψ Ψ s C A λ.. λ s λ s+ C A := h(q eα In ) λ e Λ e Λ. Λ e s C A
35 Modified Newton iterations (cont.) We assume and we use Ψ. Ψ s Z Z s z e Λ e Λ. Λ e s eq bα Q bα = T s «C A = h (( Q e bα) I nψ ) ψ.. ψ s z C A = (Q I nz ) B. z A s z s+ λ C A = h ((Q eα) I nλ ) λ s λ s+ C A C A
36 Example: The skate on an inclined plane in Cartesian coordinates L(q,v) = m 4 (v2 + v2 2 + v2 3 + v2 4 ) + mg sin(α)(q + q 3 ) 2 = ( (q3 q ) 2 + (q 4 q 2 ) 2 l 2) =: g(q) 2 = (q 3 q )(v 3 v ) + (q 4 q 2 )(v 4 v 2 ) =: g q (q)v = (q 3 q )(v 2 + v 4 ) (q 4 q 2 )(v + v 3 ) =: k(q,v) m =, g sin(α) =, l = 2 q = l ( ) T 2 v = l ( ) T 2 t = h =.
37 Example: The skate on an inclined plane in Cartesian coordinates (cont.) 3 x 3 stage Gauss Lobatto SPARK E n E t x 7 2 stage Gauss Lobatto SPARK 2 E n E t
38 Conclusions We have presented a new discrete Lagrange-d Alembert principle for Lagrangian systems with nonholonomic constraints which is consistent with the discrete Lagrange-d Alembert principle for forced Lagrangian systems of Kane, Marsden, Ortiz, and West Unlike Cortés principle the new discrete Lagrange-d Alembert principle is satisfied by both the exact solution and a large class of SPARK methods of arbitrarily high order SPARK methods have strong foundations Implementation of SPARK methods can be done efficiently but is not straightforward
39 Further issues Convergence proof for problems with mixed holonomic/nonholonomic constraints (done for Lobatto IIIA-B-C-C -D methods by Hyounkyun Oh, PhD thesis, U. of Iowa, 25) Backward error analysis Reduction theory Dissipation/friction/impact/contact
Discontinuous Collocation Methods for DAEs in Mechanics
Discontinuous Collocation Methods for DAEs in Mechanics Scott Small Laurent O. Jay The University of Iowa Iowa City, Iowa, USA July 11, 2011 Outline 1 Introduction of the DAEs 2 SPARK and EMPRK Methods
More informationRunge-Kutta type methods for differentialalgebraic equations in mechanics
University of Iowa Iowa Research Online Theses and Dissertations Spring 211 Runge-Kutta type methods for differentialalgebraic equations in mechanics Scott Joseph Small University of Iowa Copyright 211
More informationEN Nonlinear Control and Planning in Robotics Lecture 2: System Models January 28, 2015
EN53.678 Nonlinear Control and Planning in Robotics Lecture 2: System Models January 28, 25 Prof: Marin Kobilarov. Constraints The configuration space of a mechanical sysetm is denoted by Q and is assumed
More informationIterative solution of SPARK methods applied to DAEs
Numerical Algorithms 31: 171 191 2002 2002 Kluwer Academic Publishers Printed in the Netherlands Iterative solution of SPARK methods applied to DAEs Laurent O Jay Department of Mathematics 14 MacLean Hall
More informationarxiv: v1 [math.oc] 1 Feb 2015
HIGH ORDER VARIATIONAL INTEGRATORS IN THE OPTIMAL CONTROL OF MECHANICAL SYSTEMS CÉDRIC M. CAMPOS, SINA OBER-BLÖBAUM, AND EMMANUEL TRÉLAT arxiv:152.325v1 [math.oc] 1 Feb 215 Abstract. In recent years, much
More informationSolving high order nonholonomic systems using Gibbs-Appell method
Solving high order nonholonomic systems using Gibbs-Appell method Mohsen Emami, Hassan Zohoor and Saeed Sohrabpour Abstract. In this paper we present a new formulation, based on Gibbs- Appell method, for
More informationSPECIALIZED RUNGE-KUTTA METHODS FOR INDEX 2 DIFFERENTIAL-ALGEBRAIC EQUATIONS
MATHEMATICS OF COMPUTATION Volume 75, Number 254, Pages 641 654 S 0025-5718(05)01809-0 Article electronically published on December 19, 2005 SPECIALIZED RUNGE-KUTTA METHODS FOR INDEX 2 DIFFERENTIAL-ALGEBRAIC
More informationCDS 205 Final Project: Incorporating Nonholonomic Constraints in Basic Geometric Mechanics Concepts
CDS 205 Final Project: Incorporating Nonholonomic Constraints in Basic Geometric Mechanics Concepts Michael Wolf wolf@caltech.edu 6 June 2005 1 Introduction 1.1 Motivation While most mechanics and dynamics
More informationVariational Integrators for Electrical Circuits
Variational Integrators for Electrical Circuits Sina Ober-Blöbaum California Institute of Technology Joint work with Jerrold E. Marsden, Houman Owhadi, Molei Tao, and Mulin Cheng Structured Integrators
More informationPhysics 5153 Classical Mechanics. Canonical Transformations-1
1 Introduction Physics 5153 Classical Mechanics Canonical Transformations The choice of generalized coordinates used to describe a physical system is completely arbitrary, but the Lagrangian is invariant
More informationSTRUCTURE PRESERVATION FOR CONSTRAINED DYNAMICS WITH SUPER PARTITIONED ADDITIVE RUNGE KUTTA METHODS
SIAM J. SCI. COMPUT. Vol. 20, No. 2, pp. 416 446 c 1998 Society for Industrial and Applied Mathematics STRUCTURE PRESERVATION FOR CONSTRAINED DYNAMICS WITH SUPER PARTITIONED ADDITIVE RUNGE KUTTA METHODS
More informationReducing round-off errors in symmetric multistep methods
Reducing round-off errors in symmetric multistep methods Paola Console a, Ernst Hairer a a Section de Mathématiques, Université de Genève, 2-4 rue du Lièvre, CH-1211 Genève 4, Switzerland. (Paola.Console@unige.ch,
More informationSPECTRAL VARIATIONAL INTEGRATORS
SPECTRAL VARIATIONAL INTEGRATORS JAMES HALL AND MELVIN LEOK Abstract. In this paper, we present a new variational integrator for problems in Lagrangian mechanics. Using techniques from Galerkin variational
More informationDynamics of Nonholonomic Systems
Dynamics of Nonholonomic Systems 1 Example: When is the skate in equilibrium? x 2 Number of degrees of freedom n = 3 degrees of freedom m= 1 nonholonomic constraint Generalized coordinates (x 1, x 2, x
More informationChapter 3 Numerical Methods
Chapter 3 Numerical Methods Part 3 3.4 Differential Algebraic Systems 3.5 Integration of Differential Equations 1 Outline 3.4 Differential Algebraic Systems 3.4.1 Constrained Dynamics 3.4.2 First and Second
More informationSYMMETRIC PROJECTION METHODS FOR DIFFERENTIAL EQUATIONS ON MANIFOLDS
BIT 0006-3835/00/4004-0726 $15.00 2000, Vol. 40, No. 4, pp. 726 734 c Swets & Zeitlinger SYMMETRIC PROJECTION METHODS FOR DIFFERENTIAL EQUATIONS ON MANIFOLDS E. HAIRER Section de mathématiques, Université
More informationDirac Structures and the Legendre Transformation for Implicit Lagrangian and Hamiltonian Systems
Dirac Structures and the Legendre Transformation for Implicit Lagrangian and Hamiltonian Systems Hiroaki Yoshimura Mechanical Engineering, Waseda University Tokyo, Japan Joint Work with Jerrold E. Marsden
More informationREVIEW. Hamilton s principle. based on FW-18. Variational statement of mechanics: (for conservative forces) action Equivalent to Newton s laws!
Hamilton s principle Variational statement of mechanics: (for conservative forces) action Equivalent to Newton s laws! based on FW-18 REVIEW the particle takes the path that minimizes the integrated difference
More informationDiscrete Dirac Mechanics and Discrete Dirac Geometry
Discrete Dirac Mechanics and Discrete Dirac Geometry Melvin Leok Mathematics, University of California, San Diego Joint work with Anthony Bloch and Tomoki Ohsawa Geometric Numerical Integration Workshop,
More informationEXTENSIONS OF THE HHT-α METHOD TO DIFFERENTIAL-ALGEBRAIC EQUATIONS IN MECHANICS
EXTENSIONS OF THE HHT-α METHOD TO DIFFERENTIAL-ALGEBRAIC EQUATIONS IN MECHANICS LAURENT O. JAY AND DAN NEGRUT Abstract. We present second-order extensions of the Hilber-Hughes-Taylor-α HHT-α) method for
More informationVariational Collision Integrators and Optimal Control
Variational Collision Integrators and Optimal Control David Pekarek, and Jerrold Marsden 1 Abstract This paper presents a methodology for generating locally optimal control policies for mechanical systems
More informationPart of the advantage : Constraint forces do no virtual. work under a set of virtual displacements compatible
FORCES OF CONSTRAINT Lagrangian formalism : Generalized coordinate Minimum set of Eqns Part of the advantage : Constraint forces do no virtual work under a set of virtual displacements compatible with
More informationReview of Lagrangian Mechanics and Reduction
Review of Lagrangian Mechanics and Reduction Joel W. Burdick and Patricio Vela California Institute of Technology Mechanical Engineering, BioEngineering Pasadena, CA 91125, USA Verona Short Course, August
More informationOptimal Control Using Nonholonomic Integrators
7 IEEE International Conference on Robotics and Automation Roma, Italy, 1-1 April 7 ThA6.3 Optimal Control Using onholonomic Integrators Marin Kobilarov and Gaurav Suatme Robotic Embedded Systems Laboratory
More informationVideo 3.1 Vijay Kumar and Ani Hsieh
Video 3.1 Vijay Kumar and Ani Hsieh Robo3x-1.3 1 Dynamics of Robot Arms Vijay Kumar and Ani Hsieh University of Pennsylvania Robo3x-1.3 2 Lagrange s Equation of Motion Lagrangian Kinetic Energy Potential
More informationDiscretizations of Lagrangian Mechanics
Discretizations of Lagrangian Mechanics George W. Patrick Applied Mathematics and Mathematical Physics Mathematics and Statistics, University of Saskatchewan, Canada August 2007 With Charles Cuell Discretizations
More informationGeometric Mechanics and Global Nonlinear Control for Multi-Body Dynamics
Geometric Mechanics and Global Nonlinear Control for Multi-Body Dynamics Harris McClamroch Aerospace Engineering, University of Michigan Joint work with Taeyoung Lee (George Washington University) Melvin
More informationarxiv: v1 [physics.plasm-ph] 15 Sep 2013
Symplectic Integration of Magnetic Systems arxiv:309.373v [physics.plasm-ph] 5 Sep 03 Abstract Stephen D. Webb Tech-X Corporation, 56 Arapahoe Ave., Boulder, CO 80303 Dependable numerical results from
More informationVariational integrators for electric and nonsmooth systems
Variational integrators for electric and nonsmooth systems Sina Ober-Blöbaum Control Group Department of Engineering Science University of Oxford Fellow of Harris Manchester College Summerschool Applied
More informationVariational principles and Hamiltonian Mechanics
A Primer on Geometric Mechanics Variational principles and Hamiltonian Mechanics Alex L. Castro, PUC Rio de Janeiro Henry O. Jacobs, CMS, Caltech Christian Lessig, CMS, Caltech Alex L. Castro (PUC-Rio)
More informationChapter 1 Symplectic Integrator and Beam Dynamics Simulations
Chapter 1 and Beam Accelerator Physics Group, Journal Club November 9, 2010 and Beam NSLS-II Brookhaven National Laboratory 1.1 (70) Big Picture for Numerical Accelerator Physics and Beam For single particle
More informationStability of Asynchronous Variational Integrators
Stability of Asynchronous Variational Integrators Abstract The formulation of multiple time step integrators can provide substantial computational savings for mechanical systems with multiple time scales.
More informationGeneralized Forces. Hamilton Principle. Lagrange s Equations
Chapter 5 Virtual Work and Lagrangian Dynamics Overview: Virtual work can be used to derive the dynamic and static equations without considering the constraint forces as was done in the Newtonian Mechanics,
More informationNumerische Mathematik
Numer. Math. DOI 1.17/s211-14-679- Numerische Mathematik Spectral variational integrators James Hall Melvin Leok Received: 19 November 212 / Revised: 16 October 214 Springer-Verlag Berlin Heidelberg 214
More informationChapter VII. Non-Canonical Hamiltonian Systems
Chapter VII. Non-Canonical Hamiltonian Systems We discuss theoretical properties and the structure-preserving numerical treatment of Hamiltonian systems on manifolds and of the closely related class of
More informationMath 692A: Geometric Numerical Integration
Math 692A: Geometric Numerical Integration Melvin Leok Mathematics, Purdue University. mleok@math.purdue.edu http://www.math.purdue.edu/ mleok/ Mathematics, Purdue University NSF DMS-54747, 714223 and
More informationInvariant Lagrangian Systems on Lie Groups
Invariant Lagrangian Systems on Lie Groups Dennis Barrett Geometry and Geometric Control (GGC) Research Group Department of Mathematics (Pure and Applied) Rhodes University, Grahamstown 6140 Eastern Cape
More informationGeneralized Coordinates, Lagrangians
Generalized Coordinates, Lagrangians Sourendu Gupta TIFR, Mumbai, India Classical Mechanics 2012 August 10, 2012 Generalized coordinates Consider again the motion of a simple pendulum. Since it is one
More informationModified Equations for Variational Integrators
Modified Equations for Variational Integrators Mats Vermeeren Technische Universität Berlin Groningen December 18, 2018 Mats Vermeeren (TU Berlin) Modified equations for variational integrators December
More informationNonholonomic systems as restricted Euler-Lagrange systems
Nonholonomic systems as restricted Euler-Lagrange systems T. Mestdag and M. Crampin Abstract We recall the notion of a nonholonomic system by means of an example of classical mechanics, namely the vertical
More information06. Lagrangian Mechanics II
University of Rhode Island DigitalCommons@URI Classical Dynamics Physics Course Materials 2015 06. Lagrangian Mechanics II Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative Commons License
More informationRobot Control Basics CS 685
Robot Control Basics CS 685 Control basics Use some concepts from control theory to understand and learn how to control robots Control Theory general field studies control and understanding of behavior
More informationCHAPTER 10: Numerical Methods for DAEs
CHAPTER 10: Numerical Methods for DAEs Numerical approaches for the solution of DAEs divide roughly into two classes: 1. direct discretization 2. reformulation (index reduction) plus discretization Direct
More informationStochastic Variational Integrators
Stochastic Variational Integrators Jeremy Schmitt Variational integrators are numerical geometric integrator s derived from discretizing Hamilton s principle. They are symplectic integrators that exhibit
More informationPartitioned Runge-Kutta Methods for Semi-explicit Differential-Algebraic Systems of Index 2
Partitioned Runge-Kutta Methods for Semi-explicit Differential-Algebraic Systems of Index 2 A. Murua Konputazio Zientziak eta A. A. Saila, EHU/UPV, Donostia-San Sebastián email ander@si.ehu.es December
More informationMultibody simulation
Multibody simulation Dynamics of a multibody system (Euler-Lagrange formulation) Dimitar Dimitrov Örebro University June 16, 2012 Main points covered Euler-Lagrange formulation manipulator inertia matrix
More informationTransverse Linearization for Controlled Mechanical Systems with Several Passive Degrees of Freedom (Application to Orbital Stabilization)
Transverse Linearization for Controlled Mechanical Systems with Several Passive Degrees of Freedom (Application to Orbital Stabilization) Anton Shiriaev 1,2, Leonid Freidovich 1, Sergey Gusev 3 1 Department
More information1 Hamiltonian formalism
1 Hamiltonian formalism 1.1 Hamilton s principle of stationary action A dynamical system with a finite number n degrees of freedom can be described by real functions of time q i (t) (i =1, 2,..., n) which,
More information4.1 Important Notes on Notation
Chapter 4. Lagrangian Dynamics (Most of the material presented in this chapter is taken from Thornton and Marion, Chap. 7) 4.1 Important Notes on Notation In this chapter, unless otherwise stated, the
More informationIn most robotic applications the goal is to find a multi-body dynamics description formulated
Chapter 3 Dynamics Mathematical models of a robot s dynamics provide a description of why things move when forces are generated in and applied on the system. They play an important role for both simulation
More informationClassical Mechanics and Electrodynamics
Classical Mechanics and Electrodynamics Lecture notes FYS 3120 Jon Magne Leinaas Department of Physics, University of Oslo December 2009 2 Preface These notes are prepared for the physics course FYS 3120,
More informationModelling of dynamics of mechanical systems with regard for constraint stabilization
IOP Conference Series: Materials Science and Engineering PAPER OPEN ACCESS Modelling of dnamics of mechanical sstems with regard for constraint stabilization o cite this article: R G Muharlamov 018 IOP
More informationarxiv:math-ph/ v1 22 May 2003
On the global version of Euler-Lagrange equations. arxiv:math-ph/0305045v1 22 May 2003 R. E. Gamboa Saraví and J. E. Solomin Departamento de Física, Facultad de Ciencias Exactas, Universidad Nacional de
More informationAppell-Hamel dynamical system: a nonlinear test of the Chetaev and the vakonomic model
ZAMM Z. Angew. Math. Mech. 87, No. 10, 692 697 (2007) / DOI 10.1002/zamm.200710344 Appell-Hamel dynamical system: a nonlinear test of the Chetaev and the vakonomic model Shan-Shan Xu 1, Shu-Min Li 1,2,
More informationClassical Mechanics and Electrodynamics
Classical Mechanics and Electrodynamics Lecture notes FYS 3120 Jon Magne Leinaas Department of Physics, University of Oslo 2 Preface FYS 3120 is a course in classical theoretical physics, which covers
More informationDiscrete Dirac Structures and Implicit Discrete Lagrangian and Hamiltonian Systems
Discrete Dirac Structures and Implicit Discrete Lagrangian and Hamiltonian Systems Melvin Leok and Tomoki Ohsawa Department of Mathematics, University of California, San Diego, La Jolla, CA 92093-0112,
More informationConstrained motion and generalized coordinates
Constrained motion and generalized coordinates based on FW-13 Often, the motion of particles is restricted by constraints, and we want to: work only with independent degrees of freedom (coordinates) k
More informationRobotics. Dynamics. Marc Toussaint U Stuttgart
Robotics Dynamics 1D point mass, damping & oscillation, PID, dynamics of mechanical systems, Euler-Lagrange equation, Newton-Euler recursion, general robot dynamics, joint space control, reference trajectory
More informationS O A O lllllu O, W E D N E S D A Y N E X T ; D E C 2 9 t h. R E D U C T I O N S. 3 1, 3 2, 3 3, R o b e r t s o n S t., 2 & 3 Q a r e m o n t,
2» q» Z - - Qk»Q») v Q [ z -? ( 27 2 - ( q - < 5 75 7 v v 50 Q > vk 3 Q - -v- 2 3 7> 27
More informationErnst Hairer 1, Robert I. McLachlan 2 and Robert D. Skeel 3
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher ON ENERGY CONSERVATION OF THE SIMPLIFIED TAKAHASHI IMADA METHOD, Ernst Hairer 1,
More informationVariational Integrators for Interconnected Lagrange Dirac Systems
DOI 10.1007/s00332-017-9364-7 Variational Integrators for Interconnected Lagrange Dirac Systems Helen Parks 1 Melvin Leok 1 Received: 4 March 2016 / Accepted: 28 January 2017 Springer Science+Business
More informationGeneralized coordinates and constraints
Generalized coordinates and constraints Basilio Bona DAUIN Politecnico di Torino Semester 1, 2014-15 B. Bona (DAUIN) Generalized coordinates and constraints Semester 1, 2014-15 1 / 25 Coordinates A rigid
More informationarxiv: v1 [math-ph] 12 Apr 2017
Dirac structures in nonequilibrium thermodynamics arxiv:1704.03935v1 [math-ph] 12 Apr 2017 François Gay-Balmaz Hiroaki Yoshimura CNRS, LMD, IPSL School of Science and Engineering Ecole Normale Supérieure
More informationLAGRANGIAN AND HAMILTONIAN
LAGRANGIAN AND HAMILTONIAN A. Constraints and Degrees of Freedom. A constraint is a restriction on the freedom of motion of a system of particles in the form of a condition. The number of independent ways
More informationarxiv: v1 [math.ds] 18 Nov 2008
arxiv:0811.2889v1 [math.ds] 18 Nov 2008 Abstract Quaternions And Dynamics Basile Graf basile.graf@epfl.ch February, 2007 We give a simple and self contained introduction to quaternions and their practical
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 informationGeometric Numerical Integration
Geometric Numerical Integration (Ernst Hairer, TU München, winter 2009/10) Development of numerical ordinary differential equations Nonstiff differential equations (since about 1850), see [4, 2, 1] Adams
More informationDirac Structures in Lagrangian Mechanics
Dirac Structures in Lagrangian Mechanics Part II: Variational Structures Hiroaki Yoshimura Jerrold E. Marsden Department of Mechanical Engineering Control and Dynamical Systems Waseda University California
More informationThe Principle of Least Action
The Principle of Least Action Anders Svensson Abstract In this paper, the principle of least action in classical mechanics is studied. The term is used in several different contexts, mainly for Hamilton
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 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 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 informationThe Bock iteration for the ODE estimation problem
he Bock iteration for the ODE estimation problem M.R.Osborne Contents 1 Introduction 2 2 Introducing the Bock iteration 5 3 he ODE estimation problem 7 4 he Bock iteration for the smoothing problem 12
More informationThe Geometry of Euler s equation. Introduction
The Geometry of Euler s equation Introduction Part 1 Mechanical systems with constraints, symmetries flexible joint fixed length In principle can be dealt with by applying F=ma, but this can become complicated
More informationPhysics 106a, Caltech 16 October, Lecture 5: Hamilton s Principle with Constraints. Examples
Physics 106a, Caltech 16 October, 2018 Lecture 5: Hamilton s Principle with Constraints We have been avoiding forces of constraint, because in many cases they are uninteresting, and the constraints can
More informationLecture IV: Time Discretization
Lecture IV: Time Discretization Motivation Kinematics: continuous motion in continuous time Computer simulation: Discrete time steps t Discrete Space (mesh particles) Updating Position Force induces acceleration.
More informationNonholonomic Constraints Examples
Nonholonomic Constraints Examples Basilio Bona DAUIN Politecnico di Torino July 2009 B. Bona (DAUIN) Examples July 2009 1 / 34 Example 1 Given q T = [ x y ] T check that the constraint φ(q) = (2x + siny
More informationRobotics. Dynamics. University of Stuttgart Winter 2018/19
Robotics Dynamics 1D point mass, damping & oscillation, PID, dynamics of mechanical systems, Euler-Lagrange equation, Newton-Euler, joint space control, reference trajectory following, optimal operational
More informationRunge-Kutta Theory and Constraint Programming Julien Alexandre dit Sandretto Alexandre Chapoutot. Department U2IS ENSTA ParisTech SCAN Uppsala
Runge-Kutta Theory and Constraint Programming Julien Alexandre dit Sandretto Alexandre Chapoutot Department U2IS ENSTA ParisTech SCAN 2016 - Uppsala Contents Numerical integration Runge-Kutta with interval
More informationMarion and Thornton. Tyler Shendruk October 1, Hamilton s Principle - Lagrangian and Hamiltonian dynamics.
Marion and Thornton Tyler Shendruk October 1, 2010 1 Marion and Thornton Chapter 7 Hamilton s Principle - Lagrangian and Hamiltonian dynamics. 1.1 Problem 6.4 s r z θ Figure 1: Geodesic on circular cylinder
More informationCase Study: The Pelican Prototype Robot
5 Case Study: The Pelican Prototype Robot The purpose of this chapter is twofold: first, to present in detail the model of the experimental robot arm of the Robotics lab. from the CICESE Research Center,
More informationOn variational integrators for optimal control of mechanical control systems. Leonardo Colombo, David Martín de Diego and Marcela Zuccalli
RACSAM Rev. R. Acad. Cien. Serie A. Mat. VOL. Falta Falta), Falta, pp. 9 10 Falta Comunicación Preliminar / Preliminary Communication On variational integrators for optimal control of mechanical control
More informationIntroduction to the. Geometry of Classical Dynamics
Introduction to the Geometry of Classical Dynamics Renato Grassini Dipartimento di Matematica e Applicazioni Università di Napoli Federico II HIKARI LT D HIKARI LTD Hikari Ltd is a publisher of international
More informationLecture Note 12: Dynamics of Open Chains: Lagrangian Formulation
ECE5463: Introduction to Robotics Lecture Note 12: Dynamics of Open Chains: Lagrangian Formulation Prof. Wei Zhang Department of Electrical and Computer Engineering Ohio State University Columbus, Ohio,
More informationRobotics 2 Robot Interaction with the Environment
Robotics 2 Robot Interaction with the Environment Prof. Alessandro De Luca Robot-environment interaction a robot (end-effector) may interact with the environment! modifying the state of the environment
More informationAn Overview of Variational Integrators
An Overview of Variational Integrators A. Lew, M. West Stanford University J. Marsden, M. Ortiz Caltech April, 2004 Variational Integrators What? Theory to generate space and time integrators Why? Conservation
More informationEULER-LAGRANGE TO HAMILTON. The goal of these notes is to give one way of getting from the Euler-Lagrange equations to Hamilton s equations.
EULER-LAGRANGE TO HAMILTON LANCE D. DRAGER The goal of these notes is to give one way of getting from the Euler-Lagrange equations to Hamilton s equations. 1. Euler-Lagrange to Hamilton We will often write
More informationHamilton-Pontryagin Integrators on Lie Groups
Hamilton-Pontryagin Integrators on Lie Groups Thesis by Nawaf Bou-Rabee In Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy California Institute of Technology Pasadena, California
More informationPhysics 106a, Caltech 4 December, Lecture 18: Examples on Rigid Body Dynamics. Rotating rectangle. Heavy symmetric top
Physics 106a, Caltech 4 December, 2018 Lecture 18: Examples on Rigid Body Dynamics I go through a number of examples illustrating the methods of solving rigid body dynamics. In most cases, the problem
More informationKinematic and Dynamic Models
Kinematic and Dynamic Models b Advanced Robotics: Navigation and Perception 3/29/2011 1 Reasoning about mobility is important. Image: RedTeamRacing.com 2 ecture Outline Nomenclature Kinematic Constraints
More informationLinear Algebra and Robot Modeling
Linear Algebra and Robot Modeling Nathan Ratliff Abstract Linear algebra is fundamental to robot modeling, control, and optimization. This document reviews some of the basic kinematic equations and uses
More informationLecture 10 Reprise and generalized forces
Lecture 10 Reprise and generalized forces The Lagrangian Holonomic constraints Generalized coordinates Nonholonomic constraints Generalized forces we haven t done this, so let s start with it Euler Lagrange
More informationMSMS Matlab Problem 02
MSMS 2014-2015 Matlab Problem 02 Basilio Bona DAUIN PoliTo Problem formulation The planar system illustrated in Figure 1 consists of a cart C sliding with friction along the horizontal rail; the cart supports
More informationROBOTICS Laboratory Problem 02
ROBOTICS 2015-2016 Laboratory Problem 02 Basilio Bona DAUIN PoliTo Problem formulation The planar system illustrated in Figure 1 consists of a cart C sliding with or without friction along the horizontal
More informationD Alembert s principle of virtual work
PH101 Lecture 9 Review of Lagrange s equations from D Alembert s Principle, Examples of Generalized Forces a way to deal with friction, and other non-conservative forces D Alembert s principle of virtual
More informationDynamics. Basilio Bona. Semester 1, DAUIN Politecnico di Torino. B. Bona (DAUIN) Dynamics Semester 1, / 18
Dynamics Basilio Bona DAUIN Politecnico di Torino Semester 1, 2016-17 B. Bona (DAUIN) Dynamics Semester 1, 2016-17 1 / 18 Dynamics Dynamics studies the relations between the 3D space generalized forces
More informationNoether s Theorem. 4.1 Ignorable Coordinates
4 Noether s Theorem 4.1 Ignorable Coordinates A central recurring theme in mathematical physics is the connection between symmetries and conservation laws, in particular the connection between the symmetries
More informationIntroduction to SAMCEF MECANO
Introduction to SAMCEF MECANO 1 Outline Introduction Generalized coordinates Kinematic constraints Time integration Description and paramerization of finite rotations Acknowledgements Michel Géradin (ULg)
More informationUNIVERSITY OF CALIFORNIA, SAN DIEGO. Convergence of Galerkin Variational Integrators for Vector Spaces and Lie Groups
UNIVERSITY OF CALIFORNIA, SAN DIEGO Convergence of Galerkin Variational Integrators for Vector Spaces and Lie Groups A dissertation submitted in partial satisfaction of the requirements for the degree
More informationMassachusetts Institute of Technology Department of Physics. Final Examination December 17, 2004
Massachusetts Institute of Technology Department of Physics Course: 8.09 Classical Mechanics Term: Fall 004 Final Examination December 17, 004 Instructions Do not start until you are told to do so. Solve
More information