117920
117920
Stability:
Stability refers to stable working condition of a control system. Stability can be defined as:
i) A system is stable, if its output is bounded (finite) foe any bounded (finite) input.
ii) A system is asymptotically stable, if in absence of input, the output goes zero
irrespective of initial condition.
iii) A system is limitedly stable, if for a bounded input signal, output has a constant
amplitude oscillation then system may be stable or unstable under some limited
constraints.
iv) A system is absolutely stable, if a system output is stable for all variations of its
parameter.
v) A system output is stable for a limited range of variation of its parameter, then it is
conditionally stable.
vi) A system is unstable, it for a bounded disturbing input signal, and output is infinite
or oscillatory.
C ( S ) bo S m b1 S m 1 b2 S m 2 ... bm1 S m bm
M (S )
R( S ) ao S n a1 S n 1 a 2 S n 2 ... a n 1 S n a n
( S Z1 )( S Z 2 )( S Z 3 )...(S Z M )
M (S ) ….(1)
( S P1 )( S P2 )( S P3 )...(S PN )
The roots of characteristics equation are poles and the roots of numerator polynomial are zeros.
The rots or poles may lie at origin or imaginary axis or right or left half of S-plane.
From equation (1), characteristics equation: ao S n a1 S n1 a2 S n2 ... an1 S an 0
the root of the characteristics equation determines the stability of the system. the types of roots
and response is given as below
1. If the roots of the characteristics equations have negative real parts, then the system is
stable.
2. If any roots of the characteristics equation have a positive real part, then the system is
unstable.
3. If the characteristics equation has repeated roots in the jw axis, the system is unstable.
4. If one or more non repeated roots of the characteristics equation are on the jw axis, the
system is unstable.
5. If condition (1) is satisfied, except for the presence of one or more non-repeated roots on
the jw axis, the system is limitedly stable.
Routh Array
Sn ao a2 a4 a6 . .
S n1 a1 a3 a5 . . .
S n2 b1 b2 b3 . . .
S n 3 c1 c2
n4
S d1 d2
. . .
. . .
. . .
S2 e1 an
S1 f1
S 0 an
The Routh stability criterion is stated as, “For a system to be stable, it is necessary and
sufficient that each term of the first column of the Routh Array of its characteristics equation
be positive if a0>0. If this condition is not met, the system is unstable and number of sign
changes of the terms of the first column of the Routh Array corresponds to the number of
roots of the characteristics equation in the right half of the S-plane.”
EXAMPLES:
1. Find the stability of the system with characteristics equations
i) S 8S 18S 16 S 5 0
4 3 2
ii) 3S 10 S 5S 5S 2 0
4 3 2
Soln:
i) The Routh Array is given by
S4 1 18 5
0
S3 8 16 (missing
term)
8 *18 1 *16 8 * 5 1* 0
S2 16 5
8 8
16 *16 8 * 5
S1 13.5 0
16
13.5 * 5 16 * 0
S0 5
13.5
The elements of the first column are all positive and hence the system is stable.
ii) The Routh Array is given by
S4 3 5 2
S3 10 5
S3* 2 1
2 * 5 3 *1 2 * 2 3* 0
S2 3.5 2
2 2
3.5 * 1 2 * 2 0.5
S1
3.5 3.5
0
S 2
*In order to simplify the computation S3- Row is divided by 5 throughout.
The elements of the first column has both positive and negative elements, with two sign changes-
representing two poles of the characteristics equation of the right half of the S-plane. and hence
the system is unstable.
Special Cases:
Case 1: When the first term in any row in Routh Array is Zero while rest of the row has at
least one non zero term. Because of this, the term in the next row becomes infinite.
There are two methods to solve this problem
^
i) Substitute a small positive number I for the zero and proceed to evaluate the rest of
the Routh Array. Then examine the signs of the first column of Routh Array by
^
lettering I as 0
ii) Modify the original characteristics equation by replacing S by 1/Z. Apply the Routh’s
test on the modified equation in the term of Z. the number of Z roots with positive
real parts are the same as the number of S roots with positive real parts. This method
works in most but not all cases.
3. Find the Stability of the system with characteristic equation
S S 2S 2S 3S 5 0
5 4 3 2
Method 1:
The Routh Array is
S5 1 2 3
S4 1 2 5
1* 2 1* 2
S3 0 -2
1
2 2
S2 5
4 4 5 2
S1 2
2 2
S0 5
From the Routh array, it is seen that first element in the third row is 0. This is replaced by , a
2 2
small positive number. The first element in the 4th row is now which has a positive sign
4 4 5 2
as 0 . The first term in the fifth row is 2 as 0 . Therefore by
2 2
examining the first column it is noted that there is two sign changes and hence the system
becomes unstable.
Method 2:
When S is replaced by 1/Z, the characteristics equation becomes
5Z 3Z 2Z 2Z z 1 0
5 4 3 2
z3 -4/3 -2/3
z2 ½ 1
z1 2
z0 1
The elements of the first column has both positive and negative elements, with two sign changes,
hence the system is unstable.
Case 2: When all the elements in any row in the Routh Array are zero.
This condition indicates that there are symmetrically located roots (±real roots,
±conjugate roots or complex conjugate roots). The polynomial whose coefficient are the
elements of the rows just above the row of zero is called the Auxiliary polynomial. Replace the
rows of zeros with the coefficients obtained by differentiating the auxillary polynomial with
respect to S.
4. Find the Stability of the system with characteristic equation
S 2S 8S 12 S 20 S 16 S 16 0 .
6 5 4 3 2
It is noted that S3- row has all elements zero, thus considering the row above (S4 row) to
be auxiliary polynomial.
A( S ) S 4 6S 2 8
Differentiating with respect to S
dA( S )
4S 3 12 S ; replace the obtained coefficient in S3- row
dS
S6 1 8 20 16
5
S 2 12 16
S5 1 6 8
4
S 2 12 16
S4 1 6 8
S3 4 12
S3 1 3
S2 3 8
S1 1/3
S0 8
Since there is no sign change in the first column after considering the auxiliary
polynomial. The system is said to be limitedly stable because of the consideration used.
By solving the auxiliary polynomial A( S ) S 4 6S 2 8
The roots are S=±j√2 and S=±j2.
The root locus was introduced by W.R.Evans in 1948 for the analysis of control systems.
The root locus technique is a powerful tool for adjusting the location of closed loop poles to
achieve the desired system performance by varying one or more system parameters. It is a
graphical representation in s-domain and it is symmetrical about the real axis. Because the open
loop poles and zeros exist in the s-domain having the values either as real or as complex
conjugate pairs.
Consider a feedback system with characteristics equation is 1+G(S)H(S)=0
Rules for Construction of Root Locus
Follow these rules for constructing a root locus.
Rule 1 − Locate the open loop poles and zeros in the ‘s’ plane.
Rule 2 − Find the number of root locus branches.
We know that the root locus branches start at the open loop poles and end at open loop zeros.
So, the number of root locus branches N is equal to the number of finite open loop poles P or
the number of finite open loop zeros Z, whichever is greater.
Mathematically, the number of root locus branches N can be written as,
N=P if P≥Z
N=Z if P<Z
Rule 3 − Identify and draw the real axis root locus branches.
If the angle of the open loop transfer function at a point is an odd multiple of 1800, then that
point is on the root locus. If odd number of the open loop poles and zeros exist to the left side of
a point on the real axis, then that point is on the root locus branch. Therefore, the branch of
points which satisfies this condition is the real axis of the root locus branch.
If there exists a real axis root locus branch between two open loop zeros, then there will
be a break-in point in between these two open loop zeros.
Note − Break-away and break-in points exist only on the real axis root locus branches.
Follow these steps to find break-away and break-in points.
Write K in terms of ss from the characteristic equation 1+G(s)H(s)=0.
Differentiate K with respect to s and make it equal to zero. Substitute these values of s in
the above equation.
The values of s for which the K value is positive are the break points.
Rule 7 − Find the angle of departure and the angle of arrival.
The Angle of departure and the angle of arrival can be calculated at complex conjugate
open loop poles and complex conjugate open loop zeros respectively.
The formula for the angle of departure ϕd is
d 180
The formula for the angle of arrival ϕa is
a 180
Where p z
Example 1
A certain unity negative feedback control system has the following open loop Transfer function
K
GH (s) . Draw the root locus for 0 K
s( s 1)( s 3)
Solution:
1. The three root loci start from poles at 0,-1 and -3 where the value of K=0
2. There is no zero here and so the three loci terminate at infinity where k=infinity.
3. For the given open loop transfer function, P=3 and Z=0. Hence, N=3. There are three
separate loci.
4.
The asymptote of the root loci make
2n 1 angles with real axis. Hence the
P Z
5
angles are , and . The loci 1,2 and 3 make these angles respectively.
3 3
5. The point of interaction of the asymptotes is given by
0 1 3 0
x Poles of GH (s) Zeros of GH (s) 1.33
3
6. At point A there is no root locus, because two poles are lying to the right side of point
A. For even number of poles and zeros to the right side of the point there cannot be any
root locus on the real axis. On the other had consider the point B. There are three poles
on the right hand side of point B. This is an odd number. Hence, three is root locus.
Hence, on the real axis root loci exist between s=0 and s=-1 and s=-3 and s=-infinity.
7. The characteristics equation is written as
F (s) 1 G ( s) 0
s 3 4 s 2 3s K 0
Differentiating the above equation with respect to s we get,
dK
3s 2 8s 0
ds
dK
Put 0 and solving for s we get
ds
S1=-2.23 and s2=-0.43.
The breakaway point at s1=-2.23 is ruled out since there is no root locus there. Hence
s2=-0.43 is the breakaway point.
8. The critical value of K and the value of w at the imaginary axis where the root locus
crosses from LHP to RHP is obtained from the following Routh’s array which of
formed from the characteristic equation.
S3 1 3
S2 4 K
S1 12 K
0
4
S0 K
The critical value of K=12. The auxiliary equation is obtained from the s2 row as
4s 2 K 0
Substituting s=jω and K=12, we get
4 2 12 0 or 3
9. The complete root locus is shown in Fig
“ If the G(s)H(s) contour in the G(s)H(s) plane corresponding to Nyquist contour in the s-plane
encircles the point -1+j0 in the anticlockwise direction as many times as the number of right half
s-plane poles of G(s)H(s), then the closed loop system is stable.
1. No encirclement of—1+ j0 point: This implies that the system is stable if there are no
poles of G(s)H(s) in the right half s-plane. If there are poles on right half s-plane then the
system is unstable.
2. Anticlockwise encirclements of —1 + j0 point: In this case the system is stable if the
number of anticlockwise encirclements is same as the number of poles of G(s)H(s) in the
right half s-plane. If the number of anticlockwise encirclements is not equal to number of
poles on right half s-plane then the system is unstable.
3. Clockwise encirclements of the —1 + j0 point: In this case the system is always
unstable. Also in this case, if no poles of G(s)H(s) in right half s-planes then the number
of clockwise encirclement is equal to number of poles of closed loop system on right half
s-plane.
Example 1
Draw the Nyquist plot for the system whose open loop transfer function is
K
G s H s
s(s 2)( s 10)
Determine the range of K for which closed loop system is stable.
Solution:
C1-nyquist
C3 –inv of nyquist
C2
R=∞ & θ=π/2 to - π/2
C4
R=0 & θ =-π/2 to π/2
K 0.05 K
G (s) H (s)
s ( s 2)( s 10) s (1 0.5s )(1 0.1s )
s j
0.05 K
G (s) H (s)
0.6 j (1 0.05 2 )
2
At pc , imag 0, (1 0.05 2 ) 0
pc 4.472 rad / sec
0.05 K 0.05 K
G (s) H (s)
0.6 2 0.6(4.472) 2
0.00417 K
Mapping of section C1
Mapping of section C2
j
s Lt R e to
R
2 2
0.05 K 0.05 K
G(s)H (s)
s (1 0.5 s )(1 0.1s ) s (0.5 s )(0.1s )
K K
j 3
0 e j 3
Lt (R e )
3
s
R
0 3
3
2
to
3
2
Mapping of section C4
j
s Lt R e to
R0
2 2
0.05 K 0.05 K
G(s) H (s)
s (1 0.5 s )(1 0.1s ) s (1)(1)
0.05 K K
j
e j
s Lt (R e )
R0
to
2 2
Stability Analysis
When -0.00417K=-1, the contour passes through (-1+j0) point and corresponding value of K is
the limiting value of K for stability.
1
Limiting value of K 240
0.00417
When K<240
When K is less than 240, the contour crosses real axis at a point between 0 and -1+j0. On
travelling through Nyquist plot along the indicated direction it is found that the point -1+j0 is not
encircled. Also the open loop transfer function has no poles on the right half of s-plane.
Therefore the closed loop system is stable.
When K>240 greater than 240, the contour crosses real axis at a point between -1+j0 and -∞. On
travelling through Nyquist plot along the indicated direction it is found that the point -1+j0 is
encircled in clockwise direction two times.( Since there are two clockwise encirclement and no
right half open loop poles, the closed loop system has two poles on right half of s-plane).
Therefore the closed loop system is unstable.
Result
The value for stability is 0<K<240
Example 2
By Nyquist stability criterion determine the stability of the closed loop system, whose open loop
s2
transfer function is give by G(s)H(s) . Comment on the stability of open loop and closed
(s+1)(s-1)
loop system.
Mapping of C1
s2 2(1 0.5s)
G( s) H ( s)
s 1 s 1 1 s 1 s
Let s=jw
2(1 0.5 j ) 2 1 0.25 2 tan 1 (0.5 )
G ( j ) H ( j )
1 j 1 j 1 2 tan 1 ( ) 1 2 180 tan 1 ( )
2 1 0.25 2
G( j ) H ( j )
1 2
G( j ) H ( j ) 180 tan 1 (0.5 )
ω rad/sec 0 0.4 1 2 10 ∞
G ( j ) H ( j ) 2 1.76 1.12 0.57 0.1 0
G ( j ) H ( j ) -180 -168 -153 -135 -101 -90
deg
Mapping of Section C2
The mapping of section C2 from s plane to G(s)H(s) plane is obtained by letting
s lim Re j in G(s)H(s) and varying θ from +π/2 to –π/2.
R
1
G ( s) H ( s) s lim Re j
0 e j
R lim Re j
R
j
when , G ( s ) H ( s ) 0e 2
2
j
when , G ( s ) H ( s ) 0e 2
2
Mapping of Section C3
In section C3, ω varies from -∞ to 0. The mapping of section C3, is given by the locus of
G(jω)H(jω) as w is varied from -∞ to 0.This locus is the inverse polar plot of G(jω)H(jω).
Stability Analysis
Result