Stability Analysis in Control Systems
Stability Analysis in Control Systems
an
E n g i n e e r i n g
(v"Mech. - VTU) 333 Stability of Control Systems
Control.
Coefficients
for first two rows are written directly from characteristic equation.
row next rows can be obtained as follows.
these two
From
b, = b =
bË = a
F r o mg d
a n d r d row, row can be obtained as
bË a3 - a qb bË as - ajb3
C; =
bË bË
is to be continued till the coefficient for s
is obtained which will
process
This
From this array stability of system can be predicted.
be an
8.7.1 Routh's Criterion
H1x6-6 = 10
6
6
1 3 2
2 6 1
1.5
6 e-3 1
there
two sign changes, system is unstable.
are
As
Second
method :
To solve the above difficulty one more method can be used. In
by 1/z' in original equation. Taking L.C.M., rearrange characteristic
replace 's'
in descending powers of 'z'. Then coplete the Routh's array with this new
his,
quation
d e f
d f
4d 2e
iv) Complete the array interms of these new
coefficients.
[Link] Importance of an Auxiliary Equation
Auxiliary equation is always the part of the original characteristic
means the roots of the equatiot
auxiliary equation are some of the roots of on
characteristic equation. Not only this but the roots of auxiliary
dominant roots of the original characteristic equation, from equation are tie
view. the stabuity p
Key point : The stability can be predicted from the roots of A(s) = 0 rather thanth
of characteristic equation as the roots of A(s) = 0 are the most dominant fromthe st
point of view. The remaining roots of the characteristic equation are always intheletk
they do not play any significant role in the
stability equation
e.g. Let F(s) = 0 is the original characteristic analysis. of say order n =5.
Let A(s) = 0 be the auxiliary equation for the System due to occurrence
case 2 of the order mn = 2.
D m
S-plane) from the stability point of view are the 2imaginary The
willbelar
A(s) = 0.
5-2 = 3roots are not significant roots of
view as they
from the imaginary axis in the left from stability point of
half of s-plane.
E n g i n e e r i n g
(vMech. - VTU) 337 Stability of Control Systems
Control
roots
of auxiliary equation may be,
The
of real
roots of opposite sign i.e. as shown in the Fig. 8.10 (a).
pair
A
)
jo
XX repeated
XD
Control
of which saves calculation time.
determinants,
evaluation
number of roots of characteristic
No
unstable
system it gives equation having
positive r e a l p a r t .
For
i)
stability of the system can be easily judged.
i) Relative
of
sustained oscillations can be determined.
)
range of values of K for system stability.
frequency
in finding
out
helps
It out intersection points of root locus with imaginary axis.
i)
in finding
helps
vii) It
L i m i t a t i o n s of Routh's Criterion
8.11
valid only for real coefficients of the characteristic equation.
Itis
il
does not provide
exact locations of the closed loop poles in left or right half of
ii) It
splane.
does not suggest methods of stabilising an unstable system.
iü) It
Applicable only to linear systems.
iv)
Examples with Solutions
Criterion :
Examine the Stability By Routh's
0
Ex 8.4 5 + 6s + 26s + 56s + 80 =
Sol. :
1 26 80
6 56
16.66 80
27.21
80
28 +12
(-12) -15[
26 +12
15
Lim 2: +12
= 2+== +0
&’0
0-144 - 0
=- 12
0+12
1 2 3
2 15
-12 0
+ oo 15 There are two sign
-12 changes, So system is
unstable
15
C o n t r o l
1 5 2 3 1
8 30 8 6 0
1.25 1 2.25 1 0
+ 4g5. +35 -16s -64s- 48 = 0 Find the number of roots of this equation with
Ex. 8.7
s
n0sitive real part, zero real part and negative
real part.
Sol. : - 48
,6 1 3 - 16
4 64
3 48
3 ,0
dA = 12 3
A(s) = 3s -48= 0
ds
1 3 - 16 - 48
0 64
3 - 48
12
48
[e]0
576
48
Lim 576 + 0
So s = t 2 are the two roots on imaginary axis i.e. with zero real part. Root in
RHS. indicated by asign change is s = + 2as obtained by solving A(s) = 0. Total
here are 6 roots
as n = 6.
Roots with positive real part = 1
Control Engineering (V"Mech. - VTU) 342
16.8 5K
4.2
K
For system to be stable there should not be
sign change in the first
K> 0 from sO column.
and 16.8 - 5K > 0 from s1
16.8 > 5K 3.36 > K .:. K< 3.36
Range of Kis 0 <K< 3.36
Ex. 8.9 For a system with
characteristic equation
F(s) =s +355 + 4s +6s + 5s + 3s+ 2 = 0,
Sol. :
examine stability.
1 4 2
3 6 3
2 4 2
Row of zeros
A(s) = 2s +4s2 +2 = 0 i.e. s +2s +1 = 0
dA(s) = 4s + 4s
ds
1 4 5 2
3 6 3
2 4 2
4 4
2 2
change,
hence is located in [Link]. of s-plane. As row of zeros
no root
No sign unstable To
mar8naly stable or imaginary
system may De auxiliary equation.
o that find the roots of first
-2 +/4-4
2s +1 =0 =. = -1
Als) = s +
-j
_² =-1, S1,2 =tj, S3,4 = tj
$² =-1,
solve second auxiliary
he roots of A'(s) = 0are the roots of A(s) = 0. So do not
Predict the stability from the nature of roots of first auxiliary equation.
equation.
(T + 2) 5K
(T+ 2) (K+ 1)- 2Tx 5K
(T+ 2)
5K
From last row 5K > 0
Kmust be
positive.
Control Engineering (V"Mech. - VTU) 344
Unstable 2 +T Considering
values
positive
9T - 2 of K and
Stable Unstable
4.2 K
16.8 5K
4.2
K
For system to be stable there should not be sign change in the first column.
K> 0 from sO
and 16.8 5K > 0 from s'
Engineering
(vMech. - VTU) 345 Stability of Control Systems
C o n t r o l ,
K<6.5
.:. Range of values of K, 0<
Now marginal value of K is that value of K for which system becomes
marginally stable system there must be a row of zeros
narginally stable. For a
-ccurring in Routh's array. So value of K which makes any row of Routh array as
bw of zeros is called Marginal value of K. Now K 0makes row of s as row of
eros but K =0 can not be marginal value because for K =0, constant term in
naracteristic equation becomes zero i.e. one coefficient for s vanishes which makes
stem unstable instead of marginally stable.
Hence marginal value of K is a value which makes any row other than s as row
zeros.
1 + K
At marginal value of K
(2 + K) p- (1 + K) = 0
(2 + K) p = (K + 1)
K+1
p =
K+2
Now at this value,
A(s) = Ps?.+K+1= 0
(K+1) K+1
S=tj
K+1 K+1
K+2 4
4 = K +2 .. Kmar = 2
Kmar +1 = 0.75
p =
4 4
8.14 The open loop transfer function of a unity fedback system is given by,
CÀ. K
G(s)
s(1 +sT) (1+sT)
Derive an expression for gain K interms ofT and T, for the stability of the system.
(April-96)
equation is,
Sol. : The characteristic
1+ G(s) H(s) = 0
K
1+
s(1+sT)(1+ sTz)
T,I,s' +s(T, +T;)+s+K = 0
The Routh's array is
1
(4+ h)-KHh 0
(4+ 7z)
K
From s K>0
So
is the range of K for stability.
Ex. 8.15 s+285 +85 + 126+20s+16s+ 16 =0
(CE-E&C, Aug.-2003, LCS-E&E, Aug.-2003)
Gol.:
1 8 20 16
2 12 16 0
2 12 16 0
Row of zeros
Control Engineering (v"'Mech. VTU) 348
Stablity of Contro
A(s) = 2s +126 +16 = 0
dA
8s+24s = 0
ds
20
16
2 12 16
0
2 12 16
0
8 24
6 16 0
2.67 0
16
= -3+1=-2, - 4
s = -2 and s = -4
s = tN2 and s =#j2
Nonrepeated roots on imaginary axis. Hence system is marginally stable.
Ex. 8.16 Find the range of K for which the system, whose characteristic equation
s
below, is stable.
F(s) = s + (K + 0.5)s + 4 Ks +50 = 0 (Aug
Sol. : The Routh's array for the equation is,
1 4K
K+0.5 50
(K+0.5) 4K- 50
(K+0.5)
50
From row of s K+ 0.5 >0
K>-0.5
From rowof s'
(K + 0.5) 4K - 50 >0
Systems
Stability of Control
349
(v'Mech. - VTU)
Engineering
>0
4K + 2K - 50
2K' + K - 25 > 0
-1t(1 +25x4x 2
K=
2x2
can write,
3.2943, -3.7943, hence we
are
values of the roots
3.7943) > 0
= (K- 3.2943) (K +
K> 3.2943 and
K> - 3.7943
stability is,
of values of Kfor the system
the ultimate
range
ce 3.2943 < K<o
determine the range of
K
Fig, below,
cosed-loop feedback
system shoun in oscillation.
- Fx the also the frequency of sustained
stable and
for zohich the system is
(CE-E&C March-2001)
K C(s)
R(S)
s(s +s+ 1)\(s +4)
16.8-5K
S
4.2
K
ns, K>0
n s, 16.8 - 5K > 0
16.8 > 5K
E n g i n e e r i n g
(v"Mech. - VTU) 353 Stability of Control Systems
conditions
necessary
are the
What
4
0alfof s-plane2.
special Icases of Routh's criterion.
note on
Write a
an. Auxiliaryequation.
Explain the significance of
6
Routh's
criterion can be used to studythe relative stability?
How
7,
and limitation of Routh's method.
State the advantages
8.
characteristic polynomial ofa system is
9 The
s ++ 6s +15s +30s + 44s+ 24 (Ans. : Marginally stable, o= 2rad/sec)
Pls)=
Determnine the. stability ofasystem with characteristic equation.
10.
4=0
S+ 4s + 2s + 8s + s+
(Ans. : Unstable, Multiple roots on imaginary axis at t)
following system.
11. Deternmine the range of Kfor stability of the
1+ 2s
(Ans. : 0 <K<o)
12. For K= 2, determine whether the following unity feedback system is stable. Use Routh criterion.
G(s) =
K(1+ 2s) (1+ 4s)
stis +25+ 10) (Ans. : stable)
13. Apply Routh criterion to check the stability of
s+9 +20s + 12s +8s+ 16s +16= 0 (Ans. : Unstable, 2 sign changes)
[Link] the range of Kfor stability ofs+ 2s + 2s+(3+ K)s+ K= 0,for K>0
(Ans. : 0 <K<0.464)
15. Determine the range of Kfor stability +2s +3s +4s' +K(s+ 1)+ s+ 5=0
(Ans. : No real value of K makes the system stable)
16. Find the (Ans. : 0 <K<390)
range offKfor whichhsystem given below is stable.
R(s) K C(s)
s (s+3) (s + 10)
Control Engineering (v"Mech. - VTU) 354
1 ForaSwstrm having characteristic equation.
Stability of
s' 10s (50+ A)s + K=0
Find the range of Kfor stability. For what value of Ksysten vill
oscillations? (Ans. : 0 < K < 1.55, oscillate and ohat t
Keritical = 1.55, ois =frequency of
19. By means of Routh's criterion, deternine stability of 0,816 rad
a) s+ 5s + 10s + 24s + 20s + 15s+ 10 =0
b) s + s' -3s -st+2
(Ans. : (a) unstable
20. Determine the range of value of K(K> 0) (b) ung
such that the characteristic equation.
s + 3(K+1)s + (7K+5)s+ (4K+7) =
0has roots more negative that s =-1.