NivaarExam PrepOfficial exam papers ↗

23-Chem-A6 Process Dynamics and Control · May 2014

Question 6 of 8: Second-Order ODE — Standard Form and Damping

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

Notes on this paper

National Exams / EGBC — May 2014 — 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. The parts are quantitative throughout — tank and reactor modelling, transfer functions, step/impulse responses, IMC design, Nyquist and Bode stability — and every requested plot (IMC response, Nyquist locus, thermocouple response, Bode diagram) is drawn as a real figure.

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, feedback stability, frequency response and IMC design; G. Stephanopoulos, Chemical Process Control: An Introduction to Theory and Practice (Prentice Hall) — Nyquist criterion, Bode stability and dead-time systems; D. R. Coughanowr & S. E. LeBlanc, Process Systems Analysis and Control (3rd ed., McGraw-Hill) — first- and second-order tank dynamics, step/impulse response and linearisation of non-linear balances. Standard control conventions (deviation variables; unity valve/sensor gains unless stated) are used throughout.

Problem 6: Second-Order ODE — Standard Form and Damping (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. $\ddot y+k\dot y+10y=2x$, zero initial conditions (deviation variables).

Find. (a) $Y/X$ in standard form; (b) the $k$-ranges for stability / under- / over-damping; (c) $\tau$ and $\zeta$ when underdamped.

Approach. Laplace-transform the ODE, then match the denominator to the canonical $\tau^2s^2+2\zeta\tau s+1$ to read off $K,\tau,\zeta$; the sign/discriminant of the characteristic roots then classifies the response.

  1. (a) Transfer function. Transforming (zero IC): $(s^2+ks+10)Y=2X$, so $Y/X=\dfrac{2}{s^2+ks+10}$. Dividing top and bottom by $10$ to force the constant term to unity: $$\boxed{\dfrac{Y}{X}=\dfrac{0.2}{\tfrac{1}{10}s^2+\tfrac{k}{10}s+1}=\dfrac{K}{\tau^2s^2+2\zeta\tau s+1}}$$ with gain $K=0.2$, $\tau^2=\tfrac1{10}$ and $2\zeta\tau=\tfrac{k}{10}$.
  2. Read off the parameters. $\omega_n=\sqrt{10}$, so $$\tau=\dfrac{1}{\omega_n}=\dfrac{1}{\sqrt{10}}\approx0.316,\qquad \zeta=\dfrac{k}{2\sqrt{10}}\approx\dfrac{k}{6.32}.$$ The characteristic equation is $s^2+ks+10=0$ with roots $s=\dfrac{-k\pm\sqrt{k^2-40}}{2}$.
  3. (b-i) Stability. For $s^2+ks+10$ both roots have negative real part iff every coefficient is positive; since $10>0$ always, the condition is $$\boxed{k>0}.$$ At $k=0$ the poles sit on the imaginary axis ($\pm j\sqrt{10}$, sustained oscillation); for $k<0$ the response grows unbounded.
  4. (b-ii,iii) Under- vs over-damped. The roots are complex ($\zeta<1$) when the discriminant is negative: $k^2-40<0\Rightarrow|k|<2\sqrt{10}$. With stability ($k>0$): $$\underbrace{0<k<2\sqrt{10}}_{\text{underdamped }(\approx6.32)}\ ;\qquad \underbrace{k>2\sqrt{10}}_{\text{overdamped}}\ ;\qquad k=2\sqrt{10}\ (\text{critically damped}).$$
  5. (c) Underdamped parameters. Directly from the standard form (valid for any $k$, and underdamped when $0<k<2\sqrt{10}$): $$\boxed{\tau=\dfrac{1}{\sqrt{10}}\approx0.316\ \text{(time units)}},\qquad \boxed{\zeta=\dfrac{k}{2\sqrt{10}}}.$$ The time constant is fixed by the $10y$ term (independent of $k$); only the damping coefficient scales with $k$.
0 1 2 3 4 0 0.1 0.2 0.3 k=2 underdamped k=6.32 critical k=12 overdamped time t y(t) P6(b): unit-step response vs k (all ->0.2)
Problem 6(b): unit-step response for representative $k$. $k=2$ is underdamped (overshoot and ringing), $k=2\sqrt{10}\approx6.32$ is critically damped, $k=12$ is overdamped (sluggish, no overshoot); all settle at the gain $K\cdot\Delta x=0.2$.
QuantityExpression / value
Transfer function$Y/X=0.2/(\tfrac1{10}s^2+\tfrac{k}{10}s+1)$
Gain $K$$0.2$
Natural frequency / time constant$\omega_n=\sqrt{10}$, $\tau=1/\sqrt{10}\approx0.316$
Damping ratio$\zeta=k/(2\sqrt{10})\approx k/6.32$
Stable / underdamped / overdamped$k>0$ / $0<k<2\sqrt{10}$ / $k>2\sqrt{10}$