Source: http://docplayer.es/2636738-Extensibilidad-de-criterios-de-decision-para-un-empleando-caracteristicas-en-diferentes-espacios-de-representacion.html
Timestamp: 2017-02-26 17:10:45+00:00

Document:
EXTENSIBILIDAD DE CRITERIOS DE DECISIÓN PARA UN EMPLEANDO CARACTERÍSTICAS EN DIFERENTES ESPACIOS DE REPRESENTACIÓN - PDF
EXTENSIBILIDAD DE CRITERIOS DE DECISIÓN PARA UN EMPLEANDO CARACTERÍSTICAS EN DIFERENTES ESPACIOS DE REPRESENTACIÓN
Download "EXTENSIBILIDAD DE CRITERIOS DE DECISIÓN PARA UN EMPLEANDO CARACTERÍSTICAS EN DIFERENTES ESPACIOS DE REPRESENTACIÓN"
Adolfo Araya Cordero
1 EXTENSIBILIDAD DE CRITERIOS DE DECISIÓN PARA UN CLASIFICADOR DE FALLOS EN RODAMIENTOS BASADO EN HMM, SOBRE RESULTADOS EN UN CLASIFICADOR DE MENOS ESTADOS, EMPLEANDO CARACTERÍSTICAS EN DIFERENTES ESPACIOS DE REPRESENTACIÓN Juan Camilo Castaño Guzmán Carlos Hernán Agudelo Gutiérrez Proyecto de grado presentado como requisito parcial para aspirar al título de Ingeniero Eléctricista Director Ing Mauricio Holguín Londoño, M.Sc UNIVERSIDAD TECNOLÓGICA DE PEREIRA PROGRAMA DE INGENIERÍA ELÉCTRICA PEREIRA 20152 3 Nota de Aceptación Firma del Presidente del jurado Firma del jurado 1 - Evaluador Firma del jurado 2 - Evaluador Firma del jurado 3 - Director4 DEDICATORIA Dedicado a nuestros padres, quienes siempre conservaron la rme esperanza de ver a sus hijos graduados pese a todas las dicultades en modo, tiempo y lugar. A todos mil gracias por la profunda transformación que han hecho en nuestras vidas.5 AGRADECIMIENTOS Primero que todo gracias a DIOS. A nuestros padres por haber estado siempre en pie de lucha con nuestro proceso formativo, a nuestro apreciado mentor Mauricio Holguin Londoño por propiciar las viabilidades técnicas para el presente proyecto y ser apoyo constante en todas las etapas de ejecución.6 1 RESUMEN Una de las partes fundamentales más propensas a fallar en una máquina eléctrica son los rodamientos, debido a esto la detección y el diagnóstico temprano de fallos en estos elementos son necesarios para el adecuado funcionamiento de las máquinas eléctricas. Una de las mejores herramientas con las que se cuenta hoy en día para el mantenimiento predictivo centrado en los rodamientos de una máquina es la técnica de análisis de vibraciones. El presente proyecto de grado tuvo como propósito el diagnóstico de fallas tempranas en los rodamientos, esto se logró mediante el uso de la técnica de análisis de vibraciones siendo dichas señales de vibración representadas en los dominios de tiempo, frecuencia y tiempo-frecuencia. Las señales de vibración utilizadas fueron proporcionadas por el Bearing Data Center de la Western Reserve University. A cada una de estas señales le fueron extraidas un determinado número de características representadas en diferentes dominios. Para el proceso de clasicación de fallas a las cuales pertenecía las señales extraídas se hizo uso de los modelos ocultos de Markov (HMM). Con el n de vericar la extensibilidad de criterios de decisión para un clasicador de fallos basados en HMM sobre un clasicador de menos estados se realizaron siete experimentos, en algunos de ellos se combinaron los espacios de representación (Tiempo, Frecuencia, Melcepst, Wavelet y STFT), con todas sus características evaluando de qué manera inuía la combinación de espacios en el comportamiento nal del clasicador a medida que se variaban los estados ocultos de dicho clasicador entre 1 y 8. Se evaluaron además cuales eran las características que aportaban mayor relevancia al entrenamiento. Con el n de reducir la dimensionalidad de los datos con los que se trabajó, se hizo un análisis de componentes principales a las características extraídas de los diferentes espacios de representación y se realizaron las mismas observaciones que se hicieron para los datos originales. Las primeras pesquisas permiten inferir, que hay ciertas características que presentan un nivel de representatividad mayor en comparación con las demás en cada uno de los dominios de representación.7 ABSTRACT One of the fundamental parts that are prone to failure in an electric machine are the bearings. With this in mind, the diagnostic tests and early detection of aws in these pieces of equipment are necessary for their adequate functioning. One of the best tools for predictive maintenance based around the bearings of a machine, is the technique of vibration analysis. The following dissertation was based on the early detection of problems with the bearings. This was achieved using the stated technique, represented in the time domain, frequency and time-frequency. The vibration signals used were given by the Bearing Data Center at Western Reserve University. Each signal was extracted using a certain number of characteristics represented in dierent domains. The Hidden Markov Model (HMM) was used to classify the faults of the extracted signals. Seven experiments were done to verify the extensibility of criteria of decision for a classier of faults based on HMM over a classier of less states, in some of them the representation spaces were combined (Time, Frequency, Melcepst, Wavelet, and STFT), with all of their characteristics evaluating how the combination of spaces aected the behavior of the classier at the rate that the hidden states varied from 1 to 8. The characteristics that provided the most relevance to the training were also evaluated. With the purpose of reducing the dimensionality of the data which was used, the Principal Components Analysis (PCA) technique was utilized on the extracted characteristics of the varied representation spaces as it was on the original data. The rst inquiries let us infer, that there are certain characteristics that present a greater level of representation compared to the rest in every representation domain.8 CONTENIDO pág. 1 RESUMEN I PRELIMINARES 2 Introducción 17 3 Denición del Problema 19 4 Justicación 22 5 Objetivos Objetivo General Objetivos Especícos Estado del Arte 26 II MATERIALES Y MÉTODOS 28 7 Base de Datos 29 8 Mantenimiento 30 9 Vibración Análisis de Vibraciones Técnicas de Análisis de Vibración Medición de las Vibraciones9 10 Contexto Mecánico Fallas Detectadas por el Análisis de Vibraciones Desalineación Desbalanceo Caja de engranes Soltura mecánica Rodamientos Espacios de Representación Dominio del Tiempo Dominio de la Frecuencia Dominio Tiempo-Frecuencia Segmentación y Ventanas Segmentación Ventaneo Análisis de Componentes Principales PCA Cálculo de las Componentes Principales Interpretación de las Componentes Selección del número de componentes Principales Clasicación Libro de Código (Codebook) 66 III METODOLOGÍA 7110 16 Desarrollo Experimental Base De Datos Segmentación y Ventaneo Selección de Características Tiempo Frecuencia Tiempo frecuencia Clasicación Evaluación Sintonización y Validación Prueba 1: Evaluación de cada espacio de representación con todas sus características Prueba 2: Combinación de espacios de representación con todas sus características Prueba 3: Evaluación de las mejores características por espacio de representación Prueba 4: Evaluación del número de componentes principales para cada espacio de representación Prueba 5: Representatividad de cada característica Prueba 6: Evaluación de cada espacio de representación con todas sus características utilizando PCA Prueba 7: Combinación de espacios de representación con todas sus características utilizando PCA Pseudocódigos Función Frecuencia Función STFT Función Tiempo11 18.4 Función Entrenamiento y Validación Función K-Means Función Libro de Código Función Entrenamiento y Prueba Diagnóstica con PCA Función Entrenamiento y Prueba Diagnóstica Función Wavelet Función Validación Función Extracción de Características IV RESULTADOS, ANÁLISIS Y CONCLUSIONES Prueba 1: Evaluación de cada espacio de representación con todas sus características Prueba 2: Combinación de espacios de representación con todas sus características Prueba 3: Evaluación de las mejores características por espacio de representación Frecuencia Melcepst STFT Tiempo Wavelet Características más relevantes para Frecuencia Características más relevantes para Melcepst Características más relevantes para STFT12 21.9 Características más relevantes para Tiempo Características más relevantes para Wavelet Prueba 4: Evaluación del número de componentes principales para cada espacio de representación Prueba 5: Representatividad de cada característica Prueba 6: Evaluación de cada espacio de representación con todas sus características utilizando PCA Prueba 7: Combinación de espacios de representación con todas sus características utilizando PCA Conclusiones 147 V BIBLIOGRAFÍA 149 BIBLIOGRAFÍA 15013 List of Figures 1 Tipos de mantenimiento [1] Mecanismo de yugo escocés [2] Espectro de frecuencias de funcionamiento normal [3] Espectro de frecuencias de funcionamiento defectuoso [3] Onda de desplazamiento truncada por un roce [4] Desalineamiento paralelo [5] Desalineamiento angular [5] Espectro del desbalanceo y relación de fases entre soportes del eje [5] Espectro de la soltura mecánica tipo A y relación de fases en perno o placa base [6] Espectro de la soltura mecánica tipo B [6] Espectro de la soltura mecánica tipo C [6] Partes que componen un rodamiento [6] Explicación gráca de la STFT [7] Segmentación de la señal de vibración [6] Construcción del libro de código [8] Diagrama de Voronoi [8] Clustering empleando k-means [9] Algoritmo de clustering empleando K-means [9] Distribución base de datos Resultados prueba Resultados prueba Resultados prueba 3 - Frecuencia Resultados prueba 3 - Melcepst14 24 Resultados prueba 3 -STFT Resultados prueba 3 - Tiempo Resultados prueba 3 - Wavelet Mejores características Frecuencia Mejores características melcepst Mejores características STFT Mejores características Tiempo Mejores características Wavelet Resultados prueba Resultados prueba Resultados prueba Resultados prueba15 List of Tables 1 Frecuencias de vibración y posibles causas de falla [10] Funciones ventana y sus características [11] Ventanas espectrales usadas en el procesamiento digital de datos [11] Características en el dominio del tiempo Características utilizadas en el dominio de la frecuencia Características utilizadas en el dominio tiempo-frecuencia STFT Características utilizadas en el dominio tiempo-frecuencia Wavelet Características utilizadas en el dominio tiempo-frecuencia Melcepst Resultados prueba 1 Frecuencia Resultados prueba 1 Melcepst Resultados prueba 1 STFT Resultados prueba 1 Tiempo Resultados prueba 1 Wavelet Resultados prueba Resultados prueba 3 Frecuencia Resultados prueba 3 Melcepst Resultados prueba 3 STFT Resultados prueba 3 Tiempo Resultados prueba 3 Wavelet Características más relevantes Frecuencia Características más relevantes Melcepst Características más relevantes STFT Características más relevantes Tiempo16 24 Características más relevantes Wavelet Resultados prueba 4 Frecuencia Resultados prueba 4 Melcepst Resultados prueba 4 STFT Resultados prueba 4 Tiempo Resultados prueba 4 Wavelet Resultados prueba 5 Frecuencia Resultados prueba 5 Melcepst Resultados prueba 5 STFT Resultados prueba 5 Tiempo Resultados prueba 5 Wavelet Resultados prueba 6 Frecuencia Resultados prueba 6 STFT Resultados prueba 6 Tiempo Resultados prueba 6 Wavelet Resultados prueba 6 Melcepst Resultados prueba Comparación resultados con PCA y sin PCA Melcepst Comparación resultados con PCA y sin PCA Frecuencia Comparación resultados con PCA y sin PCA STFT Comparación resultados con PCA y sin PCA Tiempo Comparación resultados con PCA y sin PCA Wavelet17 Part I PRELIMINARES18 2 Introducción Con el n de mantener una operación segura y conable en una instalación industrial es fundamental evitar la ocurrencia de posibles fallas en una máquina. Para lograr esto es necesario pronósticar, prevenir y eliminar posibles anomalías en un artefacto eléctrico, así se disminuye el riesgo de que una máquina tenga que declararse fuera de servicio debido a un fallo. Una de las partes fundamentales más propensas a fallar en una máquina eléctrica son los rodaminetos. La función de los rodamientos en una máquina eléctrica es la de proporcionar el deslizamiento del rotor con respecto al estator manteniendo un espacio de aire uniforme, por este motivo la detección y el diagnóstico temprano de fallos en estos elementos son necesarios, teniendo como prioridad evitar malos funcionamientos. Una de las mejores herramientas con las que se cuenta hoy en día para el mantenimiento predictivo centrado en los rodamientos de una máquina es la técnica de análisis de vibraciones. El tema de investigación que concierne a este documento abarca el mantenimiento predictivo aplicado al diagnóstico de fallas tempranas en los rodamientos, especícamente en lo que concierne a las anomalías que se encuentran ubicadas en el anillo interior, los elementos rodantes y en la pista de rodura exterior, haciendo uso de la técnica de análisis de vibraciones siendo dichas señales de vibración representadas en los dominios de tiempo, frecuencia y tiempo-frecuencia. Con antelación se tratan ampliamente los conceptos y términos más importantes que permiten comprender la práctica del diagnóstico de fallas en rodamientos. Empleando un análisis documental se profundiza en el concepto de rodamiento, se obtienen deniciones acerca de las fallas más comunes que estos presentan; además se amplía en las deniciones y toda la información concerniente a la implementación de un sistema capaz de diagnósticar y clasicar fallas en rodamientos. Buscando simplicar la representación de la señal en otra de igual signicancia pero que facilte su análisis, se realiza la segmentación de las señales de vibración la cuales provienen de una base de datos con pruebas de rodamientos normales y defectuosos proporcionados por el Bearing Data Center de la Western Reserve University alimentando de esta manera el proceso de extracción de características de la señal. Las señales de vibración originales pasan por varias etapas de preprocesamiento con el n de 1719 realizar una óptima extracción de características, así se calcula y almacenan atributos de la señal de vibración en formato de vectores compactos los cuales utilizan los dominios de tiempo, frecuencia y tiempo-frecuencia como espacios de representación, siendo la transformada wavelets, la transformada en tiempo corto de fourier y los coecientes cepstrales las técnicas utilizadas para representar el dominio de tiempo-frecuencia. El sistema clasicador emplea como técnica de aprendizaje de máquina los Modelos Ocultos de Markov (HMM) por medio de una distribución de observación discreta la cual le permite al clasicador identicar a que clase corresponde cada una de las fallas de los rodamientos. Cabe mencionar que se hace uso de un algoritmo de validación el cual permite vericar el correcto desempeño del sistema clasicador de fallos, dicho desempeño va a depender de diversos factores entre ellos el número de estados con los que se entrena el HMM. Para observar como es el comportamiento del clasicador a medida que el número de estados va cambiando se realizan diversas puebas. Inicialmente, se evalúa el comportamiento del clasicador utilizando cada espacio de representación (Tiempo, Frecuencia, Melcepst, Wavelet y Transformada en Tiempo Corto de Fourier), con todas sus características realizando entrenamientos para 1,2,3,4,5,6,7 y 8 estados. A continuación se combinan los espacios de representación (Tiempo, Frecuencia, Melcepst, Wavelet y Transformada en Tiempo Corto de Fourier), con todas sus características y se evalúa de qué manera inuye la combinación de espacios de representación utilizado en el comportamiento nal del clasicador a medida que se varían entre 1 y 8 los estados ocultos de dicho clasicador. Una vez que se ha observado el comportamiento del clasicador a medida que cambia el número de estados ocultos cuando se utilizan todas las características de los espacios de representación, se realizan observaciones de como se comportan las características de los diferentes espacios de representación siendo entrenadas de forma individual y como se varía su capacidad de representar la señal cuando varía el número de estados ocultos. Con el n de reducir la dimensionalidad de los datos con los que se trabaja se aplica componentes principales a las características extraídas de los diferentes espacios de representación y se realizan las mismas observaciones que se realizaron para los datos originales exceptuando el comportamiento de cada una de las componentes principales de forma individual ya que esto no tenía concordancia en los objetivos propuestos para la investigación. 1820 3 Denición del Problema El desarrollo industrial de una sociedad es de vital importancia para su avance, es por esto que en la actualidad se ha suscitado un especial interés en desarrollar una política de mantenimiento dentro de un marco de diseño de estrategias empresariales orientadas a tener cada día mejores operaciones, reducción de la cantidad y frecuencia de intervenciones, disminución efectiva de la complejidad, mejora de la organización y optimización de costos en el mantenimiento y producción [12]. A nivel industrial se emplean principalmente dos tipos de mantenimiento, uno es el mantenimiento correctivo dedicado a la reparación de los equipos después de que la avería se ha detectado, con el n de evitar consecuencias inaceptables; y el otro, es el mantenimiento preventivo el cual es ejecutado en intervalos predeterminados para así reducir la probabilidades de fallo o la degradación de funcionamiento de un elemento [1]. El interés especial de las grandes industrias por desarrollar políticas claras de mantenimiento nace de la necesidad de disminuir las pérdidas que se producen cuando una máquina se descompone e interrumpe su producción, es por eso que se destina una parte considerable del presupuesto en actividades de mantenimiento predictivo, este mantenimiento realiza un seguimiento de cada una de las variables relacionadas con el funcionamiento de las máquinas para poder predecir posibles fallas y tomar las acciones correctivas más apropiadas en el momento oportuno [13]. En la actualidad el mantenimiento predictivo implica altos costos para el sector industrial, lo que ha llevado a ser el tipo de mantenimiento menos implementado en la industria, sin poder superar aún al mantenimiento preventivo. El mantenimiento predictivo posee la capacidad de poder informar en cualquier instante de tiempo el estado general de cada una de las máquinas de la planta, permitiendo un adecuado y óptimo funcionamiento, convirtiéndose así en la forma de mantenimiento ideal para la industria. El mantenimiento predictivo consta de una serie de ensayos de carácter no destructivo como análisis de aceite, análisis de desgaste de partículas, medida de vibraciones, medición de temperaturas, termografías, etc. Está orientado a realizar un seguimiento del funcionamiento de los equipos para detectar signos de advertencia que indiquen que alguna de sus partes no está trabajando de la manera correcta [13]. 1921 Las máquinas eléctricas presentan algunos niveles de vibración aun cuando estas estén operando correctamente, pero en general las vibraciones en una máquina no son buenas, pues estas pueden causar desgaste, suras por fatiga, pérdida de efectividad de sellos, rotura de aislantes, ruido, etc. Sin embargo, las vibraciones son las mejores indicadoras de la condición mecánica de una máquina y pueden ser una herramienta de predicción de la evolución de un defecto. Las vibraciones en una máquina eléctrica están directamente relacionadas con su vida útil de dos maneras: Por un lado un bajo nivel de vibraciones es una indicación de que la máquina funciona correctamente durante un largo período de tiempo, mientras que un aumento en el nivel de vibraciones es una indicación de que la máquina se encamina hacia algún tipo de falla [10]. El procesamiento digital de la señal de vibración, para diagnóstico de fallos, puede ser abordado a partir de dos contextos: tradicional y académico e investigativo. La diferencia entre ambos depende no solo de la posición geográca sino también de los equipos e instrumentación disponible. Por ejemplo, en países desarrollados, con alto nivel de tecnología, se encuentran sistemas de monitoreo de condición complejos, que contienen un alto grado de eciencia y tienden a realizar un diagnóstico automático. Contrario sucede en países en vía de desarrollo como Colombia, donde se emplea un análisis con base en sistemas simples, que tienen una capacidad de diagnóstico muy limitada ya que solamente proveen una gráca (ya sea en el dominio del tiempo o frecuencia), de donde un especialista inere el estado de la máquina. Por cuanto, el diagnóstico es netamente subjetivo a la experiencia de esta persona [6]. La extracción de características en el dominio del tiempo es la representación por naturaleza, sin embargo existen problemas para la diferenciación de las fallas, puesto que solo en un segmento del tiempo se encuentra alojado todo el aporte de vibración que ofrece cada mecanismo presente en la máquina. Otro espacio de representación muy utilizado es el dominio de la frecuencia, ya que es posible observar la inuencia de cada falla como una componente individual presente ligada a una amplitud especíca, pero esta representación es susceptible a los cambios de velocidad. Como la velocidad es variable, o las condiciones de operación son cambiantes debido a cargas no constantes, ha surgido la necesidad de emplear representaciones TF (Tiempo-Frecuencia), donde es posible observar los cambios de la frecuencia a medida que transcurre el tiempo [6]. En la etapa de clasicación se ha trabajado con diferentes técnicas o algoritmos de clasicación o reconocimiento de patrones. Los Modelos Ocultos de Markov (Hidden Markov Models, HMM) 2022 se constituyen en una de las herramientas más utilizadas y ampliamente extendidas. Los Modelos Ocultos de Markov ofrecen mejores resultados con menor despliegue de equipos y carga computacional; lo que se traduce en una minimización de costos y ahorro de tiempo [6]. HMM ha mostrado ser una herramienta ecaz para trabajar en la clasicación de fallas. El número de estados ocultos determina que tan efectivo u óptimo es el clasicador. Así por ejemplo, si se escoge un número pequeño de estados ocultos se disminuye el tiempo de cómputo pero puede no ser conable, y si por el contrario se elige un número elevado de estados, el tiempo de carga computacional será mayor pero el sistema mejora la conabilidad en la clasicación. 2123 4 Justicación Si se desea tener una operación conable y segura de un sistema de producción es necesario contar con una adecuada metodología de mantenimiento con el n de detectar anticipadamente posibles fallas en las máquinas y disminuir el tiempo en que una máquina quede por fuera de servicio. Es muy común encontrar en la industria que para las tareas de mantenimiento la táctica predominante es el mantenimiento preventivo, consistente en el reemplazo de piezas o elementos, o la reparación programada de los equipos en función del tiempo de operación, sin embargo, este tipo de mantenimiento no favorece la conabilidad que requiere en estos momentos la industria. El mantenimiento predictivo contribuye a tener una mayor conabilidad en los sistemas de producción, este mantenimiento cuenta con un grupo de técnicas no invasivas (análisis de aceite, análisis de desgaste de partículas, medida de vibraciones, medición de temperaturas, termografías, etc.), que tienen como nalidad pronosticar el punto futuro de rotura o avería de un componente de una máquina, de manera que dicha pieza pueda reemplazarse antes de que falle [14]. Una de las técnicas más utilizada es el análisis de vibraciones, que sirve para determinar el estado de cada uno de los componentes de los equipos con el n de programar las actividades de mantenimiento respectivas, sin afectar el desarrollo normal de la planta de producción. El controlar y monitorear las vibraciones de las máquinas tomando correctamente los datos y aplicando técnicas de análisis conocidas y vericadas, permite descubrir los diferentes tipos de problemas que se pueden presentar en las piezas de un equipo a pesar de que la falla apenas esté en su etapa inicial, además permite disminuir las pérdidas de producción, incrementar la vida útil de las máquinas, disminuir el lucro cesante, reducir los inventarios de repuestos y los paros por el uso del mantenimiento correctivo y lograr un consecuente aumento de la disponibilidad de los equipos; lo cual derivara en un notable aumento de la conabilidad [15]. Como técnicas para diagnosticar fallas en rodamientos a partir del análisis de vibraciones, principalmente se han utilizado aquellas que emplean características del dominio de la frecuencia, obtenidos a través de análisis espectrales. Y aunque, la demanda computacional que poseen las herramientas matemáticas utilizadas por las diferentes técnicas de transformación, extracción y clasicación de dichas características en el dominio espectral se pueden convertir en desventajas, 2224 cuando lo que se busca es un rendimiento computacional global del sistema, aún así, al procesar señales es común encontrar fenómenos oscilantes, casi monocromáticos, localizados en el tiempo, también ondas perdurables, que con amplitudes casi constantes, exhiben patrones de frecuencia variables en el tiempo. Surge entonces, en el análisis de señales de vibración no estacionarias, el problema de los fenómenos localizados en tiempo y frecuencia, es decir de posibles fallas que para su descripción requieren información conjunta de ambos dominios [16]. Debido a la no estacionariedad de las señales y la variabilidad de las condiciones de operación de la máquina, no es posible hacer un único análisis localizado en el tiempo tomando como base de caracterización el espectro. A partir de este hecho, se han empleado las representaciones tiempo frecuencia (TF) para el diagnóstico de fallas mecánicas, asumiendo como característica la concentración de energía y aunque ofrece buena diferenciación de las señales, la clasicación requiere de una plantilla o rma de referencia dado que la concentración de energía es relativa a la aplicación que se esté analizando [6]. Pero todo lo anterior no tiene sentido si no se implementa un sistema clasicador que pueda decidir a qué clase de falla corresponde cada una de las señales de vibración y que permita probar la ecacia y la abilidad de las características en diferentes dominios. Actualmente existen un gran número de técnicas o algoritmos de clasicación o reconocimiento de patrones, como los Modelos Ocultos de Markov (HMM), los cuales se constituyen en una de las herramientas más utilizadas y ampliamente extendidas, ofreciendo incluso mejores resultados que los arrojados por métodos relativamente más recientes como los basados en técnicas inteligentes como las redes neuronales, en donde al comparar el funcionamiento e implementación de ambas, se concluyó que los HMM ofrecen mejores resultados con menor despliegue de equipos y carga computacional, lo que se traduce en una minimización de costos y ahorro de tiempo y destacando además que no presenta incertidumbre con el resultado nal, al contrario de las redes neuronales, que para poder clasicar, dispone de categorías que dependen del número de unidades de procesamiento a usar [17]. Entendiendo que la demanda computacional que requiere la máquina para ser entrenada en la clasicación de fallas, y en clasicar estas mismas, debe ser lo menor posible con n de optimizar procesos, se hace necesario buscar alternativas para esta tarea, sin dejar a un lado la abilidad 2325 de la clasicación. Debido que a mayor número de estados utilizados para entrenar los modelos del clasicador HMM mayor abilidad en la clasicación pero a su vez mayor carga computacional, es apremiante reducir dicho número de estados sin que esto inuya de manera negativa en la abilidad del clasicador, para esto se pretende hacer una adecuada selección de las características de la señal en los diferentes espacios de representación (tiempo, frecuencia y tiempo-frecuencia) con el n de hacer extensibles los resultados obtenidos con un número menor de estados ocultos a los resultados obtenidos con un número mayor de estados, asegurando así una demanda computacional menor lo que se traduce en un proceso más óptimo. Con este trabajo se pretende vericar si los resultados obtenidos en el proceso de clasicación de fallas utilizando los HMM, con diferentes número de estados, son extensibles o equivalentes entre sí con una adecuada selección de características en los diferentes dominios (tiempo, frecuencia y tiempo-frecuencia), esto con el n de disminuir el tiempo de cómputo que le toma al clasicador trabajar, ya que si los resultados con un número de estados N, son extensibles desde los resultados con un número de estados L (con N>L), será deseable tomar decisiones de sintonización con el número de estados menor, pero nalmente entrenar con el número de estados mayor. Este trabajo se centrará en un análisis más académico e investigativo ya que este análisis comprende todo el proceso de reconocimiento y clasicación de patrones en las señales de vibración y es dividido en varias etapas de acuerdo a la teoría de aprendizaje de máquina: Espacio de representación, estimación de características, extracción de características relevantes y clasicación. El resultado nal del sistema de diagnóstico automático depende del funcionamiento secuencial y coordinado de todas las etapas, en otras palabras, si el conjunto de características no representa las propiedades intrínsecas de la señal ni genera buena separabilidad, los resultados de la clasicación no serán correctos. 2426 5 Objetivos 5.1 Objetivo General Vericar la extensibilidad de la validez de criterios de decisión para resultados obtenidos por un sistema de clasicación de fallos en rodamientos basado en Modelos Ocultos de Markov (HMM) de N estados, sobre resultados obtenidos por un sistema clasicador de la misma naturaleza con L estados (donde N > L), realizando una adecuada selección de características en diferentes espacios de representación. 5.2 Objetivos Especícos ˆ Indagar sobre las teorías más relevantes relacionada con el análisis de vibraciones para detección de fallas en rodamientos de máquinas rotativas utilizando un sistema de clasicación de fallos en rodamientos basado en HMM. ˆ Realizar segmentación y caracterización de señales de vibración provenientes de rodamientos mediante extracción de características en los dominios del tiempo, frecuencia y tiempofrecuencia; empleando para ello la base de datos del Bearing Data Center del Case Western Reserve University. ˆ Implementar un sistema de clasicación de degradación por estados, basado en Modelos Ocultos de Markov (HMM) con observación discreta. ˆ Evaluar el desempeño del sistema en el reconocimiento del estado de salud de los rodamientos a partir de señales de vibración, con el n de vericar la extensibilidad de decisiones en L estados sobre N estados, con N > L. 2527 6 Estado del Arte Con el surgimiento de la revolución industrial los procesos de producción tradicionalmente realizado por seres humanos ha venido siendo remplazado por el trabajo de máquinas, dichos artefactos tienen como objetivo optimizar los métodos de manufactura, para este n se hace necesario que estos elementos funcionen de manera que el tiempo de parada por concepto de falla sea lo menor posible, es por esto que el problema de diagnóstico de fallas ha suscitado particular interés en el área del mantenimiento trayendo consigo una larga lista de metodologías que se han desarrollado para darle solución al problema planteado. A nales del siglo XVIII el mantenimiento se basó en prácticas de prueba y corrección donde se esperaba por fallas que eran evidentes por sí mismas, posteriormente se entró a utilizar técnicas estadísticas las cuales proporcionaban una estimación sobre fallos [18]. El mantenimiento correctivo en el cual solo se interviene en los equipos cuando el fallo ya ha ocurrido no era la estrategia idónea para los procesos de producción en la época de la primera guerra mundial ya que era necesario disminuir el tiempo de producción, es así como en las industrias militares de la época se empiezan a realizar inspecciones un número concreto de piezas las cuales son sustituidas al alcanzar un número determinado de horas de funcionamiento[4]. Posterior a esto en los años 50 s un grupo de ingenieros japoneses desarrolla denitivamente el mantenimiento preventivo, el cual pretende disminuir en cierta medida la reparación mediante inspecciones periódicas y la renovación de elementos deteriorados [4]. No obstante, debido a la necesidad de reducir los costos de los métodos tradicionales surge el mantenimiento predictivo el cual estaba basado bajo la idea de conocer el estado de los equipos siendo posible reemplazar los elementos cuando estos no se encuentren en condiciones operativas y evitando averías imprevistas mediante la detección de cualquier anomalía [4]. Ya en los años 60 s aparecen las técnicas de vericación mecánica a través de análisis de vibraciones y ruidos, y así se desarrolla la previsión de fallos [19]. El primer analizador de vibración fue el cerebro humano, combinando con los sentidos del oído y del tacto, si bien es uno de los mejores cuando el sujeto está bien entrenado el analizador humano no tiene salida eléctrica y cuenta con una memoria deciente, fue por esto que se desarrollaron 2628 métodos mecánicos y electrónicos los cuales medían el nivel general o nivel de banda ancha de vibración de la máquina, posteriormente aparecieron los ltros análogos. En los años 70 s, con la aparición de la computadora personal y el analizador TRF (Transformada Rápida de Fourier) se posibilitó el cálculo de frecuencias muy rápido. En los 80 s, con la aparición del microprocesador dio cabida a la creación de un analizador de señales digital realmente portátil [20]. En el proceso de análisis de señales es necesario contar con ciertas técnicas tales como: El análisis temporal, en frecuencia, en tiempo-frecuencia, etc. Dichas técnicas han sido utilizadas en diferentes investigaciones, en [21] S. Poyhonen utilizó el modelo auto-regresivo en el dominio del tiempo en un motor de inducción utilizando los coecientes del modelo como características extraídas, en [22] Zhan and Jardine usaron modelos autorregresivos adaptativos para procesar las señales de vibración no estacionarias y se encontró que son capaces de proporcionar información able en el dominio de tiempo-frecuencia para monitorear condiciones, en [23] Kocur y Stanko utilizaron el dominio de la frecuencia cuando propusieron el orden del biespectro y alegando que permite la eliminación de envolventes y la modulación que a menudo se presenta en espectro y el biespectro de la energía convencional. Recientemente la transformada de Wavelet utilizada en la representación tiempo-frecuencia ha sido trabajada en el diagnóstico de fallas en máquinas, en elementos como: engranajes, cojinetes y máquinas de combustión interna [24]. Estos procesos de análisis de señal no son sucientes para obtener un diagnóstico adecuado de la máquina, por lo que posterior a dicho análisis se realiza el entrenamiento del clasicador de señales basado en técnicas de reconocimiento de patrones, entre estas técnicas se utilizan los HMM los cuales inicialmente en los años 70's fueron utilizados de manera incipiente en aplicaciones de procesamiento de voz [25]. 2729 Part II MATERIALES Y MÉTODOS 2830 7 Base de Datos Con el n de llevar un monitoreo del estado general de la máquina y así construir una línea de tendencia de los equipos se hace necesario contar con una base de datos en la cual estén organizado de forma sistemática los datos obtenidos en el análisis de vibraciones, debido a esto se hace menester la adquisición de un conjunto de señales de vibración medidas en sitios claves donde se encuentran localizados los rodamientos. Por consiguiente es fundamental que en dicha base de datos se encuentren muestras tomadas para distintos tipos de anomalías de tal manera que el clasicador se encuentre en la capacidad de identicar a qué tipo de falla pertenece la señal. La Western Reserve University realizó una detallada documentación del estado de fallo de los rodamientos para cada experimento para lo cual se utilizó un motor eléctrico de 2 hp de potencia. Este se conectó a un dinamómetro y a un sensor de par por medio de un acoplamiento de auto alineación. A los rodamientos se les instauraron fallas articiales mediante mecanizado por electro-descarga (EDM). Las características de los fallos son: ˆ Diámetro: desde pulgadas a pulgadas ˆ Profundidad: pulgadas. Cabe mencionar que los fallos se introdujeron por separado en el anillo interior, los elementos rodantes (es decir, bola) y la pista de rodadura exterior. Las especicaciones de los rodamientos empleados en la prueba fueron: ˆ Diámetro de la bola pulgadas ˆ Diámetro de paso pulgadas ˆ Número de bolas 9 ˆ Ángulo de contacto cero (0) Los rodamientos con fallas fueron instalados de nuevo en el motor de prueba y se registraron los datos de vibración para cargas del motor entre 0 a 3 HP con velocidades del motor de 1720 a 2931 797 RPM. La carga se aplicó utilizando el dinamómetro y los datos de potencia y velocidad se recolectaron por medio de un sensor. Las señales de vibración se adquirieron a una frecuencia de muestreo de 12 Khz. 8 Mantenimiento Se dene habitualmente mantenimiento como el conjunto de técnicas destinadas a conservar equipos e instalaciones en servicio durante el mayor tiempo posible, buscando la más alta disponibilidad y con el máximo rendimiento [26]. A pesar de la gran importancia que tiene el mantenimiento en la industria, fue solo en los años 30 que se comienza a hablar de mantenimiento cuando Henry Ford implementó en su empresa un área destinada exclusivamente a las actividades de reparación. En la actualidad con la globalización de los mercados y con la obligación que tienen las empresas en el mundo de cumplir con estándares de calidad internacional que les permita ser competitivas, es indispensable que las empresas cuenten con un adecuado plan de mantenimiento que permita detectar fallas antes de que estas ocurran y así poder también conservar sus equipos, herramientas e instalaciones en las mejores condiciones de funcionamiento [27]. En la actualidad se conoce dos tipos de metodologías de mantenimiento los cuales son correctivo y preventivo, el primero se caracteriza por ejecutarse después del reconocimiento de una falla, y destinado a llevar un elemento a un estado en el que pueda desarrollar una función requerida, por otro lado el preventivo es el mantenimiento que se ejecuta a intervalos predeterminados, y está destinado a reducir la probabilidad de fallo o la degradación de funcionamiento de un elemento. En la gura 1 se puede observar los principales tipos de mantenimiento según la UNE EN [1]. El mantenimiento predictivo a diferencia del preventivo se basa en la condición, ejecutado siguiendo una previsión del análisis y evaluación de los parámetros signicativos de la degradación del elemento, dicho de otra manera, se trata de vericar la condición de la máquina en caliente, para evitar detenerla y así disminuir los tiempos de parada para no afectar la producción [1]. Este tipo de mantenimiento consta de una serie de ensayos de carácter no destructivos, cuyo objetivo es realizar un seguimiento del funcionamiento de la maquinaria para detectar posibles 3032 Figura 1: Tipos de mantenimiento [1]. signos que adviertan un posible mal funcionamiento de esta. Dentro de los ensayos que más se utilizan están: Análisis de vibraciones, se basa en el estudio del funcionamiento de las máquinas rotativas a través del comportamiento de sus vibraciones; termografías, técnica que estudia el comportamiento de la temperatura de las máquinas; análisis por ultrasonido, está basado en el estudio de las ondas de sonido de alta frecuencia producidas por las máquinas cuando presentan algún tipo de problema; análisis de aceite, determina el estado de operación de las máquinas a partir del estudio de las propiedades físicas y químicas de su aceite lubricante, entre otros [13]. 9 Vibración Una vibración puede ser entendida como el movimiento de una partícula o de un cuerpo, o sistema de cuerpos conectados, desplazados desde una posición de equilibrio. Este movimiento se repite después de cierto intervalo; dicho intervalo se conoce como periodo de vibración. La mayoría de las vibraciones son indeseables ya que estas producen pérdidas de energía, desgaste extra, aumentan la carga, entre otras [28, 29]. Hoy en día los planes de mantenimiento predictivo de los sistemas mecánicos se basan en el estudio de las vibraciones de dicho sistema, utilizando sensores que permiten detectar formas de vibración anómalas en la máquina. Un ejemplo de un sistema vibrante en movimiento armónico simple es el movimiento que el sistema de yugo escocés le imparte a la masa m como se observa en la gura 2. En este sistema la manivela de radio A gira en torno al punto O, el extremo opuesto 3133 de la manivela P se mueve sobre una barra ranurada, la cual con un movimiento oscilatorio se desplaza en la guía vertical R. Cuando la manivela gira a una velocidad angular W, el extremo S del eslabón ranurado y por consiguiente la masa m del sistema de resorte y masa se desplazan de sus posiciones medias una distancia x (en un tiempo t) [2], dada por: x = Asen(wt) (1) El movimiento impartido a la masa m por el yugo escocés se muestra por medio de la curva senoidal de la gura 2. Figura 2: Mecanismo de yugo escocés [2]. 3234 Las vibraciones se pueden clasicar de diferentes maneras, las clasicaciones mas comunes son [2]: ˆ Vibraciones forzadas y vibraciones libres: Se dice que las vibraciones son forzadas cuando el sistema se somete a una fuerza externa a lo largo del tiempo. Las vibraciones son libres cuando no son forzadas, es decir, cuando se deja que el sistema vibre por sí mismo. ˆ Vibraciones no amortiguadas y vibraciones amortiguadas: Si se disipa o pierde energía en el movimiento, la vibración se conoce como vibración no amortiguada, no obstante, si se pierde energía se llama vibración amortiguada. ˆ Vibraciones lineales y vibraciones no lineales: Si todos los componentes básicos del sistema vibratorio son lineales, la vibración se conoce como vibración lineal, sin embargo, si cualquiera de los componentes básicos del sistema es no lineal, la vibración se denomina vibración no lineal. 9.1 Análisis de Vibraciones Debido a que en general todas las máquinas generan vibración como parte normal de su actividad, el análisis de vibraciones se convierte en la técnica de diagnóstico de falla más utilizado en la industria, permitiendo al técnico el diagnóstico de las posibles fallas mecánicas antes de que afecten irremediablemente la salud de la máquina. Las máquinas eléctricas poseen partes generadoras de movimiento como son el motor, biela, etc. Y otra parte que es la encargada de transportar la energía a los demás elementos de la máquina como por ejemplo el eje que es el que rige el funcionamiento y es la pieza principal, ya que este componente es quien determina la velocidad de giro de la máquina y se encarga de transmitir las vibraciones a los demás mecanismos como los rodamientos, ruedas dentadas, acoples y soportes [20]. Es sabido que entre las principales causas de fallo mecánico de una máquina se encuentran los efectos de deterioro producidos por fatiga de los materiales. Este proceso viene acompañado de vibraciones en la máquina, unas veces como fenómeno asociado, otras veces como directo 3335 responsable. Estas vibraciones son unas manifestaciones susceptibles de ser medida con relativa facilidad y muestran clara dependencia del estado del sistema vibratorio. Por lo tanto, midiendo y evaluando las vibraciones de la máquina se puede realizar una continua inspección del estado de esta con el n de elaborar el programa de mantenimiento respectivo en el cual pueden ir incluidos el calendario de reparaciones, las sustituciones de elementos, inuencia de máquinas próximas, etc [4]. 9.2 Técnicas de Análisis de Vibración ˆ Análisis Frecuencial: Es la técnica de análisis de vibración y detección de fallos más importante, consiste en la comparación de las frecuencias de ensayo con las frecuencias de los elementos con estado de salud mecánica sana, haciéndose posible tener una idea de las causas de un posible fallo o avería, y su magnitud. Las grácas y las comparaciones se realizan por medio de herramientas informáticas como se puede ver en las guras 3 y 4 [3]. Figura 3: Espectro de frecuencias de funcionamiento normal [3]. ˆ Análisis de forma de onda: Aquellos defectos tales como los roces del eje o de alguna parte del trotor con el estator, son los que mejor se pueden detectar con esta técnica de análisis ya que originan alguna deformación característica de la onda como se observa en la gura 5. Esta técnica es una de las más útiles para analizar fallos en el dominio del tiempo. Entre otras ventajas de esta técnica se encuentra que es posible observar componentes vibratorias de baja y muy baja frecuencia, hecho que resulta imposible medir por medio de las técnicas de análisis frecuenciales [4]. 3436 Figura 4: Espectro de frecuencias de funcionamiento defectuoso [3]. Figura 5: Onda de desplazamiento truncada por un roce [4]. 9.3 Medición de las Vibraciones Para capturar las señales de vibración se utilizan diferentes tipos de sensores, los cuales se colocan directamente sobre aquellos puntos que son más propensos a fallas. Por lo general, suelen colocarse sobre los ejes donde se ubican los apoyos de los rodamientos puesto que por éstos se transmiten las vibraciones y son una de las piezas que se dañan con mayor frecuencia. Se debe tener un cuidado muy especial a la hora de conectar los sensores de tal forma que haga un buen contacto con la estructura de la máquina para que puedan tomar lecturas en las tres direcciones y no afectar la calidad y veracidad de las medidas. En cualquier programa de monitoreo de máquinas, el hecho de que los datos sean recopilados de manera exactamente igual cada vez que se hace una medición es extremadamente importante [15]. 3537 10 Contexto Mecánico Debido a las vibraciones en las máquinas rotativas se dan diversos problemas tantos mecánicos como eléctricos que, por lo general, originan graves fallas como los son el desbalanceo, soltura mecánica, desalineación, rodamientos, entre otros. Sin embargo en la mayoría de los casos se pueden asociar los niveles de vibración con el tipo de problema que ocurre, y con la amplitud de las señales es posible predecir el nivel de severidad. Como se ha mencionado anteriormente, mediante el análisis de las señales de vibración es posible identicar posibles anomalías, por ello es recomendable realizar la toma y adquisición de estas señales en puntos diferentes de la máquina y en las direcciones vertical, horizontal y axial, ya que algunos problemas de las máquinas rotativas se pueden manifestar más en alguna dirección que en otra, por ejemplo, uno de los puntos más importantes para la medición de las vibraciones está ubicado sobre los cojinetes o cerca de ellos, debedio a que estos se transmiten las fuerzas de vibración [10]. Una vez obtenidos los datos de las vibraciones se procede a comparar con los valores que recomienda el fabricante de la máquina, y poder predecir el estado de salud de la máquina en estudio Fallas Detectadas por el Análisis de Vibraciones Como se ha mencionado anteriormente, basta con conocer las frecuencias de vibraciones y las amplitudes de las señales de vibración para conocer las posibles causas y el nivel de severidad de la falla en la máquina. En la tabla 1 se observa diferentes frecuencias encontradas en espectros de vibración de maquinarias. A continuación, se ilustra los tipos de fallas que por lo general son detectadas con el análisis de vibraciones, aclarando que para este trabajo se hará énfasis en la detención de fallas en los rodamientos de las máquinas rotativas Desalineación Considerado uno de los defectos más usuales en la industria. Ocurre cuando dos ejes conectados mediante acoplamientos no se encuentran perfectamente colineales uno con respecto al otro. 3638 FRECUENCIA CAUSAS PROBABLES DE FALLA 1xRPM Desbalanceo, engranajes o poleas descentradas, desalineamiento o eje exionado, resonancia, problemas eléctricos, fuerzas alternativas. 2xRPM Huelgos o juegos mecánicos, desalineamiento, fuerzas alternativas, resonancia. 3xRPM Desalineamiento, combinación de huelgos mecánicos axiales excesivos con desalineamiento. <a 1xRPM Remolino de aceite (oil whirl) en cojinetes o rotores. Frec. De línea Problemas eléctricos. Armónicas de RPM Engranajes dañados, fuerzas aerodinámicas, fuerzas hidráulicas, huelgos o juegos mecánicos, fuerzas alternativas Alta frecuencia(no armónicas) Cojinetes antifricción en mal estado. Tabla 1: Frecuencias de vibración y posibles causas de falla [10]. También se puede presentar entre los cojinetes de un eje sólido, o entre dos puntos cualesquiera de la máquina. Cuando una vibración se presenta a dos veces la velocidad de giro signica desalineamiento. Si la amplitud a dos veces la velocidad de giro es superior al 75% de la velocidad de giro, quiere decir que se avecina la ocurrencia de una falla [5].. Por consiguiente, la anomalía debe ser vigilada con cuidado y corregida a la primera oportunidad ˆ Desalineación Paralela: Ocurre cuando dos ejes están unidos en acople y están paralelos, pero no son colineales, es decir, desde un solo sitio se pueden observar dos puntos separados. Genera una alta vibración radial que se aproxima a un desfase de 180º en el acoplamiento. Cuando este tipo de desalineación presenta un alto nivel de severidad, se pueden producir elevados picos de amplitud a unos armónicos similares a los presentes en la holgura mecánica. El material de construcción del acoplamiento inuye de cierta forma en la totalidad del espectro cuando la desalineación es severa 3739 Figura 6: Desalineamiento paralelo [5]. ˆ Desalineación Angular: Ocurre cuando dos ejes están unidos en acople pero no están paralelos entre sí, es decir, entre los ejes existe un pequeño ángulo o lo que es lo mismo, desde un solo sitio se puede observar la presencia de una línea inclinada a partir del centro, como se observa en la gura 7. Este tipo de desalineación generará una elevada vibración axial, con un desfase de 180º a lo largo del acoplamiento. Figura 7: Desalineamiento angular [5] Desbalanceo Se conoce desbalanceo en una pieza cuando su eje inercial se encuentra desplazado del eje geométrico de rotación. Este tipo se fallo se ilustra en la gura 8. El origen de un desbalanceo puede darse por diferentes motivos como la deformación del material, desgaste de partes, mal ensamble, golpes, desprendimiento de partes, etc. Puede ser de dos tipos: estático y dinámico [6]. 3840 Figura 8: Espectro del desbalanceo y relación de fases entre soportes del eje [5] Caja de engranes Los engranes son otro elemento que puede ser analizado por medio de vibraciones, ya que se encuentran en movimiento permanente y al ser de carácter sólido permiten que la vibración que se experimenta en el eje se transmita a través de ellos. Una caja de engranes genera picos bastante prominentes en la frecuencia del engranaje (GMF), que es el número de dientes de los engranes por la velocidad de giro del engrane [6] Soltura mecánica Aparece cuando se pierden los ajustes de manufactura entre las partes que están en movimiento. La soltura mecánica se caracteriza principalmente por tener una larga cadena de armónicos de la frecuencia de giro con altas amplitudes anormales. Se clasica en tres categorías: Tipo A: Causada por soltura o debilitamiento estructural de las patas de amarre de la máquina o también por una cimentación deteriorada. Para este tipo de problema en particular, la frecuencia predominante es la de giro de la máquina (1x) y tiende a proyectar su amplitud en dirección vertical [6]. Una manera de vericar cual de las partes (tornillos de amarre, placa de montaje o cimentación) es el elemento que esta causando la soltura mecánica, consiste en hacer mediciones de fase en la dirección vertical en cada una de las partes de la base y un cambio de 90ž ó 180ž entre las partes indica cual es el elemento portador del problema. Tipo B: Este tipo de soltura genera vibración a dos veces la frecuencia de giro de la máquina (2x), y normalmente la ocasiona un aojamiento en los tornillos de amarre de la máquina, aumento 3941 Figura 9: Espectro de la soltura mecánica tipo A y relación de fases en perno o placa base [6]. en la holgura de los cojinetes de apoyo, o la presencia de una grieta o sura en el equipo y/o en la base de apoyo. Figura 10: Espectro de la soltura mecánica tipo B [6]. Tipo C: Esta clase normalmente es generada por un ajuste inadecuado entre las partes de los elementos (rodamiento - eje), que puede causar la presencia de numerosas frecuencias armónicas o subarmónicas de la velocidad de giro del eje (0.5x, 1x, 1.5x, 2x, 2.5x, 3x, etc.). Normalmente es causada por el aojamiento de un cojinete en su caja, por un espacio excesivo en la camisa o en los elementos rodantes del rodamiento Rodamientos Los rodamientos son dispositivos elementales en las máquinas rotativas. Utilizados para permitir el movimiento relativo entre dos elementos de una máquina, se suelen usar para guiar y reducir la fricción de una echa en una máquina rotatoria. Los problemas en rodamientos con elementos rodantes son las fallas más comunes que se diagnostican en análisis de vibraciones. En la mayoría 4042 Figura 11: Espectro de la soltura mecánica tipo C [6]. de los casos, ocurren debido a una lubricación inadecuada del elemento y a un montaje incorrecto, lo cual implica desalineación [6]. Estos dispositivos se encuentran formados a su vez por cuatro elementos: anillo interior, anillo exterior, jaula y elementos rodantes como se ilustra en la gura 12. Figura 12: Partes que componen un rodamiento [6]. Anillo interior y exterior: Son aquellos elementos que se jan a los dos elementos que acopla el rodamiento. Por lo general, el anillo exterior acopla con el alojamiento del rodamiento mientras que el anillo interior lo hace con el eje. Sin embargo, esta conguración podrá cambiar, pero siempre uno de los dos anillos deberá permanecer jo mientras el otro gira. Entre los dos anillos se forman unas ranuras denominadas caminos de rodadura, por donde se mueven los elementos rodantes [30]. Elementos rodantes: Los elementos rodantes (bolas, rodillos o agujas) son aquellos que giran sobre su propio eje, produciendo una rodadura de los rodamientos, permitiendo el giro relativo entre los dos aros y de la dos piezas que unen [30]. Jaula: La jaula o canastilla es la que contiene los elementos rodantes, permitiendo así mantener 4143 su posisción relativa, evitando que los elementos rodantes se desmonten o choquen entre ellos [30]. Fallas en rodamientos. ˆ Falla en pista Interna: Se produce por el desastillamiento o agrietamiento del material en la pista interna, producido por defectos de ensamble, esfuerzos anormales, corrosión, particulas externas o lubricación deciente. Se producen una serie de armónicos siendo los picos predominantes 1xRPM y 2xRPM de la frecuencia de la pista interna en dirección radial. ˆ Falla en pista externa: Se produce por el desastillamiento o agrietamiento del material en la pista externa, producido por defectos de ensamble, esfuerzos anormales, corrosión, particulas externas o lubricación deciente. Se producen una serie de armónicos siendo los picos predominantes 1xRPM y 2xRPM de la frecuencia de falla de la pista externa en dirección radial. ˆ Falla en elemento rodante: Se produce por el desastillamiento o agrietamiento del material en los elementos rodantes, producido por defectos de ensamble, esfuerzos anormales, corrosión, particulas externas o lubricación deciente. Se producen una serie de armónicos siendo los picos predominantes 1xRPM y 2xRPM de la frecuencia de falla de los elementos rodantes, en dirección radial. ˆ Falla en canastilla o jaula: Deformación de la jaula, caja o cubierta que mantiene en su posición a los elementos rodantes. Se producen una serie de armónicos de la frecuencia de la jaula siendo los picos predominantes 1xRPM y 2xRPM de la frecuencia de falla en jaula, en dirección radial o axial. 11 Espacios de Representación Una vez se ha obtenido la señal por medio de los transductores, se pasa a una etapa de procesamiento, con el n de monitorear esta señal, el principal objetivo es el de reducir la gran cantidad de información contenida en la señal de vibración a un pequeño número de características que 4244 contengan la información necesaria para describir el comportamiento general de la señal, a esto se le conoce como extracción de características de la señal. Los espacios de representación constituyen el desarrollo de la señal bajo una variable espacial determinada (tiempo, frecuencia, TF). La señal x en el dominio del tiempo es la señal de vibración adquirida del sistema e involucra la naturaleza del funcionamiento de la máquina. La señal de vibración puede ser analizada como una suma de funciones seno con distintas frecuencias f, amplitudes A, y fases ϕ, donde cada sinusoide representa un mecanismo en particular [6]. x(t) = i (A i sin(2πf i t + ϕ i )) (2) Donde el subíndice i denota cada componente o mecanismo presente en la máquina. Para pasar del dominio del tiempo a otro espacio de representación se hace uso de funciones ortogonales o funciones base que permitan transformar la señal x(t) a otro espacio Dominio del Tiempo Las técnicas del dominio del tiempo se basan en distintos comportamiento estadísticos de la forma de onda de la señal del tiempo. Esta técnica consiste en calcular los valores RMS de las señales y el factor de cresta, así como también los valores picos, pico a pico, kurtosis, desviación estándar, etc. Este análisis es un complemento al análisis de frecuencia, puesto que sirve para conrmar diagnósticos en aquellas fallas que poseen espectros muy parecidos. Estos problemas pueden ser: el desbalance, el desalineamiento y la holgura. También se utiliza cuando se presentan impactos, frotación y holgura, además de las máquinas de baja velocidad y cajas de cambio [15]. Estimación de características en el dominio del tiempo. Las características estadísticas que se presentan a continuación son consideradas entre las más usadas debido a su simplicidad para analizarlas y a la información que aportan para el diagnostico de fallas en las máquinas [6]. ˆ T1, Valor Medio (MED) El valor medio de la señal presenta problemas en la adquisición de la señal ya que es sensible a interferencia y presenta valores atípicos. Con este valor se identica el punto medio de la señal. 4345 n MED = (x n) N (3) ˆ T2, Desviación estándar (DES) Permite cuanticar el nivel de variabilidad y dispersión que presentan las amplitudes con respecto al valor medio. DES = n (x n MED) 2 N 1 (4) ˆ T3, Asimetría (ASI) La ventaja de este momento sobre el valor de kurtosis tradicional es su menor susceptibilidad a vibraciones espurias, la cual es considerada uno de los inconvenientes de los momentos estadísticos de alto orden. ASI = n (x n MED) 3 (N 1)DES 3 (5) ˆ T4, Kurtosis (KU R) Mide los picos o valles relativos de la distribución comparada a una distribución normal KUR = n (x n MED) 4 (N 1)DES 4 (6) ˆ T5, Pico positivo máximo (P V M) Útil para el monitoreo del cambio en la cantidad de impulsividad. P V M = max x n (7) ˆ T6, Raíz media cuadrática (RM S) 4446 Mide el contenido de potencia en la rma de vibración. Esta característica es muy efectiva en la detección de desbalanceo en maquinaria rotativa. RMS = n (x n) 2 N (8) ˆ T7, Factor de forma (F F O) Este parámetro permite centralizar el valor RMS para reducir la imprecisión ocasionada por problemas en la adquisición. F F O = 1 N RMS n x(n) (9) ˆ T8, Factor K (F AK) F AK = P V M RMS (10) ˆ T9, Factor T4 (T CU) T CU = (media( x(n) )) 2 (11) ˆ T10. Factor de Holgura (F CL) F CL = P V M T CU (12) ˆ T11, Factor de impulso (F IM) F IM = ˆ T12, Límite superior del histograma (HU B) 1 N max x(n) N n=1 x(n) (13) Proporciona un panorama generalizado de los valores demarcados en el extremo superior de un conjunto de muestras; permite evidenciar más fácilmente ciertas preferencias, tendencias, 4547 esquemas de comportamiento y pautas de variación de algún valor asociado, con una falla que puede ser dicil de captar con otros parámetros. HUB = max(x(n)) ( max(x(n)) min(x(n)) ) (14) N 1 ˆ T13, Límite inferior del histograma (HLB) De la misma manera como ocurre con el límite superior del histograma; este permite evidenciar comportamientos, variabilidades, tendencias y valores atípicos relacionados con una falla, que se encuentran acotados en el extremo inferior de un conjunto de muestras. HLB = min(x(n)) 0.3 ( max(x(n)) min(x(n)) ) (15) N 1 ˆ T14, Valor de la verosimilitud logarítmica negativa normal (N N L) Es utilizado como parámetro de entrada junto con las características anteriormente denidas. N NNL = log(normalpdf(f((x i, µ, σ)))) (16) n=1 Para calcular la función de distribución de probabilidad se emplea la siguiente ecuación: Normalpdf(f((x i, µ, σ))) = 1 σ 1 2π e 2 ( x i µ σ )2 (17) Donde µy σson la media y la desviación estándar respectivamente. ˆ T15, Entropía de Shanon (ENT ) ENT = entropia(x(n)) (18) ˆ T16, Raíz media cuadrática de picos (RMP ) RMP = mean( x(n) 2 ) (19) ˆ T17, Desviación Estándar mayor o igual que el umbral (0.7) (DEP ) DEP = std( x(n) ) (20) 4648 11.2 Dominio de la Frecuencia Por la naturaleza misma de la señal de vibración, la representación en el dominio de la frecuencia es la más ideal, sin embargo, para calcular la densidad espectral de potencia se hace necesario el uso de transformadas como la transformada rápida de Fourier (FFT, Fast Fourier Transform), donde la señal resultante, es un conjunto de componentes de frecuencia que comprende un rango desde la frecuencia 0 hasta fs/2, donde fs es la frecuencia de muestreo de la señal. Lo anterior es denido mediante el teorema de Nyquist, el cual dice que una señal discreta en tiempo será bien representada si la frecuencia de muestreo es igual o superior a la frecuencia máxima del sistema. Para llevar la señal en el tiempo al dominio espectral, es empleada la transformada rápida de Fourier [6]. ˆ Transformada rápida de Fourier (FFT) La operación de la Serie de Fourier está basada en una señal de tiempo que es periódica. Esto es una señal de tiempo cuya forma se repite en una cantidad innita de veces. Fourier demostró que una señal de este tipo es equivalente a una colección de funciones senos y cosenos cuyos frecuencias son múltiplos del recíproco del periodo de la señal de tiempo. El resultado un poco inesperado es que cualquier forma de onda, siempre y cuando no sea innita en longitud se puede representar como la suma de una serie de componentes armónicos, y la frecuencia fundamental de la serie de armónicos es 1 entre la longitud de la forma de onda. Las amplitudes de los varios armónicos se llaman los coecientes Fourier, y sus valores se pueden calcular fácilmente si se conoce la ecuación para la forma de onda [20]. La transformada integral de Fourier, o transformada de Fourier, se obtiene asumiendo las siguientes funciones base: ϕ(s, t) = ϕ(t, s) = e jst (21) Para t (, ) y s (, ) El parámetro s caracteriza la frecuencia de cada función base y, usualmente, se nota por alguna de las variables de frecuencia angular ω, o bien, de frecuencia lineal 2πf. Sustituyendo el valor de 4749 s en la función anterior se obtiene la representación integral de Fourier, la cual representa la señal x(t) como una suma continua de funciones exponenciales cuyas frecuencias están en el intervalo (, ). La amplitud relativa de los componentes a cualquier frecuencia f es proporcional a x(f) [20]. x(f) = + x(t)e j2φft dt f(x(t)) (22) En el área de procesamiento digital de señales, la transformada de Fourier no es posible implementarla por el hecho de ser una función continua, para lo cual existe un algoritmo conocido como la transformada rápida de Fourier, que es una aproximación numérica a la transformada de Fourier. Estimación de características en el dominio de la frecuencia. Las características estadísticas en el dominio de la frecuencia generalmente pueden indicar fallas en maquinaria mejor que las características de vibración en el dominio del tiempo, porque los componentes de frecuencia característicos como frecuencia de resonancia o componentes de frecuencia de defectos pueden ser más utiles de detectar y concurrentes con fallas. A partir de esto se dene un conjunto de características estadísticas donde x(k) es el k-ésimo valor del espectro obtenido con la FFT, k es el número total de líneas o componentes de frecuencia que tiene el espectro, donde k = 1, 2,..., K y f k es la k-ésima componente de frecuencia [6]. ˆ F1, Frecuencia Media (M ED) Representa el nivel de vibración medio de la señal en términos de la frecuencia, sin embargo no es una característica muy representativa, debido a que su valor es afectado por el ruido de la señal. MED = k (x k) K (23) ˆ F2, Varianza (V AR) Se dene como la variabilidad de la amplitud de la vibración en el espectro con respecto a la MF. 4850 V AR = K (x K MED) 2 K 1 (24) ˆ F3, Asimetría (ASI) Mide la estabilidad del espectro con referencia a la MED Sw = k (x k MED) 3 K( V AR) 3 (25) ˆ F4, Kurtosis (KU R) Útil para caracterizar fallas en rodamientos y engranajes. KUR = k (x k MED) 4 K( V ar) 4 (26) ˆ F5, Frecuencia central (M DI) Es útil como referente de normalización del rango de frecuencias presentes en la señal. MDI = k (f kx k xk (27) ˆ F6, Desviación estándar (DES) Describe el grado de convergencia del espectro de potencia. DES = k (f k F C) 2 x k K (28) ˆ F7, Raíz media cuadrática (RM S) Muestra los cambios de posición de las frecuencias principales. RMS = k (f k) 2 x k K (29) 4951 ˆ F8, Media cuadrática (M EC) MEC = media(x(n) 2 ) (30) ˆ F9, Inversa de la raíz cuadrática de la función por su media cuadrática (IV M) IV M = 1 media(x(n) 2 ) x(n) (31) ˆ F10. Desviación estándar sobre frecuencia central (DDI) DDI = DES MDI (32) ˆ Otros tipos de características utilizadas son: F 11 = K K=1 (f k F C) 3 x(k) Std 3 (33) F 12 = K K=1 (f k F C) 4 x(k) Std 4 (34) F 13 = K K=1 (f k F C) 1 2 x(k) K Std (35) ˆ F14, Raíz media geométrica (M GE) La media geométrica de una cantidad nita de números es la raíz n-ésima del producto de todos los números. k MGE = k x(k) (36) ˆ F15, Desviación media absoluta (DM A) k=1 La desviación absoluta de un elemento de un conjunto de datos es la diferencia absoluta entre ese elemento y un punto dado. Típicamente la desviación se calcula desde el valor central, se 5052 entienda como algún tipo de media, más a menudo la mediana o a veces la media del conjunto de datos. k DMA = 1 x(k) MDI (37) K k=1 ˆ F16, Rango intercuartílico (IQR) Se le llama rango intercuartílico o rango intercuatil a la diferencia entre el tercer y el primer cuartil de una distribución. Es una medida de la dispersión estadística. Los cuartiles son los tres valores que dividen al conjunto de datos ordenados en cuatro partes porcentualmente iguales. IQR = Q 3 Q 1 (38) 11.3 Dominio Tiempo-Frecuencia En el análisis de vibraciones es muy común encontrar que las señales se representen en el dominio del tiempo s(t) y/o en el dominio de la frecuencia s(ω) a través de su transformada de Fourier. Ambos dominios contienen la misma información sobre la señal estudiada, sin embargo, cada una de estas señales brindan enfoques distintos y complementarios. Partiendo de que la señal es no periódica y de energía nita, estas representaciones se relacionan mediante el par de fórmulas o par de Fourier: ŝ(ω) = + s(t)e jωt dt (39) s(t) = 1 2π + ŝ(ω)e jωt dω (40) Donde t representa el tiempo y ω la frecuencia angular, indicando que, en teoría, la información en uno de ellos (t ó ω) puede recuperarse a partir de la información entregada en el otro [31]. Ahora bien, existen fenómenos ubicados en tiempo y frecuencia los cuales no es suciente con abordarse con uno de los dos dominios naturales sino que por el contrario necesitan información de 5153 ambos dominios, por tal motivo, el par de Fourier se queda corto, haciendo necesario la utilización de patrones básicos capaces de ubicar conjuntamente la información de ambos dominios. Lo que plantea el problema de las representaciones tiempo-frecuencia; el desafío consiste en denir una apropiada representación acorde con las características de la señal y a los objetivos de su procesamiento. Las distribuciones o transformadas tiempo-frecuencia se pueden clasicar en lineales y cuadráticas. ˆ Distribución de Wigner-Ville. Las distribuciones cuadráticas se caracterizan porque dependen en forma cuadrática de la señal temporal permitiendo interpretarlas como una distribución de energía de la señal en el plano tiempo-frecuencia. Dentro de este tipo de distribuciones la más conocida es la Distribución de Wigner-Ville [24]. W V (t, f) = ˆ Transformada de Fourier de Tiempo Reducido (ST F T ) S(t + τ 2 )S(t τ 2 )e j2πft dt (41) Una de las distribuciones lineales extensamente utilizadas es la llamada Transformada de Fourier de Tiempo Reducido (STFT, Short Time Fourier Transform), para el análisis de señales donde sea necesario realizar un estudio tanto de las frecuencias como del tiempo, y detectar así las características espectrales y temporales [6]. El funcionamiento básico de las STFT consiste en que la función ventana inicialmente está localizada al inicio de la señal, es decir en t = 0. Si se asume que la anchura de la venta es T segundos, que están siendo escogidos de la señal. Una vez que se ha hecho esto, la nueva señal es el producto de función ventana y la señal original a la que se la aplica la transformada de Fourier [7]. El resultado de esta transformación es la transformada de Fourier de los primeros T/2 segundos de la señal original, y si esta parte de la señal es estacionaria quiere decir que los resultados obtenidos mostrarán la representación en frecuencia exacta de los primeros T/2 segundos [7]. 5254 Posteriormente el siguiente paso será desplazar esta ventana a una nueva localización hasta que toda la señal sea recorrida. Matemáticamente, se escribe como: ST F T (τ, f) = t+ τ 2 t τ 2 x(t)g(t τ)e j2πft dt (42) Donde g(t) es la función ventana, comúnmente una ventana de Hann o colina Gausiana es centrada en cero y x(t) es la señal a ser transformada. Gracamente, la ST F T se puede explicar de la siguiente manera, ver gura 13: Figura 13: Explicación gráca de la STFT [7]. En la gura 13, se muestra una función ventana de tipo Gaussiana, existen tres funciones o tres ventanas localizadas en t = t 1 en t = t 2 y t = t 3. Estas ventanas corresponden a tres transformadas de Fourier en tiempos distintos. Por lo tanto se obtendrá una buena representación tiempo-frecuencia de la señal en dichos intervalos de tiempo. El soporte de la ventana constituye un parámetro de gran importancia ya que a través de este se puede establecer el grado de resolución tanto de tiempo como de frecuencia que se desea. Si la ventana es muy angosta se analizará una porción muy pequeña de la señal lo que permitiría tener una buena resolución en tiempo pero una mala resolución en frecuencia ya que se conocería sólo una mínima fracción del espectro total existente en la señal. Por otro lado, si la ventana en muy ancha se tendría una buena resolución en frecuencia pero una mala resolución en tiempo. Por lo tanto, el problema fundamental de la STFT consiste en la incertidumbre de Heisenberg 5355 el cual establece que es imposible conocer una representación exacta tiempo- frecuencia de una señal, es decir, no se puede saber qué valor de frecuencia existe en un instante de tiempo determinado, sólo se podrá conocer que componentes de frecuencia existen dentro de un intervalo de tiempo determinado [7]. ˆ Transformada Wavelet (W T ) La Transformada Wavelet (WT, Wavelet Transform), es una distribución lineal, presenta ventajas en cierto tipo de señales donde un análisis de porcentaje de banda constante es más apropiado que el análisis de ancho de banda jo que utiliza la STFT [24], ya que este tipo de herramienta matemática es capaz de concentrarse en fenómenos transitorios y de alta frecuencia. La transformada Wavelet tiene un tamaño de ventana adaptado a las frecuencias. La transformada Wavelet, descompone la señal en una serie de pequeñas ondas u ondículas, las cuales son ondas localizadas que descienden rápidamente a cero después de unas pocas oscilaciones, y que tienen un valor promedio nulo [32]. Se dene la Wavelet respuesta al impulso como una nueva Wavelet madre. Una vez denida las funciones de comparación, se puede denir la transformada Wavelet como la suma en el tiempo del producto de la señal por la versión escalada y trasladada de la Wavelet original [32]. Matemáticamente se expresa de la siguiente manera. W T (t, a) = 1 a x(τ)w ( t τ )dτ (43) a Antes de continuar, es importante aclarar que las características estadísticas en el dominio de las frecuencia y del tiempo presentadas anteriormente solo son válidas si se asume que la señal en estudio es estacionaria es decir solo son útiles bajo régimen permanente. Estimación de características para Melcepst. Si se desea hacer un diagnóstico más acertado de la salud de la máquina, es necesario considerar a la señal de vibración en un estado no estacionario, bien sea por regímenes de carga y velocidad variable o por los transitorios causados por problemas en rodamientos. Las características estadísticas en el dominio de las frecuencia y del tiempo presentadas anteriormente solo son 5456 válidas si se asume que la señal en estudio es estacionaria es decir solo son útiles bajo régimen permanente, por tal motivo, se hace necesaria la utilización de características en el dominio TF que permitan una visualización dinámica de la señal de vibración en función de sus componentes frecuenciales a través del tiempo. Sin embargo este dominio presenta un problema, ya que los mapas TF presentan una alta dispersión de la información de la señal, por lo que denir zonas donde se encuentren información discriminante acerca de los tipos de falla es difícil [6]. ˆ Centroides de subbanda espectral (SSC) Con el n de denir los centroides de subbanda espectral, se divide la banda de frecuencia (desde0 afs/2 donde fs es la frecuencia de muestreo) en un número jo de subbandas, y se calcula el centroide de cada subbanda usando el espectro de potencia de la señal [24]. ˆ Coecientes cepstrales (LF CC) Es una técnica adelantada que consiste en tomar un espectro de un espectro. Antes de calcular el cepstro, se calcula el logaritmo natural de la amplitud del espectro [20]. Estimación de características para transformada de Fourier de tiempo corto (ST F T ) Las características usadas para para la transformada de Fourier de tiempo corto son: ˆ STFT1, Kurtosis (KU R) KUR = k k=1 (x(k) MDI)4 K (44) k=1 k( (x(k) MDI)2 k 1 ) 3 ˆ Otras Características ST F T 2 = 1 2 (45) K K=1 x(k)2 ˆ STFT3, Varianza (V AR) V AR = k k=1 (x(k) k k=1 x(k) k ) 2 k 1 (46) ˆ STFT4, Pico Máximo (M AX) 5557 MAX = max(x(k)) (47) Estimación de las características para la transformada Wavelet (T W ) Para este análisis se tiene en cuenta que la descomposición de Wavelet es de dos niveles, por lo que se analiza cada nivel por separado, para obtener más características de la misma señal con los mismos cálculos estadísticos. ˆ TW1, Media (MED) Consiste en sumar todos los datos y dividirlos por el número total de datos de las componentes de alta frecuencia. MED = K K=1 x 1(k) K (48) ˆ TW2, Pico Positivo Máximo (P V M) Valor de una función dada, que es mayor o igual que cualquier valor de la función dada. El máximo absoluto es el mayor de todos los valores. P V M = max( x 1 (k) ) (49) ˆ TW3, Kurtosis (KU R) Es una medida de la forma. Así, las medidas de Kurtosis tratan de estudiar la proporción de la varianza que se explica por la combinación de datos extremos respecto a la media en contraposición con datos poco alejados de la misma. KUR = k k=1 (x 1(k) MED) 4 k (50) k=1 k( (x 1(k) MED) 2 K 1 ) 3 ˆ TW4, Energía de Shannon (ENT ) 5658 K K=1 ENT = (x 1(k) x 1 (k)) 2 (51) Las siguientes características son las utilizadas para las componentes de baja frecuencia de la transformada Wavelet: T W 5 = k k=1 x 2(k) K (52) T W 6 = max( x 2 (k) ) (53) T W 7 = k k=1 (x 2(k) T W 1) 4 k (54) k=1 k( (x 2(k) T W 1) 2 k 1 ) 3 T W 8 = k k=1 (x 2(k) x 2 (k)) 2 (55) 12 Segmentación y Ventanas 12.1 Segmentación Una vez que se obtiene la información útil de la señal de vibración y con el n de asumir estacionalidad, se procede a realizar la segmentación de la señal, el cual es un procedimiento que consiste en dividir la señal en un número N de muestras consecutivas llamados frames. Cabe mencionar que la segmentación se obtiene multiplicando la señal x(n) por la ventana V (n) que pondera las muestras de dicho segmento [33, 34]. Un aspecto a tener en cuenta en esta etapa es el solapamiento o traslape entre frames consecutivas, las cuales deben guardar relación entre sí para que el análisis entre trama no sea aislado, así se pueden presentar tres situaciones con el traslape de los tramos: ˆ Si N N se enmascararán las muestras contiguas y el espectro resultante estará correlacionada trama a trama. 5759 ˆ Si N N la resolución en frecuencia sera demasiado baja. ˆ Si N > N no existirá traslape entre los frames adyacentes, perdiéndose parte de la señal de vibración. La gura 14 representa una segmentación en tramas para el caso en el que el traslape ( N) sea igual a la mitad del número de muestras. Figura 14: Segmentación de la señal de vibración [6] Ventaneo En la etapa de segmentación se realiza la división de la señal en tramas, esto provoca un efecto indeseable en el análisis de la evolución en el tiempo de las características de la señal. Con el n de solucionar este problema se realiza la etapa de ventaneo el cual consiste en aplicar una función ventana haciendo que los bordes de la señal tiendan a cero suavizando estos y resaltando la parte central lo cual acentúa las propiedades características de cada segmento. La expresión matemática de la aplicación de una ventana es la siguiente [11]: S i (n) = S i (n) V (n) (56) donde S i (n) es la señal segmentada, i es el índice que representa el comienzo de cada segmento, 5860 V (n) es la función ventana aplicada, y N es el tamaño en número de muestras del segmento. Características de las ventanas Existen cuatro parámetros a tener en cuenta en el momento que se dene un tipo de ventana, estos parámetros son [11]: 1. Ancho de banda a 3 db: Es la diferencia entre las frecuencias en las que su atenuación, al emplear la ventana, se mantiene igual o menor a 3 db cotejada con la frecuencia central del pico máximo. Este parámetro da testimonio acerca de la capacidad que tiene la ventana para separar componentes de similar amplitud de manera que, determina la resolución del análisis. 2. Selectividad: Es un parámetro que revela la capacidad de la ventana para discirminar componentes de niveles abundantemente disímiles. El parámetro principal para conocer la selectividad es el factor de forma, el cual surge de la correspondencia entre el ancho de banda a 60 db y a 3 db. 3. Altura del lóbulo-lateral principal: Los lóbulos laterales con altas magnitudes causan efectos de rizado en la respuesta frecuencial del frame de señal ventaneado. 4. Magnitud del rizado en la banda de paso: La banda de paso se encuentra denada como f centrada alrededor de la frecuencia central f 0. Este rizado dene la capacidad de la ventana para representar la magnitud real de las componentes frecuenciales. La tabla 2 muestra algunas funciones ventanas con sus caracteristicas, mientras que la tabla 3 muestra algunas otras ventanas con sus ecuaciones respectivas. 5961 Nombre Ancho de Banda a 3 db Factor de Forma Altura del Lóbulo- Lateral Principal Magnitud del Rizado en la Banda de Paso Rectangular 0.89Df db 3.92 db Gauss 1.37Df db 1.58 db Hanning 1.44Df db 1.40 db Chebyshev 1.85Df 3.8 Puede ser selecionado 0.87 db Kaiser-Bessel 1.71Df db 1.02 db Flat Top 3.72Df db 0.01 db Tabla 2: Funciones ventana y sus características [11]. 6062 Ventana Espectral Rectangular Ecuación w R (n) = 1 Tukkey-Hamming w T H (n) = { cos(nπ), n 1 0 P ara.otro.valor Daniel Von Hann (Hanning) w H (n) = w D = { Sen(nπ) nπ { Cos(nπ), n 1 0 P ara.otro.valor { } Bartlett-Priestley w BP (n) = 3 Sen(nπ) (nπ) 2 nπ Cos(nπ) 1 6(n) 2 + 6( n ) 3, n 0.5 Parzen w P (n) = 2(1 n ) 3, 0.5 < n 1 0, n > 1 { Cos(2nπ) Cos(4nπ), n 0.5 Blackmann w B (n) = 0, P ara.otro.valor I o(β(1 [2n] 2 1 ) 2 ) Kaiser w K (n) = I o(β), n 0.5 0, P ara.otro.valor Donde I o (x) es una función de Bessel de orden cero modicada del primer tipo y puede ser evaluada para cualquier grado de exactitud I o (x) = 1 + L k=1 [ 1 k! ( x 2 )k ] 2 Empíricamente L<25 Tabla 3: Ventanas espectrales usadas en el procesamiento digital de datos [11]. 13 Análisis de Componentes Principales PCA Cuando se recoge información de un conjunto de datos, se busca trabajar con el mayor número posible de variables, sin embargo, si se toman demasiadas variables sobre un conjunto de observaciones se tendrá que considerar un mayor número de coecientes de correlación lo que dicultaría la visualización de las relaciones entre las variables. En el análisis de componentes principales lo único que sabemos previamente es que el 100% de la 6163 variabilidad de las N variables se explica por N factores, cada uno de los cuales es combinación lineal de las variables originales. Evidentemente, no se pretende sustituir las N variables por N factores, lo cual sería impráctico, lo que interesa de esto es que los factores recogen la variabilidad de las variables originales de forma desigual, es decir, pocos factores recogen un porcentaje de variabilidad alto, por lo que, se puede explicar la mayor parte de la variabilidad a partir de ellos [35]. Dicho de otro modo, el problema que plantea resolver esta herramienta consiste en que dado el caso de disponer de los valores de p características en n ventanas de una señal de vibración dispuesto en matriz x de dimensiones nxp, donde en las columnas se encontrarán las características (variables) y en las la las ventanas (observaciones), se deseará encontrar un espacio de representación más reducido que el conjunto de datos original Cálculo de las Componentes Principales Se considera una serie de variables (x 1, x 2,..., x p ) sobre un grupo de objetos o individuos y se trata de calcular, a partir de ellas, un nuevo conjunto de variables (y 1, y 2,..., y p ) incorreladas entre sí, cuyas varianzas vayan decreciendo progresivamente. Cada y j donde j = 1,..., p es una combinación lineal de las x 1, x 2,.., x p originales, es decir y j = a j1x1 + a j2x a jpxp = a j x siendo a j = (a 1j, a 2j,..., a pj ) un vector de constantes, y x = [x 1...x p ] La primera componente se calcula eligiendo a 1 de modo que y 1 tenga la mayor varianza posible, sujeta a la restricción de que a 1 a 1 = 1. El segundo componente principal se calcula obteniendo a 2 de modo que la variable obtenida y 2 esté incorrelada con y 1. Del mismo modo se eligen y 1, y 2,..., y p, incorrelados entre sí, de manera que las variables aleatorias obtenidas vayan teniendo cada vez menor varianza [36]. En general, la matriz original de datos X y por lo tanto la matriz de varianza y covarianza S, tiene rango p, existiendo entonces tantas componentes principales como variables que se obtendrán calculando los valores propios o raíces características, λ 1,..., λ p, de la matriz de varianzas y covarianzas de las variables S, mediante: 6264 S λi = 0 (57) Y sus vectores asociados son: (S λ i I)a i = 0 (58) Los términos λ i son reales y positivos. Por ser S simétrica, y λ j y λ h son dos raíces distintas sus vectores asociados son ortogonales. a h Sa j = (a h Sa j) = a jsa h (59) a h Sa j = a jλ h a h (60) Y siλ j λ h, y a h a j = a j a h = 0 son ortogonales. Llamando Y a la matriz cuyas columnas son los valores de los p componentes en los n individuos, estas nuevas variables están relacionadas con las originales mediante: Y = XA (61) donde A A = I Interpretación de las Componentes Calcular las componentes principales equivale a aplicar una transformación ortogonal a las variables X (ejes originales) para obtener unas nuevas variables Y incorreladas entre sí. Esta operación puede interpretarse como elegir unos nuevos ejes coordenados, que coincidan con los ejes naturales de los datos [36]. Cuando existe una alta correlación positiva entre todas las variables, la primera componente principal tiene todas sus coordenadas del mismo signo y puede interpretarse como un promedio 6365 ponderado de todas las variables. Los restantes componentes se interpretan como factores de forma y típicamente tienen coordenadas positivas y negativas, que implica que contraponen unos grupos de variables frente a otros [36] Selección del número de componentes Principales Se han sugerido distintas reglas para seleccionar el número de componentes, entre ellas se encuentra la técnica de seleccionar componentes hasta cubrir determinada proporción de varianza. Esta regla debe aplicarse con cierto cuidado. Por ejemplo, es posible que una sola componente recoja el 90% de la variabilidad y sin embargo pueden existir otras componentes que expliquen adecuadamente la forma de la variable. También se puede calcular el número de componentes realizando un gráco λ i frente a i o desechar las componentes asociadas a valores propios inferiores a una cota, que suele jarse como la varianza media, λ i /P [36]. 14 Clasicación Una parte fundamental para el diagnóstico de fallas es el procesamiento de la información con el n de determinar qué tipo de información es de carácter redundante, distorsionada e imprecisa. El objetivo principal de los sistemas de clasicación es discriminar las señales de vibración, identicando los datos de entrada dentro de un número denido de clases, ya sea cuando se tiene un conjunto de datos de fallos tipicados o estandarizados llamado clasicador multiclase, los cuales están conformados por dos etapas: análisis y clasicación de características, o bien cuando, se tiene un conjunto de registros conocidos y abundantes, en el mejor de los casos registros de la máquina en un estado normal y se construye una frontera de decisión donde sea posible clasicar registros nuevos y desconocidos, a esta clasicación se le conoce como de una clase o mono clase [6]. Particularmente, en este trabajo se tratan los problemas de clasicación de fallas bajo el concepto de clasicación multiclase. En la etapa de análisis se parametriza y se extraen las características más relevantes de la señal de vibración con el n de disminuir la carga computacional, en la 6466 segunda etapa, se procesa las características extraídas y mediante un proceso de clasicación se determina el estado de salud del rodamiento, ya sea, falla en anillo interior, anillo exterior, elementos rodantes o rodamiento sin falla. De acuerdo al modelo de señal a utilizar para caracterizar las propiedades de la señal de vibración, estos se clasican en modelos deterministas y modelos estadísticos. Como consecuencia de la necesidad de caracterizar exclusivamente las propiedades estadísticas de la señal, se han desarrollado un gran número de procesos estadísticos como los procesos gaussianos, los procesos poissonianos, los clasicadores bayesianos simples (naivebayes), los discriminadores lineales, los procesos de Markov o los procesos ocultos de Markov [37]. Modelos Ocultos de Markov (HMM) Las cadenas de Markov son modelos estadísticos de señales observables representadas como un proceso aleatorio paramétrico, es decir, los parámetros pueden determinarse con precisión. Un modelo oculto de Markov corresponde a una cadena de Markov (q) en la cual la observación (O) de salida es una variable aleatoria generada de acuerdo a una función de probabilidad asociada a cada estado [25]. En los modelos ocultos de Markov, la ruta de estados no es observable, pero la salida, dependiendo del estado, si lo es. Cada estado posee una distribución de probabilidad sobre las posibles señales de salida. Por lo tanto, los modelos ocultos de Markov ofrecen indicaciones sobre el orden de los estados con ayuda de la secuencia de símbolos que estos generan [38]. Un HMM está caracterizado por lo siguiente: ˆ Matriz de Transición A = [(aij)nxn] Es una matriz (matriz de transición de estados) conformada por las probabilidades de transición donde a ij es la probabilidad de transición desde el estadoi al estado j, para todo i, j N. a ij = P (S n = j S n 1 = i) (62) ˆ Matriz de Probabilidades B = [bi(k)] 6567 Es un vector de probabilidades de emisión de símbolos (matriz de confusión), uno por cada estado, donde b j = (b j1, b j2,..., b jm ) es la probabilidad de emisión del símbolo O k del alfabeto en el estado j. b i (k) = P (x n = o k S n = i) (63) ˆ Estados Ocultos Q = [1, 2, 3,..., N] Es un conjunto nito de N estados (estados ocultos) que representan los estados (reales) del sistema que pueden ser descritos por el proceso de Markov. ˆ Estados Observables Σ = [O 1, O 2, O 3, , O M ] Es un conjunto discreto nito de M símbolos (estados observables) que representan los estados visibles del proceso, es decir, el alfabeto. ˆ Estados Iniciales π = [π 1, π 2,..., π N ] Es un vector (distribución inicial de estado) de probabilidades del estado inicial q o en Q. La notación compacta λ = (A, B, π) se emplea para representar los modelos ocultos de Markov, en otras palabras, precisar un HMM involucra la escogencia de un número proporcional de estados N, y símbolos discretos M, y especicando las probabilidades A, B, π. Las probabilidades de iniciación, transición y emisión son los parámetros del modelo. 15 Libro de Código (Codebook) El objetivo de utilizar el libro código es reducir el tamaño de la información a manejar, buscando la mejor representación de los vectores característicos que en su estado original pueden tomar innitos valores. Dicho de otro modo, la idea es encontrar y sustituir la representación de los vectores característicos de cada unos de los dominios los cuales pueden tomar innitos valores, por un conjunto nitos de vectores que constituyen lo que se denomina un libro de código, donde una señal es 6668 cotejada con cada patrón en un cierto intervalo, y la señal en cada período es proyectada por un código que índica el patrón que es más parecido a la señal [8]. La forma de medir la delidad de un cuanticador es determinar el error que éste produce al reemplazar los datos de entrada que recibe por los vectores representantes o codewords, dicho parámetro es llamado error por distorsión. La nalidad de un cuanticador es obtener un conjunto de vectores representativos llamado codebook, que presente el menor error por distorsión [8]. Durante la etapa de entrenamiento del clasicador se deben generar un conjunto de patrones de referencia apoyándose en un grupo de vectores de entrenamiento, teniendo presente que el conjunto de patrones de referencia debe ser menor que el grupo de vectores de entrenamiento. Al conjunto de patrones de referencia se los denomina libro de código y a los patrones, palabras de código. Figura 15: Construcción del libro de código [8]. La etapa de construcción del libro de código se inicia con un gran número de vectores de observación, [V 1, V 2,...V L ] que conforman el grupo de entrenamiento, posteriormente se construye el grupo de vectores de referencia o del libro de código (codebook) óptimo que representa la variabilidad observada en el grupo de entrenamiento. Inicialmente, se requiere medir la distancia entre cada par de vectores de observación con el n de agrupar el conjunto de vectores de entrenamiento como así también para clasicar vectores arbitrarios a cada entrada del libro de código. Posteriormente, un algoritmo de clasicación basado en clustering (por ejemplo, K-means), permite ubicar y calcular los centroides de los vectores de observación, agrupando estos en torno a K valores medios arrojando como resultado un libro de códigos o codebook [8]. En conclución, se puede decir que los N vectores originales se representa por M vectores, cada uno de los cuales es llamado palabra de código o codeword, luego el grupo entero de dichos vectores, forma un libro de códigos o codebook, formando entonces M regiones. En la gura 16 6769 se observa un diagrama de Voronoi con sus correspondientes regiones, conformando un codebook con sus correspondientes codewords como centroides. Figura 16: Diagrama de Voronoi [8]. Existen dos maneras de generar el libro de código, las cuales son: 1. Aprendizaje Aleatorio: Permite a los vectores ser elegidos de manera aleatoria de los datos de entrenamiento y luego ser guardados como vectores de códigos. 2. Clúster: Permite agrupar los datos de entrenamiento en grupos mutuamente separados y a cada grupo se le calculan las medias permitiendo minimizar la distorsión promedio, estas medias son archivadas como vectores de código. Clustering El proceso Clustrering consiste en agrupar datos en clases de tal modo que los elementos de un clúster posean un nivel elevado de semejanza entre ellos, y un nivel bajo con elementos de otros clúster [9]. Los algoritmos que permiten clasicar un conjunto de vectores de entrenamiento en un grupo de clases distintas (o clúster) se denominan algoritmos de clustering y se agrupan en dos categorías [9]: 1. Algoritmos de partición: Consiste en dividir el conjunto de observaciones en k conglomerados o clases (clúster), en donde k lo dene inicialmente el usuario. 6870 2. Algoritmo jerárquicos: Consiste en entregar una jerarquía de divisiones del conjunto de elementos en conglomerados. Existe una gran variedad de algoritmos de clustering como: K-means, CobWeb, Algoritmo EM (Expectation Maximization), SOM (self-organization maps), entre otros. Este trabajo está enfocado en el estudio del algoritmo de K-means. K-Means Este algoritmo de clustering es uno de los mas utilizado por su eciencia. Consiste en un procedimiento de clasicación de un grupo de objetos en una cantidad de K números de clústers. Este algortimo es usado comúnmente para situaciones en las cuales todas las variables son cuantitativas y la distancia cuadrática Euclidea es elegida como medida de diferencia. d(x i, x i ) = p (x ij x i j) 2 = x i x i 2 (64) j= Se puede observar que los pesos en la distancia Euclidea, pueden ser usados redeniendo los valores x ij. Los puntos de dispersion se denen de la siguiente manera: k W (C) = x i x i 2 (65) k=1 c(i)=k c(i )=k W (C) = Donde k k=1 N k x i x i 2 (66) c(i )=k x k = ( x 1k, x 2k,..., x pk) Es el vector de medias asociado con el k-ésimo clúster, y N k = n I(C i = K i=1 Así, el criterio es asignar las N observaciones a los K clústers de modo que dentro de cada 6971 conglomerado el promedio de las diferencias de cada observación a la media del clúster, denido por los puntos del clúster, sea mínima [39]. Cuando la representación es llevada a cabo utilizando centroides, la inmediatez en la entrega de resultados se convierten en una gran ventaja gracias a su representación de naturaleza gráca y estadística. Cada clúster por tanto es caracterizado por su centroide, que se encuentra en el medio de los elementos que componen el clúster[9], como se observa en la gura 17. Figura 17: Clustering empleando k-means [9]. El proceso de agrupamiento K-means se realiza de la siguiente manera: Como primer paso se dene el número K de conjuntos y se asume el centro de esos grupos. La denición de los centroides se puede realizar tomando k objetos de forma aleatoria como los K clústers iniciales o tomando los primeros K objetos en forma secuencial. Posterior a esto, se dene el ó los centroides iniciales con base al número de clúster esperado, luego se determina la distancia de cada objeto con relación a los centroides y nalmente se agrupan los objetos con base en la distancia mínima [9]. Los pasos del algoritmo de clustering empleando K-means puede resumirse de la siguiente manera, como se ilustra en la gura 18 7072 Figura 18: Algoritmo de clustering empleando K-means [9]. Part III METODOLOGÍA 7173 16 Desarrollo Experimental El presente trabajo profundiza en la detección de fallas en rodamientos utilizando clasicadores basados en HMM con diferentes estados, para lo cual se realiza una caracterización en tiempo, frecuencia y tiempo-frecuencia de las señales de vibración obtenidas a partir de una prueba del Bearing Data Center de la Western Reserve, el cual creó una base de datos de las señales de vibración de los rodamientos normales y defectuosos. Los rodamientos son piezas mecánicas que debido a su facilidad para ser cambiadas se diseñan de tal manera que estas fallen primero, por este motivo el continuo monitoreo del estado de dichas piezas es imperioso bajo un modelo de mantenimiento predictivo Base De Datos Las muestras tomadas fueron almacenadas en archivos con extensión.mat, agrupadas de acuerdo al estado de salud de los rodamientos: fallo en bola (Ball), fallo en la pista externa (Outer), fallo en la pista interna (Inner) y estado normal (Normal). Por su parte cada falla esta subdividida en tres niveles de severidad (N1,N2,N3), cada nivel de severidad correspondiente a cada tipo de fallo se encuentra compuesto por doce registros correspondientes a cuatro velocidades diferentes, en el cual cada registro almacena un arreglo bidimensional de 40000X Segmentación y Ventaneo Como proceso previo a la extracción de características es importante dividir la señal de vibración de ingreso en segmentos de igual tamaño (tramas o frames) lo cual permite darle estacionalidad a las señales de vibración y distinguir las partes de señal que contiene información relevante de aquella que posee la información redundante, de este modo se evita realizar cálculos inútiles que conllevan a incrementar la carga computacional y a incrementar las probabilidades de conseguir resultados erróneos en el análisis de las señales de vibraciones. Las señales de vibración se encuentran divididas en N segmentos cada uno de longitud L que se traslapan entre sí. La selección de la longitud de los segmentos se planteó de forma detallada con el n de garantizar que las tramas posean la suciente información como para obtener 7274 Figura 19: Distribución base de datos características localizadas de la señal de vibración de manera que el tiempo de cálculo sea el menor posible y que dichas características permanezcan invariantes. Dependiendo del tipo de representación se decide si se utiliza el ventaneo ya que con esto se evita que haya discontinuidades al principio y al nal de los bloques analizados, lo cual dependiendo de la representación espacial que se utilice puede afectar el análisis; en el momento de seleccionar la ventana se debe ser bastante cuidadoso ya que una ventana bastante grande puede perjudicar la captación de cambios relevantes en la señal de vibración, en este trabajo se utilizó una ventana hamming representada por la siguiente ecuación: υ(n) = ( 2πn N 1 ) (67) 16.3 Selección de Características Vale la pena señalar que la etapa de selección de características va a pesar en gran medida en el desempeño y rendimiento nal de la etapa de clasicación, así una reducción en la dimensionalidad del espacio de características que entran al clasicador sin que la información global disminuya de forma signicativa resultará benecioso en el proceso de clasicación, por este 7375 motivo se aplica métodos multivariantes (PCA) como etapa intermedia entre la extracción de características y el proceso de clasicación, de esta manera se organiza en una única matriz todas las características extraídas a todas las señales del grupo de entrenamiento y se le aplica PCA, posteriormente se calculan las componentes principales para cada trama de cada señal y se reorganiza el nuevo vector de características el cual va a estar compuesto por las componentes principales que representan las características que se tenían originalmente. Una vez se ha segmentado la señal, sobre cada una de las tramas de las señales de vibración originales se evaluaron las siguientes características estadísticas que representan el comportamigento de la señal en diferentes dominios de representación Tiempo Las siguientes características mostradas en la tabla 4 son utilizadas para el dominio temporal y están explicadas en la sección 11. T1, Valor Medio T10. Factor de holgura T2, Desviación Estándar T11, Factor de impulso T3, Asimetría T12, Limite superior del histograma T4, Kurtosis T13, Límite inferior del histograma T5, Pico positivo máximo T14, Valor de la verosimilitud logarítmica negativa normal T6, Raíz Media Cuadrática T15, Entropía de Shanon T7, Factor de forma T16, Raíz media cuadrática de picos T8, Factor K T17, Desviación estándar mayor o igual al umbral T9, Factor T4 Tabla 4: Características en el dominio del tiempo Frecuencia Las siguientes características mostradas en la tabla 5 son utilizadas para el dominio de la frecuencia y están explicadas en la sección 11. Para trabajar con estas características se debe primero realizar una transformación de su espacio a través de las transformada rápida de Fourier Tiempo frecuencia Las siguientes características son utilizadas para el dominio de tiempo frecuencia y están explicadas en la sección 11. Se dividen en tres tipos: 7476 F1, Frecuencia media F9, Inversa de la RMS F2, Varianza F10. Desviación estándar sobre frecuencia central F3, Asimetría F11 F4, Kurtosis F12 F5, Frecuencia central F13 F6, Desviación estándar F14, Raíz media geométrica F7, Raíz media cuadrática F15, Desviación media absoluta F8, Media cuadrática F16, Rango intercuartílico Tabla 5: Características utilizadas en el dominio de la frecuencia ˆ Transformada de Fourier de tiempo corto (STFT) STFT1, Kurtosis STFT3, Varianza STFT2 STFT4, Pico máximo Tabla 6: Características utilizadas en el dominio tiempo-frecuencia STFT ˆ Transformada de Wavelet Características para las componentes de alta frecuencia TW1, Media TW2, Valor Máximo TW3, Kurtosis TW4, Energía de Shanon Características para las componentes de baja frecuencia TW5, Media TW6, Valor Máximo TW7, Kurtosis TW8, Energía de Shanon Tabla 7: Características utilizadas en el dominio tiempo-frecuencia Wavelet ˆ Transformada Melcepst Centroides de sub-banda espectral Coecientes cepstrales Tabla 8: Características utilizadas en el dominio tiempo-frecuencia Melcepst 16.4 Clasicación El sistema clasicador esta basado en Modelos Ocultos de Markov (HMM) para el cual la selección de la arquitectura y de los parámetros del clasicador es fundamental. Para esto se hace necesario que la base de datos sea dividida en dos partes, la primera parte conocida como conjunto de entrenamiento es utilizada para determinar los parámetros del clasicador utilizando un 60% de los datos, la segunda parte es el conjunto de validación la cual se encuentra conformada por un 7577 40% de los datos, cabe mencionar que es posible trabajar con otros porcentajes de división de los datos muestrales. El software para cada fase del clasicador se basa en un lenguaje matemático el cual ofrece un ambiente de desarrollo completo con un lenguaje de programación integro como lo permite MATLAB, este posee toolbox los cuales son librerías dedicadas al entrenamiento y clasicación de modelos usando Modelos Ocultos de Markov (HMM), dichas librerías están conformados por pequeñas funciones y scrips que permiten vericar paso a paso en qué condiciones se encuentra el sistema. La validación cruzada es la herramienta para evaluar los resultados de los análisis realizados por el entrenamiento a partir de los HMM, de esta manera se garantiza que son independientes de la partición entre datos de entrenamiento y prueba Evaluación Los conceptos de sensibilidad y especicidad de una prueba son utilizados para determinar la validez en el sistema, la sensibilidad de una prueba responde a las preguntas de: ¾Cuántos resultados positivos se obtendrán en los elementos con falla? y ¾Cuántos casos del total de casos en la base de datos estudiada pueden identicarse por el resultado de la prueba? Mientras que la especicidad de una prueba responde a las preguntas de: ¾Cuántos resultados negativos se obtendrán en elementos sin falla? y ¾Cuántas veces el elemento sin falla se conrmará por el resultado de la prueba? [22]. Para realizar este análisis se diseña una arreglo(nxn), conocido como matriz de confusión donde el número de clase es representado por n, esta conguración permite que una vez se allá ejecutado el algoritmo se puedan visualizar los resultados. Las columnas de la matriz representan los casos en una clase real y las las representan los casos que el algoritmo predice. Los datos de la matriz de confusión son los siguientes: Tasa de Verdaderos Positivos (tvp) O Sensibilidad Es el cociente entre el número de registros con falla correctamente clasicados y el número total de datos con falla. Son los casos que pertenecen a la clase y el clasicador los denió en esa clase. 7678 tvp = vp vp + fn (68) Tasa de Falsos Negativos (tf n) Es el cociente entre el número de registros con falla incorrectamente clasicados y el número total de datos con falla. Son los casos que si pertenecen a la clase y el clasicador no los denió en esa clase. tfn = fn vp + fn (69) Tasa de Verdaderos Negativos (tvn) También llamado especicidad, es el cociente entre el número de registros normales correctamente clasicados y el número total de datos normales o sin falla. Son los casos que no pertenecen a la clase y el clasicador denió que no pertenecen a esa clase. tvn = vn vn + fp (70) Tasa de Falsos Positivos (tf p) Es el cociente entre el número de registros normales incorrectamente clasicados y el número total de datos normales o sin falla. Son los casos que no pertenecen a la clase pero el clasicador los denió en esa clase. tfp = fp vn + fp (71) Exactitud (e) Proporción del total de predicciones correctas. e = vp + vn vp + vn + fp + fn (72) 7779 17 Sintonización y Validación A partir de la base de datos proporcionados por el Bearing Data center de la Western Reserve Univsersity [40] se realizan varias pruebas para los cuatro estados de salud de la máquina: fallas en anillo interior (Inner), fallas en elementos rodantes (ball), fallas en anillo exterior (Outer) y un nivel base (normal). Con el n de comprobar la extensibilidad de criterios de decisión para un clasicador de fallos en rodamientos basado en HMM, sobre resultados en un clasicador de menos estados, donde se emplea espacios de representación tales como tiempo, frecuencia, Wavelets y la transformada de tiempo corto de Fourier; vericando el desempeño general de las fases de entrenamiento y validación de los HMM para todos los estados con los que se entrenó. Con los datos que se cuentan se realiza un fase de extracción de características y estas son divididas en dos partes, 60% de la base de datos como un conjunto de entrenamiento de tal manera que sea posible determinar los parámetros del clasicador, mientras que con el 40% restante se utiliza como conjunto de validación con el n de vericar el correcto comportamiento del modelo. Es de particular interés conocer que todas las pruebas se realizaron a una velocidad de muestreo de [Ks/seg], con un tamaño de ventana de 0.2 ms, aplicando una ventana Hamming, utilizando un traslape de 1/3, con 24 centroides, 12 coecientes cepstrales y 24 ltros para la caracterización con melcepstrum. Con este n se realiza una serie de pruebas que permitan comprobar la extensibilidad de criterios de decisión para clasicadores basados en HMM con diferentes estados; estas pruebas son: 17.1 Prueba 1: Evaluación de cada espacio de representación con todas sus características La prueba se realizó con la nalidad de observar cómo se comportaba la eciencia y el tiempo de cómputo del sistema de clasicación de fallas evaluando cada espacio de representación (Tiempo, Frecuencia, Melcepst, Wavelet y Transformada en Tiempo Corto de Fourier), con todas sus características a medida que se iban aumentando los estados de los HMM utilizados. Cada entrenamiento se realizó 30 veces y al nal se promedió los resultados para obtener valores medios de las variables de interés. 7880 17.2 Prueba 2: Combinación de espacios de representación con todas sus características La prueba se realizó con la nalidad de observar cómo se comportaba la eciencia y el tiempo de cómputo del sistema de clasicación de fallas cuando se combinaban los espacio de representación (Tiempo, Frecuencia, Melcepst, Wavelet y Transformada en Tiempo Corto de Fourier), con todas sus características a medida que se iban aumentando los estados de los HMM utilizados. Cada entrenamiento se realizó 30 veces y al nal se promedió los resultados para obtener valores medios de las variables de interés. Cabe mencionar que para los 5 espacios de representación resultaron 29 combinaciones de espacios para entrenar ya que la combinación donde se entrena un solo espacio corresponde a lo hecho en la prueba 1 y la combinación donde se entrena con ningún espacio no tiene sentido realizarla Prueba 3: Evaluación de las mejores características por espacio de representación Con el n de observar como es el comportamiento del clasicador cuando se entrena con cada característica de los diferentes espacios de representación a medida que aumenta los estados de los HMM del clasicador, se entrenó dicho clasicador con cada una de las características 30 veces, para 1,2,3,4,5,6,7 y 8 estados, promediando al nal la eciencia de los 30 entrenamientos para tener los valores medios de dicha variable y condensando la información en un gráco de barras para cada espacio de representación y para cada número de estados ocultos, el cual muestra la eciencia en la clasicación para cada característica entrenada de forma individual Prueba 4: Evaluación del número de componentes principales para cada espacio de representación Con el n de conocer cuál era el número de componentes principales que mantenía el 95% de la varianza acumulada, se realizó esta prueba para los 5 espacios de representación que se trabajaron, la cual consiste en realizar un gráco de valores propios de la matriz de covarianza corregida cumulativos vs número de componentes principales, en este se observa que cantidad 7981 de componentes permiten mantener un 95% o más de la variabilidad acumulada de las variables que originalmente se tenían Prueba 5: Representatividad de cada característica Con el n de conocer cuales característica eran las más relevante para cada espacio de representación se calculó la relevancia para cada característica y se gracaron estas en un gráco de barras donde se pudiera apreciar claramente cuáles eran sus relevancias. Este gráco se realizó para los 5 espacios de representación Prueba 6: Evaluación de cada espacio de representación con todas sus características utilizando PCA La prueba se realizó con la nalidad de observar cómo se comportaba la eciencia y el tiempo de cómputo del sistema de clasicación de fallas evaluando cada espacio de representación (Tiempo, Frecuencia, Melceps, Wavelet y Transformada en Tiempo Corto de Fourier), con todas sus componentes principales a medida que se iban aumentando los estados de los HMM utilizados. Cada entrenamiento se realizó 30 veces y al nal se promedió los resultados para obtener valores medios de las variables de interés Prueba 7: Combinación de espacios de representación con todas sus características utilizando PCA La prueba se realizó con la nalidad de observar cómo se comportaba la eciencia y el tiempo de cómputo del sistema de clasicación de fallas cuando se combinaban los espacio de representación (Tiempo, Frecuencia, Melceps, Wavelet y Transformada en Tiempo Corto de Fourier), con todas sus componentes principales a medida que se iban aumentando los estados de los HMM utilizados. Cada entrenamiento se realizó 30 veces y al nal se promedió los resultados para obtener valores medios de las variables de interés. Cabe mencionar que para los 5 espacios de representación resultaron 29 combinaciones de espacios para entrenar ya que la combinación donde se entrena un solo espacio corresponde a lo hecho 8082 en la prueba Evaluación de cada espacio de representación con todas sus PCA y la combinación donde se entrena con ningún espacio no tiene sentido realizarla. 8183 18 Pseudocódigos 18.1 Función Frecuencia Este algorítmo se encarga de recorrer la longitud total de la señal y extraer un segmento de esta de longitud N. Aplica una ventana Hamming al segmento extraído y posteriormente se aplica la transformada rápida de Fourier a la señal ventaneada, una vez realizado esto se procede a calcular las diferentes caracteristicas para el dominio de la frecuencia para nalmente ser almacenadas en un vector. Algorithm 1: Función Frecuencia Data: s, N, DeltaN, desact s: Señal; N: Número de muestras por ventana; deltan: Número de muestras en traslape; desact: Desactivación de características; begin Lt Longitud(s); while Recorre la señal desde 1 : N DeltaN : Lt do w Desde donde inicia la correspondiente trama hasta donde naliza, extraiga un segmento de señal de longitud N; Actualice el paso como punto de inicio de la trama N DeltaN; h Hamming(Longitud(w)); x w h; W abs(t(x)); MED Media(W); V AR Varianza(W); ASI Asimetria(W); KUR Kurtosis(W); MDI MedianFrecuency(W); DES DesviaciónEstándar(W); RMS RaízMediaCuadrática(W); MEC MediaCuadrática(W); IV M Inversa(MEC(W ) RMS(W )); DDI DES/MDI; F 11; F 12; F 13; MGE MediaGeométrica(W); MAR MediaArmónica(W); DMA DesviaciónMediaAbsoluta(W); Cmp [MED, VAR,ASI,KUR,MDI,DES,RMS,MEC,IVM,DDI,F11,F12,F13,MGE,MAR,DMA]; Desactive la característica indicada Result: cms 8284 18.2 Función STFT Este algorítmo se encarga de recorrer la longitud total de la señal y extraer un segmento de esta de longitud N, elimina la tendencia lineal al segmento extraído y posteriormente calcula el espectograma a la señal linealizada, una vez realizado esto se procede a calcular las diferentes caracteristicas para el dominio tiempo-frecuencia para nalmente ser almacenadas en un vector. Algorithm 2: Función STFT Data: s, N, DeltaN, desact s: Señal; N: Número de muestras por ventana; deltan: Número de muestras en traslape; desact: Desactivación de características; begin lt Longitud(s); while Recorre la señal desde 1 : N DeltaN : lt do w 0 Desde donde inicia la correspondiente trama hasta donde naliza, extraiga un segmento de señal de longitud N; Actualice el paso como punto de inicio de la trama N DeltaN; w Elimine tendencia lineal (w 0 ); W D Calcule Espectograma(w); [fil, col] Tamaño(W D ); for i = 1 : col do caract(1, i) Kurtosis(W D (i)); caract(2, i) 0.5 Suma(W D (i) 2 ); caract(3, i) Variaza(W D (i)); caract(4, i) ValorMáximo(W D (i)); Carac1 Media(carac(1)); Carac2 Media(carac(2)); Carac3 Media(carac(3)); Carac4 Media(carac(4)); cmp Concatene[Carac1,Carac2,Carac3,Carac4] Desactive la característica indicada Result: cms 8385 18.3 Función Tiempo Este algorítmo se encarga de recorrer la longitud total de la señal y extraer un segmento de esta de longitud N. Una vez realizado esto se procede a calcular las diferentes caracteristicas para el dominio del tiempo del segmento extraído para nalmente ser almacenadas en un vector. Algorithm 3: Función Tiempo Data: s, N, DeltaN, desact s: Señal; N: Número de muestras por ventana; deltan: Número de muestras en traslape; desact: Desactivación de características; begin lt longitud(s); while Recorre la señal desde 1 : N DeltaN : Lt do w Desde donde inicia la correspondiente trama hasta donde naliza, extraiga un segmento de señal de longitud N; Actualice el paso como punto de inicio de la trama N DeltaN; MED Media(w); DES DesviaciónEstándar(w); ASI Asimetría(w); KUR Kurtosis(w); P V M PicoMáximo(w); RMS RaízMediaCuadrática(w); F F O FactorForma(w); F AK FactorK(w); T CU t4; F CL ClearenceFactor(w); F IM ImpulseFactor(w); HUB LímiteSuperiordelKistograma(w); HLB LímiteInferiordelHistograma(w); NNL NormalNegativeLog-Likelihood(w); ENT EntropíaShannon(w); RMP ValorRMSPicos(w); DEP DesviaciónEstándar>=UmbralTrama; Cmp Concatene[MED,DES,ASI,KUR,PVM,RMS,FFO,FAK,TCU,FCL,FIM,HUB,HLB,NNL,ENT,RMP,DEP]; Desactive la característica indicada Result: Cms 8486 18.4 Función Entrenamiento y Validación Este algorítmo procede a cargar las características extraídas de la señal y guardarla en un vector, posteriormente divide la base de datos en dos partes, una para pruebas y la otra para entrenamiento, nalmente ejecuta las funciones para libro de códigos, para HMM y validación. Algorithm 4: Funcion Entrenamiento y Validación Data: NumCent, NumEst, P orent, rutaraiz NumCen: Número de centroides para agrupamiento; NumEst: Número de estados en los HMM; PorEnt: Porcentaje de base de datos para entrenamiento; rutaraiz: Ruta actual de directorios; begin Cargar V ibcaract; k 0; ef t 1; X V ibcaract; divisionbd(x, P orent)=[x train, X prueba ]; while ef t > do k k + 1; librocodigos2(x train, NumCen) = [V ibcodes, K, dist]; hmmtrain2(x train, V ibcodes, NumEst, 100.1e 5) = [A md, B md, pi md, logp m, ytrain]; V alid(x prueba, A md, B md, pi md, V ibcodes) = [Modelo t, Efic t, ef t]; Result: V ibcodes, ModelosF allas, V alidacin 8587 18.5 Función K-Means Este algorítmo se encarga de generar números aleatorios desde 1 hasta el numero de columnas de los datos a ser agrupados, luego calcula las distancias euclideas entre los datos a ser agrupados y los primeros números aleatorios y asigna cada vector de los datos al centroide mas cercano. Algorithm 5: Función K-Means Data: Y, K, maxiter Y: Datos a ser agrupados; K: Número máximo de clusters o grupos; maxiter: Número máximo de iteraciones; begin [M, N] Tamaño(Y ); if K > M then error, mostrar ("Mas centroides que vectores de datos"); perm Generar números aleatorios permutados de 1 hasta M; Y c Guardar los k primeros números de perm; d2y Sumar(Y 2 ); for i = 1 : maxiter do Y cold Y c; Result: Y c, c, err d2 Cálculo distancia euclidea entre Y y Y c; Asignar cada vector en Y al centroide mas cercano; Ajustar los nuevos centroides de acuerdo a la nueva asignación; err(i) Cálculo error en la iteracción i; if (Y c Y cold)< then Return; Actualizar err; 8688 18.6 Función Libro de Código Este algorítmo determina los tipos de fallas a reconocer almacenandolos en una sola matriz, la cual es presentada a la función K-Means para nalmente calcular la distorsión promedio. Algorithm 6: Función Libro Código Data: data, kmax Data: Matriz de características; Kmax: Número máximo de centroides; begin R Tipos de fallos a reconocer; for i = 1 : R do L Longitud(data(i)); for j = 1 : L do y Acceder a la carpeta i señal j; Y Concatenar todas las señales de todas la carpetas; Agrupamiento de vectores y generación de centroides; [Y c, c, err] = kmeans(y, Kmax, 5000); cb Clústers que tienen asignación; [N, K] Tamaño(cb); [N, T ] Tamaño(Y ); dist DistorsiónPromedio; Result: cb, K, dist 8789 18.7 Función Entrenamiento y Prueba Diagnóstica con PCA Este algorítmo se encarga de hacer una combinacion con los cinco espacios de representación posibles, para cada combinación recorre los 8 estados iterativamente, para el estado i-ésimo carga la base de datos, el usuario debe indicar algunos parámetros de entrada, ejecuta la función de extracción de caracteristicas y validación, se almacena el valor de la eciencia para posteriormente calcular la matriz ROC. Algorithm 7: Principal Multiple Componentes (Entrenamiento y prueba diagnostica) Data: begin for cada combinación posible de los 5 espacios do Car Crear la combinación de espacios; for NumEst = 1 : 8 do BaseDatos Se carga la base de datos; Fs Indicar velocidad de muestreo; T v Indicar tamaño de la ventana; Overlap Indicar traslape de ventana; PorEnt Indicar porcentaje de entrenamiento; ShowROC Seleccionar si se muestra curva ROC al nal; CPRIN Seleccionar PCA; Las veces que se realizan los entrenamientos es la dimensión de vector NumCoef y cada elemento del vector debe ir con el mismo valor del número de coecientes con el que se desea entrenar; NumCoef Indicar en este vector el número de coecientes y las veces que se realizarán los entrenamientos para este número de coecientes; NumF ilt Indicar en este vector el número de ltros; Las veces que se realizan los entrenamientos es la dimensión de vector NumCen y cada elemento del vector debe ir con el mismo valor del número de coecientes con el que se desea entrenar; NumCen Indicar en este vector el número de centroides; for ico = 1 : Longitud(NumCoef) do for ifi = 1 : Longitud(NumFilt) do VibExtractCaract(Fs,Tv,Overlap,NumCoef,NumFilt,rutaRaiz,BaseDatos,Car,desactC,CPRIN); for ice = 1 : Longitud(NumCen) do V ibentrev alidhmm(numcen(ice), NumEst, P orent, rutaraiz); Cargar Validación; Efic t I Se extrae de validación la eciencia del entrenamiento para cada tipo de fallo; Calcular Datos de la MatrizROC; Calcular la MatrizROC; if ShowROC then Mostrar CurvaROC; E t ime Calcular tiempo de cómputo; Eficiencia Cálculo de eciencia total de la clasicación; Eficiencia i nd Cálculo de las eciencias medias por tipo de falla; Crear carpeta con el nombre del espacio y el número de estados que se esta trabajando; Guardar MatrizROC; Guardar Efic ti ; Guardar E time ; Guardar Eficiencia; Result: MatrizROC, Efic ti, E time, Eficiencia 8890 18.8 Función Entrenamiento y Prueba Diagnóstica Este algorítmo primero carga la base de datos, el usuario debe indicar algunos parámetros de entrada, ejecuta la función de extracción de caracteristicas y validación, se almacena el valor de la eciencia para posteriormente calcular la matriz ROC. Se le da la orden de construir un gráco de barras de las eciencias de los entrenamientos trabajando con una sola característica. Algorithm 8: Principal Varios (Entrenamiento y prueba diagnostica) begin BaseDatos Se carga la base de datos; Fs Indicar velocidad de muestreo; T v Indicar tamaño de la ventana; Overlap Indicar traslape de ventana; PorEnt Indicar porcentaje de entrenamiento; ShowROC Seleccionar si se muestra curva ROC al nal; Car Indicar los dominios de representación con los que se va a entrenar; CPRIN Seleccionar PCA; Las veces que se realizan los entrenamientos es la dimensión de vector NumCoef y cada elemento del vector debe ir con el mismo valor del número de coecientes con el que se desea entrenar; NumCoef Indicar en este vector el número de coecientes y las veces que se realizarán los entrenamientos para este número de coecientes; NumF ilt Indicar en este vector el número de ltros; Las veces que se realizan los entrenamientos es la dimensión de vector NumCen y cada elemento del vector debe ir con el mismo valor del número de coecientes con el que se desea entrenar; NumCen Indicar en este vector el número de centroides; for NumEst = 1 : 8 do for i = 1 : nmerodecaractersticasdelespacio do Desactivar característica iesima; desactc Guarda característica desactivada; for ico = 1 : Longitud(NumCoef) do for ifi = 1 : Longitud(NumFilt do VibExtractCaract(Fs,Tv,Overlap,NumCoef,NumFilt,rutaRaiz,BaseDatos,Car,desactC,CPRIN); for ice = 1 : Longitud(NumCen) do V ibentrev alidhmm(numcen(ice), NumEst, P orent, rutaraiz); Cargar Validación; Efic ti Se extrae de validación la eciencia del entrenamiento para cada tipo de fallo; Calcular MatrizROC; Calcular Datos de la MatrizROC; Guardar MatrizROC; if ShowROC then Mostrar CurvaROC; E time Calcular tiempo de cómputo; Eficiencia Cálculo de eciencia total de la clasicación; Construya gráco de barras de las eciencias de los entrenamientos trabajando con una sola característica; Crear carpeta con el nombre del espacio y el número de estados que se esta trabajando; Guardar grácos de barras en las carpeta creadas; 8991 18.9 Función Wavelet Este algorítmo se encarga de recorrer la longitud total de la señal y extraer varios segmentos de esta de longitud N, elimina la tendencia lineal al segmento extraído y posteriormente normaliza cada uno de los segmentos, aplica descomposición de Wavelets de 2 niveles (Baja y Alta frecuencia) y limita las ventanas para ambos niveles, nalmente inicia el proceso de calcular las características necesarias a cada nivel. Algorithm 9: Función Wavelet Data: s, N, DeltaN, desact s: Señal; N: Número de muestras por ventana; deltan: Número de muestras en traslape; desact: Desactivación de características; begin Lt Longitud(s); while Recorre la señal desde 1 : N DeltaN : Lt do w Desde donde inicia la correspondiente trama hasta donde naliza, extraiga un segmento de señal de longitud N; Actualice el paso como punto de inicio de la trama N DeltaN ; w 0 ElimineTendenciaLineal(w 0 ); media Media(w o); w Normalice el nuevo segmento de señal; Realice descomposición de Wavelets de 2 niveles; Establecer desplazamiento para la ventana componentes alta frecuencia ; Establecer desplazamiento para la ventana componentes baja frecuencia ; Establezca el límite superior de la ventana para las componentes de alta frecuencia; Establezca el límite superior de la ventana para las componentes de baja frecuencia; for i = 1:19 do Establezca el nuevo límite superior para la ventana i,para las componentes de alta frecuencia; Establezca el nuevo límite superior para la ventana i, para las componentes de baja frecuencia; ventana Ventana i, para las componentes de alta frecuencia; ventana2 Ventana i, para las componentes de baja frecuencia; temp(1, i) PromedioNormalizado(ventana); temp(2, i) Máximo(ventana); temp(3, i) Kurtosis(ventana); temp(4, i) EnergíaShannon(ventana); temp(5, i) PromedioNormalizado(ventana2); temp(6, i) Máximo(ventana2); temp(7, i) Kurtosis(ventana2); temp(8, i) EnergíaShannon(ventana2); Calcule inicio para ventana i + 1 para las componentes de alta frecuencia; Calcule inicio para ventana i + 1 para las componentes de baja frecuencia; Carac1 Media(temp(1, :)); Carac2 Media(temp(2, :)); Carac3 Media(temp(3, :)); Carac4 Media(temp(4, :)); Carac5 Media(temp(5, :)); Carac6 Media(temp(6, :)); Carac7 Media(temp(7, :)); Carac8 Media(temp(8, :)); cmp [carac1,carac2,carac3,carac4,carac5,carac6,carac7,carac8]; Desactive la característica indicada; Result: cms 9092 18.10 Función Validación Esta función se encarga de recorrer desde 1 hasta el número de señales que se utilizan para la prueba, luego ejecuta la funcion de clasicación HMM y nalmente cálcula la eciencia de la clasicación. Algorithm 10: Funcion Validación Data: X prueba, A m, B m, P i m, cb begin nclases Longitud(X prueba ); for i = 1 : nclases do nsec Tamaño(X prueba ); Ejecutar hmmrecog2(x prueba (i), A m, B m, pi m, cb); Efic(i) Eciencia fallo i; ef t Eciencia total de la clasicación; Modelo(i) Almacena el modelo fallo i; Result: Modelo, Efic, ef t Función Extracción de Características En primer lugar se determina la base de datos a utilizar, se cálcula el tamaño de la ventana y el traslape, posteriormente se extraen las características de cada segmento para los diferentes espacios de representación las cuales son almacenadas en un vector de características. El algoritmo permite elegir si se desea trabajar con PCA, de ser así los pasos son muy similares a los descritos en el proceso anterior con la variación que los vectores característicos deben ser procesados antes de ser almacenados, para ello se construye una matriz de pesos la cual contiene en sus primeros componentes la información necesaria para describir los datos. 9193 Algorithm 11: Función Extracción de características Data: F s, T v, Overlap, NumCoef, NumF ilt, rutaraiz, BaseDatos, Opc, DesactC, CP RIN Fs: Frecuencia de muestreo; Tv: Tamaño de ventana; Overlap: Fracción de ventana en traslape; NumCoef: Número de coecientes cepstrales(12); NumFilt: Número de ltros(24); rutaraiz: ruta actual de directorios; BaseDatos: Nombre de la base de datos a emplear; Opc: Espacios de representación; desactc: Desactivación de característica; CPRIN: Selección de PCA; begin bdvibra Cargar la base de datos; V Longitud(bdvibra); V ibcaract celda(1,v ) Celda de características; N redondear(f s T v); DeltaN redondear (N overlap); for i = 1 : V do K Longitud(archivosdevibracin); for j = 1 : K do s Obtener datos archivo de vibración de ubicación (i, j); if Selecciono Tiempo then desactt desactc(1); Ct Función caracttiempo(s, N, N deltan, desactt ); [numvent, fila] Tamaño(Ct); if Selecciono Frecuencia then desactf desactc(2); Cf Función caractfrec(s, N, N deltan, desactf ); [numvent, fila] Tamaño(Cf); if Selecciono Melcepst then desactm desactc(3); Cm Función melcepst; [numvent, fila] Tamaño(Cm); if Selecciono Wavelets then desactw desactc(4); Cw Función caractwavelet(s, N, N deltan, desactw ); [numvent, fila] Tamaño(Cw); if Selecciono STFT then desacts desacts(5); Cs Función caractstft(s, N, N deltan, desacts); [numvent, fila] Tamaño(Cs); if selecciono PCA then Ctaux(i, j) Ct ; Cfaux(i, j) Cf ; Cwaux(i.j) Cw ; Cmaux(i.j) Cm ; Csaux(i.j) Cs ; else V ibcaract(i)(j) [Ct, Cf, Cm, Cw, Cs ]; Result: V ibcaract 9294 Algorithm 12: Continuación Función Extracción de características begin for i = 1 : V do if Selecciono PCA then V ibcaractprinc celda(1,v ); for i = 1 : V do V ibcaractprinc(i) celda(k(i), 1); for j = 1 : K do if Selecciono Tiempo then componentest 4; rt MatrizDePesosPCA(Ctaux); for ut = 1 : componentest do for lt = 1 : numvent do ct Componentes ut para la ventana lt; Ct Almacene(cT ) ubicación (ut, lt); if Selecciono Frecuencia then componentesf = 4; rf MatrizDePesosPCA(Cfaux); for uf = 1 : componentesf do for lf = 1 : numvent do cf Componentes uf para la ventana lf ; Cf Almacene(cF ) ubicación (uf, lf ); if Selecciono Melcepst then componentesm = 7; rm MatrizDePesosPCA(Cmaux); for um = 1 : componentesm do for lm = 1 : numvent do cm Componentes um para la ventana lm; Cm Almacene(cM) ubicación (um, lm); if Selecciono Wavelets then componentesw 3; rw MatrizDePesosPCA(Cwaux); for uw 1 : componentesw do for lw = 1 : numvent do cw Componentes uw para la ventana lw; Cw Almacene(cW ) ubicación (uw, lw); if Selecciono STFT then componentess 2; rs MatrizDePesosPCA(Csaux); for us = 1 : componentess do for ls = 1 : numvent do cs Componentes us para la ventana ls; Cs Almacene(cS) ubicación (us, ls); V ibcaractprinc(i)(j) [Ct, Cf, Cm, Cw, Cs]; V ibcaract V ibcaractprinc; Guardar V ibcaract; 9395 Part IV RESULTADOS, ANÁLISIS Y CONCLUSIONES 9496 19 Prueba 1: Evaluación de cada espacio de representación con todas sus características. En las tablas 9, 10, 11, 12 y 13 se encuentran los datos obtenidos de la prueba 1 realizada con el n de observar el comportamiento de la eciencia de la clasicación de fallas evaluando cada espacio de representación con todas sus características a medida que aumentan los estados de los HMM, la estructura de cada tabla muestra en las las el número de estados del clasicador y en las columnas la eciencia, el tiempo y la desviación estándar respectivamente. A partir de estas tablas se muestra de manera resumida los resultados en la gura 20. En cada una de las guras se gracan el tiempo de cómputo, eciencias y desviaciones vs el número de estados. La línea coloreada con azul muestra el comportamiento que tiene la eciencia, la línea de color verde corresponde al tiempo y las lineas verticales de color azul hacen referencia a la desviación estándar. La gura 20 (a) corresponde a los resultados arrojados para la Frecuencia, la gura 20 (b) muestra los resultados para Melcepst, la gura 20 (c) para STFT, la gura 20 (d) para el dominio del tiempo y la gura 20 (e) para Wavelet. Número de Estados Eciencia Tiempo Desviación Estándar PROMEDIO Tabla 9: Resultados prueba 1 Frecuencia 9597 Número de Estados Eciencia Tiempo Desviación Estándar PROMEDIO Tabla 10: Resultados prueba 1 Melcepst Número de Estados Eciencia Tiempo Desviación Estándar PROMEDIO Tabla 11: Resultados prueba 1 STFT Número de Estados Eciencia Tiempo Desviación Estándar PROMEDIO Tabla 12: Resultados prueba 1 Tiempo 9698 Número de Estados Eciencia Tiempo Desviación Estándar PROMEDIO Tabla 13: Resultados prueba 1 Wavelet (a) (b) (c) (d) 9799 (e) Figura 20: Resultados prueba 1 Se puede observar en la gura 20 (a) que la eciencia para el dominio de la frecuencia tiende a disminuir conservando una tendencia constante a medida que aumenta el número de estados. Por otro lado, el tiempo de cómputo del clasicador tiende a aumentar a medida que se incrementa el número de estados, por otro lado, la desviación estándar presenta un comportamiento ascendente a medida que aumenta el número de estados. Para la gura 20 (b) la eciencia en el dominio de Melcepst tiende a aumentar muy poco a medida que se aumentan el número de estados. El tiempo de cómputo muestra una tendencia no lineal presentado su mayor pico en el estado 5 a partir del cual comienza a disminuir nuevamente, de lo anterior se puede inferir que este espacio de representación presenta la mejor eciencia en la clasicación pero generando un tiempo de cómputo mayor. La desviación estándar tiende a aumentar a medida que se incrementa el número de estados presentándose un valor muy alto para el primer estado. El comportamiento de la eciencia para el dominio STFT el cual está representado en la gura 20 (c) presenta una relación inversamente proporcional al número de estados del clasicador, es decir, que a medida que aumenta el número de estados disminuye la eciencia. Por otro lado el tiempo de cómputo es ascendente. Cabe destacar que este espacio presenta la menor eciencia en la clasicación. En cuanto a la desviación estándar esta tiende a aumentar en cuanto se 98100 incrementa el número de estados, sin embargo para el primer estado los valores de las desviación son muy altos. Con referencia a la gura 20 (d) se puede destacar que la eciencia para el dominio del Tiempo presenta un comportamiento descendente a medida que se aumentan el número de estados presentando su mayor valor en el estado número cinco. Adicionalmente el tiempo de cómputo y la desviación estándar tienden a aumentar con el número de estados. En la gura 20 (e) se observa que la eciencia para el dominio Wavelet presenta una tendencia constante ya que a medida que aumenta el número de estados del clasicador los valores de la eciencia no varían signicativamente, por otra parte el tiempo de cómputo tiene un comportamiento exponencial el cual tiende a aumentar con el número de estados. Los valores para la desviación estándar incrementan a medida que se aumenta el número de estados presentándose un pico muy alto para el primer estado. Al evaluar cada espacio de representación con sus respectivas características se concluye que la eciencia del clasicador presenta un comportamiento descendente a medida que se incrementa el número de estados, sin embargo, se presenta un caso atípico con la transformada Melcepst en donde la eciencia aumenta levemente siendo este espacio el más eciente en la clasicación con un 93,6% en promedio, como se observa en la tabla 10,con el hecho de que el tiempo de cómputo tiende a ser mayor a medida que se aumentan el número de estados, esto se debe a que cuando se incrementa el número de estados con los demás espacios se está sobre entrenando los modelos, mientras, que cuando se incrementan los espacios para la transformada de Melcepst si se presenta aprendizaje. Por otra parte, el espacio de representación menos eciente al momento de realizar la clasicación es la transformada en tiempo corto de Fourier STFT el cual tiene una eciencia de 78,6%, como se puede ver en la tabla 11. Se destaca que el dominio del tiempo y frecuencia son los espacios de representación que ofrecen una eciencia en la clasicación más alta 91,5% y 89,5% en promedio respectivamente y su tiempo de cómputo no es tan elevado como Melcepst, ver tabla 12 y 9. Igualmente se destaca que la desviación estándar para cada uno de los espacios de representación tiende a aumentar a medida que se incrementa el número de estados, además se observa que para las transformadas Wavelet, Melcepst y STFT el valor para el primer estado de las desviaciones 99101 estándar tiende a ser mayor en comparación a los demás estados. 20 Prueba 2: Combinación de espacios de representación con todas sus características Los datos obtenidos de la prueba 2 se encuentran en la tabla 14 representando los estados en las las y en las columnas, los valores de la eciencia y desviación estándar para cada una de las combinaciones de los espacios de representación. Con estos datos se ha construido una serie de grácos donde se compara como se comporta la eciencia del sistema de clasicación a lo largo de los estados, cuando se combinan todos los espacios de representación (línea de color rojo en la gura 21) y la eciencia cuando se combinan todos los espacios menos uno (línea azul de la gura) con sus respectivas desviaciones estándar. El eje vertical representa los valores porcentuales de la eciencia y el eje horizontal representa el número de estados. En la gura 21 (a) se observa los resultados arrojados de la comparación entre TFMWS y FMWS; para la gura 21 (b) se muestra la comparación entre TFMWS y TFMS; para la gura 21 (c) se puede observar la comparación entre TFMWS y TFMW; en la gura 21 (d) la comparación se hace entre TFMWS y TFWS y por último, en la gura 21 (e) se compara TFMWS con TMWS. 100102 Número de Estados EFICIENCIA DESVIACIÓN ESTÁNDAR TFMWS FMWS TMWS TFMS TFWS TFMW TFMWS FMWS TMWS TFMS TFWS TFMW Tabla 14: Resultados prueba 2 101103 (a) (b) (c) (d) (e) Figura 21: Resultados prueba 2 102104 Como se observa, en la gura 21 (a) al tener el conjunto de todos los dominios de representación TFMWS y al quitar el dominio del tiempo T la eciencia representada por la línea de color azul tiende a permanecer por debajo de la eciencia cuando se trabaja con todos los dominios de representación manteniédonse este comportamiento a lo largo de los 8 estados. De esto se intuye que las características del dominio del tiempo pesan de manera que aumentan la eciencia en la clasicación. Se observa que la desviación estándar aumenta progresivamente junto con el número de estados para el conjunto de entrenamiento FMWS mientras que para TFMWS no arroja una tendencia clara debido a la naturaleza aleatoria del proceso. En la gura 21 (b) se observa que cuando se omite el espacio de representación de tiempofrecuencia Wavelet el comportamiento tiende a ser muy irregular, las curvas TFMWS y TFMS tienden a estar sobrepuestas la una con la otra, por lo que se puede inferir que las características de Wavelet no aumenta ni disminuye considerablemente la eciencia en la clasicación cuando se entrenan estas con las características de los demás espacios utilizados. Por su parte la desviación estándar es directamente proporcional al número de estados del clasicador. Al quitar el espacio de representacion de tiempo-frecuencia STFT tal como se observa en la gura 21 (c) se puede ver que la eciencia tiende a aumentar, se inere por lo tanto que las características de STFT disminuyen considerablemente la ecieciencia en la clasicación cuando se entrenan estas características con las características de los demás espacios. Adicionalmente la desviación estándar tiende a aumentar a medida que crece el número de estados del clasicador. Para la gura 21 (d) se ha omitido el espacio de representación de tiempo-frecuencia Melcepst, a partir de esta omision cabe anotar que la eciencia presenta una tendencia uctuante, las curvas TFMWS y TFWS tienden a estar sobrepuestas la una con la otra, por lo que se puede inferir que las características de la transformada de Melcepst no aumenta ni disminuye considerablemente la eciencia en la clasicación cuando se entrenan estas con las características de los demás espacios utilizados. El valor para la desviación es directamente proporcional al número de estados. Teniendo en cuenta los resultados representados en la gura 21 (e) al quitar el espacio de representación de la frecuencia, la eciencia representada por la línea de color azul tiende a permanecer por debajo de la eciencia cuando se entrena con todos los dominios de representación, este comportamiento se mantiene a lo largo de todos los estados, de lo que se induce que las características 103105 de frecuencia al ser combinada con las características de los demás espacios de representación, estas pesan de manera que aumentan la eciencia en la clasicación. Se observa que la desviación estándar aumenta progresivamente junto con el número de estados para el conjunto de entrenamiento TMWS mientras que para TFMWS no arroja una tendencia clara debido a que es un proceso aleatorio. En general los resultados de la tabla 14 que se obtienen en el momento de comparar el comportamiento al combinar todos los espacios de representación con el comportamiento cuando se omite un espacio de representación al entrenar con todas las características. Las caracteristicas de STFT pesan de forma negativa al momento de combinar este espacio con los demas espacios de representación ya que la eciencia presenta una disminución notable. Por otro lado, cuando las combinaciones se realizan con los dominios del tiempo o frecuencia las eciencias mejoran. Para las transformadas Wavelet y Melcepst las eciencias no presentan un cambio brusco permitiendo inferir que no realizan aportes signicativos a la clasicación. Por otro lado, la desviación estándar para el conjunto de entrenamiento TFMWS tiende a mantener un comportamiento ascendente a medida que se aumenta el número de estados con la peculiaridad que para los estados 2 y 4 las desviaciones son atípicamente altas, esto se debe a la naturaleza aleatoria del proceso. Para las demás combinaciones la desviación estándar aumenta a medida que se incrementa el número de estados. 21 Prueba 3: Evaluación de las mejores características por espacio de representación Los siguientes grácos de barras corresponden a la prueba 3 la cual tiene como nalidad mostrar el comportamiento del clasicador, cuando se entrena con cada característica individualmente en los diferentes espacios de representación a medida que aumenta el número de estados evaluando las mejores características. 104106 21.1 Frecuencia Los datos obtenidos de la prueba 3 para el dominio de la frecuencia se encuentran en la tabla 15, las las representan las características y las columnas los valores de la eciencia para cada uno de los 8 estados. A partir de los datos consignados en la tabla se construyó la gura 22 en la cual se muestra los resultados para el dominio de la frecuencia, donde cada barra indica el nivel de eciencia del clasicador cuando es entrenado con la característica individualmente. El eje vertical corresponde a los valores porcentuales de la eciencia en la clasicación mientras que el eje horizontal indica la característica utilizada para el entrenamiento. La gura 22 (a) corresponde al comportamiento para el estado 1; la gura 22 (b) para el estado 2 y asi sucesivamente hasta completar 8 estados. La línea horizontal de color rojo representa el umbral superior determinado por la característica de mas baja eciencia en la clasicación del conjunto de las más representativas mientras que la línea de color verde representa el umbral inferior determinado por la característica con menor eciencia en la clasicación en general. Características Eciencia Características Individuales Tabla 15: Resultados prueba 3 Frecuencia 105107 (a) (b) (c) (d) (e) (f) (g) (h) Figura 22: Resultados prueba 3 - Frecuencia 106108 En general se concluye que las mejores caracteristicas para entrenar individualmente en el dominio de la Frecuencia son las 1, 2, 7, 8, 9, 14, 15, 16 siendo 7, 8, 9 las de más representatividad en cuanto a eciencia se reere, mientras que la que menos representatividad tiene es la número 13, este comportamiento se mantiene a medida que se aumentan los estados, tal como se puede observar en la gura 22 y la tabla 15. Las mejores características no bajan de 81% de eciencia, esto es para la característica 8 entrenando con 8 estados mientras que la menos representativa en la caracteristica 3 que no supera el 15,6% entrenando con un estado Melcepst En la tabla 16, se puede observar los resultados de la prueba 3 para la transformada de Melcepst donde las las representan las características y las columnas los valores de la eciencia para cada uno de los 8 estados. A partir de esto se construye la gura 23 donde se muestra los resultados para Melcepst, cada barra indica el nivel de eciencia del clasicador cuando es entrenado con la característica individualmente. El eje vertical corresponde a los valores porcentuales de la eciencia en la clasicación mientras que el eje horizontal indica la característica utilizada para el entrenamiento. La gura 22 (a) corresponde al comportamiento para el estado 1; la gura 23 (b) para el estado 2 y asi sucesivamente hasta completar 8 estados. La línea horizontal de color rojo representa el umbral superior determinado por la característica de mas baja eciencia en la clasicación del conjunto de las más representativas mientras que la línea de color verde representa el umbral inferior determinado por la característica con menor eciencia en la clasicación en general. 107109 Características Eciencia Características Individuales Tabla 16: Resultados prueba 3 Melcepst 108110 (a) (b) (c) (d) (e) (f) (g) (h) Figura 23: Resultados prueba 3 - Melcepst 109111 A partir de la información brindada en la gura 23 se puede concluir que la mejor caracteristica para entrenar individualemte en el dominio Melcepst es la primera, seguido por la característica 2 y las que menos aportan son las características 5 y 10, como se observa en la tabla 16. Esta tendencia se mantiene a medida que aumentan los estados. Para este espacio la característica 1 entrenada con un estado no baja del 82,2%, mientras que el umbral inferior es de 44,7% representado por la característica 10 para 3 estados STFT En la tabla 17, se puede observar los resultados de la prueba 3 para STFT donde las las representan las características y las columnas los valores de la eciencia para cada uno de los 8 estados. En la gura 24 se muestra los resultados para STFT, cada barra indica el nivel de eciencia del clasicador cuando es entrenado con la característica individualmente. El eje vertical corresponde a los valores porcentuales de la eciencia en la clasicación mientras que el eje horizontal indica la característica utilizada para el entrenamiento. La gura 22 (a) corresponde al comportamiento para el estado 1; la gura 24 (b) para el estado 2 y asi sucesivamente hasta completar 8 estados. La línea horizontal de color rojo representa el umbral superior determinado por la característica de mas baja eciencia en la clasicación del conjunto de las más representativas mientras que la línea de color verde representa el umbral inferior determinado por la característica con menor eciencia en la clasicación en general. Características Eciencia Características Individuales Tabla 17: Resultados prueba 3 STFT 110112 (a) (b) (c) (d) (e) (f) (g) (h) Figura 24: Resultados prueba 3 -STFT. 111 Mostrar más
Análisis de vibraciones e interpretación de datos Jesús A. Royo Gloria Rabanaque Fernando Torres DIDYF Universidad de Zaragoza El análisis de vibraciones, la termografía, el análisis de lubricantes, entre Más detalles www.cdimca.com EL CENTRO DE DESARROLLO EN INGENIERÍA DE MANTENIMIENTO PRESENTA EL CURSO PROFESIONAL DE:
EL CENTRO DE DESARROLLO EN INGENIERÍA DE MANTENIMIENTO PRESENTA EL CURSO PROFESIONAL DE: www.cdimca.com Avalado por los Centros de Ingenieros de los Estados Lara, Mérida, Zulia y Colegio de Ingenieros Más detalles MANTENIMIENTO PREDICTIVO
MANTENIMIENTO PREDICTIVO CURSO sobre y sus distintas TÉCNICAS de APLICACIÓN Realizada por: Eñaut Labaien (Ingeniero Industrial) Gorka Carrasco (Ingeniero Técnico) Parque Tecnológico de San Sebastián Miramon Más detalles TEMA 8 VIBRACIONES EN MÁQUINAS. MANTENIMIENTO PREDICTIVO. Vibraciones en máquinas. Mantenimiento predictivo
Vibraciones en máquinas. Mantenimiento predictivo ELEMENTOS DE MÁQUINAS Y VIBRACIONES - 8.1 - ELEMENTOS DE MÁQUINAS Y VIBRACIONES - 8.2 - 8.1 Introducción El análisis y monitorizado de vibraciones son Más detalles Código: 101112016 Titulación: INGENIERO EN ORGANIZACIÓN INDUSTRIAL Curso: 5º
ASIGNATURA: INGENIERÍA DEL MANTENIMIENTO INDUSTRIAL Código: 10111016 Titulación: INGENIERO EN ORGANIZACIÓN INDUSTRIAL Curso: 5º Profesor(es) responsable(s): - ANICETO VALVERDE MARTÍNEZ - JOSÉ LUIS AGUIRRE Más detalles MANTENIMIENTO PREDICTIVO Y MONITOREO SEGUN CONDICION
MANTENIMIENTO PREDICTIVO Y MONITOREO SEGUN CONDICION Dr. lng. Pedro Saavedra G. Profesor Departamento de Ingeniería Mecánica Universidad de Concepción Ingeniero Consultor en Mantenimiento Predictivo. El Más detalles INGENIERÍA EN MANTENIMIENTO INDUSTRIAL
HOJA DE ASIGNATURA CON DESGLOSE DE UNIDADES TEMÁTICAS 1. Nombre de la asignatura Mantenimiento predictivo mecánico 2. Competencias Diseñar estrategias de mantenimiento mediante el análisis de factores Más detalles LA EMPRESA SERVICIOS. Reparación. Reparación de maquinaria eléctrica Sustitución de devanados de BT Mecanización de partes mecánicas
CATÁLOGO SERVICIOS 2014 LA EMPRESA Mantenimiento Industrial RAE, S.L.U., es una empresa dedicada a la reparación y mantenimiento de maquinaria eléctrica. La empresa está ubicada en Castellón de la Plana, Más detalles ALCANCES DE LA IMPLEMENTACION DE NUEVAS TECNICAS DE ANALISIS EN LOS PROGRAMAS DE MANTENIMIENTO PREDICTIVO PROACTIVO EN LA INDUSTRIA
ALCANCES DE LA IMPLEMENTACION DE NUEVAS TECNICAS DE ANALISIS EN LOS PROGRAMAS DE MANTENIMIENTO PREDICTIVO PROACTIVO EN LA INDUSTRIA Resumen Edgar A. Estupiñan P. Departamento de Ingeniería Mecánica Universidad Más detalles Índice. Capítulo 1. Introducción a las turbinas de gas 11
Índice Capítulo 1. Introducción a las turbinas de gas 11 1.1 Ciclo de las turbinas de gas en el modo de ciclo combinado o de cogeneración...12 1.2 Rendimiento de la turbina de gas...14 1.3 Consideraciones Más detalles ANALISIS DE VIBRACIONES APLICADO A LAS MAQUINAS ROTATORIAS DE BAJA VELOCIDAD
ANALISIS DE VIBRACIONES APLICADO A LAS MAQUINAS ROTATORIAS DE BAJA VELOCIDAD Edgar A. Estupiñan (eestupin@uta.cl) Profesor Dpto. de Ingeniería Mecánica Universidad de Tarapacá Casilla 6-D Arica (Chile) Más detalles OPERADORES MECANICOS
OPERADORES MECANICOS 0.- INTRODUCCION 1.- OPERADORES QUE ACUMULAN ENERGIA MECANICA 1.1.- Gomas 1.2.- Muelles 1.3.- Resortes 2.- OPERADORES QUE TRANSFORMAN Y TRANSMITEN LA ENERGIA MECANICA 2.1- Soportes Más detalles 1. CONCEPTO DE CONFIABILIDAD, HISTORIA Y CAMPOS DE APLICACIÓN
1. CONCEPTO DE CONFIABILIDAD, HISTORIA Y CAMPOS DE APLICACIÓN 1.1 CONCEPTO DE CONFIABILIDAD Desde el punto de vista de Ingeniería, se ha aceptado la definición de Confiabilidad de un elemento, como la Más detalles Prototipo de motor de combustión
2003 Emerson Process Management. Todos los derechos reservados. Vea este y otros cursos en línea en www.plantwebuniversity.com. Disponibilidad 201 Mejorando la Disponibilidad de la Planta Mediante la Instalación Más detalles Diagnóstico de Fallos por Monitoreo de Condición MONITOREO POR CONDICIÓN. Aplicación en Sistemas Expertos. Dr. Luis Amendola
Diagnóstico de Fallos por Monitoreo de Condición Dr. Luis Amendola España PUBLICADO 25/05/2003 La implementación de los sistemas de información para la gestión de activos es una nueva herramienta que contribuye Más detalles Implementando ESA como parte de su programa de mantenimiento predictivo para mejorar la Confiabilidad Eléctrica
Implementando ESA como parte de su programa de mantenimiento predictivo para mejorar la Confiabilidad Eléctrica Qué es ESA? Abstracto Este es el segundo de una serie de artículos que habla de cómo usar Más detalles Técnicas de Mantenimiento Predictivo
Entrenamiento y Capacitación Técnicas de Mantenimiento Predictivo Programa de Formación Continua www.aktiv.com.co PRESENTACIÓN El mantenimiento predictivo se compone una serie de técnicas que nos permiten Más detalles Tradicionalmente, se han distinguido 5 tipos de mantenimiento, que se diferencian entre sí por el carácter de las tareas que incluyen:
MANTENIMIENTO INDUSTRIAL 1. TIPOS DE MANTENIMIENTO Tradicionalmente, se han distinguido 5 tipos de mantenimiento, que se diferencian entre sí por el carácter de las tareas que incluyen: Mantenimiento Correctivo: Más detalles Las Técnicas de Monitoreo de Condición, como herramienta del Mantenimiento Proactivo
Las Técnicas de Monitoreo de Condición, como herramienta del Mantenimiento Proactivo Autora: Carolina Altmann El presente trabajo aborda, luego de una breve introducción a las Estrategias de Mantenimiento, Más detalles Recopilado por ForniSport
Los Amortiguadores y Resortes son parte fundamental de la Seguridad Activa del automóvil, y por ello la importancia que tienen en la Seguridad Vial, evitando accidentes provocados por falta de adherencia Más detalles Representación de señales de audio
Representación de señales de audio Emilia Gómez Gutiérrez Síntesi i Processament del So I Departament de Sonologia Escola Superior de Musica de Catalunya Curso 2009-2010 emilia.gomez@esmuc.cat 28 de septiembre Más detalles Academia de Ingeniería XVII Coloquio de Ingreso 26 de Marzo de 2015 Palacio de Minería de la Ciudad de México
Academia de Ingeniería XVII Coloquio de Ingreso 26 de Marzo de 2015 Palacio de Minería de la Ciudad de México Tecnología FPGA para el monitoreo y diagnóstico de fallas en maquinaria industrial Dr. Universidad Más detalles FAJAS. Proyectos de Ingeniería Mecánica Ing. José Carlos López Arenales
FAJAS Proyectos de Ingeniería Mecánica Ing. José Carlos López Arenales FAJAS Cintas o correas continuas que sirven para la transmisión de movimiento de dos o más ruedas que accionan o no un eje. FAJAS Más detalles Comprendiendo las estrategias de mantenimiento
2002 Emerson Process Management. Todos los derechos reservados. Vea este y otros cursos en línea en www.plantwebuniversity.com. Mantenimiento 101 Comprendiendo las estrategias de mantenimiento Generalidades Más detalles U.C: MANTENIMIENTO INDUSTRIAL. EL MANTENIMIENTO PREDICTIVO
U.C: MANTENIMIENTO INDUSTRIAL. EL MANTENIMIENTO PREDICTIVO 1. DEFINICION. TEMA Nº 3. MANTENIMIENTO PREDICTIVO. Servicios de seguimientos del desgaste de una pieza o más piezas, componente de equipos prioritarias Más detalles El monitoreo de una variable física requiere supervisión permanente de señales que
Capítulo 1 Marco Contextual 1.1. Formulación del problema 1.1.1. Definición del problema El monitoreo de una variable física requiere supervisión permanente de señales que varían con el tiempo. Tal información, Más detalles Cloud Monitoring aplicado a la supervisión y diagnóstico de activos críticos de plantas industriales. 1. Introducción
aplicado a la supervisión y diagnóstico de activos críticos de plantas industriales. David Faro Ruiz Director General Preditec/IRM Grupo Álava Ingenieros 1. Introducción El objetivo de la estrategia del Más detalles Montaje y Mantenimiento Mecánico
IES Rey Pelayo Programación 2014-15 Montaje y Mantenimiento Mecánico Instalación y Mantenimiento Electromecánico y Conducción de Líneas Formación Profesional de Grado Medio Índice Índice... 2 Introducción... Más detalles Ultrasonido: una técnica predictiva muy versátil en el mantenimiento
Ultrasonido: una técnica predictiva muy versátil en el mantenimiento Características Es una técnica bastante difundida por tres factores: fácil de usar, versatilidad y bajos costos de implementación. Pero, Más detalles Las tres áreas de SPM al servicio de la Gran Minería y sus equipos
Las tres áreas de SPM al servicio de la Gran Minería y sus equipos Maximiza la disponibilidad Aplicable a cualquier modelo de equipo Adaptable a cualquier estrategia de mantenimiento Compatible con contratos Más detalles Capítulo 1. Introducción. 1.1. Antecedentes
Capítulo 1. Introducción En este capítulo se presenta una descripción general del problema a investigar y el enfoque con el que se aborda. Se establece la necesidad de incorporar técnicas de análisis novedosas Más detalles SIMULADOR DE VIBRACIONES PARA LA REALIZACIÓN DE MANTENIMIENTO PREDICTIVO
SIMULADOR DE VIBRACIONES PARA LA REALIZACIÓN DE MANTENIMIENTO PREDICTIVO Jesús A. Royo Sánchez, Ana Brun de Buen, Maria Pilar Lambán Castillo, Daniel Pérez Soriano, José Javier Gutiérrez Aguilar. Profesores Más detalles CURSO ANALISIS DE VIBRACIONES DE MAQUINAS CATEGORÍA II, ISO 18436-2
LABORATORIO DE VIBRACIONES MECANICAS DEPARTAMENTO DE INGENIERIA MECANICA FACULTAD DE INGENIERIA http://www.dim.udec.cl/lvm/ CURSO ANALISIS DE VIBRACIONES DE MAQUINAS CATEGORÍA II, ISO 18436-2 Fechas: Junio Más detalles Comisión Nacional para el Uso Eficiente de la Energía
Comisión Nacional para el Uso Eficiente de la Energía Guía para ahorrar energía eléctrica en motores de corriente alterna, trifásicos, de inducción, tipo. Julio 2013 Contenido 1. Introducción 2. Objetivo Más detalles Unidad II: Actuadores
Unidad II: Actuadores Los actuadores tienen como misión generar el movimiento de los elementos del robot según las órdenes dadas por la unidad de control. Se clasifican en tres grandes grupos, según la Más detalles MANTENIMIENTO Y REPARACION DE MOTORES ELECTRICOS por Ing. Ciro Martínez T.
MANTENIMIENTO Y REPARACION DE MOTORES ELECTRICOS por Ing. Ciro Martínez T. ANTECEDENTES Los problemas de los motores eléctricos han sido enfocados mayormente desde el punto de vista eléctrico sin embargo Más detalles Polígono Indutrial O Rebullón s/n. 36416 - Mos - España - rodavigo@rodavigo.com
Productos de mantenimiento y lubricación Herramientas mecánicas... 1070 Llaves de gancho y de impacto... 1070 Llaves para el montaje de rodamientos de bolas a rótula y cubos axiales de ajuste de tuercas Más detalles CAPÍTULO 3 Programación en LabVIEW
CAPÍTULO 3 Programación en LabVIEW 3.1 Conexión física de los elementos Para capturar todas las señales provenientes de los sensores se utilizó una tarjeta de adquisición de datos de National Instruments, Más detalles Volante bimasa. temblores en la carrocería, con pérdidas de confort. La solución, en vehículos diésel, es el volante bimasa.
Volante bimasa LA EVOLUCIÓN TECNOLÓGICA EN LOS MOTORES HACE QUE CADA VEZ SEAN MÁS PEQUEÑOS EN CILINDRADA Y MASA, A LA PAR QUE SE REDUCEN SUS REVOLUCIONES, OFRECIENDO, NO OBSTANTE, MAYORES PRESTACIONES. Más detalles El presente documento tiene por objeto registrar y difundir internamente los casos, novedades y experiencias
1 N 20 Setiembre 2.005 El presente documento tiene por objeto registrar y difundir internamente los casos, novedades y experiencias propias dentro del Mantenimiento Predictivo /Proactivo y ponerlos a disposición Más detalles EFICIENCIA EN LOS SISTEMAS DE BOMBEO Y DE AIRE COMPRIMIDO
EFICIENCIA EN LOS SISTEMAS DE BOMBEO Y DE AIRE COMPRIMIDO 1. GENERALIDADES La sencillez en la operación, la disponibilidad, la facilidad y la seguridad en el manejo de las herramientas y elementos neumáticos Más detalles DIAGNÓSTICO DE FALLOS BASADO EN LA CONDICIÓN
DIAGNÓSTICO DE FALLOS BASADO EN LA CONDICIÓN Amendola, Luis., Ph.D Certificado Nivel B-IPMA Senior Project Manager PMM Institute for Learning Universidad Politécnica de Valencia, España Departamento de Más detalles Capítulo II. Motores
Capítulo II. Motores En la actualidad, son ampliamente utilizados en nuestra vida cotidiana, desde utensilios en el hogar como lo son licuadoras, batidoras, procesadores de alimentos, extractor de jugos, Más detalles MANTENIMIENTO PREDICTIVO.
MANTENIMIENTO PREDICTIVO. INTRODUCCIÓN. Sin dudas, el desarrollo de nuevas tecnologías ha marcado sensiblemente la actualidad industrial mundial. En los últimos años, la industria mecánica se ha visto Más detalles DIAGNÓSTICO POR VIBRACIONES: Motor de corriente continua:
DIAGNÓSTICO POR VIBRACIONES: Motor de corriente continua: Rotor de una máquina eléctrica de corriente continua: Colector de una máquina eléctrica de corriente continua: FRECUENCIAS CARACTERÍSTICAS DE UN Más detalles TECNICAS DE DIAGNOSTICO PARA EL ANALISIS DE VIBRACIONES DE RODAMIENTOS
TECNICAS DE DIAGNOSTICO PARA EL ANALISIS DE VIBRACIONES DE RODAMIENTOS Estupiñan, Edgar - Saavedra, Pedro Universidad de Concepción DEPARTAMENTO DE INGENIERÍA MECÁNICA (Casilla 160-C, Fono: 204327, Fax: Más detalles UNIVERSIDAD TECNOLÓGICA DE PEREIRA. Proyecto de desarrollo tecnológico con el fin de optar al título de Magister en Ingeniería Eléctrica.
UNIVERSIDAD TECNOLÓGICA DE PEREIRA Proyecto de desarrollo tecnológico con el fin de optar al título de Magister en Ingeniería Eléctrica. CLASIFICACIÓN DE FALLAS EN MOTORES ELÉCTRICOS UTILIZANDO SEÑALES Más detalles Contrato Certified Maintenance Partner. Introducción:
Contrato Certified Maintenance Partner. Introducción: Septiembre 2013: SKF Española y Dexis Ibérica firman el contrato de Certificación CMP (Certified Maintenance Partner). Según ello, Dexis Ibérica, es Más detalles CAUSAS DE FALLOS EN LOS COJINETES DE FRICCIÓN
CAUSAS DE FALLOS EN LOS COJINETES DE FRICCIÓN El objetivo de este manual es poder diagnosticar las causas de fallos en los motores por el aspecto que puedan presentar los cojinetes de fricción del mismo. Más detalles Fundamentos para la Representación y Análisis de Señales Mediante Series de Fourier
Fundamentos para la Representación y Análisis de Señales Mediante Series de Fourier Andrés Felipe López Lopera* Resumen. Existe una gran similitud entre vectores y las señales. Propiedades tales como la Más detalles El presente documento tiene por objeto registrar y difundir internamente los casos, novedades y experiencias
1 N 19 Agosto 2.005 El presente documento tiene por objeto registrar y difundir internamente los casos, novedades y experiencias propias dentro del Mantenimiento Predictivo /Proactivo y ponerlos a disposición Más detalles TEMA 9: Desarrollo de la metodología de Taguchi
TEMA 9: Desarrollo de la metodología de Taguchi 1 La filosofía de la calidad de Taguchi 2 Control de calidad Off Line y On Line Calidad Off Line Calidad On Line 3 Función de pérdida 4 Razones señal-ruido Más detalles GUIA DE LUBRICACION. Lubricación Basada en Condición con Ultrasonido. Productividad Maximizada a través del ULTRASONIDO
GUIA DE LUBRICACION Lubricación Basada en Condición con Ultrasonido Lubricación Basada en Condición con Ultrasonido Tradicionalmente, la programación de la lubricación ha sido basada en el tiempo. Los Más detalles ACELERÓMETROS Y SENSORES DE VIBRACIÓN
ACELERÓMETROS Y SENSORES DE VIBRACIÓN Introducción Vibración es el movimiento oscilatorio de un cuerpo respecto a su posición de equilibrio y, en consecuencia puede caracterizarse mediante tres magnitudes Más detalles Informe técnico sobre el RailScanner
Informe técnico sobre el RailScanner 1.Introducción El RailScanner es un sistema electrónico de monitoreo que verifica la condición de la vía férrea bajo una carga operativa real. Lo que hace es medir Más detalles 1º Tema.- Conceptos y definiciones en cinemática.
Universidad de Huelva ESCUELA POLITECNICA SUPERIOR Departamento de Ingeniería Minera, Mecánica y Energética Asignatura: Ingeniería de Máquinas [570004027] 5º curso de Ingenieros Industriales 1º Tema.- Más detalles Detectando Fallas Prematuras en Rodamientos
Art. 001 WAL Engineering Detectando Fallas Prematuras en Rodamientos Los Fabricantes de rodamientos durante muchos años han alertado acerca de la relación que tiene una alta temperatura con la vida útil Más detalles Desbalanceo Mecánico
LGM-09-03 2009-marzo Desbalanceo Mecánico La Guía MetAs El desbalanceo mecánico en elementos rotatorios (rotor) se ha convertido en un problema importante en el desarrollo de maquinaria moderna, especialmente Más detalles Salud de Activos Reflejo de la Estrategia de Mantenimiento
Salud de Activos Reflejo de la Estrategia de Mantenimiento Mucho se ha dicho y escrito acerca de como medir la efectividad de una estrategia de mantenimiento, sin embargo, al momento solo porciones de Más detalles QUÉ ES EL MANTENIMIENTO?
Documento elaborado por I.M. Luis Alberto Cuartas Pérez 2008 1 QUÉ ES EL MANTENIMIENTO? Es el trabajo emprendido para cuidar y restaurar hasta un nivel económico, todos y cada uno de los medios de producción Más detalles NUEVAS HERRAMIENTAS DE MANTENIMIENTO FERROVIARIO PREDICTIVO APROVECHANDO LAS TIC
NUEVAS HERRAMIENTAS DE MANTENIMIENTO FERROVIARIO PREDICTIVO APROVECHANDO LAS TIC Trilla Castelló, Alexandre Ingeniero de procesamiento de datos TLS España y Portugal, ALSTOM TRANSPORTE S.A. Coca Sola, Más detalles CAPITULO III EL MANTENIMIENTO. 3.1 Conceptos Básicos de Mantenimiento
39 CAPITULO III EL MANTENIMIENTO 3.1 Conceptos Básicos de Mantenimiento Antes que realizar cualquier tipo de clasificación o descripción de los tipos de mantenimiento, es muy importante saber realmente Más detalles MÉTODO SISTEMÁTICO DE INTERPRETACIÓN DE ANÁLISIS DE LUBRICANTE
MÉTODO SISTEMÁTICO DE INTERPRETACIÓN DE ANÁLISIS DE LUBRICANTE GERARDO TRUJILLO C CMRP, MLA, MLT. GTRUJILLO@NORIA.MX COPYRIGHT NORIA CORPORATION COACHING CONSULTORÍA ENTRENAMIENTO QUIÉNES SOMOS Somos expertos Más detalles MANTENIMIENTO PREVENTIVO DE MEDIOS DE TRANSPORTE. Programa Técnico Profesional en Transporte terrestre y de carga
MANTENIMIENTO PREVENTIVO DE MEDIOS DE TRANSPORTE Programa Técnico Profesional en Transporte terrestre y de carga Descripción del curso Si bien el mantenimiento hace parte del mapa de procesos de las compañías, Más detalles EL MANTENIMIENTO EN EL SECTOR INDUSTRIAL
I. INTRODUCCION EL MANTENIMIENTO EN EL SECTOR INDUSTRIAL Las operaciones de mantenimiento tienen lugar frente a la constante amenaza que implica la ocurrencia de una falla o error en un sistema, maquinaria, Más detalles TEMA 9 Mantenimiento en el JIT: mantenimiento productivo total (TPM)
TEMA 9 Mantenimiento en el JIT: mantenimiento productivo total (TPM) 9.1. Introducción.. 9.3. Principios básicos Anexos 1 9.1. Introducción Para resistir la competitividad, las empresas deben continuamente Más detalles Prevención de fallas en máquinas rotatorias con el monitoreo de DP en línea
Prevención de fallas en máquinas rotatorias con el monitoreo de DP en línea El monitoreo de descargas parciales (DP) en línea es una herramienta de diagnóstico esencial para un mantenimiento basado en Más detalles Seguridad y mantenimiento predictivo
Tema de tapa Seguridad y mantenimiento predictivo Beneficios de una aplicación en conjunto Por Ing. Mario Lanza Pride Argentina Introducción Las técnicas de mantenimiento han evolucionado conforme la tecnología Más detalles TUTORIAL DE VIBRACIONES PARA MANTENIMIENTO MECÁNICO REALIZADO POR A-MAQ S.A. ANÁLISIS DE MAQUINARIA ENERO 2005. Visítenos en www.a-maq.com.
TUTORIAL DE VIBRACIONES PARA MANTENIMIENTO MECÁNICO REALIZADO POR A-MAQ S.A. ANÁLISIS DE MAQUINARIA ESPECTRO Amplitud ENERO 2005 0 20 40 60 80 f [Hz] TABLA DE CONTENIDO MEMORIAS DEL SEMINARIO DE VIBRACIONES Más detalles Referencia a la Norma ISO 9001:2008 7.1, 7.2.1, 7.5.1, 7.6, 8.1, 8.2.4
Página 1 de 13 Instituto Tecnológico Superior de El Grullo Subdirección Académica Departamento de Electromecánica Instrumentación Didáctica para la Formación y el Desarrollo de Competencias Del periodo Más detalles ADQUISICIÓN DE VIBRACIONES MECÁNICAS DE UN MOTOR EN FUNCIONAMIENTO USANDO LABVIEW
ADQUISICIÓN DE VIBRACIONES MECÁNICAS DE UN MOTOR EN FUNCIONAMIENTO USANDO LABVIEW Douglas Acosta (1), Johannex Molina (2), Holger Cevallos (3) (1) (2) Facultad de Ingeniería en Electricidad y Computación Más detalles Conversor Analógico Digital (CAD)
Conversor Analógico Digital (CAD) La salida de los sensores, que permiten al equipo electrónico interaccionar con el entorno, es normalmente una señal analógica, continua en el tiempo. En consecuencia, Más detalles INTRODUCCIÓN A LOS BALASTROS ELECTRÓNICOS
1 INTRODUCCIÓN A LOS BALASTROS ELECTRÓNICOS 1.1 INTRODUCCIÓN En la actualidad existe la necesidad de controlar la potencia eléctrica de los sistemas de iluminación, tracción y motores eléctricos debido Más detalles TEORÍA DE MECANISMOS 1.- INTRODUCCIÓN TMM
TEORÍA DE MECANISMOS 1.- INTRODUCCIÓN TMM Departamento de Ingeniería Mecánica 1 Teoría de Mecanismos NIVELES 0) CONOCIMIENTOS PREVIOS Y FUNDAMENTOS TMM ELEMENTO/MIEMBRO EXPRESIÓN GRÁFICA MECÁNICA RESISTENCIAS Más detalles Operaciones y Mantenimiento 101 Estrategias de Mantenimiento y Prácticas de Trabajo para Reducir los Costos
2003 Emerson Process Management. Todos los derechos reservados. Vea este y otros cursos en línea en www.plantwebuniversity.com. Operaciones y Mantenimiento 101 Estrategias de Mantenimiento y Prácticas Más detalles Disponibilidad 301 Mejorando la Disponibilidad Mediante la Detección y el Monitoreo
2003 Emerson Process Management. Todos los derechos reservados. Vea este y otros cursos en línea en www.plantwebuniversity.com. Disponibilidad 301 Mejorando la Disponibilidad Mediante la Detección y el Más detalles 7 HERRAMIENTAS PARA EL CONTROL DE LA CALIDAD
Agencia de Cooperación Internacional del Japón Universidad de Santiago de Chile Facultad de Ingeniería Departamento de Ingeniería Industrial 7 HERRAMIENTAS PARA EL CONTROL DE LA CALIDAD Elaboración: Kiyohiro Más detalles Operaciones y Mantenimiento 201 Uso del monitoreo en líneas y de los diagnósticos predictivos para reducir los costos de mantenimiento
2003 Emerson Process Management. Todos los derechos reservados. Vea este y otros cursos en línea en www.plantwebuniversity.com. Operaciones y Mantenimiento 201 Uso del monitoreo en líneas y de los diagnósticos Más detalles RCM - Mantenimiento Centrado en Confiabilidad
RCM - Mantenimiento Centrado en Confiabilidad www.rcm-confiabilidad.com.ar 30 de Octubre, 2005 1 Qué es el mantenimiento centrado en Confiabilidad (RCM)? El mantenimiento centrado en Confiabilidad (MCC), Más detalles TEMA 6.- LA PUESTA EN MARCHA DE PROYECTOS. LOS ASPECTOS A TENER EN CUENTA
TEMA 6.- LA PUESTA EN MARCHA DE PROYECTOS. LOS ASPECTOS A TENER EN CUENTA El Programa para el Fomento de la Intraemprendeduría en Ciclos Formativos de Formación Profesional es un proyecto financiado por Más detalles 3.1 Mantenimiento preventivo
3. Mantenimiento Las estaciones de medición, ya sean fijas o móviles, representan la infraestructura remota de los SMCA, por lo que requieren de buenas y oportunas prácticas de mantenimiento preventivo, Más detalles CAPITULO I. Mantenimiento: Generalidades
3 CAPITULO I Mantenimiento: Generalidades 1.1 DEFINICIÓN DE MANTENIMIENTO El mantenimiento es un servicio que agrupa una serie de actividades cuya ejecución permite alcanzar un mayor grado de confiabilidad Más detalles Metodología de monitoreo, detección y diagnóstico de fallos en motores asíncronos de inducción
Metodología de monitoreo, detección y diagnóstico de fallos en motores asíncronos de inducción Castelli, Marcelo Universidad de Montevideo. IEEE Member Fossati, Juan Pablo Universidad de Montevideo Andrade, Más detalles Características Morfológicas. Principales características de los Robots.
Características Morfológicas Principales características de los Robots. Se describen las características más relevantes propias de los robots y se proporcionan valores concretos de las mismas, para determinados Más detalles OSCILACIONES ARMÓNICAS
APLICACIONES El campo de aplicación de las fibras ópticas es muy amplio y aumenta día a día. Algunas de las aplicaciones más importantes son: - Telecomunicaciones: En este apartado cabe incluir la red Más detalles MECÁNICA AUTOMOTRIZ. mezcla. Válvula de escape cerrada. Válvula de admisión cerrada.
MECÁNICA AUTOMOTRIZ Principio de funcionamiento La bujía inflama la mezcla. Válvula de escape cerrada. Válvula de admisión cerrada. El pistón es impulsado hacia abajo ante la expansión producida por la Más detalles EVALUACIÓN DEL IMPACTO METROLÓGICO DE TAREAS DE MANTENIMIENTO SOBRE EL MEDIDOR DE GAS TIPO TURBINA
EVALACIÓN DEL IMPACTO METROLÓGICO DE TAREAS DE MANTENIMIENTO SOBRE EL MEDIDOR DE GAS TIPO TRBINA Sergio Lupo, Juan Forastieri, Ezequiel Filipovic INTI - Instituto Nacional de Tecnología Industrial Colectora Más detalles PROGRAMA ANALÍTICO ORIENTACIÓN: SISTEMAS ELECTRÓNICOS INDUSTRIALES SISTEMAS ELECTRICOS DE POTENCIA
PROGRAMA ANALÍTICO DEPARTAMENTO: ELECTRICIDAD Y ELECTRÓNICA CARRERA: INGENIERÍA ELECTRICISTA ASIGNATURA: TÉCNICAS PARA EL MANTENIMIENTO Y DIAGNÓSTICO DE MÁQUINAS ROTATIVAS CÓDIGO: 0464 AÑO ACADÉMICO: 2014 Más detalles RESUMEN CUADRO DE MANDO
1. Objetivo Los objetivos que pueden alcanzarse, son: RESUMEN CUADRO DE MANDO Disponer eficientemente de la información indispensable y significativa, de modo sintético, conectada con los objetivos. Facilitar Más detalles SERVICIO NACIONAL DE APRENDIZAJE SENA SENA REGIONAL CALDAS CENTRO DE AUTOMATIZACIÓN INDUSTRIAL MANUAL DE PROCEDIMIENTOS PARA MANTENIMIENTO DE HARDWARE
SERVICIO NACIONAL DE APRENDIZAJE SENA SENA REGIONAL CALDAS CENTRO DE AUTOMATIZACIÓN INDUSTRIAL MANUAL DE PROCEDIMIENTOS PARA MANTENIMIENTO DE HARDWARE INTRODUCCIÓN El Servicio Nacional de Aprendizaje SENA Más detalles INSTRUMENTACIÓN GEOTÉCNICA Y ESTRUCTURAL, S.L.;
INSTRUMENTACIÓN GEOTÉCNICA Y ESTRUCTURAL, S.L.; ofrece a sus clientes un Sistema Integral de Gestión de Presas y Embalses, englobando las siguientes áreas o campos de actuación: Mantenimiento y Conservación Más detalles SESIÓN 8 TIPOS DE SISTEMAS DE MANUFACTURA
SESIÓN 8 TIPOS DE SISTEMAS DE MANUFACTURA CONTENIDO Términos en manufactura Clasificación de sistemas de manufactura Beneficios y ejemplos de sistemas de manufactura Los componentes de un sistema de manufactura Más detalles SISTEMA DE MONITOREO Y DIAGNÓSTICO AUTOMÁTICO ONLINE APLICADO AL REACTOR NUCLEAR RA6
SISTEMA DE MONITOREO Y DIAGNÓSTICO AUTOMÁTICO ONLINE APLICADO AL REACTOR NUCLEAR RA6 García Peyrano O A 1 (garciapeyrano@yahoo.com.ar); Marticorena M 1 ; Koch R G 1 ; Martínez J S 1 ; Berruti G E 2 ; Nuñez Más detalles Movimiento Armónico Simple
Movimiento Armónico Simple Introducción al Movimiento Armónico Simple En esta página se pretende que el alumno observe la representación del Movimiento Armónico Simple (en lo que sigue M.A.S.), identificando Más detalles ACERCA DEL DIAGNÓSTICO VIBROELÉCTRICO DE LOS MOTORES DE INDUCCIÓN: Experiencias del CEIM CUJAE
ACERCA DEL DIAGNÓSTICO VIBROELÉCTRICO DE LOS MOTORES DE INDUCCIÓN: Experiencias del CEIM CUJAE Dr. Evelio Palomino Marín CEIM - Centro de Estudios Innovación y Mantenimiento. Calle 114 No. 11901. Edificio Más detalles CAPITULO 1. Motores de Inducción.
CAPITULO 1. Motores de Inducción. 1.1 Introducción. Los motores asíncronos o de inducción, son prácticamente motores trifásicos. Están basados en el accionamiento de una masa metálica por la acción de Más detalles LA MEDICION Y ANALISIS DE LAS VIBRACIONES COMO TECNICA DE INSPECCION DE EQUIPOS Y COMPONENTES, APLICACIONES, NORMATIVAS Y CERTIFICACION
LA MEDICION Y ANALISIS DE LAS VIBRACIONES COMO TECNICA DE INSPECCION DE EQUIPOS Y COMPONENTES, APLICACIONES, NORMATIVAS Y CERTIFICACION Pedro Nelson Saavedra Facultad de Ingeniería - Universidad de Concepción, Más detalles Equipo de Diagnostico Mantenimiento centrado en la condición
Equipo de Diagnostico Mantenimiento centrado en la condición Por Ing. Carlos E. Flores R. ( 1 ) RESUMEN El mantenimiento centrado en la condición o mantenimiento predictivo, requiere de equipo especial Más detalles 2017 © DocPlayer.es Política de privacidad | Condiciones del servicio | Feedback

References: resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución