NivaarExam PrepOfficial exam papers ↗

22-Elec-B2 Advanced Control Systems · December 2013

Question 6 of 6: Proportional design to a specified gain margin for a plant with transport delay

Nivaar worked solution (AI-drafted; not reviewed by a licensed engineer)

Notes on this paper

Paper format. National Examinations, 07-Elec-B2 Advanced Control Systems, December 2013 — 3 hours, closed book. Six questions; the rubric states that “any four questions constitute a complete paper” and that “all questions are of equal value”, so each is worth 25 marks. Tables of Laplace and z-transforms are appended to the paper. Only a Casio FX-991 or Sharp EL-540 calculator is permitted.

Reference texts. G. F. Franklin, J. D. Powell and A. Emami-Naeini, Feedback Control of Dynamic Systems, 7th ed., Pearson (Routh–Hurwitz, stability margins, steady-state error, frequency response); K. J. Åström and R. M. Murray, Feedback Systems: An Introduction for Scientists and Engineers, 2nd ed., Princeton (sensitivity function, loops with transport delay, non-minimum-phase limitations); K. Ogata, Modern Control Engineering, 5th ed., Pearson (state-space realisations, controllability and observability, pole placement); G. F. Franklin, J. D. Powell and M. L. Workman, Digital Control of Dynamic Systems, 3rd ed., Ellis-Kagle (zero-order-hold equivalents, deadbeat design, discrete identification). These are the works listed by Engineers Canada / EGBC for the Elec-B2 syllabus.

Scope of this solution. Although the paper is marked on four questions, all six are worked here so that the set serves as a complete study resource. Every boxed number is recomputed from the question’s own data for this paper.

Check: two printing defects in the source paper. In Question 3 the leading entry of the second row of the system matrix is missing from the print — the row reads “ 0  −1” with the first column blank — and the equation is typeset as “$\dots x + B = [0\ 1\ 1]^{\mathsf T}u$” rather than $\dot x = Ax + Bu$. Reading the column alignment against rows 1 and 3, the missing entry is $a_{21}=0$; that reading is adopted throughout and makes the system controllable, so the design asked for is well posed. In Question 4(b) the paper asks for “$T(z)$ that relates $X(z)$ to $R(z)$”, but no signal $X$ is defined anywhere in the question or its figure; $X(z)$ is read as $Y(z)$.

Question 6: Proportional design to a specified gain margin for a plant with transport delay (25 marks)

Question text not reproduced: the examination questions are © Engineers and Geoscientists BC. Open the official past paper (linked at the top of this page) to read the question, then follow the worked solution below.

Given. A type-1 plant with one lag and a pure transport delay, to be closed with a proportional gain against a gain-margin specification.

Given data — Question 6
QuantityValue
Plant$P(s) = \dfrac{e^{-0.2s}}{s(2s+1)}$
Transport delay$\theta = 0.2$ s
Lag time constant$T = 2$ s (pole at $s=-0.5$)
Controller$C(s) = K$, proportional only
Specificationgain margin $= 8$ dB
Feedbackunity negative

Find. (a) the gain $K$ meeting the 8 dB gain margin; (b) the phase margin that results; (c) the steady-state output to a unit step set-point.

Approach. Locate the phase crossover from the phase equation alone — the delay contributes phase but no attenuation, so the crossover frequency does not depend on $K$ — then scale the magnitude there to the required margin, and finally find the new gain crossover to read the phase margin.

  1. Write the frequency response. A pure delay has unit magnitude at every frequency and a phase that grows linearly with frequency, so $$|K P(j\omega)| = \frac{K}{\omega\sqrt{1+4\omega^{2}}},\qquad \angle KP(j\omega) = -\,\theta\omega - 90^\circ - \arctan(2\omega),$$ with $\theta\omega$ expressed in radians before conversion. Because the delay is transcendental, the closed-loop characteristic equation $s(2s+1)+Ke^{-0.2s}=0$ is not a polynomial and the Routh test is unavailable — the frequency-domain route is the only closed-form one.
  2. Locate the phase crossover. Setting the phase to $-180^\circ$ and working in radians, $$0.2\,\omega + \arctan(2\omega) = \frac{\pi}{2}.$$ The left side rises monotonically from 0, so the root is unique; solving numerically, $$\omega_{c\varphi} = 1.5553\ \text{rad/s}.$$ Crucially this frequency is set by the phase alone and is therefore independent of $K$ — a proportional gain slides the magnitude curve up and down but never moves the $-180^\circ$ crossing.
  3. Evaluate the plant magnitude there. $$|P(j\omega_{c\varphi})| = \frac{1}{1.5553\sqrt{1+4(1.5553)^{2}}} = \frac{1}{1.5553 \times 3.2674} = 0.19679,$$ so the uncompensated loop ($K=1$) already has a gain margin of $-20\log_{10}(0.19679) = 14.12$ dB. Meeting an 8 dB specification therefore means increasing the gain.
  4. Solve for the required gain. A gain margin of 8 dB means the loop magnitude at the phase crossover must equal $10^{-8/20}$: $$K\,|P(j\omega_{c\varphi})| = 10^{-8/20} = 0.39811 \quad\Longrightarrow\quad K = \frac{0.39811}{0.19679}$$ $$\boxed{K = 2.0230}$$ Checking, $20\log_{10}(2.0230 \times 0.19679) = -8.00$ dB exactly.
  5. Find the new gain crossover. Setting $|KP(j\omega)|=1$, $$\omega\sqrt{1+4\omega^{2}} = K = 2.0230 \quad\Longrightarrow\quad \omega_{cg} = 0.9456\ \text{rad/s}.$$
  6. Read the phase margin. Evaluating the phase at that frequency, in degrees, $$\angle KP(j\omega_{cg}) = -\frac{180}{\pi}(0.2)(0.9456) - 90^\circ - \arctan\!\big(2(0.9456)\big) = -10.84^\circ - 90^\circ - 62.13^\circ = -162.97^\circ,$$ $$\boxed{\text{PM} = 180^\circ - 162.97^\circ = 17.0^\circ \ \text{at } \omega_{cg}=0.946\ \text{rad/s}}$$ The delay costs almost 11° of that phase on its own. A 17° margin is lightly damped — roughly $\zeta \approx 0.17$ by the usual $\zeta \approx \text{PM}/100$ rule of thumb — which is the price of insisting on an 8 dB rather than a larger gain margin on a delayed plant.
  7. Determine the steady-state output. The plant contains a free integrator, so the loop is type 1 and the position error constant is infinite: $$K_{pos} = \lim_{s\to0} KP(s) = \lim_{s\to0}\frac{K e^{-0.2s}}{s(2s+1)} = \infty \quad\Longrightarrow\quad e_{ss} = \frac{1}{1+K_{pos}} = 0 .$$ Therefore the output tracks the unit step exactly in steady state: $$\boxed{y(\infty) = 1, \qquad e_{ss} = 0}$$ This holds for any stable $K$, so it is a property of the plant structure rather than of the design just completed.
-60-40-2002040|L| (dB)0.1110gain crossover 0.946 rad/sGM = 8.00 dB-360-315-270-225-180-135-90frequency (rad/s), log scalephase of L (deg)0.1110phase crossover 1.555 rad/sPM = 17.0 degOpen loop KP(s) with K = 2.0230, plant e^(−0.2s)/[s(2s+1)]
Figure Q6.1 — Bode plot of the compensated loop. The delay contributes phase only, so the −180° crossing sits at ω = 1.555 rad/s regardless of K; choosing K = 2.0230 pulls the magnitude there down to −8 dB, which is the required gain margin.

Check: the delay is handled exactly, not by approximation. All numbers above use $e^{-j0.2\omega}$ directly. A first-order Padé substitution $e^{-0.2s}\approx(1-0.1s)/(1+0.1s)$ would make the characteristic equation a polynomial and allow a Routh solution, but it shifts the phase crossover by a few per cent and with it the design gain; the exact frequency-domain calculation is preferable and is what the mark scheme expects here. The $17^\circ$ phase margin also means the closed-loop step response will be markedly oscillatory: if the specification really is on transient quality rather than on gain margin, this plant wants a lead or a PI-lead compensator rather than a bare proportional gain.

Final results — Question 6
QuantitySymbolValue
Phase crossover frequency$\omega_{c\varphi}$1.5553 rad/s
Plant magnitude there$|P(j\omega_{c\varphi})|$0.19679
Uncompensated gain margin—14.12 dB
Design gain for GM = 8 dB$K$2.0230
Gain crossover frequency$\omega_{cg}$0.9456 rad/s
Phase marginPM17.0°
Steady-state output, unit step$y(\infty)$1 (zero error, type-1 loop)
Back to the paper →