UNIVERSIDAD TECNICA PARTICULAR DE LOJA

SIMULACIÓN Y CONTROL
DE PROCESOS

CONTROL DE PROCESOS PARA LA PRODUCCIÓN DE

TRINITROTOLUENO (TNT) A PARTIR DE UNA MEZCLA SULFONÍTRICA
INTEGRANTES:

JESSICA CASTILLO

ANDREA JARAMILLO

KARINA VALVERDE

PAULA VANEGAS

DOCENTE:

DR. MIGUEL A. MENESES

 CONTROL AUTOMÁTICO
1. Tema Específico
Control de proceso De Mononitrotolueno (MNT) a partir de la primera Etapa
de Nitración
2. Resumen
El presente proyecto consiste en realizar el control automático del proceso de producción
de MNT que se da en la primera nitración para la producción del TNT, para ellos se utiliza
dos equipos, un mezclador y un reactor utilizando el Software Matlab y su herramienta
Simulink.
3. Metodología

3.1 Descripción proceso

La producción de TNT en el cual se empleará un Mezclador y un Reactor. En el
mezclador existe una mezcla sulfonítrica compuesta por distintas relaciones peso/peso;
25% de HNO3, 55% de H2SO4 y 20% de H2O; que no supere una temperatura de 40ºC
que será controlada por una camisa refrigerante. Se ingresará al primer Reactor una
tonelada de tolueno al que se le adicionará una mezcla sulfonítrica, preparada en el
mezclador, donde se producirá Mononitrotolueno (MNT). La reacción que se dará es de
tipo exotérmica a una temperatura de 55ºC, es por ello que con una chaqueta de
refrigeración circulara agua para evita temperatura superior a la requerida.

3.2 Descripción equipos

Mezclador: Es un mezclador que consiste en la disolución sulfonítrica para los
compuestos de entrada. La disolución al mezclarse desprende un excesivo calor de
manera que se debe controlar la temperatura en este equipo
Reactor: En este equipo de dará la reacción de la mezcla sulfonítrica con tolueno a una
temperatura de 55ºC, además contará con un sensor de temperatura, donde al elevarse la
temperatura del reactor transmitirá una señal a la válvula de paso para que circule el agua
refrigerante.
Se requiere el mismo el mismo control del reactor en el mezclador.
3.3 Método
El método que se ha considerado para poder realizar este proyecto es del Ziegler Nichols
que permite ajustar un regulador PID de forma empírica, sin necesidad de conocer las
ecuaciones de planta o sistema controlado. Además, permite definir las ganancias
proporcional, integral y derivativa a partir de la respuesta de lazo cerrado en este caso. Es
de lazo cerrado o control de circuito cerrado ya que compara el punto de control
(referencia) con la variable controlada y determina la acción correctiva.
 3.4 Selección controladores

TIPO DE ACCIÓN
TIPO DE TIPO DE DEL TIPO DE
SENSOR CONTROLADOR CONTROLADOR VÁLVULA

Mezclador 1 Controlador de Aire para cerrar.

Sensor de Temperatura, Acción Inversa Válvula de vastago
temperatura Controlador Recíproco
(Termocupla) proporcional

Reactor 1 Controlador de
Sensor de Temperatura, Acción Directa Válvula de vastago
temperatura Controlador Recíproco
(Termocupla) proporcional

Mezclador

Nomenclatura
f3= 0.159
ft 3 H O
2
min
f1= 0.608
ft 3 H SO
2 4
min
f4=83.57
ft 3 Mezcla
min
f2= 0.24
ft 3 HN0
3
min T Ci= 68 ºF

T C=86ºF
 Cp1=¿0.56BTU/lbºF ρ1= 114.49 lb/ ft 3
Cp2=¿ 0.75 BTU/lbºF ρ2= 94.77 lb/ ft 3
Cp 3=¿ 1.81 BTU/lbºF ρ3= 62.05 lb/ ft 3
Cp4 =¿ 1.04 BTU/lbºF ρ4 = 83.57 lb/ ft 3
Cpc =¿ 1.81 BTU/lbºF ρc = 62.05 lb/ ft 3
f c= 0.55 ft 3/min Vp= 1.43 ft 3 /min

Balance de materia

d f 4 (t )
ρ1 f 1 ( t ) + ρ 2 f 2 ( t )+ ρ3 f 3 ( t ) −ρ4 f 4 ( t )=V
dt
Funciones de transferencia
K1 K K
F 4 ( s )= F ( s )+ 2 F ( s ) + 3 F ( s )
τs+1 1 τs +1 2 τs+1 3
Balance de Energía

d T 2(t )
f 1 ( t ) ρ1 C p1 T 1 ( t )+ f 2 ρ2 C p 2 T 1 ( t ) +f 3 ρ3 C p 3 T 1 ( t ) −f 4 ρ4 C p 4 T 2 ( t )−UA [ T 2 ( t )−T c ( t ) ]=V ρ4 C p 4
dt
Funciones de transferencia
K4 K
Γ 2 ( s )= Γ ( s )+ 5 F1 ( s )
τs+ 1 1 τs +1
Balance de Chaqueta

f c ( t ) ρc C pc T d T c (t )
ci +¿UA [ T 2 ( t ) −T c (t ) ] −f c ( t ) ρ c C pc T c ( t )= Vρ C p ¿
dt

Funciones de transferencia
K6 K
Γ c ( s )= F ( s) + 7 Γ 2 ( s)
τs+1 s τs+1

Parámetros del proceso

 Tabla1. Resultado de Constantes

τ s 1=0.015 min τ s 2=0.014 min τ s 3=1.296 min

K1= 1.37 K4= -0.0509 K7= 0.00276
K2= 1.13 K5= -48.78 τ v=0.20 min
K3= 0.74 K6= 32.63 = 50

Diagrama de bloques del mezclador

Figura 1. Diagrama de bloques del mezclador

Simulación de proceso
Figura 2. Simulación de Mezclador

Reactor
Nomenclatura
Ti 4 = 131 ºF ρ5= 53.98 lb/ ft 3

Ti 5= 77 ºF ρ6 = 74.63 lb/ ft 3

T c= 68 ºF ρ7 = 62.43 lb/ ft 3

Cp4 =¿1.04 BTU/lbºF ρc = 62.43 lb/ ft 3

Cp 5=¿ 0.73 BTU/lbºF C Ai4 =¿ 1.1317 lb mol/ ft 3

Cp6=¿ 0.64 BTU/lbºF C Ai5=¿ 0.2159 lb mol/ ft 3

Cp7=¿ 1.81 BTU/lbºF C Ai6 =¿ 1.347 lb mol/ ft 3

Cpc =¿ 1.81 BTU/lbºF V= 110.77 ft 3 /h

f 4= 1.2105 ft 3/min ft 3
k 0=0.044
lbmol
f 5= 0.68 ft 3/min
k final= 0.0083 lb/ ft 3
f 6= 1.84 ft 3/min
E=¿ 27.820 BTU/lb mol
f 7= 0.58 ft 3/min
R=¿1.987 BTU/lb mol
ρ4 = 83.57 lb/ ft 3
∆ H =−15384BTU/lb mol
Balance de materia de reactor

Vd C A 6 (t)
f 4 ( t ) C Ai 4 ( t ) +f 5 ( t )−f 6 ( t ) C Ai 6 ( t )−V rA ( t )=
dt
−E
RT ( t ) Vd C A 6 (t)
f 4 ( t ) C Ai 4 ( t ) +f 5 ( t )−f 6 ( t ) C Ai 6 ( t )−V ko e C 2A (t )=
dt
Funciones de transferencia

K 8 C Ai 4 ( s ) + K 9 f 4 ( s )−K 10 T ( s )+ K 11 F 5
C Ai6 ( s )=
τs+1
Balance de energía del Reactor

d T 2(t )
f 4 ( t ) ρ4 C p 4 T 0 ( t )−UA [ T ( t )−T 2 ( t ) ]−f 6 ( t ) ρ6 C p 6 T 2 ( t )−V rA (∆ H r )=Vρ C v
dt
Funciones de transferencia

K 12 T 0 ( s ) K 13 f 6 ( s ) K 14 Γ 2 ( s ) K 15 C A ( s )
Γ ( s )= + + −
τs+1 τs+1 τs+1 τs+ 1
Balance de la chaqueta reactor

f 7 ( t ) ρ7 C p 7 T d T c (t )
ci +¿UA [ T 2 ( t ) −T c (t ) ] −f 7 ( t ) ρ 7 C p7 T c ( t ) =Vρ C p ¿
dt

Funciones de transferencia

K 16 Г ci (s) K 17 F c (t) K 18 Γ ( t)
Г c ( s )= + +
τ 3 s+ 1 τ 3 s +1 τ 3 s +1

Diagrama de Bloques del reactor

Figura 3. Diagrama de bloques del reactor

Simulación de proceso del reactor
Figura 4. Simulación de Mezclador

4. Resultados

Figura 5. Respuesta de controlador del Mezclador

 Figura 6. Respuesta de controlador del Reactor

5. Conclusiones
Mediante simulink se realizó el control de los equipos para la producción de dinitrotolueno
(DNT), tanto del mezclador (figura 2) y del reactor (figura 4), el control de los equipos tiene
el mismo funcionamiento básico, siendo controlar la temperatura de los equipos con
modelo integral proporcional derivativo.
Para el mezclador, con una ganancia de 20 se observaron fluctuaciones a partir de 80; de
manera que según la respuesta del controlador de la figura 5 se puede decir que existe un
error para este control, ya que no existe una respuesta coherente, por lo que se debería
realizar nuevamente el control del proceso o evaluar otro modelo para este tipo de
sistema.
Las gráficas representadas en el reactor (figura 6) nos muestran que no existe una
estabilidad en el sistema, por lo que no fue posible ejecutarlos en controladores PI y PID,
puesto a que necesitamos el periodo de oscilación y en este caso no lo tenemos ya que
no cambia nisiquiera aumentando las ganancias, por lo que no es posible obtener debido
a que la señal del controlador varia con el tiempo.
6. Anexos

 MEZCLADOR
Balance de materia

d f 4 (t )
ρ1 f 1 ( t ) + ρ 2 f 2 ( t )+ ρ3 f 3 ( t ) −ρ4 f 4 ( t )=V
dt
Balance en estado estacionario

d f 4(t )
ρ1 f´1 + ρ2 f´2 + ρ3 f´3−ρ4 f´4 =V
dt
ρ1 ρ2 ρ3 ρ4 Vd f 4 ( t )
F1(t ) + F 2 ( t ) + F 3 (t ) - F 4 ( t ) =
ρ4 ρ4 ρ4 ρ4 ρ4 dt

Vd f 4 ( t ) ρ1 ρ ρ
+ F 4 ( t )=¿ F 1 ( t ) + 2 F 2 ( t ) + 3 F 3 (t )
ρ4 dt ρ4 ρ4 ρ4

ρ4 ρ2 ρ3 V
K 1= ; K 2= ; K 3= ; τ=
ρ4 ρ4 ρ4 ρ4

F4 (t )
τd + F 4 ( t ) =K 1 F 1 ( t )+ K 2 F 2 ( t ) + K 3 F 3 ( t )
dt
τs F 4 (s) + F 4 ( s )=K 1 F 1 ( s )+ K 2 F 2 ( s ) + K 3 F3 ( s )
K1 K K
F 4 ( s )= F ( s )+ 2 F ( s ) + 3 F ( s )
τs+1 1 τs +1 2 τs+1 3
Balance de energía

d T 2(t )
f 1 ( t ) ρ1 C p1 T 1 ( t )+ f 2 ρ2 C p 2 T 1 ( t ) +f 3 ρ3 C p 3 T 1 ( t ) −f 4 ρ4 C p 4 T 2 ( t )−UA [ T 2 ( t )−T c ( t ) ]=V ρ4 C p 4
dt
Linealización

f 1 ( t ) ρ1 C p1 T 1 ( t ) ≅ f´1 ρ1 C p 1 T´ 1+ f´1 ρ1 C p 1 Γ 1 ( t )+ ρ1 C p 1 T́ 1 F1 (t )

− f´1 ρ1 C p 1 Γ 1 (t )− ρ1 C p 1 T́ 1 F 1 ( t )+ f 1 ( t ) ρ1 C p 1 T 1 ( t ) + f 2 ρ2 C p 2 T 1 ( t ) + f 3 ρ3 C p 3 T 1 ( t )−f 4 ρ4 C p 4 T 2 ( t ) −UA [ T 2 ( t ) −T c

Balance de energía en estado estacionario 1

d Γ2 ( t )
V ρ4 C p 4 +f 4 ρ4 C p 4 Γ 2 ( t )+UA Γ 2 ( t )=− f´1 ρ1 C p 1 Γ 1 ( t )− ρ1 C p 1 T´ 1 F 1 ( t )+ f 2 ρ2 C p 2 Γ 1 ( t )+ f 3 ρ3 C p 3 Γ 1 ( t )
dt
d Γ2 ( t )
V ρ4 C p 4 +f 4 ρ4 C p 4 Γ 2 ( t )+UA Γ 2 ( t )=( f 2 ρ2 C p 2+ f 3 ρ 3 C p 3− f´1 ρ1 C p 1) Γ 1 ( t )-
dt
ρ1 C p 1 T́ 1 F 1 ( t )
K4 K
Γ 2 ( s )= Γ ( s )+ 5 F1 ( s )
τs+ 1 1 τs +1
Dónde
V ρ4C p4
τ=
f 4 ρ 4 C p 4 +UA

f 2 ρ2 C p 2+ f 3 ρ3 C p 3− f´1 ρ1 C p 1
K4=
f 4 ρ 4 C p 4 +UA

−ρ1 C p 1 T́ 1
K 5=
f 4 ρ4 C p 4 +UA
Balance de energía de la Chaqueta

f c ( t ) ρc C pc T d T c (t )
ci +¿UA [ T 2 ( t ) −T c (t ) ] −f c ( t ) ρ c C pc T c ( t )= Vρ C p ¿
dt

Balance Estacionario de la chaqueta

f c ( t ) ρc C pc T c (t)≅ f´c ρ c C pc T́ c + f´c ρc C pc Γ c ( t )+ ρc C pc T́ c Fc ( t )

d Γ c (t )
ρc C pc T ci Fc ( t ) +UA [ Γ 2 ( t ) −Γ C (t) ]−(− f´c ρc C pc Γ c ( t ) ) −ρc C pc T́ c F c ( t )=Vρ C p
dt

Vρ C p d Γ c ( t )+ UA [ Γ C ( t )− f´c ρc C pc Γ c (t ) ]= ρc C pc T ci F c ( t ) +¿ UA Γ 2 ( t ) + ρc C pc T´ c F c ( t )

Vρ C p d Γ c ( t )+ ¿UA− f´c ρ c C pc ¿ Γ c ( t )= ρc C pc F c ( t ) ( T ci + T́ c ) +UA Γ 2 ( t )

Vρ C p d Γ c ( t ) ρc C pc UA Γ 2 ( t )
+ Γ C ( t )=¿ F c ( t ) ( T ci + T´ c ) +
UA− f´c ρ c C pc UA− f´c ρ c C pc UA− f´c ρ c C pc
 Vρ C p ρ c C pc (T ci−T´ c ) UA
τ= K6= K 7=
UA + f´c ρc C pc UA + f´c ρc C pc UA + f´c ρc C pc

τd Γ c ( t ) + Γ c ( t ) =K 6 F c ( t ) + K 7 Γ 2 ( t )

K6 K
Γ c ( s )= F ( s) + 7 Γ 2 ( s)
τs+1 s τs+1

CONTROL DEL MEZCLADOR

V
τ 1= =0.015 min
ρ4

V ρ4 C p 4
τ2= =0.014 min
f 4 ρ 4 C p 4 +UA
V ρCp
τ 3= =1.296 min
UA −ρ5 C p 5 F 5

p 1 114.49
K1 = = = 1.37
p 4 83.57
p 2 94.77
K2 = = = 1.13
p 4 83.57
p 2 62.05
K3 = = = 0.74
p 4 83.57
f 2 ρ2 C p 2 + f 3 ρ3 C p 3 −f´1 ρ1 C p 1
K 4= =−0.0509
f 4 ρ4 C p 4+ UA

−ρ1 C p 1 T´ 1
K5 = = -48.78
f 4 ρ 4 C p 4 +UA

−ρ5 C p 5 (Tc−T´ c ) −62.5∗1.81(68+104)

k 6= = = 32.63
UA −ρ5 C p 5 F 5 61.94
UA
K 7= =0.00276
UA + f´c ρc C pc

τ s 1=0.015 min τ s 2=0.014 min τ s 3=1.296 min

K1= 1.37 K4= -0.0509 K7= 0.00276
K2= 1.13 K5= -48.78 τ v=0.20 min
K3= 0.74 K6= 32.63 = 50
 W ( s) Kv
Gv(s) = =
M ( s) τ vs + 1
ẃ ( ln❑ ) 0.55 ( ln 50 )
kv = = =0.0215
100 100
Sensor Transmisor

C(s) Kt
H(s) = =
Γ (s) τ Γ +1
s

100−0 %T0
KT = =1
50−150 ºF
Función de transferencia para el controlador PID

M (s) 1
Gc(s) = =kc (1+ +τ )
R ( s )−C( s) τ 1+s Ds
0.0215
G1(s)= Gv(s) Cs(s) H(s)=
0.2 s+1
1.0
G2(s)= Gc(s) H(s)=
0.75 s+ 1
kw ks kw∗ks 3.919∗(17.94)
Gs(s)= * = = =
τ 2 s+1 τ 3 s+ 1 ( τ 2 s+1 ) ( τ 3 s +1 )−k 5 ( 0.014 s +1 ) (1.296 s+1 ) +0.0509
70.306
2
0.0181 s +1.31 s +1.0509
kf (τs +1) −29.658(1.296 s +1) −38.426 s−29.65
Gf(s)= = =
( τs+1 ) ( τs+1 )−k 5 ( 0.015 s+ 1 )( 1.296 s +1 ) +48.78 0.0194 s s +1.311 s+ 49.78

REACTOR
Balance de materia
 Vd C A 6 (t)
f 4 ( t ) C Ai 4 ( t ) +f 5 ( t )−f 6 ( t ) C Ai 6 ( t )−V rA ( t )=
dt
−E
Vd C A 6 (t)
f 4 ( t ) C Ai 4 ( t ) +f 5 ( t )−f 6 ( t ) C Ai 6 ( t )−V ko e RT (t ) C 2A (t )=
dt
Linealización

f 4 ( t ) C Ai 4 ( t ) ¿ ≅ f´4 C´Ai 4 + f 4 ( t ) C´Ai 4 + f´4 C Ai 4 ( t )

f 6 ( t ) C Ai 6 ( t ) ≅ f´6 C´Ai6 + f 6 ( t ) C´Ai6 + f´6 C Ai 6 ( t )

−E
−E −E E −E
V ko e RT ( t ) 2
C ( t ) =V ko e
A
RT ( t )
C´2A + V ko e
RT
( )RT 2
Ć 2
A T(t)+ V ko e RT ( t ) 2C´ A C A ( t )

−E −E −E
E ´2
f´4 C´Ai 4 + f 4 ( t ) C´Ai 4 + f´4 C Ai 4 ( t ) + f 5 ( t )− f´6 C Ai 6−f 6 ( t ) C´Ai6 −f´6 C Ai 6−V ko e RT (t ) C´2A −V ko e RT
( )
RT 2
C A T ( t )−V ko e RT (t )

−E
f´5+ ¿ f´4 C´Ai 4 - f´ C´ −V e RT C´2 = 0
6 Ai 6 ko A

−E
−E
E ´2 Vd C A 6 (t )
f 4 ( t ) C´Ai 4 + f´4 C Ai 4 (t ) + f 5 ( t )−f 6 ( t ) C´Ai 6− f´6 C Ai 6−V ko e RT
( )
RT 2
C A T ( t ) −V ko e RT ( t) 2 C´ A C A (t )=
dt
−E
Vd C A 6 (t) ´
+ f 6 C Ai6 +V ko e RT (t ) 2 C´ A C A6 ( t )=f 4 ( t ) C´Ai 4 +
dt
−E
E ´2
f´4 C Ai 4 ( t )−f 6 ( t ) C´Ai 6+ f 5 ( t )−V ko e RT
( )RT 2
C A T (t )

Vd C A 6 ( t ) −E
+C Ai6 ( f´ +V e RT ( t ) 2 C´ )= f 4 ( t ) C´Ai 4 +
dt 6 ko A
−E
f´4 C Ai 4 ( t )−f 6 ( t ) C´Ai 6+ f 5 ( t )−V ko e RT ( RTE ) C T ( t )
2
A

V f 4 (t )
τ= −E K 8= −E
RT ( t )
f´6 +V ko e 2 C´ A f´6 +V ko e RT ( t ) 2 C´ A

K9=
C´Ai 4−C´Ai6
−E
K 10=
V
( RTE ) r 2 A

RT ( t )
f´6 +V ko e
−E
2 C´ A
f´6+ V ko e RT ( t ) 2 C´ A
f5
K 11 = −E

f´6 +V ko e RT (t ) 2 C´ A
 d C A 6 ( t)
τ + C Ai6 ( t )=K 8 C Ai 4 ( t )+¿ K 9 f 4 ( t ) −¿ K 10 T ( t ) +¿ K 11 F 5
dt
τs C A 6 ( s ) +C Ai 6 ( s )=¿ K 8 C Ai 4 ( s ) +¿ K 9 f 4 ( s )−¿ K 10 T ( s ) +¿ K 11 F 5

K 8 C Ai 4 ( s ) + K 9 f 4 ( s )−K 10 T ( s )+ K 11 F 5
C Ai6 ( s )=
τs+1
Balance de energía del Reactor

d T 2(t )
f 4 ( t ) ρ4 C p 4 T 0 ( t )−UA [ T ( t )−T 2 ( t ) ]−f 6 ( t ) ρ6 C p 6 T 2 ( t )−V rA (∆ H r )=Vρ C v
dt
Linealización

f 4 ( t ) ρ4 C p 4 T 0 ( t ) ≅ f´4 ρ 4 C p 4 T́ 0+ f´4 ρ 4 C p 4 T 0 ( t )+ f 4 ( t ) ρ4 C p 4 T´ 0

f 6 ( t ) ρ6 C p 6 T 2 ( t ) ≅ f´6 ρ6 C p 6 T́ 2 + f´6 ρ6 C p 6 T 2 (t ) + f 6 ( t ) ρ6 C p 6 T́ 2
−E
−E E −E
V rA (t) ( ∆ H r )=V ko e RT C´2 ( ∆ H r )=V ko e RT (t )
A ( )RT
C 2A T ( t ) +( ∆ H r )+V ko e RT 2 C´ C ( t )
A A

f´4 ρ4 C p 4 T 0 + f´4 ρ4 C p 4 T 0 ( t ) +¿ f 4 ( t ) ρ4 C p 4 T́ 0 −UA ¿ T 2 ( t )−¿ f´6 ρ6 C p 6 T´ 2− f´6 ρ6 C p 6 T́ 2 ( t )−¿

−E −E −E d T 2(t )
f 6 ( t ) ρ6 C p 6 T́ 2 -V ko e RT C´2 ( ∆ H r )-V ko e RT C 2A T 2 ( t ) ( ∆ H r ) - V ko e RT 2 C´ C ( t )= VρC v
A A A dt
−E
f´4 ρ4 C p 4 T´ 0−UA [T́ −T´ c] - f´6 ρ6 C p 6 T´ 2 -V ko e RT C´2 ( ∆ H r )=0
A

f´4 ρ4 C p 4 T 0 ( t ) + f 4 ( t ) ρ 4 C p 4 T´ 0−UA [ Γ ( t ) −Γ 2 ( t ) ]− f´6 ρ6 C p 6 T 2 ( t ) −¿ f 6 ( t ) ρ6 C p 6 T́ 2−V ko

−E
e RT C 2A T 2 ( t ) ( ∆ H r )

−E d T 2(t )
−V ko e RT
2 C´ A C A ( t )=¿ VρC v dt
d T 2(t ) −E
VρC v + T 2 ( t ) ¿ e RT C 2 ( ∆ H ) +UA ¿=¿ f´4 ρ4 C p 4 T 0 ( t ) + f 4 ( t ) ρ 4 C p 4 T´ 0+ Γ 2 ( t ) UA + ¿
dt A r
−E
f 6 ( t ) ρ6 C p 6 T́ 2 ( t ) +V ko e RT 2 C´ C ( t )
A A

Vρ C v ρ4 C p 4 [ T́ 0 −T́ 2 ]
τ= K 12=
E E
f´6 ρ6 C p 6 +V ( ∆ H r ) r´A
( )
R T́ 2
+UA f´6 ρ6 C p 6 +V ( ∆ H r ) r´A
( )
R T́ 2
+UA
 −E
f´4 ρ4 C p 4 2V ( ∆ H r ) ko e RT C´ A
K 13= K 14=
E
f´6 ρ6 C p 6+ V ( ∆ H r ) r´A
( )R T́ 2
+ UA f´6 ρ6 C p 6 +V ( ∆ H r ) r´A
E
( )
R T´ 2
+UA

UA
K 15=
E
f´6 ρ6 C p 6+ V ( ∆ H r ) r´A
( )R T́ 2
+ UA

dT ( t )
τ2 + T 2 ( t )=K 12 T 0 ( t ) +¿ K 13 f 6 ( t ) + K 14 Γ 2 ( t ) −K 15 C A ( t )
dt

τ 2 ( s ) +T 2 ( s )=¿ K 12 T 0 ( s )+ ¿ K 13 f 6 ( s ) + K 14 Γ 2 ( s )−K 15 C A ( s )

K 12 T 0 ( s ) K 13 f 6 ( s ) K 14 Γ 2 ( s ) K 15 C A ( s )
Γ ( s )= + + −
τs+1 τs+1 τs+1 τs+ 1
Balance de la chaqueta reactor

f 7 ( t ) ρ7 C p 7 T d T c (t )
ci +¿UA [ T 2 ( t ) −T c (t ) ] −f 7 ( t ) ρ 7 C p7 T c ( t ) =Vρ C p ¿
dt

f 7 ( t ) ρ7 C p 7 T ci (t )≅ f´7 ρ 7 C p 7 T´ci+ f´7 ρ7 C p 7 T ci ( t ) +f 7 (t) ρ7 C p 7 T´ci

f 7 ( t ) ρ7 C p 7 T c (t) ≅ f´7 ρ7 C p 7 T´ c + f´7 ρ7 C p 7 T c ( t )+ f 7 (t ) ρ7 C p7 T´ c

f´7 ρ 7 C p 7 T´ci + f´7 ρ7 C p 7 T ci ( t ) + f 7 ( t ) ρ7 C p 7 T´ci + UA [ T 2 ( t )−T c (t) ] −f´7 ρ7 C p 7 T́ c + f´7 ρ7 C p 7 T c ( t ) + f 7 (t) ρ7 C p 7 T´ c

d T c(t)
= V c ρ c Cpc
dt
f´c ρc C pc T´ci +UA [ T́ −T́ c ]−f´c ρc C pc T´ c = 0

f´7 ρ 7 C p 7 T ci + f 7 ( t ) ρ7 C p 7 T´ci +UA [ T 2 ( t )−T c ( t ) ]− f´7 ρ7 C p 7 T c ( t )−f 7 ( t ) ρ7 C p 7 T́ c =¿

d T c(t)
V c ρ c Cp c
dt
d T c (t )
Vρ C p + f́ 7 ρ 7 Cp7 Г c ( t ) +UA Г c ( t )=UA Г 2 ( t )+ f´7 ρ7 C p 7 T ci ( t )+ f 7 (t) ρ7 C p 7 T´ci
dt
d T c (t )
Vρ C pc + Г c ( t ) ( f´7 ρ7 C p 7 +UA )=UA Г 2 ( t ) + f´7 ρ7 C p 7 T ci ( t )+ f 7 (t )ρ 7 C p 7 T´ci
dt
V c ρc Cpc ρc Cpc ( T́ ci−T́ c ) f́ 7 ρ7 Cp7 T́ ci UA
τ 3= ; K 16= ; K 17 = ; K 18=
f́ 7 ρ7 Cp7 +UA f́ 7 ρ7 Cp7 +UA f́ 7 ρ7 Cp7 +UA f́ 7 ρ7 Cp7 +UA
 d T c (t)
τ3 + Г c ( t )=K 16 Г ci ( t )+ K 17 F c ( t )+ K 17 Γ ( t )
dt
K 16 Г ci (s) K 17 F c (t) K 18 Γ ( t)
Г c ( s )= + +
τ 3 s+ 1 τ 3 s +1 τ 3 s +1

Control del Reactor

(62.43)(1.81)(66.2−86)
K 16= ρc Cpc ( T́ ci −T́ c ) = =−23.80
94.17
(0.58)(62.43)(1.81)
K 17= f́ 7 ρ7 Cp7 T́ ci= =0.69
94.17
103 (0.278)
K 18=UA = =0.30
94.17

τ 4 s= 1.039 min τ 5 s= -0.052 min τ 6 s= 1.92 min

K8= 0.657 K12= 0.12 K16= -23.80
K9= 0.117 K13= -0.00566 K17= 0.69
K10= 0.067 K14= 8.42 K18= 0.30
K11= 0.369 K15= -0.016 τ v= 0.20 min

W ( s) kv
Gv(s)= =
M ( s) τ v s +1
ẃ ( ln❑ ) 0.92 ( ln50 )
kv = = =0.0359
100 100

Sensor Transmisor

C(s) Kt
H(s) = =
Γ (s) τ Γ +1
s

100−0 %T0
KT = =1
50−150 ºF
Reacción de transferencia para el controlador PID

M (s) 1
Gc(s) = =kc (1+ +τ )
R ( s )−C( s) τ 1 s Ds
0.0359
G7(s)= Gv(s) Cs(s) H(s)=
0.2 s+1
 1
G8(s)= Gk(s) H(s)=
0.75 s+ 1
kw ks kw∗ks (−23.80)∗(−0.00566)
G9(s)= * = = =
τ 5 s+ 1 τ 6 s +1 ( τ 5 s+ 1 )( τ 6 s+1 )−k 5 (−0.52 s+1 )( 1.92 s+ 1 )−(−0.00566)
0.13328
−0.9932 s 2+1.4 s+ 1.00566
kf (τ 4 s+ 1) 8.42(1.039+1) 17.17
Gf(s)= = = 2
( τ 4 s+1 ) ( τ 6 s+1 )−k 5 ( 1.039 s+1 ) ( 1.92 s+1 )−8.42 1.9949 s +2.959 s−7.42
V c ρc Cp c (1.6)(62.43)(1.81) 180.79
τ 3= = = =1.92 min
f́ 7 ρ 7 Cp 7 +UA ( 0.58 ) ( 62.43 ) ( 1.81 )+103 (0.278) 94.17