Control Systems
Control Systems
ENGINEERING COLLEGE
TIRUCHIRAPALLI-621213.
QUESTION BANK
UNIT-1
The relationship between an input variable and an output variable of a signal flow
graph is given by the net gain between input and output nodes and is known as overall gain
of the system. Masons gain formula is used to obtain the overall gain (transfer function) of
signal flow graphs.
Gain P is given by
1
P Pk k
k
is determinant of graph
k is cofactor of kth forward path determinant of graph with loops touching kth forward
path. It is obtained from by removing the loops touching the path Pk.
Closed-loop systems are also called feedback control systems. Feedback is the
property of the closed-loop systems which permits the output to be compared with the input
of the system so that appropriate control action may be formed as a function of inputs and
outputs.
Solution:
X ( s)
Transfer function =
F ( s)
Step 1:
The opposing force acting on Mass M1 are fm1, fb, fk1, fb and fk
d 2 X1
w.k.t. f m1 m1
dt 2
dX 1
f b1 B1
dt
f k1 k1 X 1
d X1 X
fb B
dt
fk k X1 X
f m1 f b1 fb f k1 fk 0 - - - - - (1)
d 2 X1 dX 1 d X1 X
M1 B1 K1 X 1 B K X1 X 0 - - - - - (2)
dt 2 dt dt
Taking Laplace transform for equation (2)
2
M 1 X 1 (s) B1 SX1 (s) BS ( X 1 (s) X (s)) K1 X 1 (s) K ( X 1 (s) X (s)) 0
X 1 ( s)[ M 1 s 2 ( B1 B) s ( K1 K )] X ( s)[ Bs K]
( Bs K )
X 1 (s) X (s) 2
- - - - - (3)
M 1s ( B1 B ) s ( K 1 K)
Step 2:
The opposing force acting on Mass M2 are fm2, fb, fk, fb2
d2X
w.k.t. f m 2 M2
dt 2
d
fb B (X X1)
dt
fk K(X X1)
dX
f b2 B2
dt
f m2 fb fk f b2 f (t )
d2X d dx
M2 B X X1 K X X1 B2 f (t ) - - - - - (4)
dt 2 dt dt
( Bs K ) 2
X ( s )[ M 2 s 2 ( B2 B) s K ] X (s) F ( s)
M 1s 2 ( B1 B) s ( K1 K)
[M 2 s 2 ( B2 B) s K ][ M 1 s 2 ( B1 B) s ( K 1 K ) ( Bs K)
X ( s) F ( s)
M 1 s 2 ( B1 B) s ( K 1 K )
X (s) M 1s 2 ( B1 B) s ( K 1 K )
2
F (s) [M 1s ( B1 B) s ( K1 K )][ M 2 s 2 ( B2 B) s K ] ( Bs k ) 2
Result:
d 2 X1 dX 1 d X1 X
1. M 1 B1 K1 X 1 B K X1 X 0
dt 2 dt dt
d2X d dx
2. M 2 B X X1 K X X1 B2 f (t )
dt 2 dt dt
X ( s) M 1s 2 ( B1 B) s ( K1 K )
1. 2
F ( s) [M 1s ( B1 B) s ( K1 K )][M 2 s 2 ( B2 B) s K ] ( Bs k ) 2
vr (t ) Vr sin c t;
vs1 (t ) KVr sin t cos( 120 ) vs1s 2 (t ) 3KVr sin c t sin( 240 )
c
When 0 , vs 3s1 (t ) 0 and maximum voltage is induced on S2 coil. This position of the rotor
is defined as the electrical zero of the transmitter and used as reference position of the
rotor.
The output of the synchro transmitter is applied to the stator winding of a synchro control
transformer. Circulating current of the same phase but of different magnitude flows through
the two sets of stator coil. The pair acts as an error detector.
The voltage induced in the control transformer rotor is proportional to the cosine of the angle
between the two rotors and is given by,
Thus the synchro transmitter-control transformer pair acts as an error detector which gives a
voltage signal at the rotor terminal of the control transformer proportional to the angular
difference between the shaft positions.
Equation (01) is represented graphically in figure below for an arbitrary time variation of (
).
It is seen that the output of the synchro error detector is a modulated signal, where the ac
signal applied to the rotor of synchro transmitter acts as carrier and the modulating signal is,
em (t ) K s ( ); K s K Vr
3.Write the rules for block diagram reduction techniques. (May12, Nov06)
Block diagram:
Lower case letters usually represent functions of time. Upper case letters
denote Laplace transformed quantities as a function of complex variable s or Fourier
transformed quantity i.e. frequency function as function of imaginary variable s j .
Take off point: In order to employ the same signal or variable as an input to more than
block or summing point, take off point is used. This permits the signal to proceed unaltered
along several different paths to several destinations.
4. (i) Derive the transfer function for Armature controlled DC motor. (May11, Nov11)
K f i f ; TM K1K f i f ia KT ia
d di
eb Kb ; La a Raia eb e
dt dt
d2 d
J f0 TM KT ia
dt 2 dt
In Laplace domain,
Eb ( s) Kb s (s)
( La s Ra ) I a ( s) E ( s ) Eb ( s ) ;
2
( Js f 0 s) ( s) KT I a ( s )
( s) KT
G( s)
E ( s) s ( Ra sLa )( Js f 0 ) KT Kb
Neglecting La ,
KT / Ra KT / Ra Km
G( s) 2
;
Js s ( f 0 KT K b / Ra ) s ( Js f ) s(s m 1)
where, f f0 KT K b / Ra and K m KT / Ra f ; m J / f.
K m and m are called the motor gain and time constant respectively. These two parameters
are usually supplied by the manufacturer.The block diagram model is,
B. Transfer Function of a Field-controlled DC Motor
Figure 1
TM K1 ia K1 K f i f ia K T' i f
di f
Lf Rf if e ;
dt
d2 d
J f TM KT i f
dt 2 dt
( L f s R f ) I f ( s ) E (s )
( Js 2 fs) ( s) TM ( s) KT I f ( s)
We obtain,
( s) KT
G( s)
E ( s) s( L f s R f )( Js f)
5. Explain the working of AC servomotor in control systems.
A.C. Servomotor
For low power application a.c. motors are preferred, because of their light weight,
ruggedness and no brush contact. Two phase induction motors are mostly used in control
system.
The motor has two stator windings displaced 90 electrical degree apart. The voltages
applied to the windings are not balanced. Under normal operating conditions a fixed voltage
from a constant voltage source is applied to one phase. The other phase, called control
phase, is energized by a voltage of variable magnitude which is 90 out of phase w.r.t. the
voltage of fixed phase.
The torque speed characteristic of the motor is different from conventional motor. X / R ratio
is low and the curve has negative slope for stabilization.
The torque-speed curve is not linear. But we assume it as linear for the derivation of transfer
function. The troque is afunction of both speed and the r.m.s. control voltage, ie.,
TM f ( , E) .
Using Tailor series expansion about the normal operating point (TM 0 , 0 , E0 ) we get,
TM TM
TM TM 0 ( E E0 ) ( 0 )
E E E0 E E0
0 0
TM TM
TM kE m J f0 ; where, k and m .
E E E0 E E0
0 0
(s) k Km k J
Or, G(s) 2
; where, K m and m .
E (s) Js ( f 0 m) s s ( m s 1) f0 m f0 m
Since m is negative the transient part is decaying as m is positive. If m would positive and
m f 0 the transient part will increase with time and the system would be unstable. k and m
are the slope of the torque-voltage and torque-speed curve.
6. Determine the overall transfer function C(S)/R(S) for the system shown in fig.
(Nov10,05)
To simplify the inner feedback loop to obtain the following block diagram.
8. Find the overall gain C(s) / R(s) for the signal flow graph shown below
(May09) G6 G7
R(s) G1 C(s)
G2 G3 G4 G5 1
X1 X2 X3 X4
X5
-H1
-H2
P2=G1G6G4G5
P3= G1G2G7
L1=-G4H1,
L2=-G2G7H2,
L3= -G6G4G5H2 ,
L4=-G2G3G4G5H2
There is one combination of Loops L1 and L2 which are non-touching with loop gain product
L1L2=G2G7H2G4H1
= 1+G4H1+G2G7H2+G6G4G5H2+G2G3G4G5H2+ G2G7H2G4H1
Therefore 1= 1, 2= 1.
Hence, 3= 1+G4H1.
10. Find the overall gain C(s) / R(s) for the signal flow graph shown below
(May09)
X6 X5 X3
Find the gains , , for the signal flow graph shown in Fig.4.
X1 X 2 X1
b -h
X1 a d X5 f X6
c e
X2 X3 X4
-g
-i
P2=abef
L1=-cg,
L2=-eh,
L3= -cdei,
L4=-bei
There is one combination of Loops L1 and L2 which are nontouching with loop gain product
L1L2=cgeh
= 1+cg+eh+cdei+bei+cgeh
Therefore 1= 1, 2= 1.
X6 P1 1 P2 2 cdef abef
T=
X1 1 cg eh cdei bei cgeh
X5
Case 2:
X2
b -h
X2 X5 X5
1 d e 1
c
X2 X3 X4
-g
-i
The transfer function can directly manipulated from case 1 as branches a and f are removed
which do not form the loops. Hence,
X5 P1 1 P2 2 cde be
T=
X2 1 cg eh cdei bei cgeh
X3
Case 3:
X1
The signal flow graph is redrawn to obtain the clarity of the functional relation as shown in
Fig.6.
-h
X1 a X2 b e X5 f 1 X3
X4 X3
-i d
-g
P2=ac
There is one combination of Loops L1 and L2 which are nontouching with loop gain product
L1L2=ehcg
= 1+eh+cg+bei+efd+befg+ehcg
X3 P1 1 P2 2 abef ac 1 eh
The transfer function T =
X1 1 eh cg bei efd befg ehcg
MAHALAKSHMI
ENGINEERING COLLEGE
TIRUCHIRAPALLI-621213.
QUESTION BANK
UNIT-2
PART A
It is the response of the system when the input changes from one state to
another
2. Name the test signals used in time response analysis. (May 12)
Step A A/s
Ramp At A / s2
Unit Ramp t 1 / s2
Parabolic At2 / 2 A / s3
Unit Parabolic t2 / 2 1 / s3
Impulse (t) 1
3. How is system classified depending on the value of damping? (May11)
Definition:
1. Undamped system = 0
2. Under damped system 0 < < 1
3. Critically damped system = 1
4. Over damped system > 1
1 1
The transient response of a practical C.S often exhibits damped oscillations before
reaching steady state. In specifying the transient response characteristic of a C.S to unit
step i/p, it is common to specify the following terms.
The steady state error of the system for a unit step input is
s 1 1
e ss lim
s 0 1 G s s 1 G0
Kp lim G s G0
s 0
Thus, the steady state error in terms of the static position error constant Kp is given
by
1
e ss
1 Kp
K T s 1 T s 1 T s 1
a b m 1
Kp lim K and e ss
s 0 1 K
Ts 1 T s 1 T s 1
1 2 p
K T s 1 T s 1 T s 1
a b m
Kp lim N 1 and e ss
s 0
N
s Ts 1 T s 1 T s 1
1 2 p
Response of a feedback control system to a step input involves a steady state error if
there is no integration in the feed forward path. If a zero steady state error for a step input is
desired, the type of the system must be one or higher.
It is the time required for the response to reach the 1st of peak of the overshoot.
It is the maximum peak value of the response curve measured from unity. The
amount of max over shoot directly indicates the relative stability of the system.
10. Consider a first order system given by G(S)=C(S)/R(S) = K/S+a .Plot its unit
step response. (May09)
C(t)
r(t) 0.95
0.865
0.632
t=0 t 0 T 2T 3T t
1. Draw the response of second order system for under damped case and when
input is unit step. (May 12)
The standard form of closed loop transfer function of second order system is
2
C ( s) n
R( s) s2 2 n s 2
n
For under damped system, 0 < < 1 and root of the denominator are complex
conjugate.
2
s n n 1
Since < 1, 2 is also less than 1
2
Therefore, the damped frequency of oscillation, d n 1
s n j d
2
n
The response in s domain, C ( s) R( s)
s2 2 n s 2
n
2
n
C (s)
s( s 2 2 n s 2
n )
2
n A Bs C
C ( s) 2 2 2 2
s( s 2 n s n ) s (s 2 ns n )
2
n A(s 2 2 n s 2
n ) (Bs c)s
1. Put s = 0
2 2
n A n
A=1
0=1+B
B = -1
3. Equating co-efficient of s
0 = 2An + C
= 2n + C
C = - 2n
1 s 2 n
C (s)
s s2 2 n s 2
n
2 2
Let us add and subtract n to denominator of the second term
1 s 2 n
C (s)
s s2 2 n s 2
n
2 2
n
2 2
n
1 s 2 n
s s2 2 n s 2
n
2 2
n
2 2
n
1 s 2 n
2 2 2 2
s (s n ) n n
1 s 2 n
2 2 2
s (s n ) n (1 )
1 s 2 n
2 2
s (s n ) d
1 s n n
2 2
s (s n ) d (s n ) d
1 s n n d
C ( s) 2 2
s (s n ) d d (s n )2 2
d
1 s
c(t ) L 1{C ( s)} L 1 2
n
2
n d
s (s n ) d d (s n )2 2
d
nt n nt
1 e cos d t e sin d t
d
nt n
1 e cos d t sin d t
d
nt n
1 e cos d t sin d t
2
n 1
nt
1 e cos d t sin d t
2
1
nt
e 2
c(t ) 1 1 cos d t sin d t
2
1
Note:
2
On constructing on right angle triangle with and 1 , we get
2
1 1
2
sin = 1
cos =
2
1-
tan
nt
e
c(t ) 1 sin cos d t cos sin d t
2
1
nt
e
c(t ) 1 sin( d t )
2
1
2
1 1
tan
nt
e
Unit step response =
sin( d t )
2
1
nt
e
Step Response =
A sin( d t )
2
1
r (t) c(t)
1 1
The response of the under damped second order system for unit step input sketched and
observed that the response oscillates before settling to a final value. The oscillation depends
on the value of damping ratio.
2. Derive the expressions for Rise time, Peak time, Peak overshoot, delay time
(May 12)
The performance characteristics of a controlled system are specified in terms of the
transient response to a unit step i/p since it is easy to generate & is sufficiently drastic.
The transient response of a practical C.S often exhibits damped oscillations before
reaching steady state. In specifying the transient response characteristic of a C.S to unit
step i/p, it is common to specify the following terms.
Response curve
It is the time required for the response to reach 50% of its final value for
st
the 1 time.
2) Rise time :- (tr)
It is the time required for the response to rise from 10% and 90% or 0%
to 100% of its final value. For under damped system, second order system the 0 to
100% rise time is commonly used. For over damped system, the 10 to 90% rise time is
commonly used.
It is the time required for the response to reach the 1st of peak of the overshoot.
It is the maximum peak value of the response curve measured from unity. The
amount of max over shoot directly indicates the relative stability of the system.
It is the time required for the response curve to reach & stay with in a range about
the final value of size specified by absolute percentage of the final value (usually 5% to
2%). The settling time is related to the largest time const., of C.S.
3. A unity feedback control system has an open loop transfer function G(S)=
10/S(S+2).Find the rise time, percentage over shoot, peak time and settling
time.
4.For a unity feedback control system the open loop transfer function
G(S) = 10(S+2)/ S2(S+1).Find (a) position, velocity and acceleration error constants.
(Nov11)
An automatic controller compares the actual value of the system output with the reference
input (desired value), determines the deviation, and produces a control signal that will
reduce the deviation to zero or a small value. The manner in which the automatic controller
produces the control signal is called the control action. The controllers may be classified
according to their control actions as
1) Electric controllers
2) Hydraulic controllers
3) Pneumatic controllers
4) Electronic controllers.
automatic controller
actuating Controller
error
error signal e(t) output u(t) Output c(t)
reference detector
Amplifier Actuator Plant
input r(t)
In a two position control action system, the actuating element has only two positions which
are generally on and off. Generally these are electric devices. These are widely used are
they are simple and inexpensive. The output of the controller is given by Eqn.1.
ut U1 et 0
.. (1)
U2 et 0
Actuating Controller
Error
Error signal e(t) Output u(t)
Reference Detector
Input r(t)
Feed back
signal b(t)
Proportional Control Action
ut K pe t ..... (2)
Us
Kp .. (3)
Es
Actuating Controller
Error
Error signal e(t) Output u(t)
Reference Detector
Kp
Input r(t)
Feed back
signal b(t)
Figure 7 Block diagram of a proportional controller
The value of the controller output u(t) is changed at a rate proportional to the actuating error
signal e(t) given by Eqn.4
du t t
K i e t or u t K i e t dt ..(4)
dt 0
Us Ki
.(5)
Es s
If the value of e(t) is doubled, then u(t) varies twice as fast. For zero actuating error,
the value of u(t) remains stationary. The integral control action is also called reset control.
Fig.4 shows the block diagram of the integral controller.
Actuating Controller
Error
Error signal e(t) Output u(t)
Reference Detector Ki
Input r(t) s
Feed back
signal b(t)
Kp t
ut K pe t e t dt ..(6)
Ti 0
Us 1
Kp 1 .....(7)
Es Ti s
Where, Kp is the proportional gain, Ti is the integral time which are adjustable.
The integral time adjusts the integral control action, while change in proportional gain
affects both the proportional and integral action. The inverse of the integral time is called
reset rate. The reset rate is the number of times per minute that a proportional part of the
control action is duplicated. Fig.5 shows the block diagram of the proportional-plus-integral
controller. For an actuating error of unit step input, the controller output is shown in Fig.6.
Actuating Controller
Error
Error signal e(t) output u(t)
Reference Detector K p 1 Ti s
Input r(t) Ti s
Feed back
signal b(t)
2Kp
Kp proportional only
0 Ti t
de t
ut K pe t K p Td ..(8)
dt
Us
K p 1 Td s ..(9)
Es
Both, Kp and Td are adjustable. The derivative control action is also called rate
control. In rate controller the output is proportional to the rate of change of actuating error
signal. The derivative time Td is the time interval by which the rate action advances the effect
of the proportional control action. The derivative controller is anticipatory in nature and
amplifies the noise effect. Fig.7 shows the block diagram of the proportional-plus-derivative.
For an actuating error of unit ramp input, the controller output is shown in Fig.8
actuating Controller
error
error signal e(t) output u(t)
reference detector
K P 1 TD S
input r(t)
feed back
signal b(t)
proportional only
Td
0 t
It is a combination of proportional control action, integral control action and derivative control
action. The equation of the controller is
Kp t de t
ut K pe t e t dt K p Td ..(10)
Ti o dt
Us 1
Kp 1 Td s ..(11)
Es Ti s
Where, Kp is the proportional gain, Ti is the integral time, and Td is the derivative time. The
block diagram of PID controller is shown in Fig.9. For an actuating error of unit ramp input,
the controller output is shown in Fig.10.
actuating
Controller
error
error signal e(t)
output u(t)
reference detector K p 1 Ti s Ti Td s 2
Td s
input r(t)
feed back
signal b(t)
proportional only
0 t
Figure 14 Response of PID controller to unit actuating error signal
11. Consider a open loop transfer function G(S) of second order system with unity
feedback system. Express damping ratio of the closed loop system in terms of
a, wn, and , where a is a constant, wn is the natural frequency of oscillations
and is the damping ratio of open loop system.(nov 11
UNIT III
PART A
The Frequency response is the steady state response (output) of a system when the
input to the system is a sinusoidal signal. The frequency response of a system is normally
obtained by varying the frequency of the input signal by keeping the magnitude of the input
signal at a constant value.
The phase cross over frequency is the frequency at which the phase of the open loop
transfer function is 180.
1
Gain Margin Kg= at pc
G( j
6. How do you calculate the gain margin from the polar plot? (May10)
Gain Margin (GM):
o The gain margin is the reciprocal of magnitude at the frequency at which the
0
phase angle is -180 .
o In terms of dB
1
GM in dB 20log10 20log10 | G( jwc) | 20log10 ( x)
| G( jwc) |
7. How do you calculate the gain margin from the polar plot? (May10)
Phase margin is that amount of additional phase lag at the gain crossover
frequency required to bring the system to the verge of instability (marginally stabile)
0
=180 +
m
Where
=G(j )
g
8. How do you find the stability of the system by using polar plot? (May10)
Stable: If critical point (-1+j0) is within the plot as shown, Both GM & PM are +ve
Unstable: If critical point (-1+j0) is outside the plot as shown, Both GM & PM are -ve
Marginally Stable System: If critical point (-1+j0) is on the plot as shown, Both GM
& PM are ZERO
9. What is cut off rate? (Nov09)
The slop of the log magnitude curve near the cut off frequency is called cut
off rate.
10. Consider the first order system G(S) = 1/1+ST .Draw its polar plot. (Nov09)
PART B
1. Define resonant peak and resonant frequency. (May09)
Peak response Mp :
The peak response Mp is defined as the maximum value of M( ) that is given in
Eq.(1.4). In general, the magnitude of Mp gives an indication of the relative stability of
a feed back control system. Normally, a large Mp corresponds to a large peak
overshoot in the step response. For most design problems it is generally accepted
that an optimum value Mp of should be somewhere between 1.1 & 1.5.
Resonant frequency p:
The resonant frequency p is defined as the frequency at which the peak resonance
Mp
2. Sketch bode plot for the following transfer function and determine the
system gain K for the gain cross over frequency to be 5rad/sec
G(S) = Ks2/(1+0.2s)(1+0.02s)
Solution:
MAGNITUDE PLOT
The various values of G(j ) are below tabulation , in the increasing order of their corner
frequencies. also slope contributed by each term and the change in slope at the corner
frequency
Choose a low frequency l such that l< c1& Choose a high frequency h such that
Let A=|G(j)| in db
=20log50/5+28=48db
=0log (100/50)+48=48db
Let points a,b,c,d be the points corresponding to frequencies l, , c1, c2, h respectively on
the magnitude plot.
PHASE PLOT
= G(j)= 180-tan-10.2-tan-10.02
5 45 5.7 129.3130
50 84.3 45 50.750
CALCULATION OF K
At =5,gain= 28db
At =5,frequency the db gain should be zero
20log K=-28db
Log K=-28/20
K=10-28/20=0.0398
NOTE: The frequency rad/sec is a corner frequency .hence in the exact plot the db gain
at =5rad/sec will be 3db less than the approximate plot.therefore for exact plot the
20logK=-25
LogK=-25/20=>K=10-25//20=0.056
Example
Problem : Construct the polar plot for the (critically damped system)
defined by : 1
T(s)
(s 1) 2
1 Im
Solution : T( j )
(1 j ) 2 Re
Limiting conditions : -0.2 0 0.2 0.4 0.6 0.8 1
(i) 0 : T( j ) 1 j0 -0.2
(ii) : T( j ) 1
2
-0.4
j 180o
i.e. T( j ) 0e
-0.6
(iii) 1 : T( j )
j
2
wc2 ; A = [Slope from wc1 to wc2 x log (wc2 / wc1 ] + Gain (A)@ wc1
= - 40 log [ 0.33 / 0.25 ] + 38
= 33 dB
Awh ; A = [Slope from wc2 to wh x log (wh / wc2 ] + Gain (A) @ wc2
= - 60 log [ 3.3 / 0.33 ] + 33
= - 27 dB
Calculations of Gain cross over frequency
The frequency at which the dB magnitude is Zero wgc = 1.1 rad / sec
Gain Margin
The gain margin in dB is given by the negative of dB magnitude of G(jw) at phase cross over
frequency
UNIT IV
STABILITY ANALAYSIS
PART A
The roots of the closed-loop characteristic equation define the system characteristic
responses Their location in the complex s-plane lead to prediction of the
characteristics of the time domain responses in terms of:damping ratio,
natural frequency, wn
damping constant, first-order modes
Consider how these roots change as the loop gain is varied from 0 to
Consider an nth-order system whose the characteristic equation (which is also the
denominator of the transfer function) is
A linear time invariant system is said to be stable if following conditions are satisfied.
1. When system is excited by a bounded input, output is also bounded & controllable.
2. In the absence of input, output must tend to zero irrespective of the initial conditions.
6. How the roots of characteristic are related to stability?
On a given section of the real axis, root axis is found in the section only if the total
number of poles and zeros of GH(s) to the right of the section is odd if the feedback is Ve
and K is +Ve, root locus is found on the real axis where even number of poles and zeros
occurred to the right of the concerned point.
10. What is absolute stability and conditional stability of a system with respect to a
parameter? (Nov09)
S 4 + 2s 3 + 0s 2 + 14s + 15 = 0
S4 1 0 15
S3 2 14 0
S2 -7 15 0
S1 18.28 0 0
0
S 15
Two sign change, there are two roots to the right of s = -2 & remaining 2 are to the left
of the line s = -2. Hence the system is unstable.
2. s6 + 4s5 +3s4 16s2- 64s 48 = 0 Find the number of roots of this equation with
positive real part, zero real part & negative real part
S5 4 0 -64 0
S4 3 0 -48 0
S3 0 0 0
dA
A(s) = 3S4 48 = 0 = 12s3
ds
S6 1 3 -16 -48
S5 4 0 -64 0
S4 3 0 -48 0
S3 12 0 0 0
S2 ( )0 -48 0 0
S1 576 0 0 0
S0 -48
Lim 576
0 +
Therefore One sign change & system is unstable. Thus there is one root in R.H.S of the s
plane i.e. with positive real part. Now
solve A(s) = 0 for the dominant roots
A(s) = 3s4 48 =0
Put S2 = Y
3Y2 = 48 Y2 =16, Y= 16 = 4
S2 = + 4 S2 = -4
S = 2 S= 2j
So S = 2j are the two parts on imaginary axis i.e. with zero real part. Root in R.H.S.
indicated by a sign change is S = 2 as obtained by solving A(s) = 0. Total there are 6 roots
as n = 6.
3. Explain the step by step procedure rules for constructing root locus
1) The root locus is symmetrical about real axis. The roots of the characteristic equation
are either real or complex conjugate or combination of both. Therefore their locus
must be symmetrical about the real axis.
2) As K increases from zero to infinity, each branch of the root locus originates from an
open loop pole (n nos.) with K= 0 and terminates either on an open loop zero (m
nos.) with K = along the asymptotes or on infinity (zero at ). The number of
branches terminating on infinity is equal to (n m).
3) Determine the root locus on the real axis. Root loci on the real axis are determined
by open loop poles and zeros lying on it. In constructing the root loci on the real axis
choose a test point on it. If the total number of real poles and real zeros to the right of
this point is odd, then the point lies on root locus. The complex conjugate poles and
zeros of the open loop transfer function have no effect on the location of the root loci
on the real axis.
4) Determine the asymptotes of root loci. The root loci for very large values of s must be
asymptotic to straight lines whose angles are given by
180 (2q 1)
Angle of asymptotes A
; q 0,1,2, n m - 1 (12)
n m
6) Find breakaway and break-in points. The breakaway and break-in points either lie on
the real axis or occur in complex conjugate pairs. On real axis, breakaway points
exist between two adjacent poles and break-in in points exist between two adjacent
dK
zeros. To calculate this polynomial 0 must be solved. The resulting roots are
ds
the breakaway / break-in points. The characteristic equation given by Eqn.(7), can be
rearranged as
dK d d
A B A B 0 (15)
ds ds ds
Note that the breakaway points and break-in points must be the roots of Eqn.(15),
but not all roots of Eqn.(15) are breakaway or break-in points. If the root is not on
the root locus portion of the real axis, then this root neither corresponds to
breakaway or break-in point. If the roots of Eqn.(15) are complex conjugate pair, to
ascertain that they lie on root loci, check the corresponding K value. If K is positive,
then root is a breakaway or break-in point.
7) Determine the angle of departure of the root locus from a complex pole
9) Find the points where the root loci may cross the imaginary axis. The points where
the root loci intersect the j axis can be found by
a) use of Rouths stability criterion or
b) letting s = j in the characteristic equation , equating both the real part and
imaginary part to zero, and solving for and K. The values of thus found
give the frequencies at which root loci cross the imaginary axis. The
corresponding K value is the gain at each crossing frequency.
10) The value of K corresponding to any point s on a root locus can be obtained using
the magnitude condition, or
Gain Margin
It is a factor by which the design value of the gain can be multiplied before the closed
loop system becomes unstable.
Find the point j 1 on the imaginary axis for which G j H j 1 for the design
value of K i.e. B j /A j K design .
4. Sketch the root locus of a unity negative feedback system whose forward path
K
transfer function is G(s) .
s
Solution:
Comments on stability:
The system is stable for all the values of K > 0. Th system is over damped.
K(s 2)
5. The open loop transfer function is G(s) . Sketch the root locus plot
(s 1) 2
Solution:
K= 0. One branch goes to open loop zero at s =-2.0 when K and other goes to
3) Root locus lies on negative real axis for s -2.0 as the number of open loop poles
plus number of open loop zeros to the right of s=-0.2 are odd in number.
180 (2q 1)
4) The asymptote angle is A = , q n m 1 0.
n m
Angle of asymptote is A = 180 .
(sum of poles) (sum of zeros)
5) Centroid of the asymptote is A
n m
( 1 1) ( 2)
0.0
1
(s 1) 2
K
(s 2)
dK
Break point is given by 0
ds
2(s 1)(s 2) (s 1) 2
0
(s 2) 2
s1 1, K 0; s 2 3, K 4
The root loci brakes out at the open loop poles at s=-1, when K =0 and breaks in onto
the real axis at s=-3, when K=4. One branch goes to open loop zero at s=-2 and
other goes to along the asymptotically.
7) The branches of the root locus at s=-1, -1 break at K=0 and are tangential to a line
s=-1+j0 hence depart at 90 .
8) The locus arrives at open loop zero at 180 .
9) The root locus does not cross the imaginary axis, hence there is no need to find the
imaginary axis cross over.
The root locus plot is shown in Fig.2.
System is stable for all values of K > 0. The system is over damped for K > 4. It is
critically damped at K = 0, 4.
K(s 4)
6. The open loop transfer function is G(s) . Sketch the root locus.
s(s 2)
Solution
7) Angle of departure from open loop pole at s =0 is 180 . Angle of departure from
pole at s=-2.0 is 0 .
9) The root locus does not cross the imaginary axis. Hence there is no imaginary cross
over.
Comments on stability:
K(s 0.2)
7. The open loop transfer function is G(s) . Sketch the root locus.
s 2 (s 3.6)
Solution:
2) There is one open loop zero at s = -0.2(m=1). There are three open loop poles at
s = 0, 0, -3.6(n=3). Three branches of root loci start from the three open loop poles
when K= 0 and one branch goes to open loop zero at s = -0.2 when K and other
two go to asymptotically when K .
3) Root locus lies on negative real axis between -3.6 to -0.2 as the number of open loop
poles plus open zeros to the right of any point on the real axis in this range is odd.
180 (2q 1)
4) The asymptote angle is A = , q n m 1 0,1
n m
Angle of asymptote are A = 90 , 270 .
6) The root locus does branch out, which are given by dK/ds =0.
(s 3 3.6s 2 )
K -
s 0.2
dK (3s 2 7.2s)(s 0.2) (s 3 3.6s 2 )
ds (s 0.2) 2
2s 3 4.8s 2 1.44s 0
s 0, 0.432, 1.67 and K 0, 2.55, 3.66 respectively.
The root loci brakeout at the open loop poles at s = 0, when K =0 and breakin onto
the real axis at s=-0.432, when K=2.55 One branch goes to open loop zero at s=-0.2
and other goes breaksout with the another locus starting from open loop ploe at s= -
3.6. The break point is at s=-1.67 with K=3.66. The loci go to infinity in the complex
plane with constant real part s= -1.67.
7) The branches of the root locus at s=0,0 break at K=0 and are tangential to imaginary
axis or depart at 90 . The locus departs from open loop pole at s=-3.6 at 0 .
8) The locus arrives at open loop zero at s=-0.2 at 180 .
9) The root locus does not cross the imaginary axis, hence there is no imaginary axis
cross over.
The root locus plot is shown in Fig.4.
Comments on stability:
System is stable for all values of K. System is critically damped at K= 2.55, 3.66. It is
under damped for 2.55 > K > 0 and K >3.66. It is over damped for 3.66 > K >2.55.
K
8. The open loop transfer function is G(s) . Sketch the root locus.
s(s 6s 25)
Solution:
2) There are no open loop zeros (m=0). There are three open loop poles at s=-0,
-3 j4(n=3). Three branches of root loci start from the open loop poles when K= 0
and all the three branches go asymptotically when K .
3) Entire negative real axis lies on the root locus as there is a single pole at s=0 on the
real axis.
180 (2q 1)
4) The asymptote angle is A = , q 0,1, n m 1 0, 1, 2.
n m
Angle of asymptote are A = 60 , 180 , 300 .
For a point to be break point, the corresponding value of K is a real number greater
than or equal to zero. Hence, S1,2 are not break points.
7) Angle of departure from the open loop pole at s=0 is 180 . Angle of departure from
complex pole s= -3+j4 is
p 180
(sum of the angles of vectors to a complex pole in question from other poles)
(sum of the angles of vectors to a complex pole inquestion from zeros)
4
p 180 (180 tan 1
90 ) 36.87
3
8) The root locus does cross the imaginary axis. The cross over point and the gain at
the cross over can be obtained by
Rouths criterion
s3 1 25
s2 6 K
150 K
s1
6
s0 6
For the system to be stable K < 150. At K=150 the auxillary equation is 6s2+150=0.
s = j5.
or
substitute s= j in the characteristic equation. Equate real and imaginary parts to zero.
Solve for and K.
s 3 6s 2 25s K 0
3 2
j 6 j 25 j K 0
( 6 K) j 25
2 2
0
0, j5 K 0, 150
The plot of root locus is shown in Fig.5.
Comments on stability:
System is stable for all values of 150 > K > 0. At K=150, it has sustained oscillation of
5rad/sec. The system is unstable for K >150.
9. Sketch the root locus of a unity negative feedback system whose forward path
K(s 2)
transfer function is G(s)H(s) . Comment on the
(s 1)(s 3 j)(s 3 j)
stability of the system.
2) There is one open loop zero at s = -2 (m = 1). There are three open loop poles at
s = -1, -3 j (n=3). All the three branches of root locus start from the open loop
3) Root locus lies on the negative real axis in the range s=-1 to s= -2 as there is one
pole to the right of any point s on the real axis in this range.
180 (2q 1)
4) The asymptote angle is A = , q n m 1 0,1.
n m
Angle of asymptote is A = 90 , 270 .
6) The root locus does not branch. Hence, there is no need to calculate break points.
7) The angle of departure at real pole at s=-1 is 180 . The angle of departure at the
complex pole at s=-3+j is 71.57 .
p 180
(sum of the angles of vectors to a complex pole in question from other poles)
(sum of the angles of vectors to a complex pole inquestion from zeros)
1
1 tan 1
26.57 or 153.43
-2
1 atan2(-2,1) 153.43
1 2
tan 1
-45 or 135 , 3 tan 1
90
-1 0
p
180 (153.43 90 ) 135 71.57
p
180 (206.57 270 ) 225 71.57
8) The root locus does not cross the imaginary axis. Hence there is no imaginary axis
cross over.
Comments on stability:
Solution:
2) There are no open loop zeros (m=0). There are four open loop poles (n=4) at s=0,
-0.5, -0.3 j3.1480. Four branches of root loci start from the four open loop poles
when K= 0 and go to (open loop zero at infinity) asymptotically when K .
3) Root locus lies on negative real axis between s = 0 to s = -0.5 as there is one pole to
the right of any point s on the real axis in this range.
180 (2q 1)
4) The asymptote angle is A = , q n m 1 0,1,2,3.
n m
Angle of asymptote is A = 45 , 135 , 225 , 315 .
K = -s(s+0.5)(s2+0.6s+10) = -(s4+1.1s3+10.3s2+5s)
dK
4s 3 3.3s 2 20.6s 5 0
ds
7) The angle of departure at real pole at s=0 is 180 and at s=-0.5 is 0 . The angle of
departure at the complex pole at s = -0.3 + j3.148 is -91.8
p 180
(sum of the angles of vectors to a complex pole in question from other poles)
(sum of the angles of vectors to a complex pole inquestion from zeros)
3.148
1 tan 1 84.6 or 95.4
- 0.3
3.148 6.296
2
tan 1 86.4 , 3 tan 1 90
0.2 0
p
180 (95.4 86.4 90 ) 91.8
p
180 (264.6 273.6 270 ) 91.8
8) The root locus does cross the imaginary axis, The cross over frequency and gain is
obtained from Rouths criterion.
s(s+0.5)(s2+0.6s+10)+K =0 or s4+1.1s3+10.3s2+5s+K=0
s4 1 10.3 K
3
s 1.1 5
2
s 5.75 K
28.75 - 1.1K
s1
5.75
s0 K
The system is stable if 0 < K < 26.13
Comments on stability:
System is stable for all values of 26.13 >K > 0. The system has sustained oscillation at
= 2.13 rad/sec at K=26.13. The system is unstable for K > 26.13.
K
11. The open loop transfer function is G(s) 2
. Sketch the root
s(s 4)(s 4s 20)
locus.
Solution:
6) The root locus does branch out, which are given by dK/ds =0.
The root loci brakeout at the open loop poles at s = -2.0, when K = 64 and breakin and
breakout at s=-2+j2.45, when K=100
7) The angle of departure at real pole at s=0 is 180 and at s=-4 is 0 . The angle of
departure at the complex pole at s = -2 + j4 is -90 .
p 180
(sum of the angles of vectors to a complex pole in question from other poles)
(sum of the angles of vectors to a complex pole inquestion from zeros)
4
1 tan 1
63.4 or 116.6
-2
1 atan2(4,-2) 116.6
4 8
2 tan 63.4 , 3 tan 1
1
90
2 0
p
180 - (116.6 63.4 90 ) 90
p
180 - (243.4 296.6 270 )
270 90
8) The root locus does cross the imaginary axis, The cross over point and gain at cross
over is obtained by either Rouths array or substitute s= j in the characteristic
equation and solve for and gain K by equating the real and imaginary parts to zero.
Rouths array
For the system to be stable K > 0 and 2080-8K > 0. The imaginary crossover is given
by 2080-8K=0 or K = 260.
At K = 260, the auxiliary equation is 26s2+260 = 0. The imaginary cross over occurs
at s= j 10.
or
s 4 8s 3 36s 2 80s K 0
put s j
4 3 2
j 8 j 36 j 80 j K 0
4 36 2 K j 83 80 0
Equate real and imaginary parts to zero
83 80 0 0, j 10 ; s j 10
4
36 2
K 0 K 260
The root locus plot is shown in Fig.3.
Figure 18 Root locus plot of K/s(s+4)(s2+4s+20)
Comments on stability:
Mapping of section C1
In section c1, varies from 0 to . The mapping of section C1 is given by the locus
G(j)H(j ) as is varied from 0 to . This polar plot of G(j)H(j)
G(S)H(S)= 0.05K/S(1+0.5S)(1+0.1S)
Let S=j
G(j)H(j S)= 0.05K/ jS(1+0.5 j )(1+0.1 j )
= 0.05K/-0.62+j(1-0.052)
when the locus of G(j) H(j ) crosses real axis the imaginary term will be zero and the
corresponding frequency is the phase crossover frequency
PC.
=+
C1
R
C4
=0
C3
=-
jv
=0 =
-0.00417K
Lt
S=R Rej in G(S)H(S) and varying from +/2 to -/2. Since S Rej and R, the
G(S)H(S) can be approximately as shown below [ ie (1+ST) ST ]
R (Rej)3
when =+/2 , G(S)H(S) = 0e-j3/2-----------------1a
j jv
C2 S=Plane
G(S)H(S)
plane
R
u
Section c2 in S-Plane
From the equation 1a and 1b we can study that section c2 in S-plane fig1a
is mapped as circular are of zero radius around origin in G(S)H(S) plane with argument
(phase varying from -3/2 and +3/2 as shown in fig 1b
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 a is varied from - to 0. This locus is the inverse polar plot of G(j)H(j)
The inverse polar plot is given by the mirror image of polar plot with respect to real axis.
The section C3 in S-Plane and its corresponding contour in G(S)H(S) plane shown in fig 2a
and fig 2b
jv
j
=0
S-Plane
=0
=-
=- -0.00417k
Mapping of section C4
From equation above we can say that sec c4 in S-Plane fig 3a is mapped as a circular are
of infinity radius with argument(phase) varying from +/2 to-/2 shown in figure
G(S)H(S)
jv plane
S-Plane
C4
R->
R->0
The entire Nyquist plot in G(S)H(S) plane can be obtained by combining the mappings
of individual sections shown
-0.00417K
G(S)H(S) = ----------------
S(S+2)(S+10)
STABILITY ANALYSIS
When -0.00417K =-1, the contour passes through (-1+j0) point and corresponding value
of K is the limiting value for K for stability
Limiting value of K = 1/0.00417 = 240
When K<240
When K is less than 240, the contour crosses real axis at appoint between 0 & -1+j0 . On
traveling through Nyquist plot along the indicated direction if it 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
When K is graeter than 240, the contour crosses real axis at a point B/W -1+j0 and -.
On traveling through Nyquist plot along the indicated direction it is found that the point -1+j0
is encircled in clockwise direction two time. [ since there are two clockwise encirclement and
no right half open loop poles, the closed loop system has 2 poles on right half of S-Plane.
Therefore the closed loop system is unstable
RESULT
The value of K for stability is 0<K<240
UNIT V
1. What are the two methods of designing a control system? (May 12)
The two methods of designing a control system are design using root locus and
design using bode plot.
In the previous lecture we discussed lead compensator design. In this lecture we would see
how to design a phase lag compensator
The essential feature of a lag compensator is to provide an increased low frequency gain,
thus decreasing the steady state error, without changing the transient response significantly.
For frequency response design it is convenient to use the following transfer function of a lag
compensator.
where,
The above expression is only the lag part of the compensator. The overall compensator is
Typical objective of lag compensator design is to provide an additional gain of in the low
frequency region and to leave the system with sufficient phase margin.
The frequency response of a lag compensator, with =4 and =3, is shown in Figure 1 where
Design a lag compensator so that the phase margin (PM) is at least 50 and steady state
where,
When
Thus,
Now let us modify the system transfer function by introducing K with the original system.
Thus the modified system becomes
Required PM is 50. Since the PM is achieved only by selecting K, it might be deviated from
this value when the other parameters are also designed. Thus we put a safety margin of 5
to the PM which makes the required PM to be 55.
To make g = 2.8 rad/sec, the gain crossover frequency of the modified system, magnitude
at g should be 1. Thus
Since the desired PM is already achieved with gain K, We should place = 1/ such that it
does not much effect the PM of the modified system with K. If we place 1/ one decade
below the gain crossover frequency, then
or,
With this compensator actual phase margin of the system becomes 52.7, as shown in
Figure 2, which meets the design criteria.
When a single lead or lag compensator cannot guarantee the specified design criteria, a lag-
lead compensator is used.
In lag-lead compensator the lag part precedes the lead part. A continuous time lag-lead
compensator is given by
where,
In a nutshell,
If it is not specified which type of compensator has to be designed, one should first
check the PM and BW of the uncompensated system with adjustable gain K.
If the BW is smaller than the acceptable BW one may go for lead compensator. If
the BW is large, lead compensator may not be useful since it provides high frequency
amplification.
One may go for a lag compensator when BW is large provided the open loop
system is stable.
If the lag compensator results in a too low BW (slow speed of response), a lag-lead
compensator may be used.
6. Lag-lead compensator design Consider the following system with transfer
function
Design a lag-lead compensator C(s) such that the phase margin of the compensated system
is at least 45 at gain crossover frequency around 10 rad/sec and the velocity error constant
Kv is 30.
where,
When
Thus K = 30 . Bode plot of the modified system KG(s) is shown in Figure 2. The gain
crossover frequency and phase margin of KG(s) are found out to be 9.77 rad/sec and -17.2
respectively.
However, we know that introduction of a lead compensator will eventually increase the gain
crossover frequency to maintain the low frequency gain.
Thus the gain crossover frequency of the system cascaded with a lead compensator is likely
to be much above the specified one, since the gain crossover frequency of the
uncompensated system with K is already 9.77 rad/sec.
From Figure 2, it is seen that at 10 rad/sec the phase angle of the system is -198.
Since the new g should be 10 rad/sec, the required additional phase at g, to maintain the
specified PM, is 45 - (180 - 198) = 63 . With safety margin 2,
And
which gives . However, introducing this compensator will actually increase the
gain crossover frequency where the phase characteristic will be different than the designed
one. This can be seen from Figure 3.
Figure 3: Frequency response of the system in Example 1 with only a lead compensator
The gain crossover frequency is increased to 23.2 rad/sec. At 10 rad/sec, the phase angle is
-134 and gain is 12.6 dB. To make this as the actual gain crossover frequency, lag part
should provide an attenuation of -12.6 dB at high frequencies.
which gives . Now, should be placed much below the new gain crossover
The frequency response of the system after introducing the above compensator is shown in
Figure 4, which shows that the desired performance criteria are met.
7. Problem Find The stability of given system
8. Find the stability of given System Using Bode plot.
9. Design compensator for given System.
10. Design a Lag Compensator of Given System.