23-Chem-A6 Process Dynamics and Control · December 2015
Question 8 of 8: Bode Plot and Gain-Margin Design for a FOPDT Process
Nivaar worked solution (AI-drafted; not reviewed by a licensed engineer)
Notes on this paper
National Exams / EGBC — December 2015 — 04-Chem-A6 Process Dynamics & Control. Three-hour open-book examination; any non-communicating calculator is permitted. Eight problems are printed and any five constitute a complete paper (each worth 20%); all eight are solved below for completeness. Most parts are quantitative (dynamic modelling, transfer functions, step/pulse/ramp responses, Routh and Nyquist stability, and IMC design); qualitative sketches are drawn as real figures where the paper asks for them.
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 and dead-time systems; D. R. Coughanowr & S. E. LeBlanc, Process Systems Analysis and Control (3rd ed., McGraw-Hill) — first-order thermal/level dynamics, step/ramp/pulse response and block-diagram algebra. Standard control conventions (deviation variables; unity valve/sensor gains unless stated) are used throughout.
Problem 8: Bode Plot and Gain-Margin Design for a FOPDT Process (20%)
Given. A first-order-plus-dead-time process: time constant $\tau=0.5$ s (corner $\omega=1/\tau=2$ rad/s), dead time $\theta=0.1$ s, unity process gain; $L(j\omega)=k_c\,e^{-0.1j\omega}/(0.5j\omega+1)$.
Find. (a) the Bode asymptotes and extremes; (b) $k_c$ for $\text{GM}=1.7$.
Approach. Read the magnitude asymptotes and phase limits directly from the FOPDT form, then locate the phase-crossover frequency ($\angle L=-180^\circ$) and invert the gain-margin definition to solve for $k_c$.
(a) Bode shape. The normalised magnitude $|G_p|=1/\sqrt{1+0.25\omega^2}$ is flat at $0$ dB (value $k_c$) below the corner $\omega=2$ rad/s, then rolls off at $-20$ dB/decade. The phase $\phi=-0.1\omega-\tan^{-1}(0.5\omega)$ starts at $0^\circ$; the first-order lag contributes $-45^\circ$ at $\omega=2$ and tends to $-90^\circ$, while the dead-time term $-0.1\omega$ keeps subtracting phase, so $\phi\to-\infty$ — the signature of a transport lag.
Problem 8(a): Bode plot of $G_p=e^{-0.1s}/(0.5s+1)$. Magnitude (blue) with $0$ dB and $-20$ dB/dec asymptotes (dashed) and corner $\omega=2$; the phase crosses $-180^\circ$ at $\omega_{co}=16.9$ rad/s and continues toward $-\infty$ because of the dead time.
(b) Phase crossover. Solve $\angle L=-0.1\omega-\tan^{-1}(0.5\omega)=-\pi$, i.e. $0.1\omega+\tan^{-1}(0.5\omega)=\pi$. Because $\tan^{-1}$ saturates near $\pi/2$, most of the phase comes from the delay, giving $\omega_{co}=16.9$ rad/s.
Amplitude ratio there. $\dfrac{|L|}{k_c}=\dfrac{1}{\sqrt{1+0.25\cdot16.9^2}}=0.1176$.
Solve for the gain. The gain margin is $\text{GM}=\dfrac{1}{|L(j\omega_{co})|}=\dfrac{1}{k_c\cdot0.1176}$. Setting $\text{GM}=1.7$: $$\boxed{k_c=\frac{1}{1.7\times0.1176}=5.0}.$$ (The ultimate gain, $\text{GM}=1$, would be $k_{cu}=1/0.1176=8.5$; a gain margin of $1.7$ backs the design off to $k_c=5.0$.)