An Introduction to Differential Flatness

Similar documents
An Introduction to Differential Flatness. Will Burgoyne May 7, 2007 ME598 Geometric Mechanics

Robotics, Geometry and Control - A Preview

Transverse Linearization for Controlled Mechanical Systems with Several Passive Degrees of Freedom (Application to Orbital Stabilization)

Dynamical Systems & Lyapunov Stability

Eigenvalues and Eigenvectors

Orbital Motion in Schwarzschild Geometry

Geometry in a Fréchet Context: A Projective Limit Approach

On the convergence of normal forms for analytic control systems

Lecture 37: Principal Axes, Translations, and Eulerian Angles

EXPERIMENTAL COMPARISON OF TRAJECTORY TRACKERS FOR A CAR WITH TRAILERS

MATH 320, WEEK 2: Slope Fields, Uniqueness of Solutions, Initial Value Problems, Separable Equations

Curves in the configuration space Q or in the velocity phase space Ω satisfying the Euler-Lagrange (EL) equations,

Classical Mechanics Review (Louisiana State University Qualifier Exam)

for changing independent variables. Most simply for a function f(x) the Legendre transformation f(x) B(s) takes the form B(s) = xs f(x) with s = df

Duality and dynamics in Hamilton-Jacobi theory for fully convex problems of control

The Complex Numbers c ). (1.1)

Physics 235 Chapter 7. Chapter 7 Hamilton's Principle - Lagrangian and Hamiltonian Dynamics

Math 164 (Lec 1): Optimization Instructor: Alpár R. Mészáros

Introduction Knot Theory Nonlinear Dynamics Topology in Chaos Open Questions Summary. Topology in Chaos

Élie Cartan s Theory of Moving Frames

Matlab-Based Tools for Analysis and Control of Inverted Pendula Systems

Parametric Equations and Polar Coordinates

How curvature shapes space

Daba Meshesha Gusu and O.Chandra Sekhara Reddy 1

Unit IV: Introduction to Vector Analysis

MATH 415, WEEKS 7 & 8: Conservative and Hamiltonian Systems, Non-linear Pendulum

1995 ASME Int l Mech Eng Congress and Expo San Francisco, CA November 12-17, 1995

Lecture 38: Equations of Rigid-Body Motion

THE RICCI BRACKET FLOW FOR SOLVABLE LIE GROUPS

Lyapunov functions on nonlinear spaces S 1 R R + V =(1 cos ) Constructing Lyapunov functions: a personal journey

Vectors. January 13, 2013

Tensor Analysis in Euclidean Space

Stabilization of a 3D Rigid Pendulum

Equivalence, Invariants, and Symmetry

Epsilon Delta proofs

HYPERBOLIC DYNAMICAL SYSTEMS AND THE NONCOMMUTATIVE INTEGRATION THEORY OF CONNES ABSTRACT

3 Parallel transport and geodesics

Lecture 38: Equations of Rigid-Body Motion

Intrinsic Vector-Valued Symmetric Form for Simple Mechanical Control Systems in the Nonzero Velocity Setting

Mathematics for Control Theory

Special Conformal Invariance

HAMILTON S PRINCIPLE

DIFFERENTIAL FLATNESS AND ABSOLUTE EQUIVALENCE OF NONLINEAR CONTROL SYSTEMS

0. Introduction 1 0. INTRODUCTION

Chapter 6. Differentially Flat Systems

Motion Planning of Discrete time Nonholonomic Systems with Difference Equation Constraints

The Liapunov Method for Determining Stability (DRAFT)

= 0. = q i., q i = E

BACKGROUND IN SYMPLECTIC GEOMETRY

EN Nonlinear Control and Planning in Robotics Lecture 2: System Models January 28, 2015

In the presence of viscous damping, a more generalized form of the Lagrange s equation of motion can be written as

Ahlswede Khachatrian Theorems: Weighted, Infinite, and Hamming

General Relativity and Cosmology Mock exam

MATH 12 CLASS 4 NOTES, SEP

MATH 353 LECTURE NOTES: WEEK 1 FIRST ORDER ODES

MATH 320, WEEK 6: Linear Systems, Gaussian Elimination, Coefficient Matrices

= w. These evolve with time yielding the

Functions and Their Graphs c 2002 Donald Kreider and Dwight Lahr

Hodge structures from differential equations

Équivalence et linéarisation des systèmes de contrôle

Math 396. Bijectivity vs. isomorphism

Finding Brake Orbits in the (Isosceles) Three- Body Problem. Rick Moeckel -- University of Minnesota

Stabilization of Relative Equilibria for Underactuated Systems on Riemannian Manifolds

Stabilization of Angular Velocity of Asymmetrical Rigid Body. Using Two Constant Torques

LOG-LINEAR ODES AND APPLICATIONS TO THE RICCI FLOW FOR HOMOGENEOUS SPACES

Quantum Mechanics without Spacetime II

Lecture 1, August 21, 2017

The Geometry of Euler s equation. Introduction

M3-4-5 A16 Notes for Geometric Mechanics: Oct Nov 2011

Chapter 1. Vectors, Matrices, and Linear Spaces

Appendix E : Note on regular curves in Euclidean spaces

Nonlinear Tracking Control of Underactuated Surface Vessel

q HYBRID CONTROL FOR BALANCE 0.5 Position: q (radian) q Time: t (seconds) q1 err (radian)

Designing Information Devices and Systems I Fall 2015 Anant Sahai, Ali Niknejad Homework 2. This homework is due September 14, 2015, at Noon.

Algebraic. techniques1

On the Linearization of Second-Order Dif ferential and Dif ference Equations

Physics 212: Statistical mechanics II Lecture XI

5 Topological defects and textures in ordered media

Lecture 13: Forces in the Lagrangian Approach

Knots and Mirror Symmetry. Mina Aganagic UC Berkeley

Energy-based Swing-up of the Acrobot and Time-optimal Motion

Symmetry and Simplicity

Chapter 23. Predicting Chaos The Shift Map and Symbolic Dynamics

Noether s Theorem: Uses and Abuses

Controllability Analysis of A Two Degree of Freedom Nonlinear Attitude Control System

M2A2 Problem Sheet 3 - Hamiltonian Mechanics

Hybrid Routhian Reduction of Lagrangian Hybrid Systems

Function approximation

Characterizing Algebraic Invariants by Differential Radical Invariants

SPECTRAL VARIATIONAL INTEGRATORS

Chapter 8. Rigid transformations

Utility of irreducible group representations in differential equations (II)

8. Prime Factorization and Primary Decompositions

Empirical Likelihood

SECTION C: CONTINUOUS OPTIMISATION LECTURE 9: FIRST ORDER OPTIMALITY CONDITIONS FOR CONSTRAINED NONLINEAR PROGRAMMING

Leaving Cert Applied Maths: Some Notes

3.C. Further comments on axioms for geometry

1 Matrices and matrix algebra

Decoherence and the Classical Limit

Multivariate Differentiation 1

Transcription:

An Introduction to Differential Flatness William Burgoyne ME598 Geometric Mechanics May 7, 007

INTRODUCTION Systems of ordinary differential equations (ODEs) can be divided into two classes, determined and underdetermined. The latter class can be defined as a system of m ODEs with n dependent variables where n > m. Systems of this sort are typically encountered by control engineers. They attempt to define n m of the dependent variables as functions of time in order force the system to act in a specified way. When the system of ODEs is linear, a number of methods exist to easily design control for the system. However, when the system of ODEs is nonlinear, control design is much more difficult and fewer design methods exist. One method that has been developed to design control somewhat easily for nonlinear systems exploits a property of such systems known as differential flatness. A differentially flat system is a system of nonlinear ODEs whose solution curves are in a smooth one-one correspondence with arbitrary curves in a space whose dimension equals n m. Hence generating solution trajectories for the system is simplified by finding curves in a lower dimension that satisfy conditions at the endpoints. In this paper I provide an example of a differentially flat system and rigorously define differential flatness. I relate differential flatness to feedback linearization. I conclude with some general results of differential flatness and some remarks about finding flat outputs for Lagrangian systems. MOTIVATIONAL EXAMPLE Suppose it is desired to control the rigid body in the vertical plane as pictured in Figure. The set of ODEs defining this system is mx && sinθ + mx && cosθ mg cosθ = F mx && cosθ mx && sinθ + mg sinθ = F I&& θ = F R which is nonlinear. The system is underdetermined by two equations. The states of the system are x, x, θ and the inputs are F, F. One method to control this system is to assign initial values to x, x, θ and their derivatives and assign F as arbitrary functions of time. The system becomes determined and evolves from a desired initial state to some final state. However, it would be difficult to determine the proper functions for F that result in the evolution of the system to a desired final value. Another option to control the system is to simplify the problem by eliminating F from the first and third ODEs. This results in a single ODE in terms of the states only mx & sinθ + mx && cosθ mg cosθ + ( I / R) & θ = 0

Figure Rigid body in the vertical plane. which is also a nonlinear ODE. The system is still underdetermined by two equations. In this case, initial conditions need to be assigned for only one of the states and its derivative. The other two states are assigned as some desired functions of time. Once the solution for this ODE is obtained, then it can be used to determine F. While this method ensures that two of the states evolve in a desired way, the third state may not. Again, it would be difficult to determine the proper functions for the two states that result in the evolution of the third state to a desired final value. One might ask if there exists a combination of two dependent variables such that defining them as arbitrary functions of time would determine the system and cause the system to evolve as desired. The answer is no. However, one can form a pair of variables that are functions of the dependent variables and that determine the system without the need of defining initial conditions. One such combination is y = x ( I / mr)cosθ + ( I / mr)sinθ where y and y are called flat outputs. These can be thought of as components of the map from the system space to the smaller dimension space. In other words, defining two new variables, which are functions of the dependent variables, as arbitrary functions of time determines the solutions for the system of five dependent variables. (Defining y and y so that the system meets control requirements is discussed in a following section.) DEFINITION OF DIFFERENTIAL FLATNESS y = x A system is said to be differentially flat or simply flat if there exist variables y,..., y n m given by an equation of the form k y = h( t, x, x, x,..., x ) such that the original variables x may be recovered from y (locally) by an equation of the form l x = g( t, y, y, y,..., y ) The variables y,..., y n m are referred to as flat outputs. When y is only a function of time and the states the system is said to be configuration flat, which is the case for the preceding example.

TRAJECTORY GENERATION Returning to the example, generating solutions that take the system from a desired initial condition to a desired final value is relatively easy. Given the initial conditions and final values for the states, their derivatives, and the inputs, the values for y, y, and their derivatives at the initial and final times can be determined. Then, one simply needs to generate functions of time that satisfy these initial and final conditions. It is interesting to note that connecting any starting and ending point in the system is possible. Hence, flat systems are controllable. CONTROL SYSTEMS AND DIFFERENTIAL FLATNESS When designing control for nonlinear systems, occasionally it is possible to transform the states in such a way that the control system can be expressed in a linear form. These systems are said to be feedback linearizable. It has been proven that a specific version of feedback linearizability is equivalent to differential flatness. Thus, there are potentially two ways to approach control design for certain nonlinear systems. One could use the feedback linearizability property of the system or to use the flatness property of the system. RESULTS FROM DIFFERENTIAL FLATNESS There do not exist general conditions that classify a system as flat or not. Also, there do not exist general conditions for determining the flat outputs of a differentially flat system. However, some specific properties do exist for flat systems. All controllable linear or state feedback linearizable systems are differentially flat. Differentially flat systems permit a transformation of the system such that the EOMs for the flat output variables become trivial. A performance index for a system can be approximated and transformed into a function of the flat outputs. Lagrangian systems have been studied in detail for flatness and catalog exists detailing which types of Lagrangian systems are flat. The results of this study have revealed properties of these systems that indicate flatness. A common feature of rigid bodies with S symmetry is that they evolve on manifolds which consist of, or include, a Lie group as part of the configuration space. Another common feature of rigid bodies with S symmetry is that the flat outputs break symmetries because if the output is invariant with respect to some action, then the motion of the system can only be retrieved up to an initial choice of the group variable. Another Lagrangian system that is flat is the system with three rigid bodies connected in a chain. This particular system and others are flat if certain parameters of the system are defined appropriately. A number of other mechanical systems are differentially flat and papers have been written explaining what the flat outputs are for those systems.

CONCLUSIONS Nonlinear systems that have the property of differential flatness are relatively easy to control. Perhaps the most difficult part or using differential flatness is determining the flat outputs. Structural information contained in Lagrange s equations can be exploited to determine is a system is flat and how to find the flat outputs. While finding the flat outputs may be difficult, geometry can be used to aid in this task. BIBLIOGRAPHY Murray, R. M., et al., Differential flatness of mechanical control systems: A catalog of prototype systems AMSE International Mechanical Engineering Congress and Expo, Nov 995. Rathinam, M., Differentially flat nonlinear control systems Thesis, 997. Rathinam, M. and Murray, R. M., Configuration flatness of Lagrangian systems underactuated by one control SIAM Journal on Control and Optimization, 36, 64, 998. Rathinam, M. and Sluis, W. M., A test for differential flatness by reduction to single input systems Technical Memorandum No. CIT-CDS 95-08, June, 995. Van Nieuwstadt, M., et al., Differential flatness and absolute equivalence of nonlinear contol systems SIAM Journal of Control and Optimization, 36, 5, 998.