ES2336842T3 - Procedimiento para la deteccion e identificacion automatica de daños en rodamientos. - Google Patents
Procedimiento para la deteccion e identificacion automatica de daños en rodamientos. Download PDFInfo
- Publication number
- ES2336842T3 ES2336842T3 ES07460032T ES07460032T ES2336842T3 ES 2336842 T3 ES2336842 T3 ES 2336842T3 ES 07460032 T ES07460032 T ES 07460032T ES 07460032 T ES07460032 T ES 07460032T ES 2336842 T3 ES2336842 T3 ES 2336842T3
- Authority
- ES
- Spain
- Prior art keywords
- test
- examined
- interval
- indicates
- signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
- G01M13/04—Bearings
- G01M13/045—Acoustic or vibration analysis
Landscapes
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- General Physics & Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Investigating Or Analyzing Materials By The Use Of Magnetic Means (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Rolling Contact Bearings (AREA)
- Feedback Control In General (AREA)
Abstract
Procedimiento para la detección y la identificación automática de daños a rodamientos, en el que las vibraciones de un cojinete en funcionamiento, medidas como una señal de onda temporal analógica se convierte en datos digitales y se filtran previamente de manera conocida, caracterizado por el hecho de que se realiza una detección de impulsos de choque con el uso de un vector binario (p) que se determina en las siguientes acciones: - en la onda temporal filtrada previamente de la señal procesada se seleccionan dos intervalos iniciales adyacentes: uno es un intervalo de referencia y el otro es un intervalo de prueba y la señal presente en los intervalos de prueba y de referencia se someten a múltiples filtrados tras lo cual después de cada filtrado, se calcula una relación de probabilidad de impulsos de choque definida como una relación entre los valores medios cuadráticos de la señal filtrada presentes en el intervalo de prueba y los valores medios cuadráticos de la señal filtrada presentes en el intervalo de referencia, - las relaciones de probabilidad calculadas se comparan con un valor límite preestablecido y se estima la probabilidad de aparición de impulsos de choque en el intervalo de prueba inicial, - la información sobre la aparición de impulsos de choque en el intervalo de prueba examinado se registra como elementos iniciales del vector binario (p), y después - todos los elementos restantes del vector binario (p) se determinan para todos los intervalos de prueba y de referencia sucesivos que son seleccionados de manera reiterada con un desplazamiento en una dirección específica y aquellos elementos del vector binario (p) adoptan valores 1 en la zona de aparición del impulso de choque y 0 en la zona en la que no se producen impulsos de choque.
Description
Procedimiento para la detección e identificación
automática de daños en rodamientos.
El objeto de la invención es un procedimiento
para la detección e identificación automática de daños a rodamientos
aplicable al diagnóstico de rodamientos y la detección de daños en
una fase muy temprana de su aparición.
En los procedimientos actuales de diagnóstico de
equipos técnicos, la detección de fallos existentes y daños que se
producen en una fase temprana de uso del equipo o de sus elementos
individuales y su localización exacta es de una importancia
primordial, y contra más pronto pueda realizarse mayor es la
importancia del procedimiento. En el diagnóstico del estado del
equipo y especialmente la detección y la identificación automática
de defectos en máquinas eléctricas giratorias y sus elementos se
utilizan procedimientos basados en la medición de señales
eléctricas o mecánicas, su tratamiento así como el análisis y la
interpretación de los resultados de las mediciones.
Los procedimientos aplicados actualmente para
controlar cojinetes durante el funcionamiento consisten en comparar
los resultados medidos y procesados digitalmente con datos obtenidos
para cojinetes que funcionan en condiciones apropiadas y que no se
encuentran en peligro de fallo. En la mayoría de los casos, dicha
comparación la realiza un experto entrenado para este fin.
El procedimiento utilizado más frecuentemente
para la determinación temprana de defectos en cojinetes consiste en
el examen de impulsos de choque combinado con el examen de las
vibraciones en los elementos analizados del cojinete. Después,
dichas formas de los impulsos son seleccionadas a partir de los
datos medidos para un análisis adicional que sugiere la posibilidad
de la presencia de un defecto. Los procedimientos informatizados de
análisis de datos permiten detectar fallos de manera temprana y
determinar el lugar de su aparición, sus tipos y el nivel de
daño.
De la descripción de la patente americana US 5
477 730 se conoce un procedimiento y un dispositivo para comprobar
el estado de un cojinete. El procedimiento consiste en detectar las
señales primarias de las vibraciones emitidas por el cojinete, la
desmodulación de estas señales a señales de salida después de la
desmodulación, la conversión de las señales de salida de señales
analógicas a señales de salida digital, el filtrado de ruido fuera
de estas señales por medio de tecnología de filtrado digital
seleccionada de numerosas de tecnologías conocidas, y emisión de
las señales digitales a un dispositivo de salida. En otra variante
de la invención, en lugar de la desmodulación de señales primarias,
se utiliza filtrado de paso de banda en un intervalo de 2000 a 15
000 hertzios, y las señales después de este filtrado se someten a
transformaciones tal como se ha descrito anteriormente. Todavía en
otra variante de la invención las señales primarias de vibración
emitidas por el cojinete se convierten de analógicas a digitales,
con lo que el ruido es filtrado de las señales utilizando
tecnología de filtrado digital seleccionada de numerosas de
tecnologías conocidas, que incluye la acción de separación de
señales de ruido de las señales de la salida digital, y después las
señales digitales se envían al dispositivo de salida.
De la solicitud de patente WO 00/04361 se conoce
un procedimiento para evaluar el estado de los elementos de un
cojinete. Este procedimiento emplea las mediciones de las
vibraciones del elemento del cojinete generando una señal de
aceleración analógica, convirtiendo esta señal en datos digitales,
muestreando la señal digital para localizar sus componentes que
vienen del defecto del cojinete. El procedimiento se caracteriza por
el hecho de que la influencia de cualquier componente de señal
indeseado determinista es eliminado de la señal de aceleración
digital a través del filtrado para así obtener una señal que
consiste esencialmente en ruido blanco proporcional a la suavidad
global de la microestructura de la superficie del cojinete y el
impulso de choque. Se seleccionan los siguientes parámetros
característicos de la señal que son sensibles a los defectos del
cojinete: impulsos agudos de vibración o impulsos de choque,
amplitud pico a pico, valor medio proporcional a la amplitud media
sin impulsos de choque, contenido máximo de la señal. La detección
de impulsos de choque de banda estrecha a intervalos regulares se
realiza utilizando un procedimiento conocido, denominado
autocorrelación de la estimación de señal de envolvente.
La esencia del procedimiento para la detección y
la identificación automática de daños a rodamientos, en el que las
vibraciones de un cojinete en funcionamiento, medidas como una señal
analógica de onda temporal se convierten a datos digitales y se
filtran previamente de manera conocida, es que la detección de los
impulsos de choque se realiza utilizando un vector binario p que se
determina en las siguientes acciones:
- -
- en la onda temporal filtrada previamente de la señal procesada se seleccionan dos intervalos iniciales adyacentes: uno es un intervalo de referencia y el otro es un intervalo de prueba y la señal presente en los intervalos de prueba y de referencia se somete a filtrado múltiple con lo cual tras cada filtrado se calcula una relación de probabilidad de impulso de choque definida como la relación entre los valores medios cuadráticos de la señal filtrada presente en el intervalo de prueba y los valores medios cuadráticos de la señal filtrada presente en el intervalo de referencia,
- -
- las relaciones de probabilidad se comparan con un valor límite preestablecido y se estima la probabilidad de aparición de impulsos de choque en el intervalo de prueba inicial,
- -
- se registra la información sobre la aparición de impulsos de choque en el intervalo de prueba examinado como elementos iniciales del vector binario p, y después
- -
- se determinan todos los elementos restantes del vector binario p para todos los intervalos de prueba y de referencia sucesivos que son seleccionados de manera reiterada con un desplazamiento en una dirección específica y esos elementos del vector binario p adoptan valores 1 en la zona de aparición de impulsos de choque y 0 en la zona donde no se produce impulso de choque.
Preferiblemente, el valor medio cuadrático de la
señal en el intervalo de prueba examinado sometido a un filtrado
L^{n} \tau_{t} [L] se determina en base siguiente
relación:
mientras que el valor medio
cuadrático de la señal en el intervalo de referencia examinado
sometido a un filtrado L^{n} \tau_{r} [L] se determina en
base siguiente
relación:
donde:
L es un número natural del intervalo 1...
N_{s} e indica el número de filtros,
K_{t}, indica la longitud del intervalo de
prueba examinado,
K_{r}, indica la longitud del intervalo de
referencia,
k_{t}, k_{r} indica números naturales
sucesivos 1 a k_{t} o k_{r},
z^{L}_{t} indica un vector que representa el
intervalo de prueba examinado de la señal sometida a un filtrado
L^{n},
z^{L}_{r}, indica un vector que representa
el intervalo de referencia examinado de la señal sometida a un
filtrado L^{n},
Preferiblemente, el valor medio cuadrático de
una señal en el intervalo de referencia examinado \tau_{r} [L]
se determina utilizando un filtrado L^{n}, que es un filtrado
adaptativo, con el uso de descomposición en valores singulares
(SVD) de una matriz de autocovarianza C_{r}, de acuerdo con la
siguiente relación:
donde la matriz \Sigma_{r} es
una matriz diagonal cuyos elementos diagonales son distintos de cero
y no son negativos e indica valores medio cuadrático \tau_{r}
[L]:
y las matrices V y U son matrices
ortogonales y al mismo tiempo la matriz de autocovarianza C_{r}
presenta la siguiente
forma:
para la
cual:
donde:
x_{r} [k_{r}] son los elementos del
intervalo de referencia,
k_{r} indica números naturales sucesivos de 1
a Kr,
i_{r} indica números naturales sucesivos de 1
a L,
Kr indica la longitud del intervalo de
referencia,
L indica el número de filtros adaptativos.
Preferiblemente, el valor medio cuadrático de
una señal en el intervalo de prueba examinado \tau_{t} [L] se
determina utilizando un filtrado L^{n}, que es un filtrado
adaptativo, en base a la siguiente relación:
donde z_{t}^{L} son los
elementos de la matriz Z_{t} presentados de la siguiente
forma:
y se calculan de acuerdo con la
siguiente
relación:
donde la matriz X_{t} se define
según:
y la matriz U es una matriz
ortogonal cuyos elementos son los coeficientes del filtro
adaptativo, mientras
que:
Kt indica la longitud del intervalo de prueba
examinado;
k_{t}, [k_{r}] son los elementos del
intervalo de prueba;
k_{t} indica los números naturales sucesivos
de 1 al K_{t},
N_{t} =
K_{t} -
L,
L indica el número de filtros adaptativos.
Preferiblemente los múltiples filtrados de
señales en los intervalos de prueba y de referencia se realizan
mediante el uso de filtros de paso de banda.
Preferiblemente, la longitud del intervalo de
prueba es preferiblemente un 15% de la longitud del intervalo de
referencia.
Preferiblemente, la identificación del defecto
en el cojinete o cojinetes se realiza automáticamente para el
vector binario p determinado para todos los intervalos de prueba y
de referencia examinados de la onda temporal X de la señal
procesada y el conjunto de frecuencias de defectos conocidos del
cojinete o cojinetes examinados F_{f} = [f_{f1},
f_{f1}, ..., f_{fK}], por medio de las siguientes acciones:
f_{f1}, ..., f_{fK}], por medio de las siguientes acciones:
- -
- para cada frecuencia de defectos conocidos sucesivos f_{fK} del conjunto F_{t} del cojinete examinado, se determina la frecuencia de prueba v_{fK} \in | (1-P_{g}) f_{fK}- (1+P_{g}) f_{fK} |) que pertenece a los alrededores P_{g} de la frecuencia de defectos conocidos f_{fK}, donde P_{g} es un parámetro de configuración que indica la gama de frecuencias de búsqueda asociadas a la frecuencia de defectos conocidos y después, utilizando principios matemáticos de probabilidad máxima conocidos, para cada frecuencia de prueba v_{fK} se determina la función de probabilidad máxima J en el dominio de tiempo utilizando el vector binario p,
- -
- después, se comprueba para qué frecuencia de prueba v_{fK} la función de probabilidad máxima J alcanza el valor máximo J_{A},
- -
- se calcula la relación F_{A} del valor máximo J_{A} respecto al nivel de ruido n_{A} y el resultado se compara con el parámetro límite preestablecido S_{A}:
- -
- si la relación F_{A} es mayor que parámetro límite preestablecido S_{A}, se genera automáticamente información sobre el defecto del cojinete examinado con frecuencia específica que pertenece al conjunto F_{t} y corresponde a la frecuencia de prueba v_{fK},
- -
- si la relación F_{A} es inferior o igual al parámetro límite preestablecido S_{A}, se supone que ninguno de los impulsos de choque presentes en el vector binario p corresponde a ninguna frecuencia de defectos conocidos f_{fk} de cojinetes examinados y se genera automáticamente información de que no se ha producido ningún defecto específico en el cojinete o cojinetes examinados, en el que el nivel de ruido (n_{A}) se calcula mediante el uso de valores de función máxima (J).
\vskip1.000000\baselineskip
Preferiblemente, la frecuencia de prueba
v_{fK}, para la cual la función de probabilidad máxima calculada
para el vector binario p en dominio de tiempo alcanza el valor
máximo J_{A}, se determina preferiblemente maximizando la suma de
cuadrados de amplitudes armónicas A_{k} de acuerdo con la
siguiente relación:
donde:
N_{h} es un número preestablecido del armónico
examinado.
A_{k} es una amplitud armónica determinada en
el dominio de tiempo utilizando el vector binario p, presentado en
forma de modelo sinusoidal.
\vskip1.000000\baselineskip
Preferiblemente los valores de las amplitudes de
armónicos A_{k} se determinan en el dominio de tiempo utilizando
el vector binario p, presentado en forma de modelo sinusoidal. El
vector binario p se presenta como una suma del modelo sinusoidal y
el vector de ruido:
el cual también puede presentarse
en forma de la siguiente
matriz:
p=
B\theta+e
donde:
\phi_{k} es una fase,
\varpi_{k} =2\pivk/f_{s}, k= 1,
2... N_{h} es una velocidad angular correspondiente a una
frecuencia de prueba (v_{fK}), y k= 1, 2... N_{h},
f_{s} es una frecuencia de muestreo
para tratamiento analógico a digital de la señal analógica
medida,
e = [\varepsilon[n]]^{T} indica el
vector de ruido,
n indica un número natural de 1 a
N-1,
N indica el número de elementos en el vector
binario (p),
B es la matriz de parámetros constantes
conocidos, presentada de la siguiente manera:
y \theta es el vector de
parámetros desconocidos presentados de la siguiente
manera:
\theta =
[A_{1} cos(\phi_{1}), A_{1} sin(\phi_{1}),
A_{2} cos(\phi_{2}), A_{2}
sin(\phi_{2}),...,A_{Nh} cos(\phi_{Nh}),
A_{Nh}
sin(_{Nh})],
que es el estimador de mínimos
cuadrados para la
relación:
\theta
(B^{T}B)^{-1}B^{T}p
donde p es un vector binario de la
siguiente
manera:
p= [p
[0], p [1],..., p
[N-1]]^{T}
Preferiblemente, el nivel de ruido n_{A} se
calcula tomando un valor medio de valores de función máxima J
encontrados en las proximidades de la frecuencia del defecto
examinado f_{fK} del conjunto F_{f}.
Alternativamente, el nivel de ruido n_{A} se
calcula tomando un valor medio de valores de función máxima J
encontrados en las proximidades de la frecuencia del defecto
examinado f_{fK} del conjunto F_{f}.
En comparación con los procedimientos conocidos
de detección e identificación automática de defectos en rodamientos,
la solución de la invención permite la detección de defectos en una
fase muy temprana de su aparición. La señal medida es procesada de
modo que los componentes de la señal asociada a impulsos de choque
se eliminan de la misma en base a un filtrado adaptativo combinado
con el cálculo de la relación de probabilidad de impulsos de choque
en una prueba específica e intervalos de referencia. Además, la
detección de cada impulso de choque potencial se realiza en el
dominio de tiempo, individualmente, uno a uno, al contrario de otras
técnicas en las que la detección de impulso de choque se basa en el
análisis de la señal desmodulada en el dominio de frecuencia. En
consecuencia, la sensibilidad de la solución de la invención es
mayor. Como resultado del algoritmo de detección de impulso de
choques, se crea el vector binario p que contiene solamente
componentes de señal asociados a impulsos de choque. Para la
identificación de un fallo en el cojinete, el vector binario p es
procesado en el dominio de tiempo mediante un procedimiento de
probabilidad máxima de alta resolución. No se utiliza en absoluto
la transformada de Fourier en el proceso de identificación de
fallo.
El procedimiento de acuerdo con la invención se
presenta como realización en un dibujo en el que la figura 1
muestra un diagrama de bloques del dispositivo utilizado para
detectar e identificar automáticamente defectos en rodamientos, la
figura 2 es una onda temporal de una señal procesada X con
intervalos de prueba y de referencia indicados, la figura 3 es una
gráfica con la aparición de impulsos de choque en forma binaria del
vector p frente al fondo con una onda temporal de la señal
procesada X, la figura 4 es un fragmento del espectro de
frecuencias en el cual se ha indicado la frecuencia f_{fault} del
defecto producido junto con la frecuencia de prueba v_{fk}, la
figura 5 es un diagrama de operaciones realizadas al detectar
impulsos de choque de acuerdo con el presente procedimiento y la
figura 6 es un diagrama de operaciones realizas al identificar
defectos en cojinetes de acuerdo con el presente procedimiento.
Un rodamiento 2, dispuesto en un eje giratorio
1, comprende unos puntos de medición, marcados en su carcasa, donde
se coloca un sensor de vibración 3 para medir las vibraciones del
cojinete y el sensor de velocidad 4 del eje 1. Los sensores están
conectados a un dispositivo de medición 5 utilizado para medir la
señal analógica de la tensión proporcional a las vibraciones del
cojinete en funcionamiento 2 y para medir la velocidad de giro del
eje. La salida del dispositivo de medición 5 se conecta a la entrada
de un procesador 6 en cuya memoria 7 puede distinguirse un módulo
de procesamiento 8 y un módulo de función 9. El módulo de
procesamiento 8 incorpora un sistema de procesamiento 10 en el cual
la señal medida se convierte a señal digital, que después se envía
como datos de entrada a un sistema de filtro de paso alto 11
incorporado también en el módulo 8, en el cual se filtran los
componentes de la señal indeseados en un análisis adicional. Los
datos de salida recibidos desde el sistema 11 se envían a un
detector de impulsos de choque 12 lo cual realiza uno de los
sistemas del módulo de función 9 del procesador 6. Los cálculos
matemáticos asociados a datos de configuración 13 realizados en el
módulo de función 9 y los datos de entrada 14 que son las
frecuencias conocidas de los defectos de los rodamientos analizados
2 se envían también al módulo de función 9. Los datos de salida del
sistema del detector de impulsos de choque 12 son transmitidos a un
identificador de fallos 15 que pertenece al módulo de función 9 y
en el cual se realiza la identificación de los defectos de los
rodamientos. La salida del procesador 6 se conecta a un dispositivo
para visualizar y/o registrar los resultados obtenidos 16, que puede
ser una pantalla de ordenador, una impresora, o cualquier soporte
de grabación de datos.
\vskip1.000000\baselineskip
El procedimiento de la invención se implementa
de acuerdo con las siguientes acciones mostradas en la figura
6.
- -
- Medición de la vibración y la velocidad de giro. Por medio del sensor de vibración 3 y el dispositivo de medición 5 las vibraciones del rodamiento operativo 2 se miden en los puntos de medición indicados 3, y por medio del sensor de velocidad 4 y el dispositivo de medición 5 se mide 1 la velocidad de giro del eje apoyado en el cojinete 1 (actividad 101).
- -
- Conversión analógico a digital. Los resultados de las mediciones se convierten de forma analógica a datos digitales en la unidad de proceso 10, o en el dispositivo de medición 5, y después estos resultados se almacenan en forma digital en la memoria 7 del procesador 6 del ordenador como onda temporal de la señal (actividad 102).
- -
- Filtrado de datos digitales. Los datos almacenados son filtrados previamente en el sistema de filtrado de paso alto 11 produciendo una onda temporal X, donde X= x [k, k+1, k+2, ..... k+K], donde k indica muestras de instantes sucesivos de una señal procesada desde el intervalo de 1 al K. Como sistema de filtrado de paso alto, puede utilizarse, por ejemplo, un filtro FIR (de respuesta finita al impulso). La característica del filtro FIR es que su respuesta a una secuencia determinada de muestras de entrada tiene una longitud finita. El objetivo del filtrado de datos es eliminar señales de una frecuencia inferior a un límite preestablecido constituyendo uno de los parámetros de configuración (actividad 103).
- -
- Filtrado adaptativo y la determinación de relaciones de probabilidad de impulsos de choque M_{S}. Los datos obtenidos como resultado del filtrado previo se transmiten al detector de impulsos de choque 12 en el cual la onda temporal de la señal experimenta un filtrado adaptativo y se determinan las relaciones de probabilidad de impulso de choques M_{S}. En el detector de impulsos de choque 12, al principio de la onda temporal filtrada previamente X, se selecciona un intervalo de referencia adyacente de una longitud K_{r} y un intervalo de prueba de una longitud K_{t}. La longitud del intervalo de referencia K_{r} es menor que el período de tiempo de un defecto de cojinete conocido con una frecuencia máxima derivada de los datos de entrada (14), y siendo la longitud del intervalo de prueba K_{t} menor o igual que al menos la mitad de la longitud del intervalo de referencia. Preferiblemente, se supone que la longitud del intervalo de prueba es un 15% de la longitud del intervalo de referencia. Después, se adopta un número definitivo de filtros N_{S}, preferiblemente N_{S}= 5 y para cada filtrado L^{n}, donde L es un número natural del intervalo de 1 a N_{S}, se realiza el filtrado adaptativo de una señal a partir de los intervalos de prueba y de referencia, tras lo cual para cada filtrado L^{n} se calcula la relación de probabilidad de impulsos de choque M_{s} de acuerdo con la siguiente la relación (actividad 104):
\vskip1.000000\baselineskip
- donde \tau_{t} [L] es el valor medio cuadrático de la señal en el intervalo de prueba examinado sometido a filtrado adaptativo L^{n}, definido de la siguiente manera:
- y \tau_{r} [L] es el valor medio cuadrático de la señal en el intervalo de referencia examinado sometido a filtrado adaptativo L^{n}, definido de la siguiente manera:
\vskip1.000000\baselineskip
\vskip1.000000\baselineskip
- donde:
- L es un número natural desde el intervalo 1... N_{S} e indica el número de filtros adaptativos,
- K es la longitud total de la onda temporal examinada de la señal procesada,
- K_{t}, K_{r} son las longitudes de los intervalos de prueba y de referencia, respectivamente,
- k_{t}, k_{r} son números naturales sucesivos de 1 a K_{t} o K_{r},
- z_{t}^{L}, indica un vector que representa el intervalo de prueba examinado de la señal sometida a un filtrado adaptativo L^{n},
- z_{r}^{L}, indica un vector que representa el intervalo de referencia examinado de la señal sometida a un filtrado adaptativo L^{n},
- -
- Verificación de la aparición del impulso de choque en el intervalo de prueba. Para cada filtrado adaptativo L^{n} el valor de la relación de probabilidad de impulso de choque M_{S} se compara con un valor límite preestablecido S_{L}, previamente enviado a la memoria del ordenador, y si la relación de probabilidad M_{S} es mayor que S_{L} en por lo menos S_{1} filtros, esto significa que se produjo un impulso de choque en el intervalo de prueba examinado, mientras que si la relación de probabilidad M_{S} es mayor que S_{L} en menos de S_{1} filtros, esto significa que no se produjo impulso de choque en el intervalo de prueba examinado, registrándose la aparición del impulso de choque en el intervalo de prueba una forma binaria como vector p_{n}, llenando los elementos iniciales del vector p_{n} con números "1" o números "0" para la no aparición de un impulso de choque (acción 105).
- -
- Determinación de elementos restantes del vector binario p_{n}. En la onda temporal X de la señal procesada se determinan intervalos de prueba y de referencia sucesivos adyacentes desplazados respecto a los determinados previamente mediante una etapa preestablecida de manera reiterada y actividades conectadas al filtrado adaptativo, determinación de relaciones de probabilidad de impulsos de choque M_{S} para todos los intervalos de prueba sucesivos, para un número de filtrados preestablecidos N_{S}, y verificación de la aparición de un impulso de choque para todos los intervalos de toda la onda temporal X de la señal procesada se repiten y de este modo se obtienen el resto de elementos del vector binario obtenido p_{n}. El vector binario P_{n} está almacenado en la memoria 7 del procesador 6 como onda temporal binaria que adopta valores 1 en la zona de aparición de impulsos de choque y 0 en la zona donde no se producen impulsos de choque (actividad 106).
- -
- Determinación del vector p_{r}. Se realiza una reversión del orden de muestras de tiempo en la onda temporal X almacenada filtrada previamente de la señal procesada de los intervalos de prueba y de referencia y para dicha onda temporal de la señal invertida se realizan todas las actividades anteriores de la 104 a la 106 como antes y la información en presencia de impulsos de choque se almacena como un vector binario p_{r}, donde el valor 1 indica la presencia de un impulso de choque para la muestra de tiempo dada, y el valor 0 indica su ausencia (actividad 107).
- -
- Determinación del vector p. El orden de los elementos en el vector binario p_{r} se invierte y se suma lógicamente con el vector p_{n} para obtener el vector binario p. El vector binario p se almacena en la memoria 7 del procesador 6 (actividad 109) como onda temporal binaria que adopta los valores 1 en la zona de aparición de impulsos de choque y 0 en la zona donde no se producen impulsos de choque (actividad 108).
- -
- Guardar el vector p en la memoria del procesador. Los resultados obtenidos se almacenan en un equipo especial 18 (actividad 109).
- -
- Después, los elementos del vector p se utilizan como datos preliminares para la identificación automática de defectos en el cojinete (actividad 110).
Para determinar el valor medio cuadrático de la
señal en el intervalo de referencia \tau_{r} [L], sometido al
filtrado adaptativo L^{n}, que es necesario para determinar la
relación de probabilidad de impulsos de choque M_{S}, se utilizan
relaciones matemáticas en las cuales se define la matriz de
autocovarianza C_{r} como la relación entre matrices ortogonales
V y U^{T} y de la matriz \Sigma_{r}. Los elementos diagonales
de la matriz \Sigma_{r}. son distintos de cero, no negativos e
indican valores medios cuadráticos \tau_{r} [L]. Esta relación
se muestra en la siguiente fórmula: (4)
(4)C_{r}= V
\Sigma_{r}
U^{T}
donde la matriz \Sigma_{r}
viene definida por la siguiente
relación:
\vskip1.000000\baselineskip
Al mismo tiempo la matriz autocovarianza C_{r}
se define como sigue:
para
cual:
donde:
x_{r} [k_{r}] son los elementos del
intervalo de referencia,
k_{r} indica números naturales sucesivos de 1
a K_{r},
i_{r} indica números naturales sucesivos de 1
a L,
K_{r} indica la longitud de la señal de
referencia.
Como resultado de la comparación de las
relaciones 4 y 6, utilizando la descomposición en valores singulares
(SVD) de la matriz de autocovarianza, se obtiene el valor
cuadrático medio de la señal en el intervalo de referencia
\tau_{r} [L].
Después, para determinar el valor cuadrático
medio de una señal en el intervalo de prueba \tau_{t} [L],
sometido al filtrado adaptativo L^{n}, cuyo valor es necesario
para determinar la probabilidad de impulso del choque M_{S}, se
utilizan las siguientes relaciones matemáticas:
donde
Z^{L}_{t} son los elementos de la matriz
Z_{t} presentados de esta manera:
y calculados de acuerdo con esta
relación:
(9)Z_{t}=
X_{t}
U
donde la matriz X_{t} se define
como:
y la matriz U es una matriz
ortogonal cuyos elementos son los coeficientes de filtrado
adaptativo, mientras
que:
K_{t} indica la longitud del intervalo de
prueba examinado;
x_{t} [K_{t}] son los elementos del
intervalo de prueba;
k_{t} indica números naturales sucesivos de 1
a K_{t},
N_{t}= K_{t} - L
\vskip1.000000\baselineskip
Entonces, la identificación del defecto en un
cojinete o cojinetes se realiza automáticamente en el identificador
de defectos 15 para un vector binario p determinado a través del
procedimiento de la invención y para la base de los datos de
entrada que son frecuencias de fallo conocidas del cojinete o
cojinetes examinados F_{f}= (f_{f1}, f_{f2}, ..... f_{fk}),
haciendo uso de las siguientes actividades representadas en la
figura 6.
- -
- Se determina un vector binario de decimación p_{d} sometiendo al vector binario p a decimación con un coeficiente de decimación preestablecido para obtener el vector binario de decimación p_{d} (actividad 111).
- -
- Se determinan las funciones de probabilidad máxima J para frecuencias sucesivas f_{fK} de defectos en el cojinete analizado. Para cada frecuencia sucesiva de defectos conocidos f_{fK} del conjunto F_{f}, del cojinete analizado, la frecuencia de prueba v_{fK} \in | (1-P_{g}) f_{fK}- (1+P_{g}) f_{fK} |) se determina de manera aproximada y después se determina de manera exacta a partir de la proximidad P_{g} de la frecuencia de defectos conocidos f_{fK}, donde P_{g} es un parámetro de configuración que indica el intervalo de búsqueda respecto a la frecuencia de defectos determinada. Después, utilizando procedimientos matemáticos conocidos, se determina la probabilidad de la función máxima J para cada la frecuencia de prueba v_{fK} utilizando el vector binario de decimación (actividad 112).
- -
- Búsqueda del impulso de choque con mayor amplitud. El valor máximo J_{A} de la función de probabilidad máxima J se determina comprobando para qué frecuencia de prueba v_{fK} la función de probabilidad máxima J alcanza el valor máximo J_{A} (actividad 113).
- -
- Generación de información sobre la aparición o ausencia de defectos en el cojinete examinado. Se calcula la relación F_{A} del valor máximo J_{A} respecto al nivel de ruido n_{A} y su resultado se compara con el parámetro límite preestablecido S_{A}. Si la relación es mayor que el parámetro límite preestablecido, entonces se genera automáticamente información de que se ha producido un defecto específico en el cojinete examinado, por ejemplo f_{f1}, correspondiente a la frecuencia de prueba v_{fK}.
Si la relación F_{A} es menor o igual que el
parámetro límite preestablecido S_{A}, significa que ninguno de
los impulsos de choque encontrados presentados en forma de vector
binario corresponde a la frecuencia examinada que sería similar a
la frecuencia de prueba de los defectos de los cojinetes analizados,
y entonces se genera automáticamente información de que no se ha
producido ningún defecto específico en el cojinete o cojinetes
analizados (actividad 114).
Para determinar el valor máximo J_{A} de la
función de probabilidad máxima J se utilizan las siguientes
relaciones matemáticas:
donde:
N_{h} es el número preestablecido de los
armónicos examinados,
y A_{k} es una amplitud armónica determinada
en el dominio de tiempo utilizando el vector binario p, presentado
en forma de modelo sinusoidal.
El vector binario p se presenta como la suma de
modelo sinusoidal y vector de ruido:
El cual puede presentarse también en forma de la
siguiente matriz:
(13)p=
B\theta +
e
donde:
\phi_{k} es una fase,
\omega_{k} = 2\pivk/f_{s}, k=1,
2... N_{h} es la velocidad angular correspondiente a una
frecuencia específica (v), k =1, 2... N_{h},
f_{s} es la frecuencia de muestreo de
la señal analógica medida cuando se realiza la conversión analógico
a digital,
e = [\varepsilon [n]]^{T} indica el vector
de ruido,
n indica un número natural de 1 a
N-1,
N indica el número de elementos en el vector
binario p,
B es la matriz de parámetros constantes
conocidos, presentada de la siguiente manera:
y \theta es el vector de
parámetros desconocidos, presentado de la siguiente
manera:
\theta =
[A_{1} cos(\phi_{1}), A_{1} sin(\phi_{1}),
A_{2} cos(\phi_{2}), A_{2}
sin(\phi_{2}),...,A_{Nh} cos(\phi_{Nh}),
A_{Nh}
sin(_{Nh})],
que es el estimador de mínimos
cuadrados para la
relación:
(15)\theta
(B^{T}B)^{-1}B^{T}p
donde p es un vector binario de la
siguiente
forma:
p= [p
[0], p [1],..., p
[N-1]]^{T}
La comparación de las relaciones 12 y 13 permite
obtener los valores de amplitudes armónicas A_{k}, necesarios
para calcular la relación J_{A} que se utiliza después para
determinar el lugar de la aparición del defecto, es decir, qué
cojinete puede ser defectuoso, y el tipo del defecto que se
produce.
El nivel de ruido n_{A} se calcula utilizando
el valor o valores de la función máxima J encontrados en las
proximidades de la frecuencia del defecto examinado f_{fK}.
El resultado final de las actividades realizadas
es la generación automática de información acerca de la presencia o
ausencia de defectos en el cojinete o cojinetes examinados. La
aparición de un defecto provoca una generación automática de
información acerca de cuál de los cojinetes examinados ha sido
dañado y cuál es el tipo de defecto.
Claims (11)
1. Procedimiento para la detección y la
identificación automática de daños a rodamientos, en el que las
vibraciones de un cojinete en funcionamiento, medidas como una
señal de onda temporal analógica se convierte en datos digitales y
se filtran previamente de manera conocida, caracterizado por
el hecho de que se realiza una detección de impulsos de choque con
el uso de un vector binario (p) que se determina en las siguientes
acciones:
- -
- en la onda temporal filtrada previamente de la señal procesada se seleccionan dos intervalos iniciales adyacentes: uno es un intervalo de referencia y el otro es un intervalo de prueba y la señal presente en los intervalos de prueba y de referencia se someten a múltiples filtrados tras lo cual después de cada filtrado, se calcula una relación de probabilidad de impulsos de choque definida como una relación entre los valores medios cuadráticos de la señal filtrada presentes en el intervalo de prueba y los valores medios cuadráticos de la señal filtrada presentes en el intervalo de referencia,
- -
- las relaciones de probabilidad calculadas se comparan con un valor límite preestablecido y se estima la probabilidad de aparición de impulsos de choque en el intervalo de prueba inicial,
- -
- la información sobre la aparición de impulsos de choque en el intervalo de prueba examinado se registra como elementos iniciales del vector binario (p), y después
- -
- todos los elementos restantes del vector binario (p) se determinan para todos los intervalos de prueba y de referencia sucesivos que son seleccionados de manera reiterada con un desplazamiento en una dirección específica y aquellos elementos del vector binario (p) adoptan valores 1 en la zona de aparición del impulso de choque y 0 en la zona en la que no se producen impulsos de choque.
\vskip1.000000\baselineskip
2. Procedimiento según la reivindicación 1,
caracterizado por el hecho de que dicho valor cuadrático
medio de la señal en el intervalo de prueba examinado sometido a
filtrado L^{n} (\tau_{t} [L]) se determina en base a la
siguiente relación:
\vskip1.000000\baselineskip
\vskip1.000000\baselineskip
mientras que el valor cuadrático
medio de la señal en el intervalo de referencia examinado sometido a
filtrado L^{n} (\tau_{r} [L]) se determina en base a la
siguiente
relación:
\vskip1.000000\baselineskip
\vskip1.000000\baselineskip
donde:
L es un número natural del intervalo 1...
N_{S} e indica el número de filtros,
K_{t} indica la longitud del intervalo de
prueba examinado,
K_{r} indica la longitud del intervalo de
referencia examinado,
k_{t}, k_{r} indica los números naturales
sucesivos 1 al Kt o al Kr,
z^{L}_{t} indica un vector que representa el
intervalo de prueba examinado de la señal sometido a filtrado
L^{n}
z^{L}_{r} indica un vector que representa el
intervalo de referencia examinado de la señal sometido a filtrado
L^{n}.
3. Procedimiento según la reivindicación 2,
caracterizado por el hecho de que dicho valor cuadrático
medio de una señal en el intervalo de referencia examinado
(\tau_{r} [L]) se determina usando filtrado L^{n}, que es un
filtrado adaptativo, con uso de descomposición en valores singulares
(SVD) de la matriz de autocovarianza (Cr) de acuerdo con la
siguiente relación:
C_{r} = V
\sum_{r}
U^{T},
donde la matriz (\sum_{r}) es
una matriz diagonal cuyos elementos diagonales son distintos de cero
y no negativos e indican valores medios cuadráticos (\tau_{r}
[L]):
y las matrices (V) y (U) son
matrices ortogonales y al mismo tiempo la matriz de autocovarianza
(Cr) presentan la siguiente
forma:
para la
cual:
donde:
x_{r} [k_{r}] son los elementos del
intervalo de referencia,
k_{r} indica números naturales sucesivos de 1
a K_{r},
i_{r} indica números naturales sucesivos de 1
a L,
K_{r} indica la longitud del intervalo de
referencia.
L indica el número de filtros adaptativos.
\vskip1.000000\baselineskip
4. Procedimiento según la reivindicación 2,
caracterizado por el hecho de que dicho valor cuadrático
medio de una señal en el intervalo de prueba examinado
(\tau_{r} [L]) se determina utilizando filtrado L^{n}, que es
un filtrado adaptativo, en base a la siguiente relación:
donde z^{L}_{t} son los
elementos de la matriz Z_{t} presentada de la siguiente
manera:
y calculada de acuerdo con la
siguiente
relación:
Z_{t} =
X_{t}
U,
\newpage
donde la matriz (X_{t}) se define
como:
y la matriz (U) es una matriz
ortogonal cuyos elementos son los coeficientes del filtro
adaptativo, mientras
que:
K_{t} indica la longitud del intervalo de
prueba examinado;
x_{t} [k_{t}], son los elementos del
intervalo de prueba;
k_{t} indica los números naturales sucesivos
de 1 a K_{t},
L indica el número de filtros adaptativos.
\vskip1.000000\baselineskip
5. Procedimiento según la reivindicación 1 o 2,
caracterizado por el hecho de que el filtrado múltiple de
señales en los intervalos de prueba y de referencia se realiza con
el uso de un filtro de paso de banda.
6. Procedimiento según cualquiera de las
reivindicaciones anteriores, caracterizado por el hecho de
que la longitud del intervalo de prueba (K_{t}) es un 15% de la
longitud del intervalo de referencia (K_{r}).
7. Procedimiento según cualquiera de las
reivindicaciones anteriores, caracterizado por el hecho de
que la identificación del defecto en el cojinete o cojinetes se
realizad automáticamente para el vector binario (p) determinado
para todos los intervalos de prueba y de referencia examinados de la
onda temporal (X) de la señal procesada y el conjunto de
frecuencias de defectos conocidos del cojinete o cojinetes
examinados ((F_{F}) = [f_{f1}, f_{f1},..., f_{fK}]), por
medio de las siguientes acciones:
- -
- para cada frecuencia sucesiva de defectos conocidos (f_{fK}) del conjunto (F_{t}) del cojinete examinado, se determina la frecuencia de prueba v_{fK} \in | (1-P_{g}) f_{fK}- (1+P_{g}) f_{fK} |) que pertenece a la proximidad (P_{g}) de la frecuencia defectos conocidos (f_{fK}), donde (P_{g}) es un parámetro de configuración que indica la gama de frecuencias de búsqueda relacionada con la frecuencia de defectos conocidos y después, utilizando principios matemáticos de probabilidad máxima conocidos, se determina para cada frecuencia de prueba (v_{fK}) la función de probabilidad máxima (J) en dominio de tiempo utilizando el vector binario (p),
- -
- después se comprueba para qué frecuencia de prueba (v_{fK}) la función máxima de la probabilidad (J) alcanza el valor máximo (J_{A}),
- -
- se calcula la relación (F_{A}) entre el valor máximo (J_{A}) y el nivel de ruido (n_{A}) y el resultado se compara con el parámetro límite preestablecido (S_{A}), donde el nivel de ruido (n_{A}) se calcula con el uso de valores de función máxima (J),
- -
- si la relación (F_{A}) es mayor que el parámetro límite preestablecido (S_{A}), se genera automáticamente información sobre el defecto del cojinete examinado con una frecuencia específica que pertenece al conjunto (F_{t}) y corresponde a la frecuencia de prueba (v_{fK}),
- -
- si la relación (F_{A}) es menor o igual que el parámetro límite preestablecido (S_{A}), se supone que ninguno de los impulsos de choque presentes en el vector binario (p) corresponde a ninguna de las frecuencias de defectos conocidos (f_{fK}) de cojinetes examinados y se genera automáticamente información de que no se produjo ningún defecto específico en el cojinete o cojinetes examinados.
\vskip1.000000\baselineskip
8. Procedimiento según la reivindicación 7,
caracterizado por el hecho de que dicha frecuencia de prueba
(v_{fK}), para la cual se calcula la función máxima de la
probabilidad para el vector binario (p) en dominio de tiempo
alcanza el valor máximo (J_{A}) se determina maximizando la suma
cuadrática de amplitudes armónicas (A_{k}) de acuerdo con esta
relación:
donde:
N_{h} es un número preestablecido de los
armónicos examinados.
A_{k} una amplitud armónica determinada en el
dominio de tiempo utilizando el vector binario p, presentado en
forma de modelo sinusoidal.
\vskip1.000000\baselineskip
9. Procedimiento según la reivindicación 8,
caracterizado por el hecho de que dichos valores de las
amplitudes de los armónicos (A_{k}) se determina en el dominio de
tiempo utilizando el vector binario (p), presentado en forma de
modelo sinusoidal en el que el vector binario (p) se presenta como
una suma del modelo sinusoidal y el vector de ruido:
qué puede presentarse también en
forma de la siguiente
matriz:
p =
B\theta +
e
donde:
\phi_{k} es una fase,
\omega_{k}=2\pivk/f_{s}, k=1,
2... N_{h} es la velocidad angular correspondiente a una
frecuencia específica (v), k =1, 2... N_{h},
f_{s} es la frecuencia de muestreo de
la señal analógica medida cuando se realiza la conversión analógico
a digital,
e = [\varepsilon [n]]^{T} indica el vector
de ruido,
n indica un número natural de 1 a
N-1,
N indica el número de elementos en el vector
binario (p),
B es la matriz de parámetros constantes
conocidos, presentada de la siguiente manera:
y \theta es el vector de
parámetros desconocidos, presentado de la siguiente
manera:
\theta =
[A_{1} cos(\phi_{1}), A_{1} sin(\phi_{1}),
A_{2} cos(\phi_{2}), A_{2}
sin(\phi_{2}),...,A_{Nh} cos(\phi_{Nh}),
A_{Nh}
sin(_{Nh})],
que es el estimador de mínimos
cuadrados para la
relación:
\theta
(B^{T}B)^{-1}B^{T}p
donde (p) es un vector binario de
la siguiente
forma:
p= [p
[0], p [1],..., p
[N-1]]^{T}
\vskip1.000000\baselineskip
10. Procedimiento según la reivindicación 7,
caracterizado por el hecho de que el nivel de ruido (n_{A})
se calcula tomando un valor medio de los valores de función máxima
(J) que se encuentran en las proximidades de la frecuencia del
defecto examinado (f_{fK}) del conjunto (F_{f}).
11. Procedimiento según la reivindicación 7,
caracterizado por el hecho de que dicho nivel de ruido
(n_{A}) se calcula tomando un valor medio de los valores de
función máxima (J) que se encuentran en las proximidades de la
frecuencia del defecto examinado (f_{fK}) del conjunto
(F_{f}).
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP07460032A EP2053375B1 (en) | 2007-10-24 | 2007-10-24 | A method for detection and automatic identification of damage to rolling bearings |
Publications (1)
Publication Number | Publication Date |
---|---|
ES2336842T3 true ES2336842T3 (es) | 2010-04-16 |
Family
ID=39154109
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
ES07460032T Active ES2336842T3 (es) | 2007-10-24 | 2007-10-24 | Procedimiento para la deteccion e identificacion automatica de daños en rodamientos. |
Country Status (11)
Country | Link |
---|---|
US (1) | US8285498B2 (es) |
EP (1) | EP2053375B1 (es) |
KR (1) | KR101436190B1 (es) |
CN (1) | CN101836100B (es) |
AT (1) | ATE453856T1 (es) |
DE (1) | DE602007004123D1 (es) |
DK (1) | DK2053375T3 (es) |
ES (1) | ES2336842T3 (es) |
PL (1) | PL2053375T3 (es) |
RU (1) | RU2470280C2 (es) |
WO (1) | WO2009053076A1 (es) |
Families Citing this family (35)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009012354A1 (en) * | 2007-07-17 | 2009-01-22 | Clemson University | System and method to assess signal similarity with applications to diagnostics and prognostics |
BR112012021738B1 (pt) | 2010-03-01 | 2020-02-04 | Siemens Ag | método e dispositivo para a detecção precoce do desenvolvimento de danos em um mancal |
JP5442553B2 (ja) * | 2010-07-21 | 2014-03-12 | 日立建機株式会社 | 軸受の損傷検出方法 |
TWI438416B (zh) * | 2011-04-14 | 2014-05-21 | Delta Electronics Inc | 使用連續遞移轉換之信號分析系統及方法 |
WO2013044973A1 (de) * | 2011-09-30 | 2013-04-04 | Siemens Aktiengesellschaft | Verfahren und anordnung zum ermitteln und/oder überwachen des zustands eines wälzlagers |
CN102426189A (zh) * | 2011-10-21 | 2012-04-25 | 沈阳黎明航空发动机(集团)有限责任公司 | 一种轴承脉冲检测方法 |
CN102519726B (zh) * | 2011-12-28 | 2015-06-03 | 昆明理工大学 | 一种针对滚动轴承复合故障的声学诊断方法 |
US8811136B2 (en) * | 2012-11-19 | 2014-08-19 | Lsi Corporation | Harmonic ratio based defect classifier |
CN104330258A (zh) * | 2014-10-23 | 2015-02-04 | 徐州隆安光电科技有限公司 | 一种基于特征参量的滚动轴承故障灰色关联度辨识方法 |
CN104596766B (zh) * | 2014-12-24 | 2017-02-22 | 中国船舶工业系统工程研究院 | 一种轴承早期故障确定方法及装置 |
CN105067248B (zh) * | 2015-08-17 | 2018-04-27 | 安徽容知日新科技股份有限公司 | 设备转速和振动数据采集方法、装置和监测系统 |
KR101667164B1 (ko) * | 2015-10-02 | 2016-10-17 | 경북대학교 산학협력단 | 회전체 속도 기반의 베어링 고장 진단 방법 |
CN105606363B (zh) * | 2016-01-29 | 2017-11-07 | 济南大学 | 一种基于域自适应的轴承故障诊断方法 |
US10082443B2 (en) | 2016-02-26 | 2018-09-25 | General Electric Technology Gmbh | System and method for monitoring bearing health in a journal assembly |
CN106017927B (zh) * | 2016-05-18 | 2018-09-11 | 唐智科技湖南发展有限公司 | 一种识别轴承滚子间距变化的保持架故障诊断方法 |
ITUA20163745A1 (it) * | 2016-05-24 | 2017-11-24 | Nuovo Pignone Tecnologie Srl | Metodo e sistema per monitorare lo stato di salute di un cuscinetto a rotolamento di una macchina, e macchina equipaggiata con tale sistema |
EP3255776A1 (en) * | 2016-06-07 | 2017-12-13 | ABB Technology AG | A method and device for determination of a torsional deflection of a rotation shaft in the electromechanical drivetrain |
CN106485070B (zh) * | 2016-09-30 | 2019-05-28 | 广州机智云物联网科技有限公司 | 一种自适应阈值调整方法 |
EP3309530A1 (en) | 2016-10-11 | 2018-04-18 | ABB Schweiz AG | Detection of bearing fault |
RU2659101C1 (ru) * | 2017-06-05 | 2018-06-28 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Владимирский Государственный Университет имени Александра Григорьевича и Николая Григорьевича Столетовых" (ВлГУ) | Способ определения момента трения в подшипниках качения |
CN107149485B (zh) * | 2017-06-07 | 2020-03-06 | 青岛海信医疗设备股份有限公司 | 基于医学的超声波信号处理方法及装置 |
JP2019045241A (ja) * | 2017-08-31 | 2019-03-22 | Ntn株式会社 | 転がり軸受の状態監視方法および状態監視装置 |
RU2680640C1 (ru) * | 2018-02-16 | 2019-02-25 | Владимир Сергеевич Давыдов | Способ вибродиагностики зарождающихся дефектов механизмов |
CN110030959B (zh) * | 2018-06-08 | 2020-12-25 | 上海微小卫星工程中心 | 用于确定转动机构的周向磨损位置的装置和方法 |
CN111553178A (zh) * | 2019-02-12 | 2020-08-18 | 斯凯孚公司 | 转动机械振动特性的智能识别方法 |
CN110174269B (zh) * | 2019-05-06 | 2021-02-26 | 江苏联能电子技术有限公司 | 变分模态分解和共振解调分析核电用泵冲击特征的方法 |
CN110307981B (zh) * | 2019-06-17 | 2021-01-26 | 西安工程大学 | 基于pnn-ifa的轴承故障诊断方法 |
CN110988130A (zh) * | 2019-12-18 | 2020-04-10 | 成都理工大学 | 一种用于岩质斜坡室内试验损伤识别的测试系统 |
CN111969639B (zh) * | 2020-08-17 | 2022-03-11 | 中国科学院电工研究所 | 电力电子化电网级联型换流器多时间尺度暂态建模方法 |
CN112561306B (zh) * | 2020-12-11 | 2023-12-08 | 领伟创新智能系统(浙江)有限公司 | 一种基于Hankel矩阵的滚动轴承健康状态评价方法 |
CN113405823B (zh) * | 2021-05-17 | 2022-05-20 | 同济大学 | 一种基于迭代扩展本征模态分解的旋转机械故障诊断方法 |
RU2769919C1 (ru) * | 2021-05-27 | 2022-04-08 | Федеральное государственное автономное образовательное учреждение высшего образования "Санкт-Петербургский государственный электротехнический университет "ЛЭТИ" им. В.И. Ульянова (Ленина)" (СПбГЭТУ "ЛЭТИ") | Способ вибродиагностики возникновения зарождающихся дефектов в отдельных узлах механизмов |
CN114441174B (zh) * | 2022-02-09 | 2024-07-09 | 上海电气集团股份有限公司 | 滚动轴承复合故障的诊断方法、系统、设备及介质 |
CN115730199B (zh) * | 2022-11-10 | 2023-07-21 | 天地(常州)自动化股份有限公司北京分公司 | 一种滚动轴承振动信号降噪和故障特征提取方法和系统 |
CN116150635B (zh) * | 2023-04-18 | 2023-07-25 | 中国海洋大学 | 基于跨域关联性表示的滚动轴承未知故障检测方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3208268A (en) * | 1961-10-05 | 1965-09-28 | Skf Ind Inc | Detection of almost periodic occurrences |
SU783620A1 (ru) * | 1979-02-26 | 1980-11-30 | Каунасский Политехнический Институт Им. Антанаса Снечкуса | Способ диагностики подшипников качени |
SU951942A1 (ru) * | 1980-05-16 | 2006-04-10 | Г.В. Морозов | Способ обнаружения повреждений подшипников качения |
US5477730A (en) | 1993-09-07 | 1995-12-26 | Carter; Duncan L. | Rolling element bearing condition testing method and apparatus |
FI112972B (fi) | 1998-07-15 | 2004-02-13 | Abb Research Ltd | Laakerin kunnon arviointi |
US6591682B1 (en) * | 2000-08-14 | 2003-07-15 | Pruftechnik Dieter Busch Ag | Device and process for signal analysis |
CN2618156Y (zh) * | 2003-05-09 | 2004-05-26 | 康乃正 | 轴承振动实时记录装置 |
JP3875981B2 (ja) * | 2004-03-17 | 2007-01-31 | 新日本製鐵株式会社 | ころがり軸受の異常診断方法および装置 |
-
2007
- 2007-10-24 DE DE602007004123T patent/DE602007004123D1/de active Active
- 2007-10-24 EP EP07460032A patent/EP2053375B1/en not_active Not-in-force
- 2007-10-24 DK DK07460032.1T patent/DK2053375T3/da active
- 2007-10-24 PL PL07460032T patent/PL2053375T3/pl unknown
- 2007-10-24 ES ES07460032T patent/ES2336842T3/es active Active
- 2007-10-24 AT AT07460032T patent/ATE453856T1/de not_active IP Right Cessation
-
2008
- 2008-10-17 CN CN2008801138280A patent/CN101836100B/zh not_active Expired - Fee Related
- 2008-10-17 RU RU2010120713/28A patent/RU2470280C2/ru not_active IP Right Cessation
- 2008-10-17 WO PCT/EP2008/008988 patent/WO2009053076A1/en active Application Filing
- 2008-10-17 US US12/738,462 patent/US8285498B2/en active Active
- 2008-10-17 KR KR1020107008456A patent/KR101436190B1/ko active IP Right Grant
Also Published As
Publication number | Publication date |
---|---|
RU2470280C2 (ru) | 2012-12-20 |
RU2010120713A (ru) | 2011-11-27 |
ATE453856T1 (de) | 2010-01-15 |
US8285498B2 (en) | 2012-10-09 |
PL2053375T3 (pl) | 2010-05-31 |
KR101436190B1 (ko) | 2014-09-01 |
DK2053375T3 (da) | 2010-04-12 |
EP2053375A1 (en) | 2009-04-29 |
KR20100094452A (ko) | 2010-08-26 |
CN101836100A (zh) | 2010-09-15 |
EP2053375B1 (en) | 2009-12-30 |
WO2009053076A1 (en) | 2009-04-30 |
US20100280772A1 (en) | 2010-11-04 |
CN101836100B (zh) | 2012-07-18 |
DE602007004123D1 (de) | 2010-02-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
ES2336842T3 (es) | Procedimiento para la deteccion e identificacion automatica de daños en rodamientos. | |
Dalpiaz et al. | Effectiveness and sensitivity of vibration processing techniques for local fault detection in gears | |
Cong et al. | Short-time matrix series based singular value decomposition for rolling bearing fault diagnosis | |
Shibata et al. | Fault diagnosis of rotating machinery through visualisation of sound signals | |
CN107209512B (zh) | 用于机器的声学或振动分析的学习阶段的方法、系统及计算机程序 | |
EP0776464B1 (en) | Method for determining resonance information | |
Li et al. | Acoustic emission analysis for bearing condition monitoring | |
BR112012010702B1 (pt) | Processo de detecção de um dano de pelo menos um rolamento de mancal | |
Hadjileontiadis et al. | Crack detection in beams using kurtosis | |
CN105067239B (zh) | 基于扫频激励振动的梁裂纹故障检测装置及方法 | |
CN108106838B (zh) | 一种机械故障精密诊断的整周期时域回归及报警方法 | |
CN107449932A (zh) | 一种航空发动机转子转速的测量方法 | |
CN109682463B (zh) | 变电站可听噪声测量方法 | |
US5031459A (en) | Turbine generator shaft torsion monitor | |
CN108414217A (zh) | 齿轮箱噪音测试系统 | |
RU2461847C2 (ru) | Способ непрерывного мониторинга физического состояния зданий и/или сооружений и устройство для его осуществления | |
RU2766130C1 (ru) | Способ диагностики подшипников роторного оборудования на основе оценки микровариаций вращения вала | |
RU2684709C1 (ru) | Способ акустико-эмиссионной диагностики динамического промышленного оборудования | |
Kumar et al. | Adaptive wavelet based signal processing scheme for detecting localized defects in rolling element of taper roller bearing | |
RU2547504C1 (ru) | Способ выявления наличия дефектов узлов и агрегатов автомобиля в реальном времени и устройство для его осуществления | |
Sheen et al. | A quantified index for bearing vibration analysis based on the resonance modes of mechanical system | |
Yan et al. | Research on fault diagnosis method of bogie bearing based on microphone arrays | |
CN117073900A (zh) | 一种电动机转子动平衡监测系统 | |
Xu et al. | Windowed Correlated Kurtosis for Isolation and Identification of Compound Faults in Locomotive Bearings | |
JPS6144016B2 (es) |