Control System Final
Control System Final
Syllabus
1.0 Introduction to Control system
1.1 Scope of Control System Engineer
1.2 Classification of Control System
1.3 Historical development of Control system
1.4 Analogues systems
1.5 Transfer function of Systems
1.6 Block diagram representation
1.7 Signal Flow Graph(SFG)
2.0 Feedback Characteristics of Control systems and sensitivity measures
IN
2.1 The Concept of Feedback and Closed loop control
2.2 Merits of using Feedback control system
2.3 Regenerative Feedback
3.0Control System Components
3.1 Potentiometers
3.2 DC and AC Servomotors
GA
3.3 Tachometers
3.4 Amplidyne
&
3.5 Hydralulic systems
3.6 Pneumatic systems
3.7 Stepper Motors
M
IN
7.3 Polar Plots and inverse Polar plots
7.4 Bode Diagrams
7.4.1 Principal factors of Transfer function
GA
7.4.2 Procedure for manual plotting of Bode Diagram
7.4.3 Relative stability Analysis
7.4.4 Minimum Phase, Non-minimum phase and All pass systems
7.5 Log Magnitude vs Phase plots.
&
7.6 Nyquist Criterion
7.6.1 Mapping Contour and Principle of Argument
7.6.2 Nyquist path and Nyquist Plot
M
MODULE#1
IN
GA
&
A M
BE
CHAPTER#1
1. Basic Concept of Control System
Control Engineering is concerned with techniques that are used to solve the following six
problems in the most efficient manner possible.
(a)The identification problem :to measure the variables and convert data for analysis.
IN
(f)The optimization problem: from a variety of design to choose the best.
The two basic approaches to solve these six problems are conventional and modern approach. The
electrical oriented conventional approach is based on complex function theory. The modern
GA
approach has mechanical orientation and based on the state variable theory.
Therefore, control engineering is not limited to any engineering discipline but is equally
applicable to aeronautical, chemical, mechanical, environmental, civil and electrical engineering.
For example, a control system often includes electrical, mechanical and chemical components.
Furthermore, as the understanding of the dynamics of business, social and political systems
&
increases; the ability to control these systems will also increase.
M
Actuator: It is the device that causes the process to provide the output. It is the device that
provides the motive power to the process.
Design: The process of conceiving or inventing the forms, parts, and details of system to
achieve a specified purpose.
Simulation: A model of a system that is used to investigate the behavior of a system by
utilizing actual input signals.
Feedback Signal: A measure of the output of the system used for feedback to control the
system.
Negative feedback: The output signal is feedback so that it subtracts from the input signal.
Block diagrams: Unidirectional, operational blocks that represent the transfer functions of
the elements of the system.
IN
Signal Flow Graph (SFG): A diagram that consists of nodes connected by several directed
branches and that is a graphical representation of a set of linear relations.
Specifications: Statements that explicitly state what the device or product is to be and to do.
GA
It is also defined as a set of prescribed performance criteria.
Open-loop control system: A system that utilizes a device to control the process without
using feedback. Thus the output has no effect upon the signal to the process.
Closed-loop feedback control system: A system that uses a measurement of the output and
&
compares it with the desired output.
Regulator: The control system where the desired values of the controlled outputs are more or
less fixed and the main problem is to reject disturbance effects.
M
Servo system: The control system where the outputs are mechanical quantities like
acceleration, velocity or position.
A
Stability: It is a notion that describes whether the system will be able to follow the input
command. In a non-rigorous sense, a system is said to be unstable if its output is out of
control or increases without bound.
BE
Multivariable Control System: A system with more than one input variable or more than
one output variable.
Trade-off: The result of making a judgment about how much compromise must be made
between conflicting criteria.
1.2. Classification
1.2.1. Natural control system and Man-made control system:
Natural control system: It is a control system that is created by nature, i.e. solar
system, digestive system of any animal, etc.
Man-made control system: It is a control system that is created by humans, i.e.
automobile, power plants etc.
1.2.2. Automatic control system and Combinational control system:
Automatic control system: It is a control system that is made by using basic theories
from mathematics and engineering. This system mainly has sensors, actuators and
responders.
Combinational control system: It is a control system that is a combination of natural
and man-made control systems, i.e. driving a car etc.
1.2.3. Time-variant control system and Time-invariant control system:
Time-variant control system: It is a control system where any one or more
parameters of the control system vary with time i.e. driving a vehicle.
Time-invariant control system: It is a control system where none of its parameters
vary with time i.e. control system made up of inductors, capacitors and resistors only.
1.2.4. Linear control system and Non-linear control system:
Linear control system: It is a control system that satisfies properties of homogeneity
and additive.
IN
Homogeneous property: f x y f x f y
Additive property: f x f x
Non-linear control system: It is a control system that does not satisfy properties of
1.2.5.
GA
homogeneity and additive, i.e. f x x3
Closed-loop control system:It is a control system where its control action depends
on both of its input signal and output response.
IN
Advantages:
Simple design and easy to construct
Economical
Easy for maintenance
Highly stable operation
Dis-advantages:
GA
Not accurate and reliable when input or system parameters are variable in
nature
Recalibration of the parameters are required time to time
&
1.3.2. Closed-loop control system:
It is a control system where its control action depends on both of its input signal and
output response as shown in Fig.1.2.
A M
BE
Examples: automatic electric iron, missile launcher, speed control of DC motor, etc.
Advantages:
More accurate operation than that of open-loop control system
Can operate efficiently when input or system parameters are variable in
nature
Less nonlinearity effect of these systems on output response
High bandwidth of operation
There is facility of automation
Time to time recalibration of the parameters are not required
Dis-advantages:
Complex design and difficult to construct
IN
4 input, provided the system parameters constant input, depending upon the
remain unaltered feedback
System output variation due to variation in
System output variation due to variation in
5 parameters of the system is greater and the
parameters of the system is less.
output very in an uncontrolled way
6
7
8
9
Small bandwidth
More stable
Affected by non-linearities
GA
Error detection is not present Error detection is present
Large bandwidth
Less stable or prone to instability
Not affected by non-linearities
10 Very sensitive in nature Less sensitive to disturbances
11 Simple design Complex design
12 Cheap Costly
&
A M
BE
1.4. Servomechanism
It is the feedback unit used in a control system. In this system, the control variable is
a mechanical signal such as position, velocity or acceleration. Here, the output signal
is directly fed to the comparator as the feedback signal, b(t) of the closed-loop control
system. This type of system is used where both the command and output signals are
mechanical in nature. A position control system as shown in Fig.1.3 is a simple
example of this type mechanism. The block diagram of the servomechanism of an
automatic steering system is shown in Fig.1.4.
IN
GA
&
M
Examples:
Missile launcher
Machine tool position control
Power steering for an automobile
Roll stabilization in ships, etc.
1.5. Regulators
It is also a feedback unit used in a control system like servomechanism. But, the
output is kept constant at its desired value. The schematic diagram of a regulating
IN
GA
Fig.1.5. Schematic diagram of a regulating system
&
M
Examples:
Temperature regulator
BE
Speed governor
Frequency regulators, etc.
CHAPTER#2
2. Control System Dynamics
2.1. Definition: It is the study of characteristics behaviour of dynamic system, i.e.
(a) Differential equation
i. First-order systems
ii. Second-order systems
(b) System transfer function: Laplace transform
2.2. Laplace Transform: Laplace transforms convert differential equations into algebraic
equations. They are related to frequency response.
L x t X s x(t )e st dt (2.1)
0
IN
L x t X s x(t )e st dt (2.2)
0
No. Function
Time-domain
x(t)=
Laplace domain
X(s)= {x(t)}
GA 7 Sine
s2 2
-1
{X(s)} 8 Cosine s
2 2
1 Delay - e - s
&
2 Unit impulse 1 9 Hyperbolic
sine
s2 2
1
3 Unit step u(t)
M
s 10 Hyperbolic s
cosine
1 s2 2
4 Ramp t
s2
A
11 Exponentiall e t
sin t
Exponential 1 y decaying (s )2 2
-
5 e sine wave
decay s
BE
12 Exponentiall e t
cos t s
Exponential y decaying
6 1 e t (s )2 2
approach s(s ) cosine wave
2.3. Solution of system dynamics in Laplace form: Laplace transforms can be solved using
partial fraction method.
A system is usually represented by following dynamic equation.
A s
N s (2.3)
B s
Example 2.1:
Expand the following equation of Laplacetransform in terms of its partial fractionsand obtain
its time-domain response.
2s
IN
Y (s )
( s 1)(s 2)
Solution:
The following equation in Laplacetransform is expandedwith its partial fractions as follows.
2s A B
( s 1)( s 2) ( s 1) ( s 2)
2s A( s 2) B( s 1)
GA
( s 1)( s 2) ( s 1)( s 2)
By equating both sides, A and B are determined as A 2, B 4 . Therefore,
2 4
&
Y ( s)
(s 1) (s 2)
Taking Laplace inverse of above equation,
t 2t
y (t ) 2e 4e
M
Example 2.2:
Expand the following equation of Laplacetransform in terms of its partial fractionsand obtain
its time-domain response.
2s
Y ( s)
(s 1) 2 (s 2)
Solution:
The following equation in Laplacetransform is expandedwith its partial fractions as follows.
2s A B C
( s 1) 2 ( s 2) ( s 1) 2 (s 1) (s 2)
By equating both sides, A and B are determined as A 2, B 4 . Therefore,
2 4 4
Y ( s)
(s 1) 2 (s 1) ( s 2)
Taking Laplace inverse of above equation,
t t 2t
y( t ) 2te 4e 4e
Example 2.3:
Expand the following equation of Laplacetransform in terms of its partial fractionsand obtain
its time-domain response.
2s 1
Y ( s)
( s 1 j)( s 1 j )
Solution:
The following equation in Laplacetransform is expandedwith its partial fractions as follows.
IN
2s 1
Y ( s)
(s 1) 1 (s 1)2 1
2
(s 1 j )(s 1 j )
Solution:
Solution of above equation,
A
lim ( s
s 1 j )( s 1 j)
2
Example 2.5:
Determine the initial value of the time-domain response of the following equation using the
initial-value theorem.
2s
Y ( s)
(s 1) 2 (s 2)
Solution:
Solution of above equation,
t t 2t
y( t ) 2te 4e 4e
Applying final value theorem,
s (2s 1)
lim ( s
s 1 j )( s 1 j)
2
IN
GA
&
A M
BE
CHAPTER#3
3. Transfer Function
3.1. Definition: It is the ratio of Laplace transform of output signal to Laplace transform of input
signal assuming all the initial conditions to be zero, i.e.
Let, there is a given system with input r(t) and output c(t) as shown in Fig.3.1 (a), then its
Laplace domain is shown in Fig.3.1 (b). Here, input and output are R(s) and C(s) respectively.
(a) (b)
IN
(c)
GA
Fig.3.1. (a) A system in time domain, (b) a system in frequency domainand (c) transfer function with differential
operator
G(s) is the transfer function of the system. It can be mathematically represented as follows.
C s
&
G s Equation Section (Next)(3.1)
R s zero initial condition
Example 3.1: Determine the transfer function of the system shown inFig.3.2.
A M
BE
Solution:
Fig.3.1 is redrawn in frequency domain as shown in Fig.3.2.
1
Vi s R Ls I s (3.2)
Cs
1
Vo s I s (3.3)
Cs
From eq (2.12),
1
I s Vo s / CsVo s (3.4)
Cs
IN
1
Vi s R Ls CsVo s
Cs
Vo s 1 1 (3.5)
1 2
Vi s GA LCs RCs 1
R Ls Cs
Cs
1
G s 2
(3.6)
LCs RCs 1
&
3.2. General Form of Transfer Function
m
s zi
K s z1 s z2 ... s zm
M
i 1
G s K n
(3.7)
s p1 s p2 ... s pn
s zj
i 1
A
Where, z1 , z2 ...zm are called zeros and p1 , p2 ... pn are called poles.
BE
Example 3.2:
Obtain the pole-zero map of the following transfer function.
( s 2)( s 2 j 4)( s 2 j 4)
G (s )
( s 3)( s 4)( s 5)( s 1 j 5)( s 1 j 5)
Solution:
The following equation in Laplacetransform is expandedwith its partial fractions as follows.
Zeros Poles
s=2 s=3
s=-2-j4 s=4
s=-2+j4 s=5
s=-1-j5
s=-1+j5
IN
GA
Fig.3.3. pole-zero map
dt
Poles and zeros can be obtained from the transfer function
Stability can be known
Impulse response can be found
3.5. Disadvantages of Transfer function:
Applicable only to linear system
Not applicable if initial condition cannot be neglected
It gives no information about the actual structure of a physical system
CHAPTER#4
4. Description of physical system
4.1. Components of a mechanical system: Mechanical systems are of two types, i.e. (i)
translational mechanical system and (ii) rotational mechanical system.
4.1.1. Translational mechanical system
There are three basic elements in a translational mechanical system, i.e. (a) mass, (b)
spring and (c) damper.
(a) Mass: A mass is denoted by M. If a force f is applied on it and it displays
d 2x
distance x, then f M 2 as shown in Fig.4.1.
dt
IN
GA
Fig.4.1. Force applied on a mass with displacement in one direction
If a force f is applied on a massM and it displays distance x1in the direction of f and
d 2 x1 d 2 x2
&
distance x2 in the opposite direction, then f M as shown in Fig.4.2.
dt 2 dt 2
X2 X1
M
M f
A
If a force f is applied on a springK and it displays distance x1in the direction of f and
distance x2 in the opposite direction, then f K x1 x2 as shown in Fig.4.4.
IN
Fig.4.5. Force applied on a damper with displacement in one direction
GA
If a force f is applied on a damperD and it displays distance x1in the direction of f and
dx1 dx2
distance x2 in the opposite direction, then f D as shown in Fig.4.6.
dt dt
&
M
There are three basic elements in a Rotational mechanical system, i.e. (a) inertia, (b)
spring and (c) damper.
BE
IN
GA
Fig.4.7. Current and voltage shown in resistor, inductor and capacitor
&
4.2.2. Current source electrical system:
v
If i is the current through a resistor and v is the voltage drop in it, then i .
R
M
If i is the current through an inductor and v is the voltage developed in it, then
1
i vdt .
L
A
If i is the current through a capacitor and v is the voltage developed in it, then
dv
i C .
dt
BE
Fig.4.8.
Solution
Voltage across resistance, eR (t ) i (t ) R
1
Voltage across capacitance, eC (t ) i(t ) dt
C
1
Total voltage drop, ei eR eC i (t ) R i(t ) dt
C
1
Laplace transform of above equation, Ei ( s) I ( s) R
Cs
System transfer function betweenvoltage drop across the capacitanceand input
EC (s ) 1 1
voltage,
IN
Ei ( s ) RCs 1 s 1
where, RC is the time-constant
GA
Q.4.2. Find system transfer function betweenfunction between the inductance
currentto the source currentin the followingRL circuit as shown in Fig.4.9.
&
M
Fig.4.9.
A
BE
e( t )
Voltage across the Resistance, e (t ) iR R iR
R
diL 1
Voltage across the Inductance, e(t ) L iL e(t ) dt
dt L
e (t ) 1
Total current, ia iR iL e(t ) dt
R L
Laplace transform of the current source,
1 1 E
Ia (s) E (s) and I L ( s )
R Ls Ls
Transfer function between the inductance current to the source current,
I L (s ) 1 1
I a (s) L s 1
s 1
R
L
where is the time-constant
R
IN
Fig.4.10.
2
EC ( s ) 1 n
E (s)
Cs R Ls
1 s2 2 n s 2
n
Cs
A
1 R
where n and
LC L
2
BE
C
[Link] the transfer function of the following Spring-mass-damperas shown
in Fig.4.11.
Fig.4.11.
Solution
X (s ) 1 1
2 2 2
F ( s) ms cs k m s 2 n s n
4.3. Analogous system: Fig.4.12 shows a translational mechanical system, a rotational control
system and a voltage-source electrical system.
IN
(a)
GA
(b)
&
A M
BE
(c)
Fig.4.12. (a) a voltage-source electrical system,(b) a translational mechanical system and (c) a rotational control
system
d 2q dq 1
L 2
R q v t
dt dt C
d2 d
J 2 D K T Equation Chapter 8 Section 0(4.1)
dt dt
d2x dx
M 2 D Kx f
dt dt
Where,
q idt (4.2)
The solutions for all the above three equations given by eq (4.2) are same. Therefore, the
above shown three figures are analogous to each other. There are two important types of
analogous systems, i.e. force-voltage (f-v) analogy and force-current analogy. From eq (4.2),
f-v analogy can be drawn as follows.
IN
Similarly, f-i analogy that can be obtainedfrom eq (4.1), can be drawn as follows.
Ra=resistance of armature
ia=armature current
BE
If=field current
Eb=back emf
f=equivalent viscous friction coefficient of motor and load referred to motor shaft
Fig.4.6. Schematic diagram of armature control type speed control system of a DC motor
IN
The air-gap flux is proportional of the field current i.e.
KfIf (4.3)
GA
The torque Tm developed by the motor is proportional to the product of armature current and
air gap flux i.e.
Tm = k1 K f I f ia (4.4)
&
In armature-controlled D.C. motor,the field current is kept constant,so that eq(4.4) can be
written as follows.
Tm =K t ia (4.5)
M
d
Eb =K b (4.6)
dt
A
dia
La Ra ia Eb Ea (4.7)
dt
d2 d
J f Tm Kt Ia (4.8)
dt 2 dt
Taking the Laplace transforms of equations (4.6), (4.7) and (4.8), assuming zero initial
conditions, we get
Eb s =sK b s (4.9)
(s 2 J sf ) ( s ) Tm ( s ) Kt Ia (4.11)
From eq(4.9) to (4.11) the transfer function of the system is obtained as,
(s ) Kt
G s (4.12)
Ea ( s ) s Ra sLa sJ f Kt K b
Kt
(s) Ra sLa sJ f 1
G s (4.13)
Ea ( s ) Kt Kb s
1
Ra sLa sJ f
1 1 1 s
IN
sLa
R a
sJ f s s
GA
Fig.4.7. Block diagram of armature control type speed control system of a DC motor
(s ) Kt Kt Kb
s2 J s f (4.14)
Ea ( s) Ra Ra
M
Kt Kb
The term f indicates that the back emf of the motor effectively increases the
Ra
viscous friction of the system. Let,
A
K t Kb
f f (4.15)
BE
Ra
Where f be the effective viscous friction coefficient. The transfer function given by eq(4.15)
may be written in the following form.
s Km
(4.16)
Ea s s s 1
Kt J
Here K m = = motor gain constant, and = motor time [Link], the motor
Ra f f
torque and back emf constant Kt, Kb are interrelated.
4.5. Mathematical model of field controlled DC motor: The field control type speed control
system of a DC motor is shown in Fig.4.8. The following components are used in this system.
Rf=Field winding resistance
If=field current
f=equivalent viscous friction coefficient of motor and load referred to motor shaft
IN
Ia
Rf (constant)
Lf Tm
ef
If GA M
&
Fig.4.8. Block diagram of field control type speed control system of a DC motor
M
In field control motor the armature current is fed from a constant current [Link] air-gap
.
A
Kf If (4.17)
The torque Tm developed by the motor is proportional to the product of armature current and
BE
Tm =k1K f I f I a Kt I f (4.18)
dI f
Lf Rf I f Ef (4.19)
dt
The torque equation is
d2 d
J f Tm Kt I f (4.20)
dt 2 dt
Taking the Laplace transforms of equations (4.19) and (4.20) assuming zero initial conditions,
we get the following equations
Lf s Rf I f s Ef s (4.21)
and
Js 2 fs s Tm s Kt I f s (4.22)
From eq(4.21) and (4.22) the transfer function of the system is obtained as
s Kt
G s (4.23)
Ef s s Rf sL f Js f
The transfer function given by eq(4.23) may be written in the following form.
s Kt Km
(4.24)
IN
Ea s s Lf s Rf Js f s s 1 s 1
Kt Lf J
Here K m = motor gain constant, and = time constant of field circuit and
Rf f Rf f
GA
= mechanical time [Link] small size motors field control is [Link] block
diagram that is constructed from eq (4.24) is shown in Fig.4.9.
1 K
&
t
Ef(s) (s)
sL f
R f
s sJ f
Fig.4.9. Block diagram of field control type speed control system of a DC motor
A M
BE
CHAPTER#5
5. Block Diagram Algebra
5.1. Basic Definition in Block Diagram model:
Block diagram: It is the pictorial representation of the cause-and-response relationship
between input and output of a physical system.
(a) (b)
Fig.5.1. (a) A block diagram representation of a system and (b) A block diagram representation with gain of a
system
IN
Output: The value of input multiplied by the gain of the system.
C s G s R s (5.1)
GA
Summing point: It is the component of a block diagram model at which two or more signals
can be added or subtracted. In Fig.15, inputs R(s) and B(s) have been given to a summing
point and its output signal is E(s). Here,
E s R s B s (5.2)
&
A M
BE
Take-off point: It is the component of a block diagram model at which a signal can be taken
directly and supplied to one or more points as shown in Fig.5.2.
Forward path: It is the direction of signal flow from input towards output.
Feedback path: It is the direction of signal flow from output towards input.
5.2. Developing Block Diagram model from mathematical model:
Let’s discuss this concept with the following example.
Example: A system is described by following mathematical equations. Find its corresponding
block diagram model.
x1 3x1 2 x2 5 x3 (5.3)
x2 x1 4 x2 3 x3 (5.4)
x3 2 x1 x2 x3 (5.5)
Example: Eq (5.3), (5.4) and (5.5) are combiningly results in the following block diagram
model.
x3(s)
5
x 1(s) + x s x 1(s)
3 ++ 1
1/s
x2(s)
2
x2(s) x s
++ 2 x2(s)
4 1/s
+
IN
x3(s)
3
GA x3(s)
x2(s) + x s
++ 3
1/s
x3(s)
x1(s)
&
2
M
1 Rule 1 Cascade
2 Rule 2 Parallel
G(s)
3 Rule 3
1 G s H s
Loop
IN
Associative
4 Rule 4
Law
GA
Move take-
5 Rule 5 off point
after a block
&
Move take-
off point
6 Rule 6
M
before a
block
A
Move
summing-
7 Rule 7
BE
point point
after a block
Move
summing-
8 Rule 8 point point
before a
block
Move take-
off point
9 Rule 9 after a
summing-
point
Move take-
off point
10 Rule 10 before a
summing-
point
IN
Fig.5.4. Rules for reduction of Block Diagram model
5.5. Procedure for finding output of Block Diagram model with multiple inputs:
Step 1: Consider one input taking rest of the inputs zero, find output using the procedure
described in section 4.3.
A
Step 2: Follow step 1 for each inputs of the given Block Diagram model and find their
corresponding outputs.
BE
CHAPTER#6
6. Signal Flow Graphs (SFGs)
It is a pictorial representation of a system that graphically displays the signal transmission in it.
IN
in between them. In Fig.6.1, there are two forward paths, i.e. path-1:‘r-x1-x2-x3-x4-x5-x6-c’ and
path-2:‘r-x1-x3-x4-x5-x6 -c’.
Feedback path: It is a path from output node or a node near output node to a node near input
node without repeating any of the nodes in between them (Fig.6.1).
GA
Loop: It is a closed path that starts from one node and reaches the same node after trading
through other nodes. In Fig.6.1, there are four loops, i.e. loop-1:‘x2 -x3-x4-x1’, loop-2:‘x5-x6-
x5’, loop-3:‘x1-x2-x3-x4-x5-x6 -x1’ and loop-4:‘x1-x3-x4-x5 -x6-x1’.
Self Loop: It is a loop that starts from one node and reaches the same node without trading
through other nodes i.e. loop in node ‘x4’ with transmittance ‘t55’ in Fig.6.1.
&
Path gain: It is the product of gains or transmittances of all branches of a forward path. In
Fig.6.1, the path gains are P1 = t1t2t3 t4t5 (for path-1) and P2 = t9t3t4t5 (for path-2).
Loop gain: It is the product of gains or transmittances of all branches of a loop In Fig.6.1,
there are four loops, i.e. L1 = -t2t3t6, L 2 = -t5t7, L 3 = -t1t2t3t4t5t8, and L4 = -t9t3t4t5t8 .
M
Dummy node: If the first node is not an input node and/or the last node is not an output node
than a node is connected before the existing first node and a node is connected after the
existing last node with unity transmittances. These nodes are called dummy nodes. In Fig.6.1,
A
Example:
IN
+ 1 1 (s) 1 s
Ea(s) KT
sL R sJ f s s
- a a
E b(s)
GA Kb
Step-3: Each ofgains is replaced by transmittances of the branches connected between two nodes
A
Step-4: Each ofgains is replaced by transmittances multiplied with (-1) of the branches connected
BE
1 1 1 s
sL a
R a
sJ f s s
(a)
1 1 1
sL a
R
a
sJ f s
(b)
IN
6.4. Mason’s gain formula:
Transfer function of a system=
G s
GA
C s
N
k 1
Pk k
(6.1)
R s
Where,
&
N= total number of forward paths
=1-( possible
two non-touching loops) - ( -touching
loops) + …
A
th
k= value of forward path
BE
Example:
Find the overall transfer function of the system given in Fig.6.1 using Mason’s gain
formula.
Solution:
In Fig.6.1,
Loop gain of individual loops: L1 t 2t3t6 , L2 t5t7 , L3 t1t2 t3t4 t5t8 and L4 t9t3t 4t5t8
1 L1 L2 L3 L4 L1 L2 0 1 L1 L2 L3 L4 L1L2
1 1 0 1and 2 1 0 1
P1 1 P2 2
G s
t1t2 t3 t4 t5 t6 t3 t4 t5
G s
1 t 2 t 3t 6 t5 t7 t1t2 t3t 4t5 t8 t9 t3t 4t5 t8 t 2 t3 t5 t6 t7
IN
GA
&
A M
BE
CHAPTER#7
7. Feedback Characteristics of Control System
7.1. Feedback and Non-feedback Control systems
Non-feedback control system: It is a control system that does not have any feedback paths.
It is also known as open-loop control system. It is shown in Fig.7.1 (a) and (b).
Feedback control system: It is a control system that has at least one feedback path. It is also
known as closed-loop control system. It is shown in Fig.7.2 (a) and (b).
(a)
(b)
IN
Fig.7.1. (a) Block diagram of a non-feedback control system and (b) SFG of a non-feedback control system
(a)
GA (b)
Fig.7.2. (a) Block diagram of a feedback control system and (b) SFG of a feedback control system
&
E s R s B s (7.1)
and
BE
G s
T1 s (7.2)
1 G s H s
7.2.2. Regenerative feedback control system: It is a control system where the feedback
signal supports or adds the input signal. Here,
Error or actuating signal = (Input signal) + (Feedback signal).
Referring Fig.7.4,
E s R s B s (7.3)
and
G s
T2 s (7.4)
1 G s H s
IN
GA
Fig.7.4. Block diagram of a regenerative feedback control system
7.3. Effect of parameter variation on overall gain of a degenerative Feedback Control system
&
The overall gain or transfer function of a degenerative feedback control system depends upon
these parameters i.e. (i) variation in parameters of plant, and (ii) variation in parameter of
feedback system and (ii) disturbance signals.
The term sensitivity is a measure of the effectiveness of feedback on reducing the influence of
M
any of the above described parameters. For an example, it is used to describe the relative
variations in the overall Transfer function of a system T(s) due to variation in G(s).
A
( )
=
( )
BE
C s G s R s
C s C s G s G s R s
C s C s G s R s G s R s
C s G s R s (7.5)
In an closed-loop system,
G s
C s R s
1 G s H s
G s G s
C s C s R s
1 G s G s H s
G s G s
C s C s R s
1 G s H s G s H s
IN
G s G s
C s C s R s
1 G s H s
G s G s
Or
C s C s GA
1 G s H s
R s
1 G s H s
R s
G s
C s R s (7.6)
&
1 G s H s
This concept can be reproved using sensitivity. Sensitivity on T(s) due to variation in G(s) is
M
given by
T T T G
SGT (7.7)
A
GG G T
T T G G
SGT 1 (7.8)
GG G G
T T 1 GH GH G 1
SGT 2
(7.9)
GG 1 GH G 1 GH 1 GH
T T T H
SHT (7.10)
H H H T
T H G H GH
S HT G 2
(7.11)
H T 1 GH G 1 GH 1 GH
For higher value of GH, sensitivity approaches unity. Therefore, change in H affects
directly the system output.
IN
Equation Chapter (Next) Section 1
GA
&
A M
BE
IN
MODULE#2
GA
&
A M
BE
CHAPTER#8
IN
GA
&
M
i.e.
Laplace transforms convert differential equations into algebraic equations. They are related to
frequency response
L x t X s x (t )e st dt (8.1)
0
1 Delay - e-
2 Unit impulse 1
4 Ramp t 1
s2
IN
5 Exponential decay e- 1
s
6 Exponential approach 1 e t
7 Sine
GA s (s )
s2 2
8 Cosine s
&
2 2
s
9 Hyperbolic sine
s2 2
10 Hyperbolic cosine s
M
2 2
s
11 Exponentially decaying sine e t
sin t
wave (s )2 2
A
K ( s z1 )( s z2 ) ( s zm )
C ( s) G ( s) R( s) R(s) (8.2)
( s p1 )( s p2 ) ( s pn )
s 1 (8.4)
IN
Dirac delta function
x (t ) xi (t a ) (8.5)
GA
Integral property of Dirac delta function
&
(t ) (t to )dt (to ) (8.6)
st sa
X (s ) e xi (t a) dt xi e (8.7)
0
0 ;t 0
u t (8.8)
BE
K ;t 0
IN
GA
Fig.7.3. Ramp signal
&
8.4.4. Parabolic Signal A step signal a(t) is mathematically defined as follows.
0
;t 0
a t Kt 2 (8.12)
M
;t 0
2
K
A s (8.13)
s3
BE
x t sin t (8.14)
st
X s e sin t dt (8.15)
0
s2 2
IN
GA
Fig.7.4. Sinusoidal signal
Here,
BE
E s R s B s (8.16)
B s C s H s (8.17)
C s E s G s (8.18)
1 G s H s E s R s (8.21)
R s
E s (8.22)
1 G s H s
Steady-state error,
ess lim e t lim sE s (8.23)
t s 0
IN
8.6. Types of input and Steady-state error:
8.6.1. Step Input
GA A
R s (8.25)
s
Using (18) in (17),
A
s
s A
ess lim lim (8.26)
&
s 01 G s H s s 01 G s H s
A A
ess (8.27)
1 lim G s H s 1 KP
s 0
M
Where,
KP lim G s H s (8.28)
A
s 0
BE
Where,
KV lim sG s H s (8.31)
s 0
IN
A A
ess
lim s 2G s H s KA
s 0
Where,
KA
s
GA
lim s 2G s H s
0
(8.34)
s 0 ess
KA
BE
Here,
KP lim G s H s K (8.37)
s 0
Therefore,
A
ess (8.38)
1 K
[Link]. Type 1
K 1 T1 s 1 T2 s ... 1 Tn s
G s H s (8.39)
s 1 Ta s 1 Tb s ... 1 Tm s
Here,
KV lim sG s H s K (8.40)
s 0
IN
Therefore,
A
ess (8.41)
K
[Link]. Type 2
GA
K 1 T1s 1 T2 s ... 1 Tn s
G s H s (8.42)
s 2 1 Ta s 1 Tb s ... 1 Tm s
Here,
&
KA lim s 2G s H s K (8.43)
s 0
Therefore,
M
A
ess (8.44)
K
A
Steady-state error and error constant for different types of input are summarized as follows.
BE
IN
GA
&
A M
BE
1
E s R s
1 G s H s
So,
E s F1 s F2 s (8.45)
1
Where, F1 and F2 s R s
1 G s H s
IN
0 0
t t t 2 t 3
e t f1 r t d r t f1 d r t f1 d r t f1 d ... (8.48)
2! 3!
&
0 0 0 0
Therefore,
t t t 2 t 3
ess lim e t lim f1 r t d r t f1 d r t f1 d r t f1 d ...
A
t t
0 0 0
2! 0
3!
2 3
BE
ess f1 r t d r t f1 d r t f1 d r t f1 d ...
0 0 0
2! 0
3!
(8.50)
Eq (44) can be rewritten as
C2 C3
ess C0 r t C1r t r t r t ... (8.51)
2! 3!
Where, C0 , C1, C2, C3, etc. are dynamic error coefficients. These are given as
C0 f1 d lim F1 s
s 0
0
dF1 s
C1 f1 d lim
0
s 0 ds
, and so on… (8.52)
2
d 2 F1 s
C2 f1 d lim
0
2! s 0 ds 2
3
d 3 F1 s
C3 f1 d lim
0
3! s 0 ds 3
IN
y y Kx (t ) (8.53)
Normalized response
In this type of response
Static components are taken out leaving only the dynamic component
The dynamic components converge to the same value for different physical systems of
the same type or order
Helps in recognizing typical factors of a system
Kxi Kxi 1
Y (s ) (8.56)
(1 s) 1
s
IN
Time-domain response
t
Kxi
y (t ) e (8.57)
The above equation is called Duhamel’s integral. Normalized response of a first-order system to
M
y (t )
Kxi
t/
Time-domain response
t
y (t ) Kxi 1 e (8.62)
IN
y (t )
Kxi GA
&
t/
M
Y (s) (8.64)
s 2 (1 s ) s2 s s
1
Time-domain response
t
y (t )
t e (8.65)
K
Normalized response of a first-order system to impulse input is shown below.
y (t )
Kxi
t/
IN
8.8.4. Sinusoidal input to a first-order system
Governing differential equation
K A s 1
Y ( s) 2 2 2 2 2 2 2
(8.67)
(1 s) s 1 s 1/ s s
&
Time-domain response
y (t ) t/ 1
2
e cos t sin t (8.68)
KA 1
M
my cy ky Kx (t ) (8.69)
Y (s) K
2 2
(8.70)
X ( s) m s 2 n s n
IN
GA
Pole-zero map
2
s1,2 n 1 (8.71)
M
s1,2 n (8.72)
2
s1,2 n j 1
(8.73)
IN
s1,2 n j d
Here, tan
2
1
(d) = 0 un-damped
Poles are:
s1,2 j n (8.75)
Solved problems:
1. A single degree of freedom spring-mass-damper system has the following data: spring stiffness 20
kN/m; mass 0.05 kg; damping coefficient 20 N-s/m. Determine
IN
(a) undamped natural frequency in rad/s and Hz
(b) damping factor
(c) damped natural frequency n rad/s and Hz.
If the above system is given an initial displacement of 0.1 m, trace the phasor of the system for three
cycles of free vibration.
Solution:
k 20 103
GA
n 632.46 rad/s
m 0.05
632.46
&
n
fn 100.66 Hz
2 2
c 20
0.32
2 km 2 20 103 0.05
M
2
d n 1 632.46 1 0.322 600 rad/s
A
d 600
fd 95.37 Hz
2 2
BE
nt 0.32 632.46 t
y (t ) Ae 0.1e
2. A second-order system has a damping factor of 0.3 (underdamped system) and an un-damped
natural frequency of 10 rad/s. Keeping the damping factor the same, if the un-damped natural
frequency is changed to 20 rad/s, locate the new poles of the system? What can you say about the
response of the new system?
Solution:
Given, n1 10 rad/s and n2 20 rad/s
2
d1 n1 1 10 1 0.32 9.54 rad/s
2
d2 n2 1 20 1 0.32 19.08 rad/s
p1,2 n1 j d1 3 j 9.54
p3,4 n2 j d2 6 j19.08
0.3
tan 17.45o
2 2
1 1 0.3
IN
8.9.1. Second-order Time Response Specifications with Impulse input
(a) O
General equation
y 2 y
GA 2
y
Kxi
(t ) (8.76)
n n
m
Laplacian of the output
&
Kxi 1
Y (s)
m s2 2 n s 2
n
(8.77)
Kxi 1 1
M
2 2 2
2m n 1 (s n n 1) (s n n 1
Time-domain response
A
Kxi nt 2
y(t ) e sinh n 1 t (8.78)
2
m 1
BE
(b) Critically
General equation
2 Kxi
y n y (t ) (8.79)
m
Laplacian of the output
Kxi 1
Y (s) 2 2
(8.80)
m s n
Time-domain response
Kxi nt
y (t ) n te (8.81)
m n
(c) Under
Time-domain response
Kxi nt
y (t ) e sin d t (8.84)
m d
IN
Normalized impulse-response of a second-order system with different damping factors are shown
graphically as follows.
GA
&
M
Solved problems:
A
3. A second-order system has an un-damped natural frequency of 100 rad/s and a damping factor of
0.3. The value of the coefficient of the second time derivative (that is m) is 5. If the static
BE
sensitivity is 10, write down the response (do not solve) for a force excitation shown in the figure
in terms of the Duhamel’s integral for the following periods of time: 0<t<t1, t1<t<t2 and t>t2.
Solution:
n =100 rad/s
2
d n 1 100 1 0.32 95.39 rad/s
Here,
t
F (t ) F ;0 t t1
t1
F
F (t ) t2 t ; t1 t t2
t2 t1
t
K
y (t ) F (t )e n
sin d d
m d 0
IN
t
10F 0.3 100
y (t ) e sin 95.39 (t )d
5 95.39t1 0
t
;0 t t1 ,
0.057 F 30
e sin 95.39 (t )d
y (t )
t1 0
0.057 F 1
e
t
30
sin 95.39 (t )d
GA
t1 0
t
; t1 t t2 and
0.057 F 30
e sin 95.39 (t2 t )d
&
t2 t1 t1
t
0.057 F 1 30
y (t ) e sin 95.39 (t )d
t1 0
M
t
;t t2
0.057 F 2 30
e sin 95.39 (t2 t )d
t2 t1 t1
A
Kxi 1
Y (s) (8.85)
m s( s n n
2
1)( s n n
2
1)
Kxi nt 2 2
y (t ) 1 e cosh n 1 t sinh n 1 t (8.86)
m n2 2
1
Kxi 1
Y (s ) (8.87)
m s (s n j d )(s n j d )
Kxi nt
y (t ) 1 e cos d t sin d t (8.88)
m n2 1 2
IN
8.10. Time Response Specifications with step-input for under-damped case
For under-damped case, the step-response of a second-order is shown as follows
GA
&
A M
BE
nt
Kxi e
y (t ) 1 sin( d t ) (8.89)
m n2 1 2
2
1 1
tan (8.90)
(ii) Rise time, tr:The time required by the system response to reach from 10% to 90% of the
final value for over-damped case, from 0% to 100% of the final value for under-damped case
IN
and from 5% to 95% of the critically value for over-damped case.
.
nt r
e
y tr 1 1 sin( d tr )
2
1
e
1
n tr GA
2
sin( d tr ) 0
d tr
tr (8.92)
wd
&
(iii) Peak time, tp:The time required by the system response to reach the first maximum value.
dy t p
0
dt
M
nt p
e
d 1 sin( dtp )
2
1
0
A
dt
nt p
e
d sin( dtp )
BE
2
1
0
dt
2
1 1
wd t p tan n ; where n 1, 2,3,...
For n=1,
wd t p n
n
tp (8.93)
wd
(iv) Peak overshoot, Mp: It is the time required to reach 50% of output.
y tp 1
M p % 100
1
n tr
e
Mp % 100 1 sin( d tr ) 1
2
1
n
nt p
e e d
2 2
1 1
e e
Mp % 100 sin( d ) 100 sin( )
2 2 2
1 1 1
2 2
1 1
e e 2
Mp % 100 sin 100 1
2 2
1 1
IN
2
1
Mp % 100 e (8.94)
4
ts (8.95)
wn
&
For 5% error band,
3
ts (8.96)
wn
M
Time Specifications
Sl. No.
Type Formula
A
1 0.7
1 Delay time td
wn
BE
2 Rise time tr
wd
3 Peak time tp
wd
4
5 Settling time ts
wn
Solved Problems:
1. Consider the system shown in Figure 1. To improve the performance of the system a feedback is
added to this system, which results in Figure 2. Determine the value of K so that the damping
ratio of the new system is 0.4. Compare the overshoot, rise time, peak time and settling time and
the nominal value of the systems shown in Figures 1 and 2.
Figure 1 Figure 2
IN
Solution:
For Figure 1,
20
c s
R s
G(s)
1 G (s ) 1
s s 1
20
s s 1
s 2
20
s 20
GA
2
Here, n 20 and 2 n 1
&
1 1
n 20 rad/s and 0.112
2 n 2 20
For Figure 2,
M
20
c s G( s) s ( s 1 20 K ) 20
2
R s 1 G( s) 20 s 1 20 K s 20
A
1
s( s 1 20 K )
2
BE
n 20 rad/s
1 20 K 1 20 K
But, given that 0.4
2 n 2 20
K 0.128
Transient characteristics of Figures 1 and 2
IN
s^2+5s+36
>>step(sys)
C s 20
Program 2: Find the step response for the following system
Solution:
>> num=[20]
GA R s s2 4 s 25
num=
20
>> den=[1 425]
&
den=
1 4 25
>>sys=tf(num,den)
Transfer function:
20
M
--------------------
s^2+4s+25
>>step(sys)
A
BE
2. Stability
2.1. Concept of stability
Stability is a very important characteristic of the transient performance of a system. Any working
system is designed considering its stability. Therefore, all instruments are stable with in a boundary of
parameter variations.
A linear time invariant (LTI) system is stable if the following two conditions are satisfied.
(i) Notion-1: When the system is excited by a bounded input, output is also bounded.
Proof:
A SISO system is given by
C s b0 s m b1s m 1
... bm
G s (9.1)
R s a0 s n a1 s n 1
... an
So,
IN
1
c t G s R s (9.2)
c t
0
GA
g r t d (9.3)
1
g G s = impulse response of the system
&
Taking absolute value in both sides,
c t g r t d (9.4)
M
Since, the absolute value of integral is not greater than the integral of absolute value of the integrand
A
c t g r t d
0
BE
c t g r t d (9.5)
0
c t g r t d
0
r t M1
(9.6)
c t M2
Then,
c t M1 g d M2 (9.7)
0
(ii) Notion-2: In the absence of the input, the output tends towards zero irrespective of initial
conditions. This type of stability is called asymptotic stability.
IN
GA
&
Under-dampedclose-loop poles
Pole-zero map Normalized response
Un-dampedclose-loop poles
Pole-zero map Normalized response
IN
GA
&
Negative Under-dampedclose-loop poles
Pole-zero map Normalized response
A M
BE
IN
Closed-loop can be located by replace the denominator of the close-loop response with s=j .
Example:
1. Determine the close-loop poles on the imaginary axis of a system given below.
GA
K
G ( s)
s ( s 1)
Solution:
Characteristics equation, B ( s ) s2 s K 0
&
Replacing s jw
B( j ) ( j )2 ( j ) K 0
M
2
(K ) j 0
Comparing real and imaginary terms of L.H.S. with real and imaginary terms of R.H.S., we get
A
K and 0
Therefore, Closed-loop poles do not cross the imaginary axis.
BE
B( j ) ( j )3 6( j )2 8 j K 0
(K 6 2 ) j(8 3
) 0
Comparing real and imaginary terms of L.H.S. with real and imaginary terms of R.H.S., we get
2
8 rad/s and K 6 48
Therefore, Close-loop poles cross the imaginary axis for K>48.
Example:
3. Consider a third order polynomial B ( s ) s 3 3s 2 16 s 130 . Although the coefficients of the
above polynomial are positive, determine the roots and hence prove that the rule about
IN
coefficients being positive is only a necessary condition for the roots to be in the left s-plane.
Solution:
Characteristics equation, B ( s ) s3 3s 2 16s 130 0
By using Newton-Raphson’s method r1
GA
5 and r2,3 1 j5
Therefore, from the above example, the condition that coefficients of a polynomial should be positive
for all its roots to be in the left s-plane is only a necessary condition.
&
2.4.2. Sufficient condition of stability:
[Link].Method I (using determinants)
The coefficients of the characteristics equation are represented by determinant form
as follows.
M
an 1 an 3 an 5
an an 2 an 4
(9.8)
A
n
0 an 1 an 3
BE
Here, the determinant decreases by two along the row by one down the
column. For stability, the following conditions must satisfy.
an 1 an 3 an 5
an 1 an 3
1 an 1 0, 2 0, 3 an an 2 an 4 0 (9.9)
an an 2
0 an 1 an 3
sn an an 2 an 4
n 1
s an 1 an 3 an 5
n 2
s bn 1 bn 3 bn 5
(9.10)
n 3
s cn 1 cn 3 cn 5
Where,
(an 1 )(an 2 ) an ( an 3 )
bn 1
an 1
IN
( an 1 )( an 4 ) an ( an 5 )
bn 3 (9.11)
an 1
(bn 1 )(an 3 ) an 1 (bn 3 )
cn 1 GA
bn 1
Solution:
A
K
G (s) s ( s 1) K
In the system, T s 2
1 G (s) H (s) K s s K
BE
1
s ( s 1)
Method-I,
Characteristics equation, B s s2 s K 0
1 1
Here, 1 0
2 K
1 K
1 0
For stability,
2 0
The system is always stable for K>0.
Method-II,
Characteristics equation, B s s2 s K 0
s2 1 K
s1 1 0
s0 K
There are no sign changes in first column elements of this [Link], the system is always stable
for K>0.
K
5. Find stability of the following system given by G ( s ) and H ( s ) 1 using
s( s 2)( s 4)
Routh-Hurwitz stability criterion.
Solution:
IN
K
C (s ) G s s ( s 2)( s 4) K
In the system, 3 2
R s 1 G s H s K s 6s 8s K
1 GA
s ( s 2)( s 4)
Method-I,
General form of characteristics equation, B s a3 s 3 a2 s 2 a1s a0 0
6 K
1 8 0, 2 48 K 0,
1 8
M
6 K 0
3 1 8 0 K 48 K 0
A
0 6 K
Therefore, the system is always stable for K 48 .
BE
Method-II,
Characteristics equation is B s s3 6s2 8s K 0
andRouth’s array
s3 1 8
6 K
s2
48 K
s1 0
6
s0 K
There are no sign changes in first column elements of this array if K 48 . Therefore, the system is
always stable for 0 K 48 .
s3 1 10
2
s 5 3
1
s 9.4 0
s0 3
There are no sign changes in first column elements of this array. Therefore, the system is always
stable.
IN
7. Find stability of the following system given by B s s3 2s2 3s 10 using Routh-Hurwitz
stability criterion.
Solution:
B s s 3 2 s 2 3s 10 0 and
GA
In this problem, given characteristics equation is
Routh’s array is
s3 1 3
&
s2 2 10
s1 2 0
s 100
M
There are two sign changes in first column elements of this array. Therefore, the system is unstable.
8. Examine stability of the following system given by s5 2s4 4s3 8s2 3s 1 using Routh-Hurwitz
stability criterion.
A
Solution:
BE
s5 1 4 3
s4 2 8 1
3
s 0 2.5
2
s
s1
s0
Here, the criterion fails. To remove the above difficulty, the following two methods can be used.
Method-1
(i) Replace 0 by very small number) and complete the array with .
(ii) Examine the sign change by taking 0
Now, Routh’s array becomes
s5 1 4 3
4
s 2 8 1
s3 2.5 0
s2 5 8
1 0
5 8
s1 2.5
5 8
s0 1
IN
4
s 2 8 1
s3 2.5 0
s2 5 8
1 0
s1 2.5
5 8
5 8
GA
s0 1
&
There are two sign changes in first column elements of this array. Therefore, the system is unstable.
Method-2
Replace s by 1 . The system characteristic equation s5 2s4 4s3 8s2 3s 1 0 becomes
M
Z
1 2 4 8 3
1 0
Z5 Z4 Z3 Z2 Z
A
Z5 3Z 4 8Z 3 4Z 2 2Z 1 0
Now, Routh’s array becomes
BE
s5 1 8 2
s4 3 4 1
s 3 6.67 1.67 0
s 2 3.25 1 0
s1 0.385 0 0
1 0 0
s0
There are two sign changes in first column elements of this array. Therefore, the system is unstable.
9. Examine stability of the following system given by s5 2s4 2s3 4s2 4s 8 using Routh-Hurwitz
stability criterion.
Solution:
In this problem, Routh’s array is
s5 1 2 4
s4 2 4 8
s3 0 0 0
s2
s1
s0
Here, the criterion fails. To remove the above difficulty, the following two methods can be used.
The auxillary equation is
As 2s4 4s2 8
dA s
8s3 8s
ds
Now, the array is rewritten as follows.
IN
s5 1 2 4
s4 2 4 8
3
s 8 8 0
s2
s1
2 8 0
24 0
8
GA
s0
There are two sign changes in first column elements of this array. Therefore, the system is
unstable.
&
10. Examine stability of the following system given by s4 5s3 2s2 3s 1 0 using Routh-Hurwitz
stability criterion. Find the number of roots in the right half of the s-plane.
Solution:
M
2
s 1 .4 2
s1 4 .1 4 0
BE
s 0 2
There are two sign changes in first column elements of this array. Therefore, the system is unstable.
There are two poles in the right half of the s-plane.
IN
GA
&
A M
BE
IN
MODULE#3 GA
Equation Chapter (Next) Section 1
&
A M
BE
CHAPTER#10
10. Root locus
10.1. Definition:
The locus of all the closed-loop poles for various values of the open-loop gain K is called root locus.
The root-locus method is developed by W.R. Evans in 1954. It helps to visualize the various
possibilities of transient response of stable systems.
Closed-loop response function
C ( s) G ( s)
(10.1)
R(s) 1 G (s )H (s)
Characteristic equation
K ( s z1 )( s z2 )...( s zm )
1 G (s) H (s) 1 0 (10.2)
( s p1 )( s p2 )...( s pn )
IN
Vector from open-loop pole to the root-locus
GA
&
K ( s z1 )( s z2 )...( s zm )
1 G (s) H (s) 1 0 (10.3)
( s p1 )( s p2 )...( s pn )
IN
For K=0,
(s p1 )( s p2 )...( s pn ) K ( s z1 )( s z2 )...( s zm )
0
(s p1 )( s
GA
( s p1 )( s p2 )...( s pn )
p2 )...( s pn ) 0 (10.4)
Open-loop poles are also closed-loop poles for K=0. A root-locus starts from every open-loop pole.
&
10.2.3. Ending points
Characteristics equation of a closed-loop system
K ( s z1 )( s z2 )...( s zm )
1 G (s) H (s) 1 0 (10.5)
M
( s p1 )( s p2 )...( s pn )
For K=
A
K ( s z1 )( s z 2 )...( s z m )
1
( s p1 )( s p2 )...( s pn )
BE
(s z1 )( s z2 )...( s z m ) 0 (10.6)
Angle criterion:
n m
Where, i angle in case of ith pole and j angle in case of jth zero
Magnitude criterion:
G (s ) H ( s) 1 (10.9)
Example:
Determine K of the characteristic equation for the root s=-0.85.
Solution:
IN
S3 6s 2 8s K 0 (10.11)
1 6
-0.85
GA
8
-4.378
K
-3.079
2. Start from open-loop poles on the real axis, extend on the real axis for increasing
values of the gain and end at an infinite value on the real axis.
3. Start from a pair of open-loop poles on the real axis, extend on the real axis for
A
increasing values of gain, meet at a point and then leave the real axis and end at a
complex open-loop zero or infinity.
BE
4. Start from a pair of open-loop poles on the real axis, extend on the real axis for
increasing values of gain, meet at a point and then leave the real axis. They may once
again enter the real axis and end at open-loop zeros or at a large value on the real axis.
5. Start from a pair of complex open-loop poles, enter the real axis and end at an open-
loop zero or an infinite value on the real axis. They could leave the real axis again and
end at a complex open-loop zero or infinity.
IN
GA
10.2.8. Determine breakaway and break-in points of the root-locus
&
A( s)
1 G ( s) H ( s) 1 K 0 (10.12)
B( s )
B(s)
K (10.14)
As )
A
f (s ) (s s1 ) r ( s s 2 )...( s sn r 1) 0 (10.15)
df ( s )
BE
df ( s )
0 (10.17)
ds s s1
B ' ( s)
K (10.19)
A' ( s )
Therefore,
dK B ' ( s ) A( s ) B ( s ) A' ( s )
0 (10.21)
ds A2 ( s )
IN
GA
Angle of asymptotes
1800 k 360
c where, k=0, 1, 2, 3..
(n m )
&
Location of asymptotes
( s p1 )( s p2 )...( s pn )
K (10.22)
( s z1 )( s z2 )...( s zm )
M
s n ( p1 p2 ... pn ) s n 1
K (10.23)
s m ( z1 z 2 ...zm ) s m 1 ...
A
K sn m
[( p1 p2 ... pn ) ( z1 z2 ... zm )]s n m 1
... (10.24)
BE
s pi s c (10.25)
( s zi ) s c (10.26)
n
(s c)
K m
sn m
( n m) cs
n m 1
... (10.27)
(s c)
( p1 p2 ... pn ) ( z1 z2 ...zm )
c (10.28)
(n m)
Angle of departure
d 180 ( 1 2) 1 (10.29)
o
d=180 - angles of vectors to the complex open-loop pole in question from other open - loop poles
+ angles of vectors to the complex open-loop pole in question from all open-loop zeros
IN
Angle of arrival
a 180 ( 1 3 ) ( 1 2 3 ) (10.30)
GA
&
A M
o
a=180- angles of vectors to the complex open-loop zero in question from other open- loop zeros
+ angles of vectors to the complex open-loop zero in question from all open-loop poles
BE
imaginary[1 G ( j ) H ( j )] 0 (10.32)
Example
Problem-1: Draw the root-locus of the feedback system whose open-loop transfer function is given
K
by G ( s) H ( s)
s( s 1)
Solution:
Step 1: Determine the number of open-loop poles and zeros
Number of open-loop poles n=2
Number of open-loop zeros m=0
Open-loop poles: s=0 and s=-1
IN
Test points in between the open-lop poles
GA
&
A M
BE
IN
c 900 k 0
0
c 270 k 1
Centroid of asymptotes
c
( p1 p2 ... pn ) ( z1
( n m)
z2 ... zm ) 0 1
2
GA
0.5
&
A M
Steps 6 & 7: Since there are no complex open-loop poles or zeros, angle of departure and arrival need
BE
not be computed
Step 8: Determine points on the root-locus crossing imaginary axis
K
1 GH 1 s2 s K 0
s ( s 1)
B( j ) ( j )2 ( j ) K (K 2
) j
2
K 0 j 0
The root-locus does not cross the imaginary axis for any value of K>0
Here,
1 1 4K
s
2
IN
GA
&
M
Problem-2: Draw the root-locus of the feedback system whose open-loop transfer function is given
K
by G ( s) H ( s)
s( s 2)( s 4)
A
Solution:
Step 1: Determine the number of open-loop poles and zeros
BE
IN
GA
&
A M
BE
1 6 8 K
IN
Step 5: Draw asymptotes of the root-locus
Angle of asymptotes:
c
1800 k 360
( n m)
180 360 k
3
GA
c 600 k 0
0
c 180 k 1
&
0
c 300 k 2
Centroid of asymptotes
( p1 p2 ... pn ) ( z1 z2 ...z m ) 0 2 4
c 2
M
( n m) 3
A
BE
Steps 6 & 7: Since there are no complex open-loop poles or zeros, angle of departure and arrival need
not be computed
1 6 8 48
1 6+j2.828 J16.97 0
1 6+j2.828 J16.97
-j2.828 -j16.97
IN
1 6 0
No.
Closed-loop pole
on the real axis
K
GA Second and third
closed- loop poles
Remarks
1 6 8 K
( s 2 1.5 s 1.25) 0
s2,3 0.75 j 0.829
IN
Problem-3: Draw the root-locus of the feedback system whose open-loop transfer function is given
K GA
by G ( s ) H ( s ) 2
s ( s 1)
Solution:
Step 1: Determine the number of open-loop poles and zeros
Number of open-loop poles n=3
Number of open-loop zeros m=0
Open-loop poles: s=0, s=0 and s=-1
&
Step 2: Mark open-loop poles and zeros on the s-plane
A M
BE
IN
Characteristic equation, K
dK
0
ds
2 s ( s 1) s 0
s 2s 3
Breakaway point as
b
0
b= -2/3and 0
GA
= -2/3is not on the root-locus and therefore not a breakaway or break-in point.
Therefore b= 0 and the two loci start from the origin and breakaway at the origin itself.
c
0
c 180 k 1
c 3000 k 2
A
Centroid of asymptotes
( p1 p2 ... pn ) ( z1 z2 ... z m ) 0 1 1
c
BE
( n m) 3 3
Steps 6 & 7: Since there are no complex open-loop poles or zeros, angle of departure and arrival need
not be computed.
2 -2.0 4 0.50±j1.32
IN
3 -2.5 9.375 0.75±j1.78
4 -3.0 18 1.00±j2.23
( s 2 1.5 s 1.25) 0
A
IN
Problem-4: Draw the root-locus of the feedback system whose open-loop transfer function is given
by G ( s ) H ( s)
Solution:
s 4
5s 3
K
8s 2 6s
GA
Step 1: Determine the number of open-loop poles and zeros
s4 5s 3 8 s 2 6s s(s2 2 s 2)( s 3) (s 1 j )( s 1 j )( s 3) s
Number of open-loop poles n=4
&
Number of open-loop zeros m=0
Open-loop poles: s=0 and s=-3, s=-1+j and s=-1-j
IN
Step 4: Determine breakaway and break-in point
Characteristic equation, K ( s 4 5s 3 8 s 2 6 s )
dK
ds
0
4 s 3 15s 2 16s 6 0
GA
s 3 3.75s 2 4s 1.5 0
f ' ( s ) 3s 2 7.5 s 4
This equation is solved using Newton-Raphson’s method
&
f ( sn )
sn 1 sn
f ' ( sn )
No.
M
sn f ( sn ) f ' ( sn ) sn 1
-5
6 -2.2886
-5.03x10
1 5 8 6 K
IN
Other closed-loop poles for K=4.3
GA
1 2.7114 1.7947 1.893
&
-2.2886 -0.9676 -1.893
1 0.4228 0.8270 0
M
s =-0.2114±j0.8814
3,4
Angle of asymptotes:
1800 k 360 180 360 k
c
( n m) 4
BE
c 450 k 0
0
c 135 k 1
c 2250 k 2
0
c 315 k 3
Centroid of asymptotes
( p1 p2 ... pn ) ( z1 z2 ...z m ) 0 3 1 j 1 j
c 1.25
( n m) 4
IN
GA
&
Step 7: As there are no complex open-loop zeros, angle of arrival need not be computed.
No. S1 S2 S3,4 K
IN
1 -0.25 -2.9217 -0.9142±0.7969j 1.0742
3
4
-0.75
-1.0
-2.8593
-2.8393
GA -0.6953±0.5938j
-0.5804±0.6063j
1.7930
2.0000
IN
GA
&
M
2
1 1
Where, tan
3. Percentage overshoot
A
/ tan
Mp e (10.34)
BE
4. Settling time
4
ts (10.35)
n
Example
Problem-1: Draw the root-locus of the feedback system whose open-loop transfer function is given
K s 2 10s 100
by G ( s ) H ( s ) ,H s 1
s 4 20 s 3 100 s 2 500 s 1500
(a) Determine the value of gain at which the system will be stable and as well have a maximum
overshoot of 5%.
(b) What is the gain margin at this point?
(c) What is the steady-state error for a unit step excitation at the above point?
Solution:
tan 1.0487
(a) ln M p
460
1
0.690 (10.36)
1 tan 2
IN
GA
&
A M
192.2
(b) GM 20 log 2.65dB
261
BE
Root locus
The locus of all the closed-loop poles for various values of the open-loop gain K is called root locus.
The root-locus method is developed by W.R. Evans in 1954. It helps to visualize the various
possibilities of transient response of stable systems.
Closed-loop response function
C ( s) G ( s)
(10.37)
R(s) 1 G (s )H (s)
Characteristic equation
K ( s z1 )( s z2 )...( s zm )
1 G (s) H (s) 1 0 (10.38)
( s p1 )( s p2 )...( s pn )
IN
GA
&
Vector from open-loop zero to the root-locus
A M
BE
Starting points
Characteristics equation of a closed-loop system
K ( s z1 )( s z2 )...( s zm )
1 G (s) H (s) 1 0 (10.39)
( s p1 )( s p2 )...( s pn )
For K=0,
IN
(s p1 )( s p2 )...( s pn ) K ( s z1 )( s z2 )...( s zm )
0
( s p1 )( s p2 )...( s pn )
(s p1 )( s GA p2 )...( s pn ) 0
Open-loop poles are also closed-loop poles for K=0. A root-locus starts from every open-loop pole.
(10.40)
Ending points
Characteristics equation of a closed-loop system
&
K ( s z1 )( s z2 )...( s zm )
1 G (s) H (s) 1 0 (10.41)
( s p1 )( s p2 )...( s pn )
M
For K=
K ( s z1 )( s z 2 )...( s z m )
1
( s p1 )( s p2 )...( s pn )
A
(s z1 )( s z2 )...( s z m ) 0 (10.42)
BE
Angle criterion:
n m
Where, i angle in case of ith pole and j angle in case of jth zero
Magnitude criterion:
G ( s ) H ( s) 1 (10.45)
Using the magnitude of vectors drawn from open-loop poles and zeros to the root-locus point, we get
n
(s pi )
i 1 | ( s p1 ) || ( s p2 ) | ... | ( s pn ) |
m
K (10.46)
| ( s z1 ) || s z 2 ) | ... | ( s zm ) |
(s z j )
j 1
Example:
Determine K of the characteristic equation for the root s=-0.85.
Solution:
S3 6s 2 8s K 0 (10.47)
IN
1 6 8 K
1 5.15
GA
3.622 K-3.079=0
again enter the real axis and end at open-loop zeros or at a large value on the real axis.
10. Start from a pair of complex open-loop poles, enter the real axis and end at an open-
loop zero or an infinite value on the real axis. They could leave the real axis again and
BE
IN
B(s)
K (10.50)
As )
df ( s )
ds
f (s )
r ( s s1 ) r 1 ( s
(s
s2 )...( s sn
GA
s1 ) r ( s s 2 )...( s
r 1)
sn r 1) 0
( s s1 ) r .( s s3 )..( s sn r 1) ...
(10.51)
(10.52)
df ( s )
0 (10.53)
&
ds s s1
B ' ( s)
M
K (10.55)
A' ( s )
Therefore,
A
IN
Angle of asymptotes
1800 k 360
c where, k=0, 1, 2, 3..
(n m) GA
Location of asymptotes
( s p1 )( s p2 )...( s pn )
K (10.58)
( s z1 )( s z2 )...( s zm )
s n ( p1 p2 ... pn ) s n 1
&
K (10.59)
s m ( z1 z 2 ...zm ) s m 1 ...
K sn m
[( p1 p2 ... pn ) ( z1 z2 ... zm )]s n m 1
... (10.60)
M
s pi s c (10.61)
( s zi ) s c (10.62)
A
n
(s c)
K m
sn m
( n m) cs
n m 1
... (10.63)
BE
(s c)
( p1 p2 ... pn ) ( z1 z2 ...zm )
c (10.64)
(n m)
Angle of departure
d 180 ( 1 2) 1 (10.65)
o
d=180 - angles of vectors to the complex open-loop pole in question from other open - loop poles
+ angles of vectors to the complex open-loop pole in question from all open-loop zeros
IN
Angle of arrival
a 180 ( 1 3 ) ( 1 2 3 ) (10.66)
GA
&
A M
o
a=180- angles of vectors to the complex open-loop zero in question from other open- loop zeros
+ angles of vectors to the complex open-loop zero in question from all open-loop poles
BE
imaginary[1 G ( j ) H ( j )] 0 (10.68)
Example
Problem-1: Draw the root-locus of the feedback system whose open-loop transfer function is given
K
by G ( s) H ( s)
s( s 1)
Solution:
Step 1: Determine the number of open-loop poles and zeros
Number of open-loop poles n=2
Number of open-loop zeros m=0
Open-loop poles: s=0 and s=-1
IN
Test points in between the open-lop poles
GA
&
A M
BE
IN
c 900 k 0
0
c 270 k 1
Centroid of asymptotes
( p1 p2 ... pn ) ( z1 z2 ... zm ) 0 1
c
( n m) 2
GA
0.5
&
A M
Steps 6 & 7: Since there are no complex open-loop poles or zeros, angle of departure and arrival need
not be computed
BE
Here,
1 1 4K
s
2
IN
GA
&
Problem-2: Draw the root-locus of the feedback system whose open-loop transfer function is given
M
K
by G ( s) H ( s)
s( s 2)( s 4)
Solution:
A
IN
GA
&
A M
BE
1 6 8 K
IN
Step 5: Draw asymptotes of the root-locus
Angle of asymptotes:
c
1800 k 360
(n m)
180 360 k
3
GA
c 600 k 0
1800 k 1
&
c
c 3000 k 2
Centroid of asymptotes
( p1 p2 ... pn ) ( z1 z 2 ...z m ) 0 2 4
2
M
c
( n m) 3
A
BE
Steps 6 & 7: Since there are no complex open-loop poles or zeros, angle of departure and arrival need
not be computed
1 6 8 48
1 6+j2.828 J16.97 0
1 6+j2.828 J16.97
-j2.828 -j16.97
IN
1 6 0
No.
Closed-loop pole
on the real axis
K
GA Second and third
closed- loop poles
Remarks
1 6 8 K
( s 2 1.5 s 1.25) 0
s2,3 0.75 j 0.829
IN
Problem-3: Draw the root-locus of the feedback system whose open-loop transfer function is given
K GA
by G ( s ) H ( s ) 2
s ( s 1)
Solution:
Step 1: Determine the number of open-loop poles and zeros
Number of open-loop poles n=3
Number of open-loop zeros m=0
Open-loop poles: s=0, s=0 and s=-1
&
Step 2: Mark open-loop poles and zeros on the s-plane
A M
BE
IN
Characteristic equation, K
dK
0
ds
2 s ( s 1) s 0
b
s 2s 3 0
b= -2/3and 0
GA
= -2/3is not on the root-locus and therefore not a breakaway or break-in point.
Therefore b = 0 and the two loci start from the origin and breakaway at the origin itself.
c
0
c 180 k 1
c 3000 k 2
A
Centroid of asymptotes
( p1 p2 ... pn ) ( z1 z2 ... z m ) 0 1 1
c
BE
( n m) 3 3
Steps 6 & 7: Since there are no complex open-loop poles or zeros, angle of departure and arrival need
not be computed.
2 -2.0 4 0.50±j1.32
IN
3 -2.5 9.375 0.75±j1.78
4 -3.0 18 1.00±j2.23
( s 2 1.5 s 1.25) 0
A
IN
Problem-4: Draw the root-locus of the feedback system whose open-loop transfer function is given
by G ( s ) H ( s)
Solution:
s 4
5s 3
K
8s 2 6s
GA
Step 1: Determine the number of open-loop poles and zeros
s4 5s 3 8 s 2 6s s(s2 2 s 2)( s 3) (s 1 j )( s 1 j )( s 3) s
Number of open-loop poles n=4
&
Number of open-loop zeros m=0
Open-loop poles: s=0 and s=-3, s=-1+j and s=-1-j
IN
Step 4: Determine breakaway and break-in point
Characteristic equation, K ( s 4 5s 3 8 s 2 6 s )
dK
ds
0
4 s 3 15s 2 16s 6 0
GA
s 3 3.75s 2 4s 1.5 0
f ' ( s ) 3s 2 7.5s 4
This equation is solved using Newton-Raphson’s method
&
f ( sn )
sn 1 sn
f ' ( sn )
No.
M
sn f ( sn ) f ' ( sn ) sn 1
b= -2.3
Gain at the breakaway point, K | 2.3 ( 3) || 2.3 0 || 2.3 ( 1 j ) || 2.3 ( 1 j ) | 4.33
1 5 8 6 K
IN
1 2.7114 1.7947 1.893
-2.2886 -0.9676 -1.893
s =-0.2114±j0.8814
1 0.4228
GA 0.8270 0
3,4
c
0
c 135 k 1
c 2250 k 2
A
0
c 315 k 3
BE
Centroid of asymptotes
( p1 p2 ... pn ) ( z1 z2 ...z m ) 0 3 1 j 1 j
c 1.25
( n m) 4
IN
GA
&
Step 7: As there are no complex open-loop zeros, angle of arrival need not be computed.
No. S1 S2 S3,4 K
IN
1 -0.25 -2.9217 -0.9142±0.7969j 1.0742
3
4
-0.75
-1.0
-2.8593
-2.8393
GA -0.6953±0.5938j
-0.5804±0.6063j
1.7930
2.0000
IN
GA
&
M
2
1 1
Where, tan
8. Percentage overshoot
A
/ tan
Mp e (10.70)
BE
9. Settling time
4
ts (10.71)
n
Example
Problem-1: Draw the root-locus of the feedback system whose open-loop transfer function is given
K s 2 10s 100
by G ( s ) H ( s ) ,H s 1
s 4 20 s 3 100 s 2 500 s 1500
(a) Determine the value of gain at which the system will be stable and as well have a maximum
overshoot of 5%.
(b) What is the gain margin at this point?
(c) What is the steady-state error for a unit step excitation at the above point?
Solution:
tan 1.0487
(b) ln M p
460
1
0.690 (10.72)
1 tan 2
IN
GA
&
A M
192.2
(b) GM 20 log 2.65dB
261
BE
num=
0 1
>>q1=[1 1];
>> q2=[1 2];
>> q3=[1 3];
>> q4=[1 4];
>>den=conv(q1,q2);
>> den=conv(den,q3);
>> den=conv(den,q4);
den=
1 12 49 78 40
>>sys=tf(num,den)
Transfer function:
1
--------------------
s^4+12s^3+49s^2+78s+40
IN
>>rlocus(sys)
GA
&
A M
BE
Here,
C s N s
G s (11.1)
R s s a s b s c ...
IN
N s R s
C s (11.2)
s a s b s c ...
R s
GA A
(11.3)
s2 2
C s ...
s a s b s c s j s j
at bt ct j t
c t A1e A2 e A3e ... B1e B2 e j t
(11.5)
BE
The term with Ai terms are decaying components. So, they tend to zero as time tends to infinity.
Then, eq (5) becomes
j t
Css t B1e B2e j t
(11.6)
Where,
A G s A j G j
B1 G j e
s j s j
2j
(11.7)
A G s A j G j
B2 G j e
s j s j
2j
Since, G j G j and G j G j
A j t A j t
c t G j e G j e (11.8)
2j 2j
j t ej e j
c t AG j e (11.9)
2j
c t AG j sin t (11.10)
c t B sin t (11.11)
Where, B AG j
Therefore, the steady-state response of the system for a sinusoidal input of magnitude A and
frequency is a sinusoidal output with a magnitude B , frequency and phase shift .
IN
The following plots are used in frequency response.
Polar plot
Bode plot
Magnitude versus phase angle plot
11.2. Definition of frequency domain specifications
GA
&
A M
BE
(v) Phase cross-over frequency: The frequency at which phase plot crosses -1800
(vi) Gain margin (GM): It is the increase in open-loop gain in dB required to drive the closed-loop
system to the verge of instability
(vii) Gain cross-over frequency: The frequency at which gain or magnitude plot crosses 0dB line
(viii) Phase margin (PM): It is the increase in open-loop phase shift in degree required to
drive the closed-loop system to the verge of instability
2
C s n
2 2
(11.12)
R s s 2 ns n
IN
Putting s j
2
C j n
2 2
R j n j2 n
C j
R j
GA
1
2
2
1
j2
(11.13)
n n
Let, u , then
&
n
C j 1
2
(11.14)
R j 1 u j2 u
M
Now,
A
M j M j M j (11.15)
Where,
BE
1
M j
2 2
1 u2 2 u
(11.16)
1 2 u
tan
1 u2
Now,
1
Mr (11.17)
2
2 1
2
r n 1 2 (11.18)
2 4 2
b n 1 2 4 4 2 (11.19)
PM 1800 (11.20)
1 2
Where, tan
2 2
4 1 2
11.4. Advantages
Good accuracy
Possible to test in lab
Can be used to obtain transfer function that is not possible with analytical techniques
Easy to design open-loop transfer function from closed-loop performance in frequency
domain
IN
It is very easy to visualize the effect of disturbance and parameter variations.
11.5. Disadvantages
Applied only to linear systems
Frequency response for existing system is possible to obtain if the time constant is up to few
minutes
Time consuming procedure
Old and back dated method
GA
Equation Chapter 12 Section 1
&
A M
BE
IN
GA
M 20log | G ( j ) H ( j ) | dB
&
Phase plot on a semi-log paper
A M
BE
IN
GA
&
1 K K Constant
2 s j Derivative factor
BE
6 s2 2 ns
2
n
2
n
2
j2 n Second order derivative factor
1 1
7 2 2 2 2 Second order integral factor
s 2 ns n n j2 n
j 900 (12.2)
2
M 20log 2 20log 1 20log dB/decade (12.3)
1
Table 12.2 Magnitude variation of a derivative factor for various multiples of the initial
frequency
IN
2
1 2 3 4 5 6 7 8 9 10
1
M dB 0 6 10 12
GA 14 16 17 18 19 20
&
A M
BE
IN
GA
&
Table 15.3Derivative factor
Frequency, rad/s
M
0.1 1 10 30 100
Magnitude, dB -20 0 20 30 40
Phase, degrees 90 90 90 90 90
A
BE
j 270 0 (12.7)
2
M 20log 2 20log 1 20log dB/decade (12.8)
1
IN
GA
Table 12.4Magnitude variation of an integral factor for various multiples of the initial frequency
2
1 2 3 4 5 6 7 8 9 10
1
&
M , dB 0 -6 -10 -12 -14 -16 -17 -18 -19 -20
Frequency, rad/s
0.1 1 10 20 100
IN
For << , M
c
For >> ,
c
M 20log dB (12.12)
Here,
c
=1/ = corner frequency
GA c
For >
c
2
M 20log 2 20log 1 20log (12.13)
&
1
Table 12.6Magnitude variation of a first-order derivative factor for various multiples of the corner
frequency
A
1 2 3 4 5 6 7 8 9 10
BE
M, 0 6 10 12 14 16 17 18 19 20
dB
wc
0 ;w
10
wc
450 1 log ; w 10 wc (12.17)
c 10
90 ; w 10 wc
Table 12.7Phase angles of a first-order derivative factor around the corner frequency
1 2 3 4 5 6 7 8 9 10
c
, deg 45 59 66 72 76 80 83 86 88 90
Frequency, rad/s
IN
0.1 1 5 10 20 100
Magnitude, dB 0 3 14 20 26 40
Phase, degrees 0
GA
45 76 90 90 90
IN
GA
&
12.10. First-order integral factor: magnitude
1 1
M 20log 20log 2
dB (12.18)
1 j 1
M
M 0, w wc
(12.19)
M 20log dB w wc
A
2
M 20log 20log 20log dB/decade (12.20)
BE
2 1
1
Table 12.9Magnitude variation of a first-order integral factor for various multiples of the corner
frequency
1 2 3 4 5 6 7 8 9 10
c
M,
0 -6 -10 -12 -14 -16 -17 -18 -19 -20
dB
Table 12.10Phase angles of a first-order integral factor around the corner frequency
1 2 3 4 5 6 7 8 9 10
c
deg 315 301 294 288 284 280 277 274 272 270
IN
c
=270, >10
GA
c
Phase, degrees 360 360 322 315 277 270 270 270
A
IN
First-order integral factor: phase
GA
&
A M
BE
M 40logwn , w wn
2
M 20log (2zwn ), w wn (12.23)
M 40 logw, w wn
For w wn
2
M 40log 2 40 log 1 40log dB/decade (12.24)
1
Magnitude variation of a second-order derivative factor for various multiples of the resonant
frequency
1 2 3 4 5 6 7 8 9 10
n
M dB 0 12 20 24 28 32 34 36 38 40
IN
2
2 2 n
| n j2 n | arctan 2
(12.27)
1 2
n
00 ,
GA w
wn
10
900 , w wn (12.28)
0
180 , w 10 wn
&
n
M
Magnitude, dB 0 0 -4 -4 18 40 80
A
IN
GA
&
A M
1 1
M 20 log 2 2
dB 20 log dB (12.29)
n j2 n 2
2 2
2
n 1 2
2
n n
M - 40log n, << n
2
M=-20log (2 n ), = n
M= - 40 log , >> n
2
M 40log 2 40log 1 40log dB / decade (12.30)
1
Magnitude variation of a second-order integral factor for various multiples of the resonant frequency
1 2 3 4 5 6 7 8 9 10
n
2
1 n
2 2
360 arctan 2
(12.32)
n j2 n
1 2
n
IN
=0, < n
=2700, = n
=1800, > GA
n
Magnitude plot
A
BE
IN
Phase plot
GA
&
A M
BE
Example 14.1
Draw the Bode magnitude and phase plot of the following open-loop transfer function and determine
gain margin, phase margin and absolute stability?
1
G (s) H (s)
s( s 1)
Solution
Applying s j ,
1
G( j )H ( j )
j (j 1)
The above frequency response function has two factors: (1) Integral factor and (2) First order integral
factor with a corner frequency of 1 rad/s
1 40 20 0 -20 -40
20log dB
j
1 0 0 -3 -20 -40
IN
20log dB
j 1
p
= 100 rad/s
GA
&
Frequency, rad/s
1
degrees 270 270 270 270 270
j
A
1
degrees 360 360 315 270 270
BE
j 1
GM=80 dB
IN
GA
&
A M
BE
Example 14.2
Draw the Bode magnitude and phase plot of the following open-loop transfer function and determine
gain margin, phase margin and absolute stability?
1
G (s) H (s)
s ( s 2) s 4)
Solution
1
G( j )H ( j )
j j
8j 1 1
2 4
The corner frequencies corresponding to first order integral factors are 2 rad/s and 4 rad/s. Minimum
frequency is chosen as 0.01 rad/s and maximum frequency 100 rad/s.
Table 14.1 Computation of Bode magnitude using asymptotic properties of the integral first-order
IN
1
term
2
x1 x2 x1 x10 x2 x1 x1 x2 x1 x10
Frequency, rad/s 2 4 2 20 20 10 20 40 10 100
Magnitude, dB 0 -6 0 -20 GA -20 -14
Table 14.2 Computation of Bode magnitude using asymptotic properties of the integral first-order
-20 -26 -14 -34
1
term
4
x1 x10 x2 x1 x2 x1 x1 x10
&
Frequency, rad/s 4 40 40 20 20 10 10 100
Magnitude, dB 0 -20 -20 -14 -14 -8 -8 -28
Frequency, rad/s
1 -18 -18 -18 -18 -18 -18 -18 -18 -18 -18 -18
20 log
8
BE
dB
Bode magnitude
IN
GA
Bode phase
&
Frequency, rad/s
1 0 0 0 0 0 0 0 0 0 0 0
M
1 270 270 270 270 270 270 270 270 270 270 270
A
j
BE
1 360 360 360 346 328 315 301 284 270 270 270
j
1
2
1 360 360 360 360 342 326 315 297 285 270 270
j
1
4
Phase 270 270 270 256 220 191 166 131 105 90 90
degrees
Phase plot
IN
GA
&
A M
BE
Bode plot
IN
GA
&
Example 12.1
Draw the Bode magnitude and phase plot of the following open-loop transfer function and determine
gain margin, phase margin and absolute stability?
M
1
G ( s) H ( s) 2
s ( s 1)
A
Solution
1
G ( j )H ( j )
BE
( j )( j )( j 1)
There are two integral factors and an integral first-order term with a corner frequency of 1 rad/s
Bode magnitude
Frequency, rad/s
1
20log dB 40 20 0 -20 -40
j
1
20log dB 40 20 0 -20 -40
j
1
20log dB 0 0 -3 -20 -40
j 1
Example 12.2
Draw the Bode magnitude and phase plot of the following open-loop transfer function and determine
gain margin, phase margin and absolute stability?
1
G (s) H (s) 4 3
s 5s 8s 2 6s
Solution
IN
1
G (s )H (s) 2
s(s 2s 2)(s 3)
1
G( j )H ( j )
j ( j )2 2( j ) 2 (( j ) 3) GA
1
G( j )H ( j ) 3
2
j (2 1) ) j2 ) 2 ( j
3
Comparing the second order term with a standard second order term,
&
2 2
n j2 n
1
M
n 2 and .
2
For > 0.5, the response at resonance is less than theresponse at frequencies less than the resonant
frequencies
BE
Table Computation of Bode magnitude using asymptotic properties of the integral second-order term
x1 x10 x1 x2 x3 x1 x1 x10 x3 x1
Frequency, rad/s 1.4 14 14 30 30 10 10 100 30 3
- -
Magnitude, dB -6 -46 -46 -58 -38 -78 -58 -18
58 38
Table Computation of Bode magnitude using asymptotic properties of the integral first-order term
x1 x3 x2 x1 x3 x1 x1 x10
Frequency, rad/s 3 30 30 14 30 10 10 100
Magnitude, dB 0 -20 -20 -14 -20 -10 -10 -30
Bode magnitude
Frequency, rad/s
n c
1 -10 -10 -10 -10 -10 -10 -10 -10 -10 -10 -10
20 log
3
IN
20log 2
(2 ) j (2 )
Bode phase
Frequency, rad/s
n c
1 0 0 0 0 0 0 0 0 0 0 0
3
1 270 270 270 270 270 270 270 270 270 270 270
degrees
j
1 360 360 360 343 297 270 221 192 180 180 180
IN
2
(2 ) j (2 )
degrees
1 360 360 360 360 336 330 315 291 285 270 270
j
3
1
, degrees GA
Bode phase, degrees 270 270 250 253 183 150 86 33 15 0 0
&
A M
BE
Nichols plot
IN
13. Polar Plots
GA
It is a graphical method of determining stability of feedback control systems by using the polar plot of
their open-loop transfer functions.
Example 13.1
Draw a polar plot of the open-loop transfer functionfor
&
K
G (s) H ( s) (14.33)
s ( s 1)
Frequency response
M
K
G( j )H ( j ) (14.34)
j (j 1)
Magnitude
A
K
G( j )H ( ) (14.35)
2
1
BE
Angle
1
G( j ) H ( j ) tan (14.36)
2
Frequency, Phase,
No. Magnitude
rad/s degrees
1 0 270
5 1 0.7071 225
6 4 0.0606 194
7 10 0.01 186
8 50 0.0004 181
IN
9 100 0.0001 181
10 200
=
A
BE
=0.2
Example 14.2
Draw a polar plot of the open-loop transfer functionfor K=1, 10, 25, 55
K
GH
s ( s 2)( s 4)
Solution
Frequency response
K
G( j )H ( j )
j (j 2)( j 4)
Magnitude
K
G( j )H ( j )
2 2
4 16
Angle
1 1
G( j ) H ( j ) tan tan
2 2 4
0
The lies in II and III quadrants as 90 270 0
G( j )H ( j )
Magnitude and phase of the open-loop frequencytransfer function (K=1)
IN
4 0.4 0.3049 253
7
1
4
GA0.1085
0.0099
229
162
8 10 0.0009 123
9 50 0 97
&
K
Polar plot of the transfer function GH for K=1, 10, 25, 55
s( s 2)( s 4)
A M
BE
Example 14.3
K
Draw a polar plot of the open-loop transfer function G ( s ) H ( s ) 2
s ( s 1)
Solution
Frequency response
K
G( j )H ( j )
( j )2 ( j 1)
Magnitude
K
G( j )H ( j )
2 2
1
Angle
G( j )H ( j ) 180 0 tan 1
The lies in II quadrant only as 90 0 G ( j ) H ( j ) 180 0
Magnitude and phase of the open-loop frequency transfer function (K=1)
Frequency, Phase,
No. Magnitude
IN
rad/s degrees
4
0.5
0.8
GA 3.5777
1.2201
153
141
5 1 0.7071 135
&
6 2 0.1118 117
7 3 0.0351 108
M
8 4 0.0152 104
9 5 0.0078 101
A
K
Polar plot of the transfer function GH for K=1, 10, 25, 55
BE
s( s 2)( s 4)
IN
GA
&
C s 40
Program 1: Sketch the bode plot for the open loop transfer function .
R s s s 1 s 4
Determine the gain margin, phase margin, gain cross over frequency and phase cross over frequency.
A
Solution:
>> num=[0 40]
num=
BE
0 40
>> q1=[1 0];
>> q2=[1 1];
>> q3=[1 4];
>> den=conv(q1,q2);
>> den=conv(den,q3);
den=
1 9 24 16 0
>>sys=tf(num,den)
Transfer function:
1
--------------------
s^4+9s^3+24s^2+16s
>>bode(sys)
>>margin(sys)
[Link]
Feedback transfer function
C ( s) G ( s)
(14.1)
R(s) 1 G (s )H (s)
(s p1 )( s p 2 )...( s pn ) K ( s z1 )( s z 2 )...( s zm )
1 G (s) H (s) (14.3)
IN
( s p1 )( s p 2 )...( s pn )
Angle
M
The s-plane to 1+GH plane mappingphase angle of the 1+G(s)H(s) vector, corresponding to a point on
the s-plane is the difference between the sum of the phase of all vectors drawn from zeros of
BE
1+GH(close loop poles) and open loops on the s plane. If this point s is moved along a closed contour
enclosing any or all of theabove zeros and poles, only the phase of the vector of each of the enclosed
zeros or open-loop poles will change by 3600. The directionwill be in the same sense of the contour
enclosing zeros and in the opposite sense for the contour enclosing open-loop poles.
IN
GA
&
M
[Link] of argument
A
BE
When a closed contour in the s-plane encloses a certain number of poles and zeros of 1+G(s)H(s) in
the clockwise direction, the number ofencirclements of the origin by the corresponding contour in the
G(s)H(s)plane will encircle the point (-1,0) a number of times given by thedifference between the
number of its zeros and poles of 1+G(s)H(s) it enclosed on the s-plane.
IN
Modified contour on the s-plane forchecking the existence of closed-looppoles
s ej GA
&
A M
BE
Magnitude of GH remains the same alongthe contourPhase of changes from 270 to 90 degrees
Gain margin
M 20log G ( j p )H ( j p ) (14.7)
Phase margin
G( j g )H ( j g ) 180o (14.8)
IN
GA
&
A M
BE
[Link]
(1) Locate open-loop poles on the s-plane
(2) Draw the closed contour and avoid open-loop poles on the imaginary axis
(3) Count the number of open-loop poles enclosed in the above contour of step 2, say P
(4) Plot G(j )H(j ) and its reflection on the GH plane and map part of the small semi-circle
detour on the s-plane around poles (if any) on the imaginary axis.
(5) Once the entire s-plane contour is mapped on to the GH plane, count the number of
encirclements of the point (-1,0) and its direction. Clockwise encirclement is considered
positive, say N.
(6) The number of closed-loop poles in the right-half s-plane is given by Z=N+P. if Z >0, the
system is unstable.
(7) Determine gain margin, phase margin, and critical value of open-loop gain.
Example 14.1
Using Nyquist criterion, determine the stability of a feedback systemwhose open-loop transfer
function is given by
K
G (s) H ( s)
s ( s 1)
Solution
Step 1Locate open-loop poles on the s-plane. Open-loop poles are at s=0 and –1. Let K=1
Step 2 Draw the closed contour on the s-plane to check the existenceof closed-loop poles in the right-
half s-plane.
Open-loop poles and s-plane contour
IN
GA
&
1
M
G ( j )H ( )
2
1
1
G( j ) H ( j ) tan
2
A
IN
17 -0.4 2.3212 112 70 280
18 -0.2 4.9029 101 80 269
GA
&
A M
BE
The above system is stable. Here, phase crossover frequency is very large (infinity) and gain
crossover frequency 0.786 rad/s. Phase angle corresponding to gain crossover frequency= 2320and
o
Phase margin is 52
Example 14.2.
Using Nyquist criterion, determine the stability of a feedback systemwhose open-loop transfer
function is given by
55
G( s) H ( s)
s ( s 2)(s 4)
Solution
Step 1Locate open-loop poles on the s-plane. Open-loop poles are at s=0, -2 and –4. Let K=1
Step 2 Draw the closed contour on the s-plane to check the existenceof closed-loop poles in the right-
half s-plane.
Open-loop poles and s-plane contour
IN
The number of open-loop pole enclosed, P is zero
K
G( j )H ( j )
2 2
4 16
G( j ) H ( j )
2
tan 1
2
tan 1
4
GA
No. Phase, , s-plane,
Frequency Magnitude degrees deg
&
1 1.5 Positive 3.4332 213 270
frequencies
2 2 2.1741 198 280
3 2.5 1.4568 187 290
M
17 -6 0.201 218 70
18 -5.5 0.2513 214 80
19 -5 0.319 210 90
20 -4.5 0.4122 204 0
21 -3.5 0.7334 191 343
22 -3 1.017 183 326
23 -2.83 1.1446 180 309
24 -2.5 1.4568 173 292
IN
25 -2 2.1741 162 276
26 -1.5 3.4332 147 259
GA
&
A M
BE
Here, Z=N+P=2.
Hence, the above system is unstable.
Again,
Phase crossover frequency 2.83 rad/s
The gain at which the system becomes marginally stable, K * 55 / 1.1446 48
Gain margin
M 20 log G ( j p )H ( j p )
20log 1.1446 1.17 dB
o
Gain crossover frequency =3 rad/s and the corresponding angle of GH=177
o
Phase margin=177-180=-3
IN
Transfer function:
1
--------------------
s^4+9s^3+24s^2+16s
>>nyquist(sys)
GA
&
A M
BE
IN
MODULE#4 GA
Equation Chapter (Next) Section 1
&
A M
BE
3 dB is considered good
G( j ) x jy (16.3)
x2 y2
IN
M(j ) (16.4)
(1 x) 2 y2
M 2 (1 x) 2 M 2 y2 x2 y2 (16.5)
GA
x 2 (1 M 2 ) (1 M 2 ) y 2
M2 M2
2M 2 x M2 (16.6)
x2 y2 2 x (16.7)
1 M2 1 M2
2
M2
&
Adding in both sides, we get
1 M2
2 2
M2 2 M
x y (16.8)
1 M2 1 M2
M
M2 M
The above equation represents a family of circles with its center at 2
,0 and radius .
1 M 1 M2
A
Family of M-circles corresponding to the closeloop magnitudes (M) of aunit feedback system
Constant M-circles for unityfeedback systems
IN
15.4. Constant N-circles
M
GA G( j )
1 G( j )
(16.9)
1 y 1 y
tan tan (16.10)
x 1 x
&
1 y 1 y
N tan tan tan (16.11)
x 1 x
tan A tan B
tan( A B ) (16.12)
M
1 tan A tan B
N 2
(16.13)
x x y2
2 2 2
BE
1 1 1 1
x y (16.14)
2 2N 4 2N
1 1
The above equation represents a family of circles with its center at ( , ) and radius
2 2N
2
1 1
4 2N
IN
Example 15.1.
GA
Determine the closed-loop magnitude ratio and bandwidth of the feedback system whose forward
10
transfer function is given by G ( s ) and H(s)=1, by (1) direct computation and (2)
s ( s 2)( s 4)
&
using M and N circles.
Solution
A M
BE
IN
Closed-loopmagnitude and phase values
Closed-loop
Frequency, Closed-loop
Magnitude
rad/s
0.1
GA
ratio
1.0
Phase angle, deg
355
0.5 1.1 335
0.8 1.2 316
&
0.9 1.2 308
1.0 1.2 300
1.1 1.3 290
M
IN
GA
&
A M
BE
C( j )
Mr 20 log dB (16.16)
R( j )
C( j )
(16.17)
R( j )
IN
GA
&
A M
Example 15.2.
Determine the closed-loop magnitude ratio and bandwidth of the feedback system whose forward
BE
10
transfer function is given by G ( s ) and H(s)=1, by (1) direct computation and (2)
s ( s 2)( s 4)
using Nichols chart.
Solution
IN
GA
Closed-loop magnitude and phase values
IN
Open-loop magnitude and phase values
0.1 22
GA
magnitude, dB angle, deg
266
0.2 16 261
&
0.3 12 257
0.4 10 253
0.5 8 249
M
0.6 6 245
A
0.7 4 241
0.8 3 237
BE
0.9 2 233
1.0 1 229
1.1 0 226
1.2 -1 222
1.3 -2 219
1.4 -3 216
1.5 -4 213
1.6 -5 210
1.7 -6 207
1.8 -7 204
1.9 -7 201
2.0 -8 198
2.1 -9 196
IN
Equation Chapter (Next) Section 1
GA
&
A M
BE
16. Controllers
16.1. Basic Control Action and response of Control systems
An automatic controller compares the actual value of the plant output with the reference input
(desired value), determines the deviation, and produces a control signal that will reduce the
deviation to zero or to a small value. The manner in which the automatic controller produces the
control signal is called the control action. Fig.1 is a block diagram of an industrial control
system, which consists of an automatic controller, an actuator, a plant and a sensor (measuring
element). The controller detects the actuating error signal, which is usually at a low power level,
and amplifies it to a sufficiently high level. The output of the controller is fed to an actuator
such as pneumatic motor or valve, hydraulic motor or electric motor. The actuator is the device
that produces the input to the plant according to the control signal so that the output signal will
approach the reference input signal.
The sensor or measuring element is device that converts the output variable into another
IN
suitable variable such as a displacement, pressure or voltage that can be used to compare the
output to the reference input signal. This element is in the feedback path of the closed-loop
system. The set point of the controller must be converted to a reference input with the same
units as feedback signal from sensor.
GA
Error Detector
Sensor
M
Present
wages
Product
Initial
Industry cost K2 Cost of
wages
K1 living
Wages
increment
Dissatisfactionf
actor
IN
16.4. POLLUTION CONTROL IN AUTO ENGINE
To meet the emission standards for automobiles, Hydrocarbons (HC), Carbon monoxide (CO), and
Nitrogen oxides (NOx) emissions can be controlled by employing a three way catalyst in conjunction
with a closed loop engine control system as shown in Fig.2. The exhaust gas sensor gives an
GA
indication of a rich or lean exhaust and compares it to a reference. The difference signal is processed
by the controller, and the output of the controller modulates the vacuum level in the carburetor to
achieve the best air-fuel ratio for proper operation of the catalytic converter.
+
&
Reference Controller Three
_ Carburetor Engine way Exhaust
catalytic
converter
Sensor
M
Fig.16.3
A
Fig.16.4
(i) P-controller
(ii) PI-controller
(iii) PD-controller
(iv) PID-controller
P-controller
(a)
IN
GA
(b)
&
Fig.16.5
M
Fig.16.6
Kp
C (s) Js 2 Kp Kp
2 2 2
(16.1)
R( s ) Kp Js Kp J (s n )
1
Js 2
Kp
C ( s) 2 2
R(s ) (16.2)
J (s n )
1
For step-input, R s
s
c ( t ) 1 cos n t (16.3)
IN
Kp
Where, n
J
GA
&
A M
Fig.16.7
BE
Solved problem
1. Consider the unity feedback system of Fig. 16.8. Let Kp=20 and J=50. Determine the equation of
response for a unit step input and determine the steady-state error.
Fig.16.8
Solution
KP
C ( s) Js 2 Kp
2
R (s ) Kp Js Kp
1
Js 2
C ( s) Kp
2 2
R (s ) J s n
IN
Kp 2
n rad/s
J 5
Kp
c(t )
J 2
n
(1 cos n t) GA
2
c(t ) 1 cos t
5
&
2 2
se ( t ) 1 1 cos t cos t
5 5
M
Fig.16.9
Solution
Kp
G1
s 1
E1 ( s ) 1 s 1
R( s) 1 G1 ( s ) s 2
1 s 1
E1 ( s ) R( s ) R(s)
1 G1 ( s ) s 2
Step response is
s 1
E1 ( s )
s s 2
U1 ( s ) K p ( s 1)
R( s ) ( s 2)
IN
C1 ( s ) 1 1
R( s ) s 1 Kp s 2
e1 (t )
1
2
1 e 2t GA
1
c1 (t ) 1 e 2t
2
&
A M
BE
Fig.16.10
I-controller
(a)
C2
R1
ei -
eo
+
IN
GA (b)
Fig.16.11
&
Ki
G2
s ( s 1)
E2 ( s ) 1 1 s( s 1)
M
R( s ) 1 G2 (s ) 1
Ki s2 s 1
s ( s 1)
A
U 2 (s) K i ( s 1)
R(s ) s2 s 1
BE
C2 ( s ) 1 1
2 2
R( s ) s s Ki s s 1
0.5 t 3 1 3
e2 (t ) e cos t sin t
2 3 2
0.5t 1 3 3
c2 (t ) 1 e sin t cos t
3 2 2
Fig.16.12
IN
PD-controller
GA
(a)
&
R2
R
R1 R
M
ei -
-
+
eo
+
A
C
BE
(b)
Fig.16.13
Fig.16.14
K p (1 Td s)
C ( s) Js 2 K p (1 Td s) K p (1 Td s)
2
(16.4)
R( s ) K p (1 Td s) Js K p (1 Td s) K pTd Kp
1 J s2 s
Js 2 J J
K p (1 Td s )
C ( s) R( s ) (16.5)
2
K pTd Kp
J s s
J J
1
For step-input, R s
s
IN
Step response becomes
Kp nt
c (t ) 2
1 e cos d t sin d t
J n GA 1 2
(16.6)
K pTd 1 nt 2
e sin n 1 t
J n 1 2
Kp
&
Where, n
J
Solved problem
M
3. Consider the unity feedback system of Figure 3. Let K p=20 and J=50. Determine the equation of
response for a unit step input and determine the steady-state error. Here, Kp =20, Td =1 and J=50.
A
BE
Fig.16.14
Solution
C ( s) 20( s 1) 20(1 s)
R( s) 50 s 2 20 s 20 2
50( s 2 n s 2
n )
Kp nt
c(t ) 2
1 e cos d t sin d t
J n 1 2
K pTd 1 nt 2
e sin n 1 t
J n 1 2
Kp
c( ) 2
1
J n
Only PD control
Transient characteristic
No system damping
IN
35.09
Maximum overshoot, %
unsatisfactory
Rise time tr, sec 3.15
Peak time tp, sec 5.24
Settling time ts , sec (5% criterion)
GA 15
PI-controller
&
M
(a)
A
BE
(b)
Fig.16.15
Fig.16.16
Ki 1 sK p Ki
G Kp (16.7)
IN
s s 1 s ( s 1)
E (s) 1 s ( s 1)
2
(16.8)
R( s ) 1 G1 ( s ) s s (1 K p ) K i
U (s)
R ( s) s 2
GA
(s 1)( sK p
s (1 K p ) Ki
Ki )
1 (16.9)
C ( s) ( sK p Ki ) 1
2
(16.10)
R(s) s s (1 K p ) K i (s 1)
&
Step response
u (t ) 1
M
t
e (t ) e
A
c (t ) 1 e t
BE
Fig.16.17
Fig.16.18
IN
GA
&
Fig.16.19
M
PID-controller
A
BE
(a)
(b)
Fig.16.20
Ki
Gc s Kp Kd s (16.11)
s
Ti
Gc s K p 1 Td s (16.12)
s
Where,
Ki
Ti
Kp
(16.13)
Kd
Td
Kp
IN
Tuning of PID-controller
Fig.16.21
A
C (s) Ke Ls
(16.14)
BE
U (s ) 1 Ts
Ke Ls
C ( s) U ( s) (16.15)
1 Ts
IN
Fig.16.22
Type of controller
P T
Kp
GA Ki
0
Kd
0
L
PI 0.9T L 0
&
L 0.3
PID T 2L 0.5L
1.2
L
M
1
Gc ( s ) Kp 1 Td s
Ti s
A
1.2T 1
1 0.5 Ls
L 2 Ls
BE
2
1
s
L
Gc ( s) 0.6T
s
Ziegler-Nichols tuning rule based on critical gain Kcr and critical period Pcr.
Type of controller Kp Ki Kd
P 0.5 Kcr 0 0
PI 0.45 Kcr 1/1.2 Pcr 0
PID 0.6 Kcr 1/0.5Pcr 0.125 Pcr
Where, Kcr proportional constant of a switched-off integral and derivative controls at which sustained
oscillations of period P cr occur.
Second Method
Ti
Gc ( s ) Kp 1 Td s
s
(16.16)
1
0.6 K cr 1 0.125 Pcr s
0.5 Pcr s
2
4
s
Pcr
Gc ( s ) 0.075 K cr Pcr (16.17)
IN
s
GA
&
A M
BE
17. Components
17.1. AC SERVOMOTORS
A two phase servomotor (Induction Motor) (A few watts to hundred watts) is commonly used
in feedback control systems. In servo applications, an induction motor is required to produce
rapid accelerations from standstill.
Schematic Diagram
IN
Constructional features GA
Squirrel Cage rotor with Cu or Al conductor
High Rotor resistance
Small diameter to length ratio to minimize inertia
Two stator windings in space quadrature(One called reference winding and the other Control
winding)
&
The two voltages to stator windings must derived from same source(Or they must be in
synchronism)
Principle of Operation
M
(i) The two applied AC voltage to stators with a phase difference produce a rotating flux.
(ii) As this moving flux sweeps over the rotor conductors, small emf is induced in rotor. Rotor being
short circuited, currents will flow and this current interacts with rotating flux to produce a torque
in the rotor. This torque causes the rotor to turn so that it chases the rotating magnetic flux.
A
BE
IN
Transfer Function
GA
The torque developed is a function of shaft angular position ( ) and control voltage Ec.
s Km Km K
G s
&
2
Ec s Js Ds s Js D s Tm s 1
Km J
Where, K = motor gain constant, Tm = motor time constant
D D
M
Merits of AC Servomotors
(i) Lower cost,
A
Demerits of AC Servomotors
(i) Nonlinear characteristics,
(ii) Used for low power applications(e.g. instrument servo),
(iii) Difficult for speed control and positioning
17.2. Synchros
It is also known as selsyn. It is a self-synchronizing device widely used in servomechanisms as
a position indicator.
IN
GA
&
IN
(c) Schematic Symbol
Fig. Synchro Transmitter
Synchro transmitter
GA
It is not a three phase machine. It is a single phase machine. Here, input is angular position of
its rotor [Link] is a set of three stator coil-to-coil voltages. Common connection
between the stator coils is not accessible.
&
Synchro system with transmitter and control transformer
A synchro error detector system may consist of synchro transmitter and synchro control
transformer.
It compares two angular displacements and the output voltage is approximately linear
M
circuit configuration, using two synchros, for maintaining synchronism between two shafts.
When rotor windings are excited, emfs are induced by transformer action in the stator windings
of transmitter and motor. If the two shafts are in similar positions (relative to that of the stator
windings), then there are two emfs of equal value are induced in the two stator [Link]
no circulating current exists and hence no torque is produced. If the two shaft positions do not
match, the emfs are unequal and result circulating current to flow. The circulating current in
conjunction with air gap magnetic field produce torque which tend to align the shafts.
Then, r t dt
If the phase sequence of stator and rotor windings of differential are reversed then
r t dt
IN
17.3.1. DC Tachometer
It is a permanent magnet DC generator. It resembles a small DC machine having a PM
stator, rotating armature, brush and commutator assembly. The rotor is connected to the
shaft to be measured. The output voltage is proportional to the angular velocity of the
GA
shaft. The direction of rotation decides the polarity of the voltage.
&
A M
BE
17.3.2. AC Tachometer
Used in AC servomechanism. It resembles 2-phase AC induction motor.
It comprises two stator windings arranged in space quadrature and a rotor which is
not conductively connected to external circuit.
One stator phase winding is excited by a suitable AC voltage of constant
magnitude and frequency. A voltage of the same frequency is generated across the
other winding known as control winding.
It is necessary that the voltage developed across the control winding is linearly
proportional to shaft speed and the phase of this voltage be fixed with respect to
voltage applied to reference winding.
The output voltage is connected to high impedance circuit of amplifier so that the
winding is considered open circuit.
An AC tacogenerator should have low inertia when rapid speed variations are
encountered. The drag cup construction gives low inertia and is used many times
IN
GA
&
A M
BE
Principle of Operation
If input x moves the pilot valve to the right, port II is uncovered, and so high pressure oil
enters the right side of the power piston. Since port I is connected to the drain port, the
oil in the left side of the power piston is returned to the drain. The oil flowing into the
power cylinder is at high pressure; the oil flowing out from the power cylinder into the
IN
drain is at low pressure. The resulting difference in pressure on both sides of the power
piston will cause it to move to the left.
GA
&
A M
Transfer Function
BE
Rate of flow of fluid Q(kg/sec) time dt(sec) is equal to the power piston displacement
dy(m) times the piston area A(sq.m) times the density of fluid P (kg/m3 ).Fluid flow rate
is proportional to pilot valve displacement x. So,
Q x
Q Kx
dy
Q AP
dt
dy
AP Kx
dt
APsY s KX s
Y s K
X s APs
Advantages
IN
Disadvantages
instrument.
ii. The actuator should have
Desired dynamic response
Desired cost, size, and weight
A
iii. Electric power is readily available, cleaner and quieter and easier to transmit.
So electric motors is mostly preferable compared to hydraulic and pneumatic actuation)
BE
Merits of DC motor
(i) Linear characteristics,
(ii) Used for large power applications,
(iii) Easier control
Demerits of DC motor
(i) Lower torque to volume and
(ii) Lower torque to inertia ratio.
Future developments
(i) Development of rare earth magnet results in DC motor high torque to volume ratio.
(ii) Advances in brush commutator technology make trouble free maintenance.
(iii) Development of brushless DC motors.
SUMMARY
IN
3.0Control System Components
3.1 Potentiometers
3.3 Tachometers
GA
3.4 Amplidyne
IN
6.5 Determination of Damping ratio, Gain Margin and Phase Margin from Root Locus
IN
GA
&
A M
BE