NivaarExam PrepOfficial exam papers ↗

20-Bio-A2 Process Dynamics and Control · Undated paper

Question 7 of 8: Internal Model Control Design for a Dead-Time-Dominant Process

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

Notes on this paper

National Exams / EGBC — May 2019 — 04-Bio-A2 Process Dynamics and Control (the paper's own header line alternates between the codes "04-BIO-A2" and "04-BIO-A3" across pages; the subject line and content are unambiguously Process Dynamics and Control, so 04-Bio-A2 is used throughout). 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 the effect of a real zero on step-response overshoot, the Nyquist criterion for an open-loop-unstable process, a two-tank linear-vs-nonlinear-valve comparison, state-space-to-transfer-function conversion, Routh–Hurwitz stability with PI control of an unstable process, inverse Laplace transforms with time delay, Internal Model Control (IMC) design for a dead-time process, and Bode/gain-margin design.

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 7: Internal Model Control Design for a Dead-Time-Dominant 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{10e^{-5s}}{200s+1}$, IMC filter time constant $\tau_c=10$, perfect model $\tilde G_p=G_p$.

Find. (a) the block diagram and $G_c^{*}(s)$; (b) the unit-step closed-loop response $C(t)$ for the perfect-model case.

Gc*(s)(200s+1)/[10(10s+1)]Gp(s)10e^-5s/(200s+1)Gp~(s)(model, =Gp)-Σysp(s)y(s)c*(s)c*(s)feedback = y(s) minus model prediction (=0, perfect model)
Problem 7(a): standard IMC structure. $G_c^{*}(s)$ drives both the true process and an internal copy of the model; with a perfect model the feedback path carries zero signal, so the closed loop reduces to open-loop $G_c^{*}(s)G_p(s)$.

Approach. Factor $G_p$ into the non-invertible dead-time part $G_p^-=e^{-5s}$ and the invertible minimum-phase part $G_p^+=10/(200s+1)$; invert $G_p^+$ and append the first-order filter $f(s)=1/(\tau_cs+1)$ to form $G_c^{*}$; then, with a perfect model, the servo transfer reduces exactly to $G_p^-f$ and inverts as a delayed first-order step.

  1. (a) Factor the model. $G_p(s)=\underbrace{e^{-5s}}_{G_p^-}\cdot\underbrace{\dfrac{10}{200s+1}}_{G_p^+}$.
  2. (a) IMC controller. With $\tau_c=10\ (\ne\tau_p=200$, unlike the convenient special case where they coincide$)$: $$\boxed{G_c^{*}(s)=\big[G_p^+\big]^{-1}f(s)=\frac{200s+1}{10}\cdot\frac{1}{10s+1}=\frac{200s+1}{10(10s+1)}}$$ — a lead-dominant (PI-like) filtered inverse, not a constant, because the chosen filter time constant is much faster than the process lag.
  3. (b) Perfect-model servo transfer. With $\tilde G_p=G_p$ the internal feedback signal is exactly zero, so the loop reduces to $Y/Y_{sp}=G_p\,G_c^{*}$: substituting and cancelling $(200s+1)$, $$\boxed{\frac{Y(s)}{Y_{sp}(s)}=\frac{10e^{-5s}}{200s+1}\times\frac{200s+1}{10(10s+1)}=\frac{e^{-5s}}{10s+1}}$$ — a pure dead time cascaded with a first-order filter of the chosen $\tau_c=10$, exactly as IMC theory predicts for a perfect model (the achievable closed loop is limited only by the non-invertible dead time).
  4. (b) Invert the unit step. Ignoring the delay momentarily, $\mathcal L^{-1}\!\left[\dfrac1{s(10s+1)}\right]=1-e^{-t/10}$; reinserting the $5\ \text{s}$ delay as a pure time shift, $$\boxed{C(t)=\begin{cases}0,&t<5\ \text{s}\\[2pt]1-e^{-(t-5)/10},&t\ge5\ \text{s}.\end{cases}}$$ At $t=15\ \text{s}$ (one time constant after the delay clears): $C=1-e^{-1}=0.632$; at $t=30\ \text{s}$: $C=1-e^{-2.5}=0.918$; as $t\to\infty$: $C\to1$ (zero offset, as expected for a step in set point with a perfect model).
ResultValue
Model split$G_p^-=e^{-5s}$, $G_p^+=10/(200s+1)$
IMC controller$G_c^{*}(s)=(200s+1)/[10(10s+1)]$
Servo transfer (perfect model)$Y/Y_{sp}=e^{-5s}/(10s+1)$
Step response$C(t)=0$ for $t<5$; $1-e^{-(t-5)/10}$ for $t\ge5$