Skip to main content

Chapter 7 · 10 hours

Performance Specifications and Compensation Design

IOE past exam questions

Past questions and answers

49 questions set from this chapter, 10 of them more than once. Most asked first.

  • Asked 4 times
  • 2076 Asoj · 2+14 marks
  • 2071 Shrawan · 16 marks
  • 2078 Bhadra · 12 marks
  • 2076 Chaitra · 16 marks

The open loop transfer function of type-II system with unity feedback is given by G(s) = K/[s²(1+0.25s)]. Design a lead compensator to meet the following specifications: (i) Acceleration error constant (Ka) = 10/sec² (ii) PM at least 35°.

Answer

A lead compensator Gc(s)=1+Ts1+αTsG_c(s)=\dfrac{1+Ts}{1+\alpha Ts} (α<1\alpha<1) adds positive phase near the new gain crossover frequency, which raises the phase margin. Its maximum phase lead ϕm\phi_m occurs at ωm=1Tα\omega_m=\dfrac{1}{T\sqrt\alpha}, with sin⁡ϕm=1−α1+α\sin\phi_m=\dfrac{1-\alpha}{1+\alpha}, and there it adds a gain of 10log⁡(1/α)10\log(1/\alpha) dB. (The design below uses the Bode plot method.)

Step 1: Gain K from the error constant

Ka=lim⁡s→0s2G(s)=lim⁡s→0s2Ks2(1+0.25s)=K=10K_a=\lim_{s\to0}s^2G(s)=\lim_{s\to0}\frac{s^2K}{s^2(1+0.25s)}=K=10 G(s)=10s2(1+0.25s)G(s)=\frac{10}{s^2(1+0.25s)}

Step 2: Uncompensated system

∣G(jω)∣=10ω21+0.0625ω2,ϕ=−180∘−tan⁡−10.25ω|G(j\omega)|=\frac{10}{\omega^2\sqrt{1+0.0625\omega^2}},\qquad \phi=-180^\circ-\tan^{-1}0.25\omega

Asymptotic Bode magnitude: −40-40 dB/dec through 20log⁡10=2020\log10=20 dB at ω=1\omega=1; corner at ω=4\omega=4, then −60-60 dB/dec.

ω\omega0.512.32.8546.3910
$G$ (dB)32.019.74.30−7.1-7.1
ϕ\phi−187∘-187^\circ−194∘-194^\circ−210∘-210^\circ−215.5∘-215.5^\circ−225∘-225^\circ−238∘-238^\circ−248∘-248^\circ

Solving ∣G∣=1|G|=1: ωgc=2.85\omega_{gc}=2.85 rad/s and

PM=180∘+(−180∘−tan⁡−10.713)=−35.5∘PM=180^\circ+(-180^\circ-\tan^{-1}0.713)=-35.5^\circ

The uncompensated system is unstable (phase is always below −180∘-180^\circ).

Step 3: Phase lead required

ϕreq=PMd−PMunc+ϵ=35∘−(−35.5∘)+ϵ≈75∘–85∘\phi_{req}=PM_d-PM_{unc}+\epsilon = 35^\circ-(-35.5^\circ)+\epsilon\approx 75^\circ\text{–}85^\circ

A single lead stage cannot do this: as the gain crossover moves right, the plant phase falls further (the (1+0.25s)(1+0.25s) lag), and a check shows one lead section can give at most about 18.6∘18.6^\circ of PM for this plant. So we use two identical lead sections in cascade, each giving about half the lead (a double-lead compensator):

Gc(s)=(1+Ts1+αTs)2G_c(s)=\left(\frac{1+Ts}{1+\alpha Ts}\right)^2

Step 4: Choose α

Trial 1: ϕm=45∘\phi_m=45^\circ per stage (α=0.172\alpha=0.172) gives PM ≈ 34.8° (just short). Trial 2: ϕm≈50∘\phi_m\approx50^\circ per stage. Take α=0.13\alpha=0.13:

sin⁡ϕm=1−0.131+0.13=0.770⇒ϕm=50.3∘ per stage (100.7∘ total)\sin\phi_m=\frac{1-0.13}{1+0.13}=0.770\Rightarrow\phi_m=50.3^\circ\ \text{per stage}\ (100.7^\circ\ \text{total})

Step 5: New gain crossover frequency

Two stages add 2×10log⁡(1/α)=20log⁡(1/0.13)=17.722\times10\log(1/\alpha)=20\log(1/0.13)=17.72 dB at ωm\omega_m. So the new crossover is where the uncompensated magnitude is −17.72-17.72 dB:

10ω21+0.0625ω2=0.13⇒ωm=ωgc′=6.39 rad/s\frac{10}{\omega^2\sqrt{1+0.0625\omega^2}}=0.13\Rightarrow\omega_m=\omega_{gc}'=6.39\ \text{rad/s}

Step 6: Time constants

T=1ωmα=16.39×0.13=0.434 sαT=0.13×0.434=0.0564 s\begin{aligned} T&=\frac{1}{\omega_m\sqrt\alpha}=\frac{1}{6.39\times\sqrt{0.13}}=0.434\ \text{s}\\ \alpha T&=0.13\times0.434=0.0564\ \text{s} \end{aligned}

Corner frequencies: zero at 1/T=2.301/T=2.30 rad/s, pole at 1/αT=17.71/\alpha T=17.7 rad/s.

Gc(s)=(1+0.434s1+0.0564s)2G_c(s)=\left(\frac{1+0.434s}{1+0.0564s}\right)^2

The attenuation of the lead network α2\alpha^2 is made up by an amplifier of gain 1/α2=59.21/\alpha^2=59.2, so that the DC gain and hence Ka=10K_a=10 are unchanged.

Step 7: Compensated system and check

Gc(s)G(s)=10(1+0.434s)2s2(1+0.25s)(1+0.0564s)2G_c(s)G(s)=\frac{10(1+0.434s)^2}{s^2(1+0.25s)(1+0.0564s)^2}

At ω=6.39\omega=6.39 rad/s:

ϕ=−180∘−tan⁡−1(1.60)+2×50.3∘=−180∘−58.0∘+100.7∘=−137.3∘PM=180∘−137.3∘=42.7∘ (≥35∘)\begin{aligned} \phi&=-180^\circ-\tan^{-1}(1.60)+2\times50.3^\circ\\ &=-180^\circ-58.0^\circ+100.7^\circ=-137.3^\circ\\ PM&=180^\circ-137.3^\circ=42.7^\circ\ (\ge35^\circ) \end{aligned}

GM of the compensated system ≈ 12.9 dB (phase crossover at 17.1 rad/s). All closed-loop poles are in the LHP (−27.0-27.0, −4.08±j7.90-4.08\pm j7.90, −2.14±j1.14-2.14\pm j1.14).

R(s)->( )->[Amp 59.2]->[(s+2.3)/(s+17.7)]^2 ->[G(s)]-+->C
       ^-                                            |
       |_____________________________________________|
ItemUncompensatedCompensated
KaK_a1010
ωgc\omega_{gc}2.85 rad/s6.39 rad/s
PM−35.5∘-35.5^\circ42.7∘42.7^\circ
Stabilityunstablestable

Answer: Gc(s)=(1+0.434s1+0.0564s)2G_c(s)=\left(\dfrac{1+0.434s}{1+0.0564s}\right)^2 (two lead stages, α = 0.13 each) with K = 10 gives Ka = 10 s⁻² and PM ≈ 43° (> 35°).

  • Asked 3 times
  • 2081 Bhadra · 4 marks
  • 2074 Asoj · 1+3 marks
  • 2068 Baisakh (old course) · 1+3 marks

State whether the statement "Derivative controllers are always used with other controllers" is true or false and justify your answer.

Answer

True. A derivative (D) controller is never used alone; it is always combined with a proportional and/or integral controller, as PD or PID.

Justification

A derivative controller gives an output proportional to the rate of change of the error:

u(t)=Kdde(t)dt,Gc(s)=Kdsu(t)=K_d\frac{de(t)}{dt},\qquad G_c(s)=K_ds
  1. No action on constant error: if the error is constant (even large), de/dt=0de/dt=0, so the controller output is zero. A steady error is never corrected, so D alone cannot hold the output at the set point.
  2. Blocks DC signals: Gc(0)=0G_c(0)=0. In the loop it acts like a zero at the origin, which can cancel the plant's integrator and reduce the system type, so the steady-state error grows.
  3. Noise amplification: ∣Gc(jω)∣=Kdω|G_c(j\omega)|=K_d\omega rises with frequency, so high-frequency noise is amplified strongly. With a P term the useful low-frequency action dominates.
  4. Anticipatory, not corrective: D action predicts the error trend and adds damping; it only improves the transient part of the response.

How it is used

With proportional action (PD):

Gc(s)=Kp+Kds=Kp(1+Tds)G_c(s)=K_p+K_ds=K_p(1+T_ds)
  • The P part removes and holds the error.
  • The D part adds a zero, increases damping, reduces overshoot and settling time, and improves stability.
ControllerSteady error correctionDamping / overshoot
D alonenoneimproves
PDyes (P)improves
PIDeliminated (I)improves

Example: in a position servo, a PD controller lets the motor slow down as it nears the target (D) while still pushing it there (P). A pure D controller would stop acting once the motor is at rest away from the target.

  • Asked 3 times
  • 2074 Asoj · 12 marks
  • 2068 Chaitra · 16 marks
  • 2081 Baisakh · 12 marks

Design a suitable lead compensator for a system having open loop TF G(s) = K/[s(1+0.1s)(1+0.001s)] such that the compensated system should have phase margin of at least 45° and static velocity error constant of at least 1000.

Answer

A lead compensator Gc(s)=1+Ts1+αTsG_c(s)=\dfrac{1+Ts}{1+\alpha Ts} (α<1\alpha<1) is designed on the Bode plot: it supplies its maximum lead ϕm\phi_m (sin⁡ϕm=1−α1+α\sin\phi_m=\frac{1-\alpha}{1+\alpha}) at the new gain crossover ωm=1Tα\omega_m=\frac{1}{T\sqrt\alpha}, where it adds 10log⁡(1/α)10\log(1/\alpha) dB of gain.

Step 1: K from the velocity error constant

Kv=lim⁡s→0sG(s)=K=1000 s−1K_v=\lim_{s\to0}sG(s)=K=1000\ \text{s}^{-1} G(s)=1000s(1+0.1s)(1+0.001s)G(s)=\frac{1000}{s(1+0.1s)(1+0.001s)}

Step 2: Uncompensated Bode plot

Corners at 10 and 1000 rad/s. Asymptotic magnitude: 60 dB at ω=1\omega=1, −20-20 dB/dec to 10 (40 dB), −40-40 dB/dec to 1000 (−40-40 dB), then −60-60 dB/dec.

ϕ(ω)=−90∘−tan⁡−10.1ω−tan⁡−10.001ω\phi(\omega)=-90^\circ-\tan^{-1}0.1\omega-\tan^{-1}0.001\omega
ω\omega11055.61001765561000
$G$ (dB)60.037.010.1−0.1-0.1−10.0-10.0
ϕ\phi−95.8∘-95.8^\circ−135.6∘-135.6^\circ−173∘-173^\circ−180∘-180^\circ−186.7∘-186.7^\circ−208∘-208^\circ−224.4∘-224.4^\circ

From ∣G∣=1|G|=1: ωgc=99.5\omega_{gc}=99.5 rad/s, ϕ=−179.94∘\phi=-179.94^\circ, so

PMunc≈0∘PM_{unc}\approx0^\circ

(and ωpc=100\omega_{pc}=100 rad/s, GM ≈ 0 dB). The system is on the verge of instability.

Step 3: Required phase lead

ϕm=PMd−PMunc+ϵ=45∘−0∘+10∘=55∘\phi_m=PM_d-PM_{unc}+\epsilon=45^\circ-0^\circ+10^\circ=55^\circ

A margin of ϵ=10∘\epsilon=10^\circ is taken because the crossover will move to a higher frequency where the plant phase is lower (a 5° margin gives only 44° and fails the check).

Step 4: α

α=1−sin⁡55∘1+sin⁡55∘=1−0.8191+0.819=0.0994≈0.1\alpha=\frac{1-\sin55^\circ}{1+\sin55^\circ}=\frac{1-0.819}{1+0.819}=0.0994\approx0.1

With α=0.1\alpha=0.1: ϕm=sin⁡−1(0.9/1.1)=54.9∘\phi_m=\sin^{-1}(0.9/1.1)=54.9^\circ.

Step 5: New gain crossover frequency

The lead adds 10log⁡(1/0.1)=1010\log(1/0.1)=10 dB at ωm\omega_m, so place ωm\omega_m where ∣G(jω)∣=−10|G(j\omega)|=-10 dB:

1000ω1+0.01ω21+10−6ω2=0.316⇒ωm≈176 rad/s\frac{1000}{\omega\sqrt{1+0.01\omega^2}\sqrt{1+10^{-6}\omega^2}}=0.316\Rightarrow\omega_m\approx176\ \text{rad/s}

Step 6: Time constants

T=1ωmα=1176×0.316=0.018 sαT=0.0018 s\begin{aligned} T&=\frac{1}{\omega_m\sqrt\alpha}=\frac{1}{176\times0.316}=0.018\ \text{s}\\ \alpha T&=0.0018\ \text{s} \end{aligned}

Zero at 1/T=55.61/T=55.6 rad/s, pole at 1/αT=5561/\alpha T=556 rad/s.

Gc(s)=1+0.018s1+0.0018sG_c(s)=\frac{1+0.018s}{1+0.0018s}

An amplifier of gain 1/α=101/\alpha=10 makes up the lead network's attenuation, so KvK_v stays 1000.

Step 7: Compensated system and check

Gc(s)G(s)=1000(1+0.018s)s(1+0.1s)(1+0.001s)(1+0.0018s)G_c(s)G(s)=\frac{1000(1+0.018s)}{s(1+0.1s)(1+0.001s)(1+0.0018s)}
ω\omega1055.6100176.8556696
$G_cG$ (dB)37.113.06.10
ϕ\phi−126∘-126^\circ−134∘-134^\circ−129∘-129^\circ−131.9∘-131.9^\circ−169∘-169^\circ−180∘-180^\circ

At ωgc′=176.8\omega_{gc}'=176.8 rad/s:

ϕ=−90∘−tan⁡−117.68−tan⁡−10.177+tan⁡−13.18−tan⁡−10.318=−90∘−86.8∘−10.0∘+72.5∘−17.6∘=−131.9∘PM=180∘−131.9∘=48.1∘ (≥45∘)\begin{aligned} \phi&=-90^\circ-\tan^{-1}17.68-\tan^{-1}0.177+\tan^{-1}3.18-\tan^{-1}0.318\\ &=-90^\circ-86.8^\circ-10.0^\circ+72.5^\circ-17.6^\circ=-131.9^\circ\\ PM&=180^\circ-131.9^\circ=48.1^\circ\ (\ge45^\circ) \end{aligned}

GM ≈ 17.5 dB at ωpc≈696\omega_{pc}\approx696 rad/s.

 dB 60 |\  uncompensated  -20
    37 |  \__ 10    compensated
     0 |------\----x----------- w
       |    100 \  176.8
       |         \__  -40/-60
 lead adds phase between 55.6 and 556 rad/s
ItemBeforeAfter
KvK_v10001000
ωgc\omega_{gc}99.5 rad/s176.8 rad/s
PM≈ 0°48.1°
GM≈ 0 dB17.5 dB

Answer: Gc(s)=1+0.018s1+0.0018sG_c(s)=\dfrac{1+0.018s}{1+0.0018s} (α = 0.1, T = 0.018 s) with K = 1000 gives Kv = 1000 s⁻¹ and PM ≈ 48° (≥ 45°).

  • Asked 2 times
  • 2080 Baisakh · 4 marks
  • 2079 Bhadra · 6 marks

What is derivative feedback controller? Draw block diagram and find the transfer function and hence show its effect on transient performance of the system.

Answer

Derivative (rate or tachometer) feedback control feeds back a signal proportional to the rate of change of the output, Kt dc/dtK_t\,dc/dt, through an inner (minor) loop, in addition to the usual unity feedback. In a position servo this signal comes from a tachogenerator on the motor shaft.

Block diagram

R(s)    E(s)        +----------------+     C(s)
--->(+)---->(+)---->| wn^2/(s(s+2z wn))|--+---->
     ^-      ^-     +----------------+  |
     |       |                          |
     |       +-------[ Kt s ]<----------+
     |                                  |
     +----------------------------------+

Transfer function

Inner loop (forward G(s)=ωn2s(s+2ζωn)G(s)=\dfrac{\omega_n^2}{s(s+2\zeta\omega_n)}, feedback KtsK_ts):

G1(s)=G1+GKts=ωn2s(s+2ζωn)+Ktωn2s=ωn2s(s+2ζωn+Ktωn2)\begin{aligned} G_1(s)&=\frac{G}{1+GK_ts}=\frac{\omega_n^2}{s(s+2\zeta\omega_n)+K_t\omega_n^2s}=\frac{\omega_n^2}{s\left(s+2\zeta\omega_n+K_t\omega_n^2\right)} \end{aligned}

Closing the unity outer loop:

C(s)R(s)=ωn2s2+(2ζωn+Ktωn2)s+ωn2\frac{C(s)}{R(s)}=\frac{\omega_n^2}{s^2+\left(2\zeta\omega_n+K_t\omega_n^2\right)s+\omega_n^2}

Comparing with s2+2ζ′ωns+ωn2s^2+2\zeta'\omega_ns+\omega_n^2:

ζ′=ζ+Ktωn2,ωn′=ωn\zeta'=\zeta+\frac{K_t\omega_n}{2},\qquad \omega_n'=\omega_n

Effect on transient performance

  • Damping increases (ζ′>ζ\zeta'>\zeta) while ωn\omega_n is unchanged.
  • Peak overshoot reduces, since Mp=e−πζ′/1−ζ′2M_p=e^{-\pi\zeta'/\sqrt{1-\zeta'^2}} falls as ζ′\zeta' rises.
  • Settling time reduces: ts=4/(ζ′ωn)t_s=4/(\zeta'\omega_n).
  • Rise time increases slightly (slower initial rise).
  • Unlike derivative error (PD) control, it adds no zero to the closed loop, so it does not cause the extra overshoot a zero can produce.

Drawback: steady-state error

Kv=lim⁡s→0sG1(s)=ωn2ζ+Ktωn,ess(ramp)=2ζ+KtωnωnK_v=\lim_{s\to0}sG_1(s)=\frac{\omega_n}{2\zeta+K_t\omega_n},\qquad e_{ss}(\text{ramp})=\frac{2\zeta+K_t\omega_n}{\omega_n}

So the ramp error increases; it is usually offset by raising the forward gain.

Example: ωn=4\omega_n=4 rad/s, ζ=0.25\zeta=0.25 (Mp=44%M_p=44\%). With Kt=0.15K_t=0.15: ζ′=0.25+0.15×4/2=0.55\zeta'=0.25+0.15\times4/2=0.55, and MpM_p falls to about 12.6%.

  • Asked 2 times
  • 2078 Kartik · 12 marks
  • 2066 Bhadra (old course) · 16 marks

The open loop transfer function of a unity feedback control system is given by G(s) = K/[s(1+0.2s)]. Design a lead compensator such that the velocity error constant Kv = 10 sec⁻¹ and phase margin = 50°. Also draw the bode diagram for compensated system.

Answer

Lead compensator: Gc(s)=1+Ts1+αTsG_c(s)=\dfrac{1+Ts}{1+\alpha Ts}, α<1\alpha<1, maximum phase lead ϕm\phi_m with sin⁡ϕm=1−α1+α\sin\phi_m=\frac{1-\alpha}{1+\alpha} at ωm=1Tα\omega_m=\frac{1}{T\sqrt\alpha}, where it adds 10log⁡(1/α)10\log(1/\alpha) dB.

Step 1: K from KvK_v

Kv=lim⁡s→0sG(s)=K=10 s−1⇒G(s)=10s(1+0.2s)K_v=\lim_{s\to0}sG(s)=K=10\ \text{s}^{-1}\Rightarrow G(s)=\frac{10}{s(1+0.2s)}

Step 2: Uncompensated system

Asymptotic magnitude: 20 dB at ω=1\omega=1 with −20-20 dB/dec, corner at 1/0.2=51/0.2=5 rad/s, then −40-40 dB/dec. Phase ϕ=−90∘−tan⁡−10.2ω\phi=-90^\circ-\tan^{-1}0.2\omega.

Gain crossover: 10ω1+0.04ω2=1⇒0.04ω4+ω2−100=0\dfrac{10}{\omega\sqrt{1+0.04\omega^2}}=1\Rightarrow0.04\omega^4+\omega^2-100=0

ω2=−1+1+160.08=39.04⇒ωgc=6.25 rad/sPM=180∘−90∘−tan⁡−1(1.25)=90∘−51.3∘=38.7∘\begin{aligned} \omega^2&=\frac{-1+\sqrt{1+16}}{0.08}=39.04\Rightarrow\omega_{gc}=6.25\ \text{rad/s}\\ PM&=180^\circ-90^\circ-\tan^{-1}(1.25)=90^\circ-51.3^\circ=38.7^\circ \end{aligned}

Step 3: Required phase lead

ϕm=50∘−38.7∘+5∘=16.3∘\phi_m=50^\circ-38.7^\circ+5^\circ=16.3^\circ

(5° added because the crossover moves right, where the plant has more lag.)

Step 4: α

α=1−sin⁡16.3∘1+sin⁡16.3∘=1−0.2811+0.281=0.561\alpha=\frac{1-\sin16.3^\circ}{1+\sin16.3^\circ}=\frac{1-0.281}{1+0.281}=0.561

Step 5: New gain crossover

Lead gain at ωm\omega_m: 10log⁡(1/0.561)=2.5110\log(1/0.561)=2.51 dB. Find where ∣G∣=−2.51|G|=-2.51 dB:

10ω1+0.04ω2=0.561=0.749⇒ωm=7.44 rad/s\frac{10}{\omega\sqrt{1+0.04\omega^2}}=\sqrt{0.561}=0.749\Rightarrow\omega_m=7.44\ \text{rad/s}

Step 6: T

T=1ωmα=17.44×0.749=0.179≈0.18 sαT=0.561×0.179=0.10 s\begin{aligned} T&=\frac{1}{\omega_m\sqrt\alpha}=\frac{1}{7.44\times0.749}=0.179\approx0.18\ \text{s}\\ \alpha T&=0.561\times0.179=0.10\ \text{s} \end{aligned} Gc(s)=1+0.18s1+0.1sG_c(s)=\frac{1+0.18s}{1+0.1s}

Zero at 5.56 rad/s, pole at 10 rad/s. An amplifier of gain 1/α≈1.81/\alpha\approx1.8 makes up the network's attenuation, keeping Kv=10K_v=10.

Step 7: Compensated system

Gc(s)G(s)=10(1+0.18s)s(1+0.2s)(1+0.1s)G_c(s)G(s)=\frac{10(1+0.18s)}{s(1+0.2s)(1+0.1s)}

Check at the new crossover ωgc′=7.47\omega_{gc}'=7.47 rad/s:

ϕ=−90∘+tan⁡−11.345−tan⁡−11.494−tan⁡−10.747=−90∘+53.4∘−56.2∘−36.8∘=−129.6∘PM=180∘−129.6∘=50.4∘\begin{aligned} \phi&=-90^\circ+\tan^{-1}1.345-\tan^{-1}1.494-\tan^{-1}0.747\\ &=-90^\circ+53.4^\circ-56.2^\circ-36.8^\circ=-129.6^\circ\\ PM&=180^\circ-129.6^\circ=50.4^\circ \end{aligned}

GM = ∞ (phase never reaches −180°).

Bode diagram of the compensated system

Asymptotes: −20-20 dB/dec from 20 dB at ω=1\omega=1; corner 5 (pole) → −40-40; corner 5.56 (zero) → −20-20; corner 10 (pole) → −40-40 dB/dec.

ω\omega0.11257.471020100
$G$ uncomp. (dB)40.019.813.33.0−2.6-2.6−7.0-7.0
$G_cG$ (dB)40.019.913.74.60−3.7-3.7
ϕ\phi uncomp.−91∘-91^\circ−101∘-101^\circ−112∘-112^\circ−135∘-135^\circ−146∘-146^\circ−153∘-153^\circ−166∘-166^\circ−177∘-177^\circ
ϕ\phi comp.−91∘-91^\circ−97∘-97^\circ−103∘-103^\circ−120∘-120^\circ−130∘-130^\circ−138∘-138^\circ−155∘-155^\circ−175∘-175^\circ
 dB 40 |\  -20 dB/dec
    20 |  \  (w=1)
       |    \     compensated
     0 |------\---x------------- w
       |     5 \ 7.47 -20 between
       |   5.56 \__   5.56 and 10
       |          \__ -40 after 10
 deg -90|__
     -130|    \___o  PM = 50.4 deg
     -180|-------------\_______

Answer: Gc(s)=1+0.18s1+0.1sG_c(s)=\dfrac{1+0.18s}{1+0.1s} (α ≈ 0.56) with K = 10 gives Kv = 10 s⁻¹ and PM ≈ 50° at ω ≈ 7.5 rad/s.

  • Asked 2 times
  • 2076 Chaitra · 4 marks
  • 2080 Bhadra · 4 marks

Discuss principle and working performance of PD controller.

Answer

A PD (proportional–derivative) controller produces a control signal proportional to the error plus the rate of change of the error:

u(t)=Kpe(t)+Kdde(t)dt,Gc(s)=Kp+Kds=Kp(1+Tds), Td=KdKpu(t)=K_pe(t)+K_d\frac{de(t)}{dt},\qquad G_c(s)=K_p+K_ds=K_p(1+T_ds),\ T_d=\frac{K_d}{K_p}
          +-->[ Kp ]-----+
 e(t) ----|              (+)---> u(t)
          +-->[ Kd d/dt ]+

Principle

  • The P part acts on the present error.
  • The D part acts on the trend of the error. When the error is falling fast (output rushing toward the set point), de/dtde/dt is negative and the D term reduces the drive before the output overshoots. So D gives anticipatory action and extra damping.
  • It adds a zero at s=−Kp/Kds=-K_p/K_d to the open-loop transfer function, i.e. phase lead.

Working performance (second-order system)

For G(s)=ωn2s(s+2ζωn)G(s)=\dfrac{\omega_n^2}{s(s+2\zeta\omega_n)} with PD control (taking Kp=1K_p=1):

C(s)R(s)=ωn2(1+Tds)s2+(2ζωn+ωn2Td)s+ωn2\frac{C(s)}{R(s)}=\frac{\omega_n^2(1+T_ds)}{s^2+(2\zeta\omega_n+\omega_n^2T_d)s+\omega_n^2} ζ′=ζ+ωnTd2\zeta'=\zeta+\frac{\omega_nT_d}{2}
PropertyEffect of PD
Damping ratioincreases
Peak overshootdecreases
Rise time, settling timedecrease (faster)
Type of systemunchanged
Steady-state errorunchanged (for given KpK_p)
Stability marginimproves
Noiseamplified (drawback)

Example: with ωn=2\omega_n=2 rad/s and ζ=0.25\zeta=0.25, choosing Td=0.3T_d=0.3 s gives ζ′=0.25+0.3=0.55\zeta'=0.25+0.3=0.55, cutting the step overshoot from about 44% to about 16% (the added zero keeps it a little above the 12.6% of a plain ζ = 0.55 system).

  • Asked 2 times
  • 2075 Chaitra · 12 marks
  • 2080 Bhadra · 12 marks

Consider a system with open loop transfer function G(s)H(s) = 2/[s(s+1)(s+2)]. It is desired to compensate the system so that the static velocity error constant Kv is 5 per second, the phase margin is at least 40° and gain margin is at least 10 dB. Determine transfer function of appropriate lag compensator.

Answer

A lag compensator Gc(s)=Kc1+Ts1+βTsG_c(s)=K_c\dfrac{1+Ts}{1+\beta Ts} (β>1\beta>1) is used here: it gives high gain at low frequency (for KvK_v) while its attenuation of 20log⁡β20\log\beta at high frequency lowers the gain crossover to a frequency where the plant already has enough phase margin.

Step 1: Gain for the required KvK_v

G(s)H(s)=2s(s+1)(s+2)=1s(1+s)(1+0.5s)⇒Kv=lim⁡s→0sGH=1G(s)H(s)=\frac{2}{s(s+1)(s+2)}=\frac{1}{s(1+s)(1+0.5s)}\Rightarrow K_v=\lim_{s\to0}sGH=1

To get Kv=5K_v=5 the compensator must give a DC gain Kc=5K_c=5:

G1(s)=5s(1+s)(1+0.5s)G_1(s)=\frac{5}{s(1+s)(1+0.5s)}

Step 2: Uncompensated (gain-adjusted) system

∣G1∣=5ω1+ω21+0.25ω2,ϕ=−90∘−tan⁡−1ω−tan⁡−10.5ω|G_1|=\frac{5}{\omega\sqrt{1+\omega^2}\sqrt{1+0.25\omega^2}},\qquad\phi=-90^\circ-\tan^{-1}\omega-\tan^{-1}0.5\omega
ω\omega0.050.10.4911.411.80
$G_1$ (dB)40.033.919.010.0
ϕ\phi−94.3∘-94.3^\circ−98.6∘-98.6^\circ−129.9∘-129.9^\circ−161.6∘-161.6^\circ−180∘-180^\circ−193∘-193^\circ
  • ωpc=2=1.414\omega_{pc}=\sqrt2=1.414 rad/s, GM =−4.4=-4.4 dB; ωgc=1.80\omega_{gc}=1.80 rad/s, PM =−13∘=-13^\circ: unstable.

Step 3: New gain crossover frequency

Required phase at ωgc′\omega_{gc}': −180∘+40∘+ϵ-180^\circ+40^\circ+\epsilon with ϵ=10∘\epsilon=10^\circ (for the lag of the compensator), i.e. −130∘-130^\circ:

−90∘−tan⁡−1ω−tan⁡−10.5ω=−130∘⇒ωgc′=0.49 rad/s-90^\circ-\tan^{-1}\omega-\tan^{-1}0.5\omega=-130^\circ\Rightarrow\omega_{gc}'=0.49\ \text{rad/s}

Step 4: β

At 0.49 rad/s, ∣G1∣=18.95|G_1|=18.95 dB =8.86=8.86. The lag network must cut this to 0 dB:

20log⁡β=18.95 dB⇒β=8.8620\log\beta=18.95\ \text{dB}\Rightarrow\beta=8.86

Step 5: T

Place the zero one decade below the new crossover so the compensator lag there is small:

1T=ωgc′10≈0.05 rad/s⇒T=20 s,βT=177 s\frac1T=\frac{\omega_{gc}'}{10}\approx0.05\ \text{rad/s}\Rightarrow T=20\ \text{s},\quad\beta T=177\ \text{s}

Pole at 1/βT=0.005651/\beta T=0.00565 rad/s.

Gc(s)=5 1+20s1+177s=0.565 s+0.05s+0.00565G_c(s)=5\,\frac{1+20s}{1+177s}=0.565\,\frac{s+0.05}{s+0.00565}

Step 6: Check

Gc(s)G(s)H(s)=5(1+20s)s(1+s)(1+0.5s)(1+177s)G_c(s)G(s)H(s)=\frac{5(1+20s)}{s(1+s)(1+0.5s)(1+177s)}
QuantityBeforeAfter
KvK_v1 (5 with gain)5 s⁻¹
ωgc\omega_{gc}1.80 rad/s0.494 rad/s
PM−13∘-13^\circ44.7∘44.7^\circ
ωpc\omega_{pc}1.414 rad/s1.37 rad/s
GM−4.4-4.4 dB13.9 dB
 dB 40 |\__ uncompensated
       |   \___
     0 |---x---\------------ w
       | 0.49   \ 1.8
       | compensated curve is
       | 19 dB lower above 0.05

All specifications are met (PM ≥ 40°, GM ≥ 10 dB, Kv=5K_v=5). The price is a lower bandwidth (slower response).

Answer: Gc(s)=5 1+20s1+177sG_c(s)=5\,\dfrac{1+20s}{1+177s} (β ≈ 8.9, T = 20 s): Kv = 5 s⁻¹, PM ≈ 45°, GM ≈ 14 dB.

  • Asked 2 times
  • 2074 Chaitra · 12 marks
  • 2078 Kartik · 12 marks

Design a suitable phase lag compensating network for G(s) = K/[s(1+0.1s)(1+0.2s)] to meet the following specifications: Kv = 30 sec⁻¹, P.M ≥ 40°.

Answer

A phase-lag network Gc(s)=1+Ts1+βTsG_c(s)=\dfrac{1+Ts}{1+\beta Ts} (β>1\beta>1) has unity DC gain and attenuates high frequencies by 20log⁡β20\log\beta dB. It lowers the gain crossover frequency to where the plant phase gives the required PM, while KvK_v is kept by the gain KK.

Step 1: K from KvK_v

Kv=lim⁡s→0sG(s)=K=30⇒G(s)=30s(1+0.1s)(1+0.2s)K_v=\lim_{s\to0}sG(s)=K=30\Rightarrow G(s)=\frac{30}{s(1+0.1s)(1+0.2s)}

Step 2: Uncompensated system

Corners 5 and 10 rad/s. Asymptotes: −20-20 dB/dec (29.5 dB at ω=1\omega=1) to 5, −40-40 to 10, −60-60 after.

ϕ(ω)=−90∘−tan⁡−10.1ω−tan⁡−10.2ω\phi(\omega)=-90^\circ-\tan^{-1}0.1\omega-\tan^{-1}0.2\omega
ω\omega0.10.2512.4657.079.77
$G$ (dB)49.541.629.320.511.6
ϕ\phi−91.7∘-91.7^\circ−94.3∘-94.3^\circ−107∘-107^\circ−130∘-130^\circ−161.6∘-161.6^\circ−180∘-180^\circ−197∘-197^\circ

ωgc=9.77\omega_{gc}=9.77 rad/s, PM =−17.2∘=-17.2^\circ; ωpc=7.07\omega_{pc}=7.07 rad/s, GM =−6=-6 dB. The system is unstable.

Step 3: New gain crossover

Required phase: −180∘+40∘+10∘=−130∘-180^\circ+40^\circ+10^\circ=-130^\circ (10∘10^\circ allowance for the lag network's own phase lag).

tan⁡−10.1ω+tan⁡−10.2ω=40∘⇒ωgc′=2.46 rad/s\tan^{-1}0.1\omega+\tan^{-1}0.2\omega=40^\circ\Rightarrow\omega_{gc}'=2.46\ \text{rad/s}

Step 4: β

At 2.46 rad/s, ∣G∣=20.53|G|=20.53 dB =10.63=10.63:

20log⁡β=20.53⇒β≈10.620\log\beta=20.53\Rightarrow\beta\approx10.6

Step 5: T

Zero about one decade below ωgc′\omega_{gc}':

1T≈2.4610≈0.25 rad/s⇒T=4 s,βT=42.5 s\frac1T\approx\frac{2.46}{10}\approx0.25\ \text{rad/s}\Rightarrow T=4\ \text{s},\quad\beta T=42.5\ \text{s}

Pole at 1/βT=0.02351/\beta T=0.0235 rad/s.

Gc(s)=1+4s1+42.5sG_c(s)=\frac{1+4s}{1+42.5s}

Network realisation (passive R-C lag): T=R2CT=R_2C, β=R1+R2R2=10.6\beta=\frac{R_1+R_2}{R_2}=10.6; e.g. C=10 μC=10\ \muF, R2=400 kΩR_2=400\ \text{k}\Omega, R1=3.84 MΩR_1=3.84\ \text{M}\Omega.

Step 6: Compensated system

Gc(s)G(s)=30(1+4s)s(1+0.1s)(1+0.2s)(1+42.5s)G_c(s)G(s)=\frac{30(1+4s)}{s(1+0.1s)(1+0.2s)(1+42.5s)}
ω\omega0.02350.2512.476.83
$G_cG$ (dB)59.224.09.1
ϕ\phi−130∘-130^\circ−134∘-134^\circ−120∘-120^\circ−135.4∘-135.4^\circ−180∘-180^\circ
PM=180∘−135.4∘=44.6∘ (≥40∘),GM=13.9 dBPM=180^\circ-135.4^\circ=44.6^\circ\ (\ge40^\circ),\qquad GM=13.9\ \text{dB}
ItemUncompensatedCompensated
KvK_v3030
ωgc\omega_{gc}9.77 rad/s2.47 rad/s
PM−17.2∘-17.2^\circ44.6∘44.6^\circ
GM−6-6 dB13.9 dB

Answer: Gc(s)=1+4s1+42.5sG_c(s)=\dfrac{1+4s}{1+42.5s} (β ≈ 10.6, T = 4 s) with K = 30 gives Kv = 30 s⁻¹ and PM ≈ 44.6°.

  • Asked 2 times
  • 2073 Shrawan · 10 marks
  • 2067 Asar (old course) · 16 marks

The open loop transfer function of a unity feedback system is G(s) = K/[s(1+0.2s)]. It is required that Kv ≥ 20 sec⁻¹ and phase margin (φm) = 44°. Design a lead compensating network to satisfy the required specifications.

Answer

Lead network: Gc(s)=1+Ts1+αTsG_c(s)=\dfrac{1+Ts}{1+\alpha Ts}, α<1\alpha<1; maximum lead sin⁡ϕm=1−α1+α\sin\phi_m=\frac{1-\alpha}{1+\alpha} at ωm=1Tα\omega_m=\frac{1}{T\sqrt\alpha}, where it adds 10log⁡(1/α)10\log(1/\alpha) dB (after the 1/α1/\alpha attenuation is made up by an amplifier).

Step 1: K

Kv=lim⁡s→0sG(s)=K=20 s−1⇒G(s)=20s(1+0.2s)K_v=\lim_{s\to0}sG(s)=K=20\ \text{s}^{-1}\Rightarrow G(s)=\frac{20}{s(1+0.2s)}

Step 2: Uncompensated system

Asymptotes: −20-20 dB/dec through 26 dB at ω=1\omega=1; corner at 5 rad/s; then −40-40 dB/dec.

20ω1+0.04ω2=1⇒0.04ω4+ω2−400=0⇒ω2=88.3, ωgc=9.40 rad/s\frac{20}{\omega\sqrt{1+0.04\omega^2}}=1\Rightarrow0.04\omega^4+\omega^2-400=0\Rightarrow\omega^2=88.3,\ \omega_{gc}=9.40\ \text{rad/s} PM=90∘−tan⁡−1(0.2×9.40)=90∘−62.0∘=28.0∘PM=90^\circ-\tan^{-1}(0.2\times9.40)=90^\circ-62.0^\circ=28.0^\circ

Step 3: Phase lead needed

ϕm=44∘−28.0∘+5∘=21.0∘\phi_m=44^\circ-28.0^\circ+5^\circ=21.0^\circ

Step 4: α

α=1−sin⁡21∘1+sin⁡21∘=1−0.3581+0.358=0.473\alpha=\frac{1-\sin21^\circ}{1+\sin21^\circ}=\frac{1-0.358}{1+0.358}=0.473

Step 5: New gain crossover

The network adds 10log⁡(1/0.473)=3.2510\log(1/0.473)=3.25 dB at ωm\omega_m. Find where ∣G∣=−3.25|G|=-3.25 dB (=α=0.688=\sqrt\alpha=0.688):

20ω1+0.04ω2=0.688⇒ωm=11.55 rad/s\frac{20}{\omega\sqrt{1+0.04\omega^2}}=0.688\Rightarrow\omega_m=11.55\ \text{rad/s}

Step 6: T

T=1ωmα=111.55×0.688=0.126 s,αT=0.0595≈0.06 sT=\frac{1}{\omega_m\sqrt\alpha}=\frac{1}{11.55\times0.688}=0.126\ \text{s},\qquad\alpha T=0.0595\approx0.06\ \text{s}

Zero at 1/T=7.941/T=7.94 rad/s, pole at 1/αT=16.81/\alpha T=16.8 rad/s.

Gc(s)=1+0.126s1+0.06sG_c(s)=\frac{1+0.126s}{1+0.06s}

Passive R-C network: EoEi=α1+Ts1+αTs\dfrac{E_o}{E_i}=\alpha\dfrac{1+Ts}{1+\alpha Ts} with T=R1CT=R_1C, α=R2R1+R2\alpha=\dfrac{R_2}{R_1+R_2}; an amplifier of gain 1/α≈2.11/\alpha\approx2.1 restores KvK_v. E.g. C=1 μC=1\ \muF, R1=126 kΩR_1=126\ \text{k}\Omega, R2=113 kΩR_2=113\ \text{k}\Omega.

Step 7: Check

Gc(s)G(s)=20(1+0.126s)s(1+0.2s)(1+0.06s)G_c(s)G(s)=\frac{20(1+0.126s)}{s(1+0.2s)(1+0.06s)}
ω\omega157.949.411.5416.8
$G$ (dB)25.99.02.60
$G_cG$ (dB)25.910.14.72.6
ϕcomp\phi_{comp}−97.6∘-97.6^\circ−119.5∘-119.5^\circ−128.3∘-128.3^\circ−131.6∘-131.6^\circ−135.8∘-135.8^\circ−143.9∘-143.9^\circ
PM=180∘−135.8∘=44.2∘≈44∘PM=180^\circ-135.8^\circ=44.2^\circ\approx44^\circ

GM = ∞ (phase stays above −180∘-180^\circ).

ItemBeforeAfter
KvK_v2020
ωgc\omega_{gc}9.40 rad/s11.54 rad/s
PM28°44.2°

Answer: Gc(s)=1+0.126s1+0.06sG_c(s)=\dfrac{1+0.126s}{1+0.06s} (α = 0.473) with K = 20 gives Kv = 20 s⁻¹ and PM ≈ 44°.

  • Asked 2 times
  • 2070 Chaitra (old course) · 6 marks
  • 2066 Jestha (old course) · 6 marks

Discuss in brief the use of PID controllers in control system.

Answer

A PID controller combines proportional, integral and derivative actions on the error. It is the most widely used controller in industry (process control of temperature, pressure, flow, level; motor speed and position control) because it can fix both transient and steady-state behaviour with three tunable gains.

u(t)=Kpe(t)+Ki∫0te dt+Kddedt,Gc(s)=Kp+Kis+Kds=Kds2+Kps+Kisu(t)=K_pe(t)+K_i\int_0^te\,dt+K_d\frac{de}{dt},\qquad G_c(s)=K_p+\frac{K_i}{s}+K_ds=\frac{K_ds^2+K_ps+K_i}{s}
        +-->[  Kp  ]---+
 e ---->+-->[ Ki/s ]--(+)--> u ---> Plant ---> c
        +-->[ Kd s ]---+

PID adds one pole at the origin and two zeros to the open loop.

Role of each term

TermActs onMain effect
Ppresent errorfaster response, reduces (not removes) error; too high → oscillation
Ipast (accumulated) errorremoves steady-state error (type +1); may raise overshoot
Dfuture trend (rate)adds damping, reduces overshoot and settling time, improves stability

Uses / advantages

  1. Zero steady-state error to step inputs (and ramp for a type-1 plant) due to the integral term.
  2. Good transient response: derivative action limits overshoot introduced by P and I.
  3. Better relative stability (PD part acts like a lead compensator, PI part like a lag; PID behaves like a lag-lead compensator).
  4. Simple tuning without an exact model, e.g. Ziegler–Nichols rules: Kp=0.6KcrK_p=0.6K_{cr}, Ti=0.5PcrT_i=0.5P_{cr}, Td=0.125PcrT_d=0.125P_{cr}.
  5. Available as standard industrial hardware and in PLC/DCS software.

Limitations

  • D action amplifies noise, so a filtered derivative is used.
  • Integral wind-up when the actuator saturates (anti-wind-up needed).

Example: in a boiler temperature loop, P gives quick heating, I removes the final offset from the set point, and D stops the temperature overshooting when the setpoint is changed.

  • 2081 Baisakh · 4 marks
  • 2081 Baisakh · 4 marks

Discuss briefly about PI controller and its effect on time response.

Answer

A PI (proportional–integral) controller produces an output proportional to the error plus the integral of the error:

u(t)=Kpe(t)+Ki∫0te(t) dt,Gc(s)=Kp+Kis=Kp(s+KiKp)su(t)=K_pe(t)+K_i\int_0^te(t)\,dt,\qquad G_c(s)=K_p+\frac{K_i}{s}=\frac{K_p\left(s+\frac{K_i}{K_p}\right)}{s}

It adds a pole at the origin and a zero at s=−Ki/Kps=-K_i/K_p to the open-loop transfer function.

          +-->[ Kp ]----+
 e(t) ----|             (+)--> u(t)
          +-->[Ki/s]----+

Effect on time response

Consider a plant G(s)=ωn2s(s+2ζωn)G(s)=\dfrac{\omega_n^2}{s(s+2\zeta\omega_n)} with PI control. The open loop becomes

Gc(s)G(s)=ωn2(Kps+Ki)s2(s+2ζωn)G_c(s)G(s)=\frac{\omega_n^2(K_ps+K_i)}{s^2(s+2\zeta\omega_n)}
  1. System type increases by one (type 1 to type 2). Step error stays zero and the ramp steady-state error becomes zero (Kv=∞K_v=\infty). This is the main benefit: steady-state accuracy improves.
  2. Order increases (2nd to 3rd order). The characteristic equation s3+2ζωns2+Kpωn2s+Kiωn2=0s^3+2\zeta\omega_ns^2+K_p\omega_n^2s+K_i\omega_n^2=0 is stable only if 2ζωnKp>Ki2\zeta\omega_nK_p>K_i; too large KiK_i causes instability.
  3. Transient response: the integrator adds phase lag, so damping decreases. Overshoot increases and settling time becomes longer; rise time usually decreases a little.
  4. The PI controller acts as a low-pass filter (like a lag compensator), so it reduces high-frequency noise.
PropertyEffect of PI
Steady-state erroreliminated / reduced
Type, orderboth increase by 1
Overshootincreases
Settling timeincreases
Relative stabilityreduces
Noiseattenuated

Example: a speed-control loop with only P control settles with a small offset; adding the integral term removes the offset, but the speed may overshoot more before settling.

  • 2082 Baisakh · 4 marks

Which type of controller do you recommend to reduce the steady state error? Explain.

Answer

To reduce (or eliminate) steady-state error, an integral-type controller, usually a PI controller (or PID), is recommended. In frequency-domain design, the equivalent is a lag compensator.

Why integral action

The steady-state error of a unity feedback system depends on the system type and error constants:

ess(step)=11+Kp,ess(ramp)=1Kve_{ss}(\text{step})=\frac{1}{1+K_p},\qquad e_{ss}(\text{ramp})=\frac{1}{K_v}

A PI controller

Gc(s)=Kp+Kis=Kp(s+Ki/Kp)sG_c(s)=K_p+\frac{K_i}{s}=\frac{K_p(s+K_i/K_p)}{s}

adds a pole at the origin, so the system type increases by one:

Plant typeWithout PIWith PI
0finite step errorzero step error
1finite ramp errorzero ramp error

Physically, the integral keeps growing while any error remains, so the controller output keeps changing until the error is exactly zero.

The zero at −Ki/Kp-K_i/K_p is placed close to the origin so that the extra phase lag (and the effect on transient response and stability) stays small.

Alternatives

  • Increasing proportional gain reduces the error but cannot remove it, and makes the response oscillatory.
  • Lag compensator 1+Ts1+βTs\frac{1+Ts}{1+\beta Ts} raises the low-frequency gain by β\beta (KvK_v multiplied by β) without changing the transient much.

Example: G(s)=10s+2G(s)=\frac{10}{s+2} with P control has step error 11+5=0.167\frac{1}{1+5}=0.167. With PI control the open loop has a pole at the origin and the step error becomes zero.

  • 2082 Baisakh · 12 marks

Design a lead compensator for unity feedback system with open loop transfer function G(s) = 2/[s(1+0.5s)] such that velocity error constant becomes greater than 20 sec⁻¹ and phase margin would be at least 50°.

Answer

Lead compensator: Gc(s)=Kc1+Ts1+αTsG_c(s)=K_c\dfrac{1+Ts}{1+\alpha Ts}, α<1\alpha<1, with maximum lead sin⁡ϕm=1−α1+α\sin\phi_m=\frac{1-\alpha}{1+\alpha} at ωm=1Tα\omega_m=\frac{1}{T\sqrt\alpha}, where it raises the magnitude by 10log⁡(1/α)10\log(1/\alpha) dB.

Step 1: Gain for KvK_v

Kv=lim⁡s→0s Gc(s)G(s)=Kc×2=20⇒Kc=10K_v=\lim_{s\to0}s\,G_c(s)G(s)=K_c\times2=20\Rightarrow K_c=10

Gain-adjusted plant:

G1(s)=20s(1+0.5s)G_1(s)=\frac{20}{s(1+0.5s)}

Step 2: Uncompensated (gain-adjusted) system

Asymptotes: −20-20 dB/dec (26 dB at ω=1\omega=1), corner at 2 rad/s, then −40-40 dB/dec.

20ω1+0.25ω2=1⇒0.25ω4+ω2−400=0⇒ωgc=6.17 rad/s\frac{20}{\omega\sqrt{1+0.25\omega^2}}=1\Rightarrow0.25\omega^4+\omega^2-400=0\Rightarrow\omega_{gc}=6.17\ \text{rad/s} PM=90∘−tan⁡−1(0.5×6.17)=90∘−72.0∘=18.0∘PM=90^\circ-\tan^{-1}(0.5\times6.17)=90^\circ-72.0^\circ=18.0^\circ

Step 3: Phase lead needed

ϕm=50∘−18∘+ϵ\phi_m=50^\circ-18^\circ+\epsilon

With ϵ=5∘\epsilon=5^\circ the check gives PM = 49.8° (just short), so take ϵ≈10∘\epsilon\approx10^\circ: ϕm≈42∘\phi_m\approx42^\circ.

Step 4: α

α=1−sin⁡42∘1+sin⁡42∘=0.198≈0.2(ϕm=41.8∘)\alpha=\frac{1-\sin42^\circ}{1+\sin42^\circ}=0.198\approx0.2\quad(\phi_m=41.8^\circ)

Step 5: New gain crossover

Lead gain at ωm\omega_m: 10log⁡(1/0.2)=6.9910\log(1/0.2)=6.99 dB. Find where ∣G1∣=−6.99|G_1|=-6.99 dB (=0.2=0.447)(=\sqrt{0.2}=0.447):

20ω1+0.25ω2=0.447⇒ωm=9.35 rad/s\frac{20}{\omega\sqrt{1+0.25\omega^2}}=0.447\Rightarrow\omega_m=9.35\ \text{rad/s}

Step 6: T

T=1ωmα=19.35×0.447=0.239≈0.24 s,αT=0.048 sT=\frac{1}{\omega_m\sqrt\alpha}=\frac{1}{9.35\times0.447}=0.239\approx0.24\ \text{s},\qquad\alpha T=0.048\ \text{s}

Zero at 4.17 rad/s, pole at 20.8 rad/s.

Gc(s)=10 1+0.24s1+0.048s=50 s+4.17s+20.8G_c(s)=10\,\frac{1+0.24s}{1+0.048s}=50\,\frac{s+4.17}{s+20.8}

Step 7: Check

Gc(s)G(s)=20(1+0.24s)s(1+0.5s)(1+0.048s)G_c(s)G(s)=\frac{20(1+0.24s)}{s(1+0.5s)(1+0.048s)}
ω\omega124.176.179.3720.8
$G_1$ (dB)25.117.06.30
$G_cG$ (dB)25.317.99.24.7
ϕcomp\phi_{comp}−106∘-106^\circ−115∘-115^\circ−121∘-121^\circ−123∘-123^\circ−126.1∘-126.1^\circ−141∘-141^\circ
PM=180∘−126.1∘=53.9∘ (≥50∘),GM=∞PM=180^\circ-126.1^\circ=53.9^\circ\ (\ge50^\circ),\quad GM=\infty
 dB 26 |\  -20
       |  \__ 2 (corner)
     0 |------\----x------- w
       |    6.17\  9.37 (new wgc)
       |         \__ -40
  lead adds phase between 4.17 and 20.8

Answer: Gc(s)=10 1+0.24s1+0.048sG_c(s)=10\,\dfrac{1+0.24s}{1+0.048s} (α = 0.2): Kv = 20 s⁻¹, PM ≈ 54° at ω ≈ 9.4 rad/s.

  • 2081 Bhadra · 12 marks

Design a compensator for unity feedback system with open loop transfer function G(s) = 10/[s(s+1)] such that the damping ratio would become 0.5 and natural frequency of oscillation 3 rad/sec.

Answer

The specification is on the closed-loop pole positions, so a lead compensator is designed by the root-locus method.

Step 1: Present system

CR=10s2+s+10⇒ωn=10=3.16 rad/s, ζ=12×3.16=0.158\frac{C}{R}=\frac{10}{s^2+s+10}\Rightarrow\omega_n=\sqrt{10}=3.16\ \text{rad/s},\ \zeta=\frac{1}{2\times3.16}=0.158

The response is very oscillatory (Mp≈60%M_p\approx60\%).

Step 2: Desired dominant poles

sd=−ζωn±jωn1−ζ2=−1.5±j2.598s_d=-\zeta\omega_n\pm j\omega_n\sqrt{1-\zeta^2}=-1.5\pm j2.598

Step 3: Angle deficiency

Angle of G(s)=10s(s+1)G(s)=\dfrac{10}{s(s+1)} at sd=−1.5+j2.598s_d=-1.5+j2.598:

∠sd=180∘−tan⁡−12.5981.5=120∘∠(sd+1)=180∘−tan⁡−12.5980.5=100.9∘∠G(sd)=−(120∘+100.9∘)=−220.9∘\begin{aligned} \angle s_d&=180^\circ-\tan^{-1}\frac{2.598}{1.5}=120^\circ\\ \angle(s_d+1)&=180^\circ-\tan^{-1}\frac{2.598}{0.5}=100.9^\circ\\ \angle G(s_d)&=-(120^\circ+100.9^\circ)=-220.9^\circ \end{aligned}

To satisfy the angle condition (−180∘-180^\circ), the compensator must add

ϕ=−180∘−(−220.9∘)=40.9∘ (lead)\phi=-180^\circ-(-220.9^\circ)=40.9^\circ\ \text{(lead)}

Step 4: Locate zero and pole

Choose the compensator zero at s=−1s=-1 to cancel the plant pole at −1-1. Then the pole −p-p must satisfy

−∠sd−∠(sd+p)=−180∘⇒−120∘−∠(sd+p)=−180∘⇒∠(sd+p)=60∘-\angle s_d-\angle(s_d+p)=-180^\circ\Rightarrow-120^\circ-\angle(s_d+p)=-180^\circ\Rightarrow\angle(s_d+p)=60^\circ tan⁡−12.598p−1.5=60∘⇒p−1.5=2.5981.732=1.5⇒p=3\tan^{-1}\frac{2.598}{p-1.5}=60^\circ\Rightarrow p-1.5=\frac{2.598}{1.732}=1.5\Rightarrow p=3

Check of angle added: zero +100.9∘+100.9^\circ, pole −60∘-60^\circ, net +40.9∘+40.9^\circ ✓.

Step 5: Gain KcK_c (magnitude condition)

Gc(s)G(s)=Kcs+1s+3⋅10s(s+1)=10Kcs(s+3)G_c(s)G(s)=K_c\frac{s+1}{s+3}\cdot\frac{10}{s(s+1)}=\frac{10K_c}{s(s+3)} ∣10Kcsd(sd+3)∣=1⇒10Kc=∣sd∣ ∣sd+3∣=3×3=9⇒Kc=0.9\left|\frac{10K_c}{s_d(s_d+3)}\right|=1\Rightarrow10K_c=|s_d|\,|s_d+3|=3\times3=9\Rightarrow K_c=0.9 Gc(s)=0.9 s+1s+3=0.3 1+s1+0.333sG_c(s)=0.9\,\frac{s+1}{s+3}=0.3\,\frac{1+s}{1+0.333s}

Step 6: Verification

C(s)R(s)=9s2+3s+9\frac{C(s)}{R(s)}=\frac{9}{s^2+3s+9}

ωn=3\omega_n=3 rad/s, 2ζωn=3⇒ζ=0.52\zeta\omega_n=3\Rightarrow\zeta=0.5 ✓. Closed-loop poles: −1.5±j2.598-1.5\pm j2.598.

         jw
   x sd  |  j2.598
    \    |
 ---x--o-x------- sigma
   -3 -1 0   (zero cancels pole at -1)
ItemBeforeAfter
ζ\zeta0.1580.5
ωn\omega_n3.16 rad/s3 rad/s
MpM_p60%16.3%
tst_s (2%)8 s2.67 s
KvK_v103

Note: KvK_v falls to 10×0.9/3=310\times0.9/3=3; if a higher KvK_v is needed, the bisector method or an extra lag section can be used.

Answer: Gc(s)=0.9 s+1s+3G_c(s)=0.9\,\dfrac{s+1}{s+3}, giving closed-loop poles at −1.5 ± j2.6 (ζ = 0.5, ωn = 3 rad/s).

  • 2081 Baisakh · 2+14 marks

Differentiate between lead and lag compensator. Design a suitable lead compensating network for G(s) = K/[s²(1+0.25s)] to meet the following specifications: Ka = 10 sec⁻² and P.M ≥ 40°.

Answer

Lead vs lag compensator

PointLead compensatorLag compensator
Transfer function1+Ts1+αTs\frac{1+Ts}{1+\alpha Ts}, α<1\alpha<11+Ts1+βTs\frac{1+Ts}{1+\beta Ts}, β>1\beta>1
Pole–zerozero nearer origin than polepole nearer origin than zero
Phasepositive (lead)negative (lag)
Filter typehigh-passlow-pass
Main useimproves PM, speed (transient)improves steady-state accuracy
Bandwidthincreasesdecreases
Noisemore sensitiveless sensitive

Design

Step 1: K. Ka=lim⁡s→0s2G(s)=K=10K_a=\lim_{s\to0}s^2G(s)=K=10, so G(s)=10s2(1+0.25s)G(s)=\dfrac{10}{s^2(1+0.25s)}.

Step 2: Uncompensated system.

∣G∣=10ω21+0.0625ω2,ϕ=−180∘−tan⁡−10.25ω|G|=\frac{10}{\omega^2\sqrt{1+0.0625\omega^2}},\qquad\phi=-180^\circ-\tan^{-1}0.25\omega
ω\omega12.8546.3910
$G$ (dB)19.70−7.1-7.1
ϕ\phi−194∘-194^\circ−215.5∘-215.5^\circ−225∘-225^\circ−238∘-238^\circ−248∘-248^\circ

ωgc=2.85\omega_{gc}=2.85 rad/s, PM=−35.5∘PM=-35.5^\circ: unstable (phase always below −180∘-180^\circ).

Step 3: Lead required. ϕ=40∘−(−35.5∘)+ϵ≈80∘\phi=40^\circ-(-35.5^\circ)+\epsilon\approx80^\circ or more. One lead stage gives at most about 60° usefully, and for this plant a single stage can reach a PM of only about 18.6°. So two identical lead stages are used:

Gc(s)=(1+Ts1+αTs)2G_c(s)=\left(\frac{1+Ts}{1+\alpha Ts}\right)^2

Step 4: α per stage. Take ϕm≈50∘\phi_m\approx50^\circ per stage: α=0.13\alpha=0.13, sin⁡ϕm=0.871.13⇒ϕm=50.3∘\sin\phi_m=\frac{0.87}{1.13}\Rightarrow\phi_m=50.3^\circ (total 100.7°).

Step 5: New crossover. Two stages add 20log⁡(1/0.13)=17.720\log(1/0.13)=17.7 dB at ωm\omega_m, so

10ω21+0.0625ω2=0.13⇒ωm=6.39 rad/s\frac{10}{\omega^2\sqrt{1+0.0625\omega^2}}=0.13\Rightarrow\omega_m=6.39\ \text{rad/s}

Step 6: T.

T=16.390.13=0.434 s,αT=0.0564 sT=\frac{1}{6.39\sqrt{0.13}}=0.434\ \text{s},\qquad\alpha T=0.0564\ \text{s} Gc(s)=(1+0.434s1+0.0564s)2G_c(s)=\left(\frac{1+0.434s}{1+0.0564s}\right)^2

with an amplifier of gain 1/α2=59.21/\alpha^2=59.2 to cancel the networks' attenuation (so KaK_a stays 10).

Step 7: Check at 6.39 rad/s.

ϕ=−180∘−tan⁡−11.60+2(50.3∘)=−137.3∘PM=42.7∘ (≥40∘)\begin{aligned} \phi&=-180^\circ-\tan^{-1}1.60+2(50.3^\circ)=-137.3^\circ\\ PM&=42.7^\circ\ (\ge40^\circ) \end{aligned}

GM ≈ 12.9 dB; all closed-loop poles in the LHP.

Answer: Gc(s)=(1+0.434s1+0.0564s)2G_c(s)=\left(\dfrac{1+0.434s}{1+0.0564s}\right)^2 with K = 10: Ka = 10 s⁻², PM ≈ 42.7°.

  • 2080 Bhadra · 4 marks

Describe working of derivative feedback controller in time response of second order system.

Answer

In derivative (rate/tachometer) feedback control, a signal proportional to the rate of change of the output, Kt c˙(t)K_t\,\dot c(t), is fed back negatively in an inner loop, besides the main unity feedback. In a motor position system, a tachogenerator on the shaft provides this velocity signal.

R --(+)--->(+)--> wn^2/(s(s+2z wn)) --+--> C
     ^-     ^-                        |
     |      +------[ Kt s ]<----------+
     +--------------------------------+

Analysis

Inner loop:

G1(s)=ωn2s(s+2ζωn+Ktωn2)G_1(s)=\frac{\omega_n^2}{s(s+2\zeta\omega_n+K_t\omega_n^2)}

Closed loop:

C(s)R(s)=ωn2s2+(2ζωn+Ktωn2)s+ωn2\frac{C(s)}{R(s)}=\frac{\omega_n^2}{s^2+(2\zeta\omega_n+K_t\omega_n^2)s+\omega_n^2} ζeff=ζ+Ktωn2,ωn unchanged\zeta_{eff}=\zeta+\frac{K_t\omega_n}{2},\qquad\omega_n\ \text{unchanged}

Working / effect on time response

  1. When the output moves fast towards the reference, the feedback signal Ktc˙K_t\dot c is large and subtracts from the actuating signal, braking the system early. This is extra damping.
  2. Overshoot falls as ζeff\zeta_{eff} rises; ts=4/(ζeffωn)t_s=4/(\zeta_{eff}\omega_n) falls.
  3. ωn\omega_n is unchanged and no zero is added, so the response is smooth (unlike PD, whose zero can add overshoot).
  4. Ramp error increases: Kv=ωn2ζ+KtωnK_v=\dfrac{\omega_n}{2\zeta+K_t\omega_n}; this is offset by raising the amplifier gain.
PropertyEffect
Dampingincreases
MpM_p, tst_sdecrease
ωn\omega_nsame
Ramp errorincreases

Example: ωn=5\omega_n=5 rad/s, ζ=0.2\zeta=0.2 (Mp=52.7%M_p=52.7\%). With Kt=0.12K_t=0.12: ζeff=0.2+0.3=0.5\zeta_{eff}=0.2+0.3=0.5, Mp=16.3%M_p=16.3\%.

  • 2080 Bhadra · 12 marks

Design a lead compensator for a unity feedback system with open loop transfer function G(s) = 5/[(s+2)(s+20)] such that velocity error constant at least 20 per second, phase margin at least 35° and gain margin at least 10 dB.

Answer

Assumption: as printed, G(s)=5(s+2)(s+20)G(s)=\frac{5}{(s+2)(s+20)} is type 0, so Kv=lim⁡s→0sG(s)=0K_v=\lim_{s\to0}sG(s)=0 for any gain and no lead network can give Kv=20K_v=20. The intended plant must have an integrator, so the design is done for G(s)=Ks(s+2)(s+20)G(s)=\dfrac{K}{s(s+2)(s+20)}, with the gain chosen to give Kv=20K_v=20.

Lead network: Gc(s)=1+Ts1+αTsG_c(s)=\dfrac{1+Ts}{1+\alpha Ts}, sin⁡ϕm=1−α1+α\sin\phi_m=\frac{1-\alpha}{1+\alpha} at ωm=1Tα\omega_m=\frac{1}{T\sqrt\alpha}, gain added there 10log⁡(1/α)10\log(1/\alpha).

Step 1: Gain

Kv=lim⁡s→0sG(s)=K2×20=20⇒K=800K_v=\lim_{s\to0}sG(s)=\frac{K}{2\times20}=20\Rightarrow K=800 G(s)=800s(s+2)(s+20)=20s(1+0.5s)(1+0.05s)G(s)=\frac{800}{s(s+2)(s+20)}=\frac{20}{s(1+0.5s)(1+0.05s)}

Step 2: Uncompensated system

Corners at 2 and 20 rad/s; asymptotes −20-20, −40-40, −60-60 dB/dec; 26 dB at ω=1\omega=1.

ϕ=−90∘−tan⁡−10.5ω−tan⁡−10.05ω\phi=-90^\circ-\tan^{-1}0.5\omega-\tan^{-1}0.05\omega
ω\omega13.336.036.3210.520
$G$ (dB)25.09.70−0.8-0.8
ϕ\phi−119∘-119^\circ−158.5∘-158.5^\circ−178.4∘-178.4^\circ−180∘-180^\circ−197∘-197^\circ−219∘-219^\circ

ωgc=6.03\omega_{gc}=6.03 rad/s, PM =1.6∘=1.6^\circ; ωpc=6.32\omega_{pc}=6.32 rad/s, GM =0.8=0.8 dB. Barely stable.

Step 3: Phase lead

ϕm=35∘−1.6∘+ϵ\phi_m=35^\circ-1.6^\circ+\epsilon

Because of the pole at 20 rad/s, the plant phase falls quickly as the crossover moves right. Trials with ϵ=5∘\epsilon=5^\circ–12∘12^\circ (ϕm=38∘\phi_m=38^\circ–45∘45^\circ) give PM of only 28°–33°. A larger margin is needed: take ϕm≈55∘\phi_m\approx55^\circ.

Step 4: α

α=1−sin⁡55∘1+sin⁡55∘≈0.1(ϕm=54.9∘)\alpha=\frac{1-\sin55^\circ}{1+\sin55^\circ}\approx0.1\quad(\phi_m=54.9^\circ)

Step 5: New crossover

The network adds 10log⁡10=1010\log10=10 dB, so ωm\omega_m is where ∣G∣=−10|G|=-10 dB:

20ω1+0.25ω21+0.0025ω2=0.316⇒ωm=10.5 rad/s\frac{20}{\omega\sqrt{1+0.25\omega^2}\sqrt{1+0.0025\omega^2}}=0.316\Rightarrow\omega_m=10.5\ \text{rad/s}

Step 6: T

T=110.5×0.1=0.30 s,αT=0.03 sT=\frac{1}{10.5\times\sqrt{0.1}}=0.30\ \text{s},\qquad\alpha T=0.03\ \text{s} Gc(s)=1+0.3s1+0.03sG_c(s)=\frac{1+0.3s}{1+0.03s}

(zero at 3.33 rad/s, pole at 33.3 rad/s; an amplifier of gain 1/α=101/\alpha=10 restores the DC gain).

Step 7: Check

Gc(s)G(s)=800(1+0.3s)s(s+2)(s+20)(1+0.03s)G_c(s)G(s)=\frac{800(1+0.3s)}{s(s+2)(s+20)(1+0.03s)}
ω\omega13.336.0310.462024.4
$G_cG$ (dB)25.412.76.20
ϕ\phi−104∘-104^\circ−119∘-119^\circ−128∘-128^\circ−141.9∘-141.9^\circ−170∘-170^\circ−180∘-180^\circ
PM=180∘−141.9∘=38.1∘ (≥35∘),GM=11.9 dB (≥10 dB)PM=180^\circ-141.9^\circ=38.1^\circ\ (\ge35^\circ),\qquad GM=11.9\ \text{dB}\ (\ge10\ \text{dB})
ItemBeforeAfter
KvK_v2020
PM1.6°38.1°
GM0.8 dB11.9 dB
ωgc\omega_{gc}6.0 rad/s10.5 rad/s

Answer (for G=K/[s(s+2)(s+20)]G=K/[s(s+2)(s+20)]): K = 800 and Gc(s)=1+0.3s1+0.03sG_c(s)=\dfrac{1+0.3s}{1+0.03s} give Kv = 20 s⁻¹, PM ≈ 38°, GM ≈ 12 dB.

  • 2080 Baisakh · 12 marks

For a unity feedback system with open loop transfer function G(s) = 4/[s(s+2)], design a Lead compensator such that settling time would become 2 seconds without change in maximum overshoot of the system.

Answer

The specification is in the time domain, so the lead compensator is designed by the root-locus method.

Step 1: Present system

C(s)R(s)=4s2+2s+4⇒ωn=2 rad/s, ζ=22×2=0.5\frac{C(s)}{R(s)}=\frac{4}{s^2+2s+4}\Rightarrow\omega_n=2\ \text{rad/s},\ \zeta=\frac{2}{2\times2}=0.5 ts=4ζωn=41=4 s (2% criterion),Mp=e−πζ/1−ζ2=16.3%t_s=\frac{4}{\zeta\omega_n}=\frac{4}{1}=4\ \text{s (2\% criterion)},\qquad M_p=e^{-\pi\zeta/\sqrt{1-\zeta^2}}=16.3\%

Step 2: Desired poles

Same overshoot means same ζ=0.5\zeta=0.5. New settling time 2 s:

ζωn=42=2⇒ωn=4 rad/s,ωd=41−0.25=3.464\zeta\omega_n=\frac{4}{2}=2\Rightarrow\omega_n=4\ \text{rad/s},\quad\omega_d=4\sqrt{1-0.25}=3.464 sd=−2±j3.464s_d=-2\pm j3.464

Step 3: Angle deficiency

∠sd=180∘−tan⁡−13.4642=120∘∠(sd+2)=90∘∠G(sd)=−(120∘+90∘)=−210∘ϕ=−180∘−(−210∘)=30∘\begin{aligned} \angle s_d&=180^\circ-\tan^{-1}\frac{3.464}{2}=120^\circ\\ \angle(s_d+2)&=90^\circ\\ \angle G(s_d)&=-(120^\circ+90^\circ)=-210^\circ\\ \phi&=-180^\circ-(-210^\circ)=30^\circ \end{aligned}

Step 4: Pole and zero by the bisector method

  1. Draw a horizontal line PA from P=sdP=s_d to the left, and line PO to the origin. Angle APO =180∘−60∘=120∘=180^\circ-60^\circ=120^\circ.
  2. Bisect it: the bisector PB makes 60∘60^\circ with PA (direction 240∘240^\circ).
  3. Draw PC and PD at ±ϕ/2=±15∘\pm\phi/2=\pm15^\circ from PB; they meet the real axis at the zero and pole.
zero: z=−2−3.464tan⁡75∘=−2−0.928=−2.93pole: p=−2−3.464tan⁡45∘=−2−3.464=−5.46\begin{aligned} \text{zero: }z&=-2-\frac{3.464}{\tan75^\circ}=-2-0.928=-2.93\\ \text{pole: }p&=-2-\frac{3.464}{\tan45^\circ}=-2-3.464=-5.46 \end{aligned}

Check: ∠(sd+2.93)=75∘\angle(s_d+2.93)=75^\circ, ∠(sd+5.46)=45∘\angle(s_d+5.46)=45^\circ, net lead =30∘=30^\circ ✓.

Step 5: Gain (magnitude condition)

Kc=∣sd∣ ∣sd+2∣ ∣sd+5.46∣4 ∣sd+2.93∣=4×3.464×4.8994×3.586=4.73K_c=\frac{|s_d|\,|s_d+2|\,|s_d+5.46|}{4\,|s_d+2.93|}=\frac{4\times3.464\times4.899}{4\times3.586}=4.73 Gc(s)=4.73 s+2.93s+5.46G_c(s)=4.73\,\frac{s+2.93}{s+5.46}

Step 6: Compensated system

Gc(s)G(s)=18.93(s+2.93)s(s+2)(s+5.46)G_c(s)G(s)=\frac{18.93(s+2.93)}{s(s+2)(s+5.46)}

Closed-loop poles: −2±j3.464-2\pm j3.464 (desired) and a third pole at −3.46-3.46, which lies close to the zero at −2.93-2.93, so its effect is small and the response is dominated by the desired pair.

          jw
    sd x  | j3.46
       \  |
 --x---o--x---x--- sigma
 -5.46 -2.93 -2  0
ItemBeforeAfter
ζ\zeta0.50.5
ωn\omega_n2 rad/s4 rad/s
tst_s4 s2 s
MpM_p16.3%≈16.3%
KvK_v25.07

Alternative: placing the zero at −2-2 (cancelling the plant pole) needs a pole at −4-4 and gives Gc=4s+2s+4G_c=4\frac{s+2}{s+4} with closed loop exactly 16s2+4s+16\frac{16}{s^2+4s+16}.

Answer: Gc(s)=4.73 s+2.93s+5.46G_c(s)=4.73\,\dfrac{s+2.93}{s+5.46}; dominant poles at −2 ± j3.46 give ts = 2 s with the same ζ = 0.5 (Mp ≈ 16%).

  • 2079 Bhadra · 12 marks

For a unity feedback system with open loop transfer function G(s) = 1.06/[s(s+1)(s+2)], design a lag compensator such that steady state error for ramp input would be less than 0.2 without significant change in transients.

Answer

Since the transient must stay almost the same while KvK_v is raised, a lag compensator by the root-locus method is used: a pole–zero pair very close to the origin raises the low-frequency gain but hardly moves the dominant poles.

Gc(s)=K^c s+1Ts+1βT,β>1G_c(s)=\hat K_c\,\frac{s+\frac1T}{s+\frac{1}{\beta T}},\qquad\beta>1

Step 1: Present system

Kv=lim⁡s→0sG(s)=1.061×2=0.53,ess=1Kv=1.89K_v=\lim_{s\to0}sG(s)=\frac{1.06}{1\times2}=0.53,\qquad e_{ss}=\frac{1}{K_v}=1.89

Characteristic equation: s3+3s2+2s+1.06=0s^3+3s^2+2s+1.06=0; roots −2.339-2.339 and −0.3307±j0.5864-0.3307\pm j0.5864.

Dominant poles: ωn=0.673\omega_n=0.673 rad/s, ζ=0.491\zeta=0.491.

Step 2: Required KvK_v and β

ess<0.2⇒Kv>5⇒β≥50.53=9.4; take β=10e_{ss}<0.2\Rightarrow K_v>5\Rightarrow\beta\ge\frac{5}{0.53}=9.4;\ \text{take}\ \beta=10

Step 3: Place the pole–zero pair

Put the zero and pole close to the origin, with ratio 10:

zero at −0.05,pole at −0.005⇒T=20 s, βT=200 s\text{zero at }-0.05,\quad\text{pole at }-0.005\Rightarrow T=20\ \text{s},\ \beta T=200\ \text{s}

Angle added at s1=−0.3307+j0.5864s_1=-0.3307+j0.5864:

∠(s1+0.05)−∠(s1+0.005)=115.6∘−119.0∘≈−3.5∘\angle(s_1+0.05)-\angle(s_1+0.005)=115.6^\circ-119.0^\circ\approx-3.5^\circ

This is less than about 5°, so the root locus near the dominant poles changes very little.

Step 4: Gain K^c\hat K_c

The new dominant poles lie on the compensated root locus with almost the same ζ (≈0.49): s1′=−0.312±j0.551s_1'=-0.312\pm j0.551. Magnitude condition:

K^c=∣s1′∣ ∣s1′+1∣ ∣s1′+2∣ ∣s1′+0.005∣1.06 ∣s1′+0.05∣=0.966\hat K_c=\frac{|s_1'|\,|s_1'+1|\,|s_1'+2|\,|s_1'+0.005|}{1.06\,|s_1'+0.05|}=0.966 Gc(s)=0.966 s+0.05s+0.005=9.66 1+20s1+200sG_c(s)=0.966\,\frac{s+0.05}{s+0.005}=9.66\,\frac{1+20s}{1+200s}

Step 5: Check

Kv=lim⁡s→0sGcG=0.966×10×0.53=5.12⇒ess=15.12=0.195<0.2K_v=\lim_{s\to0}sG_cG=0.966\times10\times0.53=5.12\Rightarrow e_{ss}=\frac{1}{5.12}=0.195<0.2

The requirement is met.

Closed-loop poles of the compensated system: −0.312±j0.551-0.312\pm j0.551, −2.326-2.326, −0.055-0.055. The pole at −0.055-0.055 is almost cancelled by the zero at −0.05-0.05, so it adds only a small, slow tail.

ItemBeforeAfter
Dominant poles−0.331±j0.586-0.331\pm j0.586−0.312±j0.551-0.312\pm j0.551
ζ\zeta0.4910.493
ωn\omega_n0.673 rad/s0.633 rad/s
KvK_v0.535.12
Ramp esse_{ss}1.890.195

Transient response is nearly unchanged (slightly slower, since ωn\omega_n dropped about 6%).

Answer: Gc(s)=0.966 s+0.05s+0.005G_c(s)=0.966\,\dfrac{s+0.05}{s+0.005} gives Kv = 5.12 s⁻¹, ramp error ≈ 0.195 (< 0.2), with almost the same damping.

  • 2078 Bhadra · 12 marks

Design a suitable phase lead compensating network for G(s) = 4/[s(s+2)] to meet the following specification: Kv = 20 sec⁻¹, P.M ≥ 50°.

Answer

Lead network Gc(s)=Kc1+Ts1+αTsG_c(s)=K_c\dfrac{1+Ts}{1+\alpha Ts} (α<1\alpha<1), maximum phase lead sin⁡ϕm=1−α1+α\sin\phi_m=\frac{1-\alpha}{1+\alpha} at ωm=1Tα\omega_m=\frac{1}{T\sqrt\alpha}.

Step 1: Gain for KvK_v

Kv=lim⁡s→0s Kc4s(s+2)=2Kc=20⇒Kc=10K_v=\lim_{s\to0}s\,K_c\frac{4}{s(s+2)}=2K_c=20\Rightarrow K_c=10 G1(s)=40s(s+2)=20s(1+0.5s)G_1(s)=\frac{40}{s(s+2)}=\frac{20}{s(1+0.5s)}

Step 2: Bode plot of G1G_1

20log20 = 26 dB at ω=1\omega=1, slope −20-20 dB/dec, corner 2 rad/s, then −40-40 dB/dec.

40ωω2+4=1⇒ω4+4ω2−1600=0⇒ωgc=6.17 rad/s\frac{40}{\omega\sqrt{\omega^2+4}}=1\Rightarrow\omega^4+4\omega^2-1600=0\Rightarrow\omega_{gc}=6.17\ \text{rad/s} PM=180∘−90∘−tan⁡−16.172=90∘−72.0∘=18.0∘PM=180^\circ-90^\circ-\tan^{-1}\frac{6.17}{2}=90^\circ-72.0^\circ=18.0^\circ

Step 3: Required lead

ϕm=50∘−18∘+8∘=40∘\phi_m=50^\circ-18^\circ+8^\circ=40^\circ

Step 4: α

α=1−sin⁡40∘1+sin⁡40∘=1−0.6431+0.643=0.217≈0.22\alpha=\frac{1-\sin40^\circ}{1+\sin40^\circ}=\frac{1-0.643}{1+0.643}=0.217\approx0.22

Step 5: New crossover

Gain added at ωm\omega_m: 10log⁡(1/0.22)=6.5810\log(1/0.22)=6.58 dB. Solve ∣G1(jω)∣=−6.58|G_1(j\omega)|=-6.58 dB =0.22=0.469=\sqrt{0.22}=0.469:

40ωω2+4=0.469⇒ωm=9.13 rad/s\frac{40}{\omega\sqrt{\omega^2+4}}=0.469\Rightarrow\omega_m=9.13\ \text{rad/s}

Step 6: T

T=1ωmα=19.13×0.469=0.233 s,αT=0.0513 sT=\frac{1}{\omega_m\sqrt\alpha}=\frac{1}{9.13\times0.469}=0.233\ \text{s},\qquad\alpha T=0.0513\ \text{s}

Zero at 1/T=4.281/T=4.28 rad/s, pole at 1/αT=19.51/\alpha T=19.5 rad/s.

Gc(s)=10 1+0.233s1+0.0513s=45.4 s+4.28s+19.5G_c(s)=10\,\frac{1+0.233s}{1+0.0513s}=45.4\,\frac{s+4.28}{s+19.5}

Step 7: Check

Gc(s)G(s)=40(1+0.233s)s(s+2)(1+0.0513s)G_c(s)G(s)=\frac{40(1+0.233s)}{s(s+2)(1+0.0513s)}

At ωgc′=9.11\omega_{gc}'=9.11 rad/s:

ϕ=−90∘−tan⁡−19.112+tan⁡−1(2.12)−tan⁡−1(0.467)=−90∘−77.62∘+64.78∘−25.05∘=−127.89∘PM=180∘−127.9∘≈52∘ (≥50∘)\begin{aligned} \phi&=-90^\circ-\tan^{-1}\frac{9.11}{2}+\tan^{-1}(2.12)-\tan^{-1}(0.467)\\ &=-90^\circ-77.62^\circ+64.78^\circ-25.05^\circ=-127.89^\circ\\ PM&=180^\circ-127.9^\circ\approx52^\circ\ (\ge50^\circ) \end{aligned}

GM = ∞.

 dB 26 |\ -20       uncompensated
       |  \__ 2   /
     0 |------\--x---x------ w
       |    6.17   9.11 compensated
       |           \__ -40
ItemBeforeAfter
KvK_v220
ωgc\omega_{gc}6.17 rad/s (with gain 10)9.11 rad/s
PM18°≈52°

Answer: Gc(s)=10 1+0.233s1+0.0513sG_c(s)=10\,\dfrac{1+0.233s}{1+0.0513s} (α ≈ 0.22) gives Kv = 20 s⁻¹ and PM ≈ 52°.

  • 2078 Kartik · 4 marks

What kind of controller would you recommend to bring changes in transient properties of a system and how?

Answer

To change (improve) the transient properties of a system (overshoot, settling time, rise time), a controller with derivative action is recommended: a PD controller (or derivative/rate feedback), or in compensator language a lead compensator.

How PD control changes the transient

Gc(s)=Kp+Kds=Kp(1+Tds)G_c(s)=K_p+K_ds=K_p(1+T_ds)

For the standard plant G(s)=ωn2s(s+2ζωn)G(s)=\dfrac{\omega_n^2}{s(s+2\zeta\omega_n)} (with Kp=1K_p=1):

C(s)R(s)=ωn2(1+Tds)s2+(2ζωn+Tdωn2)s+ωn2,ζ′=ζ+Tdωn2\frac{C(s)}{R(s)}=\frac{\omega_n^2(1+T_ds)}{s^2+(2\zeta\omega_n+T_d\omega_n^2)s+\omega_n^2},\qquad\zeta'=\zeta+\frac{T_d\omega_n}{2}
  1. Damping increases (ζ′ > ζ), so peak overshoot falls.
  2. Settling time falls (ts=4/ζ′ωnt_s=4/\zeta'\omega_n).
  3. Anticipation: D action responds to the rate of change of error, applying correction before the error becomes large, like braking a car before the stop line.
  4. The added zero gives phase lead, which raises the phase margin and bandwidth, so the response is faster and more stable.
  5. Steady-state error is not changed (type unchanged), and steady-state is handled separately by I action if needed.
Transient specEffect of PD
Overshoot MpM_pdecreases
Settling time tst_sdecreases
Rise time trt_rdecreases slightly
Stabilityimproves

Caution: derivative action amplifies high-frequency noise, so KdK_d is kept moderate or a filtered derivative is used.

Example: ωn=4\omega_n=4 rad/s, ζ=0.25\zeta=0.25 (Mp=44%M_p=44\%, ts=4t_s=4 s). With Td=0.15T_d=0.15 s: ζ′=0.25+0.3=0.55\zeta'=0.25+0.3=0.55, ts≈4/(ζ′ωn)=1.8t_s\approx4/(\zeta'\omega_n)=1.8 s, and the overshoot drops from 44% to about 16% (a little above the 12.6% of a pure ζ = 0.55 system, because of the added zero).

  • 2076 Chaitra · 12 marks

For a unity feedback system with feed forward transfer function G(s) = 10/[s(s+1)], design a lead compensator such that the settling time of the system will become 2 sec and maximum percent overshoot 5%.

Answer

Time-domain specifications are met by placing the dominant closed-loop poles with a lead compensator designed on the root locus.

Step 1: Present system

CR=10s2+s+10: ωn=3.16 rad/s, ζ=0.158, Mp=60.5%, ts=40.5=8 s\frac{C}{R}=\frac{10}{s^2+s+10}:\ \omega_n=3.16\ \text{rad/s},\ \zeta=0.158,\ M_p=60.5\%,\ t_s=\frac{4}{0.5}=8\ \text{s}

Step 2: Desired poles

From Mp=5%M_p=5\%:

ζ=−ln⁡0.05π2+(ln⁡0.05)2=2.9964.341=0.690\zeta=\frac{-\ln0.05}{\sqrt{\pi^2+(\ln0.05)^2}}=\frac{2.996}{4.341}=0.690

From ts=2t_s=2 s (2% criterion):

ζωn=42=2⇒ωn=20.690=2.898 rad/s,ωd=2.8981−0.6902=2.097\zeta\omega_n=\frac{4}{2}=2\Rightarrow\omega_n=\frac{2}{0.690}=2.898\ \text{rad/s},\quad\omega_d=2.898\sqrt{1-0.690^2}=2.097 sd=−2±j2.097s_d=-2\pm j2.097

Step 3: Angle deficiency

∠sd=180∘−tan⁡−12.0972=133.6∘∠(sd+1)=180∘−tan⁡−12.0971=115.5∘∠G(sd)=−249.1∘⇒ϕ=249.1∘−180∘=69.1∘\begin{aligned} \angle s_d&=180^\circ-\tan^{-1}\frac{2.097}{2}=133.6^\circ\\ \angle(s_d+1)&=180^\circ-\tan^{-1}\frac{2.097}{1}=115.5^\circ\\ \angle G(s_d)&=-249.1^\circ\Rightarrow\phi=249.1^\circ-180^\circ=69.1^\circ \end{aligned}

Step 4: Zero and pole

Place the zero at s=−1s=-1 to cancel the plant pole. Then the open loop becomes 10Kcs(s+p)\dfrac{10K_c}{s(s+p)} and the closed loop s2+ps+10Kc=0s^2+ps+10K_c=0. Matching with s2+2ζωns+ωn2=s2+4s+8.4s^2+2\zeta\omega_ns+\omega_n^2=s^2+4s+8.4:

p=2ζωn=4p=2\zeta\omega_n=4

Angle check: zero adds +115.5∘+115.5^\circ, pole at −4-4 adds −tan⁡−12.0972=−46.4∘-\tan^{-1}\frac{2.097}{2}=-46.4^\circ, net +69.1∘+69.1^\circ ✓.

Step 5: Gain

10Kc=ωn2=8.4⇒Kc=0.8410K_c=\omega_n^2=8.4\Rightarrow K_c=0.84

(same as magnitude condition Kc=∣sd∣∣sd+4∣/10=2.898×2.898/10K_c=|s_d||s_d+4|/10=2.898\times2.898/10).

Gc(s)=0.84 s+1s+4=0.21 1+s1+0.25sG_c(s)=0.84\,\frac{s+1}{s+4}=0.21\,\frac{1+s}{1+0.25s}

Step 6: Verification

C(s)R(s)=8.4s2+4s+8.4\frac{C(s)}{R(s)}=\frac{8.4}{s^2+4s+8.4}

Poles −2±j2.10-2\pm j2.10; ζ=428.4=0.690\zeta=\frac{4}{2\sqrt{8.4}}=0.690, Mp=5.0%M_p=5.0\%, ts=42=2t_s=\frac{4}{2}=2 s ✓.

          jw
  sd x    | j2.1
      \   |
 --x---o--x------ sigma
  -4  -1  0
 (zero at -1 cancels plant pole)
ItemBeforeAfter
ζ\zeta0.1580.69
ωn\omega_n3.16 rad/s2.90 rad/s
MpM_p60.5%5%
tst_s8 s2 s
KvK_v102.1

If a higher KvK_v is needed, a lag section can be added in cascade (lag–lead).

Answer: Gc(s)=0.84 s+1s+4G_c(s)=0.84\,\dfrac{s+1}{s+4}; closed loop 8.4s2+4s+8.4\dfrac{8.4}{s^2+4s+8.4} with ts = 2 s and Mp = 5%.

  • 2076 Asoj · 4 marks

What is derivative controller? How and why it can be useful?

Answer

A derivative controller gives an output proportional to the rate of change of the error:

u(t)=Kdde(t)dt,Gc(s)=Kdsu(t)=K_d\frac{de(t)}{dt},\qquad G_c(s)=K_ds

In practice it is used with proportional action as a PD controller, Gc(s)=Kp+KdsG_c(s)=K_p+K_ds.

How it works

 error e(t)        D output
   ^  /\              ^ +
   | /  \             |--+  (rising error)
   |/    \______      |  |   +-- (falling)
   +----------->      +--+---+--->
  • When the error is rising, D output is positive and pushes harder.
  • When the error is falling quickly (output approaching the set point fast), D output is negative and reduces the drive before overshoot occurs.
  • For a constant error, D output is zero.

So it acts on the future trend of the error: it is an anticipatory controller.

Why it is useful

  1. Adds damping: for a second-order plant with PD control, ζ′=ζ+Kdωn2\zeta'=\zeta+\frac{K_d\omega_n}{2} (with Kp=1K_p=1), so overshoot falls.
  2. Faster settling: ts=4/(ζ′ωn)t_s=4/(\zeta'\omega_n) decreases.
  3. Improves stability: the added zero gives phase lead and raises the phase margin, like a lead compensator.
  4. Quick response to sudden disturbances, because it reacts to the rate of change rather than waiting for the error to grow.

Limitations

  • Has no effect on steady-state error (zero output for constant error), so it is never used alone.
  • Amplifies high-frequency noise (∣Gc∣=Kdω|G_c|=K_d\omega).

Example: in a DC motor position servo, PD control lets the shaft reach the target angle quickly with little overshoot; P alone with high gain would oscillate.

  • 2075 Asoj · 12 marks

The open loop transfer function of a system is given by G(s) = 1/[s(s+1)(0.5s+1)]. Compensate the system such that Kv = 5 sec⁻¹ and phase margin is at least 40° and the gain margin is at least 10 dB with a lag compensator.

Answer

A lag compensator Gc(s)=KcβTs+1βTs+1G_c(s)=K_c\beta\dfrac{Ts+1}{\beta Ts+1} (β>1\beta>1) is used: it supplies the low-frequency gain for KvK_v and attenuates by 20log⁡β20\log\beta at higher frequencies, moving the gain crossover down to where the plant phase gives the required margins.

Step 1: Gain for KvK_v

Kv=lim⁡s→0s Gc(s)G(s)=Kcβ=5K_v=\lim_{s\to0}s\,G_c(s)G(s)=K_c\beta=5

Let K=Kcβ=5K=K_c\beta=5:

G1(s)=5s(s+1)(0.5s+1)G_1(s)=\frac{5}{s(s+1)(0.5s+1)}

Step 2: Bode plot of G1G_1

Corners at 1 and 2 rad/s; slopes −20-20, −40-40, −60-60 dB/dec; 14 dB at ω=1\omega=1 on the first asymptote.

ϕ=−90∘−tan⁡−1ω−tan⁡−10.5ω\phi=-90^\circ-\tan^{-1}\omega-\tan^{-1}0.5\omega
ω\omega0.010.10.20.511.412
$G_1$ (dB)54.033.927.818.810.0
ϕ\phi−90.9∘-90.9^\circ−98.6∘-98.6^\circ−107∘-107^\circ−130.6∘-130.6^\circ−161.6∘-161.6^\circ−180∘-180^\circ−198.4∘-198.4^\circ

ωgc=1.80\omega_{gc}=1.80 rad/s, PM =−13∘=-13^\circ; ωpc=1.414\omega_{pc}=1.414 rad/s, GM =−4.4=-4.4 dB: unstable.

Step 3: New gain crossover

Required phase =−180∘+40∘+12∘=−128∘=-180^\circ+40^\circ+12^\circ=-128^\circ (12° allowed for the lag network). Solving gives ω≈0.47\omega\approx0.47 rad/s; it is rounded to ωgc′=0.5\omega_{gc}'=0.5 rad/s (ϕ=−130.6∘\phi=-130.6^\circ there).

Step 4: Zero and β

  • Corner of the zero well below ωgc′\omega_{gc}': 1T=0.1\frac1T=0.1 rad/s (one-fifth of 0.5) ⇒T=10\Rightarrow T=10 s.
  • Attenuation needed at 0.5 rad/s is about 19–20 dB (18.8 dB from the table), so take 20log⁡β=2020\log\beta=20 dB: β=10\beta=10.
  • Pole at 1βT=0.01\frac{1}{\beta T}=0.01 rad/s.
Kc=Kβ=510=0.5K_c=\frac{K}{\beta}=\frac{5}{10}=0.5 Gc(s)=0.5×10 10s+1100s+1=5 10s+1100s+1=0.5 s+0.1s+0.01G_c(s)=0.5\times10\,\frac{10s+1}{100s+1}=5\,\frac{10s+1}{100s+1}=0.5\,\frac{s+0.1}{s+0.01}

Step 5: Compensated system

Gc(s)G(s)=5(10s+1)s(100s+1)(s+1)(0.5s+1)G_c(s)G(s)=\frac{5(10s+1)}{s(100s+1)(s+1)(0.5s+1)}
ω\omega0.010.10.20.4511.32
$G_cG$ (dB)51.016.98.70
ϕ\phi−130∘-130^\circ−138∘-138^\circ−131∘-131^\circ−138.4∘-138.4^\circ−167∘-167^\circ−180∘-180^\circ
PM=180∘−138.4∘=41.6∘ (≥40∘),GM=14.3 dB (≥10 dB)PM=180^\circ-138.4^\circ=41.6^\circ\ (\ge40^\circ),\qquad GM=14.3\ \text{dB}\ (\ge10\ \text{dB})
 dB 54 |\___ uncompensated
       |    \___
     0 |--x-----\----------- w
       | 0.45    \ 1.8
       | compensated: 20 dB lower
       | above 0.1 rad/s
ItemBeforeAfter
KvK_v55
ωgc\omega_{gc}1.80 rad/s0.45 rad/s
PM−13∘-13^\circ41.6∘41.6^\circ
GM−4.4-4.4 dB14.3 dB

Closed-loop poles: −0.286±j0.520-0.286\pm j0.520, −0.123-0.123, −2.32-2.32 (all in LHP).

Answer: Gc(s)=5 10s+1100s+1G_c(s)=5\,\dfrac{10s+1}{100s+1} (β = 10, T = 10 s): Kv = 5 s⁻¹, PM ≈ 42°, GM ≈ 14 dB.

  • 2074 Chaitra · 4 marks

Discuss working of PI controller.

Answer

A PI controller combines proportional action with integral action. Its output is

u(t)=Kpe(t)+Ki∫0te(τ) dτ=Kp[e(t)+1Ti∫0te dτ]u(t)=K_pe(t)+K_i\int_0^te(\tau)\,d\tau=K_p\left[e(t)+\frac{1}{T_i}\int_0^te\,d\tau\right] Gc(s)=Kp+Kis=Kp(s+1Ti)s,Ti=KpKiG_c(s)=K_p+\frac{K_i}{s}=\frac{K_p\left(s+\frac{1}{T_i}\right)}{s},\quad T_i=\frac{K_p}{K_i}
          +--->[  Kp  ]---+
 e(t) ----|               (+)---> u(t)
          +--->[ Ki/s ]---+

Working

  1. P part reacts at once to the present error, giving fast correction.
  2. I part adds up (integrates) the error over time. As long as any error remains, the integral keeps growing and keeps changing the output.
  3. Therefore the system can settle only when the error is exactly zero. A constant controller output is then held by the integrator alone; this removes the offset that a pure P controller leaves.

For a step error EE, the output is u(t)=KpE+KiEtu(t)=K_pE+K_iEt: an instant jump followed by a ramp.

Effects

  • Adds a pole at the origin (type +1) and a zero at −1/Ti-1/T_i.
  • Steady-state error is eliminated for a step (type 0 plant) or ramp (type 1 plant).
  • Adds phase lag, so overshoot increases and relative stability decreases; the zero is placed near the origin to keep this small.
  • Acts like a low-pass filter (similar to a lag compensator); noise is not amplified.
  • Problem: integral wind-up when the actuator saturates.

Example: in a water-level control with P only, the level settles 2 cm below the set point; adding integral action slowly raises the valve opening until the level reaches the set point exactly.

  • 2074 Chaitra · 4 marks

Discuss the purpose of lead and lag compensators.

Answer

A compensator is an additional network placed in a control loop to change the system's response so that it meets the specifications that adjusting gain alone cannot meet. Lead and lag compensators are the most common.

Lead compensator

Gc(s)=1+Ts1+αTs,α<1 (zero closer to origin than pole)G_c(s)=\frac{1+Ts}{1+\alpha Ts},\quad\alpha<1\ \text{(zero closer to origin than pole)}

Purpose:

  • Supplies positive phase (lead), up to ϕm=sin⁡−11−α1+α\phi_m=\sin^{-1}\frac{1-\alpha}{1+\alpha} at ωm=1/(Tα)\omega_m=1/(T\sqrt\alpha).
  • Increases phase margin and damping, so overshoot falls.
  • Increases gain crossover frequency and bandwidth, so the response is faster (smaller rise and settling time).
  • Reshapes the root locus to the left (improves transient response and stability).

Lag compensator

Gc(s)=1+Ts1+βTs,β>1 (pole closer to origin)G_c(s)=\frac{1+Ts}{1+\beta Ts},\quad\beta>1\ \text{(pole closer to origin)}

Purpose:

  • Raises low-frequency gain by β, improving steady-state accuracy (KpK_p, KvK_v, KaK_a larger) with little change to the transient.
  • Attenuates high frequencies by 20log⁡β20\log\beta, which lowers the gain crossover so that the PM increases (when the plant has enough phase at low frequency).
  • Reduces bandwidth, so it filters noise but slows the response.

Summary

AspectLeadLag
Main aimbetter transientbetter steady state
Phaseadds leadadds lag
Bandwidthincreasesdecreases
Filter behaviourhigh-passlow-pass
Analogous controllerPDPI

A lag–lead compensator combines both when both transient and steady-state improvements are needed.

  • 2073 Shrawan · 4 marks

How can a controller with transfer function Gc(s) = (1 + aTs)/(1 + Ts) be used as lead or lag compensator? Explain.

Answer

Gc(s)=1+aTs1+TsG_c(s)=\frac{1+aTs}{1+Ts}

It has a zero at s=−1aTs=-\frac{1}{aT} and a pole at s=−1Ts=-\frac{1}{T}. Whether it is lead or lag depends only on aa.

Phase of the network

∠Gc(jω)=tan⁡−1(aωT)−tan⁡−1(ωT)\angle G_c(j\omega)=\tan^{-1}(a\omega T)-\tan^{-1}(\omega T)

Case 1: a>1a>1 — lead compensator

  • aωT>ωTa\omega T>\omega T, so the phase is positive at all frequencies.
  • Zero (−1/aT-1/aT) is nearer the origin than the pole (−1/T-1/T).
  • Maximum lead at ωm=1Ta\omega_m=\frac{1}{T\sqrt a}: sin⁡ϕm=a−1a+1\sin\phi_m=\frac{a-1}{a+1}; high-frequency gain =a=a (20log⁡a20\log a dB).
  • Use: improve phase margin, damping and speed (transient response).

Case 2: a<1a<1 — lag compensator

  • aωT<ωTa\omega T<\omega T, so the phase is negative.
  • Pole (−1/T-1/T) is nearer the origin than the zero (−1/aT-1/aT).
  • High-frequency gain =a<1=a<1: attenuation of 20log⁡(1/a)20\log(1/a) dB; DC gain 1.
  • Use: improve steady-state accuracy (with the loop gain raised by 1/a1/a), or lower the crossover to gain PM.
 Lead (a>1):   -x------o-----+  sigma
               -1/T  -1/aT   0
 Lag  (a<1):   -o------x-----+  sigma
              -1/aT  -1/T    0

(a=1a=1 gives Gc=1G_c=1, no compensation.)

a>1a>1a<1a<1
Typeleadlag
Phasepositivenegative
Gain at high ω\omegaaa (boost)aa (cut)

Example: Gc=1+0.5s1+0.1sG_c=\frac{1+0.5s}{1+0.1s} (a=5a=5, T=0.1T=0.1): lead, maximum lead sin⁡−146=41.8∘\sin^{-1}\frac{4}{6}=41.8^\circ at ωm=10.15=4.47\omega_m=\frac{1}{0.1\sqrt5}=4.47 rad/s. Gc=1+2s1+20sG_c=\frac{1+2s}{1+20s} (a=0.1a=0.1, T=20T=20): lag.

  • 2072 Chaitra · 4 marks

If desired damping ratio is '1', which controller do you suggest? Explain.

Answer

To obtain a damping ratio of 1 (critical damping: fastest response with no overshoot), a controller that adds damping is needed: a PD (proportional–derivative) controller, or equivalently derivative (tachometer) output feedback. Gain adjustment alone cannot usually do it without making the response very slow.

Why PD

For a typical underdamped plant G(s)=ωn2s(s+2ζωn)G(s)=\dfrac{\omega_n^2}{s(s+2\zeta\omega_n)}, PD control Gc(s)=1+TdsG_c(s)=1+T_ds gives

C(s)R(s)=ωn2(1+Tds)s2+(2ζωn+Tdωn2)s+ωn2\frac{C(s)}{R(s)}=\frac{\omega_n^2(1+T_ds)}{s^2+(2\zeta\omega_n+T_d\omega_n^2)s+\omega_n^2} ζ′=ζ+Tdωn2\zeta'=\zeta+\frac{T_d\omega_n}{2}

The derivative term raises ζ while ωn\omega_n stays the same. For ζ′=1\zeta'=1:

Td=2(1−ζ)ωnT_d=\frac{2(1-\zeta)}{\omega_n}

With derivative output feedback (KtsK_ts in the minor loop), ζ′=ζ+Ktωn2\zeta'=\zeta+\frac{K_t\omega_n}{2}, so Kt=2(1−ζ)ωnK_t=\frac{2(1-\zeta)}{\omega_n}; this has no added zero, so the response is truly non-overshooting.

Why not others

ControllerEffect on ζ
P (lower gain)ζ rises but ωn\omega_n and speed fall; steady-state error rises
PIadds lag: ζ falls
PD / rate feedbackζ rises, speed kept

Example: G(s)=16s(s+2)G(s)=\dfrac{16}{s(s+2)}: ωn=4\omega_n=4, ζ=0.25\zeta=0.25. Rate feedback with Kt=2(1−0.25)4=0.375K_t=\dfrac{2(1-0.25)}{4}=0.375 gives s2+(2+6)s+16=(s+4)2s^2+(2+6)s+16=(s+4)^2, i.e. ζ=1\zeta=1, a critically damped response with ts≈5.8ωn=1.46t_s\approx\frac{5.8}{\omega_n}=1.46 s and no overshoot.

  • 2072 Chaitra · 4 marks

Compare the Lag and Lead compensator applications in control system.

Answer

A lead compensator Gc(s)=s+zs+pG_c(s) = \dfrac{s+z}{s+p} with p>zp > z adds positive phase near the gain crossover and mainly improves the transient response. A lag compensator with p<zp < z, both close to the origin, adds low-frequency gain and mainly improves steady-state accuracy.

PointLead compensatorLag compensator
Pole–zero positionZero nearer the origin than the polePole nearer the origin than the zero
Phase addedPositive (phase lead)Negative (phase lag), kept small
Main purposeBetter transient response and stabilityBetter steady-state error (KpK_p, KvK_v, KaK_a)
Effect on bandwidthIncreases bandwidth, faster responseReduces bandwidth, slower response
Gain crossover frequencyMoves to a higher valueMoves to a lower value
Phase marginRaised directly by the phase leadRaised by lowering the crossover frequency
Effect on noisePasses more high-frequency noiseFilters (attenuates) high-frequency noise
Root locusPulls the locus to the left (more stable)Locus nearly unchanged near dominant poles
Network typeActs like a high-pass filter (like PD)Acts like a low-pass filter (like PI)

Where each is used:

  • Lead: when the system is too slow or too oscillatory, e.g. reducing overshoot and settling time of a position servo, or adding phase margin to a type-1 or type-2 plant.
  • Lag: when the transient response is already acceptable but the steady-state error is too large, e.g. raising KvK_v of a tracking antenna or a temperature controller by about 10 times without changing the overshoot.
  • When both transient and steady-state performance must improve, a lag–lead compensator is used.
  • 2072 Chaitra · 12 marks

Design a suitable compensator for a unity feedback system with its feed forward transfer function as G(s) = 4/[s(s+2)] such that its maximum percent overshoot is 16.3% and settling time 2 sec for its step response. Also velocity error constant should not be less than 2 per sec.

Answer

The settling time must be halved while keeping ζ\zeta, so the dominant poles must move left; a lead compensator designed by the root locus method is used.

Step 1: Desired closed-loop poles

Peak overshoot:

Mp=e−πζ/1−ζ2=0.163ζ=−ln⁡0.163π2+(ln⁡0.163)2=1.8143.628=0.5\begin{aligned} M_p &= e^{-\pi\zeta/\sqrt{1-\zeta^2}} = 0.163 \\ \zeta &= \frac{-\ln 0.163}{\sqrt{\pi^2 + (\ln 0.163)^2}} = \frac{1.814}{3.628} = 0.5 \end{aligned}

Settling time (2% criterion):

ts=4ζωn=2  ⇒  ζωn=2,ωn=20.5=4 rad/st_s = \frac{4}{\zeta\omega_n} = 2 \;\Rightarrow\; \zeta\omega_n = 2,\quad \omega_n = \frac{2}{0.5} = 4\ \text{rad/s} sd=−ζωn±jωn1−ζ2=−2±j3.464s_d = -\zeta\omega_n \pm j\omega_n\sqrt{1-\zeta^2} = -2 \pm j3.464

Step 2: Uncompensated system

Closed-loop characteristic equation: s2+2s+4=0s^2 + 2s + 4 = 0, so ωn=2\omega_n = 2, ζ=0.5\zeta = 0.5.

  • Mp=16.3%M_p = 16.3\% (already correct), ts=40.5×2=4t_s = \dfrac{4}{0.5\times2} = 4 s (too slow).
  • Kv=lim⁡s→0s 4s(s+2)=2 s−1K_v = \lim_{s\to0} s\,\dfrac{4}{s(s+2)} = 2\ \text{s}^{-1}.

So the poles must move further left along the same ζ=0.5\zeta = 0.5 line, which needs a lead compensator.

Step 3: Angle deficiency

At sd=−2+j3.464s_d = -2 + j3.464:

∠sd=180∘−tan⁡−13.4642=120∘∠(sd+2)=90∘∠G(sd)=−(120∘+90∘)=−210∘ϕ=−180∘−(−210∘)=30∘\begin{aligned} \angle s_d &= 180^\circ - \tan^{-1}\frac{3.464}{2} = 120^\circ \\ \angle (s_d + 2) &= 90^\circ \\ \angle G(s_d) &= -(120^\circ + 90^\circ) = -210^\circ \\ \phi &= -180^\circ - (-210^\circ) = 30^\circ \end{aligned}

The compensator must add +30∘+30^\circ at sds_d.

Step 4: Place the zero and pole

Let Gc(s)=Kcs+zs+pG_c(s) = K_c\dfrac{s+z}{s+p}. Put the zero at s=−2s = -2 to cancel the plant pole at −2-2, so the zero gives 90∘90^\circ. The pole must then give 90∘−30∘=60∘90^\circ - 30^\circ = 60^\circ:

tan⁡60∘=3.464p−2  ⇒  p−2=2  ⇒  p=4\tan 60^\circ = \frac{3.464}{p - 2} \;\Rightarrow\; p - 2 = 2 \;\Rightarrow\; p = 4

Step 5: Gain from the magnitude condition

The compensated open loop is 4Kcs(s+4)\dfrac{4K_c}{s(s+4)}:

4Kc=∣sd∣ ∣sd+4∣=∣−2+j3.464∣ ∣2+j3.464∣=4×4=16Kc=4\begin{aligned} 4K_c &= |s_d|\,|s_d+4| = |-2+j3.464|\,|2+j3.464| = 4 \times 4 = 16 \\ K_c &= 4 \end{aligned}

Check: characteristic equation s2+4s+16=0s^2 + 4s + 16 = 0 gives ωn=4\omega_n = 4, ζ=4/(2×4)=0.5\zeta = 4/(2\times4) = 0.5. Correct.

Gc(s)=4(s+2)s+4=2(0.5s+1)0.25s+1G_c(s) = \frac{4(s+2)}{s+4} = \frac{2(0.5s+1)}{0.25s+1}

Step 6: Check velocity error constant

Kv=lim⁡s→0s Gc(s)G(s)=lim⁡s→0s⋅16s(s+4)=4 s−1  ≥  2 s−1K_v = \lim_{s\to0} s\,G_c(s)G(s) = \lim_{s\to0} s\cdot\frac{16}{s(s+4)} = 4\ \text{s}^{-1} \;\geq\; 2\ \text{s}^{-1}

Result

QuantityUncompensatedCompensated
Closed-loop poles−1±j1.732-1 \pm j1.732−2±j3.464-2 \pm j3.464
ζ\zeta, ωn\omega_n0.5, 2 rad/s0.5, 4 rad/s
MpM_p16.3%16.3%
tst_s (2%)4 s2 s
KvK_v2 s−1^{-1}4 s−1^{-1}

Answer: Gc(s)=4(s+2)s+4G_c(s) = \dfrac{4(s+2)}{s+4} (lead compensator, zero at −2-2, pole at −4-4, α=0.5\alpha = 0.5). All three specifications are met.

  • 2071 Chaitra · 12 marks

Design a suitable compensator for a unity feedback system with open loop transfer function G(s) = 4/[s(s+2)] such that the settling time will become 2 seconds without change in overshoot and velocity error constant will be 2 s⁻¹.

Answer

"Without change in overshoot" means the damping ratio stays the same, while the settling time is halved. The dominant poles must move left along the same ζ\zeta line, so a lead compensator designed by root locus is used.

Step 1: Present (uncompensated) performance

Closed loop: 4s2+2s+4\dfrac{4}{s^2 + 2s + 4}, so ωn=2\omega_n = 2 rad/s and ζ=22×2=0.5\zeta = \dfrac{2}{2\times2} = 0.5.

  • Mp=e−π(0.5)/0.75=16.3%M_p = e^{-\pi(0.5)/\sqrt{0.75}} = 16.3\%
  • ts=4ζωn=41=4t_s = \dfrac{4}{\zeta\omega_n} = \dfrac{4}{1} = 4 s (2% criterion)
  • Kv=lim⁡s→0s4s(s+2)=2 s−1K_v = \lim_{s\to0} s\dfrac{4}{s(s+2)} = 2\ \text{s}^{-1}

Step 2: Desired dominant poles

Keep ζ=0.5\zeta = 0.5 and make ts=2t_s = 2 s:

ζωn=42=2,ωn=4 rad/s,sd=−2±j41−0.25=−2±j3.464\zeta\omega_n = \frac{4}{2} = 2,\qquad \omega_n = 4\ \text{rad/s},\qquad s_d = -2 \pm j4\sqrt{1-0.25} = -2 \pm j3.464

Step 3: Angle deficiency

∠G(sd)=−[∠sd+∠(sd+2)]=−(120∘+90∘)=−210∘,ϕ=30∘\angle G(s_d) = -\left[\angle s_d + \angle(s_d+2)\right] = -(120^\circ + 90^\circ) = -210^\circ, \qquad \phi = 30^\circ

Step 4: Place the zero and pole

Let Gc(s)=Kcs+zs+pG_c(s) = K_c\dfrac{s+z}{s+p}. Put the zero at s=−2s = -2 to cancel the plant pole at −2-2, so the zero gives 90∘90^\circ. The pole must then give 90∘−30∘=60∘90^\circ - 30^\circ = 60^\circ:

tan⁡60∘=3.464p−2  ⇒  p−2=2  ⇒  p=4\tan 60^\circ = \frac{3.464}{p - 2} \;\Rightarrow\; p - 2 = 2 \;\Rightarrow\; p = 4

Step 5: Gain from the magnitude condition

The compensated open loop is 4Kcs(s+4)\dfrac{4K_c}{s(s+4)}:

4Kc=∣sd∣ ∣sd+4∣=∣−2+j3.464∣ ∣2+j3.464∣=4×4=16Kc=4\begin{aligned} 4K_c &= |s_d|\,|s_d+4| = |-2+j3.464|\,|2+j3.464| = 4 \times 4 = 16 \\ K_c &= 4 \end{aligned}

Check: characteristic equation s2+4s+16=0s^2 + 4s + 16 = 0 gives ωn=4\omega_n = 4, ζ=4/(2×4)=0.5\zeta = 4/(2\times4) = 0.5. Correct.

Gc(s)=4(s+2)s+4=2(0.5s+1)0.25s+1G_c(s) = \frac{4(s+2)}{s+4} = \frac{2(0.5s+1)}{0.25s+1}

Step 6: Velocity error constant

Kv=lim⁡s→0s⋅16s(s+4)=4 s−1K_v = \lim_{s\to0} s\cdot\frac{16}{s(s+4)} = 4\ \text{s}^{-1}

This is not less than the required 2 s−12\ \text{s}^{-1}, so steady-state accuracy is also better.

(If KvK_v must be exactly 2, the zero can be moved to −0.8-0.8 with the pole at −2.67-2.67 and Kc=3.33K_c = 3.33; the same angle condition holds and Kv=2K_v = 2, but a third closed-loop pole appears at −0.67-0.67 and slows the response, so the design above is preferred.)

Result

QuantityBeforeAfter
ζ\zeta0.50.5
MpM_p16.3%16.3% (unchanged)
ωn\omega_n2 rad/s4 rad/s
tst_s4 s2 s
KvK_v2 s−1^{-1}4 s−1^{-1}

Answer: Gc(s)=4(s+2)s+4G_c(s) = \dfrac{4(s+2)}{s+4}; closed-loop poles −2±j3.464-2 \pm j3.464, overshoot 16.3%, settling time 2 s, Kv=4 s−1K_v = 4\ \text{s}^{-1}.

  • 2071 Chaitra · 4 marks

For a compensator transfer function given by Gc(s) = (s + τ)/(s + aτ), give the condition of lead compensator. For the given value of 'a' what is the frequency that leads to maximum phase angle lead?

Answer

For Gc(s)=s+τs+aτG_c(s) = \dfrac{s+\tau}{s+a\tau} the zero is at s=−τs = -\tau and the pole is at s=−aτs = -a\tau.

Condition for a lead compensator

The network gives phase lead when the zero lies closer to the origin than the pole:

aτ>τ⇒a>1a\tau > \tau \quad \Rightarrow \quad a > 1

Then the phase angle is

ϕ(ω)=tan⁡−1ωτ−tan⁡−1ωaτ>0for all ω>0\phi(\omega) = \tan^{-1}\frac{\omega}{\tau} - \tan^{-1}\frac{\omega}{a\tau} > 0 \quad \text{for all } \omega > 0

(If a<1a < 1 the pole is nearer the origin and the network becomes a lag compensator.)

Frequency of maximum phase lead

Differentiate ϕ\phi with respect to ω\omega and set it to zero:

dϕdω=1/τ1+ω2/τ2−1/(aτ)1+ω2/(aτ)2=0ττ2+ω2=aτa2τ2+ω2a2τ2+ω2=aτ2+aω2ω2(a−1)=aτ2(a−1)ωm=τa\begin{aligned} \frac{d\phi}{d\omega} &= \frac{1/\tau}{1+\omega^2/\tau^2} - \frac{1/(a\tau)}{1+\omega^2/(a\tau)^2} = 0 \\ \frac{\tau}{\tau^2+\omega^2} &= \frac{a\tau}{a^2\tau^2+\omega^2} \\ a^2\tau^2 + \omega^2 &= a\tau^2 + a\omega^2 \\ \omega^2 (a-1) &= a\tau^2 (a-1) \\ \omega_m &= \tau\sqrt{a} \end{aligned}

So the maximum phase lead occurs at ωm=τ×aτ=τa\omega_m = \sqrt{\tau \times a\tau} = \tau\sqrt{a}, the geometric mean of the zero and pole frequencies (the midpoint on the log-frequency scale of the Bode plot).

The maximum lead is

ϕm=sin⁡−1a−1a+1\phi_m = \sin^{-1}\frac{a-1}{a+1}

Example: τ=2\tau = 2, a=4a = 4 gives ωm=24=4\omega_m = 2\sqrt{4} = 4 rad/s and ϕm=sin⁡−1(3/5)=36.9∘\phi_m = \sin^{-1}(3/5) = 36.9^\circ.

  • 2070 Chaitra (old course) · 16 marks

A system has open loop transfer function Gf(s) = 4/[s(s+2)]. It is desired to design a compensator so that the static velocity error constant Kv is 20 sec⁻¹, phase margin is at least 50° and gain margin is at least 10 dB.

Answer

The uncompensated system has a low phase margin once the gain is raised for Kv=20K_v = 20, so a phase-lead compensator is designed using the Bode plot.

Take Gc(s)=Kc α Ts+1αTs+1G_c(s) = K_c\,\alpha\,\dfrac{Ts+1}{\alpha Ts+1}, 0<α<10 < \alpha < 1.

Step 1: Gain for KvK_v

Kv=lim⁡s→0s Gc(s)Gf(s)=Kcα⋅42=2Kcα=20  ⇒  K=Kcα=10K_v = \lim_{s\to0} s\,G_c(s)G_f(s) = K_c\alpha\cdot\frac{4}{2} = 2K_c\alpha = 20 \;\Rightarrow\; K = K_c\alpha = 10

Gain-adjusted system: G1(s)=40s(s+2)=20s(0.5s+1)G_1(s) = \dfrac{40}{s(s+2)} = \dfrac{20}{s(0.5s+1)}.

Step 2: Phase margin of G1G_1

Gain crossover where ∣G1(jω)∣=1|G_1(j\omega)| = 1:

40ωω2+4=1  ⇒  ω4+4ω2−1600=0ω2=38.05,ωgc=6.17 rad/s∠G1=−90∘−tan⁡−1(6.17/2)=−162.0∘PM=180∘−162.0∘=18.0∘\begin{aligned} \frac{40}{\omega\sqrt{\omega^2+4}} &= 1 \;\Rightarrow\; \omega^4 + 4\omega^2 - 1600 = 0 \\ \omega^2 &= 38.05,\quad \omega_{gc} = 6.17\ \text{rad/s} \\ \angle G_1 &= -90^\circ - \tan^{-1}(6.17/2) = -162.0^\circ \\ PM &= 180^\circ - 162.0^\circ = 18.0^\circ \end{aligned}

The phase never reaches −180∘-180^\circ (second-order type-1 system), so GM=∞GM = \infty; the gain margin requirement of 10 dB is automatically met.

Step 3: Phase lead needed

Add about 8∘8^\circ for the extra lag caused by the shift of the crossover frequency:

ϕm=50∘−18∘+8∘=40∘\phi_m = 50^\circ - 18^\circ + 8^\circ = 40^\circ α=1−sin⁡40∘1+sin⁡40∘=0.3571.643=0.217\alpha = \frac{1-\sin 40^\circ}{1+\sin 40^\circ} = \frac{0.357}{1.643} = 0.217

Step 4: New gain crossover frequency

The compensator adds 10log⁡(1/α)=6.6310\log(1/\alpha) = 6.63 dB at ωm\omega_m. So ωm\omega_m is where ∣G1∣=−6.63|G_1| = -6.63 dB, i.e. ∣G1∣=α=0.466|G_1| = \sqrt{\alpha} = 0.466:

40ωω2+4=0.466  ⇒  ωm=9.16 rad/s\frac{40}{\omega\sqrt{\omega^2+4}} = 0.466 \;\Rightarrow\; \omega_m = 9.16\ \text{rad/s}

Step 5: Corner frequencies

T=1ωmα=19.16×0.466=0.234 szero: 1T=4.27 rad/s,pole: 1αT=19.66 rad/sKc=10α=46.1\begin{aligned} T &= \frac{1}{\omega_m\sqrt{\alpha}} = \frac{1}{9.16\times0.466} = 0.234\ \text{s} \\ \text{zero: } \frac{1}{T} &= 4.27\ \text{rad/s}, \qquad \text{pole: } \frac{1}{\alpha T} = 19.66\ \text{rad/s} \\ K_c &= \frac{10}{\alpha} = 46.1 \end{aligned}

Step 6: Compensator and check

Gc(s)=46.1 s+4.27s+19.66=10 0.234s+10.0508s+1G_c(s) = 46.1\,\frac{s+4.27}{s+19.66} = 10\,\frac{0.234s+1}{0.0508s+1} Gc(s)Gf(s)=40(0.234s+1)s(s+2)(0.0508s+1)G_c(s)G_f(s) = \frac{40(0.234s+1)}{s(s+2)(0.0508s+1)}
QuantityUncompensated (K=10K=10)Compensated
Gain crossover6.17 rad/s9.15 rad/s
Phase margin18.0°52.3°
Gain margin∞\infty∞\infty
KvK_v20 s−1^{-1}20 s−1^{-1}

Answer: Gc(s)=46.1s+4.27s+19.66G_c(s) = 46.1\dfrac{s+4.27}{s+19.66} gives Kv=20 s−1K_v = 20\ \text{s}^{-1}, PM=52.3∘ (≥50∘)PM = 52.3^\circ\ (\ge 50^\circ) and GM=∞ (≥10GM = \infty\ (\ge 10 dB)).

  • 2070 Chaitra · 4 marks

Discuss the application of a PI controller with suitable example.

Answer

A PI controller produces a control signal proportional to both the error and the integral of the error:

u(t)=Kpe(t)+Ki∫0te(t) dt,Gc(s)=Kp+Kis=Kp(s+Ki/Kp)su(t) = K_p e(t) + K_i \int_0^t e(t)\,dt, \qquad G_c(s) = K_p + \frac{K_i}{s} = \frac{K_p\left(s + K_i/K_p\right)}{s}

It adds a pole at the origin and a zero at s=−Ki/Kps = -K_i/K_p.

Why it is applied

  • The pole at the origin raises the system type by one. A type-0 plant gets zero steady-state error to a step; a type-1 plant gets zero error to a ramp.
  • It removes the steady-state offset that a pure P controller always leaves.
  • The zero, placed close to the origin, limits the extra phase lag, so stability is not badly affected.
  • It acts like a lag compensator (low-pass), so it also filters high-frequency noise.
  • Drawbacks: slightly slower response, possibly more overshoot, and integral wind-up when the actuator saturates.

Example: speed control of a DC motor

Let the motor (speed output) be G(s)=10s+2G(s) = \dfrac{10}{s+2} with unity feedback (type 0).

  • With P control Kp=4K_p = 4: KP=lim⁡s→040s+2=20K_P = \lim_{s\to0} \dfrac{40}{s+2} = 20, so the step error ess=11+20=0.048e_{ss} = \dfrac{1}{1+20} = 0.048, i.e. a 4.8% speed offset that grows when load torque is applied.
  • With PI control Gc=4+8sG_c = 4 + \dfrac{8}{s}:
GcG=4(s+2)s⋅10s+2=40sG_c G = \frac{4(s+2)}{s}\cdot\frac{10}{s+2} = \frac{40}{s}

The system becomes type 1, KP=∞K_P = \infty and ess=0e_{ss} = 0. The closed loop is 40s+40\dfrac{40}{s+40}, a fast first-order response with time constant 0.025 s and no offset, even under constant load disturbances.

Other common uses: liquid-level and flow control, temperature control of furnaces, pressure control in process plants and voltage regulators, wherever zero steady-state offset is essential.

  • 2070 Chaitra · 12 marks

Design series lag compensator for the unity feedback system with feedforward transfer function G(s) = K/[s(s+4)(s+80)]. The velocity error constant is 30 s⁻¹ and phase margin at least 33°.

Answer

A lag compensator Gc(s)=1+Ts1+βTsG_c(s) = \dfrac{1+Ts}{1+\beta Ts} (β>1\beta > 1) is designed with the Bode plot method: the gain is fixed by KvK_v, and the lag network lowers the high-frequency gain so the crossover moves to a frequency with enough phase margin.

Step 1: Gain for KvK_v

Kv=lim⁡s→0s Ks(s+4)(s+80)=K320=30  ⇒  K=9600K_v = \lim_{s\to0} s\,\frac{K}{s(s+4)(s+80)} = \frac{K}{320} = 30 \;\Rightarrow\; K = 9600 G(s)=9600s(s+4)(s+80)=30s(0.25s+1)(0.0125s+1)G(s) = \frac{9600}{s(s+4)(s+80)} = \frac{30}{s(0.25s+1)(0.0125s+1)}

Step 2: Uncompensated margins

  • Gain crossover: ∣G(jω)∣=1|G(j\omega)| = 1 at ωgc=10.55\omega_{gc} = 10.55 rad/s; ∠G=−90∘−tan⁡−110.554−tan⁡−110.5580=−166.7∘\angle G = -90^\circ - \tan^{-1}\frac{10.55}{4} - \tan^{-1}\frac{10.55}{80} = -166.7^\circ, so PM=13.3∘PM = 13.3^\circ.
  • Phase crossover: ωpc=4×80=17.9\omega_{pc} = \sqrt{4\times80} = 17.9 rad/s, GM=8.9GM = 8.9 dB.

The PM is far below 33∘33^\circ.

Step 3: New gain crossover frequency

Allow 5∘5^\circ for the lag of the compensator at the new crossover. Required plant phase:

∠G(jωgc′)=−180∘+33∘+5∘=−142∘\angle G(j\omega_{gc}') = -180^\circ + 33^\circ + 5^\circ = -142^\circ −90∘−tan⁡−1ω4−tan⁡−1ω80=−142∘  ⇒  ωgc′=4.56 rad/s-90^\circ - \tan^{-1}\frac{\omega}{4} - \tan^{-1}\frac{\omega}{80} = -142^\circ \;\Rightarrow\; \omega_{gc}' = 4.56\ \text{rad/s}

Step 4: Value of β\beta

At 4.56 rad/s:

∣G(j4.56)∣=96004.56×4.562+16×4.562+6400=4.33  (12.7 dB)|G(j4.56)| = \frac{9600}{4.56\times\sqrt{4.56^2+16}\times\sqrt{4.56^2+6400}} = 4.33 \;(12.7\ \text{dB})

The lag network must reduce the gain by this amount: β=4.33\beta = 4.33.

Step 5: Corner frequencies

Place the zero one decade below the new crossover:

1T=4.5610=0.456 rad/s,T=2.19 s1βT=0.4564.33=0.105 rad/s,βT=9.50 s\begin{aligned} \frac{1}{T} &= \frac{4.56}{10} = 0.456\ \text{rad/s},\quad T = 2.19\ \text{s} \\ \frac{1}{\beta T} &= \frac{0.456}{4.33} = 0.105\ \text{rad/s},\quad \beta T = 9.50\ \text{s} \end{aligned}

Step 6: Compensator and check

Gc(s)=1+2.19s1+9.50s=0.231 s+0.456s+0.105G_c(s) = \frac{1 + 2.19s}{1 + 9.50s} = 0.231\,\frac{s+0.456}{s+0.105} Gc(s)G(s)=9600(1+2.19s)s(s+4)(s+80)(1+9.50s)G_c(s)G(s) = \frac{9600(1+2.19s)}{s(s+4)(s+80)(1+9.50s)}
QuantityUncompensatedCompensated
Gain crossover10.55 rad/s4.57 rad/s
Phase margin13.3°33.6°
Gain margin8.9 dB20.9 dB
KvK_v30 s−1^{-1}30 s−1^{-1}

The lag network adds only −4.4∘-4.4^\circ at the new crossover, within the 5∘5^\circ allowed.

Answer: K=9600K = 9600, Gc(s)=1+2.19s1+9.50sG_c(s) = \dfrac{1+2.19s}{1+9.50s}; Kv=30 s−1K_v = 30\ \text{s}^{-1} and PM=33.6∘PM = 33.6^\circ.

  • 2069 Chaitra · 16 marks

Design a suitable cascade lag compensator network for the given system G(s) = 50K/[s(s+5)(s+10)] such that the requirement of velocity error constant of 30 sec⁻¹ and phase margin of ≥ 45° are met.

Answer

A cascade lag network Gc(s)=1+Ts1+βTsG_c(s) = \dfrac{1+Ts}{1+\beta Ts} (β>1\beta > 1, unity DC gain) is designed with the Bode plot. A lag network is suitable because the plant has enough phase at a lower frequency; we only need to lower the gain there.

Step 1: Gain for KvK_v

Kv=lim⁡s→0s 50Ks(s+5)(s+10)=50K50=K=30K_v = \lim_{s\to0} s\,\frac{50K}{s(s+5)(s+10)} = \frac{50K}{50} = K = 30 G(s)=1500s(s+5)(s+10)=30s(0.2s+1)(0.1s+1)G(s) = \frac{1500}{s(s+5)(s+10)} = \frac{30}{s(0.2s+1)(0.1s+1)}

Step 2: Uncompensated margins

  • ∣G(jω)∣=1|G(j\omega)| = 1 at ωgc=9.77\omega_{gc} = 9.77 rad/s; ∠G=−90∘−tan⁡−19.775−tan⁡−19.7710=−197.2∘\angle G = -90^\circ - \tan^{-1}\frac{9.77}{5} - \tan^{-1}\frac{9.77}{10} = -197.2^\circ, so PM=−17.2∘PM = -17.2^\circ.
  • ωpc=50=7.07\omega_{pc} = \sqrt{50} = 7.07 rad/s, ∣G(j7.07)∣=2|G(j7.07)| = 2, GM=−6.0GM = -6.0 dB.

The uncompensated system is unstable.

Step 3: New gain crossover frequency

Allow 10∘10^\circ for the lag of the network (needed for this large shift):

∠G(jωgc′)=−180∘+45∘+10∘=−125∘\angle G(j\omega_{gc}') = -180^\circ + 45^\circ + 10^\circ = -125^\circ −90∘−tan⁡−1ω5−tan⁡−1ω10=−125∘  ⇒  ωgc′=2.12 rad/s-90^\circ - \tan^{-1}\frac{\omega}{5} - \tan^{-1}\frac{\omega}{10} = -125^\circ \;\Rightarrow\; \omega_{gc}' = 2.12\ \text{rad/s}

Step 4: Value of β\beta

∣G(j2.12)∣=15002.12×2.122+25×2.122+100=12.72  (22.1 dB)  ⇒  β=12.72|G(j2.12)| = \frac{1500}{2.12\times\sqrt{2.12^2+25}\times\sqrt{2.12^2+100}} = 12.72 \;(22.1\ \text{dB}) \;\Rightarrow\; \beta = 12.72

Step 5: Corner frequencies

1T=ωgc′10=0.212 rad/s,T=4.71 s1βT=0.21212.72=0.0167 rad/s,βT=59.9 s\begin{aligned} \frac{1}{T} &= \frac{\omega_{gc}'}{10} = 0.212\ \text{rad/s},\quad T = 4.71\ \text{s} \\ \frac{1}{\beta T} &= \frac{0.212}{12.72} = 0.0167\ \text{rad/s},\quad \beta T = 59.9\ \text{s} \end{aligned}

Step 6: Compensator and check

Gc(s)=1+4.71s1+59.9s=0.0786 s+0.212s+0.0167G_c(s) = \frac{1 + 4.71s}{1 + 59.9s} = 0.0786\,\frac{s+0.212}{s+0.0167} Gc(s)G(s)=1500(1+4.71s)s(s+5)(s+10)(1+59.9s)G_c(s)G(s) = \frac{1500(1+4.71s)}{s(s+5)(s+10)(1+59.9s)}
QuantityUncompensatedCompensated
Gain crossover9.77 rad/s2.14 rad/s
Phase margin−17.2∘-17.2^\circ (unstable)49.6∘49.6^\circ
Phase crossover7.07 rad/s6.86 rad/s
Gain margin−6.0-6.0 dB15.5 dB
KvK_v30 s−1^{-1}30 s−1^{-1}

A practical RC lag network (R1R_1 in series, R2R_2–CC in shunt) gives T=R2CT = R_2C and β=(R1+R2)/R2\beta = (R_1+R_2)/R_2; e.g. C=10 μC = 10\ \muF gives R2=471R_2 = 471 kΩ\Omega and R1=5.52R_1 = 5.52 MΩ\Omega.

Answer: K=30K = 30, Gc(s)=1+4.71s1+59.9sG_c(s) = \dfrac{1+4.71s}{1+59.9s}; Kv=30 s−1K_v = 30\ \text{s}^{-1}, PM=49.6∘ (≥45∘)PM = 49.6^\circ \ (\ge 45^\circ).

  • 2068 Chaitra · 4 marks

Write a short note on PD and PI controller.

Answer

PD controller

The control signal is proportional to the error and its rate of change:

u(t)=Kpe(t)+Kdde(t)dt,Gc(s)=Kp+Kds=Kd(s+KpKd)u(t) = K_p e(t) + K_d \frac{de(t)}{dt}, \qquad G_c(s) = K_p + K_d s = K_d\left(s + \frac{K_p}{K_d}\right)
  • Adds a zero at s=−Kp/Kds = -K_p/K_d; acts like a lead compensator (high-pass).
  • Anticipates the error, so it adds damping: reduces overshoot, rise time and settling time.
  • Improves stability and increases bandwidth.
  • System type is unchanged, so steady-state error is not removed.
  • Amplifies high-frequency noise because of the derivative action.

PI controller

The control signal is proportional to the error and its integral:

u(t)=Kpe(t)+Ki∫e(t) dt,Gc(s)=Kp+Kis=Kp(s+Ki/Kp)su(t) = K_p e(t) + K_i \int e(t)\,dt, \qquad G_c(s) = K_p + \frac{K_i}{s} = \frac{K_p(s + K_i/K_p)}{s}
  • Adds a pole at the origin and a zero; acts like a lag compensator (low-pass).
  • Raises the system type by one, so the steady-state error to a step (type-0 plant) becomes zero.
  • Reduces bandwidth; the response becomes slower and may overshoot more.
  • Can reduce relative stability if KiK_i is too large; suffers integral wind-up.
  • Filters high-frequency noise.
FeaturePDPI
Added singularityZeroPole at origin + zero
ImprovesTransient responseSteady-state accuracy
System typeUnchangedIncreased by 1
NoiseAmplifiedFiltered
  • 2068 Baisakh (old course) · 1+3 marks

State whether the following statement is true or false and justify: Proportional controller makes the steady state error zero.

Answer

False. A proportional controller only reduces the steady-state error; it cannot make it zero (unless the plant already contains an integrator).

Justification. With Gc(s)=KpG_c(s) = K_p and a type-0 plant G(s)G(s) in a unity feedback loop, the step error is

ess=lim⁡s→0s⋅1s1+KpG(s)=11+KpG(0)e_{ss} = \lim_{s\to 0} \frac{s\cdot\frac{1}{s}}{1 + K_p G(s)} = \frac{1}{1 + K_p G(0)}

Example: G(s)=1s+1G(s) = \dfrac{1}{s+1}, so G(0)=1G(0) = 1.

KpK_pess=1/(1+Kp)e_{ss} = 1/(1+K_p)
10.5
90.1
990.01

The error falls as KpK_p rises but becomes zero only for Kp→∞K_p \to \infty, which is not possible. Very high gain also makes higher-order systems oscillatory or unstable and saturates the actuator.

A P controller adds no pole at the origin, so it does not change the system type. Zero steady-state error needs integral action (PI or PID controller), which adds a pole at s=0s = 0.

  • 2068 Baisakh (old course) · 16 marks

Design a lead compensator for a system having open loop transfer function 4/[s(s+2)], such that the designed system should have %Mp ≤ 16.3% and settling time (ts) ≤ 2 sec.

Answer

A lead compensator is designed by the root locus method: find the desired dominant poles, find the angle the compensator must add, place its zero and pole, then find its gain.

Step 1: Desired closed-loop poles

Peak overshoot:

Mp=e−πζ/1−ζ2=0.163ζ=−ln⁡0.163π2+(ln⁡0.163)2=1.8143.628=0.5\begin{aligned} M_p &= e^{-\pi\zeta/\sqrt{1-\zeta^2}} = 0.163 \\ \zeta &= \frac{-\ln 0.163}{\sqrt{\pi^2 + (\ln 0.163)^2}} = \frac{1.814}{3.628} = 0.5 \end{aligned}

Settling time (2% criterion):

ts=4ζωn=2  ⇒  ζωn=2,ωn=20.5=4 rad/st_s = \frac{4}{\zeta\omega_n} = 2 \;\Rightarrow\; \zeta\omega_n = 2,\quad \omega_n = \frac{2}{0.5} = 4\ \text{rad/s} sd=−ζωn±jωn1−ζ2=−2±j3.464s_d = -\zeta\omega_n \pm j\omega_n\sqrt{1-\zeta^2} = -2 \pm j3.464

Step 2: Uncompensated system

Closed-loop characteristic equation: s2+2s+4=0s^2 + 2s + 4 = 0, so ωn=2\omega_n = 2, ζ=0.5\zeta = 0.5.

  • Mp=16.3%M_p = 16.3\% (already correct), ts=40.5×2=4t_s = \dfrac{4}{0.5\times2} = 4 s (too slow).
  • Kv=lim⁡s→0s 4s(s+2)=2 s−1K_v = \lim_{s\to0} s\,\dfrac{4}{s(s+2)} = 2\ \text{s}^{-1}.

So the poles must move further left along the same ζ=0.5\zeta = 0.5 line, which needs a lead compensator.

Step 3: Angle deficiency

At sd=−2+j3.464s_d = -2 + j3.464:

∠sd=180∘−tan⁡−13.4642=120∘∠(sd+2)=90∘∠G(sd)=−(120∘+90∘)=−210∘ϕ=−180∘−(−210∘)=30∘\begin{aligned} \angle s_d &= 180^\circ - \tan^{-1}\frac{3.464}{2} = 120^\circ \\ \angle (s_d + 2) &= 90^\circ \\ \angle G(s_d) &= -(120^\circ + 90^\circ) = -210^\circ \\ \phi &= -180^\circ - (-210^\circ) = 30^\circ \end{aligned}

The compensator must add +30∘+30^\circ at sds_d.

Step 4: Locate zero and pole (bisector method)

Let Gc(s)=Kcs+zs+pG_c(s) = K_c\dfrac{s+z}{s+p}. At point P (−2,j3.464)(-2, j3.464) draw a horizontal line PA to the left and the line PO to the origin. Angle APO =120∘= 120^\circ. Draw its bisector PB, then two lines at ±ϕ/2=±15∘\pm\phi/2 = \pm15^\circ from PB; they cut the real axis at the zero and the pole.

            P (-2, j3.46)
           /|\ 
          / | \  bisector +-15 deg
         /  |  \
 -------x---o--x------x-----o---> Re
      -5.46 -2.93  -2        0
       pole  zero  plant   plant

The lines from P make angles 255∘255^\circ (zero) and 225∘225^\circ (pole) with the positive real axis:

z=2+3.464tan⁡75∘=2+0.928=2.93p=2+3.464tan⁡45∘=2+3.464=5.46\begin{aligned} z &= 2 + \frac{3.464}{\tan 75^\circ} = 2 + 0.928 = 2.93 \\ p &= 2 + \frac{3.464}{\tan 45^\circ} = 2 + 3.464 = 5.46 \end{aligned}

Check: zero angle 75∘75^\circ, pole angle 45∘45^\circ, net lead =30∘= 30^\circ.

Step 5: Gain

Kc=∣sd∣ ∣sd+2∣ ∣sd+5.46∣4 ∣sd+2.93∣=4×3.464×4.8994×3.586=4.73Gc(s)=4.73 s+2.93s+5.46\begin{aligned} K_c &= \frac{|s_d|\,|s_d+2|\,|s_d+5.46|}{4\,|s_d+2.93|} = \frac{4 \times 3.464 \times 4.899}{4 \times 3.586} = 4.73 \\ G_c(s) &= 4.73\,\frac{s+2.93}{s+5.46} \end{aligned}

Kv=lim⁡s→0s GcG=4.73×4×2.932×5.46=5.07 s−1K_v = \lim_{s\to0} s\,G_cG = \dfrac{4.73\times4\times2.93}{2\times5.46} = 5.07\ \text{s}^{-1}.

Closed-loop poles: −2±j3.464-2 \pm j3.464 and −3.46-3.46. The closed-loop zero at −2.93-2.93 is not close enough to the third pole to cancel it, so a step simulation gives Mp≈21%M_p \approx 21\%, ts≈2.04t_s \approx 2.04 s, slightly above the limits.

Step 6: Final design (pole–zero cancellation)

To meet Mp≤16.3%M_p \le 16.3\% exactly, place the zero on the plant pole at −2-2. The pole must give 90∘−30∘=60∘90^\circ - 30^\circ = 60^\circ:

p=2+3.464tan⁡60∘=4,4Kc=∣sd∣ ∣sd+4∣=4×4=16⇒Kc=4p = 2 + \frac{3.464}{\tan 60^\circ} = 4, \qquad 4K_c = |s_d|\,|s_d+4| = 4 \times 4 = 16 \Rightarrow K_c = 4 Gc(s)=4(s+2)s+4,GcG=16s(s+4)G_c(s) = \frac{4(s+2)}{s+4}, \qquad G_cG = \frac{16}{s(s+4)}

Closed loop: s2+4s+16=0s^2 + 4s + 16 = 0, so ωn=4\omega_n = 4, ζ=0.5\zeta = 0.5.

QuantityUncompensatedBisector designFinal design
Gc(s)G_c(s)14.73s+2.93s+5.464.73\frac{s+2.93}{s+5.46}4s+2s+44\frac{s+2}{s+4}
MpM_p16.3%about 21%16.3%
tst_s4 sabout 2.04 s2 s
KvK_v25.074

Answer: lead compensator Gc(s)=4(s+2)s+4G_c(s) = \dfrac{4(s+2)}{s+4} gives Mp=16.3%M_p = 16.3\% and ts=2t_s = 2 s with Kv=4 s−1K_v = 4\ \text{s}^{-1}.

  • 2066 Bhadra (old course) · 8 marks

Mention P, I and D controllers. Also explain the role of PI and PD controllers on transient and steady state performance specification.

Answer

P, I and D controllers

  • Proportional (P): u(t)=Kpe(t)u(t) = K_p e(t), Gc(s)=KpG_c(s) = K_p. Output is proportional to the present error. Raising KpK_p speeds up the response and reduces (but does not remove) steady-state error; too much gain gives overshoot and instability.
  • Integral (I): u(t)=Ki∫e dtu(t) = K_i \int e\,dt, Gc(s)=Ki/sG_c(s) = K_i/s. Output depends on the accumulated past error. It adds a pole at the origin and removes steady-state offset, but adds 90° of phase lag and slows the response.
  • Derivative (D): u(t)=Kd de/dtu(t) = K_d\, de/dt, Gc(s)=KdsG_c(s) = K_d s. Output depends on the rate of change of error (future trend). It adds damping, but is never used alone because it gives no output for a constant error and amplifies noise.

Role of the PI controller

Gc(s)=Kp+KisG_c(s) = K_p + \dfrac{K_i}{s}, a pole at the origin and a zero at −Ki/Kp-K_i/K_p.

SpecificationEffect of PI
System typeIncreased by one
Steady-state errorStep error becomes zero (type-0 plant); ramp error falls
Rise timeSlightly reduced
OvershootIncreases
Settling timeIncreases (slower)
StabilityRelative stability decreases
BandwidthDecreases, so noise is filtered

Role of the PD controller

Gc(s)=Kp+KdsG_c(s) = K_p + K_d s, a zero at −Kp/Kd-K_p/K_d. For G=ωn2s(s+2ζωn)G = \dfrac{\omega_n^2}{s(s+2\zeta\omega_n)} with PD, the characteristic equation becomes s2+(2ζωn+Kdωn2)s+Kpωn2=0s^2 + (2\zeta\omega_n + K_d\omega_n^2)s + K_p\omega_n^2 = 0, so the effective damping ratio increases.

SpecificationEffect of PD
DampingIncreased
Peak overshootReduced
Rise time and settling timeReduced (faster)
StabilityImproved
Steady-state errorType unchanged; error depends only on KpK_p
BandwidthIncreased; noise is amplified

So PI is chosen to improve steady-state performance and PD to improve transient performance. A PID controller combines both.

  • 2066 Jestha (old course) · 10 marks

The open loop transfer function of a unity feedback system is given by G(s) = K/[s²(0.2s+1)]. Design a lead compensator to meet the following specifications: Acceleration error constant = 10, Phase margin = 35°.

Answer

A type-2 plant has −180∘-180^\circ phase at low frequency, so the plant pole makes it unstable for any KK. A phase-lead network is designed on the Bode plot: Gc(s)=1+Ts1+αTsG_c(s) = \dfrac{1+Ts}{1+\alpha Ts}, α<1\alpha < 1 (unity DC gain, so KaK_a is set by KK).

Step 1: Gain for KaK_a

Ka=lim⁡s→0s2 Ks2(0.2s+1)=K=10,G(s)=10s2(0.2s+1)K_a = \lim_{s\to0} s^2\,\frac{K}{s^2(0.2s+1)} = K = 10, \qquad G(s) = \frac{10}{s^2(0.2s+1)}

Step 2: Uncompensated phase margin

∣G(jω)∣=10ω21+0.04ω2=1  ⇒  ωgc=2.94 rad/s∠G=−180∘−tan⁡−1(0.2×2.94)=−210.4∘PM=−30.4∘(unstable)\begin{aligned} |G(j\omega)| &= \frac{10}{\omega^2\sqrt{1+0.04\omega^2}} = 1 \;\Rightarrow\; \omega_{gc} = 2.94\ \text{rad/s} \\ \angle G &= -180^\circ - \tan^{-1}(0.2\times2.94) = -210.4^\circ \\ PM &= -30.4^\circ \quad (\text{unstable}) \end{aligned}

Step 3: Phase lead required

ϕm=35∘−(−30.4∘)+5∘=70.4∘\phi_m = 35^\circ - (-30.4^\circ) + 5^\circ = 70.4^\circ

This is more than about 60∘60^\circ, which a single lead network cannot give with a sensible α\alpha. Use two identical lead stages in cascade. Also, the plant magnitude falls at −40-40 to −60-60 dB/decade, so the crossover moves a lot and its phase falls quickly; a first trial with 2×35∘2 \times 35^\circ gives only PM=24.5∘PM = 24.5^\circ. So take 2×45∘2 \times 45^\circ:

α=1−sin⁡45∘1+sin⁡45∘=0.1716 per stage\alpha = \frac{1-\sin45^\circ}{1+\sin45^\circ} = 0.1716 \text{ per stage}

Step 4: New gain crossover

The two stages together add 2×10log⁡(1/α)=15.32 \times 10\log(1/\alpha) = 15.3 dB at ωm\omega_m. So ωm\omega_m is where ∣G(jω)∣=α=0.1716|G(j\omega)| = \alpha = 0.1716 (−15.3-15.3 dB):

10ω21+0.04ω2=0.1716  ⇒  ωm=6.08 rad/s\frac{10}{\omega^2\sqrt{1+0.04\omega^2}} = 0.1716 \;\Rightarrow\; \omega_m = 6.08\ \text{rad/s}

Step 5: Corner frequencies

T=1ωmα=16.08×0.4142=0.397 szero: 1/T=2.52 rad/s,pole: 1/(αT)=14.69 rad/s\begin{aligned} T &= \frac{1}{\omega_m\sqrt{\alpha}} = \frac{1}{6.08\times0.4142} = 0.397\ \text{s} \\ \text{zero: } 1/T &= 2.52\ \text{rad/s},\qquad \text{pole: } 1/(\alpha T) = 14.69\ \text{rad/s} \end{aligned}

Step 6: Compensator and check

Gc(s)=(1+0.397s1+0.0681s)2=34.0(s+2.52s+14.69)2G_c(s) = \left(\frac{1+0.397s}{1+0.0681s}\right)^2 = 34.0\left(\frac{s+2.52}{s+14.69}\right)^2 Gc(s)G(s)=10(1+0.397s)2s2(0.2s+1)(1+0.0681s)2G_c(s)G(s) = \frac{10(1+0.397s)^2}{s^2(0.2s+1)(1+0.0681s)^2}

At ω=6.08\omega = 6.08: plant phase =−180∘−tan⁡−1(1.216)=−230.6∘= -180^\circ - \tan^{-1}(1.216) = -230.6^\circ, compensator phase =+90∘= +90^\circ, so PM=180∘−230.6∘+90∘=39.4∘PM = 180^\circ - 230.6^\circ + 90^\circ = 39.4^\circ.

QuantityUncompensatedCompensated
Gain crossover2.94 rad/s6.08 rad/s
Phase margin−30.4∘-30.4^\circ39.4∘39.4^\circ
Gain marginnegative (unstable)11.5 dB at 14.5 rad/s
KaK_a10 s−2^{-2}10 s−2^{-2}

Answer: two cascaded lead stages, Gc(s)=(1+0.397s1+0.0681s)2G_c(s) = \left(\dfrac{1+0.397s}{1+0.0681s}\right)^2, give Ka=10K_a = 10 and PM=39.4∘ (≥35∘)PM = 39.4^\circ \ (\ge 35^\circ).

  • 2065 Shrawan (old course) · 16 marks

Design a lead compensator for a system with open loop transfer function G(s) = k/[s²(s+5)] for the specifications of phase margin = 30° and acceleration error constant Ka = 5 sec⁻². Also draw the bode magnitude and phase plots after compensation.

Answer

The plant is type 2 with an extra pole, so it is unstable for every kk. A phase-lead compensator Gc(s)=1+Ts1+αTsG_c(s) = \dfrac{1+Ts}{1+\alpha Ts} (unity DC gain) is designed on the Bode plot.

Step 1: Gain for KaK_a

Ka=lim⁡s→0s2 ks2(s+5)=k5=5  ⇒  k=25K_a = \lim_{s\to0} s^2\,\frac{k}{s^2(s+5)} = \frac{k}{5} = 5 \;\Rightarrow\; k = 25 G(s)=25s2(s+5)=5s2(0.2s+1)G(s) = \frac{25}{s^2(s+5)} = \frac{5}{s^2(0.2s+1)}

Step 2: Uncompensated phase margin

25ω2ω2+25=1  ⇒  ωgc=2.14 rad/s∠G=−180∘−tan⁡−1(2.14/5)=−203.2∘,PM=−23.2∘\begin{aligned} \frac{25}{\omega^2\sqrt{\omega^2+25}} &= 1 \;\Rightarrow\; \omega_{gc} = 2.14\ \text{rad/s} \\ \angle G &= -180^\circ - \tan^{-1}(2.14/5) = -203.2^\circ,\qquad PM = -23.2^\circ \end{aligned}

Step 3: Phase lead required

ϕm=30∘−(−23.2∘)+5∘=58.2∘\phi_m = 30^\circ - (-23.2^\circ) + 5^\circ = 58.2^\circ

A single stage (α=0.081\alpha = 0.081) moves the crossover to 3.75 rad/s, where the plant phase is already −216.9∘-216.9^\circ, so the PM is only 21.3∘21.3^\circ. Hence two identical lead stages of 35∘35^\circ each are used:

α=1−sin⁡35∘1+sin⁡35∘=0.271 per stage\alpha = \frac{1-\sin35^\circ}{1+\sin35^\circ} = 0.271 \text{ per stage}

Step 4: New gain crossover

Two stages add 2×10log⁡(1/α)=11.342\times10\log(1/\alpha) = 11.34 dB at ωm\omega_m, so ωm\omega_m is where ∣G∣=α=0.271|G| = \alpha = 0.271:

25ω2ω2+25=0.271  ⇒  ωm=3.83 rad/s\frac{25}{\omega^2\sqrt{\omega^2+25}} = 0.271 \;\Rightarrow\; \omega_m = 3.83\ \text{rad/s}

Step 5: Corner frequencies

T=13.830.271=0.502 s,1T=1.99 rad/s,1αT=7.35 rad/s,αT=0.136 sT = \frac{1}{3.83\sqrt{0.271}} = 0.502\ \text{s},\quad \frac{1}{T} = 1.99\ \text{rad/s},\quad \frac{1}{\alpha T} = 7.35\ \text{rad/s},\quad \alpha T = 0.136\ \text{s} Gc(s)=(1+0.502s1+0.136s)2=13.6(s+1.99s+7.35)2G_c(s) = \left(\frac{1+0.502s}{1+0.136s}\right)^2 = 13.6\left(\frac{s+1.99}{s+7.35}\right)^2 Gc(s)G(s)=25(1+0.502s)2s2(s+5)(1+0.136s)2G_c(s)G(s) = \frac{25(1+0.502s)^2}{s^2(s+5)(1+0.136s)^2}

PM check at 3.83 rad/s: 180∘−180∘−tan⁡−1(3.83/5)+2×35∘=−37.45∘+70∘=32.55∘≈32.6∘180^\circ - 180^\circ - \tan^{-1}(3.83/5) + 2\times35^\circ = -37.45^\circ + 70^\circ = 32.55^\circ \approx 32.6^\circ.

Step 6: Bode plots after compensation

ω\omega (rad/s)Magnitude (dB)Phase
0.154.0−177.0∘-177.0^\circ
0.526.5−165.3∘-165.3^\circ
115.6−153.4∘-153.4^\circ
1.996.8−142.0∘-142.0^\circ
3.830.0−147.4∘-147.4^\circ
7.35−8.4-8.4−176.1∘-176.1^\circ
10−13.7-13.7−193.3∘-193.3^\circ
20−28.7-28.7−227.0∘-227.0^\circ
50−51.5-51.5−252.1∘-252.1^\circ

Asymptotic magnitude slopes: −40-40 dB/dec up to 1.99 rad/s (14 dB at ω=1\omega = 1), 00 dB/dec from 1.99 to 5 rad/s (double zero), −20-20 dB/dec from 5 to 7.35 rad/s (plant pole), and −60-60 dB/dec beyond 7.35 rad/s (double pole).

 dB
 54 |\
    | \ -40
 7  |  \______ 0
 0  |---------\--------- 0 dB
    |          \ -20
-20 |           \
    |            \ -60
    +---+-----+--+------> w (log)
       1.99   5 7.35

 deg
-140|      .--.
-160|     /    \
-180|----'------\------- -180
-220|            \__
-250|               \___
    +---+-----+--+------> w (log)
       1.99 3.83 7.89

Phase rises from −180∘-180^\circ to a peak of about −142∘-142^\circ near 2 rad/s, is −147.4∘-147.4^\circ at the crossover 3.83 rad/s, and crosses −180∘-180^\circ at 7.89 rad/s, where the magnitude is −9.5-9.5 dB.

Answer: k=25k = 25, Gc(s)=(1+0.502s1+0.136s)2G_c(s) = \left(\dfrac{1+0.502s}{1+0.136s}\right)^2; Ka=5 s−2K_a = 5\ \text{s}^{-2}, PM=32.6∘ (≥30∘)PM = 32.6^\circ\ (\ge 30^\circ), GM=9.5GM = 9.5 dB.

  • 2081 Bhadra · 4 marks

Give a brief note on lead-lag compensator.

Answer

A lag–lead compensator is a single network that combines a lag section and a lead section, used when both the transient response and the steady-state accuracy must be improved.

Gc(s)=Kc(s+1/T1s+β/T1)⏟lead(s+1/T2s+1/(βT2))⏟lag,β>1G_c(s) = K_c \underbrace{\left(\frac{s + 1/T_1}{s + \beta/T_1}\right)}_{\text{lead}} \underbrace{\left(\frac{s + 1/T_2}{s + 1/(\beta T_2)}\right)}_{\text{lag}}, \qquad \beta > 1
  • The lag section (pole and zero very near the origin) acts at low frequency. It raises the low-frequency gain, so KvK_v or KaK_a increases and steady-state error falls.
  • The lead section acts near the gain crossover frequency. It adds positive phase, raising the phase margin, lowering overshoot and increasing bandwidth.
  • Frequencies are ordered 1βT2<1T2<1T1<βT1\dfrac{1}{\beta T_2} < \dfrac{1}{T_2} < \dfrac{1}{T_1} < \dfrac{\beta}{T_1}.

Electrical network: two resistors and two capacitors (R1∥C1R_1 \parallel C_1 in series, R2R_2 with C2C_2 in shunt) give this transfer function.

Design idea: design the lead part to meet the phase margin or the dominant-pole location, then the lag part to meet the error constant without disturbing the transient.

Its effect is similar to a PID controller: the lag part acts like PI and the lead part like PD.

  • 2081 Bhadra · 12 marks

Design a lag compensator for a unity feedback system with open loop transfer function G(s) = 0.53/[s(1+s)(1+0.5s)] such that velocity error constant would become 5 s⁻¹ without significant change in transient properties.

Answer

"Without significant change in transient properties" means the dominant closed-loop poles should stay almost where they are, so a lag compensator is designed by the root locus method, with its pole and zero close to the origin.

G(s)=0.53s(1+s)(1+0.5s)=1.06s(s+1)(s+2)G(s) = \frac{0.53}{s(1+s)(1+0.5s)} = \frac{1.06}{s(s+1)(s+2)}

Step 1: Uncompensated system

Characteristic equation: s3+3s2+2s+1.06=0s^3 + 3s^2 + 2s + 1.06 = 0.

Roots: s=−0.331±j0.586s = -0.331 \pm j0.586 (dominant) and s=−2.34s = -2.34.

ωn=0.3312+0.5862=0.673 rad/s,ζ=0.3310.673=0.491\omega_n = \sqrt{0.331^2 + 0.586^2} = 0.673\ \text{rad/s},\qquad \zeta = \frac{0.331}{0.673} = 0.491 Kv=lim⁡s→0s G(s)=0.53 s−1K_v = \lim_{s\to0} s\,G(s) = 0.53\ \text{s}^{-1}

Step 2: Required increase in KvK_v

β≥50.53=9.43⇒choose β=10\beta \ge \frac{5}{0.53} = 9.43 \quad\Rightarrow\quad \text{choose } \beta = 10

Step 3: Place the pole and zero

Gc(s)=K^c s+1/Ts+1/(βT)G_c(s) = \hat{K}_c\,\dfrac{s + 1/T}{s + 1/(\beta T)}. Put the zero close to the origin, well to the right of the dominant poles:

zero at s=−0.05,pole at s=−0.0510=−0.005\text{zero at } s = -0.05, \qquad \text{pole at } s = -\frac{0.05}{10} = -0.005

Angle added at s1=−0.331+j0.586s_1 = -0.331 + j0.586:

∠(s1+0.05)−∠(s1+0.005)=115.6∘−119.1∘=−3.5∘\angle(s_1 + 0.05) - \angle(s_1 + 0.005) = 115.6^\circ - 119.1^\circ = -3.5^\circ

This is less than 5∘5^\circ, so the root locus near the dominant poles hardly changes.

Step 4: New dominant poles and gain

Keeping ζ=0.491\zeta = 0.491, the compensated locus 1.06K^c(s+0.05)s(s+1)(s+2)(s+0.005)\dfrac{1.06\hat{K}_c(s+0.05)}{s(s+1)(s+2)(s+0.005)} crosses the ζ\zeta line at

s1′=−0.312±j0.553s_1' = -0.312 \pm j0.553

Magnitude condition at s1′s_1':

K^c=∣s1′∣ ∣s1′+1∣ ∣s1′+2∣ ∣s1′+0.005∣1.06 ∣s1′+0.05∣=0.970\hat{K}_c = \frac{|s_1'|\,|s_1'+1|\,|s_1'+2|\,|s_1'+0.005|}{1.06\,|s_1'+0.05|} = 0.970

Step 5: Compensator

Gc(s)=0.970 s+0.05s+0.005=9.70 20s+1200s+1G_c(s) = 0.970\,\frac{s+0.05}{s+0.005} = 9.70\,\frac{20s+1}{200s+1} Gc(s)G(s)=1.028(s+0.05)s(s+0.005)(s+1)(s+2)G_c(s)G(s) = \frac{1.028(s+0.05)}{s(s+0.005)(s+1)(s+2)}

Step 6: Check

Kv=lim⁡s→0s GcG=1.028×0.050.005×2=5.14 s−1  (≥5)K_v = \lim_{s\to0} s\,G_cG = \frac{1.028\times0.05}{0.005\times2} = 5.14\ \text{s}^{-1} \;(\ge 5)
QuantityUncompensatedCompensated
Dominant poles−0.331±j0.586-0.331 \pm j0.586−0.312±j0.553-0.312 \pm j0.553
ζ\zeta0.4910.491
ωn\omega_n0.673 rad/s0.635 rad/s
Other poles−2.34-2.34−2.33-2.33, −0.055-0.055
KvK_v0.53 s−1^{-1}5.14 s−1^{-1}

The extra closed-loop pole at −0.055-0.055 lies almost on the zero at −0.05-0.05, so its effect is a very small slow tail. Overshoot is unchanged and ωn\omega_n falls by only about 6%, so the transient response is practically the same while KvK_v rises nearly ten times.

Answer: Gc(s)=0.970s+0.05s+0.005G_c(s) = 0.970\dfrac{s+0.05}{s+0.005}, giving Kv=5.14 s−1K_v = 5.14\ \text{s}^{-1} with nearly unchanged transient response.

  • 2080 Baisakh · 6 marks

Explain PID controller with block diagram and its transfer function, along with its characteristics.

Answer

A PID controller produces a control signal that is the sum of three terms: proportional to the error, to the integral of the error, and to the derivative of the error.

u(t)=Kpe(t)+Ki∫0te(t) dt+Kdde(t)dtu(t) = K_p e(t) + K_i \int_0^t e(t)\,dt + K_d \frac{de(t)}{dt}

Block diagram

              +----------+
          +-->|   Kp     |----+
          |   +----------+    |
 r  +  e  |   +----------+  + v  u   +-------+   c
 --->(O)--+-->|  Ki / s  |-->(O)---->| Plant |--+-->
     -^   |   +----------+    ^      +-------+  |
      |   |   +----------+    |                 |
      |   +-->|  Kd s    |----+                 |
      |       +----------+                      |
      +-----------------------------------------+

Transfer function

Gc(s)=U(s)E(s)=Kp+Kis+Kds=Kds2+Kps+KisG_c(s) = \frac{U(s)}{E(s)} = K_p + \frac{K_i}{s} + K_d s = \frac{K_d s^2 + K_p s + K_i}{s}

or, in standard form, Gc(s)=Kp(1+1Tis+Tds)G_c(s) = K_p\left(1 + \dfrac{1}{T_i s} + T_d s\right) with integral time Ti=Kp/KiT_i = K_p/K_i and derivative time Td=Kd/KpT_d = K_d/K_p.

It adds one pole at the origin and two zeros to the open-loop transfer function.

Characteristics

ActionRise timeOvershootSettling timeSteady-state error
Increase KpK_pDecreasesIncreasesSmall changeDecreases
Increase KiK_iDecreasesIncreasesIncreasesEliminated
Increase KdK_dSmall changeDecreasesDecreasesNo change
  • Raises the system type by one, so steady-state error to a step is zero.
  • Derivative action adds damping and improves stability and transient response.
  • Gives good overall performance; used in most industrial loops (temperature, pressure, flow, motor speed).
  • Tuning is done by trial, by Ziegler–Nichols rules or by root locus/frequency methods.
  • In practice the derivative term is filtered, Kds/(1+τs)K_d s/(1+\tau s), to limit noise, and anti-windup is added to the integral term.
  • 2080 Baisakh · 12 marks

Design a suitable lag compensating network for G(s) = k/[s(s+2)(s+20)] to meet the following specification: Kv = 20 sec⁻¹, P.M ≥ 35°.

Answer

A lag network Gc(s)=1+Ts1+βTsG_c(s) = \dfrac{1+Ts}{1+\beta Ts} (β>1\beta > 1) is designed by the Bode plot method.

Step 1: Gain for KvK_v

Kv=lim⁡s→0s ks(s+2)(s+20)=k40=20  ⇒  k=800K_v = \lim_{s\to0} s\,\frac{k}{s(s+2)(s+20)} = \frac{k}{40} = 20 \;\Rightarrow\; k = 800 G(s)=800s(s+2)(s+20)=20s(0.5s+1)(0.05s+1)G(s) = \frac{800}{s(s+2)(s+20)} = \frac{20}{s(0.5s+1)(0.05s+1)}

Step 2: Uncompensated margins

  • ∣G(jω)∣=1|G(j\omega)| = 1 at ωgc=6.03\omega_{gc} = 6.03 rad/s; ∠G=−90∘−tan⁡−16.032−tan⁡−16.0320=−178.4∘\angle G = -90^\circ - \tan^{-1}\frac{6.03}{2} - \tan^{-1}\frac{6.03}{20} = -178.4^\circ, so PM=1.6∘PM = 1.6^\circ.
  • ωpc=40=6.32\omega_{pc} = \sqrt{40} = 6.32 rad/s, GM=0.8GM = 0.8 dB.

The system is on the verge of instability.

Step 3: New gain crossover frequency

Allow 10∘10^\circ for the lag of the network:

∠G(jωgc′)=−180∘+35∘+10∘=−135∘\angle G(j\omega_{gc}') = -180^\circ + 35^\circ + 10^\circ = -135^\circ −90∘−tan⁡−1ω2−tan⁡−1ω20=−135∘  ⇒  ωgc′=1.69 rad/s-90^\circ - \tan^{-1}\frac{\omega}{2} - \tan^{-1}\frac{\omega}{20} = -135^\circ \;\Rightarrow\; \omega_{gc}' = 1.69\ \text{rad/s}

Step 4: Value of β\beta

∣G(j1.69)∣=8001.69×1.692+4×1.692+400=9.02  (19.1 dB)  ⇒  β=9.02|G(j1.69)| = \frac{800}{1.69\times\sqrt{1.69^2+4}\times\sqrt{1.69^2+400}} = 9.02\;(19.1\ \text{dB}) \;\Rightarrow\; \beta = 9.02

Step 5: Corner frequencies

1T=1.6910=0.169 rad/s,T=5.92 s1βT=0.1699.02=0.0187 rad/s,βT=53.4 s\begin{aligned} \frac{1}{T} &= \frac{1.69}{10} = 0.169\ \text{rad/s},\quad T = 5.92\ \text{s} \\ \frac{1}{\beta T} &= \frac{0.169}{9.02} = 0.0187\ \text{rad/s},\quad \beta T = 53.4\ \text{s} \end{aligned}

Step 6: Compensator and check

Gc(s)=1+5.92s1+53.4s=0.111 s+0.169s+0.0187G_c(s) = \frac{1+5.92s}{1+53.4s} = 0.111\,\frac{s+0.169}{s+0.0187} Gc(s)G(s)=800(1+5.92s)s(s+2)(s+20)(1+53.4s)G_c(s)G(s) = \frac{800(1+5.92s)}{s(s+2)(s+20)(1+53.4s)}
QuantityUncompensatedCompensated
Gain crossover6.03 rad/s1.69 rad/s
Phase margin1.6°39.9°
Gain margin0.8 dB19.2 dB
KvK_v20 s−1^{-1}20 s−1^{-1}

(With only 5∘5^\circ allowance the PM comes out at 35.0∘35.0^\circ, just on the limit; the 10∘10^\circ allowance gives a safe margin.)

Answer: k=800k = 800, Gc(s)=1+5.92s1+53.4sG_c(s) = \dfrac{1+5.92s}{1+53.4s} (zero at −0.169-0.169, pole at −0.0187-0.0187); Kv=20 s−1K_v = 20\ \text{s}^{-1}, PM=39.9∘≥35∘PM = 39.9^\circ \ge 35^\circ.

  • 2079 Bhadra · 12 marks

The forward path transfer function of unity feedback system is given by G(s) = K/[s(s+10)(s+2)]. Design a suitable location of pole zero pair for a lead compensating network so that phase margin is at least 50° and velocity error constant is maintained at least 800 sec⁻¹.

Answer

A lead network 1+Ts1+αTs\dfrac{1+Ts}{1+\alpha Ts} (unity DC gain) is placed so that its maximum phase lead occurs at the new gain crossover frequency.

Step 1: Gain for KvK_v

Kv=lim⁡s→0s Ks(s+10)(s+2)=K20=800  ⇒  K=16000K_v = \lim_{s\to0} s\,\frac{K}{s(s+10)(s+2)} = \frac{K}{20} = 800 \;\Rightarrow\; K = 16000 G(s)=16000s(s+2)(s+10)=800s(0.5s+1)(0.1s+1)G(s) = \frac{16000}{s(s+2)(s+10)} = \frac{800}{s(0.5s+1)(0.1s+1)}

Step 2: Uncompensated margins

∣G(jω)∣=1 at ωgc=24.5 rad/s∠G=−90∘−tan⁡−124.52−tan⁡−124.510=−243.2∘PM=−63.2∘,GM=−36.5 dB at ωpc=4.47 rad/s\begin{aligned} |G(j\omega)| &= 1 \text{ at } \omega_{gc} = 24.5\ \text{rad/s} \\ \angle G &= -90^\circ - \tan^{-1}\frac{24.5}{2} - \tan^{-1}\frac{24.5}{10} = -243.2^\circ \\ PM &= -63.2^\circ,\quad GM = -36.5 \text{ dB at } \omega_{pc} = 4.47\ \text{rad/s} \end{aligned}

The system is badly unstable.

Step 3: Phase lead needed

ϕm=50∘−(−63.2∘)+5∘≈118∘\phi_m = 50^\circ - (-63.2^\circ) + 5^\circ \approx 118^\circ

One lead network gives at most about 60∘60^\circ in practice, so three identical lead stages are used. A first trial of 3×40∘3 \times 40^\circ gives only PM≈43∘PM \approx 43^\circ because the crossover shifts, so 3×45∘3 \times 45^\circ is chosen:

α=1−sin⁡45∘1+sin⁡45∘=0.1716 per stage\alpha = \frac{1-\sin45^\circ}{1+\sin45^\circ} = 0.1716 \text{ per stage}

Step 4: New gain crossover

Three stages add 3×10log⁡(1/α)=22.963\times10\log(1/\alpha) = 22.96 dB at ωm\omega_m, so ωm\omega_m is where ∣G∣=α1.5=0.0711|G| = \alpha^{1.5} = 0.0711:

16000ωω2+4ω2+100=0.0711  ⇒  ωm=60.6 rad/s\frac{16000}{\omega\sqrt{\omega^2+4}\sqrt{\omega^2+100}} = 0.0711 \;\Rightarrow\; \omega_m = 60.6\ \text{rad/s}

Step 5: Pole–zero locations

T=1ωmα=160.55×0.4142=0.0399 szero: s=−1T=−25.1pole: s=−1αT=−146.2\begin{aligned} T &= \frac{1}{\omega_m\sqrt{\alpha}} = \frac{1}{60.55\times0.4142} = 0.0399\ \text{s} \\ \text{zero: } s &= -\frac{1}{T} = -25.1 \\ \text{pole: } s &= -\frac{1}{\alpha T} = -146.2 \end{aligned} Gc(s)=(1+0.0399s1+0.00684s)3=197.9(s+25.1s+146.2)3G_c(s) = \left(\frac{1+0.0399s}{1+0.00684s}\right)^3 = 197.9\left(\frac{s+25.1}{s+146.2}\right)^3

Step 6: Check

GcG=16000(1+0.0399s)3s(s+2)(s+10)(1+0.00684s)3G_cG = \frac{16000(1+0.0399s)^3}{s(s+2)(s+10)(1+0.00684s)^3}

At ω=60.6\omega = 60.6: plant phase =−90∘−88.1∘−80.6∘=−258.7∘= -90^\circ - 88.1^\circ - 80.6^\circ = -258.7^\circ; lead =3×45∘=135∘= 3\times45^\circ = 135^\circ; PM=180∘−258.7∘+135∘=56.3∘PM = 180^\circ - 258.7^\circ + 135^\circ = 56.3^\circ.

QuantityUncompensatedCompensated
Gain crossover24.5 rad/s60.6 rad/s
Phase margin−63.2∘-63.2^\circ56.3∘56.3^\circ
Gain margin−36.5-36.5 dB13.6 dB
KvK_v800 s−1^{-1}800 s−1^{-1}

Answer: each stage has its zero at s=−25.1s = -25.1 and pole at s=−146.2s = -146.2 (α=0.1716\alpha = 0.1716); three such stages give PM=56.3∘≥50∘PM = 56.3^\circ \ge 50^\circ with Kv=800 s−1K_v = 800\ \text{s}^{-1}. (Such a large KvK_v needs a large lead; in practice a lag–lead compensator would give a lower bandwidth and less noise.)

  • 2078 Bhadra · 4 marks

In response to unit ramp input, discuss role of derivative feedback controller for a second order system.

Answer

Derivative (rate) feedback feeds back Kt dcdtK_t\,\dfrac{dc}{dt} in addition to the output, for example with a tachogenerator on a position servo. It increases damping but, for a ramp input, it increases the steady-state error.

Take the standard second-order plant with derivative feedback in a minor loop:

G(s)=ωn2s(s+2ζωn),H(s)=1+KtsG(s) = \frac{\omega_n^2}{s(s + 2\zeta\omega_n)}, \qquad H(s) = 1 + K_t s

Closed-loop transfer function:

C(s)R(s)=ωn2s2+(2ζωn+Ktωn2)s+ωn2\frac{C(s)}{R(s)} = \frac{\omega_n^2}{s^2 + (2\zeta\omega_n + K_t\omega_n^2)s + \omega_n^2}

Effect on damping

Comparing with s2+2ζ′ωns+ωn2s^2 + 2\zeta'\omega_n s + \omega_n^2:

ζ′=ζ+Ktωn2\zeta' = \zeta + \frac{K_t\omega_n}{2}

The natural frequency is unchanged, but damping rises, so overshoot and oscillation fall and the system becomes more stable.

Effect on the unit ramp response

Error E(s)=R(s)−C(s)E(s) = R(s) - C(s) with R(s)=1/s2R(s) = 1/s^2:

E(s)=s2+(2ζωn+Ktωn2)ss2+(2ζωn+Ktωn2)s+ωn2⋅1s2ess=lim⁡s→0sE(s)=2ζωn+Ktωn2ωn2=2ζωn+Kt\begin{aligned} E(s) &= \frac{s^2 + (2\zeta\omega_n + K_t\omega_n^2)s}{s^2 + (2\zeta\omega_n + K_t\omega_n^2)s + \omega_n^2}\cdot\frac{1}{s^2} \\ e_{ss} &= \lim_{s\to0} sE(s) = \frac{2\zeta\omega_n + K_t\omega_n^2}{\omega_n^2} = \frac{2\zeta}{\omega_n} + K_t \end{aligned}

Without rate feedback ess=2ζ/ωne_{ss} = 2\zeta/\omega_n. Rate feedback adds KtK_t to the ramp error.

Example: ωn=4\omega_n = 4, ζ=0.25\zeta = 0.25, Kt=0.15K_t = 0.15: ζ′=0.25+0.3=0.55\zeta' = 0.25 + 0.3 = 0.55 (overshoot drops from 44% to 12.6%) but esse_{ss} rises from 0.125 to 0.275.

Remedy: increase the forward gain (amplifier) to bring the ramp error back down; the extra damping from KtK_t keeps the response well damped. In contrast, a PD controller in the forward path adds similar damping while leaving KvK_v, and so the ramp error 2ζ/ωn2\zeta/\omega_n, unchanged. Derivative feedback therefore trades ramp accuracy for damping.

  • 2076 Chaitra · 3 marks

Write a short note on characteristics of PI and PD control actions.

Answer

PI control action

Gc(s)=Kp+KisG_c(s) = K_p + \dfrac{K_i}{s}; output depends on the present error and the accumulated (integrated) error.

  • Adds a pole at the origin, so the system type increases by one.
  • Removes steady-state error to a step for a type-0 plant (no offset).
  • Makes the response slower, with more overshoot; relative stability falls.
  • Reduces bandwidth, so it filters noise. Behaves like a lag compensator.

PD control action

Gc(s)=Kp+KdsG_c(s) = K_p + K_d s; output depends on the present error and its rate of change.

  • Adds a zero, giving phase lead; acts in anticipation of the error.
  • Increases damping, so overshoot, rise time and settling time decrease.
  • Improves stability and increases bandwidth.
  • Does not change the system type, so steady-state error is not removed.
  • Amplifies high-frequency noise. Behaves like a lead compensator.

Questions from Old Question Collection (EE 602) (IOE Control System exam papers (EE 602 and older course), 2065 to 2082) and Old Question Collection (BEI, EE 504) (IOE BEI Control System (EE 504) exam papers, 2076 to 2081). Answers are written for this site; check them against your class notes.

Chapter titles and hours from the IOE syllabus ↗