Dynamic Effects of Feedback Control!

Similar documents
Module 4: Time Response of discrete time systems Lecture Note 2

Some Basic Information about M-S-D Systems

LAPLACE TRANSFORM AND TRANSFER FUNCTION

Sub Module 2.6. Measurement of transient temperature

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

Math 333 Problem Set #2 Solution 14 February 2003

Lecture 13 RC/RL Circuits, Time Dependent Op Amp Circuits

Embedded Systems and Software. A Simple Introduction to Embedded Control Systems (PID Control)

CHE302 LECTURE VI DYNAMIC BEHAVIORS OF REPRESENTATIVE PROCESSES. Professor Dae Ryook Yang

Chapter 7 Response of First-order RL and RC Circuits

Chapter #1 EEE8013 EEE3001. Linear Controller Design and State Space Analysis

Dead-time Induced Oscillations in Inverter-fed Induction Motor Drives

Open loop vs Closed Loop. Example: Open Loop. Example: Feedforward Control. Advanced Control I

R.#W.#Erickson# Department#of#Electrical,#Computer,#and#Energy#Engineering# University#of#Colorado,#Boulder#

Math Linear Differential Equations

Notes 04 largely plagiarized by %khc

Electrical Circuits. 1. Circuit Laws. Tools Used in Lab 13 Series Circuits Damped Vibrations: Energy Van der Pol Circuit

Determination of the Sampling Period Required for a Fast Dynamic Response of DC-Motors

Chapter 5 Digital PID control algorithm. Hesheng Wang Department of Automation,SJTU 2016,03

Laplace transfom: t-translation rule , Haynes Miller and Jeremy Orloff

Voltage/current relationship Stored Energy. RL / RC circuits Steady State / Transient response Natural / Step response

Laplace Transforms. Examples. Is this equation differential? y 2 2y + 1 = 0, y 2 2y + 1 = 0, (y ) 2 2y + 1 = cos x,

DEPARTMENT OF ELECTRICAL AND ELECTRONIC ENGINEERING EXAMINATIONS 2008

EECS 2602 Winter Laboratory 3 Fourier series, Fourier transform and Bode Plots in MATLAB

First Order RC and RL Transient Circuits

non-linear oscillators

Chapter #1 EEE8013 EEE3001. Linear Controller Design and State Space Analysis

EECE 301 Signals & Systems Prof. Mark Fowler

Problem Set on Differential Equations

Signal and System (Chapter 3. Continuous-Time Systems)

( ) = Q 0. ( ) R = R dq. ( t) = I t

CHBE320 LECTURE IV MATHEMATICAL MODELING OF CHEMICAL PROCESS. Professor Dae Ryook Yang

Basic Circuit Elements Professor J R Lucas November 2001

Outline. Chapter 2: DC & Transient Response. Introduction to CMOS VLSI. DC Response. Transient Response Delay Estimation

2.9 Modeling: Electric Circuits

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

Lecture -14: Chopper fed DC Drives

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

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

Problemas das Aulas Práticas

( ) ( ) ( ) () () Signals And Systems Exam#1. 1. Given x(t) and y(t) below: x(t) y(t) (A) Give the expression of x(t) in terms of step functions.

Math Final Exam Solutions

Linear Time-invariant systems, Convolution, and Cross-correlation

Basic Principles of Sinusoidal Oscillators

2.4 Cuk converter example

RC, RL and RLC circuits

8. Basic RL and RC Circuits

ECE 2100 Circuit Analysis

Linear-Optimal Estimation for Continuous-Time Systems Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 2018

Math Week 15: Section 7.4, mass-spring systems. These are notes for Monday. There will also be course review notes for Tuesday, posted later.

Mon Apr 9 EP 7.6 Convolutions and Laplace transforms. Announcements: Warm-up Exercise:

10. State Space Methods

ECE 2100 Circuit Analysis

Structural Dynamics and Earthquake Engineering

Q1) [20 points] answer for the following questions (ON THIS SHEET):

Program: RFEM 5, RSTAB 8, RF-DYNAM Pro, DYNAM Pro. Category: Isotropic Linear Elasticity, Dynamics, Member

CHE302 LECTURE IV MATHEMATICAL MODELING OF CHEMICAL PROCESS. Professor Dae Ryook Yang

where the coordinate X (t) describes the system motion. X has its origin at the system static equilibrium position (SEP).

Sliding Mode Controller for Unstable Systems

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n

CHAPTER 6: FIRST-ORDER CIRCUITS

MATH 128A, SUMMER 2009, FINAL EXAM SOLUTION

Wall. x(t) f(t) x(t = 0) = x 0, t=0. which describes the motion of the mass in absence of any external forcing.

h[n] is the impulse response of the discrete-time system:

EECE 301 Signals & Systems Prof. Mark Fowler

Lab 10: RC, RL, and RLC Circuits

Chapter 4 DC converter and DC switch

Section 3.8, Mechanical and Electrical Vibrations

EECE 301 Signals & Systems Prof. Mark Fowler

Experimental Buck Converter

Chapter 4 The Fourier Series and Fourier Transform

The Laplace Transform

INDEX. Transient analysis 1 Initial Conditions 1

Introduction to AC Power, RMS RMS. ECE 2210 AC Power p1. Use RMS in power calculations. AC Power P =? DC Power P =. V I = R =. I 2 R. V p.

University of Cyprus Biomedical Imaging and Applied Optics. Appendix. DC Circuits Capacitors and Inductors AC Circuits Operational Amplifiers

Differential Equations

3. Alternating Current

Sterilization D Values

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux

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

( ) = b n ( t) n " (2.111) or a system with many states to be considered, solving these equations isn t. = k U I ( t,t 0 )! ( t 0 ) (2.

Laplace Transform and its Relation to Fourier Transform

Control Systems! Copyright 2017 by Robert Stengel. All rights reserved. For educational use only.

Chapter 2: Principles of steady-state converter analysis

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

Nature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time.

Lecture 1 Overview. course mechanics. outline & topics. what is a linear dynamical system? why study linear systems? some examples

di Bernardo, M. (1995). A purely adaptive controller to synchronize and control chaotic systems.

Load Calculations Heat Balance Method - Theory. Prof. Jeffrey D. Spitler School of Mechanical and Aerospace Engineering, Oklahoma State University

Math 334 Fall 2011 Homework 11 Solutions

System Processes input signal (excitation) and produces output signal (response)

3 General forced response

Solutions of Sample Problems for Third In-Class Exam Math 246, Spring 2011, Professor David Levermore

6.01: Introduction to EECS I Lecture 8 March 29, 2011

Design of Controller for Robot Position Control

Finite element method for structural dynamic and stability analyses

EEEB113 CIRCUIT ANALYSIS I

Dynamic Analysis of Damped Driven Pendulum using Laplace Transform Method

SOLUTIONS TO ECE 3084

Summary of shear rate kinematics (part 1)

Transcription:

Dynamic Effecs of Feedback Conrol! Rober Sengel! Roboics and Inelligen Sysems MAE 345, Princeon Universiy, 2017 Inner, Middle, and Ouer Feedback Conrol Loops Sep Response of Linear, Time- Invarian (LTI) Sysems Posiion and Rae Conrol Transien and Seady-Sae Response o Sinusoidal Inpus Copyrigh 2017 by Rober Sengel. All righs reserved. For educaional use only. hp://www.princeon.edu/~sengel/mae345.hml 1 Ouer-o-Inner-Loop Conrol Hierarchy Inner Loop! Small Ampliude! Fas Response! High Bandwidh Middle Loop! Moderae Ampliude! Medium Response! Moderae Bandwidh Ouer Loop! Large Ampliude! Slow Response! Low Bandwidh Feedback! Error beween command and feedback signal drives nex inner-mos loop 2

Naural Feedback Conrol Inner Loop Chicken Head Conrol - 1 hp://www.youube.com/wach?v=_dplkfpowcc Middle Loop Hovering Red-Tail Hawks hp://www.youube.com/wach?v=- VPVZMSEvwU Ouer Loop Osprey Diving for Fish hp://www.youube.com/wach? v=qrgpl9-n6jy 3 Ouer-o-Inner-Loop Conrol Hierarchy of an Indusrial Robo Inner Loop! Focus on conrol of individual joins Middle Loop! Focus on operaion of he robo Ouer Loop! Focus on goals for robo operaion 4

Inner-Loop Feedback Conrol Feedback conrol design mus accoun for acuaor-sysem-sensor dynamics Single-Inpu/Single-Oupu Example, wih forward and feedback conrol logic ( compensaion ) 5 Thermosaic Temperaure Conrol Dynamics! Delays! Dead Zones! Sauraion! Coupling Srucure! Layou! Insulaion! Circulaion! Muliple Spaces Exernal Effecs! Solar Radiaion! Air Temperaure! Wind! Rain, Humidiy... all conrolled by a simple (bu nonlinear) on/off swich 6

Thermosa Conrol Logic Conrol Law [i.e., logic ha drives he conrol variable, u] e =! y = u c! u b % u = % < Thermosa > 1(on), e > 0 0 (off ), e " 0 : Desired oupu variable (command) y: Acual oupu u: Conrol variable (forcing funcion) e: Conrol error 7 Thermosa Conrol Logic " u = % 1 (on), e > 0 0 (off ), e! 0...bu conrol signal would chaer wih slighes change of emperaure Soluion: Inroduce lag o slow he swiching cycle, e.g., hyseresis % u = % 1 (on), e! T > 0 0 (off ), e + T " 0 8

Thermosa Conrol Logic wih Hyseresis Hyseresis delays he response Sysem responds wih a limi cycle Cooling conrol is similar wih sign reversal 9 Speed Conrol of Direc-Curren Moor Angular Rae Linear Feedback Conrol Law (c = Conrol Gain) u = ce where e =! y How would y be measured? 10

Characerisics of he Model 1/J Simplified Dynamic Model! Roary ineria, J, is he sum of moor and load inerias! Inernal damping negleced! Oupu speed, y, rad/s, is an inegral of he conrol inpu, u! Moor conrol orque is proporional o u! Desired speed,, rad/s, is consan 11 c Angular Rae Model of Dynamics and Speed Conrol Firs-order LTI ordinary differenial equaion dy d = 1 J u = c J e = c [ J! y ], y( 0) given Inegral of he equaion, wih y(0) = 0 y = 1 J! ud = c J 0 = " c J! ed = c J 0![ y ]d + c J 0! 0! 0 [ ]d [ " y ]d Posiive inegraion of Negaive feedback of y 12

Angular Rae Soluion of he inegral Sep Response of Speed Conroller Sep inpu : " 0, < 0 y C = 1,! 0 % ( y = * 1! exp! )* " c % J ' + -,- = y () c 1! exp. +, = y ( 1! exp! / + c )*,- where!! = c/j = eigenvalue or roo of he sysem (rad/sec)! " = J/c = ime consan of he response (sec) Wha does he shaf angle response look like? 13 Angle Conrol of Direc-Curren Moor Angular Posiion Simplified Dynamic Model! Roary ineria, J, is he sum of moor and load inerias! Oupu angle, y, is a double inegral of he conrol, u! Desired angle,, is consan Feedback Conrol Law u = ce where e =! y How would y be measured? 14

Model of Dynamics and Angle Conrol Angular Posiion 2 nd -order, linear, ime-invarian ordinary differenial equaion d 2 y d 2 =!!y = 1 J u( ) = c J e( ) = c [ J y! ] Oupu angle, y, as a funcion of ime y = c J " "[! y ] d d 0 0 [ ] =! c y J " " d 2 + c J 0 0 " [ ] 0 0 " d 2 15 Model of Dynamics and Angle Conrol Corresponding se of 1 s -order equaions, wih! Angle: x 1 = y! Angular rae: x 2 = dy/d!x 2 = u J!x 1 = x 2 = c [ J! y ] = c [ J! x 1 ] Angular Posiion 16

Open-loop dynamic equaion Sae-Space Model of he DC Moor! "!x 1!x 2 % =! 0 1 " 0 0! % " x 1 x 2 % +! 0 " 1 / J u % Feedback conrol law u = c[! y 1 ]= c! x 1 [ ] Closed-loop dynamic equaion! "!x 1!x 2 % =! 0 1 " 'c / J 0! % " x 1 x 2 % +! 0 " c / J % 17 Sep Response wih Angle Feedback! "!x 1!x 2 % =! 0 1 " 'c / J 0! % " x 1 x 2 % +! 0 " c / J % % Sep Response of Undamped Angle Conrol F1 = [0 1;-1 0]; G1 = [0;1]; F2 = [0 1;-0.5 0]; G2 = [0;0.5]; F3 = [0 1;-0.25 0]; G3 = [0;0.25]; Hx = [1 0;0 1]; c/j = 1, 0.5, and 0.25 Sys1 Sys2 Sys3 = ss(f1,g1,hx,0); = ss(f2,g2,hx,0); = ss(f3,g3,hx,0); sep(sys1,sys2,sys3) 18

Conrol law wih rae feedback Wha Wen Wrong? No damping! Soluion: Add rae feedback in he conrol law u = c 1 [! y 1 ]! c 2 y 2 Closed-loop dynamic equaion! "!x 1!x 2 % =! 0 1 'c 1 / J 'c 2 / J "! % " x 1 x 2 % +! 0 c 1 / J " % 19 Alernaive Implemenaions of Rae Feedback dy u = c 1 [! y 1 ]! c 2 y 2 = c 1 [! y 1 ]! c 1 2 d One inpu, wo oupus One inpu, one oupu 20

Sep Response wih Angle and Rae Feedback c 1 /J = 1 c 2 /J = 0, 1.414, 2.828 % Sep Response of Damped Angle Conrol F1 = [0 1;-1 0]; G1 = [0;1]; F1a = [0 1;-1-1.414]; F1b = [0 1;-1-2.828]; Hx = [1 0;0 1]; Sys1 Sys2 Sys3 = ss(f1,g1,hx,0); = ss(f1a,g1,hx,0); = ss(f1b,g1,hx,0); sep(sys1,sys2,sys3) 21 LTI Model wih Feedback Conrol Command inpu, u c, has dimensions of u u = u c! Cy x = F x + Gu + Lw y = H x x + H u u 22

LTI Conrol wih Forward-Loop Gain u = C[! y ] x = F x + Gu + Lw y = H x x + H u u Wih C c = C, command inpu,, has dimensions of y 23 Effec of Feedback Conrol on he LTI Model [ ] [ ]!x = Fx + Gu = Fx + G u c! Cy { } = F open loop x + G u c! C H x x = [ F! GCH x ]x + Gu c! F closed loop x + Gu c Feedback modifies he sabiliy marix of he closed-loop sysem Convergence or divergence Envelope of ransien response 24

LTI Model wih Feedback Conrol and Command Gain Command inpu,, is shaped by C c u = u c! Cy = C c! Cy 25 Effec of Command Gain on LTI Model { } [ ]!x = Fx + Gu = Fx + G C c! Cy { } = Fx + G C c! C H x x = [ F! GCH x ]x + GC c Seady-sae response of he sysem!x = 0 x * =![ F! GCH x ]!1 GC c * Command gain adjuss he seady-sae response Has no effec on he sabiliy of he sysem 26

= sin (! ) = sin( 6.28 ), deg Response o Sine Wave Inpu wih Angle Feedback: No Damping c 1 /J = 1; c 2 /J = 0 Why are here 2 oscillaions in he oupu?! Undamped ransien response o he inpu! Long-erm dynamic response o he inpu Sysem has a naural frequency of oscillaion, n Long-erm response o a sine wave is a sine wave 27 Response o Sine Wave Inpu wih Rae Damping = sin (! ) = sin( 6.28 ), deg c 1 /J = 1; c 2 /J = 1.414 c 1 /J = 1; c 2 /J = 2.828 Wih damping, ransien response decays In his case, damping has negligible effec on long-erm response 28

Sysem Dynamics Produces Differences in Ampliude and Phase Angle of Inpu and Oupu Ampliude Raio (AR) = y Oupu Peak y Inpu Peak ( ) Phase Angle = 360! Inpu Peak Oupu Peak, deg Period of Inpu Ampliude raio and phase angle characerize he sysem model 29 Effec of Inpu Frequency on Oupu Ampliude and Phase Angle = sin( / 6.28), deg Wih low inpu frequency, inpu and oupu ampliudes are abou he same Oupu angle oscillaion lags inpu by a few degrees Rae oscillaion leads angle oscillaion by ~90 deg c 1 /J = 1; c 2 /J = 1.414 30

A Higher Frequency, Oupu Ampliude Decreases, Phase Angle Lag Increases = sin( ), deg c 1 /J = 1; c 2 /J = 1.414 31 = sin( 6.28 ), deg A Even Higher Frequency, Ampliude Raio Decreases c 1 /J = 1; c 2 /J = 1.414 32

Frequency Response of he DC Moor wih Feedback Conrol No Damping Damping More Damping!! Long-erm response o sinusoidal inpus over range of frequencies!! Deermine experimenally!! or from he ransfer funcion!! Frequency response depiced by Bode Plo!! Transfer funcion!! Laplace ransform of sysem!! TBD 33 Nex Time:! Conrol Sysems! 34