23-Chem-A6 Process Dynamics and Control · December 2014
Question 8 of 8: Bode Plot and Gain Margin of an FOPDT Loop
Nivaar worked solution (AI-drafted; not reviewed by a licensed engineer)
Notes on this paper
National Exams / EGBC — December 2014 — 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. The parts are quantitative throughout — two-capacitance thermal modelling, state-space transfer functions, IMC design with a right-half-plane zero, Nyquist stability of an open-loop-unstable plant, ultimate gain with sensor dead time, second-order damping regimes, a non-linear CSTR, and Bode gain-margin design — and every requested plot (step response, IMC servo response, Nyquist locus, damping family, Bode diagram) is drawn as a real figure.
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, feedback stability, frequency response and IMC design; G. Stephanopoulos, Chemical Process Control: An Introduction to Theory and Practice (Prentice Hall) — Nyquist and Bode stability, dead-time systems and controller design; D. R. Coughanowr & S. E. LeBlanc, Process Systems Analysis and Control (3rd ed., McGraw-Hill) — first- and second-order dynamics, linearisation of non-linear balances. Standard control conventions (deviation variables; unity valve/sensor gains unless stated) are used throughout.
Problem 8: Bode Plot and Gain Margin of an FOPDT Loop (20%)
Approach. Read the magnitude and phase asymptotes from the pole and all-pass delay, locate the phase-crossover frequency $\omega_{co}$ where $\angle L=-\pi$, then set the gain margin $=|L|^{-1}$ at $\omega_{co}$ to solve for $k_c$.
(a) Bode characteristics. $|L|=\dfrac{k_c}{\sqrt{1+0.25\omega^2}}$ (the dead time is all-pass) and $\varphi=-\arctan(0.5\omega)-0.1\omega$ (rad). Corner frequency $\omega=1/0.5=2$ rad/s. Low-frequency: AR$\to k_c$ (0 dB at $k_c=1$), slope 0, phase$\to0^\circ$. High-frequency: slope $-1$ ($-20$ dB/dec) from the single pole; phase falls without bound because the dead-time term $-0.1\omega$ keeps subtracting, crossing $-180^\circ$ at a finite frequency.
(b) Phase-crossover frequency. Set $\varphi=-\pi$: $\arctan(0.5\omega)+0.1\omega=\pi$. Solving numerically (Newton/bisection) gives $\omega_{co}\approx16.9$ rad/s. The proportional gain does not affect phase, so $\omega_{co}$ is fixed.
(b) Gain for GM = 1.7. $$\text{GM}=\frac{1}{|L(j\omega_{co})|}=\frac{\sqrt{1+0.25\omega_{co}^2}}{k_c}=1.7\ \Longrightarrow\ \boxed{k_c=\frac{\sqrt{1+0.25\omega_{co}^2}}{1.7}\approx5.0.}$$ A proportional gain of about 5 leaves a factor-of-1.7 cushion below the stability limit (the ultimate gain, GM$=1$, would be $k_{cu}\approx8.5$).
Problem 8(a): open-loop Bode plot (drawn at $k_c=1$). Magnitude is flat then breaks to $-20$ dB/dec at the corner $\omega=2$; phase starts at $0^\circ$ and, driven past $-90^\circ$ by the dead time, crosses $-180^\circ$ at $\omega_{co}\approx16.9$ — the frequency that sets the gain margin.