Soft Computing
Soft Computing
Soft Computing
TESIS DOCTORAL
Autor:
Raúl Mario del Toro Matamoros
Director:
Dr. D. Rodolfo E. Haber Guerra
¡Hemos llegado al final del largo túnel, ya podemos acariciar la luz! Aunque ahora
comienzan nuevos retos, los de verdad, los que marcarán nuestro paso por esta vida. Hay
muchas personas que han incidido directamente en la consecución de este inmenso recorrido,
ya sea desde el punto de vista profesional o personal, y que no quiero dejar de mencionar y
agradecer en este pequeño momento de gloria.
A Rodolfo, director de esta Tesis, tutor de toda la carrera investigadora y amigo. Siempre
dirás que no te debo nada, que el esfuerzo ha sido mío (lo repites con frecuencia), pero debo
agradecer enormemente tu paciencia, amistad y, sobre todo: tu confianza desde el principio, tu
conducción, tu protección y lealtad frente a todos los obstáculos que ha habido que sortear.
A los colegas (amigos) del equipo investigador de GAMHE y C4LIFE, sin este gran
equipo hubiese sido imposible llevar adelante este proyecto. Muchas gracias a Fernando,
Alfonso, Maritza (la amiguita de la eterna sonrisa) y algunos que ya no están como Pedro,
Bruno, Diego y Michael. A Agustín, por su paciencia infinita y por ir siempre despejando el
camino. Muchas gracias a Salvador Ros, Ángel Alique, José Ramón Alique, Miguel Ramírez
y Jorge Pamies por todo su apoyo. Especial mención a Rodolfo Haber Haber, por sus
enseñanzas, confianza e ideas.
Gracias a los colegas del instituto, mejor dicho “intituto”, creo que así tiene más ritmo,
que le vamos a hacer, vengo de donde vengo, a la comunidad de “La del viernes”, los que
compartimos muchos mediodías, cumpleaños y largos debates. A Miguel y Andrés por su
ánimo en las largas tardes, a Gregorio por sus chistes y cafés y a mucha gente del instituto que
hacen que el trabajo sea menos difícil.
A los amigos y amigas, no los menciono uno por uno, cuando la amistad es verdadera las
palabras sobran. Ustedes saben quienes son, les agradezco estar ahí y la comprensión ante mi
escasa comunicación durante todo este tiempo. Para los que han querido emprender también
este camino y por circunstancias de la vida han tenido que desistir (Jorge, Héctor, Garzón…)
y para los que aún lo intentan en condiciones menos propicias (Guillermo, Elsy, Frank…),
este logro también es de ustedes.
A mi madre, a ella debo mi educación y todo lo que soy, a mi abuela Benicia, mi segunda
madre, a ustedes dos está dedicado este momento. A mi padre, mis abuelos Julián y Emelina y
a toda la gran familia.
Muchas gracias a la tierra que me vio nacer y a mucha de su gente, por la educación
básica, la formación y por moldear mi actitud ante los problemas y ante la vida, por
enseñarme a sonreír inclusive en los momentos más tristes.
A la familia López-Sánchez, por acogerme como uno más. Mónica, muchas gracias por la
paciencia (han sido muchísimas horas sacrificadas), por el empuje, por el amor incondicional,
por la compañía, por aportar otros criterios y diferentes puntos de vista, este es: nuestro
momento.
RESUMEN
V
de entropía cruzada (Cross Entropy) y métodos de búsqueda directa como el algoritmo de
Nelder-Mead.
Para el ajuste óptimo de las estrategias de control inteligente, se ha utilizado, como
función objetivo a minimizar, la integral del error absoluto ponderado en el tiempo (ITAE).
Esta función penaliza las respuestas transitorias temporales de larga duración y es muy
sensible ante variaciones de parámetros, en comparación con otros criterios como la integral
del error absoluto (IAE) o la integral del error cuadrático ponderado en el tiempo (ITSE), entre
otros.
Además de las aportaciones expuestas anteriormente, otras contribuciones de esta Tesis
Doctoral constituyen los modelos computacionales eficientes que, diseñados e implementados
mediante una herramienta de alto rendimiento como Matlab/Simulink, permiten emular el
comportamiento dinámico o estático de variables representativas en procesos de fabricación.
Como parte de la metodología científica de la Tesis Doctoral, todas las estrategias
desarrolladas son validadas rigurosamente en una plataforma experimental, utilizada como
soporte tecnológico, lo cual permite el aprovechamiento mutuo de las experiencias alcanzadas
y la mejora de los resultados obtenidos. La validación de las estrategias mediante su
implementación en tiempo real, actuando directamente sobre operaciones de taladrado de alto
rendimiento, ha permitido corroborar los resultados obtenidos en simulación. Además, ha
dado la posibilidad de realizar estudios comparativos entre las estrategias basadas en técnicas
de soft-computing y las estrategias lineales para comprobar su desempeño, así como la
efectividad de los métodos metaheurísticos empleados en el ajuste de los sistemas de control
inteligente. Los estudios comparativos se han basado en la cuantificación del cumplimiento de
objetivos de control, a partir de parámetros que caracterizan el comportamiento dinámico de
sistemas de control, como la amplitud máxima de las oscilaciones (máximo sobrepaso) que
pueden aparecer en la variable controlada (ej., fuerza de corte) y los valores medios de su
error respecto a su valor deseado (ej., índices de desempeño ITAE, IAE, ITSE, entre otros).
Los estudios comparativos mostraron que, en el caso de la aplicación de las técnicas de
optimización para el ajuste de sistemas de inferencia borrosa, las respuestas temporales de los
sistemas de control, obtenidas en simulación son similares a las experimentales. Esto se debe
a que en todos los casos se han obtenido parámetros óptimos o sub-óptimos de los sistemas de
control borrosos. Pero a pesar de su similitud, el desempeño en general del sistema de control
optimizado mediante el método de entropía cruzada, es ligeramente superior al resto de
técnicas utilizadas, además de ofrecer una muy buena relación tiempo de cómputo y total de
parámetros de ajuste del método requeridos.
Por otra parte, los estudios comparativos del desempeño en tiempo real de las estrategias
de control no lineal empleadas, basadas en técnicas del soft-computing, controlando
operaciones de taladrado de alto rendimiento, han corroborado que, en general, resultan más
VI
efectivas que las estrategias lineales para hacer frente a la complejidad intrínseca de estos
procesos.
Algunos de los resultados parciales alcanzados durante las diferentes etapas de desarrollo
de esta Tesis Doctoral, han sido sometidos a la consideración de la comunidad científica
internacional. Investigadores de diferentes campos tales como las ciencias de la computación,
la inteligencia artificial, la automática, la ingeniería de fabricación, la ingeniería mecánica y
los procesos de transformación de metales, entre otros, han valorado de forma positiva varias
de las aportaciones realizadas.
VII
ABSTRACT
IX
penalizes the long time duration transient responses. Likewise, it is very sensitive to changes
in control system parameters, compared with other criteria such as integral absolute error
(IAE) or the integral of the time weighted square error (ITSE), among others.
In addition of the contributions discussed above, the design and implementation of
efficient computational models constitute other valuable contributions of this Thesis. These
models emulate the dynamic or the static relations of the most representative variables in
manufacturing processes. Being possible the implementation due to the use of a high
performance tool for the technical computing as Matlab/Simulink.
As part of the scientific methodology of this Thesis, the developed control strategies are
rigorously validated in an experimental platform used as technological support. This method
enables a mutual advantage of the acquired experience and a results improvement. The real
time control of experimental high performance drilling operations corroborates the results
obtained during the simulation stage. Moreover, these experiments give the possibility to
make comparative studies between the control strategies. Also, they allow comparing the
effectiveness of metaheuristics methods employed for optimal tuning of intelligent control
systems. The comparative studies are based on performance indices that characterize the
dynamic behaviour of the control systems. Performance indices related to the controlled
variable (i.e., cutting force) such as the maximum overshoot (i.e., maximum amplitude of
oscillations) and average values of its error with respect to the reference value (i.e., ITAE,
IAE, ITSE, among others).
Comparative studies confirm that the fuzzy control systems optimally tuned by different
methods, provide a satisfactory performance of the process. The dynamic transient responses
of control systems obtained in simulation and experimentally, are very similar. The cause is
that in all cases the adjusted controllers parameters are optimal or sub-optimal. But despite
their similarity, the overall performance of the control system optimized by the cross-entropy
method is slightly superior to the others. It also offers a good relation between the
computational time and total setting parameters that requires.
Moreover, comparative studies have confirmed that the nonlinear control strategies based
on soft-computing techniques, in general, have a better performance than linear strategies for
real-time control of high performance drilling operations. Therefore, are more effective to
deal with the inherent complexity of these types of processes.
During the development stages of this Thesis, some of the partial achieved results have
been submitted for consideration by the international scientific community. Researchers from
different fields such as computer science, artificial intelligence, automation, manufacturing
engineering, mechanical engineering and metallurgical engineering, among others, have
positively evaluated some contributions of the reported work.
X
TABLA DE CONTENIDO
INTRODUCCIÓN ............................................................................................................................. 1
Objetivos ....................................................................................................................................... 3
Metodología................................................................................................................................... 5
Estructura del trabajo ..................................................................................................................... 5
XI
2.3.5 Aplicación de los modelos a la predicción de las fuerzas de corte en materiales
de mayor dureza. Análisis de errores .......................................................................66
2.4 Diseño e implementación de un modelo dinámico múltiple entrada-simple salida,
computacionalmente eficiente de un proceso de fabricación complejo. Estudio del caso
particular del proceso de fresado de alta velocidad ..............................................................68
2.4.1 Diseño de experimento para la identificación del modelo dinámico del FAV ...........69
2.4.2 Identificación experimental del modelo dinámico del proceso de fresado de alta
velocidad ................................................................................................................71
2.4.3 Implementación del modelo en Matlab/Simulink .....................................................79
2.5 Conclusiones.......................................................................................................................82
XII
4.5.1 Algoritmo de entropía cruzada adaptado para la solución de problemas de
optimización..........................................................................................................143
4.6 Conclusiones.....................................................................................................................145
CONCLUSIONES ..........................................................................................................................199
Lista de publicaciones .................................................................................................................207
Trabajo futuro.............................................................................................................................208
BIBLIOGRAFÍA ............................................................................................................................211
XIII
LISTA DE FIGURAS
XV
Figura 2.17. Fuerza de corte resultante experimental (línea continua) y obtenida por el modelo
para diferentes condiciones de corte (línea discontinua).................................................53
Figura 2.18. Esquema del proceso de taladrado de alto rendimiento...................................................55
Figura 2.19. Parámetros geométricos utilizados en el modelo de Elhachimi et ál. ..............................57
Figura 2.20. Diagrama de flujo de la aproximación propuesta por Pirtini y Lazoglu...........................58
Figura 2.21. Diagrama de bloques teórico del proceso de taladrado de alto rendimiento.....................60
Figura 2.22. Perfil para la operación de taladrado. .............................................................................61
Figura 2.23. Señal de entrada al proceso (velocidad de avance) y fuerza obtenida, bruta y
media (filtrada) de una operación de taladrado en el material GGG40............................62
Figura 2.24. Implementación de los modelos en Matlab/Simulink. ....................................................63
Figura 2.25. Respuesta de los modelos de primer, segundo y tercer orden de la fuerza ante
cambios en la velocidad de avance y conjunto de datos de validación. ...........................64
Figura 2.26. Función de auto-correlación del error en la predicción de fuerza en el taladrado de
GGG40. ........................................................................................................................66
Figura 2.27. Aplicación de los modelos del GGG40 con ganancia unidad a la fuerza
normalizada del taladrado del material 17-4PH..............................................................67
Figura 2.28. Esquema del proceso de FAV y diagrama de bloques de una representación del
proceso mediante funciones de transferencia. ................................................................69
Figura 2.29. Probeta con perfil escalonado en su superficie con cambios de 0.5 mm. .........................69
Figura 2.30. Diseño de experimento para la identificación del proceso de FAV. Cambios
planificados para cada una de las variables de entrada al proceso...................................70
Figura 2.31. Señales registradas durante la operación de fresado con la herramienta de 12 mm..........72
Figura 2.32. Fuerza media en cada uno de los experimentos. .............................................................73
Figura 2.33. Funciones de ajuste para las ganancias de los sub-sistemas que componen el
modelo..........................................................................................................................75
Figura 2.34. Desarrollo del modelo dinámico del proceso de FAV en Matlab/Simulink. ....................80
Figura 2.35. Simulación del modelo y comparación con los datos experimentales..............................81
Figura 3.1. Diagrama general de un controlador PID en un sistema de control realimentado. .............86
Figura 3.2. Diagrama en Matlab/Simulink del controlador PID discreto.............................................88
Figura 3.3. Sub-conjunto clásico y borroso........................................................................................92
Figura 3.4. Conjunto borroso puntual o singularidad borrosa.............................................................93
Figura 3.5. Diagramas de Venn. ........................................................................................................95
Figura 3.6. Variable lingüística fuerza...............................................................................................97
Figura 3.7. Diagrama en bloques del lazo de control con DLB...........................................................99
Figura 3.8. Representación en capas de la red neuronal adaptativa...................................................115
Figura 3.9. Arquitectura del sistema ANFIS. ...................................................................................120
Figura 3.10. Ejemplo de función de distribución de tipo campana generalizada. ..............................121
Figura 4.1. Diagrama general de un sistema de control de un proceso físico. ...................................126
XVI
Figura 4.2. Diagrama de bloques del sistema de control y algoritmo de optimización.......................129
Figura 4.3. Diagrama de flujo del algoritmo de Nelder-Mead. .........................................................132
Figura 4.4. Algoritmo de Metropolis de aceptación-rechazo. ...........................................................136
Figura 4.5. Diagrama de flujo de datos del algoritmo.......................................................................137
Figura 5.1. Diagrama de simulación del sistema de control de la fuerza de corte..............................149
Figura 5.2. Gráfica del lugar de las raíces para una respuesta aperiódica. .........................................150
Figura 5.3. Simulación del sistema de control PID a lazo cerrado. ...................................................151
Figura 5.4. Diagrama del controlador PID implementado. ...............................................................152
Figura 5.5. Comportamiento en tiempo real del sistema de control. .................................................153
Figura 5.6. Estrategia CMI y modelos ANFIS. ................................................................................154
Figura 5.7. Entradas y salidas de los modelos ANFIS. .....................................................................155
Figura 5.8. Funciones de pertenencia y reglas de los modelos ANFIS..............................................157
Figura 5.9. Respuesta de los modelos ante datos de validación. .......................................................158
Figura 5.10. Aplicación de simulación de la estrategia de control por modelo interno y
modelos neuro-borrosos. .............................................................................................160
Figura 5.11. Simulación de la respuesta del sistema de control ante un escalón en la referencia .......161
Figura 5.12. Comportamiento del criterio ITAE en presencia de retardos de tiempo. ........................161
Figura 5.13. Comportamiento del sistema de control CMI-ANFIS y otras estrategias de control ......162
Figura 5.14. Comportamiento del error relativo para las estrategias de control.................................164
Figura 5.15. Particiones borrosas y funciones de pertenencia de las entradas y salidas. ....................166
Figura 5.16. Característica estática del controlador borroso. ............................................................167
Figura 5.17. Sistema de control distribuido basado en un regulador borroso para el control del
proceso de TAR ..........................................................................................................169
Figura 5.18. Aplicación de simulación del sistema de control borroso..............................................170
Figura 5.19. Resultados de la simulación del sistema de control. .....................................................171
Figura 5.20. Comportamiento del sistema de control según varía el tiempo de retardo. ....................172
Figura 5.21. Comportamiento de los índices de desempeño ITAE e IT2SE según varía el tiempo
de retardo. ...................................................................................................................173
Figura 5.22. Comportamiento del proceso según las estrategias de control PID y borrosas...............174
Figura 5.23. Comportamiento del error absoluto ponderado en el tiempo de las estrategias de
control.........................................................................................................................176
Figura 5.24. Optimización mediante temple simulado......................................................................178
Figura 5.25. Resultados tras 75 ejecuciones del algoritmo de SA. ....................................................179
Figura 5.26. Simulación del sistema de control borroso optimizado mediante SA. ...........................180
Figura 5.27. Comportamiento de la fuerza de corte según la variación del tiempo de retardo............181
Figura 5.28. Comportamiento del criterio ITAE según varía el tiempo de retardo. ............................181
XVII
Figura 5.29. Respuesta temporal del sistema de control optimizado mediante SA durante el
taladrado del material GGG40.....................................................................................182
Figura 5.30. Comportamiento del sistema de control optimizado durante una operación de
taladrado del material 17-4PH .....................................................................................183
Figura 5.31. Diagrama del sistema de control borroso y optimización mediante algoritmos
genéticos.....................................................................................................................184
Figura 5.32. Resultados tras 99 ejecuciones de la optimización mediante GA ..................................185
Figura 5.33. Simulación del sistema de control optimizado mediante GA. .......................................186
Figura 5.34. Aplicación de simulación del sistema de control con ruido aditivo. ..............................186
Figura 5.35. Diagrama del sistema de control borroso y el método de optimización de entropía
cruzada. ......................................................................................................................188
Figura 5.36. Funciones de densidad probabilísticas para los parámetros [ Ke Kce ] ..........................189
Figura 5.37. Resultados durante la optimización mediante entropía cruzada. ...................................190
Figura 5.38. Simulación del sistema de control borroso. ..................................................................190
Figura 5.39. Comportamiento de los criterios de desempeño ITAE e ISTSE en presencia de
retardos de tiempo en el sistema de control..................................................................191
Figura 5.40. Comportamiento del sistema de control durante el taladrado del material GGG40........192
Figura 5.41. Estudio comparativo sin presencia de perturbaciones. ..................................................194
Figura 5.42. Estudio comparativo en presencia de perturbaciones....................................................196
XVIII
LISTA DE TABLAS
Tabla 1.1. Aplicación del soft-computing al modelado, monitorización y control de procesos. ...........18
Tabla 1.2. Técnicas de IA en la optimización de procesos de mecanizado..........................................24
Tabla 1.3. Sistemas inteligentes en la automatización de procesos de mecanizado. ............................25
Tabla 2.1. Comparación de las fuerzas de corte obtenidas experimentalmente y mediante el
modelo. .........................................................................................................................54
Tabla 2.2. Resultados del análisis estadístico de los retardos de tiempo del bus de campo MPI. .........59
Tabla 2.3. Propiedades de los materiales de las piezas de trabajo experimentales. ..............................60
Tabla 2.4. Condiciones nominales para la operación de taladrado. .....................................................61
Tabla 2.5. Modelos del taladrado de GGG40. ....................................................................................62
Tabla 2.6. Resultados del análisis de error para los modelos del taladrado de GGG40........................65
Tabla 2.7. Análisis del error de ajuste para los modelos generados del material GGG40
empleados en el material 17-4PH...................................................................................67
Tabla 2.8. Valores que alcanzan cada una de las variables de entrada a lo largo del eje x de la
pieza, durante la operación de ranurado. ........................................................................71
Tabla 2.9. Condiciones nominales para el experimento de identificación. ..........................................71
Tabla 2.10. Parámetros de algunos de los modelos obtenidos con la herramienta de 12 mm. ..............74
Tabla 2.11. Funciones de ajuste de las ganancias de los sub-sistemas para la herramienta de 10
mm. ..............................................................................................................................76
Tabla 2.12. Análisis estadísticos de las constantes de tiempo de los modelos de los sub-
sistemas. .......................................................................................................................76
Tabla 2.13. Valores de ajuste de los coeficientes para el cálculo de las ganancias del modelo. ...........78
Tabla 2.14. Valores de ajuste de las constantes de tiempo y retardos de tiempo de los sub-
sistemas que conforman el modelo del proceso de FAV.................................................78
Tabla 2.15. Errores del modelo, réplicas de ajuste y de validación (filas sombreadas). .......................82
Tabla 3.1. Métodos de ajuste de algoritmos de control PID................................................................88
Tabla 3.2. Algunas funciones de pertenencia típicas. .......................................................................104
Tabla 4.1. Algunos tipos de representación de los individuos de una población................................138
XIX
Tabla 4.2. Ejemplos de operadores de cruce. ...................................................................................139
Tabla 5.1. Criterios de error de la respuesta temporal de la fuerza de corte. .....................................153
Tabla 5.2. Resumen de los parámetros de los modelos neuro-borrosos.............................................159
Tabla 5.3. Estudio comparativo entre el controlador CMI-ANFIS y otras estrategias de control
analizadas. ..................................................................................................................163
Tabla 5.4. Base de reglas del controlador borroso............................................................................167
Tabla 5.5. Comportamientos extremos del sistema de control para el tiempo de retardo
variable. ......................................................................................................................173
Tabla 5.6. Estudio comparativo de las estrategias de control analizadas...........................................175
Tabla 5.7. Estudio comparativo de los resultados experimentales y de simulación del
comportamiento del controlador borroso optimizado mediante el método de
entropía cruzada..........................................................................................................192
Tabla 5.8. Parámetros de las estrategias de optimización para un ajuste óptimo de un
controlador borroso de 49 reglas..................................................................................193
Tabla 5.9. Estudio comparativo del desempeño del sistema de control borroso sin presencia de
perturbaciones.............................................................................................................194
Tabla 5.10. Estudio comparativo del desempeño del sistema de control borroso en presencia
de perturbaciones. .......................................................................................................195
XX
INTRODUCCIÓN
Los sistemas físicos involucrados en la industria del siglo XXI forman parte de un sistema
social y técnico dinámico en un entorno cada vez más globalizado y turbulento. Los cambios
que ha sufrido la economía mundial en las últimas dos décadas, debido al crecimiento
económico de nuevos países emergentes, además de fenómenos como la deslocalización, han
puesto en el punto de mira a varios sectores considerados estratégicos, entre ellos un sector
muy importante y sensible para el empleo, el bienestar social y la seguridad de las naciones: el
sector de la fabricación. Los procesos de fabricación tienen, por su variedad y complejidad, un
gran impacto en el producto nacional bruto de Europa. Se estima que la actividad de
fabricación representa aproximadamente el 22% del PNB europeo, además, un 75% del PIB y
un 70% del empleo a nivel europeo están relacionados con la fabricación [1].
Existe la necesidad de desarrollar e implementar estrategias de fabricación basadas en la
investigación y la innovación, que aceleren la transformación de la industria europea,
asegurando un empleo de valor añadido y logrando que una parte importante de la producción
esté soportada en la economía del conocimiento. Esta estrategia debe garantizar una continua
competitividad de la industria de fabricación, significando una radical transformación de la
industria en todos sus niveles, desde la red de sistemas de fabricación hasta las fábricas, los
sistemas de producción, máquinas, componentes y procesos técnicos.
Entre las áreas prioritarias que asegurarán el desarrollo del sector de la fabricación en
Europa en los próximos años, se encuentra la relacionada con la transformación de la
fabricación de componentes basada en recursos (humanos, computacionales, maquinaria, etc.)
en una fabricación basada en conocimiento. Esa transformación va indisolublemente ligada a
las Ciencias de las Computación y a la Inteligencia Artificial de modo que permitan una
materialización pragmática de la fabricación adaptativa, digital y basada en conocimiento. Es
decir, hay que incorporar una respuesta automática a los cambios del entorno, involucrando
software y tecnologías convergentes en el diseño y operación de los procesos, además de la
1
Métodos clásicos y de soft-computing en la optimización de procesos complejos.
2 Introducción
Raúl M. del Toro
Objetivos
El objetivo principal de este trabajo es el diseño, desarrollo y aplicación de estrategias de
optimización basadas en soft-computing, en el contexto de los sistemas de control cognitivos
artificiales, para el control de procesos complejos [2]. Asimismo, se trata de obtener una
representación computacional eficiente y simple, de procesos de elevada complejidad, como
3
Métodos clásicos y de soft-computing en la optimización de procesos complejos.
4 Introducción
Raúl M. del Toro
Metodología
Para el desarrollo de esta Tesis Doctoral se han utilizado métodos de investigación
teóricos y experimentales, interactuando mutuamente. Las estrategias concebidas y los
algoritmos desarrollados son fruto de la aplicación de los más adecuados paradigmas elegidos
dentro del amplio universo que ponen a nuestra disposición las Ciencias de la Computación y
la Inteligencia Artificial, la Teoría de Sistemas y la Automática.
El estudio y selección de las teorías más viables para el modelado, el control y la
optimización ha sido precedido por una re-elaboración de estas estrategias y algoritmos que
ha permitido su concepción como una doctrina integradora en el diseño de sistemas
inteligentes de control óptimos para sistemas electromecánicos complejos.
Las estrategias de control y optimización desarrolladas, han sido objeto de validación
experimental rigurosa mediante la utilización de un soporte tecnológico, potenciando de esta
forma, la sinergia entre la amalgama de técnicas y procedimientos desarrollados. Esto permite
aplicar de manera iterativa, investigación teórica y experimental, pudiendo aprovechar de
forma mutua las experiencias alcanzadas y mejorando así, los resultados obtenidos. De este
modo se contribuye a una mayor completitud del trabajo de Tesis Doctoral y de los resultados
aquí plasmados.
5
Métodos clásicos y de soft-computing en la optimización de procesos complejos.
6 Introducción
Raúl M. del Toro
7
CAPÍTULO 1
ESTADO DEL ARTE
En este capítulo se presenta una breve revisión del estado de la ciencia y la técnica en los
paradigmas elegidos dentro del amplio universo que ponen a nuestra disposición la Teoría
General de Sistemas, las técnicas de Inteligencia Artificial, las Ciencias de la Computación y
la Automática, considerando específicamente la conjunción de estas grandes áreas de
conocimiento. Se presenta, además, una breve revisión de los procesos de mecanizado, como
caso de estudio de procesos electro-mecánicos complejos, presentando cómo han
evolucionado estos procesos y las tendencias actuales hacia un aprovechamiento efectivo de
la complementariedad con otras áreas científicas y tecnológicas.
El capítulo se divide en 3 secciones. En la primera sección (sección 1.1) se muestra una
revisión de las diversas técnicas provenientes, fundamentalmente, de las Ciencias de la
Computación y vinculadas con el campo de la Inteligencia Artificial, con vistas a ser
aplicadas al modelado, control y monitorización de procesos físicos. Como procesos electro-
mecánicos a estudiar y optimizar, en la sección 1.2 se realiza una introducción a los procesos
de mecanizado, presentando su evolución y tendencias actuales. Finalmente, en la sección 1.3,
se recogen varias de las aportaciones disponibles en la literatura, que demuestran la viabilidad
e importancia de la optimización de los procesos de mecanizado mediante la aplicación de
técnicas de Inteligencia Artificial.
9
Métodos clásicos y de soft-computing en la optimización de procesos complejos.
que hay que aprovechar la sinergia entre éstas, los métodos computacionales avanzados y los
avances tecnológicos en general. De este modo se pueden alcanzar diseños de sistemas de
control y/o monitorización cualitativamente superiores y aplicarlos a la mejora de procesos
productivos.
Dentro del amplio espectro de las técnicas de Inteligencia Artificial, una de las ramas con
mayor desarrollo científico-técnico en los últimos años que se nutre, principalmente, de las
sinergias entre las técnicas básicas de la IA es el “Soft Computing”. Esta abarca un amplio
conjunto de métodos para el diseño de sistemas inteligentes, capaces de manejar de forma más
provechosa información incierta, imprecisa y/o incompleta permitiendo, además, abordar
problemas reales para obtener soluciones más robustas, manejables y de menor coste que las
obtenidas mediante técnicas convencionales. El término soft-computing ha sido acuñado por
Lofti Zadeh en una publicación en 1994 [4], definiéndola como (adaptación de la traducción
realizada por Verdegay [5]):
"Básicamente, soft-computing no es un cuerpo homogéneo de conceptos y técnicas.
Más bien es una mezcla de distintos métodos que de una forma u otra cooperan desde
sus fundamentos. En este sentido, el principal objetivo del soft-computing es
aprovechar la tolerancia que conllevan la imprecisión y la incertidumbre, para
conseguir manejabilidad, robustez y soluciones de bajo costo. Los principales
ingredientes del soft-computing son la lógica borrosa, la neuro-computación y el
razonamiento probabilístico, incluyendo además como subconjunto a los algoritmos
genéticos, las redes de creencia, los sistemas caóticos y algunas partes de la teoría de
aprendizaje. En esa asociación de lógica borrosa, neuro-computación y razonamiento
probabilístico, la lógica borrosa se ocupa principalmente de la imprecisión y el
razonamiento aproximado, la neuro-computación del aprendizaje y el razonamiento
probabilístico de la incertidumbre y la propagación de las creencias".
Zadeh define como premisas básicas del soft-computing que la imprecisión y la
incertidumbre son características intrínsecas del mundo real, además que, la precisión y la
certeza conlleva un coste. Expresa, además, como principio rector del soft-computing,
explotar la tolerancia en la imprecisión, la incertidumbre y la verdad a medias, para alcanzar
tratabilidad, robustez y soluciones de bajo coste [4].
Como se ha de extraer del concepto de Zadeh, la asociación de metodologías de cómputo
que conforman el soft-computing, incluye como sus principales miembros [6] la lógica
borrosa (LB), la neuro-computación (NC), la computación evolutiva (CE) y la computación
probabilística (CP) (ver Figura 1.1). A partir de la simbiosis entre estas metodologías es
posible diseñar sistemas inteligentes híbridos para ser aplicados a procesos físicos y de esta
manera hacer frente a la imprecisión e incertidumbre inherente del mundo real.
10 Capítulo 1
Raúl M. del Toro
En un trabajo publicado en 1999, Bonissone et ál. [7] ofrecen una amplia colección de
sistemas inteligentes híbridos para realizar tareas de diagnóstico, estimación y control en
aplicaciones industriales y comerciales, tales como dispositivos lógicos borrosos (DLB)
ajustados mediante redes neuronales (NN, Neural Networks) y algoritmos de computación
evolutiva (CE), redes neuronales ajustadas mediante CE o por DLB, además de algoritmos de
CE controlados mediante dispositivos lógicos borrosos.
Entre los métodos de mayor popularidad que comparten este enfoque se encuentran el
algoritmo de ascenso a colinas (HC, Hill Climbing) [16], temple simulado (SA, Simulated
Annealing) [17, 18], búsqueda tabú (TS, Tabu Search) [19], entre otros. Otro de los
algoritmos que vale la pena mencionar, aunque menos conocido debido a que su publicación
1
La computación evolutiva y los algoritmos que la componen, como es el caso de los algoritmos evolutivos
(AE), comparten el mismo principio que estos métodos y por tanto en la literatura son considerados también de
tipo Monte Carlo (T. Weise, "Global Optimization Algorithms. Theory and Application," http://www.it-
weise.de, 2009).
12 Capítulo 1
Raúl M. del Toro
es más reciente y por ende su aplicación, es el algoritmo de entropía cruzada (CE, Cross
Entropy) [20]. Este algoritmo fue desarrollado con el objetivo de simular eventos raros, con
probabilidades muy pequeñas, pero con requerimientos de una buena precisión en su
estimación.
Existe una amplia variedad de técnicas heurísticas y sub-clases de estas, que han sido
aplicados y pueden emplearse de forma combinada con otros componentes del soft-
computing, como son, además de las técnicas mencionadas con anterioridad, los algoritmos de
búsqueda directa (ej.: algoritmo de optimización de Nelder-Mead (NM) [21]) y producto de la
hibridación de métodos heurísticos (ej.: algoritmos meméticos), por solo mencionar algunos.
La elección para la aplicación de una u otra técnica, dependerá de la eficiencia o desempeño
que demuestre ante cada problema específico. Ofrecer una taxonomía completa de los
algoritmos de búsqueda global, se escapa del ámbito de este trabajo y conllevaría una
dispersión de los objetivos específicos que se persiguen.
La Figura 1.2 muestra un diagrama que pretende abarcar, tanto los miembros principales
del soft-computing, varios de sus componentes básicos, así como sistemas híbridos que han
sido propuestos en la literatura o que pueden surgir de la combinación de los componentes
básicos. El diagrama parte de la clasificación inicial propuesta por Bonissone [7, 22], en la
que se definen varios niveles. En este caso se ha situado en el primer nivel, los cuatro
miembros principales actuales, la computación probabilística, la lógica multievaluada y
computación borrosa, la neuro-computación y tal como proponen Verdegay, Yager y
Bonissone [8], las técnicas metaheurísticas y algoritmos de búsqueda global, como cuarto
componente, ya que abarca la computación evolutiva, tal y como había propuesto Zadeh
inicialmente [6]. A continuación, para cada uno de los componentes básicos se relacionan
algunas de sus técnicas más conocidas:
• Dentro de la computación probabilística se hace énfasis en las redes bayesianas y los
sistemas Dempster-Shafer, para el desarrollo de modelos probabilísticos.
• Para el caso de la computación borrosa y multievaluada, en el diagrama se destacan por
un lado los sistemas algebraicos de lógica multievaluada y los sistemas de inferencia
borrosa, en concreto para este último, se especifican los dispositivos lógicos borrosos
(DLB), empleados para aplicaciones de modelado y control de procesos, de tipo Mamdani
[23] y de tipo Takagi-Sugeno-Kang (TSK) [24, 25].
• Dentro de la neuro-computación, algunas de sus estrategias básicas se han reagrupado de
acuerdo a su topología de red, por un lado se relacionan las redes neuronales (NN, Neural
Networks) artificiales anticipativas (FNN, Feedforward NN) y las redes neuronales
recurrentes (RNN, Recurrent NN). Dentro de la categoría de redes FNN se detallan el
perceptrón de capa simple (NN SL, NN Single Layer), el perceptrón multicapa (MLP,
Multiple Layer Perceptron), las redes neuronales de base radial (RBF, Radial Basis
Function) [26], las redes conocidas como controlador del modelo aritmético del cerebelo
(CMAC, Cerebellar Model Arithmetic Controller) [27] y con un desarrollo un poco más
reciente, las redes neuronales aumentadas (Augmented FNN) [28], las redes neuronales
dinámicas aditivas de tiempo continuo (CA–DNN, Continuous-time Additive Dynamic
NN) [29, 30], entre otras. Relacionado con la categoría de redes RNN, estas han sido, por
lo general, orientadas hacia el aprendizaje no-supervisado, como memorias asociativas y
estructuras auto-organizadas. Algunas de las estructuras más conocidas son las redes
competitivas, los mapas auto-organizados de Kohonen (Kohonen SOM, Self-Organizing
Map) [31], las redes de Hopfield [32], los modelos de red basados en la teoría de la
resonancia adaptativa (ART, Adaptive Resonance Theory) [33], las redes de
almacenamiento a corto y largo plazo (LSTM, Long Short Term Memory) propuestas por
Hochreiter y Schmidhuber [34] y las redes neuronales de conexión estocástica (SCN,
Stochastic Connection Neural networks) [35], entre muchas otras propuestas.
• Relacionado con las técnicas de búsqueda global y metaheurísticas, las diferentes
estrategias se han agrupado según las sub-clases de algoritmos mencionados con
anterioridad, la inteligencia de enjambre (SI), los algoritmos de simulación de procesos
físicos (SPF), por supuesto los algoritmos evolutivos (AE), que son considerados
miembros del soft-computing desde su definición inicial y en este caso representados por
los algoritmos genéticos (AG) [36], la programación genética (PG) [37], la programación
evolutiva (PE) [38] y las estrategias de evolución (ES, Evolution Strategies) [39],
resaltando además en el diagrama, el algoritmo de búsqueda directa de Nelder-Mead
(NM).
En un segundo nivel del diagrama, se muestran algunos sistemas híbridos fruto de la
combinación de las estrategias básicas. Por ejemplo, de la combinación de las técnicas
probabilísticas con las estrategias borrosas, se han desarrollado los denominados modelos
probabilísticos híbridos, fundamentalmente la probabilidad de eventos borrosos [40] y la
veracidad de eventos borrosos [41], entre otros.
Por otro lado, la interacción de las estrategias de neuro-computación con las metodologías
de la lógica borrosa y las técnicas de optimización, ha llevado a la elaboración de propuestas
como los llamados sistemas neuronales híbridos, particularmente la generación formal de su
estructura y los valores de los pesos en cada nodo, además del control borroso de los
parámetros de las redes neuronales artificiales. Por ejemplo, en el caso de la aplicación de los
algoritmos evolutivos, algunas propuestas han sido enfocadas hacia el entrenamiento
optimizado de las redes, como alternativa a los métodos clásicos de entrenamiento (ej.: retro-
propagación de errores) [42]. Además, las propuestas han tratado el tema de la evolución de la
topología de la red, de elementos tales como sus capas ocultas, número de nodos por capas y
enlaces entre ellos [43, 44], resaltando la idea de Maniezzo para evolucionar simultáneamente
los pesos y la topología de la red [45]. En cambio, algunas aplicaciones para el control
14 Capítulo 1
Raúl M. del Toro
borroso de las redes, han sido enfocadas a la monitorización del desempeño de las redes, al
igual que a modificar los parámetros de control de las estrategias de aprendizaje. Ejemplos
típicos han sido las propuestas de Jacobs [46] y Arabshahi et ál. [47], en las que se realiza un
aprendizaje adaptativo mediante una ley de control borrosa, que varía la tasa de aprendizaje a
partir del error obtenido durante el entrenamiento.
Más específicamente, para el caso de la idea de utilizar la lógica borrosa para mejorar el
desempeño de los algoritmos evolutivos, las propuestas de hibridación han sido más recientes,
desembocando principalmente en tres tipos de estrategias: los sistemas genéticos borrosos, los
sistemas bioinspirados y el control borroso de los parámetros de algoritmos evolutivos. Las
proposiciones realizadas van encaminadas principalmente a mejorar el procedimiento de
ajuste de los parámetros que conforman cada una de las etapas de estos algoritmos, tales como
el tamaño de la población, tasas de cruce, etc. Uno de los trabajos pioneros en este sentido,
fue la proposición de Herrera y Lozano [48] para el ajuste de algoritmos genéticos,
obteniendo una mejora en la eficiencia y la velocidad de convergencia.
Por otra parte, las iniciativas de hibridación entre estrategias borrosas con algunos de los
componentes de la neuro-computación y las estrategias de optimización, al igual que la gran
mayoría de los sistemas inteligentes híbridos, surgen de la necesidad de diseñar sistemas
inteligentes que representen el conocimiento de forma eficiente computacionalmente, además
de reducir costes de diseño. Por ejemplo, relacionado con los DLB de tipo Mamdani, en los
cuales su base de conocimientos está formada por los factores de escala (FE) de sus variables
de entrada y salida, las funciones de pertenencia (FP) también de cada una de las variables y
la base de reglas (BR), se han desarrollado estrategias para el ajuste de estos componentes, ya
sea mediante redes neuronales, algoritmos evolutivos y otras técnicas metaheurísticas. Para
este tipo de estrategias vinculadas con la lógica borrosa, se ha realizado un estudio más
exhaustivo de las propuestas claves, debido a su mayor vínculo con los objetivos específicos
de este trabajo.
En el caso de la hibridación con redes neuronales, gozan de gran popularidad los sistemas
neuro-borrosos (ver Figura 1.3) que son, en esencia, sistemas de inferencia borrosa ajustados
mediante redes neuronales o a través de sus estrategias de aprendizaje. Mediante esta
estrategia la base de reglas se representa como una red neuronal, donde los parámetros de las
funciones de pertenencia constituyen los pesos de las neuronas. Por tanto, al entrenar la red
neuronal se ajusta la base de conocimientos del dispositivo lógico borroso. Algunos ejemplos
de sistemas neuro-borrosos de tipo Mamdani disponibles en la literatura son el sistema
NEFCON [50], desarrollado por Nauck en 1994 para tareas de control y el sistema
NEFCLASS [51], propuesto por Nauck y Kruse para ejecutar labores de clasificación.
Figura 1.4. DLB de tipo TSK de primer orden y arquitectura ANFIS equivalente.
(Jang, 1993 [52])
Por otra parte, en el caso del sistema de inferencia borrosa de tipo TSK, su estructura se ha
equiparado con las redes neuronales de base radial [53, 54]. Este tipo de DLB emplea
funciones polinómicas como salida en cada regla de su base de reglas, incrementando de esta
forma, el potencial de representatividad del conocimiento. Basado en la equivalencia de
ambas estructuras, en 1993 Jang publica el sistema de inferencia borrosa basado en redes
neuronales adaptativas (ANFIS, Adaptive Network-based Fuzzy Inference System) [52] (ver
Figura 1.4), donde la base de reglas determina la topología de la red, representando los nodos
16 Capítulo 1
Raúl M. del Toro
18 Capítulo 1
Raúl M. del Toro
mecanizado, surgiendo lo que conocemos hoy como máquinas de control numérico por
computador o máquinas CNC. Estas máquinas son capaces de repetir secuencias de órdenes
con precisión, pudiendo producir por tanto, piezas de mayor complejidad.
Los procesos de mecanizado se realizan sobre piezas, por lo general tratadas previamente,
con el fin de mejorar sus propiedades mecánicas como la dureza, la resistencia, la elasticidad,
etc. La Figura 1.5 muestra una clasificación de los procesos de mecanizado, agrupados
principalmente por convencionales o no, y luego entre las operaciones tradicionales se hace
diferenciación entre los procesos por arranque de viruta y procesos de conformado. Entre las
operaciones que forman parte de los procesos de mecanizado por arranque de viruta (ver
ejemplos de la Figura 1.6) se encuentran el taladrado, torneado, fresado, roscado, grabado,
entre otras, las cuales son consideradas como operaciones de mecanizado convencionales.
Constituyendo el taladrado, el torneado y el fresado, las operaciones con mayor peso en la
industria.
a) b)
Los procesos de mecanizado tienen una gran importancia dentro de los sectores de
fabricación. Estos procesos involucran maquinaria y materiales muy costosos, así como
operadores expertos. Además, la fuerte competencia a nivel internacional en el sector del
mecanizado obliga a desarrollar e implantar sistemas de monitorización y control que mejoren
los procesos productivos en todas sus facetas, como pueden ser: incremento de las tasas de
arranque del material, disminución de tiempos de parada, mejorar la calidad superficial de los
componentes fabricados así como su calidad dimensional, conservar las herramientas de corte,
además de la durabilidad de las máquinas, entre otras [82]. Sin embargo, la complejidad
intrínseca de estos procesos hace que esta tarea no resulte nada sencilla. La fuerte interacción
en las máquinas del movimiento relativo entre la pieza y la herramienta genera problemas
como por ejemplo, vibraciones y el retemblado, la obtención de formas y de acabados
imprecisos, deformaciones de orígenes mecánicos o térmicos, etc. [83].
20 Capítulo 1
Raúl M. del Toro
a) b)
Por otra parte, el desarrollo de nuevas tecnologías ha traído consigo la implantación en los
últimos años de nuevos procesos de mecanizado, en los cuales el arranque de material se
produce de forma no convencional, empleando energía térmica, mecánica, química, eléctrica
o la combinación entre ellas, denominándose operaciones de mecanizado no tradicionales (ver
ejemplos de operaciones en la Figura 1.7). Estas operaciones suelen dividirse en cuatro
categorías [85]:
1. Químicas: la eliminación de material se produce mediante una reacción química (ej., el
mecanizado fotoquímico conocido también como fresado fotoquímico).
2. Eléctricas o electroquímicas: el material se elimina por disolución electrolítica. Ej.,
mecanizado electroquímico (ECM, Electro-Chemical Machining) y esmerilado
electroquímico.
3. Térmicas: la operación se lleva a cabo concentrando altas temperaturas en zonas muy
localizadas en las que funde y vaporiza el metal. Entre sus operaciones podemos citar el
mecanizado por electro-erosión (EDM, Electrical Discharge Machining) donde la
eliminación de material se produce mediante un arco eléctrico; el mecanizado mediante
haz de electrones (EBM, Electron-Beam Machining) en el cual electrones a muy alta
velocidad y energía inciden sobre la superficie de piezas, creando calor y evaporizando el
material en la región de contacto; el mecanizado con rayo láser y con arco de plasma,
entre otros.
4. Mecánicas: la eliminación de material se produce por erosión o corte con puntas
múltiples. Ejemplo de operaciones, mecanizado ultrasónico, chorro de agua, mecanizado
con chorro de agua abrasivo, entre otros.
22 Capítulo 1
Raúl M. del Toro
procesos, constituyen una herramienta de muy gran utilidad para la continua mejora de estos
procesos.
A lo largo de los últimos años, las aplicaciones de diferentes técnicas de Inteligencia
Artificial (IA) a las operaciones de mecanizado y de fabricación en general, han estado
enfocadas principalmente a su optimización mediante la manipulación de los parámetros que
determinan las operaciones de corte (parámetros de corte), el desarrollo de sistemas de
monitorización, supervisión y control y sistemas de diagnóstico y detección de fallos,
produciendo resultados sorprendentes. De las estrategias desarrolladas, las que se centran en
los parámetros de corte han sido las más estudiadas y las que mejores resultados han aportado.
Según la literatura a la que se ha tenido acceso, tanto las redes neuronales artificiales
(RNA) como los sistemas de inferencia borrosa han sido, de las técnicas de Inteligencia
Artificial, las que con mayor frecuencia se han empleado para la optimización de los procesos
de mecanizado. Debido a su capacidad para la aproximación de funciones, en especial de
relaciones no-lineales, las redes neuronales artificiales han resultado ser una técnica muy
atractiva para el modelado de los procesos de corte [88, 89].
En un trabajo publicado en 1998, Coit et ál. [64] consideraron los aspectos prácticos de la
construcción y validación de modelos neuronales, mientras que Viharos et ál. [90] mostraron
la aplicación de las redes neuronales artificiales al control y monitorización de procesos de
mecanizado. Diversas aplicaciones de modelado de procesos de mecanizado a través de redes
neuronales artificiales basadas en relaciones entradas-salidas han sido reportadas en la
literatura. Ejemplo de ello constituyen las aplicaciones satisfactorias del perceptrón multicapa
(MLP) para modelar el lijado de super-aleaciones propuesta por Sathyanarayanan et ál. [91],
el modelo desarrollado por Jain et ál. [92] para la predicción de la tasa de arranque de material
y el acabado superficial de un mecanizado abrasivo, además del modelo propuesto por Feng et
ál. [65] del tallado de cilindros para motores. Además, en un trabajo desarrollado por Grzesik
y Brol [93], estos autores mostraron la utilidad del modelado y el control neuronal, hibridado
con otras técnicas de análisis (ej., estadísticas), de ciertas características del acabado
superficial en procesos de mecanizado de múltiples etapas. Algunos de los inconvenientes de
las redes neuronales artificiales son que los parámetros del modelado no pueden ser
interpretados, que son dependientes de grandes conjuntos de datos de entrenamiento, que
ciertas configuraciones de las redes producen sobre-entrenamiento y que, en ocasiones, el
entrenamiento culmina sin que se hayan alcanzado soluciones óptimas.
Por otro lado, la lógica borrosa también ha jugado un papel muy importante en el
modelado de procesos de mecanizado. El modelado basado en la teoría de conjuntos borrosos
es usualmente utilizado cuando existe conocimiento subjetivo del proceso u opiniones de
operadores expertos, desempeñando un papel fundamental en la definición de funciones
objetivo y variables de decisión. Shin y Vishnupad [94] observaron que la hibridación de los
Por otra parte, en 1998 W.L.R. Ip publicó una estrategia de control basada en lógica
borrosa para la regulación de la velocidad de avance en el fresado de un acero templado,
obteniendo una mayor eficiencia de corte y un mayor aprovechamiento de la vida útil de la
herramienta [106]. Lee et ál. usaron la teoría de conjuntos borrosos para obtener un modelo
no-lineal del proceso de torneado, mostrando de este modo una herramienta más efectiva que
las técnicas de modelado matemático para el caso en el que exista cierta “vaguedad” en la
variables de control del proceso [107]. Al-Wedyan et ál., ya en el año 2001, usaron técnicas
de modelado borroso para operaciones de fresado [108]. Para una mayor ilustración de las
24 Capítulo 1
Raúl M. del Toro
Behrens y
Control borroso de la distancia entre los
Electro-erosión Ginzel, 2003
electrodos
[111]
Mecanizado Shi, Schillings y
MLP para el modelado del proceso.
químico Boyd, 2004 [112]
Micro-mecanizado Yousef et ál.,
MLP para el modelado del proceso.
con láser 2003 [113]
Pulido mecánico- Lin y Liu, 2003
ANFIS para predecir la rugosidad superficial
químico [114]
Control neuro-borroso del esfuerzo de corte Fresado de alta Haber et ál., 2005
basado en el control por modelo interno. velocidad [77]
Mecanizado
tradicional
26 Capítulo 1
Raúl M. del Toro
1.4 Conclusiones
En este capítulo se ha realizado una revisión del estado de la ciencia y la técnica en los
paradigmas elegidos dentro del amplio universo que ponen a nuestra disposición la Teoría
General de Sistemas, la Teoría de Control, las técnicas de Inteligencia Artificial y las Ciencias
de la Computación. Se ha prestado especial interés en los paradigmas con una
fundamentación teórica sólida en el campo de la Inteligencia Artificial (IA), con un enfoque
desde la perspectiva de la Automática y la Teoría de Sistemas. El estudio se ha centrado en
estrategias susceptibles de ser utilizadas en el modelado, monitorización y control de procesos
físicos. Paradigmas que aprovechan tanto las ventajas de las técnicas básicas de Inteligencia
Artificial, como las sinergias entre estas. Por este motivo, la revisión se ha enfocado hacia las
técnicas del soft-computing, analizando su evolución desde sus inicios hasta las tendencias
actuales.
Debido a que las estrategias serán aplicadas al control de procesos electro-mecánicos
complejos, como es el caso de los procesos de mecanizado, se ha realizado además, una breve
revisión del estado de la técnica de estos procesos. Presentándose de forma breve como ha
sido su evolución hasta nuestros días y las tendencias actuales en la aplicación de técnicas de
IA, desde la perspectiva del control, para hacer frente a la complejidad de estos procesos y así
mejorar su eficiencia.
En el siguiente capítulo se describe la aplicación de métodos de investigación tanto
teóricos como experimentales para el modelado, mediante algoritmos computacionales
eficientes, de algunos de estos procesos, como es el caso del fresado de alta velocidad y el
taladrado de alto rendimiento.
Uno de los aspectos claves, al cual la comunidad científica continúa prestándole atención
en la actualidad, es la obtención de modelos que permitan caracterizar los sistemas físicos
presentes en una gran variedad de procesos y con ello poder simular su comportamiento [122-
126].
Los sistemas de supervisión y control de procesos complejos, desempeñan un papel cada
vez más relevante para alcanzar un cambio cualitativo en los procesos. Para ello, es necesario
disponer de métodos de cómputo eficientes para la predicción ó estimación con cierto grado
de precisión, el posible comportamiento temporal de sus variables de salida, a partir de sus
variables de entrada y sus parámetros. La monitorización del comportamiento dinámico de
estas variables permitiría detectar fallos en el proceso, desgaste de los elementos actuadores,
predecir la presencia de perturbaciones no deseadas y aportar información acerca de la
eficiencia y calidad de los procesos. Para ello resulta clave la obtención de un modelo que
permita la simulación de procesos físicos y la predicción del comportamiento sus variables
más significativas.
Una de las herramientas más ampliamente utilizadas dentro del campo de las ingenierías,
para el análisis y modelado de sistemas y procesos físicos, es Matlab®. Esta herramienta es un
lenguaje de alto rendimiento para la computación técnica, que integra el cómputo, la
visualización y la programación en un ambiente fácil de uso donde los problemas y las
soluciones se expresan en una notación matemática familiar. En el entorno se incluyen
facilidades para el análisis numérico, cálculo matricial, procesamiento de señales,
visualización gráfica y desarrollo de aplicaciones, incluyendo el diseño de interfaz gráfica de
29
Métodos clásicos y de soft-computing en la optimización de procesos complejos.
30 Capítulo 2
Raúl M. del Toro
u (t ) PROCESO y(t )
y ( s ) b0 s m + b1s m−1 + ⋯ + bm − sτ
G( s) = = n ⋅e (2.1)
u (s ) s + a1s n−1 + ⋯ + an
1 d mu (t − τ ) d m −1u (t − τ )
y (t ) = b0 + b + ⋯ + b u ( t − τ ) −
dt m−1
1 m
an dt m
(2.2)
dy (t ) d n−1 y (t ) d n y (t )
n−1
a + ⋯ + a +
n −1
dt n
1
dt dt
Los parámetros del modelo θˆ pueden estimarse a partir de la ecuación (2.5). Una solución
directa para el sistema de ecuaciones usualmente no es posible, dado que típicamente el
número de ecuaciones, determinado por el total de muestras P , es mayor que el total de
parámetros a determinar p (sistema sobredeterminado).
−1
θˆ = XT X XT Y (2.5)
32 Capítulo 2
Raúl M. del Toro
e(t ) = y (t ) − yˆ (t ) (2.6)
con una característica de fase arbitraria y respuesta a impulso h[k ] , obteniéndose la secuencia
de datos g[k ] . Seguidamente la nueva secuencia es invertida ( g[−k ] ), aplicándose
nuevamente el mismo filtro y dando lugar a la secuencia de datos r[k ] . Por último, esta
secuencia se invierte y se obtiene s[k ] , que constituye la secuencia filtrada pero sin cambio de
fase con respecto a x[k ] .
Por otra parte, un algoritmo de filtrado muy sencillo pero de gran utilidad para el
procesamiento de señales periódicas, es el filtro de respuesta a impulso finita (FIR) de media
móvil. El principio de este filtro reside en hallar la media de las muestras de una secuencia
dada, utilizando un total de muestras igual a la ventana de datos ws . La respuesta a impulso
del filtro se describe en la ecuación (2.9).
1
h[k ] = 1 + k −1 + k −2 + ⋯ + k − ( ws −1) (2.9)
ws
La aplicación de este filtro a una señal periódica permitiría obtener el valor medio de la
señal o su componente de directa o de muy bajas frecuencias. Entre los procesos que se
encuentran caracterizados por este tipo de señal, se encuentran los procesos de fabricación,
donde señales de gran importancia como las fuerzas de corte, poseen componentes periódicas
debido a la influencia de los accionamientos presentes en estos procesos.
34 Capítulo 2
Raúl M. del Toro
El CNC abierto de la máquina, registra en tiempo real todas las señales de los motores de
los ejes de la máquina, como sus velocidades, posiciones, aceleraciones, par, etc., controlando
a su vez variables del proceso como la velocidad de giro del cabezal, la velocidad de avance
de la máquina, etc. Toda la información que maneja el CNC, se encuentra disponible a través
de un enlace de red bi-direccional denominado MPI (multi-point interface), con una velocidad
de transmisión de datos de 187.5 kbit/s, muy similar al bus de campo PROFIBUS de un uso
muy amplio en la industria en general. A través del bus MPI, también es posible cambiar los
valores de referencia de las variables controladas por el CNC, como la velocidad de giro del
cabezal y la velocidad en la dirección de avance de la máquina. Esto puede realizarse desde
un ordenador remoto o mediante un sistema de control distribuido en red basado en
comunicaciones Ethernet.
Además de las señales de fuerza y de las variables del CNC, se encuentran instalados
sensores para monitorizar otras variables del proceso como la aceleración de vibraciones y
emisión acústica. Estas señales son obtenidas mediante un sistema de adquisición de datos
(DAQ) y pre-procesadas en una aplicación desarrollada en la interfaz LabView. Tanto las
señales de fuerza como la aceleración son adquiridas hasta una frecuencia de muestreo de 50
kHz y la emisión acústica a una frecuencia de 1 MHz, mientras que las señales del CNC son
muestreadas a la misma frecuencia de giro de la herramienta.
Como se puede intuir de la descripción anterior, la plataforma experimental reproduce un
sistema de monitorización y control distribuido real pero a escala de laboratorio. La Figura
2.5 muestra una representación gráfica de todo el sistema.
36 Capítulo 2
Raúl M. del Toro
La existencia de diversos niveles de red, lleva asociado una serie de retardos de tiempo a
la hora de realizar tareas de control o monitorización, que se deben tener en cuenta en la etapa
de diseño de las diferentes estrategias. Estos retardos se adicionan a su vez al propio tiempo
muerto intrínseco en los accionamientos electromecánicos que intervienen en las operaciones
de mecanizado. Los retardos inducidos por la red usualmente no son conocidos a priori y, en
general, son de difícil medición o estimación en línea. Un tratamiento formal acerca de los
retardos inducidos por la red y el procedimiento para estimarlos puede ser encontrado en la
literatura [131, 132].
El acceso al proceso desde el primer nivel de red implica un máximo retardo τ P1 que
incluye el retardo inducido por la red de este nivel τ CA , el retardo del sistema de adquisición
τ SC y el tiempo muerto τ P del proceso (2.10). De la misma forma, el acceder al proceso
desde el segundo nivel de red, implicará un retardo adicional máximo τ12 determinado por el
tiempo de transmisión entre ambos ordenadores τ S 2 y τ A2 (2.11).
τ P1 = [τ SC + τ CA + τ P ]MAX (2.10)
38 Capítulo 2
Raúl M. del Toro
desarrollo del modelo en Matlab considerando solamente los tipos de fresa cilíndricas y de
punta esféricas [135], también en otras plataformas de desarrollo como LabView pero solo
para fresas de punta semiesférica [138]. El proceso de corte ha sido simulado, para la fresa de
punta semiesférica, en lenguajes de alto nivel como C++, interactuando con herramientas para
el modelado en 3D [139]. Además, el modelo para las fuerzas de corte en el mecanizado
convencional fue introducido en trabajos previos, considerando la geometría generalizada de
la herramienta [140] para los tipos de fresas helicoidales más utilizadas en la industria.
Asimismo, para la fresa de punta semiesférica se han desarrollado modelos teniendo en
cuenta, para el cálculo del espesor de la viruta, la influencia de la dinámica del proceso de
corte [141] y considerando el movimiento de la herramienta en superficies 3D [141, 142].
El modelo aquí presentado, cuya implementación se ha realizado en una aplicación en
Matlab/Simulink, a diferencia de trabajos previos incorpora la posibilidad de predecir las
fuerzas de corte para diferentes tipos de herramientas, caracterizadas solamente por
parámetros de su geometría. Además permite evaluar la influencia de la dinámica del proceso
de corte en las fuerzas y el par desarrollado, producto de las vibraciones de la herramienta, y
simular su comportamiento temporal considerando operaciones de mecanizado en superficies
3D.
El poder disponer de un método de cómputo eficiente para la predicción o estimación que,
a partir de las condiciones de corte prefijadas por el tipo de material a mecanizar y la
geometría de la herramienta empleada, sea capaz de predecir con cierto grado de precisión el
posible comportamiento de las fuerzas y el par de corte, resulta de gran importancia. El
comportamiento dinámico de dichas variables permitiría detectar posibles roturas en la
herramienta de corte, estimar su desgaste, predecir la presencia de vibraciones y aportar
información acerca del acabado superficial. Para ello resulta clave la obtención de un modelo
que permita la simulación del proceso y la predicción del comportamiento de estas variables.
para que se produzca el arranque de una unidad de sección de material [143]. La ecuación
(2.12) es una representación aproximada y generalizada para el cálculo de las componentes de
fuerza que se producen. El área o sección de material ( A [mm 2 ] ) que se desprende en forma
de viruta, es proporcional a la relación entre la velocidad de avance de la herramienta con
respecto a la pieza ( f [mm/min] ) y su velocidad de giro ( n [1/min o rev/min] ), es decir,
A∝ f /n.
dF = K e dS + K c A (2.12)
Por otra parte, la potencia específica o potencia requerida para que se produzca el arranque
de una unidad de sección de material, depende de la fuerza y de la velocidad de corte o
velocidad tangencial en el punto de corte, dP = dF Vc .
Figura 2.8. Diagrama del proceso de FAV visto como una caja negra.
(Haber et ál., 2004 [135])
40 Capítulo 2
Raúl M. del Toro
dFr, dFt, dFa Diferencial de fuerzas en los sentidos radial, tangencial y axial.
dz Altura diferencial de un segmento de viruta.
dS Sección infinitesimal del filo de corte de la herramienta.
El modelo matemático impuesto para el borde de corte, considera que el filo se divide en
pequeños incrementos donde los coeficientes de corte pueden variar para cada localización.
Se considera como punto de referencia inicial al primer filo de corte de la herramienta ( j = 1 ),
y su ángulo de rotación cuando z = 0 es φ . El ángulo de inmersión radial para el filo de corte j
en cierta posición axial z queda expresado como:
φ j ( z ) = φ + ( j − 1)φ p −ψ ( z ) (2.13)
42 Capítulo 2
Raúl M. del Toro
Las distancias radiales y axiales desde los puntos M y N hasta el eje de la herramienta y el
extremo de ésta, se definen a partir de las relaciones geométricas de los parámetros generales
tal y como se puede observar en la Figura 2.10. Asumiendo que la hélice es continua, es decir,
las zonas son tangentes entre sí en su intersección, las distancias para cada punto quedan
definidas de la siguiente forma:
N r = Rr + R cos β
(2.15)
N z = Rz − R sin β
Las ecuaciones anteriores son válidas sólo en el caso de que 0 ≤ α < π / 2 y β < π / 2 y su
planteamiento es más sencillo de implementar que el introducido en la referencia [140],
evitando además una posible indefinición de M r y N z para el caso de que R = 0 . Esta
reformulación constituye una contribución al modelado de este proceso.
La distancia radial r (z ) , el ángulo de inmersión axial κ (z ) , el ángulo radial de retraso
ψ (z ) , el ángulo de hélice i ( z ) y la longitud infinitesimal dS del filo en el punto de corte P,
se definen en [140], según la distancia axial z del punto para cada una de las zonas de la
geometría de la herramienta.
hr (φ j , z ) = hr _ st (φ j , z ) + hr _ din (φ j , z ) (2.16)
dr ( z )
hr _ st (φ j , z ) = stj sin Φ + stj sin φ j (t , z ) cos Φ (2.19)
dz
44 Capítulo 2
Raúl M. del Toro
Nf z2
∑ ∫
dS
Fxyz (φ ) = [ A(φ j , κ )K e + A(φ j , κ )K c hr (φ j , z ) ] dz (2.20)
j =1 dz
z1
hr (φ j , z ) K e = [ K re K ae ]
T
donde es el espesor radial de la viruta, K te y
K c = [ K rc K ac ] son los coeficientes de corte, los cuales pueden ser determinados de
T
K tc
forma experimental para un tipo dado de herramienta y material y A (φ j , κ ) es la matriz de
transformación geométrica, la cual se define como,
φ = Ω ⋅t (2.22)
A partir del diferencial de fuerza tangencial dFt se puede determinar el par de corte,
quedando expresado de la siguiente forma:
Nf N f z2
T (φ ) = ∑ T (φ j ) = ∑ ∫ r ( z )dFt (φ j , z ) =
j =1 j =1 z1
N f z2
(2.24)
dS
∑ ∫ r ( z )[ K
j =1 z1
te
dz
+ K tc hr (φ j , z )]dz
46 Capítulo 2
Raúl M. del Toro
Figura 2.12. Diagrama de bloques del modelo matemático de las fuerzas y el par de corte en SIMULINK.
El diagrama de la Figura 2.12 contiene también las entradas y salidas del modelo, las
cuales se describen a continuación:
Entradas feed Avance de la máquina [mm/min].
PHI Ángulo de elevación de la superficie que se va a
mecanizar en cada instante de tiempo.
time Tiempo de la simulación.
cutt_data Datos de corte.
ap Profundidad axial de corte [mm].
ae Profundidad radial de corte [mm].
milling_type Especifica el tipo de operación de fresado, ya sea
contorneado o planeado.
e Excentricidad entre la herramienta y la pieza [mm].
spin_speed Velocidad del cabezal [rpm].
Kc, Ke Vectores de las constantes de corte del material y la
herramienta.
Salidas Fxyz Vector de fuerzas cartesianas.
Figura 2.13. Interfaz de usuario del modelo en SIMULINK para la entrada de parámetros de la geometría
de la herramienta.
48 Capítulo 2
Raúl M. del Toro
Figura 2.14. Árbol de dependencias entre los diagramas o subsistemas que conforman el modelo
desarrollado.
A continuación se relacionan cada uno de los subsistemas con una breve descripción de su
funcionamiento:
50 Capítulo 2
Raúl M. del Toro
a)
b)
Figura 2.15. Otros diagramas de bloques que componen el modelo implementado en SIMULINK.
a) Diagrama que implementa el cálculo de los diferenciales de fuerzas cartesianas y el diferencial de par
de corte, b) diagrama de bloques para el cálculo del espesor radial estático de la viruta.
El modelo permite también, según los parámetros introducidos, la evaluación visual del
comportamiento teórico de las fuerzas de corte a través de gráficas, tanto resultante (Figura
2.16) como en cada eje cartesiano y el par resultante.
Figura 2.16. Gráfica de fuerza de corte resultante obtenida del modelo implementado en SIMULINK para
ciertos parámetros de entrada.
52 Capítulo 2
Raúl M. del Toro
b) Ktc = 496.81 N/mm2, Krc = 364.19 N/mm2, Kac = 54.08 N/ mm2, Kte= 9.84 N/mm, Kre =
12.41 N/mm, Kae = 2.08 N/mm.
c) Ktc = 593.16 N/mm2, Krc = 434.81 N/mm2, Kac = 64.57 N/ mm2, Kte= 11.75 N/mm, Kre
= 14.82 N/mm, Kae = 2.48 N/mm.
d) Ktc = 707.33 N/mm2, Krc = 518.51 N/mm2, Kac = 77 N/ mm2, Kte= 14.01 N/mm, Kre =
17.67 N/mm, Kae = 2.96 N/mm.
a) b)
c) d)
Figura 2.17. Fuerza de corte resultante experimental (línea continua) y obtenida por el modelo para
diferentes condiciones de corte (línea discontinua).
Tabla 2.1. Comparación de las fuerzas de corte obtenidas experimentalmente y mediante el modelo.
( Fe − Ft )100
Experimento Fe [N] Ft [N] AE [%] =
Fe
a 319.38 324.63 1.64
b 158.84 173.13 8.99
c 175.15 171.08 2.32
d 386.02 403.81 4.61
Para evaluar el desempeño del modelo, se utilizaron varias cifras de mérito o índices de
comportamiento para ambos experimentos. De esta manera se pueden comparar las fuerzas de
corte obtenidas a través del modelo y las experimentales. En la Tabla 2.1 se muestra para cada
experimento la fuerza resultante media obtenida experimentalmente, Fe , y la estimada
teóricamente por el modelo, Ft , para estos dos valores se determinó el error absoluto relativo
(AE).
Los coeficientes o constantes de corte (Ktc , Krc, Kte, Kre, Kac, Kae) utilizados en el
modelado fueron obtenidos mediante un procedimiento experimental para los valores de
fuerza media. Ciertamente estos coeficientes tienen un papel determinante en la precisión de
la predicción del esfuerzo de corte. La capacidad de predicción de las fuerzas de corte del
modelo presentado en este trabajo también dependen de coeficientes de corte obtenidos de
forma experimental, y este continúa siendo un punto neurálgico en los modelos matemáticos.
54 Capítulo 2
Raúl M. del Toro
Durante el proceso de TAR (ver Figura 2.18), a medida que se incrementa la profundidad
d en que penetra la herramienta en la pieza de trabajo, a una velocidad de avance f , se
dificulta el proceso de evacuación de la viruta que se produce debido al arranque de material,
lo cual conlleva un incremento de la temperatura en la punta de la herramienta, induciendo
por tanto, una disminución en su vida útil y una demanda de mayor potencia de la máquina.
Este fenómeno puede provocar además una mayor difusión de los elementos químicos de la
herramienta de trabajo hacia la pieza y viceversa [153]. Este fenómeno se agrava cuando se
realiza una operación de taladrado de corte continuo, que no permite una evacuación de la
viruta pero produce un incremento considerable en la tasa de arranque de material. Por tanto,
resulta de gran utilidad el empleo de algoritmos de control sencillos y fiables, basados en
técnicas de IA y capaces de trabajar en sistemas de monitorización y/o control distribuidos,
para reducir el desgaste de la herramienta sin afectar la productividad en el taladrado de alto
rendimiento con corte continuo.
F (s ) K
G(s) = = (2.25)
f (s ) Ts + 1
A partir del resultado anterior, Elhachimi et ál. (ver Figura 2.19) presentaron un modelo
puramente geométrico de las fuerzas de empuje y de corte para brocas, validando el modelo
mediante resultados experimentales de taladrado en piezas de acero [156, 157].
Recientemente se han aplicado las técnicas de IA para el modelado de la operación de
taladrado, mediante el empleo de redes neuronales artificiales para predecir la fuerza existente
durante el comienzo y el final de una operación de taladrado [158], como también la
predicción de la fuerza en una operación de fresado en una arquitectura de control por modelo
interno (IMC) basada en redes neuronales [159]. También, la máquina herramienta para
operaciones de taladrado ha sido evaluada mediante modelos basados en elementos finitos,
para predecir su comportamiento dinámico [160].
56 Capítulo 2
Raúl M. del Toro
Por otra parte, se han realizado estudios de modelado de la operación de taladrado para el
control de par de corte, con vista a reducir el desgaste de herramientas y empleando para ello
modelos dinámicos generalizados de segundo orden [161]. En otro trabajo, Furness et ál. han
modelado el par de corte ante cambios en los valores de referencia de la velocidad de avance
de la máquina, como un sistema de tercer orden con tiempo de retraso fijo [162]. Un modelo
similar de tercer orden pero sin retraso, ha sido implementado junto a un controlador
proporcional, integral y derivativo (PID) para el control del par de corte en operaciones de
taladrado [163].
En otro enfoque diferente para el modelado del proceso de taladrado de alto rendimiento
se obtiene la respuesta dinámica a un impulso del sistema herramienta – cabezal de la
máquina herramienta, mediante un golpe a la herramienta con un martillo calibrado [164] (ver
Figura 2.20). Estos modelos han sido verificados, en una operación de taladrado, analizando
el comportamiento de las fuerzas de corte a medida que cambia la profundidad de corte y
comparándolas con la respuesta prevista por el modelo.
Haciendo caso omiso por completo de los modelos basados en la caracterización física del
proceso de taladrado, Oh et ál. han modelado el proceso como un sistema sub-amortiguado de
segundo orden, ajustando los parámetros del modelo a partir de la respuesta temporal ante un
escalón simple. Al combinar este modelo con un modelo de primer orden del comportamiento
dinámico del cabezal de la máquina [163], se obtiene un modelo resultante de tercer orden:
T (s ) Kc Ks
G p (s ) = = (2.26)
f ( s) (τ fc1 s + 1)(τ fc 2 s + 1)(τ s s + 1)
58 Capítulo 2
Raúl M. del Toro
Tabla 2.2. Resultados del análisis estadístico de los retardos de tiempo del bus de campo MPI.
Retardos de tiempo
Muestras
Desviación
Mínimo Máximo Medio
estándar
10100 9 ms 51 ms 12.58 ms 3.58 ms
Debido a que las acciones de control (señales de referencia de las velocidades de giro y de
avance) son enviadas al CNC a través del bus MPI, los retardos de tiempo que este introduce
deben ser incluidos en el modelo dinámico del proceso. Los retrasos de tiempo del bus MPI,
en concreto el utilizado en la plataforma experimental, han sido estudiados en trabajos
previos. Los resultados del estudio se muestran en la Tabla 2.2, los cuales representan solo los
retrasos introducidos por el bus y no los del proceso físico en sí.
La dinámica de los componentes individuales del sistema no es analizada durante la
identificación, sino que se toma el proceso completo como una sola unidad y se proporciona
una función de transferencia global, que relaciona el cambio en la fuerza de corte medida
Fmed y el cambio en la velocidad de avance de referencia f ref . En la Figura 2.21 se
representa un diagrama de bloques teórico del sistema. El diagrama representa sendos
sistemas de primer orden para el sistema de accionamientos de la máquina, el proceso de
taladrado y el sistema de sensores, con retardos específicos en la dinámica de la planta,
Figura 2.21. Diagrama de bloques teórico del proceso de taladrado de alto rendimiento.
Para generar y validar el modelo del proceso de TAR, se han empleado dos tipos de
material: 17-4PH (un acero inoxidable endurecido por precipitación) y GGG40 (un hierro de
fundición nodular dúctil). Sus propiedades se representan en la Tabla 2.3.
60 Capítulo 2
Raúl M. del Toro
f
f2
f1
d t
d2
d1
t0 t1 t2 t
Los datos de fuerza se filtran mediante un filtro FIR utilizando el algoritmo de media
móvil y evitando introducir retrasos temporales durante el procesamiento (ver sección 2.1.1).
La longitud de la ventana de filtrado ws es igual al número de muestras obtenidas durante
ws = SF ⋅ 60 / n (2.27)
Para cada material, con las mismas condiciones nominales, se realizaron tres réplicas del
experimento, de tal manera que el modelo obtenido a partir de un conjunto de datos pueda
validarse con datos diferentes.
Figura 2.23. Señal de entrada al proceso (velocidad de avance) y fuerza obtenida, bruta y media (filtrada)
de una operación de taladrado en el material GGG40.
La Figura 2.23 muestra las señales de fuerza filtradas y sin filtrar, junto con la velocidad
de avance de una operación de taladro en el material GGG40. Claramente, el filtrado reduce el
ruido de alta frecuencia presente en la señal de fuerza, lo que ayuda a que la creación del
modelo sea más fácil sin pérdida resolución de datos temporales. Además, el retardo en la
salida debido tanto a la arquitectura del medio de conexión como a la propia planta física es
claramente visible en esta figura.
62 Capítulo 2
Raúl M. del Toro
A partir de la aplicación de las técnicas de minimización del error cuadrático a los datos
obtenidos, se ajustaron modelos lineales de primer, segundo y tercer orden. En la Tabla 2.5 se
relacionan las funciones de transferencia de tiempo continuo que describen los modelos
obtenidos. Las componentes dinámicas difieren, aunque la única diferencia entre los modelos
de segundo y tercer orden es la adición de un polo más y una pequeña modificación en la
ganancia. Todos ellos tienen el mismo tiempo de retardo. El uso de una función de dos
escalones en la velocidad de avance para la obtención de los datos permite obtener un modelo
más realista que incorpora una variedad de situaciones de taladrado.
a) b)
c)
Figura 2.25. Respuesta de los modelos de primer, segundo y tercer orden de la fuerza ante cambios en la
velocidad de avance y conjunto de datos de validación.
Los tres modelos han sido simulados ante el mismo cambio de doble escalón de la
velocidad de avance (entrada Feed_Rate de la aplicación) empleado durante la etapa
experimental. En la Figura 2.25 se representa la respuesta de los tres modelos junto a la fuerza
filtrada de otro conjunto de datos de validación.
n n
IAE = ∑ y (i ) − yˆ (i ) ∑ y (i ) (2.28)
i =1 i =1
n n
∑ ( y(i) − yˆ (i) ) ∑ y 2 (i)
2
ISE = (2.29)
i =1 i =1
64 Capítulo 2
Raúl M. del Toro
Los índices de comportamiento son calculados a partir del error, e(t ) , entre la salida del
modelo, yˆ(t ) , y los datos de validación, y (t ) , tal como se expresa en la ecuación (2.6). En la
Tabla 2.6 se relacionan los índices de error calculados para cada modelo. En ambos índices, el
menor error corresponde al modelo de primer orden, aunque el modelo de segundo orden se
comporta de modo similar. El modelo de tercer orden es el que peor se comporta en este análisis,
a pesar de que proporciona una representación mejor del escalón inicial.
Tabla 2.6. Resultados del análisis de error para los modelos del taladrado de GGG40.
N −( k +1)
∑ ( y(i) − yˆ (i) )( y(i + k ) − yˆ (i) )
i =k
pˆ k = N −1
(2.30)
∑ ( y(i) − yˆ (i) )
2
i =0
La Figura 2.26 muestra el análisis de correlación, de los tres modelos del taladrado del
GGG40, realizado a partir de los datos de validación. Los tres modelos se ajustan al intervalo
de confianza del 95%, por lo que todos pueden representar la dinámica del sistema de forma
adecuada. Si embargo, los modelos de segundo y tercer orden representan la dinámica del
sistema un poco mejor que el modelo de primer orden, lo cual es un resultado opuesto al
procedimiento de análisis de error. El resultado anterior es de esperarse, debido a que un
modelo de primer orden ignora dinámica de órdenes superiores. Por otra parte, aunque estos
modelos puedan realizar un buen ajuste de los datos y provean una adecuada representación
de la dinámica del proceso, puede que no sean apropiados para ser utilizados en esquemas de
control a lazo cerrado debido a su potencial inestabilidad. Por lo que resulta necesario realizar
simulaciones para determinar la aplicabilidad del modelo en cualquier esquema de control.
Figura 2.26. Función de auto-correlación del error en la predicción de fuerza en el taladrado de GGG40.
66 Capítulo 2
Raúl M. del Toro
Figura 2.27. Aplicación de los modelos del GGG40 con ganancia unidad a la fuerza normalizada del
taladrado del material 17-4PH.
La Tabla 2.7 muestra los índices de error que se obtienen al aplicar los modelos de
taladrado obtenido con los datos del material GGG40 al material 17-4PH. Los índices de
comportamiento muestran que el los modelos constituyen una representación adecuada del
proceso de taladrado de alto rendimiento del acero inoxidable 17-4PH. Nuevamente, el
modelo que realiza un mejor ajuste es el de primer orden, seguido del segundo.
Tabla 2.7. Análisis del error de ajuste para los modelos generados del material
GGG40 empleados en el material 17-4PH.
Utilizando como variable del proceso la fuerza resultante Ftotal (φ ) , la cual se define según
la ecuación (2.23) del modelo físico del FAV, y que no es más que la norma euclídea del
vector de fuerzas de tres componentes [ Fx ; Fy ; Fz ] , la fuerza media en cada revolución de la
herramienta se obtiene como:
2π
1
Fm =
2π ∫ Ftotal dφ (2.31)
0
Empleando como variable de salida del proceso de FAV, la fuerza resultante media Fm ,
este se puede representar como un proceso dinámico múltiple entrada-simple salida (MISO),
cuyas variables de entrada serían la velocidad de avance de referencia de la herramienta con
respecto a la pieza f , la velocidad de giro de referencia de la herramienta n y la profundidad
axial de corte relativa entre la herramienta y la superficie de la pieza ap . Algunos de los
parámetros que mayor influencia tienen en el comportamiento de las fuerzas son el diámetro
de la herramienta D , la profundidad radial de corte ae , la dureza del material, entre otros. La
Figura 2.28 ilustra una representación gráfica del proceso, por otra parte, si se considera el
proceso como un sistema lineal e invariante en el tiempo, en cierto intervalo de tiempo, este
68 Capítulo 2
Raúl M. del Toro
f (s) G f (s)
n( s) Gn ( s) Fm ( s)
ap( s) Gap ( s )
Figura 2.28. Esquema del proceso de FAV y diagrama de bloques de una representación del proceso
mediante funciones de transferencia.
70mm
30mm
z y 180mm
x
Figura 2.29. Probeta con perfil escalonado en su superficie con cambios de 0.5 mm.
f0
n t1 t5 t11 t
n1
n0
ap t3 t7 t10 t
70 Capítulo 2
Raúl M. del Toro
Tabla 2.8. Valores que alcanzan cada una de las variables de entrada a lo largo del eje x de la pieza,
durante la operación de ranurado.
t [seg] x [mm] f [mm/min] n [rpm] ap [mm]
0 0 f0 n0 ap0
t1 18 f1 n0 ap0
t2 36 f1 n0 ap1
t3 54 f1 n1 ap1
t4 72 f1 n1 ap2
t5 90 f1 n1 ap2
t6 108 f0 n1 ap1
t7 126 f0 n0 ap1
t8 144 f0 n0 ap0
t9 162 f0 n0 ap0
t10 180 f0 0 0
t11 198 0 0 0
n0 = 0.8nnom . En cada uno de los intervalos cambia el tamaño de la ventana de datos para el
filtrado.
En la Figura 2.31 se muestra el comportamiento de las tres componentes de fuerzas
( Fx ; Fy ; Fz ), de la fuerza resultante y de la fuerza media, en una de las réplicas de la
operación de corte con la herramienta de 12 mm. En la misma figura, se ilustra además, los
cambios en las referencias de las velocidades de giro y de avance, durante todo el corte.
a) b)
c)
Figura 2.31. Señales registradas durante la operación de fresado con la herramienta de 12 mm.
a) Fuerzas en los tres ejes, b) fuerza resultante y media y c) velocidades de giro y avance.
Tras el filtrado de los datos de fuerza en cada uno de los experimentos, se analizó el
comportamiento de la fuerza media en cada uno de ellos, para decidir de cuáles datos se
podría obtener un mejor ajuste. La Figura 2.32 muestra la fuerza media de cada una de las
réplicas de los experimentos para las herramientas de 10 y 12 mm respectivamente. Del
experimento con la herramienta de 10 mm, se escogieron la segunda y tercera réplica para la
obtención de los modelos con esta herramienta, debido a que en la primera réplica no se
aprecia claramente el cambio en la fuerza producto del primer escalón en la velocidad de
avance. De la misma forma, del experimento con la herramienta de 12 mm, se escogieron la
72 Capítulo 2
Raúl M. del Toro
a) b)
datos para cada ensayo, definiéndose en el siguiente orden y nombrados según el cambio de la
variable a que corresponden: 0 – ap0; f0 – f1; ap0 – ap1; s0 – s1; ap1 – ap2; f1 – f0; ap2 – ap1; s1
– s0; ap1 – ap0.
Tabla 2.10. Parámetros de algunos de los modelos obtenidos con la herramienta de 12 mm.
Rép. ap f n f/n TM Subconj. Aj. % K Tp1 Tp2 Td
1 0.422 238.4 954.4 0.2498 P2 0 - ap0 86.76 434.55 0.1776 0.178 0
1 0.422 357.6 954.4 0.3747 P2D f0 - f1 85.17 0.639 0.1817 0.117 0.0037
2 0.938 357.6 954.4 0.3747 P2 ap0 – ap1 97.8 504.77 0.1677 0.0520 0
2 0.938 357.6 1431.6 0.2498 P2D s0 - s1 93.439 -0.302 0.0351 0.0437 0.0262
1 1.413 357.6 1431.6 0.2498 P1 ap1 – ap2 88.01 414.02 0.161 0 0
1 1.413 238.4 1431.6 0.1665 P1D f1 - f0 84.97 1.449 0.0809 0 0.03
2 0.929 238.4 1431.6 0.1665 P2D ap2 – ap1 96.417 232.25 0.385 0.0948 0.0299
1 0.917 238.4 954.4 0.2498 P2 s1 - s0 91.06 -0.225 0.0152 0.0162 0
2 0.415 238.4 954.4 0.2498 P2 ap1 – ap0 93.235 392.91 0.512 0.115 0
74 Capítulo 2
Raúl M. del Toro
Fm ( s)
G f (s) = ⇒ K f = g (nnom , apnom ) (2.32)
f (s)
a) b)
f ( x) = p1 x + p2 f ( x) = p1 x + p2
p1 = 1463 [1408,1518] p1 = 1640 [1100, 2180]
p2 = -0.008977 [-0.05099, 0.03304] p2 = -36.19 [-183.8, 111.5]
c)
f ( x) = p1 x + p2
p1 = -0.0006879 [-0.0007297, -0.0006461]
p2 = -0.05945 [-0.07175, -0.04715]
Figura 2.33. Funciones de ajuste para las ganancias de los sub-sistemas que componen el modelo.
Caso herramienta de 12 mm: a) K f vs apnom / nnom ; b) K ap vs f nom / nnom y c) K n vs f nom .
La Figura 2.33 representa la relación de las ganancias entre los cambios en la fuerza media
y las variables de entrada de cada sub-sistema, para los ensayos con la herramienta de 12 mm.
La Figura 2.33.a describe la curva lineal de juste de la relación entre ganancia K f del sub-
sistema G f ( s ) y el coeficiente entre la profundidad de corte nominal y la velocidad de giro
nominal durante el instante de cambio de la velocidad de avance, es decir,
K f vs apnom / nnom . De la misma forma, en la Figura 2.33.b se ilustra el ajuste de la función
que relaciona la ganancia K ap con el avance por revolución nominal en el instante de cambio
de la profundidad axial de corte, es decir, K ap vs f nom / nnom . También en esa misma figura
se muestra la función ajustada que representa la relación K n vs f nom , es decir, solamente se
ha tenido en cuenta el valor nominal de la velocidad de avance. Para cada función de ajuste se
relacionan sus parámetros y sus valores para un 95% del límite de confianza.
De la misma forma se determinaron la relación entre las ganancias de los sub-sistemas
para la herramienta de 10 mm. En la Tabla 2.11 se describen las funciones de ajuste para cada
ganancia, especificando los valores de sus parámetros y su cuantía para un 95% del límite de
confianza.
Tabla 2.11. Funciones de ajuste de las ganancias de los sub-sistemas para la herramienta de 10 mm.
f ( x) = p1 x + p2 f ( x) = p1 x + p2
Tabla 2.12. Análisis estadísticos de las constantes de tiempo de los modelos de los sub-sistemas.
Desviación
Mín. Máx. Media
Sub- Standard
Parámetro
sistema 10 12 10 12 10 12 10 12
mm mm mm mm mm mm mm mm
Tp1 0.111 0.043 0.147 0.181 0.138 0.095 0.0156 0.043
Gf Tp2 0.060 0.014 0.084 0.148 0.069 0.066 0.0101 0.043
Td 0.0034 0.14 0.03 0.03 0.049
Tp1 0.0199 0.0039 0.055 0.107 0.035 0.035 0.0123 0.031
Gn Tp2 0.0197 0.0124 0.028 0.052 0.027 0.0312 0.0045 0.016
Td 0.0005 0.03 0.03 0.026 0.013
Tp1 0.149 0.115 0.615 0.447 0.343 0.298 0.1164 0.102
Tp2 0.014 0.016 0.148 0.095 0.074 0.040 0.0438 0.025
Gap Td 0.023 0.03 0.03 0.03 0.0020
Trate 2.83 3.37 8.89 8.73 4.01 5.04 1.87 1.35
76 Capítulo 2
Raúl M. del Toro
Por otra parte, en el caso de las constantes de tiempo de los diferentes tipos de modelos, se
hizo un análisis estadístico, determinando para cada uno de ellas su valor medio, desviación
estándar, el valor mínimo y el máximo. Este análisis se muestra en la Tabla 2.12.
Para el caso del comportamiento temporal del sub-sistema que relaciona la fuerza media
de corte con el cambio en la profundidad axial de corte, se hizo un análisis más exhaustivo,
planteándose la relación existente entre la mayor constante de tiempo del modelo y el tiempo
teórico Tset al cuál debe alcanzarse el estado estacionario en la fuerza
( Trate = Tset / MAX(Tp1, Tp2) ). Los resultados de este análisis se muestran también en la
Tabla 2.12. Según trabajos previos [168], se estima que, dado el diámetro de la herramienta
D expresado en milímetros y su velocidad de avance f dada en mm/min, se alcance el
estado estacionario, para el caso de una operación de rasurado, en un tiempo Tset , expresado
en segundos, según la ecuación siguiente.
60 D
Tset = (2.33)
f 2
Kf −Tdf s
Fm ( s ) = G f ( s ) f ( s ) + Gn ( s )n(s ) + Gap ( s )ap (s ) = e f ( s) +
(T f 1s + 1)(T f 2s + 1) (2.34)
Kn K ap −Tdap s
+ e −Tdn s n( s ) + e ap ( s )
(Tn1s + 1)(Tn 2 s + 1) Tap s + 1
Las ganancias del modelo estarán determinadas por las siguientes ecuaciones:
apnom
K f = Kα _ f _1 ⋅ + Kα _ f _ 2 (2.35)
nnom
f nom
Kn = − + Kα _ n _ 2 (2.36)
Kα _ n _1
f nom
K ap = Kα _ ap _1 ⋅ + Kα _ ap _ 2 (2.37)
nnom
60 D
Tap = Tα _ ap (2.38)
f 2
Tabla 2.13. Valores de ajuste de los coeficientes para el cálculo de las ganancias del modelo.
Coeficiente Unidades 10 mm 12 mm
Kα _ f _1 2 1594 1463
[N/mm ]
Kα _ f _ 2 [N/(mm/min)] -0.167 0
Kα _ n _1 [mm/(N min 2 )] 1671.68 1453.7
Kα _ n _ 2 [N min] 0.02653 -0.05945
Kα _ ap _1 [N/mm 2 ] 1014 1640
Kα _ ap _ 2 [N/mm] 80.65 -36.19
Tα _ ap 1/5
Tabla 2.14. Valores de ajuste de las constantes de tiempo y retardos de tiempo de los sub-sistemas que
conforman el modelo del proceso de FAV.
Parámetro 10 mm 12 mm
Tf 1 0.1377 0.07
Tf 2 0.0692 0.0439
La Tabla 2.13 describe los valores ajustados a los coeficientes para el cálculo de las
ganancias del modelo del proceso de FAV, para cada tipo de herramienta, empleando el
78 Capítulo 2
Raúl M. del Toro
a) b)
Figura 2.34. Desarrollo del modelo dinámico del proceso de FAV en Matlab/Simulink.
a) Entrada y salida de datos del modelo e interfaz para introducir sus parámetros y b) diagrama de
bloques de los sub-sistemas del modelo.
80 Capítulo 2
Raúl M. del Toro
a)
b)
Figura 2.35. Simulación del modelo y comparación con los datos experimentales.
a) Herramienta de 10 mm y b) herramienta de 12 mm.
unidades de fuerza (Newton) y a través de los cuales se puede comprobar el buen desempeño
del modelo ajustado e implementado en Simulink.
N
1
IAE =
N
∑ y(i) − yˆ (i) (2.39)
i =1
N
1
∑ ( y(i) − yˆ (i) )
2
ISE = (2.40)
N i =1
Tabla 2.15. Errores del modelo, réplicas de ajuste y de validación (filas sombreadas).
Diámetro Fmax Réplica IAE (IAE / Fmax) ISE (ISE / Fmax)
[mm] [N] No [N] x 100% [N] x 100%
1 11.59 5.80 13.43 6.72
10 200 2 7.57 3.79 8.82 4.41
3 5.87 2.94 6.96 3.48
1 5.09 1.70 6.86 2.29
2 7.24 2.41 8.87 2.96
12 300
3 20.47 6.82 27.31 9.10
4 14.58 4.86 19.43 6.48
Por otra parte, en la Tabla 2.15 también se muestra el valor máximo que alcanza el cambio
de la fuerza (Fmax) durante todo el rango de datos de validación y se ha calculado el valor
porcentual de los índices con respecto al valor máximo (quinta y séptima columnas),
obteniéndose valores por debajo del 10 % y corroborándose de esta manera la precisión del
modelo ajustado.
2.5 Conclusiones
En este capítulo se ha presentado un enfoque computacional en el modelado de procesos
de fabricación de alta complejidad. Han sido diseñado modelos a partir de ecuaciones físico-
matemáticas de los procesos de fabricación y también, basados en la aplicación de técnicas de
identificación experimental. Se han desarrollado modelos simple entrada-simple salida
(SISO), múltiple entrada-simple salida (MISO) y múltiple entrada-múltiple salida (MIMO).
Los modelos han sido implementados mediante aplicaciones desarrolladas en un lenguaje de
alto nivel como Matlab/Simulink.
Por otra parte, se ha mostrado la efectividad de los modelos, para la predicción de las
fuerzas que se producen en los procesos de fabricación, a través de un estudio comparativo
con datos obtenidos experimentalmente. Los modelos computacionalmente eficientes pueden
ser utilizados en el diseño, ajuste y optimización de estrategias de control, que a su vez
82 Capítulo 2
Raúl M. del Toro
85
Métodos clásicos y de soft-computing en la optimización de procesos complejos.
1
t
de ( t )
∆f ( t ) = K p ⋅ e ( t ) + ∫ e (τ ) dτ + Td (3.1)
T dt
i 0
siendo, t el tiempo, e ( t ) la señal de error (entrada) entre el valor de referencia del sistema de
control y la salida del proceso ( e(t ) = FR (t ) − F (t ) ), ∆f (t ) la señal de control (salida del
controlador) y los tres parámetros básicos del controlador dados por su ganancia proporcional
K p , la constante de tiempo integral Ti y la constante de tiempo derivativa Td .
La ecuación del controlador puede expresarse también como la suma de tres términos: la
acción proporcional, la acción integral y la derivativa
86 Capítulo 3
Raúl M. del Toro
∆f ( t ) = ∆f P ( t ) + ∆f I ( t ) + ∆f D ( t )
∆f P ( t ) = K p e ( t )
t
∆f I ( t ) = K I ∫ e (τ ) dτ (3.2)
0
de ( t )
∆f D ( t ) = K D
dt
donde K I = K p / TI y K D = K pTD .
∆f ( s ) K
GPID ( s ) = = K p + I + KD ⋅ s (3.3)
e(s) s
∆f [k ] = ∆f P [k ] + ∆f I [k ] + ∆f D [k ]
∆f P [k ] = K P e[k ]
(3.4)
∆f I [k ] = ∆f I [k − 1] + ( K I Ts ) e[k ]
∆f D [k ] = ( K D / Ts ) (e[k ] − e[k − 1])
Debido a que la acción derivativa del controlador es sensible a los cambios de la señal de
error, el ruido que puede existir en esta señal, algo muy común en los procesos físicos reales,
puede afectar el comportamiento del sistema de control. Por tanto una versión del controlador,
ampliamente utilizada, añade un filtro pasa-bajo en la componente derivativa:
Td K p ⋅ Td ⋅ N
∆f D [k ] = ∆f D [k − 1] + ⋅ (e[k ] − e[k − 1]) (3.5)
Td + N ⋅ Ts Td + N ⋅ Ts
Parámetros
Kp Ti Td
Métodos
Ziegler-Nichols – 0.6 ⋅ Ku 0.5 ⋅ Tu 0.125 ⋅ Tu
dominio frecuencial
Cohen-Coon
τ 2.5 − 2 ⋅τ 0.37 − 0.37 ⋅τ
L
a = K1 ,τ =
L (1.35 / a ) ⋅ 1 + 0.18 ⋅ L⋅ L⋅
1 −τ 1 − 0.39 ⋅τ 1 − 0.81⋅ τ
Tc L + Tc
Visioli (procesos con
retardo de transporte 1.37
L y ganancia a estado ( K1 ⋅ L ) 1.49 ⋅ L 0.59 ⋅ L
estacionario K1 )
88 Capítulo 3
Raúl M. del Toro
Por otra parte, el método de Cohen-Coon está basado en un modelo del proceso de primer
orden y con retraso de transporte (3.6), y tiene como objetivo principal el rechazo a
perturbaciones de un sistema de control. Otro de los métodos a resaltar de uso más reciente es
el de Visioli [171], basado también en un modelo lineal del proceso con retraso de transporte.
Su objetivo primario es la minimización de la integral del error ya sea ante presencia de
perturbaciones o ante cambios en la variable de referencia del sistema de control. En la Tabla
3.1 se resumen los tres métodos de ajustes, así como las ecuaciones para el ajuste de los
parámetros del controlador PID.
K1 − sL
GP = e (3.6)
1 + sTc
90 Capítulo 3
Raúl M. del Toro
Sin embargo, las aplicaciones generalizadas y el desarrollo teórico de los últimos 20 años,
se han encargado de darle a esta teoría su justo valor, probablemente menos “milagrosa” de
como nos la describiera Zadeh al principio, pero indudablemente muy útil y poderosa.
A día de hoy, es muy difícil de refutar, que la lógica borrosa ofrece una alternativa
adicional, una herramienta más, para la solución técnica al problema de toma de decisiones y
al control inteligente de procesos, es decir, que ha venido a enriquecer el arsenal a disposición
de los diseñadores de sistemas de monitorización y control, con un enfoque muy novedoso
para emular la toma de decisiones humana. La idea de incorporar al regulador la experiencia
de los operadores expertos y su forma de actuar sobre el proceso, le proporciona el
“conocimiento” necesario para hacer frente a la complejidad y severas no linealidades de
muchos sistemas.
µ A ( x) = 1 ⇔ x ∈ A
(3.7)
µ A ( x) = 0 ⇔ x ∉ A
A = {( x, µ A ( x ) ) , x ∈ U } (3.8)
a) b)
Para el caso de universos discretos U = {x1 ,..., xn } , una notación quizá más conveniente
sugerida por Zadeh es
µ A ( x1 ) µ A ( xn ) n
µ A ( xi )
A= x1 + ⋯ + xn = ∑
i =1
xi (3.9)
A = ∫ µ A ( x) (3.10)
x
U
Aα = { x ∈ U , µ A ( x ) ≥ α } (3.12)
Se dice que A está normalizado si y solo si existe al menos una x ∈ U , tal que µ A ( x) = 1 .
El conjunto vacío ∅ se define como µ∅ ( x) = 0, ∀x ∈ U , claro que µU ( x) = 1, ∀x ∈ U . Una
92 Capítulo 3
Raúl M. del Toro
µ A∪ B ( x ) = max ( µ A ( x ), µ B ( x ) ) ∀x ∈ U (3.14)
µ A∩ B ( x ) = min ( µ A ( x ), µ B ( x) ) ∀x ∈U (3.15)
A ∨ B ó µ A ∨ µ B para el máximo,
A ∧ B ó µ A ∧ µ B para el mínimo.
µ A∪ B ( x ) ≥ max ( µ A ( x ), µ B ( x) )
∀x ∈ U (3.16)
µ A∩ B ( x ) ≤ min ( µ A ( x ), µ B ( x ) )
8. ∀x ∈ U , f (1,µ A ( x ) ) = 1; g ( 0,µ A ( x ) ) = 0
µ A ( x) = 1 - µ A ( x ) ∀x ∈ U (3.18)
Un rasgo significativo que distingue a los conjuntos borrosos de los clásicos es que, para
conjuntos borrosos, en general A ∩ A ≠ ∅ ; A ∪ A ≠ U ∀A ∈ U .
Aún más, en general, si A es cierto concepto o idea, se dice que tal concepto es borroso
cuando no se distingue claramente entre A y A.
Aunque los útiles diagramas de Venn no existen para los conjuntos borrosos, en el sentido
de cómo son utilizados para sub-conjuntos ordinarios, pueden usarse los gráficos de las
funciones de pertenencia para visualizar las operaciones más arriba definidas. Sin embargo,
las definiciones para las operaciones de unión e intersección borrosas no son únicas. Otros
operadores empleados, similares a sus homólogos probabilísticos, son:
• Para la unión, conocida como suma algébrica ó probabilística:
µ A⋅B ( x) = g ′ ( µ A ( x ), µ B ( x) ) = µ A ( x) ⋅ µ B ( x) ∀x ∈ U (3.20)
94 Capítulo 3
Raúl M. del Toro
A ⋅∅ = ∅ ; A+U=U
ˆ . Para estos operadores también pueden construirse diagramas de Venn
(ver Figura 3.5).
a) b)
U ×V = { x, y , x ∈ U y y ∈ V} (3.21)
1 n
{
µ A ×⋯× A ( x1, x2 ,… , xn ) = min µ A ( x1 ),… , µ A ( xn )
1 n
} (3.22)
ó, en su lugar
T
dimensiones, donde, a cada vector µ A1 ( x1 ) ⋯ µA n ( xn ) le corresponde un punto, cuya
restante coordenada se calcula por medio de (3.22) ó de (3.23). El tema de las operaciones
lógicas borrosas se retomará más adelante, al abordar la lógica de toma de decisiones del
controlador borroso.
∗”
3.2.4 Composición “Sup-∗
Si R y S son relaciones borrosa en U×V y V×W respectivamente, la composición de R y S
es una relación borrosa denotada por R S y se define por:
R S = ( u , w ) ,sup ( µ R ( u, v ) * µS ( v, w ) ) , u ∈ U , v ∈V , w ∈ W (3.24)
v
Donde “∗” puede ser cualquier operador perteneciente a la clase de las llamadas "normas-
T" (mínimo, producto algébrico, producto acotado, producto drástico) [183].
y = xR (3.25)
96 Capítulo 3
Raúl M. del Toro
µ ALTA = 1, F ≥ 2000 N
F
M ALTA = µ ALTA = − 3, 1500 N ≤ F ≤ 2000 N
500
µ ALTA = 0, F ≤ 1500 N
La regla semántica M permitiría entonces precisar, en nuestro ejemplo, que “las fuerzas se
califican plenamente (100%) como “altas” por encima de los 2000 N, valor de fuerza en la
que la pertenencia a este conjunto decae linealmente hasta ser nula en 1500 N”. O sea
y la consecuencia
98 Capítulo 3
Raúl M. del Toro
TAMBIÉN
variable de conjunto variable de conjunto
Ri +1 : ES borroso Y ⋯ ENTONCES acción 1 ES borroso ;⋯
estado 1
TAMBIÉN
⋮
Todas las reglas del protocolo son, en principio, válidas y aplicables, lo que explica la
presencia del conectivo “TAMBIÉN”. Claro que, ante cada estado del proceso, la validez y
efectividad de cada regla será “modulada” por el grado de cumplimiento de sus respectivas
premisas. La teoría de conjuntos borrosos debe proveer entonces las herramientas
matemáticas necesarias para crear los correspondientes algoritmos de cálculo.
El DLB maneja cada una de los componentes de ambos vectores como variables
lingüísticas. De modo que, para una cualquiera de las variables de estado xi y acción uh , los
conjuntos de atributos serían, respectivamente:
{ }
T ( xi ) = A1i , A2i , …, Ami i ; i = 1, 2, …, n
(3.26)
T (u ) = {B , B , …, B } ;
h
h
1
h
2
h
qh h = 1, 2, …, ñ
donde Aij , Blh son conjuntos borrosos, mi es el número de conjuntos borrosos en que se parte
xi y qh el número de conjuntos borrosos en que se parte uh .
Si A entonces B = A ⇒ B (3.27)
A ⇒ B = A×B = R (3.28)
Ri = Ai × Bi ; i = 1, 2,..., p (3.29)
100 Capítulo 3
Raúl M. del Toro
n
Ai = ∩ Ai j (3.31)
j =1
El primero tiene, para el DLB, una esencia estrictamente numérica. En cada instante,
como resultado de la medición de las n variables de estado y su ulterior emborronado, se
tiene, para cada regla, un conjunto de n valores de pertenencia µ A j ( x 0j ) ; i = 1, 2,..., p ;
i
El segundo tiene otro cariz, pues, si bien el resultado de Ai = Ai1 ×⋯ × Ain es un número
(αi∈[0,1]), Bi está caracterizado por una curva, la función de pertenencia µ B j (u ) , y
i
p
R = ∪ Ri (3.32)
i =1
u = g ( x) (3.33)
donde, para un cierto estado del sistema x º , se obtendría una acción de control u º
u º = g ( xº ) (3.34)
y g es una función vectorial no lineal. Esta no linealidad suele ser muy marcada y responde
estrictamente a la presunta reproducción que hace el DLB de las habilidades del operador en
la dirección del proceso, por tanto la linealización no sería conceptualmente adecuada, ni
matemáticamente expedita. No obstante, no han faltado investigadores que la hayan intentado.
Si lo consideramos como parte del proceso, entonces el dispositivo es verdaderamente
borroso (aplica conjuntos borrosos sobre conjuntos borrosos). El primer punto de vista
subraya el carácter discreto del sistema y hace recordar a un relé de varios niveles. Mientras
tanto el segundo criterio pudiera ser más rico y novedoso conceptualmente, pues refleja un
dispositivo capaz de generar acciones cualitativas ante estados cualitativos.
Por tanto, de forma resumida el DLB se puede considerar como un sistema de aplicación
de un espacio de entrada Rn a un espacio de salida Qm con las siguientes propiedades:
102 Capítulo 3
Raúl M. del Toro
Emborronado
Este bloque ha de jugar el papel de interfaz entre las expresiones deterministas y las
borrosas (cualitativas) de las variables de estado del proceso. Pensando en un operador
humano, aunque que, si bien él interpreta cualitativamente el estado (la salida) del proceso, lo
cierto es que objetivamente cada variable tiene, a cada instante, un valor (determinista) dado y
tales valores son los que entran al DLB. De hecho, el emborronado sería la realización
práctica de la definición de las variables como lingüísticas y de la partición de los universos
de discurso: a cada valor (determinista) de la variable de estado, le corresponde cierto grado
de pertenencia a los conjuntos borrosos de su respectivo universo, por lo que éstas se
convierten en singularidades borrosas.
En la elección del número de particiones, o sea, el número de elementos del T(x) de la
quíntupla definida para la variable lingüística (ver Figura 3.6), influye decisivamente el
método que se emplee para la creación del algoritmo del DLB. Uno de ellos, quizá el más
popular, conocido como verbalización, es decir, la creación del algoritmo se basa
exclusivamente en descripciones verbales de los operadores, acerca de su actuación en la
dirección del proceso y resulta evidente que el número de conjuntos borrosos elegido para una
variable tiene una esencia fuertemente subjetiva. No sería extraño que, para un mismo
proceso, un operador considere que una variable lingüística dada tenga, digamos, 3
particiones cualitativas y otro considere que son 5 y como resultado, ambos dirigen
"satisfactoriamente" el proceso.
Las funciones de pertenencia pueden ser expresadas de diferente forma. En la Tabla 3.2 se
muestran algunas de las más usualmente utilizadas. Una excelente colección de estas
funciones se ofrece en un trabajo de Kaufmann [180]. Como ya se ha mencionado, la elección
del tipo de función es relativamente arbitraria y esta sujeta a factores tales como el tipo de
DLB, los requisitos de cómputo, entre otros [187].
Se podría resumir entonces que el emborronado no es más que la transformación de la
información determinista, enviada desde el entorno físico al DLB, en información cualitativa
que usa como referencia a conjuntos borrosos. También se evidencia de la explicación, que
esta transformación depende de parámetros del DLB fijados a priori, vinculados con la
caracterización de las variables medidas del proceso como variables lingüísticas: número de
conjuntos y forma de las funciones de pertenencia.
Triangular
0, x ≤ a; x ≥ c
x −a
µ ( x) = , a≤ x≤b
b − a
c− x
c − b , b ≤ x ≤ c
Trapezoidal
0, x ≤ a; x ≥ d
x−a
, a≤ x≤b
b−a
µ ( x) =
1, c ≤ x ≤ c
d − x
,c≤x≤d
d −c
Campana de Gauss
− ( x − c )2
µ ( x) = e 2σ 2
Base de Conocimiento
La base de conocimiento ha de incluir todos los datos e informaciones suministradas al
DLB, tales que posibiliten el procesamiento matemático de las variables medidas del proceso,
a fin de generar las salidas del DLB. Muchos autores sólo le adjudican las llamadas bases de
reglas y bases de datos, ambas relacionadas con informaciones muy específicas del sistema
experto borroso, como tal. Otros autores dan preferencia a completar esas informaciones con
104 Capítulo 3
Raúl M. del Toro
otras propias de cualquier controlador [186]. Éstas abarcan los valores de referencia, valores
límite de las variables, constantes varias de ajuste, parámetros necesarios para la
discretización y normalización de las variables, entre otros. El número de niveles de
cuantificación resulta de un compromiso entre la precisión y sensibilidad a que se aspira, y la
memoria dinámica y velocidad de cómputo de que se dispone. Los intervalos más usuales
para la normalización de los universos de discurso son [-1,1], [0,1], [-100%,100%] y
[0,100%].
El conjunto de datos que conforman la base de datos incluye, en primer lugar, a todos
aquellos imprescindibles para concretar la partición de todas las variables lingüísticas:
Universos de discurso, número y nombre de las particiones así como la geometría de las
funciones de pertenencia. Igualmente importantes son los relativos a las posibles alternativas
para las operaciones y cálculos que han de realizarse en las etapas de toma de decisiones y
desemborronado. Las características de la base de reglas determinan la necesidad de una
atención particular a éstas.
Base de reglas
De todo lo comentado con anterioridad se deduce claramente la importancia de las reglas
en su calidad de representación estructurada de “la política de manejo del operador”. La base
de reglas comprende pues el conjunto de reglas de producción que deben ser prefijadas en el
DLB. El número máximo posible de reglas viene dado por el producto de los números de
particiones de todas las variables de entrada al DLB.
Para hablar de una base de reglas, es preciso ante todo elegir cuáles variables se utilizarán
como de entrada y cuáles como de salida del DLB. No es ocioso subrayar que una selección
adecuada de las variables de estado y acción del proceso resulta esencial para la
caracterización y ulterior comportamiento del CLB. La experiencia y el conocimiento técnico
desempeñan un papel importante en esta etapa de selección [186]. Más que otra cosa, se
impone la necesidad de un análisis con enfoque sistémico, dónde emerjan no sólo las
variables de entrada y salida, sino también los objetivos, restricciones y parámetros. Se
pudiera decir que si esto es en general importante para cualquier sistema de control, lo es
mucho más en el caso de los llamados procesos “complejos” (a los que preferentemente se
aplica el control borroso) especialmente a la hora de concebir su base de reglas. Ejemplos
típicos de variable de entrada de un DLB son el estado, el error del estado, la integral del
error, etc.
Para crear las reglas que se prefijarán inicialmente en el DLB, existen tres vías que se
clasifican como: verbalización, emborronado, e identificación [188].
A la primera ya se ha hecho referencia con anterioridad, al describir la etapa de
emborronado del DLB. La segunda vía, denominada emborronado (no confundir con la etapa
de igual nombre del DLB) proporciona modelos borrosos a partir de expresiones matemáticas
ordinarias por medio del principio de extensión de Zadeh. Tales expresiones matemáticas no
suelen existir en los casos de procesos complejos. Por otra parte, la identificación supone la
generación de descripciones causales, a partir de datos numéricos obtenidos midiendo las
variables de estado y acción, durante la actuación de operadores eficientes.
Sobre la verbalización, la mayoría de los trabajos pioneros del control borroso de la
década de los 70 se apoyaron en modelos basados en descripciones verbales de los
operadores, tecnólogos, etc., básicamente a partir de dos métodos heurísticos: la interrogación
del operador a partir de un cuestionario riguroso y/o una verbalización de la experiencia del
operador. A pesar de estar afectadas por fuertes dosis de subjetivismo, constituye la técnica
más popular.
La identificación, por su parte, provee la objetividad que emana de mediciones reales,
procesando los archivos de datos de estas mediciones, para generar la base de conocimiento.
Existen para ello diversas técnicas, entre las cuales está el emborronado de cada medición,
convirtiéndola en un “número borroso” [189] y las técnicas de agrupamiento borroso (fuzzy
clustering) que se han empleado para diseñar tanto controladores como modelos borrosos
[190].
Takagi y Sugeno propusieron un interesante algoritmo que también se podría clasificar
entre los de identificación, pues se basa en datos procedentes de mediciones sobre el proceso
[24]. Las reglas que propusieron para su base de conocimiento eran de un nuevo tipo, con
antecedente borroso pero consecuente determinista (algébrico):
106 Capítulo 3
Raúl M. del Toro
Tipos de reglas
Si bien es cierto que las más populares y conocidas reglas de un DLB son de la forma:
La lógica borrosa pone también a disposición del diseñador del sistema una muy rica
variedad de tipos de reglas, lo que da una enorme flexibilidad en las opciones. A continuación
se relacionan varios tipos recogidos en la literatura:
• Jerarquía de reglas complejas, empleadas en dispositivos borrosos expertos [195]. Las
variables intermedias representan algún estado implicado (que no se puede medir) que
puede afectar las acciones de control a aplicar:
• De igual modo, existe la posibilidad de que las reglas sean función de objetivos y
restricciones que pudieran usarse con propósitos de adaptación, de la forma
Mención aparte merece la regla de “tipo Sugeno” mencionada con anterioridad (3.35),
cuyo rasgo distintivo radica en el carácter determinista (función ordinaria) del consecuente y
constituye uno de los primeros intentos, de reglas con consecuente no borroso. Vale la pena
insistir finalmente, citando esta vez al propio Tong, que “la construcción de la base de reglas
es el aspecto crucial y más difícil del diseño del CLB”.
108 Capítulo 3
Raúl M. del Toro
SI < estado borroso del proceso > ENTONCES < acción de control borrosa >
donde cada atributo Ai j ó Bi no son necesariamente diferente para cada regla. Los pasos
siguientes se corresponden a operaciones lógicas borrosas. Pero antes es necesario tener en
cuentas algunas consideraciones matemáticas acerca de la definición de los conectivos “Y” y
“TAMBIÉN”.
En la mayoría de los DLB el conectivo “Y” se traduce prácticamente como una
implicación borrosa en un espacio producto Cartesiano, en el que, las variables involucradas
se desenvuelven en universos de discurso diferentes. Más generalmente, la implicación
borrosa se asocia a las llamadas normas triangulares. Una norma triangular o T-norma es una
función de dos lugares de [0,1]×[0,1] en [0,1], que debe cumplir los siguientes requisitos:
1. T (0, 0) = 0 ; T (a,1) = T (1, a ) = a
a si b = 1
Tw (a, b) ≤ T (a, b) ≤ min(a, b) , siendo Tw ( a, b ) = b si a = 1
0 en otro caso
Tales requisitos los cumplen, junto a las definiciones (3.22) y (3.23), con xi ∈ [0,1] para
toda i:
intersección {
xi ∩ x j = min xi , x j } (3.41)
producto acotado {
xi ⊙ xj = max 0, xi + x j − 1 } (3.43)
xi x j = 1
producto drástico ⋅ x j = x j xi = 1
xi ∩ (3.44)
0 x i , x j < 1
Las más utilizadas en el cálculo del producto Cartesiano de los DLB constituyen la
intersección (Regla de implicación del mínimo de Mamdani) y el producto algébrico (Regla
de implicación del producto de Larsen), que aparecen reflejadas en las ecuaciones (3.22) y
(3.23).
Sobre el conectivo “TAMBIÉN” cabe apuntar, en primer lugar, que debe exhibir
propiedades de conmutatividad y asociatividad, pues el orden de las reglas es irrelevante.
Tales propiedades las poseen las antes mencionadas normas triangulares, aunque también las
tienen las llamadas co-normas triangulares o T-conorma. Las co-normas son también
funciones de dos lugares de [0,1]×[0,1] en [0,1] y las condiciones que deben cumplir son
homólogas (o coincidentes) con las de las normas, a saber:
1. S (1,1) = 1 ; S (a, 0) = S (0, a ) = a
110 Capítulo 3
Raúl M. del Toro
a si b = 0
S w (a, b) ≥ S (a, b) ≥ max(a, b) , siendo S w ( a, b ) = b si a = 0
1 en otro caso
Esas condiciones las cumplen las siguientes definiciones, válidas para xi ∈ [0,1] , para toda
i, como:
unión {
xi ∪ x j = max xi , x j } (3.45)
suma acotada {
xi ⊕ x j = min 1, xi + x j } (3.47)
xi x j = 0
. x = x
suma drástica xi ∪ j j xi = 0 (3.48)
1 x i , x j > 0
suma disjunta { }
xi ∆x j = max min( xi ,1 − x j ), min(1 − xi , x j ) (3.49)
3.2.10 Desemborronado
El bloque de desemborronado constituye la interfaz de salida del DLB, es el responsable
de entregar al proceso físico, acciones deterministas a partir de conjuntos borrosos generados
durante la toma de decisiones. La estrategia de desemborronado ha de tener como objetivo
que las salidas deterministas resultantes sean las que “mejor” representen a esos conjuntos
borrosos.
u ° = min arg sup ( µ R (u ) ) (3.50)
u∈U
u
La estrategia MDM genera como salida del DLB, los valores de u que representen la
media de todos los valores que alcanzan el mismo máximo en la distribución posibilística
resultante. Para el caso de universo discreto, la salida del DLB puede expresarse como:
1
u° = ∑ uj
m u ∈V
(3.51)
j
El método del centro de área (CDA), también conocido como del centroide, del centro de
masa o de gravedad, se basa en calcular, como salida determinista, los valores de la variable
u = u º tal que divida en partes iguales el área bajo la curva de distribución posibilística
resultante de la composición de las reglas. Analíticamente sería:
u° =
∫ µ R (u)udu (3.52)
∫ µ R (u)du
En el caso de un universo discreto se calcula como:
112 Capítulo 3
Raúl M. del Toro
∑ j u j µ R (u j )
u° = (3.53)
∑ j µ R (u j )
El método del CDA es el de mayor grado de coincidencia con el sentido común y la
intuición de los especialistas a la hora de elegir la mejor uº. La tendencia bastante frecuente,
es de salidas más drásticas del MDM respecto al CDA y aún más drásticas para el máximo.
Igualmente la composición “Sup-producto” da mayores valores de salida en igualdad de otras
condiciones.
Desde inicios de los 80, varios investigadores estudiaron comparativamente estos tres
métodos de desemborronado [197]. Es bastante generalizada la conclusión de que el método
CDA ofrece los mejores resultados para determinar el mejor uº. No obstante se plantea que,
desde el punto de vista del comportamiento dinámico de un sistema de control a lazo cerrado,
el método MDM ofrece un mejor comportamiento del sistema en régimen transitorio mientras
que el CDA da mejor comportamiento a estado de régimen estacionario. Desde el punto de
vista de la teoría de control clásica, empleando el método MDM, el DLB se comporta como
un relé de varios niveles y con el CDA se obtienen resultados parecidos a los de un
controlador proporcional integral (PI).
región de entrenamiento. Por otra parte, los sistemas de inferencia borrosa proporcionan un
entorno o marco eficiente de computación basado en la imitación del razonamiento humano
mediante el empleo de conjuntos borrosos y un modelo lingüístico constituido por un
conjunto de reglas borrosas SI-ENTONCES (IF-THEN).
La hibridación de estos dos paradigmas permite aprovechar las ventajas de ambos,
consiguiendo la sencillez de ajuste de las redes neuronales y a la vez disponiendo del
conocimiento explícito expresado mediante el modelo lingüístico de los sistemas borrosos.
A partir de la segunda mitad de los años 80, aparecen en la literatura diversos trabajos en
los cuales se introducen diferentes estructura de sistemas híbridos neuro-borrosos y que
pueden ser separados en dos grandes grupos. Uno de ellos son las redes neuronales dotadas
con la habilidad de manejar información borrosa para mejorar su capacidad de convergencia
(Fuzzy Neural Networks, FNN) [198, 199]. El otro son los sistemas borrosos combinados con
redes neuronales con el objeto de incorporar la capacidad de aprendizaje y así mejorar su
precisión (Neural-Fuzzy Systems, NFS) [52, 200]. Por otra parte, en un trabajo anterior
desarrollado por Farag et ál. se añade, además, las potencialidades de optimización de los
algoritmos de computación evolutiva (por ejemplo, algoritmos genéticos) para el ajuste de las
funciones de pertenencia de los sistema neuro-borrosos [78]. Una revisión en profundidad de
las estrategias neuro-borrosas disponibles en la literatura, puede realizarse en los trabajos de
Nauck et ál. [201] y de Li Xin Wang [202].
Referido a los sistemas borrosos que incorporan la capacidad de aprendizaje de las redes
neuronales, en los últimos 15 años es posible encontrar, disponibles en la literatura, diversos
algoritmos computacionales que emplean esta estrategia de hibridación para el ajuste de
sistemas de inferencia borrosa. Entre ellos caben destacar los siguientes: ANFIS (Adaptive
Network based Fuzzy Inference System) [52], FuNN (Fuzzy Neural Network) [203, 204],
REFuNN (Rules Extraction Fuzzy Neural Network) [204], DENFIS (Dynamic Evolving
Neural-Fuzzy Inference System) [205], HyFIS (Hybrid neural Fuzzy Inference System) [200],
TWNFI (Transductive Neuro-Fuzzy Inference model with Weighted data normalization) [206]
que facilitan tareas como el modelado y el control de sistemas no lineales y de gran
complejidad [73, 207]. Es necesario destacar que varios de estos aportes han sido realizados
por el grupo de investigación del Knowledge Engineering/Discovery Laboratory de la
Universidad de Otago de Nueva Zelanda, liderado por el Dr. Nikola K. Kasabov.
Entre los diferentes algoritmos computacionales neuro-borrosos mencionados, el conocido
como Adaptive Network based Fuzzy Inference System (ANFIS) propuesto por Jyh-Shing R.
Jang [52] es, además del trabajo pionero, uno de los más sencillos desde el punto de vista
computacional y viable para aplicaciones en tiempo real.
ANFIS está basado en redes neuronales artificiales adaptativas, introducidas también por
Jang [52] y las cuales no son más que redes neuronales artificiales multicapas con
114 Capítulo 3
Raúl M. del Toro
alimentación hacia delante, diferenciándose en que las funciones de activación pueden ser
distintas en cada nodo. La estructura y métodos de aprendizaje de la red neuronal, están
diseñados de tal forma que puedan constituir un marco unificador de casi todos los diferentes
paradigmas de redes neuronales artificiales con aprendizaje supervisado [208]. A
continuación se realiza una breve introducción a la arquitectura de la red y sus procedimientos
de aprendizaje.
En general la red puede ser heterogénea, es decir, cada nodo puede tener una función de
activación específica diferente a los otros. Los enlaces entre nodos solo especifican la
dirección de propagación de las salidas de los nodos y por lo general no existen pesos u otros
parámetros asociados, como en otros tipos de redes neuronales artificiales. La Figura 3.8
muestra una red adaptativa típica con dos entradas y dos salidas. Los nodos cuadrados
representan nodos adaptativos mientras que los circulares son fijos, es decir, no tienen
parámetros. El conjunto de parámetros de toda la red adaptativa es la unión del conjunto de
parámetros de cada nodo adaptativo.
La red puede ser entrenada de forma híbrida, empleando el gradiente descendiente para
ajustar los parámetros de las funciones de activación de los nodos y el método de mínimos
cuadrados para ajustar nodos con funciones lineales. Además, el ajuste de los parámetros
puede realizarse mediante entrenamiento por lotes (aprendizaje off-line), es decir, la
actualización de los parámetros se realiza a partir de un conjunto de datos de entrenamiento.
También el entrenamiento puede llevarse a cado de forma incremental (aprendizaje on-line),
es decir, los parámetros son actualizados para cada par de datos de entrada-salida.
A continuación se describe brevemente el método del gradiente descendiente que
constituye la estrategia básica de aprendizaje. Se parte del supuesto de una red adaptativa de
L capas, donde la capa k tiene N (k ) nodos. El nodo en la posición i en la capa k se denota
por (k , i) , su función de salida por fi k y su salida por xik . Dado que la salida de cada nodo
depende de sus entradas y de su conjunto de parámetros ( α , β , γ ,... ), esta se expresa como
(
xik = fi k x1k −1 ,..., xNk −(1k −1) , α , β , γ ,... ) (3.54)
∑ ( dm, p − xmL , p )
N ( L) 2
Ep = (3.55)
m=1
E = ∑ p =1 E p
P
(3.56)
∂E p
∂xiL, p
(
= −2 di , p − xiL, p ) (3.57)
116 Capítulo 3
Raúl M. del Toro
Para el nodo interno (k , i) , la derivada del error se calcula aplicando la regla de la cadena:
∂E p N ( k +1) ∂E p ∂f mk,+p1
∂xik, p
= ∑ ∂xmk +, 1p ∂xik, p
(3.58)
m =1
∂E p ∂E p ∂f *
∂α
= ∑ (3.59)
x*∈S ∂x * ∂α
P ∂E
∂E
=∑
p
(3.60)
∂α p =1 ∂α
∂E
∆α = −η (3.61)
∂α
siendo η el factor de aprendizaje y según una longitud k (tamaño del paso) definida de cada
transición del gradiente en el espacio del parámetro, puede determinarse así
k
η= (3.62)
2
∂E
∑α ∂α
El valor de k puede cambiarse para variar la velocidad de convergencia del algoritmo,
pudiéndose utilizar para ello estrategias heurísticas.
Como se había comentado con anterioridad, el aprendizaje de la red puede realizarse
mediante dos paradigmas. Por una parte, por medio del aprendizaje por lotes (off-line) donde
la actualización de los parámetros se realiza a partir de un gran conjunto de datos de
entrenamiento y para la actualización se utiliza la ecuación (3.60). Por otra parte, a través de
un aprendizaje incremental (on-line) donde los parámetros pueden ser actualizados mediante
un par de datos entrada-salida, empleando en este caso la ecuación (3.59). Esta última
estrategia resulta de gran utilidad para la identificación, en tiempo real, de sistemas dinámicos
con parámetros variables en el tiempo.
Aprendizaje híbrido.
El proceso de identificación del conjunto de parámetros de la red a través del gradiente
descendiente, en ocasiones, resulta lento hasta que logra una convergencia y puede quedar
atrapado en mínimos locales. En el caso de que alguna salida de la red neuronal sea lineal con
respecto a algunos de los parámetros, el algoritmo de aprendizaje puede combinarse con el
método de mínimos cuadrados y de esta forma acelerar el proceso de identificación de
parámetros.
Por ejemplo, si se asume una red adaptativa con una sola salida, a partir de un vector de
entrada a la red i y el conjunto de parámetros S , la salida puede representarse de esta forma
o = F (i,S) (3.63)
118 Capítulo 3
Raúl M. del Toro
Aθ = B (3.65)
determinadas capas de la red, los parámetros de las funciones de pertenencia de las entradas y
los coeficientes de los polinomios de salida de cada regla.
antecedentes consecuentes
O1,i = µ Ai ( x) (3.66)
120 Capítulo 3
Raúl M. del Toro
1
µ A ( x) = (3.67)
i
1 + [(( x − di ) / ai )2 ]bi
O2,i = wi = ∏ µ Ai ( x) (3.68)
wi
O3,i = wi = , i = 1, 2,...,5 (3.69)
∑ wi
i
∑wh i i
O5,1 = h = ∑ wi hi = i
, i = 1, 2,...,5 (3.71)
i ∑w i
i
Como puede observarse en la arquitectura ANFIS, dados los parámetros premisas, a partir
de las ecuaciones (3.71) y (3.70) la salida final del sistema puede expresarse como una
combinación lineal de los parámetros consecuentes
m1
c1
h = [ w1 x w1 ⋯ w5 x w5 ] ⋅ ⋮ = A ⋅ θ (3.73)
m5
c5
∂E
ai (n + 1) = ai (n) −η (3.74)
∂ai
122 Capítulo 3
Raúl M. del Toro
3.4 Conclusiones
En este capítulo se ha realizado una revisión de tres técnicas que por su consolidado
cuerpo de doctrinas e impacto científico y técnico se han considerado en esta Tesis Doctoral
con vistas a la regulación automática de procesos complejos. Se han expuesto las ideas
conceptuales y básicas de estrategias de control clásicas lineales, como el algoritmo de control
proporcional, integral y derivativo, PID, de muy amplia utilización y basadas en soft-
computing, como es el caso de los sistemas de inferencia borrosa. Por otra parte, se ha
abordado el estudio de estrategias de hibridación de los sistemas borrosos con otras técnicas
de inteligencia artificial, como las redes neuronales artificiales, combinando de esta forma el
estilo de razonamiento humano de los sistemas borrosos con la estructura conexionista y los
algoritmos de aprendizaje de las redes neuronales. Esto implicará, por tanto, una mejora del
desempeño de los sistemas borrosos en la monitorización y control de sistemas complejos.
Los algoritmos de control estudiados serán examinados en el Capítulo 5 para el diseño de
sistemas de control de procesos complejos. Se examinará, la potencialidad de los sistemas
híbridos neuro-borrosos y del control clásico PID para el control de un proceso de fabricación.
Además, como parte de la metodología de diseño de sistemas de control, se comprobará
también el desempeño de estrategias de hibridación de la lógica borrosa y diferentes métodos
metaheurísticos, para el ajuste óptimo del DLB como estrategia de control de procesos
complejos. En específico se explorarán las sinergias con los algoritmos genéticos, el
algoritmo de temple simulado, el algoritmo de Nelder-Mead y el algoritmo de entropía
cruzada, de los cuales se realizará una introducción en el siguiente capítulo.
125
Métodos clásicos y de soft-computing en la optimización de procesos complejos.
Consideremos que el controlador se define mediante una ley de control que relaciona su
variable de entrada (el error e ) con la acción de control ∆f y, siendo v un vector de
parámetros, ésta se puede expresar en el dominio temporal t de la siguiente forma:
126 Capítulo 4
Raúl M. del Toro
I = g ( e, t ) t
tF
(4.2)
0
El objetivo de control que se plantee puede cuantificarse mediante esta cifra de mérito.
Existen varias cifras de mérito o funciones de coste de un amplio uso en el diseño de sistemas
de control de procesos físicos, que son empleadas tanto en su etapa de diseño como de
evaluación. Entre las funciones de coste, más ampliamente utilizadas, se encuentra la integral
del error absoluto ponderado en el tiempo (ITAE) (4.3), que permite obtener una respuesta
transitoria del sistema de control con un pequeño tiempo de establecimiento y sin
oscilaciones.
Otros criterios también definidos son la integral del error cuadrático ponderado en el
tiempo (ITSE) (4.4), la integral del error absoluto (IAE) (4.5) y la integral del error cuadrático
ponderado por el tiempo al cuadrado (IT2SE o ISTSE) (4.6), entre otros. En específico, el
criterio ITAE penaliza las respuestas transitorias de larga duración y es muy sensible ante
variaciones de parámetros de los sistemas de control en comparación con otros criterios. Por
otra parte, el diseño basado en el criterio ITAE, por lo general, causa menores oscilaciones en
la variable controlada del proceso que los criterios IAE e ITSE.
tF
ITAE = ∫ e ( t ) ⋅ t dt (4.3)
t0
tF
ITSE = ∫ e 2 ( t ) ⋅ t dt (4.4)
t0
tF
IAE = ∫ e ( t ) dt (4.5)
t0
tF
ISTSE = ∫ e2 ( t ) ⋅ t 2 dt (4.6)
t0
Una vez definida una cifra de mérito, el procedimiento de optimización puede realizarse,
tanto en tiempo real (on-line), es decir, empleando datos reales de FR , F y ∆f como en
tiempo de simulación (off-line) utilizando algún modelo del proceso físico que relacione su
entrada y su salida:
F (t ) = G ( ∆f , t ) (4.7)
I = g ( e, t ) t = g ( FR (t ) − F (t ) ) t
tF tF
0 0
= g ( FR (t ) − G ( ∆f , t ) )
tF
(4.8)
t0
(
= g FR (t ) − G ( c ( e(t ), v ) ) )
tF
t0
tF
ITAE = ∫ FR (t ) − G ( c ( e(t ), v ) ) ⋅ t dt (4.9)
t0
128 Capítulo 4
Raúl M. del Toro
1
xc =
n
∑ xj (4.11)
j ≠h
A continuación se construye un nuevo símplex a partir del actual. Para ello se reemplaza
solamente el peor vértice xh por un mejor punto, por medio de transformaciones de reflexión,
expansión o contracción con respecto al lado opuesto a xh .
xr = xc + α ( xc − xh ) (4.12)
d rc = α d hc (4.13)
xe = xc + γ ( xr − xc ) (4.14)
130 Capítulo 4
Raúl M. del Toro
xt = xc + β ( xr − xc ) , β = dtc / d rc (4.15)
xt = xc + β ( xh − xc ) , β = dtc / d hc (4.16)
x j = xl + δ ( x j − xl ), j = 0,..., n ; j ≠ l (4.17)
En la práctica, estos coeficientes suelen ajustarse a unos valores estándares que, a la larga,
son los más utilizados en las aplicaciones del algoritmo:
1 1
α = 1, β = , γ = 2, δ = (4.19)
2 2
Entre las ventajas principales que ofrece el algoritmo de Nelder-Mead, además de ser un
método fácil de entender y aplicar, permite obtener un resultado del problema a optimizar con
132 Capítulo 4
Raúl M. del Toro
{ }
x = xi ; i = 1,..., D , refiriéndose a T como la temperatura.
e− Ek +1 / T 1
h( ∆ E ) = − Ek +1 / T − Ek / T
= ∆E / T
≈ e −∆E / T ; ∆E = Ek +1 − Ek (4.20)
e +e 1+ e
S = −∑ p ( x) ln ( p ( x ) / p( x ) ) (4.22)
x
1 − H ( x ) /T
G ( x) = e ; Z = ∑ e− H ( x ) /T (4.23)
Z x
G ( xk )h ( ∆E ( x) ) = G ( xk +1 ) (4.24)
Lo anterior es suficiente para establecer, al menos en teoría, que todos los estados del
sistema pueden ser muestreados. Sin embargo, la planificación del templado interrumpe el
134 Capítulo 4
Raúl M. del Toro
equilibrio cada vez que la temperatura cambia y de esta forma, en el mejor de los casos, esto
debe realizarse cuidadosamente y de forma gradual.
Un aspecto importante en el algoritmo de temple simulado es la selección del rango de los
parámetros a ser buscados. En la práctica, el cálculo de sistemas continuos requiere de su
discretización por lo que puede asumirse que, sin mucha pérdida de información, el espacio
de trabajo puede ser discretizado. Existen ciertas restricciones adicionales que son necesarias
cuando se utilizan funciones de coste con valores integrales. En varias aplicaciones se suelen
utilizar técnicas para reducir el espacio inicial a medida que progresa la búsqueda. Por
ejemplo, para muchos sistemas físicos, pertenecientes a la clase de procesos Gauss-Markov,
es posible escoger como función de densidad de probabilidad del espacio de estado de los
parámetros la siguiente distribución [224] de Boltzmann:
T (k ) = T0 / ln k (4.26)
T (k ) = T0 e(α −1) k
; 0 < α <1 (4.27)
T (k + 1) ≈ α T (k )
donde α es una constante de reducción cuyos valores típicos se ajustan entre 0.8 y 0.99. Por
otra parte, para una variación rápida de la temperatura T [224], también se utiliza con
frecuencia
T (k ) = T0 / k (4.28)
aleatoria, mediante una distribución uniforme, evitando de esta manera que ocurra una
convergencia prematura hacia un mínimo local y se pueda continuar la búsqueda del mínimo
global.
136 Capítulo 4
Raúl M. del Toro
138 Capítulo 4
Raúl M. del Toro
suelen utilizarse otro tipos de representaciones, como son las representaciones mediante
valores enteros y valores reales (ver Tabla 4.1).
Aleatorio: se genera un vector binario Sea el siguiente vector binario generado de forma
de genes a seleccionar de cada padre. aleatoria, donde 1 se corresponde con los genes
del padre P1 y 0 los genes del padre P2 :
1 0 1 1 0 1
Sean además, los vectores de genes de los padres:
P1 : a b c d e f
P2 : 1 2 3 4 5 6
Al combinar los genes de los padres, se obtiene el
siguiente vector de genes del individuo hijo:
H1 : a 1 c d 5 f
Población inicial: Esta puede generarse de varias formas, una de la más utilizada es generar
los individuos o sus genes, empleando una distribución uniforme. Goldberg ofrece una
colección más extensa de funciones o formas para generar la población inicial [36].
Evaluación de la población: La función de evaluación ha de medir la adaptación de cada uno
de los individuos dependiendo, fundamentalmente, del tipo de representación elegido ya que
esta función evalúa la información contenida en los genes. Si se emplea el algoritmo para
resolver problemas de optimización, la función a optimizar se convierte en la función de
evaluación.
Operación de selección: Mediante esta operación se transmiten aquellas características de las
soluciones que son consideradas valiosas a lo largo de generaciones. Los individuos mejor
evaluados tienen más probabilidades de reproducirse. Algunas de las funciones más utilizadas
para esta etapa son la selección por ruleta [36], por torneo [227] y la estocástica uniforme,
entre otras.
Cruzamiento: En esta etapa, la información de la población se combina para crear mejores
individuos. Algunos de los operadores de mayor utilización (ver Tabla 4.2) resultan el
cruzamiento aleatorio, el cruce de un punto y entre 2 puntos, entre otros.
Mutación: Mediante la mutación el algoritmo realiza pequeños cambios aleatorios en los
genes de los individuos de la población actual y así fomentar la variabilidad dentro de la
población. Entre los operadores más utilizados está la función Gaussiana que adiciona
números aleatorios, tomados de una distribución normal con media igual a 0, a cada gen del
individuo. Otro operador es la función de mutación uniforme donde, primeramente, se
selecciona una parte de los genes de un individuo asignándole a cada gen una razón de
probabilidad de ser mutado y, seguidamente, el algoritmo reemplaza cada gen seleccionado
por un número aleatorio que esté comprendido dentro del rango del gen dado.
Reemplazo de la población: En esta etapa, los nuevos individuos creados durante la etapa de
cruzamiento reemplazarán a parte de la población existente para crear la nueva generación. El
reemplazo puede realizarse mediante diferentes métodos. Algunos de los más comúnmente
empleados son el reemplazo de un número de individuos que se corresponde con los peores
evaluados de la población, cada nuevo individuo reemplaza al peor evaluado y cada nuevo
individuo sustituye a cualquier individuo elegido de forma aleatoria.
Condición de parada del algoritmo: Existen varias condiciones de paradas del algoritmo, que
pueden ser, por ejemplo, alcanzar un número máximo de generaciones, un tiempo máximo de
ejecución del algoritmo o, en el caso de optimización o ajuste de funciones, cuando se alcance
un error igual o menor a un error mínimo límite, entre otras condiciones de parada.
140 Capítulo 4
Raúl M. del Toro
Todas las etapas anteriormente descritas del algoritmo hacen que su robustez sea superior
comparada con otras técnicas clásicas de búsqueda y optimización. Se resaltan como
características principales al empleo de una población de puntos del espacio de búsqueda y no
de un único punto y el uso de mecanismos de elección aleatoria como herramienta de
búsqueda de mejores regiones del espacio. Además, se debe destacar su característica de
blindaje, ya que la única información que toman del entorno es el valor de la función objetivo
del conjunto de parámetros, es decir, no emplean información específica del problema [225].
Información más exhaustiva acerca de las diferentes etapas del algoritmo, así como los
diferentes operadores que se utilizan en cada una de ellas además de propuestas de
variaciones del algoritmo o hibridaciones con otras técnicas de optimización, pueden
encontrarse en los trabajos de Holland [216, 228], Goldberg [36, 227], Grefenstette [217,
229], entre otros.
γ * = φ ( x* ) = min φ ( x) (4.29)
x∈χ
El método CE provee una metodología para crear una secuencia x0 , x1 ,..., xN y de niveles
γ 0 , γ 1 ,..., γ N , tal que γ converge a γ * y x converge a x* .
{ }
EV = x ∈ χ φ ( x ) ≥ γ , γ ∈ ℜ+ .
1 if φ ( x ) ≤ γ
IV ( x, γ ) = I{φ ( xi ) >γ } = (4.30)
0 if φ ( x ) > γ
l ( γ ) = Pν (φ ( x ) ≥ γ )
(4.31)
= Eν ⋅ Iν ( x,ν )
N
lˆ = 1 N ⋅ ∑ I{φ ( xi ) >γ } ⋅W ( xi ) (4.32)
i =1
I φ x ≥γ ⋅ g ( x,ν )
h* ( x ) = { ( i ) } (4.33)
l
D ( g , h ) = ∫ g ( x ) ⋅ ln g ( x ) dx − ∫ g ( x ) ⋅ ln h ( x ) dx (4.34)
142 Capítulo 4
Raúl M. del Toro
cual implica:
ν v
(
max D (ν ) = max E p I{φ ( xi )≥γ } ⋅ ln g ( x,ν ) ) (4.35)
1 N p (x )
max Dˆ (ν ) = max ∑ I{φ ( x ≥γ )} ⋅ x i ⋅ ln g ( xi ,ν ) (4.36)
v N i =1 i
h ( xi )
φ ( x* ) = γ * = max φ ( x ) (4.37)
x∈χ
1. Inicializar ν t = ν 0
3. Actualizarν t mediante:
N elite
⋅ ln g ( xit ,ν t )
1
ν t +1 = arg min elite
ν N
∑ I{φ (
i =1
x ti ≥γ t )}
(4.38)
4. Repetir desde el paso 2 hasta que se produzca una convergencia o se cumpla algún
criterio de finalización del algoritmo.
N elite
xij
µɶ tj = ∑i =1
N elite
(4.39)
(x − µtj )
2
N elite
σɶ =
2
tj ∑i =1
ij
N elite
(4.40)
µˆ t = α ⋅ µɶ t + (1 − α ) ⋅ µˆt −1 (4.41)
σˆt = β t ⋅ σɶ t + (1 − β t ) ⋅ σˆ t −1 (4.42)
( )
q
βt = β − β ⋅ 1 − 1 t (4.43)
144 Capítulo 4
Raúl M. del Toro
donde 0.4 ≤ α ≤ 0.9 son valores típicos, β es una constante de suavizado (valores
recomendados entre 0.6 y 0.9) y q es un entero con valores típicos entre 2 y 7.
4.6 Conclusiones
En este capítulo se ha realizado un estudio de métodos de optimización estocásticos y
libres de gradiente, de inspiración geométrica, física y biológica, para su aplicación al ajuste
óptimo de sistemas de control de procesos complejos. Los métodos de optimización
permitirán desarrollar procedimientos de ajuste de parámetros de control. Se trata de las
estrategias de control que han sido diseñadas por técnicas de soft-computing y clásicas como
las descritas en el Capítulo 3 y que utilizarán representaciones computacionales eficientes
(modelos) descritos en el Capítulo 2. El fin último es, como veremos más adelante, utilizar
estos métodos de optimización en arquitecturas jerárquicas de control inteligente de procesos
con una relativa alta complejidad.
En el capítulo siguiente se propone un método de diseño de las estrategias de control para
ser utilizadas en el control de procesos de fabricación de alta complejidad, en específico, el
proceso de taladrado de alto rendimiento. Entre las etapas de diseño, tiene un papel muy
importante el uso de los algoritmos de optimización introducidos en este capítulo, ya que
permitirán obtener un alto desempeño de las estrategias de cara al cumplimiento de los
requisitos de control.
147
Métodos clásicos y de soft-computing en la optimización de procesos complejos.
f (s) K
G ( s) = = K p + i + Kd ⋅ s (5.1)
PID E (s) s
148 Capítulo 5
Raúl M. del Toro
Figura 5.2. Gráfica del lugar de las raíces para una respuesta aperiódica.
150 Capítulo 5
Raúl M. del Toro
a) b)
f p [k ] = K p e[k ]
fi [k ] = fi [k − 1] + ( K pTs / Ti ) e[k ]
(5.3)
f d [k ] = ( K pTd / Ts ) (e[k ] − e[k − 1])
f [k ] = f p [k ] + fi [k ] + f d [k ]
1−τ
F ( z) = (5.4)
z −τ
Con vista a probar el desempeño del controlador PID, se han realizado experimentos de
control del proceso de TAR a lazo cerrado, empleándose los siguientes parámetros:
Referencia en la fuerza: SP = 1000 N
Tiempo de muestreo del lazo: Ts = 0.069 s
Coeficiente del filtro: τ = 0.95
Velocidad de avance nominal: 100 mm/min
Velocidad de giro nominal: 870 rpm
Umbral de activación del controlador: 100 N
Ganancia de la acción proporcional: K p = 0.01086
Ganancia de la acción integral: K i = K p / Ti = 0.01348
Ganancia de la acción derivativa: K d = K pTd = 0.00219
Límites de saturación de la velocidad de avance: [10,130]mm/min
Como puede observarse, el controlador PID utiliza una versión modificada de los
parámetros ajustados inicialmente mediante Z-N. Con el objetivo de mejorar la respuesta
transitoria del sistema de control, la ganancia K p se ha ajustado a un 10% de su valor
calculado inicialmente mediante Z-N. La Figura 5.5 muestra la respuesta en tiempo real del
sistema de control, tanto la fuerza de corte controlada, como la velocidad de avance como
acción de control. En la figura se muestran tres condiciones de operación diferentes: sin el
152 Capítulo 5
Raúl M. del Toro
control PID activado, es decir, a lazo abierto, control a lazo cerrado sin filtrado en la
referencia de la fuerza y también a lazo cerrado pero con el filtro pasa bajo en la referencia.
Como puede observarse en la Figura 5.5, el filtrado del valor de referencia, mediante un
filtro discreto de primer orden, mejora el comportamiento dinámico de la fuerza de corte,
obteniéndose una respuesta del sistema de control sin sobrepaso, contrario a la respuesta
dinámica sin filtrado, en la que se obtiene un máximo sobrepaso de un 17.7 %. El sistema de
control es capaz de mantener la fuerza con muy poca variación alrededor de su valor de
referencia, durante toda la operación de taladrado, hecho que no se produce en el caso sin
control PID, en el cual la fuerza de corte se incrementa en el tiempo, a medida que aumenta la
profundidad durante el corte.
a) b)
Criterios de error
Condición de operación
IAE ITAE ITSE
(ITSE). En la Tabla 5.1 se muestran los índices de error calculados para cada una de las
condiciones de operación. Los valores de los índices relacionados en la tabla, reflejan la
incuestionable superioridad del desempeño cuando se utiliza el sistema de control a lazo
cerrado, encontrándose los índices entre un 32 y un 82 %, es decir entre 1.2 y 2.9 veces
menor, que los índices de comportamiento del sistema a lazo abierto sin el controlador PID.
Por otra parte, el tiempo de la operación de mecanizado es menor cuando se utiliza el
control a lazo cerrado. Al taladrar sin control activo, la velocidad de avance se programa en el
CNC a un valor constante, el cual se escoge de forma conservativa, de tal manera que no se
produzcan roturas de la herramienta durante el mecanizado. En el experimento realizado, el
valor máximo de la fuerza de corte debe ser 1000 N, por lo que para que no se exceda este
valor, la velocidad de avance programada debería ser entre 86 y 87 mm/min. De utilizar esta
condición nominal, el tiempo de la operación de taladrado sería aproximadamente 1.3 %
mayor que el tiempo de la operación con control a lazo cerrado, lo cual representa una
ganancia importante en el tiempo de producción. Más aún, una ventaja tecnológica adicional
como consecuencia de mantener prácticamente constante la fuerza de corte durante la
operación de taladrado, es un aumento de la vida útil de la herramienta de corte, tal y como ha
sido reportado anteriormente en la literatura por Oh et ál. [161].
154 Capítulo 5
Raúl M. del Toro
a) b)
f F F f
Una vez definida la estructura de los modelos, ambos sistemas neuro-borroso ANFIS, se
entrenan para que aprendan el comportamiento a régimen transitorio del proceso de TAR
controlado por alguna estrategia de control diseñada previamente. Para este caso, se
emplearán datos experimentales entrada-salida de operaciones de taladrado controladas
156 Capítulo 5
Raúl M. del Toro
alcanza el valor de salida deseado) y tamaño del paso de 0.01 (el aumento de este valor no
produce una mejora significativa de la salida y si un mayor cómputo de operaciones). El
tiempo de entrenamiento en los modelos, tanto para el modelo directo como para el inverso,
fue de 0.14 segundos. Estos tiempos se obtuvieron en un ordenador con procesador Intel
Core2 CPU 6400 - 2.13 GHz con sistema operativo Windows XP Profesional que no se
encontraba conectado en línea con el proceso de TAR.
a) MODELO DIRECTO
Base de reglas:
Regla 1: SI (f es “baja”) ENTONCES (F es “baja”) (1)
Regla 2: SI (f es “alta”) ENTONCES (F es “alta”) (1)
a b c
Veloc. avance baja: 3.475 2.001 93.59
Veloc. avance alta: 3.473 1.998 100.5
b) MODELO INVERSO
Base de reglas:
Regla 1: SI (F es “baja”) ENTONCES (f es “baja”) (1)
Regla 2: SI (F es “alta”) ENTONCES (f es “alta”) (1)
a b c
Fuerza media baja: 520 1.983 2.70·10-4
Fuerza media alta: 520 1.964 1040
a)
b)
1
∑
N
RMSE = i =1
[(Yi − Fi ) / Yi ]2 (5.5)
N
donde Fi es valor obtenido por el modelo, Yi es el valor deseado o real y N es el número total
de muestras del experimento.
Tanto en el caso del modelo directo como en el caso del modelo inverso, se eligieron
modelos que, a priori, tienen un error medio cuadrático mayor que otros modelos (con un
158 Capítulo 5
Raúl M. del Toro
Simple-entrada / Simple-entrada /
Entrada-Salida
Simple-salida Simple-salida
Tipo de función de
Campana generalizada Campana generalizada
pertenencia
Número de funciones de
2 2
pertenencia
Número de reglas 2 2
RMSE - datos de
1.94 % 0.39 %
entrenamiento
RMSE - datos de
5.18 % 0.68 %
validación
F ( z ) = (1 − k1 ) /( z − k2 ) (5.6)
Figura 5.10. Aplicación de simulación de la estrategia de control por modelo interno y modelos neuro-
borrosos.
siendo ω = 7.61 rad/s . Esta frecuencia se corresponde con uno de los polos del modelo de
tercer orden que representa de manera aproximada el proceso de taladrado. La amplitud de la
perturbación se estima en A = 10 para emular un ruido en la señal del 1 % de su nivel a estado
estacionario (nivel de ruido bastante común en procesos de fabricación).
Se ha simulado una operación de taladrado de 15 mm (parámetro que se utilizará durante
la validación del sistema de control en tiempo real), estimándose, previo a la simulación, una
duración aproximada de la operación de 10 segundos y fijándose este tiempo como tiempo de
simulación. Por otra parte, se ajustó como valor de referencia de la fuerza de corte el valor
FR = 1000 N . Cabe recordar que el objetivo de control es mantener la fuerza de corte
constante alrededor de su valor de referencia, una respuesta transitoria sin sobrepaso
(amortiguamiento crítico) y bajo tiempo de establecimiento.
La Figura 5.11 muestra la respuesta del sistema de control por modelo interno del proceso
de taladrado, en presencia de perturbación (5.7). Si analizamos la respuesta transitoria del
sistema a lazo cerrado se puede comprobar que el máximo sobrepaso es del 2.36 %, y el
tiempo de establecimiento es muy reducido (menor a 2 segundos, despreciando el retardo
inicial). Por otra parte a partir de la respuesta temporal del error, se ha determinado el valor de
160 Capítulo 5
Raúl M. del Toro
la integral del error absoluto ponderado en el tiempo, ITAE = 0.54 , la integral del error
absoluto, IAE = 0.96 y la integral del error cuadrático ponderado en el tiempo, ITSE = 0.32 .
a) b)
Figura 5.11. Simulación de la respuesta del sistema de control ante un escalón en la referencia
a) Fuerza de corte y b) acción de control.
Por otra parte, se analizó el comportamiento del sistema de control en presencia de retardo
de tiempo variable en el bus campo por el que se envía la acción de control al proceso de
taladrado, realizándose simulaciones con diferentes retardos comprendidos entre [0.025, 0.75]
segundos. En la Figura 5.12 se describe el comportamiento del índice de desempeño ITAE
según varía el retardo de tiempo del proceso. El valor mínimo del criterio ITAE ha sido
alcanzado por el retardo mínimo y viceversa.
a)
b)
Figura 5.13. Comportamiento del sistema de control CMI-ANFIS y otras estrategias de control
a) Fuerza de corte y b) velocidad de avance de referencia.
162 Capítulo 5
Raúl M. del Toro
Tabla 5.3. Estudio comparativo entre el controlador CMI-ANFIS y otras estrategias de control analizadas.
La Tabla 5.3 muestra los índices de error calculados a partir de la respuesta temporal de
las estrategias de control implementadas en tiempo real. Para el cálculo se ha utilizado como
tiempo inicial ( t0 = 0 ) el instante cuando se produce el contacto entre la herramienta y el
material de trabajo. La estrategia CMI-ANFIS tiene un mejor desempeño que el resto de
estrategias bajo los criterios ITAE, IAE e ITSE. Esta misma conclusión puede corroborarse a
partir de las gráficas de la Figura 5.14, donde se describe, para cada una de las estrategias de
control, el error absoluto ponderado en el tiempo e(t ) t (Figura 5.14.a) y el error absoluto
e(t ) (Figura 5.14.b) durante toda la operación de taladrado. El error está calculado relativo al
valor de referencia de la fuerza de corte. En cada una de las gráficas aparece sombreada el
área debajo de la curva correspondiente al sistema de control CMI-ANFIS y que representa el
valor de los índices ITAE e IAE, respectivamente.
a) b)
Figura 5.14. Comportamiento del error relativo para las estrategias de control
a) Error absoluto ponderado en el tiempo y b) error absoluto.
164 Capítulo 5
Raúl M. del Toro
e (k ) = F (k ) − F (k ) (5.8)
F r
∆e (k ) = e (k ) − e (k − 1) (5.9)
F F F
[e ec] = K ⋅ e K ⋅ ∆e ; u = ∆f (5.10)
e F ce F
166 Capítulo 5
Raúl M. del Toro
NM NB NB NB NM NS ZE PS
NS NB NB NM NS ZE PS PM
ZE NB NM NS ZE PS PM PB
PS NM NS ZE PS PM PB PB
PM NS ZE PS PM PB PB PB
PB ZE PS PM PB PB PB PB
49
µ (e, ec, u ) = max prod µ (e), µ (ec ), µ (u ) (5.11)
R i = 1 e ec u
i i i
∑ µ (ui ) ⋅ ui
R
u= i (5.12)
∑ µ ( ui )
R
i
f (k ) = f (k − 1) + GC ⋅ u (k ) (5.13)
168 Capítulo 5
Raúl M. del Toro
Figura 5.17. Sistema de control distribuido basado en un regulador borroso para el control del proceso de
TAR
Para el ajuste del control se ha utilizado como función objetivo a optimizar, el criterio o
función de coste ITAE, que permite obtener una respuesta transitoria del sistema de control
con pequeño tiempo de establecimiento y sin sobrepasos. Otros criterios a tener en cuenta
serán la integral del error cuadrático ponderado en el tiempo (ITSE), la integral del error
absoluto (IAE) y la integral del error cuadrático ponderado por el tiempo al cuadrado (IT2SE o
ISTSE). En específico, el criterio ITAE penaliza las respuestas transitorias de larga duración y
es muy sensible ante variaciones de parámetros en comparación con otros criterios. Por otra
parte, el criterio ITAE por lo general produce menores sobrepasos y oscilaciones que los
criterios IAE e ITSE.
El comportamiento del sistema de control ha sido examinado mediante simulación en
Matlab/Simulink empleando, como aplicación, el diagrama ilustrado en la Figura 5.18. El
criterio ITAE describiría el comportamiento del sistema ante perturbaciones externas. En este
estudio se ha considerado como perturbación un cambio de tipo escalón en el valor de
referencia de la fuerza de corte, y el criterio ITAE es una buena opción para evaluar cuan bien
es capaz el sistema de control de seguir los cambios en la referencia.
T
K OPT = [ K e K ce ]OPT = arg min ∫ e ( t ) t dt (5.14)
F
0
Valores altos de la ganancia K e del controlador borroso significan una buena capacidad
de respuesta del sistema de control, sintetizada en un error muy pequeño a estado estacionario
y bajo tiempo de establecimiento. Sin embargo esto puede conllevar un empeoramiento de la
estabilidad del sistema y provocar altos niveles de sobrepasos en su respuesta temporal. Por
otra parte, una convergencia más rápida de la variable controlada hacia su valor de referencia,
se obtiene con la ganancia K ce pequeña. Existen varios métodos para ajustar controladores
borrosos, pero no es objetivo de este trabajo hacer un estudio de estos en profundidad.
170 Capítulo 5
Raúl M. del Toro
sustituir el peor punto del símplex con otro punto reflejado en el resto de puntos, a través del
cálculo del centro de gravedad de estos M-1 puntos. Resulta típico el empleo de un símplex de
tamaño constante y pequeño, que permita desplazarse mediante gradientes locales hasta los
máximos o mínimos locales [221].
La función de coste fue evaluada mediante simulación del sistema de control, donde el
proceso de TAR se representa por su modelo dinámico de tercer orden. Al ajustar los
parámetros del controlador a los siguientes valores iniciales: [ Ke Kce ]INIT = [1.0,1.0] y
Gc = 0.5 , se obtiene un máximo sobrepaso de 4.090 % y un índice ITAE igual a 0.709.
Los resultados de la simulación del sistema de control son mostrados en la Figura 5.19. En
ella se puede observar tanto el comportamiento de la fuerza de corte ante un cambio de tipo
escalón en su valor de referencia, como en la señal de control representada por la velocidad de
avance de referencia. Como cabría de esperar, se consigue una mejora clara en la respuesta
transitoria del sistema tras el procedimiento de optimización.
a) b)
a)
b)
Figura 5.20. Comportamiento del sistema de control según varía el tiempo de retardo.
a) Fuerza de corte y b) variaciones en la velocidad de avance.
Sin embargo, estos modelos tienen una validez limitada, debido al comportamiento no
determinista de algunas redes de comunicación. Para estudiar la influencia del retardo de
tiempo en el desempeño del controlador, se ha incluido en la simulación un simple generador
de retardos. Se asume conocido el tiempo de retardo máximo y, con el objetivo de mantener
172 Capítulo 5
Raúl M. del Toro
Tabla 5.5. Comportamientos extremos del sistema de control para el tiempo de retardo variable.
Figura 5.21. Comportamiento de los índices de desempeño ITAE e IT2SE según varía el tiempo de retardo.
a)
b)
Figura 5.22. Comportamiento del proceso según las estrategias de control PID y borrosas.
a) Fuerza de corte y b) velocidad de avance de referencia.
174 Capítulo 5
Raúl M. del Toro
temporal del proceso libre de oscilaciones, a pesar de las fuerzas de perturbación que surgen
con el incremento de la profundidad de corte.
Para corroborar el buen desempeño del controlador, se ha realizado un estudio
comparativo con otras estrategias de control, diseñadas para el mismo proceso y bajo las
mismas condiciones de operación, cuyas respuestas temporales son mostradas también en la
Figura 5.22. En un trabajo previo, Haber et ál. [236] presentan el diseño e implementación de
un controlador borroso de 9 reglas en su base de conocimiento, funciones de pertenencia de
tipo trapezoidal y ajustado mediante el método iterativo de prueba y error en pruebas de
simulación. En el estudio comparativo se ha tenido en cuenta, además, el controlador clásico
PID con filtrado de referencia y optimizado mediante el método de Ziegler-Nichols,
introducido en secciones anteriores.
La Tabla 5.6 muestra los índices de error (relativo al valor de referencia de la fuerza)
calculados a partir de la respuesta temporal de las estrategias de control implementadas en
tiempo real. A primera vista, el comportamiento del controlador borroso de 49 reglas parece
muy similar al resto. Sin embargo, índices como el ITAE, el IT2SE y el IAE, indican una
mejora en la respuesta dinámica del proceso empleando la estrategia de control borrosa de 49
reglas en su base de conocimientos. Con vistas a realizar una mejor evaluación gráfica del
desempeño temporal de las estrategias de control, la Figura 5.23 ilustra el comportamiento del
error absoluto ponderado en el tiempo ( e(t ) t ) para cada una ellas. Se ha rellenado el área
debajo de la curva correspondiente a la estrategia de control borrosa de 49 reglas. El valor de
esta área representa el índice ITAE y como puede comprobarse es menor al área debajo del
resto de curvas.
Figura 5.23. Comportamiento del error absoluto ponderado en el tiempo de las estrategias de control.
176 Capítulo 5
Raúl M. del Toro
[251]. Sin embargo, si hubiera disponible una aproximación lineal del proceso, la cual
caracteriza sus componentes dinámicas principales en al menos una vecindad del punto de
operación, entonces resultaría factible realizar una prueba de estabilidad local y asintótica. La
realimentación del sistema de control, que enlaza la función de transferencia lineal del
proceso con el controlador borroso no lineal, da lugar a un lazo de control no lineal estándar
que puede ser analizado mediante el criterio de Estabilidad Absoluta con un paradigma de
dominio finito, por ejemplo, con el criterio del círculo [252].
a) b)
178 Capítulo 5
Raúl M. del Toro
el valor mínimo alcanzado, sin embargo el error es menor del 2% en el resto de las
ejecuciones. La Figura 5.25. b) muestra los valores correspondientes de los parámetros del
controlador durante cada ejecución del algoritmo. Por otra parte, el tiempo de medio de
cómputo del algoritmo de temple simulado es de 925 s, el cual se corresponde a su ejecución
en la herramienta computacional Matlab, versión 7.1, en un ordenador con un procesador
Pentium IV a 3.2 GHz de velocidad y 1.99GB de memoria RAM.
a)
b)
a)
b)
Figura 5.26. Simulación del sistema de control borroso optimizado mediante SA.
a) Comportamiento de la fuerza de corte y b) velocidad de avance de referencia.
Con la finalidad de evaluar la influencia de la variación del tiempo de retardo del proceso
en el desempeño de controlador, se ha llevado a cabo un estudio mediante simulación. Para
simplificar, se ha asumido durante la simulación que el tiempo de retardo del sistema de
control distribuido L, toma valores uniformemente distribuidos entre 0 y 0.8 s, ejecutándose
100 pruebas de simulación empleando el sistema de control borroso controlador optimizado.
La Figura 5.27 describe el comportamiento de la fuerza de corte para cada tiempo de retardo.
180 Capítulo 5
Raúl M. del Toro
Figura 5.27. Comportamiento de la fuerza de corte según la variación del tiempo de retardo.
Figura 5.28. Comportamiento del criterio ITAE según varía el tiempo de retardo.
El comportamiento del índice de desempeño ITAE, según varía el retardo de tiempo del
proceso, se muestra en la Figura 5.28. El valor mínimo del criterio ITAE se alcanza por el
retardo mínimo y viceversa. Sin embargo, se consigue un buen comportamiento del índice
ITAE aun suponiendo que el retardo máximo es mucho mayor que el utilizado para el modelo
del proceso de taladrado. El peor de los casos entre todas las pruebas de simulación, se
corresponde al máximo tiempo de retardo de 0.7954s, el cual provoca un máximo sobrepaso
de 4.95% y un máximo ITAE igual a 1.1476.
a) b)
Figura 5.29. Respuesta temporal del sistema de control optimizado mediante SA durante el taladrado del
material GGG40.
182 Capítulo 5
Raúl M. del Toro
a) b)
Figura 5.30. Comportamiento del sistema de control optimizado durante una operación de taladrado del
material 17-4PH
Figura 5.31. Diagrama del sistema de control borroso y optimización mediante algoritmos genéticos.
184 Capítulo 5
Raúl M. del Toro
a)
b)
El sistema de control fue simulado en Matlab/Simulink para las ganancias del controlador
ajustadas. Los resultados de la simulación son mostrados en la Figura 5.33, en ella se puede
observar el comportamiento de la fuerza de corte ante un cambio de tipo de escalón en el
valor de referencia y la velocidad de avance, como señal de control, la cual se expresa en
valor porcentual, referido a su valor nominal f 0 = 100 mm/min . Durante la simulación se
obtuvo un máximo sobrepaso igual a 0.5659%.
Con el fin de analizar las capacidades de rechazo ante perturbaciones del sistema de
control optimizado, se ha considerado la presencia de ruido aditivo durante la simulación (ver
diagrama en la Figura 5.34), el cual puede representar la influencia de dinámicas no
a) b)
c)
Figura 5.34. Aplicación de simulación del sistema de control con ruido aditivo.
186 Capítulo 5
Raúl M. del Toro
T
ITAE = ∫ e ( t ) ⋅ tdt (5.18)
0
T
ISTSE = ∫ e 2 ( t ) ⋅ t 2 dt (5.19)
0
Figura 5.35. Diagrama del sistema de control borroso y el método de optimización de entropía cruzada.
La Figura 5.35 muestra el diagrama de bloques del sistema de control borroso. El ajuste
óptimo de los parámetros del controlador se realiza mediante la minimización del criterio
ITAE. Tras ejecutar el método de entropía cruzada se obtiene como resultado los parámetros
del controlador óptimos. En el diagrama de bloques, las líneas discontinuas representan las
“señales” más representativas utilizadas para la optimización.
188 Capítulo 5
Raúl M. del Toro
a) b)
La optimización del sistema de control borroso basada en entropía cruzada devuelve como
parámetros óptimos del controlador [ K e , K ce ]OPT = [ 0.0165, 0.0297] , obteniéndose un índice
de desempeño ITAE de 0.45 y un máximo sobrepaso de 0.46%. La Figura 5.37 muestra el
comportamiento del índice de desempeño ITAE y la variación de los parámetros del
controlador, en cada una de las iteraciones de la ejecución del algoritmo de entropía cruzada.
En la Figura 5.38 se muestran los resultados de la simulación del sistema de control
borroso, caracterizado por el comportamiento de la fuerza corte como variable controlada y la
velocidad de avance de referencia como acción de control. El sistema de control se ha
simulado con sus parámetros ajustados con los valores iniciales [ Ke Kce ]INI y los
obtenidos tras la optimización con entropía cruzada, [ Ke , Kce ]OPT . Se puede observar
claramente en la figura que la respuesta transitoria, el tiempo de establecimiento y el máximo
sobrepaso, del sistema de control, son mejores después de la optimización con el algoritmo de
entropía cruzada.
a) b)
a) b)
Con el objetivo de evaluar la influencia de la variación del tiempo de retardo del proceso
en el desempeño del controlador, se ha llevado a cabo un estudio mediante simulación. Por
motivo de simplicidad, se ha asumido durante la simulación, que el tiempo de retardo del
sistema de control distribuido L, toma valores uniformemente distribuidos entre 0 y 0.8 s. Se
han ejecutado 100 pruebas de simulación con los parámetros del controlador borroso
ajustados a sus valores óptimos.
La Figura 5.39 describe el comportamiento de los índices de desempeños ITAE e ISTSE,
según varía el retardo de tiempo del proceso. En el estudio se ha tenido en cuenta además el
comportamiento del criterio ISTSE, para evaluar la influencia de los retardos en la
sensibilidad de este índice de desempeño. El valor mínimo del criterio ITAE ha sido alcanzado
por el retardo mínimo y viceversa. Sin embargo, los valores que alcanza el índice ITAE son
buenos aun suponiendo que el retardo máximo es mucho mayor que el modelado
experimentalmente (ver Capítulo 2, sección 2.3). La curva que describe el comportamiento
190 Capítulo 5
Raúl M. del Toro
del criterio ISTSE muestra una pendiente muy cercana a cero hasta 0,4 segundos, lo cual
simboliza que el criterio ITAE muestra mejor sensibilidad que el criterio ISTSE para este caso
de estudio. Entre todas las pruebas de simulación, el peor caso se corresponde a un tiempo de
retardo de 0.795s, donde se obtiene ITAE = 1.125, ISTSE = 0.5695 y un máximo sobrepaso de
4.454%.
Figura 5.39. Comportamiento de los criterios de desempeño ITAE e ISTSE en presencia de retardos de
tiempo en el sistema de control.
respuesta transitoria del sistema de control, está respaldado por los valores de los índices de
desempeño, tal y como se muestra en la Tabla 5.7. Estos resultados experimentales
corroboran la idoneidad del método de CE para el ajuste óptimo de los sistemas de control
borrosos.
a) b)
Figura 5.40. Comportamiento del sistema de control durante el taladrado del material GGG40.
a) Fuerza de corte y b) velocidad de avance.
Tabla 5.7. Estudio comparativo de los resultados experimentales y de simulación del comportamiento del
controlador borroso optimizado mediante el método de entropía cruzada.
Parámetros del
Sobrepaso
controlador Procedimiento ITAE ITSE ISTSE
(%)
borroso
Sin compensación
Experimental 1.376 0.243 0.552 14.28
de fuerza
K e = 0.5
Simulación 0.636 0.306 0.185 5.67
K ce = 0.5
K e = 0.0165 Simulación 0.450 0.302 0.171 0.462
K ce = 0.0297 Experimental 0.449 0.131 0.673 2.841
192 Capítulo 5
Raúl M. del Toro
Tabla 5.8. Parámetros de las estrategias de optimización para un ajuste óptimo de un controlador borroso
de 49 reglas.
En la Tabla 5.8 se resumen las técnicas de optimización empleadas, así como los
parámetros ajustados para cada una de ellas. En cada una de las estrategias se ha optimizado
el controlador borroso de 49 reglas, minimizando como función de coste la integral del valor
absoluto del error ponderado en el tiempo (ITAE) y utilizando los parámetros del controlador
[ Ke Kce ] como coeficientes de ajuste.
La Tabla 5.9 describe los parámetros de los controladores ajustados por los cuatro
métodos de optimización, así como los índices de comportamiento ITAE y la integral del error
cuadrático ponderado por el tiempo al cuadrado (ISTSE), además del máximo sobrepaso, que
caracterizan la respuesta dinámica en el dominio temporal del sistema de control borroso. Con
el objetivo de evaluar el coste computacional y el nivel de complejidad de cada método, la
tabla también muestra el tiempo de cómputo y el número de parámetros de ajuste por cada
método de optimización.
La respuesta temporal de los cuatro sistemas de control sin presencia de ruido añadido en
la simulación es mostrada en la Figura 5.41. Las respuestas transitorias son bastante similares
ya que, en los cuatro casos, se trata de parámetros óptimos o sub-óptimos de los sistemas de
control borrosos.
Tabla 5.9. Estudio comparativo del desempeño del sistema de control borroso sin presencia de
perturbaciones.
Ke = 0.0173
Temple Simulado (SA) 0.4508 0.1712 0.6705 4 925.00
Kce = 0.0333
Ke = 0.0125
Nelder-Mead (NM) 0.4604 0.1748 0.0732 5 156.03
Kce = 0.0374
Algoritmos Genéticos Ke = 0.0169
0.4513 0.1714 0.5659 11 238.60
(GA) Kce = 0.0314
a) b)
c)
194 Capítulo 5
Raúl M. del Toro
Con el fin de analizar las capacidades de rechazo ante perturbaciones de los sistemas de
control, se ha considerado durante la simulación, la presencia de ruido aditivo, además de la
influencia de dinámicas no modeladas, de tal manera que el comportamiento dinámico del
sistema se aproxime al comportamiento en tiempo real del proceso controlado. Se asume la
siguiente función como perturbación aditiva en la salida del proceso:
Tabla 5.10. Estudio comparativo del desempeño del sistema de control borroso en presencia de
perturbaciones.
Nelder-Mead Ke = 0.0125
0.4855 0.1753 0.5427
(NM) Kce = 0.0374
con una buena relación entre el tiempo de cómputo y los parámetros de ajuste, el algoritmo de
entropía cruzada es más rápido, en cuanto tiempo de cómputo se refiere, que el método de
temple simulado y menor que algoritmos genéticos y, por otra parte, tiene menos parámetros
de ajuste que el método de algoritmos genéticos. Sin embargo, se ha de destacar que el tiempo
de cómputo es un criterio un tanto engañoso, porque no se ha tenido en cuenta para la
comparación, por ejemplo, el tiempo total empleado para obtener los parámetros de ajuste
más adecuados.
a) b)
5.8 Conclusiones
En este capítulo se han diseñado sistemas de control optimizados, basados en técnicas
clásicas y de soft-computing, para su aplicación en la regulación automática de procesos
complejos, comprobándose su desempeño en una plataforma experimental compuesta por un
centro de mecanizado de alta velocidad, en específico durante el control de operaciones de
taladrado de alto rendimiento.
El diseño de las estrategias de control se ha basado en una metodología en la que, a partir
del empleo de modelos computacionales de los procesos físicos, inicialmente se realiza su
ajuste mediante simulación y seguidamente se efectúa su implementación y evaluación en
tiempo real directamente en el proceso, a través de una plataforma experimental.
Se han diseñado estrategias de control clásicas lineales, como el algoritmo de control
proporcional, integral y derivativo, PID, de muy amplia utilización en el control de procesos
196 Capítulo 5
Raúl M. del Toro
198 Capítulo 5
CONCLUSIONES
Esta Tesis Doctoral basada en una metodología científica en la que interactúan métodos de
investigación teórica y experimental, ha presentado el diseño, desarrollo e implementación de
sistemas de control inteligente para la optimización de procesos complejos. El diseño está
basado en diferentes estrategias de hibridación de técnicas del soft-computing, como son la
computación borrosa, la neuro-computación y los métodos metaheurísticos. Los sistemas de
control desarrollados han sido aplicados al proceso de taladrado de alto rendimiento, un
proceso de fabricación que tiene un gran impacto en la producción de bienes y servicios.
Como paso inicial e indispensable de la metodología científica empleada, se ha realizado
una revisión del estado de la ciencia y la técnica en los paradigmas elegidos dentro del amplio
universo que ponen a nuestra disposición la Teoría General de Sistemas, la Teoría de Control,
las técnicas de Inteligencia Artificial y las Ciencias de la Computación. Se ha prestado
principal interés en los paradigmas con una fundamentación teórica sólida en el campo de la
Inteligencia Artificial (IA), con un enfoque desde la perspectiva de la Automática y la Teoría
de Sistemas. Es decir, el estudio se ha centrado en estrategias susceptibles de ser utilizadas en
el modelado, monitorización y control de procesos físicos. Paradigmas que aprovechan tanto
las ventajas de las técnicas básicas de Inteligencia Artificial, como las sinergias entre estas.
De modo que permiten el diseño de sistemas de control cualitativamente superiores y capaces
de ofrecer altas prestaciones en la optimización de procesos físicos complejos. Por este
motivo la revisión se ha enfocado hacia las técnicas del soft-computing, analizando su
evolución desde sus inicios hasta las tendencias actuales.
Tras el estudio, de la literatura revisada resalta, principalmente, la aplicabilidad de los
sistemas de inferencia borrosa como un marco computacional eficiente para representar el
conocimiento de expertos y, por tanto, viable en el control de procesos físicos complejos. Sin
embargo, las tendencias más actuales apuntan hacia su hibridación con otras técnicas, como la
199
Métodos clásicos y de soft-computing en la optimización de procesos complejos.
200 Conclusiones
Raúl M. del Toro
201
Métodos clásicos y de soft-computing en la optimización de procesos complejos.
202 Conclusiones
Raúl M. del Toro
203
Métodos clásicos y de soft-computing en la optimización de procesos complejos.
204 Conclusiones
Raúl M. del Toro
de parada, por ejemplo que el error de salida del modelo esté por debajo del valor deseado o
se sobrepase un número de iteraciones máximas.
Una vez diseñadas y aplicadas cada una de las estrategias de control, como parte de la
metodología científica de la Tesis Doctoral y también como se ha comentado en los pasos de
la metodología de diseño, todas las estrategias desarrolladas han sido validadas rigurosamente
en una plataforma experimental utilizada como soporte tecnológico. La validación de las
estrategias mediante su implementación en tiempo real, actuando directamente sobre
operaciones de taladrado de alto rendimiento, ha permitido corroborar los resultados
obtenidos en simulación, así como la efectividad de los métodos metaheurísticos
empleados en el ajuste de los sistemas de control inteligente.
Los estudios comparativos mostraron que, en el caso de la aplicación de las técnicas
metaheurísticas para el ajuste del sistema de control borroso de 49 reglas, las respuestas
temporales de los sistemas de control, obtenidas experimentalmente son similares entre
ellas y a las obtenidas en simulación. Esto se debe a que en todos los casos se han
alcanzado parámetros óptimos o sub-óptimos de los sistemas de control borrosos. A
pesar de su similitud, el desempeño en general del sistema de control optimizado
mediante el método de entropía cruzada, es ligeramente superior al resto de técnicas
utilizadas, además de ofrecer una muy buena relación tiempo de cómputo y total de
parámetros de ajuste del método requeridos.
El algoritmo de entropía cruzada ofrece un tiempo de cómputo menor que el método
de temple simulado y menor que los algoritmos genéticos, y por otra parte, tiene menos
parámetros de ajuste que el método de algoritmos genéticos. Sin embargo, ha de
destacarse que la comparativa basada en el tiempo de cómputo puede resultar un tanto
engañosa, porque no se ha tenido en cuenta en la comparación, por ejemplo, el tiempo total
empleado, durante todo el procedimiento de ajuste mediante cada estrategia de optimización,
para obtener el mejor conjunto sub-óptimo de los parámetros del controlador (factores de
escala).
Además, la validación de las estrategias mediante su implementación en tiempo real,
también ha dado la posibilidad de realizar estudios comparativos entre las estrategias basadas
en técnicas de soft-computing, las estrategias lineales y el proceso de taladrado a lazo de
control abierto. Los estudios comparativos se han basado en la cuantificación del
cumplimiento de los objetivos de control, a partir de cifras de méritos que caracterizan el
comportamiento dinámico de sistemas de control, como el máximo sobrepaso de la fuerza de
corte (variable controlada del sistema de control) y los valores medios de su error a estado
estacionario respecto a su valor deseado (ej., índices de desempeño ITAE, IAE, ITSE, entre
otros).
205
Métodos clásicos y de soft-computing en la optimización de procesos complejos.
Las cifras de mérito han corroborado que de las estrategias diseñadas, las basadas en
técnicas del soft-computing, presentan un desempeño superior a las estrategias lineales,
al menos controlando operaciones de taladrado de alto rendimiento a lazo cerrado en
fase de simulación y en tiempo real. Esto hace que las estrategias de control no lineales,
basadas en técnicas de Inteligencia Artificial, puedan ser consideradas más efectivas
para hacer frente a la complejidad intrínseca de procesos físicos, como es el caso de los
procesos de fabricación estudiados en esta Tesis Doctoral, y conducen, por tanto, a una
mejoría sustancial de su funcionamiento y eficiencia.
Las estrategias computacionales diseñadas e implementadas, tanto basadas en
técnicas de control lineal como en técnicas de soft-computing, aplicadas a la optimización
del proceso de taladrado de alto rendimiento, constituyen, desde el punto de vista
científico y técnico, otras de las aportaciones de este trabajo. De forma resumida, las
contribuciones científico-técnicas realizadas en esta Tesis Doctoral han sido:
1. El diseño e implementación, en un lenguaje de alto rendimiento, de un modelo
computacional que emula el comportamiento de variables representativas en
procesos de fabricación. Análisis del caso particular de las fuerzas de corte en el
proceso de fresado de alta velocidad.
2. El diseño e implementación de un modelo dinámico computacionalmente eficiente de
un proceso de fabricación complejo. Aplicación al proceso de taladrado de alto
rendimiento.
3. El diseño e implementación de un modelo dinámico múltiple entrada-simple salida,
computacionalmente eficiente, de un proceso de fabricación complejo. Aplicación al
caso particular del proceso de fresado de alta velocidad.
4. El diseño e implementación de una estrategia de control lineal para la optimización
en tiempo real de un proceso de fabricación. Aplicación al caso particular del
proceso de taladrado de alto rendimiento.
5. El diseño e implementación de una estrategia de control por modelo interno, basada
en modelos híbridos neuro-borrosos, para la optimización en tiempo real de un
proceso de fabricación complejo. Aplicación al proceso de taladrado de alto
rendimiento.
6. El diseño e implementación de estrategias de control basadas en la hibridación entre
los sistemas de inferencia borrosa y métodos metaheurísticos, con vistas a optimizar
un proceso electromecánico complejo. En específico métodos de inspiración biológica
como los algoritmos genéticos, métodos basados en la simulación de procesos físicos
como el algoritmo de temple simulado y el algoritmo de entropía cruzada y métodos
de búsqueda directa como el algoritmo de Nelder-Mead. Aplicación de estas
206 Conclusiones
Raúl M. del Toro
estrategias al caso particular del control en tiempo real del proceso de taladrado de
alto rendimiento.
7. El desarrollo de un estudio comparativo entre técnicas de soft-computing y
estrategias lineales para el control en tiempo real de un proceso de fabricación
complejo.
8. El desarrollo de un estudio comparativo de los resultados obtenidos de la aplicación
de las técnicas de optimización al ajuste de un sistema de inferencia borrosa aplicado
al control de un proceso de fabricación complejo.
Lista de publicaciones
Algunas de las contribuciones científicas y técnicas realizadas durante toda la etapa de
investigación de esta Tesis Doctoral, han sido puestas a consideración por la comunidad
científica en diferentes publicaciones realizadas, ya sea en artículos de revistas indexadas
(SCI, Science Citation Index), actas de congresos científicos o en patentes solicitadas, las
cuales se relacionan a continuación.
Artículos en revistas SCI:
R. E. Haber, A. Gajate, S. Y. Liang, R. Haber-Haber, and R. M. del Toro, "An Optimal
Fuzzy Controller for a High-Performance Drilling Process Implemented over an
Industrial Network," International Journal of Innovative Computing, Information and
Control, vol. 7, March 2011 (In Press).
A. Gajate, R. Haber, R. del Toro, P. Vega, and A. Bustillo, "Tool wear monitoring
using neuro-fuzzy techniques: a comparative study in a turning process," Journal of
Intelligent Manufacturing, pp. 1-14, 2010.
R. E. Haber, R. M. del Toro, and A. Gajate, "Optimal fuzzy control system using the
cross-entropy method. A case study of a drilling process," Information Sciences, vol.
180, pp. 2777-2792, 2010.
A. Gajate, R. E. Haber, and R. M. del Toro, "Neurofuzzy Drilling Force-Based Control
in an Ethernet-Based Application," International Journal of Innovative Computing,
Information and Control, vol. 6, pp. 373–386, January 2010.
R. M. del Toro-Matamoros, R. E. Haber, J. Pamies, and A. Bustillo, "Modelling and
simulation of the high-speed milling process based on MATLAB-SIMULINK," Revista
de Metalurgia, vol. 44, pp. 176-188, 2008.
207
Métodos clásicos y de soft-computing en la optimización de procesos complejos.
R. Haber-Haber, R. Haber, M. Schmittdiel, and R. M. del Toro, "A classic solution for
the control of a high-performance drilling process," International Journal of Machine
Tools and Manufacture, vol. 47, pp. 2290-2297, 2007.
Patentes solicitadas:
R. E. Haber Guerra, F. Martínez Puente, R. M. del Toro, and B. Caballero Retamosa,
"Fuzzy logic-based control methods for drilling". Consejo Superior de Investigaciones
Científicas, Spain. PCT/ES2008/070039, 2008 (Pendiente).
R. E. Haber Guerra, F. Martínez Puente, R. M. del Toro, M. C. Schmittdiel, and D.
Martín Andrés. "Modelling method and system for drilling processes". Consejo Superior
de Investigaciones Científicas, Spain. PCT/ES2008/070038, 2008 (Pendiente).
Trabajo futuro
El desarrollo de esta Tesis Doctoral ha sido fruto de trabajos de investigación enfocados,
principalmente, en el campo de los sistemas de control inteligente y su aplicación a la
optimización de procesos complejos, basándose fundamentalmente en el desarrollo de
estrategias de hibridación de técnicas del soft-computing.
A pesar de los resultados alcanzados y las contribuciones realizadas, esta disertación sólo
constituye una etapa dentro de una línea de investigación más amplia, cuyo objetivo principal
apunta hacia el desarrollo de un método para el diseño de un sistema de control cognitivo
artificial, basado en modelos internos, para procesos de mecanizado de alto rendimiento.
Desde el punto de vista científico, el objetivo es evolucionar los controladores inteligentes a
208 Conclusiones
Raúl M. del Toro
209
BIBLIOGRAFÍA
211
Métodos clásicos y de soft-computing en la optimización de procesos complejos.
[11] M. Dorigo and M. Birattari, "Swarm intelligence," Scholarpedia, vol. 2(9), 2007.
[12] E. Bonabeau, M. Dorigo, and G. Theraulaz, Swarm Intelligence: From Natural to
Artificial System. New York: Oxford University Press, 1999.
[13] M. Dorigo, V. Maniezzo, and A. Colorni, "Ant system: Optimization by a colony of
cooperating agents," IEEE Transactions on Systems, Man, and Cybernetics, Part B:
Cybernetics, vol. 26, pp. 29-41, 1996.
[14] J. Kennedy and R. Eberhart, "Particle swarm optimization," in IEEE International
Conference on Neural Networks - Conference Proceedings, Perth, Aust, 1995, pp.
1942-1948.
[15] T. Weise, "Global Optimization Algorithms. Theory and Application," http://www.it-
weise.de, 2009.
[16] S. J. Russell and P. Norvig, Artificial Intelligence: A Modern Approach., 2nd ed.:
Prentice Hall, 2002.
[17] V. Černý, "Thermodynamical approach to the traveling salesman problem: An
efficient simulation algorithm," Journal of Optimization Theory and Applications, vol.
45, pp. 41-51, 1985.
[18] S. Kirkpatrick, C. D. Gelatt Jr, and M. P. Vecchi, "Optimization by Simulated
Annealing," Science, vol. 220, pp. 671-680, 1983.
[19] F. Glover, "Future paths for integer programming and links to artificial intelligence,"
Computers & Operations Research, vol. 13, pp. 533-549, 1986.
[20] R. Y. Rubinstein, "Optimization of computer simulation models with rare events,"
European Journal of Operational Research, vol. 99, pp. 89-112, 1997.
[21] J. Nelder and R. Mead, "A simplex method for function minimization," Computer
Journal, vol. 7, pp. 308-313, 1965.
[22] P. P. Bonissone, "Hybrid Soft Computing Systems: Where Are We Going?," in
Proceedings of the 14th European Conference on Artificial Intelligence (ECAI 2000),
Berlin, Germany, 2000, pp. 739-746.
[23] E. H. Mamdani and S. Assilian, "EXPERIMENT IN LINGUISTIC SYNTHESIS
WITH A FUZZY LOGIC CONTROLLER," International Journal of Man-Machine
Studies, vol. 7, pp. 1-13, 1975.
[24] T. Takagi and M. Sugeno, "Fuzzy identification of systems and its applications to
modeling and control," IEEE Transactions on Systems, Man and Cybernetics, vol. 15,
pp. 116-132, 1985.
[25] M. Sugeno and G. T. Kang, "Structure identification of fuzzy model," Fuzzy Sets and
Systems, vol. 28, pp. 15-33, 1988.
[26] J. Moody and C. J. Darken, "Fast Learning in Networks of Locally-Tuned Processing
Units," Neural Computation, vol. 1, pp. 281-294, 1989.
[27] W. T. Miller Iii, F. H. Glanz, and L. G. Kraft Iii, "CMAC: An associative neural
network alternative to backpropagation," Proceedings of the IEEE, vol. 78, pp. 1561-
1567, 1990.
212 Bibliografía
Raúl M. del Toro
213
Métodos clásicos y de soft-computing en la optimización de procesos complejos.
[45] V. Maniezzo, "Genetic evolution of the topology and weight distribution of neural
networks," IEEE Transactions on Neural Networks, vol. 5, pp. 39-53, 1994.
[46] R. A. Jacobs, "Increased rates of convergence through learning rate adaptation,"
Neural Networks, vol. 1, pp. 295-307, 1988.
[47] P. Arabshahi, J. J. Choi, R. J. Marks, and T. P. Caudell, "Fuzzy control of
backpropagation," in IEEE International Conference on Fuzzy Systems (FUZZ-
IEEE'92), San Diego, CA, USA, 1992, pp. 967-972.
[48] F. Herrera and M. Lozano, "Adaptive genetic algorithms based on fuzzy techniques,"
in Information Processing and Management of Uncertainly in Knowledge-Based
Systems (IPMU'96), Granada, Spain, 1996, pp. 775-780.
[49] R. Kruse, "Fuzzy neural network," Scholarpedia, vol. 3(11), 2008.
[50] D. Nauck, "A Fuzzy Perceptron as a Generic Model for Neuro-Fuzzy Approaches," in
Proceedings of Fuzzy-systeme'94, Second GI-workshop, Munich, 1994.
[51] D. Nauck and R. Kruse, "Neuro-Fuzzy Classification with NEFCLASS," in
Operations Research Proceedings 1995, Berlin, 1996, pp. 294-299.
[52] J.-S. R. Jang, "ANFIS: adaptive-network-based fuzzy inference system," IEEE
Transactions on Systems, Man and Cybernetics, vol. 23, pp. 665-685, 1993.
[53] J. S. R. Jang and C. T. Sun, "Functional equivalence between radial basis function
networks and fuzzy inference systems," Neural Networks, IEEE Transactions on, vol.
4, pp. 156-159, 1993.
[54] H. Bersini, G. Bontempi, and C. Decaestecker, "Comparing RBF and fuzzy inference
systems on theoretical and practical basis," Proc. Int. Conf. Artificial Neural
Networks. ICANN '95, vol. 1, pp. 169-174, 1995.
[55] O. Cordón, F. Herrera, and M. Lozano, "A Classified Review On The Combination
Fuzzy Logic-Genetic Algorithms Bibliography: 1989 -1995," Dept. of Computer
Science and A.I., University of Granada. Technical Report DECSAI-95129, December
1996.
[56] C. L. Karr, "Design of an adaptive fuzzy logic controller using genetic algorithms," in
Proc. Int. Conf. Genetic Algorithms (ICGA’91), San Diego, CA, 1991, pp. 450–456.
[57] M. A. Lee and H. Tagaki, "Dynamic control of genetic algorithm using fuzzy logic
techniques," in Proc. 5th Int. Conf. Genetic Algorithms., San Mateo, CA, 1992, pp.
76–83.
[58] P. P. Bonissone, P. S. Khedkar, and Y. Chen, "Genetic algorithms for automated
tuning of fuzzy controllers: a transportation application," in Proceedings of the Fifth
IEEE International Conference on Fuzzy Systems (FUZZ-IEEE’96), New Orleans,
LA, 1996, pp. 674-680.
[59] L. Zheng, "A practical guide to tune of proportional and integral (PI) like fuzzy
controllers," in IEEE International Conference on Fuzzy Systems, 1992, pp. 633-640.
[60] M. Correa, C. Bielza, M. D. J. Ramirez, and J. R. Alique, "A Bayesian network model
for surface roughness prediction in the machining process," International Journal of
Systems Science, vol. 39, pp. 1181-1192, 2008.
214 Bibliografía
Raúl M. del Toro
[61] E. B. Flynn and M. D. Todd, "Bayesian probabilistic structural modeling for optimal
sensor placement in ultrasonic guided wave-based structural health monitoring," in
Proceedings of SPIE - The International Society for Optical Engineering, San Diego,
CA, 2010.
[62] P. Elliott and B. Williams, "DNNF-based belief state estimation," in Proceedings of
the National Conference on Artificial Intelligence, Boston, MA, 2006, pp. 36-41.
[63] C. Alippi and V. Piuri, "Experimental neural networks for prediction and
identification," IEEE Transactions on Instrumentation and Measurement, vol. 45, pp.
670-676, 1996.
[64] D. W. Coit, B. T. Jackson, and A. E. Smith, "Static neural network process models:
Considerations and case studies," International Journal of Production Research, vol.
36, pp. 2953-2967, 1998.
[65] C. Feng, X. Wang, and Z. Yu, "Neural networks modelling of honing surface
roughness parameter defined by ISO 13565," SIAM Journal of Manufacturing
Systems, vol. 21, pp. 1-35, 2002.
[66] C. Barnes, S. Brown, G. Flake, R. Jones, M. O. Rourke, and Y. C. Lee, "Applications
of neural networks to process control and modelling," Artificial Neural Networks,
Proceedings of 1991 Internal. Conf. Artif. Neur. Nets, vol. 1, pp. 321-326, 1991.
[67] K. J. Hunt, D. Sbarbaro, R. Zbikowski, and P. J. Gawthrop, "Neural networks for
control systems - A survey," Automatica, vol. 28, pp. 1083-1112, 1992.
[68] L. Jin, P. N. Nikiforuk, and M. M. Gupta, "Dynamic recurrent neural networks for
control of unknown nonlinear systems," Journal of Dynamic Systems, Measurement
and Control, Transactions of the ASME, vol. 116, pp. 567-576, 1994.
[69] K. S. Narendra and K. Parthasarathy, "Identification and control of dynamical systems
using neural networks," IEEE Transactions on Neural Networks, vol. 1, pp. 4-27,
1990.
[70] R. Babuska and H. B. Verbruggen, "An overview of fuzzy modeling for control,"
Control Engineering Practice, vol. 4, pp. 1593-1606, 1996.
[71] R. E. Haber, A. Alique, J. R. Alique, R. Haber-Haber, and S. Ros, "Current trends and
future developments of new control systems based on fuzzy logic and their application
to high speed machining," Revista de Metalurgia (Madrid), vol. 38, pp. 124-133,
2002.
[72] M. Santos, R. López, and J. M. De La Cruz, "A neuro-fuzzy approach to fast ferry
vertical motion modelling," Engineering Applications of Artificial Intelligence, vol.
19, pp. 313-321, 2006.
[73] M. A. Denaï, F. Palis, and A. Zeghbib, "Modeling and control of non-linear systems
using soft computing techniques," Applied Soft Computing Journal, vol. 7, pp. 728-
738, 2007.
[74] F. C. Sun, Z. Q. Sun, and H. X. Li, "Stable adaptive controller design of robotic
manipulators via neuro-fuzzy dynamic inversion," Journal of Robotic Systems, vol. 22,
pp. 809-819, 2005.
215
Métodos clásicos y de soft-computing en la optimización de procesos complejos.
216 Bibliografía
Raúl M. del Toro
[90] Z. J. Viharos, L. Monostori, and S. Markos, "Selection of input and output variables of
ANN based modeling of cutting processes," Proceedings of the X. Workshop on
Supervising and Diagnostics of Machining Systems of CIRP, vol. 26, p. 121, 1999.
[91] G. Sathyanarayanan, I. J. Lin, and M.-K. Chen, "Neural network modelling and
multiobjective optimization of creep feed grinding of superalloys," International
Journal of Production Research, vol. 30, pp. 2421-2438, 1992.
[92] R. K. Jain, V. K. Jain, and P. K. Kalra, "Modelling of abrasive flow machining
process: A neural network approach," Wear, vol. 231, pp. 242-248, 1999.
[93] W. Grzesik and S. Brol, "Hybrid approach to surface roughness evaluation in
multistage machining processes," Journal of Materials Processing Technology, vol.
134, pp. 265-272, 2003.
[94] Y. C. Shin and P. Vishnupad, "Neuro-fuzzy control of complex manufacturing
processes," International Journal of Production Research, vol. 34, pp. 3291-3309,
1996.
[95] Y. T. Chen and S. R. T. Kumara, "Fuzzy logic and neural networks for design of
process parameters: A grinding process application," International Journal of
Production Research, vol. 36, pp. 395-415, 1998.
[96] R. J. Kuo and P. H. Cohen, "Manufacturing process control through integration of
neural networks and fuzzy model," Fuzzy Sets and Systems, vol. 98, pp. 15-31, 1998.
[97] M. K. Kamatala, E. T. Baumgartner, and K. S. Moon, "Turned surface finish
prediction based on fuzzy logic theory," Proceedings of the 20th International
Conference on Computer and Industrial Engineering, vol. 1, pp. 101-104, 1996.
[98] K. Hashmi, M. A. El Baradie, and M. Ryan, "Fuzzy logic based intelligent selection of
machining parameters," Computers and Industrial Engineering, vol. 35, pp. 571-574,
1998.
[99] M. Milfelner, J. Kopac, F. Cus, and U. Zuperl, "Genetic equation for the cutting force
in ball-end milling," Journal of Materials Processing Technology, vol. 164-165, pp.
1554-1560, 2005.
[100] T. M. A. Maksoud and M. R. Atia, "Review of intelligent grinding and dressing
operations," Machining Science and Technology, vol. 8, pp. 263-276, 2004.
[101] C. S. Rao and R. R. Srikant, "Tool wear monitoring - An intelligent approach,"
Proceedings of the Institution of Mechanical Engineers, Part B: Journal of
Engineering Manufacture, vol. 218, pp. 905-912, 2004.
[102] P. Fu and A. D. Hope, "An intelligent signal feature pattern recognition architecture
for condition monitoring of automatic machining processes," in Proceedings - 2004
International Conference on Intelligent Mechatronics and Automation, Chengdu,
2004, pp. 552-556.
[103] Y. M. Ali and L. C. Zhang, "A fuzzy model for predicting burns in surface grinding of
steel," International Journal of Machine Tools and Manufacture, vol. 44, pp. 563-571,
2004.
217
Métodos clásicos y de soft-computing en la optimización de procesos complejos.
[104] M. Liang, T. Yeap, and A. Hermansyah, "A fuzzy system for chatter suppression in
end milling," Proceedings of the Institution of Mechanical Engineers, Part B: Journal
of Engineering Manufacture, vol. 218, pp. 403-417, 2004.
[105] R. E. Haber and A. Alique, "Intelligent process supervision for predicting tool wear in
machining processes," Mechatronics, vol. 13, pp. 825-849, 2003.
[106] W. L. R. Ip, "A fuzzy basis material removal optimization strategy for sculptured
surface machining using ball-nosed cutters," International Journal of Production
Research, vol. 36, pp. 2553-2571, 1998.
[107] Y. H. Lee, B. H. Yang, and K. S. Moon, "An economic machining process model
using fuzzy non-linear programming and neural network," International Journal of
Production Research, vol. 37, pp. 835-847, 1999.
[108] H. Al-Wedyan, K. Demirli, and R. Bhat, "A technique for fuzzy logic modeling of
machining process," in Annual Conference of the North American Fuzzy Information
Processing Society - NAFIPS, Vancouver, BC, 2001, pp. 3021-3026.
[109] G. G. Wang and S. Q. Xie, "Optimal process planning for a combined punch-and-laser
cutting machine using ant colony optimization," International Journal of Production
Research, vol. 43, pp. 2195-2216, 2005.
[110] G. Skrabalak, M. Zybura-Skrabalak, and A. Ruszaj, "Building of rules base for fuzzy-
logic control of the ECDM process," Journal of Materials Processing Technology,
vol. 149, pp. 530-535, 2004.
[111] A. Behrens and J. Ginzel, "Neuro-fuzzy process control system for sinking EDM,"
Journal of Manufacturing Processes, vol. 5, pp. 33-39, 2003.
[112] X. Shi, P. Schillings, and D. Boyd, "Applying artificial neural networks and virtual
experimental design to quality improvement of two industrial processes,"
International Journal of Production Research, vol. 42, pp. 101-118, 2004.
[113] B. F. Yousef, G. K. Knopf, E. V. Bordatchev, and S. K. Nikumb, "Neural network
modeling and analysis of the material removal process during laser machining,"
International Journal of Advanced Manufacturing Technology, vol. 22, pp. 41-53,
2003.
[114] Z. C. Lin and C. Y. Liu, "Analysis and application of the adaptive neuro-fuzzy
inference system in prediction of CMP machining parameters," International Journal
of Computer Applications in Technology, vol. 17, pp. 80-89, 2003.
[115] N. Ding, L. Wang, G. Li, D. Chen, X. Li, and D. Zhang, "Human-machine intelligent
control system for grinding," Chinese Journal of Mechanical Engineering (English
Edition), vol. 17, pp. 26-29, 2004.
[116] J. Balic, "Neural-network-based numerical control for milling machine," Journal of
Intelligent and Robotic Systems: Theory and Applications, vol. 40, pp. 343-358, 2004.
[117] R. J. Lian, B. F. Lin, and J. H. Huang, "A grey prediction fuzzy controller for constant
cutting force in turning," International Journal of Machine Tools and Manufacture,
vol. 45, pp. 1047-1056, 2005.
218 Bibliografía
Raúl M. del Toro
219
Métodos clásicos y de soft-computing en la optimización de procesos complejos.
[133] M. Kupczyk and W. Misiak, "Tool life of the edges coated with the c-BN plus h-BN
coatings with different structures during hard machinable steel machining," Revista de
Metalurgia, vol. SPEC. VOL., pp. 483-487, 2005.
[134] J. M. Sánchez-Sola, M. A. Sebastián, M. S. Carrilero, J. M. González, and M. Marcos,
"Characterization of the built-up edge and the built-up layer in the machining process
of AA7050 alloy," Revista de Metalurgia, vol. SPEC. VOL., pp. 365-368, 2005.
[135] R. E. Haber, J. E. Jiménez, A. Jiménez, and J. López-Coronado, "Cutting force model
for high speed machining process," Revista de Metalurgia, vol. 40, pp. 247-258, Jul-
Aug 2004.
[136] J. L. Cantero, J. M. Sanchez-Sola, M. M. Tardio, M. A. Sebastian, M. S. Carrilero, M.
Marcos, and M. H. Miguelez, "Tool wear in the dry drilling of Ti-6Al-4V. A
SEM/EDS study," Revista de Metalurgia, vol. SPEC. VOL., pp. 355-360, 2005.
[137] M. Álvarez, M. S. Carrilero, J. E. Ares, J. M. González, and M. Marcos, "Influence of
the cutting speed-feed relation on the surface quality of drilled carbon fiber samples,"
Revista de Metalurgia, vol. SPEC. VOL., pp. 361-365, 2005.
[138] F. Cus, M. Milfelner, and J. Balic, "An intelligent system for monitoring and
optimization of ball-end milling process," Journal of Materials Processing
Technology, vol. 175, pp. 90-97, Jun 2006.
[139] B. M. Imani, M. H. Sadeghi, and M. A. Elbestawi, "An improved process simulation
system for ball-end milling of sculptured surfaces," International Journal of Machine
Tools & Manufacture, vol. 38, pp. 1089-1107, Sep 1998.
[140] S. Engin and Y. Altintas, "Mechanics and dynamics of general milling cutters. Part I:
helical end mills," International Journal of Machine Tools & Manufacture, vol. 41,
pp. 2195-2212, Dec 2001.
[141] M. H. Sadeghi, H. Haghighat, and M. A. Elbestawi, "A solid modeler based ball-end
milling process simulation," International Journal of Advanced Manufacturing
Technology, vol. 22, pp. 775-785, Dec 2003.
[142] E. E. M. Lim, H. Y. Feng, C. H. Menq, and Z. H. Lin, "The Prediction of Dimensional
Error for Sculptured Surface Productions Using the Ball-End Milling Process .1. Chip
Geometry Analysis and Cutting Force Prediction," International Journal of Machine
Tools & Manufacture, vol. 35, pp. 1149-1169, Aug 1995.
[143] Y. Altintas and P. Lee, "Mechanics and dynamics of ball end milling," Journal of
Manufacturing Science and Engineering, Transactions of the ASME, vol. 120, pp.
684-691, 1998.
[144] J. Tlusty and F. Ismail, "Special Aspects of Chatter in Milling," Journal of Vibration
Acoustics Stress and Reliability in Design-Transactions of the ASME, vol. 105, pp. 24-
32, 1983.
[145] Y. Altintas, Manufacturing Automation. Cambridge, United Kingdom: Cambridge
University Press, 2000.
[146] M. E. Martellotti, "An analysis of the milling process.," Transactions of the ASME,
vol. 63, pp. 677–700, 1941.
220 Bibliografía
Raúl M. del Toro
[147] M. E. Martellotti, "An analysis of the milling process. Part II: down milling.,"
Transactions of the ASME, vol. 67, pp. 233–251, 1945.
[148] F. Abrari, M. A. Elbestawi, and A. D. Spence, "On the dynamics of ball end milling:
Modeling of cutting forces and stability analysis," International Journal of Machine
Tools & Manufacture, vol. 38, pp. 215-237, Mar 1998.
[149] X. W. Liu, K. Cheng, A. P. Longstaff, M. H. Widiyarto, and D. Ford, "Improved
dynamic cutting force model in ball-end milling. Part I: theoretical modelling and
experimental calibration," International Journal of Advanced Manufacturing
Technology, vol. 26, pp. 457-465, Sep 2005.
[150] R. F. Hamade and F. Ismail, "A case for aggressive drilling of aluminum," Journal of
Materials Processing Technology, vol. 166, pp. 86-97, Jul 2005.
[151] L. H. Wang, P. Orban, A. Cunningham, and S. Lang, "Remote real-time CNC
machining for web-based manufacturing," Robotics and Computer-Integrated
Manufacturing, vol. 20, pp. 563-571, Dec 2004.
[152] F. R. Biglari and X. D. Fang, "Real-Time Fuzzy-Logic Control for Maximizing the
Tool Life of Small-Diameter Drills," Fuzzy Sets and Systems, vol. 72, pp. 91-101, May
1995.
[153] M. Nouari, G. List, F. Girot, and D. Coupard, "Experimental analysis and optimisation
of tool wear in dry machining of aluminium alloys," Wear, vol. 255, pp. 1359-1368,
Aug-Sep 2003.
[154] R. G. Landers, A. G. Ulsoy, and Y. H. Ma, "A comparison of model-based machining
force control approaches," International Journal of Machine Tools & Manufacture,
vol. 44, pp. 733-748, Jun 2004.
[155] J. B. Kim, S. J. Lee, and Y. P. Park, "Stable and Efficient Drilling Process by Active
Control of the Thrust Force," Mechanical Systems and Signal Processing, vol. 8, pp.
585-595, Sep 1994.
[156] M. Elhachimi, S. Torbaty, and P. Joyot, "Mechanical modelling of high speed drilling.
1. Predicting torque and thrust," International Journal of Machine Tools &
Manufacture, vol. 39, pp. 553-568, Apr 1999.
[157] M. Elhachimi, S. Torbaty, and P. Joyot, "Mechanical modelling of high speed drilling.
2: Predicted and experimental results," International Journal of Machine Tools &
Manufacture, vol. 39, pp. 569-581, Apr 1999.
[158] Y. Sheng and M. Tomizuka, "Intelligent modeling of thrust force in drilling process,"
Journal of Dynamic Systems Measurement and Control-Transactions of the Asme, vol.
128, pp. 846-855, Dec 2006.
[159] R. E. Haber and J. R. Alique, "Nonlinear internal model control using neural
networks: an application for machining processes," Neural Computing & Applications,
vol. 13, pp. 47-55, Apr 2004.
[160] Y. M. Zhang, X. L. Lin, X. D. Wang, and G. Q. Cai, "The study on the dynamic
characteristic of high speed machine tool and experiment validate," in Advances in
Materials Manufacturing Science and Technology. vol. 471-472 Zurich-Uetikon:
Trans Tech Publications Ltd, 2004, pp. 765-769.
221
Métodos clásicos y de soft-computing en la optimización de procesos complejos.
[161] Y. T. Oh, W. T. Kwon, and C. N. Chu, "Drilling torque control using spindle motor
current and its effect on tool wear," International Journal of Advanced Manufacturing
Technology, vol. 24, pp. 327-334, 2004.
[162] R. J. Furness, T. C. Tsao, J. S. Rankin, M. J. Muth, and K. W. Manes, "Torque control
for a form tool drilling operation," Ieee Transactions on Control Systems Technology,
vol. 7, pp. 22-30, Jan 1999.
[163] Y. T. Oh, G. D. Kim, and C. N. Chu, "Design of a drilling torque controller for a
machining center," International Journal of Advanced Manufacturing Technology,
vol. 22, pp. 329-335, Oct 2003.
[164] M. Pirtini and I. Lazoglu, "Forces and hole quality in drilling," International Journal
of Machine Tools & Manufacture, vol. 45, pp. 1271-1281, Sep 2005.
[165] J. Nilsson, B. Bernhardsson, and B. Wittenmark, "Stochastic analysis and control of
real-time systems with random time delays," Automatica, vol. 34, pp. 1154-1154, Sep
1998.
[166] L. A. Montestruque and P. Antsaklis, "Stability of model-based networked control
systems with time-varying transmission times," Ieee Transactions on Automatic
Control, vol. 49, pp. 1562-1572, Sep 2004.
[167] J. R. Dormand and P. J. Prince, "A family of embedded Runge-Kutta formulae,"
Journal of Computational and Applied Mathematics, vol. 6, pp. 19-26, 1980.
[168] R. M. Del Toro-Matamoros, R. E. Haber, J. Pamies, and A. Bustillo, "Modelling and
simulation of the high-speed milling process based on MATLAB-SIMULINK,"
Revista de Metalurgia, vol. 44, pp. 176-188, 2008.
[169] R. M. Del Toro, M. C. Schmittdiel, R. E. Haber-Guerra, and R. Haber-Haber, "System
identification of the high performance drilling process for network-based control," in
2007 Proceedings of the ASME International Design Engineering Technical
Conferences and Computers and Information in Engineering Conference, DETC2007,
Las Vegas, NV, 2008, pp. 827-834.
[170] K. J. Äström and T. Hagglund, PID controllers : theory, design, and tuning: Research
Triangle Park, NC : Instrument Society of America, 1995.
[171] A. Visioli, "Optimal tuning of PID controllers for integral and unstable processes,"
IEE Proceedings: Control Theory and Applications, vol. 148, pp. 180-184, 2001.
[172] B. Russell, "Vagueness," Australasian Journal of Philosophy and Psychology, vol. 1,
pp. 84-92, 1923.
[173] L. A. Zadeh, "OUTLINE OF A NEW APPROACH TO THE ANALYSIS OF
COMPLEX SYSTEMS AND DECISION PROCESSES," IEEE Transactions on
Systems, Man and Cybernetics, vol. SMC-3, pp. 28-44, 1973.
[174] L. A. Zadeh, "Fuzzy sets," Information and Control, vol. 8, pp. 338-353, 1965.
[175] R. M. Tong, "A control engineering review of fuzzy systems," Automatica, vol. 13,
pp. 559-569, 1977.
[176] M. Braae and D. A. Rutherford, "Theoretical and linguistic aspects of the fuzzy logic
controller," Automatica, vol. 15, pp. 553-577, 1979.
222 Bibliografía
Raúl M. del Toro
[177] C. C. Lee, "Fuzzy logic in control systems: Fuzzy logic controller - Part I," IEEE
Transactions on Systems, Man and Cybernetics, vol. 20, pp. 404-418, 1990.
[178] C. Elkan, "Paradoxical success of fuzzy logic," in Proceedings of the National
Conference on Artificial Intelligence, Washington, DC, USA, 1993, pp. 698-703.
[179] P. Bernhard, "Fuzzy Control: Facts, Japan and Europe," TECNews, pp. 4-5, June 1992.
[180] A. Kaufmann, "Introduction to the Theory of Fuzzy Subsets.," New York: Academic
Press, 1975.
[181] H. Zimmermann, "Fuzzy Sets Theory and Its Applications," Kluwer Academic Pub.,
1985.
[182] D. Dubois and H. Prade, "Fuzzy sets, probability and measurement," European
Journal of Operational Research, vol. 40, pp. 135-154, 1989.
[183] C. C. Lee, "Fuzzy logic in control systems: Fuzzy logic controller - Part II," IEEE
Transactions on Systems, Man and Cybernetics, vol. 20, pp. 419-435, 1990.
[184] L. A. Zadeh, "Similarity relations and fuzzy orderings," Information Sciences, vol. 3,
pp. 177-200, 1971.
[185] L. A. Zadeh, "The concept of a linguistic variable and its application to approximate
reasoning-I," Information Sciences, vol. 8, pp. 199-249, 1975.
[186] R. E. Haber Guerra, "Control Jerárquico Inteligente. Aplicación al Proceso de
Fresado," in Dpto. de Automática, Ingeniería Electrónica e Informática Industrial.
vol. Ph D. Madrid: Universidad Politécnica de Madrid, E. T. S. de Ingenieros
Industriales, 1998.
[187] F. Matía and A. Jiménez, "On optimal implementation of fuzzy controllers," Int J
Intell Contr Syst, vol. 1, pp. 407-415, 1996.
[188] R. M. Tong, "The construction and evaluation of fuzzy models," Advances in Fuzzy
Set Theory and Applications, pp. 559-576, 1979.
[189] C. A. Cabezas and R. Haber, "Objetivación de los algoritmos de control borrosos," in
IV Congreso Latinoamericano de Control Automático Puebla, México, 1990.
[190] C. Rodríguez, R. Haber, and J. R. Alique, "Modeling of Machine-Tool Cutting
Process Applying Fuzzy Logic and Clustering Analysis," in II European Congress on
Intelligent Techniques and Soft Computing, EUFIT'94. vol. II Aachen, Germany,
1994, pp. 1071-1075.
[191] T. Takagi and M. Sugeno, "DERIVATION OF FUZZY CONTROL RULES FROM
HUMAN OPERATOR'S CONTROL ACTIONS," in IFAC Proceedings Series,
Marseille, Fr, 1984, pp. 55-60.
[192] W. Pedrycz, Fuzzy Control and Fuzzy Systems. New York: John Wiley & Sons, Inc.,
1989.
[193] P. J. King and E. H. Mamdani, "The application of fuzzy control systems to industrial
processes," Automatica, vol. 13, pp. 235-242, 1977.
[194] M. Braae and D. A. Rutherford, "Selection of parameters for a fuzzy logic controller,"
Fuzzy Sets and Systems, vol. 2, pp. 185-199, 1979.
223
Métodos clásicos y de soft-computing en la optimización de procesos complejos.
224 Bibliografía
Raúl M. del Toro
[212] S. Altug, M. Y. Chow, and H. Joel Trussell, "Fuzzy inference systems implemented on
neural architectures for motor fault detection and diagnosis," IEEE Transactions on
Industrial Electronics, vol. 46, pp. 1069-1079, 1999.
[213] L. C. Chang and F. J. Chang, "Intelligent control for modelling or real-time reservoir
operation," Hydrological Processes, vol. 15, pp. 1621-1634, 2001.
[214] S. Singer and J. Nelder, "Nelder-Mead algorithm," Scholarpedia, vol. 4(7), 2009.
[215] S. Kirkpatrick, C. D. Gelatt, and M. P. Vecchi, "Optimization by Simulated
Annealing," Science, vol. 220, pp. 671-680, 1983.
[216] J. H. Holland, Adaptation in natural and artificial system: an introduction with
application to biology, control and artificial intelligence: Ann Arbor, University of
Michigan Press, 1975.
[217] J. J. Grefenstette, " Genetic Algorithms and Their Applications," in Second
International Conference on Genetic Algorithms, Cambridge, MA, 1987.
[218] P. T. De Boer, D. P. Kroese, S. Mannor, and R. Y. Rubinstein, "A tutorial on the
cross-entropy method," Annals of Operations Research, vol. 134, pp. 19-67, 2005.
[219] R. Rubinstein and D. Kroese, The Cross-Entropy Method: A Unified Approach to
Combinatorial Optimization, Monte-Carlo Simulation and Machine Learning. Berlin,
Heidelberg, New York: Springer, 2004.
[220] T. G. Kolda, R. M. Lewis, and V. Torczon, "Optimization by direct search: New
perspectives on some classical and modern methods," SIAM Review, vol. 45, pp. 385-
482, 2003.
[221] J. C. Lagarias, J. A. Reeds, M. H. Wright, and P. E. Wright, "Convergence properties
of the Nelder-Mead simplex method in low dimensions," SIAM Journal on
Optimization, vol. 9, pp. 112-147, 1999.
[222] N. Metropolis, A. W. Rosenbluth, M. N. Rosenbluth, A. H. Teller, and E. Teller,
"Equation of state calculations by fast computing machines," The Journal of Chemical
Physics, vol. 21, pp. 1087-1092, 1953.
[223] W. K. Hastings, "Monte carlo sampling methods using Markov chains and their
applications," Biometrika, vol. 57, pp. 97-109, 1970.
[224] L. Ingber, "Adaptive simulated annealing (ASA): Lessons learned," in Polish Journal
Control and Cybernetics, 1995.
[225] J. d. D. Sanz Bobi, "Control de procesos basado en eventos mediante bases de
conocimiento borrosas y algoritmos genéticos," in Automática, Ingeniería Electrónica
e Informática Industrial Madrid: Universidad Politécnica de Madrid, 2001.
[226] R. C. Fernández, J. M. Luque, M. L. Gallego, A. T. González, and T. G. Nuñez,
"Algoritmos genéticos para la resolución de problemas de Programación por Metas
Entera. Aplicación a la Economía de la Educación.," in X JORNADAS MADRID 2002.
vol. Actas_10: Rect@: Revista Electrónica de Comunicaciones y Trabajos de
ASEPUMA, 2002, p. 42.
[227] D. E. Goldberg, "A comparison of selection schemes used in genetic algorithm," in
Foundations of Genetic Algorithms, G. E. Rawlins, Ed.: Morgan Kaufmann
Publishers, 1991, pp. 69-93.
225
Métodos clásicos y de soft-computing en la optimización de procesos complejos.
226 Bibliografía
Raúl M. del Toro
227