23-Chem-A6 Process Dynamics and Control · May 2018
Question 8 of 8: Nyquist Diagram and Gain Margin of $1/(s(s+2)^2)$
Nivaar worked solution (AI-drafted; not reviewed by a licensed engineer)
Notes on this paper
Reference texts: D. E. Seborg, T. F. Edgar, D. A. Mellichamp & F. J. Doyle III, Process Dynamics and Control (4th ed., Wiley) — Laplace-domain modelling, linearisation, transfer functions, feedback stability, frequency response and IMC design; G. Stephanopoulos, Chemical Process Control: An Introduction to Theory and Practice (Prentice Hall) — Routh test, Nyquist criterion, state-space models and dead-time systems; D. R. Coughanowr & S. E. LeBlanc, Process Systems Analysis and Control (3rd ed., McGraw-Hill) — first-order tank/thermal dynamics, step/ramp response and block-diagram algebra. Standard control conventions (deviation variables; unity valve/sensor gains unless stated) are used throughout. This sitting is an open-book paper of eight problems of which any five constitute a complete exam; all eight are solved here.
Problem 8: Nyquist Diagram and Gain Margin of $1/(s(s+2)^2)$ (20%)
Given. $G(s)=\dfrac{1}{s(s+2)^2}$ — an integrator ($1/s$) in series with a repeated first-order pair at $s=-2$; unity gain assumed in the loop for the margin.
Find. (1) the qualitative Nyquist locus with its asymptotes and real-axis crossing; (2) the gain margin.
Problem 8(1): coming up from $\mathrm{Im}\to-\infty$ along the vertical asymptote $\mathrm{Re}=-\tfrac14$ (as $\omega\to0^+$), the locus crosses the negative real axis at $-\tfrac{1}{16}$ (at $\omega=2$), then curls into the origin at $\omega\to\infty$ (approaching along the $-270^\circ$ direction).
Approach. Substitute $s=j\omega$, separate the denominator into real and imaginary parts, find the phase-crossover frequency where the imaginary part vanishes (locus on the real axis), evaluate $|G|$ there, and take $\mathrm{GM}=1/|G(j\omega_{pc})|$.
(1) Frequency response. $G(j\omega)=\dfrac{1}{j\omega\,(j\omega+2)^2}$. Expand $(j\omega+2)^2=(4-\omega^2)+4j\omega$; multiplying by $j\omega$ gives the denominator $$D(j\omega)=-4\omega^2+j\,\omega(4-\omega^2).$$ So $G=1/D$.
(1) Asymptotes and axis crossing. As $\omega\to0^+$, $\mathrm{Re}\,G\to-\tfrac14$ and $\mathrm{Im}\,G\to-\infty$ (vertical asymptote at $\mathrm{Re}=-0.25$). As $\omega\to\infty$, $|G|\to0$ (phase $\to-270^\circ$). The locus meets the real axis where $\mathrm{Im}\,D=0$: $\omega(4-\omega^2)=0\Rightarrow \omega_{pc}=2$ (the finite crossover).
(2) Value at phase crossover. At $\omega=2$, $D=-4(2)^2+j\,2(4-4)=-16$, so $$G(j2)=\frac{1}{-16}=-0.0625\quad(\text{real, negative; phase }=-180^\circ).$$
(2) Gain margin. $$\boxed{\mathrm{GM}=\frac{1}{|G(j\omega_{pc})|}=\frac{1}{1/16}=16\ (=24.1\ \text{dB}).}$$ The open-loop gain can be increased $16\times$ before the crossing reaches $-1$ and the closed loop goes marginally stable (at which point $k_{cu}=16$, $\omega_u=2$).