Aa Espel Emi 0372.en - Es
Aa Espel Emi 0372.en - Es
Aa Espel Emi 0372.en - Es
[1] Departamento de Eléctrica y Electrónica, Universidad de las Fuerzas Armadas ESPE, ID:
1
II. METRO ATEMÁTICO METRO ODELING valor deseado 0,5]
Derecha - El pistón está ligeramente [0.4 1
A. Sistema de control de posición
contraído. 1.6]
El modelo matemático del control de posición de un Derecha + El pistón es muy [1 2 3]
cilindro de doble efecto se compone de un mecanismo que contratado
consiste en utilizar frases "sí-entonces" que se denominan
"reglas". Estas reglas son útiles, ya que se refieren a variables Para la salida existen dos variables lingüísticas, "Sol A +" y
y adjetivos que describen las variables. El orden en que se "Sol A-", que permiten la apertura y cierre progresivo del
utilizan las reglas no tiene relevancia. [11] solenoide hasta alcanzar el valor ingresado en la consigna.
En el diseño del modelo matemático se utiliza el método de Estas variables son las respuestas del controlador a la
inferencia tipo Mamdani, este controlador está formado por variable de entrada "error", tablas II y III.
fusor, reglas difusas y defusor, fig. 1. [12] TABLA II.
DESCRIPCIÓN DE LAS FUNCIONES DEL TRIANGULAR
TIPO DE ADHESIÓN DE LA VARIABLE "SOL A +".
Functio Descripción Paramete
norte rs
Mantener No abra el solenoide [-5 0 2]
Fig. 1 Diagrama del controlador difuso tipo Mamdani [13] en
Abierto - Abra el solenoide a un valor [2 4 6]
Este controlador se desarrolla utilizando el complemento "Fuzzy
de 4
Logic Toolbox" en el entorno MATLABTM SIMSCAPETM. [14] Para
Abrir + Abra el solenoide a un valor [5 10 15]
acceder al complemento, escriba en la ventana de comandos "fuzzy",
de 10
desde allí puede desarrollar un controlador difuso que puede ser
TABLA III.
Mamdani o Sugeno, agregar variables de entrada y salida, y luego
DESCRIPCIÓN DE LAS FUNCIONES DEL TRIANGULAR
crear un conjunto de reglas. [15]
TIPO DE ADHESIÓN DE LA VARIABLE "SOL A-".
El controlador difuso tipo Mamdani del sistema tiene una
Functio Descripción Paramete
entrada y dos salidas, fig. 2.
norte rs
Mantener No abra el solenoide [-5 0 2]
en
Abierto - Abra el solenoide a un valor [2 4 6]
de 4
Abrir + Abra el solenoide a un valor [5 10 15]
de 10
2
Los datos del subsistema se utilizarán para la
Estas reglas con sus parámetros se exportan al espacio de trabajo implementación física.
de MATLABTM, para su uso posterior en la programación del PLC
S7-1200 para su implementación.
3
Entorno de la caja de herramientas del sistema. [21] cambia la amplitud del relé, que es intrínsecamente pequeña. Para
Para importar datos es necesario implementar el circuito obtener los parámetros PID se calculan utilizando las reglas de ajuste de
hidráulico de control de presión, y obtener los valores Ziegler-Nichols [26], tabla IV.
medidos de entrada-salida del sistema generados en el Data
logging. [22], los mismos que se importan al espacio de
trabajo de MATLABTM.
En el complemento del sintonizador PID, seleccione "Identificar nueva
planta", después de esto, dado que los valores son arbitrarios en la TABLA IV.
implementación, seleccione "Datos de E / S arbitrarios" donde ingresa los PARÁMETROS DEL CONTROLADOR PID DE ZIEGLER-NICHOLS [25]
datos de entrada (punto de ajuste) y salida (sensor), el mismos que están Controlle Kp Ti Td
por separado en el espacio de trabajo. r
PAG - -
Pi -
PID
La función de transferencia obtenida del tercer grado Fig. 11. Respuesta del sistema con método de relé, SIMSCAPETM.
representa la relación que describe la dinámica del sistema
considerado. [23] Para obtener las constantes PID de la
Los parámetros en el gráfico: Tu = 9.15 y a = 2.095. El
función de transferencia, se aplica el método Relay, es una de
resultado es el gráfico de comportamiento PID, fig. 12, con los
varias formas de obtener los parámetros (Kd, Tp, Ti). Este
parámetros calculados del método de Åström y Hägglund junto
método también se conoce como método de Åström y
con la regla de Ziegler-Nichols.
Hägglund. [24], que determina el punto final utilizando una
señal de prueba de relé como entrada del sistema. [25]
4
Los valores de los parámetros calculados corresponden a
Kp = 3.6465, Ti = 0.4575, Td = 0.1143. Al obtener una
respuesta adecuada con los parámetros calculados, la
simulación del sistema físico se puede realizar en el entorno
SIMSCAPETM de MATLABTM.
Se crea el modelo de control de lazo cerrado, que consiste en
la entrada de valores de consigna mediante el bloque "Uniform
Random Number" que genera números aleatorios de manera
uniforme. En el bloque de control se encuentra el controlador
PID al que se le ha programado el control proporcional integral
derivado con sus respectivos parámetros previamente
calculados, Kp, Ti y Td. A continuación, el bloque del sistema
hidráulico, compuesto por un motor hidráulico, una válvula
limitadora de presión proporcional que permite el paso de fluido
hidráulico para aumentar o disminuir la presión al sistema, una
Fig. 15. Diagrama de bloques del sistema de control de posición.
bomba hidráulica y un sensor de presión en el que se destaca la
precisión y estabilidad de Se puede visualizar el sistema
El sistema de control de posición controla un pistón de doble
controlado por control derivativo proporcional integral PID.
efecto mediante una válvula 4/3 que se controla mediante
Finalmente, se conecta al "Scope" con las dos señales del sensor
lógica difusa. El objetivo del sistema es ingresar datos
de presión y del punto de ajuste, fig. 13.
mediante teclado (setpoint) y mediante control difuso para
entregar una respuesta que logre los valores ingresados
por teclado.
5
1 Maestro 1500 S7-1500 [CPU 1511-1 PLC_POSICION y PLC_PRESION para que ambos controlen el
PN] proceso de extrusión de plástico.
2 PLC_POSICIÓN S7-1200 [CPU 1214C Si esta usando Palabra, use el Editor de ecuaciones de
DC / DC / DC /] Microsoft o la MathType Añadir
(http://www.mathtype.com) para las ecuaciones en su artículo
3 PLC_PRESIÓN S7-1200 [CPU 1214C (Insertar | Objeto | Crear nuevo | Ecuación de Microsoft o
DC / DC / DC /]
Ecuación MathType). "Flotar sobre el texto" debería no ser
seleccionado.
4 HMI_1 [TP700 Comfort]
interfaz PROFINET se comunica con el bloques de programa ordinarios y se configuran las entradas y
salidas del bloque, fig. 20.
6
La Fig. 20 muestra el bloque "Subsystem DB_1", que ha
sido programado en MATLABTM; este bloque se activa
mediante el contacto "Crtl fuzzy" que permite encender el
controlador. En la entrada del bloque, los valores de la
consigna del sistema se introducen de forma local o remota
mediante la variable "Setpoint pistón" y los valores del sensor
de posición del sistema mediante la variable "Distancia".
Mientras que en la salida se obtiene la respuesta del
controlador difuso, con las variables "Salida Pistón" y
"Entrada Pistón", que permiten el paso y retorno del fluido
hidráulico respectivamente en el esquema físico. Fig. 21. Parámetros PID introducidos manualmente, TIA Portal.
La programación del controlador de posición contiene
entradas y salidas analógicas que se envían y reciben del La optimización fina del sistema se lleva a cabo dentro de
maestro, que luego se utilizan para realizar el sistema SCADA y la la configuración del PID_Compact, dando valores de las
configuración de la pantalla HMI. constantes, fig. 22.
7
Los resultados obtenidos en la simulación del proceso en
respuesta a números aleatorios uniformemente en el proceso de
control de posición son satisfactorios, fig. 28, un error de precisión de
Se obtiene un 0,70% admisible en un proceso industrial.
Tss [s] 36 47
Fig. 29. Curvas de ajuste para PID automático y manual
ess [%] 0 0,1 modelos de aplicación.
8
los sobrepulsos en ambos casos varían en 1%, los datos calculados y
los generados automáticamente tienen un porcentaje en rango
moderado. El tiempo de estabilización Tss correspondiente a
PLC_Manual es de 43s, mientras que el Tss de PLC_Automatic es de
73s, es decir, la estabilidad del sistema se encuentra en un rango
aceptable en respuesta a un cambio de estado. El error en estado
estable es en ambos casos 0%, lo que muestra la excelente precisión
del sistema, tabla VII.
Fig. 31. Implementación del control PID de la presión
TABLA VII. variable.
ERRORES EN EL AJUSTE PID AUTOMÁTICO Y MANUAL
CURVAS Dados los porcentajes de precisión en los controladores de
PLC_Automatic PLC_Manual posición y presión, se puede observar que el controlador PID
tiene mayor precisión, con una diferencia de 0.26%, que en los
Tss [s] 73 43
procesos industriales es mínima.
ess [%] 0 0
V. C ONCLUSIÓN
REFERENCIAS
9
[6] A. Orjuela, O. Ramos y D. Amaya., «Control PID para [22] Soporte SIO, adquisición de datos de proceso y
producto destilado y concentración de fondo en una Seguimiento, 2017.
columna de destilación binaria», IREME International [23] Bishop, RC y D. -. RJ, Modelos matemáticos de los
Review of Mechanical Engineering, vol. 11, nº 4, 2017. sistemas, Sistemas de control moderno, Madrir:
[7] C. Lapusan, R. Balan, O. Hancu y A. Plesa., PEARSON EDUCATION SA, 2007, pp. 52-62.
«Desarrollo de un edificio multisala [24] R. Alanen, control de autoajuste basado en Åström's y
modelo termodinámico utilizando SIMSCAPE Library, » Teoría de Hägglund, Universidad de Oulu, 1993.
Procedia de energía, vol. 85, págs. 320-328, 2016.
[25] S. Levy, S. Korotkin, K. Hadad, A. Ellenbogdn, M.
[8] W. Mendes, F. Meneghetti, R. Dutta y D. Heeren., «Sistema de Arad y Y. & Kadmon, «PID Autotuning using
control difuso para riego de tasa variable mediante teledetección», Sistema feedback,» Ingenieros eléctricos y electrónicos de la
experto con ubicaciones de aplicaciones, vol. 27a Convención IEEE en Israel, págs. 1-3, 2012.
124, págs.13-24, 2019.
[26] K. Ogata, de Ingeniería de control Moderna, Madrid,
[9] C. Lorenzini, AS Bazanella, LFA Pereira y GRG Pearson Education SA, 2003, págs. 683-690.
D. Silva., «El método de oscilación forzada generalizada para
[27] FA Fernandes, GS Sestillo, AL Dias, D. Brandão y
sintonizar controladores PID,» Transacciones ISA, vol. 87,
P. Ferrari., «Influencia de los parámetros de la red en el tiempo de
págs.68-87, 2019.
recuperación de una red PROFINET de topología en anillo»,
[10] W. Sánchez, Ó. Proaño, P. Salazar, E. Salazar, W. Loza Archivo de documentos de la conferencia de IFAC, vol. 49, nº 30,
y JE Sánchez., «Sistema distribuido SCADA utilizando págs.278-283, 2016.
múltiples protocolos de comunicación industrial en
[28] Soporte SIO, «Control PID con PID_ Compact»,
equipos Allen Bradley», vol. 12, nº 14, 04 2019.
2014, págs. 2-50.
[11] OP Verma y H. Gupta, «Agua basada en lógica difusa
[29] M. Rasmussen y K. Laumann, «Potential use of HMI
Sistema de control de temperatura del baño » IJARCSSE, vol. 2, nº 4,
métodos de evaluación en HRA, » FABRICACIÓN
págs. 334-335, 04 2012.
DE PROCEDIOS, págs. 1359-1363, 2015.
[12] R. León, «Diseño, programación e implementación de
[30] S. McGrady, «Introducción», de Diseñando Scada
un sistema de control difuso para el proceso de
Software de la aplicacion, Canadá, Elsiever, 2013, págs. 1-10.
producción de hidrogeno a partir de vapor
sobrecalentador de agua, »Escuela Politécnica
Nacional, Quito, 2014.
[13] S. Ramin, I. Wan y I. Wan, «Diseño del sistema de control para un
robot móvil agrícola para estudio de campo, » Conferencia
internacional sobre ingeniería agrícola y alimentaria para la
vida (Cafei 2012), 2012.
[14] Matlab, «MathWorks», 2019. [En línea]. Disponible:
https://www.mathworks.com/products/fuzzy-logic.html
? s_tid = srchtitle.
[15] D.Guzmán y VMCataño, «La lógica difusa en
ingeniería: Principios, aplicaciones y futuro., » Ciencia y
Tecnología, nº 87-107, págs. 88-106, 2006.
[16] D. Cañete, JF Saz-Orozco, P. d. Moreno-Boza y D.
Duran-Venegas, «Modelado orientado a objetos y
simulación del sistema cardiovascular de circuito cerrado
mediante SIMSCAPE» Informática en Biología y Medicina, vol.
43, nº 4, 2013.
[17] Matlab, «MathWorks», 2019. [En línea]. Disponible:
https://www.mathworks.com/help/simulink/slref/satura
tion.html? searchHighlight = saturation & s_tid = doc_srcht
itle.
[18] MATLAB, «MathWorks», 2019. [En línea]. Disponible:
http://www.mathworks.com/.
[19] MATLAB, «MathWorks», 2019. [En línea]. Disponible:
https://la.mathworks.com/discovery/pid-control.html?f
bclid = IwAR0fRvjFz1z98RclWJTBEV0AtbUxWdmlA
1z_IHML82yhgCWxKtT648msKmg.
[20] SW Sung, J. Lee y I.-B. Sotavento.,
«Control Proporcional-Integral-Derivado», de Identificación
de procesos y control PID, Singapur, John Wiley & Sons,
2009, págs. 111-147.
[21] W. Sánchez, Á. Mullo y P. Velasco., Autosintonía PID
de controladores industriales, Caracas: Universidad
Central de Venezuela, 2019.
10