23-Chem-A6 Process Dynamics and Control · May 2018
Question 6 of 8: Second-Order Transfer Function — Standard Form and Response 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, state-space models and dead-time systems; D. R. Coughanowr & S. E. LeBlanc, Process Systems Analysis and Control (3rd ed., McGraw-Hill) — first-order tank/thermal dynamics, step/ramp 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 Transfer Function — Standard Form and Response Regimes (20%)
Given. A linear second-order ODE with a variable damping coefficient $k$ and unit forcing coefficient on $x$: $\ddot y+k\dot y+2y=x$.
Find. (a) $Y(s)/X(s)$ in standard $K/(\tau^2s^2+2\zeta\tau s+1)$ form; (b) the qualitative response over $-20<k<20$ and its functional form.
Approach. Transform with zero initial conditions to read off the natural time constant and gain, express damping as a function of $k$, then classify the roots of the characteristic polynomial $s^2+ks+2$ by sign and discriminant.
(a) Transfer function. With zero initial conditions, $s^2Y+ksY+2Y=X$, so $$\frac{Y(s)}{X(s)}=\frac{1}{s^2+ks+2}.$$
(a) Standard form. Divide numerator and denominator by $2$ to make the constant term unity: $$\boxed{\frac{Y}{X}=\frac{K}{\tau^2s^2+2\zeta\tau s+1},\quad K=\tfrac12,\ \tau=\tfrac{1}{\sqrt2}=0.707,\ \zeta=\frac{k}{2\sqrt2}=0.354\,k.}$$ Matching $\tau^2=1/2$ gives $\tau=0.707$; matching $2\zeta\tau=k/2$ gives $\zeta=k/(2\sqrt2)$.
(b) Roots and stability. The poles solve $s^2+ks+2=0$: $s=\dfrac{-k\pm\sqrt{k^2-8}}{2}$. Their product is $+2$ (both signs tied) and their sum is $-k$; the real parts are negative only when $k>0$. So the response converges for $0<k<20$ and diverges for $-20<k<0$; $k=0$ gives a sustained (undamped) oscillation. The damping character switches at the discriminant $k^2-8=0$, i.e. $|k|=2\sqrt2=2.83$.
(b) Response forms (convergent side). With $\zeta=0.354k$: for $0<k<2.83$ ($\zeta<1$, underdamped) the poles are complex with negative real part, $$y(t)=y_\infty+e^{-\zeta t/\tau}\big(C_1\cos\omega_d t+C_2\sin\omega_d t\big),\quad \omega_d=\tfrac{\sqrt{8-k^2}}{2};$$ at $k=2.83$ ($\zeta=1$, critically damped) $y(t)=y_\infty+(C_1+C_2t)e^{-t/\tau}$; for $2.83<k<20$ ($\zeta>1$, overdamped) $y(t)=y_\infty+C_1e^{r_1t}+C_2e^{r_2t}$ with two negative real roots $r_{1,2}$.
(b) Response forms (divergent side). For $-2.83<k<0$ the complex poles have positive real part: $y(t)=y_\infty+e^{+|\sigma|t}(C_1\cos\omega_d t+C_2\sin\omega_d t)$ — a growing oscillation. For $-20<k<-2.83$ the two real roots are both positive: $y(t)=y_\infty+C_1e^{r_1t}+C_2e^{r_2t}$ with $r_{1,2}>0$ — monotone divergence. At $k=0$ the poles are $\pm j\sqrt2$: $y(t)=y_\infty+C_1\cos\sqrt2\,t+C_2\sin\sqrt2\,t$, an undamped oscillation of frequency $\sqrt2$.
Problem 6: the sign of $k$ sets stability (converges only for $k>0$), and $|k|=2\sqrt2$ divides oscillatory from monotone behaviour. Only the right half ($k>0$) yields a bounded, offset-free response.
Range of $k$
Poles
Response
$2.83<k<20$
real, $<0$
overdamped, converges
$0<k<2.83$
complex, $\mathrm{Re}<0$
underdamped, converges
$k=0$
$\pm j\sqrt2$
sustained oscillation
$-2.83<k<0$
complex, $\mathrm{Re}>0$
growing oscillation (diverges)
$-20<k<-2.83$
real, $>0$
monotone divergence
In standard form $K=\tfrac12$, $\tau=0.707$, $\zeta=0.354\,k$.