Rigid body simulation. Once we consider an object with spatial extent, particle system simulation is no longer sufficient

Size: px
Start display at page:

Download "Rigid body simulation. Once we consider an object with spatial extent, particle system simulation is no longer sufficient"

Transcription

1 Rigid body dynamics

2 Rigid body simulation Once we consider an object with spatial extent, particle system simulation is no longer sufficient

3 Rigid body simulation Unconstrained system no contact Constrained system collision and contact

4 Problems Performance is important!

5 Problems Control is difficult!

6 Particle simulation Y(t) = [ x(t) v(t) ] Position in phase space Ẏ(t) = [ v(t) f(t)/m ] Velocity in phase space

7 Rigid body concepts Translation Position Linear velocity Mass Linear momentum Force Rotation Orientation Angular velocity Inertia tensor Angular momentum Torque

8 Position and orientation Linear and angular velocity Mass and Inertia Force and torques Simulation

9 Position and orientation Translation of the body x(t) = x y z Rotation of the body R(t) = r xx r yx r zx r xy r yy r zy r xz r yz r zz x(t) and R(t) are called spatial variables of a rigid body

10 Body space Body space y 0 A fixed and unchanged space where the shape of a rigid body is defined z 0 r 0i x 0 The geometric center of the rigid body lies at the origin of the body space

11 Position and orientation Body space World space y0 R(t) y0 r0i z0 y x0 z0 z x(t) x r0i x0

12 Position and orientation Use x(t) and R(t) to transform the body space into world space World space y 0 What s the world coordinate of an arbitrary point r 0i on the body? R(t) z 0 r 0i x 0 y x(t) r i (t) =x(t)+r(t)r 0i z x

13 Position and orientation Assume the rigid body has uniform density, what is the physical meaning of x(t)? center of mass over time What is the physical meaning of R(t)? it s a bit tricky

14 Position and orientation Consider the x-axis in body space, (1, 0, 0), what is the direction of this vector in world space at time t? R(t) = r xx r xy r xz which is the first column of R(t) R(t) represents the directions of x, y, and z axes of the body space in world space at time t

15 Position and orientation So x(t) and R(t) define the position and the orientation of the body at time t Next we need to define how the position and orientation change over time

16 Position and orientation Linear and angular velocity Mass and Inertia Force and torques Simulation

17 Linear velocity Since x(t) is the position of the center of mass in world space, ẋ(t) is the velocity of the center of mass in world space v(t) =ẋ(t)

18 Angular velocity If we freeze the position of the COM in space then any movement is due to the body spinning about some axis that passes through the COM Otherwise, the COM would itself be moving

19 Angular velocity We describe that spin as a vector ω(t) Direction of ω(t)? Magnitude of ω(t)? Linear position and velocity are related by v(t) = d dt x(t) How are angular position (orientation) and velocity related?

20 Angular velocity How are R(t) and ω(t) related? Hint: Consider a vector c(t) at time t specified in world space, how do we represent ċ(t) in terms of ω(t) ω(t) ċ(t) = b ω(t) = ω(t) b a b c ċ ċ(t) =ω(t) b = ω(t) b + ω(t) a ċ(t) =ω(t) c(t)

21 Angular velocity Given the physical meaning of column of Ṙ(t) mean? R(t), what does each At time t, the direction of x-axis of the rigid body in world space is the first column of R(t) r xx r xy r xz At time t, what is the derivative of the first column of R(t)? r xx r xy r xz 5 =!(t) 4 r xx r xy r xz 5

22 Angular velocity Ṙ(t) = ω(t) r xx r xy r xz ω(t) r yx r yy r yz ω(t) r zx r zy r zz This is the relation between angular velocity and the orientation, but it is too cumbersome We can use a trick to simplify this expression

23 Angular velocity Consider two 3 by 1 vectors: a and b, the cross product of them is a b = a yb z b y a z a x b z + b x a z a x b y b x a y Given a, let s define a* to be a skew symmetric matrix 0 a z a y a z 0 a x a y a x 0 then a b = 0 a z a y a z 0 a x a y a x 0 b x b y b z = a b

24 Angular velocity Ṙ(t) = ω(t) r xx r xx r xy r xy ω(t) r yx r yx r yy r yy ω(t) ω(t) r zx r zy r zx r zy r xz r xz r yz r yz r zz r zz = ω(t) R(t) Vector relation: ċ(t) =ω(t) c(t) Matrix relation: Ṙ = ω(t) R(t)

25 Perspective of particles Imagine a rigid body is composed of a large number of small particles the particles are indexed from 1 to N each particle has a constant location r 0i in body space the location of i-th particle in world space at time t is r i (t) =x(t)+r(t)r 0i

26 Velocity of a particle ṙ(t) = d dt r(t) =ω R(t)r 0i + v(t) = ω (R(t)r 0i + x(t) x(t)) + v(t) = ω (r i (t) x(t)) + v(t) ṙ i (t) =ω (r i (t) x(t)) + v(t) angular component linear component

27 Velocity of a particle ṙ i (t) =ω (r i (t) x(t)) + v(t) y 0 ω(t) ω(t) (r i (t) x(t)) x(t) r i (t) ṙ i (t) v(t) y z 0 v(t) x 0 z x

28 Position and orientation Linear and angular velocity Mass and Inertia Force and torques Simulation

29 Mass The mass of the i-th particle is m i Mass M = N i=1 m i Center of mass in world space mi r i (t) M What about center of mass in body space? (0, 0, 0)

30 Center of mass Proof that the center of mass at time t in word space is x(t) mi r i (t) M = = x(t)

31 Inertia tensor Inertia tensor describes how the mass of a rigid body is distributed relative to the center of mass I = X i m i (r 0 2 iy + r 0 2 iz ) m i rix 0 r0 iy m i rix 0 r0 iz m i riy 0 r0 iz m i r 0 iy r0 ix m i (r 0 2 ix + r 0 2 iz ) m i r 0 iz r0 ix m i r 0 iz r0 iy m i (r 0 2 ix + r 0 2 iy ) r i = r i (t) x(t) I(t) depends on the orientation of a body, but not the translation For an actual implementation, we replace the finite sum with the integrals over a body s volume in world space

32 Inertia tensor Inertia tensors vary in world space over time But are constant in the body space Pre-compute the integral part in the body space to save time

33 Inertia tensor Pre-compute I body that does not vary over time I(t) = X m i r 0T i r 0 i m i r 02 ix m i r 0 ix r0 iy m i r 0 ix r0 iz m i r 0 iy r0 ix m i r 02 iy m i r 0 iy r0 iz m i r 0 iz r0 ix m i r 0 iz r0 iy m i r 02 iz 3 5 I(t) = m i (r T i r i)1 r ir T i ) = m i ((R(t)r 0i ) T (R(t)r 0i )1 (R(t)r 0i )(R(t)r 0i ) T ) = m i (R(t)(r T 0ir 0i )R(t) T 1 R(t)r 0i r T 0iR(t) T ) ( ) = R(t) mi ((r T 0ir 0i )1 r 0i r T 0i) R(t) T I(t) =R(t)I body R(t) T I body = i m i ((r T 0ir 0i )1 r 0i r T 0i)

34 Inertia tensor Pre-compute I body that does not vary over time I(t) = X m i r 0T i r 0 i m i r 02 ix m i r 0 ix r0 iy m i r 0 ix r0 iz m i r 0 iy r0 ix m i r 02 iy m i r 0 iy r0 iz m i r 0 iz r0 ix m i r 0 iz r0 iy m i r 02 iz 3 5 I(t) = m i (r T i r i)1 r ir T i ) = m i ((R(t)r 0i ) T (R(t)r 0i )1 (R(t)r 0i )(R(t)r 0i ) T ) = m i (R(t)(r T 0ir 0i )R(t) T 1 R(t)r 0i r T 0iR(t) T ) ( ) = R(t) mi ((r T 0ir 0i )1 r 0i r T 0i) R(t) T I(t) =R(t)I body R(t) T I body = i m i ((r T 0ir 0i )1 r 0i r T 0i)

35 Approximate inertia tensor Bounding boxes Pros: simple Cons: inaccurate

36 Approximate inertia tensor Point sampling Pros: simple, fairly accurate Cons: expensive, requires volume test

37 Approximate inertia tensor Green s theorem D F ds = Pros: simple, exact Cons: require boundary representation D ( F) da

38 Position and orientation Linear and angular velocity Mass and Inertia Force and torques Simulation

39 Force and torque F i (t) denotes the total force from external forces acting on the i-th particle at time t y 0 τ(t) =(r i (t) x(t)) F i (t) F(t) = i F i (t) x(t) r i (t) F i (t) τ(t) = i (r i (t) x(t)) F i (t) z 0 x 0 y z x

40 Force and torque F(t) conveys no information about where the various forces acted on the body!(t) contains the information about the distribution of the forces over the body Which one depends on the location of the particle relative to the center of mass?

41 Linear momentum P(t) = i m i ṙ i (t) = i m i v(t)+ω(t) i m i (r i (t) x(t)) = Mv(t) Total linear moment of the rigid body is the same as if the body was simply a particle with mass M and velocity v(t)

42 Angular momentum Similar to linear momentum, angular momentum is defined as L(t) =I(t)ω(t) Does L(t) depend on the translational effect x(t)? Does L(t) depend on the rotational effect R(t)? What about P(t)?

43 Derivative of momentum Change in linear momentum is equivalent to the total forces acting on the rigid body Ṗ(t) =M v(t) =F(t) The relation between angular momentum and the total torque is analogous to the linear case L(t) =τ(t)

44 Derivative of momentum Proof L(t) =τ(t) = r i F i m i r i F i = m i ( v ṙ i ω r i ω) F i = 0 r i m i ( v ṙ i ω r i ω) r i F i = 0 ( ) ( ) mi r i ṙ i ω mi r i r i ω = τ mi r i r i = m i ((r T i r i)1 r ir T i )=I(t) ( ) mi r i ṙ i ω + I(t) ω = τ İ(t) = d dt mi r i r i = m i r i ṙ i m i ṙ i r i İ(t)ω + I(t) ω = d dt (I(t)ω) = L(t) =τ

45 Position and orientation Linear and angular velocity Mass and Inertia Force and torques Simulation

46 Equation of motion Y(t) = x(t) R(t) P(t) L(t) position orientation linear momentum angular momentum d dt Y(t) = v(t) ω(t) R(t) F(t) τ(t) Constants: M and I body v(t) = P(t) M I(t) =R(t)I body R(t) T ω(t) =I(t) 1 L(t)

47 Momentum vs. velocity Why do we use momentum in the phase space instead of velocity? Because the relation of angular momentum and torque is simple Because the angular momentum is constant when there is no torques acting on the object Use linear momentum P(t) to be consistent with angular velocity and acceleration

48 Example: 1. compute the I body in body space y ( x0 2, y ) 0 2, z 0 2 I body = M 12 y z x z x y 2 0 ( x 0 2, y 0 2, z ) 0 2 z x

49 Example: 1. compute the I body in body space y 2. rotation free movement ( 3, 0, 2) (3, 0, 2) x F z F

50 Example: 1. compute the I body in body space y 2. rotation free movement 3. translation free movement ( 3, 0, 2) (3, 0, 2) x F z F

51 Force vs. torque puzzle 10 sec later F energy = 0 F energy = 1 2 MvT v Suppose a force F acts on the block at the center of mass for 10 seconds. Since there is no torque acting on the block, the body will only acquire linear velocity v after 10 seconds. The kinetic energy will be 1 2 MvT v F energy = 0 Now, consider the same force acting off-center to the body for 10 seconds. Since it is the same force, the velocity of the center of mass after 10 seconds is the same v. However, the block will also pick up some angular velocity ω. The kinetic energy will be 1 2 MvT v ωt Iω If identical forces push the block in both cases, how can the energy of the block be different?

52 Notes on implementation Using quaternion instead of transformation matrix more compact representation less numerical drift

53 Quaternion q(t) = w x y z q(t) = 1 2 apple 0!(t) q(t) quaternion multiplication

54 Equation of motion Y(t) = x(t) q(t) P(t) L(t) position orientation linear momentum angular momentum d dt Y(t) = apple v(t) 0!(t) F(t) (t) q(t) Constants: M and I body v(t) = P(t) M I(t) =R(t)I body R(t) T ω(t) =I(t) 1 L(t)

55 Exercise Consider a 3D sphere with radius 1m, mass 1kg, and inertia I body. The initial linear and angular velocity are both zero. The initial position and the initial orientation are x 0 and R 0. The forces applied on the sphere include gravity (g) and an initial push F applied at point p. Note that F is only applied for one time step at t 0. If we use Explicit Euler method with time step h to integrate, what are the position and the orientation of the sphere at t 2? Use the actual numbers defined as below to compute your solution (except for g and h).

56 What s next?

57 Collision, collision, collision Read: Unconstrained rigid body dynamics by Baraff and Witkin

Lesson Rigid Body Dynamics

Lesson Rigid Body Dynamics Lesson 8 Rigid Body Dynamics Lesson 8 Outline Problem definition and motivations Dynamics of rigid bodies The equation of unconstrained motion (ODE) User and time control Demos / tools / libs Rigid Body

More information

Rigid body dynamics. Basilio Bona. DAUIN - Politecnico di Torino. October 2013

Rigid body dynamics. Basilio Bona. DAUIN - Politecnico di Torino. October 2013 Rigid body dynamics Basilio Bona DAUIN - Politecnico di Torino October 2013 Basilio Bona (DAUIN - Politecnico di Torino) Rigid body dynamics October 2013 1 / 16 Multiple point-mass bodies Each mass is

More information

Motion in Space. MATH 311, Calculus III. J. Robert Buchanan. Fall Department of Mathematics. J. Robert Buchanan Motion in Space

Motion in Space. MATH 311, Calculus III. J. Robert Buchanan. Fall Department of Mathematics. J. Robert Buchanan Motion in Space Motion in Space MATH 311, Calculus III J. Robert Buchanan Department of Mathematics Fall 2011 Background Suppose the position vector of a moving object is given by r(t) = f (t), g(t), h(t), Background

More information

Rigid Body Dynamics and Beyond

Rigid Body Dynamics and Beyond Rigid Body Dynamics and Beyond 1 Rigid Bodies 3 A rigid body Collection of particles Distance between any two particles is always constant What types of motions preserve these constraints? Translation,

More information

Rotational & Rigid-Body Mechanics. Lectures 3+4

Rotational & Rigid-Body Mechanics. Lectures 3+4 Rotational & Rigid-Body Mechanics Lectures 3+4 Rotational Motion So far: point objects moving through a trajectory. Next: moving actual dimensional objects and rotating them. 2 Circular Motion - Definitions

More information

Artificial 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 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 information

Manipulator Dynamics 2. Instructor: Jacob Rosen Advanced Robotic - MAE 263D - Department of Mechanical & Aerospace Engineering - UCLA

Manipulator Dynamics 2. Instructor: Jacob Rosen Advanced Robotic - MAE 263D - Department of Mechanical & Aerospace Engineering - UCLA Manipulator Dynamics 2 Forward Dynamics Problem Given: Joint torques and links geometry, mass, inertia, friction Compute: Angular acceleration of the links (solve differential equations) Solution Dynamic

More information

Rigid body simulation

Rigid body simulation Rgd bod smulaton Rgd bod smulaton Once we consder an object wth spacal etent, partcle sstem smulaton s no longer suffcent Problems Problems Unconstraned sstem rotatonal moton torques and angular momentum

More information

CP1 REVISION LECTURE 3 INTRODUCTION TO CLASSICAL MECHANICS. Prof. N. Harnew University of Oxford TT 2017

CP1 REVISION LECTURE 3 INTRODUCTION TO CLASSICAL MECHANICS. Prof. N. Harnew University of Oxford TT 2017 CP1 REVISION LECTURE 3 INTRODUCTION TO CLASSICAL MECHANICS Prof. N. Harnew University of Oxford TT 2017 1 OUTLINE : CP1 REVISION LECTURE 3 : INTRODUCTION TO CLASSICAL MECHANICS 1. Angular velocity and

More information

Rotational Motion. Chapter 4. P. J. Grandinetti. Sep. 1, Chem P. J. Grandinetti (Chem. 4300) Rotational Motion Sep.

Rotational Motion. Chapter 4. P. J. Grandinetti. Sep. 1, Chem P. J. Grandinetti (Chem. 4300) Rotational Motion Sep. Rotational Motion Chapter 4 P. J. Grandinetti Chem. 4300 Sep. 1, 2017 P. J. Grandinetti (Chem. 4300) Rotational Motion Sep. 1, 2017 1 / 76 Angular Momentum The angular momentum of a particle with respect

More information

1/30. Rigid Body Rotations. Dave Frank

1/30. Rigid Body Rotations. Dave Frank . 1/3 Rigid Body Rotations Dave Frank A Point Particle and Fundamental Quantities z 2/3 m v ω r y x Angular Velocity v = dr dt = ω r Kinetic Energy K = 1 2 mv2 Momentum p = mv Rigid Bodies We treat a rigid

More information

Phys 7221 Homework # 8

Phys 7221 Homework # 8 Phys 71 Homework # 8 Gabriela González November 15, 6 Derivation 5-6: Torque free symmetric top In a torque free, symmetric top, with I x = I y = I, the angular velocity vector ω in body coordinates with

More information

General Physics I. Lecture 10: Rolling Motion and Angular Momentum.

General Physics I. Lecture 10: Rolling Motion and Angular Momentum. General Physics I Lecture 10: Rolling Motion and Angular Momentum Prof. WAN, Xin (万歆) 万歆 ) xinwan@zju.edu.cn http://zimp.zju.edu.cn/~xinwan/ Outline Rolling motion of a rigid object: center-of-mass motion

More information

6. 3D Kinematics DE2-EA 2.1: M4DE. Dr Connor Myant

6. 3D Kinematics DE2-EA 2.1: M4DE. Dr Connor Myant DE2-EA 2.1: M4DE Dr Connor Myant 6. 3D Kinematics Comments and corrections to connor.myant@imperial.ac.uk Lecture resources may be found on Blackboard and at http://connormyant.com Contents Three-Dimensional

More information

Game Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost

Game 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

Part 8: Rigid Body Dynamics

Part 8: Rigid Body Dynamics Document that contains homework problems. Comment out the solutions when printing off for students. Part 8: Rigid Body Dynamics Problem 1. Inertia review Find the moment of inertia for a thin uniform rod

More information

Multibody simulation

Multibody 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 information

Closed-Form Solution Of Absolute Orientation Using Unit Quaternions

Closed-Form Solution Of Absolute Orientation Using Unit Quaternions Closed-Form Solution Of Absolute Orientation Using Unit Berthold K. P. Horn Department of Computer and Information Sciences November 11, 2004 Outline 1 Introduction 2 3 The Problem Given: two sets of corresponding

More information

A DARK GREY P O N T, with a Switch Tail, and a small Star on the Forehead. Any

A DARK GREY P O N T, with a Switch Tail, and a small Star on the Forehead. Any Y Y Y X X «/ YY Y Y ««Y x ) & \ & & } # Y \#$& / Y Y X» \\ / X X X x & Y Y X «q «z \x» = q Y # % \ & [ & Z \ & { + % ) / / «q zy» / & / / / & x x X / % % ) Y x X Y $ Z % Y Y x x } / % «] «] # z» & Y X»

More information

Exercise 1: Inertia moment of a simple pendulum

Exercise 1: Inertia moment of a simple pendulum Exercise : Inertia moment of a simple pendulum A simple pendulum is represented in Figure. Three reference frames are introduced: R is the fixed/inertial RF, with origin in the rotation center and i along

More information

PSE Game Physics. Session (6) Angular momentum, microcollisions, damping. Oliver Meister, Roland Wittmann

PSE Game Physics. Session (6) Angular momentum, microcollisions, damping. Oliver Meister, Roland Wittmann PSE Game Physics Session (6) Angular momentum, microcollisions, damping Oliver Meister, Roland Wittmann 23.05.2014 Session (6)Angular momentum, microcollisions, damping, 23.05.2014 1 Outline Angular momentum

More information

Chapter 6: Momentum Analysis

Chapter 6: Momentum Analysis 6-1 Introduction 6-2Newton s Law and Conservation of Momentum 6-3 Choosing a Control Volume 6-4 Forces Acting on a Control Volume 6-5Linear Momentum Equation 6-6 Angular Momentum 6-7 The Second Law of

More information

In this section of notes, we look at the calculation of forces and torques for a manipulator in two settings:

In this section of notes, we look at the calculation of forces and torques for a manipulator in two settings: Introduction Up to this point we have considered only the kinematics of a manipulator. That is, only the specification of motion without regard to the forces and torques required to cause motion In this

More information

AB-267 DYNAMICS & CONTROL OF FLEXIBLE AIRCRAFT

AB-267 DYNAMICS & CONTROL OF FLEXIBLE AIRCRAFT FLÁIO SILESTRE DYNAMICS & CONTROL OF FLEXIBLE AIRCRAFT LECTURE NOTES LAGRANGIAN MECHANICS APPLIED TO RIGID-BODY DYNAMICS IMAGE CREDITS: BOEING FLÁIO SILESTRE Introduction Lagrangian Mechanics shall be

More information

Dynamics. 1 Copyright c 2015 Roderic Grupen

Dynamics. 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 information

SOLUTIONS, PROBLEM SET 11

SOLUTIONS, PROBLEM SET 11 SOLUTIONS, PROBLEM SET 11 1 In this problem we investigate the Lagrangian formulation of dynamics in a rotating frame. Consider a frame of reference which we will consider to be inertial. Suppose that

More information

Physics 106b/196b Problem Set 9 Due Jan 19, 2007

Physics 106b/196b Problem Set 9 Due Jan 19, 2007 Physics 06b/96b Problem Set 9 Due Jan 9, 2007 Version 3: January 8, 2007 This problem set focuses on dynamics in rotating coordinate systems (Section 5.2), with some additional early material on dynamics

More information

Classical Mechanics. Luis Anchordoqui

Classical Mechanics. Luis Anchordoqui 1 Rigid Body Motion Inertia Tensor Rotational Kinetic Energy Principal Axes of Rotation Steiner s Theorem Euler s Equations for a Rigid Body Eulerian Angles Review of Fundamental Equations 2 Rigid body

More information

Chapter 6: Momentum Analysis of Flow Systems

Chapter 6: Momentum Analysis of Flow Systems Chapter 6: Momentum Analysis of Flow Systems Introduction Fluid flow problems can be analyzed using one of three basic approaches: differential, experimental, and integral (or control volume). In Chap.

More information

f x f y or else = m a y

f x f y or else = m a y Rigid Body Dynamics If the aim of kinematics is to describe the body motion, the aim of dynamics is to explain it; the history of mechanics shows that the passage from description to explanation requires

More information

CIS 4930/6930: Principles of Cyber-Physical Systems

CIS 4930/6930: Principles of Cyber-Physical Systems CIS 4930/6930: Principles of Cyber-Physical Systems Chapter 2: Continuous Dynamics Hao Zheng Department of Computer Science and Engineering University of South Florida H. Zheng (CSE USF) CIS 4930/6930:

More information

Spacecraft Dynamics and Control

Spacecraft Dynamics and Control Spacecraft Dynamics and Control Matthew M. Peet Arizona State University Lecture 16: Euler s Equations Attitude Dynamics In this Lecture we will cover: The Problem of Attitude Stabilization Actuators Newton

More information

CS 480/680: GAME ENGINE PROGRAMMING

CS 480/680: GAME ENGINE PROGRAMMING CS 480/680: GAME ENGINE PROGRAMMING PHYSICS 2/14/2013 Santiago Ontañón santi@cs.drexel.edu https://www.cs.drexel.edu/~santi/teaching/2013/cs480-680/intro.html Outline Student Presentations Basics on Rigid

More information

Leonhard Euler ( September 1783)

Leonhard Euler ( September 1783) LEONHARD EULER (1707 - SEPTEMBER 1783) Leonhard Euler (1707 - September 1783) BEYOND EQUATIONS Leonhard Euler was born in Basle, Switzerland; he was in fact a born mathematician, who went on to become

More information

Lecture II: Rigid-Body Physics

Lecture II: Rigid-Body Physics Rigid-Body Motion Previously: Point dimensionless objects moving through a trajectory. Today: Objects with dimensions, moving as one piece. 2 Rigid-Body Kinematics Objects as sets of points. Relative distances

More information

Final Exam. June 10, 2008, 1:00pm

Final Exam. June 10, 2008, 1:00pm PHYSICS 101: Fundamentals of Physics Final Exam Final Exam Name TA/ Section # June 10, 2008, 1:00pm Recitation Time You have 2 hour to complete the exam. Please answer all questions clearly and completely,

More information

Angular Momentum. Objectives CONSERVATION OF ANGULAR MOMENTUM

Angular Momentum. Objectives CONSERVATION OF ANGULAR MOMENTUM Angular Momentum CONSERVATION OF ANGULAR MOMENTUM Objectives Calculate the angular momentum vector for a moving particle Calculate the angular momentum vector for a rotating rigid object where angular

More information

α = p = m v L = I ω Review: Torque Physics 201, Lecture 21 Review: Rotational Dynamics a = Στ = I α

α = p = m v L = I ω Review: Torque Physics 201, Lecture 21 Review: Rotational Dynamics a = Στ = I α Physics 1, Lecture 1 Today s Topics q Static Equilibrium of Rigid Objects(Ch. 1.1-3) Review: Rotational and Translational Motion Conditions for Translational and Rotational Equilibrium Demos and Exercises

More information

Rigid Body Rotation. Speaker: Xiaolei Chen Advisor: Prof. Xiaolin Li. Department of Applied Mathematics and Statistics Stony Brook University (SUNY)

Rigid Body Rotation. Speaker: Xiaolei Chen Advisor: Prof. Xiaolin Li. Department of Applied Mathematics and Statistics Stony Brook University (SUNY) Rigid Body Rotation Speaker: Xiaolei Chen Advisor: Prof. Xiaolin Li Department of Applied Mathematics and Statistics Stony Brook University (SUNY) Content Introduction Angular Velocity Angular Momentum

More information

The Calculus of Vec- tors

The Calculus of Vec- tors Physics 2460 Electricity and Magnetism I, Fall 2007, Lecture 3 1 The Calculus of Vec- Summary: tors 1. Calculus of Vectors: Limits and Derivatives 2. Parametric representation of Curves r(t) = [x(t), y(t),

More information

9 Kinetics of 3D rigid bodies - rotating frames

9 Kinetics of 3D rigid bodies - rotating frames 9 Kinetics of 3D rigid bodies - rotating frames 9. Consider the two gears depicted in the figure. The gear B of radius R B is fixed to the ground, while the gear A of mass m A and radius R A turns freely

More information

Symmetries 2 - Rotations in Space

Symmetries 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 information

Dynamics. Basilio Bona. Semester 1, DAUIN Politecnico di Torino. B. Bona (DAUIN) Dynamics Semester 1, / 18

Dynamics. 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 information

CS-184: Computer Graphics

CS-184: Computer Graphics CS-184: 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 information

particle p = m v F ext = d P = M d v cm dt

particle 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 information

Rotational Motion and Torque

Rotational Motion and Torque Rotational Motion and Torque Introduction to Angular Quantities Sections 8- to 8-2 Introduction Rotational motion deals with spinning objects, or objects rotating around some point. Rotational motion is

More information

z F 3 = = = m 1 F 1 m 2 F 2 m 3 - Linear Momentum dp dt F net = d P net = d p 1 dt d p n dt - Conservation of Linear Momentum Δ P = 0

z F 3 = = = m 1 F 1 m 2 F 2 m 3 - Linear Momentum dp dt F net = d P net = d p 1 dt d p n dt - Conservation of Linear Momentum Δ P = 0 F 1 m 2 F 2 x m 1 O z F 3 m 3 y Ma com = F net F F F net, x net, y net, z = = = Ma Ma Ma com, x com, y com, z p = mv - Linear Momentum F net = dp dt F net = d P dt = d p 1 dt +...+ d p n dt Δ P = 0 - Conservation

More information

Physics 201, Lecture 21

Physics 201, Lecture 21 Physics 201, Lecture 21 Today s Topics q Static Equilibrium of Rigid Objects(Ch. 12.1-3) Review: Rotational and Translational Motion Conditions for Translational and Rotational Equilibrium Demos and Exercises

More information

Rigid bodies - general theory

Rigid bodies - general theory Rigid bodies - general theory Kinetic Energy: based on FW-26 Consider a system on N particles with all their relative separations fixed: it has 3 translational and 3 rotational degrees of freedom. Motion

More information

Dynamics of Open Chains

Dynamics of Open Chains Chapter 9 Dynamics of Open Chains According to Newton s second law of motion, any change in the velocity of a rigid body is caused by external forces and torques In this chapter we study once again the

More information

Thursday Simulation & Unity

Thursday Simulation & Unity Rigid Bodies Simulation Homework Build a particle system based either on F=ma or procedural simulation Examples: Smoke, Fire, Water, Wind, Leaves, Cloth, Magnets, Flocks, Fish, Insects, Crowds, etc. Simulate

More information

PLANAR KINETIC EQUATIONS OF MOTION (Section 17.2)

PLANAR KINETIC EQUATIONS OF MOTION (Section 17.2) PLANAR KINETIC EQUATIONS OF MOTION (Section 17.2) We will limit our study of planar kinetics to rigid bodies that are symmetric with respect to a fixed reference plane. As discussed in Chapter 16, when

More information

Autonomous Underwater Vehicles: Equations of Motion

Autonomous Underwater Vehicles: Equations of Motion Autonomous Underwater Vehicles: Equations of Motion Monique Chyba - November 18, 2015 Departments of Mathematics, University of Hawai i at Mānoa Elective in Robotics 2015/2016 - Control of Unmanned Vehicles

More information

MECH 5312 Solid Mechanics II. Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso

MECH 5312 Solid Mechanics II. Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso MECH 5312 Solid Mechanics II Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso Table of Contents Preliminary Math Concept of Stress Stress Components Equilibrium

More information

Handout 7: Torque, angular momentum, rotational kinetic energy and rolling motion. Torque and angular momentum

Handout 7: Torque, angular momentum, rotational kinetic energy and rolling motion. Torque and angular momentum Handout 7: Torque, angular momentum, rotational kinetic energy and rolling motion Torque and angular momentum In Figure, in order to turn a rod about a fixed hinge at one end, a force F is applied at a

More information

Chapter 4 The Equations of Motion

Chapter 4 The Equations of Motion Chapter 4 The Equations of Motion Flight Mechanics and Control AEM 4303 Bérénice Mettler University of Minnesota Feb. 20-27, 2013 (v. 2/26/13) Bérénice Mettler (University of Minnesota) Chapter 4 The Equations

More information

Kinematics (special case) Dynamics gravity, tension, elastic, normal, friction. Energy: kinetic, potential gravity, spring + work (friction)

Kinematics (special case) Dynamics gravity, tension, elastic, normal, friction. Energy: kinetic, potential gravity, spring + work (friction) Kinematics (special case) a = constant 1D motion 2D projectile Uniform circular Dynamics gravity, tension, elastic, normal, friction Motion with a = constant Newton s Laws F = m a F 12 = F 21 Time & Position

More information

If the symmetry axes of a uniform symmetric body coincide with the coordinate axes, the products of inertia (Ixy etc.

If the symmetry axes of a uniform symmetric body coincide with the coordinate axes, the products of inertia (Ixy etc. Prof. O. B. Wright, Autumn 007 Mechanics Lecture 9 More on rigid bodies, coupled vibrations Principal axes of the inertia tensor If the symmetry axes of a uniform symmetric body coincide with the coordinate

More information

In this section, mathematical description of the motion of fluid elements moving in a flow field is

In this section, mathematical description of the motion of fluid elements moving in a flow field is Jun. 05, 015 Chapter 6. Differential Analysis of Fluid Flow 6.1 Fluid Element Kinematics In this section, mathematical description of the motion of fluid elements moving in a flow field is given. A small

More information

Chapter 8 Rotational Equilibrium and Rotational Dynamics Force vs. Torque Forces cause accelerations Torques cause angular accelerations Force and

Chapter 8 Rotational Equilibrium and Rotational Dynamics Force vs. Torque Forces cause accelerations Torques cause angular accelerations Force and Chapter 8 Rotational Equilibrium and Rotational Dynamics Force vs. Torque Forces cause accelerations Torques cause angular accelerations Force and torque are related Torque The door is free to rotate about

More information

Distance travelled time taken and if the particle is a distance s(t) along the x-axis, then its instantaneous speed is:

Distance travelled time taken and if the particle is a distance s(t) along the x-axis, 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 information

Game Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost

Game Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Essential physics for game developers Introduction The primary issues Let s move virtual objects Kinematics: description

More information

Projects. Assignments. Plan. Motivational Films. Motivational Film 1/6/09. Dynamics II. Forces, Collisions and Impulse

Projects. Assignments. Plan. Motivational Films. Motivational Film 1/6/09. Dynamics II. Forces, Collisions and Impulse Projects Dynamics II Forces, Collisions and Impulse Presentations: Dates: Week 0: Mon, Feb 6 Week 0: Wed, Feb 8 Finals Week: Tues, Feb 4 (:30 :30) room TBA 5 minutes / presentation Schedule now on Web

More information

Translational and Rotational Dynamics!

Translational 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 information

ω = 0 a = 0 = α P = constant L = constant dt = 0 = d Equilibrium when: τ i = 0 τ net τ i Static Equilibrium when: F z = 0 F net = F i = ma = d P

ω = 0 a = 0 = α P = constant L = constant dt = 0 = d Equilibrium when: τ i = 0 τ net τ i Static Equilibrium when: F z = 0 F net = F i = ma = d P Equilibrium when: F net = F i τ net = τ i a = 0 = α dp = 0 = d L = ma = d P = 0 = I α = d L = 0 P = constant L = constant F x = 0 τ i = 0 F y = 0 F z = 0 Static Equilibrium when: P = 0 L = 0 v com = 0

More information

Lecture 10 - Moment of Inertia

Lecture 10 - Moment of Inertia Lecture 10 - oment of Inertia A Puzzle... Question For any object, there are typically many ways to calculate the moment of inertia I = r 2 dm, usually by doing the integration by considering different

More information

Introduction to Haptic Systems

Introduction to Haptic Systems Introduction to Haptic Systems Félix Monasterio-Huelin & Álvaro Gutiérrez & Blanca Larraga October 8, 2018 Contents Contents 1 List of Figures 1 1 Introduction 2 2 DC Motor 3 3 1 DOF DC motor model with

More information

Kinetic Energy of Rolling

Kinetic Energy of Rolling Kinetic Energy of Rolling A solid disk and a hoop (with the same mass and radius) are released from rest and roll down a ramp from a height h. Which one is moving faster at the bottom of the ramp? A. they

More information

Chapter 8. Rotational Equilibrium and Rotational Dynamics

Chapter 8. Rotational Equilibrium and Rotational Dynamics Chapter 8 Rotational Equilibrium and Rotational Dynamics Force vs. Torque Forces cause accelerations Torques cause angular accelerations Force and torque are related Torque The door is free to rotate about

More information

Game Physics: Basic Concepts

Game Physics: Basic Concepts CMSC 498M: Chapter 8a Game Physics Reading: Game Physics, by David H Eberly, 2004 Physics for Game Developers, by David M Bourg, 2002 Overview: Basic physical quantities: Mass, center of mass, moment of

More information

Spatial Vector Algebra

Spatial Vector Algebra A Short Course on The Easy Way to do Rigid Body Dynamics Roy Featherstone Dept. Inormation Engineering, RSISE The Australian National University Spatial vector algebra is a concise vector notation or describing

More information

Applications of Eigenvalues & Eigenvectors

Applications of Eigenvalues & Eigenvectors Applications of Eigenvalues & Eigenvectors Louie L. Yaw Walla Walla University Engineering Department For Linear Algebra Class November 17, 214 Outline 1 The eigenvalue/eigenvector problem 2 Principal

More information

3-D Kinetics of Rigid Bodies

3-D Kinetics of Rigid Bodies 3-D Kinetics of Rigid Bodies Angular Momentum Generalized Newton s second law for the motion for a 3-D mass system Moment eqn for 3-D motion will be different than that obtained for plane motion Consider

More information

Chapter 9. Rotational Dynamics

Chapter 9. Rotational Dynamics Chapter 9 Rotational Dynamics 9.1 The Action of Forces and Torques on Rigid Objects In pure translational motion, all points on an object travel on parallel paths. The most general motion is a combination

More information

Sec. 1.1: Basics of Vectors

Sec. 1.1: Basics of Vectors Sec. 1.1: Basics of Vectors Notation for Euclidean space R n : all points (x 1, x 2,..., x n ) in n-dimensional space. Examples: 1. R 1 : all points on the real number line. 2. R 2 : all points (x 1, x

More information

Momentum. The way to catch a knuckleball is to wait until it stops rolling and then pick it up. -Bob Uecker

Momentum. The way to catch a knuckleball is to wait until it stops rolling and then pick it up. -Bob Uecker Chapter 11 -, Chapter 11 -, Angular The way to catch a knuckleball is to wait until it stops rolling and then pick it up. -Bob Uecker David J. Starling Penn State Hazleton PHYS 211 Chapter 11 -, motion

More information

Chapter 11. Angular Momentum

Chapter 11. Angular Momentum Chapter 11 Angular Momentum Angular Momentum Angular momentum plays a key role in rotational dynamics. There is a principle of conservation of angular momentum. In analogy to the principle of conservation

More information

Lecture 8 Analyzing the diffusion weighted signal. Room CSB 272 this week! Please install AFNI

Lecture 8 Analyzing the diffusion weighted signal. Room CSB 272 this week! Please install AFNI Lecture 8 Analyzing the diffusion weighted signal Room CSB 272 this week! Please install AFNI http://afni.nimh.nih.gov/afni/ Next lecture, DTI For this lecture, think in terms of a single voxel We re still

More information

Lecture AC-1. Aircraft Dynamics. Copy right 2003 by Jon at h an H ow

Lecture AC-1. Aircraft Dynamics. Copy right 2003 by Jon at h an H ow Lecture AC-1 Aircraft Dynamics Copy right 23 by Jon at h an H ow 1 Spring 23 16.61 AC 1 2 Aircraft Dynamics First note that it is possible to develop a very good approximation of a key motion of an aircraft

More information

13, Applications of molecular symmetry and group theory

13, Applications of molecular symmetry and group theory Subject Paper No and Title Module No and Title Module Tag Chemistry 13, Applications of molecular symmetry and group theory 27, Group theory and vibrational spectroscopy: Part-IV(Selection rules for IR

More information

General Relativity ASTR 2110 Sarazin. Einstein s Equation

General Relativity ASTR 2110 Sarazin. Einstein s Equation General Relativity ASTR 2110 Sarazin Einstein s Equation Curvature of Spacetime 1. Principle of Equvalence: gravity acceleration locally 2. Acceleration curved path in spacetime In gravitational field,

More information

OCN/ATM/ESS 587. The wind-driven ocean circulation. Friction and stress. The Ekman layer, top and bottom. Ekman pumping, Ekman suction

OCN/ATM/ESS 587. The wind-driven ocean circulation. Friction and stress. The Ekman layer, top and bottom. Ekman pumping, Ekman suction OCN/ATM/ESS 587 The wind-driven ocean circulation. Friction and stress The Ekman layer, top and bottom Ekman pumping, Ekman suction Westward intensification The wind-driven ocean. The major ocean gyres

More information

Flips are Fun: A Study in the Conservation of Angular Momentum of a Non-Rigid Human Body

Flips are Fun: A Study in the Conservation of Angular Momentum of a Non-Rigid Human Body Flips are Fun: A Study in the Conservation of Angular Momentum of a Non-Rigid Human Body Benjamin Jenkins Physics Department, The College of Wooster, Wooster, Ohio 44691, USA May 10, 2018 Abstract In this

More information

Application of Forces. Chapter Eight. Torque. Force vs. Torque. Torque, cont. Direction of Torque 4/7/2015

Application of Forces. Chapter Eight. Torque. Force vs. Torque. Torque, cont. Direction of Torque 4/7/2015 Raymond A. Serway Chris Vuille Chapter Eight Rotational Equilibrium and Rotational Dynamics Application of Forces The point of application of a force is important This was ignored in treating objects as

More information

Rotational Motion, Torque, Angular Acceleration, and Moment of Inertia. 8.01t Nov 3, 2004

Rotational Motion, Torque, Angular Acceleration, and Moment of Inertia. 8.01t Nov 3, 2004 Rotational Motion, Torque, Angular Acceleration, and Moment of Inertia 8.01t Nov 3, 2004 Rotation and Translation of Rigid Body Motion of a thrown object Translational Motion of the Center of Mass Total

More information

16.3 Conservative Vector Fields

16.3 Conservative Vector Fields 16.3 Conservative Vector Fields Lukas Geyer Montana State University M273, Fall 2011 Lukas Geyer (MSU) 16.3 Conservative Vector Fields M273, Fall 2011 1 / 23 Fundamental Theorem for Conservative Vector

More information

LOWELL WEEKLY JOURNAL

LOWELL WEEKLY JOURNAL Y -» $ 5 Y 7 Y Y -Y- Q x Q» 75»»/ q } # ]»\ - - $ { Q» / X x»»- 3 q $ 9 ) Y q - 5 5 3 3 3 7 Q q - - Q _»»/Q Y - 9 - - - )- [ X 7» -» - )»? / /? Q Y»» # X Q» - -?» Q ) Q \ Q - - - 3? 7» -? #»»» 7 - / Q

More information

Lecture 17 - Gyroscopes

Lecture 17 - Gyroscopes Lecture 17 - Gyroscopes A Puzzle... We have seen throughout class that the center of mass is a very powerful tool for evaluating systems. However, don t let yourself get carried away with how useful it

More information

Classical Mechanics Lecture 15

Classical Mechanics Lecture 15 Classical Mechanics Lecture 5 Today s Concepts: a) Parallel Axis Theorem b) Torque & Angular Acceleration Mechanics Lecture 5, Slide Unit 4 Main Points Mechanics Lecture 4, Slide Unit 4 Main Points Mechanics

More information

M E 320 Professor John M. Cimbala Lecture 10

M E 320 Professor John M. Cimbala Lecture 10 M E 320 Professor John M. Cimbala Lecture 10 Today, we will: Finish our example problem rates of motion and deformation of fluid particles Discuss the Reynolds Transport Theorem (RTT) Show how the RTT

More information

Forces of Rolling. 1) Ifobjectisrollingwith a com =0 (i.e.no netforces), then v com =ωr = constant (smooth roll)

Forces of Rolling. 1) Ifobjectisrollingwith a com =0 (i.e.no netforces), then v com =ωr = constant (smooth roll) Physics 2101 Section 3 March 12 rd : Ch. 10 Announcements: Mid-grades posted in PAW Quiz today I will be at the March APS meeting the week of 15-19 th. Prof. Rich Kurtz will help me. Class Website: http://www.phys.lsu.edu/classes/spring2010/phys2101-3/

More information

Robot Control Basics CS 685

Robot 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 information

Work - kinetic energy theorem for rotational motion *

Work - kinetic energy theorem for rotational motion * OpenStax-CNX module: m14307 1 Work - kinetic energy theorem for rotational motion * Sunil Kumar Singh This work is produced by OpenStax-CNX and licensed under the Creative Commons Attribution License 2.0

More information

Properties of surfaces II: Second moment of area

Properties of surfaces II: Second moment of area Properties of surfaces II: Second moment of area Just as we have discussing first moment of an area and its relation with problems in mechanics, we will now describe second moment and product of area of

More information

PHYSICS 44 {CLASSICAL MECHANICS 2003 SUMMER FINAL EXAM { SOLUTIONS

PHYSICS 44 {CLASSICAL MECHANICS 2003 SUMMER FINAL EXAM { SOLUTIONS PHYSICS 44 {CLASSICAL MECHANICS 3 SUMMER FINAL EXAM { SOLUTIONS Problem 1 Consider the following Atwood machine comosed of the two masses m 1 and m and attached to a stand with the hel of a sring of strength

More information

ELECTROMAGNETIC WAVES

ELECTROMAGNETIC WAVES Physics 4D ELECTROMAGNETIC WAVE Hans P. Paar 26 January 2006 i Chapter 1 Vector Calculus 1.1 Introduction Vector calculus is a branch of mathematics that allows differentiation and integration of (scalar)

More information

Classical Mechanics Comprehensive Exam Solution

Classical Mechanics Comprehensive Exam Solution Classical Mechanics Comprehensive Exam Solution January 31, 011, 1:00 pm 5:pm Solve the following six problems. In the following problems, e x, e y, and e z are unit vectors in the x, y, and z directions,

More information

Chapters 10 & 11: Rotational Dynamics Thursday March 8 th

Chapters 10 & 11: Rotational Dynamics Thursday March 8 th Chapters 10 & 11: Rotational Dynamics Thursday March 8 th Review of rotational kinematics equations Review and more on rotational inertia Rolling motion as rotation and translation Rotational kinetic energy

More information

Problem Set 2 Due Tuesday, September 27, ; p : 0. (b) Construct a representation using five d orbitals that sit on the origin as a basis: 1

Problem Set 2 Due Tuesday, September 27, ; p : 0. (b) Construct a representation using five d orbitals that sit on the origin as a basis: 1 Problem Set 2 Due Tuesday, September 27, 211 Problems from Carter: Chapter 2: 2a-d,g,h,j 2.6, 2.9; Chapter 3: 1a-d,f,g 3.3, 3.6, 3.7 Additional problems: (1) Consider the D 4 point group and use a coordinate

More information

Two-Dimensional Rotational Kinematics

Two-Dimensional Rotational Kinematics Two-Dimensional Rotational Kinematics Rigid Bodies A rigid body is an extended object in which the distance between any two points in the object is constant in time. Springs or human bodies are non-rigid

More information