NivaarExam PrepOfficial exam papers ↗

22-Elec-A2 Systems and Control · May 2017

Question 8 of 8: Series-PID design by pole placement

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

Notes on this paper

Paper format. National Exams, 16-Elec-A2 Systems & Control, May 2017, 3 hours, closed book (approved Casio/Sharp calculator plus one signed, double-sided 8.5 × 11" formula sheet; a Laplace-transform table and the standard $\zeta$–overshoot and $\zeta$–resonant-peak design plots are supplied on pages 2–3). Questions 1 and 2 are compulsory; a complete paper is five questions, so a candidate chooses three of Q3–Q8. Each question is worth 20 marks. All eight are worked below so the set is a complete study resource.

Reference texts: N. S. Nise, Control Systems Engineering (7th ed., Wiley) — Routh–Hurwitz (Ch. 6), root locus (Ch. 8), steady-state error and static error constants (Ch. 7), frequency response, Nyquist, gain/phase margins (Ch. 10), PID and compensator design (Ch. 9–11), state space, controllability/observability and Mason’s rule (Ch. 3, 5, 12); K. Ogata, Modern Control Engineering (5th ed., Prentice Hall) — dominant-poles modelling, second-order correlations and pole-placement. All block diagrams, root loci, pole–zero maps, Nyquist plots, step responses and the state-realization diagrams below are redrawn as inline figures.

Reading the supplied design charts. The percent-overshoot chart uses $PO=100\,e^{-\zeta\pi/\sqrt{1-\zeta^2}}$ and the second-order model is $G_m(s)=K_{dc}\,\dfrac{\omega_n^2}{s^2+2\zeta\omega_n s+\omega_n^2}$. Every result here is derived analytically in the $s$-domain and cross-checked; the design charts are used only to corroborate $\zeta$ from $PO$.

Question 8 — Series-PID design by pole placement [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.

−+R(s)Kp(1+Ki/s)(Kd s+1)series PIDP(s)1/(s^2+7s+8)Y(s)Figure Q8.1 - series-PID positioning loop
Figure Q8.1 — series-PID loop, controller $K_p(1+K_i/s)(K_d s+1)$ driving the process $1/(s^2+7s+8)$, unity feedback.

Item (1) — closed-loop TF and characteristic equation

  1. Open loop. $G_c(s)=\dfrac{K_p(s+K_i)(K_d s+1)}{s}$, so $L(s)=\dfrac{K_p(s+K_i)(K_d s+1)}{s(s^2+7s+8)}$ and $G_{cl}=\dfrac{L}{1+L}$ with numerator $K_p(s+K_i)(K_d s+1)$.
  2. Characteristic equation. $Q(s)=s(s^2+7s+8)+K_p(s+K_i)(K_d s+1)=0$. Expanding $(s+K_i)(K_d s+1)=K_d s^2+(1+K_iK_d)s+K_i$: $$\boxed{Q(s)=s^3+(7+K_pK_d)s^2+\big(8+K_p(1+K_iK_d)\big)s+K_pK_i=0.}$$

Item (2) — transient targets

  1. Damping. $\zeta=\dfrac{-\ln(0.15)}{\sqrt{\pi^2+\ln^2 0.15}}=\boxed{0.517}$.
  2. Natural frequency. $T_{settle(\pm2\%)}=\dfrac{4}{\zeta\omega_n}=2$ s $\Rightarrow \zeta\omega_n=2$, so $\omega_n=\dfrac{2}{0.517}=\boxed{3.87\ \text{rad/s}}$. The dominant pair is $-2\pm j\omega_d$, $\omega_d=\omega_n\sqrt{1-\zeta^2}=3.31$ rad/s.

Item (3) — pole placement and gains

Approach. the third pole is set to $-K_i$ so the closed-loop zero at $s=-K_i$ cancels it. Match $Q(s)$ to $(s^2+4s+\omega_n^2)(s+K_i)$ with $\omega_n^2=14.97$.

  1. Desired polynomial. $(s^2+4s+14.97)(s+K_i)=s^3+(K_i+4)s^2+(4K_i+14.97)s+14.97K_i$.
  2. Constant term. $K_pK_i=14.97K_i\Rightarrow \boxed{K_p=14.97}=\omega_n^2$.
  3. $s^2$ term. $7+K_pK_d=K_i+4\Rightarrow K_pK_d=K_i-3$.
  4. $s^1$ term. $8+K_p+K_pK_iK_d=4K_i+14.97$; substituting $K_pK_d=K_i-3$ gives $14.97+K_i(K_i-3)=4K_i+6.97$, i.e. $$K_i^2-7K_i+8=0\Rightarrow K_i=\frac{7\pm\sqrt{17}}{2}=5.56\ \text{or}\ 1.44.$$
  5. Derivative gains. $K_d=\dfrac{K_i-3}{14.97}$: for $K_i=5.56$, $K_d=+0.171$; for $K_i=1.44$, $K_d=-0.104$.

Item (4) — choice and commentary

  1. Select the physical set. The branch $K_i=1.44$ gives a negative derivative gain $K_d=-0.104$ (a non-minimum-phase RHP zero at $+9.6$ and a controller that is not a proper lead) — reject it. Choose $$\boxed{K_p=14.97,\quad K_i=5.56,\quad K_d=0.171.}$$
  2. Model vs actual. The third pole at $-K_i=-5.56$ is cancelled by the controller zero at $-K_i$, so it does not appear in the response; the remaining closed-loop zero at $-1/K_d=-5.85$ is far from the dominant pair, so its effect is small. The actual response therefore tracks the second-order model closely — slightly faster rise and a touch more overshoot than the ideal because a finite left-half-plane zero always adds a little lead. (The rejected set would have placed an uncancelled RHP zero, badly degrading the transient.)
Q8 - CL poles: dominant pair (zeta=0.52,wn=3.87) + pole -5.56 cancelled by zero -5.56-50-55ReIm
Q8 closed-loop poles at the chosen gains: dominant pair $-2\pm j3.31$ ($\zeta=0.52,\omega_n=3.87$) and a third pole at $-5.56$ cancelled by the controller zero; a far zero at $-5.85$ remains.
QuantityValue
$Q(s)$$s^3+(7+K_pK_d)s^2+(8+K_p(1+K_iK_d))s+K_pK_i$
$\zeta,\ \omega_n$$0.517,\ 3.87$ rad/s
Dominant pair$-2\pm j3.31$
$K_i$ roots$5.56$ or $1.44$
Chosen gains$K_p=14.97,\ K_i=5.56,\ K_d=0.171$
Rejected set$K_i=1.44,\ K_d=-0.104$ (negative $K_d$, RHP zero)
Back to the paper →