Robotics. Dynamics. Marc Toussaint U Stuttgart


 Milo Greer
 1 years ago
 Views:
Transcription
1 Robotics Dynamics 1D point mass, damping & oscillation, PID, dynamics of mechanical systems, EulerLagrange equation, NewtonEuler recursion, general robot dynamics, joint space control, reference trajectory following, operational space control Marc Toussaint U Stuttgart
2 Kinematic Dynamic instantly change joint velocities q: instantly change joint torques u: δq t! = J (y φ(q t )) u! =? accounts for kinematic coupling of joints but ignores inertia, forces, torques gears, stiff, all of industrial robots accounts for dynamic coupling of joints and full Newtonian physics future robots, compliant, few research robots 2/36
3 When velocities cannot be changed/set arbitrarily Examples: An air plane flying: You cannot command it to hold still in the air, or to move straight up. A car: you cannot command it to move sidewards. Your arm: you cannot command it to throw a ball with arbitrary speed (force limits). A torque controlled robot: You cannot command it to instantly change velocity (infinite acceleration/torque). What all examples have in comment: One can set controls u t (air plane s control stick, car s steering wheel, your muscles activations, torque/voltage/current send to a robot s motors) But these controls only indirectly influence the dynamics of state, x t+1 = f(x t, u t ) 3/36
4 Dynamics The dynamics of a system describes how the controls u t influence the changeofstate of the system x t+1 = f(x t, u t ) The notation x t refers to the dynamic state of the system: e.g., joint positions and velocities x t = (q t, q t ). f is an arbitrary function, often smooth 4/36
5 Outline We start by discussing a 1D point mass for 3 reasons: The most basic forcecontrolled system with inertia We can introduce and understand PID control The behavior of a point mass under PID control is a reference that we can also follow with arbitrary dynamic robots (if the dynamics are known) We discuss computing the dynamics of general robotic systems EulerLagrange equations EulerNewton method We derive the dynamic equivalent of inverse kinematics: operational space control 5/36
6 PID and a 1D point mass 6/36
7 The dynamics of a 1D point mass Start with simplest possible example: 1D point mass (no gravity, no friction, just a single mass) The state x(t) = (q(t), q(t)) is described by: position q(t) R velocity q(t) R The controls u(t) is the force we apply on the mass point The system dynamics is: q(t) = u(t)/m 7/36
8 1D point mass proportional feedback Assume current position is q. The goal is to move it to the position q. What can we do? 8/36
9 1D point mass proportional feedback Assume current position is q. The goal is to move it to the position q. What can we do? Idea 1: Always pull the mass towards the goal q : u = K p (q q) 8/36
10 1D point mass proportional feedback What s the effect of this control law? m q = u = K p (q q) q = q(t) is a function of time, this is a second order differential equation 9/36
11 1D point mass proportional feedback What s the effect of this control law? m q = u = K p (q q) q = q(t) is a function of time, this is a second order differential equation Solution: assume q(t) = a + be ωt (a nonimaginary alternative would be q(t) = a + b ɛ λt cos(ωt)) 9/36
12 1D point mass proportional feedback What s the effect of this control law? m q = u = K p (q q) q = q(t) is a function of time, this is a second order differential equation Solution: assume q(t) = a + be ωt (a nonimaginary alternative would be q(t) = a + b ɛ λt cos(ωt)) m b ω 2 e ωt = K p q K p a K p b e ωt (m b ω 2 + K p b) e ωt = K p (q a) (m b ω 2 + K p b) = 0 (q a) = 0 ω = i K p /m q(t) = q + b e i K p/m t This is an oscillation around q with amplitude b = q(0) q and frequency K p /m! 9/36
13 1D point mass proportional feedback m q = u = K p (q q) q(t) = q + b e i K p/m t Oscillation around q with amplitude b = q(0) q and frequency Kp /m 1 cos(x) /36
14 1D point mass derivative feedback Idea 2 Pull less, when we re heading the right direction already: Damp the system: u = K p (q q) + K d ( q q) q is a desired goal velocity For simplicity we set q = 0 in the following. 11/36
15 1D point mass derivative feedback What s the effect of this control law? m q = u = K p (q q) + K d (0 q) Solution: again assume q(t) = a + be ωt m b ω 2 e ωt = K p q K p a K p b e ωt K d b ωe ωt (m b ω 2 + K d b ω + K p b) e ωt = K p (q a) (m ω 2 + K d ω + K p ) = 0 (q a) = 0 ω = K d ± K 2 d 4mK p 2m q(t) = q + b e ω t The term K d 2m in ω is real exponential decay (damping) 12/36
16 1D point mass derivative feedback q(t) = q + b e ω t, ω = K d ± K 2 d 4mK p 2m Effect of the second term K 2 d 4mK p/2m in ω: K 2 d < 4mK p ω has imaginary part oscillating with frequency K p /m K 2 d /4m2 Kd 2 > 4mK p ω real strongly damped q(t) = q + be K d/2m t e i K p/m K 2 d /4m2 t Kd 2 = 4mK p second term zero only exponential decay 13/36
17 1D point mass derivative feedback illustration from O. Brock s lecture 14/36
18 1D point mass derivative feedback Alternative parameterization: Instead of the gains K p and K d it is sometimes more intuitive to set the wave length λ = 1 1 ω 0 =, K p = m/λ 2 Kp/m damping ratio ξ = ξ > 1: overdamped ξ = 1: critically dampled ξ < 1: oscillatorydamped K d = λk d 4mKp 2m, K d = 2mξ/λ q(t) = q + be ξ t/λ e i 1 ξ 2 t/λ 15/36
19 1D point mass integral feedback Idea 3 Pull if the position error accumulated large in the past: u = K p (q q) + K d ( q q) + K i t s=0 (q (s) q(s)) ds This is not a linear ODE w.r.t. x = (q, q). However, when we extend the state to x = (q, q, e) we have the ODE q = q q = u/m = K p/m(q q) + K d /m( q q) + K i /m e ė = q q (no explicit discussion here) 16/36
20 1D point mass PID control u = K p (q q) + K d ( q q) + K i t s=0 (q q(s)) ds PID control Proportional Control ( Position Control ) f K p (q q) Derivative Control ( Damping ) f K d ( q q) (ẋ = 0 damping) Integral Control ( Steady State Error ) f K i t s=0 (q (s) q(s)) ds 17/36
21 Controlling a 1D point mass lessons learnt Proportional and derivative feedback (PD control) are like adding a spring and damper to the point mass PD control is a linear control law (q, q) u = K p (q q) + K d ( q q) (linear in the dynamic system state x = (q, q)) With such linear control laws we can design approach trajectories (by tuning the gains) but no optimality principle behind such motions 18/36
22 Dynamics of mechanical systems 19/36
23 Two ways to derive dynamics equations for mechanical systems The EulerLagrange equation d L dt q L q = u Used when you want to derive analytic equations of motion ( on paper ) The NewtonEuler recursion f i = m v i, (and related algorithms) u i = I i ẇ + w Iw Algorithms that propagate forces through a kinematic tree and numerically compute the inverse dynamics u = NE(q, q, q) or forward dynamics q = f(q, q, u). 20/36
24 The EulerLagrange equation d L dt q L q = u L(q, q) is called Lagrangian and defined as L = T U where T =kinetic energy and U =potential energy. q is called generalized coordinate any coordinates such that (q, q) describes the state of the system. Joint angles in our case. u are external forces 21/36
25 The EulerLagrange equation How is this typically done? First, describe the kinematics and Jacobians for every link i: (q, q) {T W i(q), v i, w i} Recall T W i (q) = T W A T A A (q) T A B T B B (q) Further, we know that a link s velocity v i = J i q can be described via its position Jacobian. Similarly we can describe the link s angular velocity w i = Ji w q as linear in q. Second, formulate the kinetic energy T = i 1 2 miv2 i w i I iw i = i 1 2 q M i q, M i = J i Ji w m ii I i where I i = R i Ī ir i and Īi the inertia tensor in link coordinates Third, formulate the potential energies (typically independent of q) U = gm iheight(i) Fourth, compute the partial derivatives analytically to get something like }{{} u = d L dt q L q = }{{} M control inertia q + Ṁ q T q } {{ } Coriolis + U q }{{} gravity J i Ji w which relates accelerations q to the forces 22/36
26 Example: A pendulum Generalized coordinates: angle q = (θ) Kinematics: velocity of the mass: v = (l θ cos θ, 0, l θ sin θ) angular velocity of the mass: w = (0, θ, 0) Energies: T = 1 2 mv w Iw = 1 2 (ml2 + I 2 ) θ 2, U = mgl cos θ EulerLagrange equation: u = d L dt q L q = d dt (ml2 + I 2 ) θ + mgl sin θ = (ml 2 + I 2 ) θ + mgl sin θ 23/36
27 NewtonEuler recursion An algorithms that compute the inverse dynamics u = NE(q, q, q ) by recursively computing force balance at each joint: Newton s equation expresses the force acting at the center of mass for an accelerated body: f i = m v i Euler s equation expresses the torque (=control!) acting on a rigid body given an angular velocity and angular acceleration: u i = I iẇ + w Iw Forward recursion: ( kinematics) Compute (angular) velocities (v i, w i ) and accelerations ( v i, ẇ i ) for every link (via forward propagation; see geometry notes for details) Backward recursion: For the leaf links, we now know the desired accelerations q and can compute the necessary joint torques. Recurse backward. 24/36
28 Numeric algorithms for forward and inverse dynamics NewtonEuler recursion: very fast (O(n)) method to compute inverse dynamics u = NE(q, q, q ) Note that we can use this algorithm to also compute gravity terms: u = NE(q, 0, 0) = G(q) Coriolis terms: u = NE(q, q, 0) = C(q, q) q column of Intertia matrix: u = NE(q, 0, e i) = M(q) e i ArticulatedBodyDynamics: fast method (O(n)) to compute forward dynamics q = f(q, q, u) 25/36
29 Some last practical comments [demo] Use energy conservation to measure dynamic of physical simulation Physical simulation engines (developed for games): ODE (Open Dynamics Engine) Bullet (originally focussed on collision only) Physx (Nvidia) Differences of these engines to Lagrange, NE or ABD: Game engine can model much more: Contacts, tissues, particles, fog, etc (The way they model contacts looks ok but is somewhat fictional) On kinematic trees, NE or ABD are much more precise than game engines Game engines do not provide inverse dynamics, u = NE(q, q, q) Proper modelling of contacts is really really hard 26/36
30 Dynamic control of a robot 27/36
31 We previously learnt the effect of PID control on a 1D point mass Robots are not a 1D point mass Neither is each joint a 1D point mass Applying separate PD control in each joint neglects force coupling (Poor solution: Apply very high gains separately in each joint make joints stiff, as with gears.) However, knowing the robot dynamics we can transfer our understanding of PID control of a point mass to general systems 28/36
32 General robot dynamics Let (q, q) be the dynamic state and u R n the controls (typically joint torques in each motor) of a robot Robot dynamics can generally be written as: M(q) R n n C(q, q) R n G(q) R n M(q) q + C(q, q) q + G(q) = u is positive definite intertia matrix (can be inverted forward simulation of dynamics) are the centripetal and coriolis forces are the gravitational forces u are the joint torques (cf. to the EulerLagrange equation on slide 22) We often write more compactly: M(q) q + F (q, q) = u 29/36
33 Controlling a general robot From now on we jsut assume that we have algorithms to efficiently compute M(q) and F (q, q) for any (q, q) Inverse dynamics: If we know the desired q for each joint, gives the necessary torques u = M(q) q + F (q, q) Forward dynamics: If we know which torques u we apply, use q = M(q) 1 (u F (q, q)) to simulate the dynamics of the system (e.g., using RungeKutta) 30/36
34 Controlling a general robot joint space approach Where could we get the desired q from? Assume we have a nice smooth reference trajectory q ref 0:T (generated with some motion profile or alike), we can at each t read off the desired acceleration as q ref t := 1 τ [(q t+1 q t )/τ (q t q t1 )/τ] = (q t1 + q t+1 2q t )/τ 2 However, tiny errors in acceleration will accumulate greatly over time! This is Instable!! 31/36
35 Controlling a general robot joint space approach Where could we get the desired q from? Assume we have a nice smooth reference trajectory q ref 0:T (generated with some motion profile or alike), we can at each t read off the desired acceleration as q ref t := 1 τ [(q t+1 q t )/τ (q t q t1 )/τ] = (q t1 + q t+1 2q t )/τ 2 However, tiny errors in acceleration will accumulate greatly over time! This is Instable!! Choose a desired acceleration q t that implies a PDlike behavior around the reference trajectory! q t = q ref t + K p (q ref t q t ) + K d ( q ref t q t ) This is a standard and very convenient heuristic to track a reference trajectory when the robot dynamics are known: All joints will exactly behave like a 1D point particle around the reference trajectory! 31/36
36 Controlling a robot operational space approach Recall the inverse kinematics problem: We know the desired step δy (or velocity ẏ ) of the endeffector. Which step δq (or velocities q) should we make in the joints? Equivalent dynamic problem: We know how the desired acceleration ÿ of the endeffector. What controls u should we apply? 32/36
37 Operational space control Inverse kinematics: q = argmin q φ(q) y 2 C + q q 0 2 W Operational space control (one might call it Inverse task space dynamics ): u = argmin φ(q) ÿ 2 C + u 2 H u 33/36
38 Operational space control We can derive the optimum perfectly analogous to inverse kinematics We identify the minimum of a locally squared potential, using the local linearization (and approx. J = 0) We get φ(q) = d dt φ(q) d (J q + Jq) J q + 2J dt q = JM 1 (u F ) + 2J q u = T (ÿ 2 J q + T F ) with T = JM 1, (C T = H 1 T (T H 1 T ) 1 ) T = (T CT + H) 1 T C 34/36
39 Controlling a robot operational space approach Where could we get the desired ÿ from? Reference trajectory y0:t ref in operational space PDlike behavior in each operational space: ÿt = ÿt ref + K p (yt ref y t ) + K d (ẏt ref ẏ t ) illustration from O. Brock s lecture Operational space control: Let the system behave as if we could directly apply a 1D point mass behavior to the endeffector 35/36
40 Multiple tasks Recall trick last time: we defined a big kinematic map Φ(q) such that q = argmin q q q 0 2 W + Φ(q) 2 Works analogously in the dynamic case: u = argmin u u 2 H + Φ(q) 2 36/36
Robotics. Control Theory. Marc Toussaint U Stuttgart
Robotics Control Theory Topics in control theory, optimal control, HJB equation, infinite horizon case, LinearQuadratic optimal control, Riccati equations (differential, algebraic, discretetime), controllability,
More informationMultibody simulation
Multibody simulation Dynamics of a multibody system (EulerLagrange formulation) Dimitar Dimitrov Örebro University June 16, 2012 Main points covered EulerLagrange formulation manipulator inertia matrix
More information2007 Problem Topic Comment 1 Kinematics Positiontime equation Kinematics 7 2 Kinematics Velocitytime graph Dynamics 6 3 Kinematics Average velocity
2007 Problem Topic Comment 1 Kinematics Positiontime equation Kinematics 7 2 Kinematics Velocitytime graph Dynamics 6 3 Kinematics Average velocity Energy 7 4 Kinematics Free fall Collisions 3 5 Dynamics
More informationLecture «Robot Dynamics»: Dynamics 2
Lecture «Robot Dynamics»: Dynamics 2 151085100 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 informationDistance travelled time taken and if the particle is a distance s(t) along the xaxis, then its instantaneous speed is:
Chapter 1 Kinematics 1.1 Basic ideas r(t) is the position of a particle; r = r is the distance to the origin. If r = x i + y j + z k = (x, y, z), then r = r = x 2 + y 2 + z 2. v(t) is the velocity; v =
More informationLagrangian Dynamics: Generalized Coordinates and Forces
Lecture Outline 1 2.003J/1.053J Dynamics and Control I, Spring 2007 Professor Sanjay Sarma 4/2/2007 Lecture 13 Lagrangian Dynamics: Generalized Coordinates and Forces Lecture Outline Solve one problem
More informationRobotics. Path Planning. Marc Toussaint U Stuttgart
Robotics Path Planning Path finding vs. trajectory optimization, local vs. global, Dijkstra, Probabilistic Roadmaps, Rapidly Exploring Random Trees, nonholonomic systems, car system equation, pathfinding
More informationChaotic motion. Phys 750 Lecture 9
Chaotic motion Phys 750 Lecture 9 Finitedifference 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 information28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod)
28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod) θ + ω 2 sin θ = 0. Indicate the stable equilibrium points as well as the unstable equilibrium points.
More informationArticulated body dynamics
Articulated rigid bodies Articulated body dynamics Beyond human models How would you represent a pose? Quadraped animals Wavy hair Animal fur Plants Maximal vs. reduced coordinates How are things connected?
More informationSimple and Physical Pendulums Challenge Problem Solutions
Simple and Physical Pendulums Challenge Problem Solutions Problem 1 Solutions: For this problem, the answers to parts a) through d) will rely on an analysis of the pendulum motion. There are two conventional
More informationModeling and Experimentation: Compound Pendulum
Modeling and Experimentation: Compound Pendulum Prof. R.G. Longoria Department of Mechanical Engineering The University of Texas at Austin Fall 2014 Overview This lab focuses on developing a mathematical
More informationChapter 12. Recall that when a spring is stretched a distance x, it will pull back with a force given by: F = kx
Chapter 1 Lecture Notes Chapter 1 Oscillatory Motion Recall that when a spring is stretched a distance x, it will pull back with a force given by: F = kx When the mass is released, the spring will pull
More informationRobotics. Kinematics. Marc Toussaint University of Stuttgart Winter 2017/18
Robotics Kinematics 3D geometry, homogeneous transformations, kinematic map, Jacobian, inverse kinematics as optimization problem, motion profiles, trajectory interpolation, multiple simultaneous tasks,
More informationAdvanced Robotic Manipulation
Advanced Robotic Manipulation Handout CS37A (Spring 017 Solution Set # Problem 1  Redundant robot control The goal of this problem is to familiarize you with the control of a robot that is redundant with
More informationContents. Dynamics and control of mechanical systems. Focus on
Dynamics and control of mechanical systems Date Day 1 (01/08) Day 2 (03/08) Day 3 (05/08) Day 4 (07/08) Day 5 (09/08) Day 6 (11/08) Content Review of the basics of mechanics. Kinematics of rigid bodies
More informationSymmetries 2  Rotations in Space
Symmetries 2  Rotations in Space This symmetry is about the isotropy of space, i.e. space is the same in all orientations. Thus, if we continuously rotated an entire system in space, we expect the system
More informationAssignments VIII and IX, PHYS 301 (Classical Mechanics) Spring 2014 Due 3/21/14 at start of class
Assignments VIII and IX, PHYS 301 (Classical Mechanics) Spring 2014 Due 3/21/14 at start of class Homeworks VIII and IX both center on Lagrangian mechanics and involve many of the same skills. Therefore,
More informationColumbia University Department of Physics QUALIFYING EXAMINATION
Columbia University Department of Physics QUALIFYING EXAMINATION Monday, January 9, 2017 11:00AM to 1:00PM Classical Physics Section 1. Classical Mechanics Two hours are permitted for the completion of
More informationDamped harmonic motion
Damped harmonic motion March 3, 016 Harmonic motion is studied in the presence of a damping force proportional to the velocity. The complex method is introduced, and the different cases of underdamping,
More information7 Pendulum. Part II: More complicated situations
MATH 35, by T. Lakoba, University of Vermont 60 7 Pendulum. Part II: More complicated situations In this Lecture, we will pursue two main goals. First, we will take a glimpse at a method of Classical Mechanics
More informationLinearization problem. The simplest example
Linear Systems Lecture 3 1 problem Consider a nonlinear timeinvariant system of the form ( ẋ(t f x(t u(t y(t g ( x(t u(t (1 such that x R n u R m y R p and Slide 1 A: f(xu f(xu g(xu and g(xu exist and
More informationChapter 15  Oscillations
The pendulum of the mind oscillates between sense and nonsense, not between right and wrong. Carl Gustav Jung David J. Starling Penn State Hazleton PHYS 211 Oscillatory motion is motion that is periodic
More informationForces of Constraint & Lagrange Multipliers
Lectures 30 April 21, 2006 Written or last updated: April 21, 2006 P442 Analytical Mechanics  II Forces of Constraint & Lagrange Multipliers c Alex R. Dzierba Generalized Coordinates Revisited Consider
More informationExponential Controller for Robot Manipulators
Exponential Controller for Robot Manipulators Fernando Reyes Benemérita Universidad Autónoma de Puebla Grupo de Robótica de la Facultad de Ciencias de la Electrónica Apartado Postal 542, Puebla 7200, México
More informationControl Systems I. Lecture 2: Modeling. Suggested Readings: Åström & Murray Ch. 23, Guzzella Ch Emilio Frazzoli
Control Systems I Lecture 2: Modeling Suggested Readings: Åström & Murray Ch. 23, Guzzella Ch. 23 Emilio Frazzoli Institute for Dynamic Systems and Control DMAVT ETH Zürich September 29, 2017 E. Frazzoli
More informationRotational Kinetic Energy
Lecture 17, Chapter 10: Rotational Energy and Angular Momentum 1 Rotational Kinetic Energy Consider a rigid body rotating with an angular velocity ω about an axis. Clearly every point in the rigid body
More informationState Space Representation
ME Homework #6 State Space Representation Last Updated September 6 6. From the homework problems on the following pages 5. 5. 5.6 5.7. 5.6 Chapter 5 Homework Problems 5.6. Simulation of Linear and Nonlinear
More informationPhysics A  PHY 2048C
Physics A  PHY 2048C and 11/15/2017 My Office Hours: Thursday 2:003:00 PM 212 Keen Building Warmup Questions 1 Did you read Chapter 12 in the textbook on? 2 Must an object be rotating to have a moment
More informationOscillatory Motion SHM
Chapter 15 Oscillatory Motion SHM Dr. Armen Kocharian Periodic Motion Periodic motion is motion of an object that regularly repeats The object returns to a given position after a fixed time interval A
More informationparticle p = m v F ext = d P = M d v cm dt
Lecture 11: Momentum and Collisions; Introduction to Rotation 1 REVIEW: (Chapter 8) LINEAR MOMENTUM and COLLISIONS The first new physical quantity introduced in Chapter 8 is Linear Momentum Linear Momentum
More informationMEM04: Rotary Inverted Pendulum
MEM4: Rotary Inverted Pendulum Interdisciplinary Automatic Controls Laboratory  ME/ECE/CHE 389 April 8, 7 Contents Overview. Configure ELVIS and DC Motor................................ Goals..............................................3
More informationChapter 7 Control. Part State Space Control. Mobile Robotics  Prof Alonzo Kelly, CMU RI
Chapter 7 Control Part 2 7.2 State Space Control 1 7.2 State Space Control Outline 7.2.1 Introduction 7.2.2 State Space Feedback Control 7.2.3 Example: Robot Trajectory Following 7.2.4 Perception Based
More informationPractice Problems for Exam 2 Solutions
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Physics Physics 8.01 Fall Term 008 Practice Problems for Exam Solutions Part I Concept Questions: Circle your answer. 1) A springloaded toy dart gun
More informationRutgers University Department of Physics & Astronomy. 01:750:271 Honors Physics I Fall Lecture 20 JJ II. Home Page. Title Page.
Rutgers University Department of Physics & Astronomy 01:750:271 Honors Physics Fall 2015 Lecture 20 Page 1 of 31 1. No quizzes during Thanksgiving week. There will be recitation according to the regular
More informationCS184: Computer Graphics
CS184: Computer Graphics Lecture #25: Rigid Body Simulations Tobias Pfaff 537 Soda (Visual Computing Lab) tpfaff@berkeley.edu Reminder Final project presentations next week! Game Physics Types of Materials
More informationEXAMPLE 2: CLASSICAL MECHANICS: Worked examples. b) Position and velocity as integrals. Michaelmas Term Lectures Prof M.
CLASSICAL MECHANICS: Worked examples Michaelmas Term 2006 4 Lectures Prof M. Brouard EXAMPLE 2: b) Position and velocity as integrals Calculate the position of a particle given its time dependent acceleration:
More informationMath Assignment 5
Math 2280  Assignment 5 Dylan Zwick Fall 2013 Section 3.41, 5, 18, 21 Section 3.51, 11, 23, 28, 35, 47, 56 Section 3.61, 2, 9, 17, 24 1 Section 3.4  Mechanical Vibrations 3.4.1  Determine the period
More informationThe dynamics of a Mobile Inverted Pendulum (MIP)
The dynamics of a Mobile Inverted Pendulum (MIP) 1 Introduction Saam Ostovari, Nick Morozovsky, Thomas Bewley UCSD Coordinated Robotics Lab In this document, a Mobile Inverted Pendulum (MIP) is a robotic
More informationOscillations Simple Harmonic Motion
Oscillations Simple Harmonic Motion Lana Sheridan De Anza College Dec 1, 2017 Overview oscillations simple harmonic motion (SHM) spring systems energy in SHM pendula damped oscillations Oscillations and
More informationDynamics. 1 Copyright c 2015 Roderic Grupen
Dynamics The branch of physics that treats the action of force on bodies in motion or at rest; kinetics, kinematics, and statics, collectively. Websters dictionary Outline Conservation of Momentum Inertia
More informationParticle Systems. CSE169: Computer Animation Instructor: Steve Rotenberg UCSD, Winter 2017
Particle Systems CSE169: Computer Animation Instructor: Steve Rotenberg UCSD, Winter 2017 Particle Systems Particle systems have been used extensively in computer animation and special effects since their
More informationLecture V: The gameengine loop & Time Integration
Lecture V: The gameengine loop & Time Integration The Basic GameEngine Loop Previous state: " #, %(#) ( #, )(#) Forces (#) Integrate velocities and positions Resolve Interpenetrations Perbody change
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 information1 Lyapunov theory of stability
M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability
More informationAA242B: MECHANICAL VIBRATIONS
AA242B: MECHANICAL VIBRATIONS 1 / 50 AA242B: MECHANICAL VIBRATIONS Undamped Vibrations of ndof Systems These slides are based on the recommended textbook: M. Géradin and D. Rixen, Mechanical Vibrations:
More informationApplication Note #3413
Application Note #3413 Manual Tuning Methods Tuning the controller seems to be a difficult task to some users; however, after getting familiar with the theories and tricks behind it, one might find the
More informationModeling and Solving Constraints. Erin Catto Blizzard Entertainment
Modeling and Solving Constraints Erin Catto Blizzard Entertainment Basic Idea Constraints are used to simulate joints, contact, and collision. We need to solve the constraints to stack boxes and to keep
More informationCurves in the configuration space Q or in the velocity phase space Ω satisfying the EulerLagrange (EL) equations,
Physics 6010, Fall 2010 Hamiltonian Formalism: Hamilton s equations. Conservation laws. Reduction. Poisson Brackets. Relevant Sections in Text: 8.1 8.3, 9.5 The Hamiltonian Formalism We now return to formal
More informationHarmonic Oscillator. MassSpring Oscillator Resonance The Pendulum. Physics 109 Experiment Number 12
Harmonic Oscillator MassSpring Oscillator Resonance The Pendulum Physics 109 Experiment Number 12 Outline Simple harmonic motion The vertical massspring system Driven oscillations and resonance The pendulum
More informationImaginary. Axis. Real. Axis
Name ME6 Final. I certify that I upheld the Stanford Honor code during this exam Monday December 2, 25 3:36:3 p.m. ffl Print your name and sign the honor code statement ffl You may use your course notes,
More informationChapter 8. Dynamics II: Motion in a Plane
Chapter 8. Dynamics II: Motion in a Plane Chapter Goal: To learn how to solve problems about motion in a plane. Slide 82 Chapter 8 Preview Slide 83 Chapter 8 Preview Slide 84 Chapter 8 Preview Slide
More informationPhysics 351, Spring 2015, Homework #5. Due at start of class, Friday, February 20, 2015 Course info is at positron.hep.upenn.
Physics 351, Spring 2015, Homework #5. Due at start of class, Friday, February 20, 2015 Course info is at positron.hep.upenn.edu/p351 When you finish this homework, remember to visit the feedback page
More informationAutomatic Control 2. Nonlinear systems. Prof. Alberto Bemporad. University of Trento. Academic year
Automatic Control 2 Nonlinear systems Prof. Alberto Bemporad University of Trento Academic year 20102011 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 20102011 1 / 18
More informationPhys101 Lectures 9 and 10 Conservation of Mechanical Energy
Phys101 Lectures 9 and 10 Conservation of Mechanical Energy Key points: Conservative and Nonconservative Forces Potential Energy Generalized workenergy principle Mechanical Energy and Its Conservation
More informationAdvanced Robotic Manipulation
Lecture Notes (CS327A) Advanced Robotic Manipulation Oussama Khatib Stanford University Spring 2005 ii c 2005 by Oussama Khatib Contents 1 Spatial Descriptions 1 1.1 Rigid Body Configuration.................
More informationName: Fall 2014 CLOSED BOOK
Name: Fall 2014 1. Rod AB with weight W = 40 lb is pinned at A to a vertical axle which rotates with constant angular velocity ω =15 rad/s. The rod position is maintained by a horizontal wire BC. Determine
More informationNonlinear PD Controllers with Gravity Compensation for Robot Manipulators
BULGARIAN ACADEMY OF SCIENCES CYBERNETICS AND INFORMATION TECHNOLOGIES Volume 4, No Sofia 04 Print ISSN: 3970; Online ISSN: 34408 DOI: 0.478/cait0400 Nonlinear PD Controllers with Gravity Compensation
More informationThe Principle of Least Action
The Principle of Least Action In their neverending search for general principles, from which various laws of Physics could be derived, physicists, and most notably theoretical physicists, have often made
More informationVibrations Qualifying Exam Study Material
Vibrations Qualifying Exam Study Material The candidate is expected to have a thorough understanding of engineering vibrations topics. These topics are listed below for clarification. Not all instructors
More informationFeedback Basics. David M. Auslander Mechanical Engineering University of California at Berkeley. copyright 1998, D.M. Auslander
Feedback Basics David M. Auslander Mechanical Engineering University of California at Berkeley copyright 1998, D.M. Auslander 1 I. Feedback Control Context 2 What is Feedback Control? Measure desired behavior
More informationFor a rigid body that is constrained to rotate about a fixed axis, the gravitational torque about the axis is
Experiment 14 The Physical Pendulum The period of oscillation of a physical pendulum is found to a high degree of accuracy by two methods: theory and experiment. The values are then compared. Theory For
More informationTranslational vs Rotational. m x. Connection Δ = = = = = = Δ = = = = = = Δ =Δ = = = = = 2 / 1/2. Work
Translational vs Rotational / / 1/ Δ m x v dx dt a dv dt F ma p mv KE mv Work Fd / / 1/ θ ω θ α ω τ α ω ω τθ Δ I d dt d dt I L I KE I Work / θ ω α τ Δ Δ c t s r v r a v r a r Fr L pr Connection Translational
More informationChapter 9 Notes. x cm =
Chapter 9 Notes Chapter 8 begins the discussion of rigid bodies, a system of particles with fixed relative positions. Previously we have dealt with translation of a particle: if a rigid body does not rotate
More informationThe Acrobot and CartPole
C H A P T E R 3 The Acrobot and CartPole 3.1 INTRODUCTION A great deal of work in the control of underactuated systems has been done in the context of lowdimensional model systems. These model systems
More informationNatural and artificial constraints
FORCE CONTROL Manipulator interaction with environment Compliance control Impedance control Force control Constrained motion Natural and artificial constraints Hybrid force/motion control MANIPULATOR INTERACTION
More informationPhysics 2101 S c e t c i cti n o 3 n 3 March 31st Announcements: Quiz today about Ch. 14 Class Website:
Physics 2101 Section 3 March 31 st Announcements: Quiz today about Ch. 14 Class Website: http://www.phys.lsu.edu/classes/spring2010/phys2101 3/ http://www.phys.lsu.edu/~jzhang/teaching.html Simple Harmonic
More informationt = g = 10 m/s 2 = 2 s T = 2π g
Annotated Answers to the 1984 AP Physics C Mechanics Multiple Choice 1. D. Torque is the rotational analogue of force; F net = ma corresponds to τ net = Iα. 2. C. The horizontal speed does not affect the
More informationDr. Ian R. Manchester
Dr Ian R. Manchester Week Content Notes 1 Introduction 2 Frequency Domain Modelling 3 Transient Performance and the splane 4 Block Diagrams 5 Feedback System Characteristics Assign 1 Due 6 Root Locus
More informationEQUIVALENT SINGLEDEGREEOFFREEDOM SYSTEM AND FREE VIBRATION
1 EQUIVALENT SINGLEDEGREEOFFREEDOM SYSTEM AND FREE VIBRATION The course on Mechanical Vibration is an important part of the Mechanical Engineering undergraduate curriculum. It is necessary for the development
More informationConcept Question: Normal Force
Concept Question: Normal Force Consider a person standing in an elevator that is accelerating upward. The upward normal force N exerted by the elevator floor on the person is 1. larger than 2. identical
More informationT = 1 2 mv2. For a system of N particles, the kinetic energy of the system is obtained by adding up the individual kinetic energies: T = 2 m ivi 2.
Physics 3550, Fall 2011 Energy. Relevant Sections in Text: Chapter 4 Kinetic Energy. I probably don t need to impress upon you the importance of the notion of energy in analyzing/understanding physical
More informationHarmonic Oscillator. Outline. Oscillatory Motion or Simple Harmonic Motion. Oscillatory Motion or Simple Harmonic Motion
Harmonic Oscillator MassSpring Oscillator Resonance The Pendulum Physics 109, Class Period 13 Experiment Number 11 in the Physics 121 Lab Manual (page 65) Outline Simple harmonic motion The vertical massspring
More informationPHY218 SPRING 2016 Review for Final Exam: Week 14 Final Review: Chapters 111, 1314
Final Review: Chapters 111, 1314 These are selected problems that you are to solve independently or in a team of 23 in order to better prepare for your Final Exam 1 Problem 1: Chasing a motorist This
More informationApplications of SecondOrder Differential Equations
Applications of SecondOrder Differential Equations ymy/013 Building Intuition Even though there are an infinite number of differential equations, they all share common characteristics that allow intuition
More informationDIFFERENTIAL KINEMATICS. Geometric Jacobian. Analytical Jacobian. Kinematic singularities. Kinematic redundancy. Inverse differential kinematics
DIFFERENTIAL KINEMATICS relationship between joint velocities and endeffector velocities Geometric Jacobian Analytical Jacobian Kinematic singularities Kinematic redundancy Inverse differential kinematics
More informationExamination paper for TMA4195 Mathematical Modeling
Department of Mathematical Sciences Examination paper for TMA4195 Mathematical Modeling Academic contact during examination: Elena Celledoni Phone: 48238584, 73593541 Examination date: 11th of December
More informationLagging Pendulum. Martin Ga
02 Lagging Pendulum Martin Ga 1 Task A pendulum consists of a strong thread and a bob. When the pivot of the pendulum starts moving along a horizontal circumference, the bob starts tracing a circle which
More informationPHY 3221 Fall Homework Problems. Instructor: Yoonseok Lee. Submit only HW s. EX s are additional problems that I encourage you to work on.
PHY 3221 Fall 2012 Homework Problems Instructor: Yoonseok Lee Submit only HW s. EX s are additional problems that I encourage you to work on. Week 1: August 2224, Due August 27 (nothing to submit) EX:
More informationChapter 14. Oscillations and Resonance RESONANCE
Chapter 14 Oscillations and Resonance CHAPTER 14 RESONANCE OSCILLATIONS AND Oscillations and vibrations play a more significant role in our lives than we realize. When you strike a bell, the metal vibrates,
More informationLagrangian and Hamiltonian Mechanics (Symon Chapter Nine)
Lagrangian and Hamiltonian Mechanics (Symon Chapter Nine Physics A301 Spring 2005 Contents 1 Lagrangian Mechanics 3 1.1 Derivation of the Lagrange Equations...................... 3 1.1.1 Newton s Second
More informationECE 422/522 Power System Operations & Planning/Power Systems Analysis II : 7  Transient Stability
ECE 4/5 Power System Operations & Planning/Power Systems Analysis II : 7  Transient Stability Spring 014 Instructor: Kai Sun 1 Transient Stability The ability of the power system to maintain synchronism
More informationFinite Difference and Finite Element Methods
Finite Difference and Finite Element Methods Georgy Gimel farb COMPSCI 369 Computational Science 1 / 39 1 Finite Differences Difference Equations 3 Finite Difference Methods: Euler FDMs 4 Finite Element
More informationChapter 4. Oscillatory Motion. 4.1 The Important Stuff Simple Harmonic Motion
Chapter 4 Oscillatory Motion 4.1 The Important Stuff 4.1.1 Simple Harmonic Motion In this chapter we consider systems which have a motion which repeats itself in time, that is, it is periodic. In particular
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 informationOscillations. PHYS 101 Previous Exam Problems CHAPTER. Simple harmonic motion Massspring system Energy in SHM Pendulums
PHYS 101 Previous Exam Problems CHAPTER 15 Oscillations Simple harmonic motion Massspring system Energy in SHM Pendulums 1. The displacement of a particle oscillating along the x axis is given as a function
More informationLaboratory 11 Control Systems Laboratory ECE3557. State Feedback Controller for Position Control of a Flexible Joint
Laboratory 11 State Feedback Controller for Position Control of a Flexible Joint 11.1 Objective The objective of this laboratory is to design a full state feedback controller for endpoint position control
More informationPHY131H1S  Class 20. Preclass reading quiz on Chapter 12
PHY131H1S  Class 20 Today: Gravitational Torque Rotational Kinetic Energy Rolling without Slipping Equilibrium with Rotation Rotation Vectors Angular Momentum Preclass reading quiz on Chapter 12 1 Last
More informationMechanics II. Which of the following relations among the forces W, k, N, and F must be true?
Mechanics II 1. By applying a force F on a block, a person pulls a block along a rough surface at constant velocity v (see Figure below; directions, but not necessarily magnitudes, are indicated). Which
More informationComplex Numbers. The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition,
Complex Numbers Complex Algebra The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition, and (ii) complex multiplication, (x 1, y 1 )
More informationVideo 1.1 Vijay Kumar and Ani Hsieh
Video 1.1 Vijay Kumar and Ani Hsieh 1 Robotics: Dynamics and Control Vijay Kumar and Ani Hsieh University of Pennsylvania 2 Why? Robots live in a physical world The physical world is governed by the laws
More informationNeural Networks Lecture 10: Fault Detection and Isolation (FDI) Using Neural Networks
Neural Networks Lecture 10: Fault Detection and Isolation (FDI) Using Neural Networks H.A. Talebi Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Winter 2011.
More informationPH 120 Project # 2: Pendulum and chaos
PH 120 Project # 2: Pendulum and chaos Due: Friday, January 16, 2004 In PH109, you studied a simple pendulum, which is an effectively massless rod of length l that is fixed at one end with a small mass
More informationLecture 9  Rotational Dynamics
Lecture 9  Rotational Dynamics A Puzzle... Angular momentum is a 3D vector, and changing its direction produces a torque τ = dl. An important application in our daily lives is that bicycles don t fall
More informationComplex Numbers. The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition,
Complex Numbers Complex Algebra The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition, and (ii) complex multiplication, (x 1, y 1 )
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 informationHomework 1. Due Thursday, January 21
Homework 1. Due Thursday, January 21 Problem 1. Rising Snake A snake of length L and linear mass density ρ rises from the table. It s head is moving straight up with the constant velocity v. What force
More informationMath 312 Lecture Notes Linear Twodimensional Systems of Differential Equations
Math 2 Lecture Notes Linear Twodimensional Systems of Differential Equations Warren Weckesser Department of Mathematics Colgate University February 2005 In these notes, we consider the linear system of
More information3 DOF Scara type Robot Manipulator using Mamdani Based Fuzzy Controller
659 3 DOF Scara type Robot Manipulator using Mamdani Based Fuzzy Controller Nitesh Kumar Jaiswal *, Vijay Kumar ** *(Department of Electronics and Communication Engineering, Indian Institute of Technology,
More informationVideo 3.1 Vijay Kumar and Ani Hsieh
Video 3.1 Vijay Kumar and Ani Hsieh Robo3x1.3 1 Dynamics of Robot Arms Vijay Kumar and Ani Hsieh University of Pennsylvania Robo3x1.3 2 Lagrange s Equation of Motion Lagrangian Kinetic Energy Potential
More information