20-Bio-A2 Process Dynamics and Control · December 2019
Question 2 of 8: Transfer Function and Step Response From a State-Space Model
Nivaar worked solution (AI-drafted; not reviewed by a licensed engineer)
Notes on this paper
National Exams / EGBC — December 2019 — 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 zero-location effects on step response (overshoot/inverse response), state-space-to-transfer-function conversion, first-order sensor dynamics under a triangular forcing function, Internal Model Control (IMC) design for a dead-time process, the Nyquist stability criterion for an open-loop-unstable process, Bode/gain-margin design, a linear draining-tank model, and Routh–Hurwitz stability with a PI controller.
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 2: Transfer Function and Step Response From a State-Space Model (20%)
Given. Two-state linear ODE system with $\dot x_1=-2.4048x_1+7u$, $\dot x_2=0.8333x_1-2.2381x_2-1.117u$, $y=x_2$; zero initial conditions (deviation variables).
Find. (a) $Y(s)/U(s)$; (b) $y(t)$ for a unit step in $u$.
Approach. Laplace-transform each state equation (zero ICs), solve the first (uncoupled) equation for $X_1(s)/U(s)$, substitute into the second to eliminate $X_1$, and combine over a common denominator; for (b) multiply by $1/s$ and expand in partial fractions using the two real poles.
(a) State 1. $(s+2.4048)X_1(s)=7U(s)\Rightarrow X_1(s)=\dfrac{7}{s+2.4048}U(s)$.
(a) State 2, eliminate $X_1$. $(s+2.2381)X_2(s)=0.8333X_1(s)-1.117U(s)=\left[\dfrac{0.8333\times7}{s+2.4048}-1.117\right]U(s)$. Combining over $(s+2.4048)$: $$\frac{0.8333\times7-1.117(s+2.4048)}{s+2.4048}=\frac{5.8331-1.117s-2.6858}{s+2.4048}=\frac{3.1473-1.117s}{s+2.4048}.$$
(a) Assemble the transfer function. $$\boxed{\frac{Y(s)}{U(s)}=\frac{3.1473-1.117s}{(s+2.4048)(s+2.2381)}=\frac{-1.117(s-2.818)}{(s+2.4048)(s+2.2381)}.}$$ Both poles are in the left half-plane ($s=-2.4048,-2.2381$, stable), but the numerator has a right-half-plane zero at $s=+2.818$ — by Problem 1's criterion this predicts an inverse response to a step in $u$.
(b) Partial fractions for the unit step. $Y(s)=\dfrac{3.1473-1.117s}{s(s+2.4048)(s+2.2381)}=\dfrac As+\dfrac{B}{s+2.4048}+\dfrac{C}{s+2.2381}$, with residues evaluated at each pole: $$A=\frac{3.1473}{2.4048\times2.2381}=0.5848,\quad B=\frac{3.1473-1.117(-2.4048)}{-2.4048(-2.4048+2.2381)}=14.552,$$ $$C=\frac{3.1473-1.117(-2.2381)}{-2.2381(-2.2381+2.4048)}=-15.137.$$ (Check: $A+B+C=0.5848+14.552-15.137\approx0=y(0)$ — consistent with zero initial conditions.)
(b) Result. $$\boxed{y(t)=0.5848+14.552\,e^{-2.4048t}-15.137\,e^{-2.2381t}.}$$ At $t=0.5$: $y=0.0137$ (still near zero — the two nearly-equal-time-constant exponentials almost cancel at first, the RHP-zero footprint); by $t=5$: $y=0.5846\approx y(\infty)=0.5848$ (settled to the DC gain).