NivaarExam PrepOfficial exam papers ↗

20-Bio-A2 Process Dynamics and Control · December 2019

Question 6 of 8: Bode Plot and Gain-Margin Design for a Dead-Time-Plus-Lag Process

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

Notes on this paper

National Exams / EGBC — December 2019 — 04-Bio-A2 Process Dynamics and Control. Three-hour open-book examination; any non-communicating calculator is permitted. The cover page states that only the first five questions in the answer book are marked, yet all eight problems are printed on the paper — all eight are solved below for completeness, since the full set is a study resource. Content spans zero-location effects on step response (overshoot/inverse response), state-space-to-transfer-function conversion, first-order sensor dynamics under a triangular forcing function, Internal Model Control (IMC) design for a dead-time process, the Nyquist stability criterion for an open-loop-unstable process, Bode/gain-margin design, a linear draining-tank model, and Routh–Hurwitz stability with a PI controller.

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, transfer functions, Routh stability, Nyquist/Bode frequency response, and Internal Model Control design; G. Stephanopoulos, Chemical Process Control: An Introduction to Theory and Practice (Prentice Hall) — the Nyquist criterion for open-loop-unstable processes and dead-time systems. Standard control conventions (deviation variables; unity valve/sensor gain unless stated) are used throughout.

Problem 6: Bode Plot and Gain-Margin Design for a Dead-Time-Plus-Lag Process (20% total)

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_p(s)=\dfrac{e^{-0.1s}}{0.5s+1}$ ($\theta=0.1\ \text{s}$ dead time, $\tau=0.5\ \text{s}$ lag), proportional controller $K_c$.

Find. (a) the qualitative Bode magnitude and phase plots (slopes, corner frequency, extreme values); (b) $K_c$ such that the gain margin $=1.7$.

log10(ω) [ω in rad/s]|G(jω)| (dB)ωc=2 rad/s
Bode magnitude: flat at $0\ \text{dB}$ for $\omega\ll2\ \text{rad/s}$ (dead time never attenuates magnitude), corner frequency $\omega_c=1/\tau=2\ \text{rad/s}$, rolling off at $-20\ \text{dB/decade}$ beyond the corner.
log10(ω) [ω in rad/s]Phase (deg)ω_pc=16.89 rad/s, phase=-180°
Bode phase: starts at $0^{\circ}$ ($\omega\to0$), passes $-90^{\circ}$ near the corner frequency, and (unlike a pure lag) keeps decreasing without bound as $\omega\to\infty$ because the dead-time contributes an unbounded $-0.1\omega\ \text{rad}$ term — crossing $-180^{\circ}$ at $\omega_{pc}=16.89\ \text{rad/s}$.

Approach. Write exact magnitude/phase of $G_p(j\omega)$ (the dead time contributes phase only), read off the corner frequency and asymptotic slopes for the qualitative Bode plot, solve the transcendental phase-crossover condition $\angle G_p(j\omega_{pc})=-180^{\circ}$ numerically, then use $\text{GM}=1/[K_c|G_p(j\omega_{pc})|]$ to solve for $K_c$.

  1. (a) Frequency response. $$|G_p(j\omega)|=\frac{1}{\sqrt{(0.5\omega)^2+1}},\qquad \angle G_p(j\omega)=-0.1\omega\ (\text{rad})-\arctan(0.5\omega).$$
  2. (a) Extremes and corner. At $\omega=0$: $|G_p|=1$ ($0\ \text{dB}$), phase $=0^{\circ}$. Corner frequency (lag break point) $\omega_c=1/0.5=2\ \text{rad/s}$: magnitude flat ($0\ \text{dB/decade}$) below $\omega_c$, rolling off at $-20\ \text{dB/decade}$ above it. As $\omega\to\infty$: $|G_p|\to0$, and phase $\to-\infty$ (unbounded) — the dead-time term $-0.1\omega\ \text{rad}=-5.73\omega\ \text{deg}$ dominates and never saturates, unlike the lag's phase which saturates at $-90^{\circ}$.
  3. (b) Phase crossover $\omega_{pc}$. Solve $-0.1\omega-\arctan(0.5\omega)=-\pi$ numerically (transcendental, Newton/bisection): $$\boxed{\omega_{pc}=16.887\ \text{rad/s}}\quad(\text{check: }-0.1(16.887)-\arctan[0.5(16.887)]=-180.0^{\circ}\ \checkmark).$$
  4. (b) Magnitude at $\omega_{pc}$ (dead time contributes no attenuation). $$|G_p(j\omega_{pc})|=\frac{1}{\sqrt{(0.5\times16.887)^2+1}}=0.1176.$$
  5. (b) Solve for $K_c$. Gain margin $\text{GM}=\dfrac{1}{K_c|G_p(j\omega_{pc})|}$; setting $\text{GM}=1.7$: $$\boxed{K_c=\frac{1}{1.7\times0.1176}=5.00.}$$ (Check: $K_c|G_p(j\omega_{pc})|=5.00\times0.1176=0.588=1/1.7$ $\checkmark$.)
ResultValue
Corner frequency$\omega_c=2\ \text{rad/s}$; slope $0\to-20\ \text{dB/decade}$
DC / high-freq magnitude$0\ \text{dB}$ at $\omega=0$; $\to-\infty$ as $\omega\to\infty$
Phase behaviour$0^{\circ}\to-\infty$ (unbounded, dominated by dead time)
Phase-crossover frequency$\omega_{pc}=16.887\ \text{rad/s}$
$K_c$ for GM $=1.7$$K_c=5.00$