NivaarExam PrepOfficial exam papers ↗

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%)

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. A coupled two-state linear model with a direct feed-through of $u$ into the $x_2$ equation:

EquationExpression
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.

  1. 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.$$
  2. (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.
  3. (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}.$$
  4. 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.)
  5. 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.
0123456-0.1000.200.400.58time t (s)y(t)P3(b): unit-step response y(t) — inverse response then rise to 0.585
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.
QuantityValue
$Y/U$$(-1.117s+3.1469)/[(s+2.4048)(s+2.2381)]$
RHP zero$s=+2.817$ (inverse response)
Step response$y=0.585+14.55e^{-2.4048t}-15.14e^{-2.2381t}$
Undershoot / final value$\approx-0.094$ / $0.585$