Kinematics and kinematic functions

Similar documents
Basilio Bona ROBOTICA 03CFIOR 1

Kinematics and kinematic functions

Announcements: Warm-up Exercise:

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

From Particles to Rigid Bodies

Trajectory planning in Cartesian space

Solutions from Chapter 9.1 and 9.2

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

THE 2-BODY PROBLEM. FIGURE 1. A pair of ellipses sharing a common focus. (c,b) c+a ROBERT J. VANDERBEI

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1.

Chapter 2. First Order Scalar Equations

ME 391 Mechanical Engineering Analysis

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.

Distance Between Two Ellipses in 3D

Starting from a familiar curve

Lab #2: Kinematics in 1-Dimension

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...

Some Basic Information about M-S-D Systems

Two Coupled Oscillators / Normal Modes

The motions of the celt on a horizontal plane with viscous friction

CSE 3802 / ECE Numerical Methods in Scientific Computation. Jinbo Bi. Department of Computer Science & Engineering

Chapter 3 Boundary Value Problem

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.

EE 315 Notes. Gürdal Arslan CLASS 1. (Sections ) What is a signal?

Elements of Computer Graphics

Lecture 20: Riccati Equations and Least Squares Feedback Control

Non-uniform circular motion *

Chapter 6. Systems of First Order Linear Differential Equations

L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS. NA568 Mobile Robotics: Methods & Algorithms

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004

Traveling Waves. Chapter Introduction

10. State Space Methods

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

2. Nonlinear Conservation Law Equations

MATH 128A, SUMMER 2009, FINAL EXAM SOLUTION

This is an example to show you how SMath can calculate the movement of kinematic mechanisms.

Lecture 10: Wave equation, solution by spherical means

Multi-scale 2D acoustic full waveform inversion with high frequency impulsive source

Week 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem)

Differential Equations

ROBOTICS 01PEEQW. Basilio Bona DAUIN Politecnico di Torino

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

Solving a System of Nonlinear Functional Equations Using Revised New Iterative Method

SUMMARY. CS380: Introduction to Computer Graphics Quaternions Chapter 7. Min H. Kim KAIST School of Computing 18/04/06.

Overview. COMP14112: Artificial Intelligence Fundamentals. Lecture 0 Very Brief Overview. Structure of this course

Computation of the Effect of Space Harmonics on Starting Process of Induction Motors Using TSFEM

4.1.1 Mindlin plates: Bending theory and variational formulation

KINEMATICS IN ONE DIMENSION

STATE-SPACE MODELLING. A mass balance across the tank gives:

Theory of! Partial Differential Equations!

Theory of! Partial Differential Equations-I!

5.2 GRAPHICAL VELOCITY ANALYSIS Polygon Method

R t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t

Online Appendix to Solution Methods for Models with Rare Disasters

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

HOMEWORK # 2: MATH 211, SPRING Note: This is the last solution set where I will describe the MATLAB I used to make my pictures.

HW6: MRI Imaging Pulse Sequences (7 Problems for 100 pts)

Estimation of Poses with Particle Filters

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

Notes on Kalman Filtering

LAPLACE TRANSFORM AND TRANSFER FUNCTION

1 1 + x 2 dx. tan 1 (2) = ] ] x 3. Solution: Recall that the given integral is improper because. x 3. 1 x 3. dx = lim dx.

Translational-Rotational Motion of Earth Artificial Satellite (EAS) in Hill's Gravity Field

Vehicle Arrival Models : Headway

Class Meeting # 10: Introduction to the Wave Equation

Modal identification of structures from roving input data by means of maximum likelihood estimation of the state space model

Kinematics Vocabulary. Kinematics and One Dimensional Motion. Position. Coordinate System in One Dimension. Kinema means movement 8.

CHEMICAL KINETICS: 1. Rate Order Rate law Rate constant Half-life Temperature Dependence

Chapter Q1. We need to understand Classical wave first. 3/28/2004 H133 Spring

A DELAY-DEPENDENT STABILITY CRITERIA FOR T-S FUZZY SYSTEM WITH TIME-DELAYS

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality

System of Linear Differential Equations

Final Spring 2007

Single and Double Pendulum Models

Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi

Position, Velocity, and Acceleration

The fundamental mass balance equation is ( 1 ) where: I = inputs P = production O = outputs L = losses A = accumulation

EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE

Lecture 4 Kinetics of a particle Part 3: Impulse and Momentum

We here collect a few numerical tests, in order to put into evidence the potentialities of HBVMs [4, 6, 7].

Math Final Exam Solutions

1 Review of Zero-Sum Games

KNOWLEDGE about the mechanical state of a system is

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x

Sections 2.2 & 2.3 Limit of a Function and Limit Laws

Math 334 Fall 2011 Homework 11 Solutions

Continuous Time. Time-Domain System Analysis. Impulse Response. Impulse Response. Impulse Response. Impulse Response. ( t) + b 0.

Concourse Math Spring 2012 Worked Examples: Matrix Methods for Solving Systems of 1st Order Linear Differential Equations

Augmented Reality II - Kalman Filters - Gudrun Klinker May 25, 2004

Outline of Topics. Analysis of ODE models with MATLAB. What will we learn from this lecture. Aim of analysis: Why such analysis matters?

2. VECTORS. R Vectors are denoted by bold-face characters such as R, V, etc. The magnitude of a vector, such as R, is denoted as R, R, V

Math 315: Linear Algebra Solutions to Assignment 6

The Contradiction within Equations of Motion with Constant Acceleration

New effective moduli of isotropic viscoelastic composites. Part I. Theoretical justification

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes

Improved Approximate Solutions for Nonlinear Evolutions Equations in Mathematical Physics Using the Reduced Differential Transform Method

Design of a control system

BU Macro BU Macro Fall 2008, Lecture 4

Transcription:

Kinemaics and kinemaic funcions Kinemaics deals wih he sudy of four funcions (called kinemaic funcions or KFs) ha mahemaically ransform join variables ino caresian variables and vice versa Direc Posiion KF: from join space variables o ask space pose Inverse Posiion KF: from ask space pose o join space variables Direc Velociy KF: from join space velociies o ask space velociies Inverse Velociy KF: from ask space velociies o join space velociies

Kinemaic funcions Posiions and velociies of wha? I can be he posiion or velociy of any poin of he robo, bu, usually, is he las RF posiion and velociy LAST RF R P R BASE 0? Wha is he relaion beween hese wo RF?

Kinemaic funcions The firs sep for defining KFs is o fix a reference frame on each body In general, o move from a RF o he following one, i is necessary o define 6 parameers (hree ranslaions of he RF origin + hree angles of he RF roaion) A number of convenions were inroduced o reduce henumber of parameers and o find a common way o describe he relaive posiion of reference frames Denavi Harenberg convenions (1955) were he firs o be inroduced and are widely used in indusry (wih some minor modificaions)

A procedure o compue he KFs 1 1. Selec and idenify bodies and he relaive consrains 2. Se he RF on each body 3. Compue he homogeneous ransformaion beween he base RF and he final RF 4. Exrac he direc posiion KF from he homogeneous ransformaion 5. Compue he inverse posiion KF 6. Inverse velociy KF: analyical or geomerical approach

A procedure o compue he KFs 2 1. Selec and idenify links and joins 2. Se all RFs using DH convenions 3. Define consan geomerical parameers q () 4 q () 5 q () 2 q () 6 q () 1 R P FINAL R 0 BASE 0 T P 0 0 0 () () q () q T q = R P P P T 0 1

Join and caresian variables Join variables q1() q2() q 3() q() = q4() q 5() q6() Direc KF Task/caresian variables/pose x 1() x2() x 3() p() = α1() α2() α3() Inverse KF posiion orienaion 1 p() = f ( q () ) q() = f ( p() )

Direc posiion KF 1 q1() p1() q 2() p2() q ( ()) 3() p3() x q q() = ( ()) = = p q q 4() p4() q ( ()) 5() p5() α q q6() p6() posiion orienaion 0 0 0 0 1 2 3 4 5 6 P P P 1 ( q1) 2 ( q2) 3 ( q3) 4 ( q4) 5 ( q5) 6 ( q6) R T = T T T T T T T P = 1 0 T orienaion posiion 0 0 0 0 R = R ( q () ) = ( q () ) P P P P

Direc posiion KF 2 T 0 P 0 0 RP P = 0 T 1 Direc caresian posiion KF: easy x x 1 1 x 0 2 P 2 x 3 3 Direc caresian orienaion KF: no so easy o compue, bu no difficul We will solve he problem in he following slides

Direc posiion KF 3 ( () ) α( q() ) Rq We wan o compue angles from he roaion marix. Bu i is imporan o decide which represenaion o use α ( q() )? Euler angles RPY angles Quaernions Axis angle represenaion

q () ( ) 2 () x q q () 3 pq ( ( )) = = q( ) q () α ( ) 4 () q q () 5 q () 6 Inverse posiion KF 1 q () 1 This KF is imporan, since conrol acion are applied o he join moors, while he user wans o work wih caresian posiions and orienaions q () 3 q () 1 q () 2 ( () ) x q α( q() )

Inverse posiion KF 2 1. The problem is complex and here is no clear recipe o solve i 2. If a spherical wris is presen, hen a soluion is guaraneed, bu we mus find i... 3. There are several possibiliies Use brue force or previous soluions found for similar chains Use inverse velociy KF Usesymbolicmanipulaion symbolic manipulaion programs (compuer algebra sysems as Mahemaica, Maple, Maxima,, Lisp) Ieraively compue an approximaed numerical expression for he nonlinear equaion (Newon mehod or ohers) ( ) 1 ( ) { 1 q() f ( p() ) } q f p 1 () = () q() f p() = 0 min () ()

Direc velociy KF 1 Linear and angular direc velociy KF q () p () 1 1 q () p () 2 2 (), () x q q q () p () 3 3 q () = () q () p = p () = 4 4 (), () q q q () p () 5 5 q () p () 6 6 ( ) Linear velociy ( (), () ) vq q = (), () ω q q α ( ) ( ) Angular velociy

Direc velociy KF 2 A bi brief review of mahemaical i noaions i General rule d pq (()) d d ( () ) x q d = d ( () ) α q d d d f ( q () ) = f ( q (),, q (),, q () ) i i 1 j n d d f i f i f i = q () + + q () + + q () 1 j n q q q 1 j n

Direc velociy KF 3 q () 1 d f f f ( ()) i i i f q () T q = ( () ) () i = j fi J q q d q q q 1 j n q () n f f f 1 1 1 JACOBIAN q q q 1 j n q () 1 d ( ) f f f i i i f q() = q () ( () ) () j d = q q q J q q 1 j n q () n f f f m m m q q q 1 j n

Furher noes on he Jacobians Velociykinemaics is characerizedby Jacobians There are wo ypes of Jacobians: Geomerical Jacobian Analyical ljacobian J g J a ( ) p () = J q() q () also called Task Jacobian The firs one is relaed o Geomerical Velociies v p x = = ω J q g The second one is relaed o Analyical Velociies p x = = α α J q a

Geomerical and Analyical velociies Wha is he difference beween hese wo angular velociies? On he conrary, linear velociies do no have his problem: analyical and geomerical velociies are he same vecor, ha can be inegraed o give he caresian posiion

Derivaive of a roaion marix

Some useful formulas