PC Text Book
PC Text Book
Automatic
Process Control
Donald P. Eckman ....
Professor of Irntrurnmlaiion
Engineering
Case Institute of Technology
1, 5, 3, and/or 0, 7, 2, and 4.
1, Appendix, 4, 9, 10, 5, 7, 8.
Donald P. Eckman
Clet'eland, Ohio
April, 1958
appendix 325
index . . . 364
vii
Variables Elements
Set point A Input elements
Reference input Gx Control elements
Actuating signal G2 System elements
Manipulated variable H Feedback elements
Controlled variable N Load elements
Feedback variable
Load variable
chapter
The Science
of Automatic l.ontrol
tion of automatic control ideas to physical fiystems that only a few years
ago were impossible to analyze or control
A clear concept of the principle of automatic control benefits the
engineer because, as pointed out by J. Von Neuman^, ‘There is no point
in using exact methods where there is no clarity in concepts and the
issues towhich they are to be applied.*^ Furthermore a need for auto-
matic control has often resulted in a careful and systematic analysis of
the physical part of the process.
Automatic control devices employ the principle of feedback within
themselves in order to better their performance and are therefore a
subject of study. Although this book will discuss automatic control
from a general standpoint, most of the examples are selected from process
control. Some discussion of the principle of operation of control mech-
anisms is presented in Chapters 6 and 7.
breaks the electric circuit. This action de-energizes the solenoid valve
which turns off th5 gas to the furnace.
The most general type of automatic control system is shown in the
u ^ U (1-3)
use of automatic control, and it does seem incredible that in the ancient
water works of Home, (J recce, or Egypt no recorded use of closed-loop
water regulating devices was made. Even if there was, there is no dis-
cussion of the philosophy of automatic control. j
The first analysis of automatic control is the mathematical discussion
of the flyball governor by James Clerk Ma.xwell in 1868. A set of
historical references is provided at the end of the chapter. Further
application of the governor techni<iue.s to other engines and turbines
was made, and in the early 1900’s the application to process control
began. At the same time regulators and servomechanisms were being
studied for their application to steam-power regulators and ship
steering
The first general theory of automatic control was founded by Nyquist
in the famous article on “Hegeneration Theory.” This study provided
the basis for the determination of stability of systems without having to
solve com()lot('ly the differential equations. Further developments in
electric mnplifiers and electric servomechanisms provided the many
locus and fre(|iiency techniques in use today.
The general applications of process control did not begin until the
193()’s. The usefulness of control techniques quickly established their
value, so that by the 1940’s rather complex control networks were in
common use.
Automatic control devices are used in almost every phase of indus-
trial operations. They are commonly employed in
The Science of Automatic Control 7
Instrumentation
cybernetics
Communication Control
and and
information feedback
theory theory
upon the viewpoint of the analysis. The diagram on the left is drawn
10 Automatic Process Control
as though the outflow is the important variable being influenced by
inflow. The feedback occurs because an increase in outflow feeds back
to decrease head h. At the same time an increase iH inflow will increase
the outflow. Thus, the relationship of variables may be written as
q = /(A, t) (1-5)
h = g(m - q, t) (1-6)
g - q, 0], (1-7)
The important point is that the block diagram illustrates the relation of
these variables.
The block diagram on the right in Fig. 1-10 is for the same process
but is shown to head in the vessel
illustrate that inflow influences the
and that the outflow is the secondary variable. This diagram does not
relate one to one to the mass flow because tlie input is flow and the
output is head. Recombining equations 1-5 and 1-6, the overall
relationship is
h^gWm -f(hj)lt\ (1
~8 )
There are only a few rules for the use of block diagrams:
HISTORICAL REFERENCES
Tht* followiriK refereiK’es are selected for the contribution each makes to the
knowledpje of automatic control. Prior to 1900 there are many articles on governors,
but most of them provide consitleration of iiidi\ idual detailed problems.
1808 J. C. Maxwell, “On Cov^crriors/’ Froc Hoyal Soc. London, \ol. 16, p. 270.
1893 A. Stodola, “Uber die Itegulierung von I'uibineii," Schweiz. Bauztg., Vol. 22,
p. 117.
1907 A. M. Liapounov, “Probl^me g^n^ral de lu stabihU* du inouvement,” Ann.
fac. set. univ. Toulouse
1911 It von Misea, “Itegulierung dea Maachinemganges," Encyl. der Math. ITiss.,
PROBLEMS
1-1. Give examples of the use of automatic control devices for (a) reducing
the disturbances due to load variables in a system; (6) human safety.
» W, Name and describe briefly several automatic control devices found
in the home.
1-3. Name and
doKcnbe briefly several instruments found in the home.
1-4. Sketch a block diagram for a man in a boat with an outboard motor
performing the operation of steering for a point ashore.
1-5. Describe a device that could be employed to maintain the constant
speed of an automobile.
Process *•
Characteristics
Process Variables
The load variables of the process are all other independent \'ariables
except the controlled variable and the manipulated variable. For the
water heater of Fig. 2-1, the incoming water temperature v is the load
variable. It is expected that the automatic controller will (ujrrect for
n = Hr - Ue (2-1)
In the example of the billiard ball there are three variables of position,
one defining equation (height = constant), and therefore two degrees
of freedom.
Similarly, a process has a finite number of degrees of freedom. Con-
sider the heat exchanger of Fig. 2-1, There are four variables:
w = inlet temperature
c = outlet temperature
w = water flow rate
rn — heat input rate
• n = Ue - rip +2 (3-2)
This applies only to vhe chemical states of the process, and the number
2 in the equation atJove represents temperature and pressure. For an
isothermal process
n = Tic - rip + I (2-3)
w = w, - Tip + I (2-4)
Example 2-1. P\)r the liquid-to-liquid heat exchanger shown in the figure,
the following variables are given below the diagram:
He = 1
?i = () — ] — 5
V rir
~ 1
The number of components me two. and the number of phases are two. There-
fore
?? - 2 — 2 d- 1 - 1 (temperature v)
For the process, the number of variables are 1hre(‘ (7, r. and ?«), and the defin-
ing relations.'in* two (conservation of eiKTgy and mass), therefore the ni*niW‘r
of degrees of freedom is
n - n,. — rt, -- — 2 --
1
Process Cliaracteristios 17
Cl = overhead temperature Ui
~ bottom flow rate
cj = overhead pressure Ui feed temperature
C| = overhead coi^osition u« = feed pressure
a « overhead flow rate uj * feed composition
ui * bottom temperature u% =* feed per cent vapor
ui - bottom pressure Ui - feed flow rate
uz - bottom composition m = steam flow rate (heat input)
Note that the process may contain internal heat exchange or any other
arrangement of equipment. Also the pressures 12, and Ur are generally
the same. Calculate iFie number of degrees of freedom. What are the maxi-
mum number of automatic controllers that may l>e used?
Employing Gibbs’s phase rule at the overhead there are two com})onents
and two phases'
n - 2 - 2 -h 2 = 2
so that any two of the three variables (ci, r^, c^) are independent. The same
IS true at the bottom where any two of the three variables (ui, Ui, uz) are
independent. For the feed also any two of the three variables (ws, ur, uj)
are independent.
The degrees of freedom may now be calculated. The number of variables
are eleven
Ovei head— any two of ci, a, c^; and
—
Hottoin any two of wi, uz, and w.4
Heat input -"W
Feed --any two of 1/5, ut, uj; and Ug, uv
The basic concepts for system analysis in automatic control are llie
fundamental laws relating the behavior of the system. Statical or
steady-state behavior is important, but even more important in auto-
matic control is the dynamic behavior of mechanical, fluid, thermal,
and electric .‘systems. The important characteristics of fluid and thermal
processes are their resistance and capacitance.
Electric systems composed of pure resistance elements and jiure
capacitAiice elements in various combinations have many character-
istics in common with fluid and thermal processes Consequently a
study of electric systems serves to clarify many similar ideas in fluid
and thermal systems
The flow law for electric conductors (Fig. 2-2 ) having resistance but
no inductance or capacitance is Ohm's Law,
fi — to = iR ( 2- 5 )
Automatic Process Control
than about 4000, and laminar flow for which the Reynolds number is
lessthan about 20^. fn turbulent flow through pipes, orifices, valves,
•. I
R, resistance C, capacitance
q = KAV2g(hi - h2 ) (^9)
This law may be compared to Ohm’s law for electric conductors equa-
tion 2-5, the basic difference being that it involves the square root of
the head or potential.
The turbulent resistance is found from
resistance not constant but depends upon the flow rate and head
is
F ia the friction factor, L is the equivalent pipe length, and D is the internal diameter
of the pipe.
II
20 Automatic Process Control
Example 2^4. The flow verBUs differential head for a 1-in. orifice is
Flow, q, ft^/sec
The resistanro may be calculated at any particular head and flow and will
vary directly with flow. The resistance at any head may be calculated from
^ ^
128mL mvL
.tB-
where h — head, ft
sec/ft^ (2-12)
dq
4^ = 9
at
(2-13)
t = time, sec
'Fhis law should be compared to the charging law for electric capacitors
(equation 2 7).
Liquid capacitance is defined bv
C = (2-14)
. dh
where »
tr gas flow rate, Ib/sec
=
/iC a flow coefficient
A = area of restriction, ft*
Y - rational expansion factor
7 = gas density, Ib/ft*
p = pressure Ib/ft*
R = ^
dw
sec/ft*
**
(2-16)
R, resistance C, capacitance
Process Characteristics 23
—
V
= constant ( 2 - 19 )
'
R = gas constant for a specific gas, ft/deg
T = temperature of gas, deg
C =
V 20 X 0.134
= 0.000076 ft*
fiRT 1 X 63.3 X 660
Employing moles and the universal gas constant R,
vVi
dp"" RT
29 X 20 X 0.134 =
C = 0.000076 ft*
1646 X 660
Thermal systems involving heat transfer from one substance to another
are also characterized by resistance and capacitance. In thermal
systems it will be assumed that substances characterized by resistance
to heat flow have negligible storage of heat (as a very thin air film),
and substances characterized by heat storage have negligible resistance
to heat flow (as a block of copper or aluminum).
R, resistance C, capacitance
The fljoxv laws for heat conductors (Fig. 2-5) are of two diflerent types
corresponding to conduction or convection, and radiation. For conduc-
tion of heat through a specific conductor, the heat flow is given by the
Fourier law which reduces to
KA
TZiOi -h) (2-2\)
AX = thickness of conductor, ft
6 = temperature, deg
q^^HAiSi-do) ( 2- 22 )
„ dS AX
—
« = -- = ,
or
R= —
Jiq
= —
HA
deg sec/Btu (convection) ( 2* 24 )
dim ^
(AA’)
^
~
0.125 3600
0.65 deg sec/Btu
,lq
~ KA 12 ^ 29 X 2
7 = KAEidi^ - $2*) ( 2 - 25 )
E= emissivity
A - suiface area, ft^
0 = temperature, deg R
This is a fourth-power law as compared to a square-root relation for
turbulent fluid flow, and a linear law for electric current.
The radiation resisinnee is
(je
R ( 2 -26 )
dq
d6 1
^
dq 4KAEej * 4 X 47.4 X 1.0 X 20J10®
fl - 60 deg sec/Btu
26 Automatic Process Control
C^
at
= q (2-27)
C= TFCp (2-28)
Capacitance —
Coul.
volt
TT * fstrad Btu/deg ft*
volt
Resistance = ohm ,
deg-sec/Btu sec/ft* sec/ft*
amp
Process Characteristics 27
Mechanical systems may be composed of various masses connected
by members providing damping effects and spring effects. Although
some systems are analogous to electric systems, the analog comparison
may be made in many ways. Consequently, it is better to consider
mechanical systems on their own merits.
dt\ di)
= M dt^
(2-29)
where / = force, lb
i = time, sec
M = mass, lb«secVft
X = displacement, ft
/-«! (MO)
/ = /iTx (2-31)
Block diagram
1 Proportional element
2. Capacitance element.
3 Time-constant element
4 Oscillatory element
Process Characteristics 29
These are the '^building blocks’" of processes and almost all industrial
processes are described by arrangements of numbers of these elements
in various series and parallel arrangements.
The 'proportional element is described in Fig. 2-7. Consider, for
example, the capillary shown in the physical diagram. The flow rate
of liquid m through the capillary is regarded as the variable being
(‘hanged and thus is termed the input variable. The head c is the
!-csult variabl^Jow rate and is regarded as the output variable.
of the
A capillary constitutes a laminar resistance and the flow-head equation
may be written
c = Rm (2-32)
rhus, the output is determined from the product of the system function
;ind the input. In the example of the capillary the system function
is the resistance of the capillary. Note that the system function is a
characteristic of the system and is given in terms of the physical
characteristics
Other proportional elements are electrical resistance, gas flow resist-
dc
C— = m (2-34)
where C - capacitance
c = output variable (head)
I -- time
m= input variable (flow)
30 Automatic Process Control
c, head
,
flow
Physical diagram
1 c
i
Cs
Block diagram
(Cs)c = m (2-36)
or
(2-37)
Time, t —
Block diagram
This equation is plotted in Fig. 2-8. For a step change of input flow M,
30 Automatic Process Control
In order to obtain the Bystem function the operational notation of differ-
d
{2-35 )
Physical diagram
1 c
Ca
Block diagram
(2-37)
fit, flow
Physical diagram
R c ^
Ts+l
Block diagram
'-si'"'"-©'
This equation is plotted in Fig. 2-8. For a step change of input flow M,
32 Automatic Process Control
the output variable increases linearly with time with a rate of change
inversely proportional to capacitance.
The time-constant element is shown in Fig. 2-9 fiftid is exemplified by
the liquid tank and resistance. The input variable is the inflow rate m
and the output variable is the tank head c. For the tank capacitance C,
C- fit
= m-q (2-39)
(i
= (2-40)
R
Combining equations 2-39 and 2-40 to eliminate outflow q,
lie
KC~->tc=Rm (2-^1)
fit
c = RM {I - (2-^3)
IS given by the sum of the particular mlegral and complementary function. The
particular integral im
X =: By
if y 16 a constant. The complementary function is the solution of
Process Characteristics 33
The time (T = RC) is the time constant of the element and is found
experimentally when I = T,
M -• -B ~ - Kc + rn (2-44)
where the force m. is regarded as the input variable and the displacement
of themass is eonsidered the output variable. Rewriting equation 2-44
and in
r = Ce
0
and solving for X,
X = 1/A
The complete solution is therefore
T ^ By + Ce-"^
The arbitrary constant C is found from initial conditions.
34 Automatic Process Control
in operationnl form,
(3 -«)
Time, I
u. e, dispiacement
1
m, force-^
Physical diagram
m l/JT e
TV + 2rTs+l
Block diagram
m (2-48)
+ 2[Tt + 1 /
The system function is typical of oscillatory elements when the damp-
ing ratio { is less than one (underdamped).
The response of an oscillatory element is shown in Fig. 2~10 for a
unit step change of mput variable. Solving equation 2-48, under the
conditions that the system is initiaUy at rest,^
c-^[l--^L^8invT^l-e-f‘/»’coeVr^i;j (iM9)
is given by the sum of the particular integral and complementary function. The
particular integral is
X ^ Cy
if V is constant. The complementary function is the solution of
d*j dx
-f B d= \/jg» - 44
Xi,Xj
24
The form of the complementary function depends upon whether > 44, B* * 44,
or < 44 In the latter case the roots X are complex and the complete solution is
.
I = Cy +
Employing Euler’s relation,
c * e cos Xf t sin Xf
where K% and Ka are arbitrary constants which must be evaluated from initial con-
ditions.
:
1 l/K
(2-51)
Ms- 4- Bs -f K (Tis H- 1){T2S + 1)
J_
B
(2-52)
7'., 2M
'fherefore a second-order system with overdamped action > \KM)
may be factoied into two time-constant elements If the system is
critically damped these two time-constants are et|ual
The response* of an oscillatory elemi'iit always overshoots the final
value and comes to Hist practically at a time dependent upon the
amount of damping and the characteristic lime. The amount of over-
shoot is also greatly dependent upon the damping The frequency of
oscillation depends primarily upon the characteristic time T
Liquid Processes
Proc(‘sses involving the How ot liquuls through connected pipes and
\cssels are common cwamples ol industrial processes and are often
employed for demonstrating the operation of many different kinds of
processes.
The single-vessel process of Fig’ 2-1 1 has a constant outflow u but
a variable inflow ui The inflow m assumed to he the manipulated
is
Cc = m — u (2-53)
Process Characteristics 37
The block diagram foi tlie piocess of Fig. 2 11 may he drawni in two
ways as shown in Fig. 2 12 These two diagrams produce identical
relationships hetw'een the variables of the system.
The single vessel of Fig. 2-13 in which the outflow is through a liquid -
Ct m^ (2‘-55)
qo = outflow
9o = ^(c-u) (2-56)
Tc + c = Rm -f u (2-57)
where T => RC is the time constant of the process. The system equa-
tion is
considered constant.
The continuity relation for the upper vessel is
C 2A 2 = m -f- W2 ~ (2-59)
Process Characteristics
(2-62)
C\c = + Ml - qo (2-63)
where Ti = RiCi and is the time constant of the lower vessel. The
latter equation may be placed in operational form
the right-hand set of two blocks and circle. The process is thus seen
R, .
Ri
(2-C7)
{TiS+\)iTzs+\) (7’,.s-+1)(7’,a+1) (7VTi)
and the o[)erationul equation lor the uxeiall .system j.s obtained. The
block diagram for the overall system is shown in Fig, 2-15 where com-
parison may be made to the previous block diagram.
For two time-constant elements in series, the response is found by
solving the differential equation for the overall .system. If the load
Process Characteristics 41
variables u\ and U 2 are zero (flows turned off), equation 2-67 becomes
Fio. 2-16, Response of processes with several equal time couslante in series.
where c/ = Him is the final steady value of head in the lower vessel
and Co is the initial value of head in the lower tank. Another solution is
Example 2-9. A single capacitance process like that of Fig. 2-13 has
a normal operating head of 4 ft and a normal value of outflow of 0.125 cu ft
per sec. The cross-section area of the vessel is 6 sq ft. •Calculate the time
constant.
If the resistance to flow is parabolic, the resistance is given by
C
(-)
\ Q /normal
eft*
“
^ = 64 sec/ft*
Example 2-10. Show that the time constant of a single vessel as shown
in Fig. 2-13 is proportional to the time required to change the fluid in the
vessel.
The capacitance is the area of the vessel and the resistance of the outlet
valve depends upon head and flow. For turbulent resistance,
X A - 2 2
9o outflow
T RC -
•V.
XA , = hA— Qo
volume of fluid
outflow
Process Characteristics 43
1 (Ti + T + RiCi)
2 =b V{Ti ^Ti + RiC2y - iTiTt
2T1T2
Gas Processes
Processes involving the flow of gas through connected pipe lines and
pressure vessels are conunon in industrial practice.
|m
Cc — u\ — Wo (2-71)
“ /(r, (2-72)
"•-©'Hi)
where Mo is a constant of integration.
The inflow (ic,) is a function of three variables, the stem position m
of the inlet valve, the upstream pressure wi, and the vessel pressure c,
so that
tCa = f{m, Ui, c) (2-75)
The partial derivatives may be evaluated as for the outlet \'ul\'e witli
the addition of one more plot of flow iat(' against valve stem position
46 Automatic Process Control
Combining the equations 2-71, 2-74, and 2-77 the Tproce^s eqiwJtUm is
or in operational form
“ “+
'
(t^) ““ + (?^)
where
y
C X
\/Ra + 1/^4 *
^/Ra + V®4
1
,j. ^ i
* So/fta + ®
fil/fl. + fil/Kk
The block diagram is shown in Fig. 2-17. The process response is first
14.7 297
== 27.8, 56.0
0.53 0.53
Since 27.8 — 14.7 < 15 psig vessel pressure the outlet valve has
== 13.1 psi
critical (acoustic) Also because 56.0 — 14.7 = 41.3 < 45 psig, the
flow.
inlet valve has critical flow. Therefore /2o ® and and thereby “
Aj =* 0 in equati<|n 2-79. The time constant is
15 X 144 10
= min
RaC ^ 0.76
1.0 1.0 X 53.3 X 530
Ki RoK^
Ha
Kz - 2.0
Ri
c = Ui
is the vessel pressure c psi gage. The stem position of the outlet valve is m
the manipulated variable. The outlet valve downstream pressure is less than
critical. Derive the process equation.
Prooess Characteristics 47
The continuity relation is
Ct »
The outflow does notsdepend upon downstream pressure, so
or
Wa * Kftn -h c + Mq
The process equatioi^ is
Td -h c = -RK^m + RMo
and the operational equation is
where T = RC, The resistance R of the outlet valve is the slope of the
upstream pressure versus flow curve at constant valve stem position and
is usually obtained by test. In drawing the block diagram, constants (such
as RMo) are usually omitted because they do not directly influence process
dynamics.
Flow Processes
The automatic control of fluid flow in pipelines would at first appear
to be a simple problem of manipulating a throttling valve so as to
maintain the flow rate at the desired value. In practice, the analysis
involves a large number of factors related to flow characteristics, and
a separate consideration of flow-rate processes is necessary. The analysis
of flow-rate processes will be made by assuming no effect due to inertia
of flowing fluid. This assumption is quite close to reality, because in
practical applications the changes in variables occur relatively slowly.
In addition it is assumed that the capacitance of pipelines is negligible.
Liquid-flow problems will be discussed first, and then gas-flow problems
will be reviewed.
Liquid flow rate is assumed to be measured by some type of metering
device such as an orifice, nozzle, or Venturi as shown in Fig. 2-19. The
flow rate through the metering device is given by
c JiAiV2g{uo - wi)
'
(2-BO)
48 Automatic Process Control
(2-81)
(2-82)
^
dU{) dui 2 \Uf) — U\l 2(wn -
defined as
dc c
K„ ^ = «/2A2V2j7(t/i
/
- U2 ) - — (2-86)
dm ?7i
^^ ^
dvj <9?/2
J2A2 m /
2 \//i
2ij
— ^/2/ 2(i/i
c
— U2 )
(2-87)
— U2 ) (2-89)
where Ka
Ki -f K2
KIK 2 c
A6 =
K, -h k2 2(?<o - WlO
direct; that is, any change ui the jiosition of the throttling valve m
results immediately in the same change in flow rate modified by the
process sensitivity Ka. Changes in upstream head (uo) and down-
stream head ( 1^ 2 ) act similarly.
Gas flow-rate processes may
be analyzed by the same procetlure as
that employed for liquid flow-rate processes. However, the fl«)w rate
depends upon three variables rather than two because the weight flow
of gas through a restriction depends upon upstream temperature, up-
stream pressure, and downstream pressure. Also, the flow through the
50 Automatic Process Control
at the throttling valve is greater than critical, the flow does not depend
upon downstream pressure U 2 and the factor /ir 2 is zero.
Thermal Processes i
The heat losses and the capacitance of metal parts in the heater are
assumed to be negligible. Rewriting equation 2-92, the process equa-
tion is obtained:
|m. Btu/sec
culated from
^
W weight of fluid contained
This time is the inverse of the rate at which the fluid substance is
changed.
The thermal process of Fig. 2-21 consists of
a heating furnace such
as might be used in continuous heat treating. A gas-air mixture is
burned in the furnace and represents a specified rate of heat addition, m.
We shall suppose that the main heating load is the metal boxes con-
52 Automatic Process Control
The stack losses are assumed to be a function of the heat input and the
furnace temperature.
9, =/(m,c) (2-95)
g. = Kitn + K 2Q + M (2-97)
(2-98)
PKL
where T =
PKV + K2
1 - Ai
K, =
PKV + UA -f K 2
UA
K,
PKV +UAA-K 2
The process thus turns out to be the single time-constant type. Notice
that the amount of material on the conveyor K alters the time constant
of the process. Control of thermal processes often involves variable
parameters in this manner.
^Heat-exchanger processes of the fluid-to-fluid type require a depth
of analysis beyond the scope of this text. These processes are usually
not describable by linear differential equations with constant coefficients,
although any given heat exchanger may be tested experimentally and
its dynamic characteristics determined. Sometimes, simple approxima-
tions of time constants may be made from experimental data.
Example 2-14. For the air-heating process of Fig. 2-20, the normal air
inlet flow is 300 cu per min at atmospheric pressure. The inlet temperature
ft
is 80 F. The volume of the heater is 500 cu ft. Calculate the time constant
of the process.
T
W ^ 500 X 0.072 X 60
100 sec
g 300 X 0.072
Dead Time
Dead time is defined as any definite delay between two related
actions. Dead time is illustrated in Fig. 2 -22 When the temperature-
measuring element is installed at a distance downstream of a heat
exchanger a pure time delay occurs liefore the heated fluid passes from
the heater outlet to the pmint of location of the temperature measuring
element. For example, if the fluid flows at 10 ft per see and tHe distance
is 10 ft, the dead time is one second.
54 Automatic Process Control
Dead time may be encountered in processes involving pressure. For
example if pressure is measured at the end of a long pipe or tube a dead
time occurs because of the time for transmission*^ of the first pressure
wave down the length of the pipe. For air, the wave velocity is roughly
1000 ft per sec and appreciable dead time would occur for a pipe more
than a few hundred feet long. For liquids, the wave velocity is much
faster, and the dead time is usually negligible. Dead time may also be
found in processes involving chemical reactions when a finite time
must elapse before a reaction begins to occur.
iU it/, temperature
i IV V IV. fluid flow rate
Rt - L) (2-99)
~
f(i-L)^m -f'(0L+r(i) +/''"(«)
^+ 2- 100 )
• •
•
(
where the primes denote the derivative with respect to time. Employ-
ing the operator .s*
fit
The power
- L) = fit)
series in the
[l - Ls + ^
above equation is
+ ^ + •]
Therefore
/(f-L) = c”^*/(0 (2-102)
G = (2-103)
Time, t —
Time, t
^
Fig. 2-23. Response of process elements with dead time, L.
PROBLEMS
q = 2gh^
fj — h^ •
Controller
Characteristics
the controlling means become. This does not at all mean tnat a compli-
cated automatic controller is necessary to produce good i.utomatic
control: on the contrary, the simplest control devices are often capable
of providing a high quality of control.
The control actions discussed in this chapter may operate through
either mechanical, pneumatic, hydraulic, or electric means. The various
mechanisms of automatic controllers are discussed in Chapter 6.
c == controlled variable
ft = feedback variable
« = set point
r == reference input
e= actuating signal
m — manipulated variable
59
60 Automatic Process Control
V « r
c = 6
e = v~c-r—
Deviation is the difference between controlled variable and set point
expressed iu units of the controlled variable c. The arrangement in
Proportional Control
Proportional action is a mode of controller action in which there is a
continuous linear relation between values of the deviation and manip-
ulated variable. Thus the action of the controlled variable is repeated
and amplified in the action of the final control element. For purposes
of flexibility, an adjustment of the control action is provide^ and is
the vessel inflow m. The action of the controller is such that a rise of
level c in the vessel rotates the valve operating lever clockwise thereby
closing the valve a proportionate amount. This reduces the inflow to
the vessel and tends to prevent the level from rising. If the distance
from float bearing to connecting link is large, a small change in level
produces a large change of flow into the vessel. If this distance is small,
ControUer Characteristics 63
a large change in level results in a small change of flow into the vessel.
Thus, proportional control follows the law
+M (3-2)
c = 0 ^ < 0
e = E (3-3)
m - M KcE (3-4)
Example 3-1.For the proportional controller of Fig. 3-3, the float arm
is 14 in. long,the link to bearing distance is 2 in. and the valve stroke is 1 in.
The valve passes 10 gal per minute when full open. What is the controller
sensitivity?
Integral Control
Integral action is a mode of control action in which the value of the
manipulated variable m is changed at a rate propctional to the devia-
tion. Thus if the deviation is doubled over a previous value, the final
(3-5)
(3-6)
m = ( 3- 7 )
3-4, the units of integral time are seconds squared per square foot. A
connection of the roller link in Fig. 3-4 closer to the float produces a
small integral time, and a connection closer to the float bearing produces
a large integral time
Ei
m-M^- * (M)
Thus, the manipulated variable changes linearly with time and **inte-
grates” the area under the deviation function.For a unit step change
of deviation (E ^ 1.0), the slope of the line is inverse of the integral time.
m« sin (at dt + M
and
^”
(aJ t
(at M
Since the manipulated variable m
is given by a negative cosine function, the
Proportional-Integral Control
Integral control action is often combined additively with proportional
control action. The combination is termed proportional-integral action
and is used for the purpose of obtaining certain advantages of both
control actions.
Proportional-integral control action is defined by the following
differential equation
m= ~«
I
+ KcC (3-10)
t
iQt^rel proportional
m= M
y jedt+ Kce+ (3-11)
inteipal proportional
Controller Characteristics 67
where m» manipulated variable
Kt =“ proportional sensitivity
Ti * integral time
e ^ deviation
t =* time
M = constant of integration
These equations illustrate the simple addition of proportional and
integral control actioifs. In operational form
m= Ke (3-12)
where the system function Kc/(l\s) identifies the integral action and the
system function Kc identifies the proportional action.
Proportional-integral control action has two adjustment parameters,
Kc and as indicated in equations 3-10, 3-11, or 3-12. The propor-
tional sensitivity is defined the same as for the proportional control
action: With the integral response turned off (T* —» « ), the proportional
sensitivity is the number of units change in manipulated variable m per
(3-14)
This is the equation for a straight line. The first term, i/Tt, is the in-
tegral response, and the second term is the proportional response. The
latter is indicated by the dot ted line of Fig. 3-5. Notice that the response
68 Automatic Process Control
Time, t
e
^
Fid. 3-5 rropcjrtional-reset control action.
~
K f
~ j
siti u)t dt Kf sin wf -h Af
Controller Characteristics 69
Integrating,
m- Af -1
K, 0)^
COB + sin o)/
^
Manipulating into a function of a double angle,
m—
Phase == tan“*
The phase lag depends inversely upon integral time 'A: for large integral time,
the phase lag is small; for small integral time, the phase lag is large.
Proportional-Derivative Control
A derivative control action may be added to proportional control
action and the combination termed a proportional-derivative control
action. Other terms for derivative response are rate response and lead
component. Derivative control action may be defined as a control
action in which the magnitude of the manipulated variable is propor-
tional to the rate of change of deviation.
A proportional-derivative control action is defined by
c = Et (3-17)
where E= sl constant
t = time
70 Automatic Process Control
+ Trff)
deviatione, —5 to 5 in.
m increment = E(0A)Td
*
= A(0.1)5 = 0 017 volt
e » sinciif
or
(3
-21 )
c = Et (3-22)
where E ~ a constant
t = time
m-M = + +
/ Td] (3-24)
Controller Characteristics 73
Two^Position Control
The two-position control action, or on-off control, is undoubtedly the
most widely used type of control for both industrial and domestic service.
It is the kind of control generally employed on home-heating systems
and domestic water heaters.
Two-position control is a position type of controller action in which
the manipulated variable is quickly changed to either a maximum or
minimum value depending upon whether the controlled variable is
greater or less than the set point. The minimum value of the manipu-
lated variable is usually zero (off) This mode of control is illustrated
by the electric-level control in Fig. 3 8. A float in the vessel operates
Set
1
point
1
L_ c
h J
m ^
—
1
L_ c •-
float lever has no bearing friction and the electrical contacts draw no
arc, the action is sharp or “knife-edge” as shown in Fig. 3-9.
The equations for two-position control are t
m = Ml when c > 0
(3-25)
m = Mo when e < 0
-AT 1
m t
valve stem. When the level rises, the switch contact is made,^nd the
motor-reducer slowly closes the control valve. As soon as the level
falls, the switch contact is broken and the motor-reducer reverses its di-
m = M\ when fi > +N
m = 0 when +N > e > — JV (3-26)
PROBLEMS
.VI . Prove tliat the sum of a sine wave and a cosine wave is another sine
wave of different amplitude and phase.
Use Euler’s relation c*"" = cos x -f isin x or employ vectors in the proof.
V2. A controller has the function m= (1 Tis -f- T 2 h^)e which in-
-f-
= - / fedtdt + M
Prove that the phase of the output lags 180 degrees behind the input.
3-4. A controller has the following functions in series; {TiS d- 1),
(T 2 S -f- 1), X/Ts. What parameters determine derivative time, propor-
tional sensitivity, and integral time?
For the problems below, determine what control actions are present and
state the parameters derivative time, proportional sensitivity, and integral
time if they exist.
V5. 3-6.
vhaptpr
Closed Loop
in Automatic Control
Controller Process
u I I
2. The controller lag is zero so that the control actions are those
described in the pre^nous chapter.
The analysis of the action of the controlled system requires the use
of the differential equations describing the action of each element of the
system: the process equation as developed in Chapter 2, the controller
equation, and the deviation equation as developed in Chapter 3.
The process equation in general operational form is, from Fig. 4-1
n
c = G2m + 5Z (4-1)
It is assumed that there are n load variables at the process, whereas only
one load is indicated in Fig. 4-1.
The controller equation in general operational form is
m= G\e (4-2)
e = V — c (4-3)
c = GiGje + z
• -1
(4-4)
The deviation is the same as the actuating signal e because the measuring
lag is zero (thus h = c), and because the set point v is the same as refer-
ence input r.
'Fhe closed-loop equations 4-5 and 4-6 illustrate that there are two
direct sources of deviation:
‘
These sources are regarded as the 'disturbances,” the effect of which the
automatic controller is expected to minimize
From the standpoint of industrial process control we are generally
interested in that part of the deviation due to changes in load variable.
For a constant set point, the deviation caused by any one of the load
variables u is
=
(t —I
-f GiGoJ
(4-7)
Cr = (4-8)
(1 + W,)'
The deviation resulting from a variation of set point is sometimes
employed to establish stability 'Of performance of the system.
Proportional Control
One of the simplest closed loops in automatic control is the propor-
tional control of the single-capacitance process shown in Fig. 4-2, On
Closed Loop in Automatic Control 81
The process equation is written under the assumptions that the up-
stream head at the control valve is constant:
-(7iTTr)"' + (wVi)'
where c = controlled variable = vessel head
R = resistance of outlet valve
T = RA = vessel time constant
A = vessel capacitance = area
m “ manipulated variable = inflow
u - load variable downstream head
82 Automatic Process Control
* " "
U/fc + + i)
m* Ktfi H- M
where M is the value of inflow when deviation e is zero. The constant M is suppressed
by writing
(m - Af ) * Kce
and measuring the value of inflow m from the zero-deviation value M so that
. m » /C«e
This procedure is followed in equations in this chapter by measuring values of
all
the variables c, e, m, and u from the ‘^normal” value existing when the deviation is
zero This procedure makes it unnecessary to carry all the constants through each
step of manipulation of the equations.
Closed Loop in Automatic Control 83
- U
B. - 14
(4 )
RKc + 1
Example 4-1. Calculate the offset and stabilization time for tbe level
control system shown in the figure.
I
m = {Kc)e {Kc){v - c)
+ /
= ^
'Phe offset is
A 40
Stiibilizini; time = 4 77- ~ ~ 20 sec
Kc 2
/?!
^ " (4-16)
(7’,s + 1) (7’vs + 1) T-.s + I
+ 1/ 7’V + 2f7’s + 1
( ) -J(Tgs + l) 4-18)
V/f,/f.+ 1/ rv + 2r7’« +
T1T2
where T characteristic time (4-19)
RiK, 4- I
rj + 2fT* + .. + (4-21
86 Automatic Process Control
TH + 2m + (4-22)
and
solution depends upon the characteristic time T and the damping ratio
i. The solutions may be classified according to damping.
2 3
Time, tlTi
Example 4-3. A process can be designed so that two vessels as in Fig. 4-4
have the characteri.stics:
(a) Ti = Tz = 00 sec, and = 1/8 sec per sq ft or (b) =* 30 sec,
T2 - 120 sec, and R - 1/8 sec per sq ft. Which design provides the least
offset? Wliich design provides the fastest response (frequency)?
“ U
- ^ ^ (ri 1 )
It per unit of
This example illustrates the rule that all lags in the system except one should
be small. The process above with the small tank and the large tank provides
less offset and faster response.
Integral Control
Integral control is often employed for control of fluid flow, liquid level,
and pressure. Consider integral control of the single-capacitance process
of Fig. 4-6. The process equation is
R
c VI + “ (4-26)
7V + 1 + 1
*
^ ^It
/ Ts+\
+ 2^T,k + 1
T^/_
R \t: +
1
r damping ratio
(ti - U\
-f 2iTri +€= 0 ( 4 -29 )
90 Automatic Process Control
Several solutions are plotted in Fig. 4-7. The offset is zero as pre-
dictedby the system equation, because the integral control action forces
a complete return to the set point when there is a change in load.
T, = (4^0)
With this value of integral time a minimum deviation is obtained without
encountering excessive oscillation.
Instability is po-ssible with integral control when the process does not
have self-regulation. The arrangement of Fig. 4-8 illustrates this. The
process equation is (if the outflow does not depend upon head)
c (4-31)
Cs ,
C«
(see Chapter 2)
c = Kam + (wo - W2 ) (4-33)
92 Automatic Process Control
Ka - ^ constant
Kb = Si constant
m = position of throttling valve
Equation 4-33 is valid only in the operating region of flow rate whei*e
the partial derivative factors involved in the /C's are nearly constant.
Is KiTsii\) KbTs
^ ~
Ts +1 Ts + I 7V+ ‘I
has not been demonstrated here, a few actual tests will show that a
small dead time makes a system with integral control uncontrollable
(steady or increasing oscillation).
Proportional-Integral Control
Proportional control and integral control are often combined in order
to obtain the advantages of inherent stability of proportional control
and of elimination of offset by integral control. Proportional-integral
control is the most generally useful of all types of control. We will
consider examples of this type of control, first, for processes without
storage elements, and second, for single-capacitance processes.
Proportional-integral control of fluid-flow rate is shown in P'ig. 4-9.
The process equation is,
where the quantities are defined as for equation 4-33. For proportional-
integral control,
w = A-, 1 + e (4-37)
(
where Kc = proportional sensitivity (inverse of proportional band)
T, = integral time (inverse of reset rate)
— — - i. \ '
_ (4-38)
KaKATs+iJ ^ k„kAtx+
uo)
KqKc d" 1
where T - Ti
KaKc
The response of the system to a step change in either load variable may
be found from equation 4-38. Foi a step change in downstream head
(« 2 ),
Te +e= 0 (4-39)
R 1
c “ '
m "T“
.
u (4-40)
Ts -h 1 Ts+l
where c = controlled variable = A = vessel head
H = resistance of outlet valve
T = vessel time constant= RC
C = vessel capacitance= A = vessel area
m = manipulated variable = inflow
u = load variable = downstream head
Time, iJT
Combining the process equation 4-40 with the control equation 4-37 to
eliminate w, and the deviation equation to eliminate c, there results
T,s
’
" ~ TMTs -b RK,(l\s + 1)
'
TMTs H- + RKr{T,s -b 1)
-f 1 ) I
)
(4-41)
, j{RKc + I)*r.- , .
Time, mm
Changing integral time
Fig 4-14. Adjustment of controller sensitivity and integral time for proportional-
integral control.
98 Automatic Process Control
The particular solution for given initial conditions is sho\\Ti in Fig. 4-13.
The damping ratio in I'ig. 4-13 is maintained constant at 0.707, a
1 1
V
CT,8^ + KciTiS + 1)
or
V Tc^8^ + 2f + 1
fCTi\^
where Tc
kkJ , . 4C J
For a damping ratio of one-third
4C
KcTi
9
Use of Derivative Control
Derivative control action is useful for control of systems having a
very large number of storage elements and for control of systems having
dead time. As an example, consider the control of head in the lower
vessel of Fig. 4-15. The process equation is
R _i_
(4-44)
i',s(r2s + 1) c,«“
where c controlled variable head in lower vessel
Cl capacitance = = area of lower vessel
Ti time constant = RCi
'
T2 upper vessel time constant = RC 2
C2 capacitance = A2 = area of upper vessel
R resistance of upper outlet valve
m manipulated variable = upper vessel inflow
u load variable = lower vessel outflow
Tis(T2S + 1) + KKciT^s + 1)
RiT2S +1)
TMT2S + 1 + ) + 1)
100 Automatic ProoeBS 'Control
(4-»7)
RK^ r.V + 2iTcS + 1 (k) + 2f Trf + 1
Time, f/Tg
* “
“ + \Tis+
(Tis + DCrj* + 1) 1)
Curve
Mode
Control
of
Prop
Sensitivity.
RK,
Denv
Time.
I’d.
sec
Period
of
Cycle.
sec
Damping
e-i/r
T,
sec
Max
Error
units
B
1 Proportional^denvative 16
- 0.9 32 15 018 0.06
4
Prop int
Proportional
denv
Proportional integral
10
4
22
-
3
a
H
2.0 44
45
66
20
20
30
021
0j»
0.37
0
0.11
v{t) =0 ^ < 0
(4-60)
v{t) = F / ^ 0
_ V
(4-51)
1 -f- G1G2
(4-52)
f -^0
Offset may result when the load on the system changes, thereby causing
a steady deviation. For a step change in load u,
u(t) =0 / < 0
(4-63)
u(i) ^ U < ^ 0
-NU
(4-64)
1 + G1G2
This equation represents the deviation as a f unction of time. For steady-
state behavior we again replace s by zero and the offset is given by
NU 1
(4-65)
1 -h G1G2]
sf; = A (4-57)
fit,
= (4-68)
E. - lim
*—*0
--^1
.5 1 -h G1G2J
(4-59)
* The proof of this relationship is best found by employing the final value theorem
of the La Flare transform.
106 Automatic Process Control
Eiample 4-6. Calculate the etatio error for the eyatem in Fig. 4-4
_V
E, lim
“
=
•-»o RiK^
"
^ ifc + l
(Ti* + 1)^
(7’i. + i)(r2« + i)
E, lim „ -0
1 4 ^
^ ^8(^8 + 1)J
Thus integral control eliminates offset.
Example 4-9. Determine the velocity error for the system in Fig. 4-6.
i + A.
^ T’.AT’i + l/J
Therefore this type of system may have a steady velocity error.
Example 4-10. Determine the velocity error for the system in Fig. 4-4
with proportional-derivative control:
jF, = lim - -
•-*0 s
,
RiKc(l + Tds)
\Ti8 + 1){T2S 1)
The limit is infinite, meaning that the deviation increases without limit and
that the controlled variable will gradually fall further and further behind
the set point.
Two-Position Control
Two-position control has the widest industrial and domestic use on
processes having notmore than two energy-storage elements. In opera-
tion, two-position control is very simple, but in theory the action is
difficult to analyze because of the discontinuous nature of changes in the
manipulated variable. Nevertheless, problems of two-position control
can be solved by quantitative consideration. In the discussion that
Closed Loop in Automatic Control 107
too low, and turns it off when the level is too high. For an increase in
For constant values of manipulated variable, set point, and load, the
above two equations may be integrated. For an increase in level
Equations 4-G2 and 4-03 describe the action of the system and are
As might be expected, the oscillation consists of a
plotted in Fig. 4--19.
A (4-64)
the magnitude of the load variable. From Fig. 4-19 it is noted that a
change of load variable may cause an unsymmetrical action, but what-
ever reduction of amplitude is gained on one side is lost on the other side.
The amplitude also depends upon the size of the inlet valve Af, the dead
Closed Loop in Automatic Control 100
time L and the process capacitance C. For small amplitude the process
capacitance should be large, the dead time should be small, and the inlet
valve size should be just larger than that necessary to handle the largest
load.
The period of cycle is calculated from the amplitude divided by the
rate of change of the controlled variable:
Example 4-11. The process of Fig. 4-18 has a steady outflow of 0.5 cu ft
per sec, a vessel area of 1.0 sq ft and a dead time of 0.5 sec. The differential
gap is 1 in. Calculate the amplitude and period of oscillation.
First calculate the valve size M. This should be about 50 per cent greater
than the maximum outflow;
= 14-
0.75 X 0.5 X 12
.4 3.25 in.
The period is
For constant values of manipulated variable, set point, and load, the
above two equations can be integrated. For an increase in level
where ti - -e = - V) — deviation
(c
These two equations describe the response of the controlled variable and
are plotted in Fig. 4-21. The oscillations are composed of segments of
exponential curves Without the operation of the controller, the con-
Closed Loop in Automatic Control 111
the dead time L. The period of oscillation depends directly upon the
process time constant and upon the amplitude of oscillation. The
amplitude is reduced by a small value of dead time L and a large value
of process time constant T.
Two-position control of processes having more than one capacitance
is similar in general respects to the control of single-capacitance pro-
ceases. The additioi^l capacitances of the system act to “round-off”
the peaks of the controlled variable oscillation.
variable moves out of the neutral zone, the throttling valve begins to
close at a constant rate and the controlled variable returns to the neutral
zone as shown in Fig. 4-23. (0)hapare this result to the integral control
in Fig. 4-10 for the same process.) Therefore, single-speed floating
control is stable for processes without capacitance or energy-storage
elements as long as dead time is small. With dead time, the overshoot
into the neutral zone requires increasing of the width of neutral zone
^
N and decreasing the floating rate so that the overshoot is less than the
width of the neutral zone.
0
Time —
Single -capacitance process with self-regulation
Ziegler-Nidhols Method
The calculation of the response of the controlled system is usually
difficult because^ of the mathematical complexity of solving the high-
order differential equations for each new control problem. In addition,
one solution of these equations is generally not sufficient because it is
below the set point and end just above the set point.
3. The step changes should be made in both directions if it is su^
loop test. If the controller has a lag, the controller should be set for
proportional control action only and the response recorded at m* in
Fig. 4-24 instead.
the apparent dead time and the apparent time constant T taken in
each direction are within about 10 per cent, a simple average may be
used. If they are not within 10 per cent, there is a serious nonlinear
behavior in the system.
The lag ratio should next he calculated;
NL = L
R - — (dimensionless) (4-70)
A I
M / /2\ units of m
(4-71)
AL \ 3/ units of e
P= (4-72)
NL ( ^ \ units of c
Eo (4-73)
M \3 -f 4/?/ units of rn disturbance
m ^ 1
(4-75)
c TiS
NL^ /I 4- unit of c
(4-76)
4M \ / units of m per min
*
“Theoretical Consideration of Retarded Control" by G. H. Cohen and G. A.
Coon, ASMS Trans., Vol. 76, p. 827, July, 1963.
•“Phase Plane Analysis" by D. P. Eckman, Mech. Engrg., Vol p 682, July
1953.
118 Automatic Process Control
(4-80)
(^ 1 )
units of c
E" "I (4-^2)
M VlO
<10 H- 9r)
9R/ units of
1 m disturbance
D - l.5E^
+ (4-»3)
\ units of m
(4-84)
*
NL \10 ^ 12,) units of c
min (4-85)
*
\9 + 20R/
Closed Lwp in Automatic Conti-oi 1 19
described by
(4-86)
(4-87)
:(7f)
" (4-88)
)
,11 + 2R/
32 + 6R\ (4-89)
,13 + Br)
Note: 6.
Ti
(4-90)
K.' = (4-91)
1 + 11
T,'
Ti' = K,Ti
Kc'
K,'Ti
T,' = (4-92)
K,
- =
e
—+
s
/C. + ( 4-«8 )
120 Automatic Process Control
(4-94)
(3 + 2R)M
" GATL*
9
A = (1 - «-") + (4-96)
2M
and
Note that the amplitude A and the differential gap G are expressed in
half-amplitude —that is, twice these values give the total cycle and total
gap. The equation for amplitude is adequate for any two-position
control, but the equation for period is valid only when the two-position
cycle is relatively symmetrical. When the lag ratio R is less than
about 0.2, equation 4-95 may be written
cent/degree
Example 4-13. For ii process with R == 0.5, what is the decrease in period
of oscillation if derivative control is employed with proportional control'^
From equations 4-72 and 4-81
^ (5 + 6ie)(10-f 7P)
^ 8_
Pp (6+ l0P)(12-h3/?) 11
Example 4-15. An
open-loop transient test shows a pressure system to
have an apparent dead time L of 0.04 minutes, a reaction rate of 2.5 psi N
per minute, and a process ultimate change K
of 10 psi for a change of valve
setting M of 5 per cent. If the set point is 100 psi, and the differential gap
is 0.5 [)si, calculate the amplitude and period of on-off control.
NL 2.5X0.04
.1 = 0.5
+ [— -2^^^ - 05
.
]
0.01 = l£j
P„
4Z,^
Z„
4L/1, 1 +Z \
^ tanh g (2 " 1 - z)
or
^ ~
R R {KM„„ - VM)
O.IC 1.5X5
0.01 10 X 100 - 100 X 6
0.24 min « 14.4 sec
122 Automatic Process Control
PROBLEMS
4-1. Explain physically why offset occurs, using the mechanical float-
used and if the remainder of the system contains only time-constant elements.
4-5. Sketch the resulting oscillations for on-off control when (a) differ-
ential gap G is zero; {b) dead time L is zero; (c) capacitance C is made small;
(d) control-element size M
only 10 per cent larger than load flow U,
is
4-6. For processes with small dead time, how much should the propor-
when derivative action is added?
tional sensitivity be increased
4-7. For processes with small lag, how much should the proportional
sensitivity be decreased when integral action is added?
4-9. In Fig. 4- 1, G, = X, (l + ,
G, = ,
iV = . Calculate the
+\
M2. InFig.4-l,G, = /Ce + -,Gj = ^=r^,JV = ;=r-^. Calculate
s Is -t I i « 1 -h
the relation necessary for a damping ratio of 0.707.
^
Calculate the proportional sensitivity for critical damping. Calculate the
static error. n o
4-14. l.P«.4-l,ft - lt,(r.. + l),ft . N. .
Calculate the proportional sensitivity for a damping ratio of one third. Calcu-
late the period of oscillation.
4-15. In Fig. 4-1, Gy * G, - — T7
n
A.; . . N
Tii
1
+ l
Measuring
(Feedback) Elements
The meap
will be reviewed from
principle of operation of measuring
the standpoint of automatic control; that no emphasis is placed upon
is,
e = r — b
Gi Controiiing element
Temperature Elements
Temperature ma>’ be nieasui’ed by mechanical, fluid or electrical
means.
Measuring (Feedback) Elements 125
a liquid that partially fills the system. Since the vapor prcssui c (k'pends
solely upon the temperature at the free surface of the licjuid, th(^ vapor-
actuated thermometer indicates only the temperature existing at the
free surface. Vapor thermometers may be used for tcunppratures as
low as -50 F and as high as 600 F.
126 Automatic Process Control
Platinum, nickel, and copper resistance bulbs are used, and some of
their characteristics are shown in Table 5-2.
Pressure Elements
The measurement of pressure or vacuum for purposes of control is
generally accomplished by applying the unknown pressure to a movable,
elastic member of constant area. The force thus created is balanced by
1. Bourdon tube, spiral, and helix elements are the most common
actuating elements for the pressure gage. These types of elements do
not require a spring since the elastic modulus of the metal is utilized to
128 Automatic Proceiss Control
tial pressure by applying the two unknown pressures to the two bellows.
The force of the bellows is balanced by a spring, and the resulting dis-
placement is proportional to the applied differential pressure. Absolute
pressure may be measui’ed by evacuating one bellows to a very low ab-
solute pressure, and applying the unknown absolute pressure to the
opposite bellows.
3. The hell pressure gage is useful for the measurement of low pies-
sures near atmospheric. Thesuspended on a beam in a sealing
bell is
liquid, and the unknown pressure i.s applied underneath the bell. "Jlic
pressure- tight shaft and is utilized for operating the controlling means.
5. In the ring-balance differential pressure gage the ring is balanced
on a rotating shaft or knife-edge bearings. The differential piessiire
displaces the sealing liquid in the ring and the rotation of the ring
results from the net torque due to the different pressures acting against
the clo.sed end.s of the ring.
(). The electrical pressure gage employs an unbonded or bonded strain
gage for indicating the pressure acting against the bellows area.
Liquid-Level Measurement
Measurement of li(iuid level is accomplished by two different methods.
Head measurement by a float following the liquid-level surface i)rovide.s
.‘iccurate measurement of liquid volume independent of liquid density
1. The float and tape method is employed for open \eshels. Other
mechanisms employ a free float on a pi\'oted arm to operate an indicat-
ing means.
2. The diaphragm-box system is used for pressure measurement in
open vessels. The metal box contains a flexible diaphragm which acts
Measuring (Feedback) Elements 129
to compress the gas in the pressure receiver line when the liquid level
increases. The pressure receiver is then calibrated in terms of liquid
head. The bubbler system and the air- trap system are similar means of
liquid-level inoasiireinent.
3. The (lisplacrmenl-floal metfr may be used in either open or closed
pressure vessels The buoyant force caused by the displaced liquid is
but two factors of design serve to increase the differential head for a
given flow rate. First, the venturi section increases velocity so that
velocity head is and second, the pitot has both forward and
increased,
rearward openings to increase differential head.
The fluid passes up the tube through the annular space around the float.
buoyancy, and drag acting upward against the weight acting down-
ward. A scale etched on the glass metering tube is calibrated in flow
rate and is read by noting the position of the float against the scale.
The electromagnetic flowmeter element is an clement operating
electric
on the principle that the liquid flowing through a nonmagnetic tube in
a magnetic field induces a voltage that is proportional to flow rate and
field intensity. The field is supplied by an electromagnet. Electrodes
sealed in the tube detect a direct-current or alternating-current- voltage.
The flowing fluid should have some electric conductivity. The flow
rate is measured independently of viscosity, velocity disturbance, etc.
The by passing the fluid out one
Coriolis flowmeter element, operates
radial tube and back another radial tube. The force difference between
the two paths is caused by the Coriolis acceleration. 'J'he force is
measured by an elastic member, strain gages, and slip rings. The flow
rate is measured independently of viscosity, velocity disturbance, etc.
Pneumatic Transmission
It is often required to measure the value of a controlled variable
where the location of the point of measurement is at a considerable
distance from the location of the controller. Most measuring devices
such as a mercury thermometer, pressure gage, or flow-rate meter would
then require fluid-line connections of great length. This cannot be done
Measuring (Feedback) Elements 131
c, measured
pressure,
Receiver Transmitter
IS distnbuted along the length of the tube The resistance of the tube
IS due to fluid friction aiul is likewise distributed along the tube. Because
of these distributed paraini'lets, the calculation of the lag is extremely
difficult and it is necessary to te\y upon experimental tests to determine
the lag.^ The results of such a test, shown in Fig. 5-4, are given in
^ J C. Moise, A and Analytical InvesUyatum of Long Pneumatic TranS'
Theoretical
mission Iaucs, Doctoral thesis. Case Institute of Technology, Cleveland. Ohio, 1962.
Measuring (Feedback) Elements 133
Lag compensated
Temperature transmitter
The amplifying pilot provides a gain of one so that the output (trans-
mitted) pressure is the same as the input pressure. The pilot is simply
inserted in the transmission line. Its advantagi; derives from the large
ports of the pilot valve so that a large flow of air is allowed when a
change of pressure takes place. The operation is as follows: An increase
in input pressure acts against the diaphragm and pushes the pilot stem
to the left. This opens the supply port and allows air to flow into the
output volume. The increased output pressure acts to the right against
the second diaphragm and exactly balances the input pressure. On a
decrease in input pressure, air is exhausted from bet ween the diaphragms.
The lag-compensated transmitter provides a boosting action that
tends to overcome transmission lag. The operation is as follows: An
increase in temperature at the thermometer bulb increases the pressure
in the right inner bellows. This move.s the bellows assembly to the left,
134 Automatic Process Control
covers the exhaust nozzle, and increases the output pressure to a large
magnitude. Then, as air flows through the restriction to the outer
bellows chamber, the bellows moves gradually to the right and reduces
the output pressure to its new but larger value. The action is similar
to a derivative effect and is described more precisely by
^ “
"b 1
(5-1)
c Tlj. + 1
If the upper restriction were not present in Fig. 5-5, the compensation
Td would be zero. Thus, the right-hand (derivative) term of equa-
tion 5-1 indicates a pure derivative effect in the presence of a time-
constant lag Tl. The net effect is that of reducing the time-constant
lag of the system. The dead-time lag is also partially compensated.
Electric Transmission
Open loop
0 to 0.5 V, a-c
c, displacement
Closed-loop transmitter
(Manning, Maxwell, & Moore)
arm. Thus the output direct current of 0.5 to 5.0 milliamperes is pro-
portional to input force. The advantages of this system are that the
direct-current transmission signal is not easily interfered with and tlie
signal may be manipulated easily, giving great flexibility.
b, displacement
c, displacement
the same line source so that the primaries are in phase. The closed-
transmission-loop systems have the advantage of flexibility.
The impulse-duration transmission system is shown in Fig. 5-9. The
transmitter consists of a synchronous motor driving an electric contact
through a cam in such a way that the duration of make of the contact is
proportional to the variable being measured. The receiver is a clutch
b, displacement
e. pressure
Firsl-Order Response
Speed primary measuring element and the measur-
of response of the
ing means of a controller the most important single factor affecting
is
arm. Thus the output direct current of 0.5 to 5.0 milliamperes is pro-
portional to input force. The advantages of this system are that the
direct-current transmission signal is not easily interfered with and tlie
lb. displacement
c, displacement
Synchronous
motor
Fig. 5--9. InipulHe-dumtion tniusmiHsion.
and the signal requires only two wires. The speed is generally slow,
with a cycle of from 5 to GO seconds.
There are many other forms of electric transmission such as synchro
systems, and potentiometer systems. These are too numerous to
review here.
First-Order Response
Speed primary measuring element and the measur-
of response of the
ing means most important single factor affecting
of a controller is the
the operation of an automatic controller. Since automatic control is
a continuous, d 3 'iiainic function,* the rate of detection and the time
element in response of the measuring means form an essential part of
automatic control analysis.
If we were suddenly to immerse the bare bulb of an expansion ther-
mometer in an agitated salt bath maintamed at a constant /temperature,
the thermometer pen would rise as indicated by the curve of Fig. 5-10.
138 Automatic Process Control
12 seconds to reach 99 per cent, and another 12 seconds to reach 99.9 per
cent of the final value. The difiBculty of expressing lag in terms of time
for indicating a complete change is that theoretically the final
temper-
ature will never be reached.
Measuring (Feedback) Elements 139
The response curve for a ^ple thermal element is derived from the
laws of conservation of energy and heat transfer. As indicated in
Fig. 1^10, the increase in energy is equal to the inflow of energy, or
(WP)h - q (5-2)
q = HA{c-h) (5-3)
Tb +b= c (5-4)
= 1 - (5-6)
c - bo/
6-60
= (1 - 0.368) = 0.632
c — bo
required to reach 63.2 per cent of the final value. Note also that the
initial slope intersects the final value at a time T.
The general approach to the response of a measuring means is indi-
0 12 Time, tfT
3 4
K
(^7 )
1
Thermal-Element Lag
Temperature-measuring elements such as the expansion thermometer,
thermocouple, and resistance thermometer possess appreciable measur-
and
ing lag because of the teat capacity of the materials of construction
the heat transfer characteristics ot the surrounding media. Heat is
transferred to the thermal element by conduction, convection, and
radiation but the conduction heat transfer is usually small compared
to the convection and radiation transfer of heat.
The bare thermocouple and bare thermometer bulb installed in
rapidly moving fluid have almost first-order response as may be seen
by the actual test curves of Fig. 5-12. The element time constant can
be calculated if all the material and fluid characteristics were known.
From equation 5-4, the time constant is
WP
T
HA
Therefore a small time constant requires small weight but large surface
area, low specific heat, and high heat transfer coefficient.
The heat transfer coefficient can be calculated for forced convection
transfer of heat from a fluid flowing at right angles to a cylindrical
thermal element from
Nu = n(Nj,f\Nrr
where Nu = Nusselt number = HD/
H— convection coefficient, Btu/ft^ F~^ sec~^
D= diameter of tube or rod, ft
K - thermal conductivity of fluid, Btu/ft F~^ sec
m = a constant given below
n = a constant given below
Np = Prandtl number = gCn/K
C = specific heat of fluid, Btu/lb deg""*
=
fi absolute viscosity, lb sec/ft^ = vy/g
v = kinematic viscosity, ftVsec
7 *= fluid density Ib/ft^
Nr = Reynolds number * VD/v
V - velocity of fluid, ft/sec
142 Automatic IVooeas Control
The oonstante of the equation depoid upon R^olda numbw as follows
Nr n m
4(M000 0.68 0.47
4000-40,000 0.19 0.62
40,000-400,000 0.26 0,80
Time, sec
Fic. 5-12. Response of thermal elements. (Not necessarily typical of each type.)
inversely as the 0.6 power of velocity. The data of Fig. 5-13 show
reasonable agreement.
The characteristics of the fluid surrounding the element are also
Automatic Process Control
important. Equation 5 -9 indicates
so that high thermal conductivity, high specific heat, high density, and
low viscosity are necessary for fast response.
V o.'i X 10"' X 12
„ K .. 0.38 X 122 X 12 , „
3600 X 0.125
WP ^
ir P^Jjym P -- = 0.073
A iirDL
^ WP 0.073
^ M = IT = ^
q = HrA{c-b) (5-H)
if the coefficient of radiation transfer Hr is defined by
„ 1 dq 4KEfi
Example 5-3. A thermal well made of standard J-inch iron pipe is used
furnace (reducing atmosphere) at 1400 F. (O.D.
in a = 0.84 in., 0.85 lb per ft).
1 70 X 0.1 1
r
X 12 ^ 42
0.0 1 X X X 0.84 X 2
Prcssurc-Eli'iucul
Displacement element
measurlnp, liquTd pressure
pA - Kb (S-14)
The flow of liquid w into the bellows through the long connecting line
depends upon the pressure difference and resistance of the line and is
given by
W= (C - p) (5-15)
I
where R= resistance of line, sec/ft*
c = pressure to be measured, Ib/ft*
Tb ^c
+b=K (5-16)
T-K
(y )
The capacitance of the system is given by the parameters enclosed by
the parenthesis. Naturally small capacitance results in small element
time constant and requires a small bellows area and large spring gradient.
Therefore, small displacement leads to small element time constant.
The resistance R of the connecting line should also be small. Line
resistance is discussed later in this section.
The diaphragm pressure gage of Fig. 5-14 is considered to be a null
system which a balancing force 6, supplied by some external mech-
in
anism, maintains the diaphragm in a fixed position. Then for a fixed*
Measuring (Feedback) Elements 147
volume pressure chamber the rate of change of pressure ie given by
Cf^w (5-18)
Therefore small volume and small line rcosistance are conducive to fast
response.
The resistance of a long connecting line having bends, fittings, and
valves depends upon whether laminar or turbulent flow exists. Actu-
ally under transient conditions the flow increases from zero to some
high velocity, usually turbulent, then decreases to zero if there are no
leaks in the line. Thus the calculation of line resistance is made very
difficult and for approximation purposes it is usually assumed that
laminar flow exists. Acjcording to the Poiseuille Jlagen law, the
laminar resistance is
dp __
128 yL _ 128 /iL
(5-22)
dw irD^y
(Raa) (5-24)
tD* \npj
where To volume of element at operating pressure, ft^
n polytropic constant (n = 1 = isothermal)
Pa average absolute pressure.
A element area, ft^
Yl 1 X 0.123
=
C =• 2.0 X 10“» ft^
nPa Tra XIX 24.7 X 144
R = 23,000 8ec7'ft2
For air.
Second-Order Response
The response of many temperature and pressure elements cannot be
adequately described by a first-order differential equation. Pressure
elements, for example, may have inertia effects in addition to fluid
effects. Temperature elements often possess more than one energy
storage or capacitance. Therefore, a second-order differential equation
must be used to describe adequately many pressure- and temperature-
measuring elements. ^
The response manometers or U-tubes is sf^cond order because
of fluid
inertia of themanometer fluid is involved. Applying Newton’s second
law to the manometer of Fig. 5-15,
or
© b+ (RA)b-¥’b= -c
7
(5-26)
_RA jg
~ damping ratio =
L
2 \ (
and when the damping ratio is greater than one (f > 1 ), the ^jratem
18 overdamped.
0.6 to 0.8 (see Fig. 5~15) so that fast response is obtained without
excessive oscillation.
Many dynamic systems in engineering are described by a second-
order linear differential equation with constant coefficients. Such equa-
tions can always be expressed in the following form where the coefficient
of the zero order term is one:
Ab + Bh-\-b^ Kc (1^28)
c
“ rV + 2rr« +
where T= VZ = characteristic time
f == = damping ratio
When the system is overdamped, equation 6-28 may be written
b
^
K (5-80)
C iTi8+lKT2S+l)
Example 5-6. For the bellows gage of Fig 5-J4 with a mass of the M
spring and bellows, derive an expression for the characteristic time and damp-
ing ratio.
When measuring liquid pressure, Newton^s second law yields
MS pA - Kb
Also, the motion of the bellows depends upon the flow rate in the line:
(Ay)S = w = ~(c p)
Combining,
152 Automatic Process Control
Therefore T •
ly/UK
The natural frequency ifi
^ ^ /I.
2tT 2fr\M
Example 5^7. A pressure of 15 psi maximum is to be measured with a
mercury U-tube. What is tlic smallest characteristic time?
m /l5 X 144
C262 = ^2 — (5-31)
9* 1
t\2
(c-w (5-82)
^*
1^1 r» (^2 “ ^>l) (M3)
Ki
Measuring (Feedback) Elements 153
h 1
(5-34)
(R2C2R^C\)s^ + + R2Ci + /?2C2).^ -f I
c (Tis+ 1)(7V-|- 1)
where
1
„ R\(^\ ~t~ R 2^ 1 ~i~ ^ 2 (^2 r. 4R2^2Ri^I
^
L —
Fx'tI" m,'c,R2C7'~
[ \ {RxC\ 4 R 2 CX H- R2C2y\
A close inspection of the damping ratio will show that it cannot be less
Example 5-8, A thermal well has a time constant of 60 sec when installed
without a thermocoujile in it. The thermocouple shows a time constant of
3 3ec when suddenly inserted into the well. The weight of the thermocouple
!
junction is 1 oz and the weight of the thermal well is 0.5 lb. Assume the
specific heats to be the same. Calculate the time constants.
RiC} 4 RiCi
2
4RC =2 2 R\Ci
“T +
,
C'l E>
+
I
~
R\Ci 4RC 4RC _3
2 1 2 2 A CO j- —
'^2
2 “2 2X0.6
RiCi 4 R 2C 4 R 2 C 2
1
35.2
2
^R2C2RiCi 60 XJ
(R^Ci 4 ^22^1 4 i22Cj)* (35.2)2
154 Automatic Process Control
PROBLEMS
Controlling Elements
particularly used.
The automatic controller consists of the measuring element H of
Fig. 6-1, the actuating signal element (the circle), the control element
G. O,
Gi Controlling element
Self-Operated Controllers
Many automatic controllers operate by employing power developed
by the measuring means. These controllers are usually simple and in-
expensive, and are widely used in both industrial and domestic service.
The pressure regulator or pressure-reducing valve of Fig. 6~2 may be
either a “floAv" type or a “dead-end” type. The only difference between
these two types is that the “dead-end” type incorporates a relief valve,
indicated by the opening through the diaphragm in Fig. 6-2. The set
point determined by the adjustment of the spring compression. The
is
the nozzle opening in accordance with the actuating signal e so that the
nozzle back pressure mi is inversely proportional to the distance be-
tween nozzle opening and flapper. A
motion of the flapper of about
0.0015 in. is sufficient for nearly full This device may pro-
output range.
duce proportional control if the flapper is accurately positioned by a
Controlling Elements 159
mi = + Ml (6-1)
shown in Fig. 6-8. The system is similar to the mechanism of Fig. 6-7
except that the output nii is connected to a feedback bellows which
operates to repo.sition the flapper. The nozzle restriction is usually in-
corporated into the bcxly of the pilot and is omitted from the figure.
signal changes linearly with time. The flapper then moves slowly to
the right, covers the nozzle, and causes a slow increase in output pressure
nil. There is a flow of air through the feedback restriction resulting
from the fact that the pressure in the feedback bellows is considerably
less than it would have been without the restriction in the line. This
delays and reduces the feedback, and, since the feedback is negative, the
output pressure is higher and leads instead of lags the actuating signal e.
Thus the delayed negative feedback produces derivative response. The
operational equation •describing the proportional-derivative action is
ideally
(6-3)
nil
This arrangement thus provides the proper control action as long as the
two resistances an* not set at the same value Such controllers are
arranged so that the derivative time T 2 cannot be set to a value equal
to or greater than the integral time
If the restriction is added in the feedback line at point 1, Fig. 0-10,
the response* may be represented by
(o-r>)
c
This arrangement provides the deared control action and avoids any
Loiitations of setting derivative time.
per cent). On the other hand, if restriction Ri were wide open, restric-
tion R 2 to atmosphere being fixed, the pressure in chamber / would
always be the same as that in chamber mi. In the latter case there is
no negative feedback and the proportional sensitivity is very high. For
settings of restrictionRi between full open and closed, the pressure in
chamber / is governed by the flow through the two orifices in series (f^i
and 7^ 2 ) so that positive feedback cancels a varying portion of the
negative feedback. Therefore, the restriction Ri adjusts the proportional
sensitivity of the controller.
Integral response can be added to a proportional controller by the
addition of delayed positive feedback as shown in Fig. 6-13. The pro-
[i + + r..]
1 (6^)
Where the symbols are the same as above. This equation is easily de-
rived by combining equations 6-2 and 6-3 to eliminate the intermediate
\^ariable. This type of eontroller has certain advantuges with respect to
start-up of batch processes.
Fia. 6-15. Passive- type pneumatic control (direct and inverse derivative of the
Moore Products Co.)
When the gain ratio or is small (it is usually set at about one tenth
the effect of the control action is somewhat like a proportional-integral
action in which the time constant T adjusts the integral-like time con-
stant.
With the device on the right in Fig. 6-15, the first and second dia-
phragms are reversed, w'ith the result that the gain ratio a is also inverted
and a modified proportional-derivative action is obtained. The action
is described by etiuation (fi-S), but with the gain ratio a equal to approxi-
150-200
*
UrproducecI by permission from Ptincipla vj Indnulttal Process CorUrol, John
Wiley and Sons, Inc., New York, 1049
108 Automatic Process C.'oiitrol
t'KJ. 6-17. Air-bupply .system for small installation of pneumatic control equipment
remaining moisture, oil, or solids just before the air enters the con-
troller. The waste-type filter is most common and is quite satisfactory
under ordinaiy conditions. An individual small-capacity pressure
regulator should accompany each controller or group of controllers in
order to provide proper pressure regulation.
Moisture-free compressed air does not freeze even at subzero temper-
atures. Therefore, methods of effectively removing moisture also aid in
avoiding trouble due to freezing of air lines. Antifreeze injection sys-
tems arc not adequate, since the total amount of entrained fluid in the
Hydraulic Controllers
Hydraulically-operated controllers for industrial process control pro-
vide great power and positiveness of action since they employ a source
of high-pressure oil as auxiliary power. In certain instances high-
pressure water is used as a power source
170 Automatic Process Control
ferential link connecting the diaphragm and spring. Oil at high pressure
is supplied to the pilot valve. The operation is as follows: When the
actuating signal increases, the pilot valve stem is moved to the left
This supplies oil to the right side of the power (‘ylinder and exhausts oil
from tb(‘ left side. Consequently the power piston moves to the left.
Its rate of motion is proportional to the actuating signal. The result-
ing integral action is described by
m, =-J-e (6-9)
7 t’S
w here nii
~ position of power piston
T, ^ integral time
e actuating signal
return from the system and serves to dissipate heat. This latter
oil
Compressed
inch of gas, usually nitrogen. With suitable check valves, the accumu-
lator may run the control system for several minutes in the event of
electrical failure. The main sources of difficulty in hydraulic systems
are dirt and Extreme care must be taken to insure that the sys-
heat.
tem is initially clean. After that, particles of metal and oil deposits
are bound to occur, and the system must be periodically cleaned. Heat
causes great variations in viscosity of the oil, and cooling must be
arranged. Usually a large sump may dissipate sufficient heat by
radiation.
Feedback link
(\
— DoooU
Tx
'
0 ooo
y q
From pump f | To sump Power cylinder-^
pilot valve, and balances the power piston at its new position. Thus,
the position of the power piston mi is proportional to the actuating
signal e,
Ttii = KcC ( 6— 10 )
Propo.iional
fWvtwVi
sensibvity — 1|“
m, = KrP (6-12)
Fic; (’>-2.'^ PY'iicJliack circuits for the electric controller of Fig. 6-22.
(he (*l('ctnciil system to the pneumatic vsystem of Figs. (Ml and 6-14.)
'flu* derivative time is proportional to the electrical time constant RC
(contactor arm to the left in Fig. 6-22, away from the set point, a d-c
voltage is passed to the amplifier and drives the power motor in one
direction. The feedback voltage provided by the power-motor siidewire
Controlling Elements 175
causes a current to flow into the capacitor C, charging the capacitor and
resulting in a voltage drop across the resistor R, This is proportional
action. Now, however, in order to maintain a steady current through
must be continually charged. This can only
resistor R, the capacitor
be done by a slow and continuous motion of the power-motor arm to
the left. This is integral action. The integral time is proportional to
the electrical time constant RC and is selected by adjusting resistor R,
Proportional-integral derivative action is provided by inserting a suit-
able combination of the feedback circuits of Fig. 6-23 in the circuit
of Fig. 6-22.
The fl-c bridge proportional controller is shown in Fig. 6-24. As in the
previous electric controller the actuating-signal slide wire (at the left)
Integral amplifier
positioned by the output arm of the power motor. The feedback bridge
provides integral action as will be explained later. The operation is as
follows: When an increase in actuating signal moves the slidewire con-
tactor to the left, a-c voltage and resulting current is carried by r(‘-
Neglecting the operation of this bridge, the net efTect of the feedback
voltage is to cause a current through resistor Rh which must balance the
current through resistor Ra. This is proportional action and the power-
motor arm must ‘‘follow’* the actuating signal contactor. Integral
action is as follows: Whenever a current through resistor Ra,
there is
the balance arm moves to the right and causes the contact to close a
circuit to the power-motor. The output arm of the power-motor moves
counter clockwise and moves the piston of the feedback cylinder to
the left. This action compresses the air in the feedback system and the
feedback diaphragm provides a feedback force balancing the actuating
signal force.This is proportional action. The pressure in the feedback
system is now greater than atmospheric pressure and air flows out the
Controlling Elements 177
element with the position of the output-arm fed back in a proper manner
to provide the desired control action. The advantages of this system
arc that the actuation is electrical and therefore distance between
elements is theoretically unlimited. Another advantage is that the
actuator is enclosed within the controller loop.
The first d-c amplifier has a ‘^delayed*’ negative feedback and provides
proportional-derivative action as has been described previously for other
electric controllers as well as for pneumatic controllers. Proportional-
integral action is provided by the second d-c amplifier with ''advanced^'
negative feedback as has been described for the electric controller of
Fig. 6- 23 . Afljustable proportional sensitivity is provided by an adjust-
able gain of the a-c amplifier. Adjustable derivative time is provided
by the setting of the feedbacik resistor of the first d-c amplifier. Adjust-
able integral time is provided by the setting of thf? feedback resistor of
the second d-c amplifier. The output of the controller is a d-c current
which is carried through wires to a final control element. The com-
bination of proportional-derivative action and proportional-integral
action in cascade has an advantage with respect to start-up of control
systems. For the proportional action of the a-c amplifier and rectifier,
(6
- 13 )
^ = K2(1 + T^)
W3
(fi- 14 )
nil
(^ 15 )
m2
where e
— actuating signal = i; — c
K, = gain of a-c amplifier and rectifier
K2 =- gain of first d-c amplifier = (R2 + Rz)/Ri
Ki = gain of second d-c amplifier = R2/R\
mi = output current (or voltage) of second d-c amplifier
m2 = output voltage of first d-c amplifier
m3 --= output voltage of a-c amplifier
derivative time = R2RzC/{R2 R3)
- integral time = R2R3C/(R2 + R3)
also affect the overall proportional sensitivity. The latter is also affected
Controlling Elements 179
by the gain of each amplifier. However the gains of the d-c amplifiers
are generally fixed and the gain of the a-c amplifier made adjustable.
The **voice-coiV* motor shown in Fig. 6-27 is the primary transducing
device employed with electronic controllers in order to convert direct
current into a mechanical force. The operation is based upon the force
caused in a conductor carrying a current through a steady magnetic
field. The force is linearly proportional to the product of current and
field strength.
set point v and controlled variable c and is the actuating signal voltage.
This voltage is acted upon by the proportional-derivative amplifier and
easily shown to be
(6-17)
6- 18 )
(
e — ei = iR\ (6-19)
The current flow into the amplifier is negligible (grid current only) so
all current ^ passes through the feedback resistor 7^2-
ei — mi = iR2 (6-20)
The amplifier has a fixed gain K and in addition changes the sign of the
voltage:
wij = Kci (6—21)
mi Ro 1
~ ~
T Tf,
KRi
The amplifier gain K is usually large (10^ to 10®) so that the above
equation reduces to
- = - R2
Till
^ (6-23)
€ Hi
The negative sign shows that the output voltage (reference to ground)
is negative for a positive input voltage. As a proportional controller
the proportional sensitivity is adjusted by selecting the ratio of re-
sistances (R 2 /Ri)-
Proportional-derivatwe controllers one of two basic
employ either
principles: (1) delayed negative feedback, or (2) advanced input, as
shown in Fig. 6-31. The delayed negative feedback employs a feedback
capacitance to ground. Thus the output voltage mi must continually
change in order to charge the capacitor, 3delding a derivative effect.
The operational equation may be derived by the same procedure as
Controlling Elements 183
before. Since the amplifier gain is very large (10^ to 10®), the opera-
tional equation for the left circuit of Fig. 6-31 is
-ji:e(r<,s+ 1) (&-24)
where Kc = 2R/Ri
Td = RC/2
The derivative time,T<f is adjusted by setting either or both resistors;
the proportional sensitivity is adjusted by setting the ratio of resistors.
Ri R2
r-V\Ar-r-W\rn
i>
Advanced input
—= -KciTiS + 1) (6-26)
e
where Ke = R 2 /R 1
Td = RiC
7 +
where Kc = R4/R2
= RC/2
The integral time is set by adjusting either one or both resistors R.
The proportional sensitivity is set by ratios of feedback resistors. This
and 6-13.
The advanced negative feedback method of Fig. 6-33 is commonly
employed in electric controllers. By inspection of Fig. 6-33 it is seen
VW 11“
that unless the input voltage is zero there must exist a rate of change
mi
(6-27)
e
where Ke R/R\
Ti RC
Controlling Elements 185
where Kc = R^/2Ri*
a = 2/^2/ (2/^2 “h ^1)
T = (2/?2 -f Ri)C/2
— = — + K, + Kis (6-29)
€ S
tiometer, for example, and determines the position of the upper contact
' 7^
j©
r Furnace heater
Synchronous motor
Floating
motor
ij--
Furnace
Synchronous heater
motor
In a narrow region near the set point, neither contact is made. This
latter region is termed the neutral zone.
The HI-LOW contacts supply power to drive a slow-speed motor
in one direction or another, and to increase “time on“ when the LOW
contact is made and to decrease time on when the HI contact is made.
The floating motor usually requires several minutes to change from zero
“time on“ to full “time on.“ Single-speed floating control action re-
sults. In equation form,
1) m
when c> p H- Ar/2
The floating rate K is adjusted by selecting a gear ratio for the floating
motor. The neutral zone N is adjusted by setting the position of the
HI-LOW contacts. The period of cycle P is adjusted by selecting a gear
ratio for the synchronous motor.
A ’proportional-integral controller for electric heating is shown in
Fig. 6-36. This arrangement provides a control action in which the
slowly than they charge, the load relay is on for less than one-half
time. The fraction of time on mi is therefore proportional to the actu-
ating signal because the rates of (charge and discharge are determined
by the net voltages from the input and feedback.
PROBLEMS
&^2. Show that the pneumatic amplifier has proportional actions described
by pressure m - Ke.
6--3. Show that the pneumatic amplifier has proportional actions described
by m = Ke.
192 Automatic Process Control
6-4. The flow through the open nozzle when the velocity at the restriction
-ir
6-5. Assuming a linear relation for the nozzle, derive the output-input
relations for the hydraulic controller {A
- piston area, K — spring constant,
K 2 *= dq/dp nozzle,
= r)q/dp supply, K\ = dg/de nozzle).
6-6. For the gas-pressure regulator calculate the offset when acoustic
velocity exists at the throat.
6-7. CalculaU> the input-output relationship, and state the kind of control
action.
cl
Controlling Elements 193
M. Calculate the input-output relationship, and state the kind of oontrol
action.
6-9. Calculate the input-output relationship, and state the kind of control
*
action,
proportional-integral action; c = A+ r.
Final Control Flements
The final control element is the mechanism which alters the value of
the manipulated variable in response to the output signal from the
automatic control device. The position of the final control element in
the automatic control loop is shown in Fig. 7-1. The final control
Gj La*!" C,"*
1 I 1 I 1
^ 1
I
'
b
L-^-l
!
H '
t*
!
1
I I
Gi Controlling element
Pneumatic Actuators
Pneumatic actuators may operate directly from the pneumatic output
signal from a pneumatic controller, or they may employ a separate
source of compressed air. There are five common methods of operation
of pneumatic actuators. These are called
1. Spring actuator
2. Spring actuator with positioner
3. Springless actuator
4. Piston actuator
5. Motor actuator.
A spring achmior shown in Fig. 7-2 operates directly from the air-
pressure output of a pneumatic controller in order to provide an output
position proportional to the input air pressure. The diaphragm is
^ VK/M (7-2)
negligible.
Static friction forces must be limited to a low enough value that
excessive hysteresis does not result. For hysteresis less than one per
cent of full travel,
F/ ( 7^)
where Ft is the total thrust force acting in one direction. This may also
be a serious limitation. For example, a spring actuator with an effective
area of 100 sq in. and an initial air pressure setting of 3 psi cannot sup-
port more than 300 lb unbalanced force. In addition, thrust forces
198 Automatic Process Control
must be relatively constant otherwise the stem position will not be
directly related to the input air pressure. The performance of a spring
actuator is also influenced by the characteristics of the spring and
diaphragm. A well-designed actuator has a linear static relation be-
tween input air pressure and output stroke if the effective area of the
diaphragm and the spring gradient are constant throughout the stroke.
Hysteresis due to the stresses in the spring and diaphragm are usually
less than one or two per cent of full stroke. .
i«tein assumes a position dictated by the input air pressure. The spring
actuator becomes a power means and the characteristics of the spring
and diaphragm are relatively less important. The use of the positioner
results in several improvements in performance.
The use of a positioner with a spring actuator does not improve the
ability of the actuator to handle larger inertia or thrust forces unless
special adjustments of motor operating range are made. The only dis-
advantage in the use of a positioner is that it may require maintenance.
The 8pringle88 actuator shown in Fig. 7-4 is useful for large thrust
forces. The only difference between the springless actuator and the
spring actuator with positioner of Fig. 7-3, is that the spring of the
actuator is replaced by a pressure regulator which maintains a constant
pressure on the under side of the diaphragm. An air supply at a pres-
sure of 20 to 100 psig is required. The operation of the springless
actuator is as follows: Assume that the cushion regulator is set to pro-
vide 9 psig pressure on the under side of the diaphragm. At static
balance and with no thrust force on the actuator stem, the upper side
pressure must be 9 psig. Then if the input pressure increases, the
nozzle back pressure increases, and the upper side pressure is raised to
198 Automatic Process Control
must be relatively constant otherwise the stem position will not be
directly related to the input air pressure. The performance of a spring
actuator is also influenced by the characteristics of the spring and
diaphragm. A well-designed actuator has a linear static relation be-
tween input air pressure and output stroke if the effective area of the
diaphragm and the spring gradient are constant throughout the stroke.
Hysteresis due to the stresses in the spring and diaphragm are usually
less than one or two per cent of full stroke. ,
stem assumes a position dictated by the input air pressure. The spring
actuator becomes a power means and the characteristics of the spring
and diaphragm are relatively less important. The use of the positioner
results in several improvements in performance.
The use of a positioner with a spring actuator does not improve the
ability of the actuator tohandle larger inertia or thrust forces unless
special adjustments of motor operating range are made. The only dis-
advantage in the use of a positioner is that it may require maintenance.
The Bpringless actuaior shown in Fig. 7-4 is useful for large thrust
forces. The only difference between the springless actuator and the
spring actuator with positioner of Fig. 7-3, is that the spring of the
actuator is replaced by a pressure regulator which maintains a constant
pressure on the under side of the diaphragm. An air supply at a pres-
sure of 20 to 100 psig is required. The operation of the springleas
actuator is as follows: Assume that the cushion regulator is set to pro-
vide 9 psig pressure on the under side of the diaphragm. At static
balance and with no thrust force on the actuator stem, the upper side
pressure must be 9 psig. Then if the input pressure increases, the
nozzle back pressure increases, and the upper side pressure is raised to
200 Automatic Control
a high value. The actuator stem then moves downward and, as the
actuator stem attains the new position, the upper side pressure is
returned to 9 psig. If there is an upward thrust force on the actuator
stem, the under side pressure remains at 9 psig but the positioner
raises the upper side pressure until static balance is achieved. For a
downward thrust force the upper side pressure is reduced below 9 psig.
Thus, the springless actuator can counteract a thrust force equal to
approximately the underside pressure times the area of the diaphragm.
This is generally from three to ten times the tfirust force handled by
a spring actuator with or without a positioner.
downward, compresses the feedback spring, and returns the pilot piston
to the neutral position. Actuators of this kind are employed in sizes
from 1 to 15 hp and will handle thrusts as high as 100,000 lb.
“ _L ± iMid
2ir \M 2ir V Mx
^ Jl^ /l2 X 100 X 32.2 X 12
4.4 cps
2ir \ 300 X2
Example 7-2. A springiess actuator has a diaphragm of 100 sq in. area.
Its positioner operates from 3 to 15 psig. The cushion (under-side) pres-
sure is set to 5 psig. What range of thrust load can ke accommodated?
Electro-Pneumatic Actuators
When electric control systems are employed, it is often advantageous
to use a pneumatic actuator. a suitable air supply is available, a
If
pneumatic actuator can provide very large power output and may be
operated directly from an electric control system. This requires trans-
ducing the electrical output of the controller into an input variable for
the actuator.
The eleclro-pneurnatic pilot of Fig. 7-7 is arranged to convert an
electrical signal input to a proportional air pressure output. An elec-
tric control system employing this pilot may operate with any of the
Air supply
I
Hydraulic Actuators
Hydraulic actuators, as used for industrial process control, accept a
signal from a pneumatic controller or an electric controller and employ
hydraulic pressure to drive an output member. The hydraulic actuator
is used where high speed and large forces (or large power) are required.
The hydraulic piston actuator in Fig. 7-9 has as its input mi the
position of the vertical lever. For a pneumatic-hydraulic actuator the
input would be the position of a bellows. The balance lever pivots at
the bottom so that an increase of input (to the left) pushes the pilot
piston to the left. This action opens the left end of the piston to supply
pressure and opens the right end of the piston to drain. The large
power piston, therefore, moves to the right until, as the balance lever
rotates about the topmost end, the pilot piston is returned to center,
the motion of the output xi is therefore proportional/ to the input
motion nii. The hydraulic actuator requires a continuously running
electric motor and pump to provide a source of high-pressure oil, and a
drain or sump to collect the return. Hydraulic supply was discussed
in detail in Chapter 6.
204 Automatic Process Control
Oil
Oram supply Drain
Maximum delivery
When the cylinder body is tilted, the pistons reciprocate and the oil
motor is ph 3r8ically the same as the pump but has a fixed ^^tilt” and thus
a fixed displacement. The same arrangement may be employed for
Suction t lorain
Fia. 7-1 1 . Hydraulic motor operator with rotative output.
for
position control (output position proportional to input signal) or
speed control (output speed proportional to input signal). Hydraulic
motor actuators are used in sizes from three to three hundred horsepower.
The electro-hydraulic piston actuator is shown in Fig. 7-12. The input
device is the voice-coil motor which positions a three-land pilot spool.
The pilot controls the flow of oil to the cylinder, and the piston
motor
downward thereby pushing down the pilot spool. This action connects
side to
the lower side of the cylinder to drain and connects the upper
The piston then moves downward, and the feedback
supply pressure.
Tlie position
spring pushes back to rebalance the system at equilibrium.
current.
of tlie piston rod is therefore proportional to input direct
206 Automatic Process Control
Hydraulic actuators require a hydraulic power supply. This subject
has been discussed in Chapter 6, see Fig. 6-19. Generally, individual
power supplies are required at each actuator unless several actuators are
grouped close enough together to use a single large unit. Hydraulic
power lines should not be run for appreciable distances. The pressure
employed in most hydraulic systems usually is from 400 to 1000 psi gage.
Example 7-3. A hydraulic piston actuator with a 4-in. bore and 3-in.
stroke used with a pump having a maximum delivery of 5 gpm. Neglecting
is
acceleration and output load, calculate the time to travel the complete stroke.
_ ir4* X 3 X 60 ‘
^
Traverse t;me = ^ ^ ^ ^ ^3,
= 1-96 sec
1. The bridge type systems discussed in Chapter 6, see Figs. 6-22, 6-24
2. The relay-type reversible motor drive (see also Fig. 6-25)
3. The amplifier-t3rpe reversible motor drive
The relay-type electric actuator is shown in Fig. 7-13. A low-inertia
reversible a-c motor drives the output through gears. A mechanical
feedback is employed such that a rotation at the output rotates a lever
Feedback
signal is both motor windings are energized with the rame phase,
zero,
and the motor is stalled. There are many variations of this general
arrangement. For example, a pair of selsyns (synchros) may be used
inan a-c system. The amplifier may employ a grid-controlled thyratron,
or magnetic amplifiers may be used.
The power of the electric motor actuator is often J hp (about 200
watts) or less. One-tenth horsepower (75 watts) or less is more common.
Thus the electric motor actuator is generally employed for low power,
whereas pneumatic or hydraulic actuators are used for large power.
r.
Force = ^
0.5 X — _X
650 25 X II
j
X 10
5501b
^ I
Output
Electrical motor
Clectrii dl solenoid
That is, the input cannot be driven backward by a toixiue applied at the
output. Such actuators therefore fail in the last position by inherent
action. Failure at the open or closed position can only be a(icompli!^hed
by declutching mechanisms and like.
The thrust force on single-seat valves and butteifly valves may either
assist or prevent closing or opening on power failure. Generally it is
better to arrange that the forces due to pressure differential assist the
action desired. It is well to remember, however, that thrust forces are
not always in the direction of pressure differential.
210 Automatic Process Control
Fia. 7-16. Methods of fluid control, (o) Series throttling. (6) Bypass, (c) Variable
delivery.
= 'KiaV2g(hi (7-5)
Kia = Kx (7-6)
This equation illustrates that the flow rate m through the control valve
is directly proportional to lift;c if (1) the differential head (h\ - /12 ) is
tion formula
L
= (7-«)
The head loss therefore depends upon flow rate. Equation 7-8 may be
written in terms of flow rate:
f-Vw
nhere m= Bow rate, ftVaec. Including head loss (iJii) in equation 7-7,
A 1 ?0 psi 2 1
JS 1 18 IH 3
c 1 16 IH 5
D 1 9 1 12
2g(Ho - H 2)
(7-11)
I 1 -f aix^
Thus the flow is not directly proportional to the valve stem position
and the extent of nonlinearity depends upon line pressure difTerential.
Final Control Elements 215
This equation is plotted in Fig. 7-18 for a 1-in. sise control valve installed
in dififerent size pipelines. The pressure source Ho and pressure tenninsp
tion H2 are considered constant, and the equivalent line length L is about
100 ft. The friction factor F and overall valve coefficient K are con-
sidered constant. The effect of the larger head loss caused by smaller
pipeline size is very noticeable at high flow rates. The smaller lines do
not allow the control valve to pass the desired maximum flow of 35 gpm.
In order to obtain the desired maximum flow rate, four steps may
be taken :
m ^ Ho — H2 — X
(7-12)
M“\ Hi - H2 X
where m= flow rate at any lift
(7-13)
Ho -Hi
216 Automatic Process Control
Combining equations 7-12 and 7-13 to eliminate the head loss term and
solving for the flow rate,
ox
TIM
—a B fraction of total pressure
diff. across control valve
I
when wide open
Fraction of lift
^ Fixed resistance
ijc r
Centrifugal
pump or fan
ing flow rate. This is termed the pump or fan head characteristic and
is illustrated in Fig. 7-20. 1'he head is approximately related to flow
at constant speed by
= 1-/3 (7-15)
Ho
In general, the capacity coefficient p is about 0.2 to 0.4 for many centrif-
ugal pumps and fans.
218 Automatic Process Control
The flow through the control valve is
m jhi - H2 X
M “ \Ihi --H2X
(7-16)
4a-h (1
where a
H -H21 valve head differential at maximum flo Av
Equa-
Ho - /fa valve head differential at zero flow
tion 7-17 exactly the same as equation 7-14 for the case of series
is
E
3
E
X
19
E
:3><:
'C Bypass hand valve
or orifice
2= (7-20)
(7-21)
Example 7*5. A control valve is installed in the exit line of a heat e.\-
100 ~ 15 - 75 10
a = 0.071
lio 140
The single-seat and double-seat control valves are shown in Fig. 7-22.
The single-seat valve has only one port opening between seat and plug
and the entire flow passes through this port. It has the following
because of the differential head acting across the port and seat area.
The double-seat valve has two port openings and two seats and two
plugs. The port openings are not usually identical in size. This type has
the following features: (a) Net force acting on the valve stem is gener-
ally small (therefore ^‘pressure balanced”), (b) It cannot be shut off
tightly because of differential temperature expansion of valve plug and
valve body.
A few types of plugs for single-seat and double-seat valves are shown
in Figs. 7-23 and 7-24. The piston type plug has one or more grooves
seat ring. The parabolic plug presents an annular area to flow between
the plug and seat ring. The poppet-type plug offers a cylindrical-shaped
flow area and
used with small total lift.
is
The gate valve in Fig, 7-25 is often used for fluids containing solid
matter, because it presents an open area directly to the flow of fluid
for solid matter, and because it has no packing gland around the stem.
The flexible diaphragm of rubber or other nonmetallic material is
positioned by the plunger and stem. Fluid pressure inside the valve
body holds the diaphragm smoothly against the plunger.
1.0
2)^- in Size
double- seat valve
0.8
t 0.4
0 02 04 06 "08 10
Fraction of maximum lift, xjX
m (maximum controllable)
(7-22)
m (minimum controllable)
where R = rangeability number. Turndown is a similar concept based
on the ratio of normal maximum flow to minimum controllable flow.
m (normal maximum)
T = (7-23)
m (minimum controllable)
The minimum controllable flow of a control valve depends upon its con-
struction. Clearances must be allowed in order to prevent binding and
sticking, and the flow through these clearances constitutes the minimum
controllable flow. The minimum controllable flow for a single-seat
valve is not zero unless the throttling seat and shut-off seat are identical
and have perfect alignment. The rangeability of a sliding-stem control
valve is usually between 20 and 70.
The importance of rangeability and turndown lies in the application
of the control valve. For example, if the design of an oil burner and
furnace requires a 30 to 1 range of oil flow to accommodate various loads
on the furnace, the turndown must be at least 30 and the rangeability
must be at least 43.
Example 7-7. Flow through a linear valve (constant sensitivity type
characteristic) is given by
m
M
where m is the flow at any lift x, M
is the maximum flow at maximum lift
of 2.0 in. and the rangeability is 20, compute the valve sensitivity.
From the above equation
dm /R - \\M 19 10 gpm
4.75 gpm/in.
dx~\ R /X"20^ 2in.
^ = RMxy-u
M
where the symbols are the same as in Example 7-7. Derive an expression for
valve sensitivity.
226 Automatic Procees Control
1. Rotating-plug valves
2. Butterfly valves
3. Louvers
Final Control Elements 227
m - Kayj^^ (7-26)
Final Control Elements 229
The flow coefficient K and port area a are different for every style or
size of control valve. Consequently, it is standard practice to combine
certain terms of the above equation into a single number Cv, termed the
size coefficient:
Wl (7-26)
where niL = liquid flow rate in gallons per minute at the conditions for
which specific gravity G is taken
Cv = size coefficient
AP = pressure drop, Ib/in.’
GB specific gravity of liquid (referred to water) at either flow-
ing or standard conditions
through the valve ports. The critical pressure ratio for air and steam are
Pt/Pi
Air (70 F) 0.628
Steam (212 F) 0.539
The critical pressure ratio does not vary greatly for all common gases
and a single flow formula may be employed with reasonable accuracy^
to = 158,000 (7-28)
V K I
(h) The value of flow rate that will be required most of the time.
This is the normal flow raie.
*
See Kent, R. T., Mechanical Engtneera’ Handbook^ 12th ed., John Wil^ and Sou
Inc., New York, 19W.
Final Control Elements 231
from
data
actual
Obtain
only.
estimation
(For
manufacturer.).
valvn*.
control
fixing
for
Nomogr^h
7-30.
Fic.
232 Automatic Process Control
(c) The minimum value of flow rate required to sustain the controlled
A’ariable under any conditions of process operation. This is the normal
minimum flow rate,
3. Select the style and type of control valve to provide best operation
for the fluid to be handled. Check the rangeability to insure that the
minimum controllable flow is generously smaller than the normal
minimum flow rate desired.
4. Calculate the pressure differential at the control valve at the
maximum flow. This requires calculating line pressure losses, flow
equipment pressure losses, and determining upstream head from pump
or fan characteristic curves.
5. Determine the control-valve size from manufacturers nomographs,
charts, or slide-rules. If the fluid viscosity is high, or line velocities are
exceptionally low (Reynolds number in the pipe line less than about
10,000) the size coefficient C, will be low, and the manufacturer should
be asked to determine the valve size.
1
“
/ Y (7-30)
2500 Wmi J
where M = maximum flow rate
Mnku = normal minimum flow rate
then the control valve will, at times, be required to operate at less than
2 per cent of total lift. If this occurs, the control valve would require
accurate positioning by a powerful actuator, and the valve ports and
plug must be precisely manufactured. Otherwise it is best to inspect the
control-valve installation to determine whether the differential head can
be increased by selecting larger piunps, larger line sizes, or minimizing
losses in series flow equipment.
Example
1.
7-10. A heating furnace requires a control valve passing 10 gpm
preheated light fuel oil (sp. gr. » 0.8) at full load and only 0.2 gpm at the
emallest heating load. The pressure differential at wide o))en is 20 psi. The
source
2. pressure is constant at 50 psi gage, but there is 10 psi drop in the oil
preheater and 20 psi drop at the furnace burner nozzles. Detenninc the
valve size.
14 gpm
6.
or the turndown is 50.
The ’}-iu. valve is selected because it is the nearest size passing the desired
14 gpm.
The characteristic coefficient is (equation 7-29)
20
“ 50
so that the control-valve differential is 40 per cent of total differential and the
characteristic is satisfactory.
Electric relays for power loads must have appropriate contact life
Series rheostat
*1 - (7-31)
Af Rl +r
where m« power to load, watts
M s maximum power watts
« resistance of load, ohms
r « series control resistance, ohms
^
Re - r ^ Rctni (7-32)
m
(7-33)
M ,^Rc Rc
1
+ icl Kl
This equation is plotted in Fig. 7-32. The rangeability of this method
depends upon the value of resistance of the control rheostat,
R= 1 + (7-34)
(7-35)
12,100
Rl = 8.1 ohms
M 1500
From equation 7~34
E% 12,100 X 36
2300 watts
+ (36 + 8.1)*
PROBLEMS
1.41g ^ \
18 gpm. The line drop is 40 psi and the valve pressure drop is 10 psi. What
size plug valve is required?
Process
Instrumentation
The general block diagram is seen in Fig. 8-1 which includes a feed-
back element as well as an input element. The operational equations
for the system are
d « 1; — c deviation (8-6)
AG G2
1 N
(8-7)
1 + 0 ,(72^ 1 +G G H'*
1 2
AG1G2H NH
8-8 )
, ,
5 -S ff -1- u (
I -h G G2 H
1 1 + GiG2ff
The actual deviation is given by
1 -h GiG2(ff - A) N (8-9)
1 + G1G2H 1 G1G2//
For load changes, the change of set point is taken as zero and the part
of the deviation due to load may be calculated. For set point changes,
the change of load is taken as zero and the part of the deviation due to the
changing set point may be detennined. Notice that t^ese equations
reduce to those of Chapter 4 when the input and feedback elements
contain no time elements {A — H = 1.0).
240 Automatic Process Control
Example 8-1. Derive the general equation for actuating signal.
From equations S-1 through 8-9
A NH
”
* 1 + GiOtH 1 + GiOtH “
The actuating signal differs from deviation (see equation 8-9).
R
U2 - process time constant Tp
Tp-d- 1
R
N = effect of load
T^+l
1
H= measuring time-element constant Tm
1
^ RiTmS + 1)
“ (8-10)
+ (fp + T^)s + (RKc + 1)
^ R - 11 )
.
“ (8
+ (Tp + T„)b + {RK, + 1)
Time Time
the smallest measuring lag; that is, a measuring lag that is as small as
possible compared to the largest time constant of the process.
In pressure control, the process is often of a single energy storage
type. The magnitude of measuring lag is generally small except when
long pressure connecting lines must be used. In this ca&e also, the
measuring lag should be maintained less than a small fraction of the
largest time constant of the process
242 Automatic Process Control
In flow control, the lags of the measuring means and the controlling
means are generally the only important lags in the ^stem. For the
flow-control system of Fig. it is assumed that the measuring means
The lag of the controlling means is neglected in this first example. The
system functions are
II control valve sensitivity
integral control
The system equation for a change of load U2 with load uo and set ]K>int v
fixed is
K2T8(T^8 -f 1)
(8-12)
rr*8* + Ts + X,
The feedback variable is given by
-KjTs
(8-13)
TT„^ + r* + X,
For a sudden (step) change of load variable U 2 the solutions ,
of each of
these equations are plotted in Fig. 8-4. In each case the damping ratio
is selected as one third (therefore ^KiT^ = 9T).
Comparing the response with large measuring lag to the response
with small measuring lag it is seen that the effect is to delay and attenu-
ate the indication of changes in the controlled variable. However, a
most serious consequence of the measuring lag in flow control is that
Process Instrumentation 243
an inspection of the indicated value of the controlled variable b might
lead to the conclusion that the measuring lag of 20 sec produces a better
control than if the measuring lag Is 5 sec (compare the two solid curves
of Fig. 8-4), whereas actually the measuring lag has made the control
much less effective (compare the two dotted curves).
Time Tima
In flow control it often happens that the measuring lag and controller
lag are of approximately the same magnitude, and both are larger than
any other lags of the system. In this case, the system functions of
Fig. 8-3 are
O2 = Ki valve-flow sensitivitv
Kc
Oi proportional control with controller (or valve)
TcS-\- 1
time constant Te
- = + Dcr,* + 1)
+ iT„ + T,)s + (KiK, + 1)
The solutions for a step change in load Uf are given in Fig. S-li for
three cases; measuring lag greater than Tm » 10, equal to Tm^ 5.5,
and less than Tm = 3, the controller lag. The damping ratio is main-
tained at one third, and in addition a ratio of lags has been selected that
results in the same period of oscillation (T^Tc = 30). Thus in each of
the three cases the actual controlled variable c has precisely the same
behavior, i.e., the same damping ratio and same period of oscillation.
The upper curve in the right-hand diagram has been displaced upward
for better comparison.
It is apparent that, whereas the actual change of controlled variable
is the same in all three cases, the indicated variable shows less change as
the measuring lag Tm becomes larger. On the other hand, it is not valid
to conclude that the controller lag should be large because if the con-
troller lag is not the largest in the system, increasing it does not mate-
rially improve the control of the actual controlled variable.
Example 8-2. A flow-control system has a valve lag of 5 sec and a meas-
uring time constant of 7 sec, and proportional control is used. Would doubling
the measuring lag improve the flow control?
From equation 8-14, the damping ratio will be made 1.0 so
(T^ + T^r = miKc + DTcTm
Old KiKc - 0.029
r.r«
KiK, +
Old T = 5.8 sec
ProceM Duturbancee
A number of factors affect the quality of control, the most important
of which is the behavior of the process load. The process load may be
defined in terms of the setting of the final control element to
the controlled variable at the set point.
The location and source of process loads are illustrated in Fig. 8-6 by
the gas-fired continuofis heating furnace. The heat losses are made up
by the stack U4, radiation us, and work c. The heat supply is made up
by work U3 and fuel m.
one valve setting. For another example, the heat may be changed by
varying the temperature of the work as it enters the furnace. The most
conunon load change, however, is caused by variations in the rate at
which the work flows through the process. Other load changes are
those due to stack loss 1/4, caused by changing draft pressure, and radi-
ation or leakage loss U5, due to varying ambient conditions. The effects
of all these load changes must also be corrected by the controller.
The location of the load change has a great influence in the magnitude
of the resulting change in the controlled variable. This is illustrated
by the curves of Fig. 8-7 The curve on the left represents the result
.
able, the magnitude is large compared to the magnitude when the load
change is located far from the controlled variable. The curve on the
right shows the result of a unit change in variable uy (pressure of the
246 Automatic Process Control
fuel gas) in fig. 8-6. Notice that the magnitude of the change is about
two-thirds less.
equation. Derivatives are zero when the operator s is zero, and the
offset for a unit step change of load is given by
Dq = lim -1 (8-16)
•-»o [1 + GiGiH]
Substituting Gi = XeCl + T8)/Ta for proportional-integral control in
equation 8-16 there results
-
Do^liin (8-17)
® ^oT$ + G»H(Ta+l)
where T= integral time. The limit is zero and the offset is therefore
zero if the function G2 H/N does not contain any factors where n > 0.
Do = lim (^18)
•-*0
Thus there is no offset when the function H/N does not contain any
8^ factors.
5 ^ ^
" (Tib DTfi' Tis +l
I
Do » lim
t-HO KcR
+ IHTia l)T^
Do = 0
A disturbance of the set point or input may be required in many
process-control systeips as when a process variable is required to follow
a program, or when the system is a part of another system. The devia-
tion resulting from a disturbance of the set point; i.e., when the set
point is a function of time, is given by equation 8-9,
1 + GiGjjH - A)
(8-19)
1 O1G2H
The character of the input element A affects the performance of the
system. must be noted that the input element does not influence
First it
A ^T8 +
G\ “ 1/T|8
Gz = l/iTi8 + 1)
H= 1.0
Af = 0
T8+ I
’ (8-211
fiTi^ -f Ti8 -f 1
1 1
Setting the integral time Ti to give critioal damping {Ti 4ri) there
lesulte
r« +
(g-22)
(Tae+l)*"
1 -h - A)
Ds = lim (^23)
•-«o 1 + G1G2H
This error is usually zero in process control.
Example 8-4« The static error caused by a unit step change of input may
be written
— G162A 1
Ds = lim Tl
•-»o L 1 + OiOtHJi
For the following functions calculate the static error:
Ke
QxGt ^ = +
WfW’
Ds = lim r.1
•
— +
•-*0
I ^ (r« + i)« 1
Kc TffS
Dfi = 1 -
^
1-4-1
Controller K, = —
Ax
Ae
Am
“ “7“
Control element
Ax
•
Ac
Process
Am
II
Feedback
dm _ g(m)
(^24 )
lb KcKpKn,
or
KcKpK^
X dm (8-25)
g{m)
(gains) are
Controller
Control valve Kp
Process jR = —
Am
as shown in the operational diagram of Fig. 8-8. The system opera-
tional equation is
KcKJt(Ti8+l)
''
TMTa + \) + KcKMTi8 + l)^
RTi8
“ ^ '
TMT» + 1) + KJCMTiS + 1)
For a damping ratio of one third, the relationship among the parameters
(1 + KcKJi)^ 4T
(8-27)
KcKJl “ 9Ti
Since the loop static gain KcKJi is usually at least ten, we may write
KcKJt^^
KcK,^- (8-29)
„ dm 4C
(8-30)
dx 9KcTi
adx *= ax (8-31)
Process Instrumentation 251
The a linear relation between flow and position at the
latter requires
control valve, as shown in Fig. 8-10. If equation 8-27 is solved exactly,
the resulting flow-position characteristic does not vaiy more a few
per cent from a linear relation. It may be shown that controlling outlet
flow instead of inlet flow, and, the use of proportional control in stead of
proportional-integral control provides the same result.
Ax
Controller
Ae
Ac
Process
Am
KcKJC^{T2S + 1)
(Ts + \){T28 + 1) + KcK,K^
+
^ ^ ( ^2
8 )
(T8+l){T2S+l) + KcKJC^
252 Automatic Process Control
IjO
oa
Liquid-level controIxJ
^
0.6
-V
I
a4
/
./ Temperature control
0.2
Position
In thifi equation T2 ,
P, Wj and Kc are constants. The throughput
rate Q, however, affects the heat input because steady-state con-
servation of energy requires
m = QPD (8-3fi)
erf 0.707,
dm
^
WP r/ fzm V ,
(8-36)
dx 2KcTil\WPD/
The equation is now integrable and may he written,
dmi
xi^K (8-37)
fmi* H- 1
T'2m
where mi
WPD
X
Xl
2KcD
K = constant
Integrating,
tan ^
mi = Xi K (8-38)
or
mi = tan Xi (8-39)
magnitude. These lags are also constants. The flow equation for the
orifice or other restriction is
q = KoVmi - p (8-40)
p - Ui - Khq (8-41)
where
(8-42)
\dq/ui
q = Ml — M2 (8-43)
KJC^Ki = A (8
-47 )
L_. L_
'(S’!
1 11
(»-48 )
*
UIMI Ui
M = G262 V2 == Oie
(8^51)
^2 = t>2
“ e = V — b
62 ~ H2C2 h - HiC
The controlled variable is given by
G 2 G3 ^ N 3
(8-53)
1 + G2G3H2 1 + G2G3H2
The inner loop is generally a flow-control system and is considerably
faster than the outer loop which is generally liquid-level or temperature
control. Consequently, the type of control employed in the inner loop
(at G2 ) is not as important as the type of control in the outer loop (Gi).
The following observations are made, based upon a detailed study of
the loop equations:
1. The control action of the inner loop should be the simplest. Omit
integral or derivative responses unless speciflcally desired.
2. Set the inner loop for tight control; i.e., high loop gain.
;
3. H
Measuring lag of inner loop (at 2 ) provides a derivative action
for outer loop. Consequently this measuring lag is not always harmful
as long as it is not excessive.
applied to one side of the diaphragm and the downstream gas pressure
is applied to the other side of the diaphragm. The pressure-balanced
valve is actually a self-operated pressure regulator whose set point is
G,G2AK N
(8-54)
1 1 -hGiG2H
Several important charat^teristics of ratio control may now be shown.
Suppose, by way of example,
(8-55)
K= k
^ / TnS + l \ TtS -h 1
(8-56)
kv ” \TaS+ 1 / T.'l
Kc
In order to obtain the best “following^’' of the controlled variable to the
desired variable, the lags* Ta and the integral time T* should all be
as small as possible. The Kc should be as large
proportional sensitivity
as possible and the control system requires careful adjustment
Batch-Process Control
A process in which the materials or work are stationary at one physical
location while being treated is termed a batch process. This type of
Process Instrumentation 259
overshoot of temperature, when approaching the set point for the first
e fni fll2
K^(l + T^8) ^ 4. JLA
Continuou8-Pro<^8 Control
A process in which the materials or work flows more or less continu-
ously through a plant apparatus while being treated is termed a con-
tinuous process. The problems of continuous-process control are caused
by load changes.
Flow rate of the materials is almost always important in continuous
processing, particularly where quantitative reactions are involved as
in blending. On the other hand the supply of material is sometimes not
262 Automatic Process Control
coiistaut and a flow control is necessary as shown in Fig. 8^16. The
resulting flow of materials is nearly constant in spite of large changes of
head in the tank. It is necessary to size the tank only large enough so
that the longest charging period will not run the tank over, and the
longest off-tinie of the inlet will not allow the tank to run dry. In order
to keep load changes at a minimum, a short tank of large cross-section
area (large capacitance) should be employed.
1 1
n 1 1 1
1 1
Time-
Ac = — q (8-58)
the pressure differential at the valve is constant. The inflow, head, and
outflow are shown by the dotted curves in Fig. 8-17. Obviously if
there is no control, the tank would ultimately fill or run dry. If control
of level is added,
-q = Kci^v -- c) (8-59)
Process Instrumentation 263
where Ke is the proportional sensitivity of the proportional controller.
Combining equations 8-58 and 8-59,
1
9 ^ (8-60)
u 1
and
c 1
_ 1
(8-61)
Ic, Ts-f 1
where T = A/Kc^ The inflow, head, and outflow are also shown in
Fig. 8-17. The time constant T may be selected by adjusting the pro-
portional sensitivity of the controller so as to ''spread out" the inflow
ables constant does not allow variable inlet condiiiions or variable flow
rate. Economy can be achieved by employing a computing machine
and relating all variables, even in a transient state, so that the desired
criteria are met.
2. The process may operate in such a fast time scale that human
action is always too late.
3. The process may operate in such a slow time scale that human
attention is ineffective.
There are two difficulties t-o be overcome with this method of control;
first, the behavior of the system may 'ot be known in general analytical
form, and sen-ond, the sy.st{*m jna> have long time delays. In order to use
the continuous computing method the system behavior must be known
to some extent The and em])iri-
behavioi' can be tested experimentally
cal reliilions found either by transient or sinusoidal response methods
or l)y correlation methods.
Often it is sufficient to know only the general form of the behavior
(as the geiicial ditTerential e<iuations with unknown coefficients). In
this cose a self-checking method can be ased I\eferring to Fig. 8-19,
system variable Cm is nieasurcfl and is fed directly to servo S 2 which sets
the unknown coefficient K of the reaction. Computer C 2 takes data on
temperature, pressure, and coefficient and computes the \alue of the
system variable that should exist, C^. T1k‘ computed \alue is fed back
to the servo, ami a new coefficient determined until computed and
measured values match, d'he form of the system reaction must be
assumed and fitted into computer d'he “self-cheeking” elements,
S 2 and Co, insure that the assumed equations fit process behavior in the
1
PROBLEMS
8-1. 'I'he system block diagram ran be redrawn with GiGiH in the for-
w'ard loop when a feedback element // is present. What function is required
at the input element?
8-2. A system (see Fig. S-1 ) has the following elements.
A - 1 A - I
G, - (1/T,.9)*
G =
2 Kr, C3 = A', Hz =
7’.s- -h 1
= = A (fA - -
Gi K], G‘2 2,
-b 1
'
7Wl
I
N, - A',, Gi A 4 = 0
+
’
"1
r,.v
show that the off.set with load W3 is greater when the inner loop is broken
(Ao - 0) ,
8-8. For tlie storage vessel m Fig S Hi, sketch the outflow as a function
of time with the, inflow as shown if (a) integral control i.s used an<l (/>) propor-
For the storage vemel in Fig. 6-17 with a step inflow of U units,
vessel area A,and proportional-integral control, calculate the offset in flow rate.
6-10. For the surge vessel in Fig. 8-17 with integral control T,-, and no valve
or measuring lags, calculate the integral time for critically damped performanoe.
6-11. For the surge vessel in Fig. 8-17 with proportional control Fe, a
valve lag Tv but no measuring lag, calculate the proportional sensitivity for
critically damped performance.
cha/titT
Sinusoidal Analysis
turbance at one point in the system, holding all variables but one con-
stant, and noting the result upon the desired variable as indicated in
1 c
2^ + 1
do
di
+c= sin (9-2)
AcjT cos ut - BcjT sin col A- A sin cot + B cos cot = sin cot (9-4)
A -—1
1 +
-6
(9 )
n
1 +
The particular integral is, therefore,
c = -
—“T ;
2 ^^ ‘*^0 (9-7)
1 W*'i
c = —" T— T- sin
'-
(cot
— tan“* coT) (^8)
V\ +
Comparing the output c to the input m (equation 9-1), it is seen that
the output sine wave has lesser amplitude. The ratio of the amplitudes is
designated by the absolute value signs thus-
I
(9-9)
fn V\ +
It is also seen that the output sine wave is delayed and the phase is
—tan ^
uT (9-10)
wbere Im means *^the imaginary part of." The function above is de-
fined by Euler^s relation,
T^
at
+ c = Im e*“‘ (9-13)
I
r^
at
+c = e"‘ (9-14)
c = (g~15)
+ l)Be'* = 1 (9-17)
Solving for B,
1
B= (9-18)
(i«T + 1) c'*
Next solve for the real part in the following series of manipulations:
/ 1 - iaT\
(9-21)
Sinusoidal Analysis 273
or
f
£. - (9-25)
m “ Ta
+ 1
Replacing s by tea,
c 1
(9-26)
in iwT -f- 1
1
'
- (9-28)
- \ 1 +
The phase of a complex number is the arctangent of the imaginary part
divided by the real part,
Example 9-1. Calculate the amplitude ratio and phase of the sinusoidal
response of the system described by
c 1
m " Ts(Ts -h 1)
Substituting « = iu,
1
^
m iuiT{l + iuT)
Clearing,
^ __
^‘(1 — _ —<j)T ~ i
m“ "
o)T(l + wTd -f (^f^)
The amplitude is
1 .
1
/
\ (1+ + ^2^2)2
or
f
w I
I wT V 1 -f
The phase is
The sign of the phase must be considered carefully because the inverse func-
tions are multivalued.
z = a + t6 (9-30)
and sketched in Fig. 0-3 by marking off length u t)ii he real number
is t
axis,and marking of! the length 6 on the imaginary number axis. When
a complex number is written in polar form, only the magnitude and the
angle need be given.
The magnitude is calculated from Fig. 9-3.
= M (9-31)
Sinusoidal Analysis 275
and the angle is given by
/S = » = tan-* “ (9-32)
-h S 2 = a -h + c 4- id = (o + c) 4- i(b + d) (9-33)
(a 4" tb)(r 4' id) *= (ac — 6(f) 4- i(6c 4" ad) (9-34()
so that the phases are added. Therefore the product of two complex
numbers is obtained when the magnitudes are multiplied and the phases
are added.
The conjugate of a complex number is obtained by changing the sign
of the imaginary pari,
f = a 4- t6, z* = a — ib
276 Automatic Process Control
Example 9-2. State the magnitude and phaae of the following numben:
Phase
Number Magnitude Degrees Radians
4 -i3 5.0 -36.8 -0.643
+ v^ + v^ 2.0 +45.0 +0.785
i 1.0 +90.0 +1.67
-10 10.0 180.0 3.14
plot the amplitude and phase, since it is known that the function is
wuaoidal. Furthermore the amplitude and phase arc a function of
frequency so that sinusoidal response may be considered in the frequency
domain rather than the time domain.
Sinusoidal response may be plotted in three different ways: (1) The
rectangular plot with amplitude ratio versus frequency, and phase versus
frequency; (2) a polar plot with magnitude and phase shown in vector
form with frequency as a parameter; (3) a phaie margin plot with mag-
nitude shown versus a function of phase with frequency as a parameter.
Fr^iuency, wT Frequency, «T
Fio. 9-5. Magnitude- and phaae-frequency plot.
(9-40)
(9-41)
asymptotes:
Log|G| . -ilogd (9-42)
At high frequencies, tuT » 1.0, the magnitude ratio is a line from upper
left to lower right with a slope of minus 45 deg because as wT —> •
[G = tan“‘ - wT (^6)
and is also plotted in Fig. 9-5. This graph is made on semi-log paper.
Because of the arctangent function there is no linear approximation of
the curve but, however, the function is symmetrical about the wT = 1.0
point.
The polar plot of real versus imaginary part is also used in stability
studies. From the complex form of the system function,
1 - iu)T
G (9-46)
1 + iu>T 1 -h
The real and imaginary parts may be calculated for different values of
uiT and the locus of the end points of the vectors form a curve known as
the system locus. The plot is as easily made by plotting magnitude
and phase. The system locus is shown in Fig. 9-6 for a time-constant
elen\ent. Note that the frequency wT becomes a parameter marking
out different points along the locus. The magnitude is the distance
from a given point to the origin, and the phase is the angle taken counter-
clockwise from the positive real axis.
On the rectangular and phase margin plots, the magnitudes arc set
against a logarithmic scale, so that magnitudes are multiplied by adding
vertical distances. The calculation of phase is made from
(9-48)
r( ) = K
^
-f
• - •
-f ais + Op) (9-49)
•
H- -h •
biS bo)
where fc, w, and n are positive integers. The total phase at sufficiently
high frequencies approaches
|g| = « if ik > 0
lel-^
••-*0Oo
if - 0
|Gl=o if A: < 0
«r-^
The rate of change of magnitude may also be determined since the rate
corresponds to the units of phase lag. In terms of the slope when the sys-
tem response is plotted on log-log coordinates,
Slope of IG'I
—m -n - k (9“54)
and
Slopeof 1^1 = -ib (9-56)
4#-*0
1 1
Q= (9-66)
TsiTs -h 1) Ts -f 1
Magnitude Decibels
0.10 -20
0.50 -6.02
0.707 -3.01
1.0 0
2.0 6.02
10.0 20
282 Automatic Process (Control
(s ^
The result is plotted in the figure shown below:
215* 4-105-1-1
(75 4- 1)(35 + 1)
The response may now be plotted as two time-constant elements. The result
is plotted in the above figure.
These are plotted in i'ig. 0-9. The magnitude versus frcfiuency plot
shows the magnitude to decrease inversely Avith fre(|uency. Note that
the slope is -15 deg (1:1). The phase lag, PL is constant at 00 deg
forany frequency. The polar plot is a vertical line along the negative
end of the imaginary number axis. Th(* magnitude-phase margin dia-
gram is also a vertical line of infinite length.
The time-instant element has the operational function
«=
r/+i i+lr
The magnitude ratio and phase are,
These are plotted in Fig. 9-10. The magnitude versus frequency rela-
tion approximates l.O at low frequencies and falls off at high frequencies.
Again the slope is - 15 deg (1:1) at high frequency. The phase lag is
small at low frc(iuencics and approaches 90 deg for high frequencies.
At the frequency wT = 1,
1
<? 1 = 0.707, Z:G=_45deB
284 Autootatic Process Control
The polar plot is a semicircle as shown. The phase margin plot begins
on the right and passes downward.
The oscillatory element has the operational equation
(9-63)
TV + 2fTs + 1 (1 - «*T*) + »(2f«T)
Va -«*F)5 +'(2^r? ^ (i -
(9-64)
These are plotted in Fig. 9-11 for several values of damping ratio {.
The magnitude is greater than one for moderate frequencies and low
Sinusoidal Analysis
/^--90deg ^
0-K.(,+±)^K. (1 -h uoT)
^ + ( 9 ^)
ia>T (cT
/G = -tan-^
"(1 +ic’‘T/)
|G| = (9-70)
—“ tan ^
w7'j “f” tan ^
(9-71)
G= T8-*xoiT
The amplitude and phase are
The polar plot is a vertical line out along the positive end of the imaginary axis.
^ - 1 = Oj
Example 9-9. A thermometer and thermal \Nell have two time constants
where Ti - 100 sec, T-i — 10 sec. Plot the ainplitu<le-frequency response.
The plot of magnitude may be determined by using the approximate
(asymptotic) method: making the graplis on log- log paper, the response due
to each time constant may be added separately. The two response.** are then
combined additively to obtain the asymptotir resjionse The mterseciiuu
Automatic Prooees Control
frequencies are found from
1
0.0016 opa
2ir!r,
1
h 2rT2
0.016 ope
^ X 96 10
The amplitude is
« X 80 = 8 deg
G (9-72)
where L is the dead time and e is the base of natural logarithms. This
relation Chapter 2, by classical methods.
was derived in
The frequency response of a dead-time element is calculated in the
usual way. Replacing s by ia>,
-L$ .
^ (9-73)
These are plotted in Fig. 9-15. Dead time has the characteristic of a
pure phase shift since the magnitude is always one. It is easily seen
Sinusoidal Analysis 289
from the graphs that no other dynamic element builds up phase lag so
much as dead time. In combination with other dynamic elements the
phase lag of dead time is simply added to obtain the total phase lag.
A useful approximation for dead time L is the function from the
Fade list,^
2 -
G
U (9-^75)
2 4" IjS
The magnitude of this function is eveiywhere 1.0 and the phase lag
corresponds very closely to dead time up to frequencies given by
wL — 0.5. In addition, this function may be easier to incorporate into
complex system analysis.
Thermal process control may involve elements having distributed
parameters; that is, elements with distributed thermal resistance and
distributed thermal capacitance. Pneumatic systems having long trans-
mission lines may also be similar. As a typical example, consider the
heat transfer through a solid material such as a thick wall or tube.
The temperature distribution is governed by the equation,
d^c
(9-76)
=
^+ (a^s)c 0 (9-77)
r = (9-78)
a = ± (fr-79)
c = Ae^ + (9-80)
'
J. G. Truxal, Auio?naiic Feedback Control System Syntkme, McGraw-Hill Book
Company, Inc., New York, 1966, p. 560.
290 Autoinatie Process C'ontrol
c « C|, X » 0
ac
0 ,
X =» L
dx
A and B may be evaluated from the boundary conditions, and the result-
ing solution is o
+_e^\
c = Cl (9-81)
+ 1 /
<=1 ( 9-82 )
Ci
_ ( 9-83 )
or in hyperbolic form
C2 1
(^84)
Cl cosh V —o?L^s
This is the system function for one particular distributed parameter
system. Other boundary conditions yield a similar function, usually
one multiplied by a complicated ratio of hyperbolic functions.
In order to expand the system function it is necessary to write the
complex number a in a different form, employing the square root of t.
V V -a^LHuf = Va^Z?w V -I
The operational function is therefore
5 1
( 9-86 )
Cl cosh /3(1 — t)
where 2/5* * a*L*«. This function may be expanded and the magnitude
and phase derived. The magnitude is
C2 I 1
(f)-87)
Cl (cosh* P cos* p + sinh* p sin* p) ^
Sinusoidal Analysis 291
These are plotted in Fig. 9-16 where wT = 2/8*. The magnitude drops
off exponentially with the square root of frequency, and the phase lag
(9-89)
Tos -f 1
c%=at
-OT +? (»-90)
( 9-01 )
PL
2irfi
^0
1100 ^ 2ir'
culties. They may be used with adequate results when the element in
Bi = sin
‘ — if Xa ^ Xi (9-93)
2d4 Automatic Process Control
phase is the angle of the fundamental of the output referred to the angle
of the input. It is not very difficult to write the Fourier series for the
limiting as may be seen in Fig. 9-19 where dead zone cuts off the
“bottom’* of the sinusoidal signal. Limiting cuts off the top of the
signal. In dead-zone action, the output signal begins at an angle given by
As l>efore, the magnitude and phase may be calculated from the Fourier
Sh - sin“ (9-97)
Idle output signal continues until reversal of the velocity of input occurs.
The output signal then holds its last value while crossing the hysteresis
loop The inagiiitudo latio is
Example 9-13. A control valve has 2 per cent dead zone. The amplitude
(>1 signal generally oc(‘urriTig is o per cent. Estimate the gain.
Sinusoidal Analysis 297
From equation 9-95,
Bd * 8m“‘ I = 0.41 radian
'fhe magnitude is
not at all necessary that the input signal be purely sinusoidal. In fact,
a rectangular or triangular wave can also be emiiloyed.
The sinusoidal signal may have to have mechanical, pneumatic, or
electrical form. Thus a sine-wave generator is required. The frequency
(cycles per second) of the signal may be very approximately as follows:
There are two methods of determining phasf' lag. First, simply record
input and output and scale the difference. Second, provide a variable
jjhasing of the recording of the input so that shifting phase of the input,
to match the input and output records enables reading phase lag directly.
:
PROBLEMS
9-1. Calculate and plot the amplitude and phase of the following func-
tions
1
q = KdC
«X,= |c
Stability Analysis
( ^
+ -1 /
' fill \)
* 5,-'
the La Place transform of the equation is
/l.s'" 4- Bs + C
regardless of any initial conditions (yo, yo)* Suppose that the initial
conditions are zero {ya - 0, ya ~ 0), and tliat the system is disturbed
starting from “rest” condition.
1
£.r ( 10- 4 )
As^ 4- Bs + (
Sy 1
Stability Analysis 301
Thus there no essential diffei'ence between the system operational func-
is
= 1
( 10 -6 )
&is) = £( /) =
2 (10-7)
As- -h 5s 4- C
1
("apacitance element (10-8)
Ts
1
Time-constant element (1(^9)
Ts -f 1
Oscillatory element, i
^ 1.0 (10- 10)
Th- f + 1
1
Dead-time elemeni (10- U 1
p
“ r V + 2fn +
Fio. 10-1 . Generation circuite for time-couetaut and uacillatory elements.
The transfer function for dead time (equation 10-11) may be derived
by use of La Place transform. Consider a dynamic action such that the
disturbance (input) is
The transform is
^ = m ( 10
- 12 )
£x^£f(t) (10-13)
The only action by the element is to delay the function in time without
otherwise altering it. Thus the response (output) can be given by
»-=/(<- L) (10-14)
£y = [ d{t + L) (10-16)
Jq
or
r•
me-'dt (10-17)
Jo
£y = (10-18)
G(8) (10-19)
D l-OiG2A+GiG2n N
Closed loop U (10-23)
1 + G G2 H 1 1 + G G2H
1
(10-25)
6„a" + bis + bo •
where 0^ and &n constants and m and n are integers. In any real
physical system, furthermore, n ^ m. The denominator of the transfer
function may be factored.
„ 0**’" + • •
•
aiS + Oo (10-26)
{» - ri)(« -Ti) (« - r.)
3f)4 Automatic Process Control
whi le are the roots of the denominator polynomial. Whether or not
any roots are zero, or alike in value, the inversion of this transform can
only yield terms like
where r,i are the roots of the denominator polynomial and may be real,
G =
" + ^
.S'‘
-f - 4
Factoring the denominator,
1
^
-f - 1)
Since one of the iactorH luis h pi>Ritiv(; root — 1), the system represented (.s-
(\
a - d-
(.^ d- ,S- - f
The unknown cocthcicnt.s arc determined and substituted into the partial
fractions’
r ^ i {
9\(« +
^
2)*
2
5-h2^.s—l/
M
Using the La Place transform table to get the inverse transformations for a
unit impulse input,
c - - 2c“*' 4- 2en
Since the third term has a i>ositive exponent, c increases exponentially with
time, i.e., the output has no steady-state value, and the system is unstable.
Stability Analysis 305
-f
• •
-f k, - 0 (10-28)
1. All terms of the chariieteristic equation are not of the same sign.
2. All terms of the characteristie equation are not present in descend-
ing powers.
Either of these conditions may b(* detected by inspection and may l>e
> 1 '^ ^
Using rous two and three, tlie coedii‘i« il- ul liu‘ louitii low are tound
by cross-multiplying in the same inannei-, and now rows are formed in
this way untilno terms remain. In tin* course of development, the coeffi-
(1) all the coefficients of the characteristic erjuatiun.- are of the same
sign, and (2) all ihe derived c<x*friei(Mil'; of the first column of the anav
aie of oiu >ign
306 Automatic Process Control
There are two exceptions to the process:
(a) If the first-column term in any row is zero, but the remaining
terms are not zero, replace the zero by an arbitrarily small number e
and proceed.
(b) If all the coefficients in the second or any derived row are zero,
there are roots of equal magnitude lying radially opposite on the real-
imaginary plane. In this case form a polynomial with the coefficients
of the last row, take the coefficients of the derivative of this polynomial
for the next row and proceed.
(«* +s+
4X7-1X8 .5
1
4
2)(.< -h 2)(« -f 1)
,
«
4X4-1X0
SO
-f 4«* -f 7«*
7
+
4,
8s +4
4
0
= 0
5 X 8 - 4 X 4
^
^
5 5
The coefficients of the equation are all positive, and the left column is entirely
positive, so the system is not unstable.
3 2^' - 1
2 - * A: + 1
The third term of the first column requires that k < 2, and the last term
requires that k > Since these conditions cannot be fulfilled simultaneously,
2.
there is no value of k allowing stability.
Stability Analysis 307
C* y ^ TJ
1 4 G1G2H
“ 1 + G1G2H
As stated before, a control system is unstable when the denominator of
the closed-loop transfer function has roots with positive real parts.
Study of this condition is made easier by considering the actuating
signal rather than the controlled variable,
£= ^
R- U (10-29)
1 4" G 1 G2 H 1 4" G G2 H
1
I + Gi{s)G2{s)H(8) (10-30)
now that the tip of the vector follows the dotted line in a counterclock*
wise manner where the path of the end of the vector travels from posi-
tive infinity down the imaginary axis to negative infinity and around
a circle of infinite radius back to the starting point. Thus, the value of
complex variable s changes from
-fioi as w— 4- 00
to
- Ki) as w — — 00
and back again. In this case it is desired that th*e real part of s be zero
(o = 0) so as to make this complete encirclement of the positive half
of the 8 plane.
For a single positive root (as r2 ) the vector s - r2 makes one counter-
clockwise rotation, but for a negative root (as ri) the vector s - ri
makes no net revolutions. Therefore it is concluded that a root can
be located by the simple scheme of determining the number oa'’ rotations
of a test vector when .s == lu in the particular path encircling the right
half plane. If there happens to be a root lying squarely in the imaginary
axis or at zero, it is possible to alter the path of 5 ever so slightly around
these points and proceed in the normal manner.
It is desired that s take on values given by tw, so
and tkt resulting function can be plotted on the function ()lane as shown
in Fig. 10^-4. The problem may be transferred from the a plane of
Fig. 10-3 to the G(s) plane in Fig. 10-4 b, conformal mapping.
Suppose that the dotted ^ine in Fig. 10-4 represents G G2H
1
for dif-
on the locus r present-s the vector G 1 G 2 H. The line from the origin to
the (-1,0) j)(jint Q is the vector ( —1 ). By .subtraction the line joining
QP is the vector 1 4- G G2 G
1
Thus it is concluded that we ne(‘d only
plot the locus of the open-loop function,
Gj (lCo)G2 (fw)f^
and a vector from the point ( -1, 0) repre.sents the denominator of the
transfer function. The transfer locus is a plot of a function of variable a
convenu'm only the values from positive infinity to zero need be plotted
o,
s*(r« + 1)
310 Automatic Process Control
As shown in the figure, it is necessary to connect the ends of the locus by
making
1
Q1O2H
»(Ta -h 1)
Im
Stability Analysis 31
t (3.
and if the sum of the net number of counterclockwise rotations of the vector
1 G G2H and of the
4- 1 number of poles of the function G 1 G 2 is zero. H The
vector 1 -h G G2 H is considered
1 to follow the transfer locus G G2H from
1
1 + G1G2H
In many problems in automatic control it is interesting to note that the
number of poles and zeros is zero so that the Nyquist criterion reduces
to requiring no net rotations of the vector. The criterion then may
simplify to the statement that the transfer locus does not encircle the
(- 1, 0) point.
- rfi
GiGfH does not have a pole so P = 0. The transfer locus in the figure shows
H- 1.0 counterclockwise rotations. Thus ^
* 41.0
AT
«
P 0
312 Automatic Process Control
and
iV -h P 0
A. Plot the open-loop transfer locus GiG^H on the complex plane for
frequency varying from zero to infinity.
(a) If experimental data is iised^ plot the locus directly from the
phase and magnitude data.
Stability Analysis 3ia
(6) If the transfer function GiG^H is given, make sure that the
denominator is of order equal to or greater than that of the numer-
ator. If this is not true it is necessary to go back to the analysis and
pick up some of the higher order terms that have been neglected.
(c) If the locus does not terminate at the origin for infinite fre-
quency, then draw an arbitrary line extending the locus to the origin.
(d) Plot the locus for positive frequencies. The plot for negative
frequencies will always be the mirror image of the plot for positive
frequencies, Reflected about the real axis.
(e) If there are rf factors of s in the denominator of GiG^Hy then join
the open ends of the locus of G G2 H with
1 nr radians of arc (at infinite
radius), going counterclockwise.
B. Draw the vector from — 1 to a point on the locus. Let the point
vary from «=-!-«> to <*>
= — and count the net number of counter-
clockwise rotations of the vector. This is N.
C. Determine the number of poles in the function Gi(j 2 H. 1'hese
are the number of values of a having positive real parts (counting each
repeated value) that make 1 + G G2H infinite.
1 This is P.
D. The system is stable if and only if
AT +P- 0
The complex plane diagram can be used for more complete investiga-
tion of system performance short of solving for transient response For
example, it is often required to find the closed-loop performance know-
ing the open-loop performance. For any control system,
(open loop) (
10- 34 )
and
C
(closed loop) ( 10 -36 )
V
The analysis is made a great deal easier it the inverse functions are used
The inverse of a vector is obtained by inverting the amplitude and chang-
ing the sign of the phase. Let
(10-36)
Then
V
(10-37)
C
Ck>nBequently, as shown in Fig. 10-5, all that is necessary is to plot the
function Y and the function H/A The output-input response C/V is
.
. Automatic Process Control
rfl4
this meUi^
the final locus. Unfortunately
eaaly obtuned by inverting
measuring-element transfer function
is useful if the feedback or
of the test vector is two (iV = 2) revolutions, and since there are no poles
the system is unstable. If the proportional sensitivity is reduced to about 2.8,
the system becomes marginally stable.
B iCCeOs -f 1)
E* 60«(120s + l){30s -h 1)
This function is plotted in Ex. Fig. 10--11. The net counterclockwise rotation
of the test vector is zero and N= 0. Since there is a pole in the right half
plane at s = 1.0, therefore P= -fl O. Thus
N -{-P 9^0
Tb{T9 + 2{)
The open-loop magnitude is
i
lei
«7’V(4f* +
and the open-loop phase is
/G = tan ^
For the action within the loop, the influence of the load variable is not
important, so
= (10-39)
I
Stability Analysis 319
The maisnitude ratio and phase,
1. If the phase lag is less than 180 deg at the highest frequency for
ably over two octaves of frequency near the highest frequency at which
the magnitude is one, and if the negative slope of the magnitude fre-
320 Automatic Process Control
queiicy curve is less than 2.0 (plotted on log coordinates), the system
is probably stable. Conversely if the negative slope is greater than 2.0
C _ lU
(10-44)
E “ (16«+ 1)^
C\ 10
(10-45)
The phase is
-3tan~'a>7’ (10-46)
The magnitude and phase are plotted in Fig. lQ-6. Stability checks
1. The system is just unstable because the pha.se is slightly more than
180 deg when the magnitude i.s oiie (190 deg at 0.0191 cps)
: :
ward until the slope i^ 1.5 at the point at which the magnitude is 1.0.
c = -f 35e®
o®®'
sin (0.1 16f - 31®) (10-47)
A unit impulse disturbance was assumed and the roots of the character-
istic equation found by trial division
This solution is a stable one, as shown in Fig. 10-7. The damping ratio
of the oscillatoiy component is 0 25 and is generally conddeied
satisfactory.
Example 10-12. Determine the integral time for integral control with
two time constants and 0.625 sec. The open-loop function is
of 2.5 sec
C 1
E !r«(2.5« + 1)(0.626« -h D
The magnitude and phase are
—
wT \ (1+
/
6.25a)*) (1
^
+ 0.39a)*)
—90* — tan“^ 2.5a) — tan“^ 0.625a)
First, the phase angle does not depend upon the integral time T, The fre-
quency for 135 deg phase lag is easily calculated;
0.28
For a magnitude of 1.0,
1
ojT
v'(l + 0.49M1 + 0.0306)
and
1
2.9 sec
0.28 X 1.24
C K
""
E 6»*(20s -h 1)
A'
Vl +400u*
For 135 deg phase lag, the frequency o) must be about 0.20 radians f)er sec.
The amplitude must be one, so
PROBLEMS
IM. Check the following Bystems for stability. The function given is
1 + GiGiH.
(o) 6«» + n<* + 6« + l (c) lOOsH 210»‘ -I- 121s* + 22s + 1
A’,(r^ + 1)
as given in Chapter 4.
(?.+ !)•
and the loop closed. Does the open-loop gain depend upon the value of the
time constant?
10-4. A controller and process function are
Aj
c+^'TsfT's -I- 1)
apppndix
of reference input.
Controlled system, G2, is the body, process, or machine, a particular
quantity or Condition of which is to be controlled.
Controlled variable, c, is that quantity or condition of the controlled
system which is directly measured and controlled.
CofUrolling means are the elements of an automatic controller which
are involved in producing a corrective action.
Corner frequency of a transfer function is the frequency at which hne>
asymptotic to its log-magnitude vs. log-frequency curve intersect.
Dead time, L, is a fixed inter\^al of time between the change of an input
to an element and the beginning of response to the input.
Dead zone, Xj, is the largest range of values of the input variable to
which an element does not respond.
Derivative action is a controller action in which there is a continuou-s
linear relation between the derivative of the actuating signal and output
signal of the controller.
325
326 Automatic Process Control
variable and which is compared with the reference input to obtain the
actuating signal.
Proportional action is a controller action in which there is a continuous
linear relation between value of the controlled variable and the value of
the output signal of the controller.
Proportional sensitivity, Ke, is a proportional action adjustment. It
isthe steady-state ratio of the change of controller output variable and
the change in actuating signal.
Rangeability is the ratio of majdmum flow to minimum controllable
flow through a final-control element.
Raie time. See Derivative time.
Reset rate is the inverse of integral time.
Reference input, r, is a signal established as a standard of comparison
for a feedback control system by virtue of its relation to the command.
Response time of a system or element is the time required for the out-
put to first reach a specified value after the application of a step input
or disturbance.
Rise time of a system or element is the time required for the oui^mt to
increase from one specified percentage of the final value to another,
following the application of a step input.
Servomechanism is a feedback control system in which the controlled
variable is mechanical position.
Set point is the selected reference value of controlled variable which
it is desired to maintain.
Settling time is the time required for the absolute value of the difference
between the output and its final value to become and remain less than a
turbance.
Static error. Eg, is the steady-state deviation resulting from a change
in value of the set point.
Time constant is the time required for the output of a first-order
system to change from a given value to within 36.8 per cent of the final
value (or 63.2 per cent of the amount of total change) when a step change
of input is made.
Transfer locus is a plot of magnitude and phase angle of sinusoidal
response of a system.
FtG.
A-.
Transient
+
rntimnw
2rr..c
?r
>
do
system
-
0
to
a
step
input.
input.
impulse
unit
a do-VT’c
to
system
co»0.
second-order
with
0
=
h
of c
-f
response
2f
+
rc*c
Tranrient
A-2.
Fig.
» 'ap.-i^iuSe^
Sinusoidal Response 331
332 Automatic Process Control
La Place Transformation
The solution of linear integro-differential equations by ordinaiy
methods requires obtaining (a) a characteristic equation, (&) the roots
of the characteristic equation, (c) the complimentary solution, (d) par-
where the function z(t) is in the time domain before transformation and
.Y (s) is in the s domain after transformation. S 5mibolically this trana-
formation is written
. 2 . (0 = t
Y(,v) = /
(U
Jo
X{») = f-
p
Thus, it is possible to evaluate the transform by direct integration. In
order to save time a table of transforms, such as Table A-1 or Table A-2
may be used. The table of direct transforms is a set of functions that
may be transformed by direct integration as may be seen in Table A-1.
The inverse transform table is simply a set of solutions with arbitrary
coefficients for a corresponding set of simple differential equations.
)
can be found using the theorems of Table A-1. The disturbance u(l)
ife the step function
u(t) * 0 , / < 0
ft(i) = 1, / ^ 0
-
(a2«’ -f njji -f ao)X{s) - //o -f )
4* ni.r(0-f
+
A («)
= yo 4- (a2S^
7“tt:
«(a2«4- ai» —T—
ais)x(()+) 4- O2 «j:(0+)
4- uo)
^
X(sK'<ys (A-3)
mfi'” 4- 4- • ' •
X{8) = - -
(A-6)
(s - ri)(5 - r2)(« r3j • •
•
(6 r„)
334 Automatic Process Control
Operation or
Function Function fit) Transform F(«)
Derivative -m
d*
^/(O «*F(») - ^(0) -/'(O)
1
Unit step u(t)
8
1
Ramp t
n/
r gfi+i
OJ
Sinusoid sin u)t
«* -f
8
cos 0)1
1
g-df
Exponential
« -h a
where the roots (rn) of the denominator are found, including zero and
repeated roots. The transform can now be expanded in a series of
terms,
f'M /«)
1
/ sin \/r^«/7')
rv+2f?’s+i T’Vi-r*
«-f"’'sin(V'l-rH«/7’)-*l
1
«(T*s*+2frs+l)
where ^ - tan *
^
(7’i/7’*)e-‘"''
l-2f(7’,/7’)+(r,/7’)*
1 (7’)-‘e-'"'8in |vT^«/r)-*l
(r,«+l)(7’V+2fT<i+l) + (l-mi-2f(7’/7’,)+(7'/ri)‘l^
. . .
(T'./T’Ov/i-r’
wher«*=ten'
,
1 rj-2f(r2/7’)+(i’,/r)qx._„„
tL 1-C J
T 2 S-I-1
8in(Vl-f*(</7’)+*)
(rv+2r7’«+i)
,
(7’t/r)Vr=r
where
,
<^-tan ,
|l-2f(r2/7') + (7’2/r)»)We-f'’’
TasH-l
sin l\/rT*(«/7’)+'*’l
<i(TV+2f7'»+l)
. . ,
.(Tt/T)Vi-r ....yi-f’
where-l-tan
-f
iifl
B(8)
_di_.
S - Si
+ -A_
S - S2
+
S
^
“
(A-7)
33G Automatic Process Control
Note that a /oro root may be treated exactly like a nonsero root. In
order to solve for the A's, multiply the transform by the denominator
factor and set s equal to that particular root. Thus,
AW
® (s - Si)(s - 82)(8 - «4)
• *
•
(« - «n) «
_ Aj Az A4
.s + 2 fi -h .3 s + 4
.s-
+ 1 -1 I
(A- -h 4)(.S- -h 3) 2 X 1 2
«+ 1 -2
Az ^ -1-2
(8 -h 2)(s + 4) -1 X 1
5-h 1 -3 3
Aa
(« + 3)(s + 2) -^l X -2 2
A{&)
fi(«)
is
2 (« +
12 2) -f 3)
^
3
2
1
(8 -f 4)
2. For m repeated roots with (n - rn) distinct roots, write the re-
peated root to each power:
B{8) 8 - 8i
^
8 - Sn_m (« “ Sm)*"
Bj
-H 4. ...
(« - s^r-^ is - an.)
(A-8)
In order to solve for the A^s, multiply by the denominator factor and
set 8 equal to that particular root as in the previous case. In order to
solve for the B% a derivative is necessary,
A (a)
Bm = - -
(a ai)--*(s
^
ds (fi — ai)
• •
• (a - 8nr-^) Mm
Differential Analysers 337
Example A-t. Find the inversion for
A(s) *4- 3
B(8) a(«-|-2)(«-|-
D*
'
+ O2 ^-
«
« /
(i
• N*
+ 1)*^
'
(s-l- 1)
Ai = U + 3X I
3
(.^-h2)C.-h 1)*U 2
A, =
a + 3
I
1
.V 4- 3
- - -2
4 2)
B. =
d &“
43 -1
d« s(fi 4 2)
The desired expansion is
A is) 3 1 I
The inversion is
The
electronic analog computer employs feedback d-c amplifiers and
capacitors and resistor networks to accomplish mathematical manipula-
tion. These elements are shown in Fig. A-5. The d-c amplifier consists
Megohrnt Microfarads
6o
= -Keg (A-9)
The coefficient circuit is shown in Fig. A-6. All direct current passing
through resistor /i\ passes through resistor R; because the grid current
to the input of the amplifier is negligible. Then
ex - eg = (A- 10)
and
eg eg %R f (A>11)
(A-12)
If amplifier gain K is greater than about 1000, and if the resistance ratio
R//Ri is not greater than about 10, then the equation above reduces
to, with small error,
(A-I3)
Actually the amplifier gain is usually inuc\h higher than 1000, so that
the error is less than a small fi action of one per cent. As the last equation
show&, this circuit may be used for multiplying by a fixed coefficient.
The inUgraiing circuit is also shown in Fig. A-6. For the resistor and
capacitor
1 11 (A-14)
- «•) - * {A-16)
5
Eliminating current i and grid voltage €g by employing equation A-9,
there results
Again if amplifier gain is high enough we may write with little error,
Thus this circuit provides integration. The initial value of the output
is the constant A"' and is represented by the initial voltage charge across
the capacitor. T ntegration is always accompanied by a “scale coefficient”
of l/RC. By selecting values of the time constant RC a change of time
scale may be made. Scaling is discussed in a later section.
The adding circmi of Fig. A-'6 may be analyzed in a manner similar
to the coefficient and integrating circuits. With current ii in the
resistor Ri and current ?2 through resistor the equation becomes
Again if
K' I
the amplifier gain
i;)]
"
©
K is large and if all resistances are equal,
Co = -Cl - C2 (A-19)
Thus the summing circuit simply adds voltages while changing their
algebraic signs.
The combinations of coefficient, addition, and integration now allow
solution of ordinary linear differential equations with constant co-
efficients. A voltage input function corresponding to the desired dis-
turbance function is required. A means of indicating or recording direct
voltage is also necessary in order to display the solution. The solution
of equation proceeds as follows:
sf + *-»«)
Solve for the highest derivative:
first step
circuit. Then multiply by ^ using the coefficient circuit. This yields the
right-hand side of the equation. Close the loop by connecting the output
back to the input, thereby comjselling agreement. Note that there are three
amplifiers in the loop.
The disturbance voltage is introduced at y(t) o^nd may be a 20-volt change
in voltage representing a step change. The cutput (minus z) is recorded or
indicated and is the solution to the equation. As will be shown later this
circuit may be simplified.
Differential Analyzers 341
Example A-4. Given the simultaneous equations
* * *(0 +y
22/ -f 2y = a: +i
The equations are solved for the highest derivatives:
y = *(0 — X
—i » X — 2y — 2y
Next draw the each equation and interconnect the circuits. The
circuit!^ for
result is shownin Fig. A-8. Notice that the amplifier 6 in the lower row is
used merely to change signs. This is often necessary in computei problems.
Note that each of the three loops contains an odd number of amplifiers.
When the actual analog circuits are drawn out, a short study often
indicates that certain combination functions could be employed and
some amplifiers eliminated. For example, the circuit in tlieupper row
in Fig. A-9 could be simplified to one circuit as shown. The circuit
emplojdng fewer components is always more accurate.
Time scaling is the alteration of the time scale of the actual problem.
Consider a differential equation,
d^x dx
M-^ + B- + Kx^y(l) (A-20)
T = at (A-21)
„ d^x dx / t\
1. The actual problem solution may require so much time that speed-
ing up is necessary to save time.
2. The actual problem solution occurs so quickly that a slowing down
is necessary to ^^see^’ what is going on.
3. Better accuracy is obtained by operating a computer in the best
ranges of time events, neither too fast nor too slow.
the solution by the inverse of the integrator time constant RC. This is
handled most easily by labelling every point on the computer diagram
with the proper values of variables, including the scaling.
10
J+I = j/(t)
i
iy - x)
dr 10a
II
The circuit shown in Fig. A-10 is simplified from that of Fig. A-7. If the
computer is to speed up the problem by a factor of 60, then
T machine time !_
7 real time 60
Differential Analyzers 343
The capacitor and resistance values that may be selected are shown in the
tabl^ below.
Rt Ri c
1 megohm 1 megohm 10 mfd
5 megohms 5 megohms 2 mfd
10 megohms 10 megohms 1 mfd
All these combinations provide the same scaling and the selections of one or
the other is dictated by the particular kind of computer employed.
aHl~^ + aB'f-+KX=tiy{^
Thus the only effect of magnitude scaling is to alter the magnitude of the
disturbance signal y{t) when solving linear equations.
dh
X
dT*
dr a
Suppose that a time-scale factor of ^ is employed so that a = in order
so that
to speed up the problem, fjpppose also that 50 volts equals 10 in.
P - b volts per in. The disturbing function originally proceeded at the rate
of 1 in. per sec. The disturbing function must be a direct voltage changing
at the rate of
P I in . /OOX 5 volts
= SOOvolts/sec
a sec \1 / in.
344 Automatic Process Control
If this rate of change of voltage is too high for the available function generating
equipment then new values of time and magnitude scaling must be selected
until all parts of the problem satisfactorily fit the computing equipment.
»(«-»(;)
The input signals most commonly used in automatic control ar^:
These have two input points (accepting voltages) and one output
voltage.Whereas the product of two large numbers is even larger, it
is common to reduce the output by a factor, as
X
0 -x)^X + V
Valve Vessel
noscle or other restricting device. With the inflow set to a midvalue, and
after waiting a short time for steady-state conditions to be achieved, the
pressure or head in the vessel is constant, and owtflow equals inflow. A quick
charge of inflow of about 10 per cent magnitude will cause a pressure or head
change in the vessel. This change should be recorded on an instrument. This
is the transient response of the process. The time constant can be determiner!
from the experiment and compared to the calculated value.
:
Experiment 347
The following procedure is suggested:
^ ^
aq
(liquid) sec/ft*
dv
^ ^ sec/ft*
dq
V
C = (isothermal gas) ft*
Kq I
R
TsA' 1
^=
dv
proportional sensitivity
5.Explain why the upper and lower limits of velocity of output exist.
6.Sketch the block diagram of the apparatus and show the system function
in terms of its actual nu gibers.
" Tr.8 -h 1
TD»+lk
350 Automatic Process Control
This trsDsform esn he written as
TdEKc EKe
M(8)
T l8 1 8{T l8 -f- 1
)
and inverted
"•(0 - [i +
This curve is plotted in Fig. A-13 for Td = STl and Kc = 0.5. The ratio
of the maximum to the final value of m is the time constant ratio Td/Tl-
The 63.2 per cent time of the decay is Tl-
Td _ ynma» ~ m,
Tl
where wi, = the initial output before set point was changed.
4. From the record measure the lag time constant Tl and calculate Tji.
5. Repeat for several magnitudes of changes and in both directions.
Experiment 35I
6.Sketch the block diagram of the apparatuB and show the syitem function
in termB of its actual numbers (adjustments at mid value).
KcjX -h To 8 )
X -X- T 1,8
z Manual It
Process
control
h'
Input
element
Controller
kH
Valve Process
hO
Measuring
element
record the transient response of the process under control. The automatic
controller should have manual control means by which manual control of the
process (y~^x connection) can be obtained. Automatic control is achieved
by switching in the controller {p—* x connection).
The open-loop response, b output to x input, should be second-order, other-
wise loop gain is not easily calculated. 'J'he response is generally
R
- (7\s+l)(Tts+l)''
RRc + Tj
if* \T2 Tj
EXPERIMENT 6. PROPORTIONAl^lN1'EGR.\L-DERIVATlVE
CONTROL OF A PROCESS
Preparation: Chapters 2 and 3, or 6 and 4
The object to determine the “best" setting of a
of this experiment is
A'c * 0.6K„
and
L » 0.25/"«
The settings of integral and derivative time may now he determined by tlie
Ziegler-Nichols equations
The procedure is as follows •
1. delect operating set point v for system and check to see that all x anablcs
p, fit, c,and b are in middle of respective ranges.
2. Turn oflF integral effect and derivative effect.
Gradually increase proportional sensitivity so tlu»i a stead^> (>' '^’.dioD
3.
begins. Make sure that none of the variables p, m, c, or h is allowed to limit
during this operation. Move set point to see that marginal stability is definite
and sustained. Note proportional sensitivity and period of oscillation Pu.
Set proportional sensitivity to 60 per cent of the value causing marginal
4.
stability. Check now by moving the set point slightly to see that the transient
response is satisfactory one-ninth ratio of succeeding amplitudes is suggested.
;
T, = 2.fiL
K,' = ; A*
3
A/ = 0.4L
T,' = 2L
354 Automatic Process Control
With the selected settings check the stability by moving the set point slightly
to see that the transient repsonse is satisfactory. Trim the settings if neces-
sary» maintaining a one-ninth ratio of amplitudes, good return, and small
period of oscillation.
7. Check the settings obtained above by emplo3rmg the reaction-curve
method.
8. Compare the various values of parameters /Ce, T,, and Td obtained in
the tests and discuss any significant differences.
where RiCi is the valve time constant (usually 1 to 10 sec). The first process
time constant is represented at amplifier 6,
m R%C%s -f- 1
so that the (iroccss time (constant is set by RtCt. The second process time
constant is represented at amplifier 7,
c 1
h R^C^8 + 1
“ = (52(1 — 52)/?2C'2K + 1
^
- - +
1
^ SiRiCiS
80 that an adjustment S 4 of potentiometer
1 R4 will vary the integral time from
0 to R 4C 4 .
gated. one of the time constants is zero, simply remove the corresponding
If
capacitor in the computer circuit.
2 Time scale the computer circuit. The time scaling will depend upon the
.
25 per cent. Try several related \ alues of derivative time and loop gain until
a transient response with good stability and minimum ])eriod of oscillation is
obtained.
5. Turn off the derivative action and decrease loop gain by 35 per cent from
the last step (90 per cent of original value). Set the integral time to give a
good return on a load change.
6 Turn in derivative time to 30 per cent more than the value set in step 4.
.
Set integral time to double that value found in step 5. Set proportional
sensitivity to about 20 per cent higher than in step 5. Try several values near
these settings to see if a better response can be obtained. Use the highest
derivative time possible because this will improve loop gain and provide
smaller period of oscillation.
7. Leave all settings at their l>est values and turn out derivative action to
show its effect. Next turn out integral action and with a load change at Mi or
u? determine whether the offset is
1
"
1 -I- loop gain
as a fraction.
S. Eliminate the control-element lag by removing capacitor Ct. Find new
control settings to det(‘rmine whether there is a significant improvement when
the control-element lag is removed.
Experiment 357
1. Move the set point with the measured (feedback) variable maintained
2A
tanh^*
R K
The apparent process time constant T may now be calculated.
Estimate the amplitude and period of oscillation if differential gap 0
is zero. The following equations are used.
A = 1(1
P = 4Litonh->^
The process should be a fluid process (as in Fig. A-12) with a nonlinear
resistance at the outlet, the usual case. The controller should be set for two
responses such as proportional-derivative. Although the industrial problem
usually involves control for various load changes,it is easier in the laboratory
1 . Set the controller in operation and check all operating values of variables
to insure correct operation.
2. Beginning with the set point at a midvalue (50 per cent) determine the
best loop gain and integral or derivative time to give a certain transient
response. Record the value of flow through the control valve required to
4.
maintain the set point, and record the value of controller output x at each of
these points, l^te these values.
3. Repeat o()eration 2 for as many set-point values as are desired, at least
5. The loop gain
20, 35, 60, 65, 80 per cent. is
^
de dx
dm ^
dm
where ^ controller sensitivity = Kc
at
dm —
valve sensitivity
dx
dc
process sensitivity == Kp
dm
The set point (or controlled variable) should be plotted against valve flow.
The slope of this curve is Kp.
The required valve characteristic may now be synthesised by calculating
dm K
dx dc
^
dm (f?
dm K
Experiment 359
Uae the values of K and K„ and calculate the slope of the curve K, at each
value of flow m.
6. By cut and try, parabolic curve fitting and integration, or other methoda,
find the curve m » fix). Compare this curve to the actual plot of m veraua x
for the valve in use, and indicate whether it has the desired charaeteriatic.
7. If the prorcFS is of second order and the constants are known, uae the
analytical method given in Chapter 8 to verify these results.
Answers to Problems
1-2. Furnace thermostat, hot water heater control, water-preseure regulator, toilet
2wKL A
2-22. [T’jjrjjs* -H (I'll ri- — (rijs -f l)Asut
2-23. [RCRiCiS^ + (i2C -b RiC + ffiCOa -}- Ijc = Rim
2-24. {RCn -f l)c - -RK^m -H i?u
200
Aaswcrs to Problems 3«1
a-ts. (Ms* + Bs + Jf)i -/ + K'v + Bs»
a-27. f-Btx
a-28. mt(R + fit + JUiiCt) m mi(fii 4- RHtCt)
M9. [(Bi + A)Cs 4- l]mt » (RCt + l)mi
a-W. ZTi 4- Dm, - m,
(T*s* 4-
(B/Ay)» 4- 11c “ m/y, when y
ICB/ff)** 4- fluid deniity, K “ tube »m
»-aa. -(Ms* 4- Bs 4- K)e ~ Am, K - WRA^T/V^
V - volume under piston, R > gns constant, T ~ trmperatun
a-aa. (Ts 4-
3- D® - m, t - pDC/nr
p “ density, B =diameter, C = up. hi., B => oonv. coeff.
4-
J-i. riTj/Ti 4- 7 ’j: ?V(7’i 4- 7’,); (Ti 4- Tt)
3-5. Prop. K\ integral KT
3-5. Prop. T\IT\, integral T2
3-7. Prop. (Ti + rO/Ta; integral (Ti + T’l); derivative T^T^fiTi + P*)
3-8. Prop. JiC; derivative (Pi + Ti)/K
3-9. Prop. T\IT 2 \
derivative P2
10. Prop. 1; derivative 2P; eecond derivative P*.
6- 1 . he ^ Kv -\r (/>
—
6-.1 .
r, ,
AKx
6-5. (Pfi + Dm - A>, T
K{Ki -HR,)’ “
RCR, + Ki)
6-6. * R/AC where R spring gradient, A - diaphragm area, C
nozzle coefficient 9w/dx
K,
M. Withi;,
1 + K^Ki + KiKiKi
- Kz
Without
1 4 KxKzKi
ft-6. (a)Square wave as shown. (&) Triangular wave.
9r9. Zero
MO. Ti * 4KiA/Ki^, A vessel area, Ki *= dq/dx valve, Kz * dqfdc valve
^
-11 . K-K-
0.11
8
1 1 4 -
~ tan-‘ 26,r
— 4wZ«
0^. —uL and tan *
4 —« L*
M. <1 time constant element
RK^ a: 1 capacitance element
RKt > 1 exothermic
Answers to Problems 363
w. •V. 0.29 sec; Probably not.
60%.
10-3. /Cf *