NivaarExam PrepOfficial exam papers ↗

20-Bio-A2 Process Dynamics and Control · May 2018

Question 6 of 8: Internal Model Control of a Dead-Time-Dominant Process

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

Notes on this paper

National Exams / EGBC — May 2018 — 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 second-order Bode/phase-margin analysis, feedback stability with a non-minimum-phase sensor, dead-time Nyquist/gain-margin design, zero-location effects on step response, nonlinear radiative-heat-transfer linearization, Internal Model Control (IMC) of a dead-time process, PI-controller stability/response, and nonlinear-CSTR linearization. Note: Problem 2's printed sub-part weights (10%+20%=30%) exceed its stated 20% problem total — both sub-parts are fully answered below.

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, linearization and IMC design; G. Stephanopoulos, Chemical Process Control: An Introduction to Theory and Practice (Prentice Hall) — Nyquist criterion and dead-time systems. Standard control conventions (deviation variables; unity valve gain unless stated) are used throughout.

Problem 6: Internal Model Control of 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{5e^{-10s}}{10s+1}$, IMC filter time constant $\tau_c=20\ \text{s}$, perfect model $\tilde G_p=G_p$, no Padé approximation of the dead time.

Find. (a) IMC controller $q(s)=G_c^{*}$, the classical feedback-equivalent $G_c(s)$, and whether $G_c$ is PID form; (b) the unit-step closed-loop response $\delta C(t)$.

q(s)IMC controllerGp(s)(actual process)Gp̃(s)(process model)−Σysp(s)y(s)u(s)u(s)feedback = y(s) minus model prediction
Problem 6(a): standard IMC structure. $q(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 loop reduces to open-loop $q(s)G_p(s)$.

Approach. Factor $G_p$ into a non-invertible part $G_p^-$ (the dead time; the gain is $+5$ so there is no RHP zero, only the delay is non-invertible) and an invertible minimum-phase part $G_p^+$; invert $G_p^+$ and append a first-order filter $f(s)=1/(\tau_cs+1)$ to form $q(s)$; convert to the classical equivalent $G_c=q/(1-G_pq)$ without any Padé approximation; then, with a perfect model, the servo transfer reduces to $G_p^-f$ and can be inverted directly (delay reinserted as a pure time shift).

  1. (a) Factor the model. $G_p(s)=\underbrace{e^{-10s}}_{G_p^-,\ G_p^-(0)=1}\cdot\underbrace{\dfrac{5}{10s+1}}_{G_p^+}$; only the dead time is non-invertible (no RHP zero here).
  2. (a) IMC controller. With a first-order filter (matching the first-order $G_p^+$) and $\tau_c=20$, $$\boxed{q(s)=G_c^{*}=\big[G_p^+\big]^{-1}f(s)=\frac{10s+1}{5}\cdot\frac{1}{20s+1}=\frac{10s+1}{5(20s+1)}.}$$
  3. (a) Classical feedback-equivalent. With a perfect model, $1-G_pq=1-\dfrac{e^{-10s}}{20s+1}$, so $$\boxed{G_c(s)=\frac{q(s)}{1-G_p(s)q(s)}=\frac{10s+1}{5\big[20s+1-e^{-10s}\big]}.}$$ Because $e^{-10s}$ appears transcendentally in the denominator (no Padé approximation applied), $G_c(s)$ is not reducible to a rational $K_c(1+1/(\tau_Is)+\tau_Ds)$ — it is a Smith-predictor-type controller, not classical PID.
  4. (b) Perfect-model servo transfer. Since $\tilde G_p=G_p$, the internal feedback signal is exactly zero and the loop reduces to $Y/Y_{sp}=G_p^-f$: $$\boxed{\frac{Y(s)}{Y_{sp}(s)}=e^{-10s}\cdot\frac{1}{20s+1}.}$$
  5. (b) Invert the unit step. Ignoring the delay momentarily, $\dfrac{1}{s(20s+1)}=\dfrac1s-\dfrac{1}{s+0.05}$, giving the undelayed response $1-e^{-t/20}$. Reinserting the $10\ \text{s}$ delay as a pure time shift, $$\boxed{\delta C(t)=\begin{cases}0,&t\lt10\\[2pt]1-e^{-(t-10)/20},&t\ge10.\end{cases}}$$
ResultValue
Model split$G_p^-=e^{-10s}$, $G_p^+=5/(10s+1)$
IMC controller$q(s)=(10s+1)/[5(20s+1)]$
Classical equivalent$G_c=(10s+1)/\{5[20s+1-e^{-10s}]\}$ — not PID form
Servo transfer (perfect model)$Y/Y_{sp}=e^{-10s}/(20s+1)$
Step response$\delta C(t)=0$ for $t\lt10$; $1-e^{-(t-10)/20}$ for $t\ge10$