Constrained dynamics
|
|
- Hortense Cameron
- 5 years ago
- Views:
Transcription
1 Constrained dynamics
2 Simple particle system In principle, you can make just about anything out of spring systems In practice, you can make just about anything as long as it s jello
3 Hard constraints
4 Constraint force Single implicit constraint Multiple implicit constraint Parametric constraint Implementation
5 A simple example A bead on a wire The bead can slide freely along the wire, but cannot come off it no matter how hard you pull it. How do we simulate the motion of the bead when arbitrary forces applied to it?
6 Penalty constraints Why not use a spring to hold the bead on the wire? Problems: weak springs won t do the job strong springs give you stiff systems
7 First order world N f T In this world, f = mv ff ˆf N What is the legal velocity? What is the legal force? Add constraint force ˆf to cancel the illegal part of f f = f + ˆf ˆf = N f N N N
8 The real world N ˆf? v f a In the real world, f = ma What is the legal acceleration? depends on both N and the faster you re going, the faster you have to turn v Compute ˆf such that f only generates legal acceleration f = f + ˆf
9 Constraint force Need to compute constraint forces that cancel the illegal applied forces Which means we need to know what legal acceleration is
10 Constraint force Single implicit constraint Multiple implicit constraint Parametric constraint Implementation
11 Constraints Implicit: x C(x) = x r =0 θ r Parametric: x = r [ cos θ sin θ ]
12 Legal acceleration C =0 Ċ =0 C =0 What is the legal position? C(x) = 1 2 x x 1 2 =0 What is the legal velocity? Ċ(x) =x ẋ =0 What is the legal acceleration? C(x) =ẍ x + ẋ ẋ =0
13 Legal conditions C =0 Ċ =0 C =0 If we start with legal position and velocity C(x) =0 Ċ(x) =0 We need only ensure the legal acceleration C(x) =0
14 Constraint force Use the legal condition to compute the constraint force Rewrite the legal condition in a general form C(x) = Ċ x ẋ + C x ẍ =0 C x : constraint gradient Substitute ẍ with ẍ = f + ˆf m C = Ċ x ẋ + C x f + ˆf m =0 C x ˆf = C x f m Ċ xẋ
15 Constraint force C x ˆf = C x f m Ċ xẋ How many variables do we have? Need one more condition to solve the constraint force
16 Virtual work Constraint force is passive - no energy gain or loss Kinetic energy of the system: T = 1 2 mẋ ẋ Virtual work done by f and ˆf : T = ẋ mẍ = ẋ f + ẋ ˆf Make sure ˆf does no work for every legal velocity: ẋ ˆf =0, ẋ C x ẋ =0 Ċ(x) = C x ẋ =0
17 Constraint force ˆf ẋ ˆf =0, ẋ C x ẋ =0 must point in the direction of C x ˆf = λ C x Substituting for ˆf in C x ˆf = C x f m Ċ xẋ λ = C x f m Ċ x ẋ C x C x
18 Two conditions Legal acceleration C x ˆf = C x f m Ċ xẋ Principle of virtual work ˆf = λ C x
19 The Bead Example N v f C(x) = 1 2 x x 1 2 =0 ˆf? ˆf = λ C x λ = C x f m Ċ x ẋ C x C x
20 Feedback In principle, ensuring legal acceleration can keep the particle exactly on the circle In practice, two problems cause the particle to drift numerical errors can accumulate when the ODE is not solved exactly constraints might not be met initially
21 Feedback A feedback term handles both problems: C = k s C k d Ċ instead of C =0
22 Constraint force Single implicit constraint Multiple implicit constraint Parametric constraint Implementation
23 Tinkertoys Now we know how to simulate a bead on a wire Apply the simple idea, we can create a constrained particle system
24 Constrained particles Particles: each particle represents a point in the phase space Forces: each force affects the acceleration of certain particles Constraints: Each is a function C i (x 1, x 2,...) Legal state: C i (x 1, x 2,...)=0, i Constraint force: linear combination of constraint gradients C i x, i
25 Constraint gradients Normal of C 1 : Normal of C 2 : C 1 x C 2 x Legal states: the intersection of two planes Normal of the legal states: +
26 Implicit constraint Each constraint is represented by an implicit function What does the normal of a hypersurface mean? The direction where the particle is not allowed to move What does the intersection of the hypersurface represent? Legal state The constraint force lies in the space spanned by the constraint normals
27 Particle system notations General 3D case q Q W x 1 x 2 f 1 f 2 1 m 1 I 1 m 2 I q: 3n long position vector.... Q: 3n long force vector x n f n 1 m n I W : 3n 3n inverse mass matrix
28 Constraint notations Additional notations C λ J J C 1 C 2 λ 1 λ 2 C 1 q 1 C 1 q 2 C 1 q 3n C 2 q 1... Ċ1 q 1 Ċ2 q 1 Ċ1 q 2... Ċ1 q 3n C m λ m C m q 3n Ċm q 3n C: m long constraint vector λ: m Lagrangian multipliers J: m 3n Jacobian matrix J: time derivative of Jacobian matrix
29 Constraint equations C = Ċ x ẋ + C x ẍ C = J q + J q C = Ċ x ẋ + C x f + ˆf m =0 C = J q + JW(Q + ˆQ) =0 C x ˆf = C x f m Ċ xẋ JW ˆQ = J q JWQ C x λ C x = C x f m Ċ xẋ JWJ T λ = J q JWQ
30 Multiple force constraints To solve for force constraints, we need to solve for this linear system JWJ T λ = J q JWQ One nice property of positive definite matrix Can you prove that? JWJ T is that it is symmetric Once the linear system has been solved, the vector λ is multiplied by vector ˆQ J T to produce the global constraint force
31 Constraint force Single implicit constraint Multiple implicit constraint Parametric constraint Implementation
32 Constraints x Implicit: C(x) = x r =0 θ r Parametric: x = r [ cos θ sin θ ]
33 Parametric constraints x f x = r [ cos θ sin θ ] θ r Constraint is always met exactly 1 degree of freedom: Solve for θ
34 First order world θ x N In this world, f = mv ẋ = f + ˆf m ẋ = T θ T θ = f + ˆf m T = x θ T T θ = T f m + T ˆf m ˆf = λ C x = λn θ = 1 m T f T T
35 The real world θ x N In the real world, f = ma ẋ = T θ T = x θ ẍ = Ṫ θ + T θ = f + ˆf m T Ṫ θ + T T θ = T f m + T ˆf m θ = T f m T Ṫ θ T T
36 Particle system notations General 2D case q x 1 x 2. u θ 1 θ 2. Q f 1 f 2. m 1 I m 2 I... M q: 2n long position vector u: n long parameter vector Q: 2n long force vector x n θ n f n m n I M: 2n by 2n mass matrix
37 Constraint notations Additional notations J J q 1 q 1 u 1 u n q 2 u 1 q 1 u 1 q 2 u 1 q 1 u n J : 2n by n Jacobian matrix J : time derivative of J q 2n u 1 q 2n u n q 2n u 1 q 2n u n
38 Constraint equations ẍ = Ṫ θ + T θ = f + ˆf m T Ṫ θ + T T θ = T f m θ = T f m T Ṫ θ T T M q = M( J u + Jü) =Q + ˆQ J T M( J u + Jü) =J T Q J T MJü = J T M J u + J T Q where T = x θ where J = q u
39 Parametric vs. implicit Parametric Implicit J T MJü = J T M J u + J T Q JWJ T λ = J q JWQ where J = q u where J = C q Lagrangian dynamics
40 Parametric constraints Advantages: Fewer degrees of freedom Constraints are always met Disadvantages: Hard to formulate constraints Hard to combine constraints
41 Impress your friends The requirement that constraints not add or remove energy is called the Principle of Virtual Work The λ s are called Lagrangain Multipliers The derivative matrix J is called the Jacobian Matrix
42 Constraint force Single implicit constraint Multiple implicit constraint Parametric constraint Implementation
43 How do we implement all this? We have a global matrix equation We want to build a model on the fly Each constraint function knows how to evaluate the function itself and its various derivatives
44 How do we hook this up? We have a basic particle system which main job is to perform derivative evaluations To add constraints to the basic system, we need to modify the data structure of the basic system the deriv eval loop
45 Basic particle system system solver interface solver particles GetDim 6n n time forces Get/Set State Deriv Eval x 1 v 1 x 2 v 2 v v 1 m 1 f 1 2 f 2... m 2... x n v n v n fn m n
46 Constrained system system particles n time x 1 v 1 f 1 m 1 x 2 v 2 f 2 m 2... x n v n f n m n forces F 1 F 2... F p constraints C 1 C 2... C m
47 A Constraint C x p x q v p v q f p f q m p m q code that evaluates p j len apply_fun C Ċ C q Ċ q global structure C Ċ J J q λ ˆQ
48 Jacobian matrix n blocks Each constraint contributes one or more blocks to the Jacobian matrix m blocks C i Sparsity: many empty blocks Modularity: let each constraint compute its own blocks x j x k Constraint and particle indices determine block location
49 Deriv Eval 1. Clear force accumulators x 1 v 1 f 1 m 1 x 2 v 2 f 2 m 2... x n v n f n m n 2. Invoke apply_force functions F 1 F 2... F p 3. Return derivatives to solver [ ] [ ] ẋ v = f v m
50 Modified Deriv Eval loop 1. Clear force accumulators x 1 v 1 f 1 m 1 x 2 v 2 f 2 m 2... x n v n f n m n 2. Invoke apply_force functions F 1 F 2... F p 4. Return derivatives to solver [ ] [ ] ẋ v = f v m 3. Compute and apply constraint forces C 1 C 2... C m
51 Constraint force eval After computing ordinary forces: loop over constraints, assemble global structure call matrix solver to solve for λ, multiply by J T to get constraint force add constraint force to the force accumulator in the corresponding particle
52 What s next?
53 Beyond mass points: equations of motion for real objects Read: Constrained dynamics by Witkin and Baraff
Differential Equations
Pysics-based simulation xi Differential Equations xi+1 xi xi+1 xi + x x Pysics-based simulation xi Wat is a differential equation? Differential equations describe te relation between an unknown function
More informationConstrained motion and generalized coordinates
Constrained motion and generalized coordinates based on FW-13 Often, the motion of particles is restricted by constraints, and we want to: work only with independent degrees of freedom (coordinates) k
More informationDifferential Equations
Differential Equations Overview of differential equation! Initial value problem! Explicit numeric methods! Implicit numeric methods! Modular implementation Physics-based simulation An algorithm that
More informationPhysics 6010, Fall 2016 Constraints and Lagrange Multipliers. Relevant Sections in Text:
Physics 6010, Fall 2016 Constraints and Lagrange Multipliers. Relevant Sections in Text: 1.3 1.6 Constraints Often times we consider dynamical systems which are defined using some kind of restrictions
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 informationM2A2 Problem Sheet 3 - Hamiltonian Mechanics
MA Problem Sheet 3 - Hamiltonian Mechanics. The particle in a cone. A particle slides under gravity, inside a smooth circular cone with a vertical axis, z = k x + y. Write down its Lagrangian in a) Cartesian,
More informationPhysics 106a, Caltech 16 October, Lecture 5: Hamilton s Principle with Constraints. Examples
Physics 106a, Caltech 16 October, 2018 Lecture 5: Hamilton s Principle with Constraints We have been avoiding forces of constraint, because in many cases they are uninteresting, and the constraints can
More informationPhysics 6010, Fall Relevant Sections in Text: Introduction
Physics 6010, Fall 2016 Introduction. Configuration space. Equations of Motion. Velocity Phase Space. Relevant Sections in Text: 1.1 1.4 Introduction This course principally deals with the variational
More informationPontryagin s maximum principle
Pontryagin s maximum principle Emo Todorov Applied Mathematics and Computer Science & Engineering University of Washington Winter 2012 Emo Todorov (UW) AMATH/CSE 579, Winter 2012 Lecture 5 1 / 9 Pontryagin
More information5 Handling Constraints
5 Handling Constraints Engineering design optimization problems are very rarely unconstrained. Moreover, the constraints that appear in these problems are typically nonlinear. This motivates our interest
More informationMATHEMATICAL MODELLING, MECHANICS AND MOD- ELLING MTHA4004Y
UNIVERSITY OF EAST ANGLIA School of Mathematics Main Series UG Examination 2017 18 MATHEMATICAL MODELLING, MECHANICS AND MOD- ELLING MTHA4004Y Time allowed: 2 Hours Attempt QUESTIONS 1 and 2, and ONE other
More informationIn the presence of viscous damping, a more generalized form of the Lagrange s equation of motion can be written as
2 MODELING Once the control target is identified, which includes the state variable to be controlled (ex. speed, position, temperature, flow rate, etc), and once the system drives are identified (ex. force,
More informationUsing Simulink to analyze 2 degrees of freedom system
Using Simulink to analyze 2 degrees of freedom system Nasser M. Abbasi Spring 29 page compiled on June 29, 25 at 4:2pm Abstract A two degrees of freedom system consisting of two masses connected by springs
More informationPhysically Based Modeling: Principles and Practice Differential Equation Basics
Physically Based Modeling: Principles and Practice Differential Equation Basics Andrew Witkin and David Baraff Robotics Institute Carnegie Mellon University Please note: This document is 1997 by Andrew
More informationPhysically Based Modeling Differential Equation Basics
Physically Based Modeling Differential Equation Basics Andrew Witkin and David Baraff Pixar Animation Studios Please note: This document is 2001 by Andrew Witkin and David Baraff. This chapter may be freely
More informationCP1 REVISION LECTURE 3 INTRODUCTION TO CLASSICAL MECHANICS. Prof. N. Harnew University of Oxford TT 2017
CP1 REVISION LECTURE 3 INTRODUCTION TO CLASSICAL MECHANICS Prof. N. Harnew University of Oxford TT 2017 1 OUTLINE : CP1 REVISION LECTURE 3 : INTRODUCTION TO CLASSICAL MECHANICS 1. Angular velocity and
More 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 informationWhat make cloth hard to simulate?
Cloth Simulation What make cloth hard to simulate? Due to the thin and flexible nature of cloth, it produces detailed folds and wrinkles, which in turn can lead to complicated selfcollisions. Cloth is
More informationCoordinates, phase space, constraints
Coordinates, phase space, constraints Sourendu Gupta TIFR, Mumbai, India Classical Mechanics 2011 August 1, 2011 Mechanics of a single particle For the motion of a particle (of constant mass m and position
More informationReview: control, feedback, etc. Today s topic: state-space models of systems; linearization
Plan of the Lecture Review: control, feedback, etc Today s topic: state-space models of systems; linearization Goal: a general framework that encompasses all examples of interest Once we have mastered
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 informationDynamics 4600:203 Homework 09 Due: April 04, 2008 Name:
Dynamics 4600:03 Homework 09 Due: April 04, 008 Name: Please denote your answers clearly, i.e., box in, star, etc., and write neatly. There are no points for small, messy, unreadable work... please use
More informationN mg N Mg N Figure : Forces acting on particle m and inclined plane M. (b) The equations of motion are obtained by applying the momentum principles to
.004 MDEING DNMIS ND NTR I I Spring 00 Solutions for Problem Set 5 Problem. Particle slides down movable inclined plane. The inclined plane of mass M is constrained to move parallel to the -axis, and the
More informationLAGRANGIAN AND HAMILTONIAN
LAGRANGIAN AND HAMILTONIAN A. Constraints and Degrees of Freedom. A constraint is a restriction on the freedom of motion of a system of particles in the form of a condition. The number of independent ways
More informationCoordinates, phase space, constraints
Coordinates, phase space, constraints Sourendu Gupta TIFR, Mumbai, India Classical Mechanics 2012 6 August, 2012 Mechanics of a single particle For the motion of a particle (of constant mass m and position
More informationCS277 - Experimental Haptics Lecture 13. Six-DOF Haptic Rendering I
CS277 - Experimental Haptics Lecture 13 Six-DOF Haptic Rendering I Outline Motivation Direct rendering Proxy-based rendering - Theory - Taxonomy Motivation 3-DOF avatar The Holy Grail? Tool-Mediated Interaction
More informationClassical Mechanics Review (Louisiana State University Qualifier Exam)
Review Louisiana State University Qualifier Exam Jeff Kissel October 22, 2006 A particle of mass m. at rest initially, slides without friction on a wedge of angle θ and and mass M that can move without
More information42. Change of Variables: The Jacobian
. Change of Variables: The Jacobian It is common to change the variable(s) of integration, the main goal being to rewrite a complicated integrand into a simpler equivalent form. However, in doing so, the
More informationLecture 4. Alexey Boyarsky. October 6, 2015
Lecture 4 Alexey Boyarsky October 6, 2015 1 Conservation laws and symmetries 1.1 Ignorable Coordinates During the motion of a mechanical system, the 2s quantities q i and q i, (i = 1, 2,..., s) which specify
More informationESM 3124 Intermediate Dynamics 2012, HW6 Solutions. (1 + f (x) 2 ) We can first write the constraint y = f(x) in the form of a constraint
ESM 314 Intermediate Dynamics 01, HW6 Solutions Roller coaster. A bead of mass m can slide without friction, under the action of gravity, on a smooth rigid wire which has the form y = f(x). (a) Find the
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 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 informationThis module requires you to read a textbook such as Fowles and Cassiday on material relevant to the following topics.
Module M2 Lagrangian Mechanics and Oscillations Prerequisite: Module C1 This module requires you to read a textbook such as Fowles and Cassiday on material relevant to the following topics. Topics: Hamilton
More informationConsider a particle in 1D at position x(t), subject to a force F (x), so that mẍ = F (x). Define the kinetic energy to be.
Chapter 4 Energy and Stability 4.1 Energy in 1D Consider a particle in 1D at position x(t), subject to a force F (x), so that mẍ = F (x). Define the kinetic energy to be T = 1 2 mẋ2 and the potential energy
More informationLecture Outline Chapter 6. Physics, 4 th Edition James S. Walker. Copyright 2010 Pearson Education, Inc.
Lecture Outline Chapter 6 Physics, 4 th Edition James S. Walker Chapter 6 Applications of Newton s Laws Units of Chapter 6 Frictional Forces Strings and Springs Translational Equilibrium Connected Objects
More informationLab #2 - Two Degrees-of-Freedom Oscillator
Lab #2 - Two Degrees-of-Freedom Oscillator Last Updated: March 0, 2007 INTRODUCTION The system illustrated in Figure has two degrees-of-freedom. This means that two is the minimum number of coordinates
More informationDaba Meshesha Gusu and O.Chandra Sekhara Reddy 1
International Journal of Basic and Applied Sciences Vol. 4. No. 1 2015. Pp.22-27 Copyright by CRDEEP. All Rights Reserved. Full Length Research Paper Solutions of Non Linear Ordinary Differential Equations
More informationLecture 6, September 1, 2017
Engineering Mathematics Fall 07 Lecture 6, September, 07 Escape Velocity Suppose we have a planet (or any large near to spherical heavenly body) of radius R and acceleration of gravity at the surface of
More informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
More informationSolutions to old Exam 3 problems
Solutions to old Exam 3 problems Hi students! I am putting this version of my review for the Final exam review here on the web site, place and time to be announced. Enjoy!! Best, Bill Meeks PS. There are
More information(a) Sections 3.7 through (b) Sections 8.1 through 8.3. and Sections 8.5 through 8.6. Problem Set 4 due Monday, 10/14/02
Physics 601 Dr. Dragt Fall 2002 Reading Assignment #4: 1. Dragt (a) Sections 1.5 and 1.6 of Chapter 1 (Introductory Concepts). (b) Notes VI, Hamilton s Equations of Motion (to be found right after the
More informationLecture 10. Example: Friction and Motion
Lecture 10 Goals: Exploit Newton s 3 rd Law in problems with friction Employ Newton s Laws in 2D problems with circular motion Assignment: HW5, (Chapter 7, due 2/24, Wednesday) For Tuesday: Finish reading
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 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 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 informationCHAPTER 1 Systems of Linear Equations
CHAPTER Systems of Linear Equations Section. Introduction to Systems of Linear Equations. Because the equation is in the form a x a y b, it is linear in the variables x and y. 0. Because the equation cannot
More informationPhysics Exam 2009 University of Houston Math Contest. Name: School: There is no penalty for guessing.
Physics Exam 2009 University of Houston Math Contest Name: School: Please read the questions carefully and give a clear indication of your answer on each question. There is no penalty for guessing. Judges
More informationHow to Simulate a Trebuchet Part 1: Lagrange s Equations
Pivot Arm Counterweight Hook Base Sling Projectile How to Simulate a Trebuchet Part 1: Lagrange s Equations The trebuchet has quickly become a favorite project for physics and engineering teachers seeking
More informationSolution Set Two. 1 Problem #1: Projectile Motion Cartesian Coordinates Polar Coordinates... 3
: Solution Set Two Northwestern University, Classical Mechanics Classical Mechanics, Third Ed.- Goldstein October 7, 2015 Contents 1 Problem #1: Projectile Motion. 2 1.1 Cartesian Coordinates....................................
More informationPractice Problems for Final Exam
Math 1280 Spring 2016 Practice Problems for Final Exam Part 2 (Sections 6.6, 6.7, 6.8, and chapter 7) S o l u t i o n s 1. Show that the given system has a nonlinear center at the origin. ẋ = 9y 5y 5,
More information1 Motion of a single particle - Linear momentum, work and energy principle
1 Motion of a single particle - Linear momentum, work and energy principle 1.1 In-class problem A block of mass m slides down a frictionless incline (see Fig.). The block is released at height h above
More informationclass 25, Mon May 24 comments and summary Ph211 Spring 2010
class 25, Mon May 24 comments and summary Ph2 Spring 200 We start with energy and work, via calculations using Newton s Laws and kinematics. For motivation, we worked through a few simple setups. Generally
More informationPhysics 351, Spring 2017, Homework #2. Due at start of class, Friday, January 27, 2017
Physics 351, Spring 2017, Homework #2. Due at start of class, Friday, January 27, 2017 Course info is at positron.hep.upenn.edu/p351 When you finish this homework, remember to visit the feedback page at
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 informationPHYS 121: Work, Energy, Momentum, and Conservation Laws, Systems, and Center of Mass Review Sheet
Physics 121 Summer 2006 Work, Energy, Momentum Review Sheet page 1 of 12 PHYS 121: Work, Energy, Momentum, and Conservation Laws, Systems, and Center of Mass Review Sheet June 21, 2006 The Definition of
More informationMATHEMATICAL PROBLEM SOLVING, MECHANICS AND MODELLING MTHA4004Y
UNIVERSITY OF EAST ANGLIA School of Mathematics Main Series UG Examination 2014 2015 MATHEMATICAL PROBLEM SOLVING, MECHANICS AND MODELLING MTHA4004Y Time allowed: 2 Hours Attempt QUESTIONS 1 AND 2 and
More information10550 PRACTICE FINAL EXAM SOLUTIONS. x 2 4. x 2 x 2 5x +6 = lim x +2. x 2 x 3 = 4 1 = 4.
55 PRACTICE FINAL EXAM SOLUTIONS. First notice that x 2 4 x 2x + 2 x 2 5x +6 x 2x. This function is undefined at x 2. Since, in the it as x 2, we only care about what happens near x 2 an for x less than
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 information10.34: Numerical Methods Applied to Chemical Engineering. Lecture 19: Differential Algebraic Equations
10.34: Numerical Methods Applied to Chemical Engineering Lecture 19: Differential Algebraic Equations 1 Recap Differential algebraic equations Semi-explicit Fully implicit Simulation via backward difference
More informationENERGY. Conservative Forces Non-Conservative Forces Conservation of Mechanical Energy Power
ENERGY Conservative Forces Non-Conservative Forces Conservation of Mechanical Energy Power Conservative Forces A force is conservative if the work it does on an object moving between two points is independent
More informationBonus Section II: Solving Trigonometric Equations
Fry Texas A&M University Math 150 Spring 2017 Bonus Section II 260 Bonus Section II: Solving Trigonometric Equations (In your text this section is found hiding at the end of 9.6) For what values of x does
More informationP321(b), Assignement 1
P31(b), Assignement 1 1 Exercise 3.1 (Fetter and Walecka) a) The problem is that of a point mass rotating along a circle of radius a, rotating with a constant angular velocity Ω. Generally, 3 coordinates
More informationTHEODORE VORONOV DIFFERENTIAL GEOMETRY. Spring 2009
[under construction] 8 Parallel transport 8.1 Equation of parallel transport Consider a vector bundle E B. We would like to compare vectors belonging to fibers over different points. Recall that this was
More informationNon-Linear Response of Test Mass to External Forces and Arbitrary Motion of Suspension Point
LASER INTERFEROMETER GRAVITATIONAL WAVE OBSERVATORY -LIGO- CALIFORNIA INSTITUTE OF TECHNOLOGY MASSACHUSETTS INSTITUTE OF TECHNOLOGY Technical Note LIGO-T980005-01- D 10/28/97 Non-Linear Response of Test
More informationKinematic representation! Iterative methods! Optimization methods
Human Kinematics Kinematic representation! Iterative methods! Optimization methods Kinematics Forward kinematics! given a joint configuration, what is the position of an end point on the structure?! Inverse
More informationNewton s Laws of Motion, Energy and Oscillations
Prof. O. B. Wright, Autumn 007 Mechanics Lecture Newton s Laws of Motion, Energy and Oscillations Reference frames e.g. displaced frame x =x+a y =y x =z t =t e.g. moving frame (t=time) x =x+vt y =y x =z
More informationFinal Exam April 30, 2013
Final Exam Instructions: You have 120 minutes to complete this exam. This is a closed-book, closed-notes exam. You are allowed to use a calculator during the exam. Usage of mobile phones and other electronic
More information4.1 Important Notes on Notation
Chapter 4. Lagrangian Dynamics (Most of the material presented in this chapter is taken from Thornton and Marion, Chap. 7) 4.1 Important Notes on Notation In this chapter, unless otherwise stated, the
More informationDiscussion Question 7A P212, Week 7 RC Circuits
Discussion Question 7A P1, Week 7 RC Circuits The circuit shown initially has the acitor uncharged, and the switch connected to neither terminal. At time t = 0, the switch is thrown to position a. C a
More informationLagrangian Dynamics: Derivations of Lagrange s Equations
Constraints and Degrees of Freedom 1.003J/1.053J Dynamics and Control I, Spring 007 Professor Thomas Peacock 4/9/007 Lecture 15 Lagrangian Dynamics: Derivations of Lagrange s Equations Constraints and
More informationQuestion 1: A particle starts at rest and moves along a cycloid whose equation is. 2ay y a
Stephen Martin PHYS 10 Homework #1 Question 1: A particle starts at rest and moves along a cycloid whose equation is [ ( ) a y x = ± a cos 1 + ] ay y a There is a gravitational field of strength g in the
More informationStructural topology, singularity, and kinematic analysis. J-P. Merlet HEPHAISTOS project INRIA Sophia-Antipolis
Structural topology, singularity, and kinematic analysis J-P. Merlet HEPHAISTOS project INRIA Sophia-Antipolis 1 Parallel robots Definitions: a closed-loop mechanism whose end-effector is linked to the
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 informationIntroductory Physics. Week 2015/06/26
2015/06/26 coordinate systems Part I Multi-particle system coordinate systems So far, we have been limiting ourselfs to the problems of the motion of single particle (except a brief introduction of coupled
More informationG j dq i + G j. q i. = a jt. and
Lagrange Multipliers Wenesay, 8 September 011 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine
More informationCHAPTER 7: OSCILLATORY MOTION REQUIRES A SET OF CONDITIONS
CHAPTER 7: OSCILLATORY MOTION REQUIRES A SET OF CONDITIONS 7.1 Period and Frequency Anything that vibrates or repeats its motion regularly is said to have oscillatory motion (sometimes called harmonic
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 informationPHY321 Homework Set 10
PHY321 Homework Set 10 1. [5 pts] A small block of mass m slides without friction down a wedge-shaped block of mass M and of opening angle α. Thetriangular block itself slides along a horizontal floor,
More informationInstructions: (62 points) Answer the following questions. SHOW ALL OF YOUR WORK. A B = A x B x + A y B y + A z B z = ( 1) + ( 1) ( 4) = 5
AP Physics C Fall, 2016 Work-Energy Mock Exam Name: Answer Key Mr. Leonard Instructions: (62 points) Answer the following questions. SHOW ALL OF YOUR WORK. (12 pts ) 1. Consider the vectors A = 2 î + 3
More informationPhysically Based Modeling: Principles and Practice Implicit Methods for Differential Equations
Pysically Based Modeling: Principles and Practice Implicit Metods for Differential Equations David Baraff Robotics Institute Carnegie Mellon University Please note: Tis document is 997 by David Baraff
More informationthe gamedesigninitiative at cornell university Lecture 17 Physics and Motion
Lecture 17 The Pedagogical Problem Physics simulation is a very complex topic No way I can address this in a few lectures Could spend an entire course talking about it CS 5643: Physically Based Animation
More informationMath 302 Test 1 Review
Math Test Review. Given two points in R, x, y, z and x, y, z, show the point x + x, y + y, z + z is on the line between these two points and is the same distance from each of them. The line is rt x, y,
More informationChapter 6. Second order differential equations
Chapter 6. Second order differential equations A second order differential equation is of the form y = f(t, y, y ) where y = y(t). We shall often think of t as parametrizing time, y position. In this case
More informationConstraint Based Control Method For Precision Formation Flight of Spacecraft AAS
Constraint Based Control Method For Precision Formation Flight of Spacecraft AAS 06-122 Try Lam Jet Propulsion Laboratory California Institute of Technology Aaron Schutte Aerospace Corporation Firdaus
More informationMechanical Energy and Simple Harmonic Oscillator
Mechanical Energy and Simple Harmonic Oscillator Simple Harmonic Motion Hooke s Law Define system, choose coordinate system. Draw free-body diagram. Hooke s Law! F spring =!kx ˆi! kx = d x m dt Checkpoint
More informationANALYTISK MEKANIK I HT 2016
Karlstads Universitet Fysik ANALYTISK MEKANIK I HT 2016 Kursens kod: FYGB08 Kursansvarig lärare: Jürgen Fuchs rum 21F 316 tel. 054-700 1817 el.mail: juerfuch@kau.se FYGB08 HT 2016 Exercises 1 2016-12-14
More informationReview of course Nonlinear control Lecture 5. Lyapunov based design Torkel Glad Lyapunov design Control Lyapunov functions Control Lyapunov function
Review of corse Nonlinear control Lectre 5 Lyapnov based design Reglerteknik, ISY, Linköpings Universitet Geometric control theory inpt-otpt linearization controller canonical form observer canonical form
More informationWork and energy. 15 m. c. Find the work done by the normal force exerted by the incline on the crate.
Work and energy 1. A 10.0-kg crate is pulled 15.0 m up along a frictionless incline as shown in the figure below. The crate starts at rest and has a final speed of 6.00 m/s. motor 15 m 5 a. Draw the free-body
More informationES.1803 Topic 13 Notes Jeremy Orloff
ES.1803 Topic 13 Notes Jeremy Orloff 13 Vector Spaces, matrices and linearity 13.1 Goals 1. Know the definition of a vector space and how to show that a given set is a vector space. 2. Know the meaning
More informationCasting Physics Simplified Part Two. Frames of Reference
Casting Physics Simplified Part Two Part one of this paper discussed physics that applies to linear motion, i.e., motion in a straight line. This section of the paper will expand these concepts to angular
More informationCh 5 Work and Energy
Ch 5 Work and Energy Energy Provide a different (scalar) approach to solving some physics problems. Work Links the energy approach to the force (Newton s Laws) approach. Mechanical energy Kinetic energy
More informationDynamics and Time Integration
Lecture 17: Dynamics and Time Integration Computer Graphics CMU 15-462/15-662, Fall 2015 Last time: animation Added motion to our model Interpolate keyframes Still a lot of work! Today: physically-based
More informationParticle Systems. University of Texas at Austin CS384G - Computer Graphics
Particle Systems University of Texas at Austin CS384G - Computer Graphics Fall 2010 Don Fussell Reading Required: Witkin, Particle System Dynamics, SIGGRAPH 97 course notes on Physically Based Modeling.
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 informationEQUIVALENT SINGLE-DEGREE-OF-FREEDOM SYSTEM AND FREE VIBRATION
1 EQUIVALENT SINGLE-DEGREE-OF-FREEDOM 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 informationPre-Class. List everything you remember about circular motion...
Pre-Class List everything you remember about circular motion... Quote of the Day I'm addicted to brake fluid......but I can stop anytime I want. Table of Contents Click on the topic to go to that section
More informationCOMPLETE ALL ROUGH WORKINGS IN THE ANSWER BOOK AND CROSS THROUGH ANY WORK WHICH IS NOT TO BE ASSESSED.
BSc/MSci EXAMINATION PHY-304 Time Allowed: Physical Dynamics 2 hours 30 minutes Date: 28 th May 2009 Time: 10:00 Instructions: Answer ALL questions in section A. Answer ONLY TWO questions from section
More informationProblem Set Number 01, MIT (Winter-Spring 2018)
Problem Set Number 01, 18.377 MIT (Winter-Spring 2018) Rodolfo R. Rosales (MIT, Math. Dept., room 2-337, Cambridge, MA 02139) February 28, 2018 Due Thursday, March 8, 2018. Turn it in (by 3PM) at the Math.
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 information3.4 Application-Spring Mass Systems (Unforced and frictionless systems)
3.4. APPLICATION-SPRING MASS SYSTEMS (UNFORCED AND FRICTIONLESS SYSTEMS)73 3.4 Application-Spring Mass Systems (Unforced and frictionless systems) Second order differential equations arise naturally when
More information