HW5-Solutions

.pdf

School

Purdue University *

*We aren’t endorsed by this school

Course

375

Subject

Mechanical Engineering

Date

Apr 3, 2024

Type

pdf

Pages

13

Uploaded by ColonelHeat13546

Report
ME375 Homework 5 - Spring 2024 Due at 5:00pm on Wednesday, March 20 R ( s ) + Y ( s ) P ( s ) C ( s ) - Figure 1: Use this figure for Problem 1 1. (Lead compensator design) Given the plant P ( s ) = 1 s ( s + 1) in the block diagram of Figure 1, your objective is to design a lead compensator, C ( s ), such that the closed loop system will have about 15% maximum overshoot and 2% settling time of about 1.33 s. a. Estimate the closed loop bandwidth and phase margin corresponding to about 15% overshoot and 2% settling time of 1.33 seconds. (For simplicity, assume the closed loop bandwidth, closed loop natural frequency, and gain crossover frequency are equal.) Percent overshoot is related to damping ratio % overshoot = 100% e πζ 1 ζ 2 (1) 15% = (2) ζ = 0 . 52 (3) Damping ratio is approximately related to phase margin (in degrees) ζ PM 100 (4) PM 100 ζ (5) 100 · 0 . 52 (6) = 52 (7) 2% settling time t 2% = 4 ζω n (8) ω n = 4 ζt 2% (9) = 4 0 . 52 · 1 . 33 (10) = 5 . 78 rad/s (11) Approximate bandwidth ω BW ω n (12) 5 . 78 rad/s (13) b. Use frequency domain controller techniques to construct an initial design for a lead compensator of the form C ( s ) = K s z s p to achieve the objective. (Note, your initial design need not perfectly meet the objective.) Approach 1 (gain crossover frequency is in the ballpark) 1 Copyright 2024, Todd Lillian
ME375 Homework 5 - Spring 2024 Due at 5:00pm on Wednesday, March 20 Let’s approximate the desired gain crossover frequency with bandwidth. ω gc ω BW (14) 5 . 78 rad/s (15) Find K to satisfy gain crossover frequency K = 1 | P ( gc ) | (16) = | gc ( gc + 1) | (17) = | j 5 . 78( j 5 . 78 + 1) | (18) = | j 5 . 78 | | ( j 5 . 78 + 1) | (19) = 5 . 78 · p 5 . 78 2 + 1 2 (20) = 33 . 9 (21) How much phase angle do we need to add at ω gc to meet our margin requirement (including somewhere between 5-12 extra; lets use 7 . 8 ) additional phase = PM ̸ P ( gc ) + 7 . 8 (22) = PM (180 + ̸ P ( gc )) + 7 . 8 (23) = 52 (180 + ̸ 1 gc ( gc + 1) ) + 7 . 8 (24) = 52 (180 + ̸ 1 j 5 . 78( j 5 . 78 + 1) ) + 7 . 8 (25) = 42 . 2 + 7 . 8 (26) = 50 (27) Calculate the corresponding ratio α = z p α = 1 sin(additional phase) 1 + sin(additional phase) (28) = 1 sin(50 ) 1 + sin(50 ) (29) 0 . 133 (30) Lead compensator C ( s ) = K 1 α s + ω gc α s + ω gc α (31) = 33 . 9 1 0 . 133 s + 5 . 78 0 . 133 s + 5 . 78 1 0 . 133 (32) = 255 s + 2 . 11 s + 15 . 8 (33) Although this approach will get us sufficiently close to an acceptable controller that we can iterate, it yields a larger than required gain crossover frequency. In many applications this may be a good thing because the system will respond faster. Approach 2 (exact gain crossover frequency) Let’s approximate the desired gain crossover frequency with bandwidth. ω gc ω BW (34) 5 . 78 rad/s (35) 2 Copyright 2024, Todd Lillian
ME375 Homework 5 - Spring 2024 Due at 5:00pm on Wednesday, March 20 How much phase angle do we need to add at ω gc to meet our margin requirement additional phase = PM ̸ P ( gc ) (36) = PM (180 + ̸ P ( gc )) (37) = 52 (180 + ̸ 1 gc ( gc + 1) ) (38) = 52 (180 + ̸ 1 j 5 . 78( j 5 . 78 + 1) ) (39) = 42 . 2 (40) Calculate the corresponding ratio α = z p α = 1 sin(additional phase) 1 + sin(additional phase) (41) = 1 sin(42 . 2 ) 1 + sin(42 . 2 ) (42) 0 . 196 (43) How much gain, K , do we need to make ω gc the crossover frequency K = 1 1 α gc + ω gc α gc + ωgc α P ( gc ) (44) = 1 1 α j + α αj +1 P ( gc ) (45) = 1 1 α | P ( gc ) | (46) = α | gc ( gc + 1) | (47) = . 196 | j 5 . 78( j 5 . 78 + 1) | (48) = . 196 | j 5 . 78 | | ( j 5 . 78 + 1) | (49) = . 196 · 5 . 78 · p 5 . 78 2 + 1 2 (50) = 15 . 0 (51) Lead compensator C ( s ) = K 1 α s + ω gc α s + ω gc α (52) = 15 . 0 1 0 . 196 s + 5 . 78 0 . 196 s + 5 . 78 1 0 . 196 (53) = 76 . 6 s + 2 . 56 s + 13 . 1 (54) c. Use software to compute the phase margin, gain margin, and gain crossover frequency for your design. see Matlab script below which uses the second approach from above d. Use software to plot the closed loop step response for your system. see Matlab script below which uses the second approach from above e. Use software to determine the maximum percent overshoot and settling time for your system. (Hint, you may want to use the stepinfo(sys) command in Matlab.) see Matlab script below which uses the second approach from above 3 Copyright 2024, Todd Lillian
clear all clc close all P=tf(1,[1,1,0]); C=zpk(-2.56,-13.1,76.6); margin(P*C) figure() sys=feedback(P*C,1); step(sys) stepinfo(sys) ans = struct with fields: RiseTime: 0.2056 TransientTime: 1.0396 SettlingTime: 1.0396 SettlingMin: 0.9154 SettlingMax: 1.2093 Overshoot: 20.9291 Undershoot: 0 Peak: 1.2093 PeakTime: 0.5006 1
Published with MATLAB® R2021b 2
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help