Root Locus Design Example #4 A. Introduction The plant model represents a linearization of the heading dynamics of a 25, ton tanker ship under empty load conditions. The reference input signal R(s) is the desired heading angle for the ship ψ ref (s), and the output signal Y (s) is the actual heading (yaw) angle ψ(s). In this example, angles will be expressed in degrees. The input to the plant, U(s), is the commanded rudder angle δ r com (s) that is used to control the heading of the ship. The transfer function for the system is G p (s) = Y (s) U(s) = 3.2587 4 (s +2.55 2 ) s (s +.3333) (s +5.3288 2 )(s +6.8624 3 ) = ψ(s) δ r com (s) The gain 3.2587 4, the poles at s = 5.3288 2 and s = 6.8624 3,andthezeroat s = 2.55 2 describe the dynamics of the system between the actual rudder angle and the rate of change in heading angle. The pole at s =provides the integration from the rate of change in heading angle to the heading angle itself. The pole at s =.3333 models the hydraulic actuator dynamics between the commanded rudder angle δ r com (s) and the actual rudder angle δ r (s). The negative sign associated with the gain indicates that a negative rudder angle produces a positive rate of change in heading angle. This is from the usual convention of how a coordinate system is fixed to the ship, a step that is analogous to assigning directions of positive current flow within an electrical circuit. Because of this sign convention used in ship steering, the sign of the compensator gain must also be negative, meaning that a positive heading angle error produces a negative rudder angle. The performance specifications that are imposed on the system are: Percent overshoot to a step input must satisfy PO 2%; Settling time for a step input must satisfy T s 2 seconds; Steady-state error in the closed-loop ramp response must not exceed 2 degrees. B. Evaluating G p (s) Relative to the Specifications The first step in determining what type of compensation is needed is to evaluate the plant model relative to the specifications. Since the specifications are given in terms of percent overshoot and settling time, root locus will be the design method. Therefore, the desired location of the dominant closed-loop pole s = s must be determined. Since the plant is not second-order, it is not reasonable to assume that the second-order system equations will be valid, so a conservative approach will be used. The values used for percent overshoot and settling time will be PO design = PO spec =4%, T s design =.75 T s spec =5sec (2) 5 Using these design values and the equations for second-order systems, the dominant closed-loop pole is calculated to be ζ = s = h POdesign Ln i r ³ h =.756 (3) i 2 π 2 POdesign + Ln µ 4 +jtan cos (ζ) = 2.6667 2 + j2.626 2 (4) T s design ()
2 Figure shows the root locus and step response plots for the uncompensated system K c G p (s) with K c =. This gain of is included with the plant at this point so that the positive root locus methods (K >) can be used. The actual gains of the final compensator and the plant will be negative. K c is only used in the evaluation of the uncompensated system. Unity feedback is assumed, so H(s) =. The upper right plot in Fig. clearly shows that the root locus does not go through the point s, and the step response plot clearly shows that the overshoot and settling time specifications are not satisfied. Therefore, some form of compensation is needed. The angle of the plant and compensator at s must be computed. K c G p (s ) = tan 2.626 2 2.626 tan 2 (5) 2.6667 2 2.55 2 2.6667 2 tan 2.626 2 tan 2.626 2 2.6667 2.3333 2.6667 2 5.3288 2 tan 2.626 2 2.6667 2 6.8624 3 = 92.545 35.7 4.855 44.353 27.27 = 29.62 4.38 G c (s )=8 K c G p (s )=39.62 (6) Since the required phase shift of the compensator at s is positive, the compensator will be phase lead. C. Compensator Designs ) Overview: The compensator design technique discussed in the text which calculates both the pole and zero angles at s will be used. This requires computing the phase angle of the point s, which is 2.626 s =tan 2 =35.7 (7) 2.6667 2 The lead compensator will be designed using this method. Once that design is completed which hopefully will result in the transient performance specifications being satisfied the steady-state error of the plant/compensator combination will be checked. If the error is too large, then a special lag compensator will be designed to satisfy that specification. If the transient performance specifications are not satisfied by the lead compensator, several options exist to try and correct the problem. Some of these are shown below. Choose another value for s, using either more conservative or less conservative choices for percent overshoot and settling time. Choose different locations for the compensator zero and pole. Reconfigure the original design into the Proportional+Derivative with Derivative on Output Only (PD-DOO) version. Only the last option, the PD-DOO configuration, will be used in this example. In general, any or all of these options can be used together to try and obtain a compensator design that satisfies all specifications. K. Ogata, Modern Control Engineering, 4th Edition, Prentice Hall, Upper Saddle River, NJ, 22.
3.6 Uncompensated Root Locus. Zoomed View.4.2.5 s Imag Axis Imag Axis.2.5.4.6.5.5 Real Axis...5.5. Real Axis.5 Uncompensated Step Response.5 PO = 43.8%, T s = 45.8 sec 5 5 2 25 3 35 4 45 5 Fig.. Root locus and step response for the uncompensated system.
4 2) Design of the Lead Compensator: Using the method in the text, the phase angles from the compensator zero and pole to the point s are computed first, then the distances from the projection of s on the real axis to the zero and pole locations are computed. The angles are (s + z cd )= s + G c (s ) 2 = 35.7 + 39.62 2 =87.66 (8) (s + p cd )= s G c (s ) 35.7 39.62 = =48.4 (9) 2 2 Note that (s + z cd ) (s + p cd )=39.62 = G c (s ) as required. The distances from s to the zero and pole are d zcd = d pcd = Im [s ] tan ( (s + z cd )) = Im [s ] tan ( (s + p cd )) = 2.626 2 tan (87.66 π/8) =.635 3 () 2.626 2 tan (48.4 π/8) =2.345 2 () Since both of these distances are positive, both the pole and zero of the lead compensator are to the left of s. The zero is located at s = 2.773 2, and the pole is located at s = 5.7 2. At this point in the design, the lead compensator is G c Lead (s) = K c (s +2.773 2 ) (2) (s +5.7 2 ) Now that the lead compensator s pole and zero have been placed to satisfy the root locus phase angle criterion, the gain must be computed to satisfy the magnitude criterion at s. The gain is K c = s +5.7 2 s s +.3333 s +5.3288 2 s +6.8624 3 s +2.773 2 3.2587 4 s +2.55 2 K c = 2.299 K c = 2.299 (3) Note that the sign on the gain is negative. The forward path transfer function is now G c Lead (s)g p (s) = 7.25 4 (s +2.55 2 )(s +2.773 2 ) s (s +.3333) (s +5.3288 2 )(s +5.7 2 )(s +6.8624 3 ) (4) D. Evaluation of the Design The design will be evaluated by examining the step response of the lead-compensated system. The total forward transfer function is given in (4). The root locus and step response plots are shown in Fig. 2. The root locus plot shows that the point s is on the root locus as it should be. The closed-loop poles are located at s = 2.6667 2 ± j2.626 2,s= 2.477 2 ± j2.592 2, and s =.3477. The step response plot shows that neither of the transient performance specifications has been satisfied. The overshoot of approximately 3% is not acceptable, and the settling time of 23 seconds is too long. The various options that may be used to try and correct this problem were mentioned earlier. The only option that will be described in this example is the Proportional+Derivative (PD) form of the compensator with the Derivative on Output Only (PD-DOO) configuration. The steady-state error of this system for a ramp input is e ss = K v = lim s [s G c Lead (s)g p (s)] =.98 (5)
5 so a special lag compensator would be needed in order to satisfy that specification. However, before that is done, the transient response specifications need to be satisfied. There is no point in designing the special lag compensator until the transient performance is satisfactory. E. PD Compensator with Derivative on Output Only The gains and time constant of the Proportional+Derivative (PD) controller are τ = p cd, K p = K c zcd p cd, K d =( K c K p ) τ (6) and the values are τ =9.97 sec, K p =.2239, and K d =9.693, so if the PD compensator was to be placed in series with the plant it would be G c PD (s) =.2239 + 9.693s (7) 9.97s + where the negative sign in G c PD (s) is required since the controller gain is negative. The PD-DOO configuration is [ G p (s)] G DOO (s) = K p K d s +[ G p (s)] τs+ = (8) 3.988 4 (s +5.7 2 )(s +2.55 2 ) s (s +3.37 )(s +.482 2 )(s +4.584 2 ± j2.76 2 ) where the negative sign of the compensator is now included with the plant transfer function. The step response of the system with the PD-DOO configuration is shown in Fig. 3. Both the percent overshoot and the settling time satisfy the transient response specifications. The steady-state error specification does have to be checked to see if a special lag compensator is needed. F. Design of the Special Lag Compensator The steady-state error for a ramp input with the PD-DOO configuration / lim s [sg DOO (s)] = 28.7 (increased from.98 by the change in configurations), and the specified value is 2. Therefore, the error must be reduced by a factor of α g = e ss actual = 28.7 =4.35 = z cg (9) e ss spec 2 p cg This value for α g reduces the steady-state error to the correct value by separating the special lag s pole and zero by the same factor. Using the rule of thumb discussed in class, the compensator zero is placed to the right of s by a factor of, and as always p cg = z cg /α g, so the special lag compensator is G c Spec Lag (s) = (s +2.6667 4 ) (2) (s +.9 5 ) The step response is shown in Fig. 4. The overshoot and settling time are both very close to the values in Fig. 3, and they still satisfy the specifications. Therefore, the special lag compensator did not disturb the transient response very much. The ramp response of the final version of the compensated system is shown in Fig. 5. The graph illustrates the very long time that it might take for the ramp response to settle to essentially a constant slope. At t = seconds, the error is still larger than 3. Even though it taking a long time to reach steady-state with the ramp response, the steady-state error does have the correct value after the special lag compensator is included.
6.6 Compensated Root Locus. Zoomed View.4.2.5 s Imag Axis Imag Axis.2.5.4.6.5.5 Real Axis...5.5. Real Axis.5 Compensated Step Response.5 PO = 29.6%, T s = 23. sec 5 5 2 25 3 35 4 45 5 Fig. 2. Root locus and step response for the lead-compensated system in the normal configuration.
7.4 Compensated Step Response with PD DOO Configuration.2.8.6.4.2 PO = 7.6%, T s = 86.5 sec 5 5 2 25 3 35 4 45 5 Fig. 3. Step response for the compensated system in the PD-DOO configuration.
8.4 Compensated Step Response with PD DOO Configuration and Special Lag Compensator.2.8.6.4.2 PO = 8.7%, T s = 92.3 sec 5 5 2 25 3 35 4 45 5 Fig. 4. Step response for the compensated system in PD-DOO configuration with special lag compensator included.
9 Compensated Ramp Response 9 8 7 6 5 4 3 2 2 3 4 5 6 7 8 9 Zoomed View 9999 9998 9997 9996 9995 9994 9993 9992 999 999 999 999 9992 9993 9994 9995 9996 9997 9998 9999 Fig. 5. Ramp response for the compensated system in PD-DOO configuration with special lag compensator.