M.I.T. SERVO\251 2001 J.M. Redwood.Revised 4 December 2004.<Text-field layout="Heading 2" style="Heading 2">Synopsis</Text-field> This worksheet derives the loop transfer function of a tacho feedback compensated servo. It then plots the servo's Nyquist diagram, Bode diagrams, and closed loop responses to sinusoidal and unit step inputs with servo constants chosen to minimize noise in the input.<Text-field layout="Heading 2" style="Heading 2">Licence</Text-field> This worksheet and the package of procedures entitled "Bode2", may be used, adapted, or copied, provided that copyright is acknowledged and the whole of this Licence section is included in any copies, or adaptations, of the worksheet, or of the package. This worksheet and the "Bode2" package may be used for educational purposes free of charge. Commercial users are not required to pay a licence fee for this worksheet, or for the "Bode2" package. Instead, they are requested to make a donation to Cancer Research UK, National Office, 61 Lincoln's Inn Fields, LONDON, WC2A 3PX, UK. The donations line telephone number is +44 (0)20 7009 8820; the Fax number is +44 (0)20 7269 3100; and the E-mail address is "customerservices@cancer.org.uk". The size of donation is left to the user's discretion. This worksheet and associated "Bode2" package of procedures are offered without warranty of any kind whatsoever. <Text-field layout="Heading 2" style="Heading 2">Resources</Text-field> This worksheet uses significant computer resources. If necessary, the user should separate the worksheet into manageable portions, comprising the section on the Loop Transfer Function plus one, or two, other sections, e.g. Loop Transfer Function and Bode Diagrams. Copying the block diagram to a different file, and then deleting it from the worksheet may also help to free resources.<Text-field layout="Heading 2" style="Heading 2">Historical</Text-field> C.H. Dowker used R.S. Phillips' "RMS-Error Criterion" to carry out this classic design while they were at the M.I.T. during World War II. The design minimised the effect of noise while obtaining the best response to the signal from an aircraft-tracking radar. The details are described in Chapter 8 of the reference. Based on this design, and on another also described in Chapter 8, C.H Dowker suggested a simple design rule that produced the optimum performance, without needing to repeat the detailed calculations for minimising the RMS error.<Text-field layout="Heading 2" style="Heading 2">Introduction</Text-field> The servo described in this worksheet formed part of an automatic aircraft-tracking radar in World War II that was used by C.H. Dowker in his work on minimizing the effect of noise. The servo's input was a signal proportional to the angular error between the centre of the radar beam and the aircraft. Naturally, the signal contained a great deal of noise generated by changes in the radar signal caused by random changes in radar propagation and reflection from the aircraft's surfaces. By appropriate choice of the servo's constants, Dowker found that it was possible to minimize the effect of noise on the servo while obtaining good performance from it. Chapter 6 of the reference gives the method used by Dowker for obtaining the best choice of servo constants. This worksheet obtains the loop transfer function of the proportional control servo with tacho feedback compensation, as shown in the block diagram. The servo's constants, derived by C.H. Dowker at equations 87, 89 and 90 in Chapter 8 of the reference, are then inserted in the loop transfer function. These are the constants needed to minimize the effect of noise while producing the best response from the servo. The loop transfer function is then plotted in Nyquist diagrams and in Bode attenuation and phase margin diagrams, which yield parameters describing the servo's closed loop stability and performance. Plots showing the servo's frequency response, peak magnification, damped natural frequency and bandwidth are produced, together with transient responses to unit step functions of input position and of velocity. Important operational parameters are obtained from these plots. The nomenclature used in the reference has been preserved in this worksheet as far as possible, so that the reader can move easily from one to the other. The block diagram shown at Figure 8.4 of the reference, in slightly different form, has been inserted into this worksheet as an OLE2 object (or as an image in the Standard interface).<Text-field layout="Heading 2" style="Heading 2">Block diagram</Text-field>TUZOV3RLVWI8b2I8Uj1NRExDZE5OWmpqOnRzOkg6PHVaYXRRSj08T2BMb1xcanl5eXlRdGZ5c2hjO2VVTENUSmNEWG9YdXNWTVNFcUNEPzt2a1doQ3FXXT9kUk9nPEVxQEFYRDxoUV5uQmdjU3diZElvQVBrPXZpQV5uXl5sTWdaXFxeXUpPWk5OYWVGTmlicEFiPklCO1VGTkNGPklHTz9EXFxBVGJXVHQ9Z2RbUnVVRjxPVGc9RkFfRDtLeHBTY15pdXJ1dEZ3WVBZd2t5c3JdeU1FdFl3SENvSHVveF5ZeG1ZZXFzV1hZWVhfUmpDVEtjUz49YjtbQ1o9Qj07bTtcXEtab1pqPko6alZJdFJUb05AeWdgVEJgUGlpdG1ZUz1ITDxRdEdddEdJb1VsdGVgdVp0dENgbHBveEVXaGhpeEtHeFNodHlebEp3YHN3aD5gbVh4XVdIYmJnd2thXkZ3cnJvaXhQaW5Bc21HdmhnY25Oak9edGRWW0JGY0xPdkdnW0R3bUVXXFxBWF50aGxuVndvcHNhWG1IaF5UQFxcV0l5bHlueVFbaU5eTmdvXFxBb21oXmRfblBPdmlwaHdveE95Y0BIbHNPXVd4dlNGb2VHYFxceHRpPl1JV21BSWxzeHh1Z15LT21ZdnFlYGNWbl10X3hhQF1LUHFsRlxcUmdhTkBrU3BtQE9nbmdnVz9kZ3BrdE50R19tQWF5Pl5gd2hucGdjTlF0SHlva0FmYFFua3ZdV09dWWFpcz92YnFzPVhwYW9vPm9lQG5dSmI7WFFSVFJzaG9AaHFwPUtXQVZhaFNTUHlSZVRTcWs8dW9ScXFsXFxOW3BSaWRTZFBKXWR4eUBWVGhyWk11P0FWbmBZPmh3XFxdWE1hWGtUVVZoeVVEbDt3Pj1DQEFYeXN5WD9ic3d2ZW1maE1ZbUN0Qk1EQXVjdktZaztmVV9Sb0NXa0FYa2tyYV1yZEVGTmFjOm9XTG9XRGFUU09pU3N0VkNZaU90S0tVUmlIY1tTaXl0RU9pQU1kaE1TYFV3QGlDXXVmTD9IS0lTQ1lWRFdjW3NJcXFiSHFSXmtDTz9TRmF3RHVZXm95Q2Nldz1GXj90V1dWWV1EW3lXQHFJdWViQ0d4VVlybFNDbk9zQEFWVVVZdHdIWWVHVUtYPmFWUk9TSVFGPkdJZkFWYEFjcGNTP3FoXFxVeFdxSXdVR3BdY1dBY3V3bFFdVVdBc3BNeFxcZXFvZHJ4dW1wPVdGeFVqeVJ0UUpwSXVNQHdUXFxSUU14amlxeFRtW0hQeFltRD1NTnRtd0RraEhSUmRMQXR1Ynl4U3VyZkl3ZFR5ZklTanV0YVF4b014bWBvcUVLcVh2eHF3Q1BtcVxcdFR4bkRISmxEamVFcmZsbmRQdVtASkNESm5Ncj5ZU09US111ajxcXFJpPW12RU1fQHlKPW5IXXJreHNNeFVRZHVTaVRMaXdebXVNeUtTRFdlVG9uPFlWeWtoUXdAQHNmTFdwZWpcXHhueFVLc0RMZXlgVHFubnZcXE13cEZAc2t4bltwaUtZZkJweFVoYEVnbldRYU9Ha2phckxmZWhGeENveEdeXVlxeEZfaFlmYUlvaV1mcU5xbj5nbFFfYWtgZ09OZll4XWJPZHhHZUNeZXZfbldGYGNYbkZWbkxWZ2BOZlRAcl5uY2N4XktOZz95X1ledXhxamJHb1VRYT5pX0dBZHZYXlBOc2FHcWtWeV15ZkpAZTo/dHBhXFw8V2pTP2A7ZmVdb25DV3FScGpWYF5HSW1XXnFFYGxXV2pDSF1FX19NaXRBcWt3P2NWR2NkWW10QV5aUWlreW9EX3FTZ3Jdb1xcPXZtT3hqXVlmVmBcXEVpbTxBW0FuZ3g/cFlObV14YVFJYW1YYXBXYWo/WmxJXzplQ3VUa0NXZ0ljaD9lbmt3X3FnUj9WVkt4Y2lkXkVHRk9UW3NIP1d3bmtiT1dYa19VXFw9RF1zVjxBc2RjY2NxYk1nZnk/WUxJREZVVkY/aHhzV2o9VkRJQnVHdD1xdW5VZD1TSFZIbDxBbW9dbUtheWJhWVN0VWZMVmw8S0JFSkFFSlI8cjpsOkE6R1Zvb1NjRWJIS0NfP1RHR2JBU0hqX0hBZVlYY2dBb2VdP1hXYVRjW1RjT2Rtd1Q/W2lgQ0NTYVRmd2lvUXRRX2Rsb3VidUhtO1ZIbUN4R2k+ZXJ0Z3RxV3VTO0hpR2ZnZ0REVURyaXlAbUJGZUZKZUNmQXY8XXNJYXNkZXJpYUJMT0l2a2Q7UWNJdVl2X3dxO3JwSXJzeVZRd1c/YWRPW3VQO1Y6U2REaWZVX3hPbUJGc3U7TWduY1REVXZaX0N5SVJPa1JEO0hxVXRWT1hmdVlyRVRyc2ZgcXJCc3VGYWVnb0RqVVRgRUh0R3dcXGNSU2VzSUdodllyTm1ZRnFSbG9Dcz1ITT9ja2FZVXdFXUVTbHdHeWtGcldHYFNGQD1XX2tUPklZZj95QllGWWVkYU94VUV4SG9Vc09VPmF0dllIYz1mckFIcWlCWD9FcW9FZFtyS11Dd1FkRFN0UUdmeUdEO1V4Z3NlVE1VRnFZVFNmYU9yblFpc2ljYG93Vj9JPz9HVV9WXFxNWVdxd3FpRW1ndVphZW1NU19Zcj49WUh1VkFRc0RnaT9Pc3F5RkxbQ21ZV2dJUm9TRHNVVU8/U3ZTdFNbWHdHZj5NVFlpeFNzR1xcSVZgPUheQXVDQ0ZSd2dST0JzcUhVZURHTVRPS3hYV0c9R1RcXGtWQGdZb09kQEdkeVVCUlNFXXNyWld1ST1mb0t5bFtTRld4PmNzbndndGVzQW9HWW11Yld0UUdEU0lFaWdjU11oUHNiU0tDd0N3dXVzQ0dYWG9kcHFFdEVTW0NXc1VXTjtJc0N3Rm1lW3F2bldoSW1UQz9DZ1V0Wj1zV0FXeTtzQmN2PUFVXFxdVEhrQlBtYmthdFJTV3JnZnZRaWtvaEhxU2B5cm5FRG1VYmBzaER3R1ZNaWxHRlxcZ0hvc0VgPUJ5Q2l3P2JKUVJYaUJAXUN0QWdmPUU7cXhDcUZEaWNDbXdCaWhWaWNOO0hhP1NVb2JOWUJqc2hrbVhhW1JOXVNWWWJjS0dgQ2R5dUlTY1hKa1JkcUc7YUJxW1VfXVdERUhSP3ZAV2dGZ0dyR3ZkU1Zlc0RdT3ZbPUh5SEpGYjs8cjpidWtDPVdDO0FlP0NWREdpY29XaztlPnl5Tzt2X19nb09VUEdERGN0TEllR3VHXklUSE1YcztVZmtydD95aktSOz1FdlVHR21CanlQPTw9Pjo6WlJYdTxWeFlueEBBcGZmXUVhcV5AOjtCOk1US1dES1dnSjtlWjE6<Text-field layout="Heading 3" style="Heading 3">Definitions</Text-field>q i is the angular input to the servo (also known as R, the reference input).q o is the angular output from the servo (also known as C, the controlled output). NiMlKGVwc2lsb25H is the error between the input and output = q i - q o (= R - C). <Text-field layout="Heading 2" style="Heading 2">Loop Transfer Function</Text-field><Text-field layout="Heading 3" style="Heading 3">Initialisation</Text-field> The "Bode2" package will be required for plotting the Bode attenuation and phase margin diagrams. It is also necessary to call the "plots" package and to change the imaginary unit to j, instead of I. It is also convenient to define some settings that will be used in the various plots.restart:with(plots): with(inttrans): with(Bode2): interface(imaginaryunit=j):<Text-field layout="Heading 3" style="Heading 3">Loop Transfer Function</Text-field> By inspection, the inner loop transfer function (from amplifier output to the servo's output) isNiMvLSUjWWlHNiMlInNHKiYtJSNZZkdGJiIiIiwmRitGKyomRilGKy0lJFlmYkdGJkYrRishIiI=where Yf(s) is the forward loop transfer function = NiMqJiUjS21HIiIiKiYlInNHRiUsJiomJiUiVEc2IyUibUdGJUYnRiVGJUYlRiVGJSEiIg==, and Yfb(s) is the feedback loop transfer function = NiMqLiUjS3RHIiIiJSJzR0YlJSNLZkdGJSYlIlRHNiMlImFHRiVGJkYlLCYqJkYoRiVGJkYlRiVGJUYlISIi (also referred to as H). By definition, the loop transfer function for the servo is, NiMvLSUiWUc2IyUic0cqKCUmdGhldGFHIiIiLSUib0dGJkYqLSUoZXBzaWxvbkdGJiEiIg== = NiMqJiUjR2FHIiIiLSUjWWlHNiMlInNHRiU= It is convenient to abbreviate NiMtJSJZRzYjJSJzRw== to Y, and similarly for the other loop transfer functions.Yi := Yf/(1 + Yf*Yfb);Yf := Km/(s*(Tm*s + 1));Yfb := (Kt*s)*(Kf*Ta*s/(Ta*s + 1));Y := Ga*Yi; From the block diagram, the gain of the tacho loop is seen to be NiMvKiYlI0t0RyIiIiUjS2ZHRiYlIkFH. Substituting this in Y, and collecting terms in s, Y := collect(algsubs(Kt*Kf=A,Y),s);This is the final symbolic form of the loop transfer function, before substituting values for the gains and time constants.<Text-field layout="Heading 3" style="Heading 3">Steady State Error</Text-field> At this point, it is convenient to determine the steady state error in response to a constant velocity input, NiMlJk9tZWdhRw==. From the block diagram,epsilon(s) := 'theta[i](s)/(1+Y)'; The Laplace transform of NiMvLSYlJnRoZXRhRzYjJSJpRzYjJSJ0RyomJSZPbWVnYUciIiJGKkYt is NiMvLSYlJnRoZXRhRzYjJSJpRzYjJSJzRyomJSZPbWVnYUciIiIqJEYqIiIjISIi. Hence,epsilon(s) := eval(subs(theta[i](s)= Omega/s^2,epsilon(s))); The steady state error is obtained by using the final value theorem of the Laplace transformation, ss error =NiMtJSZMaW1pdEc2JC0lKGVwc2lsb25HNiMlInRHL0YpJSlpbmZpbml0eUc= =NiMtJSZsaW1pdEc2JComJSJzRyIiIi0lKGVwc2lsb25HNiNGJ0YoL0YnIiIh. `ss error` := limit(s*epsilon(s),s=0);<Text-field layout="Heading 3" style="Heading 3">System Constants & Loop Transfer Function</Text-field> The system constants were chosen as follows (see equations 87, 89 and 90 in Chapter 8 of the reference) A = Tacho loop gain = Kt \264 Kf = 14, where Kt = Tacho Gain and Kf = Filter Gain (Attenuation), Ta = Tacho filter time constant = 0.252 second, Tm = Motor (including load) time constant = Jm/fm = 0.316 second, where Jm = motor and load inertias and fm = viscous friction and back e.m.f. Ga = Amplifier Gain \264 Error Measuring Device Gain = 158 volts/radian, Km = Motor Gain = radian/sec/volt, Kv = Ga x Km = Velocity error coefficient for the servo = 158/second. (Note. In this case the motor gain Km is unity, so Kv is numerically equal to Ga, the amplifier gain.) Assigning these values,A := 14: Ga := 158: Km := 1: Ta := 0.252: Tm := 0.316: Kv := Ga*Km:Y := Y; The loop transfer function in the frequency domain, Y(f), is`Y(f)` := subs(s=j*omega,Y);<Text-field layout="Heading 2" style="Heading 2">Nyquist Diagram</Text-field><Text-field layout="Heading 3" style="Heading 3">Nyquist Diagram</Text-field> The loop transfer function Y(f) is plotted on the Nyquist diagram as follows,p1 := complexplot(`Y(f)`,omega=3..1000,color=BLUE):p2 := plot([[-1,0]],style=POINT,symbol=circle,symbolsize=14,color=BLACK):p3 := coordplot(polar,[0..1,Pi/2..3*Pi/2],grid=[5,13],color=[RED,RED],linestyle=[1,2]):display([p1,p2,p3],title="Nyquist Diagram");
The Nyquist point, NiM3JCwkIiIiISIiJSNqMEc=, is shown on the diagram surrounded by a small black circle, for ease of identification. The loop transfer function shows a phase margin of slightly less than 60 degrees at feedback cutoff (unity magnitude), which indicates a very stable - possibly overdamped - system.<Text-field layout="Heading 3" style="Heading 3">Closed Loop Frequency Magnification</Text-field>The peak magnification, M, of the closed loop frequency response is obtained by plotting a circle of radius NiMqJiUiTUciIiIsJiokRiQiIiNGJUYlISIiRik=, with its centre at NiMsJComJSJNRyIiIywmKiRGJUYmIiIiRikhIiJGKkYq on the real axis of the Nyquist diagram, so that it just touches the loop transfer function. The loop transfer function's phase margin of nearly 60 degrees suggests that the peak magnification of the closed loop frequency will be in the region of 1.2 - 1.3, so four M circles are plotted on the Nyquist diagram thusM := [1.2,1.25,1.3,1.35]:x := seq((-M[i]^2/(M[i]^2-1) + M[i]/(M[i]^2-1)*cos(theta)),i=1..4):y := seq((j*M[i]/(M[i]^2-1)*sin(theta)),i=1..4):p := seq(complexplot(x[i]+y[i],theta=0..2*Pi,color=BLACK),i=1..4):display(p1,p2,p, scaling=constrained, title="Nyquist Diagram with Magnification Circles\nM = 1.2, 1.25, 1.3, 1.35" ); As NiNmKjYjJSJNRzciNiQlKW9wZXJhdG9yRyUmYXJyb3dHNiIiIiJGKkYqRio=, the radius of the M circle increases and its centre approaches NiMsJCUpaW5maW5pdHlHISIi, with the circle becoming a vertical line through the Nyquist point. Noting that the loop transfer function runs between the second and third largest circles, it will touch a circle with M between 1.25 and 1.3. Choosing M = 1.27, and plotting again,M := 1.27:x := -M^2/(M^2-1) + M/(M^2-1)*cos(theta):y := j*M/(M^2-1)*sin(theta):p := complexplot(x+y,theta=0..2*Pi,color=BLACK):display( [p1,p2,p], scaling=constrained, title="Nyquist Diagram with Magnification Circle M = 1.27"); The peak magnification, NiMmJSJNRzYjJSVwZWFrRw==, of the closed loop frequency response is seen to be 1.27. This agrees with equation 109 in Section 8.8 of the reference.unassign('p1','p2','p3','p','x','y','M');<Text-field layout="Heading 2" style="Heading 2">Bode Diagrams</Text-field><Text-field layout="Heading 3" style="Heading 3">Attenuation & Phase Margin Diagrams</Text-field> The Bode diagrams, i.e. the loop gain of the transfer function, Y(f), together with the asymptotic straight-line approximations, and the phase margin plotted against normalised frequency, NiMqJiUmb21lZ2FHIiIiJSJUR0Yl, are produced by procedures in the "Bode2" package. These procedure accept the loop transfer function expressed as a function of the Laplace operator, s, or as a function of the complex frequency. The latter is used in this worksheet for ease of comparison with the transfer functions as shown in the reference. To obtain the straight-line asymptotic approximations, the denominator of the loop transfer function must be factored, thus `factored Y(f)` := numer(`Y(f)`)/factor(denom(`Y(f)`),complex); The procedure for the straight-line approximation assumes that the loop transfer function is expressed in canonical form with gain K in the numerator, NiMpKiYlImpHIiIiJSZvbWVnYUdGJiUiTkc= in the denominator, and only terms of the form NiM3IywmKiglImpHIiIiJSZvbWVnYUdGJyYlIlRHNiMlImFHRidGJ0YnRic= in the numerator and denominator. In this case, N = 1 and the only time constant in the numerator is 0.252. To reduce the denominator to the canonical form, it must be divided by the product 0.2453105557 \264 51.19129784, yielding two time constants NiMqJiIiIkYkLSUmRmxvYXRHNiQiK2RiNWBDISM1ISIi and NiMqJiIiIkYkLSUmRmxvYXRHNiQiKyV5SCI+XiEiKSEiIg==. The numerator must also be divided by this product, which results in K = 158. This is the value of Kv shown in equations 87, 89 and 90 of the reference. It is worth noting that for a Type 1 servo, provided its loop transfer function is expressed in canonical form, K is the velocity error coefficient. It is normally shown as Kv and has the units of 1/sec. The straight-line approximations (shown in red) can be superimposed on the loop gain (shown in blue) by using the "display" command.display(attensemilog(`Y(f)`,.01,100),
linesemilog(158,1,[.252],[1/.2453105557,1/51.19129784],.01,100)); The gain Kv = 158 is shown by the intercept of the initial -6dB/octave slope of the loop gain with the NiMvKiYlJm9tZWdhRyIiIiUiVEdGJkYm axis at 42dB = 158. Feedback cutoff (also known as gain crossover) occurs about 1/3 the way down a -6dB/octave slope, at about NiMvKiYlJm9tZWdhRyIiIiUiVEdGJiomJSNLdkdGJiUiQUchIiI= = 10. (Or, since the normalizing time constant, T, is 1 second, NiMvJSZvbWVnYUcqJiUjS3ZHIiIiKiYlIkFHRiclIlRHRichIiI= = 10/sec.) The phase margin is plotted thus,Y(f); The big increase in phase margin between NiMvJSZvbWVnYUcqJiIiIkYmJSNUYUchIiI= = 4/second and NiMvJSZvbWVnYUcqJiUiQUciIiIlI1RtRyEiIg== = 44/second is caused by the tacho feedback loop producing a long -6dB/octave attenuation between these frequencies.<Text-field layout="Heading 3" style="Heading 3">Improved Phase Margin Diagram</Text-field> The shapes of the Bode diagrams above with semilog scales, particularly that for phase margin, are rather distorted below about NiMqJiUmb21lZ2FHIiIiJSJUR0Yl = 0.1. (This is a known problem with Maple's semilog plots - see Robert Israel's Maple Advisor.) The diagrams can be improved by using the "attenlog" and "phaselog" procedures in the Bode2 package, which plot the loop gain and phase margin against NiMtJSZsb2cxMEc2IyomJSZvbWVnYUciIiIlIlRHRig=. For example,phaselog(`Y(f)`,.01,100);unassign('omega');<Text-field layout="Heading 2" style="Heading 2"><Font bold="true" executable="false" family="Times New Roman" foreground="[0,0,0]" italic="false" size="14" style="_cstyle275" underline="false">Closed Loop Response</Font></Text-field><Text-field layout="Heading 3" style="Heading 3">Closed Loop Transfer Function</Text-field>The closed loop transfer function, or system response, is given byNiMvKiglJnRoZXRhRyIiIi0lIm9HNiMlInNHRiYtKiZGJUYmJSJpR0YmRikhIiIqJi0lIkNHRilGJi0lIlJHRilGLg== = NiMqJi0lIllHNiMlInNHIiIiLCZGKEYoRiRGKCEiIg== Again abbreviating NiMtJSJZRzYjJSJzRw== to Y and similarly for NiMqJi0lIkNHNiMlInNHIiIiLSUiUkdGJiEiIg==, `C/R` := Y/(1 + Y):`C/R` := factor(`C/R`); The closed loop transfer function in the frequency domain, C/R(f), is`C/R(f)` := subs(s=j*omega,`C/R`);<Text-field layout="Heading 3" style="Heading 3">Closed Loop Frequency Response</Text-field> The closed loop frequency response is usually plotted as magnitude in dB against NiMtJSZsb2cxMEc2IyomJSZvbWVnYUciIiIlIlRHRig= thus,Freq_resp := 20*log10(simplify(abs(`C/R(f)`))): omega := 10^w:p := plot(Freq_resp,w=-1..1.5,-10..3,labels=["log10 (omega T )",""],title="Closed Loop Servo Frequency Response \n Magnitude dB"): q := pointplot([[-1,-3],[1.5,-3]],connect=true,colour=black,linestyle=2): display(p,q);unassign('omega'); Using the cursor, the plot shows that the peak magnification is 2.05 dB, or NiMvJiUiTUc2IyUlcGVha0cpIiM1KiYtJSZGbG9hdEc2JCIkMCMhIiMiIiIiIz8hIiI= = 1.27, occurring at the servo's damped angular natural frequency. This is given by NiMvLSUmbG9nMTBHNiMqJiYlJm9tZWdhRzYjJSJkRyIiIiUiVEdGLC0lJkZsb2F0RzYkIiN0ISIj, i.e. NiMmJSZvbWVnYUc2IyUiZEc= = 5.4/second (the normalizing constant, T, having been chosen as 1 second). The plot also shows that feedback cutoff, or gain crossover, frequency is given by NiMvLSUmbG9nMTBHNiMqJiUmb21lZ2FHIiIiJSJUR0YpLSUmRmxvYXRHNiQiJC0iISIj, i.e. w = 10.5/second. The peak magnification and feedback cutoff (or gain crossover) frequency agree well with those found from the Nyquist diagram and the Bode attenuation diagram respectively. The peak magnification and the frequency at which it occurs can be found algebraically as follows. A function C/R(a) is obtained by substituting NiMvJSZhbHBoYUclJm9tZWdhRw== in C/R(f). From this function, the maximum value of C/R(f) can then be found thus,`C/R(alpha)` := map(abs, subs(omega=alpha,`C/R(f)`) ):`C/R max` = maximize(`C/R(alpha)`,location,alpha=3..10); This shows that the maximum value of C/R(f) is 1.262 and is located at NiMvKiYlJm9tZWdhRyIiIiUiVEdGJi0lJkZsb2F0RzYkIiVOYSEiJA==, demonstrating that using the cursor is sufficiently accurate for practical purposes. As measured with the cursor, the magnitude of the servo's frequency response falls to -3 dB at NiMmJSZvbWVnYUc2IyUiY0c= = 15.5/sec. This frequency is known as the servo's cutoff frequency (not to be confused with the feedback cutoff frequency). In this case, the servo's bandwidth extends from the very lowest input frequency - an input of "zero frequency" - up to the cutoff frequency. That is, the bandwidth is 15.5/sec. The precise value can be extracted from C/R(a), thus omega[c]*T = fsolve( `C/R(alpha)`=10^(-3/20), alpha,0..40); The greater a servo's bandwidth, the faster its response, but the less the attenuation of any noise in the input.unassign('omega','alpha');<Text-field layout="Heading 3" style="Heading 3">Response to Unit Step Displacement</Text-field> The servo's response to a unit step function of displacement is given by Output(s) = NiMvKiYlJnRoZXRhRyIiIi0lIm9HNiMlInNHRiYqKC0lJEMvUkdGKUYmRiZGJkYqISIi. Abbreviating C/R(s) to C/R, the servo's response in the time domain to a unit step function of displacement, NiMvKiYlJnRoZXRhRyIiIiUib0dGJiUiQ0c= = NiMqKCUkQy9SRyIiIkYlRiUlInNHISIi, is obtained using the inverse Laplace transform. Since the input is 1 radian, the servo's normalized response is obtained by dividing by 1 radian. The time axis is normalized by dividing by T = 1 sec, so that the plot is completely free of dimensions.`C/R step disp` := invlaplace(`C/R`/s,s,t):p := plot(`C/R step disp`/1,t=0..1,labels=["t/T","C/R"],
title="Normalized Response to Unit Step Displacement", colour=blue):q := plot([[0.02,0.08],[0.13,1.0]],style=LINE, colour=red):r := plot(1,t=0..1,linestyle=DOT, colour=black):display(p,q,r);
From the transient response plot above:a. the peak of the transient response is 1.21 (i.e. 21% overshoot), which is roughly the same as the peak magnification of the frequency response (1.27).b. the response time - i.e. time to reach peak of the transient - is 0.3 second. These can be obtained analyticallyzero := diff(`C/R step disp`,t):`Max overshoot at t/T` := fsolve(zero=0,t,0.1..0.4);`Max overshoot` := eval(`C/R step disp`,t=`Max overshoot at t/T`); This demonstrates that using the cursor is accurate enough for most practical purposes. The maximum slope of the response, judged by eye, is shown by the red line on the plot. The reciprocal of its slope gives an indication of the time for the response, i.e. C/R, to reach 1.0, and is known as the build-up, or rise, time.`build-up time` := evalf((0.13 - 0.02)/(1.0 - 0.09),3); This can be obtained analytically,inflect := diff(`C/R step disp`,t,t):tim := fsolve(inflect=0,t,0..0.2):`build-up time` := 1/(eval( diff(`C/R step disp`,t), t=tim)); Again, using the cursor is accurate enough.<Text-field layout="Heading 3" style="Heading 3">Response to Unit Step of Velocity</Text-field>The servo's response to a unit step input of velocity is given by Output(s) = NiMvKiYlJnRoZXRhRyIiIi0lIm9HNiMlInNHRiYqKC0lJEMvUkdGKUYmRiZGJiokRioiIiMhIiI=. Abbreviating, NiMmJSZ0aGV0YUc2IyUib0c= = NiMvJSJDRyooJSRDL1JHIiIiRidGJyokJSJzRyIiIyEiIg==. Whence the servo's response in the time domain is obtained using the inverse Laplace transform. Again, the time axis is normalized by dividing it by T = 1 sec, while the response is normalized by dividing it by the input's (angular) position, R. That is, NiMvJiUmdGhldGFHNiMlIm9HJSJDRw== is divided by NiMvJiUmdGhldGFHNiMlImlHJSJSRw== = velocity \327 t = 1 radian/sec \327 t sec = t radian, thus,`vel step resp` := (invlaplace(`C/R`/s^2,s,t))/ t:p := plot(`vel step resp`,t=0..3/2,labels=["t/T","C/R"],
title="Normalized Response to Unit Step of Input Velocity" ):q := plot(1,t=0..3/2,linestyle=DOT,colour=black):display(p,q); The plot shows that the servo's normalized (position) response to a unit step input of velocity approaches unity smoothly and without any overshoot, reaching 0.99 at t/T = 0.8. This fast, "pursuit" response is ideal for a tracking radar's servo that has to follow targets with rapidly changing velocities. The response, NiMvJiUmdGhldGFHNiMlIm9HJSJDRw==, to a unit step input of velocity is simply the normalized response shown above, multiplied by the normalizing factor R = 1 radian/sec \327 t sec, and is shown below.p := plot(`vel step resp`*t,t=0..1, colour=blue ):q := plot(1*t,t=0..1,colour=red):display(p,q,title = "Response (blue) to Unit Step Input (red) of Velocity",labels=["t/T","C & R"]);
The input = NiMvJiUmdGhldGFHNiMlImlHJSJSRw== = 1 radian/sec \327 t is shown by the straight line on the plot. The error when responding to a unit step input of velocity is NiMvJShlcHNpbG9uRywmJiUmdGhldGFHNiMlImlHIiIiJkYnNiMlIm9HISIi = NiMsJiUiUkciIiIlIkNHISIi and is shown belowplot( 1000*(1 - `vel step resp`)*t, t=0..3/2,title="Error of Response to Unit Step Input of Velocity", labels=["t/T","1000 \327 Error"],labeldirections=[horizontal,vertical] );
The error falls to a minimum of 5.8 milliradians (mils) between t/T = 0.9 and t/T = 1.1, and then increases slowly to its steady state value. Earlier this was shown to be NiMvJSlzc35lcnJvckcqKCUmT21lZ2FHIiIiJSNHYUchIiIlI0ttR0Yp = NiMqJiUmT21lZ2FHIiIiJSNLdkchIiI=. In this case, W = 1 radian/sec. `ss error` := evalf(eval(`ss error`,Omega=1)); The design requirements included a maximum angular velocity input of 0.5 radian/second, generated by an aircraft flying at 150 yards/second (307 mph), which corresponds to the aircraft passing 300 yards away from a radar antenna controlled by the servo. (See Section 6.12 in Chapter 6 of the Reference). As noted above, the error increases slightly after t/T = 0.9, so the worst case tracking error will be obtained by using the steady state tracking error and the maximum input of 0.5 radian/sec. Since the steady state error is the response to a unit step input of 1 radian/sec, the worst case tracking error with an input of 0.5 radian/second would have been:`Tracking Error` := evalf(`ss error`/2*(300*yards),3); In fact, the maximum tracking error at any range, r(t), is given by NiMqKC0lJk9tZWdhRzYjJSJ0RyIiIi0lInJHRiZGKCUjS3ZHISIi. Noting that the maximum angular input velocity occurs when the aircraft flies on a course perpendicular to the radar's beam, NiMvJi0lJk9tZWdhRzYjJSJ0RzYjJSRtYXhHKiYtJSJWR0YnIiIiLSUickdGJyEiIg==, whence maximum tracking error = NiMqJi0lIlZHNiMlInRHIiIiJSNLdkchIiI=. <Text-field layout="Heading 2" style="Heading 2">Observations</Text-field><Text-field layout="Heading 3" style="Heading 3">Servo</Text-field><Text-field layout="Heading 3" style="Heading 3">Summary</Text-field>a. The loop transfer function plotted on the Nyquist Diagram shows that the peak magnification of the frequency response, or resonance, is 1.27. This is the same as the approximate value given by equation 109 in Section 8.8 of the reference.b. The Bode charts of the loop gain and loop phase margin show a very big phase margin of 59 degrees at feedback cutoff (gain crossover), indicating a very stable closed loop system, even with the large value of the velocity error coefficient, Kv = 158/second. Kv is readily measured as the intercept of the initial -6dB/octave slope of the loop gain with the NiMvKiYlJm9tZWdhRyIiIiUiVEdGJkYm axis at 42dB = 158.c. The servo's closed loop frequency response indicates a peak magnification of 1.27, as does the Nyquist Diagram, and a damped angular natural frequency of 5.4/second.d. The servo's closed loop response to a unit step input shows a peak transient response of 1.21 (i.e. 21% overshoot), a build-up time of 0.12 second, and a response time of 0.3 second. Although the servo's transient response might appear somewhat overdamped, this was necessary to minimise the response to noise generated by the radar, while still obtaining the best response to the radar's error signal when tracking an aircraft. (The error signal corresponds to the difference between the axis of the radar's beam and the aircraft target.)e. The peak transient overshoot of 1.21 in response to a unit step input of position, is roughly the same as the peak magnification of the frequency response, 1.27.f. The servo's error, in response to a unit step velocity input, is ideally smooth and becomes constant at 6.33 milliradians after about 1.5 second. (This corresponds to 0.95 yards when tracking a WW II aircraft flying at 307 mph past the radar antenna situated 300 yards from the flight path.)g. A root-locus is not included in this worksheet, since it adds nothing to the information given by the other plots, and could be slightly misleading. It is, perhaps, worth remembering that the root-locus method was originally introduced into servomechanism analysis and design as a method for finding the roots of complicated transfer functions, and thence establishing the closed loop response to step inputs. A PC with suitable software, such as Maple, now performs that function very much better.<Text-field layout="Heading 3" style="Heading 3">Design Rule</Text-field>h. C.H. Dowker concluded that the best performance for a servo with proportional control would always be obtained by adjusting the loop gain (Kv), so that feedback cutoff (gain crossover) occurs one third the way down a long -6 dB/octave portion of the loop transfer function. This gives the best response to the signal, while minimising the effect of noise. (For adequate phase margin, the -6 dB/octave portion should always be at least 2.5 octaves long, and in this example, it was 3.5 octaves long.)<Text-field layout="Heading 3" style="Heading 3">Maple</Text-field>a. Maple makes light work of obtaining the loop transfer function and then plotting the Nyquist diagram, Bode diagrams, closed loop frequency response, and transient responses to step inputs of position and of velocity. It is all too easy to forget that before the advent of the PC and Maple, this work used to be both demanding and time consuming.b. The cursor provides a quick and adequately accurate means for obtaining values from the plots, for example maximum overshoot from the plot of transient response to a step function of input displacement. This is a very useful feature for servo design work.c. At small values of the independent variable, the "semilogplot" command produces rather distorted plots. This has been commented upon in Robert Israel's Maple Advisor. However, the problem can be overcome as demonstrated above. Similarly, the "maximize" command is not always successful - it failed when trying to obtain the maximum overshoot of the servo's response to a position step input. This was simply overcome as shown above. Robert Israel's Maple Advisor database contains a function that can be used in place of "maximize".d. Although the "coordplot(polar)" grid was used successfully for the Nyquist diagram above, it is somewhat restrictive, and a more versatile version would be useful for servo design and analysis.<Text-field layout="Heading 2" style="Heading 2">Reference</Text-field>Theory of Servomechanisms by Hubert M. James, Nathaniel B. Nicholls, and Ralph S. Phillips, Volume 25 of the Massachusetts Institute of Technology Radiation Laboratory Series, published by the McGraw-Hill Book Company, Inc. in 1947. See particularly Sections 8.3 - 8.5, end of 8.7, and
8.8 - 8.9 in Chapter 8.Legal Notice: The copyright for this application is owned by the author(s). Neither Maplesoft nor the author are responsible for any errors contained within and are not liable for any damages resulting from the use of this material. This application is intended for non-commercial, non-profit use only. Contact the author for permission if you wish to use this application in for-profit activities.