23-Chem-A6 Process Dynamics and Control · May 2015
Question 5 of 8: IMC Design for a First-Order-Plus-Dead-Time Process
Nivaar worked solution (AI-drafted; not reviewed by a licensed engineer)
Notes on this paper
National Exams / EGBC — May 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 5: IMC Design for a First-Order-Plus-Dead-Time Process (20%)
Given. First-order-plus-dead-time (FOPDT) process $G_p=\dfrac{5e^{-2s}}{10s+1}$: gain $K_p=5$, lag $\tau_p=10\ \mathrm s$, dead time $\theta=2\ \mathrm s$; IMC filter $f=\dfrac{1}{\lambda s+1}$, $\lambda=20\ \mathrm s$.
Find. (a) $G_c^{*}$, the classical $G_c$, and whether $G_c$ is PID; (b) the perfect-model servo response.
Problem 5(a): IMC structure. Controller $G_c^{*}$ drives the real process $G_p$; the internal model $\tilde G_p$ predicts $\tilde C$, and the model error $\hat d=C-\tilde C$ is fed back. With a perfect model $\tilde G_p=G_p$, $\hat d=0$ and the servo response reduces to $C/R=G_pG_c^{*}$.
Approach. Factor the model into invertible and all-pass (dead-time) parts, invert only the invertible part and append the filter to build $G_c^{*}$, convert to the classical form $G_c=G_c^{*}/(1-\tilde G_pG_c^{*})$, then evaluate the perfect-model servo transfer function and invert the step.
Factor the model. $G_p=G_p^{-}\,G_p^{+}$ with the invertible part $G_p^{-}=\dfrac{5}{10s+1}$ and the all-pass (non-invertible) dead time $G_p^{+}=e^{-2s}$ (unity gain, $G_p^{+}(0)=1$).
(a) IMC controller. Invert the good part and add the first-order filter: $$\boxed{G_c^{*}(s)=\big(G_p^{-}\big)^{-1}f=\frac{10s+1}{5}\cdot\frac{1}{20s+1}=\frac{10s+1}{5(20s+1)}}.$$ The dead time is deliberately not inverted (that would demand prediction).
(a) Equivalent classical controller. $G_c=\dfrac{G_c^{*}}{1-\tilde G_pG_c^{*}}$; with $\tilde G_p=G_p$ and $G_pG_c^{*}=\dfrac{e^{-2s}}{20s+1}$, $$\boxed{G_c(s)=\frac{10s+1}{5\left(20s+1-e^{-2s}\right)}}.$$ Because $e^{-2s}$ appears in the denominator, $G_c$ is not a PID controller (it is a rational-plus-dead-time compensator, exact only without a Padé approximation).
(b) Perfect-model servo transfer function. With no model error the closed loop collapses to $$\frac{C(s)}{R(s)}=G_pG_c^{*}=G_p^{+}f=\frac{e^{-2s}}{20s+1}.$$ The filter $\lambda=20\ \mathrm s$ sets the closed-loop speed; the dead time $e^{-2s}$ is unavoidable.
(b) Step response. For $R=1/s$, $C(s)=\dfrac{e^{-2s}}{s(20s+1)}$. Inverting (first-order step, delayed by $2\ \mathrm s$): $$\boxed{\delta C(t)=\Big[\,1-e^{-(t-2)/20}\,\Big]\,\mathcal U(t-2)}.$$ The output stays at $0$ until $t=2\ \mathrm s$, then rises exponentially with time constant $20\ \mathrm s$ to the set point — no offset (integral-like action of IMC with $G_p^{+}(0)=1$) and no overshoot.