Numerical Methods for Rigid Multibody Dynamics
|
|
- Homer White
- 5 years ago
- Views:
Transcription
1 Numerical Methods for Rigid Multibody Dynamics Claus Führer Centre for Mathematical Sciences Lund University Lappenranta 2012
2 Unit 0: Preface These notes serve as a skeleton for the compact course. They document together with the assignments the course outline and course content. All references in the notes refer to the textbook by Eich-Soellner/Führer if not otherwise stated. The notes are a guide to read the textbook. They are no textbook. 1
3 Unit 1: Equations of Motion Unconstrained system: 2nd order ODEs Constrained system: 2nd order DAEs Linear Implicitness: O(n)-methods. 2
4 1.1 Multibody System A multibody system (MBS) consist of bodies: mass, inertia, rigid or elastic degrees of freedom (dofs) interconnections: force elements, e.g. springs, dampers, controlers joints: constraints which reduce the dofs 3
5 1.2 The MBS Simulation Loop mechanical system engineer modelling, abstraction multibody model formalism establish equations of motion M p = f(p; _p) G(p) T equations of motion 0 = g(p) numerics solution of equations of motion 25.9d d d d07 85.d0d d d d d d07 85.d0d d d d d d07 85.d0d d+03 simulation results 4
6 1.3 Basic Mathematical Tasks in Multibody Dynamics Kinematic analysis Given a robot and its hand position. How must the joint angles be chosen? Nonlinear equation system. Static equilibrium position How must the forces be chosen such that the system is in equilibrium? Nonlinear equation system. Dynamic simulation How does the system behave dynamically? Numerically Solving ODEs/DAEs. Linear system analysis Stability? Input/output behavior? Linearization, Eigenanalysis. 5
7 Design, optimal control How can the behavior of the system be optimized? Parameter Identification, System Optimization, Optimal Control. 6
8 1.4 Equations of Motion: Unconstrained System or in first order format M p = f a (t, p, ṗ). ṗ = v M v = f a (t, p, v) p positions, v velocities t time M n p n p mass matrix f a applied forces 7
9 1.5 Model of Unconstrained Truck p p p p CM5 6 CM4 p 7 loading area (body 5) cabin (body 4) chassis (body 2) p 3 p 2 CM2 p rear wheel (body 1) p 1 4 CM1 front wheel (body 3) CM3 inertial frame (body 0) 8
10 1.6 Coordinates p p p p CM5 6 CM4 p 7 loading area (body 5) cabin (body 4) chassis (body 2) p 2 p 3 CM2 p 1 rear wheel (body 1) CM1 front wheel (body 3) p 4 inertial frame (body 0) CM3 p 1 Vert motion of the rear wheel (body 1) p 2 Vert motion of the truck chassis (body 2) p 3 Rot of the truck chassis (body 2) p 4 Vert motion of the front wheel (body 3) p 5 Vert motion of the driver cabin (body 4) p 6 Rot of the driver cabin (body 4) p 7 Horiz motion of the loading area (body 5) p 8 Vert motion of the loading area (body 5) p 9 Rot of the loading area (body 5) 9
11 1.7 Equations of motion: Unconstr. Truck m 1 p 1 = f f 122 m 1 g gr m 2 p 2 = f 122 f f f f f 1d2 + f 2d2 m 2 g gr l 2 p 3 = ( a 23 f 232 a 12 f 122 h 1 (f f 121 )) cos p 3 ( a 23 f 231 a 12 f 121 h 1 (f f 122 )) sin p 3 + (a 25 f a 52 (f 1d2 + f 2d2 ) + h 2 (f f 1d1 + f 2d1 )) cos p 3 (a 25 f a 52 (f 1d1 + f 2d1 ) + h 2 (f f 1d2 + f 2d2 )) sin p 3 (a 24 f a 42 f h 2 (f f 421 )) cos p 3 (a 24 f a 42 f h 2 (f f 422 )) sin p 3 m 3 p 4 =.. f f 232 m 3 g gr 10
12 1.8 Unconstr. Truck: Geometry and Forces c d1 c 25 b 24 b 42 c d2 h 3 f 1d f 2d f 25 f 24 f 42 f 1d f 2d f 12 f 25 f 24 f 42 f 23 h 2 a 52 a 24 a 12 a 25 a 23 a 42 h 1 f 12 f 23 f 10 f 30 11
13 1.9 Constrained Multibody Systems M p = f a (t, p, ṗ) G T (p)λ 0 = g(p) with g(p) being n c constraints G(p) := d dp g(p) constraint Jacobian (n c n p ) λ Lagrange multipliers 12
14 1.10 Constrained Truck p 8 p 9 p 5 p 6 p 7 loading area (body 5) cabin (body 4) chassis (body 2) p 3 p 2 p 1 rear wheel (body 1) p 4 front wheel (body 3) reference frame (body 0) 13
15 1.11 Constrained Truck: Typical Constraint p p p p 6 p 7 loading area (body 5) cabin (body 4) p 2 chassis (body 2) p 3 p rear wheel (body 1) p 1 4 front wheel (body 3) reference frame (body 0) ρ 52 = ( ) ( ) [( p7 cc1 0 +S(p p 9 ) 8 c c2 p2) + S(p 3 ) ( ac1 a c2 )] = 0. 14
16 1.12 Different Types of constraints holonomic constraints non-holonomic constraints rhenomic constraints g(p, t) skleronomic constraints g(p). 15
17 1.13 3D-Rotation: Kardan Angles S(α, β, γ) = S(0, 0, γ)s(0, β, 0)S(α, 0, 0) with the elementary rotation matrices S, e.g. S(α, 0, 0) = cos α sin α. 0 sin α cos α α is the angle of the rotation about the x-axis, β the angle of rotation about the (new) y-axis and γ the angle about the (even newer) z-axis. 16
18 1.14 3D-Rotation: angular velocity Poisson equation: with ω: Ṡ(α, β, γ) = ωs(α, β, γ). ω := 0 ω 3 ω 2 ω 3 0 ω 1 ω 2 ω 1 0 This defines the angular velocities ω 1, ω 2, ω 3. 17
19 1.15 3D-Rotation: angular velocity vs ṗ α β = γ cos γ/ cos β sin γ/ cos β 0 sin γ cos γ 0 ω 1 ω 2. cos γ tan β sin γ tan β 1 ω 3 Singularities! 18
20 1.16 3D-Rotation: complete system ṗ = Z(p)v M v = f a (t, p, v, s) Z(p) T G(p) T λ 0 = g(p) (see also Exercise 1). 19
21 1.17 Relative vs Absolute coordinates Absolute coordinates Mass matrix constant and block-diagonal. Constraints DAEs. Relative coordinates: Mass matrix dense and position dependent. For tree structured systems the constraints can be directly eliminated. 20
22 1.18 Tree structured systems
23 1.19 Truck: relative coordinates p 9 p 8 p 6 p 5 p 7 loading area (body 5) cabin (body 4) chassis (body 2) p 2 p 3 p 1 rear wheel (body 1) p 4 front wheel (body 3) q 1 = ρ 10 q 2 = ρ 12 q 3 = p 3 q 4 = ρ 23 q 5 = ρ 24 q 6 = p 6 p 3 q 7 = p 9 p 3. 7 degrees of freedom (dofs) reference frame (body 0) 22
24 1.20 Relation between absolute/relative positions Absolute relative: n q = n dof equations 0 = g(p, q), Relative absolute: n p = n c + n q equations 0 = ( ) g(p) g(p, q) =: γ(p, q). 23
25 1.21 Relation between absolute/relative velocities Derivative assumed to be regular (Grübler condition): Γ p (p, q) = γ(p, q) p Velocities: and 0 = γ(p, q) } p {{} =:Γ p (p,q) ṗ + γ(p, q) } q {{} =:Γ q (p,q) ṗ = Γ p (p, q) 1 Γ q (p, q) }{{} =:V (p,q) q, q, 24
26 1.22 Relation between absolute/relative accelerations p = V (p, q) q + ζ(p, q, ṗ, q). with ζ collecting all terms with lower derivatives of p and q. 25
27 1.23 State Space Form in relative coordinates We note G(p)V (p, q) = 0. From this we obtain V (p, q) T MV (p, q) }{{} =: M(p,q) q = V (p, q) T (f a (p, ṗ) Mζ(p, q, ṗ, q)). Note: Mass matrix looses its structure and becomes state dependent. 26
28 1.24 Mixed Coordinate Formulation for Tree Structured Systems For chain like strucktures, mixed coordinate formulations are much more efficient. block diagonal matrices to invert. k k 1 i
29 1.25 Mixed Coordinate Formulation for Tree Structured Systems ( ) M 0 0 0) ( p q = ( ) fa (t, p, ṗ) 0 0 = γ(p, q). ( ) Γp (p, q) T Γ q (p, q) T µ We note the i-th interconnection has the form γ i (p, q) = γ i (p i, p i 1, q i ). Resorting of components gives... 28
30 1.26 Mixed Coordinate Formulation for... (Cont.) Resorting:x = (..., ( p i ) T, ( q i ) T, (µ i ) T, ( p i 1 ) T, ( q i 1 ) T, (µ i 1 ) T,...) T A k C T k C k A k 1 Ck 1 T C k 1 A k 2 Ck 2 T C k 2... C T 3 C 3 A 2 C T 2 C 2 A 1 C T 1 C 1 A 0 with x i = (( p i ) T, ( q i ) T, (µ i ) T ) T and b i = ((f i a) T, 0, (z i ) T ) T. x k x k 1 x k 2.. x 2 x 1 x 0 = b k b k 1 b k 2.. b 2 b 1 b 0 29
31 1.27 Mixed Coordinate Formulation for... (Cont.) A i := M i 0 Γ T p i,i 0 0 Γ T q i,i Γ pi,i Γ qi,i 0. (1) A i maximal dimension 18 = 3 6. C T i := Γ pi,i
32 1.28 Mixed Coordinate Formulation for... (Cont.) a.) Downward recursion Initialization: Recursion (i=k-1:-1:0)  k := A k ˆbk := b k. b.) Upward recursion: Solve for x 0 and then for x i :  i := A i C i+1  1 ˆbi := b i C i  1 i i+1 CT i+1 ˆbi+1.  0 x 0 = ˆb 0  i x i = ˆb i Ci T xi 1 with i = 1 : k. 31
33 1.29 Linearization Some assumptions f a (t, p, v) = f a (p, v, u(t)) u(t) plays in the sequel the role of a given input function Time dependency of the constraint = kinematic excitation g(t, p) = g(p) z(t). 32
34 1.30 Linearization: Nominal Solution Nominal Solutions and p N (t), v N (t) and λ N (t) Nominal input u N (t) = 0 and z N (t) = 0 33
35 1.31 Linearization: Taylor ṗ = v M(t) v = K(t) p D(t) v G(t) λ + B(t)u(t) 0 = G(t) p z(t) Here: Numerical Linearization, see Lab 1 34
36 1.32 Linearization: Matrices mass matrix M(t) := M(p N (t)) stiffness matrix K(t) := ( v T N p M(p) p f a(p, v N (t), 0) + p G(p)T λ N ) p=pn (t) damping matrix D(t) := ( v f a(p N (t), v, 0)) v=vn (t) constraint matrix G(t) := ( d dp g(p)) p=p N (t) input matrix B(t) := ( u f a(p N (t), v N (t), u)) u=0 35
Chapter 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 informationRobotics & Automation. Lecture 25. Dynamics of Constrained Systems, Dynamic Control. John T. Wen. April 26, 2007
Robotics & Automation Lecture 25 Dynamics of Constrained Systems, Dynamic Control John T. Wen April 26, 2007 Last Time Order N Forward Dynamics (3-sweep algorithm) Factorization perspective: causal-anticausal
More informationStructural Dynamics Lecture 4. Outline of Lecture 4. Multi-Degree-of-Freedom Systems. Formulation of Equations of Motions. Undamped Eigenvibrations.
Outline of Multi-Degree-of-Freedom Systems Formulation of Equations of Motions. Newton s 2 nd Law Applied to Free Masses. D Alembert s Principle. Basic Equations of Motion for Forced Vibrations of Linear
More information4.1 Introduction Issues of applied dynamics CHAPTER 4. DYNAMICS 191
Chapter 4 Dynamics Dynamics is the branch of mechanics that is concerned with the study of motion and the relation between the forces and motion. The central focus of our study is the dynamics of systems
More informationGENERAL FORMULATION OF AN EFFICIENT RECURSIVE ALGORITHM BASED ON CANONICAL MOMENTA FOR FORWARD DYNAMICS OF CLOSED-LOOP MULTIBODY SYSTEMS
Proceedings of the XI DINAME, 28th February-4th March, 2005 - Ouro Preto - MG - Brazil Edited by D.A. Rade and V. Steffen Jr GENERAL FORMULATION OF AN EFFICIENT RECURSIVE ALGORITHM BASED ON CANONICAL MOMENTA
More informationStructural Dynamics Lecture 2. Outline of Lecture 2. Single-Degree-of-Freedom Systems (cont.)
Outline of Single-Degree-of-Freedom Systems (cont.) Linear Viscous Damped Eigenvibrations. Logarithmic decrement. Response to Harmonic and Periodic Loads. 1 Single-Degreee-of-Freedom Systems (cont.). Linear
More informationMultibody simulation
Multibody simulation Dynamics of a multibody system (Newton-Euler formulation) Dimitar Dimitrov Örebro University June 8, 2012 Main points covered Newton-Euler formulation forward dynamics inverse dynamics
More informationANALYSIS OF LOAD PATTERNS IN RUBBER COMPONENTS FOR VEHICLES
ANALYSIS OF LOAD PATTERNS IN RUBBER COMPONENTS FOR VEHICLES Jerome Merel formerly at Hutchinson Corporate Research Center Israël Wander Apex Technologies Pierangelo Masarati, Marco Morandini Dipartimento
More informationDr.Vinod Hosur, Professor, Civil Engg.Dept., Gogte Institute of Technology, Belgaum
STRUCTURAL DYNAMICS Dr.Vinod Hosur, Professor, Civil Engg.Dept., Gogte Institute of Technology, Belgaum Overview of Structural Dynamics Structure Members, joints, strength, stiffness, ductility Structure
More informationArtificial Intelligence & Neuro Cognitive Systems Fakultät für Informatik. Robot Dynamics. Dr.-Ing. John Nassour J.
Artificial Intelligence & Neuro Cognitive Systems Fakultät für Informatik Robot Dynamics Dr.-Ing. John Nassour 25.1.218 J.Nassour 1 Introduction Dynamics concerns the motion of bodies Includes Kinematics
More informationDYNAMICS OF PARALLEL MANIPULATOR
DYNAMICS OF PARALLEL MANIPULATOR The 6nx6n matrices of manipulator mass M and manipulator angular velocity W are introduced below: M = diag M 1, M 2,, M n W = diag (W 1, W 2,, W n ) From this definitions
More informationA MULTI-BODY ALGORITHM FOR WAVE ENERGY CONVERTERS EMPLOYING NONLINEAR JOINT REPRESENTATION
Proceedings of the ASME 2014 33rd International Conference on Ocean, Offshore and Arctic Engineering OMAE2014 June 8-13, 2014, San Francisco, California, USA OMAE2014-23864 A MULTI-BODY ALGORITHM FOR WAVE
More information5.5 Exercises for This Chapter Two-Axle Vehicle on Cosine Track Two-Axle Vehicle on Generally Periodic Track...
Contents 1 Introduction... 1 1.1 The Basic Function of the Wheel/rail System.... 1 1.2 Significance of Dynamics on the Operation of Rail Vehicles... 2 1.3 On the History of Research in the Field of Railway
More informationTime Response of Dynamic Systems! Multi-Dimensional Trajectories Position, velocity, and acceleration are vectors
Time Response of Dynamic Systems Robert Stengel Robotics and Intelligent Systems MAE 345, Princeton University, 217 Multi-dimensional trajectories Numerical integration Linear and nonlinear systems Linearization
More information41514 Dynamics of Machinery
41514 Dynamics of Machinery Theory, Experiment, Phenomenology and Industrial Applications Ilmar Ferreira Santos 1. Recapitulation Mathematical Modeling & Steps 2. Example System of Particle 3. Example
More informationTechnische Universität Berlin
Technische Universität Berlin Institut für Mathematik M7 - A Skateboard(v1.) Andreas Steinbrecher Preprint 216/8 Preprint-Reihe des Instituts für Mathematik Technische Universität Berlin http://www.math.tu-berlin.de/preprints
More informationLecture «Robot Dynamics»: Dynamics and Control
Lecture «Robot Dynamics»: Dynamics and Control 151-0851-00 V lecture: CAB G11 Tuesday 10:15 12:00, every week exercise: HG E1.2 Wednesday 8:15 10:00, according to schedule (about every 2nd week) Marco
More informationIntroduction to structural dynamics
Introduction to structural dynamics p n m n u n p n-1 p 3... m n-1 m 3... u n-1 u 3 k 1 c 1 u 1 u 2 k 2 m p 1 1 c 2 m2 p 2 k n c n m n u n p n m 2 p 2 u 2 m 1 p 1 u 1 Static vs dynamic analysis Static
More informationLecture Notes Multibody Dynamics B, wb1413
Lecture Notes Multibody Dynamics B, wb1413 A. L. Schwab & Guido M.J. Delhaes Laboratory for Engineering Mechanics Mechanical Engineering Delft University of Technolgy The Netherlands June 9, 29 Contents
More informationApproach based on Cartesian coordinates
GraSMech course 2005-2006 Computer-aided analysis of rigid and flexible multibody systems Approach based on Cartesian coordinates Prof. O. Verlinden Faculté polytechnique de Mons Olivier.Verlinden@fpms.ac.be
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Rigid body physics Particle system Most simple instance of a physics system Each object (body) is a particle Each particle
More information, respectively to the inverse and the inverse differential problem. Check the correctness of the obtained results. Exercise 2 y P 2 P 1.
Robotics I July 8 Exercise Define the orientation of a rigid body in the 3D space through three rotations by the angles α β and γ around three fixed axes in the sequence Y X and Z and determine the associated
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 informationRigid Body Kinetics :: Virtual Work
Rigid Body Kinetics :: Virtual Work Work-energy relation for an infinitesimal displacement: du = dt + dv (du :: total work done by all active forces) For interconnected systems, differential change in
More informationNumerical Integration of Equations of Motion
GraSMech course 2009-2010 Computer-aided analysis of rigid and flexible multibody systems Numerical Integration of Equations of Motion Prof. Olivier Verlinden (FPMs) Olivier.Verlinden@fpms.ac.be Prof.
More informationSome history. F p. 1/??
Some history F 12 10 18 p. 1/?? F 12 10 18 p. 1/?? Some history 1600: Galileo Galilei 1564 1642 cf. section 7.0 Johannes Kepler 1571 1630 cf. section 3.7 1700: Isaac Newton 1643 1727 cf. section 1.1 1750
More informationVibration Dynamics and Control
Giancarlo Genta Vibration Dynamics and Control Spri ringer Contents Series Preface Preface Symbols vii ix xxi Introduction 1 I Dynamics of Linear, Time Invariant, Systems 23 1 Conservative Discrete Vibrating
More informationPlanar Multi-body Dynamics of a Tracked Vehicle using Imaginary Wheel Model for Tracks
Defence Science Journal, Vol. 67, No. 4, July 2017, pp. 460-464, DOI : 10.14429/dsj.67.11548 2017, DESIDOC Planar Multi-body Dynamics of a Tracked Vehicle using Imaginary Wheel Model for Tracks Ilango
More information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
More informationTrajectory-tracking control of a planar 3-RRR parallel manipulator
Trajectory-tracking control of a planar 3-RRR parallel manipulator Chaman Nasa and Sandipan Bandyopadhyay Department of Engineering Design Indian Institute of Technology Madras Chennai, India Abstract
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 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 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 informationAN ALGORITHM FOR TOPOLOGY OPTIMIZATION
AN ALGORITHM FOR TOPOLOGY OPTIMIZATION OF MULTIBODY SYSTEMS TOHEED GHANDRIZ Master s thesis 2014:E52 Faculty of Science Centre for Mathematical Sciences Numerical Analysis CENTRUM SCIENTIARUM MATHEMATICARUM
More informationq 1 F m d p q 2 Figure 1: An automated crane with the relevant kinematic and dynamic definitions.
Robotics II March 7, 018 Exercise 1 An automated crane can be seen as a mechanical system with two degrees of freedom that moves along a horizontal rail subject to the actuation force F, and that transports
More informationChapter 4 Statics and dynamics of rigid bodies
Chapter 4 Statics and dynamics of rigid bodies Bachelor Program in AUTOMATION ENGINEERING Prof. Rong-yong Zhao (zhaorongyong@tongji.edu.cn) First Semester,2014-2015 Content of chapter 4 4.1 Static equilibrium
More information2.003 Engineering Dynamics Problem Set 4 (Solutions)
.003 Engineering Dynamics Problem Set 4 (Solutions) Problem 1: 1. Determine the velocity of point A on the outer rim of the spool at the instant shown when the cable is pulled to the right with a velocity
More informationDynamics of Multibody Systems: Conventional and Graph-Theoretic Approaches
Dynamics of Multibody Systems: Conventional and Graph-Theoretic Approaches by Professor John J. McPhee, P.Eng. Systems Design Engineering (Cross-appointed to Mechanical Engineering) University of Waterloo,
More informationIntroduction to Robotics
J. Zhang, L. Einig 277 / 307 MIN Faculty Department of Informatics Lecture 8 Jianwei Zhang, Lasse Einig [zhang, einig]@informatik.uni-hamburg.de University of Hamburg Faculty of Mathematics, Informatics
More informationAn Implicit Runge Kutta Solver adapted to Flexible Multibody System Simulation
An Implicit Runge Kutta Solver adapted to Flexible Multibody System Simulation Johannes Gerstmayr 7. WORKSHOP ÜBER DESKRIPTORSYSTEME 15. - 18. March 2005, Liborianum, Paderborn, Germany Austrian Academy
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 34, pp. 14-19, 2008. Copyrigt 2008,. ISSN 1068-9613. ETNA A NOTE ON NUMERICALLY CONSISTENT INITIAL VALUES FOR HIGH INDEX DIFFERENTIAL-ALGEBRAIC EQUATIONS
More informationEE Homework 3 Due Date: 03 / 30 / Spring 2015
EE 476 - Homework 3 Due Date: 03 / 30 / 2015 Spring 2015 Exercise 1 (10 points). Consider the problem of two pulleys and a mass discussed in class. We solved a version of the problem where the mass was
More informationCodal Provisions IS 1893 (Part 1) 2002
Abstract Codal Provisions IS 1893 (Part 1) 00 Paresh V. Patel Assistant Professor, Civil Engineering Department, Nirma Institute of Technology, Ahmedabad 38481 In this article codal provisions of IS 1893
More informationTranslational Mechanical Systems
Translational Mechanical Systems Basic (Idealized) Modeling Elements Interconnection Relationships -Physical Laws Derive Equation of Motion (EOM) - SDOF Energy Transfer Series and Parallel Connections
More informationReduction in number of dofs
Reduction in number of dofs Reduction in the number of dof to represent a structure reduces the size of matrices and, hence, computational cost. Because a subset of the original dof represent the whole
More informationAdvanced Dynamics. - Lecture 4 Lagrange Equations. Paolo Tiso Spring Semester 2017 ETH Zürich
Advanced Dynamics - Lecture 4 Lagrange Equations Paolo Tiso Spring Semester 2017 ETH Zürich LECTURE OBJECTIVES 1. Derive the Lagrange equations of a system of particles; 2. Show that the equation of motion
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 informationReduced-order Forward Dynamics of Multi-Closed-Loop Systems
Noname manuscript No. (will be inserted by the editor) Reduced-order Forward Dynamics of Multi-Closed-Loop Systems Majid Koul Suril V Shah S K Saha M Manivannan the date of receipt and acceptance should
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 informationAA 242B/ ME 242B: Mechanical Vibrations (Spring 2016)
AA 242B/ ME 242B: Mechanical Vibrations (Spring 2016) Homework #2 Due April 17, 2016 This homework focuses on developing a simplified analytical model of the longitudinal dynamics of an aircraft during
More informationECEN 420 LINEAR CONTROL SYSTEMS. Lecture 6 Mathematical Representation of Physical Systems II 1/67
1/67 ECEN 420 LINEAR CONTROL SYSTEMS Lecture 6 Mathematical Representation of Physical Systems II State Variable Models for Dynamic Systems u 1 u 2 u ṙ. Internal Variables x 1, x 2 x n y 1 y 2. y m Figure
More informationTranslational and Rotational Dynamics!
Translational and Rotational Dynamics Robert Stengel Robotics and Intelligent Systems MAE 345, Princeton University, 217 Copyright 217 by Robert Stengel. All rights reserved. For educational use only.
More informationThe Jacobian. Jesse van den Kieboom
The Jacobian Jesse van den Kieboom jesse.vandenkieboom@epfl.ch 1 Introduction 1 1 Introduction The Jacobian is an important concept in robotics. Although the general concept of the Jacobian in robotics
More informationAN INTRODUCTION TO LAGRANGE EQUATIONS. Professor J. Kim Vandiver October 28, 2016
AN INTRODUCTION TO LAGRANGE EQUATIONS Professor J. Kim Vandiver October 28, 2016 kimv@mit.edu 1.0 INTRODUCTION This paper is intended as a minimal introduction to the application of Lagrange equations
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 information41514 Dynamics of Machinery
41514 Dynamics of Machinery Theory, Experiment, Phenomenology and Industrial Applications Ilmar (Iumár) Ferreira Santos 1. Course Structure 2. Objectives 3. Theoretical and Experimental Example 4. Industrial
More informationON THE INTEGRATION OF EQUATIONS OF MOTION: FEM AND MOLECULAR DYNAMICS PROBLEMS
8th International Congress on Computational Mechanics, Volos, 1-15 July 015 ON THE INTEGRATION OF EQUATIONS OF MOTION: FEM AND MOLECULAR DYNAMICS PROBLEMS E.G. Kakouris, V.K. Koumousis Institute of Structural
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 informationLecture «Robot Dynamics»: Dynamics 2
Lecture «Robot Dynamics»: Dynamics 2 151-0851-00 V lecture: CAB G11 Tuesday 10:15 12:00, every week exercise: HG E1.2 Wednesday 8:15 10:00, according to schedule (about every 2nd week) office hour: LEE
More informationChapter 2: Rigid Bar Supported by Two Buckled Struts under Axial, Harmonic, Displacement Excitation..14
Table of Contents Chapter 1: Research Objectives and Literature Review..1 1.1 Introduction...1 1.2 Literature Review......3 1.2.1 Describing Vibration......3 1.2.2 Vibration Isolation.....6 1.2.2.1 Overview.
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 informationAnalytical Mechanics for Relativity and Quantum Mechanics
Analytical Mechanics for Relativity and Quantum Mechanics Oliver Davis Johns San Francisco State University OXPORD UNIVERSITY PRESS CONTENTS Dedication Preface Acknowledgments v vii ix PART I INTRODUCTION:
More informationTeam-Exercises for DGC 100 Modelica Course
Team-Exercises for DGC 100 Modelica Course Hubertus Tummescheit United Technologies Research Center, East Hartford, CT 06108. November 4, 2003 Abstract This document is a preliminary version and is going
More informationMATHEMATICAL PHYSICS
MATHEMATICAL PHYSICS Third Year SEMESTER 1 015 016 Classical Mechanics MP350 Prof. S. J. Hands, Prof. D. M. Heffernan, Dr. J.-I. Skullerud and Dr. M. Fremling Time allowed: 1 1 hours Answer two questions
More informationStochastic Dynamics of SDOF Systems (cont.).
Outline of Stochastic Dynamics of SDOF Systems (cont.). Weakly Stationary Response Processes. Equivalent White Noise Approximations. Gaussian Response Processes as Conditional Normal Distributions. Stochastic
More informationMSMS Basilio Bona DAUIN PoliTo
MSMS 214-215 Basilio Bona DAUIN PoliTo Problem 2 The planar system illustrated in Figure 1 consists of a bar B and a wheel W moving (no friction, no sliding) along the bar; the bar can rotate around an
More informationPARAMETER IDENTIFICATION OF A FORK-LIFT TRUCK DYNAMIC MODEL
11 th National Congress on Theoretical and Applied Mechanics, -5 Sept. 9, Borovets, Bulgaria PARAMETER IDENTIFICATION OF A FORK-IFT TRUCK DYNAMIC MODE E. S. CHANKOV Strength of Materials Dep., Technical
More information1.053J/2.003J Dynamics and Control I Fall Final Exam 18 th December, 2007
1.053J/2.003J Dynamics and Control I Fall 2007 Final Exam 18 th December, 2007 Important Notes: 1. You are allowed to use three letter-size sheets (two-sides each) of notes. 2. There are five (5) problems
More informationSIMULATION OF OPEN-LOOP MULTIBODY SYSTEMS WITH TIME-VARYING TOPOLOGY:
Proceedings of MUSME 2005, the International Symposium on Multibody Systems and Mechatronics Uberlandia (Brazil), 6-9 March 2005 Paper n.xx-musme05 SIMULATION OF OPEN-LOOP MULTIBODY SYSTEMS WITH TIME-VARYING
More informationExercise 1b: Differential Kinematics of the ABB IRB 120
Exercise 1b: Differential Kinematics of the ABB IRB 120 Marco Hutter, Michael Blösch, Dario Bellicoso, Samuel Bachmann October 5, 2016 Abstract The aim of this exercise is to calculate the differential
More informationComputational Simulation of Dynamic Response of Vehicle Tatra T815 and the Ground
IOP Conference Series: Earth and Environmental Science PAPER OPEN ACCESS Computational Simulation of Dynamic Response of Vehicle Tatra T815 and the Ground To cite this article: Jozef Vlek and Veronika
More informationEuropean Consortium for Mathematics in Industry E. Eich-Soellner and C. FUhrer Numerical Methods in Multibody Dynamics
European Consortium for Mathematics in Industry E. Eich-Soellner and C. FUhrer Numerical Methods in Multibody Dynamics European Consortium for Mathematics in Industry Edited by Leif Arkeryd, Goteborg Heinz
More informationAppendix W. Dynamic Models. W.2 4 Complex Mechanical Systems. Translational and Rotational Systems W.2.1
Appendix W Dynamic Models W.2 4 Complex Mechanical Systems W.2.1 Translational and Rotational Systems In some cases, mechanical systems contain both translational and rotational portions. The procedure
More informationCIVL 8/7117 Chapter 12 - Structural Dynamics 1/75. To discuss the dynamics of a single-degree-of freedom springmass
CIV 8/77 Chapter - /75 Introduction To discuss the dynamics of a single-degree-of freedom springmass system. To derive the finite element equations for the time-dependent stress analysis of the one-dimensional
More informationTrajectory Planning from Multibody System Dynamics
Trajectory Planning from Multibody System Dynamics Pierangelo Masarati Politecnico di Milano Dipartimento di Ingegneria Aerospaziale Manipulators 2 Manipulator: chain of
More informationDynamics and control of mechanical systems
Dynamics and control of mechanical systems Date Day 1 (03/05) - 05/05 Day 2 (07/05) Day 3 (09/05) Day 4 (11/05) Day 5 (14/05) Day 6 (16/05) Content Review of the basics of mechanics. Kinematics of rigid
More informationLECTURE 14: DEVELOPING THE EQUATIONS OF MOTION FOR TWO-MASS VIBRATION EXAMPLES
LECTURE 14: DEVELOPING THE EQUATIONS OF MOTION FOR TWO-MASS VIBRATION EXAMPLES Figure 3.47 a. Two-mass, linear vibration system with spring connections. b. Free-body diagrams. c. Alternative free-body
More informationProblem Goldstein 2-12
Problem Goldstein -1 The Rolling Constraint: A small circular hoop of radius r and mass m hoop rolls without slipping on a stationary cylinder of radius R. The only external force is that of gravity. Let
More informationThe Modeling of Single-dof Mechanical Systems
The Modeling of Single-dof Mechanical Systems Lagrange equation for a single-dof system: where: q: is the generalized coordinate; T: is the total kinetic energy of the system; V: is the total potential
More informationIDENTIFICATION OF SEMI-RIGID JOINTS IN FRAME STRUCTURES
6th ECCOMAS Conference on Smart Structures and Materials SMART213 Politecnico di Torino, 24-26 June 213 E. Carrera, F. Miglioretti and M. Petrolo (Editors) www.smart213.com IDENTIFICATION OF SEMI-RIGID
More information(W: 12:05-1:50, 50-N202)
2016 School of Information Technology and Electrical Engineering at the University of Queensland Schedule of Events Week Date Lecture (W: 12:05-1:50, 50-N202) 1 27-Jul Introduction 2 Representing Position
More informationRobust Control of Cooperative Underactuated Manipulators
Robust Control of Cooperative Underactuated Manipulators Marcel Bergerman * Yangsheng Xu +,** Yun-Hui Liu ** * Automation Institute Informatics Technology Center Campinas SP Brazil + The Robotics Institute
More informationEuler-Lagrange's equations in several variables
Euler-Lagrange's equations in several variables So far we have studied one variable and its derivative Let us now consider L2:1 More:1 Taylor: 226-227 (This proof is slightly more general than Taylor's.)
More informationStructural Dynamics. Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma).
Structural Dynamics Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma). We will now look at free vibrations. Considering the free
More informationDynamics Algorithms for Multibody Systems
International Conference on Multi Body Dynamics 2011 Vijayawada, India. pp. 351 365 Dynamics Algorithms for Multibody Systems S. V. Shah, S. K. Saha and J. K. Dutt Department of Mechanical Engineering,
More informationPLANAR KINETICS OF A RIGID BODY FORCE AND ACCELERATION
PLANAR KINETICS OF A RIGID BODY FORCE AND ACCELERATION I. Moment of Inertia: Since a body has a definite size and shape, an applied nonconcurrent force system may cause the body to both translate and rotate.
More informationProblem 1: Find the Equation of Motion from the static equilibrium position for the following systems: 1) Assumptions
Problem 1: Find the Equation of Motion from the static equilibrium position for the following systems: 1) Assumptions k 2 Wheels roll without friction k 1 Motion will not cause block to hit the supports
More informationDynamic Modeling of Rotary Double Inverted Pendulum Using Classical Mechanics
ISBN 978-93-84468-- Proceedings of 5 International Conference on Future Computational echnologies (ICFC'5) Singapore, March 9-3, 5, pp. 96-3 Dynamic Modeling of Rotary Double Inverted Pendulum Using Classical
More informationAnalytical Dynamics: Lagrange s Equation and its Application A Brief Introduction
Analytical Dynamics: Lagrange s Equation and its Application A Brief Introduction D. S. Stutts, Ph.D. Associate Professor of Mechanical Engineering Missouri University of Science and Technology Rolla,
More informationIntroduction to Control (034040) lecture no. 2
Introduction to Control (034040) lecture no. 2 Leonid Mirkin Faculty of Mechanical Engineering Technion IIT Setup: Abstract control problem to begin with y P(s) u where P is a plant u is a control signal
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 informationSection 2.2 : Electromechanical. analogies PHILIPE HERZOG AND GUILLAUME PENELET
Section 2.2 : Electromechanical analogies PHILIPE HERZOG AND GUILLAUME PENELET Paternité - Pas d'utilisation Commerciale - Partage des Conditions Initiales à l'identique : http://creativecommons.org/licenses/by-nc-sa/2.0/fr/
More informationLinear Algebra Review (Course Notes for Math 308H - Spring 2016)
Linear Algebra Review (Course Notes for Math 308H - Spring 2016) Dr. Michael S. Pilant February 12, 2016 1 Background: We begin with one of the most fundamental notions in R 2, distance. Letting (x 1,
More informationRobot Dynamics Instantaneous Kinematiccs and Jacobians
Robot Dynamics Instantaneous Kinematiccs and Jacobians 151-0851-00 V Lecture: Tuesday 10:15 12:00 CAB G11 Exercise: Tuesday 14:15 16:00 every 2nd week Marco Hutter, Michael Blösch, Roland Siegwart, Konrad
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 informationKINEMATICS & DYNAMICS
KINEMATICS & DYNAMICS BY ADVANCED DIFFERENTIAL EQUATIONS Question (**+) In this question take g = 0 ms. A particle of mass M kg is released from rest from a height H m, and allowed to fall down through
More informationRobot Dynamics II: Trajectories & Motion
Robot Dynamics II: Trajectories & Motion Are We There Yet? METR 4202: Advanced Control & Robotics Dr Surya Singh Lecture # 5 August 23, 2013 metr4202@itee.uq.edu.au http://itee.uq.edu.au/~metr4202/ 2013
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 informationA Second Order Extension of the Generalized α Method for Constrained Systems in Mechanics
A Second Order Extension of the Generalized α Method for Constrained Systems in Mechanics Laurent O. Jay 1 and Dan Negrut 2 1 Department of Mathematics, The University of Iowa, 14 MacLean Hall, Iowa City,
More informationLecture Schedule Week Date Lecture (M: 2:05p-3:50, 50-N202)
J = x θ τ = J T F 2018 School of Information Technology and Electrical Engineering at the University of Queensland Lecture Schedule Week Date Lecture (M: 2:05p-3:50, 50-N202) 1 23-Jul Introduction + Representing
More information