Moment Gyroscope Control Experiment (MAE175a) Prof: Raymond de Callafon

Size: px
Start display at page:

Download "Moment Gyroscope Control Experiment (MAE175a) Prof: Raymond de Callafon"

Transcription

1 MAE175a Moment Gyroscope Control Experiment (MAE175a) Prof: Raymond de Callafon TAs: Jeff Narkis, Gil Collins, this lecture and lab handouts will be available on MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 1 Main Objectives of Laboratory Experiment: modeling and feedback control of a gyroscopic system rotor drum motor with encoder encoder 3 rotor (body D) rotor drum (body C) brake for axis 3 rotor motor 1 with encoder 1 inner gimbal (body B) outer gimbal (body A) encoder 4 and brake for axis 4 MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide

2 Main Objectives of Laboratory Experiment: modeling and feedback control of a gyroscopic system Overview of experiments can also be found in a short (online) movie: download at Control.WMV MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 3 Main Objectives of Laboratory Experiment: modeling and feedback control of a gyroscopic system Ingredients: modeling of dynamic behavior of moment exchange gyro estimation of model parameters using lumped mass system application of control theory for servo/positioning control design, implementation & verification of control sensitivity and error analysis Background Theory: Kinematics and Newton s Law (F = ma), Ordinary Differential Equations (derivation & solutions) Linear System Theory (Laplace transform, Transfer function, Bode plots) Proportional, Integral and Derivative (PID) control analysis and design (root-locus, Nyquist stability criterion) MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 4

3 Outline of this lecture purpose of control & aim of lab experiment hardware description - schematics - hardware in the lab theory on modeling - modeling a moment exchange gyro - step response of a 1st order system - step response of a nd order system outline of laboratory work - estimation of parameters: experiments - validation of model: simulation & experiments - design and implementation of controllers: P- & PD- & PID summary what should be in your report MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 5 Purpose of Control & Some Applications Application of automatic control: to alter dynamic behavior of a system and/or reduce effect of disturbances. industrial processes - thickness control of steel plates in a rolling-mill factory - consistency control in papermaking machines - size and thickness control in glass production processes - control of chemical, distillation or batch reactors electromechanical systems - anti-lock brakes, cruise control and emission control - position control in optical or magnetic storage media - accurate path execution for robotic systems - vibration control in high precision mechanical systems aerospace and aeronautical systems - flight control of pitch, roll and angle-of-attack in aeroplanes - reduction of sound and vibration in helicopters and planes - gyroscope for altitude control of satellites MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 6

4 Aim of Lab Experiment Focus on a (relatively simple) mechanical system of rotating masses/inertias connected in gimbals. Objective is to create a stabilizing feedback system to position mass/inertia at a specified location within a certain time and accuracy. Control is needed to reduce: stabilize rotation (progression of gyro) reduce oscillatory behavior (nutation frequency of gyro) Aim of the experiment: insight in control system principles design and implement control system evaluation of stability & performance error analysis and robustness MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 7 Schematics of Hardware Description block diagram w r K actuator u G y sensor Feedback is essential in control to address stability, disturbance rejection and robustness. For implementation of feedback: gyroscope G is equipped with sensors (to measure signals) and actuators (to activate system) For flexibility of control system K: computer control or digital control and is a combination of: ADC (analogue to digital converter) DSP (digital signal processor) DAC (digital to analogue converter) MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 8

5 Detailed Hardware Description components disturbances encoder amplifier u actuator mech. system sensor y DAC DSP ADC feed forward PC reference plant or partial gyroscope (map from input u to output y) (actuator, inertias, gimbals and encoder) real-time controller (ADC, DSP, DAC, amplifier) Personal Computer (PC) (to run ECP-software and to program DSP) MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 9 Hardware in the Lab moment gyroscope Gimbal: pivoted support that allows rotation about a single axis. Gyroscope: 4 rotating inertias/bodies mounted in gimbals. rotor drum motor with encoder encoder 3 rotor (body D) rotor drum (body C) brake for axis 3 rotor motor 1 with encoder 1 inner gimbal (body B) outer gimbal (body A) encoder 4 and brake for axis 4 MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 10

6 Hardware in the Lab moment gyroscope rotor drum motor with encoder encoder 3 rotor (body D) rotor drum (body C) brake for axis 3 rotor motor 1 with encoder 1 inner gimbal (body B) outer gimbal (body A) encoder 4 with brake for axis 4 Electro-mechanical brakes to restrict axis rotation Inputs: u = voltage V 1 to rotor motor 1 OR u = voltage V to rotor drum motor Outputs: y = angle θ from encoder (rotor drum) OR y = angle θ 3 from encoder 3 (inner gimbal) OR y = angle θ 4 of encoder 4 (outer gimbal) MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 11 Hardware in the Lab real-time control system real time controller and host PC real-time controller: To implement control algorithm and perform digital signal processing. Contains ADC, DSP, DAC & amplifiers. host-pc: To interact with DSP (run ECP software) and MatLab software. MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 1

7 Background theory: modeling a moment exchange gyro motor & encoder encoder 3 rotor (body D) rotor drum (body C) brake for axis 3 motor 1 & encoder 1 inner gimbal (body B) outer gimbal (body A) T 1 θ 1 T 3 θ3 θ T T 4 θ 4x1 x 3 encoder 4 with brake for axis 4 x T i = torque applied to bodies (inputs). θ i = angular position of bodies measured by encoders (outputs). T 1 T gyroscope θ 1 θ θ 3 θ 4 MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 13 Background theory: modeling a moment exchange gyro Important parameters (notation) for intertia modeling: body definition angular position inertia A outer gimbal θ 4 I A B inner gimbal θ 3 I B C rotor drum θ I C D rotor θ 1 I D I b = I b J b K b, b = body A,B,C or D where I b denotes the inertia along the x 1 -axis, J b denotes the inertia along the x axis and K b denotes the inertia along the x 3 axis for b = A,B,C,D. MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 14

8 Background theory: modeling a moment exchange gyro motor & encoder encoder 3 rotor (body D) rotor drum (body C) brake for axis 3 motor 1 & encoder 1 inner gimbal (body B) outer gimbal (body A) T 1 θ 1 T 3 θ3 T θ T 4 θ 4x1 x 3 encoder 4 with brake for axis 4 x Model found via (non-linear) equations of motion are found by combining: Kinematic relations Newton s Law ( T = I θ) Result: relation between angular velocity ω i = θ i and inertia I i for rotor (i = 1), rotor drum (i = ), inner gimbal (i = 3) and outer gimbal (i = 4) MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 15 Background theory: modeling a moment exchange gyro Non-linear equations of motion are of the form: T 1 = f 1 (θ,θ 3,ω,ω 3,ω 4, ω 1, ω 3, ω 4 ) (1) T = f (θ,θ 3,ω 1,ω 3,ω 4, ω 1, ω ) () 0 = f 3 (θ,θ 3,ω 1,ω,ω 3,ω 4, ω 1, ω 3, ω 4 ) (3) 0 = f 4 (θ,θ 3,ω 1,ω,ω 3,ω 4, ω 1, ω, ω 3, ω 4 )(4) and are very hard to work with! Simplification we will use during experiments in the lab: Make approximation by considering small perturbations around: ω 1 = Ω (fixed rotational speed of rotor) θ = θ (initial position of rotor drum) θ 3 = θ 3 (initial position of inner gimbal) for model approximation by linearization. MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 16

9 Background theory: modeling a moment exchange gyro Example of derivation with effects of linearization: J D ω 1 = T 1 J D cos θ ω 3 J D sin θ cos θ 3 ω 4 (1) (I C +I D ) ω = T J D Ωsin θ ω 3 +J D Ωcos θ cos θ 3 ω 4 +(I C +I D )sin θ 3 ω 4 () (J B +J C +J D (J C +J D I D K C )sin θ ) ω 3 = J D cos θ ω 1 +J D Ωsin θ ω J D Ωsin θ sin θ 3 ω 4 sin θ cos θ cos θ 3 (3) (I B +I C K B K C (J C +J D I D K C )sin θ )sin θ 3 ) ω 4 = J D sin θ cos θ 3 ω 1 J D Ωcos θ cos θ 3 ω +(I C +I D )sin θ 3 ω + J D Ωsin θ sin θ 3 ω 3 (J C +J D I D K C )sin θ cos θ cos θ 3 ω 3 (I D +K A +K B +K C +(J C +J D I D K C )sin θ (4) Looks complicated, don t expect you to remember this... However, note: all red terms are constant! Basically: 4 coupled linear equations that relate ω,ω 3,ω 4, ω 1, ω, ω 3, ω 4 to the torques T 1 and T. MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 17 Background theory: modeling a moment exchange gyro The 4 coupled linear equations that relate ω,ω 3,ω 4, ω 1, ω, ω 3, ω 4 to the torques T 1 and T can be written in the matrix notation J ω(t) = Dω(t)+QT(t) where the vector ω(t) of rotational speeds is given by ω 1 (t) [ ] ω ω(t) = (t) T1 (t) and T(t) = ω 3 (t) T (t) ω 4 (t) and where J, D and Q are all previous red terms based on a linearization around ω 1 = Ω (fixed rotational speed of rotor) θ = θ (initial position of rotor drum) θ 3 = θ 3 (initial position of inner gimbal) MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 18

10 Background theory: modeling a moment exchange gyro Choose: ω 1 = Ω (constant rotor speed) θ = θ = 0 (rotor drum perpendicular) θ 3 = θ 3 = 0 (inner gimbal perpendicular) then the 4 coupled linear differential equations reduce to J D ω 1 +J D ω 3 = T 1 (I C +I D ) ω = J D Ωω 4 +T (J B +J C +J D ) ω 3 +J D ω 1 = 0 (I D +K A +K B +K C ) ω 4 = J D Ωω and we see J ω(t) = Dω(t)+QT(t) with J D 0 J D 0 0 (I J = C +I D ) 0 0 J D 0 (J B +J C +J D ) (I D +K A +K B +K C ) J D = D Ω , Q = J D Ω MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 19 Background theory: modeling a moment exchange gyro If we extend the vector ω(t) of angular velocities to also contain (some of) the angular positions x(t) = [ θ (t) θ 3 (t) θ 4 (t) ω 1 (t) ω (t) ω 3 (t) ω 4 (t) ] T then we can write down a so-called state space model ẋ(t) = Ax(t)+Bu(t) y(t) = Cx(t) where u(t) are inputs (torques by motor 1 or motor ) and y(t) are outputs (angular rotation from encoders). In our case matrix C will be C = for θ, θ 3 and θ 4 as output With θ i (t) = ω i (t) and J ω(t) = Dω(t)+QT(t) we see [ ] [ ] 03 3 U A = J 1, B = D J 1, U Q 3 4 = MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 0

11 Background theory: modeling a moment exchange gyro Format of linearized model of gyro in state space form: ẋ(t) = Ax(t)+Bu(t), where y(t) = Cx(t) u(t) = [ T1 (t) T (t) ], y(t) = θ (t) θ 3 (t) θ 4 (t) x(t) = [ θ (t) θ 3 (t) θ 4 (t) ω 1 (t) ω (t) ω 3 (t) ω 4 (t) ] T and, A = J DΩ I D+K A+K B+K C 0 0 C = J DΩ I C+I D, B = J B+J C+J D J D(J B+J C) I C+I D J B+J C MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 1 Background theory: modeling a moment exchange gyro T 1 T gyroscope θ 1 θ θ 3 θ 4 Resulting model of gyro is a linear and multivariable model. Further assumptions for your laboratory experiment: Voltage applied to motor 1 & motor can be used to create external torque T 1 and T. Angular positions θ i = angular position for i =,3,4 can be measured via angular encoders. Further simplifications for your laboratory experiment: Lock one or more axis using brake switches for axis 3 or 4 and virtual brake for axis. Consider only ONE motor to create torgue T j and ONE measurement θ i at the time for simplified Single Input Single Output (SISO) control. MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide

12 Background theory: modeling a moment exchange gyro Impact of this simplification: Using only ONE motor to create a torgue T j and ONE measurement θ i at the time, reduces the multivariable model to a different SISO models. This SISO model can be written as a transfer function G ij (s) MAIN RESULT Let C i be the ith row of the matrix C and let B j denotes the jth column of the matrix B, then G ij (s) = C i (si A) 1 B j where i =,3,4 due to encoder, 3 or 4 and i = 1, due to motor input 1 or. Example: transfer function from motor to encoder 4 ΩJ D θ 4 (s) = G 4 (s)t (s), G 4 (s) = (I C +I D )(I D +K A +K B +K C )s +Ω JD 1 s MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 3 Background theory: modeling a moment exchange gyro motor & encoder encoder 3 rotor (body D) rotor drum (body C) brake for axis 3 motor 1 & encoder 1 inner gimbal (body B) outer gimbal (body A) T 1 θ 1 T 3 θ3 θ T T 4 θ 4x1 x 3 encoder 4 with brake for axis 4 x Configurations of the simplified gyroscope during your lab: Week 1: brake 4 on and (virtual) brake on. Measure θ 3, control motor 1: inertial control. Week : brake 3 on and ω 1 = Ω. Measure θ, control motor : nutation frequency control. Week 3: brake 3 on and ω 1 = Ω. Measure θ 4, control motor : precession control. MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 4

13 Background theory: modeling a moment exchange gyro Models during your lab: Week 1: brake 4 on and (virtual) brake on leads to 1 y 3 (s) = G 31 (s)u 1 (s), with G 31 (s) = (J B +J C )s Week : brake 3 on and ω 1 = Ω: y (s) = G (s)u (s) with I D +K A +K B +K C G (s) = (I C +I D )(I D +K A +K B +K C )s +Ω JD Week 3: brake 3 on and ω 1 = Ω: y 4 (s) = G 4 (s)u (s) with G 4 (s) = ΩJ D (I C +I D )(I D +K A +K B +K C )s 3 +Ω J D s MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 5 Background theory: reduction to 1DOF system Transfer functions G 31 (s), G (s) and G 4 (s) look complicated, but are basically of the form: G 31 (s) = G (s) = G 4 (s) = 1 (J B +J C )s = I D +K A +K B +K C (I C +I D )(I D +K A +K B +K C )s +Ω J D ΩJ D (I C +I D )(I D +K A +K B +K C )s 3 +Ω J D s = = K 0 s +β 0 s K 1 ω n s +βω n s+ω n K ω n s +βω n s+ω n 1 s where we have the (lumped) parameters: scaling K 0, (additional) damping β 0 scaling K 1, damping β and nutation frequency ω n scaling K, damping β and nutation frequency ω n MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 6

14 Background theory: reduction to 1DOF system To be able to do control design, you need to find the parameters of the transfer functions G 31 (s) = G (s) = G 4 (s) = 1 (J B +J C )s = I D +K A +K B +K C (I C +I D )(I D +K A +K B +K C )s +Ω J D ΩJ D (I C +I D )(I D +K A +K B +K C )s 3 +Ω J D s = = K 0 s +β 0 s K 1 ω n s +βω n s+ω n K ω n s +βω n s+ω n 1 s Hence, for modeling we would need to determine inertia of bodies in three different axes... Alternative to determining inertia of bodies in three different axes: determine (lumped) parameters K 0, K 1, K, β 0, β 1 and ω n via dynamic experiments (step response experiments). MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 7 Background theory: step response of a 1st order system Let us consider G 31 (s) = K 0 s +β 0 s as our model for our 1st week laboratory experiment: Instead of using the position θ 3 (t) as output, consider using the velocity ω 3 (t) as output. With ω 3 (t) = d dt θ 3(t) we have ω 3 (s) = sθ 3 (s) using Laplace. Hence, choosing ω 3 (t) as an output y 3 (t) modifies G 31 (s) s to a simple 1st order system: y 3 (s) = K 0 s+β 0 u 1 (s) You can determine (K 0,β 0 ) from a dynamic (step response) experiment... MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 8

15 Background theory: step response of a 1st order system MAIN RESULT Consider a 1st order system with a transfer function G ij (s) = K s+β then a step input u(t) = U, t 0 of size U results in the output response y(t) = K [ β U 1 e βt] Result follows directly from inverse Laplace transform of y(s) = G ij (s) U s = KU s+β 1 s MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 9 Background theory: step response of a 1st order system Typical picture of y(t) = K β U [1 e βt] for a step size of U = 1, scaling K = 5, and (inverse) time constant β = 0.: 30 5 y(t) How about the inverse problem of finding β and K from (velocity) measurement y(t)? MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 30

16 Background theory: step response of a 1st order system With d dt y(t) t=0 and y from velocity step response: 30 5 d dt y(t) t= y We see d dt y(t) t=0 = K Ue βt t=0 = K U and hence ˆK = 1 U d dt y(t) t=0 Knowing that the (velocity) lim t y(t) = y = K β ˆβ = ˆK y U U we find MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 31 Background theory: step response of a nd order system Let us consider K G (s) = 1 ω n s +βω n s+ω n as our model for our nd week laboratory experiment: This is a standard second order system, where ω n = undamped (nutation) frequency, 0 β 1 is the damping ratio, K is the gain. This is similar to the transfer function of a 1DOF mechanical system (with possibly zero damping d)! d F(t) k m x(t) mẍ(t)+dẋ(t)+kx(t) = F(t) x(s) = 1 ms +ds+k F(s) MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 3

17 Background theory: step response of a nd order system In case k > 0 we can write where G(s) = 1 ms +ds+k = Kω n s +βω n s+ω n K := 1 k (compliance), k ω n := m (resonance), β := 1 d (damping ratio) mk Note the resemblance with the nd week lab experiment model I D +K A +K B +K C K 1 ω n G (s) = (I C +I D )(I D +K A +K B +K C )s +Ω JD = s +βω n s+ω n You you determine (K 1,ω n,β) from a dynamic (step response) experiment... MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 33 Background theory: step response of a nd order system MAIN RESULT Consider a nd order system with a transfer function Kω n G ij (s) = s +βω n s+ωn then a step input u(t) = U, t 0 of size U results in the output response y(t) = K U [ 1 e βωnt sin(ω d t+φ) ] where ω d = ω n 1 β damped resonance frequency in rad/s φ = tan 1 1 β phase shift of response in rad β MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 34

18 Background theory: step response of a nd order system DERIVATION for 1DOF (second order) system Compute the dynamic response via inverse Laplace transform! Consider a step input u(t) = U, t 0 of size U. Then u(s) = U s and for the 1DOF system we have Kω n y(s) = G(s)u(s) = s +βω n s+ωn U s and the inverse Laplace transform is given by y(t) = K U [ 1 e βω nt sin(ωd t+φ) ] where ω d = ω n 1 β damped resonance frequency in rad/s φ = tan 1 1 β phase shift of response in rad β MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 35 Background theory: step response of a nd order system Typical picture of y(t) = K U [ 1 e βω nt sin(ω d t+φ) ] for a step size of U = 1, scaling K = 5, undamped resonance frequency ω n = π rad/s and damping ratio β = 0.: y(t) How about the reverse problem of finding ω n, β and K from y(t)? MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 36

19 Background theory: step response of a nd order system With the times t 0, t n and the values y 0, y n and y from (position or velocity) step response: y 0 y n y t 0 t n Allows us to estimate: n ˆω d = π (damped resonance frequency) t n t 0 ( ) βω 1 y0 y n = ln (exponential decay term) t n t 0 y n y where n = number of oscillations between t n and t 0. MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 37 Background theory: step response of a nd order system With the estimates n ˆω d = π t n t 0 ( ) βω 1 y0 y n = ln t n t 0 y n y (damped resonance frequency) (exponential decay term) we can now compute ˆω n and ˆβ: ˆω n = ˆω d +( βω n ) (undamped resonance frequency) βω n ˆβ = ˆω n and K follows from ˆK = y U (damping ratio) (compliance) MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 38

20 Background theory: step response of a nd order system Let us now consider K G 4 (s) = ωn s +βω n s+ωn 1 s as our model for our 3nd week laboratory experiment: Instead of using the position θ 4 (t) as output, consider using the velocity ω 4 (t) as output. With ω 4 (t) = d dt θ 4(t) we have ω 4 (s) = sθ 4 (s) using Laplace. Hence, choosing ω 4 (t) as an output y 4 (t) modifies G 4 (s) s back into a standard nd order system y 4 (s) = K ωn s +βω n s+ωn u (s) NOTE: we can use the same step response results and parameter estimation as we did in our nd week laboratory experiment to estimate K, β and ω n. MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 39 Outline of Lab Work: estimation of model parameters Important observation: gyroscope in 1st week lab experiment appears (dynamically) equivalent to a second order system with G ij (s) = K 0 s +β 0 s and can be reduced to a 1st order system G ij (s) s = K 0 s+β 0 when using the velocity ω 3 (t) as output y 3 (t). Main Idea: Experiment: use a step response laboratory experiment to formulate parameter estimates ˆK and ˆβ 0. Validate model, so that measured step response looks similar to simulated step response MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 40

21 Outline of Lab Work: estimation of model parameters Important observation: gyroscope in nd and 3rd week lab experiment appears (dynamically) equivalent to a single 1DOF with a transfer function where Main Idea: G ij (s) = Kω n s +βω n s+ω n K := scaling, ω n := nutation frequency, β := damping ratio Experiment: use a step response laboratory experiment to formulate parameter estimates ˆK, ˆω n and ˆβ. Validate model, so that measured step response looks similar to simulated step response MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 41 Outline of Lab Work: estimation of model parameters 8 EXAMPLE: step of U = 0.5V on motor Read from plot: t 0 = 0.5, t n = 1.5, y 0 = 7.5, y n = 5.5, y = 5 n ˆω d = π = π t n t 0 ( ) = 4π 1 y0 y βω n = ln = t n t 0 y n y ˆω n = ln y 0 y n y t 0 t n ˆω d +( βωn ) = 16π +ln(10) 1.78 ˆβ = βω n ln(10) ˆω n ˆK = y U = = ( ) = ln(10) Units of ˆK? Does it matter? MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 4

22 Outline of Lab Work: estimation of model parameters Conclusions and summary of work: Verify transfer function you are working with (week 1, or 3) It will be of the form G ij (s) = K ( ) s+β 1 s or G ij (s) = Kω ( ) n 1 s +βω n s+ω n s Estimate the parameters model parameters (K, ω n, β) via the measurement of a step response. Verify the simulation of your step response with the measurement of a step response. Substitute parameter values to create your dynamic model G(s) to be used for control design Luckely, we have Matlab and a matlab script/function called maelab.m to help you with his. Simply enter your models in a file called models.m MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 43 Outline of Lab Work: estimation of model parameters NOTES: During week 1 you DO NOT use a constant rotor speed Ω, as you will be controlling the rotor motor 1 yourself. During week & 3 you DO use a constant rotor speed Ω. For each laboratory section a different constant rotational rotor speed Ω is used. The rotational rotor speed you will be working with is given in the config.txt or info.txt file. Model parameters (K, ω n, β) will depend on physical quantities (mass, inertia, damping, etc.) and need to be determined during laboratory Important in validation of laboratory work: comparison of actual (dynamic) experiments with a (dynamic) simulation of your model G(s) before you start doing control design! With a bad model you cannot do a proper control design... MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 44

23 Outline of Lab Work: model validation w r K actuator u G y sensor Keep in mind: From parameter estimation experiments we know obtain a model specified as a transfer function model y(s) = G ij (s)u(s) G ij (s) = K s +βs or G ij(s) = Kω ( n s +βω n s+ω n where u(s) = is control input (motor Voltage or Force) and y(t) = system output (position in encoder counts) Model G(s) is created automatically for you via maelab.m script file by modifying the models.m file. You are going to use the model G(s) to design a controller K(s), so model G(s) should be validated first! MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 45 1 s ) Outline of Lab Work: model validation What is suitable for model validation: Validate the estimation of each set of model parameters by comparing actual experiments with a simulation of a 1DOF model G ij (s). This can be done with the maelab.m script file. With a bad (unvalidated) model G ij (s) you cannot do a proper model-based control K(s) design! MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 46

24 Outline of Lab Work: design of controllers Each week you construct feedback loop around G(s): r e K(s) u G(s) y Schematic view of closed loop configuration Find a feedback controller K(s) that satisfies: move a mass/inertia to a certain (angular) position as fast as possible limit overshoot during control/positioning to 5% no steady-state error e illustrate disturbance rejection when control is implemented Trade off in design specifications: high speed overshoot overshoot robustness MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 47 Outline of Lab Work: design of controllers Controller configurations that can be implemented: P-control PD-control PID control u(t) = k p e(t), e(t) = r(t) y(t) d u(t) = k p e(t)+k d e(t), e(t) = r(t) y(t) dt t u(t) = k p e(t)+k i e(τ)dτ +k d d e(t), e(t) = r(t) y(t) τ=0 dt MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 48

25 Outline of Lab Work: design of controllers All controllers are implemented in a discrete-time! Let T s = sample time, t = kt s and e(k) = r(k) y(k) we get P-control u(k) = k p e(k) PD-control PID control u(k) = k p e(k)+ k d T s [e(k) e(k 1)] u(k) = k p e(k)+k i T s k l=0 e(l)+ k d T s [e(k) e(k 1)] Templates will be available for you to enter your control parameters k p, k d and/or k i. MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 49 Outline of Lab Work: design of controllers Example of template ;***********define user variables ************** #define kp q ; proportional gain #define kd q3 ; derivative gain #define ki q4 ; integral gain #define Ts q5 ; sampling time for scaling of derivative and integral gain #define kdd q6 ; scaled derivative gain #define kii q7 ; scaled integral gain #define enc_prev q8 ; keep track of previous encoder measurement for derivative #define error q9 ; keep track of error signal for integration ;************Initialize variables**************** Ts=.0044 ; For local use only. Must ALSO be set in SETUP CONTROL ALGORITHM dialog box! kp=0. ; set proportional gain kd=0.0 ; set derivative gain ki=0. ; set integral gain kdd=kd/ts ; divide by Ts here to save real-time processing kii=ki*ts ; multiply by Ts here to save real-time processing error=0 ; initialize integrator error ;PID Controller that uses encoder and motor (nutation control) ;********* real time code which is run every servo period *** begin error=error+cmd_pos-enc_pos control_effort=kp*(cmd_pos-enc_pos)+kii*error+kdd*(enc_prev-enc_pos) enc_prev=enc_pos end MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 50

26 Outline of Lab Work: design of controllers (loop gain) Model G(s) of plant should be used for design of controller K(s)! r e K(s) u G(s) y For design of controller K(s), consider the loop gain: L(s) := K(s)G(s) r e L(s) y loop gain: series connection of K(s) and G(s) Dynamics of loop gain L(s) consists of fixed part G(s) (plant dynamics) and to-be-designed part K(s) (controller) MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 51 Outline of Lab Work: design of controllers (stability) Loop gain L(s) := K(s)G(s) important for: Stability Design specification STABILITY: With closed-loop poles found by those values of s C that satisfy 1+L(s) = 0 all closed-loop poles should have negative real values (lie in the left part of the complex plane) Stability can be checked by: Actually computing the solutions to L(s) = 1 as a function of k p, k i, k d : Root Locus Method See if Nyquist plot of L(s) encircles the point 1 as a function of k p, k i, k d : Nyquist or Frequency Domain Method MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 5

27 Outline of Lab Work: design of controllers (design specs) r e L(s) y Error rejection transfer function: E(s) = 1 1+L(s) (map from r to e) To avoid a steady state error e(t) as t, one specification for the loop gain can be found via the final value theorem. With L(s) = G(s)K(s) we have: lim e(t) = lim t s 1 s 0 1+L(s) r(s) With r(t) = (unit) step input, r(s) = 1 s and lim L(s) = s 0 is needed for zero steady-state behavior! MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 53 Outline of Lab Work: design of controllers (design specs) r e L(s) y Closed loop transfer function: T(s) = L(s) 1+L(s) (map from r to y) To make sure y follows r, we would like to make T(s) = 1 as close as possible. Notice that with Error rejection transfer function: we have E(s) = 1 1+L(s) (map from r to e) T(s)+E(s) = 1 Hence, if you can make E(s) 0 small, then T(s) 1. MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 54

28 Outline of Lab Work: design of controllers (graphical design) Computation of K(s): translate design specifications to L(s), E(S) or T(s) specifications. Use graphical analysis and design utilities (root locus or frequency domain methods) to shape loop gain L(s) and design controller K(s). Root-locus and frequency domain design method is implemented in Matlab via the rltool command. Frequency domain design method has also been implemented in a Matlab script file maelab provided during the lab. Stability via Nyquist criterion (do not encircle point -1) Phase and amplitude margin (stability and robustness) translate to shape Bode plot of loop gain L(s) = G(s)K(s): phase margin: when L(s) = 1, L(s) > π rad amplitude margin: when L(s) = π rad, L(s) < 1 MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 55 Outline of Lab Work: design of controllers (graphical design) Example of figures produced by maelab script file 10 magnitude L(s) L(s)=G(s)*C(s) 4 Nyquist contour L(s) 10 1 magnitude 10 0 imag frequency [Hz] phase L(s) real 10 1 magnitude L(s)/(1+L(s)) [r(s) > y(s)] phase [deg] frequency [Hz] magnitude closed loop stable frequency [Hz] MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 56

29 Outline of Lab Work: design of controllers (general trend) Effects of control parameters: for PD-control and a standard nd order plant model this can be analyzed as follows: ωn s +βω n s+ωn ωn T(s) = L(s) (k p +k d s) 1+L(s) = 1+(k p +k d s) s +βω n s+ωn which yields the closed-loop transfer function T(s) = ω n (k p +k d s) s + β ω n + ω n with ω n = ω n 1+k p and β = β +ω nk d / 1+kp In this case T(s) is also a second order system and with knowledge of the step response, we can conclude that the following influence of the controller parameters: k p speed of response k p damping k p steady-state error k d damping MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 57 Outline of Lab Work: design of controllers (summary) Model G(s) of plant should be used for design of your controller K(s)! Think about the controller needed (P, PD or PID) to stabilize your system Obviously, controllers K(s) for week 1, and 3 are all different and G(s) is different! Keep in mind the requirement of 5% overshoot, and no steady state error, e.g. r(t) = y(t) as t. Use graphical design tools to design your P, PD and PID control: Root-locus and frequency domain design method is implemented in Matlab via the rltool command. Frequency domain design method also implemented in a Matlab script file maelab provided during the lab. MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 58

30 Outline of Lab Work: design of controllers (summary) Consider how the frequency resp of P- and PD- and PID controller modifies the loop gain L(s) = G(s)K(s). Look at asymptotes of Bode plot of controller K(s) = k ds +k p s+k i s Phase and amplitude margin (stability and robustness) translate to shape Bode plot of loop gain L(s) = G(s)K(s): phase margin: when L(s) = 1, L(s) > π rad amplitude margin: when L(s) = π rad, L(s) < 1 Argument and motive your control design in rapport. No trial-and-error control design results are accepted. MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 59 Summary of Lab Work first week Study laboratory handout. Get familiar with the multi-input-multi-output (MIMO) moment gyroscope system, introduction to ECP software used for experiments and controller implementation. Engage axis 4 brake (lock outer gimbal) and virtual brake (lock rotor drum). Design and implement feedback using motor 1 as input and encoder 3 as output. Propose experiments to estimate (unknown) parameters in your model of the system. Design controller based on model in Matlab and implement/verify on actual system. second week Engage axis 3 brake (lock inner gimbal). Bring rotor at fixed speed Ω. Design and implement feedback using motor as input and encoder as output. Propose experiments to estimate (unknown) parameters in your model of the system. Validation of model parameters via comparison of simulation and experiments.design controller based on model in Matlab and implement/verify on actual system. third week Engage axis 3 brake (lock inner gimbal). Bring rotor at fixed speed Ω. Design and implement feedback using motor as input and encoder 4 as output. Propose experiments to estimate (unknown) parameters in your model of the system. Validation of model parameters via comparison of simulation and experiments. Validation of nutation & precession of MIMO moment gyroscope dynamics.design controller based on model in Matlab and implement/verify on actual system. MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 60

31 What should be in your report (1-) Abstract Standalone - make sure it contains clear statements w.r.t motivation, purpose of experiment, main findings (numerical) and conclusions. Introduction Motivation (why are you doing this experiment) Short description of the main engineering discipline (controls) Answer the question: what is the aim of this experiment/report? Theory Feedback system Modeling Parameter estimation Control design MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 61 What should be in your report (-) Experimental Procedure Short description of experiment How are experiments done (detailed enough so someone else could repeat them) Results Estimation of parameters K, ω n and β Validation of models G ij (s) Design of controller K(s) and Implementation Discussion Why are simulation results different from experiments? Could the model be validated? Are designed controller parameters O.K. from model? Conclusions Error Analysis Mean, standard deviation and 99% confidence intervals of estimated parameters ω n, β and K from data MAE175a Control Experiment, Winter 014 R.A. de Callafon Slide 6

Position Control Experiment MAE171a

Position Control Experiment MAE171a Position Control Experiment MAE171a January 11, 014 Prof. R.A. de Callafon, Dept. of MAE, UCSD TAs: Jeff Narkis, email: jnarkis@ucsd.edu Gil Collins, email: gwcollin@ucsd.edu Contents 1 Aim and Procedure

More information

Vibration Analysis Experiment: mode shapes and frequency response of a scaled flexible three-story building & helicopter propeller

Vibration Analysis Experiment: mode shapes and frequency response of a scaled flexible three-story building & helicopter propeller MAE175a Vibration Analysis Experiment: mode shapes and frequency response of a scaled flexible three-story building & helicopter propeller Prof: Raymond de Callafon email: callafon@ucsd.edu TAs: Jeff Narkis,

More information

Linear Circuit Experiment (MAE171a) Prof: Raymond de Callafon

Linear Circuit Experiment (MAE171a) Prof: Raymond de Callafon Linear Circuit Experiment (MAE171a) Prof: Raymond de Callafon email: callafon@ucsd.edu TA: Younghee Han tel. (858) 8221763/8223457, email: y3han@ucsd.edu class information and lab handouts will be available

More information

Control of Manufacturing Processes

Control of Manufacturing Processes Control of Manufacturing Processes Subject 2.830 Spring 2004 Lecture #18 Basic Control Loop Analysis" April 15, 2004 Revisit Temperature Control Problem τ dy dt + y = u τ = time constant = gain y ss =

More information

Dr Ian R. Manchester Dr Ian R. Manchester AMME 3500 : Review

Dr Ian R. Manchester Dr Ian R. Manchester AMME 3500 : Review Week Date Content Notes 1 6 Mar Introduction 2 13 Mar Frequency Domain Modelling 3 20 Mar Transient Performance and the s-plane 4 27 Mar Block Diagrams Assign 1 Due 5 3 Apr Feedback System Characteristics

More information

Controls Problems for Qualifying Exam - Spring 2014

Controls Problems for Qualifying Exam - Spring 2014 Controls Problems for Qualifying Exam - Spring 2014 Problem 1 Consider the system block diagram given in Figure 1. Find the overall transfer function T(s) = C(s)/R(s). Note that this transfer function

More information

Dr Ian R. Manchester

Dr Ian R. Manchester Week Content Notes 1 Introduction 2 Frequency Domain Modelling 3 Transient Performance and the s-plane 4 Block Diagrams 5 Feedback System Characteristics Assign 1 Due 6 Root Locus 7 Root Locus 2 Assign

More information

Positioning Servo Design Example

Positioning Servo Design Example Positioning Servo Design Example 1 Goal. The goal in this design example is to design a control system that will be used in a pick-and-place robot to move the link of a robot between two positions. Usually

More information

FEEDBACK CONTROL SYSTEMS

FEEDBACK CONTROL SYSTEMS FEEDBAC CONTROL SYSTEMS. Control System Design. Open and Closed-Loop Control Systems 3. Why Closed-Loop Control? 4. Case Study --- Speed Control of a DC Motor 5. Steady-State Errors in Unity Feedback Control

More information

(b) A unity feedback system is characterized by the transfer function. Design a suitable compensator to meet the following specifications:

(b) A unity feedback system is characterized by the transfer function. Design a suitable compensator to meet the following specifications: 1. (a) The open loop transfer function of a unity feedback control system is given by G(S) = K/S(1+0.1S)(1+S) (i) Determine the value of K so that the resonance peak M r of the system is equal to 1.4.

More information

Video 5.1 Vijay Kumar and Ani Hsieh

Video 5.1 Vijay Kumar and Ani Hsieh Video 5.1 Vijay Kumar and Ani Hsieh Robo3x-1.1 1 The Purpose of Control Input/Stimulus/ Disturbance System or Plant Output/ Response Understand the Black Box Evaluate the Performance Change the Behavior

More information

Course roadmap. Step response for 2nd-order system. Step response for 2nd-order system

Course roadmap. Step response for 2nd-order system. Step response for 2nd-order system ME45: Control Systems Lecture Time response of nd-order systems Prof. Clar Radcliffe and Prof. Jongeun Choi Department of Mechanical Engineering Michigan State University Modeling Laplace transform Transfer

More information

D(s) G(s) A control system design definition

D(s) G(s) A control system design definition R E Compensation D(s) U Plant G(s) Y Figure 7. A control system design definition x x x 2 x 2 U 2 s s 7 2 Y Figure 7.2 A block diagram representing Eq. (7.) in control form z U 2 s z Y 4 z 2 s z 2 3 Figure

More information

(a) Find the transfer function of the amplifier. Ans.: G(s) =

(a) Find the transfer function of the amplifier. Ans.: G(s) = 126 INTRDUCTIN T CNTR ENGINEERING 10( s 1) (a) Find the transfer function of the amplifier. Ans.: (. 02s 1)(. 001s 1) (b) Find the expected percent overshoot for a step input for the closed-loop system

More information

SRV02-Series Rotary Experiment # 1. Position Control. Student Handout

SRV02-Series Rotary Experiment # 1. Position Control. Student Handout SRV02-Series Rotary Experiment # 1 Position Control Student Handout SRV02-Series Rotary Experiment # 1 Position Control Student Handout 1. Objectives The objective in this experiment is to introduce the

More information

Coupled Drive Apparatus Modelling and Simulation

Coupled Drive Apparatus Modelling and Simulation University of Ljubljana Faculty of Electrical Engineering Victor Centellas Gil Coupled Drive Apparatus Modelling and Simulation Diploma thesis Menthor: prof. dr. Maja Atanasijević-Kunc Ljubljana, 2015

More information

Rotary Motion Servo Plant: SRV02. Rotary Experiment #11: 1-DOF Torsion. 1-DOF Torsion Position Control using QuaRC. Student Manual

Rotary Motion Servo Plant: SRV02. Rotary Experiment #11: 1-DOF Torsion. 1-DOF Torsion Position Control using QuaRC. Student Manual Rotary Motion Servo Plant: SRV02 Rotary Experiment #11: 1-DOF Torsion 1-DOF Torsion Position Control using QuaRC Student Manual Table of Contents 1. INTRODUCTION...1 2. PREREQUISITES...1 3. OVERVIEW OF

More information

Chapter 2. Classical Control System Design. Dutch Institute of Systems and Control

Chapter 2. Classical Control System Design. Dutch Institute of Systems and Control Chapter 2 Classical Control System Design Overview Ch. 2. 2. Classical control system design Introduction Introduction Steady-state Steady-state errors errors Type Type k k systems systems Integral Integral

More information

Automatic Control 2. Loop shaping. Prof. Alberto Bemporad. University of Trento. Academic year

Automatic Control 2. Loop shaping. Prof. Alberto Bemporad. University of Trento. Academic year Automatic Control 2 Loop shaping Prof. Alberto Bemporad University of Trento Academic year 21-211 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 21-211 1 / 39 Feedback

More information

Robust Control 3 The Closed Loop

Robust Control 3 The Closed Loop Robust Control 3 The Closed Loop Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University /2/2002 Outline Closed Loop Transfer Functions Traditional Performance Measures Time

More information

CDS 101/110a: Lecture 8-1 Frequency Domain Design

CDS 101/110a: Lecture 8-1 Frequency Domain Design CDS 11/11a: Lecture 8-1 Frequency Domain Design Richard M. Murray 17 November 28 Goals: Describe canonical control design problem and standard performance measures Show how to use loop shaping to achieve

More information

Prüfung Regelungstechnik I (Control Systems I) Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam!

Prüfung Regelungstechnik I (Control Systems I) Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Prüfung Regelungstechnik I (Control Systems I) Prof. Dr. Lino Guzzella 29. 8. 2 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid

More information

Quanser NI-ELVIS Trainer (QNET) Series: QNET Experiment #02: DC Motor Position Control. DC Motor Control Trainer (DCMCT) Student Manual

Quanser NI-ELVIS Trainer (QNET) Series: QNET Experiment #02: DC Motor Position Control. DC Motor Control Trainer (DCMCT) Student Manual Quanser NI-ELVIS Trainer (QNET) Series: QNET Experiment #02: DC Motor Position Control DC Motor Control Trainer (DCMCT) Student Manual Table of Contents 1 Laboratory Objectives1 2 References1 3 DCMCT Plant

More information

Topic # Feedback Control Systems

Topic # Feedback Control Systems Topic #1 16.31 Feedback Control Systems Motivation Basic Linear System Response Fall 2007 16.31 1 1 16.31: Introduction r(t) e(t) d(t) y(t) G c (s) G(s) u(t) Goal: Design a controller G c (s) so that the

More information

Performance of Feedback Control Systems

Performance of Feedback Control Systems Performance of Feedback Control Systems Design of a PID Controller Transient Response of a Closed Loop System Damping Coefficient, Natural frequency, Settling time and Steady-state Error and Type 0, Type

More information

1 x(k +1)=(Φ LH) x(k) = T 1 x 2 (k) x1 (0) 1 T x 2(0) T x 1 (0) x 2 (0) x(1) = x(2) = x(3) =

1 x(k +1)=(Φ LH) x(k) = T 1 x 2 (k) x1 (0) 1 T x 2(0) T x 1 (0) x 2 (0) x(1) = x(2) = x(3) = 567 This is often referred to as Þnite settling time or deadbeat design because the dynamics will settle in a Þnite number of sample periods. This estimator always drives the error to zero in time 2T or

More information

Outline. Classical Control. Lecture 2

Outline. Classical Control. Lecture 2 Outline Outline Outline Review of Material from Lecture 2 New Stuff - Outline Review of Lecture System Performance Effect of Poles Review of Material from Lecture System Performance Effect of Poles 2 New

More information

Intro to Frequency Domain Design

Intro to Frequency Domain Design Intro to Frequency Domain Design MEM 355 Performance Enhancement of Dynamical Systems Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University Outline Closed Loop Transfer Functions

More information

Inverted Pendulum System

Inverted Pendulum System Introduction Inverted Pendulum System This lab experiment consists of two experimental procedures, each with sub parts. Experiment 1 is used to determine the system parameters needed to implement a controller.

More information

Control of Electromechanical Systems

Control of Electromechanical Systems Control of Electromechanical Systems November 3, 27 Exercise Consider the feedback control scheme of the motor speed ω in Fig., where the torque actuation includes a time constant τ A =. s and a disturbance

More information

Lecture: Sampling. Automatic Control 2. Sampling. Prof. Alberto Bemporad. University of Trento. Academic year

Lecture: Sampling. Automatic Control 2. Sampling. Prof. Alberto Bemporad. University of Trento. Academic year Automatic Control 2 Sampling Prof. Alberto Bemporad University of rento Academic year 2010-2011 Prof. Alberto Bemporad (University of rento) Automatic Control 2 Academic year 2010-2011 1 / 31 ime-discretization

More information

INTRODUCTION TO DIGITAL CONTROL

INTRODUCTION TO DIGITAL CONTROL ECE4540/5540: Digital Control Systems INTRODUCTION TO DIGITAL CONTROL.: Introduction In ECE450/ECE550 Feedback Control Systems, welearnedhow to make an analog controller D(s) to control a linear-time-invariant

More information

Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam!

Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Prüfung Regelungstechnik I (Control Systems I) Prof. Dr. Lino Guzzella 3. 8. 24 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid

More information

Control Systems Design

Control Systems Design ELEC4410 Control Systems Design Lecture 18: State Feedback Tracking and State Estimation Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science Lecture 18:

More information

Mechatronics Engineering. Li Wen

Mechatronics Engineering. Li Wen Mechatronics Engineering Li Wen Bio-inspired robot-dc motor drive Unstable system Mirko Kovac,EPFL Modeling and simulation of the control system Problems 1. Why we establish mathematical model of the control

More information

Analysis and Design of Control Systems in the Time Domain

Analysis and Design of Control Systems in the Time Domain Chapter 6 Analysis and Design of Control Systems in the Time Domain 6. Concepts of feedback control Given a system, we can classify it as an open loop or a closed loop depends on the usage of the feedback.

More information

Lecture 5 Classical Control Overview III. Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore

Lecture 5 Classical Control Overview III. Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Lecture 5 Classical Control Overview III Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore A Fundamental Problem in Control Systems Poles of open

More information

Richiami di Controlli Automatici

Richiami di Controlli Automatici Richiami di Controlli Automatici Gianmaria De Tommasi 1 1 Università degli Studi di Napoli Federico II detommas@unina.it Ottobre 2012 Corsi AnsaldoBreda G. De Tommasi (UNINA) Richiami di Controlli Automatici

More information

Lecture 6 Classical Control Overview IV. Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore

Lecture 6 Classical Control Overview IV. Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Lecture 6 Classical Control Overview IV Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Lead Lag Compensator Design Dr. Radhakant Padhi Asst.

More information

EXPERIMENTALLY DETERMINING THE TRANSFER FUNCTION OF A SPRING- MASS SYSTEM

EXPERIMENTALLY DETERMINING THE TRANSFER FUNCTION OF A SPRING- MASS SYSTEM EXPERIMENTALLY DETERMINING THE TRANSFER FUNCTION OF A SPRING- MASS SYSTEM Lab 8 OBJECTIVES At the conclusion of this experiment, students should be able to: Experimentally determine the best fourth order

More information

Laboratory 11 Control Systems Laboratory ECE3557. State Feedback Controller for Position Control of a Flexible Joint

Laboratory 11 Control Systems Laboratory ECE3557. State Feedback Controller for Position Control of a Flexible Joint Laboratory 11 State Feedback Controller for Position Control of a Flexible Joint 11.1 Objective The objective of this laboratory is to design a full state feedback controller for endpoint position control

More information

Radar Dish. Armature controlled dc motor. Inside. θ r input. Outside. θ D output. θ m. Gearbox. Control Transmitter. Control. θ D.

Radar Dish. Armature controlled dc motor. Inside. θ r input. Outside. θ D output. θ m. Gearbox. Control Transmitter. Control. θ D. Radar Dish ME 304 CONTROL SYSTEMS Mechanical Engineering Department, Middle East Technical University Armature controlled dc motor Outside θ D output Inside θ r input r θ m Gearbox Control Transmitter

More information

AMME3500: System Dynamics & Control

AMME3500: System Dynamics & Control Stefan B. Williams May, 211 AMME35: System Dynamics & Control Assignment 4 Note: This assignment contributes 15% towards your final mark. This assignment is due at 4pm on Monday, May 3 th during Week 13

More information

Time Response Analysis (Part II)

Time Response Analysis (Part II) Time Response Analysis (Part II). A critically damped, continuous-time, second order system, when sampled, will have (in Z domain) (a) A simple pole (b) Double pole on real axis (c) Double pole on imaginary

More information

Part IB Paper 6: Information Engineering LINEAR SYSTEMS AND CONTROL. Glenn Vinnicombe HANDOUT 5. An Introduction to Feedback Control Systems

Part IB Paper 6: Information Engineering LINEAR SYSTEMS AND CONTROL. Glenn Vinnicombe HANDOUT 5. An Introduction to Feedback Control Systems Part IB Paper 6: Information Engineering LINEAR SYSTEMS AND CONTROL Glenn Vinnicombe HANDOUT 5 An Introduction to Feedback Control Systems ē(s) ȳ(s) Σ K(s) G(s) z(s) H(s) z(s) = H(s)G(s)K(s) L(s) ē(s)=

More information

VALLIAMMAI ENGINEERING COLLEGE SRM Nagar, Kattankulathur

VALLIAMMAI ENGINEERING COLLEGE SRM Nagar, Kattankulathur VALLIAMMAI ENGINEERING COLLEGE SRM Nagar, Kattankulathur 603 203. DEPARTMENT OF ELECTRONICS & COMMUNICATION ENGINEERING SUBJECT QUESTION BANK : EC6405 CONTROL SYSTEM ENGINEERING SEM / YEAR: IV / II year

More information

Lecture 5: Frequency domain analysis: Nyquist, Bode Diagrams, second order systems, system types

Lecture 5: Frequency domain analysis: Nyquist, Bode Diagrams, second order systems, system types Lecture 5: Frequency domain analysis: Nyquist, Bode Diagrams, second order systems, system types Venkata Sonti Department of Mechanical Engineering Indian Institute of Science Bangalore, India, 562 This

More information

DC-motor PID control

DC-motor PID control DC-motor PID control This version: November 1, 2017 REGLERTEKNIK Name: P-number: AUTOMATIC LINKÖPING CONTROL Date: Passed: Chapter 1 Introduction The purpose of this lab is to give an introduction to

More information

Laboratory handouts, ME 340

Laboratory handouts, ME 340 Laboratory handouts, ME 340 This document contains summary theory, solved exercises, prelab assignments, lab instructions, and report assignments for Lab 4. 2014-2016 Harry Dankowicz, unless otherwise

More information

Mechatronics Assignment # 1

Mechatronics Assignment # 1 Problem # 1 Consider a closed-loop, rotary, speed-control system with a proportional controller K p, as shown below. The inertia of the rotor is J. The damping coefficient B in mechanical systems is usually

More information

YTÜ Mechanical Engineering Department

YTÜ Mechanical Engineering Department YTÜ Mechanical Engineering Department Lecture of Special Laboratory of Machine Theory, System Dynamics and Control Division Coupled Tank 1 Level Control with using Feedforward PI Controller Lab Date: Lab

More information

Modeling and Experimentation: Compound Pendulum

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

Computer Aided Control Design

Computer Aided Control Design Computer Aided Control Design Project-Lab 3 Automatic Control Basic Course, EL1000/EL1100/EL1120 Revised August 18, 2008 Modified version of laboration developed by Håkan Fortell and Svante Gunnarsson

More information

Control of Manufacturing Processes

Control of Manufacturing Processes Control of Manufacturing Processes Subject 2.830 Spring 2004 Lecture #19 Position Control and Root Locus Analysis" April 22, 2004 The Position Servo Problem, reference position NC Control Robots Injection

More information

Massachusetts Institute of Technology Department of Mechanical Engineering Dynamics and Control II Design Project

Massachusetts Institute of Technology Department of Mechanical Engineering Dynamics and Control II Design Project Massachusetts Institute of Technology Department of Mechanical Engineering.4 Dynamics and Control II Design Project ACTIVE DAMPING OF TALL BUILDING VIBRATIONS: CONTINUED Franz Hover, 5 November 7 Review

More information

Control for. Maarten Steinbuch Dept. Mechanical Engineering Control Systems Technology Group TU/e

Control for. Maarten Steinbuch Dept. Mechanical Engineering Control Systems Technology Group TU/e Control for Maarten Steinbuch Dept. Mechanical Engineering Control Systems Technology Group TU/e Motion Systems m F Introduction Timedomain tuning Frequency domain & stability Filters Feedforward Servo-oriented

More information

CHAPTER 1 Basic Concepts of Control System. CHAPTER 6 Hydraulic Control System

CHAPTER 1 Basic Concepts of Control System. CHAPTER 6 Hydraulic Control System CHAPTER 1 Basic Concepts of Control System 1. What is open loop control systems and closed loop control systems? Compare open loop control system with closed loop control system. Write down major advantages

More information

General procedure for formulation of robot dynamics STEP 1 STEP 3. Module 9 : Robot Dynamics & controls

General procedure for formulation of robot dynamics STEP 1 STEP 3. Module 9 : Robot Dynamics & controls Module 9 : Robot Dynamics & controls Lecture 32 : General procedure for dynamics equation forming and introduction to control Objectives In this course you will learn the following Lagrangian Formulation

More information

MEM 355 Performance Enhancement of Dynamical Systems

MEM 355 Performance Enhancement of Dynamical Systems MEM 355 Performance Enhancement of Dynamical Systems Frequency Domain Design Intro Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University /5/27 Outline Closed Loop Transfer

More information

Singular Value Decomposition Analysis

Singular Value Decomposition Analysis Singular Value Decomposition Analysis Singular Value Decomposition Analysis Introduction Introduce a linear algebra tool: singular values of a matrix Motivation Why do we need singular values in MIMO control

More information

Raktim Bhattacharya. . AERO 632: Design of Advance Flight Control System. Preliminaries

Raktim Bhattacharya. . AERO 632: Design of Advance Flight Control System. Preliminaries . AERO 632: of Advance Flight Control System. Preliminaries Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. Preliminaries Signals & Systems Laplace

More information

Homework 7 - Solutions

Homework 7 - Solutions Homework 7 - Solutions Note: This homework is worth a total of 48 points. 1. Compensators (9 points) For a unity feedback system given below, with G(s) = K s(s + 5)(s + 11) do the following: (c) Find the

More information

Example: DC Motor Speed Modeling

Example: DC Motor Speed Modeling Page 1 of 5 Example: DC Motor Speed Modeling Physical setup and system equations Design requirements MATLAB representation and open-loop response Physical setup and system equations A common actuator in

More information

Dr Ian R. Manchester Dr Ian R. Manchester AMME 3500 : Root Locus

Dr Ian R. Manchester Dr Ian R. Manchester AMME 3500 : Root Locus Week Content Notes 1 Introduction 2 Frequency Domain Modelling 3 Transient Performance and the s-plane 4 Block Diagrams 5 Feedback System Characteristics Assign 1 Due 6 Root Locus 7 Root Locus 2 Assign

More information

Collocated versus non-collocated control [H04Q7]

Collocated versus non-collocated control [H04Q7] Collocated versus non-collocated control [H04Q7] Jan Swevers September 2008 0-0 Contents Some concepts of structural dynamics Collocated versus non-collocated control Summary This lecture is based on parts

More information

Pole placement control: state space and polynomial approaches Lecture 2

Pole placement control: state space and polynomial approaches Lecture 2 : state space and polynomial approaches Lecture 2 : a state O. Sename 1 1 Gipsa-lab, CNRS-INPG, FRANCE Olivier.Sename@gipsa-lab.fr www.gipsa-lab.fr/ o.sename -based November 21, 2017 Outline : a state

More information

Course Summary. The course cannot be summarized in one lecture.

Course Summary. The course cannot be summarized in one lecture. Course Summary Unit 1: Introduction Unit 2: Modeling in the Frequency Domain Unit 3: Time Response Unit 4: Block Diagram Reduction Unit 5: Stability Unit 6: Steady-State Error Unit 7: Root Locus Techniques

More information

ECEn 483 / ME 431 Case Studies. Randal W. Beard Brigham Young University

ECEn 483 / ME 431 Case Studies. Randal W. Beard Brigham Young University ECEn 483 / ME 431 Case Studies Randal W. Beard Brigham Young University Updated: December 2, 2014 ii Contents 1 Single Link Robot Arm 1 2 Pendulum on a Cart 9 3 Satellite Attitude Control 17 4 UUV Roll

More information

(Refer Slide Time: 00:01:30 min)

(Refer Slide Time: 00:01:30 min) Control Engineering Prof. M. Gopal Department of Electrical Engineering Indian Institute of Technology, Delhi Lecture - 3 Introduction to Control Problem (Contd.) Well friends, I have been giving you various

More information

Control Systems I Lecture 10: System Specifications

Control Systems I Lecture 10: System Specifications Control Systems I Lecture 10: System Specifications Readings: Guzzella, Chapter 10 Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich November 24, 2017 E. Frazzoli (ETH) Lecture

More information

Course roadmap. ME451: Control Systems. Example of Laplace transform. Lecture 2 Laplace transform. Laplace transform

Course roadmap. ME451: Control Systems. Example of Laplace transform. Lecture 2 Laplace transform. Laplace transform ME45: Control Systems Lecture 2 Prof. Jongeun Choi Department of Mechanical Engineering Michigan State University Modeling Transfer function Models for systems electrical mechanical electromechanical Block

More information

Lab 6a: Pole Placement for the Inverted Pendulum

Lab 6a: Pole Placement for the Inverted Pendulum Lab 6a: Pole Placement for the Inverted Pendulum Idiot. Above her head was the only stable place in the cosmos, the only refuge from the damnation of the Panta Rei, and she guessed it was the Pendulum

More information

Example: Modeling DC Motor Position Physical Setup System Equations Design Requirements MATLAB Representation and Open-Loop Response

Example: Modeling DC Motor Position Physical Setup System Equations Design Requirements MATLAB Representation and Open-Loop Response Page 1 of 5 Example: Modeling DC Motor Position Physical Setup System Equations Design Requirements MATLAB Representation and Open-Loop Response Physical Setup A common actuator in control systems is the

More information

Lecture 25: Tue Nov 27, 2018

Lecture 25: Tue Nov 27, 2018 Lecture 25: Tue Nov 27, 2018 Reminder: Lab 3 moved to Tuesday Dec 4 Lecture: review time-domain characteristics of 2nd-order systems intro to control: feedback open-loop vs closed-loop control intro to

More information

Table of Laplacetransform

Table of Laplacetransform Appendix Table of Laplacetransform pairs 1(t) f(s) oct), unit impulse at t = 0 a, a constant or step of magnitude a at t = 0 a s t, a ramp function e- at, an exponential function s + a sin wt, a sine fun

More information

APPLICATIONS FOR ROBOTICS

APPLICATIONS FOR ROBOTICS Version: 1 CONTROL APPLICATIONS FOR ROBOTICS TEX d: Feb. 17, 214 PREVIEW We show that the transfer function and conditions of stability for linear systems can be studied using Laplace transforms. Table

More information

Alireza Mousavi Brunel University

Alireza Mousavi Brunel University Alireza Mousavi Brunel University 1 » Control Process» Control Systems Design & Analysis 2 Open-Loop Control: Is normally a simple switch on and switch off process, for example a light in a room is switched

More information

Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam!

Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Prüfung Regelungstechnik I (Control Systems I) Prof. Dr. Lino Guzzella 9. 8. 2 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid -

More information

A brief introduction to robust H control

A brief introduction to robust H control A brief introduction to robust H control Jean-Marc Biannic System Control and Flight Dynamics Department ONERA, Toulouse. http://www.onera.fr/staff/jean-marc-biannic/ http://jm.biannic.free.fr/ European

More information

FATIMA MICHAEL COLLEGE OF ENGINEERING & TECHNOLOGY

FATIMA MICHAEL COLLEGE OF ENGINEERING & TECHNOLOGY FATIMA MICHAEL COLLEGE OF ENGINEERING & TECHNOLOGY Senkottai Village, Madurai Sivagangai Main Road, Madurai - 625 020. An ISO 9001:2008 Certified Institution DEPARTMENT OF ELECTRONICS AND COMMUNICATION

More information

Systems Analysis and Control

Systems Analysis and Control Systems Analysis and Control Matthew M. Peet Illinois Institute of Technology Lecture 8: Response Characteristics Overview In this Lecture, you will learn: Characteristics of the Response Stability Real

More information

ECSE 4962 Control Systems Design. A Brief Tutorial on Control Design

ECSE 4962 Control Systems Design. A Brief Tutorial on Control Design ECSE 4962 Control Systems Design A Brief Tutorial on Control Design Instructor: Professor John T. Wen TA: Ben Potsaid http://www.cat.rpi.edu/~wen/ecse4962s04/ Don t Wait Until The Last Minute! You got

More information

Lecture 12. Upcoming labs: Final Exam on 12/21/2015 (Monday)10:30-12:30

Lecture 12. Upcoming labs: Final Exam on 12/21/2015 (Monday)10:30-12:30 289 Upcoming labs: Lecture 12 Lab 20: Internal model control (finish up) Lab 22: Force or Torque control experiments [Integrative] (2-3 sessions) Final Exam on 12/21/2015 (Monday)10:30-12:30 Today: Recap

More information

School of Engineering Faculty of Built Environment, Engineering, Technology & Design

School of Engineering Faculty of Built Environment, Engineering, Technology & Design Module Name and Code : ENG60803 Real Time Instrumentation Semester and Year : Semester 5/6, Year 3 Lecture Number/ Week : Lecture 3, Week 3 Learning Outcome (s) : LO5 Module Co-ordinator/Tutor : Dr. Phang

More information

The control of a gantry

The control of a gantry The control of a gantry AAE 364L In this experiment we will design a controller for a gantry or crane. Without a controller the pendulum of crane will swing for a long time. The idea is to use control

More information

MEM 355 Performance Enhancement of Dynamical Systems

MEM 355 Performance Enhancement of Dynamical Systems MEM 355 Performance Enhancement of Dynamical Systems Frequency Domain Design Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University 5/8/25 Outline Closed Loop Transfer Functions

More information

EE 468 Control System Laboratory Supporting Materials and Experiments

EE 468 Control System Laboratory Supporting Materials and Experiments EE 468 Control System Laboratory Supporting Materials and Experiments Khosrow Rad Helen Boussalis DEPARTMENT OF ELECTRICAL & COMPUTER ENGINEERING CALIFORNIA STATE UNIVERSITY, LOS ANGELES ECP Systems, Inc

More information

Step Response Analysis. Frequency Response, Relation Between Model Descriptions

Step Response Analysis. Frequency Response, Relation Between Model Descriptions Step Response Analysis. Frequency Response, Relation Between Model Descriptions Automatic Control, Basic Course, Lecture 3 November 9, 27 Lund University, Department of Automatic Control Content. Step

More information

Appendix A: Exercise Problems on Classical Feedback Control Theory (Chaps. 1 and 2)

Appendix A: Exercise Problems on Classical Feedback Control Theory (Chaps. 1 and 2) Appendix A: Exercise Problems on Classical Feedback Control Theory (Chaps. 1 and 2) For all calculations in this book, you can use the MathCad software or any other mathematical software that you are familiar

More information

Aircraft Stability & Control

Aircraft Stability & Control Aircraft Stability & Control Textbook Automatic control of Aircraft and missiles 2 nd Edition by John H Blakelock References Aircraft Dynamics and Automatic Control - McRuler & Ashkenas Aerodynamics, Aeronautics

More information

MAS107 Control Theory Exam Solutions 2008

MAS107 Control Theory Exam Solutions 2008 MAS07 CONTROL THEORY. HOVLAND: EXAM SOLUTION 2008 MAS07 Control Theory Exam Solutions 2008 Geir Hovland, Mechatronics Group, Grimstad, Norway June 30, 2008 C. Repeat question B, but plot the phase curve

More information

LABORATORY INSTRUCTION MANUAL CONTROL SYSTEM I LAB EE 593

LABORATORY INSTRUCTION MANUAL CONTROL SYSTEM I LAB EE 593 LABORATORY INSTRUCTION MANUAL CONTROL SYSTEM I LAB EE 593 ELECTRICAL ENGINEERING DEPARTMENT JIS COLLEGE OF ENGINEERING (AN AUTONOMOUS INSTITUTE) KALYANI, NADIA CONTROL SYSTEM I LAB. MANUAL EE 593 EXPERIMENT

More information

Belt Tension Clamp. Drive Motor. Friction Brake. Load. Encoder 2. Drive. (4000 lines/rev incremental) Encoder 1. (4000 lines/rev incremental)

Belt Tension Clamp. Drive Motor. Friction Brake. Load. Encoder 2. Drive. (4000 lines/rev incremental) Encoder 1. (4000 lines/rev incremental) Industrial Servo System Introduction The first part this lab is to investigate how the dynamic response of a closed-loop system can be used to determine the mass moment of inertia of a model industrial

More information

DC Motor Position: System Modeling

DC Motor Position: System Modeling 1 of 7 01/03/2014 22:07 Tips Effects TIPS ABOUT BASICS INDEX NEXT INTRODUCTION CRUISE CONTROL MOTOR SPEED MOTOR POSITION SUSPENSION INVERTED PENDULUM SYSTEM MODELING ANALYSIS DC Motor Position: System

More information

Introduction to Feedback Control

Introduction to Feedback Control Introduction to Feedback Control Control System Design Why Control? Open-Loop vs Closed-Loop (Feedback) Why Use Feedback Control? Closed-Loop Control System Structure Elements of a Feedback Control System

More information

Frequency methods for the analysis of feedback systems. Lecture 6. Loop analysis of feedback systems. Nyquist approach to study stability

Frequency methods for the analysis of feedback systems. Lecture 6. Loop analysis of feedback systems. Nyquist approach to study stability Lecture 6. Loop analysis of feedback systems 1. Motivation 2. Graphical representation of frequency response: Bode and Nyquist curves 3. Nyquist stability theorem 4. Stability margins Frequency methods

More information

Root Locus. Motivation Sketching Root Locus Examples. School of Mechanical Engineering Purdue University. ME375 Root Locus - 1

Root Locus. Motivation Sketching Root Locus Examples. School of Mechanical Engineering Purdue University. ME375 Root Locus - 1 Root Locus Motivation Sketching Root Locus Examples ME375 Root Locus - 1 Servo Table Example DC Motor Position Control The block diagram for position control of the servo table is given by: D 0.09 Position

More information

2.004 Dynamics and Control II Spring 2008

2.004 Dynamics and Control II Spring 2008 MIT OpenCourseWare http://ocw.mit.edu 2.004 Dynamics and Control II Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Massachusetts Institute

More information

State Feedback Controller for Position Control of a Flexible Link

State Feedback Controller for Position Control of a Flexible Link Laboratory 12 Control Systems Laboratory ECE3557 Laboratory 12 State Feedback Controller for Position Control of a Flexible Link 12.1 Objective The objective of this laboratory is to design a full state

More information

Industrial Servo System

Industrial Servo System Industrial Servo System Introduction The goal of this lab is to investigate how the dynamic response of a closed-loop system can be used to estimate the mass moment of inertia. The investigation will require

More information