NivaarExam PrepOfficial exam papers ↗

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%)

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. $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.

P8: Nyquist of $G=1/(s(s+2)^2)$ ($\omega>0$ branch) asymptote Re$=-\tfrac14$ -0.25-0.125ImRe $\omega=2:\ -1/16$ $\omega\to0^+$ (Im$\to-\infty$) $\omega\to\infty:0$ plot symmetric about the real axis ($\omega<0$ is the mirror image)
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. (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$.
  2. (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).
  3. (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).$$
  4. (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$).
QuantityValue
Low-frequency asymptote$\mathrm{Re}=-\tfrac14$, $\mathrm{Im}\to-\infty$
Phase-crossover frequency$\omega_{pc}=2$ rad/time
Real-axis crossing$G=-1/16=-0.0625$
Gain margin$\mathrm{GM}=16$ ($24.1$ dB)
Ultimate gain / frequency$k_{cu}=16,\ \omega_u=2$
Back to the paper →