NivaarExam PrepOfficial exam papers ↗

22-Elec-A2 Systems and Control · May 2017

Question 2 of 8: System type, error constants, dominant 2nd-order model (compulsory)

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 2 — System type, error constants, dominant 2nd-order model (compulsory) [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.

Part A — open loop, type and error constants

Given. the Q1 loop at $K=K_{op}=20$. Find. $G_{open}(s)$, the type, and $K_{pos},K_v,K_a$.

  1. Equivalent unity-feedback open loop. Since the actuating error is $e_1=R-Y$, the equivalent forward transfer function is $G_{open}=G_{cl}/(1-G_{cl})$: $$G_{open}(s)=\frac{10K(s+0.5)}{s\,(s^3+13s^2+30s+10)}.$$ One free integrator ($s$ in the denominator) $\Rightarrow$ Type 1.
  2. Error constants at $K_{op}=20$. $K_{pos}=\lim_{s\to0}G_{open}=\infty$ (Type 1); $K_v=\lim_{s\to0}sG_{open}=\dfrac{10K(0.5)}{10}=0.5K=\boxed{10}$; $K_a=\lim_{s\to0}s^2G_{open}=\boxed{0}$.
  3. Consistency. $e_{ss(ramp)}=1/K_v=1/10=0.1$ V/V and $e_{ss(step)}=1/(1+K_{pos})=0$ — exactly the Q1 specs.

Part B — dominant-pole model and step specs

Given. factored $G_{cl}$: poles $-11.9,\,-0.505,\,-0.296\pm j4.068$ and a zero at $-0.5$. Find. $G_{cl}(0)$, the reduced model, and the step specs. Approach. the zero at $-0.5$ nearly cancels the slow real pole at $-0.505$ and the pole at $-11.9$ is far away, so the complex pair dominates.

Q2 - Gcl poles: dominant pair -0.30+/-j4.07 (zero -0.5 ~cancels pole -0.505)-10-50-55ReIm
Q2 pole–zero map. The zero $-0.5$ all but cancels the real pole $-0.505$; the far pole $-11.9$ decays quickly, leaving the lightly-damped pair $-0.30\pm j4.07$ as the dominant modes.
  1. DC gain. $G_{cl}(0)=\dfrac{200(0.5)}{(11.9)(0.505)(16.64)}=\dfrac{100}{100.0}=\boxed{1}.$ (A unity DC gain is expected from the zero step-error result.)
  2. Cancel and drop. The near-cancelling factor $(s+0.5)/(s+0.505)\approx1$, and the far pole contributes its DC value $s+11.9\to11.9$. What remains is $\dfrac{200(s+0.5)}{(s+11.9)(s+0.505)}\Big|_{s=0}=16.64$ over the quadratic, so $K_{dc}\omega_n^2=16.64$.
  3. Model parameters. From $s^2+0.5925s+16.64$: $\omega_n=\sqrt{16.64}=\boxed{4.08\ \text{rad/s}}$, $\zeta=\dfrac{0.5925}{2\omega_n}=\boxed{0.0726}$, and $K_{dc}=\dfrac{16.64}{\omega_n^2}=\boxed{1.0}$. Hence $$\boxed{G_m(s)=\frac{16.64}{s^2+0.5925\,s+16.64}.}$$
  4. Overshoot. $PO=100\,e^{-\zeta\pi/\sqrt{1-\zeta^2}}=100\,e^{-0.2288}=\boxed{79.6\%}$ — a very lightly damped response.
  5. Rise / settle / period. With $\omega_d=\omega_n\sqrt{1-\zeta^2}=4.068$ rad/s and $\sigma=\zeta\omega_n=0.296$: the numerically evaluated 10–90% rise time is $T_{rise}=\boxed{0.265\ \text{s}}$; $T_{settle(\pm2\%)}=\dfrac{4}{\sigma}=\boxed{13.5\ \text{s}}$; $T_{period}=\dfrac{2\pi}{\omega_d}=\boxed{1.54\ \text{s}}$.
  6. Step error. $K_{dc}=1\Rightarrow e_{ss(step)}=\boxed{0\%}$.
Q2 - dominant 2nd-order step (PO=79.6%, Ts=13.5 s)0.00.20.40.60.81.01.21.41.61.82.001234567891011121314151617181920y_ss=1.0time (s)
Unit-step response of the dominant model $G_m$: first peak $\approx1.80$ (79.6% overshoot), ringing at $T_{period}=1.54$ s inside the $\pm2\%$ band reached at 13.5 s.
QuantityValue
System type / $K_v$Type 1 / $K_v=10$ ($K_{pos}=\infty,K_a=0$)
$G_{cl}(0)$$1$
$K_{dc},\ \zeta,\ \omega_n$$1.0,\ 0.0726,\ 4.08$ rad/s
$PO$$79.6\%$
$T_{rise(10-90)}$ / $T_{settle(2\%)}$$0.265$ s / $13.5$ s
$T_{period}$ / $e_{ss(step)}$$1.54$ s / $0\%$