Process Control Engineering
Process Control Engineering
P. Sai Krishna
Dept. of Instrumentation Engineering
M.J. College of Engineering & Technology
Hyderabad (A.P).
©Copyright 2020 I.K. International Pvt. Ltd., New Delhi-110002.
This book may not be duplicated in any way without the express written consent of the publisher,
except in the form of brief excerpts or quotations for the purposes of review. The information
contained herein is for the personal use of the reader and may not be incorporated in any commercial
programs, other books, databases, or any kind of software without written consent of the publisher.
Making copies of this book or any portion for any purpose other than your own is a violation of
copyright laws.
Limits of Liability/disclaimer of Warranty: The author and publisher have used their best efforts in
preparing this book. The author make no representation or warranties with respect to the accuracy or
completeness of the contents of this book, and specifically disclaim any implied warranties of
merchantability or fitness of any particular purpose. There are no warranties which extend beyond the
descriptions contained in this paragraph. No warranty may be created or extended by sales
representatives or written sales materials. The accuracy and completeness of the information provided
herein and the opinions stated herein are not guaranteed or warranted to produce any particulars
results, and the advice and strategies contained herein may not be suitable for every individual.
Neither Dreamtech Press nor author shall be liable for any loss of profit or any other commercial
damages, including but not limited to special, incidental, consequential, or other damages.
Trademarks: All brand names and product names used in this book are trademarks, registered
trademarks, or trade names of their respective holders. Dreamtech Press is not associated with any
product or vendor mentioned in this book.
ISBN: 978-93-89633-42-9
EISBN: 978-93-89976-69-4
Edition: 2020
I dedicate this work at the divine feet of my Master and Mentor
Bhagavan Sri Satya Sai Baba
Preface
P. Sai Krishna
Acknowledgements
Author
Contents
Preface vii
Acknowledgements xi
1. Introduction 1
1.1 Introduction to Process Control 1
1.2 Process Characteristics 2
1.3 Elements of Process Dynamics 9
1.4 Interacting and Non-interacting Systems 30
1.5 Process Degrees of Freedom 37
1.6 Piping and Instrumentation Drawings 38
2. Modes of Controller Operation 58
2.1 Reverse and Direct Action of Controllers 59
2.2 Modes of Control Actions 59
2.3 Discontinuous Controller Modes 60
2.4 Continuous Controller Modes 64
2.5 Composite Controller Modes 71
3. Effect of Control Actions on Processes 95
3.1 Concept of Closed Loop Automatic Control 95
3.2 Error Equation and Error Terminology 97
3.3 Effect of Two Position Controller 100
3.4 Effect of Proportional Controller Mode 103
3.5 Effect of Integral Controller Mode 114
3.6 Effect of Derivative Controller Mode 120
3.7 Effect of Composite Controller Modes 123
4. Tuning of Controllers 152
4.1 Introduction 152
4.2 Process Reaction Curve Method 154
4.3 Quarter Amplitude Criteria [Cohen-coon Corrections] 156
4.4 Ziegler-Nicol Tuning Method 158
xiv 2HAB=?A
speed of steam engines was invented by James Watt in the year 1770. After a
long gap of nearly ten decades, Maxwell analysed the dynamics of flyball governor.
During the world war I, positioning of heavy masses like guns and ships
became an absolute necessity.
Minorsky performed the classic work on automatic steering of ships and
positioning of guns on the ship boards.
Servo mechanisms were developed by Hazen during the year 1934. During
world war II, many mathematical, analytical and practical methods were
developed in the field of Automatic Control.
qi (Inlet)
Sight tube
Human operator
h
qo
(Control valve) (Outlet)
FIG. 1.1
Consider a tank having inlet with inlet flow rate of liquid (qi) and outlet
with outflow rate (qo). The objective is to maintain the level of the liquid in the
tank at a prescribed height h measured from the bottom of the tank.
A control valve is placed at the outlet, so that the outlet flow rate may be
altered and monitored.
To achieve the control objective or to convert it into a control system a
human operator is placed near the control valve as shown. This human
operator has the desired value of level in his mind, also known as set point. He
constantly monitors the sight tube and measures the level of the liquid in the
tank. If the level is above the set point or desired value he performs an
operation on the control valve by opening it till the level reaches the set point
and fixes the control valve at particular position which makes outlet flow rate
equal to inlet flow rate, thus maintaining the level at desired value.
Hence, any changes in inlet flow rate causing the level to change, is
continuously monitored by the human operator, who is performing certain
functions like measuring, comparing, evaluating and operating on the control
valve to ensure that the control objective is met.
5TH PROOF 3-9-09
COM -6\D-DRIVE\I.K. INTERNATIONAL\PROCESS CONTROLL ENGINEERING
Introduction 3
Thus the control system shown in Fig 1.1 is manual control system. In this
context it is important to define process as applied to Figure 1.1.
The process obviously refers to ensemble of tank, liquid in the tank, the
inlet and outlet. It is an environment where the level of the liquid is measured,
monitored and controlled.
The human operator who is performing the functions of measuring the
level, comparing with the set point, evaluating and opening or closing of
control valve based on results of evaluation is called as Controller.
The control valve is known as Final Control Element.
To change the manual control system into an automatic process control
system, the human operator is replaced with automatic controller.
Let us study and analyse an automatic example of liquid level
measurement and control as shown.
Set point
Inlet (qi)
Automatic
Controller
LT
h Final control element
qo
(outlet)
FIG. 1.2
Controller block
Comparator
Input or
Automatic Final Control Output
Measurand or Process
– Controller Element (Level)
Set point
Measurement
FIG. 1.3
Introduction 5
Process Load
The set of all load variables that cause a change or deviations of controlled
variable from set point is known as process load. Process load excludes the
controlled variable. Any one of the load variable may be made a manipulating
or controlling variable.
Process Equation
In the example shown in Fig. 1.2, a process equation may be written which
describes a process mathematically. The liquid level h is a function of changes
in inlet flow rate and outlet flow rate. A simple process equation may be written
as
level (h) = f(qi, qo)
Measurement
It is the determination of the existence (or) the magnitude of a variable i.e.
controlled variable. A transducer is normally employed for the measurement.
Controller
It is an automatic device having an output that varies to regulate a controlled
variable in a specified manner. A controller performs operations such as
comparison, evaluation and operation upon the next stage to meet the
demands of control objectives.
A controller may be a self-contained analog or digital instrument, or it may
be the equivalent of such an instrument in a shared control room.
FIG. 1.4
5TH PROOF 4-9-09
COM-6\D-DRIVE\I.K. INTERNATIONAL\PROCESS CONTROLL ENGINEERING
Figure 1.4 shows the physical diagram of Flow Control System. The control
objective is to control the fluid flow rate and obtain constant outlet flow rate,
depending on the set point value. The different load variables of flow control
loop are inlet flow rate (qi) and outlet flow rate (qo).
Since there are only two variables in which outlet flow rate is monitored
and controlled, a simple process equation can be written as
qo = f [qi]
Operation
A orifice plate inserted in the pipe line causes a pressure differential. This
differencial pressure is measured by a differential pressure transmitter (D.P.
Cell) which is a pneumatic device. Since its output is pressure proportional to
the fluid flow rate, it has to be converted into proportional electrical signal if
the controller is an electronic device. Hence, a pressure to current converter is
employed, which converts, pressure into current. The output of the controller,
after evaluation of the measured value by comparing it with set point, is further
converted into a pressure signal so that it can be applied to control valve
actuator, which positions the valve stem to regulate flow rate at the outlet of the
pipe line. It is to be noted that the standard range of 4-20 mA, and the
standard pneumatic signal transmission is 3-15 psi.
Control system evaluation and analysis can be done by considering the
block diagram of above flow control system, in addition to knowing about
elements of flow control loop.
Block Diagram
Controller block
Comparator
Set point Process Flow rate
Automatic I/P Control (Pipe line inlet
– Controller Converter Valve & outlet fluid)
P/I Measurement
Converter (Orifice)
FIG. 1.5
Process constitutes the pipe line, fluid, inlet and outlet. Measurement is
done by orifice plates and electrical signal conversion is done by P to I
converter.
TEMPERATURE CONTROL SYSTEM
A simple temperature control system is showing in Fig.1.6. The control
objective is to maintain the temperature of the liquid at a desired value. The
controlled process variable is therefore TL i.e. temperature of liquid in tank.
5TH PROOF 3-9-09
COM -6\D-DRIVE\I.K. INTERNATIONAL\PROCESS CONTROLL ENGINEERING
Introduction 7
TA
Hot steam TL
qs , T s qo
(outlet)
Control (T T)
valve
Automatic
Controller
Set point
FIG. 1.6
Operation
The temperature transmitter (TT) which is a temperature transducer measures
the temperature of the liquid in the tank. The measured value is sent to the
controller, which compares it with the set point. After proper evaluation, it
sends a signal to the control valve which then regulates the steam flow rate
depending upon the requirement. The standard signal transmission should be
maintained in all process control applications. The output of the temperature
transducer and that of the controller should always be in standard signal
transmission range of 4-20 mA and 3-15 psi respectively.
5TH PROOF 4-9-09
COM-6\D-DRIVE\I.K. INTERNATIONAL\PROCESS CONTROLL ENGINEERING
Block Diagram
Controller block
Comparator
Set point Process Output TL
Automatic Control Tank, liquid, inlet
– Controller Valve outlet, steam flow
Temperature
Transmitter
(Transducer)
FIG. 1.7
Set point
Controller
Pressure Current to
Sensor I/P
Pressure Converter
Pressure
Vessel Control valve
Gas inflow (qi)
Gas outflow (qo)
FIG. 1.8
Measurement
(Pressure
Transducer)
FIG. 1.9
Introduction 9
h
Capillary (1) Capillary (2)
Flow rate q fi
of fluid Pipe line
FIG. 1.10
* Pressure drop through pipes and other equipment is the most obvious illustration of
resistance type element.
H(s) = RQ(s)
Q(S) H(S)
R
q = KA 2 g( h1 - h2 ) ...1(a)
3
q = liquid flow rate (m /s)
K = Flow constant
A = Area of restriction (m2)
g = Acceleration due to gravity (m/s2)
h = head of fluid (m)
Squaring equation 1(a)
q2 = K 2 A 2 2g(h1 h2)
Þ q.q = K 2 A 2 2g(h1 h2)
q 2 ( h1 - h2 )
2 2
= ...1(b)
gK A q
h dh
From equation (1) Resistance R = (or) R = [Rate of change of head
q dq
with flow]
5TH PROOF 3-9-09
COM -6\D-DRIVE\I.K. INTERNATIONAL\PROCESS CONTROLL ENGINEERING
Introduction 11
dh 2 ( h1 - h2 ) q
R= = =
dq q gK 2 A2
128 NL 128 VL
h1 h2 = 4
×q = ×q
πrD g π D4
where
h = head (m)
v = kinematic viscocity (m2/s)
D = Inside distance of pipe (m)
q = Inlet flow rate (m)
m = Absolute viscosity
g = Fluid density (kg/m3)
laminar flow is generally not encountered in industrial applications.
However, the resistance element for laminar flow may be defined as
h dh 128 VL
R= = = ...1(c)
q dq g π D4
Capacitance Element
Capacitance is the ability to store. An illustration of
capacitance element is a tank with inlet as shown in
Figure (1.11). h
The flow of the fluid into the tank is the input and q
the level of the liquid in the tank (h) is the output.
FIG. 1.11
The ability of the tank to store liquid is the
capacitance.
The changes in the level of liquid in the tank (dh/dt) is directly proportional
to inlet flow rate
dh
¥q
dt
dh
Þ C =q ...(2)
dt
where C = capacitance
Applying L.T. to equation (2)
CS H (S) = Q(S)
5TH PROOF 4-9-09
COM-6\D-DRIVE\I.K. INTERNATIONAL\PROCESS CONTROLL ENGINEERING
Transfer function
H (S) 1
=
Q ( S) CS
Q(S) 1 H(S)
CS
z dh
dd
dt = z q
c
dt
Level
(h)
(q/c) t
q
h= t ...(3)
c
C R
h
qo
(outlet)
Control
Valve
FIG. 1.12
5TH PROOF 3-9-09
COM -6\D-DRIVE\I.K. INTERNATIONAL\PROCESS CONTROLL ENGINEERING
Introduction 13
Mathematically,
dh
q i qo µ
dt
(net flow rate) µ (changes in level)
(liquid stored in the tank)
dh
Þ qi qo = C
dt
dh
qi = C + qo
dt
qo = outlet flow rate = h/R (from equation(1))
dh h
\ qo = C +
dt R
Pre multiplying both sides by R
dh
Rqi = RC +h ...(4)
dt
Applying L.T. to equation (4)
R.Q i (s) = RCS H(S) + H (S)
R.Q i ( s) = H(S) (RCS + 1)
Defining time constant element = RC = T
R.Q i (s) = (TS + 1) H (s)
H(S) R
\ = ...(5)
Q i (S) TS + 1
Equation (5) represents the transfer function of Time Constant Process R
in the numerator constitutes proportional element.
Qi(S) R H(S)
(TS + 1)
R
H(S) = Q i (s)
TS + 1
F 1 - T I (partial fractions)
H(S) =
R
S (TS + 1) GH S 1 + TS JK
=R
F1
H(S) = R G -
I
H S T (S + 1 / T ) JK
T
\
F1 1 I
H(S) = R G -
H S S + 1 / T JK
Applying Inverse L.T.
h(t) = R(1 e t /T ) ...(6a)
R(1.0)
R(0.98)
R(0.86)
R(0.63)
t
t = 0 t = 1t t = 2t t = 4t
FIG. 1.13
Introduction 15
C(S) K
=
R(S) 1 + TS
Oscillatory Element
Another important element, which may not be generally encountered in
ordinary liquid, gas and thermal system is oscillatory element.
This element is a typical characteristic of higher order systems. It can be
proved that the response of second order or higher order processes or systems
exhibit oscillations about the steady state value of input. By using appropriate
damping methods these oscillations can be damped to steady state value of
input. A typical example of second order system is PMMC galvanometer
(permanent magnet moving coil). Its principle of operation is based on
Principle of suspension galvanometer which states that a coil of fine wire
suspended between the poles of a permanent magnet with current passing
through it deflects in the magnetic field because of the fundamental property
of electromagnetic force or torque. It continues to deflect until the
electromagnetic torque also known as deflecting torque is counter balanced by
mechanical controlling torque, which is produced due to control springs. At
the equilibrium position when both the torques are equal and opposite in
magnitude, they cancel and the pointer attached to the moving coil, moving
over a calibrated scale, should indicate the magnitude of unknown current
passing through it. But the pointer continues to oscillate about the final value
on account of insufficient friction inherent in the moving parts of the system.
Therefore additional friction is generated by using different damping
methods and this torque is known as damping torque. It is on account of this
damping torque that the pointer rests at the final value indicating the
magnitude of current at the equilibrium position.
Therefore, we conclude that a second order system is characterised in
terms of an important factor called as Damping Factor.
To illustrate the behaviour of a oscillatory element consider a mechanical
translational system as shown.
When a force (F) is applied to this mechanical system, the mass, friction or
damper and spring oppose the applied force. From the principle of
translational mechanics
F = F1 + F2 + F3
dx2 dx
Since F1 = M 2
; F2 = f ; F3 = K x
dt dt
d2 x dx
F=M 2
+f +K x ...(6)
dt dt
5TH PROOF 4-9-09
COM-6\D-DRIVE\I.K. INTERNATIONAL\PROCESS CONTROLL ENGINEERING
F3 F2
Force (F)
FIG. 1.14
X ( S) output 1
\ Transfer function = = ...(7)
F ( S) input MS2 + f S + K
We can compare equations (7) and (8) and obtain expressions for damping
factor xand undamped natural frequency W n in terms of mass (M) spring (K)
and damper (f).
Rewriting equation (7) as
X ( S) 1/ M
= 2 ...(9)
F ( S) S + ( f / m) S + K / M
Comparing the characteristic equation of equation (8) and (9)
f K
S2 + 2 x W n S + W n2 = S2 + S+
M M
K
Wn2 =
M
K FG rad IJ
Þ Wn =
M H sec K
5TH PROOF 3-9-09
COM -6\D-DRIVE\I.K. INTERNATIONAL\PROCESS CONTROLL ENGINEERING
Introduction 17
f
2xWn =
M
K f
\ 2x =
M M
Squaring both sides
4 ξ2 K f2
=
M M2
f2
x2 =
4K M
f
Þ x=
2 KM
R LS
V(S) I(S) 1
CS
F LCS 2
+ RCS + 1 I
V(S) = I(S) GH CS JK
I(S) CS
\ Transfer function = 2
V (S) LCS + RCS + 1
Rewriting it as
I( S) (1 / L)S
= 2 ...(10)
V ( S) S + ( R / L)S + (1 / LC)
Comparing characteristic equation of (8) and (10)
S2 + 2 x W n S + W n2 = S2 + (R/L)S + (1/LC)
5TH PROOF 4-9-09
COM-6\D-DRIVE\I.K. INTERNATIONAL\PROCESS CONTROLL ENGINEERING
1
Wn2 =
LC
1 FG rad IJ
Þ Wn =
LC H sec K
R
2xW nS =
L
Squaring both sides
4 ξ2 R2 R2C
= 2 Þ x2 =
LC L 4L
R C
x=
2 L
- b ± b2 - 4 ac
S1, S2 =
2a
- 2 ξ Wn ± 4 ξ2 Wn2 - 4Wn2
\ S1, S2 =
2a
= x Wn ± Wn ξ2 - 1 ...(11)
Introduction 19
t
jwn
Wn(r/s)
value of input
Steady state
–jwn
–Wn
Oscillatory Response
Consider a second order system with transfer function.
C(S) Wn2
= 2
R(S) S + 2 ξ Wn S + Wn2
Let a sudden input be applied to this system
1
Hence R(S) =
S
1 (S + 2 ξ Wn )
C(S) = 2
S S + 2 ξ Wn S + Wn2
Rewriting (S + 2xW n) as
S x Wn
+ x Wn
S 2x Wn
and the denominator
S2 + 2 x W nS + Wn2 = (S + x W n)2 x2 W n2 + W n2
= (S + x W n)2 + Wn2(1 x2)
Introduction 21
(S + a )
L 1 = e a t cos wt
(S + a )2 + W 2
w
L 1 = e a t sin wt
( s + a )2 + w2
1 (S + ξ Wn ) Wd ξ Wn
\ C(S) = 2 2
2 2
´
S (S + ξ Wn ) + Wd (S + ξ Wn ) + Wd Wd
1 (S + ξ Wn ) Wd ξ Wn
C(S) = ´
S (S + ξ Wn )2 + Wd2 (S + ξ Wn )2 + Wd2 Wn 1 - ξ2
1 (S + ξ Wn ) Wd ξ
C(S) = 2 2
´
S (S + ξ Wn ) + Wd (S + ξ Wn )2 + Wd2 1 - ξ2
e - ξ Wn t FG x IJ
C(t) = 1
1- ξ 2 H 1 - x 2 cos Wd t +
A
sin Wd t
K
We know that Sum of sine wave and cosine wave is another sine wave with
different magnitude and phase angle
FG B IJ
A Sin wt + B Cos wt =
H
A2 + B2 Sin wt + Tan -1
A K
where A=x
B = 1 - ξ 2 as shown
e - ξ Wn t F 1 - ξ2 I
\ C(t) = 1 GG
sin W t + Tan -1 JJ ...(12)
H K
d
1 - ξ2 ξ
C(t)
1.0
Frequency of
Damped
Oscillations
} = wd r/s
Temperature indicator
T1 T2
Tank
Pipe line
Heating element
FIG. 1.16
Let x(t) be the temperature indicated by T1 and y(t) be the temperature
indicated by T2.Under steady state conditions
x(t) = y(t)
Suppose, if the temperature of water at T 1 is suddenly changed. This
change would be indicated at T2 only after an elapse of T unit of Time by the
fluid to flow from T1 end to T 2 end. It is known as transportation lag. Since the
temperature indicated by T2 is delayed by T units of time
y(t) = x(t T) ...(13)
5TH PROOF 4-9-09
Introduction 23
Graphically
Curve - 1
t
T
Y ( S) X(S) Y(S)
= e T S Þ e –TS
X ( S)
1
Writing e T S =
eTS
Expanding using Taylor Series
1
e T S =
FT S I + FT S I
2 2 3 3
(1 + TS) + GH 2! JK GH 3! JK
Keeping only the first two terms
1
e T S = ...(16)
1 + TS
Equation (16) is a crude approximation of dead time. An improvement in
equation (16) is known as Pade approximation.
e -TS /2 1 - TS / 2
e TS = TS /2
= ...(17)
e 1 + TS / 2
Equation (17) is known as 1st order Pade approximation.
SECOND ORDER PADE APPROXIMATION
Introduction 25
F Imaginary part I
|F(jw)| = tan1 GH real part JK
\ |e j w T | = (cos wT )2 + ( - sin wT )2 = 1
F - sin wT I
|e - jwt = + tan1 GH cos wT JK
= tan1( tan wT)
= wT rad
Since p radians ® 180°
wT radians ® ?
-wT
\ ´ 180°
π
Þ 57.3 wT degrees
Therefore we conclude that the magnitude of dead time is unity. It effects
only the phase angle by contributing 57.3 wT degrees
Inlet (qi)
C R u
h
qo
(outlet)
FIG. 1.17
R = fluid resistance
C = fluid capacitance or hydraulic capacitance
u = down stream head
dh
qi qo ¥
dt
dh
qi qo = C
dt
dh
qi = C + qo
dt
Since down stream head is given
h-u
qo =
R
dh h-u
\ qi = C +
dt R
Multiplying both sides by R
dh
Rqi = RC +hu
dt
dh
Rqi + u = RC +h
dt
Let RC = Time constant = T
dh
Rqi + u = T +h ...(21)
dt
Applying L.T. to equation (21)
Rqi(S) + U(S) = (TS + 1) H(S)
R 1
Þ H(S) = Qi(S) + U(S) ...(22)
TS + 1 TS + 1
U(S) 1
TS + 1
Qi(S) R + H(S)
+
1 + TS
Introduction 27
dh
qi + u qo = C
dt
C R
h
dh qo
qi + u = C + qo
dt
FIG. 1.18
h
3 qo =
R
dh h
\ qi + u = C +
dt R
Multiplying both sides by R
dh
Rqi + R.u = RC +h ...(23)
dt
Applying L.T. to equation (23)
RQi(S) + RU(S) = (RCS + 1) H(S)
T = time constant = RC
RQi(S) +RU(S) = (TS + 1) H(S)
R R
H(S) = Q(S) + U(S)
TS + 1 TS + 1
U(S) R
TS + 1
Qi(S) R + H(S)
+
TS + 1
Thermal Systems
Let hi be heat input rate and q be the temperature of the liquid at the
outlet.
The mathematical description of the above system can be obtained by
expressing heat input rate hi in terms of two quantities.
hi = h1 + h2
h1 = Amount of heat stored in tank
h2 = Amount of heat responsible for temperature q
5TH PROOF 4-9-09
Inlet
Stirrer
Outlet
hi
(Temperature of
liquid at outlet)
FIG. 1.19
dq
h1 = MS
dt
where M = Mass of liquid (kg)
S = Specific heat J/kg-°C
Defining thermal capacitance C = MS
dθ
\ h1 = C and h2 = q s q.
dt
q = outlet flow rate (kg/min)
q = temperature of outflowing liquid
Defining Thermal Resistance
1
R=
qs
q
\ h2 =
R
dq θ
hi = c
+
dt R
Multiplying both sides by R
dq
Rhi = Rc +q ...(24)
dt
Applying L.T. to equation (23)
5TH PROOF 4-9-09
Introduction 29
Hi(S) R q (S)
TS + 1
R C C = Gas capacitance
Pi Po R = Gas resistance
Gas inflow
Pressure vessel
FIG. 1.20
Pi(S) 1 Po(S)
TS + 1
5TH PROOF 4-9-09
Non-Interacting System
qi
C1 R1
h1
q1
Tank-1
C2 R2
h2
q0
Tank-2
FIG. 1.21
Consider two tanks (1) and(2) connected as shown in Fig. (1.21). This type
of connection is known as non-interacting system. Any changes in level h, of
tank-1, will change the level h2 of tank-2, but the converse, is not possible i.e.
any changes in level h2 does not affect the level in tank-1. This type of system
is referred to as non-interacting system.
Analysis:
Consider tank-1
dh1
qi q1 = C 1
dt
dh1
qi = C1 + q1 ...(26)
dt
h1
q1 = ...(27)
R1
dh1 h
\ qi = C1 + 1
dt R1
5TH PROOF 4-9-09
Introduction 31
dh1
R i q i = R1C1 + h1 ...(28)
dt
Consider tank-2
dh2
q1 qo = C 2
dt
dh2
q1 = C 2 + qo
dt
h2
qo =
R2
dh2 h
\ q1 = C 2 + 2
dt R2
dh2
R 2q 1 = R 2C 2 + h2
dt
Let T2 = R 2 C 2
dh2
R 2 q 1 = T2 + h2
dt
h1
Since q1 = (from equation (27))
R1
R2 dh
h1 = T 2 2 + h2 ...(29)
R1 dt
Applying L.T. to equations (28) and (29)
R 1 Q i (S) = (T1S + 1) H1(S)
R2
H1(S) = (T2S + 1) H2(S)
R1
The overall input = Qi(S)
The overall output = H2(S)
From equation (30), the overall transfer function is the product of transfer
functions of individual tanks.
Interacting System
Inlet
qi
R1 R2
h1 C 1 h2 C 2
q0
Tank-1 q1 Tank-2
FIG. 1.22
Introduction 33
H2(S)
Q1(S) Q0(S)
1
R2
Feed back
1 1
C2S C2S
Þ
1 R2C2S + 1
1+
R2C2 S R2C2S
R2
Þ
R2C2S + 1
Feed back
Qi(S) 1 – R2 H2(S)
– C1S R1–R2C2S + 1
R2C2S + 1
R2
R2
R1 R2 C2S + R1
R2
1+
R1 R2C2 S + R1
5TH PROOF 4-9-09
R2
Þ
R1 R2C2S + R1 + + R2
Qi(S) 1 R2 H2(S)
C1S R1R2C2S + –R1+ R2
R2C2S + 1
R2
R2
Þ 2
R1C1 R2C2 S + R1C1S + R2C1S + R2C2S + 1
Therefore the transfer function of interacting system is
H2 (S) R2
= 2
...(39)
Q i ( S) R1C1 R2C2S + S( R1C1 + R2C2 + R2C1 ) + 1
As seen from equation (39), unlike the non-interacting system, the transfer
function of interacting system is not the product of transfer transfer functions
of individual tanks.
1
Let Qi(S) = (unit step input)
S
5TH PROOF 4-9-09
Introduction 35
1
H2(S) =
S(TS + 1)2
Considering the partial fractions
1 A B C
2
= + +
S(TS + 1) S TS + 1 (TS + 1)2
A=1 B = T C = T
1 T T
\ H2(S) =
S TS + 1 (TS + 1)2
1 1 T
H2(S) = 2
S S +1/T T (S + 1 / T )2
Applying Inverse Laplace Transform
t t/T
h2(t) = 1 e t/T e ...(40)
T
Now consider the transfer function of interacting system
let R1 = R2 = R
C1 = C2 = C and T= RC
H2 (S) R
\ = 2 2
Q i ( S) T S + 3TS + 1
Let the proportional element in the numerator R = 1
H2 (S) 1
= 2 2
Q i ( S) T S + 3TS + 1
The factor of denominator polynomial are
- 3T ± 9T 2 - 4T 2
=
2T 2
- 3T ± 5T 2
=
2T 2
- 3T ± 2.23T
=
2T 2
- 3T ± 2.23T - 3T ± 2.23T
= ,
2T 2 2T 2
5TH PROOF 4-9-09
- 0.777 - 5.23T
= ,
2T 2 2T 2
- 0.385 2.61
= ,
T T
Therefore the factors are
FG S + 0.385IJ FG S + 2.61IJ = 0
H T KH T K
= (TS + 0.385) (TS + 2.61) = 0
= (2.61TS + 1) (0.38TS + 1) = 0
H2 (S) 1
=
Q i ( S) (2.61TS + 1) (0.38TS + 1)
Subjecting the system to step input
1
H2(S) =
S (2.61TS + 1) (0.38TS + 1)
A B C
= + +
S 2.61TS + 1 0.38TS + 1
A=1 B = 1.17 C = 0.17
1 ( -117)
. 0.17
H2(S) = + +
S 2.61TS + 1 0.38TS + 1
Applying Inverse Laplace Transform
h2(t)
Steady state value of input
1.0
Non-interacting
system response
(faster)
Interacting system
response (sluggish)
Introduction 37
h = h c, h p + 2
=12+2
n=1
Either temperature of steam pressure may be considered an independent
variable.
IDENTIFICATION OF LETTERS
FIRST LETTER SUCCEEDING LETTER
Measured Modifier Readout or Output Modifier
Variable Passive function function
A Analysis Alarm
B Burner UC UC UC
C UC Control
D UC Differential
E Voltage Sensor
F Flow rate (fraction) Ratio
G UC Glass
H Hand High
I Current Indicate
J Power Scan
K Time of change Time rate Station Control
L level Light Low
M UC Momentary Middle
Intermediate
N UC UC UC UC
O UC Orifice
P Pressure Point Cnnection
Q Quantity Integrate
R Radiation Record
S Speed, Frequency Safety Switch
5TH PROOF 4-9-09
Introduction 39
T Temperature Transmitter
U Multi Variable Multi function Multi function Multi function
V Vabration Valve
W Weight, Force Well
X Unclassified X-axis Unclassified Unclassified Unclassified
Y Event (or) State Y-axis Relay
Z Position Z-axis Driver Actuator,
FCE
Note: UC = Users Choice.
INSTRUMENT LINE SYMBOLS
1. Instrument Supply 7. EM or Sonic Signal (Guided)
2. Undefined signal 8. EM or Sonic Signal (Unguided)
3. Pneumatic Signal 9. Internal System Link
4. Electrical Signal (or) 10. Mechanical Link
5. Hydraulic Signal 11. Peneumatic Binary Signal
6. Capillary tube 12. Electrical Binary Signal
(or)
ABBREVIATIONS
AS ® Air Supply HS ® Hydraulic Supply
IA ® Instrument Air NS ® Nitrogen Supply
PA ® Plant Air SS ® Steam supply
ES ® Electric Supply WS ® Water Supply
GS ® Gas Supply
GENERAL INSTRUMENT SYMBOLS
Primary Location Field Mounted Auxiliary Location Not
(Accessible to Operator) not accessible to operator
Discrete
Instruments
Shared Display or
Shared Control
Computer
Function
PLC
Control Valve
Symbol of
with Actuator Two way valve Two way valve
(Fail open) (Fail closed)
PIC
Example: 204
Introduction 41
SOLVED EXAMPLES
1. Calculate the resistance of a flow metering device whose flow
characteristics are related by q = hn.
Solution: We know that
head ( h )
Resistance (R) =
flow rate ( q )
h Dh dh
R= = = . ...(1)
q Dq dq
Taking Log on both sides for the equation q = hn
Log q = Log [hn]
log q = n log h.
Differentiating w.r.t. q
d d
Log q = [n Log h]
dq dq
1 1 dh
= n. .
q h dq
dh h
Þ = .
dq nq
Comparing with eqn (1)
h
R= Ans.
nq
q = KA 2 g( h1 - h2 )
where, K = flow constant, A = Area of restriction (M2).
Squaring b.s.
q2 = K2A2[2g(h1 h2)].
q.q = K2A22g(h1 h2).
q 2( h1 - h2 )
2 2
=
K A g q
5TH PROOF 4-9-09
Let h1 h2 = h
q 2h
\ 2 2
= .
K A .g q
2h
Since R=
q
q
\ R= Ans.
K A2 . g
2
q = KA 2 g( h1 h2 )
q2 = K2A22g(h1 h2).
q.q = K2A22g(h1 h2).
q 2( h1 h2 )
2
2
= .
K A g q
Let h1 h2 = H.
Capacitance C a Area (A)
q
=
2H
= R3 R=
h
.
LM OP
2
2
K A g q q N Q
2H
\ R= [Expressing q = Q].
q
2H
R= .
Q
Multiplying b.s. by C
2HC
RC =
Q
2HA t = RC
t= 3 C = A [ Area ] .
Q
5TH PROOF 4-9-09
Introduction 43
LM L OP OP
2
N MN QQ
Log q2 = Log Cv2 2 gh 5
dh
Since R= .
dq
2h
R= Ans.
5q
100°C
0°C
t
5 secs
T = 5 secs t = 5 secs.
c(t) = 100°C [1 e5/5]
c(t) = 100°C [1 0.37] [e 1 = 0.37]
c(t) = 63°C.
7. For the process shown in figure calculate the time constant for the
following specifications.
qi
R1
h c
qo
2h
R= = given h = 5 m. q = 0.5 m3/s
q
2´5
\ R= = 20 s/m2.
0.5
Capacitance of the tank = Area = 10 m2.
c = 10 m2
\ Time constant element
20s
T = Rc = ´ 10 m2
2
m
T = 200 secs.
5TH PROOF 4-9-09
Introduction 45
2h 2´3
R= = \ h=3m q = 25 m3/s.
q 25
R = 0.24 s/m2.
Since C = A = 5 m2 .
T = Rc = 0.24 s/m2 ´ 5 m2
T = 1.2 secs
T = 1.2 secs.
9. A liquid level process having first order dynamics is suddenly subjected
to a change in level from 2 m 10 m at time t = 0. If the time constant
of the process is 10 secs, estimate the rise in level after 20 secs.
Solution: From the first order dynamics o/p c(t) = K[1 e t/T] when the initial
value of the output is zero.
If the initial value of the output is not zero, then the above equation may be
modified as
c(t) = Ko + [Ki ko]e t/T
where, Ki = Initial value.
Ko = Final value.
For liquid level process
output h(t) = ho + [hi ho] et/T.
ho = Final value of level = 10 m
hi = Initial value of level = 2 m.
After t = 20 secs.
h(t) = 10 + [2 10] e20/10
h(t) = 10 + [ 8]e2
h(t) = 10 + [ 8][0.13].
h(t) = 10 1.04 = 8.96 m
\ h = 8.96 m.
10. A temperature sensor having time constant characteristics is subjected
to a sudden change in temperature of 25°C 200°C. If it has a time
constant of 5 secs, what temperature will be indicated after 8 secs?
5TH PROOF 4-9-09
k 1000
Wn = = = 100 ´ 103
m 10 ´ 10 3
W n = 316.22 r/s.
(b) Given x = 0.6.
W n = 316.22 r/s.
\ Damped natural frequency
W d = Wn 1 x 2
W d = 316.22 1 (0.6)2
W d = 252.97 r/s.
12. Obtain the time domain specifications of second order under damped
response subjected to step input.
Solution: Mp Tolerance
band
The time domain specifications are: 1.0 2.5%
1. Delay time (td): It is defined as the 0.5
time taken by the system response to
reach 50% of the final value. td
tr
tp
1 + 0.7x ts
td = secs.
Wn
5TH PROOF 4-9-09
Introduction 47
2. Rise time (tr): It is defined as the time taken by the system response to
reach 100% of the final value.
pf
LM
tan 1 1 x 2 OP
tr =
Wd
f=
MN
x.
(radians) .
PQ
3. Peak time (tp): It is defined as time taken by the system response to
reach maximum peck value.
p
tp = (secs).
Wd
4. Setting time (ts): It is defined as time taken by the system response to
reach and settle down within the limits of tolerance band about steady state
value of the input.
4
ts = for 2% Tolerance band
xwn
3
ts = for 5% Tolerance band.
xwn
5. Maximum peak overshoot (mp): It is the maximum value of the
response measured from steady state value.
xp / 1 x 2
MP = e
xp / 1 x 2
% MP = e ´ 100 .
13. Derive the expression for rise time (tr) for the second order under
damped system subjected to step input.
Solution: We know that under damped response
e xwn t
c(t) = 1 sin [wd t + f].
1 x2
tan 1 1 x 2
where, f= .
x
At t = tr c(t) = 1
e xwn tr
1 =1 sin (wd tr + f).
1 x2
5TH PROOF 4-9-09
e xwn tr
sin [wd tr + f] = 0.
1 x2
xw n tr
e
Since is finite
1 x2
sin [wd tr + f] = 0.
This is true if
wdtr + f = p
pf 1 x2
Þ tr = where, f = tan 1 .
wd x
14. Derive an expression for peak time and maximum peak overshoot
shown in time domain specifications.
Solution: To obtain an expression for peak time
d
c(t) = 0.
dt
d e xwn t
1 sin [wd t + f] = 0
dt 1 x2
LM e xw n t
e xwn t
OP
=0
MN 1 x 2
cos ( w d t + f). wd + sin ( w d t + f)
1 x2
[ xw n ]
PQ
e xwn t
Þ [ xwn(sin wdt + f) wd cos (wdt + f)] = 0
1 x2
e xwn t
Since is finite
1 x2
sin ( wd t + f) w wn 1 - x2
= d =
cos ( wd t + f) ywn x wn
1 x2
tan [wd t + f] =
x
5TH PROOF 4-9-09
Introduction 49
1 x2
Since f = tan 1
x
LM 1 x2 OP 1 x2
MN
tan wd t + tan 1
x PQ = x
tan [wd t] = 0
Þ wdt = p
t = tP = p/wd
Substituting t = tP = p/wd in the expression for response
c( t ) t = t = c(t)|max.
p
p
e xwn ´
wd é p ù
c(t)|max = 1 sin ê wd + fú .
2
1 x êë wd úû
xp / 1 x2
e
c(t)|max = 1 sin [p + f].
1 x2
From the under damped pole locations
sin f = 1 x 2
wn
wn 1 – x2
xp / 1 x2
e
\ c(t)|max = 1 [ sin f] xwn
1 x2
-xp / 1- x 2
e é ù
c(t)|max = 1 - ´ ê- 1 - x2 ú
2 ë û
1- x
xp / 1 x 2
c(t)max = 1 + e
MP = c(t)|max 1
xp / 1 x 2
MP = e .
15. A certain oscillatory process has a transfer function.
C( s) K ( s + 1)
= 2 .
R( s) ( s + 2 s + 5)
5TH PROOF 4-9-09
If the steady state value of the response when subjected to unit step
input is 2, find the system gain K.
RS K ( s + 1) UV
T| ( s + 2 s + 5) W|
Solution: C(s) = R(s) 2
K ( s + 1)
C(s) = .
S( s2 + 2 s + 5)
Applying final value theorem.
Lt S. C(s) = 2
s®0
K ( s + 1)
Lt s . =2
s®0 s ( s2 + 2 s + 5)
K
\ = 2 Þ K = 10.
5
16. Obtain the dynamic equations for the mechanical system shown in
figure.
x2
M2
B1 x1
K2 K1 x
M1 F
B2
F M1 M2 K2
5TH PROOF 4-9-09
Introduction 51
F = k1(x x1).
At node (x1)
d2 x1 d
0 = M1 2
+ B1 (x1 x2) + K1(x1 x),
dt dt
d2 x1 d
K1(x x1) = M1 + B1 (x1 x2)
dt 2 dt
d2 x1 d
F = M1 + B1 (x1 x2). ...(1)
dt 2 dt
At node x2
d2 x2 dx d
0 = M2 + B2 2 + B1 (x2 x1) + k2x2. ...(2)
dt 2 dt dt
17. Find the transfer function of lead compensator network.
R1
I(s)
1
Vi(s) CS R2 V0(s)
Vi(s) = I(s)
LM R + R OP 1
N R cs + 1 Q
1
2
Vi(s) = I(s)
LM R + R + R R CS OP
1 2 1 2
N R CS + 1 Q 1
Vo(s) = I(s)[R2].
Vo ( s) R2 [ R1CS + 1]
\ =
Vi ( s) R1 + R2 + R1 R2CS
R2 ( R1CS + 1)
=
L
R + R M1 +
R R CS O
.
N R + R PQ
1 2
1 2
1 2
5TH PROOF 4-9-09
R2
Let T = R1c ; a = (a < 1)
R1 + R2
Vo ( s) a(1 + Ts)
= .
Vi ( s) 1 + aTs
I(s) R1
R2
Vi(s) Vo(s)
1
cs
LM R cs + 1OP 2
N
Vi(s) = I(s) R1 +
cs Q
LM R cs + R cs + 1OP
1 2
Vi(s) = I(s)
N cs Q
LM R cs + 1OP .
2
Vo(s) = I(s)
N cs Q
Vo ( s) R2 cs + 1
=
Vi ( s) R1 cs + R2 cs + 1
R1 + R2
Let b= & T = R2 c
R2
Vo ( s) R2 cs + 1 1 + Ts
= = .
Vi ( s) cs( R1 + R2 ) 1 + bTs
+1
R2
R1
I(s)
1 R2
Vi(s) c1s Vo(s)
1
c2s
5TH PROOF 4-9-09
Introduction 53
NR c s +1 c s Q
1 1 2
Vo(s) = I(s)
LM R c s + 1OP
2 2
N cs Q
2
Vo ( s) ( R1c1 s + 1) ( R2 c2 s + 1)
\ =
Vi ( s) ( R1c1 s + 1) ( R2 c2 s + 1) + R1c2 s
Vo ( s) ( R1 c1 s + 1) ( R2 c2 s + 1)
=
Vi ( s) R1c1 R2 c2 s 2 + s[ R1c1 + R2 c2 + R1 c2 ] + 1
Vo ( s) a(1 + T1 s) (1 + T2 s)
=
Vi ( s) (1 + aT1 s) (1 + bT2 s)
X1( s)
20. Find the transfer function .
F ( s)
K2 x2
x1
K1
M2 M1 F
B2
x1 K1 x2
B2
F K2
M1 M2
5TH PROOF 4-9-09
At node x1
d2 x1
F= + K1 (x1 x2).
dt 2
At node x2
M2 d2 x2 B2 dx2
0= + + K2x2 + K1(x2 x1).
dt 2 dt
Applying Laplace transform
F(s) = X1(s) [M1s2 + K1] K1X2(s)
0 = X2(s) [M2s2 + B2s + K2 + K1] K1X 1(s)
K1
\ X 2(s) = . X1(s).
M2 s2 + B2 s + K 2 + K1
2
K . X1 ( s)
1
\ F(s) = X1(s)[m1s2 + K1]
M2 s2 + B2 s + K2 + K1
LM( M s 2
+ K1 )( M2 s2 + B2 s + K2 + K1 ) K
2 OP
F(s) = X (s) N Q
1
1
1
M2 s2 + B2 s + K2 + K1
X1( s) M2 s 2 + B2 s + K1 + K2
= .
F ( s) 2
( M1 s 2 + K1 )( M2 s 2 + B2 s + K1 + K2 ) K
1
QUESTION SET
1. With a neat diagram explain the temperature control system indicating
all the elements of process control loop.
2. Explain the importance of elements of process dynamics in modelling
of processes.
3. What is the significance of time constant in a first order process?
4. Obtain the different approximations of dead time element.
5. Explain clearly how damping affects the nature of response of
oscillatory process.
6. Derive an expression of oscillatory process (under damped) for sudden
changes in input. Mark all the time domain specifications and explain
clearly each one of them?
5TH PROOF 4-9-09
Introduction 55
7. Obtain the mathematical models for (a) Liquid level system (b)
Thermal system (c) Pneumatic system.
t1( s)
8. A time constant process is shown in fig. Find the transfer function .
Qi( s)
C = Hydraulic
qi
capacitance
c R1
h
qo1
R2
qo2
9. For an interacting system, draw the block diagram schematic and hence
write the mathematical model.
qi
c1 R1 c2 c3 R3
R2
h1 h2 h3
qo
10. For the above interacting system (Prob 9) draw equivalent electrical N/
W and compare the mathematical models.
t13 ( s)
11. Find the transfer function for three tank non-interacting system .
Q i ( s)
qi
c1 R1
h1
c2 R2
h2
R3
h3 c3
qo
5TH PROOF 4-9-09
(a) (b)
(c) (d)
(e)
LC SP FC
104 100
(a)
FT
(b)
TDI
TT TT
103
(c) TIC
TT
101
(d)
PI PI PI
(A) 102 (B)
5TH PROOF 4-9-09
Introduction 57
Gas Flow
PT
FT TT
108
(e)
FY FIC
111
R L
f
K
+ C
V i
–
F x
R(s) 16
C(s)
– – s(s + 4)
as
What should be the value of a such that the response for unit step input
will give an overshoot of 1.5%?
19. A thermometer is subjected to a sudden temperature change of 35°C
® 115°C. What would be the time constant, if it indicates 87.63°C
after 6 secs?
20. For a second order system having transfer function
C( s) 8
= 2
R( s) s + 8 s + 25
Find the time domain specifications.
+0)26-4
D p - Pmin
%P= ´ 100
Pmax - Pmin
D p = 12 mA
Pmin = 4 mA, Pmax = 20 mA
12 - 4
\ %P = ´ 100
20 - 4
8
= ´ 100 Þ 50 %
16
Similarly, the standard pneumatic signal transmission range is 315 psi.
Input e
controller P
output
CLASSIFICATION
Modes of
Control Actions
Discontinuous Continuous
Modes of Modes of
Action Action
OFF
0%
error
(e)
Neutral Zone
FIG. 2.1
Applications
The main application of two position control are:
(1) Room heating or air conditioning systems
(2) Liquid bath temperature control
(3) Level control in large volume tanks.
R|100 % e > e2
P= S| 50 % - e1 < e < e2
T 0% e < - e1
5TH PROOF 4-9-09
As long as the error is between e2 and e1 of the set point, the controller stays
at a nominal setting of 50% if the error exceeds the set point by e1 or more, the
output is increased by 100 %. If it is less than the set point by e1 or more the
controller output is reduced to zero.
Controller
0/P
P 100%
50%
0
error
– e1 0 + e1 (e)
FIG. 2.2
floats at whatever setting it was when the error went to zero. And whenever
error occurs the controller output again begins to change.
dp
= ± KF |e|> Neutral Zone
dt
dp
where = rate of change of controller output,
dt
KF = Rate constant [% / sec]
The specific controller output is obtained by integrating the above
equation.
P = ± KF t + P (0)
where P (0) refers to initial setting of the controller when the error is zero. The
characteristics of single speed floating controller mode is
dp
dt KF
– KF
error
0 (e)
FIG. 2.3
Between the instants 0 t1 the error is within the limits of Neutral Zone.
Hence the controller output floats at whatever setting it was when the error
went to zero i.e. P (0). From t1 to t2 error exceeds the neutral zone and the
controller output changes at a fixed rate, ± KF t, depending upon Nature of
Process, i.e. reverse or direct action.
Neutral
Zone
Error t
(e)
Floating
Action
+KFt –KFt
Controller Floating Floating
O/P P(0)
Action Action
t
t1 t2 t3 t4
FIG. 2.4
dp
= ± KF i e > ei
dt
Application
The single speed or multiple speed floating action can be employed with
Controller
O/p
+ K2
+ K1
– K1
– K2
error
– e2 – e1 0 e1 e2
FIG. 2.5
general flow control system or pressure control system with small capacities
and having self-regulation to some extent with small dead time or Lag.
\ E(S) Kp P(S)
Controller
o/p
100 % Saturation
P 100%
0% Saturation
0%
error
P.B (e)
FIG. 2.6
error
(e)
t
t = t1
Controller
o/p
P
t
t = t1
FIG. 2.7
10 0
as P . B = . Thus by proper choosing of Kp the P . B can be changed
Kp
In Figure 2.7, assume a step change in deviation,
Let e = A [step change].
P = Kp e,
Þ P = K p A.
5TH PROOF 4-9-09
Offset
One of the main disadvantage of proportional control action is that it produces
a permanent residual error in the operating point of the controlled variable
when a process load change results in error. This residual error is known as
OFFSET.
Illustration of offset
Consider a first order process having process equation or transfer function
C ( s) 1
= .
R ( s) 1+ TS
FIG. 2.8
S . R ( s)
ess = Lt
s®0
1 + G ( s) H ( s)
A
Let R (s) = [step change]
S
A
s/
s/
ess = Lt
s®0 Kp
1+
1 + TS
A
e ss = ...(2)
1+ Kp
5TH PROOF 4-9-09
This steady state error is not zero but inversely proportional to gain Kp.
Thus, this control action fails to make the error zero whenever the o /p changes
due to changes in input or Load (process) variables. This is known as offset.
From eqn. (2),
1
offset µ ...(3)
Kp
Controller
dp
o/p ae
dt
kF
P
0
– kF
error
o
(e)
FIG. 2.9
dp
µe
dt
dp
= KI . e ...(4)
dt
For particular error, the output will begin to increase or decrease at a rate
KI % /second for every 1% of error.
The specific controller output is obtained by integrating eqn. (4).
P = KI z e d t + P(0) ...(5)
1
Ti = ...(6)
KI
P=
1
Ti z e d t + P (0) ...(7)
1
P (s) = . E (s).
Ti s
E(s) 1 P(s)
Ti s
P=
1
Ti z A. d t
Þ P=
LM A OP t
NT Qi
From Figure (2.10), the main disadvantage of integral action is that, its
response to error is slow or sluggish. But however, due to its basic nature of rate
of change of output it eliminates the error and brings back the controlled
variable to the set point value.
5TH PROOF 4-9-09
error
e
A
A
Controller t
Ti
o/p P
t
t = t1
FIG. 2.10
Since the rate of change of controller output can be reset at any instant of
time t by changing Ti. This mode is often referred to as RESET
CONTROLLER MODE. The integral time Ti is also known as Reset time. Let
us see the effect of integral control action on error by taking a first order
process.
1 P(s) 1
R(s) C(s)
Tis 1 + Ts
FIG. 2.11
S . R ( s)
ess = Lt
s®0
1 + G ( s) H ( s)
A
Let R(s) = [causing step change in deviation],
S
A
s/
ess [error] = Lt /s
s®0
1
1+
b
Ti s 1 + T s g
A
ess [error] = = 0.
1+ ¥
Thus it can be established that integral control action successfully
eliminates the error, but in the process its response is sluggish. The integral
action, due to its basic nature, might introduce oscillations in the controlled
variable about the set point value, whenever it deviates from the set point
value.
5TH PROOF 4-9-09
de
P = KD ...(8)
dt
Defining T d derivative time as
Td = KD [derivative gain]
de
P = Td ...(9a)
dt
Applying L .T to equation (9a)
P (s) = Td s . E (s)
E(s) P(s)
Td s
Characteristics
® If the error is step change i.e. sudden, the controller output will be zero
because it responds only to rate of change of error. When e = A (step change)
d
P = Td (A) = 0.
dt
Considering the error to be rate of change i.e. linearly changing w.r.t. time,
e = At [Ramp change in deviation]
d
P = Td [A t]
dt
Þ P = Td A [Controller output]
As seen from Fig. (2.12), the derivative action is anticipating the nature of
error and giving out its output much before the error has occurred. Because of
this action, it is also called as Anticipatory Controller Mode.
Further, since the controller output is proportional to rate of change of
error it is also called as Rate Controller Mode and Td is called as rate time.
Unlike integral action, the derivative actions response is much faster. For a
very rapid change of error there could be a large sudden changes of controller
output, which could lead to un-stable condition. This mode is suitable if the
process reaction rate itself is very sluggish, as it considerably improves the
speed of response.
5TH PROOF 4-9-09
Controller
o/p
P TdA { t
t = t1
FIG. 2.12
P = Kp e +
Kp
Ti z e d t + P (0) ...(9)
E(S) 1 P(S)
Kp 1 +
Tis
error
(e)
A
t
t = t1
P-only
t
t = t1
I-only
t
t = t1
P + I
n
only A ctio AKp
P +I P = KpA + t
Ti
{
t
t = t1
FIG. 2.13
e = A [step change]
P = K pe +
Kp
Ti z edt.
P = Kp A +
Kp
Ti z Adt.
5TH PROOF 4-9-09
F AK I t
P = Kp A +
" ! GH T JK i
p
Therefore, it may be stated that the integral action eliminates the inherent
offset of proportional action. To Illustrate consider a first order process whose
transfer function is
C ( s) 1
=
R ( s) 1 + Ts
P + I controller (Process)
{
{
1
R(s) Kp 1 + 1 C(s)
Tis 1 + Ts
S . R( s)
ess = Lt
S®0
1 + G ( s) H ( s)
A
Since R (s) =
S
LM F 1 I OP LMF 1 I OP
MN GH JK PQ MNGH1 + T JK PQ
G (s)H (s) = K p 1 +
Ti s s
A
s/
Lt /s
ess = S®0
F IF 1 I
1+ Kp GH 1+
1
Ti s JK GH 1 + Ts JK
A A
ess = = =0
e
1 + K p + ¥ (1) j ¥
P = Kp sin w t +
Kp
Ti z sin w t. d t .
F - K I cos w t
P = K p sin w t + GH w T JK p
i
...(10)
We know that the sum of sine wave and co-sine wave is another sine wave
having different magnitude and phase angle
A2 + B2 sin w t + Tan - 1
LM B OP
A sin w t + B cos w t =
N A Q
In equation (10)
- Kp
A = Kp and B =
wTi
therefore,
LM F - K/ I OP
F - K I sin Mw t + Tan GG w T JJ P
p
2
P = eK j + G
H w T JK MM
2
JJ PP
p -1 i
p
i GG K/
H K PQ
p
MN
e K j + FGH w T IJK sin LMNw t - Tan w1T OPQ .
2
2 K p -1
P = p
i i
As seen from eqn. (11) when error (e) is sin w t the phase of the P +I
F 1 I . This implies the overall effect on
controller output lags by Tan1 GH w T JK i
transient state or speed of response in the process is slowed down due to the
rate of change of controller output of integral mode. Hence this mode can be
optimally used for the improvisation of steady state response characteristics
only in a process.
The characteristics of P + I control action are also comparable with Lag
compensator (or) low pass filter.
de
P = K p e + K p Td + P0 ...(11)
dt
5TH PROOF 4-9-09
Characteristics of P + D controller
We cannot consider a step change in deviation because the derivative action
does not respond to such errors. Hence considering rate of change of error.
At an instant of time t = t1, the error is changing at a continuous rate. As
shown in the figure, the derivative action first sends its output in its in imitable
anticipatory mode which is followed by proportional action. Such a response
would lead to large instability in a process which may result in error increasing
beyond proportions rather than becoming zero. If we consider this control
action for level control system discussed in earlier chapters, when the liquid
level in the tank starts rising beyond the actual set point resulting in an error
shown in Fig (2.14), the derivative mode opens the final control element to its
maximum capacity causing large quantity of outflow, followed by proportional
action further opening the control valve resulting in outflow greater than
inflow. This condition results in large error and High instability. Since this
mode is incapable of improving the steady state behaviour i.e. eliminate steady
state error and inherently exhibits good speed of response, we can use it to
processes having slow response characteristics, to effectively increase their
speed of response. Referring to the figure below
error At
(e)
t
t = t1
P-only
t
t = t1
D-only
t
t = t1
P+D
Action
P+
D { Ac
tion
t
t = t1
FIG. 2.14
5TH PROOF 4-9-09
Controller Process
A
let R (s) = [ Ramp change in deviation].
S2
A
/s 2
/s
ess = Lt
S®0
K p 1 + Td s
1+
1 + T s
A
= Lt
S®0
s K p 1 + Td s
s+
1 + T s
A
\ ess = =¥
o
Therefore this mode has no effect on the steady state response. Whether or
not it has any effect on transient state behaviour can be established by
considering sinusoidal deviation and respective analysis
Let e = sin w t
d
P = Kp sin w t + Kp Td sin w t.
dt
P = K p sin w t + w K p Td cos w t.
2
P= K p2 + wK p Td . sin [wt + Tan1wTd]
P= e
K 2p + wK p Td j2
. sin w t + Tan - 1 w Td ...(12)
As seen from equation (12), for sinusoidal deviation the phase of the
controlled variable leads by Tan-1 wTd. Hence a P+D control action can be
employed with processes having slow response characteristics, which can be
conveniently improvised by this mode of control action.
Since the phase of controller output is leading by Tan1 wTd for sinusoidal
deviation its characteristics are similar to Lead Compensator and High Pass
Filter.
5TH PROOF 4-9-09
P = Kpe +
Kp
Ti z e dt + Kp Td + P(0) ...(13)
E(S) 1 P(S)
Kp 1 + + TdS
TiS
P = Kp At +
Kp
Ti z At dt + Kp T d
d
dt
(At)
P = (KpA) t +
LM A K OP t p 2
+ [A Kp Td] ...(14)
N 2T Q i
t
t = t1
I-only
P-only
P
D-only
t
t = t1
tio .
Ac I. D
}
n
P
P.
t
t = t1
FIG. 2.15
S × R( s)
ess = Lt
S®0 1 + G( s) H( s)
A
Let R(s) =
S2
A
S
ess = Lt
S®0 FG 1 IJ FG 1 IJ
1+ Kp 1+
H Ti S
+ Td S
K H 1 + Ts K
A
ess =
L K
0 + M0 + i
+ 0 [1]
OP
N T i Q
(PID controller) Process
1 1
R(s) Kp 1 + + Tds C(s)
Tis 1 + Ts
5TH PROOF 4-9-09
ATi
ess = since Ti << Kp
Kp
ess = 0
let e = sin wt
P = Kp sin wt +
Kp
Ti z sin wt dt + Kp Td
d
dt
sin wt
FG - K IJ cos wt +w K T
p
P = Kp sin wt +
H wT K i
p d cos wt
LM
P = Kp sin wt + w K p Td -
Kp OP cos wt
N wTi Q
2
æ Kp ö é æ 1 öù
P= K p2 + ç wK p Td - ÷ sin ê w t + Tan-1 ç wTd ú ...(15)
è wTi ø êë è wTi ø÷ úû
Therefore from eqn. (15), that phase of the controller output leads by
FG
Tan1 wTd -
1 IJ
H wTi K
, which implies it is capable of improving the transient state
SOLVED EXAMPLES
Type-1 Problems
1. Indicate the type of control action and obtain the controller parameters
Kp, T d and Ti for the controller represented as Shown in figure.
E(s)
K + P(s)
+
1
Ts
LM 1 OP E(s)
Solution:
N Ts Q
P(s) = k +
R L 1 OUV E(s).
P(s) = Sk M1 +
T N kTs PQW
It is P + I controller comparing with
P(s) = KP 1 +
LM 1 OP
N Ti s Q
KP = k Ti = kT
E(s) 1
1 + T2s
2. T1s P(s)
RS
P(s) = (1 + T2 s).
1 UV
E(s)
T T1 s W
ì 1 T sü
P(s) = í + 2 ý E(s)
î T1 s T1 s þ
ìï T é T1 1 ù üï
P(s) = í 2 ê1 + . ú ý E(s)
ï T ëê T2 T1 s úû ï
î 1 þ
It is P + I controller.
T2
Kp = Ti = T2
T1
5TH PROOF 4-9-09
é T T s2 s [T1 + T2 ] 1 ù
P(s) = ê 1 2 + + ú E(s)
ë Ts Ts Ts û
éT T T +T 1ù
P(s) = ê 1 2 . s + 1 2 + ú E(s)
ë T T Ts û
ìï T + T é T TT T 1 ù üï
E( s). í 1 2 ê . 1 2 .s + 1+ . úý
îï T ë T1 + T2 T T1 + T2 T s û ï
þ
ìï T + T é æ T1T2 ö 1 ù üï
\ P(s) = í 1 2 ê1 + ç T + T ÷ s + (T + T ) s ú ý E(s)
T ë è 1 2ø
îï 1 2 ûþ ï
It is a PID controller comparing with
ì é 1 ùü
P(s) = í K p ê1 + + Td s ú ý E(s)
î ë Ti s ûþ
T1 + T2 T1T2
Kp = Ti = T 1 + T 2 Td =
T T1 + T2
E(s) 1 + T1s 1
4. 1+
T2s
P(s)
é æ 1 öù
P(s) = ê(1 + T1s) ç1 + ÷ ú E(s)
ë è T2sø û
é 1 T sù
P(s) = ê1 + + T1s + 1 ú E(s)
ë T2 s T2 s û
éæ Tö 1 ù
P(s) = êç1 + 1 ÷ + + T1 s ú E(s)
ëè T2 ø T2 s û
éT + T 1 ù
P(s) = ê 1 2 + + T1s ú E(s)
ë T2 T2 s û
ìï T + T é T2 1 T2 ùüï
P(s) = í 1 2
ê1 + . + .T1 s úý E(s)
ï T2 ëê T1 + T2 T2 s T1 + T2 ûúï
î þ
5TH PROOF 4-9-09
ìï T + T é 1 æ T1T2 ö ù üï
\ P(s) = í 1 2 ê1 + (T + T ) s + ç T + T ÷ s ú ý E(s)
ï T2 ë è 1 2 ø ûï
î 1 2 þ
It is a PID controller.
T1 + T2 T1T2
Kp = Ti = T 1 + T 2 Td =
T2 T1 + T2
E(s) T
+
+
T1s P(s)
5. +
T2s
P(s) = [T + T 1s + T 2s]E(s)
= [T + (T1 + T 2]E(s)
ì é æ T + T2 ö ù ü
P(s) = íT ê1 + ç 1 ÷ s ú ý E(s)
î ë è T ø ûþ
It is a P + D controller comparing with
P(s) = [Kp(1 + T d s)]E(s)
T1 + T2
Kp = T Td =
T
E(s) T1s 1
6. 1+
T2s
P(s)
é T sù
P(s) = êT1s + 1 ú E(s)
ë T2 s û
ìT é T2 ùïü
P(s) = ïí 1 ê1 + T . T1 s ú ý E(s)
ï T ë ûï
î 2 1 þ
ìT ü
P(s) = í 1 [1 + T2 s]ý E(s)
T
î 2 þ
It is a P + D controller.
Kp = T1/T2 Td = T 2
5TH PROOF 4-9-09
7.
2
E(s) (1 + Ts) P(s)
d2 e de
P = T2 + 2T + e.
dt 2 dt
de
1et = eP.
dt
deP
P = T2
dt
+ 2T . eP + ò e dt .
P
é 1ù
EP (s) êT 2 . s + 2T + ú .
ë sû
ì é T 1 ùü
P(s) = í2T ê1 + . s + ý EP (s)
î ë 2 2Ts úûþ
It is a PID controller.
\ Kp = 2T Ti = 2T
T
Td =
2
1
8. E(s) T1s P(s)
T2s
éT s ù é T1 ù
P(s) = ê 1 ú E(s) Þ P(s) = ê T ú E(s)
ë T2 s û ë 2û
It is proportional controller
Kp = T1/T2
5TH PROOF 4-9-09
T1
E(s) +
P(s)
9.
+
T2s
ì é T ùü
P(s) = íT1 ê1 + 2 . s ú ý E(s)
îï ë T1 û ï þ
It is P + D controller
Kp = T 1 Td = T2 /T1
d
Put s=
dt
de de p
P = e + T1 + T 22
dt dt
de
let = ep.
dt
de p
P= ò e dt + T e
p 1 p + T 22
dt
It is PID controller.
Applying laplace transform.
é 1 ù
P(s) = êT1 + + T22 s ú E(s)
ë s û
ìï é 1 T 2 ù üï
P(s) = íT1 ê1 + + 2 . s ú ý E(s)
îï ë T1s T1 û þï
T22
Kp = T 1 Ti = T 1 Td =
T1
5TH PROOF 4-9-09
Type-2 Problems
11. In an integral controller with reset time set to 0.6 minutes, what will be the
phase of controller output for sinusoidal deviation?
1
Solution: P=
Ti ò e dt.
Since Ti = 0.6 min. e = sin wt
1
P=
0.6 ò sin wt. dt
- cos wt
P=
w ´ 0.6
Since cos q = sin [90 q]
cos q = sin [q 90°]
sin [ wt - p /2] 1
\ P= Þ sin (wt p/2)]
0.6w 0.6w
The phase of the controller o/p Lays by p/2 degrees.
12. Determine the phase of the controller output of derivative controller for
sinusoidal deviation?
de
Solution: P = Td Since e = sin wt
dt
d
p = Td (sin wt) = wTd cos wt
dt
P = wTd sin (wt + p/2). [3 sin (90° + q) = cos q]
The phase of the controller o/p leads by p/2.
13. The input error signal for P + I controller is sinusoidal in nature. Prove
that the phase lag is a function of reset time?
Kp
Solution: P = Kpe +
Ti ò e dt
e = sin wt
Kp
p = Kp sin wt +
Ti ò sin wt. dt
é-Kp ù
p = Kp sin wt + ê ú cos wt
ë wTi û
5TH PROOF 4-9-09
é æ - Kp öù
ê 2
ç ú
æ -Kp ö ê wt + tan -1 ç wTi ÷ ú
\ p= K p2 +ç sin ÷
è wTi ÷ø ê èç Kp ø÷ úû
ë
2
æ Kp ö é æ 1 öù
\ P= K p2 +ç ÷ sin ê wt - tan - 1 ç ú
è wTi ø ë è wTi ø÷ û
Thus the phase lag q = tan 1 (1/wTi) is a function of reset time (T i).
14. A P + I controller subject to sinusoidal deviation, is tuned to T i = 3 min
and Kp = 10. Determine its output when f = 50 Hz.
Kp
Solution: P = Kpe +
Ti ò e dt.
Given e = sin wt; kp = 10; Ti = 3 min
10
P = 10e +
3 ò e dt
( - 3.33)
Þ P = 10 sin wt + cos wt.
w
2
æ 3.33 ö é æ - 3.33/ w ö ù
P = 102 + ç .sin ê wt + tan -1 ç
è w ÷ø ë è 10 ÷ø úû
æ 11 ö é 3.33 ù
P = 102 + ç 2 ÷ .sin ê wt - tan -1
èw ø ë 10 w ûú
w = 2pf = 2 ´ p ´ 50 = 314.15 radians.
11 é 3.33 ù
P= 100 + .sin ê314.15t - tan -1 ú
(314.15)2 ë 10 ´ 314.15 û
15. What is the phase of the controller output for sinusoidal error for P+D
controller?
de
Solution: P = Kpe + KpT d
dt
e = sin wt
d
P = Kp sin wt + KpTd sin wt
dt
5TH PROOF 4-9-09
1 E ( s)
P(s) =
T 2 s2
Put s = jw
P ( jw) 1 1 1 1
F( jw) = = = . .
E( jw) ( jw)2 . T 2 T 2 jw jw
1 1 + jw 1 + j 0
F( jw) = .
T 2 0 + jw 0 + jw
1 12 + 02 12 + 02 1
|F( jw)| = . =
T 2 02 + w2 2
0 +w 2 w T22
0° 0°
Þ -1 = = 180°.
2 T an ¥ 2 ´ 90°
Therefore the phase of the controller 0/p lags by 180°.
17. Show that for a PID controller the phase of its output is a function of rate
time and reset time for sinusoidal deviation?
Kp de
Solution: P = Kpe +
Ti ò e dt + K T p d
dt
e = sin wt
Kp d sin wt
P = Kp sin wt +
Ti ò sin wt + K T p d
dt
5TH PROOF 4-9-09
æ - Kp ö
P = Kp sin wt + ç cos wt + wKpTd cos wt.
è wTi ÷ø
æ Kp ö
P = Kp sin wt + ç wK pTd - cos wt
è wTi ÷ø
é æ Kp öù
2 ê ç wK T - ú
æ Kp ö
ê wt + tan - 1 ç
p d
wTi ÷ ú
P= K p2 + ç wK pTd - .sin ÷
è wTi ÷ø ê çè Kp ÷ø ú
ë û
2
æ Kp ö é æ 1 öù
P= K p2 + ç wK pTd - .sin ê wt + tan -1 ç wTd -
è ÷
wTi ø ë è wTi ø÷ ûú
æ 1 ö
f = tan 1 ç wTd -
è wTi ø÷
The phase of the controller output is a function of rate time (Td) and reset
time (Ti).
18. A certain P + D controller has Td = Kp = 1. Find the phase of its 0/p
when w = p/3?
Solution: From prob (15). The phase of P + D controller leads by tan 1(wTd).
\ Given w = p/3 Td = 1
tan1 (p/3) = 89°.
Type-3 Problems
19. A pneumatic controller is employed to control the valve plug position
between 0 10 mm. It the range of the controller output is 3 15 psi.
What pressure will bring the valve\plug to a position of 8.5 mm.?
Solution: Since controller final control element operation is linearly
dependent
Y = MX + C Y = valve plug position
0 = M.3 + C X = controller 0/p
10 = M.15 + C
Solving for M & C
M = 0.833 C = 2.5.
\ 8.5 = 0.833x 2.5
8.5 + 2.5 = 0.833x x = 13.20 psi
5TH PROOF 4-9-09
20. In the above problem express 13.20 psi as percentage controller output in
the range of 3 15 psi?
Solution: Using the formula
DP - Pmin
%p = ´ 100
Pmax - Pmin
13.20 - 3
= ´ 100
15 - 3
% P = 85%
42 P-4
=
100 20 - 4
42 ´ 16
+ 4 = P.
100
\ P = 10.72 mA
\ Y = 28.75(I ) + 25
5TH PROOF 4-9-09
\ I = 9.91 mA
DP - Pmin
(b) %P = ´ 100
Pmax - Pmin
9.91 - 4
\ %P = ´ 100 = 36.9%
20 - 4
\ % P = 36.9%
54.78 DP - Pmin
Solution: =
100 Pmax - Pmin
54.78 DP - 3
\ =
100 15 - 3
\ DP = 9.57 psi
24. In a level control system, the controller output 3 15 psi signal to the
final control element to monitor liquid level between 2.5 m 22.5 m.
What should be the output of controller to maintain level at 19.75 m?
Solution: Using linear dependence equation
Y[level] = M[X (pressure)] + C
2.5 = 3M + C
22.5 = 15M + C
Solving for M and C
M = 1.67 m/psi C = 2.5 m
\ 19.75 = 1.67x 2.5
19.75 + 2.5
\ x= = 13.32 psi
1.67
Therefore the rewired pressure is 13.32 psi
13.32 - 3
%P = ´ 100
15 - 3
% P = 86%
5TH PROOF 4-9-09
26. In a pressure control system a floating controller with single speed mode
is employed to control pressure in a gas vessel. The controller output
when the process error is within the limits of neutral zone is 35% If the
rate constant KF = +2% per second,
(a) What will be the controller output after 10 seconds.
(b) Find the time when the output saturates
Solution: (a) For floating controller
p = KFt + P(0).
At t = 10 seconds.
P = +2 ´ 10 + 35% = 55%
(b) At p = 100% [output saturation]
100 = +2 ´ t + 55
100 - 55
t= = 22.5 seconds.
2
5TH PROOF 4-9-09
27. In the level control system shown in figure valve A is linear one with valve
coefficient of 20 m3/h per percent controller output. In the initial
controller output for inlet flow rate of 300 m2/h is 15%. By how much
percent the valve A should be opened/closed if the inlet f low rate
suddenly becomes 600 m3/h. Also calculate offset for setting of
proportional gain (kp) at 10% per %?
Inlet
[Valve 'A']
h
Outlet
m3
Solution: Let x % be the new controller output; 20 / % ´ x % = 600 m3/h.
h
600
\ x= = 30%
20
x = 30%
Since only proportional gain setting is given the controller is proportional
controller
P = K p e + P0
30% = 10 e + 15%
30 - 15
= e = 1.5%
10
28. Find the response for PI controller with Kp = KI = 1 for the error as
shown in figure for P(0) = 10%
t
0 1 3
P = Kpe + KpKI ò e dt
For 1>t>0 ; e=t
3>t>1 ; e=1
t >3 ; e=0
5TH PROOF 4-9-09
Consider 1>t>0
p = 1.t + 1 ò t + 10%
t2
p=t+ + 10%
2
At t = 1 sec.
1
P=1+ + 10 + 11.5%
2
For 3>t>1
P=1+1 ò 1 dt + 11.5%
P = 1 + t|31 + 11.5%
P = 1 + (3 1) + 11.5%
P = 14.5%
For t > 0e = 0
\ P = 14.5%
29. Find the P + D controller output with Kp = KD = 1 for the error shown in
Pwb (28)?
de
Solution: P = Kpe + KpKD
dt
Considering the same initial controller output as in Problem. (28)
For 1>t>0; e=t
d
P=t+ (t) + 10%
dt
At t = 1 sec
P = 1 + 1 + 10% = 12%
For 3>t>1
P = 1 + 0 + 12% = 13%
For t > 3e = 0
\ P = 13%
5TH PROOF 4-9-09
30. Find P + I + D controller output for same settings as in pwb (28) & (29)?
Solution: For 1>t>0 ; e=t
1
t2
P=t+ + 1 + 10%
2 0
At t = 1 sec
P = 1 + 0.5 + 1 + 10% = 12.5%
For 3>t>1; e=1
P = 1 + t°|13 + 0 + 12.5%
1 + (3 1) + 12.5%
P = 15.5%
For t >3; e=0
\ P = 15.5%
QUESTION SET
1. What are the types of discontinuous controller modes?
2. Explain the characteristics of two position controller mode?
3. What is floating action? Illustrate with an example?
4. What are the different types of continuous controller modes?
5. Explain the characteristics of
(a) P controller
(b) I controller
(c) D controller
6. Compare and contrast the relative merits and de-merits of above control
actions?
7. For sinusoidal deviation obtain the magnitude and phase angle for P, I
and D actions?
8. What are composite control actions? Explain each one of then?
9. Compare the relative advantages and disadvantages of composite control
actions?
10. Obtain the phase relation for sinusoidal deviation for composite modes?
!
+0)26-4
U(s)
N(s)
R(s) E(s) +
Gc(s) Gp(s) + C(s)
–
Controller Process
Hm(s) = 1
FIG. 3.1
C( s) GP ( s)GC ( s)
=
R( s) 1 + GP ( s)GC ( s)
C(s) = R(s)
LM G ( s)G ( s) OP .
P C ...(1)
N1 + G ( s)G ( s) Q
P C
Similarly, the transfer function with respect to load variable u(s) can be
obtained by setting R(s) = 0.
u(s)
N(s) C(s)
+
– [GP(s) GC(s)]
5TH PROOF 4-9-09
u(s) 1
N(s) 1 + GP(s) GC(s) C(s)
Therefore,
C ( s) N ( s)
=
u( s) 1 + GP ( s)GC ( s)
C(s) = u(s)
LM N ( s) OP . ...(2)
N1 + G ( s). G ( s) Q
P C
C(s) = R(s)
LM G ( s)G ( s) OP + u(s) LM N ( s) OP .
P C
...(3)
N1 + G ( s)G ( s) Q
P C N1 + G ( s)G ( s) QP C
S. R( s)
ess[Static error] = Lt . ...(8)
s ® o 1 + GC ( s)GP ( s)
In Eq. 8, if the change in set point is sudden in nature [Step change] i.e.,
A
R(s) = [A = Step magnitude]
S
A
s.
s
ess = Lt
s ® o 1 + G ( S )G ( S )
C P
A
ess|static error = . ...(9)
1 + Lt GC ( s)GP ( s)
s®o
5TH PROOF 4-9-09
A
Let R(s) = [Steadily changing w.r.t time or Ramp change]
s2
A
s.
s2
ess|Velocity error Þ Lt
s ® o 1 + GC ( s)GP ( s)
A
ess|Velocity error Þ . ...(10)
Lt s + Lt s. GC ( s)GP ( s)
s® o s® o
From Eq. (10) it is clear that the magnitude of error for ramp change in
input predominantly depends on.
Lt S. GC(s)GP(s). This error constant is termed as Velocity error Constant
s® o
abbreviated as Kv.
Similarly, by considering non-linear (or) acceleration type of input, the
error due to such input may be defined as Steady State acceleration error.
The non-linear/Acceleration type of input may be considered as a
A
parabolic input where R(s) = .
s3
A
s/ . 3/
\ ess|Acceleration error = Lt s
s® 0 1 + GC ( s)GP ( s)
5TH PROOF 4-9-09
A
ess|Acceleration error = .
Lt s + Lt s2GC ( s)GP ( s)
2
s ®0 s ®0
KV = Lt SGC(s) GP(s).
s®0
And these error constants are further the function of controller i.e., GC(s).
A controller, properly chosen and tuned definitely improves the transient
and steady state responses of the process.
3.2.3 OFFSET
In many of the process industries, the changes in load variables or
disturbances affect the performance of process or plant. This error which
results due to changes in load variables is termed as offset. This offset error
should not be confused with that of offset of proportional controller. Here, it is
a special type of error which occurs mainly because of changes in load
variables.
In eqn. 7
ess|offset = Lt
LM S.U ( s). N ( s) OP . ...(11)
s ®0
N1 + G ( s)G ( s) Q
C P
Control
qi valve
Inlet LC Level Controller
h LT (Level transducer)
u
Outlet
FIG. 3.2
V(s) 1
cs
Qi(s) 1 –
cs + H(s)
FIG. 3.3
5TH PROOF 4-9-09
U(s) 1
Q(s)
cs
Set point ON –
1
+ H(s)
R(s) – OFF cs
Controller Process
FIG. 3.4
In Fig. 3.2, when the valve is open i.e., [H(s)(h) < (r)(R(s)] [level less than
set point value], where
h = o/P (level) and r = set point.
dh
qi u = c [h < r + N]
dt
N = [Positive Limit of Neutral Zone]
When the valve is closed i.e.,
h > r, qi = 0
dh
u = c [h > r N] [N = ve limit of neutral zone]
dt
When the level is exactly at the set point value i.e., when h = r, the controller
oscillates continuously between the ON and OFF states, explained analytically
by integrating the above two equations.
z dh
dt
dt = z qi u
c
dt
qi u
h(t) = t + h01 [h01 = integral constant.] ...(13)
c
Similarly,
z dh
dt
. dt = z u
c
. dt
u
h(t) = . t + h02 [h02 = Integral constant.] ...(14)
c
From equations 13 and 14, the level rises and falls continuously. within the
limits of neutral zone ±N when the error is zero, i.e., when c = r.
5TH PROOF 4-9-09
Level
h(t) Limits of
Neutral
+N –u t zone.
qi – u t c
O
c
–N
Time
FIG. 3.5
Controller Process
Qi(s) 1
Kp + H(s)
– cs –
U(s) 1
cs
FIG. 3.6
Qi(s) 1
Kp H(s)
– cs
FIG. 3.7
Let U(s) = 0.
The overall transfer function
KP KP
H( s) cs K P C
= = = . ...(15)
Q ( s) KP cs + K P K
1+ S+ P
cs C
The nature of response with proportional controller can be obtained by
applying inverse laplace transform to eqn. 15.
æ KP ö
K ç ÷t
\ h(t) = P . e è C ø
. ...(16)
C
Comparing e (KP/c)t » e t/T.
1 K C
= P Þ T=
T C KP
The time constant element (T), which indicates the speed of response is
inversely proportional to proportional gain. Thus, as KP , time constant (T)
¯ and speed of response increases.
h(t)
KP
[ K P1 > K P2 ]
c
– KP t
KP
h(t) = e c
c
K
KP = KP2
P
K
=
P1
FIG. 3.8
S. R( s)
ess|Static error = Lt
s ®0 1 + GC ( s)GP ( s)
1
Gc(s) = KP and GP(s) = .
cs
Since u(s) i.e., outlet flow rate is assumed to be zero.
R(s) = Qi(s) = [continuously changing quantity with respect to time]
A
Thus, Qi(s) = [Ramp change].
s2
A
s.
s2
ess|Static error Þ Lt
s®0 K
1+ P
Cs
A A
Þ Lt = = C .
s®0 s . KP KP
s+
Cs
Thus, the steady state error also called as offset of proportional controller
is inversely proportional to KP.
An increase in KP will improve the speed of response as well as reduces the
error.
But a very large increase in proportional gain KP implies less
proportional band. This converts the operation of proportional mode into an
ON/OFF (two position mode) action.
To summarize, the proportional gain KP, no doubt improves both transient
and steady state responses but at the same time its value should not be large
because the proportional band may become narrow resulting in ON/OFF
characteristics.
Set point
qi
Inlet
LC Level controller
LT = level transmitter
or LT
level transducer h (Control valve)
q0
Outlet
FIG. 3.9
The process equation/transfer function from chapter (1)
Qi(s) R
H (s )
T = RC
Ts + 1
R = Control valve characteristics
C = Area/capacity of tank (m2).
The control loop with level controller (proportional mode) will be
Qi(s) R H(s)
KP
– Ts + 1
RK P
H ( s) 1 + Ts RK P
= =
Q i ( s) RK P Ts + ( RK P + 1)
1+
1 + Ts
RK P
H ( s) T
= . ...(17)
Q i ( s) (1 + RK P )
S+
T
The transient response is obtained by applying inverse Laplace transform
to eqn. 17.
æ RK P + 1ö
RK P e èç T ÷ø t
h(t) = . .
T
Comparing
æ RK P + 1ö
ç
è T ÷ø t
e » et/T
5TH PROOF 4-9-09
1 RK P + 1 RK P + 1
= =
T T RC
RC
T= .
RK P + 1
4 RC
Stabilization Time [ts] = 4T = . ...(18)
RK P + 1
S. R( s)
ess|Static error = Lt
s ® 0 1 + GC ( s)GP ( s)
R
where, Gc (s) = KP ; Gp (s) = and R(s) = Qi(s).
1 + Ts
Let there be sudden change in inlet flow rate (qi) causing the level to
change from the set point value, i.e., Qi(s) = A/s (step change).
A
s.
s
ess|Static error = Lt .
s®0 RK P
1+
1 + TS
A
ess|Static error = . ...(19)
1 + RK P
qi
m[Load variable]
LC T = Time constant = RC
LT
h 'c' 'R'
q0
FIG. 3.10
U(s) R
Ts + 1
+ H(s)
+
Qi(s) R
Ts + 1
qi
h2 c2 R2
q0
Tank (2)
FIG. 3.11
5TH PROOF 4-9-09
dh1 h1
qi = c1 + .
dt R1
dh1
R1qi = R1C1 + h1 .
dt
dh1
R1qi = T1 + h1 . ...(20)
dt
For tank (2)
dh2
qi + u = c2 + qo .
dt
q1 + u = c2
dh2 h2
+
hLM
qo = 2 .
OP
dt R2 R2 N Q
dh2
R2[q1 + u] = R2c2 + h2 .
dt
h1
Since q1 =
R1
R2
LM h
1 OP
+ u = R2 c2
dh2
+ h2 .
NR1 Q dt
R2 dh
h1 + R2u = T 2 2 + h2. ...(21)
R1 dt
Applying L.T to eqn. 20 and 21
R1Qi(s) = [T1s + 1]H1(s)
R1
H1(s) = .Qi (s).
T1 s + 1
R2
.H1(s) + R2u(s) = [T 2s + 1]H2(s)
R1
5TH PROOF 4-9-09
u(s) R2
T2s + 1
FIG. 3.12
H2(s) = o/p and Qi(s) and u(s) are input and load variable, respectively.
Since, the overall output of Fig. 3.11 is considered as level in tank (2) i.e.,
h2 a control loop may be formulated either by considering the manipulating
variables to be
qi = inlet flow rate to tank (1)
u = load variable to tank (2)
q0 = outlet flow rate of tank (2) and/or.
q1 = outlet flow rate of tank (1).
The various control schemes are shown in the following figures, 3.13 to
3.15.
qi
h1 R1
q1 Set point
LC
LT
h2 R2
q0
FIG. 3.13
qi
h1 c1 R1
u LC Set point
q1
LT R2
h2 c2
q0
FIG. 3.14
5TH PROOF 4-9-09
qi
h1 c1 R1
u LC Set point
q1
LT
h2 R2
q0
FIG. 3.15
u(s) R2
T2s + 1
Qi(s) R1 R2/R1 +
KP + H (s )
– T1s + 1 T 2s + 1
Controller Process
FIG. 3.16
Let u(s) = 0.
Qi(s) KPR2
H(s)
– (T1s + 1)(T2s + 1)
K P R2
H ( s) ( 1 1)(T2 s + 1)
T s + K P R2
= =
Q i ( s) K P R2 (T1 s + 1)(T2 s + 1) + K P R2
1+
(T1 s + 1)(T2 s + 1)
H ( s) K P R2
= 2
Q i ( s) T1T2 s + s[T1 + T2 ] + ( K P R2 + 1)
5TH PROOF 4-9-09
K P R2
H ( s) T1T2
Þ ...(22)
Q i ( s) 2 s[T1 + T2 ] K P R2 + 1
s + +
T1T2 T1T2
Comparing the characteristic equation of 22 with std. characteristic
equation of second order system s2 + 2x w ns + w n2 = 0,
T1 + T2 K P R2 + 1
2xwn = ; wn 2 =
T1T2 T1T2
K P R2 + 1
wn = (r/s).
T1T2
K P R2 + 1 T + T2
2x = 1
T1T2 T1T2
F R K + 1I = F T + T I 2
4x 2 GH T T JK GH T T JK
2 P
1 2
1
1 2
2
T1T2 (T1 + T2 ) 2 1
x2 = . .
R2 K P + 1 (T1T2 ) 2 4
(T1 + T2 )2
x= . ...(23)
4T1T2 ( R2 K P + 1)
As seen from the above analysis, the response of the two time constant
process with proportional controller depend on two factors.
(1) Characteristic time, defined as reciprocal of undamped natural
frequency (wn).
1 T1T2
T= = .
wn K P R2 + 1
(T1 + T2 )2
(x) = .
4T1T2 ( R2 K P + 1)
The proper choice of proportional gain (KP) affects the damping ratio of
the system which may give rise to different types of responses such as over
damped, under damped, and critically damped responses.
5TH PROOF 4-9-09
1
Since, x a , for smaller values of K P the response may be faster
KP
(x = 1) or sluggish (x > 1). As KP is increased, the response exhibits damped
oscillations (x < 1).
Hence, a compromise value of KP such that the damping ratio is 0.3 is
generally considered to be satisfactory.
From eqn. 23,
4x 2 ( R2 K P + 1) (T1 + T2 ) 2
=
T1T2 (T1T2 ) 2
(T1 + T2 ) 2
R2KP + 1 =
T1T2 . x 2 .4
1 é T1 T2 ù
R2KP + 1 = 2 ê
+ + 2ú
4x ë T2 T1 û
ìï 1 é T1 T2 ù üï
Þ R2KP = í 2 ê + + 2ú ý 1. ...(24)
îï 4x ë T2 T1 û ïþ
The optimum value of proportional gain may be calculated from eqn. 24
for this particular two time constant process.
Error Computation
Error may be computed by considering sudden changes either in the inlet flow
rate qi or the load variable u.
In Fig. 3.26, let there be a sudden change in inlet flow rate Qi(s)
S . R( s)
ess|Static error = Lt
s®0 1 + Gc ( s)GP ( s)
(offset of proportionnal controller)
S . Q i ( s)
ess|Static error = Lt
s ®0 K P . R2
1+
(T1 s + 1)(T2 s + 1)
1
let Qi(s) = (Sudden change)
S.
5TH PROOF 4-9-09
1s.
s
ess|Static error = Lt
s ®0 K P R2
1+
(T1s + 1) (T2 s + 1)
1
ess|Static error = . ...(25)
1 + K P R2
R2 KP = 0
Error
esc
R2 KP = 1
R2 KP = 10
offset
t
FIG. 3.17
h "c" R u
FIG. 3.18
5TH PROOF 4-9-09
V(s) 1
1 + Ts
Qi(s) 1 R + H(s)
Tis 1 + Ts +
–
Controller Process
FIG. 3.19
Qi(s) 1 R
1 + Ts H(s)
– Tis
R
H( s ) Ti s(1 + Ts) R
= Þ
Q i ( s) R 2
TiTs + Ti s + R
1+
Ti s(1 + Ts)
R
H( s ) i .
TT
= . ...(26)
Q i ( s) 1 R
s2 + . s2 +
T TiT .
Comparing the characteristic equation of eqn. 3.26 with standard second
order characteristic equation i.e., s2 + 2x wns + w n2.
1 R
s2 + .s + » s 2 + 2 x wn s + wn 2
T TiT
R
wn = r/s.
TiT .
1 R 1
2x w n = Þ 2x = ...(27)
T TiT T.
4x 2 R 1
Þ = Þ 4x 2RT = Ti
Ti T T 2
5TH PROOF 4-9-09
1
Ti = 4 ´ ´R´T
(3)2
4 RT
Ti = . ...(28)
9
This setting of Ti results in minimum deviation of the response from the
set point or input value and also avoids encountering excessive oscillations in
the response.
From the error terminology, the integral controller effect can be analyzed
in terms of state error and offset (error due to load variables). The expression
for state error is
S . R( s)
ess|Static error = Lt .
s ® 0 1 + Gc ( s)GP ( s)
S . u( s). N ( s)
And ess|offset = Lt .
s®0 1 + Gc ( s)GP ( s)
S . Q i ( s)
ess|Static error = Lt .
s®0 1 + Gc ( s)GP ( s)
1
Qi(s) = [unit step input].
s
1
s.
s
ess|static error = Lt .
s®0 R
1+
Ti s(1 + Ts)
R
where, Gc(s)GP(s) =
Ti s(1 + Ts)
5TH PROOF 4-9-09
1
ess|static error = Lt
s®0 R
1+
Ti s(1 + Ts)
1
ess|Static error = = 0.
1+ ¥
Similary, let there be a sudden change in load variable u(s).
1
u(s) = [Unit step change].
s
S . u( s). N ( s)
ess|offset = Lt .
s®0 1 + Gc ( s)GP ( s)
From the block diagram representation shown in Fig. 3.19.
1 R
N(s) = GP(s) =
1 + Ts 1 + Ts
1
Gc(s) = .
Ti s .
1 1
- s..
s 1 + Ts
ess|offset = Lt
s ®0 R
1+
Ti s(1 + Ts)
1
ess|offset = = 0.
1+ ¥
The integral controller is very effective in improving the steady state
response characteristics as it completely eliminates the error between input or
set point and controlled variable whether due to variations in the input or set
point or due to load variable changes.
Eqn. 28 may be utilized for obtaining properly tuned value of integral
time Ti depending on the dynamics and characteristics of process with which
it is employed.
V(s) R2
T2s + 1
Qi(s) 1 R1 R2/R1 +
+ H (s )
– Tis T1s + 1 T 2s + 1
Controller Process
FIG. 3.20
Qi(s) 1 R2
Tis H(s)
– (T1s + 1)(T2s + 1)
R2
(Ti s)(T1 s + 1)(T2 s + 1)
H( s ) R2
= 1+ 2
Q i ( s) (Ti s)[T1T2 s + s(T1 + T2 ) + 1]
H ( s) R2
= 2
Q i ( s) Ti s[T1T2 s + s[T1 + T2 ] + 1] + R2
H( s ) R2
= 3 2
...(29)
Q i ( s) TiT1T2 s + s Ti (T1 + T2 ) + Ti s + R2
with integral controller, the characteristic equation of the overall system
function has increased to third order polynomial.
Any higher order systems nature of response is obtained by
approximating it to second order system with respect to dominant pole
analysis. Even in this case also, the integral time T i is a function of damping
ratio (y), which is now called as relative damping factor or damping ratio.
The effect of integral action, whether single time constant process or
multiple time constant process is that it induces oscillations in the nature of
response of the system. The number of oscillations and their peak values can
be monitored by appropriately choosing the value of damping factor, which in
turn depends on the optimum setting of integral or reset time Ti. However,
as the order of system is increased with this mode of action, there is definitely
a possibility for the response to exhibit sluggishness. The response of the
system with integral controller for different values of integral or reset time Ti
settings are shown below.
5TH PROOF 4-9-09
Case 1: Case 2:
Ti = 0.75 Ti = 0.25
h (t)
o/p
MP1 MP2
h (t)
o/p
Steady
state
input
tr1 t t
t= o t = o tr2
ts1 ts2
FIG. 3.21
As seen from the above two cases as the value of Ti is reduced, the rise time
decreases and, the number of damped oscillations in the response increases with
maximum peak overshoot shooting up resulting in more settling time.
Note: MP = Maximum peak overshoot
t s = Settling time.
tr = Rise time.
S . Q i ( s)
ess|Static error = Lt
s®0 1 + GC ( s)GP ( s)
1
s.
s
= Lt
s ®0 1 R2
1+ ´
Ti s (T1s + 1) (T2 s + 1)
1
ess|Static error = = 0.
1+ ¥
S . u( s). N ( s)
ess|offset = Lt
s®0 1 + GC ( s)GP ( s)
1 R2
.
-s.
s T2 s + 1
= Lt
s ®0 1 R2
1+ ´
Ti s (T1s + 1) (T2 s + 1)
5TH PROOF 4-9-09
R2
ess|offset = = 0.
1+ ¥
Thus, both the error i.e., due to changes in input and load variables are
successfully eliminated by integral action, even for two time constant process.
E(s) T dS P(s)
Let us consider the effect of this mode on single time constant process as
shown in Fig. 3.22.
qi
LC
h "c" 'R' u
FIG. 3.22
The block diagram showing the controller and process is presented below.
V(s) 1
1 + Ts
Qi(s) R + H(s)
Tds +
– 1 + Ts
Controller Process
FIG. 3.23
Qi(s) R
Td s H(s)
– 1 + Ts
RTd s
H ( s) 1 + Ts RTd s
= =
Q i ( s) RTd s 1 + Ts + RTd s
1+
1 + Ts
H ( s) RTd s
=
Q i ( s) s(T + Td R) + 1
RTd s
H ( s) T + Td R
= . ...(30)
Q i ( s) [ s + 1/ T + Td R ]
Consider a sudden change in inlet flow rate qi (step change).
1
Let Qi(s) =
s
RTd s
T + Td R 1
H(s) = ´ .
1 s
s+
(T + Td R )
RTd
T + Td R
\ H(s) = .
1
S+
(T + Td R)
Applying final value theorem
Lt h(t) = Lt S . H ( s)
t®¥ s ®0
RTd
S.
T + Td R
\ Lt = 0.
s®0 1
S+
(T + Td R)
Þ h(t) = 0
From the above analysis it is evident that for step change in deviation, the
derivative controller does not respond i.e., its output to the final control
element is zero and hence there is no effect on the output which is level of the
liquid in the tank [h(t) = 0].
5TH PROOF 4-9-09
F RT s I
H( s)
GH T + T . RJK
d
d
= .
Q i ( s) 1
S+
(T + Td R)
1
Considering Qi =
s2
RTd s
T + Td R 1
H(s) = . .
1 s 2
s +
T + Td R
Again applying final value theorem
RTd
T + Td R
Lt S.H(s) = Lt s . .
s®0 s®0 s
s2 +
T + Td R
RTd
T + Td R.
Þ Lt = RT d.
s®0 1
S+
T + Td R.
This implies that whenever there is a continuous constant change in inlet,
disturbance i.e., load variable or even the set point causing the output i.e.,
liquid level in this example, to also change at a particular rate, the derivative
controller sends a signal to the final control element initiating a sudden action
before the error begins to occur as is illustrated graphically in figure 3.24.
Error
Due to
changes in
inlet flow
rate t
t = t1
Derivative
controller
o/p
} RTd
t
FIG. 3.24
5TH PROOF 4-9-09
This type of response results in large errors rather than reducing it to zero
and bringing the controlled variable to the required set point value.
Evaluating error
S . Q i ( s)
e ss|Inlet flow changes = Lt
s ® 0 1 + GC ( s)GP ( s)
1
s.
s2
Þ Lt
s ®0 T s. R
1+ d
1 + Ts
1
= Lt = ¥.
s ®0 s2 RTd
S+
1 + Ts
Infinity error implies, very large error may result due to derivative action.
P ( s) é 1 ù
= Kp ê1 +
E( s) ë Ti s úû
1
1 + Ts
Qi(s) R H(s)
KP 1 + 1 +
+
– Tis 1 + Ts
Controller Process
FIG. 3.25
let u(s) = 0, since analysis is done by considering only one input and output at
active for a transfer function
æ H ( s) ö
The overall transfer function ç can be obtained by considering the
è Q i ( s) ÷ø
block diagram shown in Fig. (3.25).
Qi(s) R H(s)
KP 1 + 1
– Ts i 1 + Ts
é K pTi s + K p ù é R ù
ê úê ú
æ H ( s) ö ë Ti s û ë1 + Ts û
\ çè Q ( s) ÷ø =
i
( K pTi s + K p )( R)
1+
Ti s(1 + Ts)
RK pTi s + RK p
Þ
Ti s(1 + Ts) + RK pTi s + RK p
5TH PROOF 4-9-09
RK p [Ti s + 1]
Þ 2
TTs
i + Ti s + RK pTi s + RK p
RK p (1 + Ti s)
TT
i
Þ
s( Ti + RK p Ti ) RK p
s2 + +
Ti T TiT
RK p (1 + Ti s)
H ( s) TiT
Þ s(1 + RK p ) RK p ...(30a)
Q i ( s) 2
s + +
T Ti T
RK p RK p
w2n = Þ wn = ...(32)
TiT TiT
1 TiT
\ Tc = = ...(33)
wn RK p
1 + RK p
Now, 2xwn =
T
RK p 1 + RK p
2x =
TiT T
5TH PROOF 4-9-09
4 x2 RK p (1 + RK p )2
=
Ti T T2
Ti (1 + RK p )2
4x2 =
RK p
Ti (1 + RK p )2
x2 =
4 RK p
Ti (1 + RK p )
Þ x= ...(34)
4( RK p )
é 1 ù
Gc(s) = Kp ê1 + ú
ë Tis û
R
GP(s) =
1 + Ts
Assuming step change in input Q i(s).
1
s.
s
ess|static error = Lt
s®0 æ 1 öæ R ö
1 + K p ç1 +
è Ti s ø÷ çè 1 + Ts ÷ø
1
= =0
1 + K p (1 + ¥)( R)
é - s. u( s). N ( s) ù
ess|offset = Lt ê
s® 0 1 + Gc ( s)GP ( s) ú
ë û
5TH PROOF 4-9-09
1
N(s) =
1 + Ts
æ 1 ö
Gc(s) = Kp ç1 +
è Ti s ø÷
R
Gp(s) =
1 + Ts
1 1
- s.
.
s 1 + Ts
ess|offset = Lt
s®0 æ 1 öæ R ö
1 + K P ç1 +
è Ti s ÷ø èç 1 + Ts ø÷
-1
ess|offset = =0
1 + K P (1 + ¥)( R)
Hence it may be concluded that in both the cases the error is eliminated
completely. Thus a P + I controller is effective in reducing the error to a
minimum valve let us see the effect of proportional gain and integral time on
the response. For a constant value of integral/Reset time Ti and changing the
value of proportional gain, for higher values of proportional gain the response
exhibits more number of damped oscillations as a result more peak overshoots
which clearly indicates that the damping factor value is between 0 and 1.
Ti = Constant
Curve - 1
Error e
Curve - 2
Curve - 3
O
Curve - 4
Time (t)
FIG. 3.26
5TH PROOF 4-9-09
Curve1 Indicates the response for a valve of Kp between 1 and 10. The
response does not exhibit any oscillation and is a smooth curve decaying to
zero error position.
Curve2 Is for the valves of proportional gain Kp chosen between 10 and
20. Very small oscillations can be observed.
Curve3 Is the response characteristics for 50 > Kp > 20. The response for
this setting exhibits oscillations in response before reaching the steady state
value.
For higher values of proportional sensitivity there are considerable
overshoots and undershoots in the response before decaying to zero error
position or steady state valve of the input or set point.
Similarly, the effect of reset time on the characteristics of response also can
be analyzed by keeping the proportional gain constant. Figure 3.26 indicates
the nature of response for different valves of Ti.
Error e
Curve 1
Curve 2
Curve 3
O
Curve 4
Time (t)
FIG. 3.27
For very larger values of Ti (Ti = ¥), the number of oscillations in the
response are very less and the error is also considerable between input and
output as shown in Curve1.
The response shown in curve (2) is for valves of Ti between 10 and 15.
However when Ti valve is chosen between 0 and 1 [0 < Ti < 1].
Controller
Process
Qi(s) R1 R2/R1 H(s)
KP[1 + Td s]
– T1s + 1 T2s + 1
FIG. 3.28
R2 K P (1 + Td s)
H ( s) ( 1s + 1)(T2 s + 1)
T
=
Q i ( s) R K (1 + Td s)
1+ 2 P
(T1 s + 1)(T2 s + 1)
H ( s) R2 K P (1 + Td s)
\ = ...(35)
Q i ( s) (T1 s + 1)(T2 s + 1) + R2 K P (1 + Td s)
s[T1 + T2 + R2 K P Td ] 1 + R2 K P
Þ s2 + + =0
T1T2 T1T2
1 + R2 K P
wn =
T1T2
T1 + T2 + R2 K P Td
2xwn =
T1T2
1 + R2 K P T + T2 + R2 K P Td
2x = 1
T1T2 T1T2
4x2 (1 + R2 K P ) (T + T2 + R2 K P Td )
= 1
T1T2 (T1T2 )2
(T1 + T2 + R2 K P Td )2
x2 =
4(1 + R2 K P )T1T2
5TH PROOF 4-9-09
(T1 + T2 + R2 K PTd )2
x= ...(36)
4T1T2 (1 + R2 K P )
(2 + K P Td )2
Þ x= [Assuming T1 = T2 = R2 = 1]
4(1 + K P )
(2 + Td )2
x= ...(37)
8
Td Damping (y)
0 0.707
0.1 0.74
. .
. .
0.5 0.88
. .
. .
1.0 1.16
The table above indicates that an increase in derivate time Td increases
the damping ratio ( x), thereby reducing the oscillations in the response.
Hence the addition of a properly tuned derivative action makes the
response of the system faster, in approaching the steady state value of input
Let us see its effect on error
s . R( s)
ess = Lt
s®0 1 + Gc ( s).GP ( s)
1
Assuming R(s) = (Linearly changing input & hence error)
s2
1
s.
s2 1
ess = Lt = =¥
s®0 K P (1 + Td s). R2 0
1+
(T1 s + 1)(T2 s + 1)
The steady state error and/or offset remains when P + D action is used
with processes for obtaining stable operation. However a large proportional
gain may reduce or minimise offset, but practically Large proportional gain
should not be employed because of the reasons already studied in previous
chapters. Therefore it is clear that this system cannot eliminate offset of
proportional controller or any type of steady state error in a process. It is
capable of handling fast process load changes as long as the load change offset
error is acceptable. The transient response for slow processes is also
considerably improved.
5TH PROOF 4-9-09
ì é 1 ùü
P(s) = í K P ê1 + + Td sú ý E(s).
î ë Ti s ûþ
Consider a single time constant process controlled by PID action
Qi(s) H(s)
1 R
KP(1 + + T d s)
– Tis 1 + Ts
FIG. 3.29
1
= Lt
s®0 é K ùé R ù
s + ê sK P + P + K PTd s2 ú ê ú
ë Ti û ë1 + Ts û
Ti
ess Þ (Since Ti << KP)
RK P
ess » 0.
Thus PID action successfully eliminates the error between input and
output with regard to speed of response Fig. 3.30 indicates the different
characteristics for different values of KP, Ti, Td.
5TH PROOF 4-9-09
KP = 5; Ti = 3 min;
Td = 0 min.
KP = 5; Ti = 3 min;
Td = 0.8 min.
FIG. 3.30
From the above two curves, with addition of derivative gain the systems
stabilization time is considerable reduced thus indicating faster response. The
above curves are based on experimental trials done for processes having single
or multiple time constant characteristics as mathematical analysis of such
systems with PID control is very complicated as far as calculation of responses
for different values of KP, Td and Ti are concerned.
Summary
The choice of a particular mode of controller depends on the characteristics of
the process, i.e. whether single capacitance type, single or multiple time
constant type etc.
For processes with more number of storage elements i.e. more time
constants the choice of composite controller modes is best suited. We have
studied in this chapter, the effects of controller modes on different types of
processes. Depending on the mode of controller operation particular process
was chosen for better understanding of the effect of controller mode.
The overall analysis may be summarized as follows.
(1) Proportional controller produces a permanent residual error known as
offset.
(2) Integral controller successfully eliminates offset or error but its overall
speed of response is much slower resulting in large stabilization time.
(3) Proportional + integral controller has no offset but the unstabilizing
influence of integral response is reflected in large stabilization time.
(4) Proportional + Derivative controller cannot eliminate offset or error
how ever due to the presence of derivative action the stabilization time
is considerably reduced thus making the speed of response much faster.
The offset of proportional controller can be made only half of what
would have been without derivative action.
5TH PROOF 4-9-09
SOLVED EXAMPLES
1. For the process shown in Fig. Find
(a) offset
(b) velocity error
1 10 C(s)
R(s) s (s + 1)
+
– +
U(s) 5
(s + 1)
1
s.
s2 1
e ss|velocity error = Lt Þ Lt
s®0 10 s®0 s .10
1+ s+
s( s + 1) s ( s + 1)
e 1
ss|velocity error =
= 0.1 units.
10
2. A two vessel process shown in figure has T1 = T2 = 60 secs; R =
1
sec/m2. And another design with T1 = 30 sec s; T2 = 120 sec s; R =
6
1
sec/m2. Which design provides less offset for setting of proportional
6
gain at 10?
5TH PROOF 4-9-09
qi
R1
(C1) q1
h1 u(s) [Disturbance]
Tank 1
R2
(C2)
h2
q0
Tank 2
dh1
\ R1qi = R1C1 + h1
dt
R1Q1(s) = [T1s + 1]H 1(s) [T 1 = R1C1] ...(1)
For tank (2)
q1 q1 + u = c2 dh2
dt
q1 = c2 dh2 + q0 u.
dt
h1 dh
= c2 2 + h2 R2u.
R1 dt
R2 dh
. h = h1 = R2C2 2 + h2 R2u.
R1 1 dt
R2
H (s) = [T2s + 1]H2(s) R2u(s). [T2 = R2C2]
R1 1
R2
H (s) + R2u(s) = [T 2s + 1]H2(s)
R1 1
R2
R1 R2
H2(s) = . H1(s) + u(s)
(T2 s + 1) (T2 s + 1)
5TH PROOF 4-9-09
R2
u(s) (T2s + 1)
Process
R2 R2
GP(s) = ; N(s) =
(T1 s + 1)(T2 s + 1) (T2 s + 1)
Case 1: T1 = T 2 = 60 sec s
1
R= sec/m2
6
s . u( s). N ( s)
ess|offset = Lt
s®0 1 + Gc ( s). GP ( s)
1 1/6
. s.
s (60 s + 1)
ess|offset = Lt 1/6
s®0
1 + 10.
(60 s + 1)2
-1/6 - 0.166
ess|offset = =
10 2.66
1+
6
ess|offset = 0.062 units
1 1/6
. s.
s (120 s + 1)
ess|offset = Lt
s®0 1/6
1 + 10.
(30 s + 1)(120 s + 1)
5TH PROOF 4-9-09
-1/6
ess|offset = = 0.062 units.
1 + 10/6
qi
Set point
LC
h
qo u t
Design specifications
LC = level controller (proportional mode)
KP = 5
A = Area of the tank
= 10 m2
qout = Load variable
= 0.1 step change
qi = Inlet flow rate
= 1 m3/s
Solution: to obtain the Process equation
dh
qi qout = c [c = capacity of tank].
dt
Applying laplace transform
Q i(s) Q out(s) = CS. H(S)
1 1
H(s) = . Q i(s) . Q out ( s)
cs cs
It is a single capacitance process.
5TH PROOF 4-9-09
u(s) 1
cs
Qi(s) 1 – H (s )
KP +
cs
0.1
u(s) = Q out(s) =
s
0.1 1
s. .
s 10 s
ess /offset = Lt
s®0 5
1+
10 s
0.1 0.1
= Lt Þ .
s®0 10 s + 5 5
Qi(s) 1 H(s)
KP
– cs
H ( s) K P / cs KP KP
= = =
Q i ( s) 1 + K P / cs cs + K P ( s + K P / c)
Applying inverse L.T.
K P (KP/c)t
h(t) = e .
c
e[KP /c]t » e t/T
5TH PROOF 4-9-09
1 K
= P Þ Time Constant T = c/KP
T c
4 c 4 ´ 10
Stabilization time = 4T = =
KP 5
= 8 secs
4. Determine static error and velocity error for the block diagram shown in
figure.
u
N
R G1 G2 +
+ C
–
1 R 1
G1(s) = ; G2(s) = ; N(s) =
Ti s Ts + 1 (Ts + 1)
s . R( s)
Solution: ess|static error Þ Lt
s®0 1 + Gc ( s) GP ( s)
Gc(s) = G1(s) = R/(Ts + 1)
\ R(s) = 1/s (unit step). (Assume)
1
s.
s 1
ess|static error = Lt = =0
s®0 1 R ¥
1+ .
Ti s (Ts + 1)
1
s.
s2
ess|velocity error = Lt
s®0 1 R
1
Ti s 1 + Ts
(R(s) = 1/s2 for velocity error).
1
ess/velocity errors = Lt
s®0 s R
s+ .
Ti s 1 + Ts
Ti
ess/velocity error = .
R
5TH PROOF 4-9-09
4 RC
5. For single time constant process prove that stabilization time is ?
1 + RK P
Solution: (Refer text)
6. For the process shown in figure obtain the setting of integral time for
damping ratio 1/3rd?
Design specification
qi
R = 5 s/m2
sp
C = 10 m2
LC
h
q0
dh
Rqi = Rc +h
dt
RQ i(s) = [Ts + 1] + H(s)
Gc(s) Gp(s)
Qi(s) 1 R H(s)
– Ti s Ts + 1
H ( s)
Solving for
Q i ( s)
5TH PROOF 4-9-09
R
Ti s(Ts + 1) R
Þ
R TT s2
+ Ti s + R
1+ i
Ti s(Ts + 1)
H ( s) R /TTi
\ =
Q i ( s) 1 R
s2 + s +
T TTi
Comparing the character equation with
s2 + 2xw ns + w2n = 0
R
wn =
TTi
1
2xwn =
T
R 1
2x =
TTi T
4x2 R 1
=
T Ti T2
4x2 R 1
= .
Ti T
Ti = 4x2 R.T.
Given in design specifications
x = 1/3; R = 5; T = RC = 50 sec s.
1
Ti = 4 ´ ´ 5 ´ 50 = 111.11 sec s.
9
Ti = 111.11 sec s. (or)
Ti = 1.85 mins
Solution: The block diagram with controller for single capacitance process is
(Controller) (Process)
Q(s) 1 1 H(s)
KP 1 + +
– Tis cs –
1
u(s) cs
Q(s) 1 H(s)
KP 1 + 1
– Tis cs
é æ 1 öù 1
ê K P ç1 + T s ÷ ú cs
H ( s) è i øû
= ë
Q ( s) æ 1 ö1
1 + K P ç1 +
è Ti s ø÷ cs
KP K
1+ + P2 = 0.
cs Ti cs
Tics2 + KPTis + KP = 0
K P Ti s KP
s2 + + =0
Ti c Ti c
KP K
s2 + .s + P = 0
c Ti c
KP KP
wn = and 2xwn =
Ti c c
KP K
2x = P
Ti ´ c c
5TH PROOF 4-9-09
4x2 K P K P2
= c2
Ti c
4x 2 K
= P
Ti c
Ti K P TK
4x2 = Þ x2 = i P
c 4c
Ti K P
x= or 4cx2 = TiKP
4c
qi u
R1
(C1) q1
h1
R2
C2
h2
q0
dh1 h
qi + u = c 1 + 1
dt R1
dh1
R1qi + R1u = R1c1 + h1.
dt
5TH PROOF 4-9-09
h1 dh h
= c2 2 + 2
R1 dt R2
R2 dh
h = R2c2 2 + h2
R1 1 dt
R2
H (s) = [T 2 S + 1] [T2 = R2C2] ...(2)
R1 1
Substituting (1) in (2).
R2 ì é R1 ù é R1 ù ü
R1 í ê T s + 1 ú Q i ( s) + ê (T s + 1) ú 4( s)ý = (T2 s + 1)H2(s)
ïî ë 1 û ë 1 û ïþ
é R2 ù é R2 ù
\ H2(s) = ê
(T s + 1)(T s + 1) ú Q i ( s) + ê (T s + 1)(T s + 1) ú u(s)
ë 1 2 û ë 1 2 û
Qi(s) R2
(T1s + 1)(T2s + 1)
+
+ H2(s)
V(s) R2
(T1s + 1)(T2s + 1)
u(s) R2
(T1s + 1)(T2s + 1)
Qi(s) R2 +
+ H2(s)
– (T1s + 1)(T2s + 1)
5TH PROOF 4-9-09
s . R( s)
® ess|static error = Lt
s®0 1 + Gc ( s)GP ( s)
1
let R(s) = (unit step)
s
R2
Gc(s) = KP; GP(s) =
(T1 s + 1)(T2 s + 1)
1
s.
s
ess|static error = Lt
s®0 K P . R2
1+
(T1 s + 1)(T2 s + 1)
1
ess|static error =
1 + K P . R2
s . u( s) N ( s)
® ess|offset = Lt
s®0 1 + Gc ( s). GP ( s)
1
let u(s) = (unit step)
s
R2
N(s) = .
(T1 s + 1)(T2 s + 1)
1 R2
. s.
s (T1s + 1)(T2 s + 1)
ess|offset = Lt
s®0 K P R2
1+
(T1 s + 1)(T2 s + 1)
- R2
ess|offset =
1 + K P R2
9. For the process shown in 6 prob (8) show that addition of integral
controller eliminate offset.
æ 1 ö
Solution: Gc(s) = KP ç1 + ÷ Þ P + I controller
è Tisø
R2
GP(s) =
(T1 s + 1)(T2 s + 1)
1 R2
- s. .
s (T1 s + 1)(T2 s + 1)
\ ess|offset =
æ 1 ö R2
1 + K P ç1 + ÷
è Ti s ø (T1s + 1)(T2 s + 1)
R2
ess|offset = =0
1 + K P (1 + ¥). R2
10. For a closed loop automatic control from the operational diagram
R1 R1
G1 = Kc; G2 = and N= .
T1 s(T2 s + 1) T1 s(T2 s + 1)
(a) Calculate proportional sensitivity for critical damping?
(b) Calculate static error and offset.
Solution: Given
G1 = Gc(s) = Kc = proportional gain
R1
G2 = GP(s) =
T1 s(T2 s + 1)
R1
N = N(s) = .
T1 s(T2 s + 1)
The block diagram may be formulated as.
u(s) R1
T1s(T2s + 1)
Qi(s) R1 +
KC + H(s)
– (T1s(T2s + 1)
let u(s) = 0.
R1K P
H ( s) T1s(T2 s + 1)
=
Q i ( s) R1K P
1+
T1 s(T2 s + 1)
5TH PROOF 4-9-09
R1K c
1+ =0
T1 s(T2 s + 1)
1 RK
s2 + .s + 1 c = 0
T2 T1T2
R1K c
wn = (r/s)
T1T2
R1K c 1
2x ´ =
T1T2 T2
4x2 R1 K c 1
= 2
T1 T2 T2
T1
4x2R1Kc =
T2
T1
Kc =
4 R1T2
s . R( s)
(b) ess|static error = Lt
s®0 1 + Gc ( s) GP ( s)
1
s.
s 1
it Lt = =0
s®0 K c R1 1 + ¥
1+
T1 s(1 + T2 s)
5TH PROOF 4-9-09
- s . u( s). N ( s)
® ess|offset = Lt
s ®0 1 + Gc ( s) GP ( s)
1 R1
- s. .
s T1 s(1 + T2 s)
= Lt
s ®0 K c R1
1+
T1 s(1 + T2 s)
- R1
=
T1 s(1 + T2 s) + K c R1
-1
ess|offset =
Kc
R1
N(s) =
T1 s(T2 s + 1)
(a) Calculate KP for damping ratio 1/3rd?
(b) Calculate the period of oscillations?
Solution: Constructing the characteristic equation
1 + Gc(s) GP(s) = 0
K P (Td s + 1). R1
1+ =0
T1 s(T2 s + 1)
s2 +
[T1 + K P R1Td ] + K P R1
=0
T1T2 T1T2
K P R1
wn = r/s.
T1T2
5TH PROOF 4-9-09
K P R1
4p2 f 2n =
T1T2
K P R1
f 2n =
4p2T1T2
Period of oscillation (undamped oscillations)
1 1 4p2T1T2
T= ; T2 = =
fn f n2 K P R1
4 p2T1T2
T= sec s
K P R1
K P R1 T + K P R1Td
(a) 2x = 1
T1T2 T1T2
4x2 . K P R1 (T + K P R1Td )2
= 1
T1T2 (T1T2 ) 2
1
For x=
3
4 (T + K P Td )2
KPR1 = 1
9 T1T2
é 1 R1 ù
ê - s . s . T s(T s + 1) ú
1 2
Lt ê ú
s®0 ê K P (Td s + 1). R1 ú
1+
êë T1 s(T2 s + 1) úû
é ( - R1) ù
Lt ê ú
s®0 T s
ë 1 2(T s + 1) + K (
P dT s + 1). R1û
- R1 -1
ess|offset Þ =
K P R1 K P
5TH PROOF 4-9-09
13. Prove that for a single capacitance process, integral action is not required
to eliminate static error, when only P-controller is used?
Solution: let Gc(s) = KP = proportional controller
1
GP(s) = [single capacitance process]
cs
s . R( s)
ess|static error = Lt
s®0 1 + Gc ( s) GP ( s)
1
let R(s) =
s
1
s.
s 1
Lt = =0
s ®0 KP 1+ ¥
1+
cs
ess|static error = 0
1
GP(s) =
cs
Calculate the system time constant?
Solution: The block diagram representation is
Qi(s) 1 H(s)
KP(1 + Td s)
– cs
H ( s) K P (1 + Td s)/ cs K P (1 + Td s)
= Þ
Q i ( s) K P (1 + Td s) cs + K P (1 + Td s)
1+
cs
K P (1 + Td s) K P (1 + Td s)
Þ
cs + K P + K P Td s cs + K P Td s + K P
K P (1 + Td s)
H ( s) ( c + K P Td )
Þ
Q i ( s) æ KP ö
çè s + c + K T ÷ø
P d
5TH PROOF 4-9-09
K P (1 + Td s)
c + K P Td A
=
æ KP ö æ KP ö
çè s + c + K T ÷ø ç s + c + K T ø÷
è
P d P d
Put s = 0
KP
c + K P Td A
=
KP KP
c + K P Td c + K P Td
KP
\ A=
c + K P Td
æ KP ö
H ( s) çè c + K T ÷ø
P d
\ =
Q i ( s) ( s + K P / c + K P Td )
KP F I
h(t) =
c + K P Td
.e GHKP
c + K PTd
t JK
Comparing with e t/T
1 KP
=
T c + K P Td
c + K P Td
T=
KP
c
T= + Td
KP
15. Calculate the velocity error for single time constant process with P + I
controller?
Design specifications are
R = 5 s/m2; c = 5 m2; KP = 10
Ti = 2 mins.
s . R( s)
Solution: Ess|velocity error = Lt
s®0 1 + Gc ( s) GP ( s)
5TH PROOF 4-9-09
æ 1 ö é 1 ù
Gc(s) = KP ç1 + = 10 ê1 +
è ÷
Ti s ø ë 120s úû
R 5
GP(s) = = [T = Rc]
1 + Ts 1 + 25s
1
s.
ess|velocity error = Lt s2
s®0 æ 1 öæ 5 ö
1 + 10 ç1 + ÷
è 120 s ø çè 1 + 25s ÷ø
1
= Lt
s®0 æ 10 ö 5
s + ç10 s + ÷
è 120 ø 1 + 25s
1
=
é 10 ù
0 + ê0 + [5]
ë 120 úû
12 0
ess|velocity error = = 2.4 units
50
QUESTION SET
1. How are errors classified in process control? Obtain the expressions of
various error representations?
2. Discuss the effect of proportional controller on
(a) Single capacitance process
(b) Time constant process.
3. Why derivative action is not suitable in single mode configuration?
Explain?
4. Discuss the effect of PI action on time constant process?
5. Compare the effects of PI, PD and PID actions on time constant process?
"
+0)26-4
Tuning of Controllers
4.1 INTRODUCTION
The optimum values of controller parameters Kp, Ti and Td i.e. proportional
gain, integral time and derivative time are intimately related to the dynamics
and characteristics of the process, where a particular dynamic variable is
measured monitored and controlled.
Tuning may be defined as the optimal adjustment of controller
parameters to achieve satisfactory control.
The procedure of tuning may depend on elaborate optimization
calculation based on process dynamics and model to trial and error attempt in
achieving the required control objective. Some methods involve semi-
empirical equations, in that they depend on measurements made on the
control system with all its complexities to determine factors used in adjustment
formulas.
The frequency response analysis approach is more analytical and is based
on model or transfer function of the process and the control loop. From the
time response characteristics of second order system subjected to step change
in input or load or set point, the typical criteria for good control is that, there
should be minimum overshoot, rise time and settling time and one-quarter
decay ratio.
One criteria that is often used to evaluate a response of a control system is
the integral of the square of the error with respect to time (ISE)
z
¥
ISE = e2 dt
0
z
¥
Integral of the Absolute value of Error (IAE) = |e|dt
0
5TH PROOF 4-9-09
z
¥
Integral of time Weighted Absolute Error (ITAE) = |e|t × dt
0
o/p
I TAE
IAE
IsE
t
FIG. 4.1
The various tuning sequences may thus be summarised as follows.
Tuning Methods
LT
Outlet
FIG. 4.2
Level Transducer
Hm (s)
FIG. 4.3
For the measurement of process reaction curve, the process control loop is
opened as shown in Fig. 4.4, where the output of the controller is
disconnected from the control valve.
Hm (s)
To Recorder
for Measurement
of Process Reaction Curve
FIG. 4.4
A step change of magnitude A is introduced at the control valve and the
values of output as time t ® ¥ are recorded at the outup of the measuring
element as shown in Fig. 4.4.
As maintained earlier, the transfer function between control valve, process
and measuring element is approximated to be first order system with dead
time.
K e - td s
Gv(s) Gp(s) Hm(s) = ...(1)
1 + τs
o/p
B
Tangential
Line
Inflection
Point
t
td t
I/p
FIG. 4.5
obtained from the process reaction curve i.e. eqns. 2 to 4. They derived
expressions for the best controller settings using load changes and various
1
performance criteria such as decay ratio (Quarter amplitude criteria)
4
minimum offset and minimum ISE. Thus, when quarter amplitude criteria,
which is the most favourable performance specification by practising process
control technologists is desired, these c-c-corrections are an absolute necessity.
1
decay ratio implies the ratio between second peak overshoot and first
4
1
peak overshoot from the steady state value of the input should be , in the
4
step response of standard second order system.
y 1
x y x = 4
FIG. 4.6
y 1
A decay ratio = is a reasonable trade off between fast rise time and
x 4
reasonable settling time.
Kp =
A LM1 + S × t OP
d
...(5)
S × td N 3× B Q
(2) Proportional + Integral mode [P + I]
A LM0.9 + R OP ; LM 30 + 3 R OP t
Kp =
S × td N 12 Q Ti =
N9 + 20 R Q d ...(6)
5TH PROOF 4-9-09
Td =t M
L 4 OP ...(7)
d
N11 + 2 R Q
S × td
where R = log ratio (unit less) = .
B
Tuning procedure
(1) The integral and derivative actions are reduced to the minimum effect
after the process reaches steady state at the normal level of operation.
5TH PROOF 4-9-09
F Imaginary Part I
|F ( jw ) = Tan1 GH Real Part JK
Bode plots are asymptotic plots that give magnitude and phase
characteristics versus frequency (w) calibrated on log scale. To obtain
asymptotic characteristics of magnitude, its decibel values (db) are considered
where db value is defined as 20 log [F(jw)].
The stability criteria from bode plots is derived from the characteristic
equation 1 + G(s)H(s) = 0 expressed as
G(s)H(s) = 1 (or) G( jw) H( jw) = 1+ jo.
We define the following terms to determine the actual stability of the
system.
(1) Gain cross over frequency [MC?]
It is defined as the frequency at which the magnitude of open loop transfer
function G(jw) H( jw) is unity or 0 db.
1
G.M. =
X
1
G.M. [db] = 20 log
X
(4) Phase Margin [P ×M.]
It is the allowable phase lag i.e. the amount of additional phase lag that
cab be introduced into a system till the system reaches the verge of instability.
It is obtained at gain cross over frequency (wgc).
Let |G( jw ) H( jw )|w = w gc = f
P×M. = 180° + f.
0db
G.m. = o – (– x)
– Xdb = xdb
f
P.M. P.M. = f – (–180)
– 180°
= 180° + f = + VE
Stable system
G.M. = P.M. = +VE
FIG. 4.7
5TH PROOF 4-9-09
db 0 – (x) = – x db
0 db
– 180°
180° + f = – ve
f
log w
wpc wgc
Unstable system
G.M. = P.M. = – ve
FIG. 4.8
0 db
– 180°
log w
wgc = wpc
G.M. = P.M. = 0
Marginally stable system
FIG. 4.9
To find gain margin, extend the phase cross over frequency line till it
meets the magnitude characteristics. The distance between 0 db line to the
value in db where wpc line intersects the magnitude plot is the G.M.
Similarly to find phase margin the distance between the point where the
gain cross over frequency intersects the resultant phase angle plot to 180°
reference line gives P.M.
(2) From the Bode diagram find gain margin and phase cross over
frequency.
Let G.M. = y db
FG 1 IJ = y db
i.e. 20 log
H xK
1 y
log =
x 20
1
= log1 (y/20)
x
Thus is the value of critical (or) ultimate gain Kc
1
Therefore Kc = = log1 (y/20).
x
The critical period is,
2p
Tc = [sec/min]
w pc
where wpc = phase cross over frequency
w pc = 2p fc,
1 1
Tc = Þ fc = ,
fc Tc
2p 2p
\ w pc = Þ Tc =
Tc w pc
(3) Use Zeigler Nicol Settings to find Kp, Ti and Td values.
Rooth Array =
S4 a0 a2 a4
3
S a1 a3 0
2
S b1 b2 0
1
S c1 0 0
0
S d1 0 0
The elements of other rows and columns b1, b2 etc. are obtained as follows
a1a2 - a0 a3 a1a 4 - a0 ´ 0
b1 = b2 =
a1 a1
b1a3 - a1b2
c1 = c2 = 0
b1
c1b2 - b1 ´ 0
d1 =
c1
After obtaining all the value of other elements, the condition for stability is
there should not be any sign change in first column of Routh array i.e. all the
first column elements should be greater then zero.
Procedure
(1) Construct the characteristic equation
1 + Gc(s) Gv(s) Gp(s) Hm(s) = 0.
With Gc(s) = controller transfer function with only K p [proportional
gain] with derivative and integral actions reduced to minimum effect.
(2) Construct Routh array and evaluate all the elements of the array in
terms of Kp.
(3) Apply, the condition of stability, where all the first column elements
should be greater than zero.
(4) The value of proportional gain Kp which gives marginally stable case,
thus, corresponds to critical or ultimate gain Kc.
(5) For this value of proportional gain Kp construct the auxiliary equation
A(s).
(6) Compare the roots of the auxiliary equation A(s) with S = jw and find
w.
2π
Therefore w = wc [critical frequency], wc =
Tc
2π
Þ Tc = (sec/min)
wc
(7) Apply Ziegler-Nicol tuning equations for finding the optimum values of
Kp, T i and Td.
Table 4.1
A A s.t
Proportional only Kp = Kp = 1+ d Kp = 0.5Kc
LM OP
S .t d S.t d 3B
N Q
A = Step change in input Kc = Critical/utlimate gain
B = Ultimate value of output Tc = Critical/ultimate period
td = Dead time
R = Log ratio
0.9A A R
Proportional + Integral [P + I] only Kp = Kp = 0.9 + Kp = 0.45Kc
LM OP
S ×t d S ×t d N 12 Q
5TH PROOF 4-9-09
Tc
Ti = 3.33 td Ti = d Ti =
LM 30 + 3R OP t
N 9 + 20R Q 1.2
1.2A A L RO
Proportional + Integral + Derivative only [P + I + D] Kp = K =
p 1.33 + P Kp = 0.6Kc
S ×t d d
S ×t MN 4Q
Tc
Ti = 2td i
T =M d Ti =
L 32 + 6R OP t
N 13 + 8R Q 2
Tc
Td = 0.5td Td = M d Td =
L 4 OP t
N11 + 2R Q 8
Tuning of Controllers
165
5TH PROOF 4-9-09
SOLVED EXAMPLES
1. A Process exhibits oscillations for a critical gain setting of 2 in a 30
minute period. Find the settings of Ziegler-Nicol for.
(a) P + I Controller mode
(b) P + I + D Controller mode
Solution:
(a) For P + I mode
Kp = 0.45 Kc. Given Kc = 2
Kp = 0.45 ´ 2 = 0.9 Tc = 30 min
Tc 30
Ti = = = 25 mins.
1.2 1.2
(b) For P + I + D mode
Kp = 0.6 Kc = 0.6 ´ 2 = 1.2
Tc 30
Ti = = = 15 min.
2 2
Tc 30
Td = = = 3.75 min.
8 8
2. In the application of Ziegler-Nicol method oscillations are observed in
the process with proportional band set to 400/01 in the time period of
10 minutes, Find the settings of three mode controller?
Solution: A three mode controller is P I D controller.
100
Since, P.B =
Kp
100 100
Kc = = = 2.5
P. B 40
The critical period Tc is given as T c = 10 minutes
\ Kp = 0.6 Kc = 0.6 ´ 2.5
= 1.5
Tc 10
Ti = = = 5 min.
2 2
Tc 10
Td = = = 1.25 min.
8 8
5TH PROOF 4-9-09
3. For the process shown in figure, determine the optimum settings for
PID controller. (Use Routh Array)
Kp 1 c(s)
s(s + 2) (s + 4)
R(s)
Kp
1+ =0
s( s + 2) ( s + 4)
S(S + 2) (S + 4) + Kp = 0
S(S2 + 6S + 8) + Kp = 0
S3 + 6S2 + 8S + Kp = 0
S3 1 8
S2 6 Kp
48 - K p
S1
6
S0 Kp
48 - K p
= 0 Þ Kp = 48
6
\ The critical Gain Kc = 48
Auxillary equation
A(S) = 6 S2 + Kp = 0
6S2 + 48 = 0
S2 = 8
S = ± j 2.8
Comparing with S = jw
S = ± j2.8 jw
w = 2.8 r/s
The critical period Tc will be
2p 2p
Tc = = = 2.24 secs.
w 2.8
5TH PROOF 4-9-09
Tc 2.24
Ti = = = 1.12 secs
2 2
Tc 2.24
Td = = = 0.285 secs.
8 8
A first order process has transfer function
e -2 s
G(S) =
1 + 4s
A Test run is conducted by subjecting to a sudden disturbance of
magnitude.
4. The ultimate value of output stabilized at 3.5. Find the settings of P + I
controller?
Solution: Comparing the transfer function
e -2 s e - td s
=
1 + 4s 1 + Ts
td = 2 secs.
T = 4 secs.
Given A =5
B = 3.5
B 3.5
Slope [S] = = = 0.875
T 4
S td 0.875 ´ 2
Log ratio [R] = = = 0.5
B 3.5
The settings of P + I mode are
0.9 A
Kp = ; Ti = 3.33 td
S td
0.9 ´ 5
Kp = = 10.3
0.875 ´ 2
6. What is the critical gain and critical period of the system whose
characteristic equation is
S3 + 4S2 + 6S + (4 + K) = 0
Solution:
S3 1 6
S2 4 (4 + K)
24 (4 + K )
S1 0
4
S0 (4 + K) 0
24 (4 + K )
=0 Þ K = Kc = 20
4
® The Auxillary equation is
A(S) = 4S2 + (4 + K) = 0
4S2 + (4 + 20) = 0
S2 = 6
S = ± j 6 jw
w = wc = 6 r/s.
2p 2p
Tc = Þ = 2.61 secs
wc 6
5
+0)26-4
Implementation of
Control Actions
Inverting V1
input –
A
Non-inverting V + V0
2
input
FIG. 5.1
5TH PROOF 4-9-09
The differential inputs are designated by the signs (+) and () notations.
These polarity signs refer to the phase relationship between individual input
signals and the output signal.
The output signal is in phase with the signal applied to the input pin
marked + and hence is known as noninverting input and anti phase to
the signal applied to the pin marked and is known as inverting input.
When the power supply connections are also shown, they are represented
as V + and V for positive and negative voltages, respectively.
The characteristics of ideal op-Amp are
Infinite input impedence.
An ideal amplifier can be connected to any signal source with no
loading effects.
Infinite gain A.
Zero output impedence
An ideal op-Amp is capable of driving any load.
Infinite bandwidth:
An ideal op .Amp can be made to amplify any signal frequency from 0 to
¥ Hz. without attenuation.
Infinite common mode rejection ratio so that the output common mode
noise voltage is zero.
Infinite slew rate so that output voltage changes occur simultaneously
with input voltage changes.
However, practical op. Amps can be made to exhibit some of these
characteristics using a negative feedback arrangement.
The input impedence, output impedence and bandwidth can be brought
close to ideal values by incorporating negative feedback arrangement.
When operated in negative feedback configuration the output of an op.
Amp does not switch between positive and negative saturation levels.
In the implementation of control actions initially the ideal characteristics
are assumed and Laplace transformation of the inputs and impedences is
considered for the case of understanding and analysis.
Vin ( s) - V1( s) V ( s) - V0 ( s)
= 1
Z1( s) Z2 ( s)
Z2(s)
I(s)
Vin(s) Z1(s) –
V1(s)
I(s) A
+
V0(s)
FIG. 5.2
Vin ( s) - V0 ( s)
=
Z1 ( s) Z2 ( s)
\ V0(s) =
LM Z ( s) OP .V (s)
2
...(1)
N Z ( s) Q
1
in
Z2 ( s)
The amplifier gain A = .
Z1( s)
By properly choosing the values of z1(s) and z2(s), the op. Amp can be used
for various applications including the controllers variety.
On/off Controller
An on/off controller or two position controller can be implemented in two
designs, namely,
Electro mechanical design.
Electronic system design.
They find extensive use in many household applications like air
conditioning systems and heating systems in particular.
The electro mechanical design involves the use of bimetallic strip and
mercury switch.
The inertia of the mercury tends to keep the system either in ON or OFF
position until the temperature changes from the set point value.
This provides the required neutral zone the prevent excessive cycling of
the system.
The electronic implementation of ON/OFF controller with adjustable
neutral zone using op. Amps is shown in Figure 5.3.
R2
Vi R1 R1
R
– R Vo
– Vout =
– 0
+
+
+
+ Vs
Comparator
Summer Inverter
Vref
FIG. 5.3
R1
VL = Vref . VO ...(3)
R2
The width of the neutral between VL and VH can be adjusted by varying the
resistor R2 in Fig. 5.3. The extent of width of neutral zone can be calculated
from the difference of equations (2) and (3). The typical value of inverter
resistance may be in the range of 1 W to 100 K W range.
Vout
(ON) V0
(OFF) O
Vi
VL VH
FIG. 5.4
Multiposition Controller
The implementation of three position controller (multiposition) using two
comparators, a summing amplifier and an inverter is shown in figure 5.5.
Vout
Vout 2
Vout 1
O Vi
Vref (1) Vref (2)
(Set point 1) (Set point 2)
FIG. 5.5
5TH PROOF 4-9-09
Figure 5.5 shows the response of three position controller with bias voltage
zero.
The bias voltage is the minimum output voltage when the input voltage Vi
is less than Vref (1) or setpoint 1.
In the circuit shown in figure 5.6, the comparator output is assumed to be
zero for OFF state and VO volts for ON stage. The bias voltage input is
included in the summing amplifier represented as VB.
VB (Bias Voltage)
R3
Vi – R2 R3
1
Vref (2) +
–
2 – R
–
+ R1
Vref (1) + 0
+ Vout =
Vo
Inverter
Comparators Summer
Amplifier
FIG. 5.6
V1 - V V - Vout
\ = [(V = 0) virtual GND]
R 1/ CS
J
- V1 -8
Vout = Vout = J
RcS 4+
5TH PROOF 4-9-09
1
Cs
+ V1
Vi Vref 1 R
O –
Vref 2 Vi V
+
– V1
Vout
Integrator
FIG. 5.7
During the interval Vref 1 < Vi < Vref 2 the output will float at previous
setting, on whatever magnitude of the Vout.
V1
- = - Vout . CS
R
t2
V1 V
\ Vout = = 1 t
RC S RC 0
Times t1 and t2 indicate the intervals when the error is in the limits of
neutral zone.
R2 R2
Vo R2
R1
Ve – R
–
+ Vout
+ + Vout
Summer Inverter
FIG. 5.8
R2
Vout = .Ve + V0.
R1
and the output of inverter will yield
R2
+Vout = .Ve + V0. ...(6)
R1
Integral Controller
The analytic equation of integral controller from previous chapters is
reproduced below.
P=
1
Ti z e dt + P(0). ...(7)
1
Cs
R
R
Ve – R
V –
+ Vout
+ + Vout
Integrator
Inverter
FIG. 5.9
Ve - V V - Vout
= (Since V = 0)
R 1/ CS
Ve
Vout =
RCS
V
Þ +Vout = e (o/p of inverter).
RCS
1
= + Vout =
RC z
Ve dt + V(0)
Derivative Controller
Practically, Derivative controller mode is not used alone for the reasons
already stated in chapter 2. However, since this controller is used in
combination with other modes its electronic implementation is discussed here.
The analytic equation of derivative controller is
de
P = Td ...(9)
dt
Thus, to get this kind of response a differentiator along with an inverter is
employed.
The analysis of differentiator yields the equation given by (9).
Ve - V V - Vout
=
1/ CS R
-Vout
Ve.CS =
R
5TH PROOF 4-9-09
R
R
1
Cs
Ve – R
V –
+ Vout
+ + Vout
Differentiatior Inverter
FIG. 5.10
Þ Vout = RCSVe.
dV
Vout = RC e
dt
dVe
O/p of Inverter will be + Vout = RC ...(10)
dt
Comparing equations (9) and (10), the derivative time Td = RC
Vout = controller output
Ve = Error voltage (Input).
From practical perspective the Fig. 5.10 cannot be used because it tends
to be unstable, i.e. it may exhibit oscillations in the output voltage. The
occurrence of this instability is because the magnitude of the output voltage
increases linearly with frequency. Therefore, when the frequency is infinity, the
output tends to be infinity.
Consider the equation
Vout = RCSVe
The o/p of the inverter will be
+Vout = RCSVe
Put S = jw
+Vout = jwRCSVe
|Vout| = 02 + ( wRCVe )2
|Vout| = wRCSVe
Since w = 2p f
|Vout| = 2p f RCSVe. ...(11)
Thus, |Vout|af which implies a little high frequency noise will cause large
excursions in the output voltage which is not acceptable in control system
design.
5TH PROOF 4-9-09
R2
1
R1 Cs
Ve –
V
+ Vout
FIG. 5.11
Ve - V V - Vout
= (Since V = 0).
R1CS + 1 R2
CS
VeCS -Vout
=
R1CS + 1 R2
R2CSVe
= Vout.
R1CS + 1
Put s = jw
jwR2CVe
Vout =
jwR1C + 1
wR2C|Ve|
|Vout| =
1 + ( wR1C)2
Since w = 2p f.
2pf R2C|Ve|
|Vout| = ...(12)
1 + (2pf R1C)2
For frequencies in practical range.
2p f R1C << 1
|Vout| = 2p f R2C |Ve|. The circuit acts as a derivative controller.
5TH PROOF 4-9-09
2p f R2 C |Ve|
|Vout| =
2p f R1 C
R2
|Vout| = |V |. It acts as a fixed gain amplifier. Thus, it can be stated
R1 e
that this modification essentially clamps the gain above some frequency to a
constant value. And it is rest assured that the clamped frequencies are well
above anything that could occur in the actual control system.
Composite Controllers
Proportional + Integral Controller
It is basically the combination of proportional and integral actions,
respectively. The analytic equation of this combination is
P = Kpe +
Kp
Ti z e dt + P(0) ...(13)
1
R2 Cs
R
R1
Ve – R
V –
+ Vout
+ + Vout
PI Controller Inverter
FIG. 5.12
Ve - V V - Vout
=
R1 R2CS + 1
CS
Since V = 0 (Virtual Ground)
Ve ( R2 CS + 1)
= Vout
R1CS
5TH PROOF 4-9-09
Ve . R2 CS Ve
Vout = +
R1 CS R1CS
R2 Ve
Vout = .Ve +
R1 R1CS
The output of the inverter will be
R2Ve Ve
+Vout = +
R1 R1CS
Þ +Vout =
R2
R1
R
.Ve + 2 ×
1
R1 R2C z Ve dt + V(0) ...(14)
Ve – R
V –
+ Vout
+ + Vout
1
Cs
P + D Controller Inverter
FIG. 5.13
Ve - V V - Vout
=
R1 R2
R1CS + 1
Since V = 0 (Virtual ground)
5TH PROOF 4-9-09
Ve [ R1CS + 1] -Vout
=
R1 R2
R2 R
Vout = .Ve + 2 .R2CSVe
R1 R1
o/p of inverter will be
R2 R dV
+Vout = .Ve + 2 .R2C e ...(16)
R1 R1 dt
\ Comparing equations (15) and (16)
R2
Kp = ; Td = R2C
R1
P = Kpe +
Kp
Ti z e dt + KpTd
de
dt
+ P(0) ...(17)
where P(0) = Initial value of controller o/p. There are many modes of
realisation of this control action but with minimum number of op-Amps its
implementation is shown in Fig. 5.14.
1
R2 C2s
R
R1
Ve – R
V –
+ Vout
+ + Vout
1
C1s
FIG. 5.14
Ve - V V - Vout
=
R1 R2C2S + 1
R1C1S + 1 C2S
(Since V = 0)
Ve [ R1C1S + 1][ R2C2S + 1]
= Vout
R1C2S
5TH PROOF 4-9-09
Ve
\ Vout = {[R1C1R2C2S2 + [R1C1 + R2C2]S + 1]}
R1C2 S
\ Vout = R1C1SVe +
LM C
1
+
R2
V +
VeOP
NC2 R1 e R1C2SQ
C1 R
Since the ratio < < 2 and at the output terminals of inverter, the
C2 R1
above equation reduces to
+Vout =
R2
R1
dV
.Ve + R2C1 e +
dt
1
R1C2 z Vedt + V(0).
Þ +Vout =
R2
R1
R
Ve + 2 .
1
R1 R2 C2 z Vedt +
R2
R1
dV
.R1C1 e + V(0)
dt
...(18)
Flapper
Pivot Flapper motion
Gap
Supply
Pressure
(20 psi) (Restriction)
Signal pressure
FIG. 5.15
Signal Pressure
20 psi
(Supply pressure)
15 psi
Linear characteristics
3 psi
Gap Distance
FIG. 5.16
Spring
Assembly
Flapper
Pivot
Gap
4-20
N (mA)
o
z
Supply z
Pressure l Coil
(20 psi) e (Force)
3-15 psi
signal
FIG. 5.17
Pneumatic Pilots
These pneumatic elements are basically used for improvisation of the speed of
pneumatic transmission. One such method is pneumatic amplifier as shown in
Fig. 5.18.
Air Supply
Output
Pressure
Leak
Signal Pressure
FIG. 5.18
The pneumatic amplifying pilot is also known as booster or relay. Its main
function is to raise the air pressure and/or air flow volume by some linearly
proportional amount from the input signal. For the input of 3 15 psi, if the booster has
a pressure gain of 100 the output would be 300 1500 psi.
In the figure 5.18 if the signal pressure is such that the diaphragm motion moves
the plug in the body block of the booster down the gas leak is reduced and the pressure
in the output line is increased.
Similarly, if the diaphragm moves the plug upwards blocking the air supply then,
the pressures in the input transmission line decreases as there will be considerable air
leak.
A part from power amplification there are also there pilots which are used
for overcoming transmission lags. They are made of double bellows assembly
and are directly inserted in the transmission line.
5TH PROOF 4-9-09
Kp La Feedback bellows
rge
all
K p Sm
Flapper
e (Pilot)
Nozzle ‘P’ Output
(Error) Pressure
Air Supply
FIG. 5.19
link is connected in the upper position as shown in Fig. 5.19. The feedback
pressure is small making the proportional gain large (small P.B) thus resulting
in large output pressure change.
Similarly when the feedback link is connected in the lower position, Kp is
small i.e. very large P.B. (Proportional band). By proper design we can choose
the appropriate location for connecting the feedback link to the flapper, so
that the output pressure (P).
P = Kp e.
The effect of variations in temperature, any leakage and supply pressure
causing deviation of the output pressure from the required value is reduced to
minimum because of the negative feedback action of feedback bellows. Also
the P.B and proportional gain may be varied over a wide range of values
(1000 to 1).
)
mall
K p (S
Flapper
Pilot
Nozzle Assembly
e Output
Error Pressure
signal ‘P’
Supply
Pressure
FIG. 5.20
5TH PROOF 4-9-09
Kp (H
igh )
Negative
feedback
bellows
all)
K p (Sm
R
Flapper
Pilot
Nozzle Assembly
error
Output
e Air gap Pressure
‘P’
Supply Pressure
(Air)
FIG. 5.21
5TH PROOF 4-9-09
Positive Negative
f.b. bellows f.b. bellows
Negative f.b.
Restriction (R1)
Flapper
Pilot
Nozzle Assembly
Output
error Pressure
(e) (P)
Supply Pressure
FIG. 5.22
)
C2
2
, R I e
p P+ ons
(K sp
(P) re
Lag
PID
response
R1C1 (Derivative)
t
(Lead)
FIG. 5.23
5TH PROOF 4-9-09
From
Hydraulic
System
Filter
Sump
FIG. 5.24
The electric motor which is usually of gear type pumps oil from the sump
through filter and passes it to the hydraulic system with certain pressure. A
relief valve is provided in the supply line which when opened sends oil back
into the sump during any untoward emergencies.
The sump collects the return oil from the system and serves to dissipate
heat, which is very important in hydraulic system. Sometime oil coolers are
also used for dissipating the heat.
5TH PROOF 4-9-09
To Hydrautic System
Oil
Compressed Diaphragm
Nitrogen Gas
Accumulator
FIG. 5.25
Oil
Error P
e
Power Power
Cylinder Piston
From To
Pump Sump
FIG. 5.26
i.e.
dp
dt
µ e ® P = KI z e dt.
Unlike the proportional action, the power cylinder outlets are connected
simultaneously to pump and sump depending upon the magnitude of the
error signal.
dP/dt
Oil Oil
Error
Signal
(e)
From To
Pump Sump
FIG. 5.27
5TH PROOF 4-9-09
B
a
l a
a r
n m
c
e
Diaphragm
Jet
Pipe
Power Power
Cylinder Piston
From To
Controlled Pump Sump
Variable
“C”
FIG. 5.28
5TH PROOF 4-9-09
QUESTION SET
1. With next diagrams obtain the implementation of composite control
actions electronically using op-Amps?
2. What are the different types of pneumatic elements used in process
industries? Explain with next diagrams?
3. With a next diagram explain the operation of P to I converter?
4. Explain with next diagrams:
(a) Pneumatic P + I controller
(b) Pneumatic P + I + D controller.
5. Obtain the mechanism of hydraulic implementation of controller
modes? Explain with next diagrams?
$
+0)26-4
6.1 INTRODUCTION
One of the most important components of any process control system is the final
control element. The output of the automatic controller is fed to the final control
element which directly affects the manipulating variable to achieve the desired
control objective. Thus, a final control element may be defined as
The mechanism which alters the value of the manipulated variable in
response to the output signal from the automatic controller.
The final control element, abbreviated as FCE essentially consists of two
parts
(1) An ACTUATOR, which is primarily used to translate the output signal of
the automatic controller into a position of a member exerting the required
force or displacement or power.
(2) A control valve to adjust the value of the manipulated variable.
Some of the other actuators are operated by compressed air, hydraulics and
also by electrical means. Thus, the actuator designs depend on the particular
signal with which it is operated upon as enlisted above.
Actuator Types
Electromechanical Actuators
Electro-hydraulic Actuators
Electrical Types
Pneumatic Actuators
Hammer
blow
Gear drive
Toroue De-clutching
Drive
control handwheel
slfeve
FIG. 6.1
5TH PROOF 4-9-09
Reversible
Low-Inertia
AC Motor
Continously
Connected
Handwheel
Motor Drive
Pinion Gears
Load Connection
Linear Motion
FIG. 6.2
In Fig. 6.2, a linear output actuator which uses a worm and a rack and Pinion
arrangement to translate horizontal shaft motor output to vertical linear motion is
shown.
This actuator is designed with a conventional globe valve bonnet for ease of
mounting.
A continuously connected hand wheel, which must rotate the rotor of the
motor, can be used when there is a short stem travel and relatively low force
output. However these electrical actoators employing output gear system for
reduction in speed employ AC motors and/or D.C. motors respectively.
The maximum force output in this case is 680 kgf at about 2.5 mm/minute.
Since the subject matter of electric motors is extensive, the reader is advised
to refer the related texts for more details. In our present context, the application
of these motors is stressed which enables them to be used for actuation purposes.
Apart from worm gears showed in Fig. 6.1, of another most common actuator
employed in industrial application is the spur gear actuator.
The spur gear actuator shown in Fig. 6.3 in contrast to worm gears has very
low power loss. The advantage of spur gear is that it prevents the back driving
torque on the actuator even if the load possesses frictional characteristics.
The spur gear actuator consists of adjustable cams to operate limit and
auxillary switches. A potentiometric or tachogenerator feedback calibrated with
respect to the rotation of the actuator is required for formulating a control circuit.
Actuators of this type must be adaptable to mounting on and operating a variety
5TH PROOF 4-9-09
Output
Shaft Adjustable
Cams
(Optional)
Travel Potentiometer
Limit
Switches Auxiliary
Switches
FIG. 6.3
of quarter turn valves. Hence, the opening of control valve can be controlled by a
limit switch whereas for shutting it is advisable to incorporate a torque limiting
device to sense the closure of the valve.
Both the spring and diaphragm actuator and the piston actuator produce
linear motion to move the control valve. They are ideal for use on valves whose
stem motion is linear, i.e., globe valves. For rotary valves such as butterfly valves
some mechanism of linear to rotary motion conversion must be employed.
Þ PA = KX + Pv Av ...(3)
p A
Flow Flow
At PV
AV
(a) (b)
Direct acting actuator Reverse acting actuator
FIG. 6.4
In the reverse acting mode the valve stem is pulled upwards on an increase of
air pressure. Equation 3 is modified as
PA + KX PvAv = 0
Þ KX = PA + PvAv ...(4)
The reverse and direct acting spring and diaphragm actuators are shown in
Fig. 6.4 (a) and 6.4 (b) respectively.
Figure 6.5 indicates plot of diaphragm pressure and stem travel neglecting
friction and fluid forces on the valve plug for equation 3.
-1
ve
100% Cu
r e-2
u rv
Closed C
Stem es
travel f orc
plug ce
s
No for
ug
pl
ith
W
0%
Open
Diaphragm
0 3 6 9 12 15 pressure (Psi)
FIG. 6.5
5TH PROOF 4-9-09
100%
Ideal
Closed
characteristics
Stem Actual
travel characteristics
0%
Open
Diaphragm
0 3 15 pressure (Psi)
FIG. 6.6
Piston
Cylinder
Operating air Operating air
pressure inlet pressure inlet
Cylinder
Piston
Spring
The double acting cylinders are equipped with two inlets for the operating
air pressure on both front and rear ends of the cylinder. Thus, the piston motion
is a two way motion. The movement of the piston depends on the quantity of the
compressed air supplied at the front and reverse ends, and the pressure
difference makes the piston travel possible. The throttling of the control valve in
both single and double acting cylinders is achieved by the piston movement or
displacement.
Cylinder
Piston
Inlet Inlet
supply as against compressed air supply of pneumatic actuation. The single and
double acting cylinders discussed in Section 6.13 are also applicable for
hydraulic actuation.
However, the principle of hydraulic actuation is discussed in this section. A
hydraulic actuator converts a small force Fx into an amplified force called as
operating force Fw
Fx Fw
PH
(A1)
PH(A2) Hydraulic fluid
FIG. 6.10
Fx
PH = ...(5)
A1
A2
\ Fw = Fx × ...(6)
A1
The working piston force can be made very large by changing the area of
forcing piston A1 as seen by equation (6).
Valve stem
Valve plug
Þ Flow Þ Flow
Valve seat
FIG. 6.11
The pressure on the diaphragm of actuator moves the valve stem which in
turn moves the plug. The amount of resistance offered to the fluid flow depends
on the extent to which the plug covers the valve seat. As the pressure on the valve
stem varies over the normal pneumatic range of 3 15 psi, the valve stem motion
also varies from a fraction of inch to several inches, depending on the design of
the valve. Both the direct and reverse acting actuators discussed in Section 6.13,
Fig. 6.4 (a) and (b) can be used for operating the control valve. The direct acting
actuator is also called as air-to-close (or) pressure to close actuator where an
increase in signal pressure above the diaphragm causes the valve stem to move
down. The reverse acting actuator also called as air-to-open (or) pressure to open,
has the air inlet below the diaphragm. An increase in air pressure causes the valve
stem to move upwards. Based on the design of actuators, the valves are also
referred to as AIR - TO - OPEN or AIR - TO - CLOSE valves.
The valve shown in Fig. 6.11 is known as THE SINGLE SEATED
CONTROL VALVE. It implies the valve contains one plug with one seating
surface only.
However, to reduce the effort needed to open the valve against full pressure
drop across it, valves are constructed with double seating as shown in Fig. 6.12.
5TH PROOF 4-9-09
In double seating, there will be two valve plugs and two seatings. The two
valve plugs are attached to the single valve stem.
The flow pattern through the valve is designed in such a way that pressure
drop across seat A tends to open the plug and the pressure drop across seat B
tends to close the plug. This kind of action tends to reduce the effort needed to
open the valve and consequently less expensive actuator is needed.
Valve stem
A
Þ Flow Þ Flow
B
Flow
profile
Air supply
Fd = PA
Fu = Kx Spring
1
B.dx
Fu = dt Valve stem
2
Valve plug
fi Flow fi Flow
Valve body
FIG. 6.13
The air pressure exerted on the top of diaphragm having cross-sectional area
of A m2 causes the diaphragm to deflect downwards. Let p be the pressure
signal that opens or closes the valve. This force acts downwards.
Thus, downward force Fd is expressed as Fd = P ´ A.
The force exerted by the spring attached to the valve stem is upwards and is
expressed as
5TH PROOF 4-9-09
Fu1 = Kx
where, K = spring constant
x = displacement of the stem.
The downward force Fd is further balanced by frictional force exerted
upwards resulting from the proximity of valve stem with valve packing. This
frictional force is expressed as
dx
Fu2 = B ×
dt
where, B = damping coefficient
The third force acting upwards results from the mass of the stem and the
plug and is expressed
m d2 x
Fu3 =
g dt2
where,
M = mass of the stem/plug
g = acceleration due to gravity m2/s
From the principles of translational mechanics
Algebraic sum of downward forces = Algebraic sum of upward forces
Fd = Fu1 + Fu2 + Fu3
dx M d 2 x
PA = Kx + B + ...(7)
dt g dt 2
Equation 7 indicates that the valve stem follows second order dynamics with
pneumatic spring and diaphragm actuator.
Practically the mass of the valve stem and plug are neglected and equation (7)
is approximated to that of first order dynamics.
dx
B + Kx = PA
dt
B dx PA
+x=
K dt K
dx
T× + x = Kp × P. ...(8)
dt
Thus, the stem position x for the applied pressure P follows first order
dynamics. The presence of time constant element indicates sluggishness of the
response for a given input. But generally most of pneumatic valves respond very
rapidly to changes in air pressure. The valve plug opening or closing is so fast
that even the first order dynamics also may be neglected.
In such a case only a constant gain term will remain which relates the output
from the controller, i.e. air pressure signal to the fluid flow through the valve.
This approximation further makes the control valve to exhibit zero order
dynamics. Equuation 6.8 further reduces to
x = KP × P ...(9)
xµ P ...(10)
(b) Linear
In this type of control valve the flow rate varies linearly with stem motion.
Let Q = Flow rate (m2/s)
Qmax = Maximum flow rate (m3/s)
X = Stem position (m)
Xmax = Maximum stem position (m)
Q X
= ...(11)
Q max X max
Equation (11) describes the operation of linear valve, for a constant pressure
drop.
Q
Let Q¢ = where, Q ¢ = fraction of maximum flow.
Q max
X
X¢ = X ¢ = fraction of maximum lift.
X max
We can write a general equation for flow through the valve for constant
pressure drop as
Q¢ = f(X¢)
Sensitivity is defined as
dQ ¢
S= ...(12)
dX ¢
dQ ¢
= K; where K =
For linear valve, equation (12) must be constant. Thus,
dX ¢
any constant, the relationship between flow and valve opening can be
represented as straight line.
Q max
R= ...(13)
Q min
The flow rate through equal percentage valve is given in terms of R as
Q(Actual flow rate) = Qmin × R(X/Xmax) ...(14)
The valve characteristics for the three types listed above is shown in
Fig. 6.15.
Quick opening
(c)
(1)
1 ® Quick opening
Flow rate through valve (Q)
2 ® Linear
3 ® Equal percentage
(2)
(3)
In practical applications the quick opening valve may be used as shut off
valve and for throttling purposes equal percentage is best suited.
Linear valve, whose sensitivity is constant is suitable for both complete shut
off and throttle applications.
5TH PROOF 4-9-09
FC FO
Three way Fail Closed Fail Open
Angle Value
FIG. 6.16
FIG. 6.17
5TH PROOF 4-9-09
Inline Globe
FIG. 6.18
5TH PROOF 4-9-09
FIG. 6.20 Top & bottom guided invertible Single Seated Globe Valve
The top entry, top guided single seated globe valve has only one body
opening for the bonnet and valve plug mass is reduced, which thereby increases
the trim throttling precisely when compared to the top and bottom guided valve.
When dealing with fluids containing solids, gummy or highly viscous fluids,
and highly corrosive fluids stem guided valve is always a better choice than top-
guided valve because it minimises cavities and ensures trouble free operation.
Cage Valves
One of the most popular variant of single seated globe valve is the cage valve. The
valve plug design and top entry bonnet makes the job of trouble shooting work
very easy, especially when changing the trim or plug. A wide variety of trim or
valve plug designs such as anti-cavitation and reduced-noise trims can be
conveniently installed.
5TH PROOF 4-9-09
There are two basic design configuration available for cage valves.
(a) cage valve with clamped-in seat ring and characterized plug shown in Fig.
6.21
(b) cage valve with unbalanced plug and flow characterized ports shown in
Fig. 6.22.
Very Heavy Guiding:
Two widely-spaced
guides on heavy
stem. plug does not
guide in retainer.
Interchange Ability:
class 600 body used
for class 150, 300 Clamped-in seat ring:
and 600 Ratings valve may be dis-
assembled quickly,
easily by removing
four bonnet bolts
FIG. 6.21
FIG. 6.22
The cage valve design shown in Fig. 6.21 uses cage solely to clamp the seat
ring into the valve body.
5TH PROOF 4-9-09
This configuration is usually stem or top guided design and the valve plug
does not contact the cage.
In Fig. 6.22, the cage is used to guide the plug as well as clamp the seat ring
into the body. The cage openings are shaped in such a way so as to provide
desired flow characteristics and the valve plug acts as a curtain to expose the
ports.
Angle Valves
These valves are generally used in piping arrangements which require high
pressure drop services and applications such as in controlling molten metals,
polymer materials, cryogenic fluids and liquid slurries.
The Y-shaped angle valve is a popular design for the above mentioned
applications and is shown in Fig. 6.23.
The valve can be installed in any of the vertical, horizontal or angled piping.
They are usually fitted along with a vacuum jacket, which provides the required
thermal insulation when used in dealing with flow of cryogenic fluids such as
liquid hydrogen. The vacuum jacket is designed with metal bellows to allow for
thermal changes or mechanical tolerances and also has a provision for welding to
the similar jacket around the adjacent piping.
The cooling rate depends on valve design, thickness etc. Since these valves
are used in high pressure applications they are highly susceptible to cavitation
problems (discussed later) and erosion problems.
5TH PROOF 4-9-09
Ball Valves
A typical ball valve consists of a spherical valve through the valve body. They are
available in three types
(a) conventional or quarter turn pierced ball type
(b) characterized type
(c) Cage type.
FIG. 6.24
5TH PROOF 4-9-09
The structural details of top entry quarter turn pierced ball valve is shown in
Fig. 6.25.
Stem
Body
Ball
Seal
Tea Ports
Angle Ports
Porting arrangements of multiport ball valves.
The size of conventional ball valve varies between ½ inch to 42 inches (12.5
mm to 1.06m) [ANSI class 150 standards]. The design temperature varies
between 150°C to 315°C, with special designs available up to 1020°C also.
They are usually made of cast or bar stock brass or bronze, carbon steel,
stainless steel, aluminium, monel etc for the body.
For balls, forged naval bronze, carbon steel, plastics, glass, ceramics etc. The
seats are usually made up of teflon Delrin, neoprene etc and related synthetic
materials.
5TH PROOF 4-9-09
Bearings
Seal
Rings
Body
It consists of a venturi-ported body with two seat rings, a ball that acts as a
plug (valve plug) or trim. A stem is used for positioning the cage which in turn
position the ball to obtain fully closed or open and/or throttling flow
characteristics. The structural details and various positions enlisted above are
shown in Fig. 6.28.
The cage rolls the ball out of the seat as it is tilted by the stem supporting it and
positions the ball firmly during throttling and lifts it out completely during full opening
position from the flow stream. The cage is designed with four inclined surfaces. The two
surfaces next to the downstream seat lifts the ball out of the seat and roll it over the top
edge of the seat ring thus opening the valve.
During fully open position, the ball rolls down to the centre of the cage to
rest on all four inclined surfaces. The Bernoullis effect of the flowing stream
holds it cradled in this position during the rest of the operation. They are
available in sizes ranging from ¼ inches to 14 inches for design pressures limited
to 17 mpa. They can withstand temperatures in the range of 254°C to 980°C.
The ball is usually made up of stainless steel.
Butterfly Valves
The butterfly valve is perhaps the most popular and oldest designs available in
process industries. The earlier designs were particularly used in the low pressure
water works applications. However the modern butterfly valves are used in a wide
variety of fluid applications, those with high pressure drop, tight shut off and
also involving corrosive & erosive characteristics.
The butterfly valve is defined as A damper or vane or throttle valve in a fluid
pipe line consisting of a disc turning on a diametrical axis.
The operation of butterfly valve therefore is based on the above definition of
it where it involves only rotating the vane, disc, a flapper or a louver by means of
a shaft to which it is fastened. The rotation is accomplished using Pneumatic,
electric or hydraulic motor drives attached to the shaft by various methods. The
5TH PROOF 4-9-09
fully closed to fully open position is achieved when the disc or vane rotates or
moves through 90° angle. The free area developed by the disc as it moves towards
the full open position provides the throttling operation.
Based on the construction and flow characteristics the butterfly valves are
available in two designs.
(a) Swing through butterfly valves
(b) Shut-off butterfly valves (or) High performance butterfly valves [HPBV].
The vane positions are shown in Fig. 6.29 for fully open, closed and
throttling.
FIG. 6.29
Packing
Bearings
Disc
DP
q =CV ...(15)
G
CV = factor associated with capacity of valve
(control valve coefficient)
q =design flow rate (gpm)
G =Specific gravity relative to water.
DP =Pressure drop across the valve (Psi)
The size of the control valve depends upon the factor CV i.e. valve coefficient
value. The valve coefficient CV is defined as the flow (gpm) of a fluid of unit
specific gravity through a fully open valve, across which a pressure differential of
1.0 psi exists. In equation 15 by putting q = DP Þ G = 1; CV = 1 i.e. the
definition just stated above.
The value of CV varies with the design of the value i.e. shape, size, reynolds
number etc.
Typical values of CV for different valve are shown in the table.
5TH PROOF 4-9-09
Table 6.2
Valve Size (Inches) CV
1
1. 0.3
4
1
2. 3
2
3. 1 14
1
4. 1 35
2
5. 2 55
6. 3 08
7. 4 174
8. 6 400
9. 8 725
Applications
used for processes where instant large flow is required.
used for frequent on/off service.
Best suited control for globe and ball valves.
Now, let us see a quick comparison between different types of valves.
Globe Valves
The recommended uses are for
1. Throttle service/Flow regulation.
2. Frequent operation.
It can be used for fluids that are liquids, vapours, gases corrosive substances,
and slurries.
Advantages
1. Efficient throttling
2. Accurate flow control
3. Available in multiple ports
Disadvantages
1. High pressure drop
2. Expensive compared to other valves.
Ball Valves
The recommended uses are for Fully open or fully closed, limited throttling
applications.
They can be used for high temperature fluids, most liquids and slurries.
Advantages
1. Low cost
2. High capacity
3. Relatively lower leakages & maintenance.
4. Tight sealing
Disadvantages
1. Poor throttling characteristics
2. Prone to cavitation
Butterfly Valves
The recommended uses are for fully open or fully closed or throttling services,
minimal fluid trapping in line and frequent operation.
They can be used for fluids such as liquids, gases, slurries, and liquids with
suspended solids.
Advantages
1. Low cost and maintenance
2. Low pressure drop
3. Good flow control
5TH PROOF 4-9-09
Disadvantages
1. High pressure required for valve control
2. Prone to cavitation at lower flows
QUESTION SET
1. Discuss about the spring and diaphragm pneumatic actuator in reverse
and direct acting modes?
2. Explain the principle of hydraulic actuator?
3. Obtain the dynamics of the control valves and the related
approximations?
4. Discuss about the control valve characteristics?
5. Draw the types of valve plugs for different control valves characteristics?
6. Discuss about the different types of sliding stem valves?
7. What are the basic designs available in rotary valves? List their advantages
over conventional globe valves?
8. What is meant by quarter turn in ball valve design? Write about the ball
valves generally used in industries?
9. Explain with neat diagrams the vane positions for butterfly valves during
fully open, closed and throttling?
10. Write about
control valve selection
sizing of control valves
cavitation in control valves.
+0)26-4
%
Process Control System
Configurations
INTRODUCTION
The main objective of a process control system is to monitor and control the
controlled variable in a process and keep it at set point, the desired value, in spite
of disturbances occurring within the process or through any external agencies.
The are various configurations in setting up of process control systems to regulate
and control the dynamic variable in a process. This chapter outlines the
structures and schemes of developing different types of process control loop
configurations. However, the emphasis is more on structural representation than
on mathematical designs and considerations. Also the selection of a particular
type of configuration depends on characteristics and dynamics of the process,
and is usually made by process design engineers.
Manipulated Variable
Controller
Set Point
Load
Variables .
or . Process
. Controlled
Disturbances Variable
FIG. 7.1
In Fig. 7.1, the deviations of controlled variable from set point value due to
variations in load variable or disturbances are first measured through feedback
elements, compared with input or set point at the controller and the extent of
deviation or error is eliminated by converting it into changes in the manipulated
variable and sending back to the process to restore the balance. Consider the
liquid level control system employing feedback control, shown in Fig. 7.2.
Set Point
Inlet Controller
LT
h
Outlet
FIG. 7.2
Whenever, the inlet flow changes inducing error in the plant or process, the
liquid level is brought back to the set point value by the controller after evaluating
upon the magnitude of error.
In contrast to this control loop configuration, FEED FORWARD provides a
more direct solution to control than finding the correct value of the manipulated
variable by trial and error, as observed in the feed back control.
In feed forward control system, the major components of the load variables
or disturbances are entered into the dynamic model to calculate the value of the
manipulated variable required to maintain the controlled variable at the set point
value.
Fig. 7.3 illustrates the structural scheme of feed forward control loop.
Fig. 7.3 illustrates how information flows forward from the load to the
manipulated variable input of the process. When a change in load is sensed, the
manipulated variable is automatically adjusted to the correct value at a rate that
keeps the process continually in balance.
In feed forward control loop the effect of disturbances or load variables are
not seen in the process since the control action starts immediately after a change
in the disturbance has been detected.
5TH PROOF 4-9-09
Set Point
Controller Manipulated
Variable
Load
Variables
or . Process
. Controlled
Disturbances .
Variable
FIG. 7.3
Let us now consider the same liquid level control system employing feed for-
ward control. The disturbance is the inlet flow rate changes and the manipulated
variable is the outlet flow rate.
The changes in the inlet flow rate are measured, before they affect the
controlled variable (level) and the controller adjusts the outlet flow rate
accordingly to maintain the control objective. As seen in Fig. 7.4 the changes in
the load variable or disturbances are never felt in the process.
P/I
Controller Set Point
Converter
D.P. Cell
Inlet
flow (orifice)
rate
Outlet
FIG. 7.4
(1) It acts before the effect of disturbance (1) It acts after the disturbance affecting
has been felt in the process. the control objective of the process.
(2) It requires identification of all possible (2) It does not require any identification
disturbances and their measurement. and measurement of load variables or
disturbances.
(3) It is more stable and does not introduce (3) It is relatively less stable and induces
instability when evaluating the closed instability in closed loop response.
loop response.
5TH PROOF 4-9-09
(4) It is a better control configuration for (4) It is unsatisfactory for slow processes
multi capacity or slow processes. as there could be a scope of large
inherent lags.
(5) Feed forward control is a suitable option (5) It is unsuitable for processes having large
for processes having significant dead or significant dead time.
time effects.
(6) Any changes in the process parameters (6) It is insensitive to process parameter
excluding disturbances and manipulated changes.
variables cannot be compensated by
feed forward controller because their
effect cannot be detected.
(7) It requires good knowledge of the (7) It is insensitive to modelling errors.
process model.
(8) It is suitable for processes where load (8) It is not suitable for such processes
or disturbances to the process change as it cannot cope up with rapidly changing
rapidly with time. quantities.
(9) The use of feedback in a feed (9) With feed forward, the feedback cont-
forward system does not detract from roller must only change its output by an
the performance improvement which amount equal to what the feed forward
was gained by feed forward control. system fails to correct.
(10) The feed forward system is more (10) They are comparatively less costly and
costly and requires more engineering require less complicated design methodo-
effort. logies.
Inlet (Stirrer)
TT
(Control valve)
P/I
Converter
Outlet
Set point Feed forward
controller
(Control valve)
Steam flow
Let us now see the generalized block diagram for feedback-feed forward
control loop configuration. This representation is helpful in mathematical
analysis of these control loop configurations.
Fig. 7.7 is the generalised block diagram showing both control actions.
R2(s) D(s)
– Hm2(s)
Gc2(s)
R1(s) + +
Gc1(s) + Gp(s) C(s)
–
Hm1(s)
FIG. 7.7
Consider the liquid level control system employing feedback control action
shown in Fig. 7.2. The manipulating variable is the outlet flow rate directly
affecting the level of the liquid in the tank, which is the controlled variable in the
process.
This inherent interaction between the outlet flow rate and level of the liquid
in the tank is exploited to provide better overall control through cascade
configuration. A cascade control system consists of two control loops, namely
primary and secondary control loops respectively.
The loop that measures the controlled variable i.e., level of the liquid in the
tank (h) is the dominant or primary or master control loop and has its setpoint
supplied by the operator.
The loop that measures the outlet flow rate with the output of the primary or
master controller as its set point is called secondary control loop or slave loop.
Together with these loops exploiting the inherent interaction occurring between
the two control systems is known as cascade control loop.
The major advantage of this configuration is if at all any disturbances
occurring within the secondary loop are corrected by secondary or slave
controller before they can effect the value of the primary controlled output.
5TH PROOF 4-9-09
Inlet
Primary (or) Set Point
Master
Controller
LT
Set Point
Secondary (or)
P/I
Slave
Converter
Controller
h
(level) Control
D.P. Valve
Cell
Orifice Outlet
flow rate
FIG. 7.8
Disturbance Disturbance
236 Process Control Engineering
Primary or Secondary + +
Process Process
Master or Slave + +
Set –2 –1 Output
– Controller – Controller
Point
Measurement
5TH PROOF 4-9-09
Secondary Loop
Measurement
Process-2 is the outlet pipe line with outlet flow rate as the output i.e.
manipulating variable. This cascade configuration ensures that the stability of the
overall system is never affected due to disturbances in process-1 or process-2
thereby providing better overall control.
Inlet Stirrer
Set Point
Primary (or)
Master Controller TT
Set Point
Secondary (or)
Slave Controller
FT
Steam
flow rate (Steam flow line)
Outlet
FIG. 7.10
Cascade loops consist of two or more controllers in series, but have only a
single, independently adjustable set point by the operator to the master or
primary controller.
The use of secondary loop and slave controller whose set point is
determined from the output of master controller is to act as the first line of
defense against disturbances, preventing these upsets from entering the
primary process.
The secondary loop must have faster dynamics than the primary loop for
the cascading effect to yield better results.
Some rules of thumb suggest that slaves time constant should be under 1/
4th to 1/10th that of the master loop and the slaves period of oscillations
should be below 1/2 to 1/3 that of the master control loop. Thus, an even
distribution of time constants between master and slave loops is obtained.
The even distribution of time constants between the two loops eliminates the
problems related to instability and saturation in the overall cascade control
configuration.
X
Ratio (R) =
Y
The most common method manipulates a flow loop whose set point is
calculated as
X = YR
The set point for the flow controller is generated by an adjustable gain device
known as the ratio station.
A common example is when the ratio of two reactants must be controlled. A
ratio control system should maintain a constant ratio of controlled flow i.e., X to
wild flow i.e. Y.
Consider the Figure 7.11 where reactant Y is measured but not regulated
and reactant 'X' is both measured and regulated to provide specified constant
ratio of
5TH PROOF 4-9-09
X
R= .
Y
The flow rate of reactant Y is measured and added, with relevant scaling, to
the measurement of flow rate X. The controller reacts to the resulting input
signal by adjustment of the control valve in the reactant X input line.
Reactant 'Y'
Flow sensor FT
for 'Y' Ratio X
station (or)
¸
Controller or
Set point +
FT
Reactant 'X'
FIG. 7.11
1 B2
R2 =
3 A2
5TH PROOF 4-9-09
For example, let the true flow ratio of the additive to gasoline be 2.0 cc/
gallon.
If the additive flow scale is 0 to 1200 cc/min and gasoline flow scale is 0 to
500 gal/min then, the setting of ratio R would be
Set point
Controller
Pressure vessel
Gas in V1 Pressure
vessel V2 Gas out
FIG. 7.12
controllers output, which is split into two parts affecting both inlet and outlet
valves simultaneously as shown in Fig. 7.12. Suppose the controller output
increases from 6 psi to 9 psi, valve V2 continues to open while the valve V1
remains completely open, thus resulting in reduction of pressure in the vessel.
For very large increase in the pressure vessel, say if the controller output
increases beyond 9 psi the control valve V2 is completely opened while valve V1
starts closing, till the pressure in the vessel reaches the desired value. This
operation of split range control can be clearly understood from Fig. 7.13.
Valve stem
position
Open
V2 V1
Closed Psi
3 6 9 15 Controller output
FIG. 7.13
Controllers output signal (psi) Valve V1 stem position Valve V2 stem position
3 psi Open Closed
6 psi Open Partially open
9 psi Open Open
12 psi Partially closed Open
15 psi Closed Open
5TH PROOF 4-9-09
Soft constraint
Over ride set point
Variable Normal
operating range
time (t)
FIG. 7.14
rated capacity. The normal control loop is the pressure control loop on the steam
discharge line. The over ride control loop is formulated at the level side of the
boiler drum and the outputs of normal controller and over ride controller are fed
to low signal selector (LSS).
Whenever the liquid level falls below the allowable limit, the LSS switches the
control action from pressure control to level control and closes the valve on the
discharge line.
PT
LT LC LSS PC
Water
flow
Normal
Steam drum Over ride controller
control loop loop
FIG. 7.15
QUESTION SET
1. Discuss about the merits and demerits of feed forward control.
2. What are the relative advantages and disadvantages of feed forward
control over feed back control configuration?
3. Explain how the inherent interaction existing in a process can be used for
better control using cascade loops.
4. Explain with suitable examples about ratio control systems.
5. What are over rides? Explain the concept of over ride control
configuration.
6. Explain with an example the split range control configuration.
7. Which control configuration is suitable for a process that has more than
one controlled variable and single manipulated variable? Explain.
8. Justify how ratio control can be useful in blending applications.
9. What is selective control loop configuration?
10. It is required to monitor one controlled variable in a process with more
than one manipulated variables. Explain the relevant control loop
configuration.
Index