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 design of mechanical systems 3
. Time Domain Tuning 4
Servo force? F s Mass M Disturbance F d 5
Mechanical solution: k d F Disturbance F d Mass M Forcespring-damper: F = k d & Eigenfrequ ency: f = π k M 6
Servo analogon: F servo Mass M ServoForce: F s = k p k v & Eigenfrequ ency: f = π k p M k k p v : : servo stiffness servo damping 7
Eample: Disturbance F d F s Mass M Slide: mass = 5 kg Required accuracy 0 µm at all times Disturbance (f.e. friction) = 3 N. Required servo stiffness?. Eigenfrequency? 8
9 h or s ) ( ) ( : -damper Spring h d h k F & & + = How to move to / follow a setpoint: ) ( ) ( : Controller k k F s v s p s & & + =
s F d F disturbance controller process s + - e k p + k v d dt + + F servo F Mass K p /k v -controller or PD-controller 0
s+ - e controller k k p v =? =? F servo process k k p v error Trade off stability
Concluding remarks time domain tuning A control system, consisting of only a single mass m and a k p /k v controller (as depicted below), is always stable. k p will act as a spring; k v will act as a damper As a result of this: when a control system is unstable, it cannot be a pure single mass + k p /k v controller (With positive parameters m, k p and k v ) s k p M k v
Setpoints: s What should s look like as a function of time, when moving the mass? (first order, second order, third order,.?) 3
Apply a force F (step profile): F ( t) = M & ( t) F M (t) is second order, when F constant Second order profile requires following information: - maimum acceleration - maimum velocity - travel distance 4
Eample Pos Vel Acc ma ma = = = 0π 500π π 6.3rad 63rad.6e3rad / sec / sec 5
3 Frequency domain Time domain: Monday and Thursday at :0 Frequency domain: twice a week 6
F M amplitude (log) H = F - phase H -80 o frequency (log) 7
going from Time-domain to the Frequency-domain weak spring (f =.5 Hz) M M = 5 kg F 0. 0.5 0. 0.05 ecitation force (offset 0. &scaling e-4) Ft ( ) = 400sin( π 7t) 0-0.05-0. response H( 7Hz) 0045. / 400= e 4 m/ N H( 7Hz) 80 0-0.5 0 0.5.5.5 time in sec 8
finding a solution of the equation of motion: choose input: then: solution: F = M F = F sin( ωt) = sin( ω t+ ϕ ) =?; ϕ =? t ( ) = F sin( ωt) + ct+ c Mω H = = F Mω log( H) = log logω M F M H = = F Mω H = ϕ = 80 9
measurement mechanics stage amplitude in db -0-40 -60-80 0 0 3 frequency in Hz phase in deg 00 0-00 0 0 3 0
Derivation of transfer function - make a model of the dynamics: differential equations - substitute s=d./dt - rearrange the equations and get the transfer function e.g. H(s) - for sinusoids make a Bode plot using s=jω
Transfer function: H ( s) ( s) = = F( s) Ms + ds + k F Ms + ds + k
consider sinusoidal signals ('Euler notation'): ( t) = ˆ(cosωt + j sin ωt) = e ˆ jωt & ( t) = ωˆ( sin ωt + j cosωt) = jωe ˆ jωt apparently: s = jω for sinusoidal signals Frequency Response Function: s jω H( jω) = Mω + jdω + k 3
e k p + k v d dt F F = k p e + k v e& F() s = ( k + k s)() e s transfer function: F Cs () = () e s = ( k k s p + v ) frequency response: p C = k + jk ω p v v 4
Amplitude: C = + k p k v ω ω 0 C k p C 0 ω C k v ω C 90 ω ( C ) = log k logω log + v break point: log k = log k + log ω p ω= k k p v v 5
Bode plot of the PD-controller: amplitude in db 00 kp = 500 N/m; kv = 0 Ns/m 80 +0 60 0 0 0 0 frequency in Hz phase in deg 00 + 50 0 0 0 0 0 6
Block manipulation e C(s) H(s) s () es () = CsHs () () C + H - k M p k p k v frequency CH - - frequency 7
s + e - C(s) H(s) H c = s = CH + CH 8
Four important transfer functions. open loop: Ho( s) = C() s H() s s + F d e + C(s) H(s) - F s +. closed loop: H c ( s) = s ( s) = C( s) H ( s) + C( s) H ( s) 3. sensitivity: 4. process sensitivity: e S() s = () s = + CsHs () () H ps s s F s Hs () () = () = + CsHs () () d 9
s + F d e + C(s) H(s) - F s + Derivation of closed-loop transfer functions: start with the output variable of interest go back in the loop, against the signal flow write down the relations, using intermediate variables stop when arrived at the relevant input variable eliminate the intermediate variables 30
amplitude in db 0 closed loop 0-0 -40 phase in deg 00 0 0 3 0-00 amplitude in db 0 0-0 -40 phase in deg 00 0 0 3 open loop 0 0 3 frequency in Hz Eperimental results: stage servo 0-00 0 0 3 3
bandwidth: 0 db crossing open loop (cross-over frequency) 3
The Nyquist curve litude in db 50 Bode plot.5 Nyquist plot 0-50 -00 0 0 hase in deg 0 frequency 0.5 0-0.5 amplitude (appr 0.7) phase (appr. -75 deg) -00 - -00 -.5 0 0 -.5 - -0.5 0 0.5.5 frequency 33
Stability: The open-loop FRF CH(jω) should have the (-,0) point at left side.5 Nyquist plot Im 0.5 0 Re -0.5 ω ω - -.5 -.5 - -0.5 0 0.5.5 34
4. Filters Mimo Integral action Differential action Low-pass High-pass Band-pass Notch ( sper ) filter PeeDee PeeEye 35
Integral action : X(t) τ i s Y(t) - τ I integral time constant τ I =/k i 0 ω=πf -90 36
Differential action + u u H = ks = ; s = jω; = ε ε kω +90 ω ε ks u 0 u ε tamme differentiator : = τ d ks s + +90 τ d = ω d = πf d ω 0 37
38 s s s s u H d d γ τ τ τ τ ε + + = + + = = lead filter γ> 0 +90 ω τ ω= τ ω = 5 τ τ γ τ τ ω ω ω = = c
39 +90 0-90 ω i i τ ω = d d τ ω = - + P+I+D + + + = = γ τ τ τ ε s s s k u d d i
ε(t) ω - s nd order filter ω s u(t) H = u ε = s ω + k β s ω + β β Top:. ω o = ω β 0 - ε -90-80 ω=ω 40
4 + + + + = = ω β ω ω β ω ε s s s s u H General nd order filters + 0 +80 ω ω General: ω ω
ω ω + ω ω - ω 0-80 ω ω 4
Notch -filter :ω = ω ampl. β β fase 0-80 43
W.B.E. 44
Loop shaping procedure. stabilize the plant: add lead/lag with zero = bandwidth/3 and pole = bandwidth*3, adjust gain to get stability; or add a pure PD with break point at the bandwidth. add low-pass filter: choose poles = bandwidth*6 3. add notch if necessary, or apply any other kind of first or second order filter and shape the loop 4. add integral action: choose zero = bandwidth/5 5. increase bandwidth: increase gain and zero/poles of integral action, lead/lag and other filters during steps -5: check all relevant transfer functions, and relate to disturbance spectrum 45
Implementation issues. sampling = delay: linear phase lag for eample: sampling at 4 khz gives phase lag due to Zero-Order-Hold of: 80º @ 4 khz 8º @ 400 Hz 9º @ 00 Hz. Delay due to calculations 3. Quantization (sensors, digital representation) 46
5 Feedforward design 47
Why feedforward? Consider the simple motion system m F Control problem: track setpoint 0 0 s Setpoint t [s] s.8 Is this possible with a PD-controller? 48
Analysis (IV) -3 0 error [m].5 0.5 0-0.5 K p = 65000 [N/m] K p = 6500 [N/m] m = 5 [kg] K v = 60 [Ns/m] - -.5-0 0. 0.4 0.6 0.8..4.6.8 t [s] 49
Feedforward based on inverse model ms s K p + K v s ms 50
Eample: m=5 [kg], b= [Ns/m], nd degree setpoint s [m] 0.5 0.5 vs [ms-] 0.5 0 4 as [ms-] 0 - -4 0 0. 0.4 0.6 0.8..4.6 t [s] 5
Eample: tracking error, no feedforward -3 0.5 viscous damping effect error [m] 0.5 0-0.5 - -.5-0 0. 0.4 0.6 0.8..4.6.8 t [s] 5
Eample: tracking error, with feedforward -3 0.5 K fv = 0.9, K fa = 0 error [m] 0.5 0-0.5 = 0.9, = 4.5 K fa K fv - -.5-0 0. 0.4 0.6 0.8..4.6.8 t [s] 53
feedforward structure sign( & s ) K fc & & s K fa & s K fv s C(s) H(s) 54
3rd degree setpoint trajectory.5 s [m] 0.5 0 vs [ms-] 0.5 0 4 as [ms-] 0 - -4 0 0. 0.4 0.6 0.8..4.6.8 t [s] 55
6. Servo-oriented design of mechanical systems 56
Eample of measurement: mechanical system (force to position) modelling understanding the dynamical behaviour 57
Three Types of Dynamic Effects - Actuator fleibility - Guidance fleibility - Limited mass and stiffness of frame 58
. Actuator fleibility F s Motor k Sensor d 59
. Guidance fleibility F s M, J k 60
3. Limited mass and stiffness of frame F s Motor Frame 6
M M Positioning the load M (while using for feedback): Rule of thumb: Optimal bandwidth with 0 db crossing of open loop between the antiresonance and resonance frequency of the mechanical system. 6
Concluding Remarks bit of control into mechanical design bit of mechanics into control design same language ( mechatronics ) 63