EE2253 - Control Systems PDF
EE2253 - Control Systems PDF
com powers Engineers, Doctors, Managers & Lawyers in India by providing 'free'
resources for aspiring students of these courses as well as students in colleges.
You get FREE Lecture notes, Seminar presentations, guides, major and minor projects.
Also, discuss your career prospects and other queries with an ever-growing community.
ALL FILES ON STUDYNAMA.COM ARE UPLOADED BY RESPECTIVE USERS WHO MAY OR MAY NOT BE THE OWNERS OF THESE FILES. FOR ANY SUGGESTIONS OR FEEDBACK, EMAIL US AT [email protected]
L T PC
EE2253 CONTROL SYSTEMS 3104
(Common to EEE, EIE & ICE)
2
UNIT I
SYSTEMS AND THEIR REPRESENTATION
Basic elements in control systems – Open and closed loop systems – Electrical
analogy of Mechanical and thermal systems – Transfer function – Synchros – AC
and DC servomotors –Block diagram reduction techniques – Signal flow graphs.
Control System
Definitions:
Controller: Controller means measuring the value of the controlled variable of the system &
applying the manipulated variable to the system to correct or to limit the deviation of
the measured value to the desired value.
Disturbance: A disturbance is a signal that tends to affect the value of the output of a
system. If disturbances created inside the system, it is called internal While an
external disturbance is generated outside the system.
Open loop System: In an Open loop System, the control action is independent of the
desired output. OR When the output quantity of the control system is not fed back to
the input quantity, the control system is called an Open loop System.
Closed loop System: In the Closed loop Control System the control action is
dependent on the desired output, where the output quantity is considerably controlled
by sending a command signal to input quantity
3
Classification of Control Systems
(a) Method of analysis & design, as Linear & Non- Linear Systems.
(b)The type of the signal, as Time Varying, Time Invariant, Continuous data, Discrete data
systems etc.
In a linear system, the principle of superposition can be applied. In non- linear system,
principle & homogeneity.
While operating a control system, if the parameters are unaffected by the time, then
the system is called Time Invariant Control System.
Most physical systems have parameters changing with time. If this variation is
measurable during the system operation then the system is called Time Varying
System.
If there is no non-linearity in the time varying system, then the system may be called
as Linear Time varying System.
If the signal is not continuously varying with time but it is in the form of pulses.
Thenthe control system is called Discrete Data Control System.
If the signal is in the form of pulse data, then the system is called Sampled Data
Control System.
Here the information supplied intermittently at specific instants of time. This has the
advantage of Time sharing system.
On the other hand, if the signal is in the form of digital code, the system is called
Digital Coded System.
Here use of Digital computers, µp, µc is made use of such systems are analyzed by
the Z- transform theory.
If the signal obtained at various parts of the system are varying continuously with time,then
the system is called Continuous Data Control Systems.
4
Adaptive Control systems:
In some control systems, certain parameters are either not constant or vary in an
unknown manner.
If the parameter variations are large or rapid, it may be desirable to design for the
capability of continuously measuring them & changing the compensation, so that the
system performance criteria can always satisfied.
In simple control system there will be One input & One output such systems are calledSingle
variable System (SISO Single Input & Single Output).
In Multivariable control system there will be more than one input & correspondingly more
output - Multiple Inputs & Multiple Outputs)
5
CONTROL Output/
Input/
SYSTEM Result
Objective
Open-loop systems:
The open-loop system is also called the non-feedback system. This is the simpler of
the two systems. A simple example is illustrated by the speed control of an automobile as
shown in Figure. In this open-loop system, there is no way to ensure the actual speed is
close to the desired speed automatically. The actual speed might be way off the desired
speed because of the wind speed and/or road conditions, such as uphill or downhill etc.
Closed-loop systems:
The closed-loop system is also called the feedback system. A simple
closed-system is shown in Figure. It has a mechanism to ensure the actual speed is close
to the desired speed automatically
+
Desired Controller Engine Actual
Speed Speed
_
Speed
Sensor
An open loop system has the ability A closed loop system has got the
to perform accurately, if its ability to perform accurately because
calibration is good. If the calibration of the feed back
is not perfect its performance will go
down
6
In general it is more stable as the Less Stable Comparatively
feedback is absent
If non- present; the system operation Even under the presence of non loop
is not good system
Transfer Function
A simpler system or element maybe governed by first order or second order
differential equation
When several elements are connected in sequence, say ―n‖ elements, each one with
first order, the total order of the system will be nth order
In general, a collection of components or system shall be represented by nth order
differential equation
7
It is applicable to Linear Time Invariant system.
It is the ratio between the Laplace Transform of the o/p variable to the
Laplace Transform of the i/p variable.
An equation describing the physical system has integrals & differentials, the step involved in
obtaining the transfer function are;
d
(2) Replace the terms
dt
by ‗S‘ & dt by 1/S.
(3) Eliminate all the variables except the desired variables.
Once a physical model is obtained, the next step is to obtain Mathematical model.
When a mathematical model is solved for various i/p conditions, the result represents
the dynamic behavior of the system.
Analogous System:
The concept of analogous system is very useful in practice. Since one type of system may be
easier to handle experimentally than another. A given electrical system consisting of
resistance, inductance & capacitances may be analogous to the mechanical system consisting
of suitable combination of Dash pot, Mass & Spring.
An Electrical engineer familiar with electrical systems can easily analyze the system
under study & can predict the behavior of the system.
The model of mechanical translational systems can obtain by using three basic elements
mass, spring and dash-pot. When a force is applied to a translational mechanical system, it is
opposed by opposing forces due to mass, friction and elasticity of the system. The force
acting on a mechanical body is governed by Newton‘s second law of motion. For
translational systems it states that the sum of forces acting on a body is zero.
Translational System:
9
Force balance equations of idealized elements
Consider an ideal mass element shown in fig. which has negligible friction and elasticity.
Let a force be applied on it. The mass will offer an opposing force which is proportional to
acceleration of a body.
d 2x
Here fm M
dt 2
d 2x
By Newton‘s second law, f = f m= M
dt 2
2. Damping Force [Viscous Damping]: Due to viscous damping, it is proportional to
velocity & is given by,
Consider an ideal frictional element dash-pot shown in fig. which has negligible mass and
elasticity. Let a force be applied on it. The dashpot will be offer an opposing force which is
proportional to velocity of the body
10
Let f = applied force
f b = opposing force due to friction
dx
Here, f b B
dt
dx
By Newton‘s second law, f = f b B
dt
3. Spring Force: Spring force is proportional to displacement.
Consider an ideal elastic element spring shown in fig. which has negligible mass and friction.
The model of rotational mechanical systems can be obtained by using three elements,
11
moment of inertia [J] of mass, dash pot with rotational frictional coefficient [B] and torsional
spring with stiffness[k].When a torque is applied to a rotational mechanical system, it is
opposed by opposing torques due to moment of inertia, friction and elasticity of the system.
The torque acting on rotational mechanical bodies is governed by Newton‘s second law of
motion for rotational systems.
Consider an ideal mass element shown in fig. which has negligible friction and
elasticity. The opposing torque due to moment of inertia is proportional to the angular
acceleration.
12
d
T= Tb = B (1 2 )
dt
Consider an ideal elastic element, torsional spring as shown in fig. which has negligible
moment of inertia and friction. Let a torque be applied on it. The torsional spring will offer an
opposing torque which is proportional to angular displacement of the body.
Mechanical Mechanical
Sl.No. F.VAnalogy
Translational System Rotational System F.IAnalogy
13
Modeling of electrical system
Electrical circuits involving resistors, capacitors and inductors are considered. The
behaviour of such systems is governed by Ohm‘s law and Kirchhoff‘s laws
Resistor: Consider a resistance of ‗R‘ carrying current ‗i‘ Amps as shown in Fig
(a), then the voltage drop across it is v = R I
Inductor: Consider an inductor ―L‘ H carrying current ‗i‘ Amps as shown in Fig (a),
then the voltage drop across it can be written as v = L di/dt
Capacitor: Consider a capacitor ―C‘ F carrying current ‗i‘ Amps as shown in Fig (a),
then the voltage drop across it can be written as v = (1/C) i dt
Apply Kirchhoff‘s voltage law or Kirchhoff‘s current law to form the differential
equations describing electrical circuits comprising of resistors, capacitors, and
inductors
Form Transfer Functions from the describing differential equations
Then simulate the model
Example
14
t
1
c 0
R1i(t ) R2 i (t ) i (t )dt v1 (t )
t
1
R2 i (t ) i (t )dt v2 (t )
c0
Electrical systems
LRC circuit. Applying Kirchhoff‘s voltage law to the system shown. We obtain the
following equation;
Resistance circuit
Equation (1) & (2) give a mathematical model of the circuit. Taking the L.T. of equations
(1)&(2), assuming zero initial conditions, we obtain
Armature-Controlled dc motors
15
They are either armature-controlled with fixed field or field-controlled with fixed
armature current.
For example, dc motors used in instruments employ a fixed permanent-magnet field,
and the controlled signal is applied to the armature terminals.
kg.m2
f = equivalent viscous-friction coefficient of the motor and load referred to the motor
shaft. Newton*m/rad/s
T = k1 ia ψ where ψ is the air gap flux, ψ = kf if, k1 is constant
For the constant flux
16
The differential equation for the armature circuit
The armature current produces the torque which is applied to the inertia and friction; hence
Assuming that all initial conditions are condition are zero/and taking the L.T. of
equations (1), (2) & (3), we obtain
Kpsθ(s) = Eb(s)
17
The T.F can be obtained is
Field-Controlled dc motors
kg.m2
f = equivalent viscous-friction coefficient of the motor and load referred to the motor
shaft. Newton*m/rad/s
18
T = k1 ia ψ where ψ is the air gap flux, ψ = kf if , k1 is constant
In the field controlled DC motor, the armature current is fed from a constant current source.
On Laplace Transform,
On Laplace Transform,
19
The complete block diagram is as shown below,
Analogous Systems
Let us consider a mechanical (both translational and rotational) and electrical system as
shown in the fig.
20
From the fig (b)
d 2 d
We get J 2
D K T -------------- (2)
dt dt
From the fig (c)
d 2q dq 1
We get L 2
R q v(t ) ------------ (3)
dt dt C
Where q i dt
21
PROBLEM
1. Find the system equation for system shown in the fig. And also determine f-v and f-i
analogies.
22
From eq (2) we get
Force–current analogy
23
The system can be represented in two forms:
Block diagram
A pictorial representation of the functions performed by each component and of the
flow of signals.
It is a representation of the control system giving the inter-relation between the
transfer function of various components.
The block diagram is obtained after obtaining the differential equation & Transfer
function of all components of a control system.
A control system may consist of a number of components.
A block diagram of a system is a pictorial representation of the functions performed
by each component and of the flow of signals.
The elements of a block diagram are block, branch point and summing point.
The arrow head pointing towards the block indicates the i/p
& pointing away from the block indicates the o/p.
After obtaining the block diagram for each & every component, all blocks are combined to
obtain a complete representation. It is then reduced to a simple form with the help of block
diagram algebra.
24
Basic elements of a block diagram
Blocks
Transfer functions of elements inside the blocks
Summing points
Take off points
Arrow
Block diagram
+ G(s)
―
-
Block
In a block diagram all system variables are linked to each other through functional
blocks. The functional block or simply block is a symbol for the mathematical operation on
the input signal to the block that produces the output.
Summing point
25
A summing point has only one output and is equal to the algebraic sum of the
inputs.
A takeoff point is used to allow a signal to be used by more than one block or summing
point.
The transfer function is given inside the block
• The input in this case is E(s)
• The output in this case is C(s)
• C(s) = G(s) E(s)
Functional block – each element of the practical system represented by block with its
T.F.
Branches – lines showing the connection between the blocks
Arrow – associated with each branch to indicate the direction of flow of signal
Closed loop system
Summing point – comparing the different signals
Take off point – point from which signal is taken for feed back
Advantages of Block Diagram Representation
Very simple to construct block diagram for a complicated system
Function of individual element can be visualized
Individual & Overall performance can be studied
Over all transfer function can be calculated easily
26
Disadvantages of Block Diagram Representation
No information about the physical construction
Source of energy is not shown
Because of their simplicity and versatility, block diagrams are often used by control
engineers to describe all types of systems.
Also, it can be used, together with transfer functions, to represent the cause-and-
effect relationships throughout the system.
Transfer Function is defined as the relationship between an input signal and an output
signal to a device.
27
Cascaded blocks
28
Eliminating a feedback loop
Cascaded Subsystems
Parallel Subsystems
29
Procedure to solve Block Diagram Reduction Problems
Step 1: Reduce the blocks connected in series
Step 4: Try to shift take off points towards right and Summing point towards left
PROBLEM
Solution:
By eliminating the feed-back paths, we get
30
Combining the blocks in series, we get
Transfer function
31
Solution:
32
Transfer function
3. Obtain the transfer function for the system shown in the fig
4. Obtain the transfer function C/R for the block diagram shown in the fig
33
Solution
Transfer function
34
Signal Flow Graph Representation
Signal Flow Graph Representation of a system obtained from the equations, which
shows the flow of the signal
For complicated systems, Block diagram reduction method becomes tedious &
time consuming. An alternate method is that signal flow graphs developed by S.J.
Mason.
In these graphs, each node represents a system variable & each branch connected
between two nodes acts as Signal Multiplier.
Definitions
Branch: A branch is a line joining two nodes. The signal travels along a branch.
35
Input node [Source]: It is a node which has only out going signals.
Mixed node: It is a node which has both incoming & outgoing branches (signals).
Path: It is the traversal of connected branches in the direction of branch arrows. Such
that no node is traversed more than once.
Non-Touching Loops: Loops are Non-Touching, if they do not possess any common
node.
Forward Path: It is a path from i/p node to the o/p node w hich doesn‘t cross any
node more than once.
PROBLEM
36
Solution:
No of forward paths
37
(2) Obtain the closed loop TF, by using Mason‘s gain formula.
Solution:
Forward Paths:
38
3. Obtain C(s)/R(s) by Block Diagram Reduction method & verify the result by signal flow
graph
Solution:
39
Signal flow graphs:
A feedback control system usually consists of several components in addition to the actual
process. These are error detectors, power amplifiers, actuators, sensors etc. Let us now
discuss the physical characteristics of some of these and obtain their mathematical models.
DC Servo Motor
40
DC servo motors have a high ratio of starting torque to inertia and therefore they have
a faster dynamic response.
DC motors are constructed using rare earth permanent magnets which have high
residual flux density and high coercively.
As no field winding is used, the field copper losses am zero and hence, the overall
efficiency of the motor is high.
The speed torque characteristic of this motor is flat over a wide range, as the armature
reaction is negligible.
Moreover speed in directly proportional to the armature voltage for a given torque.
Armature of a DC servo motor is specially designed to have low inertia.
In some application DC servo motors are used with magnetic flux produced by field
windings.
The speed of PMDC motors can be controlled by applying variable armature voltage.
These are called armature voltage controlled DC servo motors.
Wound field DC motors can be controlled by either controlling the armature voltage
or controlling rho field current. Let us now consider modelling of these two types or
DC servo motors.
The field is produced either by a permanent magnet or the field winding is separately
excited and supplied with constant voltage so that the field current If is a constant.
When the armature is supplied with a DC voltage of e a volts, the armature rotates and
produces a back e.m.f eb.
41
The armature current ia depends on the difference of eb and en. The armature has a
remnant of inertia J, frictional coefficient B0
The back emf is proportional to the speed of the motor and hence
42
The block diagram representation of the armature controlled DC servo motor is developed in
steps
Usually the inductance of the armature winding is small and hence neglected
Where
43
The electrical circuit is modeled as
Where
44
AC Servo Motors
It has two windings displaced by 90 on the stator One winding, called as reference
winding, is supplied with a constant sinusoidal voltage.
The second winding, called control winding, is supplied with a variable control
voltage which is displaced by -- 90 out of phase from the reference voltage.
The major differences between the normal induction motor and an AC servo motor are
The rotor winding of an ac servo motor has high resistance (R) compared to its
X
inductive reactance (X) so that its R ratio is very low.
X
For a normal induction motor, R ratio is high so that the maximum torque is obtained
in normal operating region which is around 5% of slip.
The torque speed characteristics of a normal induction motor and an ac servo motor
are shown in fig
The Torque speed characteristic of a normal induction motor is highly nonlinear and
has a positive slope for some portion of the curve.
This is not desirable for control applications. as the positive slope makes the systems
unstable. The torque speed characteristic of an ac servo motor is fairly linear and has
45
negative slope throughout.
The rotor construction is usually squirrel cage or drag cup type for an ac servo motor.
The diameter is small compared to the length of the rotor which reduces inertia of the
moving parts.
The supplies to the two windings of ac servo motor are not balanced as in the case of
a normal induction motor.
The control voltage varies both in magnitude and phase with respect to the constant
reference vulture applied to the reference winding.
The direction of rotation of the motor depends on the phase (± 90°) of the control
voltage with respect to the reference voltage.
For different rms values of control voltage the torque speed characteristics are shown
in Fig.
The torque varies approximately linearly with respect to speed and also controls
voltage.
The torque speed characteristics can be linearised at the operating point and the
transfer function of the motor can be obtained.
Synchros
Synchro transmitter
46
The principle of operation of these two devices is sarne but they differ slightly in their
construction.
The stator consists of a balanced three phase winding and is star connected.
The rotor is of dumbbell type construction and is wound with a coil to produce a
magnetic field.
When a no voltage is applied to the winding of the rotor, a magnetic field is produced.
The coils in the stator link with this sinusoidal distributed magnetic flux and voltages
are induced in the three coils due to transformer action.
Than the three voltages are in time phase with each other and the rotor voltage.
The magnitudes of the voltages are proportional to the cosine of the angle between the
rotor position and the respective coil axis.
The position of the rotor and the coils are shown in Fig.
When 90 the axis of the magnetic field coincides with the axis of coil S2 and
maximum voltage is induced in it as seen.
47
For this position of the rotor, the voltage c, is zero, this position of the rotor is known
as the 'Electrical Zero' of die transmitter and is taken as reference for specifying the
rotor position.
In summary, it can be seen that the input to the transmitter is the angular position of
the rotor and the set of three single phase voltages is the output.
The magnitudes of these voltages depend on the angular position of the rotor as given
Hence
Now consider these three voltages to he applied to the stator of a similar device called
control transformer or synchro receiver.
Since the rotor is cylindrical, the air gap is uniform and the reluctance of the magnetic
path is constant.
Usually the rotor winding of control transformer is connected teas amplifier which
requires signal with constant impedance for better performance.
Since the some currents flow through the stators of the synchro transmitter and
receiver, the same pattern of flux distribution will be produced in the air gap of the
control transformer.
48
The control transformer flux axis is in the same position as that of the synchro
transmitter.
Thus the voltage induced in the rotor coil of control transformer is proportional to the
cosine of the angle between the two rotors.
**************************
UNIT I
49
PART – A
2. Determine the transfer function Y2(S)/F(S) of the system shown in fig. (16)
50
3. Write the differential equations governing the Mechanical rotational system shown in
fig. Draw the Torque-voltage and Torque-current electrical analogous circuits. (16)
4. Determine the overall transfer function C(S)/R(S) for the system shown in fig. (16)
5. For the system represented by the block diagram shown in fig. Determine C1/R1 and
C2/R1. (16)
6. Find the overall gain of the system whose signal flow graph is shown in fig. (16)
51
7. Draw a signal flow graph and evaluate the closed loop transfer function of a system
Whose block is shown in fig?
8. Derive the transfer function for Armature controlled DC servo motor. (16)
9. Derive the transfer function for Field controlled DC servo motor. (16)
UNIT II
TIME RESPONSE
Time response – Time domain specifications – Types of test input – I and II order
system response – Error coefficients – Generalized error series – Steady state error –
P, PI, PID modes of feedback control
52
Introduction
The next step would be, to obtain its response, both transient and steady state, to a
specific input.
The input can be a time varying function which may be described by known
mathematical functions or it may be a random signal.
Thus it is customary to subject the control system to some standard input test signals
which strain the system very severely.
These standard input signals arc: an impulse. a step, a ramp and a parabolic input.
Analysis and design of control systems are carried out, defining cream performance
measures for the system, using these standard test signals.
It is also pertinent to mention that any arbitrary time function can be expressed in
terms of linear combinations of these test signals and hence, if the system is linear, the
output of the system can be obtained easily by using superstition principle.
Further, convolution integral can also be used to determine the response of a linear
system for any given input, if the response is known for a step or an impulse input.
Normally use the standard input signals to identify the characteristics of system‘s
response
Step function
Ramp function
Impulse function
53
Parabolic function
Sinusoidal function
Time response analysis
It is an equation or a plot that describes the behavior of a system and contains much
information about it with respect to time response specification as overshooting, settling time,
peak time, rise time and steady state error. Time response is formed by the transient response
and the steady state response.
Transient time response (Natural response) describes the behavior of the system in its
first short time until arrives the steady state value and this response will be our study
focus. If the input is step function then the output or the response is called step time
response and if the input is ramp, the response is called ramp time response ... etc.
Classification of Time Response
Transient response
Steady state response
y(t) = yt(t) + yss(t)
Transient Response
The transient response is defined as the part of the time response that goes to zero as
time becomes very large. Thus yt(t) has the property
Lim yt(t) = 0
t -->
The time required to achieve the final value is called transient period. The transient
response may be exponential or oscillatory in nature. Output response consists of the sum of
forced response (form the input) and natural response (from the nature of the system).The
transient response is the change in output response from the beginning of the response to the
final state of the response and the steady state response is the output response as time is
approaching infinity (or no more changes at the output).
54
Steady State Response
The steady state response is the part of the total response that remains after the
transient has died out. For a position control system, the steady state response when
compared to with the desired reference position gives an indication of the final accuracy of
the system. If the steady state response of the output does not agree with the desired reference
exactly, the system is said to have steady state error.
Typical Input Signals
Impulse Signal
Step Signal
Ramp Signal
Parabolic Signal
Impulse Signal
The impulse function is zero for all t 0 and its infinity at t 0 .It rises to infinity at
t 0 and comes uback to zero at t 0 enclosing a finite area. If this area is A it is called as
impulse function of strength A. If A 1 it is called a unit impulse function. Thus an impluse
signal is denoted by f (t ) A (t ) .
step Signal
55
It is zero for t 0 and suddenly rises to a value A at t=0 and remains at this value for t>0: it is
denoted by f (t ) A u(t ) . If A=1, it is called a unit step function.
Ramp Signal
It is zero for t<0 and uniformly increases with a slope equal to A. It is denoted by
f (t ) A t .If the slope is unity then it is called a unit ramp signal.
Parabolic Signal
At 2
A parabolic signal is denoted by f (t ) If A is equal to unity it is known as a unit
2
parabolic signal.
56
Step response K K Step, exponential
-
s s 1/ T
Order of a system
The Order of a system is given by the order of the differential equation governing the
system
Alternatively, order can be obtained from the transfer function
In the transfer function, the maximum power of s in the denominator polynomial
gives the order of the system
Dynamic Order of Systems
Order of the system is the order of the differential equation that governs the dynamic
behaviour
Working interpretation: Number of the dynamic elements / capacitances or holdup
elements between a manipulated variable and a controlled variable
Higher order system responses are usually very difficult to resolve from one another
The response generally becomes sluggish as the order increases
1
Consider a feedback system with G( s) as shown in fig
s
57
---------------------------(1)
1
For a unit step input R ( s ) and the output is given by
s
---------------------------(2)
---------------------- (3)
This is 63.2 percent of the steady value. This time is known as the here constant of
the system. One of the characteristics which we would like to know about the system
is its speed of response or how fast the response is approaching the final value. The
time constant is indicative of this measure and the speed of response is inversely
proportional to the time constant of the system.
Another important characteristic of the system is the error between the desired value
and the actual value under steady state conditions. This quantity is known as the
steady state error of the system and is denoted by e ss
E(s) =R(s)-C(s)
The response of the system in Fig.for a unit ramp input, for which.
is given by,
which is the response of the system to a step input. Thus no additional infommtion
about the speed of response is obtained by considering a ramp input. But let us see
the effect on the steadystate error As before,
Thus the steady state error is equal to the time constant of the system. The first order system,
therefore, cannot track the ramp input without a finite steady state error. If the time constant
is reduced not only the speed of response increases but also the steady state error for ramp
59
input decreases. Hence the ramp input is important to the extent that it produces a finite
steady state error. Instead of finding the entire response, it is sufficient to estimate the steady
state value by using the final value theorem. Thus
The response of a first order system for unit ramp input is plotted in Fig.
The response of a first order system to a unit parabolic input, for which
which is the response of the first order system to unit velocity input. Thus subjecting the first
order system to a unit parabolic input does not give any additional information regarding
transient behavior of the system. But, the steady state error, for a parabolic input is given by.
60
Thus a first order system has infinite state error for a parabolic input .The steady state error
can be easily obtained by using the final value theorem as:
Summarizing the analysis of first order system. we can say that the step input yields the
desired information about the speed of transient response. Isis observed that the speed of
response is inversely proportional lathe time constant r of the system. The ramp and parabolic
inputs do not give any additional information regarding the speed of response. However, the
steady statesman is different for these three different inputs. For a step input, the steady state
error e0 is zero, for a velocity input there is a finite error equal to the time constant r of the
system and for an acceleration input the steady state error is infinity.
It is clear from the discussion above, that it is sufficient to study the behavior of any system
to a unit step input for understanding its transient response and use the velocity input and
acceleration input fur understanding the steady state behavior of the system.
61
Second order system
Consider a Type second order system as shown in Fig. Since G(s) has Pole at the origin,
The transient response of any system depends on the poles of the transfer function T(s).the
roots of the denominator polynomial in s of T(s) are the poles of the transfer function. Thus
the denominator polynomial of T(s) is given by
Is known as the characteristic polynomial of the system and D(S)=0is known as the
characteristic equation of the system.
--------------(1)
Where
62
K
n = natural frequency
1
= damping factor
2 K
s 2 2n s n 0
2
2n s 4 2 n 4 n
2 2
s1, 2
2
Are given by n j n 1 2
n j d
Where d n 1 2 is known as the damped natural frequency of the system. If >1, the
roots s1, s2 are real and we have an over damped system. If =1 the system is known as a
critically damped system. The more common case of <1 is known as the under damped
system.
For a unit step input R(s) =1/s and equation (1) can be written as
C(s) =T(s).R(s)
n 2 1
. --------------- (2)
s 2n s n s
2 2
K1 K 2 s K3
2
s s 2n s n 2
1 s 2n
C ( s)
s ( s n ) 2 n 2 (1 2 )
1 s n n (1 2 )
s ( s n ) 2 n 2 (1 2 ) ( s n ) 2
63
Second order system responses
Overdamped response:
Poles: Two real at
1, 2
Natural response: Two exponentials with time constants equal to the reciprocal of the pole location
Underdamped response:
d j d
Natural response: Damped sinusoid with an exponential envelope whose time constant is equal to the
reciprocal of the pole‘s real part. The radian frequency of the sinusoid, the damped frequency of oscillation,
is equal to the imaginary part of the poles
c(t ) Ae d t cos(d t )
Undamped response:
j1
64
Natural response: Undamped sinusoid with radian frequency equal to the imaginary part of the poles
c(t ) A cos(1t )
Critically damped responses:
Natural response: One term is an exponential whose time constant is equal to the reciprocal of the pole
location. Another term is the product of time and an exponential with time constant equal to the reciprocal of
the pole location
It is observed that, for an underdamped system, there are two complex conjugate poles.
Usually, even if a system is of higher order, the two complex conjugate poles nearest to the
j - axis (called dominant poles) are considered and the system is approximated by a second
order system. Thus, in designing any system, certain design specifications are given based on
the typical underdamped step response shown as Fig.
Delay time t d : It is the time required for the response to reach 50% of the steady state value
for the first time.
Rise time t r : It is the time required for the response to mach 100% of the steady state value
for under damped systems. However, for over damped systems, it is taken as the time
required for the response to rise from 10% to 90% of the steady state value.
Peak time t p : It is the time required for the response to reach the maximum or Peak value of
the response.
Peak overshoot M p : It is defined as the difference between the peak value of the response
and the steady state value. Iris usually expressed in percent of the steady state value. If the
time for the peak is t p , percent peak overshoot is given by,
66
Steady State Error
E(s) 1 R(s)
or E(s)
R(s) 1 G(s) 1 G(s)
sR(s)
ess lime(t) limsE(s) lim
t s 0 s 0 1 G(s)
Proportional controllers
– pure gain or attenuation
Integral controllers
67
– integrate error
Derivative controllers
– differentiate error
Proportional Controller
Integral Controller
Integral of error with a constant gain
Increase system type by 1
Infinity steady-state gain
Eliminate steady-state error for a unit step input
Integral Controller
Derivative Control
68
PI controller, PD controller
PID controller
Controller Performance
P controller
PI controller
PD controller
PID controller
Design of PID Controllers
Based on the knowledge of P, I and D
– trial and error
– manual tuning
– simulation
Design of PID Controllers
Time response measurements are particularly simple.
A step input to a system is simply a suddenly applied input - often just a constant
voltage applied through a switch.
The system output is usually a voltage, or a voltage output from a transducer
measuring the output.
A voltage output can usually be captured in a file using a C program or a Visual Basic
program.
You can use responses in the time domain to help you determine the transfer function of
a system.
First we will examine a simple situation. Here is the step response of a system. This is
an example of really "clean" data, better than you might have from measurements. The input
to the system is a step of height 0.4. The goal is to determine the transfer function of the
system.
Impulse Response of A First Order System
The impulse response of a system is an important response. The impulse response is
the response to a unit impulse. The unit impulse has a Laplace transform of unity (1).That
gives the unit impulse a unique stature. If a system has a unit impulse input, the output
transform is G(s), where G(s) is the transfer function of the system. The unit impulse
response is therefore the inverse transform of G(s), i.e. g(t), the time function you get by
inverse transforming G(s). If you haven't begun to study Laplace transforms yet, you can just
file these last statements away until you begin to learn about Laplace transforms. Still there
is an important fact buried in all of this.
69
Knowing that the impulse response is the inverse transform of the transfer function of
a system can be useful in identifying systems (getting system parameters from measured
responses).
In this section we will examine the shapes/forms of several impulse responses.
We will start with simple first order systems, and give you links to modules that discuss
other, higher order responses.
A general first order system satisfies a differential equation with this general
form.
If the input, u(t), is a unit impulse, then for a short instant around t = 0 the input is
infinite.
Let us assume that the state, x(t), is initially zero, i.e. x(0) = 0. We will integrate both
sides of the differential equation from a small time, , before t = 0, to a small time, after t =
0. We are just taking advantage of one of the properties of the unit impulse.
The right hand side of the equation is just Gdc since the impulse is assumed to be a unit
impulse - one with unit area. Thus, we have:
We can also note that x(0) = 0, so the second integral on the right hand side is zero. In
other words, what the impulse does is it produces a calculable change in the state, x(t), and
this change occurs in a negligibly short time (the duration of the impulse) after t = 0 That
leads us to a simple strategy for getting the impulse response. Calculate the new initial
condition after the impulse passes. Solve the differential equation - with zero input - starting
from the newly calculated initial condition.
************************
70
UNITII
IIII
TIME RESPONSE ANALYSIS
PART-A
1. What is Proportional controller and what are its advantages?
2. What is the drawback in P-controller?
3. What is integral control action?
4. What is the advantage and disadvantage in integral controller?
5. What is PI controller?
6. What is PD controller?
7. What is PID controller?
8. What is time response?
9. What is transient and steady state response?
10.Name the test signals used in control system.
11. Define Step signal:
12. Define Ramp signal:
13. Define parabolic signal:
14. What is an impulse signal?
15. What is the order of a system?
16. Define Damping ratio.
17. Give the expression for damping ratio of mechanical and electrical system.
18. How the system is classified depending on the value of damping?
29. What will be the nature of response of a second order system with different types of
damping?
20. Sketch the response of a second order under damped system.
21. What is damped frequency of oscillation?
22. List the time domain specifications:
71
24. Define Delay time.
25. Define rise time.
26. Define Peak time.
27. Define Peak overshoot.
28. Define settling time.
29. What is type number of a system? What is its significance?
30. Distinguish between type and order of a system:
31. What is steady state error?
32. Define acceleration error constant:
33. What are generalized error coefficients?
34. Give the relation between generalized and static error coefficients:
35. Mention two advantages of generalized error constants over static error constants
PART-B
1. (a) Derive the expressions and draw the response of first order system for unit step
input. (8)
(b) Draw the response of second order system for critically damp case and when input is
unit step. (8)
2. Derive the expressions for Rise time, Peak time, and Peak overshoot. (16)
3. A potential control system with velocity feedback is shown in fig. What is response of the
system for unit step input?
72
7. A closed loop servo is represented by the differential equation, where c is the displacement
of the output shaft, r is the displacement of the input shaft and e= r-c. Determine undamped
natural frequency, damping ratio and percentage maximum overshoot for unit step input. (16)
8. 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.
(b) The steady state error when the input is R(S) where R(S) =3/S –2/S2 +1/3S3 (16)
9. The open loop transfer function of a servo system with unity feedback system is G(S) = 10/
S(0.1S+1).Evaluate the static error constants of the system. Obtain the steady state error of
the system when subjected to an input given by the polynomial r(t) = a0+a1t +a2 /2 t2 . (16)
UNIT III
FREQUENCY RESPONSE
Frequency response
Introduction
If a time signal like step and ramp are used to excite the system and its time response
is obtained, we call it a time response analysis or time domain analysis.
On the other hand, if a sinusoidal signal of variable frequency is used to excite the
system and the magnitude and phase of the steady state output from the system is
measured, we call it frequency response analysis or frequency domain analysis.
73
The time domain response of second order system to a unit step unit. In this section
the frequency response of a second order system is obtained and a correlation between
time domain at frequency domain response.
We have
The time response for a unit sinusoidal input with frequency is given by
The magnitude and phase of steadystate sinusoidal response for variable frequency can be
plotted. It is to be noted that when
u0 M 1 0
1
u 1 M 90
2
u M 0 and 180
74
Monotonically decrease s from 1 to 0 or it attains a maximum value and then decreases to
zero. If it attains the maximum value at some frequency ,its derivative should be zero at that
frequency.
This frequency where the magnitudes become maximum is known as the resonance
frequency. Substituting u u r ,we get the maximum value of the magnitude response. This
value of M M r is known as the resonance peak. Hence
75
Fig.(a) Frequency response :Magnitude Vs Normalised frequency
One important observation can be made about the resonant peak. It is clear that the resonant
peak depends only on the damping factor . For a given , the resonance frequency is
indicative of the natural frequency of the system. In other words, resonant frequency is a
4
measure of the speed of response of the system since settling time t s
n
Just as peak overshoot and settling time are used as performance measures of a control
system in time domain, the resonant peak and resonance frequency can be used as
performance measures for a control system in frequency domain.
As M=1 at u=0 control system are considered as lowpass filters and the frequency at which
the magnitude falls to -3db is known as the bandwidth of the system. This shown in the fig .
frequencies beyond u=ub are greatly attenuated.
76
Squaring on both sides of above equation and simplifying, we have
Only positive sign is considered in above equation because ub must be positive and real
The band width is indicative of the noise characteristics of the system. If the bandwidth is
more the system is more susceptible to noise signals. Also for a given , the bandwidth b is
a measured n and hence the speed of response. If the bandwidth is more, the speed of
response is high.
Bode plot
• Plots of the magnitude and phase characteristics are used to fully describe the
frequency response
• A Bode plot is a (semilog) plot of the transfer function magnitude and phase angle as
a function of frequency
dB = 20 log10 A
Gain factor, K
(jω)±p factor: pole and zero at origin
(1+jωT)±q factor
77
Quadratic factor
r
2
1 j 2
n n 2
Gain margin:
The gain margin is the number of dB that is below 0 dB at the phase crossover
frequency (ø=-180º). It can also be increased before the closed-loop system becomes unstable
Phase margin:
The phase margin is the number of degrees the phase of that is above -180º at the gain
crossover frequency
For the following T.F draw the Bode plot and obtain Gain cross over frequency (w gc) ,
Phase cross over frequency , Gain Margin and Phase Margin.
Solution:
78
Term Corner Frequency Slope Change
db/dB in slope
= 38.06 dB
A @ wc2 ; A = [Slope from wc1 to wc2 x log (wc2 / wc1 ] + Gain (A)@ wc1
= 33 dB
= - 27 dB
79
Calculation of Phase angle for different values of frequencies [PHASE PLOT]
When
80
GM = - { 20 log [G( jwpc )] = - { 32 } = -32 dB
Phase Margin
Ґ = 180o+ Øgc= 180o + (- 240o) = -60o
Conclusion
For this system GM and PM are negative in Values. Therefore the system is unstable
in nature.
Polar plot
To sketch the polar plot of G(jω) for the entire range of frequency ω, i.e., from 0 to
infinity, there are four key points that usually need to be known:
(1) the start of plot where ω = 0,
(3) where the plot crosses the real axis, i.e., Im(G(jω)) = 0, and
(4) where the plot crosses the imaginary axis, i.e., Re(G(jω)) = 0.
81
Write the frequency corresponding to each of the point of the plot.
Minimum Phase Systems
Systems with all poles & zeros in the Left half of the s-plane – Minimum Phase
Systems.
For Minimum Phase Systems with only poles
Type No. determines at what quadrant the polar plot starts.
Order determines at what quadrant the polar plot ends.
Type No. → No. of poles lying at the origin
Order → Max power of ‗s‘ in the denominator polynomial of the transfer function.
Gain Margin
Gain Margin is defined as ―the factor by which the system gain can be increased to drive the
system to the verge of instability‖.
For stable systems,
ωgc ωpc
G(j)H(j) at ω=ωpc 1
GM = in positive dB
More positive the GM, more stable is the system.
For marginally stable systems,
ωgc = ωpc
G(j)H(j) at ω=ωpc = 1
GM = 0 dB
For Unstable systems,
ωgc ωpc
G(j)H(j) at ω=ωpc > 1
GM = in negative dB
Gain is to be reduced to make the system stable
Note:
If the gain is high, the GM is low and the system‘s step response shows high
overshoots and long settling time.
On the contrary, very low gains give high GM and PM, but also causes higher ess,
higher values of rise time and settling time and in general give sluggish response.
Thus we should keep the gain as high as possible to reduce ess and obtain acceptable
response speed and yet maintain adequate GM and PM.
An adequate GM of 2 ( ie 6 dB ) and a PM of about 30 is generally considered good
enough as a thumb rule.
82
At = pc , G(j)H(j) = -180
Let G(j)H(j) at = pc be taken as B
If the gain of the system is increased by a factor 1/B, then the G(j)H(j) at =
pc becomes B(1/B) = 1 and hence the G(j)H(j) locus pass through -1+j0 point
driving the system to the verge of instability.
GM is defined as the reciprocal of the magnitude of the OLTF evaluated at the phase
cross over frequency.
GM = Kg = 1 / G(j)H(j) = pc
GM in dB = 20 log (1/B) = - 20 log B
Phase Margin
Phase Margin is defined as ― the additional phase lag that can be introduced before the system
becomes unstable‖.
Let ‗A‘ be the point of intersection of G(j)H(j) plot and a unit circle centered
at the origin.
Draw a line connecting the points ‗O‘ & ‗A‘ and measure the phase angle
between the line OA and +ve real axis.
This angle is the phase angle of the system at the gain cross over frequency.
G(jgc)H(jgc) = gc
If an additional phase lag of PM is introduced at this frequency, then the phase
angle G(jgc)H(jgc) will become 180 and the point ‗A‘ coincides with (-
1+j0) driving the system to the verge of instability.
This additional phase lag is known as the Phase Margin.
= 180 + G(jgc)H(jgc)
= 180 + gc
83
A polar plot may be constructed from experimental data or from a system transfer
function
If values of w are marked along the contour, a polar plot has the same information
as a Bode plot
Usually, the shape of a polar plot is of most interest
In the following a comparison of time domain and frequency domain analysis is given.
Variable frequency, sinusoidal signal generators are readily available and precision
measuring instruments are available for measurement of magnitude and phase angle.
The time response for a step input is more difficult to measure with accuracy.
If the system has large time constants, it makes more time to reach steady state at each
frequency of the sinusoidal input. Hence time domain method is preferred over
frequency domain method in such systems.
In order to do a frequency response test on a system, the system has to be isolated and
the sinusoidal signal has to be applied to the system. This may not be possible in
systems which cannot be interrupted. In such cases, a step signal or an impulse signal
may be given to the system to find its transfer function. Hence for systems which
cannot be interrupted, time domain method is more suitable.
The design of a controller is easily done in the frequency domain method than in time
domain method. For a given set of performance measures in frequency domain.
The effect of noise signals can be assessed easily in frequency domain rather than
time domain.
The most important advantage of frequency domain analysis is the ability to obtain
the relative stability of feedback control systems. The Routh Hurwitz criterion is
essentially a time domain method which determines the absolute stability of a system.
The determination of relative stability by Routh Hurwitz criterion is cumbersome.
Nyquist criterion, will not only give stability but also relative stability of the system
without actually finding the roots of the characteristic equation.
Since the time response and frequency response of a system are related through Fourier
transform, the time response can be easily obtained from the frequency response. The
correlation between time and frequency response can be easily established so that the time
domain performance measures can be obtained from the frequency domain specifications and
vice versa.
The time domain specifications are obtained by subjecting the second order system to a unit
step input.
84
Peak overshoot
M p e / 1 2
For 0 1
d n 1 2
Settling time
4
ts
n
In the frequency domain, the second order system is subjected to constant amplitude, variable
frequency; sinusoidal input and the magnitude and phase response are obtained. The
important frequency domain specifications are
Resonant peak
1
Mr For <0.707
2 1 2
Resonance frequency
r n 1 2 2
Band width
1
b n [1 2 2 4 4 ]
2 2 4 2
Both M p and M r are dependent on the damping factor only and hence they are both
inductive of damping in the system.
Given M p , the resonant peak M r can be evaluated provided is less than 0.707. This
condition is usually satisfied by many practical control system as is seldom greater than
0.707. Thus the resonant peak M r and peak overshoot are well correlated. A plot of M p and
M r with respect to damping factor is given in the fig.
85
The correlation between the time domain specifications of damped natural frequency domain
specification of resonant frequency. For a given damping factor , the ratio r is fixed and
d
given the frequency domain specification, the corresponding time domain specification and
vice versa, can be easily obtained.
The speed of response is indicated by settling time in time domain. The bandwidth a
frequency domain concept is also indicative of speed of response. Thus we can see that there
is a perfect correlation between time domain and frequency domain performance measures
given one, the other can be obtained easily. Of course one should keep in mind that the
correlation is valid only for < 0.707, which is usually satisfied in many practical control
systems.
*************************
UNIT III
IVIV IV
FREQUENCY RESPONSE
PART -A
1. What is frequency response?
2. What are advantages of frequency response analysis?
3. What are frequency domain specifications?
86
4. Define Resonant Peak.
5. What is resonant frequency?
6. Define Bandwidth.
7. What is cut-off rate?
8. Define gain margin.
9. Define phase margin.
10. What is phase and Gain cross-over frequency?
11. What is Bode plot?
12. Define corner frequency.
13. What are the advantages of Bode Plot?
14. What is a Nichols plot?
15. What are M and N circles?
16. What is Nichols chart?
17. What are the advantages of Nichols chart?
18. What is polar plot?
19. What is minimum phase system?
20. What are All-Pass systems?
PART-B
1. Plot 1. Plot the Bode diagram for the following transfer function and obtain the gain and
phase cross over frequencies: G(S) = 10/ S(1+0.4S) (1+0.1S) (16)
2. The open loop transfer function of a unity feedback system is G(S) = 1/ S (1+S) (1+2S)
Sketch the Polar plot and determine the Gain margin and Phase margin (16)
3. Sketch the Bode plot and hence find Gain cross over frequency, Phase cross over
frequency, Gain margin and Phase margin. G(S) = 0.75(1+0.2S)/ S(1+0.5S) (1+0.1S)(16)
4. Sketch the Bode plot and hence find Gain cross over frequency, Phase cross over
frequency, Gain margin and Phase margin. G(S) = 10(S+3) / S(S+2) (S2+4S+100) (16)
5. Sketch the polar plot for the following transfer function and find Gain cross over
frequencies, Phase cross over frequency, Gain margin and Phase margi G(S) = 10(S+2)
(S+4)/ S (S2 -3S+10) (16)
6. Construct the polar plot for the function GH(S) =2(S+1)/ S2. Find Gain cross over
frequency, Phase cross over frequency, Gain margin and Phase margin. (16)
7. Plot the Bode diagram for the following transfer function and obtain the gain and phase
cross over frequencies. G(S) =KS2 / (1+0.2S) (1+0.02S). Determine the value of K for a gain
cross over frequency of 20 rad/sec. (16)
87
8. Sketch the polar plot for the following transfer function and find Gain cross over
frequency, Phase cross over frequency, Gain margin and Phase margin.G(S) = 400/ S (S+2)
(S+10) (16)
9. A unity feedback system has open loop transfer function G(S) = 20/ S(S+2) (S+5).Using
Nichol‘s chart determine the closed loop frequency Response and estimate all thefrequency
domain specifications. (16)
10. Sketch the Bode plot and hence find Gain cross over frequency, Phase cross over
frequency, Gain margin and Phase margin: G(S) = 10(1+0.1S)/ S (1+0.01S) (1+S). (16)
11. Write short notes on correlation between the time and frequency response? (16)
12. What is compensation? Why it is needed for control system? Explain the types of
compensation (16)
13. Realize the basic compensators using electrical network and obtain the transfer function.
(16)
14. Design a suitable lead compensators for a system with unity feedback and having open
loop transfer function G(S)= K/ S(S+1) (S+4) to meet the specifications.
(i)Damping ratio=0.5
(ii) Undamped natural frequency _n =2 rad/sec. (16)
15. A unity feedback system has an open loop transfer function G(S) = K/ S(S+1)
(0.2S+1).Design a suitable phase lag compensators to achieve the following specifications
Kv= 8 and Phase margin 40 deg with usual notation. (16)
16. Explain the procedure for lead compensation and lag compensation. (16)
17. Explain the design procedure for lag-lead compensation. (16)
18. Consider a type 1 unity feedback system with an OLTF G(S) =K/S (S+1) (S+4).
The system is to be compensated to meet the following specifications Kv > 5sec
and PM>43 deg. Design suitable lag compensators.
UNIT IV
STABILITY ANALYSIS
88
Stability
A system is stable if any bounded input produces a bounded output for all bounded
initial conditions.
Characteristic Equation
Consider an nth-order system whose the characteristic equation (which is also the denominator
of the transfer function) is
a( s ) s n a1 s n 1 a 2 s n 2 a n 1 s 1 a 0 s 0
89
Real negative i.e in Absolutely
1.
L.H.S of s-plane stable
Complex conjugate
with negative real Absolutely
2.
part i.e in L.H.S of stable
s-plane
Complex conjugate
with positive real
4. unstable
part i.e in R.H.S of
S-plane
Repeated pair on
imaginary axis
6. without any pole Unstable
in R.H.S of s- plane
A necessary condition for stability of the system is that all of the roots of its
characteristic equation have negative real parts, which in turn requires that all the
coefficients be positive.
90
A necessary (but not sufficient) condition for stability is that all the coefficients of the
polynomial characteristic equation are positive & none of the co-efficient vanishes.
A system is stable if and only if all the elements of the first column of the Routh array
are positive
a ( s ) 1 s n a1s n 1 a2 s n 2 a3 s n 3 an 1s1 a0 s 0
Routh array method
sn : 1 a2 a4
s n 1 : a1 a3 a5
n2
s : b1 b2 b3
s n 3 : c1 c2 c3
s2 : * *
s1 : *
s0 : *
1 a3 a2 a1
b1 ,
a1
1 a5 a4 a1
b2 ,
a1
1 a7 a6 a1
b3
a1
a( s) s 6 4s 5 3s 4 2s 3 s 2 4s 4
Is the system described by this characteristic equation stable?
Answer:
91
s6 : 1 3 1 4
5
s : 4 2 4 0
4
s : 52 0 4
3
s : 2 12 5 0
2
s : ? ?
1
s : ? ?
0
s : ?
s6 : 1 3 1 4
5
s : 4 2 4 0
s4 : 52 0 4
3
s : 2 12 5 0
s2 : 3 4
s : 76 15
1
0
s0 : 4
The elements of the 1st column are not all positive. Then the system is unstable.
Form Auxiliary equation by using the co-efficient of the row which is just above the
row of zeros
Zero and the same remaining row contains atleast one non-zero element
92
Substitute a small positive no. ‗ε‘ in place of zero and complete the array.
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:
o damping ratio,
o natural frequency, wn
o damping constant, first-order modes
o Consider how these roots change as the loop gain is varied from 0 to
93
Find BA and BI points
AOD = 180º- (sum of angles of vectors to the complex pole from all other poles) +
AOA = 180º- (sum of angles of vectors to the complex zero from all other zeros) +
(sum of angles of vectors to the complex zero from poles)
1 F ( s) 0
Step 2: Rewrite preceding equation into the form of poles and zeros as follows
m
(s z
j 1
j )
1 K n
0
( s pi )
Step 3: i 1
Locate the poles and zeros with specific symbols, the root locus begins at the open-
loop poles and ends at the open-loop zeros as K increases from 0 to infinity
If open-loop system has n-m zeros at infinity, there will be n-m branches of the root
locus approaching the n-m zeros at infinity
Step 4:
The root locus on the real axis lies in a section of the real axis to the left of an odd
number of real poles and zeros
Step 5:
Step 6:
The root loci must be continuous and symmetrical with respect to the horizontal real
axis
Step 7:
The loci proceed to zeros at infinity along asymptotes centered at centroid and with
angles
n m
pi z j
i 1 j 1
a
n94 m
(2k 1)
a (k 0,1, 2, n m 1)
nm
Step 8:
The actual point at which the root locus crosses the imaginary axis is readily evaluated
by using Routh‘s criterion
m n
1 1
d z
j 1
i 1 d pi
j
Step 9:
Step 10:
Determine the parameter value K1 at a specific root using the magnitude criterion
n
(s p ) i
K1 i 1
m
1. Two necessary but not sufficient conditions that all the roots have negative real parts are
a) All the polynomial coefficients must have the same sign.
b) All the polynomial coefficients must be nonzero.
2. If condition (1) is satisfied, then compute the Routh-Hurwitz array as follows
95
where the ai's are the polynomial coefficients, and the coefficients in the rest of the
table are computed using the following pattern
96
3. The necessary condition that all roots have negative real parts is that all the elements of the
first column of the array have the same sign. The number of changes of sign equals the
number of roots with positive real parts.
4. Special Case 1: The first element of a row is zero, but some other elements in that row are
nonzero. In this case, simply replace the zero elements by "", complete the table
development, and then interpret the results assuming that "" is a small number of the same
sign as the element above it. The results must be interpreted in the limit as 0.
5. Special Case 2: All the elements of a particular row are zero. In this case, some of the roots
of the polynomial are located symmetrically about the origin of the s-plane, e.g., a pair of
purely imaginary roots. The zero rows will always occur in a row associated with an odd
power of s. The row just above the zero rows holds the coefficients of the auxiliary
polynomial. The roots of the auxiliary polynomial are the symmetrically placed
roots. Be careful to remember that the coefficients in the array skip powers of s from one
coefficient to the next.
Let P = no. of poles of q(s)-plane lying on Right Half of s-plane and encircled by s-plane
contour.
Let Z = no. of zeros of q(s)-plane lying on Right Half of s-plane and encircled by s-
plane contour.
For the CL system to be stable, the no. of zeros of q(s) which are the CL poles that lie
in the right half of s-plane should be zero. That is Z = 0, which gives N = -P.
Therefore, for a stable system the no. of ACW encirclements of the origin in the q(s)-
plane by the contour Cq must be equal to P.
St
ab
le
Un 97
sta
Problem
Sketch the Nyquist stability plot for a feedback system with the following open-loop transfer
function
1
G( s) H ( s)
s ( s s 1)
2
Solution
1
G ( j ) H ( j )
j (1 2 j )
On section bcd, s = RejR ; therefore i.e. section bcd maps onto the origin of the G(s)H(s)-
plane
1
G( s) H ( s) 0
R3
98
Relative stability
The main disadvantage of a Bode plot is that we have to draw and consider two different
curves at a time, namely, magnitude plot and phase plot. Information contained in these two
plots can be combined into one named polar plot.
The polar plot is for a frequency range of 0 . while the Nyquist plot is in the
frequency range of The information on the negative frequency is redundant
because the magnitude and real part of G ( jw) an are even functions.
. In this section. We consider how to evaluate the system performance in terms of relative
stability using a Nyquist plot.
The open-loop system represented by this plot will become unstable beyond a certain value.
As shown in the Nyquist plot of Fig. the intercept of magnitude 'a on the negative real axis
corresponds lost phase shift of - 180° and - 1 represents the amount of increase in gain that
can be tolerated before closed-loop system tends toward instability. As 'a' approaches
(1 j 0) point the relative stability is reduced, The gain and phase margins are represented
as follows in the Nyquist plot.
Gain margin
As system gain is increased by a factor 1/a, the open loop G( jw) H ( jw) will increase by
1
a factor a ( ) 1and the system would be driven to instability.Thus, the gain margin is the
a
reciprocal of the gain at the frequency at which the phase angle of the Nyquist plot is - 1800.
The gain rnargin, usually measured in dB, is a positive quantity given by
GM = 20 log adB
Phase Margin m
Importance of the phase margin has already in the content of Bode. Phase margin is defined as
99
the change in open-loop phase shift required al unity gain to make a closed loop system
unstable. A closed-loop system will be unstable if the Nyquist plot encircles 1 j 0 point.
Therefore, the angle required to make this system marginally stable in a closed loop is the
phase margin .In order to measure this angle, we draw a circle with a radius of 1, and find the
point of intersection of the Nyquist plot with this circle, and measure the phase shift needed for
this point to be at an angle of 1800. If may be appreciated that the system having plot of Fig
with larger PM is more stable than the one with plot of Fig.
100
UNITIV
IVIV IV
STABILITY ANALYSIS
PART-A
1. Define BIBO Stability.
2. What is impulse response?
3. What is characteristic equation?
4. How the roots of characteristic equation are related to stability?
5. What is the necessary condition for stability?
6. What is the relation between stability and coefficient of characteristic polynomial?
7. What will be the nature of impulse response when the roots of characteristic equation
8. What will be the nature of impulse response if the roots of characteristic equation are
9. What is principle of argument?
10. What is the necessary and sufficient condition for stability?
11. What is routh stability condition?
12. What is auxiliary polynomial?
13. What is quadrantal symmetry?
14. In routh array what conclusion can you make when there is a row of all zeros?
15. What is limitedly stable system?
16. What is Nyquist stability criterion?
17. What is root locus?
18. How will you find root locus on real axis?
19. What are asymptotes?
20. What is centroid, how it is calculated?
21. What is breakaway point and break in point?
22. What is dominant pole?
PART-B
1. Using Routh criterion determine the stability of the system whose characteristics
equation is S4+8S3+18S2+16S+5 =0 . (16)
2. F(S) =S6 +S5-2S4-3S3-7S2-4S1-4 =0.Find the number of roots falling in the RHS plane
and LHS plane. (16)
101
3. Draw the Nyquist plot for the system whose open loop transfer function is
G(S) H(S) =K/S (S+2) (S+10).Determine the range of K for which closed loop system is
stable. (16)
4. Construct Nyquist plot for a feedback control system whose open loop transfer
function is given by G(S)H(S) =5/ S(1-S).comment on the stability of open loop and
closed loop transfer function. (16)
5. Sketch the Nyquist plot for a system with the open loop transfer function
G(S)H(S) =K (1+0.5S) (1+S) / (1+10S) (S-1).determine the range of values of K for
which the system is stable. (16)
6. Sketch the root locus for the open loop transfer function of unity feedback control system
given below: G(S) H(S)=K/S(S+2)(S+4). (16)
7. Sketch the root locus for the open loop transfer function of unity feedback control system
given below: G(S) H(S) =K/S(S+1) (S+2).Also find K of breakaway point. (16)
UNIT V
COMPENSATOR DESIGN
102
Performance criteria – Lag, lead and lag-lead networks – Compensator design using
bode Plots.
The three types of compensators, viz, lead, lag and lag-lead compensators can be
realized by electrical, mechanical.
Pneumatic, hydraulic and other components. We will discuss the realization of these
compensators by electrical RC networks.
We will also obtain their frequency response characteristics which are useful in their
design.
Lead Compensator
We have seen in section that the pole zero configuration of a lead compensator is as
given in fig
103
The transfer function of the network is given by
Since 1 it is a lead compensator. Given the time constant and the d.c gain of the
network the three components of RC network namely R1, R2 and C have to determined.
This means that the network produces an attenuation of , and this d.c attenuation has to be
1
cancelled by using an amplifier of gain A .thus the compensating network takes the
form as shown in fig.
104
Under the steady conditions the output of this network leads the input and hence this network
is known as a lead network.
G( j ) is given by
The angle is a function of and it will be a maximum when tan is a maximum with
respect to
105
1 1
m is the geometric mean of the two corner frequencies and
We have
1
Substituting m we have
106
The lead compensator is required to provide the necessary angle at the gain crossover
frequency, to obtain the required phase margin. Knowledge the value of required angle m to
be provided by the network parameter can be obtained .the value of required to get a
maximum phase lead m is tabulated
m
0 1
20 0.49
40 0.217
60 0.072
80 0.0076
From the table it can be seen that for obtaining larger phase angles, has to be reduced to
very low value. For very low values of . the pole of the lead network is located for away
from the j ,axis and from the Bode magnitude diagram it can be seen that the network has a
1
large gain of 20 log db at high frequencies. Than the high frequency noise signals are
amplified and signal to noise ratio becomes very poor. Hence it snot desirable to have low
values . A value of = 0.1 is considered to be suitable. When large phase angles are
desired, two lead networks in cascade may be used.
Lag Compensator
1 ,Ensures that the pole is to the right of the zero. The pole zero configurations is given
in fig
The d.c gain of this network is unity and hence no d.c amplifier is required as in the case of a
lead network.
108
Under steady state conditions the output of the network lags the input and therefore this
network is known as a lag network. The maximum phase lag m is again obtained from the
equation and frequency at which it occurs, replacing by .At high frequencies the network
produces an attenuation of 20 log and hence signal to noise ratio is improved by this
network.
Lag-Lead Compensator
109
Comparing the equation
And
1
110
Cascade Compensation in Frequency Domain Using Bode Plots
In order to design a compensator using Bode plot, the specifications are to be given in
frequency domain.
If specifications are given in time domain, they can be easily converted into frequency
domain using time and frequency domain correlations.
The design is carried out in the frequency domain using Bode plots; the time response
of the compensated system must be obtained to check the time response
specifications.
This is necessary because the time and frequency domain specifications arc correlated
under the assumption that the compensated system is of second order, or it has a pair
of dominant closed loop poles.
Based on the time response characteristics. The design may be altered suitably.
The frequency domain method of design is easy to apply but the time response
specifications cannot be directly controlled in this method.
The procedure for the design of a lead compensator is developed in the following.
1. The system gain K is adjusted to satisfy the steady stale error criterion. as specified by
the appropriate error constants K p, K v or K a or the steady state error ess.
111
2. The Bode plot of the system with the desired K is plotted. The phase margin of the
uncompensated system is read from the graph ( 1 ). If this is satisfactory no
compensator is necessary. If the phase margin falls short of the desired phase margin,
additional phase lead m must be provided by a lead compensator at the gain cross
over frequency.
The frequency m ,at which maximum phase lead occurs in a phase lead network is
made to coincide with the gain cross over frequency gc1 , of the uncompensated
system. The effect of this is twofold.
The gain cross over frequency moves to gc 2 which is higher than gc1 . Hence the
phase margin obtained after the compensation is not equal to 1 m as desired
The new phase margin has to be read at the new gain cross over frequency gc 2 .and is
equal to 2 3 as shown in the fig. Obviously 2 1 and 3 m and hence the
phase margin obtained will fall short of the desired value.
It is clear from that if m of lead network is made to coincide with gc1 ,desired phase
margin is not obtained for the compensated network. Thus it is desirable to make
m coincide with gc 2 so that maximum phase lead occurs at the new gain crossover
frequency. But this new gain crossover frequency is not known before hand.
The above difficulty is overcome as follows. The phase lead network provides a gain
1
of 10 log at m .if m is to be coincide with the new gain cross over frequency
gc 2 , gc 2 must be chosen such that the uncompensated system has a gain of
1
10 log at this frequency
112
The maximum phase angle to be provided by the lead network is obtained as
follows:
At the gain crossover frequency of the uncompensated system. The system has a
phase margin of 1 which is less than the desired phase margin d . Hence an
additional phase lead of d 1 , has to be provided by the lead network. But
thesystem will not be providing a phase lead of 1 after compensation because
the gain crossover is going to shift to a new higher value gc 2 . at which the
phase angle provided by the system will be 2 . as shown in Fig. which is less
than 1 . Also gc 2 and hence 2 is not known to start with. In order to
compensate this unknown shortfall, an additional phase angle ranging from
5to 20
is to be provided. Thus the maximum phase lead m to be provided by
the lead network is taken to be,
113
Knowing , we can locate gc 2 from the magnitude plot.
From these values, the corner frequencies of lead network are obtained as
It is observed that the gain cross over frequency is increased after compensation.
Gain cross over frequency can be taken as a measure of the bandwidth of the
system. Since the bandwidth is increased by the lead compensation, Speed of
response is improved. The actual bandwidth can be obtained by transferring the
data from Bode plot to Nichols chart.
1
Since the gain at higher frequencies is increased by 20 log .the noise
frequencies' are amplified and hence the signal to noise ratio deteorates.
Since one of the three elements R1, R2and C can be chosen arbitrarily, as there are only two
equations, impedance level required or the cost of the compensator maybe used to select this
element.
(ii) Using this value of K draw the Bode magnitude and phase plots of the uncompensated
system. Obtain the phase margin of the uncompensated system. Let this be 1
Where
114
d - Desired phase margin
1 -Available phase margin
-Additional phase angle required to compensate the reduction in
phase angle due to increase in gain crossover frequency after
compensation.
1
(v)Find the gain of the phase lead network at m , which is 10 log , Locate the
1
frequency at which the gain of the uncompensated system is 10 log .Let this
frequency be gc 2 .Make gc 2 = m .Read the phase margin 2 , provided by the
system at this frequency m .If 1 2 is ,increase and recalculate
m , and m until 1 2 .
(vii) Draw the bode plot of the compensated system and check the phase margin.
If it falls short of the desired phase margin increase and redesign.
115
1. The gain constant is set as per the steady state error requirements.
2. Bode plot is drawn. If the phase margin obtained is adequate no other compensation is
required. If phase margin is not adequate, locate a frequency gc 2 , at which the desired phase
margin is available.
3. The gain of the uncompensated system at this frequency is measured. This gain will be
positive for type 1 or type 0 systems.
4. The attenuation characteristic of phase lag network at high frequencies is used to reduce
the gain of the uncompensated system to zero at gc 2 , so that this frequency becomes the gain
crossover frequency of the compensated system, as shown in Fig. The phase lag introduced
by the lag network is detrimental to the system. Hence the frequency at which maximum
phase lag occurs for the phase lag network must be located for away from the gain crossover
5. The corner frequencies of the lag network are located far to the left of the new gain
cross over frequency gc 2 ' so that the phase lag provided by the lag network does not
appreciably affect the phase of the system at this frequency. However a small negative
angle will be added to the system phase at gc 2 , as shown in Fig. To compensate this, the
phase margin to be provided is taken to be slightly higher than the desired phase margin.
Thus the phase margin is taken to be
m d
is taken to be in the range of 50o to 15o. If the corner frequencies of lag network are
closer to the new gain cross over frequency. higher value of is chosen. The upper
116
corner frequency 2 , of the lag network is placed, usually, one octave to 1 decade
below the new gain cross over frequency gc 2 .
2. Draw the Bode plot with this value of K. If the phase margin of this system is
inadequate, design a lag compensator.
4. Locate the frequency gc 2 , at which the phase margin is available on the gain plot
of G(s). Find the magnitude of G(s) at this frequency.
gc 2
.Thus
10
7. With and known, the phase lag network can be designed. The transfer function
of the Lag network is,
8. Draw the Bode plot of the compensated system to check the specifications.
117
A lead compensator increases the gain cross over frequency and hence the bandwidth
is also increased. For higher order systems with large error constants, the phase lead
required may be large, which in turn results in larger bandwidth.
This may be undesirable due to noise considerations. On the other hand a lag
compensator results in a lower bandwidth and generally a sluggish system.
Hence a lag-lead compensator is used. The procedure for the design of lag-lead
compensator is as follows:
I. The required error constant is satisfied by choosing a proper gain. The phase margin
and the bandwidth are obtained from the Nichols chart. If the phase margin falls short and
the bandwidth is smaller than the desired value, a lead compensator is designed.
2. If the bandwidth is larger than the desired value, lag compensator is attempted, if the
phase margin desired is available at any frequency. If this results in a lower bandwidth
than the desired bandwidth, a lag lead compensator is to be designed.
3. To design a lag-lead compensator, we start with the design of a lag compensator. A lag
compensator is designed to partially satisfy the phase margin requirement. It means that
the gain cross over frequency gc 2 is chosen to be higher than that to be used, if full lag
compensation is designed. This ensures that the bandwidth is not reduced excessively by
the lag compensator.
4. The and gc 2 are known and the comer frequencies of the lag compensator can be
1
obtained. Since the frequency at which 20 log is available on the magnitude
plot is the new gain crossover frequency gc 3 . The maximum phase lead available is,
With these parameters, a lead compensator is designed. Here, for the lead network,
The log magnitude vs phase angle curve of the lag-lead compensated network is drawn
on the Nichols chart, from which the bandwidth can be obtained.
118
***********************
UNIT V
COMPENSATOR DESIGN
PART-A
119
10. Using Electrical lag-lead network derive the transfer function
11. Design a lead compensator for a unity feedback system with open loop transfer function
G(S) = K/ S(S+1) (S+5) to satisfy the following specifications (i) Kv > 50 (ii) Phase
Margin
is 20 .
12. Design a lead compensator for G(S) =K / S2 (0.2S+1) to meet the following
Specifications (i) Acceleration ka=10; (ii) P.M=35.
13. Design a Lag compensator for the unity feedback system whose closed loop transfer
function C(s) / R(s) = K / (s (s+4) (s+80) + K) is to meet the following specifications
P.M 33 º. And Kv 30.
14. Consider the unity feed back system whose OLTF is G(s) = K / s (s+3)(s+6). Design a
lag-lead compensator to meet the following specifications.
(i) Velocity error constant, Kv = 80, (ii) P.M _ 35º (16)
15. A unity feedback system has an OLTF G(s) = K / s(s+2)(s+60). Design a Lead-Lag
compensator is to meet the following specifications.
(i) P.M is atleast 40º , (ii) Steady state error for ramp input 0.04 rad.
*********************************END*********************************
120
121
122
123
124