Instrumentation & Process Control
Instrumentation & Process Control
PROCESS CONTROL
For
CHEMICAL ENGINEERING
INSTRUMENTATION AND PROCESS
CONTROL
SYLLABUS
Measurement of process variables; sensors, transducers and their dynamics, process
modeling and linearization, transfer functions and dynamic responses of various
systems, systems with inverse response, process reaction curve, controller modes (P,
PI, and PID); control valves; analysis of closed loop systems including stability,
frequency response, controller tuning, cascade and feed forward control.
ANALYSIS OF GATE PAPERS
3. MATHEMATICAL MODELING
4. TRANSFER FUNCTION
5. FEEDBACK CONTROLLER
7.1 Introduction 71
7.2 Piping and Instrumentation Diagram 71
7.3 Process Flow Diagram 75
7.4 Actuators 76
7.5 Sensors 81
8. ADVANCE CONTROL
What is a Process?
The word ‘Process’ used popularly in process control and the process industry refers to the
ways and methods applied in changing or refining raw materials into end products suitable
for mankind. The raw materials which can either be in a liquid, gaseous or a mixture of solid
and liquid(slurry) are during processing transferred, measured, mixed, heated or cooled,
filtered, stored, or handled in some other way to produce the end product.
Process industries include the chemical industry, the oil and gas industry, the food and
beverage industry, the pharmaceutical industry, the water treatment industry, and the
power industry.
Process Control refers to the methods used to control process variables during the
manufacture of a product. The manufacturing process for a product is controlled for the
following reasons:
(i) To reduce product variability – Reducing product variability leads to increase in product
quality.
(ii) Increase efficiency – Increase in the efficiency of a manufacturing process will ensure
that minimum resources are required to produce a product which consequently saves
money for the manufacturers.
(iii) Ensure process safety – Every process must be controlled to operate safely. If this does
not happen, the consequences of a run-away process can be catastrophic as witnessed in the
1986 nuclear plant disaster in Chernobyl near pripyat, about 130 km (about 80 mi) north of
Kyiv (now in Ukraine) in the then USSR which exploded and burned.
Process Variable:
The setpoint is the target value of the process variable that is desired to be maintained. For
example, if a process temperature needs to be kept within 5 °C of 100 °C, then the setpointis
100 °C.
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 1
Introduction
Controller:
Manipulated Variable:
The variable to be manipulated, in order to have control over the process variable (PV), is
called the manipulated variable or MV. If we control a particular flow for instance, we
manipulate a valve to control the flow. Here, the valve position is called the manipulated
variable and the measured flow becomes the process variable. The diagram below illustrates
the relationship between setpoint, measured variable and a controller:
Error
Error is the difference between the measured variable and the setpoint and can be either
positive or negative. The objective of any control scheme is to minimize or eliminate error.
Therefore, it is imperative that the concept of error be well understood. Any error can be
seen as having three major components. These include – the magnitude, duration and rate of
change of the error:
(i) Error magnitude - The magnitude of the error is simply the deviation between the values
of the setpoint and the process variable. The magnitude of error at anypoint in time
compared to the previous error provides the basis for determining the change in error. The
change in error is also an important value.
(ii) Duration of error - Duration refers to the length of time that an error condition has
existed.
(iii) Rate of change of error - The rate of change of error is shown by the slope of the error
plot shown below:
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 2
Introduction
Offset
Offset is a sustained deviation of the process variable from the setpoint. For example, if in a
pressure control loop, the control system held the process pressure at 10.5 bar consistently,
even though the setpoint is 10 bar, then an offset of 0.5 bar exists.
Load Disturbance
A load disturbance is an undesired change in one of the factors that can affect the process
variable. For example, in a temperature control loop adding cold process fluid to a vessel we
are maintaining the temperature at a constant value would be a load disturbance because it
would lower the temperature of the process fluid.
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 3
2 LAPLACE TRANSFORM
Most control system analysis and design techniques are based on linear systems theory.
Although we could develop these procedures using the state space models, it is generally
easier to work with transfer functions. Basically, transfer functions allow us to make
algebraic manipulations rather than working directly with linear differential equations
(state space models). To create transfer functions, we need the notion of the Laplace
transform.
The Laplace transform of a time-domain function, f(t), is represented by L[f(t)] and is defined
as
∞
L f ( t ) = F ( s ) ∫ f ( t ) e − st dt
0
The Laplace transform is a linear operation, so the Laplace transform of a constant (C)
multiplying a time-domain function is just that constant times the Laplace transform of the
function,
∞
f ( t ) C=
L C= f ( s ) C ∫ f ( t ) e − st dt
0
The Laplace transforms of a few common time-domain functions are shown next.
Exponential Function
Exponential functions appear often in the solution of linear differential equations. Here
f ( t ) = e − at
The transform is defined for t > 0 (we also use the identity that ex + y
= ex ey )
∞ ∞
L f (=
t ) L e=
− at
(s)
F= ∫e
− at − st
e=
dt ∫e
−( s + a ) t
dt
0 0
1 −( s + a ) t ∞ 1 1
L f ( t ) =− e =− [0 − 1] =
s+a 0 s+a s+a
So we now have the following relationship:
1
L e − at =
s+a
Derivatives
This will be important in transforming the derivative term in a dynamic equation to the
Laplace domain (using integration by parts),
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 4
Laplace Transform
df ( t ) ∞ df ( t ) − st ∞
∞
=
L ∫ = e dt
e − st
f ( t ) 0 ∫ f ( t ) se dt
+ − st
dt 0 dt 0
df ( t ) ∞
L =
0 − f ( 0 ) ∫ f ( t ) e dt =
+ s − st
−f ( 0 ) + sF ( s )
dt 0
so we can write
df ( t )
L = sF ( s ) − f ( 0 )
dt
For an nth derivative, we can derive
dnf ( t ) n
s F ( s ) − s f ( 0 ) − s f ' ( 0 ) − .... − f ( 0)
n −1 n −2 ( n −1)
L =
dt
n initial conditions are needed: f ( 0 ) ,..., f ( 0)
( n – 1)
One reason for using deviation variables is that all of the initial condition terms in Equation
are 0, if the system is initially at steady-state.
Time delays often occur owing to fluid transport through pipes, or measurement sample
delays. Here we use q to represent the time delay. If f(t) represents a particular function of
time, then f(t – q) represents the value of the function q time units in the past.
∞ ∞ ∞
L f ( t − =
θ ) ∫ f ( t − θ) e =
dt ∫ f ( t − θ ) e
− st − s( t −θ+θ )
=
dt ∫ f ( t − θ) e
− s( t −θ ) − sθ
e dt
0 0 0
∞
t − θ ) e − sθ ∫ f ( t − θ ) e − s( t −θ) d ( t − θ )
L f (=
0
We can use a change of variables, t* = t – q, to integrate the function. Notice that the lower
limit of integration does not change, because the function is defined as f(t) = 0 for t < 0.
∞
L f ( t − θ ) e −=
= sθ
∫ f ( t *) e dt * e F ( s )
− st* − sθ
So the Laplace transform of a function with a time delay (q) is simply e –qs times the Laplace
transform of the non-delayed function.
Step Functions
Step functions are used to simulate the sudden change in an input variable (say a flow rate
being rapidly changed from one value to another). A step function is discontinuous at t = 0.
A "unit" step function is defined as
𝑓𝑓(𝑡𝑡) = 1 𝑓𝑓𝑓𝑓𝑓𝑓 𝑡𝑡 ≥ 0
𝑓𝑓(𝑡𝑡) = 1 𝑓𝑓𝑓𝑓𝑓𝑓 𝑡𝑡 ≥ 0
and using the definition of the Laplace transform,
∞ 1 1 1
𝐿𝐿[𝑓𝑓(𝑡𝑡)] = ∫0 𝑒𝑒 −𝑠𝑠𝑠𝑠 𝑑𝑑𝑑𝑑 = [𝑒𝑒 −𝑠𝑠𝑠𝑠 ]∞
0 = − [ 0 − 1] =
𝑠𝑠 𝑠𝑠 𝑠𝑠
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 5
Laplace Transform
so
1
𝐿𝐿[1] =
𝑠𝑠
Similarly, the Laplace transform of a constant, C, is
𝐶𝐶
𝐿𝐿[𝐶𝐶 ] =
𝑠𝑠
Pulse
Impulse
An impulse function can viewed as a pulse function, where the pulse period is decreased
while maintaining the pulse area, as shown in Figure. In the limit, as tp approaches 0, the
pulse function becomes (using L'Hopital's rule)
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 6
Laplace Transform
1 1
𝐿𝐿[𝑒𝑒 −𝑎𝑎𝑎𝑎 ] = 𝑎𝑎𝑎𝑎𝑎𝑎 𝐿𝐿−1 � � = 𝑒𝑒 −𝑎𝑎𝑎𝑎
𝑠𝑠 + 𝑎𝑎 𝑠𝑠 + 𝑎𝑎
The following theorems are very useful for determining limiting values in dynamics and
control studies. The long-term behavior of a time-domain function can be found by analyzing
the Laplace domain behavior in the limit as the s variable approaches zero. The initial value
of a time-domain function can be found by analyzing the Laplace domain behavior in the limit
as s approaches infinity.
The final-value theorem is,
lim 𝑓𝑓(𝑡𝑡) = lim[𝑠𝑠𝑠𝑠(𝑠𝑠)]
𝑡𝑡→∞ 𝑠𝑠→0
The initial-value theorem is,
lim 𝑓𝑓(𝑡𝑡) = lim [𝑠𝑠𝑠𝑠(𝑠𝑠)]
𝑡𝑡→0 𝑠𝑠→∞
It should be noted that these theorems only hold for stable systems.
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 7
Laplace Transform
Question: Determine the initial value of the time-domain response of the following equation
using initial-value theorem
2𝑠𝑠 + 1
𝑌𝑌(𝑠𝑠) =
(𝑠𝑠 + 1 + 𝑗𝑗)(𝑠𝑠 + 1 − 𝑗𝑗)
Solution: Solution of the above equation is,
𝑦𝑦(𝑡𝑡) = 2𝑒𝑒 −𝑡𝑡 cos 𝑡𝑡 + 𝑒𝑒 −𝑡𝑡 sin 𝑡𝑡
Applying initial value theorem,
𝑠𝑠(2𝑠𝑠 + 1)
lim =2
𝑠𝑠→∞ (𝑠𝑠 + 1 + 𝑗𝑗)(𝑠𝑠 + 1 − 𝑗𝑗)
Question: Solve the following linear differential equation using Laplace Transform:
d2 y dy dy
2
+3 − y = 5t ; = 0 ; y t =0 = 2
dt dt dt t =0
Solution:
d2 y dy
2
+3 − y = 5t
dt dt
Taking Laplace Transform of the equation we obtain,
Hence,
or,
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 8
Laplace Transform
Hence,
or,
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 9
3 MATHEMATICAL MODELLING
The state variables are a set of fundamental dependent quantities whose values describe the
natural state of a given process. The state equations are set of differential equations which
describe the progression of the states with time. The state variables are primarily the
fundamental quantities of a process viz., mass, energy and momentum, whereas the state
equations are generated out of conservation principle involving these fundamental
quantities.
Let S be the fundamental quantity. The state equation can be written as:
Rate of accumulation of S within the system
= Flow rate of S coming to the system
− Flow rate of S going of the system
+ Rate of s generated within the system
− Rate of S consumed within the system
If there are N no. of streams going inside a process and M no. of streams coming out of a
process, the mass balance equations for the process can be written as:
(1)
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 10
Mathematical Modeling
Where ρ and V represent the density and volume of the material inside the system
and ρ(i/ j) and F(i/ j) are the density and flow rate of the incoming/outgoing streams
respectively. In case of a chemical reaction, let CA and CA(i/ j) be the concentrations of
component A of material inside the process and the streams respectively and r be the rate of
the reaction. The component balance equations can be written as
(2)
Let h and h (i/ j) be the enthalpies of material inside the process and incoming/outgoing
streams respectively, Q be the heat supplied to the system (or heat removed from the
system), W is the work done on (or by) the system. Then the energy (E) balance equation for
such system is written as
(3)
As the chemical process plants are usually static, momentum balance equations are usually
not required for such cases.
In addition to the above state equations, a few algebraic equations are also useful for
modeling of a chemical process. Examples of such equations are:
Transport rate equation: (4a)
Consider the stirred tank heater in the Fig 3.1 as shown below. The question is what would
change in case a change is occurred in the input condition (either in the manipulated variable
or the disturbance). It is evident that inlet flow rate and its temperature are the input
condition which can undergo a change and in such situation the mass and energy content
(state variables) of the tank would show a progression. In normal situation, flow rate or
temperature of an inlet flow does not have a potential to displace the tank physically from
its normal position. Hence, there is no scope of progression of momentum of the tank. In
other words, one need not carry out momentum balance operation on this process, rather
mass and energy balance operation would suffice.
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 11
Mathematical Modeling
Let us now apply the material balance and energy balance operation on this process that
would yield the following two equations:
Rate of accumulation of heat = rate of heat in - rate of heat out + rate of heat supplied
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 12
Mathematical Modeling
If we assume the density and specific heat of water to be constant and the reference
temperature to be zero, then the energy balance equation would take the form as
(6)
Equation 6 and 7 represent the mathematical model of the stirred tank heater.
The degree of freedom is defined by the total number of independent variables that must be
specified in order to define the system completely. In other words, it is the minimum number
of independent variables required to be specified so that the model equation(s) can be
solved.
Let us analyze the case of the stirred tank heater:
Total no. of variables : 5 ( h, T, F , T , Q ) i i
The model consists of two ordinary differential equations. Simulation of the model
(i.e. solving these model equations) with an “input condition” yields the data along with some
percentage of error, which are otherwise expected, from the real-time process operation at
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 13
Mathematical Modeling
same “input condition”. Although equations are simple first order ODEs, the nonlinear terms
of those model equations may be hindrance against getting analytical solution of the model.
Numerical solution, on the other hand, is a popular approach for simulation of such model
equations. Various simulation software are available for easy use, e.g. MATLAB® and
SIMULINK® to name a few. Following subsection describes the model of stirred tank heater
developed in SIMULINK® domain.
In order to solve the model equation one needs to specify the model parameters and satisfy
the degree of freedom. Let us consider the following numerical values for the said
parameters and variables: A = 10, c = 1;c p = 1, ρ = 1, Fi = 1, Ti = 29, Q = 1 . At steady state, left
hand side of equations (5) and (6) would be zero.
(7)
(8)
Solution of equations (7) and (8) yields the steady state values of the state
variables: h s = 1 and Ts = 30 . The values of h and T would stay at 1 and 30 respectively as
long as there is no change in the input condition.
Let us analyze what would happen if the temperature of inlet flow increases from 29 to 30.
Solution of equations (7) and (8) with Ti = 30 yields a new steady state condition of state
variables as h s = 1 and Ts = 31 . That means a disturbance in Ti would not change the height
of the liquid inside the tank however it would increase the temperature of the liquid to an
increased new steady state. The transition of Ts from 30 to 31 can be observed (Fig 3.3) by
simulating the model in SIMULINK® domain.
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 14
Mathematical Modeling
It is observed from the figure that the process takes 70 s to reach the new steady state of the
liquid temperature inside the stirred tank heater.
From Eq (7), we observe that the steady state value of h = 1 when the inlet flow rate is
maintained at Fi = 1 .
Case I: Inlet flow rate undergoes a step change and new value is Fi = 2 . As a result
the h would attain a new steady state value of 4. The gain of the system is defined as the ratio
∆h s 4 − 1
of change in output to change in input. Hence, gain (K)= = = 3.
∆Fis 2 − 1
Case II: Inlet flow rate undergoes a step change and new value is Fi = 3 . As a result
∆h s 9 − 1
the h would attain a new steady state value of 9. Hence, gain (K)= = = 4.
∆Fis 3 − 1
It is observed that the gain of the system is not constant. The gain changes with various input
conditions of the system. This is the identity of a nonlinear process. In general, if a system
follows the principle of superposition then it is termed as linear process. Usually, a nonlinear
process is identified with existence of a nonlinear term in its model such as logarithmic,
power, exponential, product etc. Existence of square root term with h yields nonlinearity in
the stirred tank heater process.
Most of significant developments in the control theories have occurred for linear processes
as easy analytical solutions are available for them. Unfortunately, most of the chemical
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 15
Mathematical Modeling
processes are nonlinear in nature. Fig II.4 pictorially describes the effect of linearization of a
nonlinear process.
Suppose the blue line indicates the true dynamics of a nonlinear process. A steady state
point Fis , h s is identified on this line that represents the nominal point of operation of the
process. If the process is linearized in the neighborhood of the nominal operating point then
the resulting gradient can be indicated by the red line. It is evident from the figure that there
lies a considerable portion of domain in the neighborhood of the nominal point of operation
over which the dynamics of both nonlinear and linearized process match. Hence, it can be
argued that the basic aim of a control system is to maintain the operation of the process at a
predefined nominal point. Even if a disturbance takes the process away from this nominal
point, an ideal control system should ensure that the departure of the process from its
nominal point of operation remains within the limit of such neighborhood and the process
is eventually brought back at its nominal operating point within reasonable time. A
linearized form of the model is likely to cover sufficient domain around the nominal
operating point so that its dynamics reasonably matches with that of the actual nonlinear
model. In such situation, it is customary to design a linear controller on the basis of linearized
model and apply it on the nonlinear process.
The linearization of model is carried out by using Taylor series expansion. The expansion of
a function f ( x ) around x 0 is given as
(9)
When the value of x is very close to x 0 , then the power terms of ( x − x 0 ) are very small and
hence can be neglected. The resulting function will have the form:
(10a)
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 16
Mathematical Modeling
Consider the equation (5) of the stirred tank heater. Only nonlinear term in the equation
is h . Linearization of this nonlinear term would yield the following:
(11)
Or
(12)
As the aim of the control system is to maintain the process at a nominal operating point, the
deviation or dislocation of the state from its nominal point is a crucial variable to observe.
The deviation variable of a state describes the departure of the state from its nominal point
of operation.
Let us consider the linearized model of the stirred tank heater as given by Eq.12. At steady
state,
(13)
Or,
(14)
Eq. (14) represents the deviation form of the model. The variables h and Fi represent the
deviation variables of the height of liquid in the stirred tank heater and that of the inlet flow
rate.
Question:Consider two tanks in series as shown in the figure. F, h and Arepresent the liquid
flow rate, height of liquid and cross sectional area of tank respectively. Assume that the flow
rate of an effluent stream from a tank is proportional to the cubic root of the liquid height.
(a) Starting from usual material balance, develop a time domain mathematical model for this
process.
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 17
Mathematical Modeling
(b) Identify all the state(s), manipulated input(s), disturbance input(s) and process
output(s).
(c) It is intended to control height of liquid in the second tank. Develop a Laplace transfer
function based model relating input-output variables.
Solution:It is understood that the outlet flow of a tank is proportional to the cubit root of the
height of the liquid in that tank. Hence the outlet flows of the tanks are written as:
(a) The material balance on the tanks would yield the following equations:
Hence the time domain model for the whole process can be written as:
(b) The process has two states viz., h1 and h 2 , two inputs viz. F1 and F4 and two
outputs viz., h1 and h 2 . One control valve is attached with input F1 and hence it can be
manipulated as and when desired. Whereas, F4 is not associated with any such tool and
hence it is a source of potential disturbance in the process. Thus F1 is the manipulated input
and F4 is the disturbance input.
(c) As h 2 is the output to be controlled and F1 is the only manipulated variable available
h2 (s)
with the process, the desired transfer function is .
F1 ( s )
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 18
Mathematical Modeling
The only nonlinear part in the above model is with F2 and F3 . Linearizing these portions, we
obtain
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 19
Mathematical Modeling
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 20
4 TRANSFER FUNCTION
A transfer function relates inputs to outputs in the Laplace domain. In particular, the Laplace
domain relationship between a manipulated input and output is called a process transfer
function.
Consider an nth order differential equation,
dn y d n −1 y d n −1u
+ a n −1 + =
... + a y b n −1 + ... + b 0 u
dt n −1 dt n −1
0
dt n
Since we are assuming that the model is based on deviation variables, and that the system is
initially at steady state, the initial conditions are
d n −1 y d n −2 u ( 0 )
= =
... y ( 0 )
= = ...= u ( 0 )= 0
dt n −1 t =0 dt n − 2 t =0
Taking the Laplace transform of each term,
dn y n d n −1 y
L = s Y ()
s − s n −1
y ( )
0 − ... − n −1
= sn Y ( s )
dt n dt t =0
.
.
.
L [ y] = Y ( s )
d n −1u n −1 d n −2 u
L = n −1
s U ( s ) − s u ( 0 ) − ... − n −=
n −2
2
s n −1U ( s )
dt dt t =0
.
.
.
L [u] = U (s)
b n −1s n −1 + ... + b 0
Y (s) = n .U ( s )
s + a n −1s n −1 + ... + a 0
Solving for Y(s), we find
The ratio of polynomials is called the transfer function. When it relates a manipulated input
to an output it is commonly called a process transfer function. In general, we will use gp(s) to
represent the process transfer function.
Y (s) = gp (s) U (s)
In this case the process transfer function is
b s n −1 + ... + b 0
g p ( s ) = n n −1 n −1 .U ( s )
s + a n −1s + ... + a 0
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 21
Transfer Function
The roots of the numerator polynomial are known as zeros, and the roots of the denominator
polynomial are call poles. We have used capital letters to distinguish Laplace domain
variables from the time-domain variables. In the rest of this text we generally use lowercase
letters for all input and output variables. If the argument is s, then we assume that we are
referring to the Laplace domain.
Transfer functions are often used in block diagrams. For example, the relationship between
an input and output is shown as
In the rest of this chapter we study the dynamic behavior of some commonly used transfer
functions. Our focus is on Step Responses, since process engineers often apply step changes
to understand dynamic behavior.
Many chemical processes can be modeled as first-order systems. The differential equation
for a linear first-order process is often written in the following form:
dy
τp +y= kpu
dt
This can also be written as
dy 1 kp
=− y+ u
dt τp τp
where the parameters (tp and kp) and variables (y and u) have the following names: tp is the
process time constant (units of time), kp the process gain (units of output/input), y the
output variable, and u the input variable. Taking the Laplace transform of each term (notice
that we are now using lower-case variables to represent the Laplace domain input and
output), and assuming that the initial condition is y (0) = 0,
dy dy
L τp = τp L = τp sy ( s ) − y ( 0 ) =
τp .sy ( s )
dt dt
L [ y] = y ( s )
=
L k p u ( t ) k= u ( t ) k p .u ( s )
p .L
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 22
Transfer Function
Step Response
Consider the case where the output is initially zero (steady state in deviation variable form),
and the input is suddenly step changed by an amount Δu. The Laplace transform of the input
is
∆u
u (s) =
s
So Equation can be written
k p ∆u
y (s ) = .
τps + 1 s
Using a partial fraction expansion and inverting to the time domain,
y ( t ) =k p ∆u 1 − e p
−t/τ
( )
Here the notion of a process gain is clear. After a substantial amount of time (t>> tp), we find,
from Equation
y ( t → ∞ )= k p ∆u
That is,
y ( t → ∞)
kp =
∆u
and, since y(0) = 0, we can think of y ( t → ∞ ) as Dy, so
∆y long-term output change
=
kp =
∆u input change
We can think of the process time constant as the amount of time it takes for 63.2% of the
= tp ,
ultimate output change to occur, since when t
y ( t =τp ) =k p ∆u 1 − e ( −τp / τp
) =k ∆u (1 − e ) =0.632k ∆u =0.632∆y
p
−1
p
Impulse Response
Consider now an impulse input of magnitude P, which has units of the input*time; if the input
is a volumetric flow rate (volume/time), then the impulse input is a volume. The output
response is
kp
y (s) = .P
τps + 1
You should find that the time domain solution is
Pk p − t / τp
y(t) = e
τp
which has an immediate response of Pk p / t p followed by a first-order decay with time.
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 23
Transfer Function
Integrating System
We found in the previous chapter that material balances on liquid surge vessels or gas drums
often yielded models with the following form:
dy
= ku
dt
In the Laplace domain, this is
k
y (s) = u (s)
s
Consider an integrating process initially at steady state, with y ( 0 ) = 0 .
Step Response
Impulse Response
Second-order models arise from systems that are modeled with two differential equations
(two states). In this section we separately consider transfer functions that do not have
"numerator" dynamics and those that do.
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 24
Transfer Function
where (obviously a 0 ≠ 0 )
τ2 a 2 / a 0 , 2ζτ = α1 / a 0 , k = b 0 / a 0
and the parameters are as follows: kis the gain (units of output/input), zthe damping factor
(dimensionless), and t the natural period (units of time).
The second-order model shown in Equation generally arises by changing a set of two first-
order equations (state-space model) to a single second-order equation. For a given second-
order ODE, there are an infinite number of sets of two first-order (state-space) models that
are equivalent.
Taking the Laplace transform
τ2 s 2 y ( s ) − sy ( 0 ) − y ( 0 ) + 2ζτ sy ( s ) − y ( 0 ) + y ( s ) =
k.u ( s )
where y(s) indicates the Laplace transformed variable.
Assuming initial conditions are zero, that is dy / dt=
t =0
( 0 ) 0 , we find
y=
k
y (s) =
.u ( s )
τ s + 2ζτs + 1 2 2
The characteristic equation of the second-order transfer function is τ2s 2 + 2ζτs + 1 . We can
find the roots (known as the poles) by using the quadratic formula
−b ± ∆
α, β = , where ∆= b 2 − 4ac
2a
which yields the following values for the roots:
ζ ζ2 −1 ζ ζ2 −1
p1 =− + , p 2 =− −
τ τ τ τ
The following analysis assumes that z >0 and t > 0. This implies that the real portions
of p1 and p2 are negative and, therefore, the system is stable.
Step Responses
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 25
Transfer Function
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 26
Transfer Function
polynomials in the transfer function. If the relative order is 1, then output response has a
nonzero slope at the time of the step input; the step response of a system with a relative
order greater than 1 has a zero slope at the time of the step input.
A number of insights can be obtained from Fig 4.2 and from an analysis of the step response
equations. For z < 1 , the ratio of the imaginary portion to the real portion of the pole is
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 27
Transfer Function
Im 1− ζ2
=
Re ζ
As the imaginary/real ratio gets larger, the response becomes more oscillatory (z becomes
smaller). We also notice that a decreasing ratio corresponds to a larger negative value for
the real portion. As the real portion becomes larger in magnitude (more negative), the
response becomes faster.
The following common measures of underdamped second-order step responses are shown
in Fig. 4.3, and defined below: (1) rise time, (2) time to first peak, (3) overshoot, (4) decay
ratio, and (5) period of oscillation.
Rise time is the amount of time it takes to first reach the new steady-state value. Time to first
peak is the time required to reach the first peak. Overshoot is the distance between the first
peak and the new steady state. This is usually expressed as the overshoot ratio, as shown
in Fig. 4.3. Decay ratio is a measure of how rapidly the oscillations are decreasing. A b/a ratio
of 1/4 is commonly called "quarter wave damping." Period of oscillation is the time between
successive peaks.
The previous discussion involved pure second-order systems, where the relative
order (difference between the denominator and numerator polynomial orders) was two.
Consider now a second-order system with numerator dynamics with the gain/time constant
form
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 28
Transfer Function
k p ( τn s + 1)
gp (s) =
( τ1s + 1)( τ2s + 1)
which is relative order one. The reader should show that the pole-zero form is
k pz ( s − z1 )
gp (s) =
( s − p1 )( s − p 2 )
where the parameters are
kτn 1 1 −1
k pz = , p1 =
− , p2 = − , z1 =
τ1τ2 τ1 τ2 τn
The gain/time constant form has the following time-domain response to a step input:
τ −τ τ −τ
y ( t ) =k p ∆u. 1 + n 1 e − t / τ1 + n 1 e − t / τ2
τ1 − τ2 τ1 − τ2
The reader should show that, if τn =τ2 , the response is the same as a first-order process.
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 29
Transfer Function
Notice also that a numerator time constant that is greater than the denominator time
constant causes overshoot before settling to the final steady state. Also notice that the
inverse response becomes "deeper" as the numerator zero ( −1/ τn ) approaches a value of 0
from the positive side.
Lead-Lag Behavior
Lead-lag transfer functions have the same order numerator polynomial as denominator
polynomial. This occurs when the input has a direct effect on the output variable. In terms of
the state space model, this means that D ≠ 0.
Consider a lead-lag transfer function where the numerator and denominator polynomials
are first order:
τ s +1
gp (s) = k p . n
τp s + 1
For a step input of magnitude ∆u , the output response is
τ −t/τ
y ( t ) = k p ∆u 1 − 1 − n e p
τp
A dimensionless output can be defined by dividing Equation by k p ∆u , and τ / τp is a natural
dimensionless time. The responses to a step input at t = 0 are shown in figure. Notice that
there is an immediate response that is equal to the τn / τp ratio. This can also be found by
applying the initial value theorem to Equation for a step input change.
There are a number of different ways to represent process transfer functions. The
"polynomial" form is
b s m + b m −1s m −1 + ... + b1s + b 0
g (s) = m n
a n s + a n −1s n −1 + ... + a1s + a 0
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 30
Transfer Function
The values of s that cause the numerator to equal zero are known as the "zeros" of the
transfer function. The values of s that cause the denominator of to equal zero are known as
the "poles" of the transfer function.
The "pole-zero" form is
k pz ( s − z1 )( s − z 2 ) ... ( s − z m )
g (s) =
( s − p1 )( s − p 2 ) ... ( s − p n )
and complex poles (or zeros) must occur in complex conjugate pairs.
The "gain-time constant" form is the one that we use most often for control system design.
k p ( τn1s + 1)( τn 2s + 1) ... ( τnms + 1)
g (s) =
( τp1s + 1)( τp2s + 1) ... ( τpns + 1)
where τni is a numerator time constant and τpi is a denominator time constant. This form is
normally used when the roots (poles) of the denominator polynomial are real.
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 31
Transfer Function
As poles move further to the left they yield a faster response, and increasing the magnitude
of the imaginary portion makes the response more oscillatory. This behavior is summarized
in Fig. 4.7. Recall also that a process with a pole at the origin (and none in the RHP) is known
as an integrating system; that is, the system never settles to a steady state when a step input
change is made.
Figure: 4.7 Effect of pole-zero location on dynamic behavior (x, poles; o, zero). As poles
become more negative, the response is faster. As the imaginary/real ratio increases, the
response becomes more oscillatory.
Multiple RHP zeros cause multiple "changes in direction"; for example, with two RHP zeros,
the step response, initially going in one direction, switches direction, then switches back to
the initial direction.
Many processes have a delayed response to a process input, either due to transport lags
(such as flow through pipes) or high-order effects. In Section 3.3.4 we found that a time-delay
of q time units had a transfer function of e −θs .
A first-order process combined with a time-delay has the transfer function
k p e −θs
gp (s) =
τps + 1
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 32
Transfer Function
Consider a process with k p = 1 , τp =10 , and q = 5 . A unit step input change at t = 0 yields
the response shown in Figure 4.8. We see that the time-delay shifts the response of the
output.
Question: The initial value ( t = 0+ ) of the unit step response of the transfer fucntion
( s + 1) / ( 2s + 1) is?
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 33
Transfer Function
Solution:
1 𝑠𝑠 + 1
lim 𝑓𝑓(𝑡𝑡) = lim [𝑠𝑠𝑠𝑠(𝑠𝑠)] = lim 𝑠𝑠
𝑡𝑡→0 𝑠𝑠→∞ 𝑠𝑠→∞ 𝑠𝑠 2𝑠𝑠 + 1
By L-Hospital rule:
(𝑠𝑠 + 1)
lim [𝑠𝑠𝑠𝑠(𝑠𝑠)] = = 1/2
𝑠𝑠→∞ (2𝑠𝑠 + 1)
Question: For an input forcing function, X ( t ) = 2t 2 , the Laplace transform of this function
is ?
Solution:
𝑛𝑛!
𝐿𝐿(𝑡𝑡 𝑛𝑛 ) = 𝑛𝑛+1
𝑠𝑠
4
= 3
𝑠𝑠
Question: A process is initially at steady state with its output y = 1 for an input u = 1 . The
input is suddenly changed to 2 at t = 0 . The output response is y ( t )= 1 + 2t . The transfer
function of the process is?
Solution:
𝑌𝑌(𝑠𝑠)
𝐺𝐺𝑝𝑝 (𝑠𝑠) =
𝑈𝑈(𝑠𝑠)
Given: y ( t )= 1 + 2t and u ( t ) = 2 − 1 = 1
(s)
Y= L y (=
t ) =
1 / s + 2 / s2 ; and U (s) u ( t )
L= 1 / s . Therefore
1 2
+ 2 2
𝐺𝐺𝑝𝑝 (𝑠𝑠) = 𝑠𝑠 𝑠𝑠 = 1 +
1 𝑠𝑠
𝑠𝑠
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 34
Transfer Function
Question: The transfer function of a pure dead time system with dead time τd is?
1
(a) (b) 𝜏𝜏𝑑𝑑 𝑠𝑠 + 1 (c) 𝑒𝑒 −𝜏𝜏𝑑𝑑𝑠𝑠 (d) 𝑒𝑒 𝜏𝜏𝑑𝑑𝑠𝑠
𝜏𝜏𝑑𝑑 𝑠𝑠+1
Answer: (c)
Question: A rectangular tank is filled with a valve at the bottom and is used for storing a
liquid. The area of cross-section of the tank is 10 m2 and the flow resistance of the valve
(assumed constant) is 0.1 s / m2 . The time constant of the tank will be?
(a) 1 (b) 100 (c) 10.1 (d) 9.9
Question: A certain thermocouple has a specific time constant of 2 sec. If the process
temperature changes abruptly from 800°C to 900°C , the temperature reading in an indicator
attached to the thermocouple after 6 sec will be approximately?
(a) 800°C (b) 900°C (c) 890°C (d) 895°C
Answer: (a)
Question: The second order system with the transfer function 4 / ( s2 + 2s + 4 ) has a damping
ratio of
(a) 2.0 (b) 0.5 (c) 1.0 (d) 4.0
Solution: (b)
𝑠𝑠 2 + 2𝑠𝑠 + 4 = 𝜏𝜏 2 𝑠𝑠 2 + 2𝜁𝜁𝜁𝜁 + 1
1 1 2
𝜏𝜏 = ; 𝜏𝜏 = ; 2𝜁𝜁𝜁𝜁 = = 1/2
4 2 4
2 1
𝜁𝜁 = = = 0.5
2×2 2
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 35
5 FEEDBACK CONTROLLER
Consider a chemical process plant shown in the Fig 5.1. It has an output y that needs to be
maintained at a desired setpoint ysp . The process is subjected to unavoidable disturbance
(load) d . One input is available in the process that can be manipulated in order to steer the
process to its desired setpoint . Hence, in the event of process disturbance drifting the
process output away from its desired setpoint, the objective of feedback control is to steer
the process output at its setpoint by manipulating the process input.
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 36
Feedback Controller
manipulation of u . However, a setpoint ( ysp ) for the process output is available for manual
manipulation in case the process objective demands so.
Fig. 5.2 demonstrates some typical feedback control systems which are often encountered in
a chemical process plant. This is an example of a distillation column. There are six control
loops in this example, one each for pressure, temperature, flow, level, reflux and
composition. The measurements are done through an appropriate measuring device cum
transmitter (denoted by PT, TT, FT, LT, RT and CT). The measured output is compared with
the respective setpoints (denoted by Psp , Tsp , Fsp , Lsp , R sp and Csp ) and errors are fed to
the respective controllers (denoted by PC, TC, FC, LC , RC and CC) and control actions are
subsequently initiated by their corresponding valves.
The feed is pumped into the column whose flow rate and composition are controlled. Feed
mainly comes through a major supply line where the feed composition may vary due to
various factors which are beyond the scope of the operator of the distillation column. An
auxiliary supply line with guaranteed purity of material can be mixed with the main supply
through a regulated dose in order to maintain the composition of the feed entering the
column. A flow controller at the entry point of the column maintains the feed flow rate. The
pressure of the top section of the column is maintained by regulating the flow of vapor going
out of the column. The vapor is condensed and collected in a tank. The outlet temperature of
distillate coming out the condenser (heat exchanger) is controlled by regulating the flow rate
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 37
Feedback Controller
of cooling water used in the heat exchanger. The liquid level inside the storage tank in the
downstream of condenser is controlled with a level controller by regulating the tank outlet.
The L/D ratio of a distillation column is very important and it needs to be controlled by a
ratio controller. The ratio transmitter is essentially a flow measuring device that sends the
ratio of the two flow rates rather than the flow rates themselves. In all the above examples,
only one type of final control element has been demonstrated viz .control valve.
Nevertheless, there are various other final control elements available such as variable speed
metering pump, heating regulators etc.
Note: The composition transmitter/controller is usually a high value delicate instrument
which is not used in a rugged process plant environment. On-line gas chromatograph is an
example of such sophisticated instrument. They are used only in a few critical operations
where tighter regulation of composition is mandatory, such as automated drug delivery.
The above example indicates that the basic hardware components of a feedback control
configuration are the following:
• Process
• Measuring devices
• Transmission lines
• Controllers
• Final Control Elements
5.1.1 Process
Any equipment that serves the targeted physical/chemical operation of the plant is termed
as a process. Reactors, separators, exchangers, pressure vessels, tanks, etc. are examples of
a process. Typically these processes are connected in a logical fashion and the output of one
process becomes input to the other. Any disturbance/malfunction of one process may affect
other processes in the downstream side (and upstream too, in case recycle streams are
used). Detailed discussions on these processes are not within the scope of this course,
however, the modeling techniques and related issues have already been discussed before.
Process variables are primarily pressure, temperature, flow rate, level, composition, etc.
From the process control perspective, it is crucial to study how the changes in one process
variable affect the other, so that an educated measure of control action on one variable can
be taken in order to maintain the other.
The success of any feedback control operation depends largely on accurate measurement of
process variables through appropriate sensors. There are a large number of commercial
sensors available in the market. They differ in their measuring principle(s) and/or their
construction characteristics. Module VII accounts of a few of such measuring instruments.
Details of such devices may be found in technical booklets dedicated for those individual
items.
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 38
Feedback Controller
5.1.3 Controllers
A controller is basically a mathematical function block that reads the error between desired
set point and the measured output and then computes the corrective action for the
manipulated input that would steer process towards the desired set point. There are three
basic types of feedback controllers which are widely used in the industry.
• Proportional (P) controller
• Proportional Integral (PI) controller
• Proportional Integral Derivative (PID) controller
The actuating output of a P controller is proportional to the error between the set point and
process output. Higher the error, higher will be the control action. The control law is given
as:
c ( t ) = K c ε ( t ) + cs (1)
where K c is called the gain of the controller and cs is the bias signal. When error
signal ε ( t ) is zero (i.e., the process output reaches its desired setpoint), the control
signal c ( t ) stabilizes at its bias value cs . The deviation form of actuating signal is
c ( t ) = c ( t ) − cs = K c ε ( t ) (2)
Hence the transfer function of the proportional controller is
c (s) (3)
=Gc = Kc
ε (s)
In industrial lingo, the proportional controller is also termed as “Gain” controller. Equivalent
representation of proportional gain is proportional band. It is the amount of change in error
that will cause the control action to go from full OFF to full ON. The amount of change in error
is calculated as a percentage of full-scale error,
100 (4)
PB ( in% ) =
Kc
e.g., consider a level controller acting on a tank where we measure the level from bottom to
top as 0 to 100%. A control valve on the outlet of the tank maintains the level in the tank. The
PB is defined as the range of level over which the control valve will go from fully
closed to fully open. Let us take the example of a tank whose maximum level is 5 m. Suppose
we decide that if the tank level should fall to 20% (1 m) we want the control valve fully closed
(0% open) and if the tank level rises to 60% (3 m) we want the control valve to be fully open
(100% open). If the tank level is between 20% and 60% we want the control valve to be open
in pro rata basis. This controller would have a PB of 40% (60% - 20%). So if the tank level
were to rise to 2.5 m or 50% of the full tank (75% up the PB), the control valve should be set
to 75% open).
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 39
Feedback Controller
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 40
Feedback Controller
At this point, it is worth explaining the significance of the term reset . Suppose the error
between desired setpoint and process output changes by a constant step of magnitude ε .
The effect of integral term of eqn. (5.6) after every τ1 minutes is given as
1 t ε τ1 ε (8)
K c ∫ ε ( t ) dt = K c ∫ dt = K c τ1 = K c ε
τ1 0 τ1 0 τ1
In other words, the integral action repeats the response of the proportional action every τ1
minutes and “resets” itself for an integral action. Sometimes the controllers are calibrated in
terms of reciprocal of reset time, 1/ τ1 (repeats per minute). This is known as reset rate.
The reset term causes the control action changing as long as there exists a non-zero error in
the system. Often this error cannot be eliminated quickly and given enough time, they
produce larger values for integral terms. Such situation is often observed when the system
undergoes a large change in set point (say a positive change) and the integral term
accumulates a significant error during the rise. This condition is termed as Integral Windup.
The control action in turn keeps on increasing until it reaches the control valve saturation
(i.e. control valve fully open or fully closed). Even if the error changes its sign (as the process
output overshoots the desired set point), this accumulated error has to unwind completely
before control action is reversed. Various measures can be taken to address the issue of
integral windup such as:
• Re-initializing the integral action to a desired value
• Increasing the set point in a suitable ramp (rather than a single step jump)
• Disabling the integral action until the process output enters the controllable region
• Preventing the integral term from accumulating above or below pre-determined bounds
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 41
Feedback Controller
Measurement and/or control signals are carried through various transmission lines. Various
process piping, connection and transmission lines, as per the standard set by International
Society of Automation (ISA), are listed in the following figure.
A heavy solid line represents piping, a thin solid line represents process connections to
instruments, a dashed line represents electrical signals (e.g., 4–20 mA connections), a
slashed line represents pneumatic signal tubes, a line with circles on it represents data links.
Other connection symbols include capillary tubing for filled systems, (e.g., remote diaphragm
seals), hydraulic signal lines, and guided/unguided electromagnetic or sonic signals.
Electric/electromagnetic signals are instantaneous. Unless the process changes very fast
and/or the transmission lines are too long, the dynamic behaviors of
electric/electromagnetic transmissions are also usually ignored.
The hardware component of a control loop that resides between the process and the
controller and implements the control action is called the Final Control Element (FCE). It
receives the control signal from the controller and regulates the value of the manipulated
variable accordingly. An ideal FCE should be an instantaneously operating hardware that
does not induce any time lag between process and controller. However, in reality it is
impossible to find such an instantaneously operating FCE. Nevertheless any FCE, whose time
constant is very small, is considered to be a good hardware for a control system.
Various types of FCEs are available in process industries that are widely used in control
applications. They can be largely classified based on their energy source, viz. Pneumatic,
Hydraulic and Electric, etc. .
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 42
Feedback Controller
Consider the generalized closed loop process given in the Fig 5.5.
Fig. 5.5: Generalized block diagram of feedback control loop of a chemical process
If we assume that the transmission line does not affect the signal flow, dynamics of the
transmission lines can be completely ignored. Hence the following sub-processes will
constitute the overall dynamics of the process:
Process (11a)
Measurement (11b)
Controller (11c)
(11d)
Final Control Element (11e)
Hence using eqs. (5.11 b-e) on eq. (5.11a), following algebraic manipulation will be useful.
(12)
The eq. (5.13) represents the closed loop response of the process. The pictorial
representation is given in the Fig 5.6
Fig. 5.6: Simplified block diagram of feedback control loop of a chemical process
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 43
Feedback Controller
Where
(14)
(15)
G sp and G Load are the closed loop transfer functions (CLTF) of the process. G sp maps the
effect of change of setpoint on the process output whereas G Load maps the effect of change
of load (disturbance) on the process output.
Two types of control problems are encountered with the feedback control systems, viz .servo
and regulatory.
When setpoint of a process undergoes a change while the disturbance affecting the process
remains constant, i.e. d ( s ) = 0 , the objective of the control system would be to steer the
output as close as the setpoint trajectory. In such situation,
y ( s ) = G sp ysp (16)
However, when the setpoint remains constant, i.e. ysp = 0 , while the disturbance forces the
process output to move out of the track of the setpoint, the objective of the control system
would be to reject the effect of disturbance as soon as possible and steer the output back to
the setpoint trajectory. In such situation,
ysp = G Load d ( s ) (17)
Note that the CLTFs depend not only on the process transfer functions, but also on the
transfer functions of measuring element, controller and final control element.
To expedite the construction of overall closed loop transfer function of any feedback loop,
following rules may be applied:
1. The denominator of the overall transfer function is “the product of all transfer functions
in the feedback loop PLUS one” i.e ., {1 + G p ( s ) G f ( s ) G c ( s ) G m ( s )}
2. The numerator of the overall transfer function is “the product of all transfer functions in
the forward path between setpoint and the controlled output (for servo problem) or the load
and the controlled output (for regulatory problem)” i.e., {G p ( s ) G f ( s ) G c ( s )} or {G d ( s )}
5.2.1 Closed loop response of liquid level in a storage tank: A case study
Fig 5.7 represents the closed loop block diagram of a storage tank whose liquid level needs
to be controlled at a predefined setpoint.
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 44
Feedback Controller
Process: The process has two inputs and one output. Input Fi can be manipulated whereas
input Fd is the source of disturbance. Output Fo varies proportionally with the square root
of the height of the liquid in the tank as, Fo ∝ h The material balance around the tank gives
the following model,
dh
A = Fi + Fd − Fo = Fi + Fd − c h (18)
dt
Where, c is the constant of proportion. Linearization of the model around its nominal
operating point {h s , Fi,s , Fd,s } , expressing it in its deviation form and converting it in Laplace
domain yields the following open loop transfer function:
(19)
Measuring device: Differential Pressure (DP) cell is a measuring device that measures
differential pressure (ΔP) between two ends. Output of DP cell maps a second order behavior
with the differential pressure. In this example, DP cell measures the height of liquid level in
the tank by comparing pressure exerted by liquid at one end with that of atmospheric
pressure at the other. The pressure exerted by the liquid on the DP cell is linearly
proportional to height of the liquid. Hence,
K (20)
= hm (s) + h (s)
τ s + 2ξτs + 1
2 2
where h m = measured value of the height of the liquid, h = actual height of the
liquid, K, τ, ξ =constants.
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 45
Feedback Controller
Control Valve: Let us assume that the control valve for this system follows a first order
dynamics, hence
Kv (23)
Fi ( s ) = c (s)
τvs + 1
The following figure is a representation of closed loop control configuration of liquid level in
a tank.
Fig. 5.8: Schematic of closed loop control configuration of liquid level in a storage tank
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 46
Feedback Controller
(27)
The ultimate closed loop response at t → ∞ under P controller never reaches its desired
setpoint. There always remains a discrepancy between the final value of response and the
setpoint. This discrepancy is called the offset. Following derivation gives an account of the
offset shown by first order and second order processes.
The final values of the processes on unit step change in setpoint are
(28)
(29)
Hence,
(30)
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 47
Feedback Controller
The following figures show the offsets exhibited by a first and second order processes. In
both the cases, offset is {1/ (1 + 10 ×1)} =
0.0909
Note that both natural period of oscillation and damping factor decreased in closed loop
response; more the value of controller gain less are the values of natural period of oscillation
and damping factor. As a result the closed loop response becomes more oscillatory than its
open loop counterpart.
We shall find closed loop response of first order process under the influence of an “integral
only” controller.
(31)
It indicates that the order of the dynamics for the closed loop response is increased when
integral control action is used. As a result, the response becomes sluggish.
Following derivation gives an account of the offset shown by first order closed loop process
with a PI controller. The final value of the process on unit step change in set point is
(32)
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 48
Feedback Controller
The offset = 1 - 1 = 0, indicates that the integral control action eliminates any offset. It also
indicates that the increase in controller gain K c and decrease in reset value τ1 reduce the
value of natural period of oscillation and damping factor of closed loop response. Hence, we
can improve the speed of the closed loop response at the expense of higher deviations and
longer oscillations. In other words, high K c and low τ1 makes the closed loop response more
sensitive and proportional action will be stronger. On the other hand, low K c and
high τ1 makes system sluggish but integral action would be stronger (no offset, no
overshoot).
We shall find closed loop response of first order and second order processes under the
influence of a “derivative only” controller.
(33)
(34)
It indicates that the derivative control action does not change the order of the response. It
also indicates that time constant of the first order process is increased by a positive quantity
of K p K c τD , i.e ., the derivative control action makes a process sluggish. On the other hand,
natural period of oscillation of second order process remains same. Damping factor
increases with increase in K c and/or τ1 Increase in damping factor means less sensitivity of
process response and more robustness in the process behavior.
The existence of offset in a proportionally controlled process indicates that the input
manipulation has reached some sort of saturation before the controlled output could even
reach the setpoint. As a result, manipulated input no longer changes even if the error
between setpoint and controlled output is non-zero. The above fact can be demonstrated
mathematically as follows.
The equations of P and PI controllers are:
u= K c ε (35)
t (36)
=u K c ε + ∫ εdt
0
The gradient of the dynamics of above control actions can be derived by simple
differentiation,
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 49
Feedback Controller
du dε (37)
= Kc
dt dt
du dε (38)
= Kc + ε
dt dt
The above equations show that in case of P controller gradient of control action becomes
zero if there is no change in the value of error. In other words, when the error dynamics
reaches a saturation point, the control action reaches saturation too. Mere existence of non-
zero error does not yield any change in control action. On the other hand, gradient of control
action in a PI controller not only depends upon the gradient of error dynamics but also upon
the absolute value of the error at that moment of time. Hence, until the error and its gradient
both are zero, the gradient of control action will not be zero. This fact ensures offset-free
control action by PI controller.
The stability of a process is realized by the location of its poles. Existence of a positive pole
leads to the instability of the process. One of the main purpose of the controller is to ensure
stable closed loop process for the otherwise unstable open loop process.
It is evident from the above equation that one needs to set 3K c > 1 in order to make the
closed loop transfer function to be stable.
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 50
Feedback Controller
(42)
The roots of the denominator of the closed loop transfer function are {−3.93, 0.46 ± j7.96} .
The closed loop system has complex poles with positive real parts. Hence the open loop
stable system becomes unstable due to inappropriate choice of controller parameters. On
the other hand, if we choose a PI controller with gain 10 and reset 0.2, then the closed loop
transfer function will be
(43)
The roots of the denominator of the closed loop transfer function are {−2.63, −0.19 ± j4.36} .
The closed loop system has complex poles with negative real parts. Hence the system
remains stable.
Following observations are made from the above two exercises:
• Stability characteristics are determined by the poles of (1 + G c G p G f G m )
• An otherwise stable system can be destabilized by wrong choice of controller parameters
The Routh Hurwitz criterion does not calculate the actual values of closed loop poles, rather
it calculates whether any of the poles is on the left hand side of the imaginary axis of complex
plane. In other words, the Routh Hurwitz criterion find whether a closed loop transfer
function is stable. In this process, it also finds a limiting condition for controller parameters
which would ensure stability for closed loop system. The denominator of a closed loop
process transfer function is termed as characteristic equation. The Routh Hurwitz criterion
works on the characteristic equation as follows:
Step 1. Expand the characteristic equation into a polynomial
1 + G c G p G f G m= a 0s n + a1s n −1 + a 2s n − 2 + ... + a n −1s + a n (44)
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 51
Feedback Controller
The first column of the Routh array is {0.2, 0.6, −4.667,50} . Hence the system has unstable
pole(s). The sign is changed twice in the column (positive to negative and then negative to
positive). Hence the closed loop system has two unstable poles {0.46 ± j7.96} .
The Routh Array can also be used for finding the stability limits of the controller parameters.
(45)
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 52
Feedback Controller
Only the third element of the first column is capable of being negative due to inappropriate
choice of controller parameters. The stability of the closed loop system will be ensured if
K
3 × (10 + K c ) > c
τ1
The root locus analysis is a useful tool for determining the stability characteristics through
graphical means. Root loci are plotted in complex plane as the gain varies from zero to
infinity.
The Fig. 5.10 demonstrates the root locus diagram of the said example. The loci of the
complex pair of poles cross the imaginary axis at K c = 15 .
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 53
Feedback Controller
The following table shows the values of roots of the characteristic equation as Kc changes
from 0 to ∞.
Frequency response analysis is a useful tool for designing feedback controllers. It enables
the designer to study the stability characteristics of a closed loop system using Bode or
Nyquist plots and also to select the appropriate design values of controller parameters.
It is also worth to study the Bode diagram of PID controller in the following figure.
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 54
Feedback Controller
And
It is left to the reader to exercise how to arrive at above eqs. Note that the PID controllers
ideally have three asymptotes.
• At ω → 0, AR → ∞ with a slope ( −1)
• At ω → 0, AR → ∞ with a slope ( +1)
• At =
ω 1/ τ1τD , AR / K=
c 1 with a zero slope.
Let us have a control loop with two components, viz ., one PI controller and one first order
plus dead time model.
5e − s 1
G= =
;G 3 1 +
10s + 1
p c
20s
The open loop transfer function is
1 5e −2s 1 1 −2s
G OL =G pGc = 3 1 + = 15 1 + e
20s 10s + 1 20s 10s + 1
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 55
Feedback Controller
The above open loop transfer function is a combination of four individual processes in series,
viz., pure gain, pure dead time, first order system and PI controller. Two time constants are
observed in the series that would yield the location of corner
frequencies viz ., ωco1
= 1/ 20; ωco2
= 1/10 or ωco1= 0.05; ωco2= 0.1
The above figure shows the Bode plot of the open loop process indicated in this example. AR
and phase shift for all individual transfer functions as well as the overall transfer function
have been indicated along with the location of corner frequencies.
Consider a simple first order plus dead time process to be controlled by a proportional
controller:
Fig. 5.13: Example of first order system for studying Bode stability criterion
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 56
Feedback Controller
The Bode plot of the above open loop transfer function is given by the following figure.
We are interested to know the frequency where the phase shift is −180 . Numerically it can
be solved by the equation
= −π tan −1 ( −4ωco ) + ( −0.3ωco ) (49)
The frequency is ωco =5.39 and the value of AR / 2K c at ωco is observed to be 0.0463 which
can also be found numerically by,
(50)
The above exercise indicates that in order to obtain AR = 1 at this frequency ( ωco ) one
needs to set the value of K c as,
(51)
At this juncture, one needs to perform a thought experiment in order to understand the Bode
stability criterion. Let us set the value of controller gain, K c = 10.8 and let us “open up” the
feedback loop as indicated in the figure before. Suppose, we vary the setpoint as a sinusoidal
function y= sp ( t ) A sin ( ωco t ) . As the loop is open, the error will be equal to the
=
setpoint {ε A sin ( ωco t )} and thereby yield an output,
y ( t ) = ( AR × A ) sin {ωco t + ( −π )} = −A sin ( ωco t ) (52)
Now, suppose two events occur simultaneously
• Thesetpoint perturbation is stopped
• The feedback loop is reconnected
Then, the error signal will remain as ε ( t ) = 0 − {−A sin ( ωco t )} = A sin ( ωco t ) . In other words,
the response of the system will continue to oscillate with constant amplitude even when the
setpoint signal is withdrawn.
Alternatively, if we choose the value of controller gain less than 10.8 , (say K c = 8 ) then
AR = 2 × Κ c × 0.0463 =
2 × 8 × 0.0463 =0.74 (53)
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 57
Feedback Controller
If we repeat the above thought experiment, the output signal will take the form
y ( t ) = ( AR × A ) sin {ωco t + ( −π )} = −0.74A sin ( ωco ) t (54)
Upon closing the loop and withdrawing the setpoint perturbation, the new value for the error
=
for the next cycle will be ε 0.74A sin ( ωco t ) that will eventually yield an output response
of y ( t )
= ( 0.74 ) A sin ( ωco t ) and so on. It is evident that the amplitude of the error signal
2
The Bode stability criterion states that the maximum value of the controller gain that can be
chosen for stable closed loop response is called the ultimate gain Ku. In other words, the
value of controller gain must always be less than K u in order to ensure stability. The gain
margin (GM) is a design parameter such that
K (56)
Kc = u
GM
Gain margin should always be chosen as greater than one (GM>1) to ensure stability.
Gain margin acts as a safety factor for model uncertainty. Since process parameters such as
gain, time constant and dead time can never be estimated exactly, a safety factor of
magnitude more than one is necessary for stable operation. For relatively well modeled
processes, a low safety factor will be acceptable whereas poorly modeled processes need
higher safety factors. For an example, let us choose GM=2 for the process we have discussed
above (eq.48), the design value of the controller gain = =
is K c 10.8 / 2 5.4 ; suppose there
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 58
Feedback Controller
exists a modeling error of 50% in estimating the dead time of the process and the true value
of the dead time is 0.45 instead of 0.3, then the revised value of crossover frequency is
= −π tan −1 ( −4ωco ) + ( −0.45ωco ) (57)
or, ωco =3.64 , and the corresponding K u = 7.3 which is still higher than the designed value
of K c = 5.4 . The system is still stable despite the error by 50% we made in estimation of dead
time of the process.
Phase margin is another safety factor which is used for controller design. Here we are
interested to compute a frequency ωpm that satisfies the following expression,
−π
= + φ tan −1 ( −4ωpm ) + ( −0.3ωpm ) (58)
φ is called phase margin (PM) and it is the extra phase lag needed to destabilize a system.
For an example, let us choose PM = 30 ωpm . ωpm can be calculated from the following
expression
−150
=
tan −1 ( −4ωpm ) + ( −0.3ωpm ) (59)
or, ωpm =
3.72 . The gain is designed from the expression
(60)
or, K c = 7.44 Suppose there exists a modeling error of 50% in estimating the dead time of the
process and the true value of the dead time is 0.45 instead of 0.3,then the phase lag encountered
by the process would be
tan −1 ( −4ωpm ) + ( 0.45ωpm ) =−182 (61)
which is 2 more than the safety limit for stability. Hence, the phase margin of 30 is not sufficient
for handling 50% error in dead time estimation. It is left to the reader to verify that a phase margin
of 45 will suffice for handling 50% error in dead time estimation.
Unlike process reaction curve method which uses open loop response data, Ziegler Nichols tuning
technique uses closed loop response data. The following settings are given by this technique for
feedback controllers:
Table IV.6: Controller settings through Ziegler-Nichols tuning technique
Controller
P
PI
PID
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 59
Feedback Controller
The Bode stability criterion is valid for systems where amplitude ratio and phase shift decreases
monotonically with ω Nyquist stability criterion does not have any such limitation and is
applicable in more general sense. The criterion states that:
If open-loop Nyquist plot of a feedback system encircles the point (-1,0) as the frequency ω varies
from −∞ to +∞ the closed loop response is unstable.
Question: An ideal PID controller has the transfer function 1 + 1 / ( 0.5s ) + 0.2s . The
frequency at which the magnitude ratio of the controller is 1, is
(a) 0.5/0.2 (b) 0.2/0.5 (c) 0.2×0.5 (d) 1/√0.2 × 0.5
Question: Which of the systems having the following transfer fucntions is stable?
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 60
Feedback Controller
(a) 1 / ( s2 + 2) (b) 1 / ( s2 – 2s + 3)
(c) 1 / ( s2 + 2s + 2) (d) e−20s / ( s2 + 2s − 1)
Solution: (c)
Linked data type questions 1 & 2: The response of a first order system with a
proportional controller K c has an offset to a step input. The offset can be eliminated by
which one of the following actions?
(a) Decreasing K c
(b) Adding an integral mode
(c) Adding a derivative mode
(d) No change needs to be done because offset cannot be eliminated
Answer: (a)
Q.2 By this action, what happens to the overall order of the system?
(a) Increased
(b) Decreased
(c) Does not change
Answer: (a) Decreasing K c increases offset.
Linked data type questions 3 & 4: Oscillations of a system may be reduced by one of the
following controllers
(a) Proportional controller
(b) Proportional-integral controller
(c) Proportional-derivative controller
(d) No controller can reduce
Answer: (c)
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 61
Feedback Controller
Question: Bode diagrams are generated from output response of the system subjected to
which of the following input?
(a) Impulse (b) step (c) ramp (d) sinusoidal
Answer: (d)
Question: According to Bode stability criterion, a system is unstable if the open loop
frequency response exhibits amplitude ratio exceeding unity at frequency for which phase
lag is:
(a) 0° (b) 45° (c) 90° (d) 180°
Answer: (d) Bode stability criterion: A feedback control system is unstable if the
Amplitude Ratio of the corresponding open loop transfer function is larger than 1 at the
cross-over frequency (frequency at which phase lag is 180° )
Question:Design a PID controller for the following process using Ziegler-Nichols tuning
4e − s
=Gp = ; G =
0.25; Gv 2
( 5s + 1) m
Solution:
The overall transfer function is
4e − s 2e − s
=
G OL G p G m=
Gv × 0.25=
×2
( 5s + 1) ( 5s + 1)
If the crossover frequency is denoted by ωco , then we have
−π tan −1 ( −5ωco ) + ( −ωco )
=
or
ωco =
1.68841
Hence the amplitude ratio can be given as
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 62
Feedback Controller
For a unit step input in its setpoint, find the offset, overshoot, decay ratio, period of
oscillation, natural period of oscillation, damping coeffcient.
Solution:The closed loop transfer function is
(a) The final steady state value of step input is 0.8889, hence the offset is (1 -0.8889) =
0.1111
(b) The natural period of oscillation is, √0.2222 = 0.4714
0.3333
(c) The damping coefficient is = 0.3536
2×0.4714
(d) The overshoot is
𝜋𝜋×0.3536
− =0.305
𝑒𝑒 √1−0.35362
(e) The decay ratio is, 0.3052 = 0.093
(f) The period of oscillation is,
2 × 𝜋𝜋 × 0.4714
= 3.1664
√1 − 0.35362
Question: Consider the characteristic equation of a fourth-order system
s 4 + 8s3 + 18s 2 + 16s + 5 =0
Check the stability of the system using Routh - Hurwitz criteria
Solution:
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 63
Feedback Controller
Routh Table
The elements of the first column of Routh array are positive and hence the system the above
characteristics equation is stable.
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 64
6 FREQUENCY RESPONSE ANALYSIS
Let the sinusoidal input u(t) = A sin ωt perturb the system. Then the output of the process
will be
(2)
Computing the constants C1 , C2 and C3 and taking inverse Laplace Transform of the above
equation we obtain,
(3)
After sufficiently long time ( t → ∞ ) , the first term disappears as Hence, using the identity
eq. we obtain,
(4)
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 65
Frequency Response Analysis
Fig.6.1: Input and Output wave profile for a frequency response analysis
The modulus (or absolute value or magnitude) of W is a 2 + b 2 and the argument (or phase
angle) is tan−1 𝑏𝑏�𝑎𝑎. Let us put 𝑠𝑠 = 𝑗𝑗𝑗𝑗 in the transfer function of the first order process as
(5)
As 𝐺𝐺(𝑗𝑗𝑗𝑗) is now a complex number, the modulus and argument can be calculated as,
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 66
Frequency Response Analysis
(6)
The last two relationships indicate the amplitude ratio and phase lag for the ultimate
response of the first order process. Hence the observations in the last subsection can also be
stated in the light of the above results as follows:
• Sinusoidal output wave has the same frequency as that of input sinusoid
• Amplitude Ratio between the output wave and input wave is AR = G ( jω )
• Output wave lags behind the input wave with a phase difference of =
arg {G ( jω)}
6.3 Example of frequency response of a second order system
The process is
(7)
Then,
(8)
and,
(9)
Question: A first order system with a time constant of 1 min is subjected to frequency
response analysis. At an input frequency of 1 radian/min, the phase shift is?
(a) 45° (b) −90° (c) −180° (d) −45°
lag φ tan −1 (=
Answer: (d) Phase = −τω) tan −1=
( −1) 45° i.e, phase shift = − 45°
Question: Derive analytical expression for the amplitude ratio and phase shift of the
ultimate frequency response of sinusoidal input for two first order systems in series.
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 67
Frequency Response Analysis
Solution: Two first order processes in series will have an overall transfer function block:
Similarly,
Now,
=
Similarly, φ2 tan −1 ( −τ2 ω)
Hence the overall phase shift is
Thus, the amplitude ratio and phase lag are given as:
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 68
Frequency Response Analysis
This is the response of a process under unit step change in its input. Find the transfer
function of the process and its response time.
Solution:It is evident from the transient behavior of the process that this is a second order
plus time delay underdamped system. Hence it has a transfer function of the form,
Hence, ξ =0.15
The period of oscillation, i.e. the time elapsed between which the transient attains two
successive peaks, is 12.7, hence recall the formula
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 69
Frequency Response Analysis
The response time is defined as the time when the transient reaches within ±5% of its final
steady state value and stays within that band. The visual observation shows the response
time to be 45 sec.
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 70
7 INSTRUMENTATION
7.1 Introduction
The instrumentation part refers mainly to the hardware of a control systems that includes
various measuring instruments, transmitters, valves etc . In this module, some basic
knowledge on instrumentation symbols as well as process drawings will be discussed.
The P&ID refers to the detailed drawing of plant layout that includes pictorial representation
of entire piping and instrumentation blocks used in a plant. It has been standardized by
American National Standards Institute (ANSI) and Instrument Society of America (ISA).
Equivalent Indian Standard is also available, e.g. “IEC/PAS 62424 Ed. 1.0 en - Representation
of process control engineering requests in P&I diagrams and data exchange between P&ID
tools and PCE-CAE tools”. It describes and specifies how process control engineering
requests are represented in a P&I diagram. It also defines the exchange of process control
engineering request relevant data between a process control engineering tool and a P&I tool
by means of a data transfer language (called CAEX). These provisions apply to the
export/import applications of such tools. Following are the components of P&ID:
• Plant equipment and vessels showing location, capacity, pressure, liquid level operating
range, usage and so on
• All interconnection lines distinguishing between the types of interconnection, i.e. gas or
electrical and operating range of line
• All motors giving voltage and power and other relevant information
• Instrumentation showing location of instrument, its major function, process control loop
number, and range
• Control valves giving type of control, type of valve, type of valve action, fail save features,
and flow plus pressure information
• The ranges for all safety valves, pressure regulators, temperatures, and operating ranges
• All sensing devices, recorders, and transmitters with control loop numbers
P&ID can run into multiple sheets as it is not always possible to code them all in a single
sheet. Each drawing should have a part list and that should be properly numbered. It should
also have an area for revisions, notes and approval signature. More importantly they should
always be up-to-date.
The symbols of transmission lines have already been discussed in the lecture 3 of Module 4
(Introduction to Feedback Controller). Symbols of standardized instruments are given in the
following figure.
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 71
Instrumentation
Each instrument is identified with a letter of English alphabet. Table. 7.1 refers to the
identification letters
Table 7.1: Instruments Identification Letters
Fig 7.2 represents a few examples of the letters and numbering code that can be used to refer
the instruments. Consider the example TY178 which has two letters T & Y. Table 7.1 refers
to the first letter T as temperature and second letter Y as converter. Hence it is a temperature
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 72
Instrumentation
transducer which converts a 4-20 mA current signal into a 3-15 psi pressure signal. The
number 178 refers to the location of the transducer such as “zone 1, equipment number. 7,
transducer number. 8”. From Fig 7.1, we can further infer that transducer is a discrete
instrument located in the field itself.
Symbols can be broadly classified into four types depending on their functions:
• Actuator
• Primary elements
• Regulators and safety valves
• Math function
Following are the examples of those categories
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 73
Instrumentation
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 74
Instrumentation
Process Flow Diagram, on the other hand, shows the relationships between the major
components in the system. PFD also tabulates process design values for the components in
different operating modes, typical minimum, normal and maximum. Following figure shows
a typical PFD.
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 75
Instrumentation
7.4 Actuators
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 76
Instrumentation
Pneumatic valves are air-operated valves that are used for regulating fluid flow through a
pipeline (Fig 7.9). Varying air pressure (3 to 15 psig) is used as an actuating signal for the
pneumatic valves. Usually control signal is generated and transmitted in electrical form (4 to
20 mA) from the controller. Hence, a current to pressure converter (I/P) converts an analog
signal (4-20 mA) to a proportional linear pneumatic output (3-15 psig). Its purpose is to
translate the analog output from a controller into a precise, repeatable pressure value to
control pneumatic actuators/operators, pneumatic valves, dampers, vanes, etc. .
The bottom portion of the pneumatic valve has an orifice that separates the upstream and
downstream flows. A tapered plug, capable of blocking the orifice in varying proportion, is
attached to a stem that is connected to a diaphragm in the top portion of the valve. A spring
restricts the movement of the stem. When air pressure above the diaphragm forces the stem
to move downwards, the plug starts reducing the aperture of the orifice and eventually
blocks the orifice at high air pressure. As a result the flow of fluid through the orifice is
gradually decreased from “FULL FLOW” to “NO FLOW” condition. As the air pressure at the
top of the diaphragm is released, the plug moves back to its original position resulting in full
flow of fluid. This is called “FAIL OPEN” valve (Fig. 7.9 (a)) because when the control signal
fails to provide enough air pressure, the valve remains in fully open condition. Similarly the
“FAIL CLOSED” valve is shown in Fig 7.9 (b) where the shape of the plug is opposite.
The dynamics of a typical pneumatic valve is usually second order in nature. The position of
stem (or plug) determines the size of the aperture of the orifice that consequently
determines the fluid flow rate. The position of the stem is determined by balancing all the
forces acting on it. These forces are:
• Force exerted by the compressed air at the top of the diaphragm
• Force exerted by the spring
• Force exerted due to friction between stem and the valve packing
Hence,
(1)
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 77
Instrumentation
(2)
The transfer function between stem position and the actuating pressure is
(3)
which is indeed a second order dynamics. However, the movement of stem is always aided
by proper lubrication which eventually reduces its coefficient of friction. Moreover, mass of
the stem and its attachments are very small compared to the Hooke's constant
( M << K and C << K ) . Hence the dynamics of pneumatic valve can be neglected and
movement of stem x can be regarded as instantaneous in response to change of air pressure.
(4)
Where, F =flow rate of liquid, ∆P = pressure drop across the upstream and downstream side
of the valve, ρ = density of liquid, Cv =flow coefficient that depends on valve size, f(x) = flow
characteristic curve of the valve which is a function of the stem position x.
The shape of the flow characteristic curve depends on the shape of plug of the valve. There
are five major types of flow characteristics used for designing the valve:
• Linear, f ( x ) = x
• Square root, f ( x ) = x
• Equal percentage, f ( x ) = α x −1
{α − ( α − 1) x}
−1
• Hyperbolic, f ( x ) =
• Quick opening, f ( x ) = x1/ α
Where, α =design constant. The types of plug are given in Figs. (7.10).
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 78
Instrumentation
Hydraulic actuators employ hydraulic pressure to move a target device. These are used
where high speed and large forces are required. Pressure applied to a confined
incompressible fluid at any point is transmitted throughout the fluid in all directions and acts
upon every part of the confining vessel at its interior surfaces. Fig. 7.12 represents the
schematic of a hydraulic actuator.
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 79
Instrumentation
Electric actuators can be of the following forms: Electric current/voltage, Solenoid, Stepping
Motor and DC/AC Motor.Current or voltage can easily be regulated to adjust the flow of
energy into the process, e.g. heater, fan speed regulator etc. Energy supplied by the heater
element is W = I2Rt , where I =current, R =resistance, t =time of heating. The current/voltage
can be regulated usingpotentiometer (or rheostat), amplifier or a switch.
A rheostat is a device that has variable resistance to current flow (Fig. 7.13). The current
=
flowing through the circuit is I V / ( R1 + R 2 ) , where R1 can be varied by suitably sliding
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 80
Instrumentation
the pointer. The power transmitted to the heater would be P = I 2 R 2 . As the pointer slides
towards a lesser value of R1 , heater receives more power for heating.
7.5 Sensors
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 81
Instrumentation
Most common devices for relatively low temperature measurement are thermocouples and
resistance thermometers. Radiation pyrometers are used for high temperature
measurements. Other temperature measurement devices are filled system thermometers,
bimetal thermometers, oscillating quartz crystals, etc.
7.5.1.1 Thermocouple
Thermocouples are the most popular temperature sensors. They are inexpensive,
interchangeable, have standard connectors and can measure a wide range of temperatures.
Their main limitation is accuracy as the system errors of less than 1°C can be difficult to
achieve. Following figure represents internal construction of thermocouple and its circuitry.
(b) Circuit
(a) Internal Construction
Fig. 7.16: Thermocouple
A thermocouple is constructed of two dissimilar metal wires joined at one end. It works on
the principle of “Seabeck effect” whereby electromagnetic force is generated when two
dissimilar metals are joined at two different temperature ends. When one end of each wire
is connected to a measuring instrument, the thermocouple becomes a sensitive and highly
accurate measuring device. Heating the measuring junction of the thermocouple produces a
voltage which is greater than the voltage across the reference junction. The difference
between the two voltages is proportional to the difference in temperature and can be
measured on the voltmeter (in mV). Thermocouples may be constructed of several different
combinations of materials. The most important factor to be considered when selecting a pair
of materials is the "thermoelectric difference" between the two materials. A significant
difference between the two materials will result in better thermocouple performance.
The RTD incorporates pure metals or certain alloys that increase in resistance as
temperature increases and, conversely, decrease in resistance as temperature decreases.
RTDs act somewhat like an electrical transducer, converting changes in temperature to
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 82
Instrumentation
voltage signals by the measurement of resistance. The metals that are best suited for use as
RTD sensors are pure, of uniform quality, stable within a given range of temperature, and
able to give reproducible resistance-temperature readings. RTD elements are normally
constructed of platinum, copper, or nickel. These metals are best suited for RTD applications
because of their linear resistance-temperature characteristics.
(6)
The above equation represents the Resistance vs. Temperature relationship where R is the
resistance and temperature T, with reference Resistance R 0 and coefficient of resistance
as α The coefficient of resistance is the change in resistance per degree change in
temperature, usually expressed as a percentage per degree of temperature.
Most common devices for pressure measurement are differential pressure transducer.
The differential pressure cell (DPC) measures the difference between two or more pressures
introduced as inputs to its sensing unit. The sensing unit consists of a diaphragm and a
pressure cavity to create a variable capacitor which detects strain due to applied pressure.
The following figure presents a schematic of the sensing unit of the DPC.
The cell contains two compartments separated by a diaphragm. The differential pressure
acting on the diaphragm, due to pressures exerted on its two sides, results in displacement
of the diaphragm. This displacement of diaphragm is transduced to electrical signal which is
transmitted and recorded through appropriate instrumentation. The displacement of
diaphragm is usually linearly calibrated to the differential pressure acting on it. The ends of
the DPC are capillary tubes. Compartment is filled with some liquid ( e.g . silicone oil) of
mass m. The ends are connected to processing units of which the differential pressure needs
to be measured. Often the End 2 is left open to the atmosphere when the pressure at only
one end needs to be measured. The following force balance will be in place:
Force due to pressure p1 at End 1- Force due to pressure p2 at End 2 - Resistance exerted by
the diaphragm - Viscous frictional force exerted by the capillary liquid = Mass of capillary
liquid x Acceleration of capillary liquid towards the direction of motion
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 83
Instrumentation
where A= cross-sectional area of the capillary tube, K= Hooke's constant of diaphragm which
acts like a spring, x= displacement of diaphragm/capillary liquid, C= damping coefficient of
capillary liquid.
Rearranging the above, we obtain
(8)
In Laplace domain,
(9)
The value of x is linearly transduced to electrical signal that indicates the actual
measurement variable. Hence the output of DP cell maps a second order behavior with the
differential pressure.
Bellows Resistance Transducer is resistance type transducers which has moving contacts
(slide wire variable resistors). The pressure change causes a resistance change due to the
distortion of the wire. The value of the pressure can be found by measuring the change in
resistance of the wire grid. Eq.(10) shows the resistance relationship.
(10)
Where R is resistance of the wire grid (in ohms) K is resistivity constant for the particular
type of wire grid, L is the length of wire grid and A is the cross sectional area of wire grid.
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 84
Instrumentation
Fig. 7.18 shows the schematic of a Bellows Resistance Transducer. As pressure changes, the
bellows will either expand or contract. This expansion and contraction causes the attached
slider to move along the slide wire, increasing or decreasing the resistance, and thereby
indicating an increase or decrease in pressure.
Various flow measurement devices are already taught in the relevant course of Fluid
Mechanics such as Orificemeter, Venturimeter, pitot tube, rotameter etc. . A few other types
of flow meters are discussed here.
The hot-wire anemometer, consists of an electrically heated, fine platinum wire which is
immersed into the flow. As the fluid velocity increases, the rate of heat flow from the heated
wire to the flow stream increases. Thus, a cooling effect on the wire electrode occurs, causing
its electrical resistance to change.
In a displacement flowmeter, all of the fluid passes through the meter in almost completely
isolated quantities. The number of these quantities is counted and indicated in terms of
volume or weight units by a register. The most common type of displacement flowmeter is
the nutating disk, or wobble plate meter. A typical nutating disk is shown in Fig.7.20.
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 85
Instrumentation
The movable element is a circular disk which is attached to a central ball. A shaft is fastened
to the ball and held in an inclined position by a cam or roller. The disk is mounted in a
chamber which has spherical side walls and conical top and bottom surfaces. The fluid enters
an opening in the spherical wall on one side of the partition and leaves through the other
side. As the fluid flows through the chamber, the disk wobbles, or executes a nutating motion.
Since the volume of fluid required to make the disc complete one revolution is known, the
total flow through a nutating disc can be calculated by multiplying the number of disc
rotations by the known volume of fluid. To measure this flow, the motion of the shaft
generates a cone with the point, or apex, down. The top of the shaft operates a revolution
counter, through a crank and set of gears, which is calibrated to indicate total system flow.
The differential pressure (P) detector method is also used for liquid level measurement
which uses a P detector connected to the bottom of the tank being monitored. The higher
pressure, caused by the fluid in the tank, is compared to a lower reference pressure (usually
atmospheric). This comparison takes place in the P detector. In addition to this a couple of
other level detectors are discussed below.
The ball float method is a direct reading liquid level mechanism. The most practical design
for the float is a hollow metal ball or sphere. However, there are no restrictions to the size,
shape, or material used. The design consists of a ball float attached to a rod, which in turn is
connected to a rotating shaft which indicates level on a calibrated scale (Fig.7.21).
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 86
Instrumentation
The operation of the ball float is simple. The ball floats on top of the liquid in the tank. If the
liquid level changes, the float will follow and change the position of the pointer attached to
the rotating shaft.
The magnetic bond mechanism consists of a magnetic float which rises and falls with changes
in level. The float travels outside of an on-magnetic tube which houses an inner magnet
connected to a level indicator.
When the float rises and falls, the outer magnet will attract the inner magnet, causing the
inner magnet to follow the level within the vessel (Fig.7.22).
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 87
8 ADVANCE CONTROL
A feedback controller responds only after it detects a deviation in the value of the controlled
output from its desired set point. On the other hand, a feedforward controller detects the
disturbance directly and takes an appropriate control action in order to eliminate its effect
on the process output.
Consider the distillation column shown in Fig (8.1) The control objective is to keep the
distillate concentration at a desired set point despite any changes in the inlet feed stream.
Fig 8.2: The comparative schematic of feedback and feedforward control structure
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 88
Advance Control
Fig 8.2 shows the general form of a feedforward control system. It directly measures the
disturbance to the process and anticipates its effect on the process output. Eventually it
alters the manipulated input in such a way that the impact of the disturbance on the process
output gets eliminated. In other words, where the feedback control action starts after the
disturbance is “felt” through the changes in process output, the feedforward control action
starts immediately after the disturbance is “measured” directly. Hence, feedback controller
acts in a compensatory manner whereas the feedforward controller acts in
an anticipatory manner.
Let us consider the block diagram of a process shown in Fig 8.3. The Fig 8.3(a) presents the
open-loop diagram of the process. The process and disturbance transfer functions are
represented by G p and G d respectively. The controlled output, manipulated input and the
disturbance variable are indicated as y, u and d respectively.
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 89
Advance Control
(1)
The control objective is to maintain y at the desired setpoint ysp . Hence the eq (1) can be
rewritten as
(2)
The eq. (2) can be rearranged in the following manner:
or
.....................................................(3)
For the sake of simplicity, measuring element and final control element were not considered
as parts of the feedforward control configuration as shown in Fig 8.3(b). In a more
generalized case, when such elements are added in the controller configuration, the resulting
control structure takes the form of Fig 8.3(c). A generalized form of controller equation can
be written as
..............................................................................................(4)
And
In case of regulatory problem (disturbance rejection) i.e. when ysp = 0 , the controller should
be able to reject the effect of disturbance and ensure no deviation in the output, i.e. y = 0 . In
other words,
G d − G pG cG f G m =0 (5)
Or
Gd
Gc = (6)
G mG f G p
In case of servo problem (setpoint tracking), i.e. when d = 0 , the controller should be able to
ensure that output tracks the setpoint, i.e. y = ysp . In other words,
G p G c G f G sp = 1 (7)
Or
1 1Gm
=G sp = = (8)
G pG cG f Gd Gd
Gp Gf
G G G
m f p
Consider an overflow type continuous stirred tank heater shown in Fig 8.4. The fluid inside
the tank is heated with steam whose flow rate is Fst and supplying heat at a rate of Q to the
fluid. Temperatures of the inlet and outlet streams are Ti and T respectively. V is the volume
of liquid which is practically constant in an overflow type reactor. A control valve in the
steam line indicates that the steam flow rate can be manipulated in order to keep the liquid
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 90
Advance Control
temperature at a desired setpoint. Temperature of the inlet stream flow is the source of
disturbance (change in Ti ) to the process.
Fig 8.4: Feedforward control configuration of an overflow type continuous stirred tank heater
A simple energy balance exercise will yield the model equation of the above process as:
(9)
All the variables are assumed to be in the deviation form. Hence, taking Laplace transform
on both sides we obtain:
(10)
or, (11)
or, (12)
or,
(13)
Or
(15)
Hence, one needs to set Fst in such a way that Q amount of heat as given in eq. (15) is
transferred to the process. Fig 8.4(b) represents the feedforward structure of the controller.
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 91
Advance Control
Remarks:
• The feedforward controller ideally does not get any feedback from the process output.
Hence, it solely works on the merit of the model(s). The better a model represents the
behavior of a process,the better would be the performance of a feedforward controller
designed on the basis of that model. Perfect control necessitates perfect knowledge of
process and disturbance models and this is practically impossible. This inturn is the main
drawback of a feedforward controller.
• The feedforward control configuration can be developed for more than one disturbance in
multi-controller configuration. Any controller in that configuration would act according to
the disturbance for which it is designed.
• External characteristics of a feedforward loop are same as that of a feedback loop. The
primary measurement (disturbance in case of feedforward control and process output in
case of feedback control) is compared to a setpoint and the result of the comparison is used
as the actuating signal for the controller. Except the controller, all other hardware elements
of the feedforward control configuration such as sensor, transducer, transmitter, valves are
same as that of an equivalent feedback control configuration.
• Feedforward controller cannot be expressed in the feedback form such as P, PI and PID
controllers. It is regarded as a special purpose computing machine.
• Let us consider a system where process delay is higher than disturbance delay, eg.
5e −0.1s 5e −0.02s Gd
Gp = and G d = =
; in such case, G = e0.08s . That means one needs to know
s +1 s +1
c
Gp
the future values of disturbance in order to decide present control action. This is physically
unrealizable controller.
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 92
Advance Control
Insensitive to mild errors in modeling Bad for sluggish systems and/or system
with large deadtime
Insensitive to mild changes in process May affect the stability of the process
parameters
Let us now explore how a combination of feedforward and feedback controller would
perform when they are designed to act simultaneously. The schematic of a feedforward-
feedback controller is shown in Fig 8.5.
Without losing the generality we shall ignore the transfer functions of the measuring element
and the final control element.
Now the closed loop transfer function of feedforward-feedback controller can be derived in
the following manner:
(16)
It is observed that the stability of the closed loop response is determined by the roots of the
characteristic equation: 1 + G p G cFB =0 . Hence, the stability characteristics of a process does
not change with the addition of a feedforward loop.
The following numerical example demonstrates the efficacy of a feedforward-feedback
controller. Consider a process having process and disturbance transfer functions as
(18)
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 93
Advance Control
A feedback PID controller,with K c= 10; τ1= 1 , is used to control the process for disturbance
rejection purpose. A feedforward controller has also been designed for the process however
it has been assumed that the time constant of the process has been measured erroneously as
2.1 instead of 2. A first order filter with time constant 0.1 has been augmented to the transfer
function of G sp in order to make it causal.
A Simulink code (Fig 8.6) has been generated to simulate the process under the three types
of controllers as said above.
Fig 8.6: Block diagram of Simulink code for closed-loop process under three controllers
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 94
Advance Control
The primary disadvantage of conventional feedback control is that the corrective action for
disturbances does not begin until after the controlled variable deviates from the setpoint. In
other words, the disturbance must be “felt” by the process before the control system
responds. Feedforward control offers large improvements over feedback control for
processes that have large time constant and/or delay. However, feedforward control
requires that the disturbances be measured explicitly and that a model be available to
calculate the controller output. Cascade control is an alternative approach that can
significantly improve the dynamic response to disturbances by employing a secondary
measurement and a secondary feedback controller. The secondary measurement point is
located so that it recognizes the upset condition sooner than the controlled variable, but the
disturbance is not necessarily measured.
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 95
Advance Control
= + d 2 G p2 G c2 ( Ysp2 − Y2 ) + d 2
Y2 G p2 G c2 e 2= 19
Simplifying
G p2 G c2 1 20
=Y2 Ysp2 + d
(1 + G p2 G c2 ) (1 + G p2G c2 ) 2
Similarly
=
Y1 G p1Y2 + d1
G G 1
=
Y1 G p1 p2 c2 Ysp2 + d 2 + d1
1 + G p2 G c2 1 + G p2 G c2
G G 1
=
Y1 G p1 p2 c2 G c1e1 + d 2 + d1
1 + G p2 G c2 1 + G p2 G c2
G G
Y1 G p1 p2 c2 G c1 ( Ysp1 − Y1 ) +
1
= d 2 + d1
1 + G p2 G c2 1 + G p2 G c2
Again simplifying the above eqn:
G p1G p2 G c2 G c1 G p1 1 + G p2 G c2 (21)
Y1 = Ysp1 + d2 + d1
1 + G p2 G c2 + G p1G p2 G c2 G c1 1 + G p2 G c2 + G p1G p2 G c2 G c1 1 + G p2 G c2 + G p1G p2 G c2 G c1
Now see what happens if the secondary loop is absent. In that case:
Y2 = G p2 Ysp2 (22)
And
=
Y1 G p1Y2 + d1
(23)
=Y1 G p1 G p2 Ysp2 + d 2 + d1
=Y1 G p1 G p2 G c1e1 + d 2 + d1
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 96
Advance Control
K p2 K p2 K c2
K c2 (25)
Y2 G p2 G c2 τp2s + 1 K p2 K c2 1 + K p2 K c2
= = =
Ysp2 1 + G p2 G c2 K p2 τp2s + 1 + K p2 K c2 τp2
1+ K c2 s +1
τp2s + 1 1 + K K
p2 c2
τp2
Time constant =
1 + K p2 K c2
Clearly there is a decrease in time constant in presence of slave loop. This will guarantee that
the time constant of the overall process decreases when a secondary loop is introduced to
the system. Eventually that decreases the overall response time.
The stability part can also be explained in general derivation, but it will become too
cumbersome. So let us take a simple example.
1 1
=
Let us take, G p1 = ;G p2 . Thus overall system is third order.
( s + 1) s +1
2
In absence of the secondary loop the characteristic equation of the overall process is:
1 1 ( s + 1) + K
3 (26)
1 + G p1G p2 G c1 =
1+ K c1 =3 c1
( s + 1) s + 1 ( s + 1)
2
In presence of the secondary loop the characteristic equation of the overall process is:
( s + 1) + ( s + 1) K c2 + K c1K c2 (27)
3 2
1 1 1
1 + G p2 G c2 + G p1G p2 G c1G c2 =
1+ K c2 + K c1K c2 =
s +1 ( s + 1) s + 1 ( s + 1)
2 3
Now apply Routh-Hurwitz criteria to find the limits of both these cases:
Routh table:
1 3
3 1 + K c1 .
9 − (1 + K c1 )
3
For a stable closed loop process, 9 > 1 + K c1 . Thus the upper limit of K c1 is 8.
= s3 + ( 3 + K c2 ) s 2 + ( 3 + 2K c2 ) s + (1 + K c2 + K c1K c2 )
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 97
Advance Control
Routh table:
1 3 + 2K c2
3 + K c2 1 + K c2 + K c1K c2 .
( 3 + K c2 )( 3 + 2K c2 ) − (1 + K c2 + K c1K c2 )
3 + K c2
For a stable closed loop process, ( 3 + K c2 )( 3 + 2K c2 ) > (1 + K c2 + K c1K c2 ) . With further
8
simplification, we get: 8 + K c2 + > K c1
K c2
Thus the upper limit of
8
K c1 increases by a positive quantity of K c2 + . In other words, the stability limits of the
K c2
process increases by introducing a secondary loop.
One may put some reasonable values to the controller gain and simulate the system.
However it will generate some offset for obvious reasons. To get an offset free control action,
one needs to assign the master controller with a Gain-Reset setting.
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 98
Advance Control
The chemical process industries have various applications for ratio controllers. Following
are a few such examples:
• Reflux ratio and reboiler feed ratio in a distillation column
• Maintaining the stoichiometric ratio of reactants in a reactor
• Keeping air/fuel ratio in a combustion process
During the operation of a process plant it is possible that a dangerous situation may arise
due to unacceptable process conditions which may destruct the process or its personnel. In
such case the normal operation should temporarily be stopped and preventive measures
should be initiated to avert the unacceptable situation. In order to facilitate such measures,
a single-purpose “switch” can be used that can take preferential instruction from one
controller over the others to manipulate the final control element in such a way that the
dangerous situation can be averted. This is called override control. The technique can be
illustrated with the following example.
Consider a boiler shown in the Fig 8.10. It has one water inlet and one steam outlet. The
steam outlet is regulated by the valve in the discharge line that takes the control signal from
the control mechanism in Loop1 (pressure transducer and pressure controller). In other
words, the discharge of steam is regulated on the basis of its pressure desired in the supply
line elsewhere. However, the water is boiled using a heating coil that needs to be always
submerged below the water level so that the heating coil does not burn out. Hence, in order
to ensure a certain minimum level of water inside the boiler, the control Loop 2 is set in place
that contains a level transducer and a level controller. Both level controller and pressure
controller give the control signal to the valve through an intermediate switch LSS (Low
Selection Switch) that takes the preferential signal from the level controller. In other words,
Loop 2 remains inactive during the normal operation and the Loop 1 regulates the process.
Nevertheless, at critical situation when the water level drops below the minimum allowable
limit, the Loop 2 takes over and takes corrective measures.
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 99
Advance Control
There are conditions in process plant where multiple process measurements are available
for a particular variable that needs to be regulated through a single control action. Thus it is
evident that the said control action should be given based on the most critical measurement
condition for the process variable. This is termed as Auctioneering Control. The technique
can be illustrated with the following example.
When there are more than one manipulated inputs available for one single controlled output
then Split Range Control scheme can be implemented. The technique can be illustrated with
the following example.
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 100
Advance Control
Consider a reactor shown in the Fig 8.12. The pressure inside the reactor needs to be
controlled. The control valve is available in both inlet and outlet lines. Hence the control
action can be coordination between two valves. When the pressure is low, inlet valve should
be FULLY OPEN and when the pressure is high, outlet vale should be FULLY OPEN. In either
care the other valve needs to be PARTIALLY OPEN/CLOSED depending upon the
requirement of control action.
It is understood in the previous chapters that task of controller design for a process is very
much domain specific. First the model of the process is linearized around a certain nominal
point and the controller is designed on the basis of that linearized model and finally
implemented on the process. Hence, the controller is applicable for certain domain around
the nominal operating point around which the model has been linearized. However, if the
process deviates from the nominal point of operation, controller will lose its efficiency. In
such cases, the parameters of the controller need to be re-tuned in order to retain the
efficiency of the controller. When such retuning of controller is done through some
“automatic updating scheme”, the controller is termed as adaptive controller. The technique
can be illustrated with the following figure.
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 101
Advance Control
One of the most popular adaptive control techniques is gain scheduling technique. The
overall gain of an open loop process is usually given as
K overall = K c K p K m K f (29)
It is customary to keep the overall gain constant. In case of changes in the process (or valve
characteristics or measuring element), K c should be tuned in such a manner that overall
gain remains constant.
K overall (30)
Kc =
KpKmKf
The above is called the gain scheduling control law.
When the process is poorly known, one cannot rely much on the correctness of the value
of K p In such cases, the self-adaptive control may be helpful. A self-adaptive controller
optimizes the value of certain objective function (criterion) in order to obtain updated
controller parameters. Two examples of self-adaptive controllers are Model Reference
Adaptive Control (MRAC) and Self-Tuning Regulator (STR)
The following figure shows the schematic of a Model Reference Adaptive Controller.
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 102
Advance Control
Self-Tuning Regulator on the other hand estimates the model parameters by measuring
process inputs and outputs. The re-tuned model eventually guides the controller parameter
adjustment mechanism. Figure 8.15 shows the schematic of Self Tuning Regulator.
Often the process plant has certain variables that cannot be measured on-line, however,
needs to be controlled on-line. In such cases, the unmeasured variables to be controlled can
be estimated by using other measurements available from the process. Consider the
following example:
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 103
Advance Control
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 104
GATE QUESTIONS
Q.1 The transfer function of a PD Q.6 The response of two tanks of same
controller is size and resistance in series is
(a) 𝐾𝐾𝑐𝑐 (1 + 𝜏𝜏𝐷𝐷 𝑆𝑆) (b) 𝐾𝐾𝑐𝑐 (1 + 1�𝜏𝜏 𝑆𝑆) (a) Under damped
𝐷𝐷
(b) Critically damped
(c) 𝐾𝐾𝑐𝑐 (𝜏𝜏𝐷𝐷 𝑆𝑆) (d) 𝐾𝐾𝑐𝑐 /𝜏𝜏𝐷𝐷 𝑆𝑆
(c) Over damped
[GATE-1995]
(d) None of the above
[GATE-1995]
Q.2 The offset introduced by proportional
controller with gain K c in response of first
Q.7 The transfer function of a pure
order system can be reduced by dead time system 𝜏𝜏𝐷𝐷 is
(a) Reducing value of K c 1
(a) (b) 1 + 𝜏𝜏𝐷𝐷 𝑆𝑆
(b) Introducing integral control 1+𝜏𝜏𝐷𝐷 𝑆𝑆
(c) Introducing derivative control (c) 𝑒𝑒 −𝜏𝜏𝐷𝐷 𝑆𝑆 (d) 𝑒𝑒 𝜏𝜏𝐷𝐷 𝑆𝑆
(d) None of the above [GATE-1995]
[GATE-1995]
Q.8 The time constant of a first order
Q.3 Bode diagrams are generated from process with resistance r and capacitance
output response of the system subjected to C is
which of the following input (a) 𝑅𝑅 + 𝐶𝐶 (b) 𝑅𝑅 − 𝐶𝐶
(a) Impulse (b) step 1
(c) ramp (d) sinusoidal (c) 𝑅𝑅𝑅𝑅 (d)
𝑅𝑅𝑅𝑅
[GATE-1995] [GATE-1995]
Q.4 The root locus method, a pole of a Q.9 When a thermometer at 30°C is
transfer function G(s) is the value of s for placed in a water bath at 90°C , the initial
which G(s) approaches rate of rise in thermometer temperature is
(a) -1 (b) 0 found to be 2°C / sec . What is the time
(c) 1 (d) ∞ constant of the temperature, assuming it is
[GATE-1995] a first order device with unity steady state
gain? What will thermometer read after
Q.5 According to Bode stability one minute?
criterion, a system is unstable if the open [GATE-1995]
loop frequency response exhibits an
amplitude ratio exceeding unity at Q.10 A first order system with transfer
frequency for which phase lag is function 𝐺𝐺𝑃𝑃 = 𝑌𝑌(𝑆𝑆)⁄𝑋𝑋(𝑆𝑆) = 1(𝑆𝑆 + 1) is
(a) 0° (b) 45° subjected to input X(t) = t. Derive the
(c) 90° (d) 180° expression for change in output Y(t) as a
[GATE-1995] function of time. What is the maximum and
minimum difference between input and
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 105
Gate Questions
output? At what time does these Q.15 The transfer function for an ideal
differences occur? proportional plus reset controller (Reset
[GATE-1995] time T) is
1
(a) 𝐾𝐾𝑐𝑐 �1 + � (b) 𝐾𝐾𝑐𝑐 [1 + 𝑇𝑇𝑇𝑇]
Q.11 Amplitude ratio corresponding to 𝑇𝑇𝑇𝑇
crossover frequency of 10 rad/sec for a 𝐾𝐾𝑐𝑐 𝐾𝐾𝑐𝑐
(c) (d) 𝑠𝑠
process is found to be 0.01. (a) What is the [1+𝑇𝑇𝑇𝑇] �1+𝑇𝑇�
ultimate gain? (b) What is the period of [GATE-1996]
sustained oscillations corresponding to
ultimate gain? (c) If a proportional Q.16 The transfer function of a PID
controller is to be used, what is the value controller is
of proportional gain you would (a) 𝐾𝐾𝑐𝑐 (1 + 𝜏𝜏1 𝑠𝑠 + 𝜏𝜏𝐷𝐷 𝑠𝑠)
recommend? 1
[GATE-1995] (b) 𝐾𝐾𝑐𝑐 �1 + + 𝜏𝜏𝐷𝐷 𝑠𝑠�
𝜏𝜏1 𝑠𝑠
1
Q.12 In a feed-back control system G (c) 𝐾𝐾𝑐𝑐 �1 + 𝜏𝜏1 𝑠𝑠 + �
𝜏𝜏𝐷𝐷 𝑠𝑠
and H denote open loop and closed loop 1 1
transfer functions respectively. The (d) 𝐾𝐾𝑐𝑐 �1 + + �
𝜏𝜏1 𝑠𝑠 𝜏𝜏𝐷𝐷 𝑠𝑠
output-input relationship is [GATE-1997]
𝐺𝐺 𝐻𝐻
(a) (b)
1+𝐻𝐻 1+𝐺𝐺 Q.17 The transfer function for a first-
𝐺𝐺 𝐻𝐻 order process with time delay is
(c) (d)
𝐻𝐻 𝐺𝐺 𝑒𝑒 𝜏𝜏 𝑑𝑑 𝑠𝑠 𝑒𝑒 −𝜏𝜏 𝑑𝑑𝑠𝑠
[GATE-1996] (a) (b)
𝜏𝜏𝜏𝜏+1 𝜏𝜏𝜏𝜏+1
1 𝜏𝜏𝑑𝑑 𝑠𝑠
Q.13 The open loop transfer function of (c) (d)
(𝜏𝜏𝜏𝜏+1)(𝜏𝜏𝑑𝑑 𝑠𝑠+1) 𝜏𝜏𝜏𝜏+1
KR [GATE-1997]
a control system is this represents
1 + τs '
(a) A first order system Q.18 The open-loop transfer function of
(b) Dead time system (𝑠𝑠+1)(𝑠𝑠+4)
(c) A first order time lag a process is 𝐾𝐾 . In the root
(𝑠𝑠+2)(𝑠𝑠+3)
(d) A second order system locus diagram, the poles will be at
[GATE-1996] (a) -1, -4 (b) 1, 4
(c) -2,-3 (d) 2, 3
Q.14 The Laplace transform of a [GATE-1997]
staircase function s (b, t) formed by
successive addition of unit step functions Q.19 For an input forcing function,
at 0, b, 2b, 3b, etc.is
1 1 X ( t ) = 2t 2 , the Laplace transform of this
(a) (b) function is
𝑠𝑠 𝑏𝑏𝑏𝑏
1 2 4
(c) 𝑠𝑠(1 − 𝑒𝑒 −𝑏𝑏𝑏𝑏 ) (d) (a) 2 (b)
𝑠𝑠(1−𝑒𝑒 −𝑏𝑏𝑏𝑏) 𝑠𝑠 𝑠𝑠 2
2 4
[GATE-1996] (c) 3 (d)
𝑠𝑠 𝑠𝑠 3
[GATE-1997]
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 106
Gate Questions
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 107
Gate Questions
Q.28 The second order system with the (a) becomes more oscillatory (frequency
4 increases)
transfer function 2 has a damping
𝑆𝑆 +2𝑆𝑆+4 (b) becomes less oscillatory (frequency
ratio of decreases)
(a) 2.0 (b) 0.5 (c) gets amplified (magnitude increases)
(c) 1.0 (d) 4.0 (d) gets attenuated (magnitude decreases)
[GATE-1998] [GATE-1999]
Q.29 In a closed loop system, the process
to be controlled is an integrating process Q.33 A typical example of a physical
1 system with under-damped
with transfer function . The controller
2𝑠𝑠 characteristics is a
proposed to be used is an integral (a) U-tube manometer
1 (b) spring-loaded diaphragm valve
controller with transfer function .
𝑇𝑇1 𝑠𝑠 (c) CSTR with first-order reaction
When a step change in set point is applied (d) thermocouple kept immersed in a
to such a closed loop system, the liquid-filled thermowell
controlled variable will exhibit [GATE-1999]
(a) overdamped response
(b) underdamped response Q.34 A control system has the following
(c) undamped system transfer function,
(d) unstable response (𝑠𝑠 − 1)(𝑠𝑠 + 1)
[GATE-1998] 𝐹𝐹(𝑆𝑆) =
𝑠𝑠(𝑠𝑠 − 2)(𝑠𝑠 + 4)
The initial value of the corresponding time
Q.30 The Laplace transform integral of function is
exp(at), where a >0, is defined only for the (a) 1 (b) 1/8
Laplace parameter, s>a, since (c) 7/8 (d) -1
(a) the function is exponential [GATE-1999]
(b) the Laplace-transform integral of
exp(at) has finite values only for s>a Q.35 Which of the systems having the
(c) the Laplace-transform integral of following transfer functions is stable?
exp(at) has initial values only for s>a 1 1
(d) the function exp(at) is piece-wise (a) 2 (b) 2
𝑠𝑠 +2 𝑠𝑠 −2𝑠𝑠+3
continuous only for s>a 1 exp(−20𝑠𝑠)
[GATE-1999] (c) 2 (d) 2
𝑠𝑠 +2𝑠𝑠+2 𝑠𝑠 +2𝑠𝑠−1
[GATE-1999]
Q.31 A system with a double pole at the
origin is unstable since the corresponding Q.36 A proportional controller is used
term in the time for the control of a first-order system. If
(a) is a constant the dynamics of all the other units in the
(b) grows exponentially with time control loop are negligible and their
(c) grows linearly with time steady-state gains are all equal to unity,
(d) decays linearly with time show that
[GATE-1999] (a) The response of the controlled process
is faster than that of the uncontrolled
Q.32 A sinusoidal variation in the input process, and
passing through a linear first-order system
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 108
Gate Questions
(b) Then offset, both for changes in the set (c) The amplitude ratio of a pure capacity
point and in the load variable, decreases as process is inversely proportional to the
the parameter of the controller is frequency
increased. (d) The amplitude ratio of a pure time
[GATE-1999] delay system increases with frequency
[GATE-2000]
Q.37 A second- order process with the
transfer function Q.41 For a feedback control system to be
1 stable, the
𝐾𝐾𝑃𝑃 𝐺𝐺𝑃𝑃 (𝑠𝑠) = 2
𝑠𝑠 + 2𝑠𝑠 + 3 (a) Roots of the characteristics equation
Is controlled with a proportional-integral should be real
controller. If all the other lags in the (b) Poles of the closed loop transfer
control loop are negligible and their function should lie in the left half of the
steady-state gains are all equal to unity, complex plane
find the relation between the parameters (c) Bode plots of the corresponding open
to the controller that should be satisfied loop transfer function should
for the control system to be stable. monotonically decrease
[GATE-1999] (d) Poles of the closed loop transfer
function should lie in the right half of the
Q.38 The unit step response of the complex plane
2𝑠𝑠−1 [GATE-2000]
transfer function reaches
(3𝑠𝑠+1)(4𝑠𝑠+1)
its final steady state asymptotically after Q.42 The initial value (t = 0+ ) of the
(a) A monotonic increase
unit step response of the transfer function
(b) A monotonic decrease
(c) Initially increasing and then decreasing ( s + 1) / ( 2s + 1) is
(d)Initially decreasing and then increasing (a) 0 (b) ½
[GATE-2000] (c) 1 (d) 2
[GATE-2000]
Q.39 The unit step response of the
1 Q.43 The time constant of a unity gain,
transfer function 2
𝑠𝑠 +2𝑠𝑠+3 first order plus time delay process is 5 min.
(a) Has a non-zero slope at the origin If the phase lag at frequency of 0.2 rad/min
(b) Has a damped oscillatory is 60° , then the dead time (in minutes) is
characteristic (a) 5π /12 (b) π / 6
(c) Is overdamped (c) π /12 (d) π / 3
(d) Is unstable [GATE-2000]
[GATE-2000]
Q.44 An ideal PID controller has the
Q.40 Select the correct statement from transfer function 1 + 1 / ( 0.5s ) + 0.2s .
the following
(a) The frequency response of a pure The frequency at which the magnitude
capacity process is unbounded ratio of the controller is 1 is
(b) The phase lag of a pure time delay (a) 0.5 / 0.2 (b) 0.2 / 0.5
system decreases with increasing (c) 0.2 × 0.5 (d) 1/ 0.2 × 0.5
frequency [GATE-2001]
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 109
Gate Questions
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 110
Gate Questions
(a) 1 + 𝑒𝑒 𝑡𝑡 (b) 1 − 𝑒𝑒 𝑡𝑡 1
(c) (1 + 𝑒𝑒 +2𝑠𝑠 )
(c) 1 + 𝑒𝑒 −𝑡𝑡 (d) 1 − 𝑒𝑒 −𝑡𝑡 𝑠𝑠 2
1
[GATE-2004] (d) (1 − 𝑒𝑒 −2𝑠𝑠 )
𝑠𝑠 2
[GATE-2006]
Q.53 For the time domain function
f ( t ) = t , the Laplace transform of Q.56 A liquid level control system is
𝑡𝑡 configured as shown in the figure below. If
∫0 𝑓𝑓(𝑡𝑡)𝑑𝑑𝑑𝑑 is given by
3 the level transmitter (LT) is direct acting
(a) 1 / 2s (b) 2 / s3 and the pneumatic control valve is air-to-
(c) 1 / s3 (d) 2 / s2 open, what kind of control action should
[GATE-2004] the controller (LC) have and why?
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 111
Gate Questions
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 112
Gate Questions
Column-A
P- Control Strategy
Q- Primary Control Variable
R- Slowest Controller
S- Fastest Controller
Column-B
1-Feed Forward Control
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 113
Gate Questions
Q.66 The first two rows of Routh’s Q.67 The types of action for the two
tabulation of a third order equation are controllers are
s3 2 2 (a) direct acting for the pressure control
s2 4 4 and direct acting for the flow control
Select the correct answer from the (b) indirect acting for the pressure control
following choices and indirect acting for the flow control
(a) The equation has one root in the right (c) direct acting for the pressure control
half s-plane and indirect acting for the flow control
(b) The equation has two roots on the axis (d) indirect acting for the pressure
at s = j and ‒j. The third root is in the left control and direct acting for the flow
half plane. control
(c) The equation has two roots on the axis [GATE-2007]
at s = 2j and s = ‒2j. The third root is in the
right half plane. Q.68 The bias values for the two
(d) The equation two roots on the j-axis at controllers, so that no offset occurs in
s = 2j and s = ‒2j. The third root is in the either controller are
right half plane. (a) Pressure controller: 40%; Flow
[GATE-2007] controller: 60%
(b) Pressure controller: 33%; Flow
Linked Questions 67 & 68 & 69 controller: 67%
A cascade control system for pressure (c) Pressure controller: 67%; Flow
control is shown in the figure given below. controller: 33%
Thepressure transmitter has a range of 0 (d) Pressure controller: 60%; Flow
to 6 bar (g) and the flow transmitter range controller: 30%
is 0 to 81 nm3 / h . The normal flow rate [GATE-2007]
through the valve is 32.4 nm3 / h
corresponding to the value of set point for Q.69 Given that the actual tank pressure
pressure = 1 bar (g) and to give the flow, is 4 bar(g) and a proportional controller is
the valve must be 40% opened. The employed for pressure control, the
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 114
Gate Questions
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 115
Gate Questions
𝑡𝑡
−
(b) 2𝑒𝑒 2 − 𝑒𝑒 −𝑡𝑡
⁄
(c) 𝑒𝑒 −𝑡𝑡 − 2𝑒𝑒 −𝑡𝑡 2
(d) 𝑒𝑒 −𝑡𝑡 − 𝑒𝑒 −𝑡𝑡⁄2
[GATE-2009]
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 116
Gate Questions
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 117
Gate Questions
Q.88 The Bode stability criterion is oscillation 1 second and decay ratio 0.25.
applicable when The transfer function of the process is
(a) Gain and phase curves decreases 1
(a)
continuously with frequency 2
0.024𝑠𝑠 +0.067𝑠𝑠+1
(b) Gain curve increases and phase curve 1
(b) 2
decreases with frequency 0.67𝑠𝑠 +0.024𝑠𝑠+1
(c) Gain curve and phase curve both 1
(c)
increase with frequency 0.021𝑠𝑠 2+0.1176𝑠𝑠+1
(d) Gain curve decreases and phase curve 1
(d)
increases with frequency 0.1176𝑠𝑠 2 +0.021𝑠𝑠+1
[GATE-2012] [GATE-2012]
Q.89 The characteristic equation for the Q.92 A control valve, with a turndown
system is ratio of 50, follows equal percentage
𝑠𝑠 3 + 9𝑠𝑠 2 + 26𝑠𝑠 + 12(2 + 𝐾𝐾𝑐𝑐 ) = 0 characteristics. The flow rate of a liquid
Using Routh test, the value of K c that will through the valve at 40% stem position is
1 m3 / h . What will be the flow rate in
keep the system on the verge of instability
is m3 / h at 50% stem position, if the
(a) 20.9 (b) 18.4 pressure drop across the valve remains
(c) 17.5 (d) 15.3 unchanged? (Up to 2 digits after the
[GATE-2012] decimal point) _______?
[GATE-2012]
Q.90 The block diagram of a system
with proportional controller is shown Q.93 Consider the following transfer
below function
5
𝐺𝐺𝑝𝑝 (𝑠𝑠) =
(2𝑠𝑠 + 1)4
(Note: The unit of the process time
constant is in seconds.)
The crossover frequency (in rad/s) of the
process is
(a) 20 (b) 0.1
(c) 0.5 (d) 0.05
[GATE-2012]
A unit step input is introduced in the set Q.94 For the computation of Zeigler-
point. The value of K c to provide a Nichols settings, the ultimate period (in
critically damped response for U = 0 , s/cycle) and the ultimate gain are
τp =8 and τm = 1 is (a) π and 0.8, respectively
(b) 4π and 0.8, respectively
(a) 3.34 (b) 2.58 (c) 4π and 1.25, respectively
(c) 1.53 (d) 1.12 (d) π and 1.25, respectively
[GATE-2012] [GATE-2012]
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 118
Gate Questions
Q.95 For the time domain function Q.99 A step change of magnitude 2 is
f ( t ) = t 2 , which one of the following is introduced into a system having the
𝑡𝑡 following transfer function
the Laplace transform of ∫0 𝑓𝑓(𝑡𝑡)𝑑𝑑𝑑𝑑 ? 2
3 1 𝐺𝐺𝑝𝑝 = 2
(a) 4 (b) 2 𝑠𝑠 + 2𝑠𝑠 + 4
𝑠𝑠 4𝑠𝑠 The percent overshoot is __________?
2 2 [GATE-2014]
(c) 3 (d) 4
𝑠𝑠 𝑠𝑠
[GATE-2014] Q.100 Given below is a simplified block
diagram of a feed forward control system.
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 119
Gate Questions
[GATE-2015]
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 120
Gate Questions
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 121
Gate Questions
Assume the transfer functions of the Q.115 The number of positive roots of
measuring element and the final control the function f(x) shown below in the range
(Gm 1=
element are both unity= , Gf 1) . of 0 < x < 6 is ________?
The gain (rounded off to the first decimal
place) of the controller that will constitute
the critical condition for stability of the PI
feedback control system is _______?
[GATE-2016]
Q.114 For a unit step input, the response
of a second order system is
𝑦𝑦(𝑡𝑡)
= 𝐾𝐾𝑝𝑝 �1
1 �1 − 𝜁𝜁 2
𝜁𝜁𝜁𝜁
− 𝑒𝑒 − 𝜏𝜏 sin �
𝑡𝑡 + 𝜙𝜙��
�1 − 𝜁𝜁 2 𝜏𝜏
where, is the steady state gain, is the
damping coefficient, is the natural period [GATE-2017]
of oscillation and is the phase lag. The
πζ Q.116 The transfer function of a system
overshoot of the system is exp − is
1− ζ2
1
. For a unit step input, the response of the 4𝑠𝑠 2 + 1.2𝑠𝑠 + 1
system from an initial steady state For a unit step increase in the input, the
condition at t = 0 is shown in the figure fractional overshoot, rounded to 2
below. decimal places is_________?
[GATE-2017]
What is the natural period of oscillation (in Q.118 The Laplace Transform of a
seconds) of the system? function is
(a) 15.9 (b) 50 𝑠𝑠+1
(c) 63.2 (d) 100 𝑠𝑠(𝑠𝑠+2)
[GATE-2016] The initial and final values respectively, of
the function are
1
(a) 0 and 1 (b) 1 and
2
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 122
Gate Questions
1 1
(c) and 1 (d) and 0
2 2
[GATE-2017]
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 123
Gate Questions
[GATE-2019]
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 124
ANSWER KEY:
1 2 3 4 5 6 7 8 9 10
a b d d d a c c 81.70- 1, 0
81.90
11 12 13 14 15 16 17 18 19 20
* b a d a b b c a a
21 22 23 24 25 26 27 28 29 30
a * a d b a a b c b
31 32 33 34 35 36 37 38 39 40
c d a b a * * b b c
41 42 43 44 45 46 47 48 49 50
b b a d a * c c d d
51 52 53 54 55 56 57 58 59 60
d d c c a a d c a d
61 62 63 64 65 66 67 68 69 70
b a d d a b d d b a
71 72 73 74 75 76 77 78 79 80
a c d a d b a a a b
81 82 83 84 85 86 87 88 89 90
b c b d c b b a c c
91 92 93 94 95 96 97 98 99 100
a 1.47- c b d c d c 16.0- c
1.49 16.8
101 102 103 104 105 106 107 108 109 110
0.49- d c d 44-47 a 0.85- 3.5-3.7 a c
0.51 0.9
111 112 113 114 115 116 117 118 119 120
d 10.00- 2.49- a 3 0.3725 10 b 0.785 c
10.25 2.51
121 122 123 124 125 126
0.2 d 0.4137 a 43.375 b
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 126
Gate Questions
Q.34 (b)
Q.22 (*)
Open loop transfer function for the Q.35 (a)
1
process is
4(3𝑠𝑠+1)4 Q.36 (*)
Amplitude ratio for this system AR (a) The response of the uncontrolled first
1⁄4
= order process
�(3𝑤𝑤)2 + 1�(3𝑤𝑤)2 + 1�(3𝑤𝑤)2 + 1�(3𝑤𝑤)2 + 1 𝑦𝑦�(𝑠𝑠) 𝐾𝐾𝑃𝑃
1⁄4 │uc =
(a) AR = (3𝑤𝑤)2+12 𝑥𝑥̅ (𝑠𝑠) 𝜏𝜏1 𝑠𝑠 + 1
Phase lag, ϕ = −4 tan−1 (3𝑤𝑤) Now, the response of the controlled
Putting φ = −180º , we have process
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 127
Gate Questions
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 128
Gate Questions
Q.54 (c)
Q.43 (a)
Q.55 (a)
Q.44 (d)
Q.56 (a)
Q.45 (a)
Q.57 (d)
Q.46 (*)
(a) Yes Q.58 (c)
(b) Plot the open loop bode diagram Q.59 (a)
individually for all components in control
loop. Q.60 (d)
Then draw am overall cross frequency 𝑋𝑋1 = (𝑌𝑌𝑠𝑠𝑠𝑠 + 𝑦𝑦)𝐾𝐾𝑐𝑐
(𝜔𝜔𝑐𝑐 ) and its corresponding amplitude (𝐺𝐺1 𝑋𝑋1 + 𝐷𝐷)𝐺𝐺3 + 𝑋𝑋1 𝐺𝐺2 = 𝑦𝑦
ratio. 𝑋𝑋1 (𝐺𝐺1 𝐺𝐺3 + 𝐺𝐺2 ) + 𝐷𝐷𝐺𝐺3 = 𝑦𝑦
If α is the amplitude ratio then �𝑌𝑌𝑠𝑠𝑠𝑠 − 𝑦𝑦�𝐾𝐾𝑐𝑐 (𝐺𝐺1 𝐺𝐺3 + 𝐺𝐺2 ) + 𝐷𝐷𝐺𝐺3 = 𝑦𝑦
1 𝑌𝑌𝑠𝑠𝑠𝑠 (𝐺𝐺1 𝐺𝐺3 + 𝐺𝐺2 )𝐾𝐾𝑐𝑐 + 𝐷𝐷𝐺𝐺3
𝐾𝐾𝑐𝑐 max =
𝛼𝛼 = 𝑦𝑦(1
Based on Ziegler Nicholas controller + 𝐾𝐾𝑐𝑐 (𝐺𝐺1 𝐺𝐺3 + 𝐺𝐺2 ))
setting for a PID controller 𝑦𝑦(𝑠𝑠) 𝐺𝐺3
𝑇𝑇𝑇𝑇 = 2𝜋𝜋⁄𝜔𝜔𝑐𝑐 =
𝐷𝐷(𝑠𝑠) (1 + 𝐾𝐾𝑐𝑐 (𝐺𝐺1 𝐺𝐺3 + 𝐺𝐺2 ))
𝐾𝐾𝑐𝑐 = 0.6𝑘𝑘𝑐𝑐,𝑚𝑚𝑚𝑚𝑚𝑚 ; 𝜏𝜏1 = 𝑇𝑇𝑇𝑇⁄2
𝜏𝜏𝐷𝐷 = 𝑇𝑇𝑇𝑇⁄8
Q.61 (b)
Q.47 (c)
Q.62 (a)
1
Q.48 (c) 𝐺𝐺1 (𝑠𝑠) = , 𝐺𝐺2 (𝑠𝑠)
𝜏𝜏1 𝑠𝑠 + 1
1
Q.49 (d) =
𝜏𝜏2 𝑠𝑠 + 1
Q.50 (d) 𝜏𝜏2 = 2𝜏𝜏1
1
𝐺𝐺2 (𝑠𝑠) =
Q.51 (d) 2𝜏𝜏1 𝑠𝑠 + 1
The closed loop transfer function is X (=
t) sinωt and
𝐾𝐾𝑐𝑐
𝐺𝐺 = 1 + =Y A sin ( ωt + φ )
(4𝑠𝑠 + 1)(2𝑠𝑠 + 1)(𝑠𝑠 + 1)
ϕ = tan−1 (−𝜔𝜔𝜏𝜏1 ), 𝜏𝜏1 = 1
The characteristic equation is
(4𝑠𝑠 + 1)(2𝑠𝑠 + 1)(𝑠𝑠 + 1) + 𝐾𝐾𝑐𝑐 = 0 ϕ = tan−1 (−𝜔𝜔)
ϕ = 45, ω=1
Put 𝑠𝑠 = 𝑗𝑗𝑗𝑗 and equate real and imaginary
part to solve for 𝜔𝜔 and𝐾𝐾𝑐𝑐 .
Q.63 (d)
1 1
Q.52 (d) 𝐴𝐴𝐴𝐴1 = =
�1 + 𝜏𝜏12 𝜔𝜔 2 √1 + 12 12
Q.53 (c) 1
=
√2
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 129
Gate Questions
1 1
𝐴𝐴𝐴𝐴2 = =
�1 + 𝜏𝜏22 𝜔𝜔 2 �1 + (2𝜏𝜏1 )2 𝜔𝜔 2 Q.74 (a)
1
= Q.75 (d)
�1 + (2.1)2 + 12
1
= Q.76 (b)
√5
1
Amplitude ratio for L = AR1. AR2= Q.77 (a)
√2
1 1 1 1
. = = 0.32 2
=
√5 √10 2𝑠𝑠 + 3𝑠𝑠 + 1 (2𝑠𝑠 + 1)(𝑠𝑠 + 1)
2 1
Q.64 (d) = −
2𝑠𝑠 + 1 𝑠𝑠 + 1
1
Q.65 (a) 𝐿𝐿 � 2 �
2𝑠𝑠 + 3𝑠𝑠 + 1
2 1
Q.66 (b) = 𝐿𝐿 � � − 𝐿𝐿 � �
2𝑠𝑠 + 1 𝑠𝑠 + 1
Characteristic equation correspond to the = 𝑒𝑒 −𝑡𝑡 ⁄2 − 𝑒𝑒 −𝑡𝑡
above criteria is 2s2 + 4s2 + 2s + 4 =
0
(2s + 2) ( s + 2) =
2
0 Q.78 (a)
s2 + 1 = 0, s + 2 = 0 Q.79 (a)
s = ± j, s = -2 R ( s ) = Y ( s ) = E ( s ) ; 2 × E ( s ) = 2U ( s )
So the equation has two roots on the axis
2
at s = j and ‒j. The third root is in the left 𝑈𝑈(𝑠𝑠) ∙
half plane. (2𝑠𝑠 + 1)(3𝑠𝑠 + 1)
𝑈𝑈 (𝑠𝑠)
= 𝑌𝑌 (𝑠𝑠), 𝑅𝑅(𝑠𝑠) − 𝑌𝑌 (𝑠𝑠) =
Q.67 (d) 2
𝑌𝑌(𝑠𝑠)
=
Q.68 (d) 2
2×(
2𝑠𝑠 + 1 3𝑠𝑠 + 1)
)(
Q.69 (b) 𝑄𝑄(𝑥𝑥 ) = 𝑆𝑆𝑆𝑆𝑆𝑆 + 𝐶𝐶
𝑄𝑄(𝑥𝑥 ) = 0 𝑎𝑎𝑎𝑎 𝑥𝑥 = 0 ⇒ 𝐶𝐶 = 0
Q.70 (a) 𝑄𝑄(𝑥𝑥 ) = 𝑆𝑆𝑆𝑆𝑆𝑆
𝑄𝑄(𝑥𝑥)
� = 𝑆𝑆𝑆𝑆
Q.71 (a) 𝐴𝐴 𝑥𝑥=𝐿𝐿
𝑦𝑦(𝑠𝑠) 𝑘𝑘𝑝𝑝
= , 𝑋𝑋 (𝑠𝑠) = 1.0, 𝑌𝑌(𝑡𝑡) Q.80 (b)
𝑥𝑥(𝑠𝑠) 𝜏𝜏𝜏𝜏 + 1
𝑘𝑘𝑝𝑝 −𝑡𝑡
= 𝑒𝑒 𝜏𝜏 , 𝑌𝑌(𝑡𝑡) Q.81 (b)
𝑧𝑧
The given function can be written in the
= 2𝑒𝑒 −0.5𝑡𝑡 ,
𝜏𝜏 = 2, 𝑘𝑘𝑝𝑝 = 4 unit step functions
F ( t=) V (U(t − a) − U(t − b))
Q.72 (c) 1
𝐿𝐿�𝑈𝑈(𝑡𝑡 − 𝑎𝑎)� = 𝑒𝑒 −𝑎𝑎𝑎𝑎
𝑠𝑠
Q.73 (d) 𝑉𝑉 −𝑏𝑏𝑏𝑏
𝐿𝐿�𝑓𝑓 (𝑡𝑡)� = (𝑒𝑒 − 𝑒𝑒 −𝑎𝑎𝑎𝑎 )
𝑠𝑠
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 130
Gate Questions
𝐾𝐾𝑐𝑐,2 + 1.5
� , For the system to be stable,
Q.82 (c) 2𝐾𝐾𝑐𝑐,2 − 6
From the figure, the transfer function is of 𝐾𝐾𝑐𝑐,2 = 3⁄2
the type 𝐴𝐴(𝑇𝑇𝑇𝑇 + 1)
For the above, the amplitude ratio is Q.85 (c)
𝑙𝑙𝑙𝑙𝑙𝑙(𝐴𝐴𝐴𝐴/𝐴𝐴) = 𝑙𝑙𝑙𝑙𝑙𝑙√𝜏𝜏 2 𝜔𝜔 2 + 1 𝑝𝑝(𝑡𝑡) = 30 + 5𝑒𝑒(𝑡𝑡)
𝑡𝑡
We will get the values of A and T by
plotting Ar vs ω curve for τω < 1 and > + 1.25 � 𝑒𝑒(𝑡𝑡)𝑑𝑑𝑑𝑑
0
1. 𝑑𝑑𝑑𝑑(𝑡𝑡)
𝐴𝐴𝐴𝐴 + 15
When τω < 1 , 𝑙𝑙𝑙𝑙𝑙𝑙 � � = 𝑙𝑙𝑙𝑙𝑙𝑙√1, 𝐴𝐴𝐴𝐴 = 𝐴𝐴 𝑑𝑑𝑑𝑑
𝐴𝐴
𝑝𝑝(𝑡𝑡) − 30 = 𝑂𝑂(𝑡𝑡)
constant part of the figure which is 10.
→ 𝑝𝑝(𝑡𝑡) = 5𝑒𝑒(𝑡𝑡)
τ =1 / 10 because till ω = 10 it is having 𝑡𝑡
Ar =10. + 1.25 � 𝑒𝑒(𝑡𝑡)𝑑𝑑𝑑𝑑
0
Q.83 (b) 𝑑𝑑𝑑𝑑(𝑡𝑡)
+ 15
Given range is on the scale, 𝑈𝑈(𝑠𝑠) = 𝑑𝑑𝑑𝑑
2(1+0.2𝑠𝑠) 1.25
𝐸𝐸(𝑠𝑠) 𝑝𝑝(𝑠𝑠) = 5𝑒𝑒(𝑠𝑠) + 𝑒𝑒(𝑠𝑠)15𝑠𝑠. 𝑒𝑒(𝑠𝑠)
1+0.02𝑠𝑠 5
𝑌𝑌(𝑠𝑠) 𝑃𝑃(𝑠𝑠) 5 12𝑠𝑠 2 + 4𝑠𝑠 + 1
2 𝐺𝐺(𝑐𝑐) = = � �
= 2 𝑈𝑈(𝑠𝑠), at the junction, 𝑅𝑅 (𝑠𝑠) 𝑒𝑒(𝑠𝑠) 4 𝑠𝑠
𝑠𝑠 + 11𝑠𝑠 + 1
∗ 0.2 − 0.2𝑌𝑌(𝑠𝑠) = 𝐸𝐸(𝑠𝑠)
Q.86 (b)
𝑅𝑅(𝑠𝑠) ∗ 0.2 − 0.2𝑌𝑌(𝑠𝑠)
Characteristic equation
2
=� 2 1 + 𝐺𝐺(𝑐𝑐) . 𝐺𝐺(𝑠𝑠)
𝑠𝑠 + 11𝑠𝑠 + 1
−1
10 5 12𝑠𝑠 2 + 45𝑠𝑠 + 1
2(1 + 0.2𝑠𝑠) 1+ . � �
∗ � 𝑌𝑌(𝑠𝑠) (200𝑠𝑠 + 1) 4 𝑠𝑠
1 + 0.02𝑠𝑠 = 700𝑠𝑠 2 + 102𝑠𝑠 + 25 = 0
𝑌𝑌(𝑠𝑠) 0.8(1+0.2𝑠𝑠)
On solving → = =
𝑅𝑅(𝑠𝑠) 0.2𝑠𝑠 3+3.2𝑠𝑠 2+11.36𝑠𝑠+1.8
0.8⁄1.8(1+0.2𝑠𝑠) Q.87 (b)
0.2⁄1.8𝑠𝑠 3+3.2⁄1.8𝑠𝑠 2+11.36⁄1.8𝑠𝑠+1 𝑓𝑓(𝑡𝑡) = 𝑒𝑒 −𝑡𝑡 ⁄𝜏𝜏
Steady state gain = 0.8/1.8 = 4/9 By Taylor’s expansion, 𝑒𝑒 𝑥𝑥 = 1 + 𝑥𝑥 +
𝑥𝑥 2
+
2!
𝑥𝑥 3
Q.84 (d) 3!
+⋯
𝑅𝑅(𝑠𝑠) − 𝑌𝑌(𝑠𝑠) = 0.5𝑈𝑈1 (𝑠𝑠), 𝑡𝑡 𝑡𝑡 2 𝑡𝑡 3
𝑈𝑈1 (𝑠𝑠) − 𝑌𝑌(𝑠𝑠) 𝑒𝑒 −𝑡𝑡 ⁄𝜏𝜏 = 1 − + 2 − 3 + ⋯ 𝑡𝑡 ≪ 𝜏𝜏
𝜏𝜏 2𝜏𝜏 6𝜏𝜏
= 𝐾𝐾𝑐𝑐,2 𝑈𝑈2 (𝑠𝑠) 𝑎𝑎𝑎𝑎𝑎𝑎 𝑈𝑈2 (𝑠𝑠) 𝑒𝑒 −𝑡𝑡 ⁄𝜏𝜏
𝑆𝑆 + 2 =1
= 𝑌𝑌(𝑠𝑠)
𝑆𝑆 − 3 𝑡𝑡
Substituting U2(s) in equation 2 and − , as all higher expressions of τ can be ignored
𝜏𝜏
U1 ( s ) in equation 1
𝑅𝑅(𝑠𝑠) 0.5 ∗ (𝑠𝑠 − 4) Q.88 (a)
=
𝑌𝑌(𝑠𝑠) (𝐾𝐾𝑐𝑐,2 + 1.5)𝑠𝑠 + (2𝐾𝐾𝑐𝑐,2 − 6)
Q.89 (c)
Routh Stability Criteria
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 131
Gate Questions
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 132
Gate Questions
Q.116 0.3725 4 1
1 𝑌𝑌(𝑠𝑠) = ×
7𝑠𝑠 + 5 𝑠𝑠
⇒ 𝜏𝜏 2 = 4 ⇒ 𝜏𝜏 = 2 4
4𝑠𝑠 2 + 1.2𝑠𝑠 + 1 𝑌𝑌(∞) = lim 𝑌𝑌(𝑡𝑡) = lim 𝑠𝑠𝑠𝑠 (𝑠𝑠) =
1.2 1.2 𝑡𝑡→∞ 𝑠𝑠→0 5
2𝜏𝜏𝜏𝜏 = 1.2 ⇒ 𝜀𝜀 = = = 0.3 𝑋𝑋(∞) = 1
2𝜏𝜏 2 × 2
−𝜋𝜋𝜋𝜋 4 1
Fractional overshoot, = exp � � = 1 − = = 0.2
√1 − 𝜀𝜀 2 5 5
−𝜋𝜋(0.3)
= exp � � = 0.3725 Q.122 (d)
√1 − 0.32
K = 10 = 0.4137
© Copyright Reserved Gateflix. No part of this material should be copied or reproduced without permission 133