Mansol Dalpro 1-14
Mansol Dalpro 1-14
Parameters : ρ, V
a) Mass Balances:
Variables : w1, w2, T1, T2, T3
NE = 1 dh1
2.1 ρA1 = w1 − w2 − w3 (1)
NV = 5 dt
Thus, NF = 5 – 1 = 4 dh
a) Overall mass balance: ρA2 2 = w2 (2)
dt
Because w1, w2, T1 and T2 are determined by upstream units, we assume
d (ρV )
= w1 + w2 − w3 (1) they are known functions of time: Flow relations:
dt
w1 = w1(t) Let P1 be the pressure at the bottom of tank 1.
Energy balance: w2 = w2 (t)
T1 = T1(t) Let P2 be the pressure at the bottom of tank 2.
d ρV (T3 − Tref ) T2 = T2(t)
C = w1C (T1 − Tref ) + w2C (T2 − Tref )
dt (2) Let Pa be the ambient pressure.
− w3C (T3 − Tref )
Thus, NF is reduced to 0.
P1 − P2 ρg
Then w2 = = (h1 − h2 ) (3)
Because ρ = constant and V = V = constant, Eq. 1 becomes: R2 g c R2
w3 = w1 + w2 (3) 2.2
P1 − Pa ρg
w3 = = h1 (4)
R3 g c R3
b) From Eq. 2, substituting Eq. 3 Energy balance:
d (T3 − Tref ) b) Seven parameters: ρ, A1, A2, g, gc, R2, R3
dT d ρV (T − Tref )
ρCV = ρCV 3= w1C (T1 − Tref ) + w2C (T2 − Tref )
dt dt (4) Cp = wC p (Ti − Tref ) − wC p (T − Tref ) − UAs (T − Ta ) + Q Five variables : h1, h2, w1, w2, w3
dt
− ( w1 + w2 ) C (T3 − Tref )
Simplifying Four equations
Constants C and Tref can be cancelled:
dT Thus NF = 5 – 4 = 1
ρVC p = wC p Ti − wC p T − UAs (T − Ta ) + Q
dT3 dt
ρV = w1T1 + w2T2 − ( w1 + w2 )T3 (5) 1 input = w1 (specified function of time)
dt dT
ρVC p = wC p (Ti − T ) − UAs (T − Ta ) + Q 4 outputs = h1, h2, w2, w3
dt
The simplified model now consists only of Eq. 5.
1/ γ
2.4 2.5 V2 γ
1 dP2
P2
(1− γ ) / γ
= wb − wc
γ C dt
Assume constant liquid density, ρ . The mass balance for the tank is a) For linear valve flow characteristics,
as the new dynamic model. If the ideal gas law were not valid, one would
d (ρAh + m g ) Pd − P1 P − P2 P2 − Pf use an appropriate equation of state instead of Eq. 3.
= ρ(qi − q ) wa = , wb = 1 , wc = (1)
dt Ra Rb Rc
Mass balances for the surge tanks
Because ρ, A, and mg are constant, this equation becomes
dm1 dm2
dh = wa − wb , = wb − wc (2)
A = qi − q (1) dt dt
dt 2.6
where m1 and m2 are the masses of gas in surge tanks 1 and 2,
The square-root relationship for flow through the control valve is respectively.
a) Assumptions:
1/ 2
ρgh If the ideal gas law holds, then
q = C v Pg + − Pa (2) 1. Each compartment is perfectly mixed.
gc
P1V1 =
m1
RT1 , P2V2 =
m2
RT2 (3) 2. ρ and C are constant.
M M 3. No heat losses to ambient.
From the ideal gas law,
where M is the molecular weight of the gas Compartment 1:
(m g / M ) RT T1 and T2 are the temperatures in the surge tanks.
Pg = (3) Overall balance (No accumulation of mass):
A( H − h)
Substituting for m1 and m2 from Eq. 3 into Eq. 2, and noticing that V1, T1, 0 = ρq − ρq1 thus q1 = q (1)
where T is the absolute temperature of the gas. V2, and T2 are constant,
Energy balance (No change in volume):
Equation 1 gives the unsteady-state model upon substitution of q from Eq. 2 and V1M dP1 V2 M dP2
= wa − wb and = wb − wc (4)
of Pg from Eq. 3: RT1 dt RT2 dt dT1
V1ρC = ρqC (Ti − T1 ) − UA(T1 − T2 ) (2)
dt
1/ 2
dh (mg / M ) RT ρ gh The dynamic model consists of Eqs. 1 and 4.
A = qi − Cv + − Pa (4) Compartment 2:
dt A( H − h) gc b) For adiabatic operation, Eq. 3 is replaced by
Overall balance:
Because the model contains Pa, operation of the system is not independent of Pa. γ γ
V V
For an open system Pg = Pa and Eq. 2 shows that the system is independent of Pa. P1 1 = P2 2 = C , a constant (5) 0 = ρq1 − ρq2 thus q2 = q1= q (3)
m1 m2
1/ γ 1/ γ
PV γ PV γ Energy balance:
or m1 = 1 1 and m2 = 2 2 (6)
C C dT2
V2ρC = ρqC (T1 − T2 ) + UA(T1 − T2 ) − U c Ac (T2 − Tc ) (4)
dt
Substituting Eq. 6 into Eq. 2 gives,
1/ γ
b) Eight parameters: ρ, V1, V2, C, U, A, Uc, Ac
1 V1 γ dP1 Five variables: Ti, T1, T2, q, Tc
P1
(1− γ ) / γ
= wa − wb
γ C dt
Two equations: (2) and (4)
2 outputs = T1, T2
(b) In this part, two controllers have been added to the system. Each controller
3 inputs = Ti, Tc, q (specify as functions of t) provides an additional equation. Also, the flow out of the tank is now a
manipulated variable being adjusted by the controller. So, we have 2.9
c) Three new variables: ci, c1, c2 (concentration of species A). 4 parameters: C , ρ , Tsp , Vsp
Two new equations: Component material balances on each compartment.
6 variables: V , T , wi , Ti , Q, w Assume that the feed contains only A and B, and no C. Component balances for
c1 and c2 are new outputs. ci must be a known function of time.
4 equations A, B, C over the reactor give.
2.10 2.11
Note that the only conservation equation required to find h is an overall Rearrange Eq. 5 to get the feedforward (FF) controller relation,
mass balance:
dx3 1
= [120(0.6 − x3 ) + 100(0.5 − x3 )]
dt (800)(1.75)π
d) In this case, both h and x3 will be changing functions of time. Therefore, d) At steady state, the dynamic model given by Eqs. 1, 2 and 3 becomes:
both Eqs. 1 and 3 will have to be solved simultaneously. Since
concentration does not appear in Eq. 1, we would anticipate no effect on 0 = rg - DX DX = rg (5)
the h response.
0 = rp - DP DP = rp (6)
1
0 = 𝐷𝐷�𝑆𝑆𝑓𝑓 − 𝑆𝑆� − 𝑌𝑌 𝑟𝑟𝑔𝑔 (7)
2.14 𝑋𝑋/𝑆𝑆
dS 1
Substrate: V = F (S f − S ) − Vrg Substituting Eq. 9 into Eq. 8,
dt YX / S
or DX = Y X / S ( S f − S ) D (10)
dS F 1 1
= ( S f − S ) − rg − rP (3)
dt V YX / S YP / S
From Eq. 4
b) At steady state, DK S
S=
mmax − D
r kcs t kc s dS SX 1 0.2
∫ r dr =
o
−
ρ ∫0
dt ∴ r (t ) = ro −
r
t (4)
dt
= 0.2 − −
1 + S 0.5 0.1
By using MATLAB, this system of differential equations can be solved. The time
Finally,
to achieve a 90% conversion of S is t = 22.15 h.
M = rV = rπhr 2
then (b) the drain is opened for 15 mins; assume a time constant in a linear
2
kc transfer function of 3 mins, so a steady state is essentially reached. (3 < t < 18).
M (t ) = ρπh ro − s t Assume that the process will return to its previous steady state in an exponential
ρ
manner, reaching 63.2% of the response in three minutes.
b) The time required for the pill radius r to be reduced by 90% is given by
Eq. 4:
kc s 0.9ρo ρ (0.9)(0.4)(1.2)
0.1ro = ro − t ∴ t= = = 54 min
r kc s (0.016)(0.5)
Therefore, t = 54 min .
2.18
The graph should show an exponential decrease to the previous steady state of 0.5
ft. The initial value should coincide with the final value from part (c).
(c) the inflow rate is doubled for 6 minutes (18 < t < 24) Degrees of freedom = (11-6) (number of variables) – 2 (number of equations) = 3
The height should rise exponentially towards a new steady state value double that
of the steady state value in part b), but it should be apparent that the height does
not reach this new steady state value at t = 24 min.. The new steady state would The three input variables (w, Ti, Q) are assigned and the resulting system has zero
be 1 ft. degrees of freedom.
(d) the inflow rate is returned to its original value for 16 minutes (24 < t < 40)
(a) First we simulate a step change in the Fuel Gas Purity (FG_pur) from 1 to
(a) First we simulate a step change in the vapor flow rate from 0.033 to 0.045 0.95.
m3/s. The resulting plots of Oxygen Exit Concentration (C_O2) and Hydrocarbon
The resulting plots of xD and xB are shown below. Outlet Temperature (T_HC) are shown below.
Figure: Plot of xD, xB, and zF versus time for a step change in zF from 0.5 to 0.55
By examining the resulting data, we can find the steady-state values of xD and xB
before and after the step change in zF.
By examining the resulting data, we can find the steady-state values of xD and xB xB 0.0050 0.0068 +0.0018 By examining the resulting data, we can find the steady-state values of C_O2 and
before and after the step change in V. T_HC before and after the step change in FG_pur.
(c) Increasing V causes xD and xB to decrease, while increasing zF causes both
Start End Change xD and xB to increase. The magnitude of the effect is greater for changing V Start End Change
than for changing zF. When changing V, xB changes more quickly than xD.
xD 0.85 0.73 -0.12 C_O2 0.92 1.06 0.14
xB 0.15 0.0050 -0.145 T_HC 609 595 -14
(b) Next we simulate a step change in the feed composition (zF) from 0.5 to 0.55. (b) Next we simulate a step change in the Hydrocarbon Flow Rate (F_HC_sp)
Note that the vapor flow rate, V, is still set at 0.045 m3/s. from 0.035 to 0.0385. Note that the Fuel Gas Purity, FG_pur, is still set at
0.95.
( 4 ) =
1
2.21 τ AR
= = 2 min
2
Part a
The key to this problem is solving the mass balance of the tank in each part.
dh
Mass balance: A = qi − qo (Mass Balance)
dt
d
( ρ Ah
= ) ρ qi − ρ qo dh
dt 4 =2 (Separable ODE)
dt
- ρ (density) and A (tank cross-sectional area) are constants, therefore:
1
dh
∫ dh = ∫ 2 dt
A = qi − qo
dt
1
h(t ) = t +C h(0) =
1
- The problem specifies qo is linearly related to the tank height 2
1
qo =
1
h h(t ) = t +1 0≤t <3
R 2
dh 1
A = qi − h
Figure: Plot of C_O2, T_HC, and F_HC_sp versus time for a step change in F_HC_sp dt R
from 0.035 to 0.0385.
By examining the resulting data, we can find the steady-state values of C_O2 and - Next, we can obtain R (valve constant) from the steady state information in the
T_HC before and after the step change in F_HC_sp. problem
dh
Start End Change = 0 at steady state
dt
C_O2 1.06 1.06 0
1
T_HC 595 572 -23 0= qi − h
R
dh 1
+ h=
1 (Solution by integrating factor) Where Ci is the proportionality constant for each tank.
dt 2
∫ d (e h )
= ∫ 1e dt h(18) 1
=
t/2 t/2
Results:
- Method is same as part b. a) The final equations for the height of liquid in each tank are shown below.
γ
dh1
dt
=
C C
− 1 h1 + 3 h3 + 1 F1
A1 A1 A1
(1) Chapter 3 © 3.2
dh2 C C γ
=
− 2 h2 + 4 h4 + 2 F2 (2) Break the pulse into three step functions. First, a step up to 10 at t=0. Then, a step
dt A2 A2 A2 down by 8 at t=1. Finally, a step down by 2 at t=3:
dh3 C (1 − γ 2 ) 3.1
=
− 3 h3 + F2 (3) f(t) = 10 S(t) – 8 S(t-1) – 2 S(t- 3)
dt A3 A3
dh4 C ( 1 − γ1 ) (a)
dt
=
− 4 h4 +
A4 A4
F1 (4)
f (t ) 5 e 3t
te 4t F (s ) =
1
s
10 - 8e-s - 2 e-3s
Transform each term using rules 2, 5, and 7 from Table 3.1, respectively.
b) Now we can substitute γ= γ= 0.5 5 1 1
1 2
F ( s)
s s 3 ( s 4)2 3.3
dh1 C C 0.5
=
− 1 h1 + 3 h3 + F1
dt A1 A1 A1 (b)
5 a) Pulse width is obtained when x(t) = 0. Since x(t) = h – at
dh2 C C 0.5 f (t ) sin(4t ) (t 3) S (t 3) e (t 3) S (t 3)
=
− 2 h2 + 4 h4 + F2 t
dt A2 A2 A2 t : h at = 0 or t = h/a
To transform sin(4t), use rule 14 from Table 3.1 b)
dh3 C 0.5 To transform (t-3)S(t-3) use rules 3 and 26 together. To use rule 26,
=
− 3 h3 + F2
dt A3 A3 set f(t) = t and t0 =3.
dh4 C 0.5 To transform e-(t-3)S(t-3) use rules 5 and 26 together. To use rule h
=
− 4 h4 + F1 26, set f(t) = e-t and t0 =3. slope = -a
x(t)
slope = a
dt A4 A4 x(t)
Note that there is no Laplace transform for 1/t.
The differential equations for the tank heights are coupled, so the heights 4 e 3 s e 3 s 1
F (s) 2 5
cannot be solved for or controlled independently. F1 and F2 can be used to s 16 s 2 s 1 t
slope = -a
s
In Eq. 1, substitute any s-1 to determine 1. Arbitrarily using s=0, Eq. 1 Then
a3
gives
3c a a a4 1 12 3 / 4 2 3
Y ( s) 1 2 Xˆ ( s )
s ( s 2)( s 4) s s 2 ( s 2) ( s 4)
2
2 1 1 s s4 s3
or 1 1
12 1 12
Regardless of the numerical values of a1, a2, a3, and a4, the inverse Laplace 1 3 4t 2 3t
transform indicates y(t) includes a1 , t , e2t , and e4t . xˆ(t ) e e
1 1 12 4 3
X ( s) and x( t ) et tet
s 1 ( s 1)2
1 1 1
3.5 c) X ( s) 2 Using the Real Translation Theorem,
s s 1 1 2 3 s b 2 2
s
2 4 1 3 4(t 0.5) 2 3(t 0.5)
55 55 x(t ) xˆ (t 0.5) e e
T(t) = 20 S(t) + t S(t) – (t-30) S(t-30) 1 3 12 4 3
30 30 where b and
2 2
T ( s)
20 55 1
s 30 s 2
55 1 30s 20 55 1
30 s 2
e
s 30 s 2
1 e 30s t
for t 0.5
1 bt 2 2 3
x(t ) e sin t e sin t
3 2
3.6 3.7
s 1
d) X(s) = e0.5 s
s( s 1) s( s 4)( s 3) 6( s 1) 6 1 2
a) X ( s) 1 2 3 Y ( s) 2
( s 2)( s 3)( s 4) s 2 s 3 s 4 a) s 2 ( s 1) s 2 s s
To invert, first ignore the time delay term. Using the Heaviside expansion
s( s 1) with the partial fraction expansion, 6
1 1 2 s2 6 1 0
( s 3)( s 4) s2 s 0
s 2 s 1 A B C
Xˆ ( s) 6
s( s 4)( s 3) s s 4 s 3 Y (s) 2
s( s 1) s
2 6
( s 2)( s 4) s 3 Multiply by s and let s 0
12( s 2) s 3 1 s 2 3s 4 2
b) Y ( s ) 1 22 = 5 or (s2 + 4s + 3) X(s) =
s ( s 2 9) s s 9 s 2 2 s 2 ( s 2 2 s 2) 2 s 2 s( s 1)
2
X(s) =
Multiplying both sides by s(s2+9) Multiplying both sides by (s 2 2s 2) 2 (s 2) gives s( s 1)2 ( s 3)
12( s 2) 1 ( s 9) ( 2 s 3 )( s)
2
or 1 = 1s4 + 41s3 + 61s2 +41s + 2s3 +42s2 +62s +42 + 3s2 +23s + Performing partial fraction expansion and taking the inverse Laplace
12s 24 (1 2 ) s 2 3 s 91 4s + 24 + 5s4 + 45s3 + 85s2 + 85s + 45 transform (either manually or using a symbolic software program), we get:
s1: 3 = 12 b) Applying the Final Value Theorem (note that the theorem is applicable here)
s0: 91 = 24 s3 : 41 + 2 + 45 = 0
2 2
s2 : 61 + 42 + 3 + 85 = 0 lim x(t ) lim sX ( s) lim
Solving simultaneously, t s 0 s 0 ( s 1)2 ( s 3) 3
( s 2)( s 3) 2 2 A B C
2 6 iii) Y ( s)
( s 4)( s 6) s 5
a) From Eq. 3-66 s( s 2 4s 4) s( s 2) 2 s ( s 2) 2 s 2
t 1
L f (t * )dt * F ( s)
( s 2)( s 3) 0 s y(t) will contain terms of form: constant, e-2t , te-2t
3 6
( s 4)( s 5) s 6 t 1 1
we know that L e d = L e t iv) Y (s)
2
0 s s( s 1) s( s 2 4s 8)
1 6 6
Y (s)
s4 s5 s6
Laplace transforming yields s 2 4s 8 (s 2 4s 4) (8 4) (s 2) 2 2 2
1 1 2
Y ( s) Y ( s)
(s 1)
d) 2
2 2
1 ( s 2) ( s 2s 2) 2 ( s 2)
2 2
s X(s) + 4X(s) + 3X(s) = s[( s 2) 2 2 2 ]
s( s 1)
y(t) will contain terms of form: constant, e-2t sin2t, e-2tcos2t
2( s 1) 1
A = lim 2 1
s 3 X (s) X (s)
500
s 0 ( s 4) 2 c)
s 1
2
400
1 1
2
2(s+1) = A(s +4) + Bs(s) + Cs X (s) 2
x(t)
( s 1)( s 3 1) 1 3 1 3
( s j )( s j )( s 1)( s j )( s
300
j)
2s+2 = As2 + 4A + Bs2 + Cs 2 2 2 2
200
Equating coefficients on like powers of s The denominator contains complex factors; x(t) is oscillatory. The 100
s1 : 4 time
2= C C=2 d) s 2 X ( s) sX ( s)
1 s Figure S3.10b. Simulation of X(s) for case b)
s: 0
2 = 4A A= 4 4
2 X (s)
s ( s 2 s ) s 2 ( s 1)
80
1 2 (1 2) s 2
Y(s) 2
s s 22 s 2 22 60
The denominator of [sX(s)] contains no complex factors; x(t) is not
1 1 2 oscillatory. The denominator of [sX(s)] vanishes at s = 0; x(t) is not
y(t) = cos 2t sin 2t convergent. See Fig. S3.10d.
40
2 2 2
x(t)
3.5 20
1
y(t) = (1 cos 2t ) sin 2t 3
2 0
2.5
-20
2
3.10
x(t)
1.5 -40
0 1 2 3 4 5 6 7 8 9 10
time
3
s 3 X ( s ) 2s 2 X ( s ) 2sX ( s ) X ( s )
s 0
3 3
X (s) -0.5
0 1 2 3 4 5 6 7 8 9 10
s ( s 3 2s 2 2s 1) 1 3 1 3 time
s ( s 1)( s j )( s j)
2 2 2 2 Figure S3.10a. Simulation of X(s) for case a)
The denominator of [sX(s)] contains complex factors so that x(t) is
oscillatory, and the denominator vanishes at real values of s= 1 and -½
which are all <0; thus x(t) is converges. See Fig. S3.10a.
2
b) s 2 X ( s) X ( s)
s 1
2 2
X ( s)
( s 1)( s 2 1) ( s 1) 2 (s 1)
b 1 1
18 b) If u(t) = be-t/ then U(s) = 1
s 1 ( s 3 4) 5
16 s 1
1 1
14 Kb 2
X b ( s) 1 2 ( s 1)( s 0.79 1.37 j )( s 0.79 1.37 j ) 19.6
12 (s 1)( 1 s 1)( 2 s 1) s 1.59
1
10 3 j 3 0.74 0.59 j
xb(t) = fb(e-t/ , e-t/1, e– t/2) ( s 1)( s 1.59)( s 0.79 1.37 j )
x(t)
s 0.791.37 j
8
c
If u(t) =ce-t/ where = 1 , then U(s) =
6
c) 1 1
1 s 1
0.074 0.059 j 0.074 0.059 j
X(s) 5 19.6
4
Kd 3
X d ( s) sX (s) 12 X(s)
Since the time function in the solution is not a function of initial conditions, take ( s 2 2 )( 1 s 1)( 2 s 1) s2 9
the Laplace transform with:
xd(t) = fd(e– t/1, e– t/2, sin t, cos t)
time
3 3
dx(0) X ( s)
x(0) 0 ( s 2 9)( s 12) ( s 3 j )( s 3 j )( s 12)
dt
3.12 1 j1 1 j1 3
12s2X(s) + (1+2)sX(s) + X(s) = KU(s)
s 3j s 3j s 12
K dx3 d 2 x(0) dx(0)
X ( s) U (s) a) 4 x et with x(0) 0 3 3 1 4
1 2 s 2 (1 2 ) s 1 dt 3 dt 2 dt 1 j1 j
( s 3 j )( s 12) s 3 j
18 72 j 102 102
Factoring the denominator Take the Laplace transform of the equation:
3 1
K 1 3
X (s) U ( s) s 3 X(s) 4 X(s) ( s 2 9 ) s 12 51
(1 s 1)( 2 s 1) s 1
1 4 1 4 1
a 1 1 j j
a) If u(t) = a S(t) then U(s)= X ( s) X (s) 102 102 102 102 51
s ( s 1)( s 3 4) ( s 1)( s 1.59)( s 0.79 1.37 j )( s 0.79 1.37 j ) s3j s 3j s 12
Ka 1 2 3 j 3 3 j 3
X a (s) 1 2 x(t )
1 1
(cos 3t 4 sin 3t ) e12 t
s(1 s 1)( 2 s 1) s 1 s 1.59 s 0.79 1.37 j s 0.79 1.37 j 51 51
20 20 40
Figure S3.13. Output variable for a step change in x of magnitude 1.5 3.15
3.13
3.14 f (t ) hS (t ) hS (t 1 / h)
4 x hS (t ) S (t 1 / h)
dx
d 2 y(t ) dy(t ) d ( x 2) , x(0)=0
3 y(t ) 4 x(t 2) dt
dt 2 dt dt First, take the Laplace transform of each term in the equation
Take the Laplace transform assuming zero initial conditions: Take the Laplace transform,
d2y 2 1 e s / h
2 s Y ( s ) sy (0) y '(0) s Y ( s ) s
2
sX ( s ) 4 X ( s ) h
s2Y(s) + 3sY(s) + Y(s) = 4 e-2ssX(s) e-2sX(s)
dt s s
Rearranging, dy 1
5 5( sY ( s ) y (0)) 5sY ( s ) 5 X ( s) h(1 e s / h ) h(1 e s / h ) 1 2
dt s( s 4) s s 4
Y (s) (1 4 s )e 2 s 6 y 6Y ( s) 1 1 1 1
G (s) 2 (1) 1 , 2
X (s) s 3s 1 s 4 s 0 4 s s 4 4
7
7
s
a) The standard form of the denominator is : 2s2 + 2s + 1 h 1 1
X ( s) (1 e s / h )
The final transformed equation is: 4 s s 4
From (1) , = 1 , = 1.5
h 1 e s / h 1 e s / h
Thus the system will exhibit overdamped and non-oscillatory responses.
4 s s s 4 s 4
b) Steady-state gain
s 3 8s 2 2 s 8
0 t <0 Y(s) =
( s 3) 2 ( s 2 1) 1 / 8 1 / 8( s 2) 3/ 4
Y ( s)
x(t )
h
(1 e 4t ) 0 < t < 1/h s ( s 2) 2 2 2 ( s 2) 2 2 2
4 To find y(t) we have to expand Y(s) into its partial fractions
4
h 4 ( t 1 / h )
e e 4t t > 1/h 1 1 3
y(t) = e-2t cos 2t + e-2t sin 2t
Y ( s)
A
B
Cs
D 8 8 8
( s 3) 2 s 3 s 2 1 s 2 1
1
h=1
y(t) = Ate-3t + Be-3t + C cost + D sin t 3.17
h=10
0.9 h=100
s 1
0.8
b) Y(s) = Laplace transform of the system of ODEs gives:
0.7
s( s 2 4s 8)
dy1
y2 e
0.6
t
2
4
Since 8 , there are complex factors. dt
x(t)
0.5
4
dy2
3 y2 2 y1
0.4
0.3
complete the square in denominator dt
0.2
s2 + 4s + 8 = s2 + 4s + 4 + 84 1
0.1
sY1 Y2 (1)
0
s 1
= s2 + 4s + 4 + 4 = (s+2)2 + (2)2 { b = 2 , =2}
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
time sY2 3Y2 2Y1 (2)
Figure S3.15. Solution for values h= 1, 10 and 100
Partial fraction expansion gives Next solve Equation 2 for Y2 in terms of Y1
3.16 A B( s 2) C s 1 Y2 ( s 3) 2Y1
Y ( s)
s s 2 4s 8 s 2 4s 8 s( s 2 4s 8) 2Y
Y2 1 (3)
a) Take the Laplace transform: s3
Multiply by s and let s0 Substitute equation 3 into equation 1 and solve for Y1
Multiply by s(s2+4s+8) s 3 s 1
s 2 1
(s + 6s + 9)Y(s) s(1) 2 –(6)(1)= 2
2
Y1 s
s 1 A(s2+4s+8) + B(s+2)s + Cs = s + 1 s 3 s 1
1
s As2 + 4As + 8A + Bs2 + 2Bs + Cs = s + 1 Y1
(s2 + 6s + 9)Y(s) = 2 2
+s+8 ( s 1)( s )
s 1 s3
1
s2 : A+B=0 B = A = Expand using partial fractions:
s s 3 s 8s 2 8 8
(s2 + 6s + 9)Y(s) = 1 1 3
s2 1 1
C = 1 + 2 4 = s3 1 1 2
s: 4A + 2B + C = 1 Y1
8 8 4 ( s 1)2 ( s 2) s 2 s 1 ( s 1)2
1
s0 : 8A = 1 A= (This checks with above result)
8
41
2Y1 2 2 2 2 Plot[g[t],{t,0,100},AxesLabel{time,Y},PlotRange{{0,100},{0 ,2}}]
Y2 Using MATLAB, the concentration response is shown in Fig. S3.18. (Consider
s 3 ( s 1)2 ( s 2) s 2 s 1 ( s 1) 2 V = 2 m3, ci=50 Kg/m3 and q = 0.4 m3/min) y
4
Finally, get both time-domain solutions using the inverse Laplace transform:
50
y1 (t ) e 2t e t 2te t 45
y2 (t ) 2(e 2t e t te t ) 40
3
35
30
3.18 2
c(t)
25
dc
V qc qci 20
dt
15 1
5
sVC(s) + qC(s) = q Ci(s)
0
0 t
0 5 10 15 20 25 30 0 2 4 6 8 10
Note that c(t = 0) = 0 Time
Figure S3.19a: Tank level response to a unit step change in flow rate.
Figure S3.18. Concentration response of the reactor effluent stream.
Also, ci(t) = 0 , t 0 (b) Define the time when y(t) reaches its maximum as tmax. This time occurs
ci(t) = c i , t>0 when y′(t)=0. Solve for this time using Mathematica and find that tmax=0.877
3.19 and y(tmax)=1.558. Therefore, the tank will not overflow.
Taking Laplace transform of the input function, a constant, gives
(c) Now find the general solution for any input step size, M (the solution is
ci (a) Take the Laplace transform of each term, taking into account that all initial denoted in this case as YM(s) and yM(t) for clarity). The input U(s) is M/s.
Ci ( s ) conditions are zero:
s
M
so that U (s)
ci qci s 2Y sy (0) y '(0) 5sY 5 y (0) Y 8sU 8u (0) U s
sVC(s) + qC(s) = q or C(s) =
s 2Y 5sY Y U (8s 1) M (8s 1)
s ( sV q)s YM ( s 5s 1)
2
1 s
U ( s) M (8s 1)
Dividing numerator and denominator by q s YM MY
8s 1 s ( s 2 5s 1)
ci Y ( s 2 5s 1)
C(s) = s
V 8s 1
YM is the previous Y, multiplied by the size of the step, M. Since M is a constant,
q s 1 s Y taking the inverse Laplace transform gives:
s ( s 2 5s 1)
yM (t ) My(t )
Use Transform pair #3 in Table 3.1 to invert ( =V/q) Now use symbolic mathematical software (ex. Mathematica) to solve for y(t).
Now solve the equation:
InverseLaplaceTransform[(8*s+1)/(s*(s^2+5*s+1)),s,t]
where ci(t) is a pulse of magnitude A and width tw. A pulse can be described by the
yM (tmax ) 2.5 My (tmax ) M (1.558) b) Taking Laplace transform of above equations and eliminating sum of two step functions. The first will be a step function of magnitude A at time
2.5 C1 (s) and C2 (s) gives 0. The second will be a step function of –A at t=tw.
M 1.605
1.558 q q q
ci (t ) AS (t ) AS (t t w )
C3 ( s ) V1 V2 V3
The maximum step change in the flow rate into the tank that will not overflow the Ci ( s ) A A tw s 1 e tw s
tank is 1.605. s q / V1 s q / V2 s q / V3 Ci ( s )
e A (2)
s s s
Now substitute Equation (2) into Equation (1). For simplicity, define a new variable
y
Since ci (t ) (t ), Ci (s) 1
4 f=q/V.
V1 V2 V3 V
A f (1 etw s )
3
1.
( q / V )3 1 2 3 C3 ( s)
ss f
3
3 C3 ( s )
( s q / V )3 ( s q / V ) ( s q / V ) 2 ( s q / V )3
Now use a symbolic mathematics software to find the inverse Laplace transform,
giving c3(t). The solution is formulated as a function of t, f, A, and tw. Then as an
c3 (t ) 1e(q /V )t 2te(q /V )t 3t 2e(q /V )t example, we plot c3(t) for f=1/20, A=10, and tw=1.
2
2.
V1 V2 V3 V1 In Mathematica, take the inverse Laplace transform:
If desired, the fractions can be approximated as decimals: From the Final Value Theorem, y(t) = 12 when t
y
y
12 4.4 4.5
10 0.25
dy1
8 0.20
From Exercise 4.2, a) 2 = -2y1 – 3y2 + 2u1 (1)
6 0.15
dt
s dy 2
4 0.10 Y (s) 3e = 4y1 – 6y2 + 2u1 + 4u2 (2)
dt
2 0.05 U ( s ) 10s 1
time time Rearrange, Taking Laplace transform of the above equations and rearranging,
10 20 30 40 50 10 20 30 40 50
Fig. S4.2a. Output for parts c) and d). Fig. S4.2b. Output for part e).
Y ( s )[10s 1] 3e sU ( s ) (2) (2s+2)Y1(s) + 3Y2(s) = 2U1(s) (3)
y
0.30 y
Take L-1 of (2), -4 Y1(s) + (s+6)Y2(s)=2U1(s) + 4U2(s) (4)
0.25
2 dy
0.20 10 y 3 u (t 1) (3) Solving Eqs. (3) and (4) simultaneously for Y1(s) and Y2(s),
0.15 dt
1
(2s 6)U1 ( s) 12U 2 ( s) 2( s 3)U1 ( s) 12U 2 ( s)
0.10 Take L of (3) for y (0)=4, Y1(s) =
time 2s 2 14s 24 2( s 3)( s 4)
0.05 10 20 30 40 50
10[ sY ( s ) 4] Y ( s ) 3e sU ( s )
time 1
10 20 30 40 50 (4s 12)U1 ( s) (8s 8)U 2 ( s) 4( s 3)U1 ( s) 8( s 1)U 2 ( s)
Substitute U ( s) 2 / s and rearrange to give, Y2(s) =
2s 2 14s 24 2( s 3)( s 4)
Fig. S4.2c. Output for part f). Fig. S4.2d. Output for part g).
6e s
10sY-40+Y= Therefore,
s
4.3
6e s Y1 ( s ) 1 Y1 ( s ) 6
Y (10 s 1) 40 ,
s U1 ( s ) s 4 U 2 ( s ) ( s 3)( s 4)
The transfer function for the pressure transmitter is given by, Partial fraction expansion:
Pm ( s) 1 Y2 ( s ) 2 Y2 ( s ) 4( s 1)
(1) ,
P(s) 10s 1 Y (s) e s
6
40 U1 ( s ) s 4 U 2 ( s ) ( s 3)( s 4)
s (10 s 1) (10 s 1)
and P( s) 15 / s for the step change from 35 to 50 psi. Substituting (1) and
6 a a2
rearranging gives: 1 4.6
1 15 s (10 s 1) s 10 s 1
Pm ( s)
10s 1 s Find 1: Multiply by s and set s 0 1 6 a)
Taking the L-1 gives,
From item #13 in Table 3.1, the step response is given by: Find 2 : Multiply by 10s 1 and set s 0.1 2 60
Pm (t ) 15 (1 et /10 ) (2) x (t ) 0.09e t /10 and x(t ) x x (t ) 0.3 0.09e t /10
6 6 4
Y (s) e s
s s 0.1 ( s 0.1) The intial values are x (0) 0.09 and
Let ta be the time that the alarm sounds. Then,
Pm (ta ) 45 35 10 psi (3) x(0) x (0) x 0.09 0.3 0.39.
Take L-1 ,
Substituting (3) and t=ta into (2) and solving gives ta = 11s. Thus, the alarm will y (t ) 6 S (t 1)(1 e ( t 1)/10 ) 4e t /10 The plot of the concentration response is shown in Fig. S4.6.
sound 11 seconds after 1:30PM. Check: At t =0, y (0)=4.
a) The model for the system is given by Additional assumptions a) Mass balances on the surge tanks:
The last partial derivative is zero, because qi q from the steady-state relation, and K2 q 3.153 10 4.364 10
6 7
q
Taking the Laplace transforms and rearranging, 13.84
q 13840
the derivative term in brackets is finite for all 0<h<D. Consequently, the linearized 0.15766q 1000 T2 T2
model of the process, after substitution of deviation variables, is
Vs q Vk (T )C (s) qC (s) Vc
A Ai A k (T )
20000
T ( s) (5) 8.87 10 3
dh ' 1 T2
(q 'i q ')
dt 2 L ( D h )h 20000
VCs qC (H )VcA k (T ) T 2 T ( s) (H )Vk (T )C A ( s) (6)
4.15
Recall that the term 2L ( D h)h in the previous equation represents the variable
Substituting for CA (s) from Eq. 5 into Eq. 6 and rearranging,
surface area of the tank. The linearized model treats this quantity as a constant that
Assumptions:
depends on the nominal (steady-state) operating level. Consequently, operation of
T ( s ) HVk (T )q (7) 1. Constant physical properties
the horizontal cylindrical tank for small variations in level around the stead-state
(s)
C Ai 20000 20000 2. Perfect mixing
Vs q Vk (T ) VCs qC (H )Vc A k (T ) T 2 (H )c AV k (T ) T 2
2 2
value would be much like that of any tank with equivalent but constant liquid
surface. For example, a vertical cylindrical tank with diameter D’ has a surface are
c A is obtained from the steady-state version of Eq.1, Dynamic model: Balances on cell mass and substrate concentration
of liquid in the tank equal to ( D ')2 / 4 2 L ( D h )h . Note that the coefficient
qcAi dX
1 cA 0.001155 mol/cu.ft. ( S ) X DX f1 ( S , X , D) (1)
L ( D h )h is infinite for h 0 or for h D and is a minimum at h D / 2 . q Vk (T ) dt
2 dS
Thus, for large variations in level, this equation would not be a good approximation, Substituting the numerical values of T , , C, ( H), q, V, c A into Eq. 7 and ( S ) X / YX / S D( S f S ) f 2 ( S , X , D, S f ) (2)
simplifying gives, dt
because dh/dt is independent of h in the linearized model. In these cases, the
where:
horizontal and vertical tanks would operate very differently.
T ( s) 11.38 S
( S ) X is defined as m X ,
C Ai ( s) (0.0722s 1)(50s 1) Ks S
4.14
T ( s) D is defined as
F
b) The gain K of the above transfer function is equal to ,
C Ai ( s) s 0
V
Assumptions
Linearization of (1) about the nominal steady state gives a linearized model of the
1. Perfectly mixed reactor 0.15766 q
2. Constant fluid properties and heat of reaction. K (8) form:
q c q c
3.153 106 A2 13.84 4.364.107 A2
a) Component balance for A, 1000 T 1000 T dX f1 f1 f1
S X D
dc dt S X D
V A q (c Ai c A ) Vk (T )c A (1) It is obtained by setting s=0 in Eq. 7 and substituting numerical values for ,
ss ss ss
dt
Energy balance for the tank, C, ( H), V. Evaluating sensitivities gives,
dX m ( K s S ) m S S
dT X S ' m D X ' X D (3)
VC qC (Ti T ) ( H )Vk (T )c A (2) dK K K2 q c dt ( K s S )2 Ks S
dt 2 0.01384 3153 A2 6.50 10 4
Since a transfer function with respect to cAi is desired, assume the other inputs, dq q 0.15766q 10 6 T
namely q and Ti, to be constant. Linearization of (2) about the nominal steady state:
dc A dc A dT dT dK K 2 q 3.153 10 c A 2 2 4.364 10 c A
6 7
13.84 dS
Linearize (1) and (2) and note that , , f f 2 f 2 f 2
dt dt dt dt dT 3.153 1000 T3 T3 2 S' X ' D' S f
dt S ss X ss D ss S f
dc 20000 2.57 10 5 ss
dX '
0.113S ' 2.25D ' 5.1
dt
No, the time required for the output Y(s) to reach steady state does not depend on
dS '
0.326 S ' 0.2 X ' 9 D ' 0.1S f the magnitude of the step input in U, it only depends on the time constant 𝜏𝑖 and 1000
dt delay θ . Since the Laplace transform of a step change is M/s, we have:
Taking Laplace transforms, assuming steady state initially: KM
Y s G s U s e s
s 1s 1 2 s 1 3 s 1 4 s 1
T'
sX ' ( s) 0.113S ' ( s) 2.25D' ( s )
The inverse Laplace transform takes the following form:
sS '( s) 0.326 S '( s) 0.2 X ( s) 9 D( s) 0.1 S f ( s)
500
12 22
e t / 1 e t / 2
1 2 1 3 1 4 1 2 2 3 2 4
In order to derive the transfer function between X and D, assume that Sf is Y t KMu t
32 42
constant at its nominal steady-state value, S f (t ) S f ; thus Sf 0. Rearranging e t / 3
e t / 4
1 3 3 2 3 4 1 4 4 2 4 3
gives, 0
As shown in above equation, the settling time is not related to the magnitude of 0 1 2 3 4 5 6 7 8 9 10
0.113 2.25
X ' ( s) S ' (s) D' ( s ) (5) input signal M. Time(/min)
s s
Figure S5.2. Temperature response for a ramp input of magnitude 0.5 Kw/min.
0.2 9
S '( s) X '( s) D '( s) (6) 5.2 5.3
s 0.326 s 0.126
(a) For a step change in input of magnitude M: The contaminant concentration c increases according to this expression:
Substitute (6) into (5) and rearrange gives,
y t KM 1 e t / y 0 c(t) = 5 + 0.2t
X '( s) (2.25s 1.7)
(7)
D '( s) s 2 0.326s 0.0226 We note that KM y y 0 500 100 400 C Using deviation variables and Laplace transforming,
400 C 0.2
Rearrange (7) to a standard form: Then K 400 C / Kw c(t ) 0.2t or C ( s)
2 1 Kw s2
X '( s ) K ( s 1) 400 100 Hence
2 2 a At time t 4 , y 4 400 C ; thus, 1 e4/ , or 2.89 min 1 0.2
D '( s ) s 2 s 1 500 100 C m ( s)
T ' s 400 10s 1 s 2
where: C / Kw
4-18
5-1 5-2
The time at which the actual concentration exceeds the limit (t = 10 s) is 5.5
subtracted from the previous result to obtain the requested t . 6
dT
2.5
mC hA(Ts T )
4
dt
1
1 Substituting numerical values in (1) and noting that
dT dT
0.5
0
0 2 4 6 8 10 12 14 16 18 20 Ts T and ,
time
dt dt
Figure S5.4. Exit concentration response for a rectangular input. dT
0
0 2 4 6 8 10 12 14 16 18 20
15 Ts T
dt
time
b) By inspection of Eq. 1, the time at which this function will reach its
Figure S5.3. Concentration response for a ramp input of magnitude 0.2 Kw/min. T ( s ) 1
maximum value is 2, so maximum value of the output is given by Taking Laplace transform,
c (2) 8(1 e 1 ) 8(1 e 0 / 2 ) S (0) Ts( s) 15s 1
5.4
and since the second term is zero, c(2) 5.057 b) Ts(t) = 23 + (80 23) S(t)
a) Using deviation variables, the rectangular pulse is
c) By inspection, the steady state value of c(t ) will be zero, since this is a Ts T 23
0 t<0 first-order system with no integrating poles and the input returns to zero.
cF = 2 0t<2 To obtain c() , simplify the function derived in a) for all time greater From t = 0 to t = 20,
0 2t than 2, yielding
57
Ts(t ) 57 S(t) , Ts( s )
Laplace transforming this input yields
c (t ) 8(e (t 2) / 2 e t / 2 ) s
1 57
T ( s ) Ts( s )
CF' ( s)
2
s
1 e2s which will obviously converge to zero. 15s 1 s (15s 1)
Substituting c(t ) 0.05 in the previous equation and solving for t gives Applying inverse Laplace Transform,
The input is then given by
t = 9.233 T (t ) 57(1 e t / 15 )
8 8e2 s
C ' ( s)
s(2s 1) s(2s 1) Then
and from Table 3.1 the time domain function is T (t ) T (t ) T 23 57(1 e t / 15 )
c ' (t ) 8(1 et /2 ) 8(1 e (t 2)/2 ) S (t 2) Since T(t) increases monotonically with time, maximum T = T(20).
25
Partial fraction expansion:
T'
20
a1 = 25, a2 = -15.
25 15 5.7
Y ( s)
15
5s 1 3s 1
10
Inverse Laplace: Assume that at steady state the temperature indicated by the sensor Tm is equal to
5 y (t ) 5e t / 5 5e t / 3 the actual temperature at the measurement point T. Then,
0
Tm ( s ) K 1
Then,
0 5 10 15 20 25 30 35 40 45 50
time
lim y(t ) lim 5e t / 5 5e t / 3 0 T ( s ) s 1 1.5s 1
t t
Figure S5.5. Thermocouple output for parts b) and c)
Tm T 350 C
Or, final value theorem from Chapter 3 applies:
5.6
10
lim y (t ) lim sY ( s ) lim sG ( s )U ( s ) lim s 0 Tm (t ) 15sin t
t s o s o s o (5s 1)(3s 1)
(a)
10 M (c) where =2 0.1 rad/min = 0.628 rad/min
Y ( s) G ( s )U ( s ) 10 1
(5s 1)(3s 1) s Y ( s ) G ( s )U ( s )
(5s 1)(3s 1) s 2 1 At large times when t/ >>1, Eq. 5-26 shows that the amplitude of the sensor
a a a signal is
Y ( s) ( 1 2 3 ) M a1 a2 a3 b3 j a3 b3 j
5s 1 3s 1 s Y (s) A
5s 1 3s 1 s j s j Am
2 2 1
Partial fraction expansion:
a1 = 125, a2 = -45, a3 = 10. Partial fraction expansion: where A is the amplitude of the actual temperature at the measurement
125 45 10 a1 = 625/26, a2 = -27/2, a3 = -2/13, b3 = 7/26.
Y ( s) ( )M point.
5s 1 3s 1 s 625 / 26 27 / 2 2 / 13 7 / 26 j 2 / 13 7 / 26 j
Y (s)
5s 1 3s 1 s j s j A 15 (0.628) 2 (1.5) 2 1 = 20.6C
Inverse Laplace: Therefore
y (t ) (25e t / 5 15e t / 3 10) M Inverse Laplace:
Maximum T T A =350 + 20.6 = 370.6
y (t ) 125 / 26e t / 5 9 / 2e t / 3 7 / 13sin t 4 / 13 cos t
Then,
Maximum Tcenter = 3 (max T) – 2 Twall
Therefore, the catalyst will not sinter instantaneously, but will sinter if operated 2
10
8
h'(t)
h(t)/ft
1
5.8 7
0.5 6
a) Assume that q is constant. Material balance over the tank,
5
dh
q1 q 2 q
0
A 0 5 10 15 20 25 30 35 40 45 50
dt time
4
0 5 10 15 20 25
Figure S5.8a. Liquid level response for part b) Time/min
Writing in deviation variables and taking Laplace transform
Figure S5.8b. Liquid level response for part d)
As H ( s ) Q1 ( s ) Q2 ( s ) c) h 6.122 ft at the new steady state t 12
H ( s ) 5.9
1 d) q1(t ) 10S t 5S t 12
Q1 ( s ) As
a) Material balance over tank 1.
b) q1 (t ) 5 S(t) – 5S(t-12) 10 5 12 s
Q1( s) e
s s dh
A C (qi 8.33h)
5 5 12 s 10 / A 5 / A dt
Q1 ( s) e H ( s) 2 2 e12 s
s s s s
1 5/ A 5/ A 10 where A = (4)2/4 = 12.6 ft2
H ( s) Q1( s) 2 2 e12 s 4 t 4 0.354t 0 t 12
As s s h t A
6.122 0.177t t 12 ft 3 /min
C = 0.1337
5 5 USGPM
h(t ) t S(t) (t 12) S(t-12) The liquid level will keep increasing and there will be no steady-state
A A value of liquid level h .
AsH ( s) CQi( s) (C 8.33) H ( s)
4+
5
t 4 0.177t 0 t 12 H ( s ) 0.12
A Qi( s ) 11.28s 1
h(t) =
5
4 + 12 6.122 12 < t For tank 2,
A
dh
A C (qi q )
dt
Tank 2 overflows first, at 9.4 min. K That way we can get K and
e) The red line (h’(t)=2 ft, or h(t)= 8 ft) shows that tank 2 overflows first at H ( s ) P ( s)
2 s 2 2s 1
9.4 min.
y(t)
1
Hence K = C/B =
9
2g Slope = KM
1/ 2
8 1 2L
2 then 1 / B
7 B 3g
1/ 2
A 3 2 L
6
2 then
B R 2 3g
5
h'(t)
5.12 t t
5.13 Let E = ln1 s s 3
y Ky 4 y x t
a) and find value of s that makes E 0 by trial-and-error.
a) The solution of a critically-damped second-order process to a step change
Assuming y(0) = y (0) 0 of magnitude M is given by Eq. 5-50 in text.
ts/ E
4 0.6094
Y (s) 1 0.25 t
y(t) = KM 1 1 e t / 5 -0.2082
X ( s) s 2 Ks 4 0.25s 2 0.25Ks 1 4.5 0.2047
4.75 -0.0008
b) Characteristic equation is Rearranging
a value of t = 4.75 is ts, the settling time.
s2 + Ks + 4 = 0 y t
1 1 e t / Ka a a a a4
KM b) Y(s) = 1 22 3
K K 2 16 s 2 (s 1) 2 s s s 1 (s 1) 2
The roots are s =
2
t t / y We know that the a3 and a4 terms are exponentials that go to zero for large
1 e 1
-10 K < -4 Roots : positive real, distinct KM values of time, leaving a linear response.
Response : A + B e t / 1 + C et / 2
Ka
When y/KM = 0.95, the response is 0.05 KM below the steady-state value. a2 = lim Ka
K = -4 Roots : positive real, repeated s 0 (s 1) 2
Response : A + Bet/ + C et/
KM Ka
Define Q(s) =
-4 < K < 0 Roots: complex with positive real part. (s 1) 2
t t 0.95KM
Response: A + et/ (B cos 1 2 + C sin 1 2 )
dQ 2 Ka
y
ds (s 1) 3
K=0 Roots: imaginary, zero real part.
Response: A + B cos t/ + C sin t/
1 2 Ka
Then a1 = lim
1! s 0 (s 1) 3
0<K<4 Roots: complex with negative real part.
0 ts
Response: A + e-t/ (B cos 1 2
t t
+ C sin 1 2 )
time a1 = 2 Ka
the long-time response (after transients have died out) is
K=4 Roots: negative real, repeated.
t s t /
Response: A + Be-t/ + C t e-t/ 1 e 1 0.95 0.05 y (t ) Kat 2 Ka Ka (t 2)
a(t 2) for K = 1
4 < K 10 Roots: negative real, distinct t t
ln1 s s ln(0.05) 3.00 and we see that the output lags the input by a time equal to 2.
Response: A + B et / 1 + C et / 2
Response will converge in region 0 < K 10, and will not converge in
region –10 K 0
5.15
2
Y
1.5
P ( s ) 3
C / kW
1 T ( s ) (3) 2 s 2 2(0.7)(3) s 1
0 actual response time
Note that the input change p(t ) 26 20 6 kw
0.5
(c) . 0
0 0.5 1 1.5 2
Time
2.5 3 3.5 4
Since K is 3 C/kW, the output change in going to the new steady state
1
Output
will be
0.9
Figure S5.13b Computer simulation results on part (b)
95% threshold
T (3 C / kW )6 kW 18 C
0.8 t
0.7 5.14
a) Therefore the expression for T(t) is Eq. 5-51
0.6
0.1
Overshoot = exp 0.47 , = 0.234
0 1 2
0 1 2 3 4 5 6 20
Time
2
Figure S5.13a Computer simulation results on part (a) Period = 2.3 sec
1 2 15
T'(t)
1 0.2342
= 2.3 sec 0.356 sec 10
2
R ( s ) 0.2
(1)
P ( s ) 0.127s 2 0.167s 1 5
time
0.167 R R 0.2 P
0.127 R Figure S5.15. Process temperature response for a step input
b) The overshoot can obtained from Eq. 5-53 or Fig. 5.11. From Figure 5.11 20
Time to the first peak, tp ts
we see that OS 0.05 for =0.7. This means that maximum temperature is 1 2
Therefore, ln
1 2
Tmax 70 + (18)(1.05) = 70 + 18.9 = 88.9 y (t p ) KM
b) Overshoot, OS =
KM
From Fig S5.15 we obtain a more accurate value. t
5.17
The time at which this maximum occurs can be calculated by taking OS = exp cos() sin( )
derivative of Eq. 5-51 or by inspection of Fig. 5.8. From the figure we see 1 2
that t / = 3.8 at the point where an (interpolated) =0.7 line would be.
exp exp a) Assume underdamped second-order model (exhibits overshoot)
1 1
2 2
Δoutput 15 10 ft ft
tmax 3.8 (3 min) = 11.4 minutes K 1/ 6
y (t 3 p ) KM Δinput 210 180 gal/min gal/min
c) Decay ratio, DR =
y (t p ) KM
5.16 16.5 15 1.5
3 Fraction overshoot = 0.3
where y (t 3 p ) is the time to the third peak. 15 10 5
1 2
For underdamped responses, From Fig 5.11, this corresponds (approx) to = 0.35
KM e
t 3 p /
2
1 2 1 2
DR = exp (t 3 p t p ) exp From Fig. 5.8 , = 0.35 , we note that tp/ 3.5
y (t ) KM 1 e t / cos t sin t (5-51) KM e
t p /
1 2
1 2
Since tp = 4 minutes (from problem statement, assuming first peak),
2
exp (OS) 2 t p 1 ζ2
1.19 min
1
a) At the response peaks, 2
dy
1 2 1 2 d) Consider the trigonometric identity 1/ 6 0.17
KM e t / cos t sin t G p (s)
1 2 (1.19) 2 s 2 2(0.35)(1.19)s 1 1.42s 2 0.83s 1
dt
sin (A+B) = sin A cos B + cos A sin B
b) 4 minutes might not be the first peak (as shown in Figure 5.8); thus, the
1 2 1 2 1 2
1 2 solution may be not unique.
et / sin t cos t 0 Let B = t , sin A = 1 2 , cos A =
5.18
1 1 1 (a)
2 2 2 2
0 cos t sin t KM 1
e t /
sin( A B) τ=1,ζ=0.5 .
1 2
1 2
Roots of denominator are: s 2 s 1 0 , s = -0.50+0.87j and -0.50-0.87j.
Imaginary roots suggest oscillation.
1 2 Hence for t t s , the settling time,
0 sin t sin n , t n
2 2 4
lim y (t ) lim sY ( s) lim s 2 lim 4.
1 2 t s 0 s 0 s s 1 s s 0 s 2 s 1
where n is the number of the peak.
e t /
1 2
0.05 , or
t ln 0.05 1 2 Time to first peak: πτ/ 1-ζ 2 =3.6
Overshoot: 2 2 exp(-πζ/ 1-ζ 2 )=0.652
dh h
q (t ) 0.134M 1 e t / 2.64
5.19 A Cqi
dt R
Original system provides better damping since q2 (t ) < q(t ) for t < 3.4.
c) Computer simulation result
a) For the original system, A (4) 2 / 4 12.6 ft 2
V = V1 + V2 = 2 7.07ft 2 5ft = 70.7ft3
dh h hmax = V/A = 5.62 ft
A1 1 Cqi 1
dt R1
dh h h h 0.5 5.62 ft
A2 2 1 2 R= 0.21 3
dt R1 R2 Cqi 0.1337 100 ft /min
3 dc A
0.14 b) C1 ( s) V q(c A i c A ) Vk (T )c A (1)
Original s dt
Modified 1 .5
0.12 C m ( s )
s (0.1s 1)(49s 1) Energy balance for the tank,
0.1 1
cm (t ) 1.51 (0.1e t / 0.1 49e t / 49 ) dT
(49 0.1) VC qC (Ti T ) (H R )Vk (T )c A (2)
0.08
dt
q/M
0.5 3 1.5 Since a transfer function with respect to cAi is desired, assume the other
0.06 c) C m ( s )
(49s 1) s s (49s 1) inputs, namely q and Ti, are constant. Linearize (1) and (2) and note that
0.04
cm (t ) 1.5 1 e t / 49 dc A dc A dT dT
, ,
0.02 dt dt dt dt
d) The responses in b) and c) are nearly the same. Hence the dynamics of the dcA 20000
0
0 1 2 3 4 5 6 7 8 9 10 V qcA i (q Vk (T ))cA Vc A k (T ) T (3)
Time/min
conductivity cell are negligible. dt T2
Figure S5.19 Computer simulation results on part (b) 1.5
dT 20000
VC qC H RVc A k (T ) T H RVk (T )c A (4)
dt T2
5.20
Taking Laplace transforms and rearranging
1
a) Caustic balance for the tank,
Vs q Vk (T )C (s) qC (s) Vc
A Ai A k (T )
20000
T ( s) (5)
T2
Cm'(t)
dC 20000
V w1c1 w2 c 2 wc VCs qC ( H )Vc k ( T ) T ( s ) ( H )Vk (T )C ( s ) (6)
dt T 2
R A R A
0.5
Since V is constant, w = w1 + w2 = 10 lb/min
Substituting C A (s ) from Eq. 5 into Eq. 6 and rearranging,
For constant flows,
T ( s) (H R )Vk (T )q
Part b)
Vs C ( s ) w1C1 ( s ) w2 C 2 ( s ) wC ( s ) Part c) C Ai ( s) 20000 20000
Vs q Vk (T ) VCs qC (H R )VcA k (T ) (H R )V 2c Ak 2 (T )
T 2
0
0 20 40 60 80 100
time
120 140 160 180 200
T2
C ( s) w1 5 0.5
Figure S5.20 Step responses for parts b) and c) (7)
C1 ( s) Vs w (70)(7) s 10 49s 1
c A is obtained from Eq. 1 at steady state,
C m ( s) K
, K = (3-0)/3 = 1 , 6 sec = 0.1 min qc Ai
C ( s) s 1 5.21 cA = 0.01159 lb mol/cu.ft.
(from the graph) q Vk (T )
C m ( s) 1 0.5 0.5 Assumptions: 1) Perfectly mixed reactor Substituting the numerical values of T , , C, –HR, q, V, c A into Eq. 7
C1 ( s) (0.1s 1) (49s 1) (0.1s 1)(49s 1) 2) Constant fluid properties and heat of reaction and simplifying,
(c)
For part (a), the heating rate returns to steady state after time 1, the tank
temperature will gradually return to the steady state value once the hotter fluid is
passed out of the stirred tank heater.
Figure S5.21 Step responses for the ODE system, 2nd order t.f and 1st order For part (b), the heating rate rises linearly with time and so does the outlet
approx. temperature. Physical limitations include element burnout, boiling of the liquid,
and constraints on the amount of electrical power available. However, there
should not be any short term physical limitations to the ramp, but it is an unsafe b) T (t ) t T0 (1)
situation.
Define deviation variables: Tm’ = Tm – Tss = Tm – T0 ; T’ = T – Tss = T – T0,
Step Response
1.4 substitute these into Eq. 1 :
5.24 1.2
1
T ' (t ) t and LT: T ' ( s) (2)
1 s2
0.8 As known, thermometer can be modeled by a first order system with time
Amplitude
a)
From block algebra, constant 0.1 and gain 1:
0.6
Tm ' ( s ) G 1
𝑌(𝑠) = 𝐺1 (𝑠)𝑈(𝑠) + 𝐺2 (𝑠)𝑈(𝑠) + 𝐺3 (𝑠)𝑈(𝑠) 0.4
(3) Eq. 5-19 in the book
T ' ( s ) s 1 0.1s 1
0.2
or 𝑌(𝑠) = [𝐺1 (𝑠) + 𝐺2 (𝑠) + 𝐺3 (𝑠)]𝑈(𝑠) Combine Eqs. 2 and 3 :
0
0 2 4 6 8 10 12
After some simple operations, and by account that U(s) = 1, then 1 1 1
Time (seconds)
Tm ' ( s) T ' (s) (4)
0.1s 1 0.1s 1 s 2
1 4 −3 4𝑠+1 Apply PFE to Eq. 4:
𝑌(𝑠) = [𝑠 + 2𝑠+1 + 𝑠+1] 𝑈(𝑠) = 𝑠(2𝑠+1)(𝑠+1) Figure S5.24 Step response for part (a)
5.25 a1 a a
1 1 1
or
4𝑠+1 1 4𝑠+1
𝑌(𝑠) = (2𝑠+1)(𝑠+1) × 𝑠 = 2𝑠2 +3𝑠+1 × 𝑠
1 Tm ' ( s ) T ' (s) 2 3
0.1s 1 0.1s 1 s 2 0.1s 1 s s 2
For such an integrating process at steady state, any positive/negative step 1
Notice that this system is equivalent to a step input response of an overdamped a1 2 0.01
s s 10
change in inlet flow will cause the tank level to increase/decrease with time.
( = 1.06) second-order transfer function with numerator dynamics (see 1
a3 1
Example 6.2 in your textbook). Thus, no new steady state will be attained, unless the tank overflows or 0.1s 1 s 0
1 1 a1 a a
For this example, a > 1 (e.g., 4 > 2), so the response will exhibit some empties. set s 10, 0.005; 2 3 0.015 0.1a 2
0.1s 1 s 2 0.1s 1 s s 2
overshoot.
Integrating processes do not have a steady-state gain in the usual sense. Note a 2 0.1
The system poles (-0.5, -1) lie in the LHP, so y(t) will be bounded.
that G(0) is undefined because of dividing by zero. As a result:
Finally,
K 0.01 0.1 1
K lim G ( s) lim Undefined. Tm ' ( s) 2 (5)
s 0 s 0 s 0.1s 1 s s
𝑦(0) = 𝑙𝑖𝑚 𝑠𝑌(𝑠) = 0 𝑦(∞) = 𝑙𝑖𝑚 𝑠𝑌(𝑠) = 1
s⟶∞ s⟶0 Use inverse LT to time domain:
5.26
1 t
Tm ' (t ) 0.01 exp( ) 0.1 t
0.1 0.1 Eq. 5-21 in the book
(a) Tm ' (t ) 0.1(exp( 10t ) 1) t
At time 0, Tm (0) and T(0) are the same which is t0. Then T (bath temperature) At t = 0.1 min and t = 1.0 min after the change in T(t), the difference would be:
follows a ramp:
By looking at Figure 5.5, the maximum difference occurs when t→∞ and the 1 a a a a4 a5
Y s 2
Combine Eqs. 2 and 3:
corresponding difference is: lim{0.1(exp( 10t ) 1)} 0.1 1 2 3
1 1 20 40 10 s s (2 s 1)3 s s 2 s 1 ( s 1) 2 ( s 1)3
t
Tm ' ( s ) T ' ( s) ( e )
0.1s 1 s 1 s s We know that the a3 , a4 , a5 terms are exponentials that go to zero for large values
(c) 20 20 40 40 10 s of time, leaving a linear response.
For large time, exp(-10t) approaches zero and: Tm ' ( s ) ( )( )e
s 1 s s 1 s 1 a a 1
Tm (t ) Tm (0) 0.1(exp( 10t ) 1) t Tm (0) (t 0.1) 2 1 22 a1s a2
s 2s 1 2s 1
3 3
s s
which indicates there is a 0.1 min time delay between measurement and true value Use inverse LT to time domain:
1
after a long time. Tm ' (t ) 20(1 exp( t )) 40(1 exp( (t 10)))S (t 10) (4) a2 lim 1
2s 1
s 0 3
a1 = 6
the long-time response (after transients have died out) is
y (t ) t 6
We see that the output lags the input by a time equal to 6.
Figure S5.26 T(t) and Tm(t)
5.29
Imaginary Part
(b) Since the second stage involves a phase change with a constant temperature,
thus, the time spent on phase change can be calculated based on the following 0
equation:
ρV UA T TA t -0.5
-1
1000kg / m3 0.52 1m3 334 103 J / kg 120W / m2 K 0.5 1m2 15 K t
4
-1.5
5.30
b) The process output will be bounded because there is no pole in the right
(a) From the results after 15 hr, we can see: half plane, but oscillations will be shown because of pure imagine roots.
It is a first order system, the gain K is:
c) Simulink results:
0 1 K
K 103 K / kW
0 1000 kW Solution Manual for Process Dynamics and Control, 4th edition
Copyright © 2016 by Dale E. Seborg, Thomas F. Edgar, Duncan A. Mellichamp,
and Francis J. Doyle III
6-1
5-39 5-40
6.3 a) Extremum y (t ) 0
Step response
0.5
1 1 t / 1 1 a 2 t / 2
Y ( s) K ( a s 1) KM 0 a e e 0
1 1 2 2 1 2
M
0.45
, X ( s)
0.4
X ( s) (1 s 1) s
From Eq. 6-13 1 1
0.35
1 a / 2 t
e 1 2 1 since 1>2
0.3 1 a / 1
y(t) = KM 1 1 a e t /1 KM 1 a 1 e t /1
Amplitude 0.25 1 1
0.2
1 a b) Overshoot y(t ) > KM
a) y (0 ) KM 1 a KM
0.15
1 1
1 t / 1 a 2 t / 2
0.1
b) Overshoot y(t) > KM KM 1 a e e KM
0.05 1 2 1 2
1 t / 1 1 1
KM 1 a e KM a 1 t
0
0 20 40 60 80 100 120 140 1 e 2 1 0 , therefore a>1
Time a 2
Figure S6.1c. Response of the output to a unit step input. or when, a 1 > 0 , that is, a > 1
As shown in Fig. S6.1c, the system is stable but oscillations show up c) Inverse response y (t ) 0 at t = 0+
( a 1 )
because of pure imaginary roots. y KM et /1 0 for KM 0
2
1
1 1 t / 1 1 a 2 t / 2
KM 0 a e e 0 at t = 0
+
c) Inverse response y(t) < 0
1 1 2 2 1 2
6.2
4 s 2 8 0.5s 1 1 1 1 a 2
(a) Standard form: G s e 5 s e 5 s KM 1 a 1 et /1 0 a 0
0.5s 1 2s 1 0.5s 1 2s 1 1 1 1 2 2 1 2
(b) Apply zero-pole cancellation: a 1 a
e t /1 or 1 e t /1 0 at t = 0.
8 1 1 1 1
G s e5 s a
2s 1 Therefore, a < 0. 2 1 0
Gain =8; Pole=-0.5; Zeros=None 1 2
(c)
Since 1 > 2, a < 0.
6.4
1 5 / 2s
1/1 Padéapproximation: e5 s d) If an extremum in y exists, then from (a):
1 5 / 2s Y ( s) K ( a s 1)
, 1>2, X(s) = M/s
The transfer function becomes X ( s) (1 s 1)( 2 s 1) 1 1
t 1 a 2
8 (1 5 / 2s) e 1 2
G( s) From Eq. 6-15 1 a 1
2s 1 (1 5 / 2s)
Gain = 8; poles = -0.5, -0.4; zero = 0.4
1 t / 1 a 2 t / 2 1 2 1 a 2
y (t ) KM 1 a e e t ln
1 2 1 2 1 2 1 a 1
K s
2 2
( K 2 K ) s ( K K )
6.6 ( 2 s 2 2s 1)(s 1)
K1 K K K 6.7
Y ( s) U ( s) 2 U ( s) 1 2 U ( s) Pm ( s ) %
s s 1 s s 1 d) Overall process gain = K K 6 3 3
P ( s ) s 0
psi
2
Y ( s) K1 s K1 K 2 s ( K1 K 2 ) s K1 a) p(t ) (4 2)S (t ) , P ( s)
s
U ( s) s (s 1) s(s 1) 3 3 2
Q( s) P ( s)
20s 1 20s 1 s
a) Transfer Function for the blending tank: 1.2 V is the volume of each tank.
bt s 1
C5 ( s) C ( s) 1
5 5
i
3
where K bt
qin
1 and bt
2m
2 min
0.8
,
C0 ( s) i 1 Ci1 ( s) 6s 1
q i 1m 3 / min
0.6
Then, by partial fraction expansion,
Output/Kbt
Transfer Function for the transfer line
t 1 t 2 1 t 3 1 t 4
0.4 c5 (t ) 0.60 0.15 1 et / 6 1
Gtl ( s)
K tl
6 2! 6 3! 6 4! 6
tl s 15 0.2
where:
K tl 1 0
Exact model
0.1m 3 Approximate model
tl 0.02 min -0.2
5 1m 3 / min 0 5 10 15
Time
20 25 30
Thus,
Figure S6.8. Unit step responses for exact and approximate models.
( s)
C out K bt
Cin ( s) (2s 1)(0.02s 1) 5
6.9
which is a 6th-order transfer function.
320 1 4s e3s 80 1 4s e 3s
Since bt >> tl [ 2 >> 0.02], we can approximate
1
by e-s (a) G s
b)
(0.02s 1) 5 24s 2 28s 4 6s 1 s 1
5 Gain= 80; time delay = 3; time constants 1 6, 2 1 ; poles = -1, -1/6; zeros =
where (0.02) 0.1 0.25
i 1
(b) Since a 4 0; it will show an inverse response.
( s) K bt e 0.1s
Cout
6.10
Cin ( s) 2s 1
c) Because bt 100 tl, we anticipate that this approximate TF will yield a) The transfer function for each tank is
results very close to those from the original TF (part (a)). This
approximate TF is exactly the same as would have been obtained using a C i( s ) 1
plug flow assumption for the transfer line. Thus we conclude that ,i = 1, 2, …, 5
Ci1 ( s ) V
q s 1
investing a lot of effort into obtaining an accurate dynamic model for the
transfer line is not worthwhile in this case.
Note: if bt tl , this conclusion would not be valid. where i represents the ith tank.
Denote the undelayed response (for =0) by cm (t ). Then, Let ta denote the time that the alarm goes on for the undelayed system; thus, the alarm lights up when
b) Using Simulink, c m (ta ) 25 min; i.e., when c m (ta ) 25 5 20 min . .Substituting into (7) and solving for ta by trial
cm (t ) cm (t ) (1) and error gives
0.6
c5
c4 Taking the Laplace transforms give, ta 9.24 min
0.58
c3
c2 Let ta denote the time that the alarm goes off for the system with time delay. It follows from the
0.56
c1
Cm ( s) e Cm ( s)
s
(2)
definition of a time delay that,
0.54
Concentration
The transfer functions for the delayed and undelayed systems are:
ta ta 9.24 2.00 11.24 min
0.52
Cm ( s ) e s
(3)
0.5 C ( s ) s 1
6.12
0.48 Cm ( s ) 1
(4)
C ( s ) s 1
0.46 a) Using Skogestad’s method
For the ramp input, c(t ) 2t; from Table 3.1: 5e (210.2) s 5e3.2 s
0.44 G s
0 5 10 15 20 25 30 35 40 45 50
(12s 1) 12s 1
time
2
C ( s) (5)
Figure S6.10. Concentration step responses of the stirred tank. s2
Using Simulink,
The value of the expression for c5(t) verifies the simulation results above:
Substituting (5) into (4) and rearranging gives: 5
5 5 5 2 3 4
c5 (30) 0.60 0.15 1 e 5 1 5 0.5161
2! 3! 4! 1 2
Cm ( s)
4
2 (6)
s 1 s
3
The corresponding response to the ramp input is given by Eq. 5-19 with K =1, a =2, and =10:
6.11 2
True
Approximate
-1
0 5 10 15 20 25 30 35 40 45 50
Figure S6.12 Unit step responses for the exact and approximate models.
1 Ka Kc 1 Ka Kc
6.13
Substituting for P2 ( s ) from Eq. 5 into 4, We want an approximate model of the form,
From the solution to Exercise 2.5(a) , the dynamic model for isothermal
operation is Kb Ke s
(2 s 1) Gapprox ( s)
V1 M dP1 Pd P1 P1 P2
P1( s )
1 Ka Kc (6)
s 1
(1) Pd ( s ) 12 2 1 2
RT1 dt Ra Rb s s 1
1 Ka Kc 1 Ka Kc In order for the approximate model and the original models to have the
V2 M dP2 P1 P2 P2 Pf same steady-state gain, we set K 4.
(2) To determine whether the system is overdamped or underdamped, consider the
RT2 dt Rb Rc
denominator of the transfer functions in Eqs. 5 and 6.
The largest time constant in G (s) to neglect is 1. Thus,
Taking Laplace transforms, and noting that
1
2 1 2.5
Pf ( s) 0 1 2 1 2 2
2 , 2 1 K K
1 Ka Kc a c
Approximate the smallest time constant by:
since Pf is constant, 1
e 0.4 s
0.4 s 1
K P ( s) K a P2( s) Therefore, Thus,
P1( s) b d (3)
1 s 1 1
1 (1 2 ) (1 K a K c ) 1 1 2 1 ( ) 2 2 (0.4) 1.8
K P ( s) 2
P2( s) c 1 (4) 2 (1 K a K c ) 1 2 2 2 1 (1 K a K c )
2 s 1
where
Since x + 1/x 2 for all positive x,
6.15
K a Ra /( Ra Rb )
1
K b Rb /( Ra Rb ) (1 K a K c )
From Eqs. 6-71 and 6-72,
K c Rc /( Rb Rc ) Since KaKc 0,
R2 A2 R1 A1 R2 A1 1 R1 A1 R2 A2 1
R2 A1
1
V1 M Ra Rb 2 R1 R2 A1 A2 2 R2 A2 R1 A1 2 R1 A2
1
RT1 ( Ra Rb ) Hence the system is overdamped.
1
Since x 2 for all positive x and since R1, R2, A1, A2 are positive
V2 M Rb Rc x
2
RT2 ( Rb Rc )
1
2 1 R2 A1 1
Substituting for P1 ( s ) from Eq. 3 into 4, 2 2 R1 A2
Q2(s) 1 algebraic equations to obtain the coefficients of the four partial fractions. The
6.16 Q0(s) 2 s 1 second method requires that the numerator and denominator be defined as
coefficients of descending powers of s prior to calling the MATLAB residue
function:
a) Mass balance on Tank 2: b) Mass balances on the two tanks yield (after dividing by , which is constant)
dh2 dh dh MATLAB Commands:
A2 q0 q2 A1 1 q1 A2 2 q0 q1 q2
dt dt dt
>> b = [ 36*0.262 72*0.262]
Dividing by , Valve resistance relations: b=
dh2 9.4320 18.8640
1 1
A2 q0 q2 q1 (h1 h2 ) q2 h2 >> a = conv([10 1], conv([5 1], [1 0 0.262^2]))
dt R1 R2
b=
For a linear resistance, (cf. Eq. 4-50), c) These equations clearly describe an interacting second-order system; one 50.0000 15.0000 4.4322 1.0297 0.0686
1 or more transfer functions may contain a single zero (cf. Section 6.4). For >> [r,p,k] = residue(b,a)
q2 h2 the Q2/Q0 transfer function
we know that the steady-state gain must be r=
R2
equal to one by physical arguments (the steady-state material balance 6.0865 4.9668j
Substitute, around the two tank system is q2 q0 ). 6.0865 + 4.9668j
38.1989
dh2 1
A2 q0 h2 d) The response for Case (b) will be slower because this interacting system is 50.3718
dt R2
second order, instead of first order. p=
or
dh2 0.0000 0.2620j
A2 R2 R2q0 h2 0.0000 + 0.2620j
dt
0.2000
6.17
0.1000
Introducing deviation variables and Laplace transforming yields
k=
The input is Ti(t ) 12 sin t where []
H2(s) R2
2 radians
Q0(s) A2R2 s 1 0.262 hr 1 Note: the residue function re-computes all the poles (listed under p). They
24 hours are, in reverse order: p1 = 0.1( 1 10) , p = 0.2( 2 5) , and the two purely
Because
imaginary poles corresponding to the sine and cosine functions. The
1 The Laplace transform of the input is from Table 3.1,
Q2(s) H 2(s)
residues (listed under r) are exactly the coefficients of the corresponding
R2 12 poles; in other words, the coefficients that would have been obtained via a
Ti( s) manual partial fraction expansion. In this case, we are not interested in the
s 2 2
real poles since both of them yield exponential functions that go to 0 as
we obtain,
Multiplying the transfer function by the input transform yields t .
Q2(s) 1 R2 1
Q0(s) R2 A2R2s 1 A2R2s 1 (72 36s) The two complex poles are interpreted as the sine/cosine terms using
Ti( s)
(10s 1)(5s 1)( s 2 2 ) Appendix L.
Letting 2 A2 R2 To invert, either (i) make a partial fraction expansion manually, or (ii) use the The coefficients of the periodic terms:
MATLAB residue function. The first method requires solution of a system of
b= 0, thus the exponential terms = 1. Using (L-13) and (L-15), 0.264 . (a) The mathematical model is derived based on material balance: 1
F02 1
dc R C0' s
V 1 F0 c0 RF2 c2 F1c1 Vkc1
y(t ) 12.136 cos t 9.9336sin t ... dt V 2 s 2 2 F0 F2
Vs 2 F0Vs F02 0
dc R R R
The amplitude of the composite output sinusoidal signal, for large values, V 2 F1c1 1 R F2c2 Vkc2 F0
dt When R , C2' s C0' s , equivalent to a single tank with a
of t is given by
Subtracting the steady-state equation and substituting deviation variables 2Vs F0
yields: volume = 2V.
A (12.136) 2 (9.9336) 2 15.7 The gain of above transfer function is 1.
dc '
V 1 F0c0' RF2c2' F1c1' Vkc1'
dt
Thus the amplitude of the output is 15.7 for the specified 12 amplitude
dc ' 6.20
input. V 2 F1c1' 1 R F2c2' Vkc2'
dt
The dynamic model for the process is given by Eqs. 2-45 and 2-46,
which can be written as
(b) The transfer function model can be derived based on Laplace transform:
6.18
dh 1
VsC1' s F0C0' s RF2C2' s F1C1' s VkC1' s ( wi w) (1)
dt A
150 VsC2' s F1C1' s 1 R F2C2' s VkC2' s
dT w Q
148 Solve above equations, we have: i (Ti T ) (2)
dt Ah AhC
146 F0 F1
y1 C s 2 2
'
C0' s where h is the liquid-level
V s 2Vk F1 F2 RF2 Vs F1F2 Vk 1 R F2 FVk V 2k 2
2
y2 A is the constant cross-sectional area
144 1
System outputs: h , T
142 (c) When R 0 , we have:
System inputs : w, Q
F0 F1
C2' s C0' s
y 140
V 2 s 2 2Vk F1 F2 Vs F1F2 VkF2 FVk
1 V 2k 2 Assume that wi and Ti are constant. In Eq. 2, note that the nonlinear term
138
F0 F1 dT
C' s h can be linearized as
136 Vs Vk F1 Vs Vk F2 0 dt
134 which is equivalent to the transfer function of the two tanks connected in
dT dT
series. h h
132 dt dt
130 (c) When k=0, Equation in (b) becomes:
0 1 2 3 4 5 6 7 8 dT dT
or h since 0
t F0 F1 dt dt
C2' s C0' s
V 2 s 2 F1 F2 RF2 Vs F1F2
Figure S6.18 Comparison between y1 and y2 Then the linearized deviation variable form of (1) and (2) is
Since F1 RF2 F2 , F0 F2
dh 1
w
dt A
dV CAT h
dT wi 1
Eq. 2 is simplified by substituting for from Eq. 1. Also, replace V by AT h Let
T Q dt CqF UpT h
dt Ah Ah C (where AT is the tank area) and replace A by pT h (where pT is the perimeter of
the tank). Then, Then from Eq. 7
Taking Laplace transforms and rearranging,
dh H ( s) 1 H ( s) 1 H ( s )
H ( s ) K1 H ( s ) T ( s ) T ( s) K2 AT qF q (3) , , 0
, 0 , 0 , dt QF ( s) AT s Q( s ) AT s Ts ( s )
W ( s ) s Q ( s ) W ( s ) Q ( s) 2 s 1
dT
CAT h qF C (TF T ) UpT h(Ts T ) (4) And from Eq. 9
1 1 Ah dt
where K 1 ; and K 2 , 2
A wi C wi UpT (Ts T ) C (TF T ) AT
Then, Eqs. 3 and 4 are the dynamic model for the system. s 1
T ( s ) CqF UpT h UpT (Ts T )
For an unit-step change in Q: h(t ) h T (t ) T K 2 (1 e t / 2 )
, a) Making a Taylor series expansion of nonlinear terms in (4) and introducing QF ( s ) AT s s 1
deviation variables, Eqs. 3 and 4 become:
For an unit step change in w: h(t ) h K 1t , T (t ) T
dh UpT (Ts T )
AT q F q (5)
dt T ( s )
CqF UpT h
6.21 Q( s ) AT s s 1
dT
Additional assumptions: CAT h C (TF T )qF (CqF UpT h )T UpT hTs UpT (Ts T )h (6)
dt UpT h
(i) The density and specific heat C of the liquid are constant.
T ( s ) CqF UpT h
Taking Laplace transforms,
Ts( s ) s 1
(ii) The temperature of steam, Ts, is uniform over the entire heat transfer
1 1
area. H ( s) QF ( s) Q( s) (7)
AT s AT s
(iii) The feed temperature TF is constant (not needed in the solution).
CAT h C (TF T )
Mass balance for the tank is s 1 T ( s ) QF ( s )
CqF UpT h CqF UpT h
dV UpT h UpT (Ts T )
qF q (1) Ts( s) H ( s) (8)
dt CqF UpT h CqF UpT h
Energy balance for the tank is Substituting for H ( s) from (7) into (8) and rearranging gives
d [V (T Tref )]
C q F C (TF Tref ) qC (T Tref ) UA(Ts T ) (2) CAT h C (TF T )
dt AT s s 1 T ( s ) AT s QF ( s )
CqF UpT h CqF UpT h
where Tref is a constant reference temperature and A is the heat transfer area
UpT hAT s UpT (Ts T )
Ts( s) QF ( s) Q( s) (9)
CqF UpT h CqF UpT h
UpT h (Ts T ) h q transfer function is an integrator with a negative gain. h accumulates dh2' 1 1 ' '
or C (TF T ) (h1 h2 (2)
qF changes in q, in the opposite direction, decreasing as q increases and vice dt A2 R
Substituting versa.
hAT Laplace transforming
2 h Ts transfer function is zero. Liquid level is independent of Ts and
qF
steam pressure Ps . A1RsH1' ( s) RWi ' ( s) H1' ( s) H 2' ( s) (3)
V
Let V hAT so that 2 = (initial residence time of tank)
qF T q transfer function is second-order due to the interaction with liquid A2 RsH ( s) H ( s) H ( s)
' ' '
(4)
T ( s )
2 1 2
T ( s ) level; it is the product of an integrator and a first-order process.
For and the “gain” in each transfer function is
QF ( s ) Q( s ) From (4)
T qF transfer function is second-order due to the interaction with liquid
Up (T T ) level; it has numerator dynamics since qF affects T directly as well if (A2 Rs 1) H 2' ( s) H1' ( s) (5)
K T s
AT CqF UpT h
TF T .
or
T Ts transfer function is first-order because there is no interaction with H 2' ( s) 1 1
and must have the units of temperature/volume . (The integrator s has units of (6)
liquid level. H1' ( s) A2 Rs 1 2 s 1
t-1).
To simplify the transfer function gain, we can substitute c) h qF : h increases continuously at a constant rate. where 2 A2 R
(12 ) s 2 (12 ) s H1' ( s ) R(2 s 1)Wi ' ( s ) (9) dh' 1 ' - However, the two-tank TF’s contain a pole (3 s 1) that will “filter
wi (16)
H1' ( s ) R(2 s 1) dt A out” changes in level caused by changing wi(t).
(10)
W1' ( s) s 12 s (1 2 ) - On the other hand, for this special case, we see that the zero in the first
Note that , which is constant, subtracts out.
tank transfer function ( H i' ( s) / Wi ' ( s)) is larger than the pole:
Dividing numerator and denominator by (1 2 ) to put into standard form
Laplace transforming and rearranging: 2 3 3
H1' ( s ) [ R /(1 2 )](2 s 1) Thus we should make sure that amplification of changes in h1(t)
(11)
W1' ( s) H ( s) 1/ A
'
caused by the zero do not more than cancel the beneficial filtering of
s 1 2 s 1 (17)
Wi ' ( s)
1 2
s the pole so as to cause the first compartment to overflow easily.
Again Now look at more general situations of the two-tank case:
Note that 1
K
A
H1' ( s) K (A2 Rs 1) K (2 s 1)
K
R
R
1
1
(12) H ' ( s) K (20)
1 2 A1R A2 R ( A1 A2 ) A (18) Wi ' ( s) RA1 A2 s(3 s 1)
Wi ' ( s) s s s 1
A
since A A1 A2 H 2' ( s) K
which is the expected integral relationship with no zero. (21)
Wi ' ( s) s(3 s 1)
Also, let
For either A1 0 or A2 0 ,
12 2 R 2 A1 A2 RA1 A2 b) For A1 A2 A / 2
s (13)
1 2 R ( A1 A2 ) A RA1 A2
3 0
2 AR / 2 A
(19)
so that 3 AR / 4 Thus the beneficial effect of the pole is lost as the process tends to
look more like the first-order process.
H1' ( s) K (2 s 1) Thus 2 23
(14)
Wi ' ( s) s(3 s 1) c) The optimum filtering can be found by maximizing 3 with respect
We have two sets of transfer functions: to A1 (or A2)
and
One-Tank Process Two-Tank Process RA1 A2 RA1 ( A A1 )
3
H 2' ( s) H 2' ( s) H1' (s) 1 K (2 s 1) A A
Wi ' ( s) H1' ( s) Wi ' ( s) (2 s 1) s(3 s 1) H ' ( s) K H i' ( s) K (23 s 1)
3 R
K Wi ' ( s) s Wi ' ( s) s(3 s 1) Find max 3 : ( A A1 ) A1 (1)
(15) A1 A
s (3 s 1)
H 2' ( s) K
Set to 0: A A1 A1 0
Wi ' ( s) s(3 s 1)
Transfer functions (6), (14) and (15) define the operation of the two-tank
process. Remarks: 2 A1 A
The single-tank process is described by the following equation in - The gain ( K 1/ A) is the same for all TFs. A1 A / 2
deviation variables:
- Each TF contains an integrating element. Thus the maximum filtering action is obtained when A1 A2 A / 2.
2 A2 R A order process element (1 = 0.1) will have little effect.
3 RA1 A2 A1
2.5
Output
As A1 goes to 0, 2 Ke 0.2 s
1.5
3 G ( s)
4s 2 2s 1 1
Summing up: 2
Exact model
1
Approximate model
-0.5
0 5 10 15 20 25 30
Output
3.5 1.5
6.23
3
1
2.5
The process transfer function is 0.5
Y ( s) K 0
G( s)
Output
1.5
(0.1s 1) 2 (4s 2 2s 1)
Exact model
U ( s) Approximate model
-0.5
1 0 5 10 15 20 25 30
Time
where K = K1K2. 0.5
-0.5
Approximate model
As is apparent from the plots, the smaller 1 is, the better the quality of the
2 4 2 approximation. For large values of 1 (on the order of the underdamped
0 5 10 15 20 25 30
Time
6.26
6.24 6.25
dx3
-0.8
Using a state space representation, 0.634x2 1.173x3 0.539x f
dt
-1
x Ax Bu y i 0.72xi
y = Cx + Du
-1.2
0.72 0 0 0
C 0 0.72 0 D 0
0 0 0.72 0
Y2( s) 0.4022s 0.4717 In the absence of more accurate data, use a first-order transfer function: Using Eq. 7-15,
Y f ( s) s 3 3.5190s 2 3.443s 0.8123
T '( s) Kes KM 0.339 (30.1 0.1)
0.82 min
Qi '( s) s 1 dh 1.25
Y2( s) 0.2550
dt t 0
Y f ( s) s 3 3.5190s 2 3.443s 0.8123 T () T (0) (124.7 120) F o
K 0.235 h(t ) h(0)
qi 520 500 gal/min The slope of the linear relation between ti and z i ln 1 i
c) gives
= 3:08 am – 3:05 am = 3 min h() h(0)
an approximation of (-1/), according to Eq. 7-13.
Assuming that the operator logs a 99% complete system response as “no change
after 3:34 am”, five time constants elapse between 3:08 and 3:34 am. Using h() = h(5.0) = 6 .52, the values of zi are
5 = 3:34 min 3:08 min = 26 min
= 26/5 min = 5.2 min ti zi ti zi
Therefore, 0.0 0.00 1.4 -1.92
T '( s) 0.235e3s 0.2 -0.28 1.6 -2.14
0.4 -0.55 1.8 -2.43
Qi '( s) 5.2s 1 0.6 -0.82 2.0 -2.68
To obtain a better estimate of the transfer function, the operator should log more 0.8 -1.10 3.0 -3.93
data between the first change in T and the new steady state. 1.0 -1.37 4.0 -4.62
1.2 -1.63 5.0 -
7.2
Then the slope of the least squares fit, using Eq. 7-6 is
Process gain,
1 13Stz St S z
slope (1)
h(5.0) h(0) 6.52 5.50 min 13Stt ( St )
2
K 0.339 2
qi 30.1 0.1 ft
where the datum at t = 5.0 has been ignored.
Let z1 and z2 be the natural log of the fraction incomplete response for T1
6.8 and T2, respectively. Then, 28
26
T (50) T1 (t ) 18 T1 (t )
6.6
z1 (t ) ln 1 ln
T1 (50) T1 (0)
24
8
6.4 22
T (50) T2 (t ) 26 T2 (t )
z2 (t ) ln 2 ln
20
2
T2 (50) T2 (0)
, T
6.2 6
1
T
18
T2 ' ( s) 0.75
Figure S7.2 Comparison between models a), b) and c) and the step response data. (3) 7.4
T1 ' ( s) 1.15s 1
2 1.5
0.0 Y ( s) G ( s) X ( s)
0 5 10 15 20 (5s 1)(3s 1)( s 1) s
7.3 -1.0
-2.0
Taking the inverse Laplace transform,
-3.0
a)
y(t ) (-75/8) exp(-t /5) + (27/4) exp (-t /3) - (3/8) exp(-t ) + 3
z 1,z 2
(1)
T1( s) K1 T2( s ) K2 -4.0
Q( s) 1s 1 T1( s ) 2 s 1 -5.0
a) Fraction incomplete response
-6.0
where the approximation follows from Eq. 6-58 and the fact that 1>2, as
revealed by an inspection of the data. Figure S7.3a z1 and z2 as a function of t
T (50) T1 (0) 18.0 10.0 b) Using Simulink-MATLAB, the following results are obtained:
K1 1 2.667
q 85 82
z(t)
-5.0 True data
0.9 Approximate data
-6.0
1.5
0.8
-7.0
-8.0 z(t) = -0.1791 t + 0.5734 0.7
Third-order model
y(t)
First order model
-9.0 Second order model 0.6
Response
time,t 1
0.5
0.4
Figure S7.4a Fraction incomplete response; linear regression
0.5 0.3
0.2
From the plot: slope = - 0.179 and intercept 3.2 0.1
0
0 5 10 15 20 25 30 35 40
0
Hence, time,t
0 5 10 15
Time/s
Figure S7.4b Comparison of three models for a step input
-1/ = -0.179 and = 5.6
Figure S7.5b Comparison of true data and approximate model
= 3.2 7.5
2e 3.2 s Sum of squared error = 0.0232
G(s) 1
5.6s 1 For a first-order plus time-delay model G e s , assume 1; 0.1,1,10 ,
s 1 b) / 1 :
b) In order to use Smith’s method, find t20 and t60: we have:
a) / 0.1 : 1
X: 2.918
y(t20)= 0.2 3 =0.6 0.9 Y: 0.8531
1
0.8
y(t60)= 0.6 3 =1.8 0.9
0.7
0.8 X: 2.017
Y: 0.853 0.6
Using either Eq. 1 or the plot of this equation, t20 = 4.2 , t60 = 9.0
Response
0.7
0.5
0.6
Using Fig. 7.7 for t20/ t60 = 0.47 X: 1.436
Response
0.4 Y: 0.3537
0.5
0
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
Time/s Figure S7.5c Plot of the true data; / 1
5.0 0.3
0.3 0.5 7.0 0.4
8.0 0.5
0.2
0.4 9.0 0.6
11.0 0.7
0.1
0.3 14.0 0.8
0 17.5 0.9
0 5 10 15 30.0 1.0
Time/s 0.2
0.1
For the step response of Eq. 5-48, the time constants were calculated so as
Figure S7.5d Comparison of true data and approximate model to minimize the sum of the squares of the errors between data and model
0 predictions. Use Excel Solver for this Optimization problem:
0 5 10 15
Sum of squared error = 0.1050 Time/s
1 =6.76 min and 2 = 6.95 min
Figure S7.5f Comparison of true data and approximate model
1
c) / 10 : G ( s)
Sum of squared error = 4.3070 (6.95s 1)(6.76s 1))
1
The models are compared in Fig.S7.6:
0.9
0.8 X: 11.92
Y: 0.853
7.6
0.7
0.6
Response
0
0 5 10 15
b) Smith’s method
Time/s
From the plot, t20 = 3.9 , t60 = 9.6 ; using Fig 7.7 for t20/ t60 = 0.41
Figure S7.5e Plot of the true data; / 10
= 1.0 , t60/= 2.0 , hence = 4.8 and 1 = 2 = = 4.8
1 The inlet flow rate is quickly changed from 1.5 gallon/min to 4.8 gallon/min so it
GTC ( s ) since 2 1
0.6
is a step change, Q’(s) 3.3/s:
2 s 1
Output
0.5
H ' s
K
Q' s
K 4.8 1.5 0.1337 K 0.4412
0.4 Then s 1 s 1 s s 1 s
0.3
3 s
10e Apply the inverse Laplace transform:
GHE ( s) ,
1 s 1
h' t 0.4412 K 1 et /
0.2
1s 1 2 s 1
K
46.31
Figure S7.6 Comparison of three models for unit step input Q' s s 1 2.65s 1
7.8 A comparison of the data and the two models is shown in Fig.S7.8.
7.7
(a) 63% response method
a) From the plot, time delay = 4.0 min From inspection of the data, it is obvious that there is no time delay in the system
(=0).
Using Smith’s method,
Time constant is estimated by the 63% response method:
from the graph, t20 5.6 , t60 9.1
h' 0.63 hss' h0' 0.63* 20.3 12.78 ft
t20 1.6 , t60 5.1 , t20 / t60 1.6 / 5.1 0.314 h 12.78 10.4 23.18 ft
From Fig.7.7 , 1.63 , t60 / 3.10 , 1.645 From inspection at the data, ≈ 270 min.
The process gain is calculated as:
Using Eqs. 5-45 and 5-46, 1 4.81 min , 2 0.56 min
hss' h0' 20.3 0 20.3
K 46 min/ ft 2 Figure S7.8 A comparison of the step responses of the data and the two models.
q q0 4.8 1.5 0.1337 0.4412
b) Overall transfer function
The estimated process model is:
10e4 s The Sum of Squared Errors for the two models are:
G( s) , 1 2
(1s 1)(2 s 1) H ' s K 46
SSE (63.2%) = 0.75
Q' s s 1 2.70s 1
7.12 In the integrated results tabulated below for t = 0.1, the values are shown y(k 1)
n 1
2
G s
Ke s Table S7.13 Integrated results for the first order differential equation
x(k 1)y(k ) y(k 1) y(k 1)x(k 1) y(k 1) y(k ) 2
1.5
9
10
2.356
1.929
2.180
1.744
2.340
1.912
x(k 1)
n 1
2
14 , y(k 1)
n 1
2
198.112
10
Thus t = 0.1 does improve the finite difference model making it a more
n 1
0.5
accurate approximation of the exact model.
To find a1 and b1 , use the given first order model to minimize a1 = 0.8187 , b1 = 1.0876
-0.5
0 5 10 15 20 25 30
t,sec
10
The fitted model is y(k 1) 0.8187y(k ) 1.0876x(k )
Figure S7.12 The response of the derived FOPTD model J ( y (k ) a1 y (k 1) b1 x(k 1)) 2
n 1
or y(k ) 0.8187y(k 1) 1.0876x(k 1) (1)
(c) The inverse response at the initial state is caused by a right-half plan zero and where y(k) denotes the data.
is not captured by FOPTD model. Let the first-order continuous transfer function be,
J 10
2( y (k ) a1 y (k 1) b1 x(k 1))( y (k 1) 0 Y (s) K
a1 n1
7.13
J 10 X ( s ) s 1
2( y (k ) a1 y (k 1) b1 x(k 1))( x(k 1)) 0
b1 n1
For Eq. 7-34, the discrete model is
a) Replacing by 5, and K by 6 in Eq. 7-25 Solving simultaneously for a1 and b1 gives
t / 5 t / 5
y (k ) e t / y (k 1) [1 et / ]Kx(k 1) (2)
y (k ) e y (k 1) [1 e ]6u (k 1)
b) Replacing by 5, and K by 6 in Eq. 7-22 Comparing Eqs. 1 and 2, for t=1, gives
7 0.7 Next, we generate model predictions for the calibration data (dataset basal1)
0.6 using past inputs and past model predictions, but not past output data.
Response
6 0.5 Figure S7.16a compares the calibration data and the model predictions,
0.4 where y = y - y(0). Metric S denotes the corresponding sum of squared
0.3 errors,
5
y(t)
0.2 N
S y (k ) yˆ (k )
2
0.1 (2)
4 1
0
0 1 2 3 4 5 6 7 8 9 10 ARX2 Cal
Time 50
3 basal1
2nd Order Discrete Time Model Prediction
Figure S7.15. Comparison of true data and model responses. 0 S = 1921
The result obtained using the ARX model is different from that obtained using an -50
0 1 2 3 4 5 6 7 8 9 10 FOPTD model, because the extra constraint “K=1” is not used. In other words,
y (mg/dL)
time,t
the discrete time data do not include the final steady-state value, so the calculation
Figure S7.14 Responses of the fitted model and the data gives a different gain. If more data points are added on steady-state values, the -100
y (k ) a1 y (k 1) a2 y (k 2) b1u (k 1) b2u (k 2)
-200
a) FOPTD model: (1)
Since K=1, using linear interpolation to find times corresponding to the 35.3%
-250
and 85.3% of response: a) For the model in (1), the least squares parameter estimates are given by 0 25 50 75 100 125 150 175 200 225 250
Time (min)
t35.3% 2.89; t85.3% 8.66
βˆ = X X X Y
-1 Figure S7.16a Comparison of model predictions and calibration data for the 2nd
1.3t35.3% 0.29t85.3% 1.24
T T
(2)
order discrete-time model (S is the sum of squared errors in Eq. 2).
0.67 t85.3% t35.3% 3.87
For the basal1 dataset: b) The comparison of the validation data (dataset basal2) and the corresponding
model predictions is shown in Figure S7.16b.
b) Discrete-time ARX model:
y k 0.911y k 1 0.1329u k 1 (since u k u k 1 1 )
Thus:
y (mg/dL)
-60 0 S = 2428 (basal1)
Validation data 60,282 48,194
-80
-50 (basal2)
-100
y (mg/dL)
-100
The discrete-time model is more accurate than the transfer function model for
-120
the calibration data, which is not surprising because the former has more
-140 model parameters. Although, the transfer function model is more accurate for
-150
the validation dataset, neither model is very accurate for this dataset.
-160
0 50 100 150 200 250 300 350
Time (min) -200
Figure S7.16b Comparison of model predictions and validation data for the 2nd
order discrete time model (S is the sum of squared errors in Eq. 2). -250 7.17
0 25 50 75 100 125 150 175 200 225 250
Time (min)
b) Now, consider the first-order transfer function model
Figure S7.16c Comparison of the model predictions and calibration data for the a) Fit a first-order model:
Y (s) K 1st order transfer function (S is the sum of squared errors in Eq. 2).
(3)
U (s) s 1 20
Let y = hydrocarbon exit temperature, THC
basal2 u = air flow rate, FA
1st Order Transfer Function Model Prediction
First, determine the steady-state gain, K = Δy/Δu. The output finally reaches a 0
S = 48194
Note: There is a typo in the 1st printing. The step change in u should start at 17.9
new steady state of about 250 mg/dL. For dataset basal1, the input change is -20 m3/min, not 17.0 m3/min.
u=2.5 units/day. Thus,
-40 The step response data is shown in Fig. S7.17a. The step change in u from 17.9 to
250 mg day
y (mg/dL)
K 100 21 m3/min occurs at t = 14 min. By inspection of the noisy y data, the time delay
-60
2.5 dL units is approximately = 4 min.
-80
To identify time constant τ, determine the time at which 63.2% of the total
change has occurred. This corresponds to the time at which the output, Δy has a -100
value of - 250 × 63.2% = - 158. For inspection of the data, τ = 134 min when
y= 158 mg/dL. -120
The model predictions for the model in (3) can be calculated from the step -140
0 50 100 150 200 250 300 350
response for a first-order transfer function in (5-18) Time (min)
Figure S7.16d Comparison of the model predictions and validation data for the
y (t ) KM (1 e t / ) (5-18)) 1st order transfer function (S is the sum of squared errors in Eq. 2).
The input step sizes are M = 2.5 units/day for basal1 and 1.5 units/day for c) Discussion of results
basal2.
Table S7.16 lists the calculated values of S for the two models.
t20 1.5
0.3 The model comparisons in Fig. S7.17b indicate that the two models are very
t60 5
similar and reasonably accurate. However, the low-order transfer function models
From Fig. 7.7: fail to capture the higher order dynamics of the physical furnace model that was
used to generate the step response data. The first-order model has a lower value of
2.2 the least squares index, S:
and
t60 t 5 First order model: S = 1.71 x 104
4 = 60 1.25 min
4 4 Second-order model: S = 2.01 x 104
Thus the second-order transfer function can be written in standard form as:
From the step response data, the following information can be obtained: (a) Fit a FOPTD model to the column step response data:
THC ( s ) 11.9e 4 s
y 25 K K FA ( s ) (1.25) 2 s 2 2(2.2)(1.25) s 1 Let y = distillate MeOH composition, xD
K 11.9 3
u 2.1 m3 /min m /min u = reflux ratio, R
From (5-45) and (5-46) or by factoring (e.g., using MATLAB command roots )
The step response data is shown in Fig. S7.18a with the step change in u from
y(0) = 609.5 K, y(∞) = 584.5 K; thus y = 609.5 – 584.5 = - 25 K gives:
1.75 to 2.0 occurring at t = 3950 s. By inspection of the noisy data, the time delay
THC 11.9e 4 s is ≈ 50 s.
y63.2 = 609.5 + (0.632)(-25) = 593.7 K
FA (5.2 s 1) (0.3 s 1)
From the figure, t63.2 = 19.5 – 14 - 4 = 5.5 min. Thus, the transfer function model The following information can be obtained from the step response data:
is: c) Simulations
Y ( s ) 11.9 xD(0) = 0.85, xD(∞) = 0.88;
U ( s ) 5.5s 1 thus
xD = 0.88 – 0.85 = 0.03
b) Fit a second-order model: and
y 0.03
K 0.12
Use Smith’s method to find 1 and 2. u 0.25
Also,
y20 = y(0) + (0.2)(y) = 609.5 + (0.2)(-25) = 604.5 K
y63.2 = 0.85 + (0.632)(0.03) = 0.869
y60 = y(0) + (0.6) (y) = 609.5 + (0.6)(-25) = 594.6 K
xD ( s) 0.12 e50 s
R( s ) 981s 1
K 2
Chapter 8 © 8.1
or 1
K1 K 2
The response of a PI controller to a unit step change in set point at t = 0 is shown in
Many of the problems in this chapter require determining whether a controller should be Fig. 8.6. The instantaneous change at t = 0 is Kc and the slope of the response is K1 K 2 K 2
direct-acting or reverse-acting. The following chart can help guide the thinking process Kc/I. Now consider a more general step change in the set point of magnitude M.
for these problems when also considering the style of the valve. Note that the chart P( s) 1 K K1 K 2 K 2 K 2 ( 1)
assumes all unmentioned gains are positive (measurement, I/P, etc.). K c (1 ) Kc c
E (s) Is Is
Substituting,
Table S8.1: Chart for determining if controller should be direct-acting or reverse-acting. M K1 ( K c K1 )( 1) ( 1) K c ( 1) K1
E (s)
Then the controller s
If the process gain is: And the valve is: KM KM
should be: P( s) c c 2 Then,
s Is 1
K1 K c
Positive
Fail Close (Air-to-Open) Reverse-Acting
p (t ) K c M
Kc M
t for t 0
Kv + Kc +
MV ↑, CV ↑ I
For MV ↓, want p ↓ For CV ↑, want p ↓ c) If Kc = 3 , D = 2 , = 0.1 then,
Kp + The instantaneous change at t = 0 is KcM and the slope of the response is KcM/I
Fail Open (Air-to-Close) Direct Acting From the data given in the table, the initial instantaneous change is -1.3 mA and the
For control, if CV ↑, 0.9
Kv – Kc – slope is -0.0335 mA/s for a step change of M = 2.5 mA. Thus, K1 3 27
want MV ↓ 0.1
For MV ↓, want p ↑ For CV ↑, want p ↑ 1.3 mA 1.3mA
KC 0.52
Fail Close (Air-to-Open) Direct Acting
M 2.5 mA K 2 3 (27) 30
Negative KC M
Kv + Kc – 0.0335 mA/s
MV ↑, CV ↓ I 1 = 0.1 2 = 0.2
For MV ↑, want p ↑ For CV ↑, want p ↑
Kp –
KC M 0.52(2.5mA)
Fail Open (Air-to-Close) Reverse-Acting I = = 26 s Hence
For control, if CV ↑, 0.0335 mA/s 0.05 mA/s K1 + K2 = -27 + 30 = 3
Kv – Kc +
want MV ↑
For MV ↑, want p ↓ For CV ↑, want p ↓ Because KC is negative, we classify this controller as direct acting.
K 2 1 30 0.2
2
K1 K 2 3
8.2
2s 1
Gc ( s) 3
K 2 1 0.2s 1
P ( s ) K1 K K 2 1 s K 2 K K s 1
a) K2 1 ( K1 K 2 ) 1 2
8.3
E ( s ) 1 s 1 1 s 1 1 s 1
b) Kc = K1 + K2 K2 = Kc K1 a) From Eq. 8-14, the parallel form of the PID controller is :
1 D 1
Gi ( s ) K c 1 D s
I s
K 2 1 K
D 2 D
K1 K 2 K1 K 2
t
8.7 p k i p K c r (1 i ) r , i = 1, 2, …
I
e(t ) 0.5t
0.5
E ( s)
s2
D
Kc r
The PID controller transfer function is given by (Eq. 8-14): pk t
P' ( s) 1
K C 1 D s
E (s) Is
t
Substituting gives the controller output: K c r Kc r
I
p
0.5 K C 1
P '( s ) 1 Ds
s2 I s k-1 k k+1 k+2 k+3
1 2
p '(t ) 0.5 K C t t D S (t )
2 I Figure S8.7: PID controller output response a) To eliminate derivative kick, replace (ek – ek-1) in Eq. 8-25 by - (yk-yk-1).
(Note the minus sign.)
Substituting numerical values and adding p 12 mA gives:
8.8
1 8.9
pPID (t ) 12 t 2 t 0.5S (t ) D
3 From inspection of Eq. 8-25, the derivative kick = K c r
t
(b) The equation for a PI controller is obtained by setting τD to zero. a) Proportional kick = K c r a) Let the constant set point be denoted by y sp . The digital velocity P
algorithm is obtained by setting 1/I = D = 0 in Eq. 8-27:
1
pPI (t ) 12 t 2 t b) e1 = e2 = e3 = …. = ek-2 = ek-1 = 0
3 pk = Kc(ek – ek-1)
ek = ek+1 = ek+2 = …= r
(c) The plot of the controller response for both controllers is shown in Fig. S8.7. = K c ( y sp y k ) y sp y k 1
p k 1 p
The two controllers have similarly-shaped responses. The difference is the
t = K c y k 1 y k
sudden jump at t=0 that occurs with the PID controller as a result of the p k p K c r r D r
derivative term. When the set point begins to change with a constant slope, there I t
is a step change in the error derivative from 0 to 0.5. The derivative term in the The digital velocity PD algorithm is obtained by setting 1/I = 0 in Eq. 8-
controller gives it a jumpstart right when the setpoint begins to change that the 27:
PID controller does not have. pk = Kc [(ek – ek-1) + D (ek – 2ek-1 + ek-2)]
t
= Kc [ (-yk + yk-1) + D (-yk – 2yk-1 + yk-2) ]
t
Kc
I s(D s 1) D s 1 D sI s 8
8.14
I s(D s 1) 6
1 ( I D ) s (1 ) I D s 2
4
Kc First consider qualitatively how solute mass fraction x responds to a change in steam
I s ( D s 1) 2 flow rate, S. From physical considerations, it is clear that if S increases, x will also
0 2 4 6 8 10
Time increase. Thus, if x is increasing, we want S to decrease, and vice versa. For a fail-
Cross- multiplying open (air-to-close) control valve, the controller output p should increase in order to
Figure S8.10. Step responses for both parallel and series PID controllers have S decrease. In summary, if x increases we want S to decrease, which requires
(I D s 2 I s) P ( s) K c 1 ( I D ) s (1 ) I D s 2 E ( s) with a derivative filter. an increase in controller output p; thus a direct-acting controller is required.
I s 2 (D s 1) P ( s) K c ( I s 1)( D s 1) E ( s) For an air-to-close control valve, the temperature controller output p should
increase in order to have wc decrease. In summary, if Th2 decreases we want wc to
decrease, which requires controller output p to increase; thus a reverse-acting
controller is required.
8.16 Chapter 9 ©
a) The safest conditions are achieved by the lowest temperatures and pressures
in the flash vessel.
Two pieces of information are needed to specify controller action:
i) Is the control valve fail open or fail close 9.1 VALVE 1.- Fail close (air-to-open)
ii) Is x1 > x2 or x1 < x2 VALVE 2.- Fail open (air-to-close)
VALVE 3.- Fail open (air-to-close)
If x1 > x2, then the mass balance is: a) Flow rate transmitter: VALVE 4.- Fail open (air-to-close)
15 psig - 3 psig VALVE 5.- Fail close (air-to-open)
x1w1 x2 w2 xw x ( w1 w2 ) qm(psig)= (q gpm - 0 gpm) 3 psig
400 gpm-0 gpm
x1 x2 Setting valve 1 as fail close prevents more heat from going to flash drum and
psig setting valve 3 as fail open to allow the steam chest to drain. Setting valve 3
x1 x2 = 0.03 q(gpm) 3 psig
gpm as fail open prevents pressure build up in the vessel. Valve 4 should be fail-
( x2 ) w1 x2 w2 x( w1 w2 ) Pressure transmitter: open to evacuate the system and help keep pressure low. Valve 5 should be
x2 w1 w1 x2 w2 x ( w1 w2 ) fail-close to prevent any additional pressure build-up.
x2 ( w1 w2 ) w1 x ( w1 w2 ) 20 mA - 4 mA
Pm(mA)= ( p in.Hg 10 in.Hg) 4 mA b) Vapor flow to downstream equipment can cause a hazardous situation
w1 30 in.Hg - 10 in.Hg
x x2
( w1 w2 ) mA VALVE 1.- Fail close (air-to-open)
= 0.8 p(in.Hg) 4 mA
in.Hg VALVE 2.- Fail open (air-to-close)
Since all the variables in the equation are positive, then x > x2 . The only way to Level transmitter: VALVE 3.- Fail close (air-to-open)
decrease x is to increase w2 (but x can never be less than x2). Therefore, w2 should be VALVE 4.- Fail open (air-to-close)
increased when x increases, in order to have x decrease. If the control valve is fail 5 VDC - 1 VDC VALVE 5.- Fail close (air-to-open)
open (air-to-close), then the composition controller output signal p should decrease. hm(VDC)= (h(m) - 0.5m) 1 VDC
Thus a reverse-acting controller should be selected. Conversely, for a fail close (air- 10 m - 0.5 m
Setting valve 1 as fail close (air-to-open) prevents more heat from entering
to-open) control valve, a direct-acting controller should be used. VDC flash drum and minimizes future vapor production. Setting valve 2 as fail
= 0.421 h(m) 0.789 VDC
m open (air-to-close) will allow the steam chest to be evacuated, setting valve 3
If x1 < x2, then as fail close (air-to-open) prevents vapor from escaping the vessel. Setting
Concentration transmitter:
x1w1 x2 w2 xw x ( w1 w2 ) 10 VDC - 1 VDC valve 4 as fail open (air-to-close) allows liquid to leave, preventing vapor
Cm(VDC)= (C (g/L)-3 g/L)+1 VDC build up. Setting valve 4 as fail close (air-to-open) prevents pressure buildup.
x1 x2 20 g/L - 3 g/L
x1 x2 VDC c) Liquid flow to downstream equipment can cause a hazardous situation
= 0.529 C (g/L) 0.59VDC
( x2 ) w1 x2 w2 x( w1 w2 ) g/L
VALVE 1.- Fail close (air-to-open)
x2 w1 w1 x2 w2 x( w1 w2 ) VALVE 2.- Fail open (air-to-close)
b) The gains, zeros and spans are:
x2 ( w1 w2 ) w1 x ( w1 w2 ) VALVE 3.- Fail open (air-to-close)
w1 VALVE 4.- Fail close (air-to-open)
Flow Pressure Level Concentration
x x2 VALVE 5.- Fail close (air-to-open)
( w1 w2 ) Gain 0.03 psig/gpm 0.8 mA/in.Hg 0.421 VDC/m 0.529 VDC/g/L
Zero 0 gpm 10 in.Hg 0.5 m 3 g/L Set valve 1 as fail close to prevent all the liquid from being vaporized (This
Since all the variables are positive, then x < x2 . If x increases, the controller will Span 400 gpm 20 in.Hg 9.5 m 17 g/L would cause the flash drum to overheat). Setting valve 2 as fail open will
need to decrease it to bring it back to the set point. The only way to decrease x is to
allow the steam chest to be evacuated. Setting valve 3 as fail open prevents
decrease w2 (although x can never be smaller than x1). If the control valve is fail
pressure buildup in drum. Setting valve 4 as fail close prevents liquid from
open (air-to-close), then the composition controller output signal p should increase
escaping. Setting valve 5 as fail close prevents liquid build-up in drum
in order to reduce w2. Thus the composition controller should be direct acting.
Conversely, for a fail close control valve, a reverse acting controller should be used. 9.2
qmax Cv q v
Cv gs
P Kqmax 2
1
9.4 N (1 2.44 106 3502 )106.4 (40 1.953 104 3502 ) 2
gs Cv 350
Now that all of the variables on the right hand side of the equation are 0.9
Starting from Eq. 9-7: known, plug in to solve for Cv.
q kPa m3 Then using Eq. 9-27,
Cv (1) P 450 kPa, K 0.154 3 , N 0.0865 ,
Pv (m / h) 2 h(kPa)1/2
Nf (l ) q 66.4 4.55 10 4 q 2 1 / 2
gs g s 1.2, qmax 54 m /h 3
ln
m3 101.6 0.9
The pressure drop in the valve is: 54 l 1
h ln 50
Pv P Ps (2) Cv
kPa
where 450kPa 0.154 54 2 (m 3 /h) 2 The plot of the valve characteristic is shown in Figure S9.5. From the plot of the
m3 (m3 /h) 2
Ps Kq 2 (3) 0.0865 valve characteristic for the rated Cv of 101.6, it is evident that the characteristic is
h(kPa)1/2 1.2
Solve for K by plugging in the nominal values of q and Ps . First, convert reasonably linear in the operating region 250 q 350.
3 3
m m
the nominal value of q into units of m3/h to match the metric units version 54 54
h h The pumping cost could be further reduced by lowering PDE to a value that would
of N (the parameter N = 0.0865 m3/h(Kpa)1/2 when q has units of m3/h and
m3 make Pv/Ps = 0.25 at q 320 gpm. Then PDE = 100 and for qd = 320 gpm, the
pressure has units of KPa). 0.0865
h(kPa)1/2
0.88(kPa)1/2
0.076
m3
h rated Cv = 133.5. However, as the plot shows, the valve characteristic for this design
Cv 710.5 is only slightly more nonlinear in the operating region. Hence, the selected valve
coefficient is Cv = 133.5.
ii. Level controller: As the liquid level h increases, we want the product 9.8
400 flow rate B to increase. Since the control valve is air-to-open, this
means that the controller output signal to the control valve (via the
350
I/P) should increase. Thus, the controller should be direct-acting. Configuration I: This series configuration will not be very effective because a
iii. Pressure controller: As the pressure P increases, we want the solvent large flow rate has to pass through a small control valve. Thus, the pressure drop
300
flow rate D to increase. Since the control valve is air-to-close, this will be very large and flow control will be ineffective.
means that the controller output signal to the control valve (via the
250
I/P) should decrease. Thus, the controller should be reverse-acting. Configuration II: This parallel configuration will be effective because the large
q (gpm) control valve can be adjusted to provide the nominal flow rate, while the small
200
control valve can be used to regulate the flow rate. If the small valve reaches its
maximum or minimum value, the large valve can be adjusted slightly so that the
150
9.7 small valve is about half open, thus allowing it to regulate flow again.
------- Cv = 101.6
100
- - - - Cv = 133.5 Because the system dynamic behavior would be described using deviation 9.9
50
variables, the dynamic characteristic can be analyzed by considering that the input
terms (not involving x) can be considered to be constant, and thus deviations are
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 zero. The starting form is the linear homogeneous ODE: First write down the time-domain step response for a step change of 10°C. Then
l (valve lift) M d 2x dx solve the equation to find when y(t) is equal to 5 (since the variables are in deviation
R Kx 0 variables, this represents when TM will reach 30°C).
Figure S9.5. Control valve characteristics. gc dt 2 dt
Taking the Laplace transform gives,
ym (t ) KM (1 e t / )
M 2
9.6 X ( s) s Rs K 0 where M 10o C, K 1, and 10s
gc
ym (t ) 10(1 e t /10 )
M 2 R
a) There are three control valves. The selection of air-to-close vs. air-to-open is X ( s) s s 1 0 5 10(1 e ta /10
)
based on safety considerations: Kgc K
ta 6.93s
Calculate and ζ by comparing this equation to the standard form of the second-
i. Steam control valve: Air-to-open to prevent overheating of the order model in (5-39) (keeping in mind that gc = 32.174 lbm ft/(lbf s2)).
evaporator. Therefore, the alarm will sound 6.93 seconds after 1:10PM.
ii. Level control valve (that adjusts liquid flow rate B): Air-to-open to M
0.00965s
prevent the steam coils from being exposed to the vapor space, Kg c
which could lead the coils to being burned out. 9.10
M R
iii. Pressure control valve (that adjusts solvent flow rate D): Air-to-close 2 2
Kg c K 0.1 psig
to prevent over-pressurization of the evaporator. precision = 0.5% of full scale
20 psig
R gc
155.3 accuracy is unknown since the "true" pressure in the tank is unknown
b) For the three controllers: 2 KM
0.1 psig
resolution = 0.5% of full scale
The valve characteristics are highly overdamped and can be accurately 20 psig
i. Concentration controller: As the product concentration xB increases, approximated by a first-order model obtained by neglecting the d2x/dt2 term.
we want the steam pressure, Ps to increase. Since the steam valve is ±0.1 psig
repeatability = =±0.5% of full scale
air-to-open, this means that the controller output signal to the control 20 psig
valve (via the I/P) should increase. Thus, the controller should be
direct-acting.
Assume that the gain of the sensor/transmitter is unity (i.e. there is no steady-state
measurement error). Then, 10.2
10.1
Tm ( s ) 1 Assumptions: (a) Start with a mass balance on the tank. Then solve the equation to find how much
T ( s) ( s 1)(0.1s 1) 1. Incompressible flow. time it takes for the height to decrease from 1 m to 0.25 m.
2. Chlorine concentration does not affect the air sample density.
where T is the temperature being measured and Tm is the measured value. For the 3. T and P are approximately constant. dV (t )
C h(t )
ramp temperature change: dt
The detection time, td, depends on the transportation time delay, , and the response
0.3 Adh(t )
T (t) = 0.3t (C/s) , T (s) = 2 time of the analyzer, tr = 10 s: C h(t )
s dt
Tm ( s)
1
0.3 td = + tr (1) dh C
h0.5
( s 1)(0.1s 1) s 2 dt A
Time delay can be calculated as the ratio of the volume of the tubing V divided C
h 0.5 dh dt
Tm (t ) 0.00333e10t 0.333et 0.3t 0.33 by the volumetric flow rate of chlorine q: A
0.25 tf
C
h dh
V 0.5
The maximum error occurs as t : θ (2) dt
q 1 0
A
Maximum error = |0.3t (0.3t 0.33)| = 0.33 C 2( 0.25[m] 1[m])
C
(t f 0)
where q = 10 cm3/s and, A
If the smaller time constant is neglected, the time domain response is slightly π Di 2 L C
V (3) 1m0.5 tf
different for small values of t, although the maximum error (t) does not change. 4 A
where the inside diameter Di is: A
T °C tf 1[m0.5 ]
Di = 6.35 mm – 2(0.762 mm) = 4.83 mm = 4.83 x 10-3 m C
6 (0.5) 2 [m 2 ]
tf [m0.5 ]
Substitute Di and L = 60 m into (3): 0.065[m 2.5 / min]
5
t f 12.1[min]
V = 1.10 x 10-3 m3
4
Substitute Di into (2): Therefore, the alarm will sound at 5:12:06AM
3
(b) To find how much liquid has leaked out of the tank, calculate the difference in
V 1.10 x10 m 100cm
3 3
2 θ
3
110 s
volume between the starting level and the alarm level.
q 10 cm3 /s 1 m
1 2
1
Substitute into (1): V Vh 1m Vh 0.25 m [m] 1[m] 0.25[m] 0.59m3
0 Time s 2
0 5 10 15 20 td = + tr = 110 + 10 = 120 s = 2 min
Figure S9.11. Response for process temperature sensor/transmitter. Orange solid 0.59m3 of liquid has leaked out when the alarm sounds.
Carbon monoxide (CO) is one of the most widely occurring toxic gases, especially
line is T’(t), and purple dashed line is T’m(t). for confined spaces. High concentrations of carbon monoxide can saturate a
10.3
person’s blood in matter of minutes and quickly lead to respiratory problems or
9-9 10-2
[Type here] 10-1 [Type here]
This SIS system is shown in Fig. S10.3 with conventional control loops for
The solenoid-operated valve is normally closed. But if the pressure in the Level switch: LS = 0.45
column exceeds a specified limit, the high pressure switch (PSH) activates
pressure and liquid level. Level alarm: A = 0.3
an alarm (PAH) and causes the valve to open fully, thus reducing the
Vapor pressure in the tank. Notation:
PS = the probability that the solenoid switch fails
PSH
10.5 PLS = the probability that the level switch fails
PA = the probability that the level alarm fails
PT PC
LSH PI = the probability that the interlock system (solenoid & level switch fails)
Define k as the number of sensors that are working properly. We wish to calculate
the probability that k 2 , P (k 2) . We wish to determine,
S P
P = the probability that both safety systems fail (i.e., the original system
Feed and the additional level alarm)
Because k = 2 and k = 3 are mutually exclusive events (cf. Appendix F),
Because the interlock and level alarm systems are independent, it follows that (cf.
P(k 2) P(k 2) P(k 3)
LT LC
(1) Appendix F):
LSL P = PI PA (1)
These probabilities can be calculated from the binomial distribution 1
From the failure rates, the following table can be constructed, in analogy with
S
3 Example 10.4:
P(k 2) 0.05 (0.95)2 0.135
1
2
Liquid
Component R P=1–R
3
P(k 3) 0.05 (0.95)3 0.857
0
Solenoid:
Figure S10.3: SIS system 3 Level switch: 0.45 0.638 0.362
Level alarm 0.3 0.741 0.259
10.4
Mean time between failures, MTBF: To calculate the overall system reliability, substitute R3 = 0.449 for the reliability D2
value for two D/P flowmeters (R2=0.33) in the R calculation from Exercise 10.7: Gd2
X2 X1
From (10-6) through (10-8): 5
R Ri (0.449)(0.95)(0.61)(0.55)(0.64)
~
R = 1 – P = 1 – 0.095 = 0.905 Ysp Ysp E P U X3 + +
i 1 Km +- Gc + + Y
Gv Gp
= - ln (0.905) = 0.0998 R 0.092
1
MTBF 10.0 years
Thus, the addition of the third D/P flowmeter has increased the overall system
reliability from 0.067 (for Exercise 10.7) to 0.092.
B
10.7 Gm
10.9
Let P2 = the probability that neither D/P flowmeter is working properly. Then P2 11.2
and the related reliability, R2, can be calculated as (cf. Appendix F):
Assume that the switch and solenoid are independent. From the failure rate data, 1
P2 = (0.82)2 = 0.672 the following table can be constructed, in analogy with Example 10.4: Gc ( s) K c 1
I s
R2 =1 - P2 = 1 - 0.672 = 0.33 The closed-loop transfer function for set-point changes is given by Eq. 11-36
Component R P=1–R
1
To calculate the overall system reliability, substitute R2 = 0.33 for the reliability Pressure switch 0.34 0.712 0.288 with Kc replaced by K c 1 ,
value for a single D/P flowmeter, 0.18, in the R calculation of Example 10.4: I s
Solenoid switch/valve:
1 1
K c K IP K v K p K m 1
5 H ( s) I s (s 1)
R Ri (0.33)(0.95)(0.61)(0.55)(0.64) Assume that the switch and solenoid are independent. Then, the overall reliability
H sp ( s ) 1 1
i 1
of the interlock system is, 1 K c K IP K v K p K m 1
R 0.067 I s (s 1)
R = (0.712)(0.657) = 0.468
Thus, the addition of the second D/P flowmeter has increased the overall system
where Kp = R = 1.0 min/ft2 ,
reliability from 0.037 (for Example 10.4) to 0.067. = - ln (0.468) = 0.760
and = RA = 3.0 min. Note also that 𝜏𝐼 = 𝜏 = 3.0 min.
1
MTBF 1.32 years
Solution Manual for Process Dynamics and Control, 4th edition
Copyright © 2016 by Dale E. Seborg, Thomas F. Edgar, Duncan A. Mellichamp,
and Francis J. Doyle III
10-6 10-7
11-1
s(5s 1 1.440K 2 )
C'3m(s) C'3(s)
(3 2) 1 T () lim sT ( s)
7.20K 2 Gm
( s)
For H sp ma lb sol
s s s 0 (1 1.440K 2 ) ft3
h(t ) 2.5 ft h(t ) 0.5 ft t 0.65 min Pt() K v K 2 T () , Pt() 1.03 psi (20 4)ma 2 s ma 2 s
Gm ( s ) e 2.67 e
lb sol lb sol/ft 3
(9 3) 3
h(t ) 3.0 ft h(t ) 1.0 ft t and PtK v K 2 Ti K1 T , Pt 3.74 psi ft
G IP ( s ) K IP 0.3 psi/ma
h(t ) 3.0 ft
11.4
(10 20) USGPM USGPM
Gv ( s ) K v 1.67
11.3 (12 6) psi psi
(a) Controlled variable: c3
Manipulated variable: q2 For process and disturbance transfer function:
Disturbance variable: c2 (note: q1 and c1 are kept constant.)
Gc ( s ) K c 5 ma/ma If c2 changes, then q2 must be adjusted to keep c3 at the set point. Overall material balance for the tank,
Assume m = 0, v = 0, and K1 = 1, in Fig 11.7.
USgallons dh
7.481 A q1 q 2 q3
ft 3 dt
a) Offset = Tsp () T () 5 F 4.14 F 0.86 F
K As h is held constant at 4 ft by the overflow pipe:
K m K c K IP K v 2
T ( s ) s 1
b)
Tsp ( s ) K 0 10 15 q3 (1)
1 K m K c K IP K v 2
s 1
Thus q3 25
Using the standard current range of 4-20 ma,
Component balance for the solute,
q1c1 q1c1 q1c'1 q'1 c1 ii. Km changes. The close loop transfer functions changes, hence Gc(s) needs to be
q 2 c 2 q 2 c 2 q 2 c' 2 q' 2 c 2 (3) adjusted to compensate for changes in Gm and Km. The PI controller should be
retuned.
q 3 c3 q 3 c3 q 3 c ' 3 q ' 3 c3
iii. Km remains unchanged when zero is adjusted. The controller does not need to
At steady state: be retuned.
0 q1c1 q 2 c 2 q3 c3 (4)
Put (2) into deviation variable by considering (3) and (4): To verify the linearization results, the nonlinear model is used:
d (c 3 )
7.481 Ah q1c1 q 2 c2 q3 c3
dc'3 dt
7.481 Ah q 2 c ' 2 q ' 2 c 2 q 3 c ' 3 q ' 3 c3
dt q1 q 2 q3 The results agree with linearization.
As q1 is constant, q’3 = q’2: Step response of c3 to q2: (Gain 0.077 compared with linearized gain (Kp) 0.08 in 11.5
Eq. 6)
dc'3
7.481 Ah q 2 c ' 2 q ' 2 c 2 q 3 c ' 3 q ' 2 c3
dt (a)
dc' From Eq. 11-26 we get the closed loop transfer function for set point changes
7.481 Ah 3 q 2 c' 2 q' 2 (c2 c3 ) q3 c'3 (5)
dt
𝑌 𝐾𝑚 𝐺𝑐 𝐺𝑣 𝐺𝑝
=
Taking Laplace transform and rearranging gives 𝑌𝑠𝑝 1 + 𝐺𝑐 𝐺𝑣 𝐺𝑝 𝐺𝑚
4 The block diagram for this control system is the same as in Fig.11.8.
𝑌(𝑠) = Hence Eqs. 11-26 and 11-29 can be used for closed-loop responses to
𝑠(𝑠 2 + 4𝑠 + 4)
setpoint and load changes, respectively. 11.7
Using the Final Value Theorem we get
The transfer functions G p (s ) and G d (s ) are as given in Eqs. 11-66 and
4 4 11-67, respectively. Using block diagram algebra
𝑙𝑖𝑚 𝑠𝑌(𝑠) = 2 = =1
𝑠⟶0 𝑠 + 4𝑠 + 4 4
a) Substituting for Gc, Gm, Gv, and Gp into Eq. 11-26 gives Y G d D G pU (1)
Therefore y(∞) = 1 1
Km Kc Kv
~
U Gc Ysp Y G pU (2)
Y
As 1
(c) As the step change is a unit step change, and we have shown that y(∞) = 1, s 1 GcYsp GcY
Ysp 1
we can say that offset = 0. This is consistent with the fact that the gain of 1 Kc Kv K m From (2), U ~
the overall transfer function is 1, so no offset will occur. Normally As 1 Gc G p
proportional control does not eliminate offset, but it does for this A
where (1)
integrating process. Kc Kv Km Substituting for U in Eq. 1
(d) Using Eq. 5-50 or taking the inverse Laplace transform of the response
given above we get
For a step change in the setpoint, Ysp ( s) M / s 1 G (G
c p
~
~
G p ) Y Gd (1 Gc G p ) D G p Gc Ysp
M /s Therefore,
y t 1 1 2t e 2t
Y (t ) lim sY ( s) lim s M
s 0 s 0
s 1
Y G p Gc
Substituting the value of 0.5 for t gives ~
Offset = Ysp (t ) Y (t ) M M 0 Ysp 1 Gc (G p G p )
y t 0.264
and
b) Substituting for Gc, Gm, Gv, Gp , and Gd into (11-29) gives Gd (1 Gc G p )
Y
1 1
(e) We can tell from the response derived above that the response will not be K K K D 1 Gc (G p G p )
Y ( s)
As c v m
oscillatory, since 1 . D( s) 1 s 1
1 Kc Kv K m
As
11.8
where is given by Eq. 1.
11.6 For a step change in the disturbance, D(s) M / s
The available information can be translated as follows
M /( K c K v K m ) M
For proportional controller, Gc ( s ) K c Y (t ) lim sY ( s) lim s K K K 1. The outlets of both the tanks have flow rate q0 at all times.
s 0 s 0
s(s 1) c v m
Assume that the level transmitter and the control valve have negligible 2. To ( s ) 0
M
dynamics. Then, Offset = Ysp (t ) Y (t ) 0 0
Kc Kv Km 3. Since an energy balance would indicate a first-order transfer function
Gm ( s) K m between T1 and Q0 ,
Gv ( s ) K v Hence, offset is not eliminated for a step change in disturbance.
T (t ) 2
1 e t / 1 or 1 e 12 / 1 , 1 = 10.9 min
T () 3
T1 ( s) 3 F / (0.75 gpm) 4
The given block diagram is equivalent to
Q0 ( s) 10.9s 1 10.9s 1
D
Gd
T3 ( s) (5 3) F /(0.75gpm)
2.67
for T2(s) = 0 Ysp E E'
Q0 ( s) 2s 1 2s 1 +-
P + Y
+- Gc + G
~ ~
T ( s) (78 70) F /(12 10)V 4 Y1- Y2 ~
4. 1 ~
G*(1-e-s)
V1 ( s) 10s 1 10s 1
11.10
a) Derive CLTF:
Y Y3 Y2 G3 Z G2 P
Y G3 ( D Y1 ) G2 K c E a)
CF'(s) GD ( s) K c D s
GL
Y G3 D G3 G1 K c E G2 K c E Kg/m3
with Kc >0 as the controller should be reverse-acting, since P(t) should
X1
Y G3 D (G3G1 K c G2 K c ) E E K mY ~ increase when Cm(t) decreases.
C'sp(s) C'sp(s) E P' Pv' Q'A X2
+ C'(s)
Km +- GPI + KIP Gv Gp + I/P transducer
Y G3 D K c (G3G1 G2 ) K mY Kg/m3 ma ma +
ma psig m3/min
a) The transmitter gain, Km, changes as the span changes. Thus Gm(s) Set s j : Re: −14𝜔2 + 4𝐾𝑐 − 50 = 0 (1)
changes and the characteristic equation is affected. Stability limits would Im: (𝑗) − 𝜔3 + 35𝜔 = 0 (2)
be expected to change.
5 1 K j 6
3
c
2
Kc 0 0 : Kcm 135
b) The zero on the transmitter does not affect its gain Km. Hence Gm(s) 0 : Kcm 12.5
remains unchanged and stability limits do not change. 12.5 K c 135
Re: −6𝜔2 + 𝐾𝑐 = 0 (1)
3
−5𝜔 + (1 + 𝐾𝑐 )𝜔 = 0
Changing the control valve trim changes G s . This affects the
Im: (2)
c)
0 : K cm 6
characteristic equation and the stability limits would be expected to
change as a result. To have a stable system, we have:
Kc 0
1 1.3 1
1 K c 1 (4) 0 6
Stability region
I s 0.167s 1 22.4s
5
Step response of a closed loop function Im: (𝑗) − 3.73𝜏𝐼 𝜔3 + 5.2𝐾𝑐 𝜏𝐼 𝜔 = 0 (2)
0 : cm 0.167
0.7
1
K c 0, I 0.167
0 1 2 3 4 5 6 7
K
c
0.5
c) Find I as K c
0.4
Response
11.17 25K c 25
lim 5
0.3
lim
s 1 5 N (s) 1 5K c Kc 1/ K c 5
K c
GOL ( s ) K c I
2
I s (10s 1) D( s )
0.2
Kc = 139.7 = Kcm
c) Because Kc can be much higher without the RHP zero being present, With a different initial guess ( e.g., Kc = 5, = 5), the solution is:
the process can be made to respond faster.
Frequency of oscillation is 0.579 rad/sec Kc = 11.78, = 0.5805
Again, c = 0.5805 and the stability is ensured when
b) Substituting the Pade approximation into (1) gives: 0 < Kc < 11.78
(c). 11.21 0
Kc =-1, unstable
Kc =15.33, unstable -2
150
a) To approximate GOL(s) by a FOPTD model, the Skogestad approximation -4
100
technique in Chapter 6 is used.
50 Initially, -6
response
3K c e (1.50.30.2) s 3K c e 2 s
0
GOL ( s) -8
-50
-10
-100
Skogestad approximation method to obtain a FOPTD model:
Time constant 60 + (5/2) -12
-150
Time delay 2 +(5/2) + 3 + 2 =9.5 -14
0 100 200 300 400 500 600 700 800 900 1000
-200 Then time
1 0 (1) 0.6
50
62.5s 1
response
0.5
0 Substituting s = j in (1) and using Euler's identity. 0.4
response
-50 0.3
9.5 j
e cos(9.5 ) j sin(9.5 ) 0.2
-100
gives 0.1
1 3K c cos(9.5 ) 0
Kc =15, unstable
-250
0 50 100 150 200 250 300 350 400 450 500 (2) 250
time
150 150
response
0
0
-50
Eq. 4 has infinite number of solutions. The solution for the range
response
/2 < 9.5 < 3/2 (to make sure 𝐾𝑐 is positive) is found by trial and error
-50 -100
-100
to be = 0.1749. -150
(c) Substituting the values from above and in part (a), we get, A sketch would look like this
Figure S11.26b Step response to unit step change with integral control
Kc E
1 Kc E As evident in the sketch, there is no offset for this controller.
11.27
3
1 Kc=7.9
For a step change D , applying the final value theorem:
s 2.5 Kc=8
1 Kc=8.1
8 Is
2
Offset: lim y t lim s s
0
I
c I
t s 0 3
s s 2 K s 1 1.5
1
Y
So there is no offset for PI controller.
0.5
0
11.28
-0.5
The closed loop transfer function for set point changes is given by -1
0 50 100 150
Time (min)
𝑌 𝐾𝑚 𝐺𝑐 𝐺𝑣 𝐺𝑝 Figure S11.28 y(t) responses with different Kc
=
𝑌𝑠𝑝 1 + 𝐺𝑐 𝐺𝑣 𝐺𝑝 𝐺𝑚 Figure S11.29a: System response to a unit step setpoint change. Note that the
11.29 system is stable at Kc=7.9, marginally stable at Kc=8, and unstable at Kc=8.1.
Substituting the information in the problem gives
b) PD controller:
𝑌 𝐾𝑐 (𝑠 + 3) 𝐾𝑐 (𝑠 + 3) We derive the transfer function as follows:
= = a) Proportional controller:
𝑌𝑠𝑝 (𝑠 + 1)(0.5𝑠 + 1)(𝑠 + 3) + 3𝐾𝑐 0.5𝑠 3 + 3𝑠 2 + 5.5𝑠 + 3 + 3𝐾𝑐 We derive the transfer function as follows:
Gc K c 1 D s
Y K mGc GvG p
So the characteristic equation is Y K mGc Gv G p
YSP 1 GmGc Gv G p
YSP 1 GmGc Gv G p
0.5𝑠 3 + 3𝑠 2 + 5.5𝑠 + 3 + 3𝐾𝑐 = 0 1
Kc (1)
s 1 K c 1 D s
1
3
Y Kc Kc
Substituting s=j in above equation, we have: s 1
3
s 1 K c s3 3s 2 3s 1 K c Y
3
YSP 1 K 1
3 2 3 3K c 5.5 0.5 3 j 0
s 1 c
YSP 1 K 1 s
D
1
c 3
s 1
3
Thus, we have:
The characteristic equation of (1) is the following:
3 2 3 3K c 0 11 Kc D s 1
, (4)
5.5 0.5 0 K c ,max 10 s3 3s 2 3s 1 K c 0
3
(2) s3 3s 2 3s 1 K c D K c s
Substituting s=j in above equation, we have: where Kc=10.
Chapter 12 ©
The characteristic equation of (4) is the following:
12.2
s3 3s 2 10 D 3 s 11 0 (5)
4(1 s )
Substituting s=j in above equation, we have: G Gv G p Gm
s
3 2 11 10 D 3 3 j 0
a) Let G G . Factor the model as
Thus, we have: 12.1
3 11 0 G GG
2
1
D ,min 0.0667
10 D 3 0
3
15 with:
For K = 1, 1=10, 2=5, and =0, the PID controller settings are obtained using
4
1.8
\tau_d=0.1, stable
Eq.12-14 as G 1 s, G
s
1.6 The controller design equation in (12-20) is:
1 1 2 15 1 2
1.4 Kc , I = 1+2=15 , D 3.33 Gc
* 1
f
1.2
K c c 1 2 G
with a given first-order “filter”,
response
1
The characteristic equation for the closed-loop system is
1
f
0.8
cs 1
0.6
1 1.0
1 K c 1 D s 0 Substitute,
I s (10s 1)(5s 1)
0.4
1 s
0.2
G *c
0 Substituting for Kc, I, and D, and simplifying gives 4 cs 1
0 10 20 30 40 50 60 70 80 90 100
time b) The equivalent controller in the classical feedback control configuration in
(i) 𝜏𝐷 > 𝜏𝐷,𝑚𝑖𝑛 c s (1 ) 0 Fig. 12.6(a) is:
G *c
\tau_d=0, unstable
In order for the closed loop system to be stable, the coefficients of this first-order Gc
1500
1 G *c G
polynomial in s must be positive. Thus,
1000
Substitute to give,
c > 0 Gc
1
500 and 4( c 1)
(1+) > 0 > 1. Thus Gc is a proportional-only controller.
response
Results:
-500
12.3
a) The closed loop system is stable for > 1
-1000
-1500
b) Choose c > 0 For the FOPTD model, K = 2, = 1, and = 0.2.
0 10 20 30 40 50 60 70 80 90 100
17
Apply Euler’s identity, eθj cos(θ) j sin(θ) :
16
e3u j cos(3u ) j sin(3u )
Output
15
12
0 2 4 6 8 10 12
Collect terms for the real and imaginary parts:
Time
4Kcu cos(3u ) 0 (3)
Figure S12.5. Step response data and tangent line at the inflection point.
u 4Kcu sin(3u ) 0 (4)
This estimated model parameters are:
For (3), because Kcu 0 , it follows that:
psi psi 16.9 12.0 mA
K = KIP Kv(KpKm) = 0.75 0.9 = 1.65
cos ( 3u ) = 0 3u
π
(5) mA psi 20 18 psi
2 = 1.7 min
π
u 0.5236 rad/min (6) + = 7.2 min = 5.5 min
6 1
a) Since / > 0.25, a conservative choice of c = is used. Thus, c = 2.75
From (4) – (6), 2
min. From Case H in Table 12.1:
0.5236
Figure S12.3. Comparison of PI controllers for a unit step disturbance. 0.5236 4 K cu sin ( ) 0 K cu 0.130
2 4
θ 3
1
4 12.7
Kc 2 1.76 Kc 2 2 0.24
K θ 3
c K ( c ) 5(3 ) a.i) The model reduction approach of Skogestad gives the following
2 2 2
approximate model:
3
θ θ I 4
5.5
I 6.35 min, D 0.736 min 2 2
2 2 θ e 0.028s
4(3) G(s)
D 1.1 ( s 1)(0.22s 1)
b) From Table 12.6, the AMIGO tuning parameters are: 2 2(4) 3
Since / < 0.25, an aggressive choice of c = = 0.028 is made. From
1 1 5.5 The simulated process for a step change in the set point is plotted below for both
Case I in Table 12.1 with 3 = 0, the IMC settings are:
Kc 0.2 0.45 0.2 0.45 1
K 1.65 1.7 the PI and PID controllers. Note that the PID controller was implemented in the
proper form to eliminate derivative kick (see chapter 8). 1 1 2
0.4 0.8 0.4(1.7) 0.8(5.5) Kc 21.8
I (1.7) 3.8 min K c θ
0.1 1.7 0.1(5.5)
1 2
0.5 0.5(1.7)(5.5) I 1 2 1.22, D 0.180
D 0.78 min 1 2
0.3 0.3(1.7) 5.5
c) From Table 12.4, the ITAE PID settings for a step disturbance are a.ii) To use the AMIGO tuning relations in Table 12.6, the model reduction
method of Skogestad can be used to reduce the model to a FOPDT model.
KKc = 1.357(/)-0.947 or Kc = 2.50 The time constant in the resulting FOPDT model is the largest time constant
/I = 0.842 (/)-0.738 or I = 2.75 min in the full-order model plus one half of the next biggest time constant, 1 +
D/ = 0.381 (/)0.995 or D = 0.65 min 0.5(0.2) = 1.1. The time delay in the resulting FOPDT model is half of the
second-biggest time constant in the full-order model, 0.5(0.2) = 0.1. The
other smaller time constants are neglected.
d) The most aggressive controller is the one from part c, which has the
highest value of Kc and smallest value of τI
e0.1s
G( s)
1.1s 1
12.6
The AMIGO rules for a PID controller in Table 12.6 give:
The model for this process has K=5, τ=4, and θ=3. The PI controller parameters 1 1.1
Kc 0.2 0.45 0.2 0.45 5.15
for an FOPDT model using IMC tuning are given by entry G in Table 12.1: K 0.1
Figure S12.6: Responses to a step change in the set point at t = 1for PI and PID 0.4 0.8 0.4(0.1) 0.8(1.1)
controllers. I (0.1) 0.44
Kc
4
0.13 0.1 0.1 0.1(1.1)
K ( c ) 5(3 3) The PID controller allows the controlled variable to reach the new set point more 0.5 0.5(0.1)(1.1)
D 0.049
I 4 quickly than the PI controller, due to its larger Kc value. This large Kc allows an 0.3 0.3(0.1) 1.1
initially larger response from the controller during times from 1 to 4 minutes. The
The parameters for a PID controller are given by entry H in Table 12.1: reason that the Kc can be larger is that, after the controlled variable begins to change b) The simulation results shown in Figure S12.7 indicate that the IMC
and move toward the set point, the derivative term can “put on the brakes” and slow controller is superior for a step disturbance due to its smaller maximum
down the aggressive action so the controlled variable lands nicely at the set point. deviation and lack of oscillations. This result makes sense, given that we
made an aggressive choice for τC for the IMC controller.
3
Parallel form
Series form
2.5
y(t) 1.5
Figure S12.7. Closed-loop responses to a unit step disturbance at t=1. As shown in Figure S12.8, as β increases the set-point response becomes faster
but exhibits more overshoot. The value of β =0.5 seems to be a good choice. The 0
0 50 100 150 200 250 300
disturbance response is independent of the value of β. Time
12.8
Figure S12.9. Closed-loop responses for parallel and series controller forms.
12.9
From Eq. 12-40 (with γ=0):
1 dy 12.10
p (t ) p K c ysp (t ) ym (t ) K c e(t*)dt * D m
I dt a) From Table 12.2, the controller settings for the series form are:
This control law can be implemented with Simulink as follows: a) Block diagram
: K c K c 1 D 0.971 X1'
I Gd
CONTROLLER
I I D 26.52
X'sp X'sp E P' W'2 X'
+
I D Km +- GC Gv Gp +
WEIGHTING FACTOR
D 2.753 (mA) (mA) (mA) (Kg/min)
SET-POINT PROPORTIONAL CONTROLLER
I D
+
K
ß +- ACTION C
OUTPUT
+
Closed-loop responses generated from Simulink are shown in Fig. S12.9. X'm
The series form results in more oscillatory responses; thus, it produces more Gm
INTEGRAL
+ ACTION
aggressive control action for this example. (mA)
-
b) By changing the derivative term in the controller block, the Simulink results
show that the system becomes more oscillatory as D increases. For the
CONTROLLER
INPUT
b) Process and disturbance transfer functions: 15 3 300 / 1.2 187.5 d) Figure S12.10d indicates that c 1.57 min gives very good results.
Gv ( s )
20 4 0.0833s 1 0.0833s 1
Overall material balance: 0
w1 w2 w 0 (1) Controller:
Solute balance:
-0.02
w1 x1 w2 x 2 wx V
dx 187.5 2.6 10 4
(2) Let G Gv G p G m 32e s
dt 0.0833s 1 1 4.71s
Substituting (1) into (2) and putting into deviation variables: -0.04
1.56e s
dx then G
w1 x1 w2 x2 w1 x w2 x w2 x V (4.71s 1)(0.0833s 1) y(t) -0.06
dt
For process with a dominant time constant, c dom / 3 is recommended.
Taking the Laplace transform: -0.08
Hence . c 1.57 min. From Table 12.1
w1 X 1(s) ( x2 x )W2(s) (w1 w2 Vs) X (s)
Kc = 1.92 and I = 4.71 min -0.1
Finally:
x2 x c) Simulink results: -0.12
0 2 4 6 8 10 12 14 16 18 20
X ( s) x2 x w w2
G p (s) 1
Time
0.04
W2( s) w1 w2 Vs 1 s Figure S12.10d. Closed-loop response for set-point change.
0.035
w1
X ( s) w1 w w2 e) Improved control can be obtained by adding derivative action: D 0.4 min.
Gd ( s) 1 0.03
X 1 ( s ) w1 w2 Vs 1 s
0
0.025
V
where
w1 w2 y(t) 0.02 -0.02
2.6 10 4 0.01
G p ( s)
y(t)
1 4.71s -0.06
0.005
0.65
Gd ( s) -0.08
1 4.71s 0
0 5 10 15 20 25 30 35
Time
Composition measurement transfer function: -0.1
Final control element transfer function: Figure S12.10e. Closed-loop response after adding derivative action.
1
y(t)
Figure S12.10f. Closed-loop response for = 3 min. c) The Simulink simulation results are shown in Fig. S12.12. Both controllers Substitute the Padé approximation,
are rather aggressive and produce oscillatory responses. The IMC controller
is less aggressive (that is, more conservative). 1 0.85s
e s
12.11 1.8 1 0.85s
The PID controller may have to be retuned if any of the transfer functions, GI/P, 0.8 4.675s2 + (6.35 –1.403Kc)s + 1 + 1.65Kc = 0
Gv, Gp or Gm, change.
0.6
Substitute s=ju at Kc = Kcu:
a) Gm changes. The controller may have to be retuned. 0.4
b) The zero does not affect Gm. Hence the controller does not require 0.2 4.675 u2 + j(6.35 1.403Kcu) u + 1 +1.65Kcu = 0 + j0
retuning. 0
0 10 20 30 40 50 60 70 Equate real and imaginary coefficients,
c) Gv changes. Retuning may be necessary. Time
c) Simulation results
Ignoring u= 0, the approximate values are: 12.16
The closed-loop responses for the IMC and RAT controller settings and a step
Kcu = 4.53 and u= 1.346 rad/min change in feed composition from 0.5 to 0.55 are shown in Figs. S12.16a and
a) IMC design: S12.16b, respectively.
2
Pu 4.67 min From Table 12.1, Controller H with c = /2 =3.28 min is:
u
0.8
Note that the values of Kcu and u are approximate due to the Padé
K c 0.00802 600
approximation. By using Simulink, more accurate values can be obtained 0.7
by trial and error. For this case, no Padé approximation is needed and: θ 2 590
τI τ 6.5 7.5 min 0.6
Kcu = 3.76 Pu = 5.9 min 2 2
580 0.5
0 20 40 60 80 0 20 40 60 80
The Z-N PI settings from Table 12.7 are: τθ (6.5)(2)
τD 0.867 min t (min) t (min)
2τ+θ 2(6.5) 2
Kc = 1.69 I = 4.92 min (more accurate) 1.6 19
The sight glass has confirmed that the liquid level is rising. Because the controller K c 0.6 K cu 0.0071
output is saturated, the controller is working fine. Hence, either the feed flow is
higher than recorded, or the liquid flow is lower than recorded, or both. Because Pu Pu
the flow transmitters consist of orifice plates and differential pressure transmitters, I 7 min, D 1.75 min
2 8
a plugged orifice plate could lead to a higher recorded flow. Hence, the liquid-flow-
transmitter orifice plate would be the prime suspect.
a) IMC design:
Pu P
I 425 s, D u 119 s
From Table 12.1, Controller H with c = /2 =381 is: 2 8
K c 14.7
θ 138
τI τ 762 831min
2 2
τθ (762)(138)
τD 63.3min
2τ+θ 2(762) 138
Fig. S12.16b. Performance of the RAT controller for a step change in The distillation column model includes an RAT option for the xB control loop,
hydrocarbon flow rate from 0.035 to 0.040 m3/min. but not the xD control loop. Thus, the Simulink diagram must be modified by
copying the RAT loop for xB and adding it to the xD portion of the diagram.
The RAT controller is superior due to its smaller maximum deviation and shorter
settling time. Also, the parameters for the relay block must be changed. The new Simulink
diagram and appropriate relay settings are shown in Fig. S12.17a. The results
d) Due to the high noise level for the xD response, it is difficult to obtain from the RAT are shown in Fig.S12.17b.
improved controller settings. The RAT settings are considered to be
satisfactory. From the documentation from the RAT results, it follows that:
a = 5.55 x 10-3, d = 0.2 Fig. S12.17a. Modified RAT Simulink diagram and relay settings.
From (12-46),
4d 4(0.2)
K cu 45.9
a (5.55 x103 )
Pu 950 s
Fig. S12.17c. Performance of the IMC-PID controller for a step change in feed
composition from 0.5 to 0.55.
c) Simulation results Fig. S12.17d. Performance of the RAT controller for a step change in feed
composition from 0.5 to 0.55.
The closed-loop responses for the IMC and RAT controller settings and a step
change in feed composition from 0.5 to 0.55 are shown in Figs. S12.17c and
S12.17d, respectively.
d) Due to the high noise level for the xD response, it is difficult to obtain
improved controller settings. The RAT settings are considered to be
satisfactory.
13.2 Conclude: NE = 5
(a) The larger the tank, the more effective it will be in “damping out”
disturbances in the reactor exit stream. A large tank capacity also provides Degrees of freedom: = NF = NV – NE = 8 5 = 3
Exit flow rate w4 has no effect on x3 or x4 because it does not change the a large feed inventory for the distillation column, which is desirable for
relative amounts of materials that are blended. The bypass fraction f has a periods where the reactor is shut down. Thus a large tank is preferred from
dynamic effect on x4 but has no steady-state effect because it also does not a process control perspective. However a large tank has a high capital cost, Disturbance variable: q 6 ND = 1
change the relative amounts of materials that are blended. Thus, w2 is the so a small tank is appealing from a steady-state, design perspective. Thus,
best choice. the choice of the storage tank size involves a tradeoff of control and design NF = NFC + ND
objectives. NFC = 3 1 = 2
13.3
(b) After a set-point change in reactor exit composition occurs, it would be 13.7
desirable to have the exit compositions for both the reactor and the storage
Both the steady-state and dynamic behaviors need to be considered. From a tank change to the new values as soon as possible. But concentration in the
steady-state perspective, the reflux stream temperature TR would be a poor storage tank will change gradually due to its liquid inventory. The time Consider the following energy balance assuming a reference temperature
choice because it is insensitive to changes in xD, due to the small nominal constant for the storage tank is proportional to the mass of liquid in the tank of Tref = 0:
value of 5 ppm. For example, even a 100% change in from 5 to 10 ppm (cf. blending system models in Chapters 2 and 4). Thus, a large storage tank
would result in a negligible change in TR. Similarly, the temperature of the will result in sluggish responses in its exit composition, which is not
top tray would be a poor choice. An intermediate tray temperature would be Heat exchanger:
desirable when frequent set-point changes are required. In this situation, the
more sensitive to changes in the tray composition but may not be storage tank size should be smaller than for case (a).
representative of xD. Ideally, the tray location should be selected to be the Cc (1 f )wc (Tc 0 Tc1 ) Ch wh (Th1 Th 2 ) (1)
highest tray in the column that still has the desired degree of sensitivity to
composition changes. Overall:
13.6
The choice of an intermediate tray temperature offers the advantage of early Cc wc (Tc 2 Tc1 ) Ch wh (Th1 Th 2 ) (2)
detection of feed disturbances and disturbances that originate in the
stripping (bottom) section of the column. However, it would be slow to Variables : q1, q2,…. q6, h1, h2 Nv= 8 Mixing point:
respond to disturbances originating in the condenser or in the reflux drum.
But on balance, an intermediate tray temperature is the best choice. wc (1 f )wc fwc (3)
Equations :
13.4 Thus,
Since three flow rates ( fw1 , w2 and w3 ) can be independently adjusted, it be demonstrated by a more complicated theoretical analysis or by
NE =3 , NV = 8 ( f , wc , wh , Tc1 , Tc 2 , Tc 0 , Th1 , Th 2 ) would appear that there are three control degrees of freedom. But the bypass simulation studies.
flow rate, fw1 , has no steady-state effect on x4. To confirm this assertion,
NF =NV NE = 8 3 = 5 (b) Consider the steady-state model in (1) and (2). Substituting (1) into
consider the overall steady-state component balance for the tank and the (2) and solving for T gives:
mixing point:
NFC = 2 (f, wh)
also x1 w1 x2 w2 x4 w4 wcTc whTh
(5)
T (5)
wc wh
ND = NF NFC = 3 (wc, Tc1, Tc2)
This balance does not depend on the fraction bypassed, f, either directly or
indirectly, Since w does not appear in (5), it has no steady-state effect on T.
The degree of freedom analysis is identical for both cocurrent and Consequently, the proposed multiloop control strategy is not
countercurrent flow because the mass and energy balances are the same for Conclusion: NFC = 2 (w2 and w4) feasible.
both cases.
13.9 13.10
13.8
(a) In order to analyze this situation, consider a steady-state analysis. (a) Model degrees of freedom, NF
The dynamic model consists of the following balances:
NF = NV - NE (13-1)
Assumptions:
Mass balance on the tank:
1. Steady-state conditions with w, Th, and Tc at their nominal NV = 11 (xF, TF, F, wL, L, wV, V, T, P, h, VT)
dh values. where TF is the feed temperature and VT is the volume of the
A (1 f ) w1 w2 w3 (1) 2. Constant heat capacities
dt flash separator.
3. No heat losses
Component balance on the tank: 4. Perfect mixing NV = 7:
Mass balance
d (hx3 ) Steady-state balances: Component balance
A (1 f ) x1 w1 x 2 w2 x3 w3 (2) Energy balance
dt
wc wh w (1) Vapor-liquid equilibrium relation
Valve relations (2)
Mixing point balances: wcTc whTh wT (2) Ideal gas law
w4 = w3 + fw1 (3)
Assume that T=Tsp, where Tsp is the set point. Thus, NF = 11 – 7 = 4
x4w4 = x3w3 + fx1w1 (4)
wc wh w (3) (b) Control degrees of freedom, NFC
Thus,
wcTc whTh wTsp (4) NF = NFC + ND (13-2)
NE = 4 (Eqs. 1- 4)
Equations (3) and (4) are two independent equations with two Typically, some knowledge of the feed conditions would be
NV = 9 (h, f , w2 , w3 , w4 , x1 , x2 , x3 , x4 ) unknown variables, wh and wc. For any arbitrary value of Tsp, these available. We consider two cases:
equations have a unique solution. Thus the proposed multiloop
control strategy is feasible. Case 1: xF and TF are disturbance variables
NF = NV NE = 5
This simple analysis does not prove that the liquid level h can also Here ND =2 and:
be controlled to an arbitrary set point hsp. However, this result can
Substituting the numerical value of the frequency: Since the angle of e j is negative and becomes unbounded as ,
we see that the Pade representation also provides the better approximation
3 16 1 2
to the time delay element's phase angle, matching of the pure time delay
Aˆ AR A 2 0.12 2 0.24
4 64 2 1 element to a higher frequency than the two-term representation.
Thus the amplitude of the resulting temperature oscillation is 0.24 degrees.
14.3
14.2
First approximate the exponential term as the first two terms in a truncated 128 F 120 F
Nominal temperature T 124 F
Taylor series 2
1
Aˆ (128 F 120 F) 4 F
e s 1 s 2
5sec., 2(1.8 / 60 sec) 0.189 rad/s
Then G( j) 1 j
and ARtwo term 1 () 2 1 2 2 Using Eq. 13-2 with K=1,
Solution Manual for Process Dynamics and Control, 4th edition Actual maximum air temperature = T A 129.5 F
Copyright © 2016 by Dale E. Seborg, Thomas F. Edgar, Duncan A. Mellichamp,
and Francis J. Doyle III Actual minimum air temperature = T A 118.5 F
14.4 14.7
Bode Diagram
G(s)
Tm ( s ) 1 Using MATLAB G(s) w ith Pade approx.
Magnitude (abs)
T ( s ) 0.1s 1 0
10
5
Bode Diagram
T ( s) (0.1s 1)Tm ( s) 10
0
phase angle of T = + tan-1(0.1) = + 0.02
-1
10 10
Since only the maximum error is required, set = 0 for the comparison of T and -5
10
0
Phase (deg)
-100
-10
Error = Tm T =3.464 sin (0.2t) – 3.465sin(0.2t + 0.02) 10
-150
-45
= 3.464 sin(0.2t) –3.465[sin(0.2t) cos 0.02 + cos(0.2t)sin 0.02]
-200
= 0.000 sin(0.2t) 0.0693 cos(0.2t) -90
-250
Phase (deg)
-135
Since the maximum absolute value of cos(0.2t) is 1, -2 -1 0 1
10 10 10 10
-180
maximum absolute error = 0.0693 Frequency (rad/sec)
-225 Figure S14.8. Bode diagram for G(s) and G(s) with Pade approximation.
As we can see from the figures, the accuracy of Pade approximation does not
-270
-1 0 1 2 change as frequency increases in magnitude plot, but it will be compromised in
10 10 10 10
14.5 the phase plot as frequency goes higher.
Frequency (rad/sec)
(b) No, cannot make 2nd order overdamped system unstable for closed-loop. The value of that yields a -180 phase angle and the value of AR at that =2f where f is in cycles/min
frequency are: For the standard thermocouple, using Eq. 14-13b
(c) Yes, 3rd order system can be made unstable.
= 0.807 rad/sec 1 = -tan-1(1) = tan-1(0.15)
(d) Yes, anything with time delay can be made unstable. AR = 0.202
Phase difference = 1 – 2
14.8
Thus, the phase angle for the unknown unit is
14.6
Using MATLAB,
2 = 1
Engineer A is correct. and the time constant for the unknown unit is
Second order overdamped process cannot become unstable with a proportional 1
controller. 2 = tan( 2 )
FOPTD model can become unstable with a large Kc due to the time delay.
overflow is not a concern, the two-tank system will provide a smaller amplitude in Case i
H 1 ( s) / h1max 0.01 K the output flow for those values of that satisfy Case ii(a)
Magnitude (abs)
0 Case ii(b)
Q1i ( s) 1.32s 1 s 1 10
Case iii
H 2 ( s) / h2 max Aˆ q 2 Aˆ q
0.01 K
Q1i ( s) (1.32s 1) 2 (s 1) 2
-1
10
0.1337A 0.1337A
or
Q2 ( s) 0.1337 0.1337 ( 1)
2 2 2
4 1
2 2
-2
10
Q1i ( s) (1.32s 1) 2 (s 1) 2 90
or 2 / = 1.07
and for the one-tank system 0
Therefore, the two-tank system provides better damping of a sinusoidal
Phase (deg)
H ( s) / hmax
0.01 K disturbance for 1.07 if and only if
-90
Q1i ( s) 2.64s 1 2s 1
1.32 2 2 1
Aˆ h1 / h1max 1 , that is, A -180
Q ( s) 0.1337 0.1337 0.01
Q1i ( s) 2.64s 1 2s 1 -270
-2 -1 0 1
10 10 10 10
For a sinusoidal input q1i (t ) A sin t , the amplitudes of the heights and flow
Frequency (rad/sec)
rates are
Figure S14.11. Bode plot for each of the four cases of numerator dynamics.
Gc ( j) K c
D
4 D 4 D
c
For a third order process, a Kc can always be chosen to make the process unstable.
A stability analysis would verify this but was not required.
b) From Eq. 8-15 with I = 4D and = 0.1
Substitute s = jω into Eq. (1), we have:
(4 D s 1) D s 1 ( j 1)(50 j 1)(2 j 1) K c (4)(2)(0.6) 0
Gc ( s) K c
4 D s0.1 D s 1 For τ = 1, we have:
(−100𝜔3 + 53𝜔)𝑗 + (1 + 4.8𝐾𝑐 − 152𝜔2 ) = 0 (2)
16 2 2 1 2 2 1 Thus, we have ωc = 0.53 and Kcu = 16.58.
Gc ( j) K c D
D
(b) False. If the open loop system is first order, increasing Kc cannot result in
c
AR/K
oscillation.
0
(c) Generally true. Increasing the controller gain can cause real part of the roots of
10 the characteristic polynomial to turn positive. However, for first or second order
processes, increasing Kc will not cause instability.
Parallel controller (asymptote)
(d) Always true. Increasing the controller gain will decrease offset. However, if
Parallel controller (actual) the gain is increased too much, oscillations may occur. Even with the oscillations
Series controller (actual)
10
-1 the offset will continue to decrease until the system becomes unstable.
-2 -1 0 1 2
10 10 10 10 10
D
Figure S14.12. Nominal amplitude ratio for parallel and series controllers
Cannot become unstable – max phase angle (GOL) is -180 degree while at Can become unstable due to time delay at high frequency.
low frequency the integrator has -90 degree phase angle.
(c)
14.16
By using MATLAB,
2
Bode Diagram
10
1
10
Magnitude (abs)
0
10
Figure S14.15d Bode plot of GOL. (Kc = 10)
-1
10
1 s (1 s) K c
GOL GGc Kc
(4s 1)(2s 1) (4s 1)(2s 1) -2
10
0
Can become unstable – max phase angle (GOL) is -270 degree. Parallel
-45
Series w ith filter
Figure S14.15c Bode plot of GOL. (Kc = 10) -90
Phase (deg)
(e)
s 1 2s 1 ( s 1) K c -135
GOL GGc Kc
(4s 1)(2s 1) s s(4s 1) -180
-225
Cannot become unstable – lead lag unit has phase lag larger than -90,
-270
integrator contributes -90 degree; the total phase angel is larger than -180. 10
-2
10
-1
10
0
Frequency (rad/sec)
Figure S14.16 Bode plot for Exercise 13.8 Transfer Function multiplied by
PID Controller Transfer Function. Two cases: a)Parallel b) Series with Deriv.
Filter (=0.2).
.
Amplitude ratios:
0
a) Method discussed in Section 6.3: 10
0
Magnitude (abs)
Magnitude (dB)
-50
12e 0.3s -1
Gˆ 1 ( s)
10
G(s) G(s)
(8s 1)(2.2s 1) -100
G1(s) G'(s)
-2
-150 G2(s) 10
Method discussed in Section 7.2.1: G3(s)
4
-200 x 10 0
0
Step response of G(s) -360
12
-720
-1.152
Phase (deg)
-1080
Phase (deg)
10 -1440
-2.304
-1800
X: 10.09
8 Y: 7.2 -3.456 -2160
-2520
Value
6 -4.608 -2880
-3 -2 -1 0 1 2 3 -2 -1 0 1
10 10 10 10 10 10 10 10 10 10 10
Frequency (rad/s) Frequency (rad/sec)
4
X: 4.034
Inflection Figure S14.17b Bode plots for the exact and approximate models. Figure S14.18 Bode plots for the exact and approximate models.
Y: 2.4
point
2
14.21
1 5 1
AR OL K c 1
(1.0)
25 2 100 1 1
2 2
Figure S14.19b Step response of closed loop system with different Kc.
Figure S14.19c Bode plot of FOPTD model.
(b)
Use Skogestad’s half rule
14.20
60 0.5 5 62.5
2.5 3 2 2 9.5
Using the Bode plot, at a phase angle of -180°, we require that K c KV K p K m 1
The approximated FOPTD model:
G p ( s) e s GV 0.5 Gm 1.0
0
10 10
1
0.5
Gp(s) = (1)
8s 2 7 s 1
AR/Kc
0
-1 10
10
0.5
-1 For R2 = 0.5, Gp(s) = (2)
-2
10
2s 2 5.8s 1
10
10
-2 a) For R2=2
-100
-150
7 c
|Gp| =
-150 -180
0.5
Gp= tan-1 2
,
Phase (deg)
2
Phase (deg)
1 8 c
-200
(1 8 c ) (7 c )
2 2
-200
-250
-250 -300
Kcu and c are obtained using Eqs. 14-7 and 14-8:
-300 -350
-2 -1 0 1 2
10 10 10 10 10
-350
Frequency (rad/sec)
7 c
10
-2
10
-1
10
0 1
10
2
10 Figure S14.21b Solution for part b) using Bode plot -180 = 0 + 0 + tan-1 2
tan-1(0.5c)
Frequency (rad/sec) 1 8 c
Bode Diagram
Figure S14.21a Bode plot 10
2
b) Set = 180 and solve for to obtain c = 0.4695 Solving, c = 1.369 rad/min
1
10
AR/Kc
1 ( K cu )(2.5)
Then AR OL c
= 1 = Kcu(1.025) 10 0.5 1.5
(1 8 c ) (7 c ) (0.5c ) 1
-1 2 2 2 2
10
For R2=0.5
Phase (deg)
-200
c) For Kc = 0.2, set AROL = 1 and solve for to obtain g = 0.1404 -250
5.8 c
|Gp| =
-300 0.5
Gp = tan-1 2
,
Then g = = -133.6 2
1 2 c (1 2c ) (5.8c )
-350 2 2
g -2 -1 0 1 2
10 10 10 10 10
Frequency (rad/sec)
From Eq. 14-61, PM = 180 + g = 46.4 Figure S14.21c Solution for part c) using Bode plot For Gv = Kv = 2.5, v=0, |Gv| = 2.5
1.5 1.5
d) From Eq. 14-60 14.22 For Gm = , m= -tan (0.5) ,
-1
|Gm| =
0.5s 1 (0.5c ) 2 1
1 1
GM = 1.7 = = From modifying the solution to the two tanks in Section 6.4, which have a
Ac AR OL
c slightly different configurations,
Substituting c = 2.51 rad/min, Kcu = 15.93, cKcu = 40.0 Substituting for g gives Kc = 1.34. Then from Eq. 14-8
Substituting into Eq. 14-57 gives g=|=g = 172.7
0.5c
b) From part a), for R2=2,
Kc = 0.45 Kcu = 4.932 , I= Pu/1.2 = 3.825 min Using Eq. 12-11, the PI settings are
From Eq. 14-60, GM = 1/Ac = 2.61
1
Using Eqs. 13.63 and 13-62 , Kc 1 , I = = 1 min,
K c
c) By using Simulink-MATLAB, these two control systems are compared for
c= -tan-1(-1/3.825) a unit step change in the set point.
Using Eq. 14-58 ,
2 1.4
1 part a)
|Gc| = 4.932 1 1 part b)
c
1
Then, from Eq. 14-56
/ 2
or c = = 7.85 rad/min 0.8
-1
1 -1 7 c
0.2
Output
-180 = tan + 0 + tan 2
tan-1(0.5c)
3.825c 1 8 c
0.6
c
c
2
c 1
Using Eq. 14-57 2
c
0
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
time
Ac = AROL|=c =
From Eq. 14-60, GM = 1/Ac = 3.93 Figure S14.23 Close-loop response for a unit step change in set point.
2
1
1 (2.5)
0.5 1.5
= 4.932
3.825c 2 b) Using Eq. 14-61, The controller designed in part a) (Direct Synthesis) provides better
(1 8c ) (7c ) (0.5c ) 1
2 2 2
performance giving a first-order response. Part b) controller yields a large
g = PM 180 = 140 = tan-1(-1/0.5g) 0.2g tan-1(g) overshoot.
= 0.7362
Solving, g = 3.04 rad/min
AR OL AR OL
14.24 From the plot of vs , = 0.144 Bode Diagram
Kc Kc g 10
2
1
a) Using Eqs. 14-56 and 14-57 Since AR OL g
=1 , Kc = = 6.94 0
10
0.144
AR/Kc
Ysp 4 2 1 2 0.4 c) From the plot of vs. , = -180 at c = 4.05 rad/min -2
Kc
10
AR OL (1.0)
Ym 0.01 2 1 0.25 2 1 25 2 1
AR OL AR OL -4
vs ,
10
From the plot of = 0.0326 -90
= tan-1(2) tan-1(0.1) tan-1(0.5) – (/2) tan-1(5) Kc Kc c
-135
Ac = AR OL c
= 0.326
Phase (deg)
Bode Diagram From Eq. 14-60, GM = 1/Ac = 3.07 -180
2
10
-225
2
Bode Diagram
10
0 -270
-2 -1 0 1 2
10 10 10 10 10 10
0
Frequency (rad/sec)
AR/Kc
10
AR/Kc
-4
-4 10
10 -90
14.25
-90
-135
Phase (deg)
-135 -180
Phase (deg)
a) Schematic diagram:
-180 -225
-270
-2 -1 0 1 2 TC
-225 10 10 10 10 10
Frequency (rad/sec)
Hot fluid
-270
Figure S14.24b Solution for part b) using Bode plot
10
-2 -1
10 10
0
10
1 2
10 TT
Frequency (rad/sec)
Figure S14.24a Bode plot
Cold fluid
b) Using Eq.14-61 Mixing Point Sensor
Kc = 0.45 Kcu = (0.45)(0.167) = 0.075 Kcu and c are obtained using Eqs. 14-56 and 14-57. Including the filter GF into -180= 0 tan-1(0.083c) tan-1(0.432c) tan-1(0.017c)
I = Pu/1.2 = 16 s/1.2 = 13.33 s these equations gives tan-1(0.024c)
Using Eq. 14-60, GM = 1/Ac = 1.54 By using MATLAB, the Nyquist diagram for this open-loop system is
1 ( K cu )
5.624 2
(0.083 ) 2 1 (0.432 ) 2 1 (0.017 ) 2 1
c c c Solving Eq. 14-56 for g gives Nyquist Diagram
1
0.12
AR OL =1 at g = 11.78 rad/min 0.5
(0.024 ) 2 1 g
c Substituting into Eq. 14-57 gives 0
Imaginary Axis
g
-1.5
Using Table 12.4, the Ziegler-Nichols PI settings are tan-1(0.017g) tan-1(0.024g) = -166.8
-2
Kc = 0.45 Kcu = 5.84 , I=Pu/1.2 = 0.288 min Using Eq. 14-61 -2.5
-3
PM = 180 + g = 13.2
-3.5
b) Using Eqs.14-39 and 14-40
-4
-3 -2.5 -2 -1.5 -1 -0.5 0
c = Gc = tan-1(-1/0.288)= -(/2) + tan-1(0.288) Real Axis
14.28 Figure S14.28a The Nyquist diagram for the open-loop system.
2
1
|Gc| = 5.84 1
0.288 a) From Exercise 14.28, b) Gain margin = GM =
1
AR c
Then, from Eq. 14-57, 5.264
Gv ( s )
0.083s 1 where ARc is the value of the open-loop amplitude ratio at the critical
2 frequency c. By using the Nyquist plot
- = (/2) + tan-1(0.288c) tan-1(0.083c) tan-1(0.432c) G p ( s)
(0.432s 1)(0.017s 1)
Nyquist Diagram
0.12
tan-1(0.017c) tan-1(0.024c) Gm ( s) 1
(0.024s 1) 0.5
1
The PI controller is Gc ( s ) 51 0
Solving, c = 15.11 rad/min. 0.3s
-0.5
Hence the closed-loop transfer function is
Using Eq. 14-56 -1
Imaginary Axis
GOL Gc Gv G p Gm -1.5
2
1 5.264
Ac AR OL 5.84 1 Rearranging,
-2
c
0.288c (0.083 ) 2 1
c -2.5
6.317s 21.06
2 0.12 GOL -3
1.46 10 5 s 5 0.00168s 4 0.05738s 3 0.556s 2 s
(0.432c ) 2 1 (0.017c ) 2 1 (0.024c ) 2 1 -3.5
-4
-3 -2.5 -2 -1.5 -1 -0.5 0
Real Axis