EJERCICIOS

1. Se tiene un tanque en donde se almacenará un líquido. El tanque tiene

un flujo a la entrada y un flujo a la salida. Manipularemos la válvula de
la salida para mantenerla dentro de control.

Fo

h́=4 m
D́=2m

h −2 m3
F́ o=1 x 10
s
F

Vamos a establecer las variables

Variable controlada VC → h

Variable manipulada VM → F

Variable de perturbación VP → F o

Realizaremos un balance global

Acum=entrada−salida
dm
=m entrada−m salida ; m= ρv ; asumiento T es constante( ρ constante)
dt
ṁ=ρF
dρV
=ρ F o− ρF
dt
dV
ρ =ρ ( Fo −F )
dt
dV
=F o−F
dt
Como el: V = Ah
Adh
=F o−F
dt

La corriente de salida presenta una válvula, que la usaremos para mantener el

nivel de líquido deseado

∆P
F=C V
√ ρ

( gρ h1+ P atm ) −(gρ h2 + Patm ) gρ h1−gρ h2

CV
√ ρ
=C v
√ ρ
=C v
gρh
ρ √
F=∝ √ h
Sustituyendo la ecuación de la válvula en la ecuación obtenida del balance:
dh
Á =F o−ά √h
dt

Observamos que el ultimo termino no está lineal, por lo que lo linealizaremos

haciendo uso de la serie de Taylor para una variable.
dF
f ( x)=F ( x́ ) + |
dx x́
( x− x́ )

ά
f (∝ √h )=2 √ h́+ ( h−h́ )
2 √ h́
Sustituyendo esta la anterior, tenemos
dh ά
Á =F o−ά √h́− ( h−h́ )
dt 2 √ h́
Ahora procedemos a llevar las variables que estén en el tiempo a condiciones en
desviación.

X d= X− X́

d (hd + h́) ά ( d )
Á =(Fod + F́ o )−ά √ h́− h
dt 2 √ h́
Haremos cero todo lo que este en estado estacionario que no esté multiplicando
una desviación.
 d hd ά ( d )
Á =F do− h
dt 2 √ h́
Trataremos de darle la forma de un modelo matemático de primer orden, primera
derivada más la variable.

d hd ά ( d )
( Á +
dt 2 √ h́
h =F do
)
2 √ h́
Para que la variable este sola, multiplicaremos toda la expresión por:
ά
Obteniendo:

2 √ h́ Á d hd d 2 √ h́ d
+h = Fo
ά dt ά
Donde la constante de tiempo será:

2 √ h́ Á
T p=
ά
Como:
2
π D2 π ( 2 m)
A= = =3.1416 m2
4 4

F 1 x 10−2
F=∝ √ h , despejando ∝nos queda ,∝= = =5 x 10−3
√ h́ √ 4
Entonces:

2 √ h́ Á 2 √ 4(3.1416)
T p= = =2 513.28 s
ά 5 x 10−3
La ganancia o sensibilidad será:

2 √ h́ d 2 4 m
K p= F o= √ −3 =800 3
ά 5 x 10 m
s
El proceso se estabilizará para:
T =5 Tp=5 ( 2 513.28 )=12566.4 s
Finalmente podemos decir que el modelo matemático es:

d hd d d
2 513.28 +h =800 F o
dt
el diagrama de bloque, será:
F od hd
d hd d d
2 513.28 +h =800 F o
dt

La Función Transferencia de:

d hd d d
2 513.28 +h =800 F o
dt
Aplicando transformada de Laplace
2 513.28 S h(s) +h(s )=800 Fo(s )

Sacando factor común la variable de control

h( s ) ( 2513.28 S+1 ) =800 Fo(s)

Podemos tener la función transferencia:

h(s) 800
=
Fo(s) 2 513.28 S+1

El Diagrama de bloque con la función transferencia será:

Fo(s) h(s)
800
2513.28 S+1

¿Cuál será el valor final de la variable de control si F́o varía a 3 x 10−2 m3 /s ?

Por ser un proceso de primer orden aplicaremos la solución general.
−T
=Kp Q ( 1−e )
d d Tp
X G

3 −5 Tp
d m ( −2
h =800 3 3 x 10 −1 x 10−2 m
)
s
(
1−e Tp )=31.99 m
m
s
Y como la pregunta es el valor final, entonces:

h=hd + h́=31.99+ 4=35.99 m

2.- Se tiene un tanque en donde entra un componente A, como se desea
mantener en un valor desea, entonces usaremos el tanque superior en
donde tenemos al componente puro.

A
Puro
C AO C´ A=0.89lbmolA / f t 3
Fo
C´AO=0.1 lbmolA /f t 3

F́ O=80 f t 3 /min

f t3
F´AP=1
min
CA V́ =3000 f t 3
F

Variable:

VC →C A

VM → F AP

VP →C AO

Balance global.
Acum=entrada−salida
dm
=m entrada−m salida
dt
Donde m=ρV ; ( asumiendo T es constante ρ constante )
ṁ=ρF
dρV dV dv
=ρ F o− ρF + ρ F AP =ρ = ρ ( F o−F + F AP )= =F o−F+ F AP
dt dt dt

dV dv
ρ =ρ ( Fo −F+ F AP )= =F o−F + F AP
dt dt
dV
= F́O − F́+ F AP
dt
Asumiendo que el volumen es constante, y como la derivada de una constante es
cero, se tiene:
Se tiene:
F́= F́ O+ F AP

Haciendo un balance por componente

dn
=n entrada−n salída
dt
Donde:
n=CV
Tenemos:
dCA
V́ =C AO F O +C AP F AP−C A F
dt
Sustituyendo el balance global dentro del balance por componente, se tiene:
dCA
V́ =C AO F́ O + C´AP F AP−C A ( F́O + F AP )
dt

dCA
V́ =C F́ O + C´AP F AP−C A F´O−C A F AP
dt AO

Podemos observar que el último término no está lineal, entonces debemos

linealizar usando la serie de Taylor para dos variables.
dF dF
f ( x , x )=f ( X́ 1 , X́ 2 ) +
1 2
d X1
(|X 1− X́ 1 ) + ( X − X́ 2)
d X2 2

F(C A
=C A ´F AP + F´AP ( C A −C´ A ) + C´ A (F AP− F´AP )
, F AP )

Sustituyéndolo en el balance,
dCA
V́ = F́ O C AO + C´AP F AP− F́O C A− C A ´F AP− F´AP ( C A −C´ A )−C´ A ( F AP− F´AP)
dt
Llevando a condiciones en desviación,

X d= X− X́ , donde , X=X d + X́
 (CdA +C ¿¿ A) d d d d
V́ d = F́O ( C AO + C´AO ) + C´AP ( F AP + F´AP ) − F́ O ( C A +C A )−C A ´F AP− F´AP C A ¿
dt

−C´ A F dAP

Se hacen cero las variables estacionarias que no multiplican a una desviación.

d C dA d d d d d
V́ = F́ O C AO + C´AP F AP− F́O C A− F´AP C A −C´ A F AP
dt
Agrupando términos:

d C dA d d d
V́ = F́ O C AO + ( C´AP−C´ A ) F AP −( F´O + F´AP) C A
dt
Llevando a la forma de la ecuación de un proceso de primer orden.

d C dA 1
( V́
dt )
+ ( F́ O + F´AP ) C dA= F´O C dAO + ( C´AP −C´ A ) F dAP x
( F´O + F´AP )

V́ d C dA d F´O d C´AP−C´ A d
+C A = C AO + F
F´O + F´AP dt F́O + F´AP F́O + F´AP AP

Donde:

V́ 3000 f t 3
T P= = =37.04 min
F́ O + F´AP f t3 f t3
80 +1
min min

f t3 lbmol
80
F´O min f t3
Kp 1= = 3
=0.9876
F́O + F AP
´ ft f t3 lbmol
80 +1
min min f t3

Para el cálculo de Kp2 se necesita tener el valor de C AP del cual se obtendrá

mediante un balance de componentes.
dV
=C AO F O +C AP F AP−C A F
dt
en condición estacionaria
C A F−C AO F O ( 80 ) ( 0.89 )−( 80)( 0.1) lbmol
C AP= = =63.2
F AP 1 f t3
lbmol lbmol lbmol
63.2 −0.89
C −C A
´ ´ ft 3
f t3 f t3
Kp2= AP = =0.7802
F́ O + F´AP f t3 f t3 f t3
80 +1
min min min
Se estabilizará para
T =5 Tp=5 ( 37.04 )=185.2 min
Finalmente, el modelo matemático es:

d C dA d d d
37.04 +C A=0.9876 C AO +0.7802 F AP
dt
Diagrama de bloque del modelo matemático.

C dAO d C dA d d
37.04 +C A=0.9876 C AO
dt
+
C dA
+
F dAP d C dA d d
37.04 +C A=0.7802 F AP
dt

La Función transferencia.
37.04 S C A (s )+ C A (s )=0.9876 C A 0(S )+ 0.7802 F AP (S )

( 37.04 S+1 ) C A (s )=0.9876 C A 0(S ) +0.7802 F AP (S )

C A ( s) 0.9876 C () 0.7802
= ; As =
C A 0(S ) 37.04 S+1 F AP( S ) 37.04 S+1

El Diagrama de bloque con función transferencia.

C A 0 (S ) 0.9876
37.04 S +1
+
C A ( s)
+
F AP( S ) 0.7802
37.04 S +1
¿Cuál será la desviación de la variable de control si C AO varía 0.05 lbmolA /f t 3 y
F AP varía a 1.1 f t 3 /min?
−T
X d
G =Kp∗Q ( 1−e )
d Tp

lbmolA
C dA= [0.9876 ( 0.05 )+ 0.7802 ( 1.1 ) ] ( 1−e−5 ) =0.1265
f t3

3.- el mismo ejercicio 2, pero ahora cambiamos la variable perturbadora

A
C AO Puro C´ A=0.89lbmolA / f t 3
Fo
C´AO=0.1 lbmolA /f t 3

F́ O=80 f t 3 /min

F´AP=1 f t 3

V́ =3000 f t 3
CA
F

Variables
VC →C A

VM → F AP

VP → F O

Balance global.
Acum=entrada−salida
dm
=m entrada−m salida + m AP ;
dt
Donde m=ρv ; ( asumiendo T es constante )
dρV
=ρ F o− ρF + ρ F AP
dt
 dV
ρ =ρ ( Fo −F+ F AP )
dt
dV
=F o−F+ F AP
dt
Asumiendo que el volumen no varía,
F=F o+ F AP

Balance por componentes.

dn
=n entrada−n salída
dt
n=CV
dCA
V́ =C AO F O +C AP F AP−C A F
dt
Sustituyendo el balance global dentro del balance por componente.
dCA
V́ =C´AO F O + C´AP F AP−C A (F o+ F AP )
dt

dCA
V́ =C´AO F O + C´AP F AP−C A F o−C A F AP
dt

Linealizando. Serie de Taylor para dos variables, los dos últimos términos.
dF dF
f ( x , x )=f ( X́ 1 , X́ 2 ) +
1 2
d X1 1|
( X − X́ 1 ) + d X ( X 2− X́ 2)
2

F(C A , F O)= C A´F O + C´ A ( FO − F́O ) + F́O (C A −C´ A )

F(C A , F AP ) =C A ´F AP + F´AP ( C A −C´ A ) + C´ A (F AP− F´AP )

Sustituyendo:

dCA
V́ =C´AO F O + C´AP F AP−C´A FO −C´ A ( FO − F́O )− F́ O ( C A −C´ A ) −C A´F AP− F´AP ( C A −C´ A ) −C´ A (F AP − F´AP )
dt
Se aplica lacondición de desviación

X d= X− X́ , donde , X=X d + X́
 (CdA + Ć ¿¿ A) ´ d d d d d d
V́ d =C AO ( F O + F́ O )+ C´AP (F AP+ F´AP )−C´A FO −C´ A F O − F́O C A−C A ´F AP − F´AP C A −C´ A F AP ¿
dt

Se hacen cero las condiciones estacionarias

Aplicando condición de desviación y haciendo ‘’0’’ cero lo estacionario que no
multiplica lo que está en desviación nos queda:

d C dA d d d d d d
V́ =C´AO F O + C´AP F AP−C´ A F AP + F´AP C A + C´ A F O + F́ O C A
dt
Agrupando:
¿

V́ d C dA d ( C´AO−C´ A ) d C´AP −C´ A d

+C A = FO+ ´ F AP
F´AP + F́ O dt F´AP + F́ O F AP+ F́ O

Calculando:

V́ 3000 f t 3
Tp= = =37.04 min ,
F´AP+ F́ O 80 f t 3 1 f t 3
+
min min

( C´AO−C´ A ) ( 0.1−0.89 ) lbmolA / f t 3 llbmolA / f t 3

Kp 1= = =−0.0098
F´AP + F́ O 80 f t 3 1 f t 3 f t 3 /min
+
min min
C´AP −C´ A
Kp 2= ; Del balance por componente del ejercico2 obtuvimos C AP resultando :
F´AP+ F́ O

lbmol lbmol
63.02 3
−0.89
ft f t3 llbmolA/ f t 3
Kp2= =0.7802
80 f t 3 /min +1 f t 3 f t 3 /min

T =5 Tp=5 ( 37.04 )=185.2 min

Modelo matemático.
 d C dA d d d
37.04 +C A=−0.0098 F O +0.7802 F AP
dt

Diagrama de bloque del modelo matemático.

F dO d C dA d d
37.04 +C A=−0.0098 F O
dt
+
C dA
+
F dAP d C dA d d
37.04 +C A=0.7802 F AP
dt

Función transferencia.
37.04 S C A (s )+ C A (s )=−0.0098 F0 (S) +0.7802 F AP(S )

( 37.04 S+1 ) C A (s )=−0.0098 F 0(S )+ 0.7802 F AP(S )

C A ( s) −0.0098 C A ( s) 0.7802
= ; =
F 0 (S) 37.04 S +1 FdAP ( s) 37.04 S +1

Diagrama de bloque con la función transferencia.

F 0(S ) −0.0098
37.04 S +1
+
C dA
+
F AP(S ) 0.7802
37.04 S +1

¿Cuál será la desviación de la variable de control si F 0(S )varía 0.05 ft 3 /min F AP

varía 1.1 f t 3 /min?
−T
(
X Gd =Kp∗Q d 1−e Tp )
llbmolA / f t 3 llbmolA / f t 3
[
C dA= −0.0098 3
f t /min
( 0.05 ) f t 3
/min+ 0.7802 3
f t /min ]
( 1.1 ) f t 3 /min ( 1−e−5 )
 lbmolA
C dA=0.8519
f t3
4.- repetir el problema 2, pero ahora las dos variables de entrada son de
perturbación.
A
Puro

C AO
Fo
C´ A=0.89lbmolA / f t 3

C´AO=0.1 lbmolA /f t 3

F́ O=80 f t 3 /min

F´AP=1 f t 3
CA
V́ =3000 f t 3
F

Variables
VC →C A

VM → F AP

VP → F O ,C AO

Balance global.
Acum=entrada−salida
dm
=m entrada−m salida + m AP ;
dt
m=ρv ; ( asumiendo T es constante )
dρV
=ρ F o− ρF + ρ F AP
dt
dV
ρ =ρ ( Fo −F+ F AP )
dt
dV
=F o−F+ F AP
dt
Resultando :F =F o + F AP

Balance por componentes.

dn
=n entrada−n salída
dt
n=CV

dCA
=C AO F O +C AP F AP−C A F
dt
sustituyendo balance global dentro del balance por componente
dCA
V́ =C AO F o + C´AP F AP −C A (F o+ F AP )
dt

dCA
V́ =C AO F́ o + C´AP F AP −C A F o−C A F AP
dt

Linealizando. Serie de Taylor

dF dF
f ( x , x )=f ( X́ 1 , X́ 2 ) +
1 2
d X1 1|
( X − X́ 1 ) + d X ( X 2− X́ 2)
2

F(C AO , FO ) =C´AO F́O + C´AO ( FO − F́ O ) + F́ O ( C AO −C´AO )

F(C A , F O)= C A´F O + C´ A ( FO − F́O ) + F́O (C A −C´ A )

F(C A , F AP ) =C A ´F AP + F´AP ( C A −C´ A ) + C´ A (F AP− F´AP )

Sustituyendo:
dCA
V́ =C´AO F́ O + C´AO ( F O− F´O ) + F´O ( C AO−C´AO ) + C´ A F AP−C A´F O−C´ A ( F O− F́ O ) −F´O ( C A− C´ A )−C A ´F AP −C´ A ( F
dt

Se aplica la condición de desviación : X d =X − X́ , donde , X =X d + X́

(CdA + Ć ¿¿ A) ´
=C AO F́O ¿+ C´AO FO + F O ´C AO + C´ A (F AP + F AP )−C A´FO −C´ A F O − F́O C A−C A ´F AP−C´ A F AP− F´AP
d d d d d d
V́ d
dt

Hacer cero lo estacionario

 d C dA
=+ C´AO F O + FO ´C AO + C´ A F AP−C´ A FO − F́ O C A −C´ A F AP− F´AP C A
d d d d d d d
V́
dt

V́ d C dA d C AO −C A d C AP−C´ A d FO d
+C A = FO + F AP + C AO
F O + F AP
´ dt F́O + F AP F O + F AP
´ F́ O + F AP

V́ 3000 f t 3
Tp= = =37.03 min
F́ O + F AP ( 80+1 ) f t 3 /min

C −C A ( 0.89−0.1 ) lbmolA /f t 3 lbmolA /f t 3

Kp 1= ´AO = =0.0097
F O + F AP ( 80+1 ) f t 3 /min f t 3 / min
lbmolA lbmolA
( 63.2−0.89 )
C AP −C´ A ft 3
f t3
Kp2= = =0.7693 3
F́ O + F AP ( 80+1 ) f t 3 /min f t /min
FO 80 f t 3 /min lbmolA / f t 3
Kp 3= = =0.9877
F́O + F AP ( 80+1 ) f t 3 /min lbmolA / f t 3

Tiempo de estabilización
T =5 Tp=5 ( 37.04 )=185.2 min

Modelo matemático.

d C dA d d d d
37.03 + C A =0.0097 F O+ 0.7693 F AP +0.9877 C AO
dt

Diagrama de bloque del modelo matemático.

F dO
d C dA d d
37.03 + C A =0.0097 F O
dt

F dAP d C dA d d
+
C dA
37.04 +C A=0.7693 F AP +
dt +

C dAO
 d C dA d d
37.04 +C A=0.9877 C AO
dt

Función transferencia.

37.04 S C A (s )+ C A (s )=0.0097 F0 (S) + 0.7693 F AP(S )+ 0.9877 C AO (S)

( 37.04 S+1 ) C A (s )=0.0097 F 0( S )+ 0.7693 F AP (S ) +0.9877 C AO ( S )

CA(s) 0.9877 C () 0.7693 C A (s ) 0.0097

= ; As = ; =
C AO ( S ) 37.04 S+1 F AP ( S ) 37.04 S+1 F 0 ( S) 37.04 S+1

Diagrama de bloque de la función transferencia.

F 0 (S ) 0.0097
37.04 S +1

F AP( S ) 0.7693 +
C dA (s)
37.04 S +1 +
+

C AO (S ) 0.9877
37.04 S +1

¿Cuál será la desviación de la variable de control si F O varía a 53 f t 3 /min y F AP

varía a 0.05 f t 3 /min?
−T
(
X Gd =Kp∗Q d 1−e Tp )
lbmolA
C dA= [0.0097 ( 53−80 ) +0.7693 ( 0.05−1.1 ) ] ( 1−e−5 )=−1.0625
f t3
5.- SE TIENE UN SISTEMA DONDE SE DESEA CONTROLAR LA
TEMPERATURA, CON LA YUDA DE UN FLUJO DE VAPOR QUE CIRCULA
POR UN SERPENTÍN
.F
O

F́ O=50 lb/ min

TO
T´O=50 ℉
mV
F T́ =120 ℉
mv =3.5 lb/min

Ḿ =2 000 lb
λ=1000 btu/lb
Variables

VC →T
VP →T O

VM → mV

Balance de materia
Acum=entrada−salida
dm
=m entrada−m salida ; m= ρv ;( asumir T es constante)
dt

dρV
=ρ F o− ρF
dt
dV
ρ =ρ ( Fo −F )
dt
dV
=F o−F
dt

Asumir volumen constante

Resultando :F =F o
Balance de energía.
dE ^ ^
= E E + E S +Q
dt

^ H
E=H donde , E= ^

dH ^ ^
= H entrada + H salida +Q ; H
^ =CpF ( T −T referencia )
dt

Q= λ mv , donde , H =CpMT , sustutiyendo nos queda :

dCpMT
=Cp FO ( T O −T referencia ) −CpF ( T −T referencia ) + λm v
dt
dT 1
[ CpM
dt
=Cp FO T O−CpFT + λm v x
Cp ]
dT λ́ m v
Ḿ = F́O T O−FT + ,
dt Cp
Sustituyendo el balance global dentro del balance de energía

dT λ́ m v
Ḿ = F́O T O + F́ O T +
dt Cp

Se aplica lacondición de desviación : X d =X − X́

d (T d −T́ ) λ́
Ḿ − F́O ( T dO −T́ O )− F´O ( T d−T́ ) + (mdV + ḿV )
dt Cp

haciendo ‘cero lo estacionario que no multiplica lo que está en desviación nos

queda:

d Td d d λ́ d
Ḿ = F́O T O− F´O T + m
dt Cp V
 d T d ´ d ´ d λ́ d 1
( Ḿ
dt
+ FO T = FO TO+
Cp
mV x)F´O

¿ Ḿ d T d d λ́
+T =T dO + mdV
F´O dt Cp F́O

Ḿ 2000 lb
Tp= = =40 min
F́ O 50 lb/ min
℉
Kp 1=1
℉
λ́ 1000 btu/lb ℉
Kp2= = =20
Cp F́ O 1 btu ∗50℉ lb /min
lb

Tiempo de estabilización
T=5Tp=5∗40 min=200 min

Modelo Matemático.

d Td d d d
40 +T =T O + 20 mV
dt

Diagrama de bloque del modelo matemático.

T dO d Td d d
40 +T =T O
dt
+
Td
+
mdV d Td d d
40 +T =20 mV
dt

Función Transferencia.
 d Td d d d
40 +T =T O + 20 mV
dt

40 S T (s )+T ( s)=¿ (s )+ 20 m v ( s)

T (s ) ( 40 S+1 ) =¿( s) +20 m v (s )

T ( s) 1 T ( s) 20
= ; =
¿ (s ) 40 S+1 mv (s ) 40 S+1

Diagrama de la función transferencia.

¿( s ) 1
40 S +1
+
T (s )
+
m v ( s) 20
40 S +1

6.- el mismo ejercicio 5 pero ahora la variable perturvadora el el flujo.

FO
F́ O=50 lb/ min

TO T´O=50 ℉

mV T́ =120 ℉
F mv =3.5 lb/min

Ḿ =2 000 lb
λ=1000 btu/lb
Variables
VC →T
VP → F O

VM → m V
Balance de materia
Acum=entrada−salida
dm
=m entrada−m salida
dt
m=ρv ; ¿
dρV
=ρ F o− ρF
dt
dV
ρ =ρ ( Fo −F )
dt
dV
=F o−F
dt
Resultando :F =F o

Balance de energía.
dE ^ ^
= E E + E S +Q ; E=H donde , ^
E=^
H
dt
dH ^ ^
= H entrada + H salida +Q ; H
^ =CpF ( T −T referencia ) ; Q=λ m v ,donde , H =CpMT , sustutiyendo nos queda:
dt
dCpMT
=Cp FO ( T O −T referencia ) −CpF ( T −T referencia ) + λm v
dt
dT 1
[ CpM
dt
=Cp FO T O−CpFT + λm v x
Cp ]
dT λ́ m v
donde nos queda , Ḿ =F O T´O −F O T +
dt Cp

Linealizando.

F(F o
= F o´T o + F́ o ( T −T´ o ) + T´ o (F o− F́ o)
,T )

Sustituyendo

dT ´ λ́ m v
Ḿ =T o F o− F o´T o − F́o ( T −T́ )−T́ (F o− F́ o)+
dt Cp
Se aplica lacondición de desviación : X d =X − X́

d (T d −T́ ) ´ λ́
Ḿ =T O ( F− F́ do )− F´o T − F́ o T d −T́ F do + (mdV + ḿV )
dt Cp
haciendo cero lo estacionario que no multiplica lo que está en desviación nos
queda:

d Td ´ d d d λ́ d
Ḿ =T O F o −F́ o T −T́ F o + m
dt Cp V

d Td d d λ́ d
Ḿ =( T´O−T́ ) Fo − F́ o T + m
dt Cp V

dTd λ́ d 1
[ Ḿ
dt
+ F́ o T d =( T´O−T́ ) Fdo +
Cp ]
mV x
F́ o

Ḿ d T d d T´O −T́ d λ́
+T = Fo+ mVd
F́ o dt F́ o Cp F́ o

Ḿ 2000 lb
Tp= = =40 min
F́ o 50 lb/min

T´O −T́ ( 50−120)℉ ℉

Kp1= = =−1.4
F́ o 50 lb /min lb/min

λ́ 1000 btu /lb

Kp 2= =
Cp F́ o (1 btu )¿ ¿
lb∗℉
Modelo Matemático.

d Td d d d
40 +T =−1.4 F O+ 20 mV
dt
Diagrama de bloque del modelo matemático.

F dO d Td d d
40 +T =−1.4 F O
dt
+
Td
+
m dV d Td d d
40 +T =20 mV
dt
Función Transferencia.

d Td d d d
40 +T =−1.4 F O+ 20 mV
dt
40 S T (s )+T ( s)=−1.4 Fo (s )+20 m v (s)

( 40 S+1 ) T (s )=−1.4 Fo( s) +20 m v(s)

T(s) −1.4 T( ) 20
= ; s =
Fo(s) 40 S+1 m v (s) 40 S +1

Diagrama de la función transferencia.

Fod(s) −1.4
40 S +1
+
T (s )
+
m v (s) 20
40 S +1

¿Cuál será el valor final de la variable de control si ¿ varía a 45 ℉ y mv varía 5

lb/min?
℉ ℉
[
T d= −1.4
lb /min
( 45−50 ) lb /min+20
lb/ min ]
( 5 ) lb /min ∗( 1−e−5 )=106.99 ℉

T =T d + T́ =( 106.99+120 ) ℉=226.99 ℉

7.- el mismo ejercicio 5, pero ahora ambas variables de entrada son

perturbadoras.

FO
F́ O=50 lb/ min

T´O=50 ℉

T́ =120 ℉
m v =3.5 lb /min

Ḿ =2 000 lb
λ=1000 btu/lb
 TO

mV
F

variables

VC →T
VP → F O ,T O

VM → mV

Balance de materia
Acum=entrada−salida
dm
=m entrada−m salida
dt
m=ρv ; ¿
dρV
=ρ F o− ρF
dt
dV
ρ =ρ ( Fo −F )
dt
dV
=F o−F
dt
Resultando :F =F o

Balance de energía.
dE ^ ^
= E E + E S +Q ; E=H donde , ^
E=^
H
dt
dH ^ ^
= H entrada + H salida +Q ; H
^ =CpF ( T −T referencia ) ; Q=λ m v ,donde , H =CpMT , sustutiyendo nos queda:
dt
dCpMT
=Cp FO ( T O −T referencia ) −CpF ( T −T referencia ) + λm v
dt
 dT 1
[ CpM
dt
=Cp FO T O−CpFT + λm v x
Cp ]
dT λ́ m v
donde nos queda , Ḿ =F O T 0−FO T +
dt Cp

Linealizando.

F(F o ,T o) = F o´T o+ F́ o ( T 0−T´ o ) + T´ o ( F o− F́ o)

F(F o ,T )= F ´o T + F́ o ( T −T́ ) + T́ (F o− F́ o)

Sustituyendo

dT λ́ m v
Ḿ = F o´T o + F́ o ( T 0−T´ o ) + T´ o ( F o− F́ o ) −F ´o T − F́ o ( T −T́ )−T́ ( F o− F́ o ) +
dt Cp

Se aplica la condición de desviación : X d =X − X́

d ( T d−T́ ) λ́
Ḿ = F o´T o + F́ o T do + T´ o F do− F ´o T −F́ o T d−T́ F do + (mdV + ḿv )
dt Cp
Haciendo’ cero lo estacionario que no multiplica lo que está en desviación nos
queda:

d Td d d d d λ́ d
Ḿ = F́ o T o + T´ o F o − F́ o T −T́ F o + m
dt Cp V

dTd λ́ d 1
( Ḿ
dt
+ F́ o T d= F́ o T do + ( T´ o−T́ ) Fdo +
Cp
mV x
F́o )
Ḿ d T d d ( T´ −T́ ) d λ́ d
+T =T do + o Fo+ mV
F́ o dt F́ o F́o Cp

Ḿ 2000 lb
Tp= = =40 min
F́ o 50 lb/min
℉
Kp 1=T do =1
℉
 ( T´ o−T́ ) (50−120) ℉ ℉
Kp2= = =−1.4
F́ o 50 lb/min lb/min

λ́ 1000 btu/lb
Kp3= =
F́ o Cp (1 btu )¿ ¿
lb∗℉

Modelo matemático.

d Td d d d d
40 +T =T o −1.4 F o +20 m V
dt

Diagrama de bloque del modelo matemático.

T dO d Td d d
40 +T =T o
dt

F do d Td d d
+
Td
40 +T =−1.4 F o +
dt +

m dv d Td d d
40 +T =20 mV
dt

Función Transferencia.
40 S T (s )+T (s)=T o (s)−1.4 Fo(s) +20 mV (S)

T (s ) ( 40 S+1 ) =T o(s) −1.4 F0 (S) +20 mV (S)

T ( s) 1 T( ) −1.4 T( ) 20
= ; s = ; s =
T o (s) 40 S+1 F0 (S) 40 S+1 m V (S) 40 S+1

Diagrama de bloques de la función transferencia.

T o(s ) 1
40 S +1

F 0(S ) −1.4 +
Td
40 S +1 +
 +

m V (S) 20
40 S +1

8.- RESOLVER EL MODELO MATEMATICO DE SEGUNDO ORDEN.

CONSIDERANDO DOS ECUACIONES DE PRIMER ORDEN
Modelo (1): Modelo (2):

d Td d d d d d m dv d d
40 +T =−1.4 F O+ T O +20 m V 1.2 +m v =0.8 P v
dt dt

d Td d d d d
Del modelo (1): 40 +T =−1.4 F O+ T O +20 m V , nos queda :
dt

d Td d d d d
40 +T +1.4 F O−T O=20 m V
dt

dTd 1
multiplicando por la inversa : 40 (
dt
+T d +1.4 FOd −T dO=20 m dV x
20 )
40 d T d 1 d 1.4 d 1 d d
nos queda: + T + F O− T O=mV
20 dt 20 20 20

sustituyendo mdV del modelo ( 1 ) en el modelo ( 2 ) nos queda :

40 d T d 1 d 1.4 d 1 d
1.2d ( 20 dt 20
+ T + F − T )
20 O 20 O 40 d T d 1 d 1.4 d 1 d
+ + T + F O− T O=0.8 P dv
dt 20 dt 20 20 20
ordenando la ecuación , nos queda :
d d
(1.2)(40) d 2 T d (1.2)(1) d T d (1.2)(1.4) d F o ( 1.2 ) ( 1 ) d T o 40 d T d 1 d 1.4 d 1 d d
2
+ + − + + T + F o − T o=0.8 P v
20 dt 20 dt 20 dt 20 dt 20 dt 20 20 20
multiplicando por 20 nos quedala siguiente ecuación :
d d
( 1.2 )( 40 ) d 2 T d (1.2 ) ( 1 ) d T d ( 1.2 )( 1.4 ) d F o ( 1.2 )( 1 ) d T o 40 d T d 1 d 1.4 d 1 d
( 20 dt 2
+
20 dt
+
20 dt
−
20
+
dt 20 dt 20
+ T +
20 20 )
F o− T o =0.8 Pdv x 20
 d d
d2 T d d Td d Fo d To d Td d d d d
( 1.2 ) ( 40 ) 2
+ (1.2 )( 1 ) + (1.2 ) ( 1.4 ) −( 1.2 ) ( 1 ) + 40 +T + 1.4 Fo −T o =(0.8 Pv )
dt dt dt dt dt
d d
d2T d d Td d Fo d To d Td d d d d
48 2
+1.2 + 1.68 −1.2 + 40 +T +1.4 F o −T o=0.8 P v
dt dt dt dt dt

Agrupandotérminos nos queda finalmente el siguiente modelo matemático :

d d
d2T d dTd d d Fo dTo d d d
48 2
+ 41.2 +T =−1.68 +1.2 −1.4 F o +T o + 0.8 P v
dt dt dt dt

T p2=48 , donde :Tp= √ 48=6.92 min

41.2 41.2
2 εTp=41.2 , donde despejando ε nos queda :ε = = =2.98( proceso sobreamortiguado)
Tp∗2 6.92∗2

Diagrama de bloque del modelo matemático.

F do d2T d dTd d d Fo d
d
48 2
+ 41.2 +T =−1.68 −1.4 F o
dt dt dt

+
T dO d2T d dTd d dTo d
d
+ Td
48 2
+ 41.2 +T =1.2 +T o
dt dt dt +

Pdv d2T d dTd d d

48 2
+ 41.2 +T =0.8 P v
dt dt
Función transferencia.
d d
d2T d dTd d d Fo dTo d d d
48 2
+ 41.2 +T =−1.68 +1.2 −1.4 F o +T o + 0.8 P v
dt dt dt dt
2
48 S T ( s )+ 41.2 S T (s )+T ( s)=−1.68 S F o (s )−1.4 F o (s )+ 1.2 S T o (s )+T o ( s) +0.8 Pv( s )

Sacando factor comúnT ( s) nos queda :

T (s ) ( 48 S 2+ 41.2 S+ 1 )=−1.68 S Fo ( s )−1.4 Fo ( s )+1.2 S T o ( s) +T o (s )+ 0.8 Pv ( s)

T (s ) T ( s) 1.2 S+1 T ( s) 0.8

−1.68 S−1.4 = =
= T o ( s) 48 S2 +41.2 S+1 2
Pv ( s) 48 S +41.2 S +1
F o ( s) 48 S2 +41.2 S +1
Diagrama de bloque de la función transferencia.
F o (s ) −1.68 S−1.4
48 S 2 +41.2 S +1

T o (s ) 1.2 S+1 + T (s )
48 S 2 +41.2 S +1 +
+

Pv (s ) 0.8
2
48 S +41.2 S +1