Process Control
Process Control
Process Control
Advanced
Process Control
Prof. Shi-Shang Jang
Chemical Engineering Department
National Tsing-Hua University
Hsin Chu, Taiwan
Course Description
Course: CHE517 Advanced Process Control
Instructor: Professor Shi-Shang Jang
Text: Seborg, D.E., Process Dynamics and
Control, 2nd Ed., Wiley, USA, 2003.
Course Objective: To study the application of
advanced control methods to chemical and
electronic manufacturing processes
Course Policies: One Exam(40%), a final
project (30%) and biweekly homework(30%)
Course Outline
1. Review of Feedback Control System
2. Dynamic Simulation Using MATLAB
and Simu-link
3. Feedforward Control and Cascade
Control
4. Selective Control System
5. Time Delay Compensation
6. Multivariable Control
Chapter 1 Review of
Feedback Control Systems
Feedback Control
Terminology
Modeling
Transfer Functions
P, PI, PID Controllers
Block Diagram Analysis
Stability
Frequency Response
Stability in Frequency Domain
Feedback Control
Transmitter
Controller
Set point
Temp
sensor
Heat loss
stream
condensate
Examples:
Room temperature control
Automatic cruise control
Steering an automobile
Supply and demand of chemical engineers
Feedback Control-block
diagram
error
+
Set point
Manipulated
variable
Controller
disturbance
process
Controlled variable
Measured value
Sensor +
transmitter
Terminology:
Set point
Manipulated variable (MV)
Controlled variable (CV)
Disturbance or load (DV)
Process
controller
Instrumentation
Transmitter
Controller
Set point
Temp
sensor
Heat loss
stream
condensate
Modeling
Q
Mass M
Cp
T
Q=UA(T-T0)
d
( MC PT ) Q UA(T T0 )
dt
Transfer Functions
Laplace Transforming:
M Cp S Td(S) = qd(S) - U A (Td(S) Tod(S))
Or
q d S
UATod S
Td S
SMC p UA SMC p UA
qd(S)
UA
MsC p UA
1
Tod(S)
MsC p UA
Td(S)
Non-isothermal CSTR
F0
0
CA0
T0
steam
AB
rA = - KCA mol/ft3
T V C A CB
condensate
CA
T
K = e-E/RT
Mass balance :
d
(V C A ) F0 C A0 F C A ( KC A )V
dt
d
Energy balance : (V C PT ) F0C PT0 F C PT (Hr )( KC A ) UA(Ts T )
dt
Initial conditions :
Input variables :
Linearization of a Function
F(X)
aX+b
X0 -
X0
X0+
Linearization
dx
f
f x, u
dt
x
x x0
u u0
x x0
f ( x0 , u0 )
dxd
axd bud 0
dt
Laplace Transform
sxd s axd s bud s
or
xd ( s )
b
K
ud s s a s 1
f
u
x x0
u u0
u u0
C
A
,
d
dt
Td
21
a31
0
a22
a32
0
0
a 33
Vd
b11 b12
C
b b
A
,
d
21 22
Td
b31 b32
0
b23
b33
F0,d
Vd
y 0 0 1 C A,d 0 0 0 Fd CX DU
Ts ,d
Td
Td s C ( sI A) 1 B D U (s )
G p s Tsd s GL s F0,d (s ) GL ' s Fd s
F0,d
F
d
AX ( s ) BU ( s )
Ts ,d
CV ( s )
K
MV s s 1
CV ( s )
K
2 2
MV s s 2s 1
CV ( s )
K Ds
e
MV s s 1
Transfer Functions of
Controllers
Proportional Control (P)
m(s) = Kc[ e(s) ]
e(s)
e = Tspt - T
m(s)
Kc
1
m ( t ) K c e( t )
I
e(t )dt
0
1
m(s) K c e(s)
e(s)
s
I
e(s)
K c (1
1
)
Is
m(s)
1
m ( t ) K c e( t )
I
m(s) K c
de
0 e(t )dt D dt
t
1
e(s) 1 Ds
Is
e(s)
K c (1
1
D s)
Is
m(s)
Purdy oscillatory
Exponential growth
with oscillatory
Fast Exponential
growth
Exponential Decay
Fast Decay
Slow Decay
Slow growth
Purdy oscillatory
Stable (LHP)
Unstable (RHP)
Re
1/
U1(S)
X2(S)
G1(S)
X1(S)
G2(S)
U2(S)
X2(S)
X (S)
Kc
QS
control
1
MC P S UA
process
Measuring device
1
Proportional control
No measurement lags
Td
Tset
1
MC P S UA
1
1 KC
MC P S UA
KC
Td
Gc(S)
L(S)
GL(S)
GP(S)
Xm
X1
+
+
Gm(S)
e = Xs Xm
m = Gc (S) e(s) = Gc e
X1 = Gp m = Gp Gc e
X = GL L + X1 = GL L + Gp Gc e
Xm = Gm X = Gm GL L + Gp Gc e
X = GL L + Gp Gc[Xs Xm]
= GL L + Gp Gc [Xs] Gp Gc [Xm]
=GL L + Gp Gc Xs Gp Gc Gm X
G pG c
GL
X
L
Xs
1 G pG cG m
1 G pG cG m
X(S)
Level System
dh
Fin Fout Fin k h
dt
Given a reference point Fin 0 , h0
A
dhd
f
Fin Fin 0 f h h0 Fin,d k hd
dt
Fin
h
2 h0
Laplacing
or
hd s
1
1/ a
K
Fin.d s As a A / a s 1 s 1
TRC
FC
2A B
Tc
Wc
T, Ca
dC A W (C A f C A )
2
kC A
dt
V
(ii) Energy balance of the jacket
dt
M cCP '
Mc
(iii) Energy balance for the reactor
2
dT W (T f T ) A(T Tc ) HkC A
dt
V
C PV
C P
Q
k A0 exp(
)
T 273
Linearization of
Nonisothermal CSTR
CV=T(t)
MV=Wc(t)
It can be shown that
Td s
K
3
Wc ,d s as bs 2 cs 1
Maximum slope
C
D
Dead time
Time constant
time
PI
PID
Kc
/DKp
0.9/DKp 1.2/DKp
n.a.
D/0.3
D/0.5
n.a.
n.a.
0.5D
D= 1
=13
k = -0.0202
Kc= -579.2079
i =3.33
setpoint
Ku/2
Ku/2.2
Ku/1.7
n.a.
Pu/1.2
Pu/2
n.a.
n.a.
Pu/8
Measuring Controller
Performance
IAE ys y t dt e t dt
ISE y s y t dt e t dt
2
ITAE t y s y dt t e t dt
Approaches
Direct Substitution for Kc
Root Locus method for Kc
Frequency Analysis for all
parameters
An Example
Kc
1
( s 1)( s 2)( s 3)
6+Ku-611=0 Ku=60
3. Frequency Domain
Analysis
Definitions: Given a transfer function
G(s)=y(s)/x(s); Given x(t)=Asint;
we have y(t) Bsin(t+)
We denote Amplitude Ratio=AR()
=B/A; Phase Angle=()
Both AR and are function of
frequency ; we hence define AR
and is the frequency response of
system G(s)
An Example
A sin(t)
1
s 1 s 2 s 3
B = sin(t+)
Frequency Response of a
first order system
y(s)
K
A
G( s)
; x(t ) A sin t x( s )
x( s )
s 1
s
A
K
y(s)
s s 1
KA
y (t )
sin(t ); tan
1
K
AR
1
tan
2
Basic Theorem
Given a process with transfer function
G(s);
AR()= G(i)
()= G(i)
Basically, G(i)=a+ib
AR a 2 b 2
tan b / a
1
i
a ib
2 2
2 2
2 2
i 1 1
1
1
1
2
2
AR a b
1 2 2
b
tan 1 tan 1
a
Note that
lim AR 0
lim 90
Corollary
If G(s)=G1(s)G2(s)G3(s)
Then AR(G)=AR(G1) AR(G2) AR(G3)
(G)= (G1) + (G2)+ (G3)
Proof: Omitted
Example
K1
K2
G ( s)
G1 s G2 s
1 1 2 1
AR AR1 AR2
K1
1 1
2
K2
2
2 1
1 2 tan 1 1 tan 1 2
Nyquist Plot
sys=tf(num,den)
NYQUIST(sys,{wmin,wmax}))
Simulink Example
Response
0.165 0.5 s
GP
e
2.5s 1
D1.4
3.7-1.4=2.3
time
Transfer function:
1
---------------------s^3 + 6 s^2 + 11 s + 6
u=3.5
>> bode(sys)
ARu=-38db
=10-38/20
=0.0162
Ku=1/ARu=80
2
1.8
1.6
1.4
1.2
1
0.8
0.6
0.4
0.2
0
10