NivaarExam PrepOfficial exam papers ↗

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

Question 7 of 8: Stability and Step Response of a PI-Controlled First-Order 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 7: Stability and Step Response of a PI-Controlled First-Order 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)=\dfrac1{s+5}$, $G_c(s)=k_c\big(1+\tfrac1s\big)$ (the exam text labels this "proportional-derivative (PI)" — the acronym and the given equation are both for a Proportional-Integral controller; the equation, not the mismatched label, is what is solved here).

Find. (a) the range of $k_c$ for closed-loop stability; (b) the unit-step closed-loop response with $k_c=1$.

Approach. Form the open-loop $L(s)=G_cG_p$, clear fractions to a quadratic characteristic equation, apply the (trivial, for a quadratic) Routh positive-coefficient test, then invert the closed-loop step response by partial fractions at $k_c=1$.

  1. (a) Characteristic equation. $L(s)=\dfrac{k_c(s+1)}{s(s+5)}$, so $1+L=0\Rightarrow s(s+5)+k_c(s+1)=0$, i.e. $$\boxed{s^2+(5+k_c)s+k_c=0.}$$
  2. (a) Stability. A quadratic is stable iff both coefficients are positive: $5+k_c\gt0$ and $k_c\gt0$; the second condition already implies the first, so $$\boxed{k_c\gt0.}$$
  3. (b) Poles at $k_c=1$. $s^2+6s+1=0\Rightarrow s=-3\pm2\sqrt2=-0.1716,\,-5.8284$ (both real and negative — an overdamped, non-oscillatory closed loop).
  4. (b) Closed-loop TF and step response. $Y/Y_{sp}=\dfrac{k_c(s+1)}{s^2+(5+k_c)s+k_c}\Big|_{k_c=1}=\dfrac{s+1}{s^2+6s+1}$. With $Y_{sp}(s)=1/s$, partial fractions give residues $A=1$, $B=-0.8536$, $C=-0.1464$, so $$\boxed{y(t)=1-0.8536\,e^{-0.1716t}-0.1464\,e^{-5.8284t}.}$$
  5. Check. $y(0)=1-0.8536-0.1464=0$ (correct start) and $y(\infty)=1$ (zero steady-state offset), as expected from the controller's integral action acting on a type-0 process.
ResultValue
Stability range$k_c\gt0$
Closed-loop poles ($k_c=1$)$s=-0.1716,\,-5.8284$
Step response$y(t)=1-0.8536e^{-0.1716t}-0.1464e^{-5.8284t}$
Steady-state value$y(\infty)=1$ (zero offset)