Aspen Hysys - Simulacion Dinamica
Aspen Hysys - Simulacion Dinamica
Aspen Hysys - Simulacion Dinamica
PROCESOS QUMICOS
Control, Modelamiento y Simulacin dinmica
Facultad de Ingeniera
Programa de Ingeniera Qumica
Universidad del Atlntico
AUTORES:
Melanio Coronado Hurtado.
Jennifer Juliao Paula
Alberto Fuentes Enamorado
2012
PREFACIO
En un proceso qumico siempre estn ocurriendo cambios, y si no se responde con
las acciones apropiadas, las variables de proceso ms importantes (es decir aquellas
relacionadas con la seguridad, la calidad del producto, y el volumen de produccin)
no alcanzaran las condiciones de diseo adecuadas. Adems, si se consideran los
desafos que posee la industria qumica en el siglo XXI como son: el ahorro
energtico, el cuidado del medio ambiente, la seguridad de los procesos y el
desarrollo de mejores ambientes de trabajo, el control automtico es un
complemento que sirve a los procesos, tanto para compensar esos cambios
inherentes como para hacer ms eficientes todos los sistemas de gestin.
Sin embargo, la teora de control es muy extensa y en los ltimos tiempos ha
avanzado mucho gracias a la creacin de herramientas informticas ms
sofisticadas, hecho que ha contribuido a que el control de procesos se pueda
comprender mucho ms fcilmente con la ayuda de lenguajes de programacin y
simuladores de procesos qumicos. Es precisamente gracias a este avance, que en la
teora y la prctica del control automtico se han aportado los medios necesarios
para obtener un desempeo ptimo de los sistemas dinmicos, mejorar la
productividad y aligerar la carga de las operaciones manuales repetitivas y rutinarias
antes realizadas por el personal operador.
Para llegar a un mejor entendimiento de los procesos y de la mejor forma como
controlarlos se hace necesario realizar un estudio riguroso de la teora de control y
para ello se requiere un conocimiento matemtico avanzado y complejo. Pero esta
dificultad se puede disminuir, sin detrimento de su entendimiento, a travs de
herramientas informticas como: Simulink, LabView, Aspen Dynamics, y
Aspen-HYSYS, entre otras, que permiten hacer simulaciones de los procesos
reales con la realizacin interna de una gran cantidad de algoritmos para la solucin
de los diferentes problemas. Sin embargo, este ltimo es uno de los ms utilizados,
porque posee una de las interfaces ms amigables con el usuario, sumado al hecho
de que puede realizar simulaciones tanto en estado estacionario como en estado
dinmico, en una misma plataforma.
El desarrollo de un curso de control, modelamiento y simulacin asistido con
Aspen-HYSYS como ste es, por lo tanto, una oportunidad de multiplicar el
conocimiento acadmico de estudiantes de pre grado de Ingeniera Qumica y buscar
unir las herramientas informticas actuales con los procesos qumicos.
Autores
i
Ingeniera Qumica | Uniatlntico | 2012
TABLA DE CONTENIDO
2.
3.
4.
5.
6.
7.
8.
9.
iii
Ingeniera Qumica | Uniatlntico | 2012
MDULO I:
CONTROL DE PROCESOS
Diseo y Estrategias
1
Ingeniera Qumica | Uniatlntico | 2012
Los instrumentos que se requieren para hacer este lazo de control son: un
sensor/transmisor que es el dispositivo que mide el valor de la variable de proceso
a controlar y la transforma en seales entendibles para el siguiente elemento del
lazo. Un controlador que es un equipo capaz de interpretar la desviacin o
diferencia que hay entre el SP y la medida hecha por el sensor (de ah el nombre del
lazo), y genera una seal correctora que busca alcanzar el valor deseado; y por
2
Ingeniera Qumica | Uniatlntico | 2012
Las variables de salida del proceso son aquellas que dan informacin sobre el
estado de ste y usualmente estn asociadas con corrientes de salida o con
mediciones hechas dentro del proceso. A menudo se refiere a estas como
variables controladas, variables de proceso o PV, estas son medidas y reguladas
hasta su valor deseado o SP. Las variables comnmente controladas son: flujo,
por ejemplo alimentacin a un equipo, nivel de lquido, presin y temperatura de
un equipo (separador de fases, tanques de mezclado, tambores de reflujo,
intercambiadores de calor, etc.).
Las variables de entrada del proceso son variables independientes que afectan
las variables de salida. Estas, a su vez, pueden ser clasificadas en:
Variables manipuladas: (tambin llamadas variables de control) son ajustadas
libremente por un operador o un mecanismo de control a travs de un
elemento de control final como una vlvula de control. Un ejemplo de estas
puede ser el flujo msico del fluido caliente en un intercambiador de calor, el
flujo de una corriente de agua de proceso, etc.
Perturbaciones (o variables definidas externamente) que estn sujetas al
medio ambiente externo y por lo tanto no pueden ser controladas.
Generalmente estas se relacionan con corrientes de entrada y de salida del
proceso. Un ejemplo de estas es la temperatura de vapor de alta como fluido
caliente en un intercambiador o la composicin de la corriente de alimento de
un equipo.
Tericamente, en un proceso deberan controlarse todas las variables de salida, pero
en la realidad esto no es posible por diversas razones como:
3
Ingeniera Qumica | Uniatlntico | 2012
4
Ingeniera Qumica | Uniatlntico | 2012
Cuando dos de estas directrices parezcan entrar en conflicto, adopte aquella que sea
ms importante. Por otro lado se debe tener en cuenta que antes de seleccionar
cuales son las variables controladas y manipuladas para el sistema de control, por lo
cual se debe determinar el nmero de variables manipuladas permisibles. Esto se
hace a travs del anlisis de grados de libertad, que se determina usando un modelo
del proceso y de acuerdo a:
(. )
(. )
5
Ingeniera Qumica | Uniatlntico | 2012
1.3.2 Controlador
El controlador se considera como la mente maestra detrs del lazo de control, debido
a que este elemento toma la accin correctiva frente a un cambio en la variable de
salida ocasionado por una perturbacin en el proceso. Para que esta vuelva al valor
deseado, esta accin es enviada al elemento de control final. Generalmente los
controladores tienen dos modos de trabajar: uno es el manual donde el operario
decide que correccin debe hacerse, y el otro es el modo automtico en el que el
controlador es el que decide cmo actuar frente a la perturbacin.
Casi todos los controladores industriales emplean como fuente de energa la
electricidad o un fluido presurizado, tal como el aceite o el aire. Las seales que le
enva al elemento de control final pueden ser del tipo neumtico, hidrulico o
electrnico. Por ltimo, la escogencia de un controlador adecuado para una industria
en particular debe hacerse con base en la naturaleza de la planta y las condiciones
operacionales, incluyendo consideraciones tales como: seguridad, costo,
disponibilidad, confiabilidad, precisin, peso y tamao.
1.3.2.1
Los controladores poseen dos formas de accin y esto depende del tipo de naturaleza
de la variable de salida y del elemento de control final; hay que aclarar que si esto
no se escoge correctamente, el lazo de control nunca va a cumplir su funcin. Los
dos tipos de accin son la accin directa y la accin inversa.
Accin directa
Cuando la variable a controlar o de salida aumenta su valor por encima del valor
deseado o SP, es necesario que el controlador aumente el valor de la variable
manipulable o de entrada, por medio del elemento de control final, y as regresar la
variable de salida hasta su valor deseado y viceversa, es decir, cuando disminuye la
variable de salida por debajo del SP, el controlador por medio del elemento de
control final, disminuye el valor de la variable manipulada.
7
Ingeniera Qumica | Uniatlntico | 2012
Este puede ser el caso del control de nivel de lquido de un tanque, que cuando este
nivel aumenta por encima de su valor deseado, debe aumentar el flujo de salida de
lquido de tanque, para que regrese a su nivel normal.
Accin Inversa
Cuando el valor de la variable de salida aumenta por encima de su SP, se hace
necesario que el controlador actu sobre el elemento de control final, para que
disminuya el valor de la variable a manipular, y as poder disminuir el valor de la
variable de salida hasta su valor deseado y viceversa.
Un tpico ejemplo seria el control de flujo de un fluido caliente que pasa a travs de
los tubos de un intercambiador de calor. En caso que aumente demasiado la
temperatura del fluido fro, se hace necesario disminuir el caudal del fluido caliente,
para que este regrese a su temperatura deseada.
1.3.2.2
Tipos de controladores
8
Ingeniera Qumica | Uniatlntico | 2012
Siendo u(t) seal de salida del controlador que va hacia el elemento de control final
y e(t) es la seal del error definida como la deviacin entre el valor del SP y el valor
de la variable a controlar; para un controlador con accin de control proporcional la
relacin entre estos estara dada por:
( ) = ( )
(. )
( )
=
()
(. )
1.8
Cambio paso en el SP
Kp= 2
Kp= 10
Kp=100
1.6
1.4
1.2
1
0.8
0.6
0.4
0.2
0
-1
4
5
Tiempo (s)
10
( ) = ( ) + ( )
0
(. )
( )
1
= (1 + )
( )
(. )
10
Ingeniera Qumica | Uniatlntico | 2012
1.8
Cambio paso en el SP
Accin P
Accin PI
1.6
1.4
1.2
1
0.8
0.6
0.4
0.2
0
10
20
30
40
Tiempo (s)
50
60
70
11
Ingeniera Qumica | Uniatlntico | 2012
( ) = ( ) +
()
(. )
( )
= (1 + )
( )
(. )
1.8
Cambio paso en el SP
Accin P
Accin PI
Accin PD
1.6
1.4
1.2
1
0.8
0.6
0.4
0.2
0
10
20
30
40
Tiempo (s)
50
60
70
Figura 1.4. Comparacin entre la accin solo proporcional, proporcionalintegral y proporcional-derivativa para un cambio paso unitario en el SP
()
( ) = ( ) + ( ) +
0
(. )
( )
1
= (1 +
+ )
( )
(. )
1.8
Cambio paso en el SP
Accin P
Accin PI
Accin PID
1.6
1.4
1.2
1
0.8
0.6
0.4
0.2
0
10
20
30
40
Tiempo (s)
50
60
70
14
Ingeniera Qumica | Uniatlntico | 2012
( )
1
+ 1
)
= (1 + ) (
( )
+ 1
(. )
( )
1
)
= (1 +
+
( )
+ 1
(. )
Serie a Paralelo
Paralelo a Serie
Proporcional
= (1 )
4 0,5
) ]
= 0,5 [1 + (1
Integral
= +
4 0,5
) ]
= 0,5 [1 + (1
Derivativa
=
+
4 0,5
) ]
= 2 [1 + (1
15
Ingeniera Qumica | Uniatlntico | 2012
Es importante destacar que las ecuaciones para pasar de paralelo a serie solo son
Presin
Diafragma
Resorte
Vstago
Y
Plomada
Asiento
1.3.3.2
El otro tipo de falla, seria de falla abierta o fail-open (FO), la cual se abre, cuando
sucede algunos de los dos casos mencionados anteriormente y esta requiere de
energa para cerrarse o air-to-close (AC), un ejemplo de esta seria el lazo de control
de presin de un tanque, que es regulado por vapor que se encuentra en su interior,
cuando se va la energa o posee una falla interna el FCE, es necesario liberar todo el
vapor para evitar que el tanque explote por aumento de presin a causa del vapor
retenido, por lo cual esta vlvula debe ser de tipo FO o AC. Para explicar con mayor
claridad el tipo de accin de falla con que opera una vlvula de control se muestra en
la Figura 1.7 un modelo de cada una de ellas.
a)
b)
actuador hacia arriba, por lo cual esta quedar cerrada, lo que conlleva a decir que
esta vlvula es del tipo FC o AO.
Otro punto a considerar es la accin de falla de las vlvulas de control y como est
como afecta la accin de controlador, por lo cual se usara el siguiente ejemplo:
Suponga que se piensa controlar el nivel de un tanque con la corriente de salida del
mismo, en teora lo que hace el lazo es que si el nivel aumenta por encima de su SP,
debe salir ms liquido del recipiente para que el nivel vuelva a la normalidad. Ahora
si a esta situacin se le agrega una vlvula de tipo FC, la cual abrir al aumentar la
seal del controlador, eso significa que si el nivel aumenta el controlador mandara la
seal y la vlvula abrir y realizar la accin correcta para el lazo, es decir que el
controlador es de accin directa. Sin embargo, si se usa una vlvula de tipo FO, esta
al llegarle seal del controlador lo que har es cerrarse y esta accin ser contraria a
lo que el lazo necesita hacer, por eso la accin del controlador debe ser inversa.
Con respecto un controlador on/off se puede realizar este mismo anlisis y generara
los mismos resultados, es decir si la vlvula es de FC y el nivel de lquido aumenta,
el controlador debe encenderse cuando el valor de la variable controlada supera el
lmite superior de tolerancia ingresada al controlador y por ende abrir
completamente la vlvula hasta que el nivel alcance el lmite inferior de tolerancia,
momento en que se apagar el controlador y por ende se cerrara el FCE. Si la
vlvula de control es de FO el controlador se encender cuando el valor del nivel de
lquido alcance el lmite inferior de tolerancia y por consiguiente cerrar por
completo la vlvula y este controlador se apagar cuando el nivel alcance el lmite
superior de tolerancia y la vlvula se abrir por completo. En conclusin y para que
el lector no se confunda, en este documento la accin del controlador se estimara
como si la vlvula fuese de falla cerrada, porque es el mtodo ms sencillo de
asimilar y adems esta es la accin de falla que tienen por defecto las vlvulas en
Aspen-HYSYS y si el lector quiere estimar la accin del controlador usando un
FCE de falla abierta es la accin contraria a cuando se usa uno de falla
cerrada.
Fisher (10)
Ahora como parte de Emerson Electric Co., es el lder mundial en el suministro de
confiables vlvulas de control, reguladores y servicios de instrumentacin y
desempeo a las industrias de control de procesos. Fisher lleva a cabo
innovadoras soluciones diseadas para ayudar a los clientes a reducir los costos de
mantenimiento de planta, reducir los requerimientos de capital, reducir los costos de
cumplimientos legales y aumentar la disponibilidad del proceso.
1.3.3.4
(. )
(. )
(. )
= 500
(. )
Las ecuaciones de 1.12 hasta 1.15, son las ecuaciones bsicas para el diseo de
vlvulas, cabe decir que de estas ecuaciones, los creadores de vlvulas de control,
las modifican para el diseo de sus propias vlvulas y la creacin de sus catlogos y
tablas de seleccin. Adems del clculo de flujo de lquido, a travs de las vlvulas
de control, tambin existen modelos matemticos, generados por los fabricantes de
vlvulas, para predecir el flujo de corrientes gaseosas, fluidos compresibles, flujo
bifsico, flujo de vapor, que pasan a travs de ellas. Generalmente estas ecuaciones
realizan sus clculos en SCFH (standard cubic feet per hour) o clculo de flujo
volumtrico a condiciones estndar a 14,7 psia (1 atm) y 60 F (15,556 C) y otros
parmetros especficos, dados por el fabricante.
1.3.3.6
23
Ingeniera Qumica | Uniatlntico | 2012
24
Ingeniera Qumica | Uniatlntico | 2012
(. )
(. )
Debido a que existen tres tipos de caractersticas, para que una vlvula regule el
flujo, deben existir criterios que permitan la seleccin de mejor opcin. Una forma
prctica de hacerlo serian por medio de reglas de dedo gordo o heursticas que son
fruto de la experiencia, estas dictan las siguientes reglas:
Use vlvulas con caracterstica lineal e en procesos lineales cuando la cada de
presin a travs de ellas no cambia con el flujo, es decir que los controles de
flujos y nivel de lquido se ajustan a esta situacin.
Use vlvulas con caracterstica de igual porcentaje se utilizan, generalmente,
cuando la cada de presin a travs de la vlvula vara con el flujo, y cuando se
encuentre frente a procesos no lineales como control de presin de un equipo
(tanques, separadores de fases, etc.) o controles hechos a estructuras complejas
como reactores o columnas de destilacin, por eso esta es la caracterstica de las
vlvulas de control ms usada en la industria.
Use vlvulas con caractersticas de abertura rpida, solo cuando el lazo de
control requiera una gran cantidad de flujo, con solo un pequeo porcentaje de
abertura, como los casos de control de flujo, aunque es menos usado, que l con
caracterstica lineal.
1.3.3.7
Debido a que la vlvula de control, es el nico elemento del lazo de control que
interacta directamente con el proceso, es importante saber la ubicacin estratgica
que debe poseer este elemento de control final, por lo cual se usaran
recomendaciones o heursticas para su proceso de instalacin (15):
Instale las vlvulas de control en la corriente de descarga y no en la de succin
una bomba centrifuga.
Use una sola vlvula de control en una lnea de lquido.
Nunca instale una vlvula de control en la corriente de descarga de un
compresor.
26
Ingeniera Qumica | Uniatlntico | 2012
1.4
Controladores de flujo
La dinmica de los lazos en los que se involucran controladores de flujo suelen ser
rpidas, estos incluyen generalmente un sensor tipo orificio, un transmisor de
diferencial de presin, un controlador y una vlvula de control. El tiempo que
transcurre para obtener la respuesta en el actuador de las vlvulas es bastante corto,
a menos que se tengan vlvulas muy grandes. Por otro lado, la ganancia que se
obtiene en la respuesta ha de ser no muy grande dado que cambios muy bruscos
pueden ocasionar flujos an ms turbulentos, llevando as a tener mediciones poco
27
Ingeniera Qumica | Uniatlntico | 2012
precisas en el medidor que se use. Los valores tpicos de cada uno de los parmetros
seran:
Ganancia Kp = 0,5
Tiempo integral = 0,3 minutos
Tiempo derivativo = no suele ser necesario
Controladores de nivel
Para este tipo de control, generalmente se usa un controlador proporcional, a pesar
de que esto implique tener siempre un off-set, pues tener cierto nivel de lquido no
es ms crtico que alcanzar el estado estacionario, siempre y cuando el nivel de
lquido se encuentre dentro de los lmites tolerables para el sistema. El valor tpico
para este controlador seria:
Ganancia Kp = 2
Para estos casos se requieren que el nivel de lquido posea el menor off-set posible,
por ende recomienda aumentar el valor Kp de 2 a 10.
Controladores de presin
De igual manera, que los otros 2 controladores, estos son fcilmente sintonizados.
La variable de tiempo del sistema se calcula haciendo una razn entre el volumen de
gas del sistema y el flujo volumtrico de gas a travs del sistema. Donde entonces el
tiempo integral es de 2 a 4 veces el tiempo del sistema y junto con una ganancia
razonable se obtienen resultados satisfactorios. Los valores recomendados son los
siguientes:
28
Ingeniera Qumica | Uniatlntico | 2012
Ganancia Kp = 2
Tiempo integral = 2 - 10 minutos
Tiempo derivativo = no suele ser necesario
Referencias bibliogrficas
[Cited:
12,
2012.]
12. . Principles and Practice of Automatic Process Control. Tercera edicin. s.l. :
Jhon Wiley & Sons, Inc, 2006. pp. 164-170.
13. Fisher. Control Valve Handbook. Cuarta edicin. Marshalltown : Emerson
Process Management, 2005. pp. 58-60 .
14. . Control Valve Handbook. Cuarta edicin. Marshalltown : Emerson Process
Management, 2005. pp. 109-111.
15. Willian L. Luyben. Plantwide Dynamic Simulators in Chemical Processing and
Control. New York : Marcel Dekker,Inc, 2002. pp. 17-20.
16. . Plantwide Dynamic Simulatos in Chemical Processing and Control. New
York : Marcel Dekker, Inc, 2002. pp. 26-28.
30
Ingeniera Qumica | Uniatlntico | 2012
Otra variable
de entrada
Seal de salida
del
controlador
Seal del
error
SP
Controlador
Elemento de
control final
(vlvula de
control)
Otras
Perturbaciones al
proceso
Variable de
entrada
Proceso
Variable de
salida
Sensor/Transmisor
31
Ingeniera Qumica | Uniatlntico | 2012
0.2
a)
0.1
variable de salida
variable de salida
10
20
1.5
b)
1
0.5
0
30
10
tiempo (s)
20
30
tiempo (s)
variable de salida
2
Pu
1.5
c)
1
0.5
0
10
15
tiempo (s)
20
25
30
20
25
30
variable de salida
100
d)
50
0
-50
10
15
tiempo (s)
Las figuras 2.2a y 2.2b, se pueden considerar estables porque el valor de la variable
de salida se vuelve constante despus de haber pasado un tiempo considerable, o en
otras palabras estas respuestas pudieron alcanzar el estado estacionario; por su parte,
la Figura 2.2d es una respuesta inestable, porque no es capaz de alcanzar el estado
estacionario. Por ltimo, la Figura 2.2c se puede considerar un caso intermedio
entre la estabilidad y la no estabilidad del lazo, debido a que esta respuesta es
oscilatoria de amplitud contante y si se aumenta un poco la ganancia del controlador,
la respuesta se vuelve inestable, mientras que si se disminuye un poco este valor la
respuesta se vuelve estable. A esta ganancia, que limita la estabilidad del lazo de
control, se le conoce como ganancia ltima Kpu, y el tiempo entre las crestas en
cada oscilacin se le conoce como periodo ltimo Pu (dado en min o en segundos)
o tambin se puede expresar como frecuencia ltima wu, dado por la siguiente
frmula.
2
)
(
(. )
32
Ensayo y error
Mtodo de sustitucin directa
Prueba de Routh
Lugar de las races (Root Locus)
Diagrama de Bode y Nyquist
ATV (Auto Tuning Variation)
De todos los mtodos aqu presentados, solo se explicar el mtodo ATV, porque los
otros requieren del planteamiento de un modelamiento matemtico, como la
dinmica del proceso estudiado afecta a la variable de salida, lo cual en la mayora
de procesos industriales es difcil de realizar por la complejidad que ellos poseen,
mientras que con este mtodo no es necesario. Otra razn es porque es un mtodo no
necesita mucha fundamentacin matemtica para su comprensin.
2.2.1.1
Mtodo ATV
El mtodo de ATV (18) (19) (20) o Auto Tuning Variation, es un procedimiento sencillo
para la estimacin de la ganancia ltima y el periodo ltimo, el cual viene como
algoritmo interno de controladores PID comerciales y es usado cuando el
controlador se encuentra modo manual. Este mtodo parte del mismo principio que
posee el controlador on/off o relay, que consiste que la seal de salida del
controlador solo posee dos valores uno mximo y otro mnimo, representando con
las siguientes ecuaciones:
< 0
(. )
> 0
(. )
Variable de salida
estacionario (ideal), el cual se designa como h. Esto significa, que cuando la seal
del controlador es mnima, posee un valor igual a h, y cuando es mxima es igual a
+h. La respuesta grafica para un cambio paso en el SP de este mtodo, con un valor
de h igual a la unidad, se puede apreciar en la Figura 2.3.
10
5
0
-5
Set-Point de 0.96
-10
10
20
30
40
50
60
70
80
90
100
70
80
90
100
Variable de entrada
Tiempo (s)
1
0.5
0
-0.5
-1
0
10
20
30
40
50
60
Tiempo (s)
34
Ingeniera Qumica | Uniatlntico | 2012
(. )
[ ()]
=
2
(. )
@ (2.4a) =
@ (2.4b) =
4(1)
= ,
(5,8524)
(1)
= ,
(2 )( 5,8524)
= 5 4
= 52,45 39,66 = ,
35
Ingeniera Qumica | Uniatlntico | 2012
2
= ,
12,79
Root Locus
1.5
System: untitled1
Gain: 0.207
Pole: -0.000922 + 0.497i
Damping: 0.00186
Overshoot (%): 99.4
Frequency (rad/sec): 0.497
Imaginary Axis
0.5
-0.5
-1
-1.5
-2
-1.5
-1
-0.5
0.5
Real Axis
2.2.1.2
valor del ruido del proceso. Al agregar este nuevo trmino a las ecuaciones 2.2 y
2.3, cambia el criterio en el cual la seal del controlador es mxima o mnima de la
siguiente manera:
.
<
(. )
> +
(. )
= 0,05(2) = 0,1
Variable de salida
La representacin grfica del mtodo ATV sin banda muerta y con banda muerta, se
presenta en la Figura 2.5
10
5
0
Set-Point
Sin banda muerta
Con banda muerta
-5
-10
10
20
30
40
50
60
70
80
90
100
Variable de entrada
Tiempo (s)
1
0.5
0
-0.5
-1
0
10
20
30
40
50
60
70
80
90
100
Tiempo (s)
Figura 2.5 Oscilaciones usando el mtodo ATV, con y sin banda muerta
37
Ingeniera Qumica | Uniatlntico | 2012
Como se muestra en la Figura 2.5, las oscilaciones con banda muerta, estn
desfasadas con respecto a la que no poseen, creando un efecto de histresis
rectangular. Cabe anotar acerca de las lneas que poseen banda muerta, el cambio de
la variable de entrada de 1 a -1 se hace en 0,86 (0,96-0,1= 0,86) y de -1 a 1 en 1,06
(0,96+0,1=1,06), lo cual confirma el desfase de las dos respuestas.
A pesar que esta mejora al mtodo ATV soluciona el problema con respecto al ruido
de proceso, trae como consecuencia, y es debido al desfase producido por la banda
muerta, que ocurre un aumento en la amplitud entre las oscilaciones del valor de la
variable de salida, con respecto al mtodo sin banda muerta, y por consiguiente
genera un clculo no tan exacto de la ganancia y periodo ultimo usados en la seccin
2.2.1.1.
Una posible aproximacin, tomando como referencia la ecuacin 2.4 b sera:
( 2
+ 4 2 )0,5
(. )
Al usar esta ecuacin, la ganancia ultima seria 0,2683, valor casi igual al calculado
con la ecuacin 2.4b, lo que no generara ningn error usar las ecuaciones 2.4a y
2.4b si el valor de la banda muerta es muy pequea.
Aspen-HYSYS posee este mtodo, para la obtencin de los parmetros
dinmicos para controladores PI o PID (22). Este mtodo posee las siguientes cinco
especificaciones:
Alpha ()
Phi ()
Es el ngulo de fase y es entendido, como cuantos radianes o grados est atrasada
(valor negativo) o adelantada (valor positivo) la respuesta de la variable de salida,
con respecto a la de entrada (esto tiene ms significado viendo la teora del control
desde el dominio de la frecuencia). Este valor se toma como referencia y calculara la
ganancia a este valor de . El uso de esta ganancia y no la ganancia ltima, es
debido a que internamente el simulador tiene es el algoritmo de un controlador PID,
38
Ingeniera Qumica | Uniatlntico | 2012
Beta ()
Tambin es conocido como margen de ganancia, el cual es definido como o el
cociente entre la ganancia obtenida al valor del ngulo de fase seleccionado y la
ganancia real del controlador ( ( ) ). Los lmites de este
parmetro son de 0,1 hasta 1 y el valor que viene por defecto es 0,25.
Ganancia
Proporcional
Kp
Tiempo
Integral
Tiempo
Derivativo
Kpu/2
PI
Kpu/2,2
Pu/1,2
PID serie
Kpu/1,7
Pu /2
Pu /8
PID paralelo
0,75Kpu
0,625Pu
Pu /10
40
Ingeniera Qumica | Uniatlntico | 2012
Ganancia
Proporcional
Kp
Tiempo
Integral
Tiempo
Derivativo
PI
Kpu/3,2
2,2Pu
PID paralelo
Kpu/2,2
2,2Pu
Pu /6,3
PI
PID serie
PI
PID paralelo
Ganancia
Proporcional
Kp
Ziegler - Nichols
0,0989
0,1280
Tyreus - Luyben
0,0680
0,0989
Tiempo
Integral
Tiempo
Derivativo
10,6583
6,3950
1,5987
28,1380
28,1380
2,0302
41
Ingeniera Qumica | Uniatlntico | 2012
Variable de salida
PI con ZN
PI con TL
a)
1.5
1
0.5
0
10
20
30
40
50
60
70
Variable de salida
Tiempo (s)
2
PID con ZN
PID con TL
b)
1.5
1
0.5
0
10
20
30
40
50
60
70
Tiempo (s)
Este mdulo solo trata tcnicas de sintonizacin en lazo cerrado, sin embargo
tambin existen las tcnicas de sintonizacin en lazo abierto, y aunque el simulador
Aspen-HYSYS, permite hacer ambas metodologas, la sintonizacin en lazo
cerrado, es ms amigable con la interfaz que provee el simulador., Sin embargo, en
un mdulo posterior se trataran las tcnicas de sintonizacin en lazo abierto.
Referencias bibliogrficas
42
Ingeniera Qumica | Uniatlntico | 2012
43
Ingeniera Qumica | Uniatlntico | 2012
Variable de Entrada
Variable de Salida
Funcin de
Transferencia
Medicin
f, m^3/s
Ti(t), C
Medicin
f, m^3/s
T(t), C
E-1
44
Ingeniera Qumica | Uniatlntico | 2012
=
( ) ( ) =
( )
(. )
Puesto que en control, lo que interesa saber en cuanto esta desviado la variable de
salida con respecto de su valor deseado, se usa las variables desviacin que es
simplemente la resta entre la variable en modo dinmico y la misma en estado
estacionario (la cual se representa con una barra sobre la variable). Entonces la
ecuacin en estado estacionario es la siguiente:
=0
(. )
) (( ) ) =
( ( )
(( ) )
(. )
( )
(. )
( ) ( ) =
Donde
)
( ) = ( ( )
(. )
( ) = (( ) )
(. )
45
Ingeniera Qumica | Uniatlntico | 2012
( )
+ ( ) = ( )
(. )
Donde
=
[3 ][3 ][ ]
=
[3 ]3
(. )
3 [3 ]
=
=
=
3 [3 ]
(. )
() (0) + () = ()
(. )
Puesto que se est usando las variables desviacin, la condicin inicial es igual a
cero ((0) = 0) porque las variables en estado dinmico son iguales a las del
estado estacionario. Entonces la ecuacin 3.10 queda:
() + () = ()
()
=
() +
(. )
(. )
( ) = (( ) )( ) = () = ()()
(. )
() =
1
=
(. )
() = ()( )
(. )
Usando los siguientes datos para obtener el valor del atraso dinmico.
Cp = Cpi = 4,184 Kj/Kg K
Cv= 3,7221 Kj/Kg K
= =1000 Kg/ m3
V = 10 m3
f = 0,5 m3/s
10 1000 3,7221 [3 ][3 ][ ]
=
= ,
0.5 1000 4,184 [3 ]3
47
Ingeniera Qumica | Uniatlntico | 2012
Con el valor del atraso dinmico, el valor de K (es igual a uno) y la ecuacin 3.15 se
hace una grfica de los valores de las variables de salida y entrada con respecto al
tiempo, el cual se puede apreciar en la Figura 3.3
0.9
Respuesta ltima
0.8
0.6321 de la respuesta ltima
Respuesta
0.7
0.6
0.5
0.4
0.3
0.2
Constante de tiempo
0.1
0
10
20
30
40
50
60
Tiempo (s)
80
90
100
Como se observa en la Figura 3.3, mientras que a aplicarle un cambio paso unitario
a la variable de entrada es instantneo, la variable no responde instantneamente
debido al atraso dinmico que posee la funcin de transferencia. Esto se debe a que
si en la ecuacin 3.15, cuando el tiempo es igual al atraso dinmico, este apenas ha
alcanzado el 63,21 % de la respuesta ultima ( () = (1)(1 1 ) = (1)0,6321)
que es igual a la unidad y que la respuesta ltima se alcanza cuando el tiempo es
aproximadamente 5 veces el valor de ( () = (1)(1 5 ) = (1)0,9932).
Otra apreciacin que se puede hacer al atraso dinmico, mencionada al principio de
esta seccin, es que este parmetro depende de la dinmica del sistema, es decir que
si se manipula alguna de las variables de proceso, estas pueden favorecer o no, a que
una perturbacin en la variable de entrada afecte ms o menos rpido la variable de
salida. Un ejemplo seria aumentando de valor de flujo de alimentacin del tanque a
0.7 m3/s (= 12.708 s) o disminuyndolo a 0.3 m3/s (= 29.653 s) y se usa la
conclusin, ya deducida, de que la respuesta ltima se alcanza cuando el tiempo es 5
veces el valor de entonces:
48
Ingeniera Qumica | Uniatlntico | 2012
3
@ 0.7
= 12,708 (5) = ,
3
@ 0.5
= 17,792 (5) = ,
@ 0.3
3
= 29,653 (5) = ,
Esto nos lleva a la conclusin, de entre ms pequeo el valor del atraso dinmico
sea, ms instantnea ser la respuesta de la variable de salida frente a una
perturbacin hecha de la variable de entrada y viceversa.
El ejemplo utilizado anteriormente, presenta una dinmica de primer orden, sin
embargo, es posible tener funciones de transferencia de rdenes superiores, que
pueden definirse mediante las siguientes funciones de transferencia donde Y(s)
representa la variable de salida y X(s) la variable de entrada (todo se encuentra en
variables desviacin):
Segundo orden
Tercer orden
N orden
()
()
()
()
()
()
1
(1 +1)(2 +1)
1
(1 +1)(2 +1)(3 +1)
= =1
1
( +1)
(. )
(. )
(. )
49
Ingeniera Qumica | Uniatlntico | 2012
1
0.9
0.8
Respuesta
0.7
0.6
0.5
0.4
0.3
Variable de entrada
V.S. con una dinmica de primer orden
V.S. con una dinmica de segundo orden
V.S. con una dinmica de tercer orden
0.2
0.1
0
10
15
Tiempo (s)
Medicin
f, m^3/s
Ti(t), C
Medicin
f, m^3/s
T(t), C
Con esta condicin, debe existir un tiempo en variable de entrada, el cual al hacerle
perturbacin sobre esta, no afectara la variable de salida, sino despus que este haya
transcurrido (representado por el tramo de tubera). Matemticamente se aprecia de
la siguiente forma:
( ) = ( )
(. )
() = [ ( )] = 0 ()
(. )
()
=
() +
(. )
() = ()( ) (
( )
)
(. )
51
Ingeniera Qumica | Uniatlntico | 2012
1
0.9
0.8
Respuesta
0.7
0.6
0.5
0.4
0.3
0.2
0.1
tiempo muerto to
0
0
10
20
30
50
60
Tiempo (s)
70
80
90
100
( + 1)
()
=
()
( + 1)
(. )
52
Ingeniera Qumica | Uniatlntico | 2012
funcin de transferencia con dos atrasos y otra que adems posee un adelanto
dinmico.
()
(0.8 + 1)
=
() ( + 1)( + 1)
(. )
()
1
=
() ( + 1)( + 1)
(. )
1
0.9
0.8
Respuesta
0.7
0.6
0.5
0.4
0.3
0.2
Variable de entrada
Variable de salida con adelanto dinmico
Variable de salida sin adelanto dinmico
0.1
0
0
10
15
Tiempo (s)
53
Ingeniera Qumica | Uniatlntico | 2012
()
=
()
=1( + 1)
=1( + 1)
(. )
para las mediciones de estos lazos de control, como se muestra en la Tabla 3.1. Con
respecto a los elementos de control final, no hay problema con el uso de vlvula con
accin instantnea, porque este tipo de accin es normalmente representada por el
uso de vlvulas neumticas, las cuales de uso comn en la industria.
Tipo
de
Medicin
Ganancia
Constante
de
Tiempo
(min)
Tipo
de
Atraso
Lquido
0.5
Atraso
de
Primer orden
Vapor
Atraso
de
Primer orden
3 a 10
Tiempo
muerto
Temperatura
Composicin Cromatografa
Referencias bibliogrficas
55
Ingeniera Qumica | Uniatlntico | 2012
Otra variable
de entrada
Seal del
error
SP
+
-
Controlador
Seal de salida
Variable de
del
entrada
controlador Elemento de control
final (vlvula de
control)
Otras
Perturbaciones al
proceso
+
Proceso
Variable de
salida
Sensor/Transmisor
Para que la Figura 4.1 se convierta en lazo abierto literalmente se debe hacer un
corte de las seales de salida del controlador, con el resto de los elementos que
posee el lazo y la seal del sensor/transmisor al controlador, tal como se muestra en
la Figura 4.2.
La ventaja de trabajar con la Figura 4.2 para la sintonizacin de controladores, es
que no se hace necesario llevar a la variable de salida (controlada) hasta su valor
ltimo (respuesta oscilatoria de amplitud constante), situacin que no es posible
realizar en muchos procesos por cuestiones de seguridad, puesto que conociendo las
dinmicas del elemento de control final y del sensor/transmisor, se pueden estimar
los parmetros dinmicos que requiere el controlador para actuar frente a la
dinmicas de ellos.
56
Ingeniera Qumica | Uniatlntico | 2012
Otra variable
de entrada
Cambio
Paso
Seal del
error
SP
+
-
Controlador
Elemento de control
final (vlvula de
control)
Otras
Perturbaciones al
proceso
Variable de
entrada
+
Proceso
Variable de
salida
Sensor/Transmisor
Seal del
Sensor/transmisor
() = () () () =
+
(. )
57
Ingeniera Qumica | Uniatlntico | 2012
/
=
(. )
= 1.5(2 1 )
(. )
= 2
(. )
() = () () =
0.9
(3 + 1)( + 1)(0.5 + 1)
Valor ltimo
0.8
Respuesta
0.89961
0.6
(t2=4.6098 , 0.623ST=0.56046)
0.4
(t1=2.4143 , 0.283ST=0.25459)
0.2
0
0
10
15
20
25
Tiempo (s)
0,8996 0
=
= ,
10
= 1.5(4,6098 2,3142) = .
= (4,6098 3.2933) = ,
() =
, ,
. +
0.9
0.8
0.7
Respuesta
0.6
exp(-1.3165to)
FOPDT= 0.89961
-------------------------------3.2933s + 1
0.5
0.4
0.3
0.2
0.1
0
10
15
20
25
Tiempo (s)
valores de dentro del siguiente rango: 0,1 0,5 y para el control PID se usar la
PI
Ganancia
Proporcional
Kp
Tiempo
Integral
Tiempo
Derivativo
1 1
( )
0.9 1
( )
3.33
PID serie
1.2 1
( )
0.5
PID paralelo
1.5 1
( )
2.5
0.4
Una comparacin grafica entre estas dos tcnicas para controladores PI y PID (con
un valor de de 0.1), y usando el ejemplo de la seccin 4.2 se puede apreciar en la
Figura 4.5 y los parmetros dinmicos en la Tabla 4.3 para un cambio paso unitario
en el SP. Como se puede apreciar en la figura y en la tabla, las dependencias de las
acciones integral y derivativa, con el atraso dinmico del FOPDT, son las que
genera las discrepancias entre las tcnicas de ZN y CC, generando esta ltima,
valores con ms sobresalto y por consiguiente respuestas potencialmente inestables.
61
Ingeniera Qumica | Uniatlntico | 2012
Ganancia
Proporcional
Kp
Tiempo
Integral
Tiempo
Derivativo
(1 + )
PI
(0.9 +
)
12
3
)
20
(9 +
)
PD
1 5
( + )
4 6
2
)
3
(22 + )
PID paralelo
1 4
( + )
3 4
6
)
8
(13 + )
4
2
(11 + )
(30 +
(6
(32 +
Ganancia
Proporcional
Kp
Tiempo
Integral
Tiempo
Derivativo
Ziegler - Nichols
PI
2,5024
4,3839
PID serie
3,3366
2,6330
0,6583
PI
PID paralelo
Cohen - Coon
2,5951
2,4167
3,9852
2,7957
0,3251
62
Variable de salida
1.5
a)
0.5
PI con ZN
PI con CC
0
10
Tiempo (s)
15
20
25
Variable de salida
1.5
b)
0.5
PID con ZN
PID con CC
0
10
Tiempo (s)
15
20
25
63
Ingeniera Qumica | Uniatlntico | 2012
= Proceso ( () = () () ())
U(s)
Y(s)
R(s)
Rm(s)
-
Q(s)
D(s)
U(s)
G(s)
P-32
+
Y(s)
Proceso
Controlador IMC
Gm(s)
Ym(s)
Dm(s)
Proceso realizado internamente por el
controlador IMC
Figura 5.1 Estructura de un IMC (todas las variables estn como variables
desviacin)
() = () () = ()( () ()) + ()
() = () () = () () + ()
(. )
(. )
De la Figura 5.1 se le puede hacer el arreglo mostrado en la Figura 5.2. Cabe notar
que las Figuras 5.1 y 5.2 generan iguales resultados, aunque su diagrama de bloques
sea diferente, sin embargo en la Figura 5.2 se puede hacer el arreglo mostrado en la
Figura 5.3, si lo delimitado por las lneas punteadas de la figura se convierte en un
solo bloque, eso significara que ese bloque fuese equivalente la funcin de
65
Ingeniera Qumica | Uniatlntico | 2012
Controlador Convencional
Gc(s)
D(m)
R(s)-Y(s)+Ym(s)
E(s)=R(s)-Y(s)
R(s)
Q(s)
U(s)
G(s)
Ym(s)
Y(s)
Gm(s)
E(s)=R(s)-Y(s)
R(s)
+
D(m)
Gc(s)
Y(s)
U(s)
G(s)
Controlador basado en
IMC
Para deducir Gc(s), se usar el lgebra del diagrama de bloques y la Figura 5.2
() = ()[ () + () ()]
(. )
Despejando la ecuacin (5.4) para obtener el cociente U(s)/E(s) que es igual a Gm(s)
()
()
= () =
()
() ()
(. )
66
Ingeniera Qumica | Uniatlntico | 2012
(. )
(. )
( ) =
1
( + 1)
(. )
1 ( )
() = ()() =
()
(. )
67
Ingeniera Qumica | Uniatlntico | 2012
2 + + 1
1
)(
)
() = (
+ 1
() =
+ 1
(. )
0,5 + 1
0,5 + 1
(. )
() =
(0,5 + 1)
=
+ 1 (0,5 + 1)( + 1)
(. )
68
Ingeniera Qumica | Uniatlntico | 2012
() =
(0,5 + 1)( + 1)
(. )
+ () = (0,5 + 1)
(. )
1 ( )
() =
=
(0,5 + 1)( + 1)
(. )
Y el controlador de modo interno ideal + el filtro de bajo paso, con n igual uno, sera
igual a
() = () () =
(0,5 + 1)( + 1)
1
(
)
+ 1
(. )
() =
( )
() ()
=
1 () () 1 () () ()
1 ( ) ( )
() =
1 () ()
1 ()+ ()
(. )
1 ( ) ( )
1 ( )
() =
=
1 + ()() 1
( )
() + +
1
() = ( )
(. )
(0,5 +1)(+1)
+1(0,5 +1)
69
Ingeniera Qumica | Uniatlntico | 2012
(, +)(+)
(, +)
() = ( )
(. )
1 [0,5 2 + ( + 0,5 ) + 1]
() = ( )
(0,5 + )
(. )
() = [
+ 0,5
0,5 2 + ( + 0,5 ) + 1
][
]
(0,5 + )
( + 0,5 )
(. )
() = (1 +
+ ) = (
2 + +1
+ ,
(, + )
= + ,
= 0,5
70
Ingeniera Qumica | Uniatlntico | 2012
Entonces,
,
+,
Para este caso en la bibliografa citada (74), se recomienda usar > 0.8 . El lector
debe notar que para este caso se opt el uso de una aproximacin de Pad de primer
orden, pero si se hubiese usado series de Taylor de primer orden ( + 1),
el resultado hubiese sido un controlador PI con los siguientes parmetros.
= ( )
+
Lo que lleva a pensar que entre ms o menos robusto se quiere que sea el
controlador as sern las ecuaciones generadas para los parmetros dinmicos.
Generalmente eso se logra usando una ecuacin de filtro mejorado y que n se ajuste
de tal forma que el orden del numerador sea igual que el denominador (funcin semi
impropia) y no un grado mayor, usar el artificio matemtico de factorizacin all
pass o aumentado el orden de aproximacin de Pad, aproximacin de Taylor o
inclusive no usar ninguna de las dos, sin embargo en este documento no se hondara
sobre esto, pero el lector puede revisar las bibliografa de este mdulo para poder
profundizar en este tema. (Nota: Cuando se refiere a un control robusto hace
referencia a un control que no se ve menos afectado por las perturbaciones del lazo).
Por ltimo en la Tabla 5.1 se resume el diseo de los ms representativos casos de
modelo de proceso hay que mencionar que para los modelos de proceso sin tiempo
muerto, uno puede aumentar o disminuir a gusto del lector, pero debe tener en
cuenta que entre ms grande es el valor ms robusto ser el controlador, sin embargo
cuando este posee tiempo muerto se debe aumentar desde un lmite inferior
recomendado.
De la tabla presentada, solo el ejemplo que intencionalmente se demostr en este
mdulo y est en negrilla de la tabla puede ser aplicado directamente en AspenHYSYS conociendo los parmetros de ajuste de FOPDT, sin embargo el lector
puede con otras tcnicas de ajuste pero tiene que realizar los clculos de los
parmetros dinmicos manualmente.
71
Ingeniera Qumica | Uniatlntico | 2012
Kp
( + )( + )
1 + 2
1 + 2
1 2
1 + 2
( + + )
( + + )
( + )
+ 0,5
(0,5 + )
+ 0,5
2 +
> 0.8
Pad 1er orden
No se us Pad
> 1,7
( )
+
Serie de Taylor de
primer orden
+ 0,5
( + )
+ 0,5
2 +
2( + )
+ 2
( + 2)2
+ 2
2
( + 0,5 )
+ 2
+ 0,25
+ 2
> 0,2
( + 2)
( + 2)
2
6
+ 4
22
Recomendaciones
() = 2
() = 2
Factorizacin
> 0,25
Serie de Taylor de
primer orden
> 0,2
Referencias bibliogrficas
67. Carlos A. Smith, Armando B. Corripio. Principles and Practice of Automatic
Process Control. Tercera edicin. s.l. : Jhon Wiley & Sons, Inc, 2006. pgs. 234239.
68. B. Wayne Bequette. Process Control: Modeling, Design and Simulation. Upper
Saddle River,NJ : Prentice Hall, 2003. pgs. 129-135.
72
Ingeniera Qumica | Uniatlntico | 2012
73
Ingeniera Qumica | Uniatlntico | 2012
SP
TR
TC1
TH SP
TT1
TH
TT2
Horno
TC2
TH
Agua de Enfriamiento
Alimento
Reactor
Producto
FC
Aire
Combustible
R(s)
E(s)
CO1(s)
CO2(s)
Taire (s) C
Dinmica del
flujo de aire
gpm
Controlador
primario
Controlador
Secundario
Elemento de
control final
TH (s) C
T(s) C
Dinmica del
flujo de
combustible
Dinmica del
Reactor
%TO2(s)
%TO1(s)
Dnde:
R(s) = Set-Point del lazo de control primario.
E(s) = Seal de error que entra al controlador maestro.
CO1(s) = Seal de salida del controlador primario (maestro), que sera el Set-Point
del lazo de control secundario.
CO2(s) = Seal de salida del controlador secundario (esclavo).
TO1(s) = Seal de salida del transmisor del lazo de control primario.
TO2(s) = Seal de salida del transmisor del lazo de control secundario.
TH(s) = Temperatura del fluido antes de entrar al reactor.
76
Ingeniera Qumica | Uniatlntico | 2012
primario disminuye el SP del control secundario, entonces Tc2 quedara con un valor
por encima de su SP y como su accin es tambin inversa, obliga que la abertura de
la vlvula que regula el flujo de combustible disminuya. Todo lo mencionado
anteriormente se puede resumir de esta forma: el control maestro usa el sistema de
control regulatorio (SP constante perturbaciones variable), mientras que el control o
controles esclavos usa el sistema de servo control (SP variable, perturbaciones
constante), que para ambos casos la accin de controlador depende de gran medida
del elemento de control final que posee lazo de control en cascada.
En el ejemplo presentado, se mostr un control en cascada de segundo nivel, pero
tambin pueden existir de niveles superiores, sin embargo solo es recomendable
hasta el tercer nivel, puesto que si se usan rdenes superiores, la configuracin de
los lazos de control se vuelve ms compleja, que la dinmica del proceso que obligo
el uso de esta estrategia de control.
(. )
Referencias bibliogrficas
37. Carlos A. Smith, Armando B. Corripio. Principles and Practice of Automatic Process Control.
Tercera edicin. s.l. : Jhon Wiley & Sons, Inc, 2006. pp. 310-315.
38. B. Wayne Bequette. Process Control: Modeling, Design and Simulation. Upper Saddle River,
NJ : Prentice Hall, 2003. pp. 314-317.
78
Ingeniera Qumica | Uniatlntico | 2012
39. AspenTech. Aspen HYSYS Dynamics, user guide. Burlington, Ma : Aspen Technology, Inc., 2009.
pp. 3-28 - 3.29.
40. Willian Y. Svrcek, Donald P. Mahoney, Brent R. Young. A Real Time Approach to Process
Control. Segunda edicin. Chichester : John Wiley & Sons. Ltd, 2006. p. 134.
79
Ingeniera Qumica | Uniatlntico | 2012
SP
Perturbacin
Controlador
Feedforward
Proceso
Salida del
controlador
CO%
a)
Perturbacin
SP
Controlador
Feedback
Salida del
controlador
CO%
Proceso
Variable
de salida
b)
LT
SP
LC
Controlador
Feedback
Vapor
Tanque
Hervidor
Agua de
alimento
Gas
Caliente
81
Ingeniera Qumica | Uniatlntico | 2012
FFC
Controlador
Feedforward
FT
Vapor
Tanque
Hervidor
Agua de
alimento
Gas
Caliente
FEEDBACK
Ventajas
Desventajas
las
4. Simple de implementar
82
Ingeniera Qumica | Uniatlntico | 2012
FEEDFORWARD
Ventajas
Desventajas
Controlador
Feedback
+
LC
SP
FFC
LT
Controlador
Feedforward
FT
Vapor
Tanque
Hervidor
Agua de
alimento
Gas
Caliente
D(s)
FFC(s)
%TOD
GD(s)
UFF(s)
%CO FF
R(s)
%TO
E(s)
-
%TO
GC(s)
UFB(s)
%CO FB
U(s)
%CO
FV(s)
Lb/h
HD(s)
GM(s)
+
+
C(s)
%TO
Siendo:
Gc(s) = Controlador del lazo feedback
84
Ingeniera Qumica | Uniatlntico | 2012
() = () () + () ()() ()
(. )
() =
()
() ()
(. )
85
Ingeniera Qumica | Uniatlntico | 2012
Debido a que generalmente la deduccin de GD(s) y GM(s), puede ser engorrosa, una
forma sencilla de deducir estos parmetros es por medio del ajuste de FOPDT, con
el mtodo enseado en la leccin 4, y las ecuaciones son las siguientes:
%
+ 1 /
(. )
%
+ 1 %
(. )
() =
() =
%
/
(. )
()
+ ( )
=
(
)
()
()
(. )
De la ecuacin (7.4) es importante destacar que se puede dividir en tres partes: una
parte que contiene solo ganancias y es conocida como controlador feedforward en
estado estacionario (FFCSS), un compensador que posee la forma de una funcin de
transferencia de adelanto/atraso y un compensador de tiempo muerto. El producto de
los tres elementos de la ecuacin (7.4) tambin es conocido como control
feedforward dinmico (FFCDYN), los cuales van a hacer estudiados
individualmente a continuacin. Si se observa primer elemento, correspondientes al
producto y divisin de las ganancias (FFCSS), este posee las siguientes unidades:
%
()
%
/
=
=
() () % %
%
/ %
(. )
86
Ingeniera Qumica | Uniatlntico | 2012
()
=
=
(+)(+)
(. )
87
Ingeniera Qumica | Uniatlntico | 2012
Cuando la relacin
< 1,3, no
Referencias bibliogrficas
77. Dale E. Seborg, Thomas F. Edgar, Duncan A. Mellichamp. Process Dynamics and Control.
Segunda edicin. s.l. : Jhon Wiley & Sons, Inc, 2004. pgs. 389-391.
78. Lyle Albright. Albrights chemical engineering handbook. s.l. : CRC Press, 2009. pgs. 12301233.
79. Carlos A. Smith, Armando B. Corripio. Principles and Practice of Automatic Process Control.
Tercera edicin. s.l. : Jhon Wiley & Sons, Inc, 2006. pgs. 375-389.
80. Steven E. LeBlanc, Donald R. Coughanowr. Process Systems Analysis and Control. Tercera
edicin. s.l. : McGrawHill, 2009. pgs. 361-370.
81. Jonathan Love. Process Automation Handbook: a guide to theory and practice. s.l. : Springer,
2007. pgs. 183-186.
82. B. Wayne Bequette. Process Control: Modeling, Design and Simulation. Upper Saddle River,NJ :
Prentice Hall, 2003. pgs. 324-328.
89
Ingeniera Qumica | Uniatlntico | 2012
90
Ingeniera Qumica | Uniatlntico | 2012
(. )
Siendo:
: Es el volumen del reactor.
: Es el flujo de alimentacin al reactor del componente j.
: Es el flujo de salida del reactor del componente j.
: Es la velocidad de formacin del componente j por unidad de volumen. (El signo
negativo representa que el reactivo se est consumiendo en el interior del reactor).
P(t)
T(t)
XA
P(t)
T(t)
XA
91
Ingeniera Qumica | Uniatlntico | 2012
Ffrio
Reactivos
E-1
Productos
Ffrio
(. )
92
Ingeniera Qumica | Uniatlntico | 2012
Referencias bibliogrficas
41. H. Scott Fogler. Elementos de Ingeniera de las Reacciones Qumicas. Cuarta edicin.
Naucalpan de Juarz : Prentice Hall, 2008. pp. 12-14.
42. . Elementos de Ingeniera de las Reacciones qumicas. Cuarta edicin. Naucalpan de Juarz :
Prentice Hall, 2008. pp. 14-17.
93
Ingeniera Qumica | Uniatlntico | 2012
TT1
TC1
SP1
SP2
TT2
TC2
Falla abierta
Fluido X frio
0-50% CO
Reactor Batch
Falla cerrada
Fluido X recirculado
Fluido X caliente
50-100% CO
(. )
[%] = 0
(. )
(. )
95
Ingeniera Qumica | Uniatlntico | 2012
[%] = 2( 50)
(. )
100
90
80
70
60
50
40
30
20
10
0
-10
10
20
30
40
50
60
70
80
90
100
Temperatura C
Recipiente
Chaqueta
4
2
0
0
10
12
14
16
18
20
tiempo min
50
CO %
-50
10
12
14
16
18
20
tiempo min
vp %
100
0
0
10
12
14
16
18
20
tiempo min
110
100
90
80
70
60
50
40
30
20
10
0
10
20
30
40
50
60
70
80
90
100
Referencias bibliogrficas
46. B. Wayne Bequette. Process Control: Modeling, Design and Simulation. Upper Saddle River,NJ :
Prentice Hall, 2003. pgs. 374-378.
47. Myke King. Process Control, A Practical Approach. The Atrium, Southern Gate, Chichester, UK :
John Wiley & Sons Ltd, 2011. pgs. 140-145.
48. Cecil L. Smith. Practical Process Control, Tuning and Troubleshooting. Hoboken,NJ : John Wiley
& Sons, Inc, 2009. pgs. 387-389.
98
Ingeniera Qumica | Uniatlntico | 2012
FT-A
Estacin de Relacin
R=B/A
SP= B=AxR
FT-B
PV- B
FC
Flujo B
(. )
B
A
(. 2)
Flujo A
FT-A
PV-A
Estacin de Relacin
FT-B
PV-R =B/A
RC
SP= R deseado
PV-B
Flujo B
100
Ingeniera Qumica | Uniatlntico | 2012
= 2=
(. )
(. )
101
Ingeniera Qumica | Uniatlntico | 2012
(. )
Flujo A
PV-A
FT-A
Estacin de Relacin
FC
SP
R=B/A
SP= B=A*R
FT-B
PV- B
FC
Flujo B
Como se observa de la figura, la seal del flujo A, esta es mandada a dos sitios
distintos: uno se dirige como un seal de entrada hacia el controlador que
manipulada la vlvula que regula el flujo A y el segundo camino es hacia la estacin
de relacin donde esta seal se convierte en el SP del controlador que manipulara la
vlvula del flujo B. La ventaja practica de usar esta estrategia de control que usar los
dos lazos de control flujos individuales, es que manteniendo la relacin constante,
solo modificando el SP del controlador que manipula el flujo A este
automticamente ajusta el flujo B, mientras que si tuvieran los lazos
independientemente, habra que modificar dos SP y no uno.
102
Ingeniera Qumica | Uniatlntico | 2012
Referencias bibliogrficas
55. Carlos A. Smith, Armando B. Corripio. Principles and Practice of Automatic Process Control.
Tercera edicin. s.l. : Jhon Wiley & Sons, Inc, 2006. pgs. 333-336.
56. B. Wayne Bequette. Process Control: Modeling, Design and Simulation. Upper Saddle River,NJ :
Prentice Hall, 2003. pgs. 372-373.
57. Dale E. Seborg, Thomas F. Edgar, Duncan A. Mellichamp. Process Dynamics and Control.
Segunda edicin. s.l. : Jhon Wiley & Sons, Inc, 2004. pgs. 391-394.
58. Willian Y. Svrcek, Donald P. Mahoney, Brent R. Young. A Real Time Approach to Process
Control. Segunda edicin. Chichester : John Wiley & Sons. Ltd, 2006. pgs. 138-140.
103
Ingeniera Qumica | Uniatlntico | 2012
Alimento
PV1
SP
TIC
PV
>
PV2
%CO
PV3
TT
TT
TT
Flujo de refrigerante
Reactor
Producto
SP
>
CC
PV1
%CO
PV2
Analizador
Reactivo limite
Reactivo en
exceso
Analizador
Producto
Reactor
Mediana
PV1
FT
PV2
FT
PV
PV3
FIC
FC
SP
%CO
FT
106
Ingeniera Qumica | Uniatlntico | 2012
Otra de los grandes usos de este tipo de control es que mientras los interlocks
deciden tomar una accin sbita de parada de un equipo cuando se alcanza
alguna condicin peligrosa para la seguridad del proceso, los equipos o para el
personal en el rea, esta estrategia de control puede evitar este tipo de situaciones
haciendo que el sistema tome medidas correctivas antes de que se alcance una
condicin accionante de un interlock de seguridad, logrando de esta forma que el
sistema contine operando aunque no sea a su ptimo nivel, es por esto que esta
estrategia es mtodo sencillo de tener un control robusto.
Flujo A
h1
SP
LT
LT
LC
%CO
<
%CO
FC
F
%CO
FT
h2
Al proceso
En conclusin las diferentes variantes del control selectivo son en esencia un lazo de
control feedback sencillo, que posee la adicin de un bloque selector, el cual le
puede seleccionar el valor mximo, mnimo, mediana, promedio etc. de ciertas
mediciones que ingrese al l (estn deben tener las mismas unidades) y dependiendo
de las necesidades del proceso y esta seal puede ser enviada ya sea al controlador o
al elemento de control final dependiendo del caso. Todo esto con la finalidad de
obtener un control efectivo, robusto y sin usar excesivamente el elemento de control
final para el control override.
Referencias bibliogrficas
63. Willian Y. Svrcek, Donald P. Mahoney, Brent R. Young. A Real Time Approach to Process
Control. Segunda edicin. Chichester : John Wiley & Sons. Ltd, 2006. pgs. 140-145.
64. Techmation. Protuner. [En
http://www.protuner.com/Select.PDF.
lnea]
[Citado
el:
20
de
de
2012.]
65. Carlos A. Smith, Armando B. Corripio. Principles and Practice of Automatic Process Control.
Tercera edicin. s.l. : Jhon Wiley & Sons, Inc, 2006. pgs. 340-344.
107
Ingeniera Qumica | Uniatlntico | 2012
MDULO II:
DINMICA DE PROCESOS
Modelamiento y Anlisis de variables
108
Ingeniera Qumica | Uniatlntico | 2012
V(t)
Yi(t)
Fo(t)
Zi(t)
To(t)
P(t)
T(t)
Fw(t)
Tw(t)
h(t)
L(t)
Xi(t)
[ ()]
= () + () () (1 ) ()
( )
Balance de componentes
Un balance de componente i en estado estacionario se plantea con la siguiente
ecuacin:
[ () ()]
= () () + () () () () (1 ) () ()
( )
Relaciones de equilibrio
Una relacin de equilibrio para un componente i es dada por la siguiente ecuacin
() = ()
( )
Balance de energa
Suponiendo que todo el calor cedido por el fluido de calentamiento o enfriamiento
sea igual al calor que entra al separador de fases entonces:
() () = () ()
110
Ingeniera Qumica | Uniatlntico | 2012
[ () ()]
= () () + () () () () (1 ) () ()
( )
= 1
( )
=1
= 1
( )
=1
= 1
( )
=1
Ecuaciones: = + ( ) + + + = +
Variables:
1
111
Ingeniera Qumica | Uniatlntico | 2012
Total (NVariables)
Grados de libertad:
3n + 9
=
= (3 + 9) (2 + 4)
= +
n-1
n+1
=
= = + 5 ( + 1) =
113
Ingeniera Qumica | Uniatlntico | 2012
V(t)
Yi(t)
F3(t)
Z3i(t)
T(t)
F2(t)
Z2i(t)
T(t)
F1(t)
Z1i(t)
T(t)
P(t)
T(t)
H(t)
L(t)
Xi(t)
poder plantear las ecuaciones que rigen este proceso, y seguidamente ser usadas para
la determinacin de los lazos de control que necesita el proceso estudiado. Otra
consideracin e suponer que las densidades (1 , 2 , 3 , , ) y las capacidades
calricas (Cp y Cv) se pueden considerar constante en el tiempo.
Balance global de materia
Un balance global de materia en estado no estacionario se plantea con la siguiente
ecuacin:
[ ()]
= 1 1 () + 2 2 () + 3 3 () () ()
( )
Balance de componentes
Un balance de componente i en estado estacionario se plantea con la siguiente
ecuacin:
[ () ()]
= 1 1 ()1 () + 2 2 ()2 () + 3 3 ()3 () () () () () ( )
Relaciones de equilibrio
Una relacin de equilibrio para un componente i es dada por la siguiente ecuacin
() = ()
( )
Balance de energa
En este caso sera cero ecuaciones porque el proceso de mezclado se realiza a
temperatura constante.
115
Ingeniera Qumica | Uniatlntico | 2012
( )
= 1
=1
= 1
( )
=1
( )
= 1
=1
Las fracciones en los alimentos son tres ecuaciones porque son 3 alimentos que
entran al tanque.
Ecuaciones:
= + ( ) + + = +
Variables:
Nmero de variables
Altura del lquido (h)
3n
Total (N variables)
5n + 7
116
Ingeniera Qumica | Uniatlntico | 2012
Grados de libertad:
=
= (5 + 7) (2 + 5)
= +
Parmetros especificados
Flujo de la corriente de inerte
3(n - 1)
3n - 2
117
Ingeniera Qumica | Uniatlntico | 2012
118
Ingeniera Qumica | Uniatlntico | 2012
QC(t)
PD(t)
Etapa m-1
hD(t)
Lm-1(t)
Xim-1(t)
Pm-1(t)
Tm-1(t)
Etapa m
Vm(t)
Yim(t)
Pm(t)
Tm(t)
D(t)
XDi(t)
R(t)
hm (t)
F(t)
XFi(t)
Q transferido m(t)
Vm+1(t)
Yim+1(t)
Pm+1(t)
Tm+1(t)
Lm(t)
Xim(t)
Pm(t)
Tm(t)
V(t)
Yi(t)
hR(t)
QR(t)
Etapa m+1
B(t)
XBi(t)
a)
b)
Unidades
Balance Global
de
Materia
Balance
de
Componentes
Balance
de
Energa
Equilibrio
Lquido Vapor
ELV
xFi
xi
yi
Platos
Tericos
Np
Np(n - 1)
Np
nNp
Np
Np
Condensador
Total
n-1
Rehervidor
Parcial
n-1
Subtotal
Np + 2
Np(n - 1) + 2(n - 1)
Np + 2
nNp + n
Np + 1
Np + 1
= + 2 + + 2 2 + + 2 + + + 1 + + 1 + + 1
= + + +
Variables
De manera anloga a como se dedujeron las ecuaciones del modelo matemtico de
la columna de destilacin, se organizan las variables de la siguiente forma: variables
correspondientes a los platos tericos, a la parte superior de la Columna y a la parte
inferior de la misma.
Se aclara que los siguientes datos son constantes o dependientes de otras variables y
por tal razn se omiten en el conteo de las variables independientes:
La cada de presin en cada plato, eso es porque estn hechos del mismo
material y las mismas dimensiones (ya sea zona de rectificacin o
agotamiento) y es lgico pensar que la cada de presin en cada plato debera
ser la misma y no vara con respecto al tiempo.
121
Ingeniera Qumica | Uniatlntico | 2012
Nmero de variables
Platos Tericos
Flujo de lquido (Lm) en cada plato
NP
NP
NP
2n NP
1
122
Ingeniera Qumica | Uniatlntico | 2012
Total (Nvariables)
Grados de libertad:
2n NP + 3NP + 4n + 10
=
= (2 + 3 + 4 + 10) (2 + 3 + 3 + 5)
= +
n-1
n -1
123
Ingeniera Qumica | Uniatlntico | 2012
Variables a controlar
=
=
= = + 5 ( 1) =
125
Ingeniera Qumica | Uniatlntico | 2012
126
Ingeniera Qumica | Uniatlntico | 2012
F(t)
Zi(t)
T1(t)
F(t)
Zi(t)
T1(t)
f(t)
Zi(t)
t1(t)
f(t)
Zi(t)
t2(t)
F(t)
Zi(t)
T2(t)
a)
f(t)
Zi(t)
t2(t)
(1-D(t))f(t)
t2'(t)
D(t)f(t)
t1(t)
f(t)
Zi(t)
t1(t)
F(t)
Zi(t)
T2(t)
b)
Son cero ecuaciones debido a la suposicin que la misma cantidad de fluido caliente
o frio que ingresa al intercambiador de calor es igual a su flujo de salida
correspondiente, sin embargo estos flujos son variables en el tiempo, para suplir las
necesidades de recuperacin de calor en el proceso. Este anlisis es igual para la
configuracin con y sin derivacin.
Balance de energa
De forma general, el balance de energa viene dado por la siguiente expresin
() (1 () 2 ()) = () (2 () 1 ())
( )
( )
Ecuacin de diseo
Para la Figura 4.1a la ecuacin viene dada por la siguiente expresin:
(1 () 2 ()) (2 () 1 ())
() (1 () 2 ()) = (
) ( )
1 () 2 ()
ln (
)
2 () 1 ()
(1 () 2 ()) (2 () 1 ())
() (1 () 2 ()) = (
) ( )
1 () 2 ()
ln (
)
2 () 1 ()
()() 1 () + (1 ())() 2 () = () 2 ()
129
Ingeniera Qumica | Uniatlntico | 2012
Simplificando
()1 () + (1 ())2 () = 2 ()
( )
= 1
( )
=1
= 1
( )
=1
Variables
Figura 15.1a
Figura 15.1b
130
Ingeniera Qumica | Uniatlntico | 2012
2n + 6
2n + 8
Total (N variables)
Grado de libertad
=
,
4.1
= (2 + 6) 4 = +
4.1
= (2 + 8) 5 = +
Figura 15.1b
2(n - 1)
2(n - 1)
2n + 1
2n + 1
Variables a controlar
=
131
Ingeniera Qumica | Uniatlntico | 2012
=
.5 ) = 2 + 2 (2 + 1) =
.5 ) = 2 + 3 (2 + 1) =
Fluido caliente
SP
TC
Fluido caliente
Fluido frio
TT
Fluido frio
TT
a)
TC
SP
b)
132
Ingeniera Qumica | Uniatlntico | 2012
Fluido caliente
Fluido frio
TT
TC
SP
Fluido frio
SP
TC
TT
a)
b)
Fluido caliente
Fluido caliente
Fluido frio
TC
Fluido frio
TT
TT
SP
c)
TC
SP
d)
Para escoger la mejor o mejores opciones para la realizacin del lazo de control, de
recirculado en intercambiadores de calor, del punto de vista del diseo y
controlabilidad, se usa las siguientes heursticas (49):
Para un buen control (relaciones ms lineales entre variables de controladas e
manipuladas y actuacin ms rpida del lazo de control), generalmente se
deben derivar la corriente, donde se va controlar la temperatura.
Por cuestiones econmicas de instalacin del sensor/transmisor y vlvula es
recomendable realizar el control en la lnea del fluido frio.
La corriente derivada, tiene que ser de 5 - 10% del flujo total derivado, para
poder manejar las perturbaciones al sistema y ejecutar un buen control.
Siguiendo las heursticas solo las Figuras 15.3a y 15.3c, pueden son buenas y
viables opciones para la aplicacin del lazo de control y tambin la aplicacin del
control de rango dividido, sin embargo, solo la Figura 15.3a cumple con la
133
Ingeniera Qumica | Uniatlntico | 2012
134
Ingeniera Qumica | Uniatlntico | 2012
F1(t)
Z1i(t)
To(t)
V(t)
Yi(t)
F2(t)
Z2i(t)
To(t)
P(t)
T(t)
L(t)
Xi(t)
h(t)
Fc(t)
Tco(t)
Fc(t)
Tc(t)
[ ()]
= 1 1 () + 2 2 () () ()
( )
[ () ()]
= 1 1 () 1 () + 1 1 () 2 () () () () () + () ()
( )
136
Ingeniera Qumica | Uniatlntico | 2012
Relaciones de equilibrio
Una relacin de equilibrio para un componente i es dada por la siguiente ecuacin
() = ()
( )
[ () ()]
= 1 1 ()1 () + 2 2 ()2 () + () ()
Balance de energa
En este caso, el balance de energa queda de la siguiente forma:
[ ()]
= () () + () () (() ())
( )
137
Ingeniera Qumica | Uniatlntico | 2012
( )
= 1
=1
= 1
( )
=1
( )
= 1
=1
Las fracciones en los alimentos son dos ecuaciones porque son 2 alimentos que
entran al tanque.
Variables:
Nmero de variables
Flujos de alimento (F1, F2)
1
138
Ingeniera Qumica | Uniatlntico | 2012
2n
Total (N variables)
4n + 12
Grados de libertad:
=
= (4 + 12) (2 + 6)
= +
Parmetros especificados
Temperaturas de entrada al reactor (To)
2(n - 1)
139
Ingeniera Qumica | Uniatlntico | 2012
2n + 1
Variables a controlar
=
=
= = 2 + 6 (2 + 1) =
140
Ingeniera Qumica | Uniatlntico | 2012
141
Ingeniera Qumica | Uniatlntico | 2012
160
P = 8 bar
B2
D2
(0.6443,138.0027)
140
T C
120
Fase lquida
Fase vapor
B1
100
80
P = 1 bar
60
(0.5,0.5)
0
0.1
0.2
0.3
0.4
0.5
F1
D1
(0.8272,63.3158)
0.6
0.7
0.8
0.9
x1 y1
Figura 17.4. Diagrama Txy para la mezcla THF (1) y agua (2) a presiones
de 1 bar y 8 bar
142
Ingeniera Qumica | Uniatlntico | 2012
P=(P1-P2)
P=(P2-P1)
F1
D1
0,8272 THF
0,1728 Agua
P= 1bar
P= 8bar
D2
0,6443 THF
0,3557 Agua
0,5 THF
0,5 Agua
B1
Agua pura
B2
THF puro
Estos son los pasos a seguir para poder deducir la Figura 17.5 a partir de la Figura
17.4
Si se posee una mezcla equimolar y se destila a una presin de 1 bar la
mxima composicin del destilado ser la del azetropo (D1), y la de fondo es
agua (B1).
Es necesario recordar que en la destilacin, salen los vapores que posean menor
punto de ebullicin de la mezcla, y por el fondo los productos que posean mayor
punto de ebullicin.
El destilado obtenido en primera columna se transporta hacia la columna por
una bomba, la cual tiene la finalidad de aumentar la presin de la corriente D1
hasta la del plato que entra en la segunda columna de destilacin.
De manera anloga a la primera columna, si a la mezcla de la corriente D 1 se
alimenta a una columna de destilacin, la mxima composicin del destilado
seria la azeotropa (D2) y como el producto de fondo se obtiene THF puro
(B2).
Debido a que no se puede desperdiciar THF y agua el destilado azeotrpico
de la columna 2 (D2), primero se le disminuye la presin igual a F, se
recircula y se mezcla con la corriente de alimento a la primera columna
formando la corriente F1. Otra opcin hubiese sido que la corriente D2 se
recircula y se alimenta a un plato de la primera columna, por lo cual se le
disminuye la presin hasta la del plato correspondiente.
143
Ingeniera Qumica | Uniatlntico | 2012
flujo de calor que sale condensador puesto que afecta directamente esta
variable (directriz 8).
La composicin en el fondo de la columna se debe controlar porque esta
afecta la calidad de los productos (directriz 3), pero, como este tipo de
controlador por su misma dinmica, es mejor escoger un control, ms rpido,
y ms sensible, que de manera indirecta controle la composicin (directrices
11 y 12), la cual sera la temperatura en el plato de la columna donde la
variacin de composicin del componente clave ligero sea mayor. Esta ser
manipulada con el flujo calrico del rehervidor, porque esta variable afecta
directamente la variable controlada (directriz 8).
La composicin de destilado se controlara porque afecta la calidad del
producto (directriz 3), pero como este tipo de controlador por su misma
dinmica, es mejor escoger un control, ms rpido, y ms sensible, que de
manera indirecta controle la composicin, por eso se usara controlar el flujo
de reflujo (directrices 11 y 12), y este ser manipulado con el actuador de la
vlvula que manipula el flujo, porque esta afecta esta variable de entrada
(directriz 8).
El nivel de lquido en el tambor de reflujo, se controlara, porque su dinmica,
no es auto reguladora (directriz 1) esta variable ha de ser controlada con el
flujo de destilado, puesto a que esta afecta directamente esta variable de
salida (directriz 8). A este lazo de control se le acopla el control por relacin
por multiplicacin.
Columna de alta presin (segunda columna) configuracin LV
Flujo de alimento a la columna (es el mismo flujo de destilado) no se
realizara, porque este se hace indirectamente en el control de tambor de
reflujo de la primera columna.
El flujo de reflujo ser controlado porque afecta directamente la calidad del
destilado y es una variable ms sensible y presenta una respuesta ms rpida
que la medicin y control que la composicin del destilado (directrices 3, 11
y 12). Esta variable ser manipulado por medio del actuador vlvula que
regula el mismo flujo (directriz 8).
El nivel de lquido en la base de la columna ser controlado por su dinmica
no auto reguladora (directriz 1) y ser manipulado por el flujo de fondo que
sale de la segunda columna (directriz 8).
145
Ingeniera Qumica | Uniatlntico | 2012
146
Ingeniera Qumica | Uniatlntico | 2012
MDULO III
SIMULACIN DINMICA DE PROCESOS
Diseo y Control asistido con Aspen-HYSYS
147
Ingeniera Qumica | Uniatlntico | 2012
Paquete Fluido
Componentes
Ecuacin
Peng Robinson
Reacciones
Sistema de unidades
Euro Si
149
Ingeniera Qumica | Uniatlntico | 2012
Aqu aparecer una ventana, en la cual hay que buscar un panel de nombre
Variables, en este aparecer un push button de nombre Add, el cual se
le hace clic (Figura 18.3).
= 2492
(. )
3
1
3
= 0,0758
547,6 60
151
Ingeniera Qumica | Uniatlntico | 2012
(. )
0,61
11.08
= 0.1832
13
1
3
= 0,06919
11,08 3600
(. )
0,06919
=
=
0,1832
= ,
152
Ingeniera Qumica | Uniatlntico | 2012
Esta rea seccional debe ser mnima, puesto que fue calculada con velocidad
mxima de vapor permisible, entonces el dimetro que equivale a esta rea
viene da dado por la siguiente expresin:
2
=
4
(. )
4
=
= 0,6935
=
3
2 2 (2) 3
=
=
4
4
2
2
(. )
= .
= 2 = ,
Deben introducirse 2 datos (altura-dimetro, altura-volumen o dimetrovolumen) y le calcula la especificacin que falta, tambin solo puede ingresar el
volumen de separador de fases y este le calcula el dimetro y la altura con una
relacin H/D =1.5.
Observe que hay tres opciones para el mtodo de diseo: Cv para lquidos, Cg
para gases (el simulador dependiendo de las condiciones de flujo, escoge Cv o
Cg) y el mtodo k que se basa en ecuaciones de resistencia del tipo k.
Tambin se puede escoger entre tres tipos de caractersticas de la vlvula:
Linear, Quick Opening (abertura rpida) y Equal Percentage (igual
porcentaje), recuerde que tenga en cuenta que el porcentaje de abertura
realmente debe entenderse como la posicin del vstago de la vlvula que
154
Ingeniera Qumica | Uniatlntico | 2012
que colocar la accin del controlador inversa, mientras que si lo colocas como
falla cerrada, no sufre ese problema. Para no confundir al lector con las acciones
del controlador, se deja todo, como vienen por defecto.
Siguiendo dimensionando las otras 2 vlvulas, la vlvula V3, se puede
dimensionar con cualquiera de los 2 mtodos aqu mencionados, pero la vlvula
V2, se dimensiona, con el primer mtodo, porque este elemento de control final,
va a ser usado en control de presin, y la caracterstica con la que debe trabajar
debe ser de igual porcentaje y no lineal, que es la que viene por defecto.
a)
b)
Figura 18.8. Asistente Dinmico a) cuando se abre por primera vez para hacer
el cambio de modo y b) despus que todos los cambios se han hecho de forma
satisfactoria
156
Ingeniera Qumica | Uniatlntico | 2012
157
Ingeniera Qumica | Uniatlntico | 2012
158
Ingeniera Qumica | Uniatlntico | 2012
= (
(. )
Hasta este punto el controlador est listo para ponerlo a funcionar como se muestra
en la Figura 18.12. Si se aprecia bien el SP es la mitad del rango de la variable a
controlar (PV), lo ltimo que hace falta es comprobar si el controlador realmente
sirve, entonces:
a)
b)
Figura 18.13. Face Plate del controlador de flujo a) en modo manual sin
activar el integrador y b) en modo automtico, activado el integrador
164
Ingeniera Qumica | Uniatlntico | 2012
168
Ingeniera Qumica | Uniatlntico | 2012
Cabe notar que como el valor de la variable a controlar o PV es igual al valor del SP
en cualquier intervalo de tiempo, la lnea correspondiente al PV queda sobrepuesta a
la del SP y no ve; otra anotacin es cuando hace clic a alguna de las lneas, el eje de
las ordenadas cambia de las unidades en las que se encuentra, a la de la lnea
169
Ingeniera Qumica | Uniatlntico | 2012
Referencias bibliogrficas
17. . Plantwide Dynamic Simulators in Chemical Processing and Control. New York : Marcel
Dekker, Inc, 2002. pp. 9-15.
170
Ingeniera Qumica | Uniatlntico | 2012
En esta prctica se har uso de unas de las herramientas que posee Aspen-HYSYS
con respecto a la simulacin dinmica, que es la estimacin de los parmetros
dinmicos, por medio del mtodo de sintonizacin ATV. Por esta razn, se simulara
un proceso de mezclado (26), en cual se agrega etanol (F1), con un alto grado de
pureza (95% en fraccin msica), que ha de ser diluido con agua (F2) en un tanque
de 1000 gal U.S hasta una composicin de ms o menos 27% en fraccin msica.
Tambin se agrega un flujo pequeo y constante de aire comprimido (F3), para
regular la presin del tanque y una bomba para poder transportar la mezcla de etanol
diluida.
El sistema, adems del tanque y la bomba, cuenta con 5 vlvulas, 3 en las corrientes
de alimentacin, una de la corriente de venteo, y en la lnea de descarga del tanque
(esta siempre va despus de la bomba para evitar la cavitacin).
Ecuaciones
Reacciones
Sistema de unidades
Field
Etanol: Instale una corriente de materia de nombre F1_etanol, que posee 100
Lbmol/h de 95% de etanol y 5% de agua en fraccin msica, una presin de 150 psia
y una temperatura de 90F. Despus instale una vlvula de nombre V1, la cual posee
una corriente de entrada F1_etanol y de salida Sale_V1. Para que la vlvula quede
totalmente especificada asigne una cada de presin de 40 psia.
Agua: Instale una corriente de materia de nombre F2_agua, que posee 600 Lbmol/h
de 100% de agua, una presin de 150 psia y una temperatura de 90F. Despus
instale una vlvula de nombre V2, la cual posee una corriente de entrada F2_agua y
de salida Sale_V2. Para que la vlvula quede totalmente especificada asgnele una
cada de presin de 40 psia.
Aire: Instale una corriente de materia de nombre F3_Aire, el cual posee 7 Lbmol/h
de 79% de nitrgeno y 21% de oxgeno en fraccin molar, una presin de 150 psia y
una temperatura de 90F. Despus instale una vlvula de nombre V3, la cual posee
una corriente de entrada F3_Aire y de salida Sale_V3. Para que la vlvula quede
totalmente especificada asigne una cada de presin de 40 psia.
Tanque de mezclado: Instale un tanque de la paleta de objetos
(tambin se
puede de esta forma: instalar un separador de fases, en la pestaa Design, en la
opcin Parameters, aparecern tres radio botones y escoge el de nombre Tank)
172
Ingeniera Qumica | Uniatlntico | 2012
2 2 (2) 3
=
=
4
4
2
3 2
= ,
= 2 = ,
El volumen y dimetro (o la altura) se insertan en la pestaa de Dynamics
del tanque de mezclado.
Diseo de las vlvulas de control
Las cinco vlvulas de control estn dimensionadas, con un porcentaje de
abertura del 50%. Las vlvulas V1 y V4, van a hacer diseadas, con
caractersticas de igual porcentaje, debido que las PV que regulan estos FCE,
no poseen comportamiento lineal. Mientras que las otras tres se dimensionan
con caracterstica lineal.
Para dimensionar una vlvula con caracterstica lineal, se hace en pestaa
Dynamics de la ventana de la vlvula, y se hace clic en Size Valve y
hacer check en prevenir el retroflujo Check Valve (Prevents Backflow).
Para dimensionar una vlvula con caracterstica de igual porcentaje se hace
en la pestaa Rating de la ventana de la vlvula, se escoge el radio botn
Equal Percentage y se hace clic en Size Valve.
174
Ingeniera Qumica | Uniatlntico | 2012
176
Ingeniera Qumica | Uniatlntico | 2012
179
Ingeniera Qumica | Uniatlntico | 2012
a)
b)
a)
b)
181
Ingeniera Qumica | Uniatlntico | 2012
100
0,2797 0,269
100 = 7.64%
0,339 0,199
4(5%)
= .
(7,64%)
182
Ingeniera Qumica | Uniatlntico | 2012
= 4 3
= (155,99 145,48) = .
Para un controlador PID paralelo, con igual a cero, usando las heursticas de ZN y
TL, los parmetros dinmicos son los siguientes:
Kp
min
min
PID segn ZN
0,6225
6,562
1,05
PID segn TL
0,377
23,1
1,667
Como se muestra en la Tabla 2.3 los valores obtenidos por el criterio de la ganancia
ltima y el periodo ltimo son diferentes a los calculados por Aspen-HYSYS
(Figura 19.8b)), sin embargo los parmetros calculados por el simulador poseen ms
similitud con los parmetros dinmicos de TL.
Por ltimo active el Face Plate de este controlador, active el integrador,
dejar pasar unos segundos, pasar el controlador de manual a automtico, dejar
pasar un par de segundos ms y apague el integrador
A partir de este momento active todas los Face Plates y active el integrador y la
simulacin ha debido converger satisfactoriamente, como se muestra en la Figura
19.10.
183
Ingeniera Qumica | Uniatlntico | 2012
Referencias bibliogrficas
26. Willian L. Luyben. Plantwide Dynamic Simulators in Chemical Processing and Control. New
York : Marcel Dekker, Inc, 2002. pp. 87-94.
27. AspenTech. Aspen HYSYS, operations guide. Burlington, Ma : Aspen Technology, Inc, 2009. pp.
5-179 - 5-188.
184
Ingeniera Qumica | Uniatlntico | 2012
Ecuacin
Peng Robinson
Reacciones
Sistema de unidades
Field
= +
= + (
(. )
(. )
5
1
33,48
1 2
1
= 201,7 + 30 (
)(
)(
)
)
(
)(
3
2
1
12
144
2
= ,
Presin
Cada de Presin
Condensador total
Rehervidor
201,7 psia
204,6 psia
0 psia
0 psia
186
Ingeniera Qumica | Uniatlntico | 2012
La Figura 20.1 muestra la ventana final del asistente con las corrientes conectadas a
la columna, el nmero de platos, el plato de alimentacin, las presiones y cadas de
presin en condensador y rehervidor.
Para que el sistema quede especificado por completo, hacen falta especificar 2
parmetros: las composiciones de cabeza y fondo de propano. Esto se asigna en la
opcin Spec de la pestaa Design. Primero se borra las especificaciones que
viene por defecto, excepto la relacin de reflujo Reflux Ratio. Seguidamente se
agrega 2 veces Column Component Fraction y se especifican como se muestra en
las Figuras 20.2a y 20.2b.
a)
b)
187
Ingeniera Qumica | Uniatlntico | 2012
De la figura se puede deducir que la mejor configuracin del control es LV, porque
el valor del reflujo externo es igual a 3.58.
Despus instale una bomba de nombre B1, el cual tiene como corriente de succin
Destilado y de descarga Sale_B1. Para especificar la bomba totalmente asigne un
aumento de presin de 30 psi. Seguidamente instale una vlvula de nombre V2, que
tiene como corriente de entrada Sale_B1 y como corriente de salida Sale_V2,
asigne una cada de presin de 20 psia para que quede especificado.
De forma anloga, instale una bomba de nombre B2, el cual tiene como corriente
de succin Fondo y de descarga Sale_B2, para especificar la bomba, asigne un
aumento de presin de 30 psia. Despus instale una vlvula de nombre V3, que
tiene como corriente de entrada Sale_B2 y como corriente de salida Sale_V2, asigne
una cada de presin de 22,9 psia para que la simulacin converja
satisfactoriamente como se muestra en la Figura 20.4. Para que en el PFD aparezca
el Workbook en forma de tabla se debe hacer clic derecho en el PFD y se escoge
la opcin Add WorkBook Table.
Por ltimo, se define en qu plato se va a instalar el control de temperatura, en el
cual la variacin de composicin del componente ms liviano (propano) sea mayor.
Para hacer esto, despliegue la ventana de propiedades de la columna, seleccione en
la pestaa Performance, en la opcin Plots. En esta ventana se selecciona
Composition dentro del cuadro de ttulo Tray by Tray Properties. Por ltimo
188
Ingeniera Qumica | Uniatlntico | 2012
haga clic en View Table y se mostrar la tabla que muestra variacin de los tres
componentes en toda la columna (Figura 20.4). En la figura se aprecia que los
mayores cambios de composicin de propano ocurren entre los platos 5 y 6 (0,76950,694 =0,0755) y entre los platos 6 y 7 (0,694-0,6174=0,0766). Puesto que el plato
en comn de estos 2 rangos es el plato 6, este sera el plato en el cual habra que
controlar la temperatura (cuya temperatura es 123,8 F) usando el flujo de calor del
rehervidor. Hasta este punto guarde esta prctica de nombre COLUMNA DE
DESTILACIN ESTADO ESTACIONARIO.
189
Ingeniera Qumica | Uniatlntico | 2012
( 1) =
( 30) =
1
1,926
=
= 0,7205
1
2,305
= ,
Los clculos demuestran que la columna debe ser diseada con respeto al
plato 30, esto tambin se puede corroborar calculando el rea mnima para
evitar el arrastre usando la ecuacin de continuidad.
( 1):
(. )
1
)
(
3600
=
= 1,2154 2
(0,7205 ) ( 1,926 3 )
( 30): =
6072
1
(
)
3600
(0,6586 )( 2,305 3 )
9693
= ,
@ 30 4
1,7736 2 4
=
= ,
191
Ingeniera Qumica | Uniatlntico | 2012
=
= 1
: 6072
3
1
(
)(
) 10min = ,
29,17
60
: 13760
3
1
(
)(
) 10min = ,
29,33
60
Otra forma de insertar los valores para el diseo de los platos, el tambor de reflujo y
acumulador de la base, sera en el ambiente normal de simulacin, en la ventana de
la Columna Despropanizadora, en la pestaa Rating en la opcin Tray Sections,
se puede agregar las dimensiones y el tipo de plato, y en la opcin Vessels, se
pueden insertar las dimensiones del tanque de reflujo y acumulador en la base.
a)
b)
196
Ingeniera Qumica | Uniatlntico | 2012
198
Ingeniera Qumica | Uniatlntico | 2012
de la paleta de objetos.
Asigne como PV el nivel del lquido del rehervidor, Vessel Liq Percent
Level (en la ventana donde se selecciona el PV se escoge como Flowsheet
Columna Despropanizadora, en "Object Reboiler y en Variable Vessel
Liq Percent Level), y como valor de OP no se selecciona nada, porque en
Aspen-HYSYS solo se puede conectar la seal OP de la funcin de
transferencia, desde el PV del controlador y no al revs.
En la pestaa Parameters en la opcin Configuration, se insertan los
rangos de la variable a controlar (PV). Debido a que la finalidad de la
funcin de transferencia, es atrasar o alentar el valor de PV, es lgico pensar
que los rangos de PV y OP son iguales. Como el valor en estado estacionario
es un nivel de lquido de 50%, se supone una tolerancia ms o menos flexible
de 25 %, es decir que los lmites de PV a insertar serian 25% y 75% (Figura
20.13a).
En la misma pestaa, pero en la opcin Delay, se hace check Delay, para
que la funcin de transferencia posea un atraso del tipo de tiempo muerto. En
Delay parameters, se le asigna a este atraso una ganancia de 1 y un tiempo
muerto de 30 s, como se muestra en la Figura 20.13b.
a)
b)
201
Ingeniera Qumica | Uniatlntico | 2012
de la paleta de objetos
203
Ingeniera Qumica | Uniatlntico | 2012
Referencias bibliogrficas
33. Warren D. Seider, J.D. Seader, Daniel R.Lewin, Soemantri Widagdo. Product and Process
Design Principles: Synthesis, Analysis and Evaluation. Tercera edicin. Westford : John Wiley &
Sons. Inc, 2009. pp. 329-331.
34. Willian L. Luyben. Plantwide Dynamic Simulators in Chemical Processing and Control. New
York : Marcel Dekker, Inc, 2002. pp. 155-172.
35. Donald Q. Kern. Procesos de Transferencia de Calor. Colonia San Juan Tlihuaca, Delegacin
Azcapotzalco : Cecsa, 1965. p. 118.
36. Willian L. Luyben. Plantwide Dynamic Simulators in Chemical Processing and Control. New
York : Marcel Dekker, Inc, 2002. pp. 11-15.
204
Ingeniera Qumica | Uniatlntico | 2012
Paquete Fluido
Componentes
Ecuacin
Reacciones
Sistema de unidades
cual posee de nombre Q_separador y con corrientes de salida Vapor y Lquido, para
las corrientes de vapor y liquido respectivamente; para que el separador de fases
converja, asigne una temperatura a la corriente Vapor de 30C. Despus instale la
vlvula V2, la cual tiene como corriente de entrada, Vapor, y como corriente de
salida Sale_V2, a esta vlvula se asigna una cada de presin de 2 bar; Instale una
bomba de nombre B1, el cual posee como corriente de entrada la corriente Lquido,
de salida Sale_B1 y corriente de energa W_B1, para que esta bomba converja,
asigne aumento de presin de 2 bar. Por ltimo instale una vlvula de nombre V3, la
cual posee como entrada Sale_B1, y de salida, Sale_V3, a esta vlvula se le asigna
una cada de presin de 4 bar, despus de hacer esto, el sistema en estado
estacionario ha debido converger satisfactoriamente, como se muestra en la Figura
21.1.
Usando los datos obtenidos en la seccin 18.3, se estimaron las dimensiones del
separador de fases, las cuales fueron Volumen = 1,516 m3 y Dimetro = 0.9882 m,
datos que se insertan en la pestaa Dynamics del separador de fases.
206
Ingeniera Qumica | Uniatlntico | 2012
Con respecto a las vlvulas todas se disean con un porcentaje de abertura de 50%,
con caracterstica lineal las vlvulas V1 y V3, mientras que la vlvula V2 se disea
con caracterstica de igual porcentaje.
Como el flujo de calor que ingresa al separador de fases es 8,4537*104 Kcal/h, se
usa una tolerancia de 4*104 Kcal/h, entonces los rangos para dimensionar la
vlvula que regula este flujo de calor son 4,4537*104 y 1,24537*105 Kcal/h.
En este momento guarde esta simulacin con el nombre SEPARADOR DE
FASES LAZO ABIERTO.
209
Ingeniera Qumica | Uniatlntico | 2012
Abra la ventana del integrador Ctrl + I y asigne que este se detenga End
Time en 50 min.
Active el integrador y vera que este se detendr automticamente en 50 min.
El lector debe notar que el registrador debi aparecerle, despus de hacerle unos
cambios como cambiar los lmites de los ejes, color, grosor de las lneas (esto se
hace en clic derecho en Graph Control), como la Figura 21.2.
1 = 0,283(0,2714 ) = 0,0768
2 = 0,632(0,2714 ) = 0,1715
1 = (30 + 0,0768) = ,
2 = (30 + 0,1715) = ,
1 @ 1 1295 = ,
2 @ 2 1415 = ,
Antes de calcular los parmetros del FOPDT, es importante destacar que la ganancia
K, el valor del denominador esta en porcentaje y por consiguiente el
numerador , tambin debe estar en porcentaje, y que los tiempos 1 y 2 , tienen
que ser restados por el tiempo inicial (20 min), haciendo alusin de las variables
desviacin, concepto fundamental para el desarrollo de este mtodo.
Entonces los parmetros de ajuste del FOPDT seran igual a:
0,2714
100%
=
= ,
(60 50)% (40 20)
derivativo para un controlador PID en paralelo con un igual a cero, los cuales son
presentados en la Tabla 21.1.
Ganancia
Proporcional
Kp
Tiempo Integral
min
Tiempo Derivativo
min
ZN
PID paralelo
56,8806
1,4575
0,2332
1,32849
0,20476
CC
PID paralelo
1,97902
Con respecto al IMC, se supone un valor de lambda sea 50% ms grande que el
valor mnimo o:
= 1,5(0.8 ) = ,
[3 + 0,5(0,583)]
+ 0,5
=
= ,
(0,5 + )
0,1357[( 0,5(0,583) + 0,6996 )]
= + 0,5 = [3 + 0,5(0,583)] = ,
213
Ingeniera Qumica | Uniatlntico | 2012
(3)(0,583)2
=
=
= ,
2 + [2(3) + 0,583]
Este mismo proceso se repite, pero usando los parmetros dinmicos de CC, y el
registrador le presentara el comportamiento de PV, OP y SP del controlador de
temperatura, como se muestra en la Figura 21.5.
214
Ingeniera Qumica | Uniatlntico | 2012
Con respecto al mtodo IMC, Aspen-HYSYS posee una forma de calcular los
parmetros dinmicos, solo ingresando, los valores del ajuste de FOPDT de la curva
de reaccin. Por tal motivo, se abre la ventana del control de temperatura y en la
pestaa Parameters, en la opcin IMC Design, se insertan los valores de , ,
y , como se muestra en la Figura 21.6.
En esencia los tres mtodos de sintonizacin son efectivos, puesto que llevan la
variable controlada, temperatura del separador de fases, hacia su nuevo estado
estacionario de 31 C, sin embargo, en las figuras I.12 y I.15 se puede apreciar que
el porcentaje de abertura llega a 100% por un breve periodo de tiempo, pero el valor
de la variable controlada segua aumentando, a este fenmeno se le conoce como
Windup y debe ser evitado, en especial en procesos donde esta situacin es crtica y
no deseada, como es el caso del control de temperatura de un reactor CSTR,
entonces en este caso las heursticas de CC son ms convenientes que los otros 2
mtodos. Una comparacin de estos mtodos con el ATV (usando los parmetros
que Aspen-HYSYS trae por defecto), este genera una ganancia de 14, un tiempo
integral de 1,26 min y un tiempo derivativo de 0,28 min, valores que son muy
parecidos al mtodo de sintonizacin en lazo abierto.
Otra consideracin importante que a pesar que el mtodo IMC produjo Windup, este
problema se puede solucionar aumentado el parmetro , el lector puede intentar
aumentar este parmetro a su gusto, sin embargo se le recomienda el uso de =10, y
podr apreciar en el registrador (Figura 21.8), tendr una figura parecida que usando
las heursticas de CC, pero con la ventaja que est el OP sube progresivamente,
semejante a una respuesta monotnica estable hasta llegar a su OP de 88,31 %.
216
Ingeniera Qumica | Uniatlntico | 2012
Referencias bibliogrficas
76. AspenTech. Aspen HYSYS, operations guide. Burlington, Ma : Aspen Technology, Inc, 2009.
pgs. 5-121 - 5-122 .
218
Ingeniera Qumica | Uniatlntico | 2012
Ecuacin
Peng Robinson
Reacciones
Sistema de unidades
Euro Si
Usando los datos obtenidos en la seccin 3.2, se estimaron las dimensiones del
separador de fases, las cuales fueron Volumen=1,516 m3 y Dimetro =0,9882 m,
datos que se insertan en la pestaa Dynamics del separador de fases.
220
Ingeniera Qumica | Uniatlntico | 2012
Con respecto a las vlvulas todas se disean con un porcentaje de abertura de 50%,
con caracterstica lineal las vlvulas V1 y V3, mientras que la vlvula V2 se disea
con caracterstica de igual porcentaje.
En este momento guarde esta simulacin con el nombre SEPARADOR DE
FASES EN CASCADA- ESTADO ESTACIONARIO y cree 2 copias de
nombre SEPARADOR DE FASES 2 CONTROLES EN CASCADA- MODO
DINMICO y SEPARADOR DE FASES 3 CONTROLES EN CASCADAMODO DINMICO.
control de flujo de
221
Ingeniera Qumica | Uniatlntico | 2012
222
Ingeniera Qumica | Uniatlntico | 2012
a)
b)
a)
b)
226
Ingeniera Qumica | Uniatlntico | 2012
Figura 22.4. Diagrama de flujo del Separador de fases con cascadas de segundo
nivel en modo dinmico
230
Ingeniera Qumica | Uniatlntico | 2012
231
Ingeniera Qumica | Uniatlntico | 2012
De la figura se puede apreciar que dice que el SP tiene las opciones de Track PV,
Track Remote y Use %. La primera opcin hace que referencia que si el SP
fuese tomado dentro del controlador, como un lazo sencillo de feedback, se mueva
con el valor de PV. La segunda opcin es igual que la primera, slo que es para el
SP remoto (control en cascada) y la tercera opcin es que SP remoto no es recibido
como unidades sino de forma porcentual, lo cual permite que no exista problemas de
unidades entre la seal de salida del controlador maestro y el SP del controlador
esclavo, la otra opcin de Use PV Units, tiene mas significado cuando los
controladores son del mismo tipo como flujo-flujo, nivel-nivel, situacin que es ms
comn en la estrategia de control por relacin (mdulo G).
232
Ingeniera Qumica | Uniatlntico | 2012
Agua
Lquido: UNIQUAC - Vapor: Ideal
No hay. No es un proceso con reaccin
Field
calor Exchanger Design (End point), debido a que en este sistema no hay cambios
de fase y el Cp es relativamente constante y como tipo de coraza segn Tubular
Exchanger Manufacturers Association (TEMA) de tipo E, lo quiere decir que la
coraza es de un solo paso.
Para las cadas de presin en el lado de los tubos y la coraza, se asignaran 7,5psia y
7,7 psia, respectivamente. Asumiendo coeficiente global de transferencia de calor
(U) de 259 BTU/ft2 F h y un rea de transferencia (A) de 502 ft2, asignar como UA,
el producto entre estos valores, 130018 BTU/F h y observara que la simulacin ha
convergido satisfactoriamente como se muestra en la Figura 23.1.
Para cambiar el icono del intercambiador al mostrado en la figura, se hace clic
derecho sobre icono del intercambiador y escoge la opcin Change Icon y el
simulador le presentara las diferentes opciones que tiene el simulador para cambiar
el icono, para este caso se escogi la segunda.
En este momento guarde la simulacin con el nombre de INTERCAMBIADOR
DE CALOR AGUA-AGUA.
23.3
234
Ingeniera Qumica | Uniatlntico | 2012
Tubos
Coraza
Dimetro = 15,25 in
OD = 0,75 in
Pitch = 0,9375 in
Di = 0,652 in
Arreglo Triangular
L =16 in
Coeficiente de ensuciamiento =
0,0005 ft2 F h /BTU
Coeficiente de ensuciamiento =
0,0015 ft2 F h /BTU
Espaciado de bafles = 12 in
23.4
la ventana del asistente dinmico, en la cual, primero debe hacer check en Save
Steady State Case. Despus de eso, hacer clic en Make Changes y los cambios se
habrn cumplido satisfactoriamente. Por ltimo se hace clic en Dynamic Mode de
nuevo, y aparecer un cuadro de dialogo para confirmar el cambio de modo, el cual
se da clic en s.
23.5
primero se instala el control feedforward, por tal motivo en este punto haga 2 copias
de este archivo: uno con el nombre INTERCAMBIADOR DE CALOR AGUAAGUA FOPDT TEMPERATURA en el cual se estimaran los parmetros del
FOPDT para el control de temperatura en cascada con un controlador de flujo y el
otro archivo INTERCAMBIADOR DE CALOR AGUA-AGUA FOPDT
PERTURBACIN donde se estimaran los parmetros de FOPDT de la
perturbacin del flujo de Agua Destilada y como esta afecta la temperatura de la
corriente Sale_coraza. .
En el primer archivo realice los siguientes pasos:
En la pestaa Stripchart del controlador de temperatura, cree un registrador
que solo muestre el SP, OP y PV de este lazo de control.
Abra la ventana del integrador Ctrl + I y reinicie el integrador y adems
asigne que este se detenga End Time en 50 min.
Para que el registrador muestre una curva ms suave, se debe abrir el
Databook Ctrl +D y en la pestaa Stripchart se hace clic en Setup All y el
Logger Size se cambia el valor de 300 a 3000 y en Sample Time cambie
de 20 s a 5 s.
Active el integrador y vera que este se detendr automticamente en 50 min.
Abra la ventana del controlador de temperatura en la pestaa Parameters en
la opcin Configuration, cambie el OP de 50% a 70%.
Abra la ventana del integrador Ctrl + I y asigne que este se detenga End
Time en 100 min.
Active el integrador y vera que este se detendr automticamente en 100 min.
1 = 0,283(0,087 ) = 0,0246
2 = 0,632(0,087 ) = 0,0549
1 = (85.058 0,0246) = 85,0334
2 = (85.058 0,0549) = 85,0031
1 @ 1 3035 = 50,583
2 @ 2 3080 = 51.333
0,087
100%
%
=
= ,
(70 50)% (90 80)
%
, , %
() =
, + %
239
Ingeniera Qumica | Uniatlntico | 2012
Los pasos para obtener la curva de reaccin de esta perturbacin con respecto a la
variable controlada son los siguientes (este controlador ficticio debe estar en modo
manual).
240
Ingeniera Qumica | Uniatlntico | 2012
241
Ingeniera Qumica | Uniatlntico | 2012
1 = 0,283(0.6827F ) = 0,1932
2 = 0,632(0.6827F ) = 0,4314
1 = (85.058 + 0,1932) = 85,2512
2 = (85.058 + 0,4314) = 85,4894
1 @ 1 3025 = 50,4167
2 @ 2 3060 = 51
Entonces los parmetros de ajuste del FOPDT seran igual a (se recuerda que
es la diferencia del flujo cuando el OP es 70% y cuando era 50%)
0.6827F
100%
%
=
,
() =
, , %
, +
100 %
%
5105
()
+ 1 ( )
(
)
=
() ()
+ 1
242
Ingeniera Qumica | Uniatlntico | 2012
, +
%
) ,
= , (
, +
%
(. )
funcione debe hacer check en Enable FeedForward, y que la seal de salida del
controlador feedforward es 0%, la razn es que como la variable controlada no est
siendo afectada por la perturbacin, entonces este controlador no enviara alguna
seal correctiva.
Seguidamente en la misma pestaa, pero en la opcin Autotuner, escoger la
opcin de sintonizar un controlador PID, hacer check en Automatically Accept,
abrir la caratula Face Plate, hacer clic en Start Autotuner, encender el
integrador, espere a que le aparezcan los parmetros dinmicos en la ventana del
controlador y apague el integrador.
El lector se ha percatado el OP de este controlador no quedo en aproximadamente
50% cuando se usa el mtodo ATV, sino en que el OP queda aproximadamente en
28% y si se dejara activado el integrador este comenzara a oscilar entre 28% y
aproximadamente 73%. Esta es la situacin se debe a la mala conmutacin de las
seales del control feedback y feedforward. Para solucionar existen dos formas, la
primera es simplemente pasando el controlador a modo automtico y despus se
pasa a automtico el controlador feedforward, o active el control feedforward
despus que haya sintonizado. Para este caso se usara la primera opcin, aunque el
lector puede intentar la segunda, pero con el problema es que si usted retira el check
en Enable FeedForward se borrara la informacin que haya agregado.
Una vez activado el integrador y pasado el controlador PID y el feedforward la
simulacin ha debido convergerle satisfactoriamente como se muestra en la Figura
23.6.
244
Ingeniera Qumica | Uniatlntico | 2012
lazo
usando
el
control
Referencias bibliogrficas
83. Donald Q. Kern. Procesos de Transferencia de Calor. Colonia San Juan Tlihuaca, Delegacin
Azcapotzalco : Cecsa, 1965. pgs. 191-194.
84. AspenTech. Aspen HYSYS, operations guide. Burlington, Ma : Aspen Technology, Inc, 2009.
pgs. 5-130 - 5-131.
247
Ingeniera Qumica | Uniatlntico | 2012
6 7 () + 32() 6 13 ()
(. )
O de forma ms general
(. )
[ 3 ]
(. )
Por otro lado la constante de velocidad de reaccin especfica viene dada por la
ecuacin de Arrhenius.
20000 3
) [ ]
= 1108 exp (
(. )
248
Ingeniera Qumica | Uniatlntico | 2012
El reactor posee un volumen de 1200 ft3, con nivel de lquido del 80%. Este opera a
una presin de 600 psia, una temperatura de 250 F y una conversin de la anilina en
aproximadamente del 74%. Debido a que la reaccin es exotrmica (entalpia de
formacin estndar -8.2x104 BTU/Lbmol), se utiliza un sistema de enfriamiento.
Adems del reactor, el sistema cuenta con 5 vlvulas, 4 instalada manualmente y
una ltima, como un algoritmo interno en la corriente de energa. Es importante
anotar que en este caso la corriente de energa ser tomada como proveniente de una
corriente de utilidad externa, (agua de una torre de enfriamiento, etc.), utilizando el
mtodo denominado en Aspen-HYSYS como From Utility Fluid.
Ecuaciones
Reaccin
6 7 () + 32() 6 13 ()
Sistema de unidades
Field
249
Ingeniera Qumica | Uniatlntico | 2012
a)
b)
c)
24.3
Dimensionamiento de equipos
2
= ,
= 2 = ,
Alguno de estos 2 valores se inserta en la ventana del reactor en la pestaa
Dynamics.
(. )
=
=
(. )
4,9343106
(150
)(420,063 2 )
2
= ,
(. )
252
Ingeniera Qumica | Uniatlntico | 2012
( )
= ,
4,9343106
=
(1
) (171,678 90)
= ,
(. )
(. )
= (
420,063 2
1
1
) (6 ) (
) (62,5 3 ) (
)
0.8
12
18
= ,
Para asignar estos parmetros, se hace doble clic en la corriente de energa y
en la pestaa Dynamics, se escoge el radio button de Utility Fluid.
Seguidamente se hace clic en Utility Valve y debe aparece la ventana
donde se le asignan los parmetros calculados como se muestra en la Figura
24.3 (el flujo mximo permisible del fluido frio es el doble del valor
calculado).
253
Ingeniera Qumica | Uniatlntico | 2012
24.4
24.5
refrigerante para disminuir esta temperatura. Los pasos de su instalacin son los
siguientes:
Instale una funcin de transferencia de la paleta de objetos y asigne como PV la
temperatura del reactor, Vessel Temperature.
En la pestaa Parameters, en la opcin Configuration, se asigna los lmites
permisibles de la temperatura del reactor. Para escoger los lmites de esta
variable controlada, se tiene en cuenta que el controlador que posee este lazo de
control, ser sintonizado por medio de heursticas, por esta razn se le asigna una
tolerancia de 100F del valor en estado estacionario que es 250F. Entonces los
lmites de PV, que son iguales a los del OP, son 150F y 350F
En la misma pestaa, pero en la opcin 2nd Orden, se insertaran 2 atrasos en
serie, a razn de las mediciones de temperatura del reactor y la chaqueta de
enfriamiento. Para realizar esto, se hace check en 2nd Orden, para active la
funcin de transferencia, y se le asigna un valor al atraso dinmico Time
constant de 1 min.
Instale un controlador PID, que posee como variable de proceso a controlar (PV),
la seal de salida de la funcin de transferencia OP Value y seal de salida del
controlador (OP), el porcentaje de abertura de la vlvula que regula la corriente
de energa Control Valve. Como se puede notar, esta vlvula no regula el flujo
de calor de esta corriente, sino el flujo de agua de enfriamiento, que ocasiona
esta transferencia de calor.
En la pestaa Parameters, colocar como accin del controlador directa e
insertar los lmites de PV son los mismos que el OP de la funcin de
transferencia, es decir 150F y 350 F.
Este controlador ser sintonizado como si fuese uno de accin proporcional, por
lo cual se asigna, en esta misma pestaa, una ganancia proporcional de 10 (al
final del mdulo se explica porque no se us un PID envs de un P).
Abrir la caratula de este controlador y cambiar el modo del controlador de
manual Man a automtico Auto
Como el integrador solo se pondr a funcionar, cuando sea instalado el control de
nivel, se contina con el proceso de instalacin del mismo.
255
Ingeniera Qumica | Uniatlntico | 2012
257
Ingeniera Qumica | Uniatlntico | 2012
258
Ingeniera Qumica | Uniatlntico | 2012
la simulacin End Time (aqu debe salir Non-Stop, es decir que el integrador
solo se detendr manualmente), entre otras opciones.
259
Ingeniera Qumica | Uniatlntico | 2012
Una vez hecho esto, se hace clic en Start y tiempo de simulacin empieza de cero
hasta 500 min, momento en el cual todas las variables controladas han alcanzado su
estado estacionario.
260
Ingeniera Qumica | Uniatlntico | 2012
(. )
100
250,94 250
100 = 0,47%
350 150
(. )
4(5%)
= ,
(0,47%)
261
Ingeniera Qumica | Uniatlntico | 2012
< <
(. )
En conclusin, los procesos con reaccin poseen una dinmica desafiante para su
control, puesto que su lazo de control de temperatura, es condicionalmente estable,
situacin que fcil superarla con el uso de un controlador de accin proporcional y
no uno de accin proporcional-integral- derivativo, adems la situacin presentada
con la sintonizacin del control de temperatura indica que para procesos complejos
como este se debe hacer un estudio ms profundo sobre la escogencia de los
parmetros dinmicos.
262
Ingeniera Qumica | Uniatlntico | 2012
Referencias bibliogrficas
43. Willian L. Luyben. Plantwide Dynamic Simulators in Chemical Processing and Control. New
York : Marcel Dekker, Inc, 2002. pp. 105-116.
44. Donald Q. Kern. Procesos de Transferencia de Calor. Colonia San Juan Tlihuaca, Delegacin
Azcapotzalco : Cesca, 1965. p. 945.
45. Michael L. Luyben, William L. Luyben. Essencials of Process Control. McGraw-Hill. pp. 320-323.
263
Ingeniera Qumica | Uniatlntico | 2012
25.2
Paquete Fluido
Componentes
Metanol y Agua
Ecuacin
Reacciones
Sistema de unidades
SI
posee como corriente de entrada Agua y de corriente de salida Sale_V4, este equipo
se especifica con una cada de presin de 68 Kpa.
Intercambiador: Instale un intercambiador de calor de la paleta de objetos y asigne
en la pestaa Design, opcin Connections como corriente de entrada en el lado
de los tubos a Sale_V4 y de salida a Sale_tubos. Paralelamente como corriente
entrada en el lado de la coraza a Sale_V2 y de salida a Sale_coraza. En la misma
pestaa, pero en la opcin Parameters, asigne como modelo de intercambiador de
calor Exchanger Design (End point), debido a que en este sistema no hay cambios
de fase y el Cp es relativamente constante (52) y como tipo de coraza segn Tubular
Exchanger Manufacturers Association (TEMA) de tipo E, lo quiere decir que la
coraza es de un solo paso (53).
Para las cadas de presin en el lado de los tubos y la coraza, se asignaran 7,2 Kpa y
68 Kpa, respectivamente. Asumiendo coeficiente global de transferencia de calor
(U) de 615 W/m2 C y un rea de transferencia (A) de 271,42m2, asignar como UA,
el producto entre estos valores, 166923,3 W/C.
En la corriente Sale_coraza se instala una vlvula de nombre V5, con una cada de
presin de 68 Kpa y una corriente de salida denominada Sale_V5.
Punto de Mezcla: se instala un punto de mezcla (MIX), en el que las corrientes
Sale_V3 y Sale_V5 se mezclan, formando una nica corriente llamada
Metanol_subenfriado a la salida del mismo. Para que esta mezcla pudiera darse se
haca necesario que ambas corrientes de entrada tuvieran la misma presin, razn
por la que en la vlvula V3 se asign una cada de presin de 204Kpa, producto del
paso del Metanol por la vlvula V2, P=68Kpa, el intercambiador, P=68Kpa y la
vlvula V5, P=68Kpa, (Ptotal=68Kpa x 3 = 204Kpa).
Despus de realizado lo anterior, el sistema en estado estacionario ha debido
converger satisfactoriamente, como se muestra en la Figura 25.1.
(a)
265
Ingeniera Qumica | Uniatlntico | 2012
(b)
Figura 25.1. Simulacin de un intercambiador de calor con derivacin en
estado estacionario: (a) Diagrama de flujo y (b) Especificaciones de las
corrientes
Tubos
Numero de tubos= 981
OD = 20 mm
Di = 16 mm
L = 4,830 m
Coraza
Dimetro = 894 mm
Pitch = 25 mm
Arreglo Triangular
Coeficiente de ensuciamiento =
(1/5000) m2 - C/W
Ktuberia = 50 W/m - C
Coeficiente de ensuciamiento =
(1/3000) m2 - C/W
Tabla 25.1 Datos del intercambiador de calor
266
Ingeniera Qumica | Uniatlntico | 2012
a)
b)
Figura 25.2. Parmetros de diseo del intercambiador a) Coraza y b)
Tubos
25.4
25.5
269
Ingeniera Qumica | Uniatlntico | 2012
El lector debe notar el valor de PV que sale en la caratula esta vacio (empty),
por tal motivo se enciende el integrador y en el momento que se alcance un
valor estable, que es aproximadamente 51,364 C se hacer clic en Start
Autotuner. En este momento se espera que hasta que el simulador obtenga
los parmetros dinmicos. Una vez se hayan obtenido, se detiene el
integrador. Los parmetros obtenidos son una ganancia proporcional de
1,05, un tiempo integral de 0,359 min y un tiempo derivativo de 7,98x10 -2
min
Active la caratula Face Plate de este controlador, active el integrador, dejar
pasar unos segundos, pasar el controlador de manual a automtico, dejar
pasar un par de segundos ms y apague el integrador.
Con respecto al controlador que manipula la vlvula V5 se repiten todos los pasos
anteriores solo con las diferentes modificaciones:
Para generar la seal del OP de la funcin de transferencia, solo se debe
permitir que el integrador aumente en 20 min del tiempo de simulacin actual
que posee, por ejemplo si el tiempo de simulacin est en 101,20 min, ingrese
en End Time, 121,20 min.
La accin del controlador debe ser directa.
271
Ingeniera Qumica | Uniatlntico | 2012
272
Ingeniera Qumica | Uniatlntico | 2012
(a)
(b)
observa que los valores de la caratula comienzan a oscilar, hasta obtener los
parmetros dinmicos. Una vez se hayan obtenido, se detiene el integrador.
En la Figura 25.9 muestra como son las variaciones del OP de las vlvulas V3 (lnea
morada) y V5 (lnea azul agua marina), y la temperatura de la corriente
Metanol_subenfriado (lnea verde), si se reinicia el integrador con este controlador
en modo automtico. En la figura se muestra cmo acta el controlador a los
elementos de control final, es decir, accin inversa a V3 y directa para la V5.
Tambin cabe destacar que la respuesta del controlador fue buena y rpida como es
de esperarse de esta estrategia de control.
Referencias bibliogrficas
49. Willian L. Luyben, Bjrn D. Tyreus, Michael L. Luyben. Plantwide Process Control. s.l. :
McGraw-Hill, 1998. pgs. 152-154.
50. Jonathan Love. Process Automation Handbook: a guide to theory and practice. 225-229. s.l. :
Springer, 2007.
51. R. K. Sinnott. Coulson & Richardson's Chemical Engineering. Cuarta edicin. Linacre House,
Jordan Hill, Oxford OX2 8DP 30 Corporate Drive, MA 01803 : Elsevier Butterworth-Heinemann,
2005. pgs. 675-679.
52. AspenTech. Aspen HYSYS, operations guide. Burlington, Ma : Aspen Technology, Inc, 2009.
pgs. 4-97 - 4-98.
53. Richard L. Shilling, Patrick M. Bernhagen, Victor M. Goldschmidt, Predrag S. Hrnjak, David
Johnson, Klaus D. Timmerhaus,. Perry's Chemical Engineers' Handbook. Octava edicin. s.l. :
McGraw-Hill, 2008. pgs. 11-34. Vol. XI.
54. Willian L. Luyben. Plantwide Dynamic Simulators in Chemical Processing and Control. New
York : Marcel Dekker, Inc, 2002. pgs. 36-37.
277
Ingeniera Qumica | Uniatlntico | 2012
En esta simulacin se estudiara como se separara un mezcla equimolar de THFagua por medio de dos columnas de destilacin, cada una con 10 platos un
condensador total y un rehervidor parcial, donde la columna de baja presin operara
a un presin de 1 bar y la columna de alta presin opera a 8 bar (60).
La simulacin contara adems de las dos columnas de destilacin, once vlvulas a
razn de los once lazos de control a instalar, de las cuales cinco son instalados
manualmente y las seis restantes son dimensionadas con los algoritmos internos de
las corrientes de energa y de reflujo. Tambin la simulacin contar con tres
bombas, una cada corriente de fondo de las columnas y una en el flujo de destilado
para aumentar la presin de esa corriente hasta la presin de la columna de alta
presin.
THF(C4H8O) y agua
Fase Vapor: Ideal Fase Lquida: Wilson
No hay. No es un proceso con reaccin
EuroSI
)
= + (
(. )
278
Sin embargo esta ecuacin es solo aplicable para el sistema ingls, porque esta
ecuacin parte del hecho que numricamente son iguales el peso especfico ( =
) y la densidad del sistema. Para ser ms claro se usara el siguiente ejemplo con
la densidad del agua @ 4C que es igual a 62,4Lbm/ft3.
= = 62,4
32,2 2
3
= ,
32,2
2
)
= + (
(. )
= 915,3
9,8 2
= ,
3
= 1 + 10 (
5
0,02540
1
1
)(
) (8,97103 3 ) (
)( 5 )
1
10
2
. = ,
279
Ingeniera Qumica | Uniatlntico | 2012
= 849
9,8 2
= ,
3
1 2
= 8 + 10 (
5
0,02540
1
1
)(
) (8,32103 3 ) (
)( 5 )
1
10
2
. = ,
(a)
281
Ingeniera Qumica | Uniatlntico | 2012
(b)
Figura 26.1. Simulacin de una destilacin de una mezcla azeotrpica
homognea en estado estacionario: a) Diagrama de flujo y b) Propiedades
de corrientes
prctica
de
nombre
DESTILACIN
282
Ingeniera Qumica | Uniatlntico | 2012
a)
b)
283
Ingeniera Qumica | Uniatlntico | 2012
Columna
de
Baja presin
Columna
de
Alta presin
2,193
12,6
1,897
16,42
1,436x104
7230
9751
1,537x104
Propiedades
Para calcular la velocidad mxima del vapor, para evitar el arrastre se aplica
la ecuacin 18.2 con F = 1,22 y el flujo de vapor es calculado con la ecuacin
18.1 (flujo volumtrico = flujo msico*densidad-1), el rea mnima con la
ecuacin 18.3 y el dimetro mnimo con la ecuacin 18.4. Los resultados
estn resumidos en la Tabla 26.2.
Columna
de
Baja presin
Columna
de
Alta presin
0,8238
0,3437
0,8858
0,3011
2,2079
0,4638
1,6120
0,8636
1,6766
0,7684
1,4326
1,0486
Propiedades
Columna
De
Baja presin
Columna
De
Alta presin
848,5
769,1
879,8
732,6
1,436x104
7230
1,068x104
1,893x104
Propiedades
Columna
De
baja presin
Columna
De
Alta presin
2,8207
1,5668
2,0232
4,3066
Propiedades
Vlvula
Valor
En
estado estacionario
Kcal/h x106
Tolerancia
Kcal/h x106
Rango
Kcal/h x106
Qc_baja
Qc_alta
Qr_baja
Qr_alta
1,8014
0,97902
1,7309
1,3775
1
0,4
1
0,5
0,8014 - 2,8014
0,47902 - 1,47902
0,7309 - 2,7309
0,8775 - 1,8775
286
Ingeniera Qumica | Uniatlntico | 2012
287
Ingeniera Qumica | Uniatlntico | 2012
a)
b)
289
Ingeniera Qumica | Uniatlntico | 2012
290
Ingeniera Qumica | Uniatlntico | 2012
a)
b)
Es necesario recalcar que ingresando de esta forma los datos la relacin ser
calculada como Flujo de reflujo/ flujo de destilado, pero si se hubiese ingresado al
revs el lector debe percatarse que la accin del controlador debe ser directa (62).
En la parte inferior izquierda se hace clic en Control Valve y agrege los
rangos del flujo del destilado. Como su flujo en estado estacionario es 67,06
Kgmol/h, se escoge una tolerancia de 20 Kgmol/h, es decir que los rangos a
insertar serian 47,06 y 87,07 Kgmol/h.
Rellene en la pestaa Parametes en las opciones Operation y
Configuration, como se muestra en las Figuras 26.6a y 26.6b
respectivamente. Es importante hacer check en Enable Ratio control,
porque esto le permite manipular la relacin entre el flujo de destilado y el
reflujo.
291
Ingeniera Qumica | Uniatlntico | 2012
a)
b)
La razn de que no uso el mtodo ATV para este controlador, es porque si este se
sintoniza como un PI, solo generara accin integral, y no se puede usar un PID
porque no es recomendable para este tipo de lazo de control, por tal motivo se us
las heursticas de Luyben para este lazo de control.
Active el integrador y cuando aparezcan valores de PV y OP en la caratula,
pase el controlador de modo manual a automtico, deja pasar un par de
segundos ms y detiene el integrador.
26.5.4 Control de nivel del tambor de reflujo de la columna de baja presin
Este controlador posee una accin directa, puesto que un aumento en el nivel de
lquido de tambor de reflujo, implicara un aumento en la en la abertura de la
vlvula V2, para disminuir el nivel de lquido. Los pasos de su instalacin son los
siguientes:
Instale un controlador PID, el cual posee como PV, el nivel del lquido del
condensador, Liquid Percent Level y como OP, el porcentaje de apertura de
la vlvula V2 Actuactor Desired Position.
En la pestaa Parameters, colocar como accin del controlador directa,
como el valor ideal es que el tambor este lleno en un 50%, por eso se toma
una tolerancia de 25 %, es decir que los lmites de PV a insertar serian 25 %
y 75%.
292
Ingeniera Qumica | Uniatlntico | 2012
Instale un controlador PID, el cual posee como PV, el nivel del lquido del
condensador, Liquid Percent Level y como OP, el porcentaje de apertura de
la vlvula V5 Actuactor Desired Position.
En la pestaa Parameters, colocar como accin del controlador directa,
como el valor ideal es que el tambor este lleno en un 50%, por eso se toma
una tolerancia de 25 %, es decir que los lmites de PV a insertar serian 25 %
y 75%.
En la misma pestaa Parameters, en la opcin Autotuner, escoger la
opcin de sintonizar un controlador PI, hacer clic en Start Autotuner y
encender el integrador. En este momento se observa que los valores de la
caratula comienzan a oscilar, hasta obtener los parmetros dinmicos. Una
vez se hayan obtenido, se detiene el integrador.
Active el Face Plate de este controlador, active el integrador, dejar pasar
unos segundos, pasar el controlador de manual a automtico, dejar pasar un
par de segundos ms y apague el integrador.
26.5.7 Control de nivel de lquido en base de la columna de alta presin
Este controlador posee una accin directa, puesto que un aumento en el nivel de
lquido en la base de la columna, implicara un aumento en la en la abertura de la
vlvula V4, para disminuir el nivel de lquido. Los pasos de su instalacin son los
siguientes:
Instale un controlador PID, el cual posee como PV, el nivel del lquido del
rehervido, Vessel Liq Percent Level y como OP, el porcentaje de apertura
de la vlvula V4 Actuactor Desired Position.
En la pestaa Parameters, colocar como accin del controlador directa,
como el valor ideal es que el tambor este lleno en un 50%, por eso se toma
una tolerancia de 25 %, es decir que los lmites de PV a insertar serian 25 %
y 75%.
En la misma pestaa Parameters, en la opcin Autotuner, escoger la
opcin de sintonizar un controlador PI, hacer clic en Start Autotuner y
encender el integrador. En este momento se observa que los valores de la
caratula comienzan a oscilar, hasta obtener los parmetros dinmicos. Una
vez se hayan obtenido, se detiene el integrador.
294
Ingeniera Qumica | Uniatlntico | 2012
298
Ingeniera Qumica | Uniatlntico | 2012
Referencias bibliogrficas
59. Warren D. Seider, J.D. Seader, Daniel R.Lewin, Soemantri Widagdo. Product and Process
Design Principles: Synthesis, Analysis and Evaluation. Tercera edicin. Westford : John Wiley &
Sons. Inc, 2009. pgs. 233-236.
60. Willian L. Luyben. Plantwide Dynamic Simulators in Chemical Processing and Control. New
York : Marcel Dekker, Inc, 2002. pgs. 227-242.
61. Page S. Buckley, Willian L. Luyben, Joseph P. Shunta. Design of Distillation Column Control
Systems. s.l. : Edward Arnold, 1985. pgs. 229-231.
62. AspenTech. Aspen HYSYS, operations guide. Burlington, Ma : Aspen Technology, Inc, 2009.
pgs. 5-81 - 5-102.
300
Ingeniera Qumica | Uniatlntico | 2012
301
Ingeniera Qumica | Uniatlntico | 2012
Agua: Instale una corriente de materia de nombre F2_agua, el cual posee 600
Lbmol/h de 100% de agua, una presin de 150 psia, 90F. Despus instale una
vlvula de nombre V2, la cual posee una corriente de entrada F2_agua y de salida
Sale_V2. Para que la vlvula quede totalmente especificada asgnele una cada de
presin de 40 psia.
Aire: Instale una corriente de materia de nombre F3_Aire, el cual posee 7 Lbmol/h
de 79% de nitrgeno y 21% de oxgeno en fraccin molar, una presin de 150 psia,
90F. Despus instale una vlvula de nombre V3, la cual posee una corriente de
entrada F3_Aire y de salida Sale_V3. Para que la vlvula quede totalmente
especificada asgnele una cada de presin de 40 psia.
Tanque de mezclado: Instale un tanque de la paleta de objetos y asigne de nombre
Mezclador, conctelo con las corrientes de entrada Sale_V1, Sale_2 y Sale_V3 y
como corrientes de salida Venteo y Lquido.
Instale una vlvula de nombre V4, como corriente de entrada Venteo, y de corriente
de salida Sale_V4, a esta vlvula se le ingresa una cada de presin de 40 psia.
Despus instale una bomba de nombre B1 la cual tiene corriente de entrada Lquido,
de salida Sale_B1, y como corriente de energa W_B1. Para que la bomba quede
totalmente especificada, asigne un aumento de presin de 50 psia. Por ultimo instale
una vlvula V5, la cual conecta como corriente de entrada Sale_B1 y de salida
Sale_V5, asigne una cada de presin de 90 psia.
A partir de este momento la simulacin ha podido converger satisfactoriamente
como se muestra en la Figura 27.1, por eso guarde esta simulacin de nombre
TANQUE DE MEZCLADO.
(a)
302
Ingeniera Qumica | Uniatlntico | 2012
(b)
Figura 27.1. Simulacin estacionaria (a) Diagrama de flujo y (b)
Propiedades de las corrientes
27.3
Dimensionamiento de equipos
2 2 (2) 3
=
=
4
4
2
2
= ,
= 2 = ,
3
303
Ingeniera Qumica | Uniatlntico | 2012
304
Ingeniera Qumica | Uniatlntico | 2012
305
Ingeniera Qumica | Uniatlntico | 2012
a)
b)
307
Ingeniera Qumica | Uniatlntico | 2012
Instale un controlador PID, el cual posee como PV, el nivel de lquido del
tanque Liquid Percent Level y como OP el porcentaje de abertura de la
vlvula V5 Actuactor Desired Position.
(. )
308
Ingeniera Qumica | Uniatlntico | 2012
a)
b)
Figura 27.4. Bloque selector en la pestaa:
a) Connections y b) Parameters
El lector debe recordar que se han diseado todos los procesos con un
porcentaje de abertura del 50% y como en Aspen-HYSYS por defecto la
dinmica de la vlvula de control es instantnea (ganancia de 1) es decir que
los controladores en estado estacionario generaran una seal de salida del
50%, pero como estos no estn conectados al FCE su seal seria cero. Por tal
motivo las entradas deben ser transformadas en entradas escaladas usando la
ecuacin 27.1.
= 1(0%) + 50% = 50%
= 1(0%) + 50% = 50%
La Figura 27.5 muestra cmo se deben ingresar los datos presentados en los
clculos previos:
309
Ingeniera Qumica | Uniatlntico | 2012
311
Ingeniera Qumica | Uniatlntico | 2012
En esta simulacin solo se mostr una de las variaciones del control selectivo, sin
embargo el bloque selector que presenta Aspen-HYSYS, usado los mismos
conceptos explicados en este mdulo, puede usarse como instrumentacin
redundante o auctioneering, por ejemplo se tiene una columna de destilacin con el
selector captura las temperaturas del plato donde quiere implantar un control de
temperatura y los platos superior e inferior del mismo y a esos valores el selector le
puede aplicar un promedio ponderado o la mediana y ese sera el PV del controlador
de temperatura.
Referencias bibliogrficas
66. Aspentech. Aspen HYSYS, operations guide. Burlington, Ma : Aspen Technology, Inc, 2009.
pgs. 5-218 - 5-224.
312
Ingeniera Qumica | Uniatlntico | 2012
313
Ingeniera Qumica | Uniatlntico | 2012
314
Ingeniera Qumica | Uniatlntico | 2012