Nonlinear Regulation for Motorcycle Maneuvering

Size: px
Start display at page:

Download "Nonlinear Regulation for Motorcycle Maneuvering"

Transcription

1 Nonlinear Regulation for Motorcycle Maneuvering John Hauser Univ of Colorado in collaboration with Alessandro Saccon* & Ruggero Frezza, Univ Padova * asaccon@dei.unipd.it for dissertation

2 aggressive maneuvering we seek to understand dynamics and control issues of aggressively maneuvering systems an opinion: maneuvering is one of the most common and interesting ways that nonlinear effects are seen in control systems examples include aircraft, motorcycles, skiers

3 motorcycles motorcycles possess unstable nonlinear dynamics coupling of inputs control vector field sign changes nonminimum phase response broad range of operation: mph, lateral g s rapidly changing trajectories: turn-in, chicane, accel, braking just plain fun! Note: we do not intend to replace rider

4 motorcycles: engineering objectives provide strategies to test-drive various virtual prototypes: human rider is not able to evaluate virtual needed: a virtual rider (a control system) to enable complex maneuvering near the limits of performance (max roll, max lateral accel) and that can exploit input coupling better understand performance tradeoffs: what setup (bike geometry, tires, suspension, ) is best for different circuits.

5 aggressive Moto maneuvers are desired! Loris Capirossi

6 Circuit Catalunya

7 max acceleration and braking Loris Capirossi Valentino Rossi

8 complex Moto behaviors are possible! Isle of Man 1999

9 Hierarchy of models: motorcycle specifics - nonholonomic motorcycle infinitely sticky tires, simplified geometry - sliding plane motorcycle more realistic contact forces, simplified geometry... - articulated motorcycle include suspension, chain, flexible frame, semi-empirical tire models, art / magic!

10 planning maneuvering objectives - track specification inner and outer track boundaries go fast stay on track - path or race line specification arc length parametrized curve go fast on this line - ground trajectory specification time parametrized curve leads to a desired maneuvering objective

11 test track

12 velocity profile

13 velocity and accel trajectory

14 maneuvers and maneuver regulation Given ẋ = f(x, u) and a trajectory (x(t),u(t)), t R, with ẋ(t) and ẍ(t) bdd and ẋ(t) bdd away from zero, the corresponding maneuver is the curve swept out by (x( ),u( )) together with local temporal separation. The maneuver has unique projection within a tube prop In practice, a maneuver is specified using a parametrized curve ( x(θ), ū(θ)), θ R The param θ could be time-like or arc-length s.

15 transverse dynamics Around a maneuver, choose transverse coordinates θ = 1+g 1 (ρ,u ū(θ)) ρ = A(θ)ρ + B(θ)(u ū(θ)) + g 2 (ρ,u ū(θ)) locally, we may eliminate time d dθ ρ = A(θ)ρ + B(θ)(u ū(θ)) + f 2(ρ,u ū(θ)) key: study stability, control, robustness of time-varying nonlinear control systems discuss

16 nonholonomic motorcycle model. nonholonomic car model ẋ = v cos ψ ẏ = v sinψ v = u 1 ψ = vσ σ = u 2 coupled roll dynamics R =1/σ h ϕ = g sinϕ ((1 hσ sinϕ)σv 2 + b ψ)cosϕ ψ ϕ (x, y) b h p δ

17 to get a trajectory path and velocity profile directly provide a nonholonomic car trajectory the desired motorcycle maneuver is determined by lifting the car trajectory to a moto traj, adding a roll traj in this fashion, the class of motorcycle trajectories is parametrized by the family of smooth curves in the plane

18 lifting to an executable Moto trajectory given the desired flatland traj, find a roll trajectory consistent with, roughly, h ϕ = g sinϕ a lat (t)cosϕ + u hog after dynamic embedding, we optimize away the hand of God for now, we do the whole trajectory

19 quasi-static static roll trajectory when the desired flatland traj is a constant speed, constant radius circle, there is a static roll trajectory given by for more dynamic flatland trajectories, we define the quasi-static roll trajectory according to we expect (hope) that the desired roll traj is close to this!

20 achievable motorcycle trajectories problem: given a smooth velocity-curvature profile, find, if possible, an upright roll trajectory satisfying with h ϕ = g sinϕ a lat (t)cosϕ a lat (t) =[σv 2 + b( vσ + v σ)](t) in fact, such inverted pendulum dynamics is always a part of the dynamics of every motorcycle also, the lateral acceleration will, in general, be much more complicated and may not be smooth

21 wanted: an upright soln of the geometric story ~Thm: if ϕ( ) is an upright soln, the phase traj lies in 6 phase plane pi/2 -pi/4 0 pi/4 pi/2

22 Thm: existence of an upright roll traj with a bdd that is const before some t 0 possesses an upright soln 6 phase plane pi/2 -pi/4 0 pi/4 pi/2

23 dynamics w.r.t.. quasi-static static roll traj defining the quasi-static roll angle and total acceleration the roll dynamics is given by inverted pendulum dynamics with gravity that varies in strength and direction we seek a bounded traj of the driven unstable system.

24 bounded solutions: dichotomy when will a system like have a bounded solution? [and with upright roll] the unique bounded solution of the LTI system is given by.

25 bounded solutions: dichotomy can we find a bounded solution for the time-varying linear system the LTI system is hyperbolic for time-varying systems, we seek a dichotomy?. [this will be used to show the TV nonlinear sys has a soln]

26 bounded solutions: dichotomy Thm: the unique soln of is given by the noncausal bounded operator where c(.) and d(.) are nonl filtered versions of α(.)

27 solution algorithm Fact: under some conditions, the unique soln of can be computed by the algo and, furthermore, is small.. (note: above optimization can also be used) ³ h(t) α/2 e α t

28 maneuver regulationg with an executable trajectory in hand (reparametrized by arclength), we may write the system dynamics in transverse maneuver coordinates so that the transverse dynamics are given by

29 maneuver regulation MP maneuver regulation may then be implemented using possibly subject to some constraints (e.g., lateral accel) a first order controller may be obtain by solving a TV Riccati equation (where time is arclength)

30 cost function design how should we choose Q and R? the many heuristics suggested in the literature did not seem effective to us performance requires a certain speed of response physical motion requires a restricted speed of response nonlinearities (seem to) require a certain uniformity of response under aggressive maneuvering plus all the usual control performance expectations...

31 Q = I, R = I not too interesting σ root locus too fast desired region

32 another heuristic for Q & R design get a desired lateral response first for SS system (e.g., place poles for driving in a high g circle) solve, if able, an inverse optimal control problem (must satisfy return difference ineq ) requiring Q, R > 0 (resulting 5x5 Q is far from diagonal) [can be done as a convex problem---we use SeDuMi] augment the lateral Q, R with a choice of Q, R for the (scalar) longitudinal subsystem evaluate over a range of velocity and lateral accel and iterate reasonable results have been obtained for nonholonomic motorcycle

33 Q, R obtained by inverse opt heuristic σ root locus

34 Q, R obtained by inverse opt heuristic v root locus v root locus

35 example performance eval

36 remarks robustness: we have applied maneuver regulation (based on simple moto model) to regulation of high fidelity motorcycle model (multi-body)---with great success! Ale Saccon for details (in his dissertation).

37

Aircraft Maneuver Regulation: a Receding Horizon Backstepping Approach

Aircraft Maneuver Regulation: a Receding Horizon Backstepping Approach Aircraft Maneuver Regulation: a Receding Horizon Backstepping Approach Giuseppe Notarstefano and Ruggero Frezza Abstract Coordinated flight is a nonholonomic constraint that implies no sideslip of an aircraft.

More information

Robotics, Geometry and Control - A Preview

Robotics, Geometry and Control - A Preview Robotics, Geometry and Control - A Preview Ravi Banavar 1 1 Systems and Control Engineering IIT Bombay HYCON-EECI Graduate School - Spring 2008 Broad areas Types of manipulators - articulated mechanisms,

More information

Single-track models of an A-double heavy vehicle combination

Single-track models of an A-double heavy vehicle combination Single-track models of an A-double heavy vehicle combination PETER NILSSON KRISTOFFER TAGESSON Department of Applied Mechanics Division of Vehicle Engineering and Autonomous Systems Vehicle Dynamics Group

More information

An homotopy method for exact tracking of nonlinear nonminimum phase systems: the example of the spherical inverted pendulum

An homotopy method for exact tracking of nonlinear nonminimum phase systems: the example of the spherical inverted pendulum 9 American Control Conference Hyatt Regency Riverfront, St. Louis, MO, USA June -, 9 FrA.5 An homotopy method for exact tracking of nonlinear nonminimum phase systems: the example of the spherical inverted

More information

EQUATIONS OF MOTION: NORMAL AND TANGENTIAL COORDINATES (Section 13.5)

EQUATIONS OF MOTION: NORMAL AND TANGENTIAL COORDINATES (Section 13.5) EQUATIONS OF MOTION: NORMAL AND TANGENTIAL COORDINATES (Section 13.5) Today s Objectives: Students will be able to apply the equation of motion using normal and tangential coordinates. APPLICATIONS Race

More information

Physics 12. Unit 5 Circular Motion and Gravitation Part 1

Physics 12. Unit 5 Circular Motion and Gravitation Part 1 Physics 12 Unit 5 Circular Motion and Gravitation Part 1 1. Nonlinear motions According to the Newton s first law, an object remains its tendency of motion as long as there is no external force acting

More information

Nonholonomic Constraints Examples

Nonholonomic Constraints Examples Nonholonomic Constraints Examples Basilio Bona DAUIN Politecnico di Torino July 2009 B. Bona (DAUIN) Examples July 2009 1 / 34 Example 1 Given q T = [ x y ] T check that the constraint φ(q) = (2x + siny

More information

Time-Optimal Automobile Test Drives with Gear Shifts

Time-Optimal Automobile Test Drives with Gear Shifts Time-Optimal Control of Automobile Test Drives with Gear Shifts Christian Kirches Interdisciplinary Center for Scientific Computing (IWR) Ruprecht-Karls-University of Heidelberg, Germany joint work with

More information

Chapter 6. Force and Motion-II (Friction, Drag, Circular Motion)

Chapter 6. Force and Motion-II (Friction, Drag, Circular Motion) Chapter 6 Force and Motion-II (Friction, Drag, Circular Motion) 6.2 Frictional Force: Motion of a crate with applied forces There is no attempt at sliding. Thus, no friction and no motion. NO FRICTION

More information

Control Systems I. Lecture 2: Modeling. Suggested Readings: Åström & Murray Ch. 2-3, Guzzella Ch Emilio Frazzoli

Control Systems I. Lecture 2: Modeling. Suggested Readings: Åström & Murray Ch. 2-3, Guzzella Ch Emilio Frazzoli Control Systems I Lecture 2: Modeling Suggested Readings: Åström & Murray Ch. 2-3, Guzzella Ch. 2-3 Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich September 29, 2017 E. Frazzoli

More information

The Role of Zero Dynamics in Aerospace Systems

The Role of Zero Dynamics in Aerospace Systems The Role of Zero Dynamics in Aerospace Systems A Case Study in Control of Hypersonic Vehicles Andrea Serrani Department of Electrical and Computer Engineering The Ohio State University Outline q Issues

More information

Lecture 10. Example: Friction and Motion

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

Control Systems I. Lecture 2: Modeling and Linearization. Suggested Readings: Åström & Murray Ch Jacopo Tani

Control Systems I. Lecture 2: Modeling and Linearization. Suggested Readings: Åström & Murray Ch Jacopo Tani Control Systems I Lecture 2: Modeling and Linearization Suggested Readings: Åström & Murray Ch. 2-3 Jacopo Tani Institute for Dynamic Systems and Control D-MAVT ETH Zürich September 28, 2018 J. Tani, E.

More information

13 Path Planning Cubic Path P 2 P 1. θ 2

13 Path Planning Cubic Path P 2 P 1. θ 2 13 Path Planning Path planning includes three tasks: 1 Defining a geometric curve for the end-effector between two points. 2 Defining a rotational motion between two orientations. 3 Defining a time function

More information

Control of a Car-Like Vehicle with a Reference Model and Particularization

Control of a Car-Like Vehicle with a Reference Model and Particularization Control of a Car-Like Vehicle with a Reference Model and Particularization Luis Gracia Josep Tornero Department of Systems and Control Engineering Polytechnic University of Valencia Camino de Vera s/n,

More information

Examination paper for TMA4195 Mathematical Modeling

Examination paper for TMA4195 Mathematical Modeling Department of Mathematical Sciences Examination paper for TMA4195 Mathematical Modeling Academic contact during examination: Elena Celledoni Phone: 48238584, 73593541 Examination date: 11th of December

More information

Chapter 8: Newton s Laws Applied to Circular Motion

Chapter 8: Newton s Laws Applied to Circular Motion Chapter 8: Newton s Laws Applied to Circular Motion Circular Motion Milky Way Galaxy Orbital Speed of Solar System: 220 km/s Orbital Period: 225 Million Years Mercury: 48 km/s Venus: 35 km/s Earth: 30

More information

Control of Mobile Robots Prof. Luca Bascetta

Control of Mobile Robots Prof. Luca Bascetta Control of Mobile Robots Prof. Luca Bascetta EXERCISE 1 1. Consider a wheel rolling without slipping on the horizontal plane, keeping the sagittal plane in the vertical direction. Write the expression

More information

Lecture D4 - Intrinsic Coordinates

Lecture D4 - Intrinsic Coordinates J. Peraire 16.07 Dynamics Fall 2004 Version 1.1 Lecture D4 - Intrinsic Coordinates In lecture D2 we introduced the position, velocity and acceleration vectors and referred them to a fixed cartesian coordinate

More information

Polynomial-Time Verification of PCTL Properties of MDPs with Convex Uncertainties and its Application to Cyber-Physical Systems

Polynomial-Time Verification of PCTL Properties of MDPs with Convex Uncertainties and its Application to Cyber-Physical Systems Polynomial-Time Verification of PCTL Properties of MDPs with Convex Uncertainties and its Application to Cyber-Physical Systems Alberto Puggelli DREAM Seminar - November 26, 2013 Collaborators and PIs:

More information

EML5311 Lyapunov Stability & Robust Control Design

EML5311 Lyapunov Stability & Robust Control Design EML5311 Lyapunov Stability & Robust Control Design 1 Lyapunov Stability criterion In Robust control design of nonlinear uncertain systems, stability theory plays an important role in engineering systems.

More information

Problem 1: Ship Path-Following Control System (35%)

Problem 1: Ship Path-Following Control System (35%) Problem 1: Ship Path-Following Control System (35%) Consider the kinematic equations: Figure 1: NTNU s research vessel, R/V Gunnerus, and Nomoto model: T ṙ + r = Kδ (1) with T = 22.0 s and K = 0.1 s 1.

More information

Unit Speed Curves. Recall that a curve Α is said to be a unit speed curve if

Unit Speed Curves. Recall that a curve Α is said to be a unit speed curve if Unit Speed Curves Recall that a curve Α is said to be a unit speed curve if The reason that we like unit speed curves that the parameter t is equal to arc length; i.e. the value of t tells us how far along

More information

CIRCULAR MOTION AND GRAVITATION

CIRCULAR MOTION AND GRAVITATION CIRCULAR MOTION AND GRAVITATION An object moves in a straight line if the net force on it acts in the direction of motion, or is zero. If the net force acts at an angle to the direction of motion at any

More information

Curvilinear Motion: Normal and Tangential Components

Curvilinear Motion: Normal and Tangential Components Curvilinear Motion: Normal and Tangential Components Coordinate System Provided the path of the particle is known, we can establish a set of n and t coordinates having a fixed origin, which is coincident

More information

NEWTON S LAWS OF MOTION, EQUATIONS OF MOTION, & EQUATIONS OF MOTION FOR A SYSTEM OF PARTICLES

NEWTON S LAWS OF MOTION, EQUATIONS OF MOTION, & EQUATIONS OF MOTION FOR A SYSTEM OF PARTICLES NEWTON S LAWS OF MOTION, EQUATIONS OF MOTION, & EQUATIONS OF MOTION FOR A SYSTEM OF PARTICLES Objectives: Students will be able to: 1. Write the equation of motion for an accelerating body. 2. Draw the

More information

Real-time trajectory generation technique for dynamic soaring UAVs

Real-time trajectory generation technique for dynamic soaring UAVs Real-time trajectory generation technique for dynamic soaring UAVs Naseem Akhtar James F Whidborne Alastair K Cooke Department of Aerospace Sciences, Cranfield University, Bedfordshire MK45 AL, UK. email:n.akhtar@cranfield.ac.uk

More information

Name: School: Class: Teacher: Date:

Name: School: Class: Teacher: Date: ame: School: Class: Teacher: Date: Materials needed: Pencil, stopwatch, and scientific calculator d v λ f λ λ Wave Pool Side View During wave cycles, waves crash along the shore every few seconds. The

More information

AP PHYSICS Chapter 5. Friction Inclines Circular Motion

AP PHYSICS Chapter 5. Friction Inclines Circular Motion AP PHYSICS Chapter 5 Friction Inclines Circular Motion Friction Force that opposes motion due to contact between surfaces. Depends on: Composition and Qualities of the two surfaces in contact (μ) Roughness,

More information

Analysis of Critical Speed Yaw Scuffs Using Spiral Curves

Analysis of Critical Speed Yaw Scuffs Using Spiral Curves Analysis of Critical Speed Yaw Scuffs Using Spiral Curves Jeremy Daily Department of Mechanical Engineering University of Tulsa PAPER #2012-01-0606 Presentation Overview Review of Critical Speed Yaw Analysis

More information

The PVTOL Aircraft. 2.1 Introduction

The PVTOL Aircraft. 2.1 Introduction 2 The PVTOL Aircraft 2.1 Introduction We introduce in this chapter the well-known Planar Vertical Take-Off and Landing (PVTOL) aircraft problem. The PVTOL represents a challenging nonlinear systems control

More information

First Year Physics: Prelims CP1. Classical Mechanics: Prof. Neville Harnew. Problem Set III : Projectiles, rocket motion and motion in E & B fields

First Year Physics: Prelims CP1. Classical Mechanics: Prof. Neville Harnew. Problem Set III : Projectiles, rocket motion and motion in E & B fields HT017 First Year Physics: Prelims CP1 Classical Mechanics: Prof Neville Harnew Problem Set III : Projectiles, rocket motion and motion in E & B fields Questions 1-10 are standard examples Questions 11-1

More information

A Model-Free Control System Based on the Sliding Mode Control Method with Applications to Multi-Input-Multi-Output Systems

A Model-Free Control System Based on the Sliding Mode Control Method with Applications to Multi-Input-Multi-Output Systems Proceedings of the 4 th International Conference of Control, Dynamic Systems, and Robotics (CDSR'17) Toronto, Canada August 21 23, 2017 Paper No. 119 DOI: 10.11159/cdsr17.119 A Model-Free Control System

More information

Decentralized Stabilization of Heterogeneous Linear Multi-Agent Systems

Decentralized Stabilization of Heterogeneous Linear Multi-Agent Systems 1 Decentralized Stabilization of Heterogeneous Linear Multi-Agent Systems Mauro Franceschelli, Andrea Gasparri, Alessandro Giua, and Giovanni Ulivi Abstract In this paper the formation stabilization problem

More information

Integration of Constraint Equations in Problems of a Disc and a Ball Rolling on a Horizontal Plane

Integration of Constraint Equations in Problems of a Disc and a Ball Rolling on a Horizontal Plane Integration of Constraint Equations in Problems of a Disc and a Ball Rolling on a Horizontal Plane Eugeny A. Mityushov Ural Federal University Department of Theoretical Mechanics Prospect Mira 19 62 2

More information

Interior-Point Methods for Linear Optimization

Interior-Point Methods for Linear Optimization Interior-Point Methods for Linear Optimization Robert M. Freund and Jorge Vera March, 204 c 204 Robert M. Freund and Jorge Vera. All rights reserved. Linear Optimization with a Logarithmic Barrier Function

More information

Drag Force. Drag is a mechanical force generated when a solid moves through a fluid. Is Air fluid?

Drag Force. Drag is a mechanical force generated when a solid moves through a fluid. Is Air fluid? Feline Pesematology Drag Force Drag is a mechanical force generated when a solid moves through a fluid. Is Air fluid? Drag factors Does drag increase/decrease with 1. Density of fluid? 2. Velocity of the

More information

3 Space curvilinear motion, motion in non-inertial frames

3 Space curvilinear motion, motion in non-inertial frames 3 Space curvilinear motion, motion in non-inertial frames 3.1 In-class problem A rocket of initial mass m i is fired vertically up from earth and accelerates until its fuel is exhausted. The residual mass

More information

28. 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) 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 information

Feedback Linearization based Arc Length Control for Gas Metal Arc Welding

Feedback Linearization based Arc Length Control for Gas Metal Arc Welding 5 American Control Conference June 8-, 5. Portland, OR, USA FrA5.6 Feedback Linearization based Arc Length Control for Gas Metal Arc Welding Jesper S. Thomsen Abstract In this paper a feedback linearization

More information

Robotics & Automation. Lecture 25. Dynamics of Constrained Systems, Dynamic Control. John T. Wen. April 26, 2007

Robotics & Automation. Lecture 25. Dynamics of Constrained Systems, Dynamic Control. John T. Wen. April 26, 2007 Robotics & Automation Lecture 25 Dynamics of Constrained Systems, Dynamic Control John T. Wen April 26, 2007 Last Time Order N Forward Dynamics (3-sweep algorithm) Factorization perspective: causal-anticausal

More information

Practice Exam 2. Name: Date: ID: A. Multiple Choice Identify the choice that best completes the statement or answers the question.

Practice Exam 2. Name: Date: ID: A. Multiple Choice Identify the choice that best completes the statement or answers the question. Name: Date: _ Practice Exam 2 Multiple Choice Identify the choice that best completes the statement or answers the question. 1. A roller-coaster car has a mass of 500 kg when fully loaded with passengers.

More information

DYNAMICS ME HOMEWORK PROBLEM SETS

DYNAMICS ME HOMEWORK PROBLEM SETS DYNAMICS ME 34010 HOMEWORK PROBLEM SETS Mahmoud M. Safadi 1, M.B. Rubin 2 1 safadi@technion.ac.il, 2 mbrubin@technion.ac.il Faculty of Mechanical Engineering Technion Israel Institute of Technology Spring

More information

The Inverted Pendulum

The Inverted Pendulum Lab 1 The Inverted Pendulum Lab Objective: We will set up the LQR optimal control problem for the inverted pendulum and compute the solution numerically. Think back to your childhood days when, for entertainment

More information

Chapter 2 Review of Linear and Nonlinear Controller Designs

Chapter 2 Review of Linear and Nonlinear Controller Designs Chapter 2 Review of Linear and Nonlinear Controller Designs This Chapter reviews several flight controller designs for unmanned rotorcraft. 1 Flight control systems have been proposed and tested on a wide

More information

Physics 207 Lecture 10. Lecture 10. Employ Newton s Laws in 2D problems with circular motion

Physics 207 Lecture 10. Lecture 10. Employ Newton s Laws in 2D problems with circular motion Lecture 10 Goals: Employ Newton s Laws in 2D problems with circular motion Assignment: HW5, (Chapters 8 & 9, due 3/4, Wednesday) For Tuesday: Finish reading Chapter 8, start Chapter 9. Physics 207: Lecture

More information

Upon successful completion of MATH 220, the student will be able to:

Upon successful completion of MATH 220, the student will be able to: MATH 220 Matrices Upon successful completion of MATH 220, the student will be able to: 1. Identify a system of linear equations (or linear system) and describe its solution set 2. Write down the coefficient

More information

Announcements. Introduction and Rectilinear Kinematics: Continuous Motion - Sections

Announcements. Introduction and Rectilinear Kinematics: Continuous Motion - Sections Announcements Week-of-prayer schedule (10:45-11:30) Introduction and Rectilinear Kinematics: Continuous Motion - Sections 12.1-2 Today s Objectives: Students will be able to find the kinematic quantities

More information

EE16B, Spring 2018 UC Berkeley EECS. Maharbiz and Roychowdhury. Lectures 4B & 5A: Overview Slides. Linearization and Stability

EE16B, Spring 2018 UC Berkeley EECS. Maharbiz and Roychowdhury. Lectures 4B & 5A: Overview Slides. Linearization and Stability EE16B, Spring 2018 UC Berkeley EECS Maharbiz and Roychowdhury Lectures 4B & 5A: Overview Slides Linearization and Stability Slide 1 Linearization Approximate a nonlinear system by a linear one (unless

More information

Study Questions/Problems Week 4

Study Questions/Problems Week 4 Study Questions/Problems Week 4 Chapter 6 treats many topics. I have selected on average less than three problems from each topic. I suggest you do them all. Likewise for the Conceptual Questions and exercises,

More information

Calculating Mechanical Transfer Functions with COMSOL Multiphysics. Yoichi Aso Department of Physics, University of Tokyo

Calculating Mechanical Transfer Functions with COMSOL Multiphysics. Yoichi Aso Department of Physics, University of Tokyo Calculating Mechanical Transfer Functions with COMSOL Multiphysics Yoichi Aso Department of Physics, University of Tokyo Objective Suspension Point You have a pendulum like the one shown on the right.

More information

Estimation of Lateral Dynamics and Road Curvature for Two-Wheeled Vehicles: A HOSM Observer approach

Estimation of Lateral Dynamics and Road Curvature for Two-Wheeled Vehicles: A HOSM Observer approach Preprints of the 19th World Congress The International Federation of Automatic Control Cape Town, South Africa. August 24-29, 214 Estimation of Lateral Dynamics and Road Curvature for Two-Wheeled Vehicles:

More information

CEE 271: Applied Mechanics II, Dynamics Lecture 9: Ch.13, Sec.4-5

CEE 271: Applied Mechanics II, Dynamics Lecture 9: Ch.13, Sec.4-5 1 / 40 CEE 271: Applied Mechanics II, Dynamics Lecture 9: Ch.13, Sec.4-5 Prof. Albert S. Kim Civil and Environmental Engineering, University of Hawaii at Manoa 2 / 40 EQUATIONS OF MOTION:RECTANGULAR COORDINATES

More information

AP Physics 1 Lesson 10.a Law of Universal Gravitation Homework Outcomes

AP Physics 1 Lesson 10.a Law of Universal Gravitation Homework Outcomes AP Physics 1 Lesson 10.a Law of Universal Gravitation Homework Outcomes 1. Use Law of Universal Gravitation to solve problems involving different masses. 2. Determine changes in gravitational and kinetic

More information

Linearization problem. The simplest example

Linearization problem. The simplest example Linear Systems Lecture 3 1 problem Consider a non-linear time-invariant system of the form ( ẋ(t f x(t u(t y(t g ( x(t u(t (1 such that x R n u R m y R p and Slide 1 A: f(xu f(xu g(xu and g(xu exist and

More information

Linear and Nonlinear Oscillators (Lecture 2)

Linear and Nonlinear Oscillators (Lecture 2) Linear and Nonlinear Oscillators (Lecture 2) January 25, 2016 7/441 Lecture outline A simple model of a linear oscillator lies in the foundation of many physical phenomena in accelerator dynamics. A typical

More information

Curves - A lengthy story

Curves - A lengthy story MATH 2401 - Harrell Curves - A lengthy story Lecture 4 Copyright 2007 by Evans M. Harrell II. Reminder What a lonely archive! Who in the cast of characters might show up on the test? Curves r(t), velocity

More information

Just what is curvature, anyway?

Just what is curvature, anyway? MATH 2401 - Harrell Just what is curvature, anyway? Lecture 5 Copyright 2007 by Evans M. Harrell II. The osculating plane Bits of curve have a best plane. stickies on wire. Each stickie contains T and

More information

In the absence of an external force, the momentum of an object remains unchanged conservation of momentum. In this. rotating objects tend to

In the absence of an external force, the momentum of an object remains unchanged conservation of momentum. In this. rotating objects tend to Rotating objects tend to keep rotating while non- rotating objects tend to remain non-rotating. In the absence of an external force, the momentum of an object remains unchanged conservation of momentum.

More information

The single track model

The single track model The single track model Dr. M. Gerdts Uniersität Bayreuth, SS 2003 Contents 1 Single track model 1 1.1 Geometry.................................... 1 1.2 Computation of slip angles...........................

More information

Types of Forces. Pressure Buoyant Force Friction Normal Force

Types of Forces. Pressure Buoyant Force Friction Normal Force Types of Forces Pressure Buoyant Force Friction Normal Force Pressure Ratio of Force Per Unit Area p = F A P = N/m 2 = 1 pascal (very small) P= lbs/in 2 = psi = pounds per square inch Example: Snow Shoes

More information

DISTURBANCES MONITORING FROM CONTROLLER STATES

DISTURBANCES MONITORING FROM CONTROLLER STATES DISTURBANCES MONITORING FROM CONTROLLER STATES Daniel Alazard Pierre Apkarian SUPAERO, av. Edouard Belin, 3 Toulouse, France - Email : alazard@supaero.fr Mathmatiques pour l Industrie et la Physique, Université

More information

Comparison of two non-linear model-based control strategies for autonomous vehicles

Comparison of two non-linear model-based control strategies for autonomous vehicles Comparison of two non-linear model-based control strategies for autonomous vehicles E. Alcala*, L. Sellart**, V. Puig*, J. Quevedo*, J. Saludes*, D. Vázquez** and A. López** * Supervision & Security of

More information

Control of Mobile Robots

Control of Mobile Robots Control of Mobile Robots Regulation and trajectory tracking Prof. Luca Bascetta (luca.bascetta@polimi.it) Politecnico di Milano Dipartimento di Elettronica, Informazione e Bioingegneria Organization and

More information

11.6. Parametric Differentiation. Introduction. Prerequisites. Learning Outcomes

11.6. Parametric Differentiation. Introduction. Prerequisites. Learning Outcomes Parametric Differentiation 11.6 Introduction Often, the equation of a curve may not be given in Cartesian form y f(x) but in parametric form: x h(t), y g(t). In this section we see how to calculate the

More information

Uniform Circular Motion

Uniform Circular Motion Uniform Circular Motion Motion in a circle at constant angular speed. ω: angular velocity (rad/s) Rotation Angle The rotation angle is the ratio of arc length to radius of curvature. For a given angle,

More information

Differential Geometry of Curves

Differential Geometry of Curves Differential Geometry of Curves Cartesian coordinate system René Descartes (1596-165) (lat. Renatus Cartesius) French philosopher, mathematician, and scientist. Rationalism y Ego cogito, ergo sum (I think,

More information

MULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question.

MULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question. Exam Name MULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question. 1) Two men, Joel and Jerry, push against a wall. Jerry stops after 10 min, while Joel is

More information

Components of a Vector

Components of a Vector Vectors (Ch. 1) A vector is a quantity that has a magnitude and a direction. Examples: velocity, displacement, force, acceleration, momentum Examples of scalars: speed, temperature, mass, length, time.

More information

11.6. Parametric Differentiation. Introduction. Prerequisites. Learning Outcomes

11.6. Parametric Differentiation. Introduction. Prerequisites. Learning Outcomes Parametric Differentiation 11.6 Introduction Sometimes the equation of a curve is not be given in Cartesian form y f(x) but in parametric form: x h(t), y g(t). In this Section we see how to calculate the

More information

Linearize a non-linear system at an appropriately chosen point to derive an LTI system with A, B,C, D matrices

Linearize a non-linear system at an appropriately chosen point to derive an LTI system with A, B,C, D matrices Dr. J. Tani, Prof. Dr. E. Frazzoli 151-0591-00 Control Systems I (HS 2018) Exercise Set 2 Topic: Modeling, Linearization Discussion: 5. 10. 2018 Learning objectives: The student can mousavis@ethz.ch, 4th

More information

This is example 3 on page 44 of BGH and example (b) on page 66 of Troutman.

This is example 3 on page 44 of BGH and example (b) on page 66 of Troutman. Chapter 4 The brachistochrone This is example 3 on page 44 of BGH and example (b) on page 66 of Troutman. We seek the shape of a frictionless wire starting at the origin and ending at some point (, d)

More information

Topic # /31 Feedback Control Systems. Analysis of Nonlinear Systems Lyapunov Stability Analysis

Topic # /31 Feedback Control Systems. Analysis of Nonlinear Systems Lyapunov Stability Analysis Topic # 16.30/31 Feedback Control Systems Analysis of Nonlinear Systems Lyapunov Stability Analysis Fall 010 16.30/31 Lyapunov Stability Analysis Very general method to prove (or disprove) stability 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

On my honor, I have neither given nor received unauthorized aid on this examination.

On my honor, I have neither given nor received unauthorized aid on this examination. Instructor(s): Field/inzler PHYSICS DEPATMENT PHY 2053 Final Exam April 27, 2013 Name (print, last first): Signature: On my honor, I have neither given nor received unauthorized aid on this examination.

More information

Cooperative Control and Mobile Sensor Networks

Cooperative Control and Mobile Sensor Networks Cooperative Control and Mobile Sensor Networks Cooperative Control, Part I, A-C Naomi Ehrich Leonard Mechanical and Aerospace Engineering Princeton University and Electrical Systems and Automation University

More information

Math 302 Outcome Statements Winter 2013

Math 302 Outcome Statements Winter 2013 Math 302 Outcome Statements Winter 2013 1 Rectangular Space Coordinates; Vectors in the Three-Dimensional Space (a) Cartesian coordinates of a point (b) sphere (c) symmetry about a point, a line, and a

More information

First-Year Engineering Program. Physics RC Reading Module

First-Year Engineering Program. Physics RC Reading Module Physics RC Reading Module Frictional Force: A Contact Force Friction is caused by the microscopic interactions between the two surfaces. Direction is parallel to the contact surfaces and proportional to

More information

INC 693, 481 Dynamics System and Modelling: Introduction to Modelling Dr.-Ing. Sudchai Boonto Assistant Professor

INC 693, 481 Dynamics System and Modelling: Introduction to Modelling Dr.-Ing. Sudchai Boonto Assistant Professor INC 693, 481 Dynamics System and Modelling: Introduction to Modelling Dr.-Ing. Sudchai Boonto Assistant Professor Department of Control System and Instrumentation Engineering King Mongkut s Unniversity

More information

Module 02 Control Systems Preliminaries, Intro to State Space

Module 02 Control Systems Preliminaries, Intro to State Space Module 02 Control Systems Preliminaries, Intro to State Space Ahmad F. Taha EE 5143: Linear Systems and Control Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/ taha August 28, 2017 Ahmad

More information

Chapter 6. Force and Motion II

Chapter 6. Force and Motion II Chapter 6 Force and Motion II 6 Force and Motion II 2 Announcement: Sample Answer Key 3 4 6-2 Friction Force Question: If the friction were absent, what would happen? Answer: You could not stop without

More information

a = v2 R where R is the curvature radius and v is the car s speed. To provide this acceleration, the car needs static friction force f = ma = mv2

a = v2 R where R is the curvature radius and v is the car s speed. To provide this acceleration, the car needs static friction force f = ma = mv2 PHY 30 K. Solutions for mid-term test #. Problem 1: The forces acting on the car comprise its weight mg, the normal force N from the road that cancels it, and the static friction force f that provides

More information

Simple Car Dynamics. Outline. Claude Lacoursière HPC2N/VRlab, Umeå Universitet, Sweden, May 18, 2005

Simple Car Dynamics. Outline. Claude Lacoursière HPC2N/VRlab, Umeå Universitet, Sweden, May 18, 2005 Simple Car Dynamics Claude Lacoursière HPC2N/VRlab, Umeå Universitet, Sweden, and CMLabs Simulations, Montréal, Canada May 18, 2005 Typeset by FoilTEX May 16th 2005 Outline basics of vehicle dynamics different

More information

Chapter 6 Dynamics I: Motion Along a Line

Chapter 6 Dynamics I: Motion Along a Line Chapter 6 Dynamics I: Motion Along a Line Chapter Goal: To learn how to solve linear force-and-motion problems. Slide 6-2 Chapter 6 Preview Slide 6-3 Chapter 6 Preview Slide 6-4 Chapter 6 Preview Slide

More information

8. THE FARY-MILNOR THEOREM

8. THE FARY-MILNOR THEOREM Math 501 - Differential Geometry Herman Gluck Tuesday April 17, 2012 8. THE FARY-MILNOR THEOREM The curvature of a smooth curve in 3-space is 0 by definition, and its integral w.r.t. arc length, (s) ds,

More information

Alberto Bressan. Department of Mathematics, Penn State University

Alberto Bressan. Department of Mathematics, Penn State University Non-cooperative Differential Games A Homotopy Approach Alberto Bressan Department of Mathematics, Penn State University 1 Differential Games d dt x(t) = G(x(t), u 1(t), u 2 (t)), x(0) = y, u i (t) U i

More information

Delay Coordinate Embedding

Delay Coordinate Embedding Chapter 7 Delay Coordinate Embedding Up to this point, we have known our state space explicitly. But what if we do not know it? How can we then study the dynamics is phase space? A typical case is when

More information

A Robust Controller for Scalar Autonomous Optimal Control Problems

A Robust Controller for Scalar Autonomous Optimal Control Problems A Robust Controller for Scalar Autonomous Optimal Control Problems S. H. Lam 1 Department of Mechanical and Aerospace Engineering Princeton University, Princeton, NJ 08544 lam@princeton.edu Abstract Is

More information

ENGINEERING MECHANICS

ENGINEERING MECHANICS ENGINEERING MECHANICS BASUDEB BHATTACHARYYA Assistant Professor Department of Applied Mechanics Bengal Engineering and Science University Shibpur, Howrah OXJFORD UNIVERSITY PRESS Contents Foreword Preface

More information

EQUATIONS OF MOTION: NORMAL AND TANGENTIAL COORDINATES

EQUATIONS OF MOTION: NORMAL AND TANGENTIAL COORDINATES EQUATIONS OF MOTION: NORMAL AND TANGENTIAL COORDINATES Today s Objectives: Students will be able to: 1. Apply the equation of motion using normal and tangential coordinates. In-Class Activities: Check

More information

Automatic Control II Computer exercise 3. LQG Design

Automatic Control II Computer exercise 3. LQG Design Uppsala University Information Technology Systems and Control HN,FS,KN 2000-10 Last revised by HR August 16, 2017 Automatic Control II Computer exercise 3 LQG Design Preparations: Read Chapters 5 and 9

More information

Lectures 25 & 26: Consensus and vehicular formation problems

Lectures 25 & 26: Consensus and vehicular formation problems EE 8235: Lectures 25 & 26 Lectures 25 & 26: Consensus and vehicular formation problems Consensus Make subsystems (agents, nodes) reach agreement Distributed decision making Vehicular formations How does

More information

Assignment 4.2 Frictional Forces CONCEPTUAL QUESTIONS: 1. What is the SI unit of the coefficient of friction (μ s or μ k )?

Assignment 4.2 Frictional Forces CONCEPTUAL QUESTIONS: 1. What is the SI unit of the coefficient of friction (μ s or μ k )? CONCEPTUAL QUESTIONS: 1. What is the SI unit of the coefficient of friction (μ s or μ k )? 2. Tennis is played on clay, grass, and hard surfaces. Please explain why you think tennis players have or don

More information

A conjecture on sustained oscillations for a closed-loop heat equation

A conjecture on sustained oscillations for a closed-loop heat equation A conjecture on sustained oscillations for a closed-loop heat equation C.I. Byrnes, D.S. Gilliam Abstract The conjecture in this paper represents an initial step aimed toward understanding and shaping

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department 8.01 Physics I Fall Term 2009 Review Module on Solving N equations in N unknowns

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department 8.01 Physics I Fall Term 2009 Review Module on Solving N equations in N unknowns MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department 8.01 Physics I Fall Term 2009 Review Module on Solving N equations in N unknowns Most students first exposure to solving N linear equations in N

More information

P - f = m a x. Now, if the box is already moving, for the frictional force, we use

P - f = m a x. Now, if the box is already moving, for the frictional force, we use Chapter 5 Class Notes This week, we return to forces, and consider forces pointing in different directions. Previously, in Chapter 3, the forces were parallel, but in this chapter the forces can be pointing

More information

Nonlinear Programming Models: A Survey

Nonlinear Programming Models: A Survey Nonlinear Programming Models: A Survey Robert Vanderbei and David Shanno NJ-INFORMS Nov 16, 2000 Rutgers University Operations Research and Financial Engineering, Princeton University http://www.princeton.edu/

More information

CDS 101 Precourse Phase Plane Analysis and Stability

CDS 101 Precourse Phase Plane Analysis and Stability CDS 101 Precourse Phase Plane Analysis and Stability Melvin Leok Control and Dynamical Systems California Institute of Technology Pasadena, CA, 26 September, 2002. mleok@cds.caltech.edu http://www.cds.caltech.edu/

More information

Modeling and Analysis of Dynamic Systems

Modeling and Analysis of Dynamic Systems Modeling and Analysis of Dynamic Systems by Dr. Guillaume Ducard Fall 2016 Institute for Dynamic Systems and Control ETH Zurich, Switzerland based on script from: Prof. Dr. Lino Guzzella 1/33 Outline 1

More information

Homework Solution # 3

Homework Solution # 3 ECSE 644 Optimal Control Feb, 4 Due: Feb 17, 4 (Tuesday) Homework Solution # 3 1 (5%) Consider the discrete nonlinear control system in Homework # For the optimal control and trajectory that you have found

More information