NivaarExam PrepOfficial exam papers ↗

22-Mec-A3 System Analysis and Control · May 2014

Question 4 of 6: Transient-response terms under proportional control

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

Notes on this paper

Paper format. National Exams, May 2014 — 07-Mec-A3 System Analysis and Control. Closed book, 3 hours; a Casio or Sharp approved calculator and semi-log graph paper are the only aids. Six questions are printed; any four constitute a complete paper and all questions are of equal value (25 marks each). All six questions are solved in full. A Laplace-transform table is appended to the exam (pages 5–6) and two sheets of blank four-cycle semi-log paper are provided for the Bode work.

Reference texts. K. Ogata, Modern Control Engineering, 5th ed. (Pearson) — Chapters 5–7; N. S. Nise, Control Systems Engineering, 8th ed. (Wiley) — Chapters 6–10; R. C. Dorf and R. H. Bishop, Modern Control Systems, 13th ed. (Pearson); G. F. Franklin, J. D. Powell and A. Emami-Naeini, Feedback Control of Dynamic Systems, 8th ed. (Pearson). These are the standard Canadian undergraduate control texts and cover every technique used below (Routh–Hurwitz, root locus, Bode margins, static error constants, frequency response).

Notation. Throughout, K denotes the root-locus (pole–zero form) gain that multiplies the open-loop transfer function as printed, and ω is in rad/s. Magnitudes in decibels are 20 log10 of the magnitude ratio.

Question 4: Transient-response terms under proportional control (25 marks)

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. The Figure-1 loop: reference $R(s)$, error into a proportional controller $G_c(s)=K_p$, plant $200/[(s+1)(s+2)]$, and a feedback element of gain 0.02 (a non-unity, sensor-scaling feedback).

ElementTransfer function
Controller$G_c(s)=K_p$
Plant$G_p(s)=\dfrac{200}{(s+1)(s+2)}$
Feedback (sensor)$H(s)=0.02$
Cases required$K_p = 0.05$; $K_p = 0.5$; then the minimum $K_p$ for fastest decay

Find. The transient (natural-response) terms — i.e. the exponential and damped-sinusoidal modes fixed by the closed-loop poles — for each of the two gains, and the smallest proportional gain that makes the transient decay as fast as the loop allows.

[Figure not reproduced: Figure 1 redrawn: proportional controller, plant 200/[(s+1)(s+2)] and a 0.02 sensor gain in the feedback path. See the official exam paper.]

Approach. Reduce the loop to a single second-order transfer function, so that its characteristic polynomial gives the closed-loop poles as a function of $K_p$; the transient terms are then read directly from the pole pattern, and the "fastest decay" question becomes a question about where the two real poles coalesce.

  1. Reduce the block diagram. The forward path is $G_cG_p$ and the feedback path is $H = 0.02$, so $$\frac{C(s)}{R(s)}=\frac{G_cG_p}{1+G_cG_pH}=\frac{200K_p}{(s+1)(s+2)+200K_p(0.02)} =\frac{200K_p}{s^{2}+3s+2+4K_p}$$ Note the useful simplification: the loop gain is $200\times0.02 = 4$ times $K_p$, so the characteristic equation depends on $K_p$ only through the product $4K_p$.
  2. Write the characteristic equation. The transient modes are governed by $$\boxed{\;s^{2}+3s+(2+4K_p)=0\;}\qquad\Longrightarrow\qquad s_{1,2}=\frac{-3\pm\sqrt{9-4(2+4K_p)}}{2}=\frac{-3\pm\sqrt{1-16K_p}}{2}$$ The damping term $-3/2$ is fixed by the plant; only the discriminant moves with gain.
  3. Part (a): $K_p = 0.05$. The constant term is $2+4(0.05)=2.2$ and the discriminant is $9-8.8=0.2>0$, so the poles are real and distinct: $$s_{1,2}=\frac{-3\pm 0.4472}{2}=-1.2764,\ -1.7236$$ The transient response therefore consists of two decaying exponentials, $$\boxed{\;c_{tr}(t)=A_1e^{-1.2764t}+A_2e^{-1.7236t}\;}$$ with time constants $\tau_1 = 1/1.2764 = 0.783$ s and $\tau_2 = 1/1.7236 = 0.580$ s. The response is over-damped and settles in roughly $4\tau_1 \approx 3.1$ s, dominated by the slower mode. For reference, the closed-loop dc gain is $200(0.05)/2.2 = 4.545$.
  4. Part (b): $K_p = 0.5$. Now the constant term is $2+4(0.5)=4$ and the discriminant is $9-16=-7<0$, so the poles form a complex-conjugate pair: $$s_{1,2}=\frac{-3\pm j\sqrt{7}}{2}=-1.5\pm j1.3229$$ and the transient terms are a decaying sine and cosine at the damped natural frequency, $$\boxed{\;c_{tr}(t)=e^{-1.5t}\left(B_1\cos 1.3229t+B_2\sin 1.3229t\right)\;}$$ Comparing with the standard form $s^{2}+2\zeta\omega_n s+\omega_n^{2}$ gives $\omega_n=\sqrt{4}=2$ rad/s and $\zeta=1.5/2=0.75$, i.e. a well-damped oscillatory response with about 2.8 % overshoot and a settling time near $4/1.5 = 2.7$ s. The dc gain has risen to $200(0.5)/4 = 25$.
  5. Part (c): the gain that gives the fastest decay. The decay rate of the transient is set by the least negative real part among the closed-loop poles. While the poles are real (small $K_p$) the slower pole is $s_1=(-3+\sqrt{1-16K_p})/2$, which becomes more negative as $K_p$ rises. Once the discriminant turns negative both poles sit at $\sigma=-1.5$ and the decay rate stops improving — further gain only adds oscillation. The transition is at the repeated root: $$9-4(2+4K_p)=0\quad\Longrightarrow\quad 2+4K_p=2.25\quad\Longrightarrow\quad \boxed{\;K_{p,\min}=0.0625\;}$$ at which the characteristic equation is $(s+1.5)^{2}=0$, giving a critically damped response with transient terms $c_{tr}(t)=(C_1+C_2t)e^{-1.5t}$.
  6. Confirm the claim numerically. Sweeping the gain: at $K_p=0.03$ the slower pole is $-1.139$; at $K_p=0.05$ it is $-1.276$; at $K_p=0.0625$ it is exactly $-1.5$; and for every larger gain (0.1, 1, 10) the real part stays pinned at $-1.5$. Hence 0.0625 is the minimum gain achieving the fastest available decay envelope $e^{-1.5t}$ — which is precisely why the question asks for a minimum rather than an optimum.
Normalised closed-loop step response of the Figure-1 loop0.01.02.03.04.05.06.00.00.20.40.60.81.01.2Kp = 0.05 (over-damped, two real modes)Kp = 0.5 (under-damped, ζ = 0.75)time t (s)c(t)/c(∞)
Normalised unit-step responses for the two gains, showing the over-damped (Kp = 0.05) and under-damped (Kp = 0.5) characters.
Check: part (c) is printed with the symbol $K$ while parts (a) and (b) use $K_p$; since the only adjustable gain in Figure 1 is the proportional controller gain, $K$ is read here as $K_p$. Under this reading the three parts form a coherent progression (over-damped → under-damped → the critically damped boundary between them), and the answer 0.0625 lies between the two gains examined in (a) and (b).
CaseClosed-loop polesTransient termsCharacter
(a) Kp = 0.05−1.2764, −1.7236 $A_1e^{-1.2764t}+A_2e^{-1.7236t}$over-damped, τ = 0.783 s and 0.580 s
(b) Kp = 0.5−1.5 ± j1.3229 $e^{-1.5t}(B_1\cos 1.3229t+B_2\sin 1.3229t)$ under-damped, ωn = 2 rad/s, ζ = 0.75
(c) Kp,min = 0.0625−1.5 (double) $(C_1+C_2t)e^{-1.5t}$critically damped — fastest decay envelope