Design and Implementation of A Temperature Predictive Controller For A Dryer Plant of Cocoa Beans
Design and Implementation of A Temperature Predictive Controller For A Dryer Plant of Cocoa Beans
Design and Implementation of A Temperature Predictive Controller For A Dryer Plant of Cocoa Beans
Conferencia de Asia y el Pacífico de 2015 sobre ingeniería de sistemas asistida por computadora
Licencia de uso autorizada limitada a: UNIVERSIDAD POLITECNICA SALESIANA. Descargado el 13 de diciembre de 2022 a las 15:16:48 UTC de IEEE Xplore. Se aplican restricciones.
Machine Translated by Google
se ensayó un horno a gas con charolas perforadas dispuestas verticalmente y misma dirección en la que se mueve el cacao. La figura 1 muestra el sistema
se utilizó una muestra de cacao con una masa aproximada de 13kg. de secado.
Para el segundo experimento se utilizó un horno con flujo de aire a la entrada y
un extractor de aire a una temperatura de 55°C. Luego de analizar
computacionalmente las simulaciones y el comportamiento de las variables de
temperatura y humedad, los autores concluyeron que el cacao debe ser agitado.
en el horno para evitar quemaduras por el contacto entre los granos y las
bandejas de metal. Recomendaron no eliminar grandes cantidades de vapor de
agua en el aire ya que el aire seco deshidrata la superficie de los granos,
provocando quemaduras.
Los autores CL Hii et al. [10] realizaron una investigación para estudiar la
cinética de secado de las capas individuales del grano de cacao; estos son la
cubierta o cáscara de la semilla y el cotiledón. Una fina capa de granos de
Figura 1. Planta de secado de cacao.
cacao se secó utilizando una bomba de calor y aire seco a temperaturas de
28,2 °C, 40,4 °C y 56 °C, registrando durante todo el proceso el contenido de En el sistema de dosificación, que se describe en la TABLA I, se deposita
humedad de la cubierta de la semilla y el cotiledón. Observaron que la reducción cacao fermentado con una humedad de 0.35 en base seca, equivalente a 35%
del contenido de humedad en la cubierta de la semilla fue más rápida que en el de humedad. En la parte inferior del dispensador, la gravedad y la vibración se
cotiledón durante las dos primeras horas del proceso. Luego, el paso de secado utilizan para medir el flujo de granos de cacao de forma continua hacia la
intermedio registró mayor contenido de humedad en la cabeza con respecto al cámara de secado.
cotiledón. Finalmente se encontró mayor humedad en la cubierta seminal en el TABLA I. CARACTERÍSTICAS DEL DISPENSADOR.
cotiledón, demostrando la migración de humedad desde el cotiledón hacia la el menor, el mayor Tipo
testa. En sus conclusiones, los autores argumentan que el comportamiento del
Capacidad de dosificación 0,002 kg/s 0,020 kg/s Vibrador
contenido de humedad registrado en la cubierta de la semilla y el cotiledón
Almacenamiento 100 kg. capacidad del conducto
explicado al postular la transferencia de masa de Luikov en este caso mostró
que el potencial de transferencia de masa para la testa es siempre menor que
en el cotiledón, por lo tanto, la transferencia de masa a la la cabeza del cotiledón La generación de aire caliente (gas de secado) se realiza con una cámara
se ve afectada, aunque hay menos contenido de humedad del cotiledón. de combustión que toma aire del ambiente y lo
se calienta por la combustión de gas licuado de petróleo (GLP - C4H10), el
aire calentado se envía al cilindro rotativo. Ver las características del generador
de aire caliente en la TABLA II.
En este trabajo mostramos los resultados de la identificación de una planta TABLA II. GENERADOR DE CARACTERÍSTICAS.
secadora de cacao, la misma que se encuentra instalada en el campus mínimo máximo Tipo 0,28 m
centenario de la Universidad Politécnica Salesiana del Ecuador. Diámetro - -
0,53 m. 31,2 kilovatios
En los siguientes párrafos del artículo se hace una descripción de la planta - -
Longitud
a controlar, los resultados del modelado y simulación de la planta mostrada en Poder 12,4kW -
lazo abierto y finalmente el diseño e implementación de un controlador predictivo Flujo de aire -
0,1 kg/s 0,8 kg/s 0,0001 kg/s
detallado para secado de gas a temperatura variable. Además se realiza un El 0,0007 kg/s BPL
análisis comparativo de los resultados obtenidos utilizando el controlador consumo de combustible
predictivo versus los resultados obtenidos utilizando un controlador PID. 40 ÿC 120 ÿC -
Temperatura de salida
En este trabajo se presenta un aporte significativo al secado del cacao. Se El cilindro giratorio (cámara de secado), recibe flujos de cacao y aire
demuestra la factibilidad de implementar una técnica de secado artificial caliente. Internamente el efecto de transferencia de calor se da por convección,
utilizando un cilindro rotatorio y la transferencia de calor por convección como el grano se mueve debido al giro del cilindro ya la velocidad del gas de secado.
alternativa al proceso de secado natural, basado en técnicas artesanales. El cacao permanece por un período de aproximadamente 10 horas (tiempo de
residencia) dentro del cilindro para lograr un contenido de humedad de 0,08 en
Las técnicas de secado artesanal no están controladas y el éxito de este base seca; valor que permite el almacenamiento, transporte y procesamiento.
proceso depende de simples observaciones y análisis empíricos. La
investigación desarrollada en torno a la implementación de un algoritmo de
La figura 2 muestra esquemáticamente la sección transversal del cilindro
control para la planta piloto de secado artificial garantiza que los granos de
que se muestra en la TABLA III y se describen las características de la cámara
cacao obtengan un nivel de humedad óptimo, sin dependencia directa de las
de secado.
condiciones climáticas de la época de cosecha.
CUADRO III. Características de la cámara de secado
mínimo máximo Tipo
II. DESCRIPCIÓN DE LA PLANTA
Velocidad de 0,5 rpm cilíndrico de 4,0 rpm
La planta piloto consta de un sistema de dosificación de granos, un rotación
ventilador de cámara de combustión de empuje con una cámara de secado, el Longitud 3,30 m.
cual consta de un cilindro rotativo por donde circula aire caliente en el diámetro 0,6 metros
203
Licencia de uso autorizada limitada a: UNIVERSIDAD POLITECNICA SALESIANA. Descargado el 13 de diciembre de 2022 a las 15:16:48 UTC de IEEE Xplore. Se aplican restricciones.
Machine Translated by Google
3.1 Identificación.
tercero DISEÑO DEL CONTROLADOR El punto de funcionamiento nominal de la planta es de 0,00015 kg/s de
caudal másico de combustible (GLP - C4H10), valor que se alcanza
La estrategia de control PID es la más utilizada en los procesos industriales
estableciendo una temperatura del gas de secado en el interior de la cámara
en general, sin embargo, los controladores predictivos son los más utilizados
de 54,85 °C en condiciones ambientales de 70% de humedad relativa y una
en secciones económicamente importantes dentro del área industrial [14]. El
temperatura de 30° C. Durante el día la temperatura ambiente varía unos
nivel de aceptación de los controladores predictivos actualmente es bastante
pocos grados centígrados. Si la planta opera en circuito abierto, tal cambio
alto, una muestra de esta aceptación es la caja de herramientas desarrollada
para MATLAB® por Bemporad et al [15], llamada Model Predictive Control provoca cambios en el nivel de la temperatura del gas de secado dentro de la
204
Licencia de uso autorizada limitada a: UNIVERSIDAD POLITECNICA SALESIANA. Descargado el 13 de diciembre de 2022 a las 15:16:48 UTC de IEEE Xplore. Se aplican restricciones.
Machine Translated by Google
la respuesta en lazo abierto. En la TABLA IV se tabulan los datos para la ruido blanco con valor medio cero. La función de costo para este
generación de los datos. la estrategia de control viene dada por la ecuación (5).
R=1
Figura 4. Datos de salida. Modelo Identificado y Planta.
S = [1.7020 -1.5221 0.0767]
(2) y (3) representan la estructura matemática para el modelo de
tiempo discreto identificado por ARX. T = 0,2566.
() = 1 ÿ 1,038 + 0,05276 (2) Donde R, S y T son los valores que puede implementar el controlador
predictivo generalizado. Otra opción para calcular el controlador, utilizando
() = 0.00578 (3)
el mismo algoritmo, es obtener la ley de control, como se muestra en la
3.2 Diseño GPC. ecuación (6).
El algoritmo básico de control predictivo generalizado (GPC) fue () = (ÿ1) ÿ 1,7020() + 1,5221(ÿ1) ÿ 0,0767(ÿ2) + 0,0072(+1)
propuesto por Clarke et al. [12] en el año 1987. Actualmente este método + 0,0147(+2) + 0,0221 (+3) + 0,0289(+4) + 0,0358(+5) +
0,0427(+6) + 0,0494(+7) + 0,056(+8) (6)
es muy utilizado en la academia y la industria. Con esta estrategia predictiva
se calcula una secuencia de acciones de control futuras minimizando una
función de costo. Una función cuadrática, por un lado, mide la distancia 3.3 Diseño del controlador PID.
entre la salida del sistema pronosticada y la ruta de referencia dentro del
horizonte de pronóstico y, por otro lado, se minimiza el esfuerzo de control Para el cálculo del controlador PID se eligió el método de la curva de
necesario para obtener esta salida. [13]. A diferencia de otras estrategias reacción de Ziegler-Nichols. En la Figura 5, el gráfico muestra los datos
de control predictivo, el GPC puede proporcionar soluciones explícitas en obtenidos experimentalmente en la planta de cacao.
ausencia de restricciones y puede trabajar con procesos inestables o de Se procede a calcular las constantes k0, ÿ0 y ÿ0, con
fase no mínima. ecuaciones (7), (8) y (9).
= = ..
= 1,2 (7)
Para los procesos SISO (una entrada y una salida), al estar linealizados
alrededor de un punto de operación dado, se puede utilizar un modelo de = ÿ = 13 ÿ 8 = 5 (8)
control autorregresivo integrado con medias móviles (CARIMA), como se
describe en la ecuación (4). ÿ = ÿ = 164 ÿ 13 = 151 (9)
205
Licencia de uso autorizada limitada a: UNIVERSIDAD POLITECNICA SALESIANA. Descargado el 13 de diciembre de 2022 a las 15:16:48 UTC de IEEE Xplore. Se aplican restricciones.
Machine Translated by Google
= 2 = 2(5) = 10 (11)
206
Licencia de uso autorizada limitada a: UNIVERSIDAD POLITECNICA SALESIANA. Descargado el 13 de diciembre de 2022 a las 15:16:48 UTC de IEEE Xplore. Se aplican restricciones.
Machine Translated by Google
[9] Tinoco H., Ospina D. Análisis del proceso de deshidratación del Cacao
para la disminución del tiempo de secado. Revista EIA, ISSN 1794-1237
Número 13, p. 53-63. (Julio 2010). Escuela de Ingeniería de Antioquia,
Medellín-Colombia.
[10] CL Hii, CL Law, S. Suzannah. Cinética de secado de la capa individual
de granos de cacao durante el secado con bomba de calor.
Revista de Ingeniería de Alimentos. ELSEVIER. (2011).
[11] Ogata, K. Ingeniería de Control Moderna (5ª Edición). Respuesta al
escalan. pag. 136. 2007 [12] Clark, DW., Mohtadi, C., Tuffs, PS. Control
predictivo generalizado. Parte I. El algoritmo básico. Automática 23 (2),
Figura 10. Acciones de control GPC y PID. 137-148, 1987. 3869, 1987.
207
Licencia de uso autorizada limitada a: UNIVERSIDAD POLITECNICA SALESIANA. Descargado el 13 de diciembre de 2022 a las 15:16:48 UTC de IEEE Xplore. Se aplican restricciones.