NivaarExam PrepOfficial exam papers ↗

22-Elec-B2 Advanced Control Systems · December 2016

Question 3 of 5: Cascade-compensator design for a double-integrator plant

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

Notes on this paper

Paper format. 07-Elec-B2 Advanced Control Systems, December 2016 — open book, three hours, five questions. The cover page states that “any four questions constitute a complete paper” and that “all questions are of equal value (25%)”, so each question below carries 25 marks and a candidate answers four. Tables of Laplace and z-transforms are appended to the paper. All five questions are worked here, because this set is a study resource rather than a three-hour sitting. Question 1 is a nineteen-item short-answer block whose per-item mark weights are printed in the margin as [1] or [2] and total 25.

Reference texts. N. S. Nise, Control Systems Engineering, 8th ed. (Ch. 7 steady-state error and system type, Ch. 8 root-locus sketching rules, Ch. 9 root-locus design of cascade compensators, Ch. 10 frequency response and stability margins); R. C. Dorf and R. H. Bishop, Modern Control Systems, 13th ed. (Ch. 5, 7, 8); K. Ogata, Modern Control Engineering, 5th ed. (Ch. 5 transient response, Ch. 6 root-locus design, Ch. 7 frequency-response design); G. F. Franklin, J. D. Powell and A. Emami-Naeini, Feedback Control of Dynamic Systems, 8th ed. (Ch. 3, 5, 6). These are the references listed for exam code 07-Elec-B2 / 16-Elec-B2 in the Engineers Canada syllabus; the vocabulary of this particular paper (percent overshoot, settling time, “compensated system gain”, asymptote intercept) follows Nise closely.

Question 3: Cascade-compensator design for a double-integrator plant (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. A double-integrator plant inside unity feedback, with a cascade compensator whose zero is fixed by the question and whose pole and gain are the free variables.

Given data
SymbolValueMeaning
$G(s)$$K/s^{2}$plant: two poles at the origin, no zero
$T_{s}$1.667 sspecified 2% settling time
$\%OS$16.3%specified percent overshoot
$z_{c}$$-1$compensator zero, fixed by the question
$p_{c}$, $K$to be foundcompensator pole and compensated gain

Find. The dominant closed-loop pole pair, the compensator pole and gain that place the locus through it, all remaining closed-loop poles, an honest assessment of the second-order approximation, and the steady-state errors to step, ramp and parabolic commands.

Approach. Translate the two transient specifications into a single design point, use the angle condition to solve for the one unknown compensator pole that puts the locus through it, use the magnitude condition for the gain, then factor the resulting cubic to find the third pole and simulate the true response to test the approximation. The steady-state part is a system-type argument on the finished loop.

R(s)+−Gc(s)K / s²Y(s)unity feedback
Question 3: the cascade structure. The compensator $G_{c}(s)$ supplies one zero at $-1$ plus one pole and one gain; the plant $K/s^{2}$ contributes the double pole at the origin.
  1. Part (a) — convert the specifications into a design point. Inverting the overshoot relation, $$\zeta=\frac{-\ln(0.163)}{\sqrt{\pi^{2}+\ln^{2}(0.163)}}=0.5001 ,$$ a value the paper has clearly chosen so that $\zeta=0.5$ exactly. The settling time fixes the real part, $\sigma=4/T_{s}=4/1.667=2.400$, so $\omega_{n}=\sigma/\zeta=4.798$ rad/s and $\omega_{d}=\omega_{n}\sqrt{1-\zeta^{2}}=4.155$ rad/s: $$\boxed{s_{d}=-2.400\pm4.155j}$$
  2. Part (b) — the angle condition leaves the compensator pole as the only unknown. With one zero at $-1$, the double plant pole at the origin and an unknown pole at $-p_{c}$, $$\angle(s_{d}+1)-2\angle s_{d}-\angle(s_{d}+p_{c})=-180^\circ .$$ Substituting $\angle(s_{d}+1)=108.61^\circ$ and $2\angle s_{d}=2\times120.00^\circ=240.01^\circ$ gives the required pole angle $\angle(s_{d}+p_{c})=108.61^\circ-240.01^\circ+180^\circ=48.61^\circ$, and therefore $$p_{c}=\sigma+\frac{\omega_{d}}{\tan48.61^\circ} =2.400+\frac{4.155}{1.1346}=\boxed{6.062} ,$$ i.e. a compensator pole at $s=-6.062$. Re-evaluating the full angle sum at $s_{d}$ returns $-180.000^\circ$, so the locus does pass through the design point.
  3. Part (c) — the magnitude condition gives the gain. $$K=\frac{|s_{d}|^{2}\,|s_{d}+p_{c}|}{|s_{d}+1|} =\frac{(4.7986)^{2}\times5.5394}{4.3850}=\boxed{29.09}$$ so the compensated forward path is $G_{c}(s)G(s)=29.09(s+1)/[s^{2}(s+6.062)]$.
  4. Part (d) — factor the closed-loop cubic. The characteristic polynomial is $$s^{2}(s+6.062)+29.09(s+1)=s^{3}+6.062s^{2}+29.09s+29.09 ,$$ and dividing out the designed pair $s^{2}+4.799s+23.03$ leaves a single real root: $$\boxed{s_{3}=-1.263}$$ The other two roots come back as $-2.400\pm4.155j$ to six figures, confirming the design. There are no other closed-loop poles — the third-order loop has exactly three.
  5. Part (e) — the second-order approximation is poor, and the reason is geometric. Two tests both fail. The third pole is at $-1.263$, only $1.263/2.400=0.53$ times the real part of the dominant pair, where dominance conventionally needs a factor of at least five; and it is not cancelled by the zero at $-1$ either, since the pole-zero separation is $0.263$, a full 21% of the pole's own distance from the origin. Its step-response residue is $+0.327$, far from negligible. Simulating the true third-order step response against the ideal second-order model, $$\%OS_{\text{actual}}=34.7\%\ \text{vs}\ 16.3\%\ \text{designed}, \qquad T_{s,\text{actual}}=2.41\ \text{s vs}\ 1.667\ \text{s},$$ with the actual peak arriving early at $t=0.711$ s. Both specifications are missed, and missed on the pessimistic side, so the approximation is not merely inaccurate — it is optimistic, which is the dangerous direction.
  6. Part (f) — steady-state errors follow from the system type. The compensated loop $L(s)=29.09(s+1)/[s^{2}(s+6.062)]$ retains both poles at the origin, so it is Type 2 and $$K_{p}=\infty,\qquad K_{v}=\infty,\qquad K_{a}=\lim_{s\to0}s^{2}L(s)=\frac{29.09}{6.062}=4.798 .$$ Hence $$\boxed{e_{\text{step}}=0,\qquad e_{\text{ramp}}=0,\qquad e_{\text{parabola}}=\frac{1}{K_{a}}=0.208}$$ for the unit parabola $r=t^{2}/2$. If the parabola is written $r=t^{2}$ instead, the error doubles to 0.417 — worth stating explicitly, because the two conventions differ by a factor of two and the question does not say which it means.

The design is textbook-correct and still unsatisfactory, which is the real lesson of part (e). Forcing a plant with two poles at the origin through a specified point using only one zero and one pole leaves no freedom to push the third closed-loop pole out of the way: the same angle condition that places the dominant pair also determines where the leftover root lands, and here it lands close in. A designer who needed the specifications actually met would either move the compensator zero further left, accept a larger $p_{c}$, or add a second zero — and would check the simulated response rather than trusting the dominant-pole formulas.

Compensated locus: K(s+1) / [s²(s+6.062)]−8−7−6−5−4−3−2−112−6−5−4−3−2−1123456Re sIm sζ = 0.5−2.400 + 4.156jthird pole −1.263open-loop poleopen-loop zerodouble pole at the origincompensator zero −1,compensator pole−6.062locus traced for K from 0 upward
Parts (a) to (d): the compensated locus of $K(s+1)/[s^{2}(s+6.062)]$. The design point sits on the $\zeta=0.5$ ray; the third closed-loop pole at $-1.263$ is the root that spoils the second-order approximation.
Part (e): the second-order approximation against the true response0.000.240.490.730.971.211.460123456final value 1actual third-order responsesecond-order approximation34.7 % overshoot16.3 %time t (s)y(t)
Part (e): the true third-order step response against the ideal second-order model. The actual overshoot is 34.7% where 16.3% was specified, and the actual settling time is 2.41 s where 1.667 s was specified.
Final results
QuantityValue
(a) dominant poles$s=-2.400\pm4.155j$ ($\zeta=0.500$, $\omega_{n}=4.798$ rad/s)
(b) compensator pole$s=-6.062$ (required angle $48.61^\circ$)
(c) compensated gain$K=29.09$
(d) nondominant poleone real pole at $s=-1.263$; no others
(e) accuracy of the approximationpoor: actual $\%OS=34.7\%$, $T_{s}=2.41$ s, $T_{p}=0.711$ s
(f) $e_{ss}$, unit step0 (Type 2 loop)
(f) $e_{ss}$, unit ramp0 (Type 2 loop)
(f) $e_{ss}$, unit parabola $t^{2}/2$$1/K_{a}=0.208$ ($K_{a}=4.798$); $0.417$ if $r=t^{2}$