23-Chem-A6 Process Dynamics and Control · May 2015
Question 3 of 8: State-Space Model → Transfer Function and Step Response
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 3: State-Space Model → Transfer Function and Step Response (20%)
Given. A coupled two-state linear model with a direct feed-through of $u$ into the $x_2$ equation:
Equation
Expression
State 1
$\dot x_1=-2.4048\,x_1+7u$
State 2
$\dot x_2=0.8333\,x_1-2.2381\,x_2-1.117u$
Output
$y=x_2$
Find. (a) $Y/U$; (b) $y(t)$ for a unit step in $u$.
Approach. Transform each state equation, eliminate $X_1$ to get $Y/U$, note the sign of the numerator zero, then invert the step-forced output by partial fractions.
Transform the states. With zero initial conditions, $X_1=\dfrac{7U}{s+2.4048}$ and $(s+2.2381)X_2=0.8333X_1-1.117U$. Substituting $X_1$: $$(s+2.2381)X_2=\Big[\tfrac{0.8333\times7}{s+2.4048}-1.117\Big]U.$$
(a) Transfer function. Clearing the denominator, $$\boxed{\frac{Y(s)}{U(s)}=\frac{-1.117\,s+3.1469}{(s+2.4048)(s+2.2381)}}.$$ (Numerator: $0.8333\times7-1.117\times2.4048=5.8331-2.6862=3.1469$.) The numerator has a right-half-plane zero at $s=+3.1469/1.117=+2.817$, which produces an inverse response.
(b) Step-forced output. For $U=1/s$, $$Y(s)=\frac{-1.117s+3.1469}{s\,(s+2.4048)(s+2.2381)}=\frac{A}{s}+\frac{B}{s+2.4048}+\frac{C}{s+2.2381}.$$
Residues. Evaluating by cover-up: $A=\dfrac{3.1469}{(2.4048)(2.2381)}=0.5847$; $B=\dfrac{-1.117(-2.4048)+3.1469}{(-2.4048)(-2.4048+2.2381)}=+14.55$; $C=\dfrac{-1.117(-2.2381)+3.1469}{(-2.2381)(-2.2381+2.4048)}=-15.14$. (Check: $A+B+C=0$, so $y(0)=0$ — consistent with a strictly proper transfer function.)
Response. Inverting, $$\boxed{y(t)=0.5847+14.55\,e^{-2.4048t}-15.14\,e^{-2.2381t}}.$$ Because the two exponentials nearly cancel at $t=0$ but the faster ($e^{-2.4048t}$) term with the large positive coefficient fades first, $y$ first dips negative (to $\approx-0.094$ near $t\approx0.3\,\mathrm s$) before climbing to its steady value $0.585$ — the classic inverse response of the RHP zero.
Problem 3(b): unit-step response. The right-half-plane zero at $s=+2.82$ makes the output first dip to $\approx-0.09$ before rising monotonically to its steady value $0.585$ — an inverse response.