23-Chem-A6 Process Dynamics and Control · December 2018
Question 6 of 8: Second-Order ODE — Standard Form and Damping Regimes
Nivaar worked solution (AI-drafted; not reviewed by a licensed engineer)
Notes on this paper
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, dead-time systems and state-space models; D. R. Coughanowr & S. E. LeBlanc, Process Systems Analysis and Control (3rd ed., McGraw-Hill) — first-order tank/thermal dynamics, step/impulse response and block-diagram algebra. Standard control conventions (deviation variables; unity valve/sensor gains unless stated) are used throughout. This sitting is an open-book paper of eight problems of which any five constitute a complete exam; all eight are solved here.
Problem 6: Second-Order ODE — Standard Form and Damping Regimes (20%)
Given. A linear second-order ODE with a variable damping term:
Quantity
Value
ODE
$\ddot y+k\dot y+10y=2x$
Static gain factor
output/input $=2/10$
Natural coefficient
$10$ ($=\omega_n^2$)
Damping parameter
$k$ (free)
Find. (a) $Y/X$ in standard form; (b) $k$-ranges for stable / underdamped / overdamped; (c) $\tau(k)$ and $\zeta(k)$ when underdamped.
Problem 6: the character of the response is set by where $k$ falls. $\omega_n=\sqrt{10}$ is fixed; increasing $k$ moves the system from unstable, through underdamped, to overdamped, crossing critical damping at $k=2\sqrt{10}$.
Approach. Laplace-transform to get $Y/X$, normalise to $K/(\tau^2s^2+2\zeta\tau s+1)$ to read $K,\tau,\zeta$, then classify by the sign of the discriminant and the coefficient positivity.
(a) Transfer function. With zero initial conditions, $(s^2+ks+10)Y=2X$, so $Y/X=\dfrac{2}{s^2+ks+10}$. Dividing numerator and denominator by $10$: $$\boxed{\frac{Y}{X}=\frac{0.2}{\tfrac{1}{10}s^{2}+\tfrac{k}{10}s+1},\qquad K=0.2,\ \tau=\frac{1}{\sqrt{10}}=0.316,\ \zeta=\frac{k}{2\sqrt{10}}.}$$ Matching $\tau^2=1/10$ gives $\tau=1/\sqrt{10}$; matching $2\zeta\tau=k/10$ gives $\zeta=k/(2\sqrt{10})$.
(b-i) Stable. Roots $s=\dfrac{-k\pm\sqrt{k^2-40}}{2}$ have negative real part iff the coefficients are positive: $k>0$ (with $10>0$ automatic). So stable for $k>0$ ($k=0$ is marginally stable, $k<0$ unstable).
(b-ii) Underdamped. Complex roots require $k^2-40<0\Rightarrow|k|<2\sqrt{10}$; combined with stability, underdamped for $0 ($\zeta<1$).
(b-iii) Overdamped. Real distinct roots require $k^2-40>0$; taking the stable branch, overdamped for $k>2\sqrt{10}\approx6.32$ ($\zeta>1$). Critical damping ($\zeta=1$) occurs exactly at $k=2\sqrt{10}$.
(c) Underdamped expressions. The standard-form parameters are independent of the input, so directly from step 1: $$\boxed{\tau=\frac{1}{\sqrt{10}}=0.316\ (\text{constant}),\qquad \zeta=\frac{k}{2\sqrt{10}}=0.158\,k.}$$ Equivalently $\omega_n=1/\tau=\sqrt{10}$ is fixed and only the damping ratio $\zeta$ scales with $k$; for $0