Colorado School of Mines CHEN403: FT Yt Fs Ys Gs
Colorado School of Mines CHEN403: FT Yt Fs Ys Gs
f (t )
Dynamic
Process
Input
y (t )
f (s)
Output
G(s)
y(s)
Consider the simple process with one input & one output. The describing n-th order ODE
is:
an
dn y
d n1 y
d2 y
dy
+
a
+
+
a
+ a1
+ a0 y = bf ( t )
n 1
2
n
n 1
2
dt
dt
dt
dt
d n1 y
d2 y
dy
=
=
=
= 0.
n 1
2
dt t =0
dt t =0 dt t =0
Taking the Laplace transform of this gives:
an s n y ( s ) + an1 s n1 y ( s ) + + a2 s 2 y ( s ) + a1 sy ( s ) + a0 y ( s ) = bf ( s )
y(s)
f (s)
an s + an1 s
n
Transfer Functions
n 1
b
G(s)
+ + a2 s 2 + a1 s + a0
-1-
All resistance to heat transfer is in a thin film around the bulb i.e., neglect thermal
resistance of glass & mercury.
All thermal capacity is in the mercury.
Mercury always has uniform temperature.
The glass wall does not expand or contract.
= hA (Ta T )
dt
dt
dt dt
dT
mC p
= hA (Ta T ) for constant C p
dt
mC p dT
= Ta T
hA dt
dT
T = Ta
dt
where the time constant is:
mC p
hA
At steady state:
T * = Ta*
dT
T = Ta where T ( 0) = 0
dt
Transfer Functions
-2-
( s + 1)T = Ta
so the transfer function is:
G(s) =
T
1
=
Ta s + 1
So, we would expect the heat transfer resistance around a thermometer to be a 1st order
system.
dT
+ T = Ta ( sT T ( 0 ) ) + T = Ta
dt
( s + 1)T = Ta + T (0)
1
Ta +
T ( 0)
s + 1
s + 1
1
=
Ta +
Ta*
s + 1
s + 1
T=
Now there are two inputs & two transfer functions: one for the driving function ( Ta ( t ) or
Ta ( s ) ) and one for the initial condition ( Ta* ).
Transfer Function for Process with Multiple Inputs and/or Multiple Outputs
What if there are multiple inputs and/or multiple outputs? We would associate a transfer
function with each pairing of an input & output. The block diagram for 2 inputs & 1 output
is:
f1 ( s )
G(s)
f1 ( s )
G1 ( s )
f2 ( s )
G2 ( s )
y(s)
f2 ( s )
+
+
y(s)
Transfer Functions
-3-
y ( s ) = G1 ( s ) f1 ( s ) + G2 ( s ) f2 ( s ) + G3 ( s ) f3 ( s ) + + Gn ( s ) fn ( s )
n
y ( s ) = Gi ( s ) fi ( s )
i =1
f1 ( s )
G1,1 ( s )
+
+
y1 ( s )
G2,1 ( s )
f2 ( s )
G2,1 ( s )
G2,2 ( s )
+
+
y2 ( s )
y1 ( s ) = G1,1 ( s ) f1 ( s ) + G1,2 ( s ) f2 ( s )
y2 ( s ) = G2,1 ( s ) f1 ( s ) + G2,2 ( s ) f2 ( s )
For n inputs and m outputs, then:
yi ( s ) = Gi ,1 ( s ) f1 ( s ) + Gi ,2 ( s ) f2 ( s ) + Gi ,3 ( s ) f3 ( s ) + + Gi ,n ( s ) fn ( s )
n
yi ( s ) = Gi , j ( s ) f j ( s ) for i = 1,2,3, , m .
j =1
y ( s) = G(s ) f ( s )
Transfer Functions
-4-
h, A, T,
Fs, Ts
F1, T1, 1
dt
dt
dt dt
So:
dH d VH
=
= F00 H 0 F1H + Q
dt
dt
If we assume that the enthalpy can be expressed as:
Transfer Functions
-5-
d
VC pT = F00C pT0 F1C pT + Q
dt
d
C p ( VT ) = F00C pT0 F1C pT + Q
dt
d
Q
( VT ) = F00T0 F1T +
dt
Cp
dh
= F0 F1
dt
and:
A
d
Q
( hT ) = F0T0 F1T +
dt
Cp
F
F
d
Q
( hT ) = 0 T0 1 T +
dt
A
A
AC p
TA
dh
dT
Q
+ hA
= F0T0 F1T +
dt
dt
C p
dT
Q
= F0T0 F1T +
dt
C p
dT
Q
= F0 (T0 T ) +
where h = h( t ) .
dt
C p
dT
Q
= F0 (T0 T ) +
dt
C p
If we are using steam for the heating medium, then we could relate the rate of heat added,
Q , to the steam temperature, Ts , as:
Transfer Functions
-6-
Q = UA (Ts T ) .
So:
UA (Ts T )
dT
= F0 (T0 T ) +
dt
C p
dT
UA
UA
+ F0 + T = F0T0 +
Ts
dt
C p
C p
F
dT F0
UA
UA
Ts
+ +
T = 0 T0 +
dt V VC p
V
VC p
dT 1
1
+ + K T = T0 + KTs
F
dt F
dT
1
+ aT = T0 + KTs
dt
F
where:
1 F0
UA
1
,K
, and a + K .
F V
F
VC p
At steady state:
aT * =
1 *
T0 + KTs*
F
so:
dT
1
+ aT = T0 + KTs
dt
F
where the deviation variables are defined as:
Transfer Functions
-7-
sT T ( 0) + aT =
1
T0 + KTs
F
1
T0 + KTs
F
1
( s + a )T = T0 + KTs
F
1 F
K
T =
T0 +
Ts
s +a
s +a
sT + aT =
T = G0 ( s )T0 + Gs ( s )Ts
where:
G0 ( s )
1/ F
K
and Gs ( s )
s +a
s +a
A block diagram for the stirred tank heater can be drawn as follows.
1/ F
s+a
T0 ( s )
+
+
T ( s )
K
s+a
Ts ( s )
G1(s)
y2(s)
G2(s)
Gn(s)
yn(s)
Transfer Functions
-8-
yn ( s )
f (s)
= Gi ( s )
i =1
y(s)
f (s)
= G(s)
where:
G(s) =
Q(s)
P (s)
and where Q ( s ) and P ( s ) are usually polynomials in s (time delays will introduce
exponential terms, however). In general, the order of Q ( s ) will be less than that of P ( s ) .
The roots of the numerator Q ( s ) are referred to as the zeros of the transfer function. At
the zeros, G ( s ) becomes zero. The roots of the denominator P ( s ) are referred to as the
poles of the transfer function. At the poles, G ( s ) becomes infinite.
We can get a qualitative sense of the response of a system by knowing the poles. Let:
f (s) =
r (s)
q(s)
Since:
G(s) =
Q(s)
P( s )
then:
Transfer Functions
-9-
Q(s) r (s)
P (s) q(s)
i =1
i =1
P ( s ) = ( s pi ) ( s pi )
mi
and:
y(s) =
Q(s)
N
r (s)
(s p ) (s p )
i
i =1
q(s)
mi
i =1
When split into partial fractions, each of the factors in the denominator will lead to a
separate term. Splitting up the factors of the transfer function (while leaving the
denominator from the input function aside for now) gives:
mi 1
N
C
y(s) = i +
i =1 s pi
i =1
D (s p )
j ,i
j =0
( s pi )
mi
M mi 1
D j ,i
Ci
y(s) =
+
mi j
i =1 s pi
i =1 j =0 ( s pi )
r* (s)
q( s)
r* (s)
+
.
q(s)
Note that for the repeated root, the numerator can be a polynomial of order up to one less
the order of denominator. Also, each repeated root can have a different order. The only
requirement on the number of roots is that they have to add up to n , i.e.:
M
n = N + mi .
i =1
N C N
L1 i = C i exp ( pi t ) .
i =1 s pi i =1
Transfer Functions
- 10 -
i =1 j =0 ( s pi )
M mi 1
D j ,i
= L1
mi j
i =1 j =0
s
p
(
)
i
M
mi 1
D
= exp ( pi t ) L1 mij,i j
i =1
s
j =0
mi 1
= exp ( pi t )
j =0 ( mi
i =1
M
mi 1
= exp ( pi t )
j =0 ( m
i =1
i
( m j 1)!
L1 i mi j
j 1 )!
s
D j ,i
t mi j 1
j 1 )!
D j ,i
Note that the roots pi are important for the long-time characteristics of the solution. For
the real non-repeating roots:
For the complex non-repeating roots (which will occur as a complex conjugate pair), then
pi can be expressed as i i i . These roots will give rise to terms of the form
exp ( i t ) sin ( i t + i ) . Now, the important term with regards to stability is the real portion
of the root, i :
For the repeating roots, the situation is similar. The polynomial term will always grow
towards infinity as t , so the behavior of the exponential term will dictate the overall
behavior.
If pi < 0 or i < 0 , then the exponential term will go to zero as t and the entire
term will also go to zero. This exponential decay leads to a zero contribution from
this pole.
Transfer Functions
- 11 -
If pi > 0 or i > 0 , then the exponential term will grow to infinity as t and the
entire term will also grow to infinity. This exponential growth leads to an explosive
contribution from this pole.
If pi = 0 or i = 0 , then the polynomial term will dictate the behavior for t .
This polynomial term will lead to an explosive contribution from this pole..
So, in general:
G(s) =
Q(s)
P (s)
Q(s)
s + 3s + 5s 2 + 4s + 2
4
Transfer Functions
- 12 -
P (s ) 5
4
3
2
1
0
-2.5
-2
-1.5
-1
-0.5
0.5
1.5
2.5
P ( s ) = s 4 + 3s3 + 5s 2 + 4s + 2 = ( s 2 + s + 1 )( s2 + 2s + 2)
From this, we find that the roots are:
r=
1 1 4 1 1
1
3
i
=
2
2 2
r=
2 22 4 1 2
= 1 i
2
Since the real portion of the roots are all negative, the system is stable.
Transfer Functions
- 13 -