Análisis de
vibraciones
en maquinaria
Edwin Espinel Blanco
Andrés Camilo Melgarejo Ramos
Hugo Alberto Ramírez Galvis
A N Á L I S I S D E V I B R AC I O N E S
E N M AQ U I N A R I A
E D W I N E S P I N E L B L A N CO
A N D R É S C A M I LO M E LG A R E J O R A M O S
H U G O A L B E R TO R A M Í R E Z G A LV I S
U N I V E R S I D A D F R A N C I S CO
D E PAU L A S A N TA N D E R O C A Ñ A
FA C U LTA D D E I N G E N I E R Í A
PLAN DE ESTUDIOS DE
INGENIERÍA MECÁNICA
O C A Ñ A , CO LO M B I A
Catalogación en la publicación - Biblioteca Nacional de Colombia
Espinel Blanco, Edwin
Análisis de vibraciones en maquinaria / Edwin Espinel Blanco, Andrés Camilo Mel-
garejo Ramos, Hugo Alberto Ramírez Galvis. -- 1a. ed. -- Bogotá : Ecoe Ediciones : Universidad
Francisco de Paula Santander, 2019.
61 p.
Contiene: complemento virtual SIL (Sistema de Información en Línea) [Link].
-- “Incluye: Teoría de las vibraciones mecánicas aplicadas en el diagnóstico de maquinaria rotativa.
Introducción en el manejo de equipos utilizados en la toma y medición de vibraciones. Formas
de realizar la colección de datos explicando los tipos de espectros asociados a cada falla”. -- Incluye
datos biográficos de los autores. -- Contiene bibliografía.
ISBN 978-958-771-690-0
1. Vibración 2. Vibración - Mediciones 3. Mecánica I. Melgarejo Ramos, An-
drés Camilo II. Ramírez Galvis, Hugo Alberto III. Título IV. Serie
CDD: 620.3 ed. 23 CO-BoBN– a1038818
Colección: Ingeniería y salud en el trabajo
Área: Ingeniería
▶ Edwin Espinel Blanco Primera edición: Bogotá, diciembre de 2018
▶ Andrés Camilo Melgarejo Ramos
▶ Hugo Alberto Ramírez Galvis ISBN: 978-958-771-690-0
© Ecoe Ediciones Limitada Coordinación editorial: Angélica García Reyes
[Link] Corrección de estilo: Osmar Peña
Carrera 19 # 63C 32 Diagramación: Denise Rodríguez
Bogotá, Colombia Carátula: Wilson Marulanda Muñoz
Impresión: La Imprenta Editores
© Universidad Francisco Calle 77 # 27 A - 39
de Paula Santander
Vía Acolsure, Sede el Algodonal
Ocaña Norte de Santander -
Colombia
Teléfono (057)(7) 5690088
Prohibida la reproducción total o parcial por cualquier medio
sin la autorización escrita del titular de los derechos patrimoniales.
Impreso y hecho en Colombia - Todos los derechos reservados
CO N T E N I D O
Introducción.............................................................................................. IX
Capítulo 1: Modelos matemáticos y contexto
de las vibraciones..................................................................................... 1
1.1. Naturaleza de las vibraciones......................................................................1
1.2. Modelos matemáticos de las vibraciones...................................................3
1.2.1. Vibración libre sin amortiguamiento...............................................3
1.2.2. Vibración libre con amortiguamiento.............................................5
1.2.3. Vibración forzada sin amortiguamiento.........................................8
1.2.4. Vibración forzada con amortiguamiento........................................10
1.3. Caracterización del movimiento vibratorio..............................................13
1.3.1. Unidades de medición de la vibración............................................16
1.4. Límites de severidad en la vibración..........................................................18
1.4.1. Carta de Rathbone..............................................................................18
1.5. Norma ISO 2372...........................................................................................20
1.6. Norma ISO 10816.........................................................................................21
Capítulo 2: Toma de datos de vibraciones.......................................... 25
2.1. Pasos previos a la medición........................................................................25
2.2. Puntos de medición......................................................................................26
2.3. Problemas al inicio de la medición............................................................27
2.4. Proceso de medición....................................................................................29
VI Análisis de vibraciones en maquinaria
Capítulo 3: Análisis de vibraciones en maquinaria rotativa........ 31
3.1. ¿Por qué es importante un análisis de vibraciones?.................................31
3.2. El espectro de vibraciones...........................................................................32
3.3. Identificar la frecuencia de giro..................................................................34
3.3.1. Fase.......................................................................................................35
3.3.2. Armónicos...........................................................................................35
3.3.3. Bandas laterales...................................................................................35
3.3.4. Ruido de piso......................................................................................35
3.4. Comparación con medidas de referencia..................................................36
3.5. Condiciones de falla en máquinas rotativas.............................................38
3.5.1. Desbalance...........................................................................................38
3.5.2. Desalineamiento.................................................................................41
3.5.3. Soltura estructural..............................................................................42
3.5.4. Falla en rodamientos..........................................................................43
3.5.5. Caja de engranajes..............................................................................46
3.5.6. Recomendaciones básicas para el analista de vibraciones............46
Capítulo 4: Casos de estudio: expansión fourier.............................. 47
4.1. Observaciones...............................................................................................52
4.2. Expansión de Fourier para describir el espectro de un motor
desbalanceado...............................................................................................53
Anexo ............................................................................................................ 59
Bibliografía................................................................................................ 61
Contenido VII
ÍNDICE DE FIGURAS
Figura 1. Oscilación de una masa a un resorte....................................................2
Figura 2. Trayectoria sinusoidal............................................................................2
Figura 3. Sistema libre amortiguado.....................................................................5
Figura 4. Diagrama de cuerpo libre del sistema..................................................6
Figura 5. Sistema con vibración forzada amortiguada.......................................11
Figura 6. Comportamiento de la velocidad.........................................................14
Figura 7. Comportamiento de la aceleración......................................................15
Figura 8. Comparación de las fases.......................................................................15
Figura 9. Descripción de la escala RMS...............................................................16
Figura 10. Descripción del movimiento periódico...............................................17
Figura 11. Carta Rathbone........................................................................................19
Figura 12. Puntos de medición en motor...............................................................26
Figura 13. Medición defectuosa con efecto ski......................................................28
Figura 14. Posicionamiento correcto para evitar impacto...................................28
Figura 15. Representación de una función periódica en el dominio
de la frecuencia........................................................................................33
Figura 16. Espectro de vibración.............................................................................34
Figura 17. Identificación de la frecuencia de giro.................................................34
Figura 18. Gráficos de tendencia.............................................................................36
Figura 19. Espectros desplegados apilados............................................................37
Figura 20. Visualización en cascada........................................................................38
Figura 21. Representación del espectro por desbalanceo estático,
realizada por equipo a-predictor..........................................................39
Figura 22. Representación del espectro de falla por desbalanceo
dinámico..................................................................................................41
Figura 23. Representación del espectro de falla por soltura estructural............43
Figura 24. Representación del espectro de falla en rodamiento con fisura.......45
Figura 25. Señal de vibración en el tiempo............................................................47
Figura 26. Fracción de la señal de vibración correspondiente
a un período de 0,1 s..............................................................................48
Figura 27. Representación de la onda en la expansión de Fourier
para n 5 3................................................................................................51
Figura 28. Representación de la onda en la expansión de Fourier
para n 5 6................................................................................................51
Figura 29. Representación de la onda en la expansión de Fourier
para n 5 10..............................................................................................52
Figura 30. Gráfico armónicos an del dominio de la frecuencia...........................57
Figura 31. Gráfico armónicos bn del dominio de la frecuencia...........................58
VIII Análisis de vibraciones en maquinaria
Í N D I C E D E I M ÁG E N E S
Imagen 1. Verificar el buen estado del equipo.........................................................27
Imagen 2. Equipo de pruebas para análisis de vibración
con disco desbalanceado...........................................................................39
Imagen 3. Equipo de pruebas y análisis con desbalanceo dinámico....................40
Imagen 4. Equipo de mediciones para prueba de vibración..................................59
Imagen 5. Simulación de falla por desbalanceo dinámico.....................................59
Imagen 6. Simulación de falla por desbalanceo estático........................................60
I N T R O D U CC I Ó N
Las vibraciones mecánicas son en general un fenómeno complejo y su estudio
busca describir el comportamiento de cuerpos sujetos a una excitación como
puede ser una fuerza. El movimiento o los movimientos que resultan como conse
cuencia de esta fuerza corresponden al análisis de vibraciones, gracias al cual se
descifran las características de dicha fuerza que ocasionó el movimiento o excita-
ción, así como la cantidad de energía liberada durante el mismo, que bien puede
ser un desplazamiento o una velocidad resultante del cuerpo afectado.
Cuando observamos todo lo que nos rodea, por ejemplo, un automóvil que es
sacudido al pasar sobre un camino destapado de manera que empieza a oscilar,
este movimiento de vaivén es conocido como vibración y su estudio es muy im-
portante para entender cómo el automóvil puede resistir (de mejor manera) las
fuerzas generadas sobre sus neumáticos y transmitidas a su carrocería por la sus
pensión, mediante el uso de materiales adecuados que resistan las deformaciones
internas. Del mismo modo, el análisis de vibraciones se utiliza para identificar fa
llas en maquinarias rotativas, tuberías y en diagnósticos de enfermedades cardia
cas, gracias al espectro de vibración que relaciona los movimientos de distintas
velocidades presentes en un cuerpo y su intensidad o severidad. A lo largo de este
trabajo se desarrollará la teoría básica que permite entender y modelar el compor-
tamiento de estos sistemas vibratorios.
X Análisis de vibraciones en maquinaria
El análisis de vibraciones es una importante herramienta utilizada para el diseño,
el diagnóstico y la evaluación del comportamiento de casi cualquier sistema me-
cánico, por eso hace parte de la formación del ingeniero mecánico y es fundamen-
tal que este pueda dominar estas teorías para su aplicación.
Este documento presenta de forma concreta los fundamentos de vibraciones y su
aplicación en análisis de máquinas rotativas. Se busca informar a los estudiantes
de Ingeniería mecánica sobre el empleo del equipo de vibraciones y la interpreta-
ción de los datos obtenidos desde las mediciones para que puedan aprovechar los
conocimientos para diagnósticos de fallas en maquinaria rotativa, a través de la
comprensión de los conceptos básicos sobre vibraciones, el equipo de mediciones
A-MAQ, las mediciones y la aplicación de procesos matemáticos más complejos
que son de aplicación generalizada y el análisis de vibraciones de todo tipo de es
tructuras, pero que no se cubren en el curso de vibraciones.
C APÍTULO 1
M O D E LO S M AT E M ÁT I CO S
Y CO N T E X TO D E L A S V I B R AC I O N E S
1.1. Naturaleza de las vibraciones
La vibración se define como un movimiento que genera oscilaciones o perturbacio
nes en una posición de equilibrio de un cuerpo en cualquier dirección, estas son
causadas por fuerzas externas o internas y en general se repiten después de un
intervalo de tiempo, como en los casos de maquinaria rotativa. Debido a que los
niveles altos de vibración evidencian la exposición de un equipo a una fuerza des-
controlada que ocasionará un fallo en cualquier momento, es importante conocer
su naturaleza.
La forma más simple de vibración es el movimiento armónico simple y el modelo
fundamental del análisis de vibraciones se obtiene cuando una masa “M”, unida
a un resorte que está fijo a una pared, se suelta desde una posición x1 respecto a
una posición de equilibrio x0, hasta que llega a una posición x2 (figura 1). Este per
manecería en ese movimiento de vaivén indefinidamente, pero eventualmente la
fricción hará que retorne al punto de equilibrio.
2 Análisis de vibraciones en maquinaria
Figura 1. Oscilación de una masa a un resorte.
x2
M
x0
M
x1 x
M
Fuente: Autores, 2018.
Si observamos el movimiento de la masa en el tiempo, notamos que esta oscila libre
mente alrededor del punto de equilibrio y el comportamiento de su trayectoria en
el tiempo es sinusoidal (figura 2). En las máquinas, la mayoría de sus componentes
individuales vibran de esta manera, aunque la vibración de toda la máquina es una
superposición de múltiples movimientos, por lo que, en general, el movimiento
resultante no es sinusoidal.
Figura 2. Trayectoria sinusoidal.
M
x0
Tiempo
Fuente: Autores, 2018.
Capítulo 1: Modelos matemáticos y contexto de las vibraciones 3
1.2. Modelos matemáticos de las vibraciones
El modelo matemático consiste en hacer idealizaciones para destacar aquellos fac-
tores que intervienen en el proceso estudiado y en el que se extraen aquellos pará
metros que influyen decisivamente y se le atribuyen propiedades características, al
mismo tiempo que se descartan los parámetros poco determinantes. Un sistema
vibratorio se compone de un medio almacenador de energía potencial (como un
resorte), un medio para almacenar energía cinética o masa y un medio para disipar
la energía o amortiguador. A continuación, se hará un resumen de los modelos
fundamentales.
Podemos describir dos tipos de modelos: los correspondientes a vibraciones li-
bres que se dan cuando no existen fuerzas externas aplicadas sobre algún sistema
vibratorio estudiado en un tiempo específico y las vibraciones forzadas, las cuales
aparecen cuando hay fuerzas externas aplicadas al sistema durante el transcurso
del tiempo en cuestión, y en momentos internos (Roca Vila y León, 1981).
A su vez, los modelos de las vibraciones forzadas y de las no forzadas podrían o no
ser afectados por fuerzas que amortigüen su movimiento y, con base en esto, divi-
dirse en:
»» Sin amortiguamiento: no existe fuerza o efecto que absorba la energía
del movimiento del sistema.
»» Con amortiguamiento: hay fuerzas o momentos que disipan el movi-
miento vibratorio del sistema.
1.2.1. Vibración libre sin amortiguamiento
Al considerar el sistema vibratorio más sencillo descrito en la figura 1, este está
formado por un resorte lineal (K) y una masa (m). Como sabemos, el resorte ge-
nera una fuerza del tipo:
F 5 2K ? x
Al aplicar la ley de Newton en la dirección vertical:
2K ? x 5 m ? x
Es decir,
m?x1K?x50 (1)
Esta ecuación representa un movimiento simple no amortiguado y su solución
es la correspondiente a una ecuación diferencial lineal de coeficientes constantes,
4 Análisis de vibraciones en maquinaria
por lo tanto, este modelo aplica como solución de todos los modelos de vibración
lineal.
La ecuación característica de (1) es:
K
s2 1 50
m
Al ser sus raíces imaginarias conjugadas:
K
s56 i
m
La solución general de (1) está dada por:
x(t) 5 A1 cos vn ? t 5 A2 sin vn ? t (2)
Esta ecuación se puede escribir de forma más compacta si se considera:
x 5 A2
( A1
A2
cos vn ? t 1 sin vn ? t
)
Al introducir el ángulo definido por:
A1
tan w 5
A2
A2
x5 (sin w ? cos vn t 1 sin vnt ? cos w)
cos w
Al aplicar la relación de la suma de ángulos:
x(t) 5 A12 1 A22 ? sin(vnt 1 w)
Se obtiene:
x(t) 5 A sin(vnt 1 w) (3)
La ecuación 3 representa la vibración de un sistema libre no amortiguado, co-
múnmente llamada ecuación del movimiento armónico simple.
K
Donde A 5 A12 1 A22 es la amplitud del movimiento y vn 5 es la frecuencia
m
natural del sistema.
Capítulo 1: Modelos matemáticos y contexto de las vibraciones 5
Para un sistema en movimiento oscilatorio la frecuencia natural será:
vn 1 K
fn 5 5 ?
2p 2pp m
Y el período o duración será:
1 m
T5 5 2p
fn K
1.2.2. Vibración libre con amortiguamiento
Si ahora imaginamos un sistema mecánico que consta de un mecanismo de amorti
guamiento (figura 3), se tiene en cuenta que un mecanismo real presenta siempre
una amortiguación de su movimiento, ya sea debido a la fricción o porque existe
un elemento destinado para tal fin, en todo caso, siempre deberíamos tener esto
en cuenta.
Figura 3. Sistema libre amortiguado.
K C
m
x
Fuente: Autores, 2018.
Donde:
»» m 5 masa
»» c 5 coeficiente de amortiguación
»» x 5 desplazamiento
»» K 5 constante de rigidez elástica
6 Análisis de vibraciones en maquinaria
El principio de funcionamiento de la amortiguación lineal es la estrangulación de
un flujo y es lineal porque el efecto de las fuerzas cortantes generadas en el flujo de
fluidos newtonianos es proporcional a la velocidad. Así que tenemos:
F 5 2cx
La constante c depende las propiedades del fluido y de la geometría del sistema.
Al aplicar la ecuación de equilibrio al sistema se tiene que 2c ? x 2 K ? x 5 m ? x
(figura 4).
Figura 4. Diagrama de cuerpo libre del sistema.
Kx
cx
Fuente: Autores, 2018.
Al reacomodar:
m?x1c?x1K?x50 (4)
La ecuación diferencial (4) describe el movimiento libre amortiguado y su ecua-
ción característica es:
ms2 1 cs 1 k 5 0
Cuyas raíces se pueden escribir como:
c K c2
s1,2 5 2 56 ? 21
2m m 4mK
Donde:
c
z5
2 Km
Capítulo 1: Modelos matemáticos y contexto de las vibraciones 7
z se llama al factor de amortiguación.
K
Si introducimos nuevamente vn 5 , se tiene que:
m
s1,2 5 2zvn 6 vn ? z2 21
La solución formal de (4) es:
x(t ) 5 A1e s1t 1 A2 e s2t (5)
Esta ecuación puede tomar varias formas según si los valores de s1,2 son reales o
complejos, lo cual corresponde físicamente a una amortiguación grande (z . 1) o
a una amortiguación débil (z , 1), aunque también está el caso extremo de z 5 1,
correspondiente a la amortiguación crítica.
Por lo que se presentan los siguientes casos posibles:
»» Si z . 1 el movimiento será sobre amortiguado.
Entonces las raíces s1,2 son reales y distintas a la solución de esta ecuación con amor
tiguamiento y será:
x(t) 5 A1e s1t 1 A2e s2t
Donde A1 y A2 son constantes que surgen de integrar.
El sistema permanece sin movimiento y vuelve al punto de equilibrio, si el amorti
guamiento se incrementa, el sistema tardará más tiempo para regresar al equilibrio.
»» Si z , 1 el movimiento es sub-amortiguado, las raíces s1,2 son imaginarias
conjugadas e iguales a:
s1,2 5 2zvn 6 ivd
Siendo:
vd 5 vn 1 2 z2
O bien llamada: frecuencia amortiguada.
Al tener en cuenta estas anotaciones, la solución de 5 se resuelve como:
x 5 e2zvnt(A1 cos vd ? t 1 A2 sin vd ? t)
8 Análisis de vibraciones en maquinaria
Al reescribir:
A ? e 2zvnt ? sin(vdt 1 w) (6)
Donde:
2 2
A 5 A1 1 A2
A1
tan w 5
A2
En esta solución, el movimiento tiene características periódicas y los intervalos
establecidos por t9, correspondientes a una vibración amortiguada, se expresan
en función del período t, así:
2p
t9 5
vd
»» Si z 5 1, la amortiguación es crítica.
c 5 2 Km 5 cc
Las raíces s1,2 son iguales y están definidas por:
s1,2 5 2vn
En este caso, la solución de la ecuación diferencial es del tipo:
x(t) 5 (A1 1 A2t)e2vnt
Es decir, no contiene funciones periódicas, por lo que este sistema retorna al equi-
librio sin oscilar y en el menor tiempo. Este caso representa la menor cantidad de
amortiguamiento para que el sistema no oscile, por ello resulta útil su análisis en
aplicaciones donde se necesita evitar la vibración y que el sistema alcance el equi
librio rápidamente.
Al volver a la ecuación del factor de amortiguación z, se puede observar que para
una amortiguación c cualquiera se puede escribir como:
c c
z5 5
2 Km cc
1.2.3. Vibración forzada sin amortiguamiento
Para lograr que un sistema oscile se debe suministrar energía a este. Cuando se
realiza esta operación de manera continua, como puede ser a través de una fuerza,
Capítulo 1: Modelos matemáticos y contexto de las vibraciones 9
entonces la vibración es forzada y cuando se agrega energía en el sistema con una
tasa mayor de la disipada, esta aumenta con el tiempo, lo que se manifiesta por un
incremento de la amplitud del movimiento (Crawford, 1977). Cuando se sumi
nistra energía a la misma tasa de disipación, la amplitud permanece constante
en el tiempo y si la fuerza que potencia el movimiento es periódica, la ecuación
diferencial del movimiento estará definida así:
m ? x 1 K ? x 5 F 5 F0 sin vt (7)
Donde F0 es la amplitud y v es la frecuencia de la fuerza excitadora.
La solución general de la ecuación diferencial se obtiene al añadir a la solución ge
neral de la homogénea una solución particular x 5 xh 1 xp.
La ecuación característica será:
ms2 1 k 5 0
K
Las raíces de esta ecuación son imaginarias y conjugadas: s 5 6 i
m
La solución general de la homogénea es:
xh 5 A sin(vnt 1 w)
La solución particular es:
xp 5 A cos(vnt)
Y la solución general se expresa como:
F0
x 5 A cos(vnt 1 w) 1 K 2 cos vt (8)
v
12 2
vn
Para los sistemas no amortiguados y forzados armónicamente, su movimiento es
la suma de dos armónicos: uno a frecuencia natural vn y otro a frecuencia de la
fuerza exterior v.
La amplitud del primero depende de las condiciones iniciales y se anula para algu-
nos valores particulares, la amplitud del segundo se establece con la proximidad
de ambas frecuencias a través de la expresión llamada factor de resonancia:
1 A
r5 2 5
v x est
12 2
vn
10 Análisis de vibraciones en maquinaria
Batimiento: se produce cuando la frecuencia natural del sistema vn se aproxima
a la frecuencia exterior v, es decir, para el caso en que:
vn 5 v 1 Dv. Donde sin perturbación inicial (x0 5 x 09 5 0) se obtiene:
x5
F0vn
kDv
? sin
(
Dv
2
? t sin vt
)
Esta ecuación representa un movimiento con frecuencia vn y amplitud armónica,
la cual crece hasta un máximo y disminuye hasta que se anula, este ciclo se repite
de forma periódica (Beer & Russell, 2010).
Resonancia: cuando una fuerza excitadora de las vibraciones tiene frecuencias
coincidentes con el sistema original, esta provoca alteración de la configuración
del sistema mecánico al alcanzar amplitudes que generalmente ocasionan un fallo
estructural del material sometido a esfuerzos, incluso con fuerzas iniciadoras de
muy bajo nivel.
Cuando la frecuencia de la fuerza exterior es igual a la frecuencia natural del sis-
tema (v 5 vn), es decir, cuando Dv → 0, se produce la resonancia y la ecuación
que rige este fenómeno es:
F0v
x5 ? t sin(vn t)
2K
Expresión que corresponde a un movimiento armónico y cuya amplitud tiende a
infinito cuando t → `.
1.2.4. Vibración forzada con amortiguamiento
Analizaremos el caso de la vibración forzada al considerar el sistema masa-amor-
tiguador-resorte, sometido a una fuerza exterior F(t), (figura 5).
Se ha visto que la energía mecánica de un oscilador amortiguado disminuye en el
tiempo como resultado de esa fuerza resistiva y es posible compensar esta disminu
ción de energía al aplicar una fuerza externa que haga trabajo positivo al sistema.
Por ejemplo, un niño en un columpio se puede mantener en movimiento mediante
“empujones” adecuadamente cronometrados. La amplitud del movimiento perma-
nece constante si la entrada de energía por cada ciclo iguala exactamente a la dis-
minución de energía mecánica en cada ciclo, la cual resulta de las fuerzas resistivas
(Serway & Jewett, 2008).
Capítulo 1: Modelos matemáticos y contexto de las vibraciones 11
Figura 5. Sistema con vibración forzada amortiguada.
K C
F(t ) m
Fuente: Autores, 2018.
Se consideran las siguientes hipótesis para el sistema descrito:
a) La masa se mueve únicamente en dirección vertical y ausente de fricciones
y movimiento giratorio.
b) La masa del resorte es despreciable, a la vez que su fuerza recuperadora es
proporcional a su deformación.
c) El amortiguador tiene un comportamiento de rozamiento de tipo viscoso.
d) Los elementos actúan en vacío.
Se consideran F, la fuerza aplicada sobre el sistema yc ? x, la fuerza amortiguadora
de tipo viscoso. Con las condiciones de: m . 0; c . 0; K . 0, la ecuación de equili
brio será:
m ? x 1 c ? x 1 K ? x 5 F(t)
En general, esta fuerza es una función cualquiera del tiempo, de tipo periódico,
siendo el caso más sencillo F(t) 5 F0 cos vt, donde nuevamente F0 se refiere a la
amplitud de la fuerza aplicada y v es su frecuencia. Para tal caso, la ecuación se
puede escribir:
m ? x 1 c ? x 1 K ? x 5 F0 cos vt (9)
12 Análisis de vibraciones en maquinaria
Donde la ecuación característica de (9) es:
ms2 1 cs 1 k 5 0
Como sabemos, la solución de esta ecuación se puede descomponer en:
x 5 xh 1 xp
La solución de la ecuación homogénea es:
m?x1c?x1K?x50
Se toma la forma previamente estudiada y descrita en la ecuación (6), donde he-
mos supuesto un amortiguamiento sub-crítico para que resulte una vibración y
como todas las soluciones de esta ecuación tienden a cero cuando crece t, esta
representa la parte transitoria del movimiento vibratorio:
xh 5 A ? e2zvnt ? sin(vdt 1 w)
Pasemos a determinar la solución particular. Al proceder por el método de identi
ficación de constantes, obtenemos la solución:
xp 5 B cos(vt 2 w)
Donde B es:
F0
B5
( K 2 mv ) 1 (cv)
2 2 2
Y consideramos el ángulo:
cv
tan w 5
K 2 mv 2
Por lo que la solución de la ecuación 8 es:
2zvn t
x5A? e ? sin(vdt 1 w) 1 xp 5 B cos(vt 2 w) (10)
La primera parte generalmente desaparece al cabo de un corto tiempo. La se-
gunda parte es la componente estacionaria que oscila con amplitud B constante y
desfase w.
Capítulo 1: Modelos matemáticos y contexto de las vibraciones 13
1.3. Caracterización del movimiento vibratorio
El movimiento vibratorio cuenta con características que hacen posible determinar
su origen. Para lograr un manejo práctico del fenómeno, utilizaremos el modelo
matemático descrito en la ecuación 3, el cual representa el movimiento armónico
simple (Serway & Jewett, 2008).
El desplazamiento de un movimiento armónico simple (figura 1) está determina
do por la ecuación:
x(t) 5 A sin(v ? t 1 w) ec. 3
Donde
»» x(t) 5 Desplazamiento
»» A 5 Amplitud
»» v 5 Frecuencia
»» t 5 Tiempo o período
»» w 5 Fase
El desplazamiento x(t) es la cantidad de movimiento que se produce en un ciclo,
por ser una vibración que presenta mayor grado de amplitud a bajas frecuencias,
se pueden identificar un rango de anomalías muy limitado en la máquina.
La amplitud (A) hace referencia a la intensidad del movimiento o a la energía
involucrada en este, por lo que entre mayor sea la amplitud más severo será el
efecto del movimiento.
El período (t) del movimiento es el tiempo que tarda la masa en ir al punto x2
y volver hasta el punto x1 (ilustrados en la figura 1) y representa la duración del
movimiento cíclico del vaivén.
La frecuencia (v) es la velocidad a la que se realiza el movimiento cíclico y es el
inverso del período, es decir, la frecuencia es la cantidad de ciclos en un período
de tiempo.
La fase (w) es la medida de referencia para identificar la diferencia entre dos se-
ñales sinusoidales en el tiempo, por lo que la fase relaciona el período entre dos
señales, aunque en realidad esta es una medida de tiempo y usualmente se repre-
senta en grados, esto se hace al considerar que un ciclo consta de 360°, de manera
que la diferencia entre dos señales se mide por el retraso o adelanto de una onda
respecto al ciclo de otra. Por ejemplo, un retraso en la fase de 180° en una señal
significa que esta apenas ha completado la mitad de su ciclo cuando la señal de
referencia tenga un ciclo completo.
14 Análisis de vibraciones en maquinaria
Estos comportamientos se pueden apreciar en la onda si se comparan la velocidad
y la aceleración con desplazamiento y se identifica que estas están retrasadas o
desfasadas en 90° y 180° respecto a este (figuras 5 y 6). Estos conceptos se pueden
observar como se muestra en la figura 7.
Al derivar una y dos veces la ecuación 3, hallaremos la velocidad y la aceleración:
v(t) 5 Av cos(v ? t 1 w)
a(t) 5 2Av2 sin(v ? t 1 w)
Los valores de velocidad y de aceleración son comúnmente usados en vez del valor
de desplazamiento para medir las vibraciones de la maquinaria rotativa. En gene-
ral, el valor de velocidad ofrece información muy buena a casi todos los niveles de
frecuencia, esto es debido a que la velocidad es un claro indicador de la severidad
de la vibración por ser directamente proporcional a la frecuencia, por otro lado,
la aceleración será especialmente útil a alta frecuencia porque al ser proporcional
a la frecuencia al cuadrado se harán más notorias las vibraciones ubicadas en ese
rango, de manera que se facilita el encontrar todos los elementos que vibran en
dicho nivel, como los de rodamientos cuando tienen una falla de tipo metal-metal.
Por último, la gráfica obtenida a partir del desplazamiento será importante a muy
bajas frecuencias, donde el movimiento es amplio y fácilmente identificable.
Figura 6. Comportamiento de la velocidad.
Desplazamiento - velocidad
Desplazamiento
Velocidad máxima
Tiempo
Velocidad mínima
Fuente: Autores, 2018.
Capítulo 1: Modelos matemáticos y contexto de las vibraciones 15
Figura 7. Comportamiento de la aceleración.
Desplazamiento - aceleración
Aceleración mínima
Tiempo
Aceleración máxima
Fuente: Autores, 2018.
Si superponemos las gráficas anteriores y se comparan con el desplazamiento, no-
tamos que la velocidad se encuentra retrasada o desfasada por 90°, respecto al
desplazamiento, y la aceleración también desfasada, pero en 180°.
Figura 8. Comparación de las fases.
Amplitud Velocidad
Aceleración
Desplazamiento
180˚
90˚ Fase
Fuente: Autores, 2018.
16 Análisis de vibraciones en maquinaria
1.3.1. Unidades de medición de la vibración
La amplitud (A) es la medida de la intensidad del movimiento, comúnmente se
toma en tres escalas: pico-pico, cero-pico y RMS. Es importante conocer la escala
en la que se toma la medición, ello con el fin de compararla con niveles estableci-
dos en la norma, ya que estos aparecen en la escala RMS, además de que mantener
una escala única para cada punto de medición permite el monitoreo de un mismo
rango de frecuencia, lo cual hace posible un análisis preciso.
Pico-pico (pk-pk): representa la distancia desde la cresta inferior hasta la cresta
superior de la onda.
Cero-pico (pk): es la distancia máxima medida desde el punto de equilibrio hasta
la cresta de la onda, si la onda es sinusoidal cero-pico es la mitad de pk-pk.
RMS (Root Mean Square): es la raíz de los promedios de amplitud elevados al
cuadrado. Se considera esta escala de medida muy útil porque representa el pro-
medio de la intensidad global de la vibración. El procedimiento consiste en elevar
al cuadrado los picos de amplitud positivos y negativos, sumarlos y dividirlos entre
el número total, de forma que el área entre las crestas sea igual al área debajo de
cada cresta (A1 5 A1’ 5 A2 5 A2’), después se saca la raíz del promedio obtenido.
Figura 9. Descripción de la escala RMS.
RMS A1’
A1 A2
A2’
Fuente: Autores, 2018.
La frecuencia (f ): las veces que se realiza un movimiento en un determinado tiem
po. Sus unidades son:
Hz (Hertz): ciclos por segundo
CPM: ciclos por minuto
RPM: revoluciones por minuto
Capítulo 1: Modelos matemáticos y contexto de las vibraciones 17
El período (t): es el tiempo que tarda en completarse un movimiento. Se mide en
segundos (s) o milisegundos (ms) y es el inverso de la frecuencia. t 5 1/f.
El desplazamiento (D): representa el cambio de posición respecto a un punto de
referencia, se mide en mm o mils (milésimas de pulgada) y comúnmente se utiliza
la escala pk-pk para obtener la distancia total recorrida por un eje.
La velocidad (V): es la rapidez con que varía la posición de un punto a otro, se mi
de en pulg/s (pulgadas por segundo) o mm/s (milímetros por segundo) y se utiliza
la escala 0-pk o RMS.
La aceleración (a): es la variación de la velocidad, se mide en mm/s2, pulg/s2 o co-
múnmente respecto a la aceleración de la gravedad G (1G 5 386 pulg/s2 5 9.800
mm/s2). Se utiliza la escala RMS, dado que es una medida bastante útil a grandes
velocidades . 10.000 RPM.
Como ilustra la figura 10, dependiendo de la característica del movimiento que
se mida (desplazamiento, velocidad o aceleración) será más adecuada una escala
u otra.
Figura 10. Descripción del movimiento periódico.
Parámetros descriptores de un movimiento armónico simple
T 5 1/f (1 ciclo)
A0-pico RMS
Apico-pico
Promedio
T: Período 5 [segundos] 5 [s] f: Frecuencia 5 [ciclos/segundos] 5 [Hz]
Fuente: Autores, 2018.
18 Análisis de vibraciones en maquinaria
Conversión entre escalas
Sistema internacional
D 5 diámetro (mm)
V 5 velocidad (mm/s)
a 5 aceleración (mm/s2)
g 5 gravedad (9,8 m/s2)
DPK 2PK 1.000 gaPK
DRMS 5 DPK2PK 5
2 ( pf )
2
2? 2
VPK
VRMS 5 V 5 p ? f ? DPK2PK
2 PK
aPK 2pf ? VPK
aRMS 5 aPK 5
2 1.000 g
1.000 ? g ? aPK
VPK 5
2pf
1.4. Límites de severidad en la vibración
1.4.1. Carta de Rathbone
Esta constituyó el primer establecimiento de límites de vibración aceptados y usa-
dos en el sector industrial, la cual fue realizada en los años treinta y posteriormente
mejorada. La carta está compuesta por una escala logarítmica de frecuencias (her-
cios o RPM) y de una logarítmica de amplitudes en desplazamiento (Pico o Pico-
Pico) y en velocidad, las cuales permitían establecer un límite de severidad de
forma rápida al ubicar las características de la máquina analizada en ocho zonas.
Capítulo 1: Modelos matemáticos y contexto de las vibraciones 19
Figura 11. Carta Rathbone.
100.000
150.000
1.000
200
400
600
800
RPM
1.000
100
Desplazamiento pico a pico (micrones)
16
10
8
1 4
Velocidad mm/s
0.8
0.6
0.4 2
0.2 1
0.1
0.1 0.2 0.5
25 5
Fuente: Rathbone, 1930.
La carta presenta las siguientes limitaciones:
»» No se consideran las características de la máquina como rigidez y poten-
cia ni tampoco su tipo.
»» Esta guía aplica solo a equipo rotativo y no a los alternativos u otros.
»» Si la frecuencia es más alta, la amplitud de vibración en desplazamiento
tiene que ser menor para mantener el mismo nivel de severidad, por ejem-
plo, cuando una máquina vibra a 200 CPM con 200 micras de amplitud
pk-pk, este es un nivel “bueno” de severidad, pero cuando esta amplitud
corresponde a una frecuencia de 2.000 CPM, el nivel de severidad sería
“grave”. La vibración a baja frecuencia es menos peligrosa que la vibración
a alta frecuencia, de ahí que las averías de engranajes y rodamientos que
se producen generalmente a alta frecuencia sean muy peligrosas. Por esta
razón, las amplitudes de baja frecuencia se miden en desplazamientos y
las de alta frecuencia en velocidad o aceleración. La carta de Rathbone fue
creada para máquinas de bajas RPM y hoy está fuera de uso.
20 Análisis de vibraciones en maquinaria
1.5. Norma ISO 2372
Vibración mecánica de máquinas con velocidades de operación
entre 100 y 200 rev/s
Consideraciones para aplicar ISO 2372:
»» Se aplican a equipos con velocidades de giro entre 600 y 12.000 RPM.
»» Los niveles de vibración expresan el nivel global de vibración en velocidad
RMS para frecuencias en un rango de 10 y 1.000 Hz. Inicialmente, los
equipos se clasifican mediante la tabla 1, luego para determinar el nivel
de vibración en ISO 2372 se debe ubicar la máquina dentro de la clase
correspondiente, después colocar el valor global de vibración obtenido
(en la tabla 2 aparecen identificados para frecuencias entre 600 y 60.000
CPM) y finalmente identificar la zona.
Tabla 1. Clasificación de equipos.
Clase Descripción
Clase I Equipos pequeños de hasta 15 Kw
Clase II Equipos medios, de 15 a 75 kW o hasta 300 kW con cimentación especial
Equipos grandes, por encima de 75 kW con cimentación rígida o de 300
Clase III
kW con cimentación especial
Clase IV Turbomaquinaria (equipos con RPM . velocidad crítica)
Fuente: ISO 2372, 1974.
Capítulo 1: Modelos matemáticos y contexto de las vibraciones 21
Tabla 2. Severidad de la vibración.
Velocidad Tipos de máquinas
(mm/s, RMS) Clase I Clase II Clase III Clase IV
0,18 a 0,28
0,28 a 0,45 A
0,45 a 0,71
0,71 a 1,12
1,12 a 1,8 B
1,8 a 2,8
2,8 a 4,5 C
4,5 a 7,1
7,1 a 11,2 D
11,2 a 18
18 a 28
Buena A Insatisfactoria C
Satisfactoria B Inaceptable D
Fuente: ISO 2372, 1974.
1.6. Norma ISO 10816
La norma ISO 10816 remplaza desde 1995 a la ISO 2372. Actualmente esta es la
normativa vigente para establecer niveles de vibración y funciona al realizar una
clasificación previa de las máquinas de acuerdo a sus potencias, apoyos, velocida-
des de operación y aplicaciones. Para complementar esta visión, se ilustra una tabla
con los niveles de vibración globales dependiendo del tipo y de las características
de la maquinaria, esto resulta práctico para una evaluación del estado de la máqui
na como se muestra en la tabla 3. La ISO 10816 consta de las siguientes partes, bajo
el título general de vibración mecánica y de la evaluación de la vibración de la má
quina mediante mediciones en piezas no giratorias:
22 Análisis de vibraciones en maquinaria
Parte 1: indicaciones generales.
Parte 2: turbinas de vapor y generadores que superen los 50 MW con velocidades
típicas de trabajo de 1.500, 1.800, 3.000 y 3.600 RPM.
Parte 3: maquinaria industrial con potencia nominal por encima de 15 kW y velo
cidades entre 120 y 15.000 RPM.
Parte 4: conjuntos movidos por turbinas de gas que excluyen a las empleadas en
aeronáutica.
Parte 5: conjuntos de máquinas en plantas de hidrogeneración y bombeo.
La condición de la vibración se clasifica conforme a los siguientes parámetros: tipo
de máquina, potencia o altura de eje y flexibilidad del soporte.
Debido a las características típicas de construcción de las máquinas, se establecen
grupos para la selección de cada tipo de ellas: grupo 1, grupo 2, grupo 3 y grupo 4,
tal como se muestra en la tabla 3.
La magnitud de vibración se define dentro de este grupo de estándares como el
valor máximo de velocidad RMS de banda ancha, en el rango de frecuencias espe
cificado (típicamente de 10 a 1.000 Hz) evaluado sobre la estructura en puntos
prescritos. Se debe tener en cuenta que otras cantidades como el desplazamiento,
los valores de aceleración y que toman valores pico, en lugar de valores RMS, es-
tán permitidos, pero pueden no correlacionarse fácilmente a criterios basados en
valores RMS.
Los criterios para evaluar la severidad de la vibración incluyen magnitud de la
vibración y cambios en la magnitud de la vibración (Robichaud & Eng, 2009). A
continuación, las zonas de evaluación se definen para permitir una apreciación de
la vibración y para proporcionar directrices sobre posibles acciones.
Zona A: valores de vibración de máquinas recién puestas en funcionamiento o
reacondicionadas.
Zona B: máquinas que pueden funcionar indefinidamente sin restricciones.
Zona C: la condición de la máquina no es adecuada para una operación continua,
sino solamente para un período de tiempo limitado. Se deberían llevar a cabo me
didas correctivas en la siguiente parada programada.
Zona D: los valores de vibración son peligrosos, la máquina puede sufrir daños.
Capítulo 1: Modelos matemáticos y contexto de las vibraciones 23
Tabla 3. Niveles de vibración.
in/s mm/s
RMS RMS
0,43 11 C
0,28 7,1 B
0,18 4,5
0,14 3,5 A
0,11 2,8
0,09 2,3
0,06 1,4
0,03 0,71
Base Rígida Flexible Rígida Flexible Rígida Flexible Rígida Flexible
Grandes
Tamaño medio
Bombas . 15 kW máquinas
15 kW , P
Flujo radial, axial o mixto 300 kW , P
Tipo de , 300 kW
, 50 MW
máquina
Motores
Motor Motor Motores
160 mm , H
integrado separado 315 mm , H
, 315 mm
Grupo Grupo 4 Grupo 3 Grupo 2 Grupo 1
Máquina en buenas La máquina no podrá
A C
condiciones operar en forma continua
La máquina puede operar La vibración está
B D
indefinidamente provocando daños
Fuente: (Standardization, ISO 10816-8:2014, 2014).
24 Análisis de vibraciones en maquinaria
Tabla 4. Organismos de normalización de vibraciones.
Principales organismos normalización
Organismo Ámbito Nombre Dirección web
Internacional International
IEC [Link]
(sede Suiza) Electrotechnical Comission
Internacional International Organization [Link]
ISO
(sede Suiza) For Standarization [Link]
Asociación Española de
Aenor España [Link]
Normalización y Certificación
American National
ANSI USA [Link]
Standards Institute
[Link]
BSI Gran Bretaña Britis Standards Institution
com/en-GB/
Deutsches Institut
DIN Alemania [Link]
Für Normung
Japanese Standards
JSA Japón [Link]
Association
Association of German
VDI Alemania [Link]
Engineers
American Petroleum
API USA [Link]
Institute
Fuente: Autores, 2018.
C APÍTULO 2
TO M A D E D ATO S
D E V I B R AC I O N E S
Medir adecuadamente las vibraciones de una máquina es el paso más importante
para el análisis de vibración, se necesitan datos confiables con el fin de lograr resul
tados seguros para que quien revise pueda realizar diagnósticos exactos, para lo
grarlo, los datos deben tomarse cada vez con el mismo sensor, el mismo montaje
y con la máquina en las mismas condiciones de operación. A continuación, se
describirán las recomendaciones más importantes en la toma de datos.
Los datos no serán considerados como óptimos si la máquina no está funcionando
bien, si se está usando un sensor incorrecto o si el sensor ha sido instalado inco-
rrectamente. Hay dos partes que se deben tener en cuenta para recolectar buenos
datos:
»» Asegurarse que el transductor correcto está montado y está siendo usado
correctamente.
»» Asegurarse que las medidas se hacen de la misma manera cada vez.
2.1. Pasos previos a la medición
»» Seleccione el tipo correcto de transductor.
»» Seleccione entre desplazamiento, velocidad o aceleración.
»» Seleccione el modelo con la estructura que se acomode al ambiente de me
dición.
26 Análisis de vibraciones en maquinaria
»» Considere el tipo de máquina, la naturaleza de su superficie y la accesibili
dad, todo esto para identificar la mejor posición del transductor.
»» Identificar las condiciones de carga o velocidad de la máquina, ya que se
espera que tales condiciones sean las mismas en cada medición para au
mentar la confiablidad de los datos.
2.2. Puntos de medición
»» Los puntos de medición deben estar en sitios de acceso seguro.
»» Las mediciones deben realizarse en diferentes puntos de la máquina, las
mejores partes son los rodamientos. Estos puntos y su dirección deben
identificarse en la figura 12.
»» El sensor mide en una sola dirección y para obtener un panorama comple
to del movimiento se debe hacer una medición vertical, una axial y una
horizontal.
En general, estos puntos pueden nombrarse de la forma que convenga, pero una
gran parte de analistas utilizan la siguiente designación:
Rodamiento 1 dirección medida 1 máquina (o parte de ella) medida
Al utilizar un número para designar la posición del rodamiento, ‘1’ representa el
primer rodamiento, ‘2’ el segundo, etc.
Las direcciones se presentan como “VER”, “HOR” y “AXI” para expresar una me-
dición vertical, horizontal y axial, respectivamente. Una medición común sería:
‘1 VER MOTOR’ equivalente a decir: medición del rodamiento
uno del motor, en la dirección vertical.
Figura 12. Puntos de medición en motor.
Sensor
Fuente: Autores, 2018.
Capítulo 2: Toma de datos de vibraciones 27
»» De la forma más directa, se debe ubicar el sensor hacia el rodamiento para
que la vibración viaje sin intermitencias, ya que en la máquina existen
muchas partes que vibran y algunas de estas pueden resonar sin que repre-
senten fuerzas en los rodamientos, lo que los hace adecuados para medir.
»» Marcar o etiquetar claramente los sitios de medición, esto con el fin de au
mentar la repetitividad de los datos.
2.3. Problemas al inicio de la medición
Es importante hacer una revisión de los datos, justo como fueron obtenidos antes
de ser cargados a la computadora, con el fin de identificar si estos son razonables.
Un cable en mal estado, una falla en la conexión de montaje o el movimiento du-
rante la medición pueden ocasionar fallas en la medición que pueden detectarse,
un ejemplo se puede ver en la imagen 1.
Imagen 1. Verificar el buen estado del equipo.
Fuente: Toll Morris Inc, 2018.
Problemas de estabilización: cuando se aplica energía al transductor la señal ten-
drá un sobresalto, en el caso de un acelerómetro la señal se vuelve inestable mien-
tras el circuito interno amplificador se establece. Cualquier dato que se recolecte en
ese período de tiempo tendrá una onda que presentará un desbalance y el espectro
tendrá una caída pronunciada con levantamiento (caída ski), (figura 13).
28 Análisis de vibraciones en maquinaria
Figura 13. Medición defectuosa con efecto ski.
Amplitud
Caída descendente del espectro (ski)
Frecuencia
Fuente: Autores, 2018.
Para que este fenómeno no ocurra, el tiempo de estabilización establecido en software
deberá incrementarse.
»» Impactos: se debe tener cuidado al usar imanes como montura porque
un mal manejo ocasiona que el sensor sea golpeado. Esto producirá un
voltaje de salida más alto, el espectro tendrá una característica de caída
y la onda presentará una oleada desbalanceada, por tanto, se debe ubicar
suavemente el sensor, ayudándose con un punto de apoyo.
Figura 14. Posicionamiento correcto para evitar impacto.
Posicionar cuidadosamente
acelerómetros con imanes
Fuente: Autores, 2018.
Capítulo 2: Toma de datos de vibraciones 29
»» Problemas térmicos: cuando el sensor se mueve de una, entre superficies,
a una diferencia importante de temperatura, en ese momento se presenta
un choque térmico que sobreexcita la salida de datos del sensor, por lo que
se debe esperar a que se equilibre la temperatura del sensor para ejecutar
la medición.
2.4. Proceso de medición
Cuando se monta un transductor en una máquina, la salida eléctrica es una señal
analógica continua que contiene toda la información del desplazamiento, la velo-
cidad y la aceleración de la vibración. Esta medida debe ser digitalizada para que
un dispositivo colector registre y almacene la vibración, para ello, la señal analó-
gica continua es convertida en un conjunto de números discretos, denominado
registro de tiempo, con el cual se reconstruye la señal analógica.
La frecuencia de muestreo afecta el rango de la onda que se puede reconstruir, es
necesario que esta sea muestreada a una frecuencia mayor que el doble de la fre-
cuencia de la onda, la cual se conoce como frecuencia Nyquist; la mayoría de ana-
lizadores toman esta frecuencia como 2,56 veces la frecuencia máxima de interés.
Una vez la señal ha sido digitalizada mediante el proceso descrito anteriormente,
se puede reconstruir y mostrar, como también exponerla en datos de frecuencia
mediante la transformada rápida de Fourier y así desplegarla como un espectro.
El número de muestras de la onda de tiempo determina el espesor o la resolución
del espectro. Este espesor o resolución lo determina el usuario mediante la introduc-
ción del número de líneas llamadas líneas de resolución (LOR) para una frecuencia
máxima. Un ejemplo de ello es un espectro con una frecuencia máxima de 1.000 Hz
y una selección de 400 líneas, donde tendrá una resolución de 1.000/400 5 2,5 Hz.
Esto quiere decir que dos frecuencias aparecerán como una sola si están a menos
de 2,5 Hz.
Se debe tener en cuenta que la frecuencia de muestro es 2,56 la frecuencia máxima,
entonces las muestras necesarias serán 2,56 veces el número de LOR. Esto se tra-
duce en que, a mayor resolución, mayor tiempo del proceso.
Luego, la señal pasa por un proceso que ayuda a su identificación llamado: filtrado
de señal:
Los filtros principales son:
»» Filtros de paso bajo: permiten el paso de bajas frecuencias.
»» Filtros de paso de banda: permiten el paso de frecuencias dentro de una
banda.
30 Análisis de vibraciones en maquinaria
»» Filtros de paro de banda: bloquean el paso de frecuencias en una banda
específica.
»» Filtros de paso alto: permiten el paso de altas frecuencias.
C APÍTULO 3
A N Á L I S I S D E V I B R AC I O N E S
E N M AQ U I N A R I A R OTAT I VA
El análisis de vibración es una técnica importante en el mantenimiento predictivo,
estudia el funcionamiento de maquinaria rotativa y proporciona la identificación
de una variedad de problemas que pueden afectar el correcto funcionamiento de
esta. Aunque todas las máquinas presentan niveles de vibración aceptables en
condiciones de funcionamiento, se debe priorizar un seguimiento correcto, pues
cuando estos niveles exceden estados vibratorios elevados aparece el mal funcio-
namiento de las partes o, en el peor de los casos, el daño total de la máquina. Esto
ocurre a medida que las máquinas se deterioran, lo que acarrea que haya desgaste
y soltura en las partes y puede haber fuerzas de excitación producidas debido a al-
gunos fallos, tales como desbalances y desalineación. Una correcta aplicación del
análisis de vibraciones impide que estas anomalías se presenten esporádicamente,
lo que lleva un control adecuado en su implementación.
3.1. ¿Por qué es importante un análisis de vibraciones?
Los componentes de una estructura o máquinas sometidos a vibración pueden fa
llar debido a la fatiga del material producida por la vibración cíclica inducida por
los esfuerzos, además, la vibración provoca un desgaste más rápido de las partes de
la máquina como cojinetes y engranes que generan ruido excesivo. En las máqui-
nas, la vibración puede aflojar las sujeciones como tuercas y separar la soldadura
o los puntos de anclaje. En procesos de maquinado, la vibración puede provocar
rechinidos lo que conduce a un acabado superficial pobre.
32 Análisis de vibraciones en maquinaria
También en las máquinas, cuando la frecuencia natural de la vibración de su estruc
tura coincide con la frecuencia de la excitación externa, se presenta un fenómeno
conocido como resonancia, el cual conduce a continuas deflexiones y fallas. Por
lo tanto, uno de los propósitos del estudio de la vibración es reducirla mediante el
diseño o el mantenimiento, de modo que el estado de desequilibrio en la máquina
sea insignificante.
Las mediciones de vibración pueden usarse para monitorear partes de una má-
quina y relacionar las características del patrón de vibración con los parámetros
asociados a una falla. Por lo cual, analizar los problemas y su causa desde la pers-
pectiva de la condición del equipo, antes de que estas fallas colapsen la máquina
y que al final se conviertan en costos de mantenimiento y producción, constitu
ye una herramienta sumamente útil en la planeación de estrategias de manteni-
miento asociadas a fallas específicas.
Concretamente, la frecuencia y la velocidad angular de las máquinas rotativas se
relacionan mediante un valor único para cada elemento, lo cual hace posible su
identificación en un espectro, y como se verá más adelante, varios problemas se
presentan a la frecuencia de velocidad de giro o en sus múltiplos.
“El analista y el software de vibración deben buscar cambios en el patrón de
vibración mediante la toma de mediciones a intervalos regulares”.
3.2. El espectro de vibraciones
Cada máquina en funcionamiento tiene una forma de movimiento específico que
se puede evidenciar en sus vibraciones, para analizar este movimiento se deben
colectar estas, las cuales son de varios tipos. Cuando se utiliza un colector de vibra
ciones estas aparecen mezcladas, por lo que no se puede identificar su origen.
Afortunadamente, existen procesos como la transformada de Fourier y la transfor
mada rápida de Fourier, las cuales permiten separar cada onda a su frecuencia es
pecífica, estas tienen origen en la serie de Fourier, por la que cualquier función
periódica x(t) continua en el tiempo puede expresarse como la suma de términos
senos y cosenos que son funciones armónicas (James, 1994).
a0 `
x(t) 5 5 ∑ (an cos nvt 1 bn sin nvt )
2 n51
Los términos an cos nvt y bn sin nvt son llamados “armónicos”, estos armónicos de
orden n se grafican con líneas verticales de amplitud an o bn contra la frecuencia vt.
Capítulo 3: Análisis de vibraciones en maquinaria rotativa 33
Figura 15. Representación de una función periódica en el dominio de la frecuencia.
x(t) 5 A sin(vt 1 w)
an
v 2v 3v 4v 5v Frecuencia nv
Fuente: Autores, 2018.
Por ejemplo, la figura 15 muestra la representación en el dominio de la frecuencia
de la función armónica x(t) 5 A sin(vt 1 w), para el caso A 5 an.
En general, el resultado de este procedimiento, llamado espectro de vibración (fi-
gura 16), trata de obtener las diferentes frecuencias a las que se están produciendo
movimientos, de esta manera y conociendo las partes de la máquina se puede iden-
tificar el origen de cada vibración, como también su importancia según la ampli-
tud. Este es el principio de funcionamiento del software analizador de vibraciones
(White, 2010).
34 Análisis de vibraciones en maquinaria
Figura 16. Espectro de vibración.
La vibración compleja 5 suma de Espectro de vibración
vibraciones simples
Dominio de
la frecuencia
Dominio del tiempo
Señales en el dominio del tiempo
y en el dominio de la frecuencia (espectro)
Fuente: Tutorial de vibraciones A-MAQ, 2005.
3.3. Identificar la frecuencia de giro
El análisis de espectros en maquinaria rotatoria se hace con base en la frecuencia de
giro del motor, por lo que identificarla es el inicio del proceso. Comúnmente a esta
frecuencia se le llama “1x” (uno equis) y el resto de frecuencias de interés pueden
ser fácilmente expresadas relacionadas a esta. De manera que si la frecuencia de
velocidad de giro es 1x 5 30Hz, una frecuencia de 90 Hz será 3x y una de 270 Hz
corresponderá a 9x (figura 17).
Figura 17. Identificación de la frecuencia de giro.
Amplitud
Frecuencia de giro para un motor a 1.800 RPM
1x
3x
2x
4x
30 Hz 60 Hz 90 Hz 120 Hz Frecuencia
Fuente: Autores, 2018.
Capítulo 3: Análisis de vibraciones en maquinaria rotativa 35
3.3.1. Fase
Una adecuada lectura de fase permite hacer un diagnóstico de averías claro, donde
se identifican y se diferencian los problemas mecánicos que en el espectro lucen
similares, como pueden ser: el desequilibrio, la desalineación, el eje deformado, la
resonancia, las holguras, la falta de rigidez en la base y la excentricidad.
3.3.2. Armónicos
Los armónicos son una serie de picos espaciados equitativamente que son múlti-
plos del primer pico en la serie, donde el primer pico es la frecuencia fundamental.
Los armónicos más comunes son múltiplos de la velocidad de giro del eje; sin em
bargo, pueden ser múltiplos de cualquier frecuencia.
3.3.3. Bandas laterales
Las bandas laterales son el resultado de una modulación de amplitud entre dos se-
ñales. Se presentan comúnmente al estudiar elementos giratorios del eje, la caja de
engranaje y el sistema eléctrico cuando hay resonancia. En el espectro, las bandas
laterales se asemejan a picos espaciados equitativamente, centrados en otro pico
llamado “frecuencia central”.
La frecuencia central puede ser ocasionada por la frecuencia del engranaje, la de
fallo de pista interna y las de paso de falla del rotor. Las frecuencias de bandas la
terales incluyen la velocidad de giro, la de jaula y la de paso de polos.
3.3.4. Ruido de piso
El ruido de piso representa la vibración del ambiente y comúnmente sus niveles
están cerca del cero con picos que se levantan levemente; sin embargo, habrá casos
en los que estos picos parecieran levantarse de forma pronunciada sobre el nivel
de referencia de amplitud cero, sobre todo en ciertas áreas. Esto puede deberse
generalmente a dos tipos de falla:
»» Puede tratarse de un desgaste de rodamiento serio, lo que incrementa el
levantamiento de ruido de piso y hace que el problema sea más serio.
»» Puede haber resonancia, con frecuencia los picos del espectro tienen una
base muy ancha y se reflejan como acumulación debido a resonancia, si
no hay otras.
36 Análisis de vibraciones en maquinaria
3.4. Comparación con medidas de referencia
Las comparaciones entre datos actuales y datos antiguos es una ayuda real para es
tablecer la condición de una máquina al determinar cuándo y cuánto han cambia
do los datos.
»» Superposición de gráficas: la mejor manera de ver los datos es superpo-
niendo las gráficas. Entonces, se puede ver claramente cuánto han cam-
biado los datos en cada pico del espectro. Tres o cuatro gráficas se pueden
desplegar de esta manera y aún así identificar de qué espectro se derivó
cada dato.
»» Tendencia: los gráficos de tendencia proveen una rápida mirada a los cam
bios que están ocurriendo. La mayoría del software provee una gráfica de
nivel global de vibración (figura 18).
En el análisis de vibración es necesario aplicar tendencias en el área de espectro
que nos interesen, por ejemplo, a la frecuencia de giro, frecuencia de aspas, engra-
naje y frecuencia de rodamiento.
Figura 18. Gráficos de tendencia.
NBF O-PK(mm/s)
30
20
10
0
0 50 200 150 200 250 300 350 400
Fuente: Tutorial de vibraciones A-MAQ, 2005.
Cuando una de estas áreas muestra un incremento de valor, entonces se comienza
a analizar el espectro y la onda para identificar dónde ha iniciado el incremento.
»» Gráficos apilados: algunos softwares permiten ver los gráficos desplegados
verticalmente. La diferencia en cada nivel de cada pico principal se puede
ver de manera clara al tiempo que se observa el espectro completo (figura
19), se tiene la desventaja de solo poder observar los datos que quepan en
la pantalla.
Capítulo 3: Análisis de vibraciones en maquinaria rotativa 37
Figura 19. Espectros desplegados apilados.
Lin g 5kHz gSE Peak-Peak Mag
1.
1xBPIR
1.
8. 1xBFS 2xBPIR
5.
6.
0. 20.000 40.000 60.000.
X: 8.925 CPM 1xBPIR 30/10/2002 [Link] p.m Lin
Y: 678.115
4. 2xBPIR
1xBFS
2.
0. 20.000 40.000 60.000.
CPM 30/10/2002 [Link] p.m
Lin
X: 24.150 Y: 0378393
Fuente: Ingeniería de mantenimiento, 2016.
»» Gráficos de cascada: una manera común y usada para observar los datos es
desplegar una gran cantidad de espectros en forma de cascada, este plano
permite ver cómo han cambiado los picos de las frecuencias desplegadas a
lo largo del tiempo (figura 20), por lo que ofrecen información importante
de las tendencias. Esto es una herramienta importante para revisión sobre
cómo fue el desarrollo de la falla y qué niveles se alcanzaron después de la
reparación.
38 Análisis de vibraciones en maquinaria
Figura 20. Visualización en cascada.
Preview
Signal FFT
Configuration
Zoom X Zoom Y
22 1.0
✓ 30 View ✓ Names 2013-04
✓ Severty
Parameter
2013-04-01
mm/s
Print 2013-03-15
✓ Signal ✓ FFT
2013-03-01
2013-02-15
12.9
✓ 2013-04-27
✓ 2013-04-01 2013-02-01
✓ 2013-03-15
✓ 2013-03-01 2013-01-15
✓ 2013-02-15
✓ 2013-02-01
✓ 2013-01-15 -23828 CPM Generate report
Fuente: Erbessd-Instruments, 2018.
»» Visualización logarítmica: la escala de gráficos logarítmicos se usa con
frecuencia para resaltar armónicos y bandas laterales. Esto se debe a que la
escala logarítmica despliega amplitudes pequeñas en contraste con ampli
tudes grandes, también es muy útil para hacer comparaciones gráficas con
datos de referencia.
3.5. Condiciones de falla en máquinas rotativas
Existen diversos documentos que registran las características de espectros asocia
dos a fallas comunes, unos de estos documentos son las tablas de Charlotte, des-
cripciones que se incluyen a continuación.
3.5.1. Desbalance
Es una condición donde el eje geométrico no pasa por el centro de masa de un
cuerpo en rotación, donde el centro de masa es el punto en el que toda esta se distri
buye equitativamente, lo que significa que hay un punto pesado a lo largo del eje.
El punto más pesado en el rotor produce una fuerza centrífuga en los rodamientos
al girar. Los principales tipos de desbalanceo son:
Capítulo 3: Análisis de vibraciones en maquinaria rotativa 39
»» Desbalance estático: ocurre en elementos rotativos para los cuales su largo
es despreciable en comparación con su diámetro, ya sea por acumulación
de material o pérdida del mismo, en todo caso, debido al cambio en la dis-
tribución de su masa. El espectro presenta vibración dominante en 1x.
En la ilustración se observa espectro de vibración para un motor eléctrico a 30 Hz,
acompañado de un disco para la simulación del desbalanceo estático, las medicio-
nes se toman en base a las escalas principales RMS, 0-pk y pk-pk. Esta prueba se
realizó con el equipo de vibraciones de la Universidad Francisco de Paula Santander
Ocaña.
Imagen 2. Equipo de pruebas para análisis de vibración con disco desbalanceado.
Fuente: Autores, 2018.
Figura 21. Representación del espectro por desbalanceo estático, realizada por equipo a-predictor.
Punto: 1 HOR
1
0,8
0,6
0,4
0,2
0
0 200 400 600 800 1.000 1x dominante en
Espectro ISO RMS (mm/s) vs RPS ambas direcciones,
con niveles altos
de vibración en
Punto: 1 VER
dirección vertical
4 por desbalanceo.
3
2
1
0
0 200 400 600 800 1.000
Espectro ISO RMS (mm/s) vs RPS
Fuente: Autores, 2018.
40 Análisis de vibraciones en maquinaria
Imagen 3. Equipo de pruebas y análisis con desbalanceo dinámico.
Discos Acelerómetro
Acople
Rodamientos
Fuente: Autores, 2018.
»» Desbalance dinámico: 1x domina el espectro, la diferencia de fase radial
entre los rodamientos externos e internos puede estar en cualquier lugar
en el rango de 0-180°; sin embargo, la diferencia de fase horizontal debe
cuadrar con la diferencia de fase vertical.
Diagnóstico: un rotor desbalanceado generará una vibración en la frecuencia de
la velocidad de giro debido a la fuerza centrífuga del desbalance de masa. Por lo
tanto, se espera que una máquina con condición de desbalance genere una onda
sinusoidal y un pico dominante en el espectro a la frecuencia de la velocidad de
giro 1x. Se tiene que determinar esta condición al apoyarse en los niveles históri-
cos de medición.
Capítulo 3: Análisis de vibraciones en maquinaria rotativa 41
Figura 22. Representación del espectro de falla por desbalanceo dinámico.
Punto: 1 HOR
0,8
0,6
0,4
0,2
0
0 200 400 600 800 1.000
Espectro ISO RMS (mm/s) vs RPS
1x dominante en
dirección vertical
Punto: 1 VER
y horizontal, el
1,5 desbalanceo
1 dinámico se
0,5 evidencia con
0
0 200 400 600 800 1.000 la excitación
Espectro ISO RMS (mm/s) vs RPS de fuera de fase
de los niveles
Punto: 1 AXI en la dirección
axial.
0,8
0,6
0,4
0,2
0
0 200 400 600 800 1.000
Espectro ISO RMS (mm/s) vs RPS
Punto: 2 HOR
1,5
1
0,5
0
0 200 400 600 800 1.000
Espectro ISO RMS (mm/s) vs RPS
Fuente: Autores, 2018.
3.5.2. Desalineamiento
Ocurre cuando las líneas de rotación de dos ejes que se unen mediante un acople
no son colineales.
42 Análisis de vibraciones en maquinaria
Desalineamiento angular: tiene lugar cuando las líneas de rotación de dos ejes no
colineales se proyectan formando un ángulo entre ellas.
Desalineamiento paralelo: ocurre cuando las líneas de rotación de dos ejes no co
lineales están paralelamente separadas.
Diagnóstico
Desalineamiento angular: produce una fuerte excitación en 1x y algo de excita
ción en 2x de la dirección axial en ambos rodamientos. Habrá niveles radiales
fuertes (horizontal y vertical) en 1x y 2x; sin embargo, ambos estarán en fase.
Desalineamiento paralelo: altos niveles de vibración en 1x y 2x en la dirección
radial (horizontal y vertical) en los rodamientos de cada lado del acoplamiento.
Con una frecuencia componente 2x que será mayor a 1x y que estará fuera de fase
en la dirección axial y radial.
3.5.3. Soltura estructural
Ocurre cuando hay debilidad en la base, en la fundición o en las patas, haciendo que
la máquina se balancee de un lado a otro.
Diagnóstico: la soltura entre una máquina y sus bases incrementará el componen
te de vibración 1x de menor rigidez, usualmente la horizontal, siendo esta con
dición más evidente si 1x horizontal es mayor a 1x vertical.
Capítulo 3: Análisis de vibraciones en maquinaria rotativa 43
Figura 23. Representación del espectro de falla por soltura estructural.
Punto: 1 HOR
0,5
0,4
0,3
0,2
0,1
0
0 200 400 600 800 1.000
Espectro ISO RMS (mm/s) vs RPS
Punto: 1 VER
Aparecen
2 múltiples
1,5
arnómicos y
1
0,5 subarmónicos
0 de 1x en todas
0 200 400 600 800 1.000 las direcciones.
Espectro ISO RMS (mm/s) vs RPS
Punto: 1 AXI
0,5
0,4
0,3
0,2
0,1
0
0 200 400 600 800 1.000
Espectro ISO RMS (mm/s) vs RPS
Punto: 2 HOR
1
0,5
0
0 200 400 600 800 1.000
Espectro ISO RMS (mm/s) vs RPS
Fuente: Autores, 2018.
3.5.4. Falla en rodamientos
Los rodamientos generan frecuencias diferentes basadas en su construcción física.
Hay cuatro frecuencias de fallo de interés:
44 Análisis de vibraciones en maquinaria
Frecuencia de fallo de pista interna (BPFI):
BPFI: v ?
[ (
1
2
d
n 1 1 cos a
D )]
Frecuencia de fallo de pista externa (BPFO):
BPFO: v ?
[(1
2
d
n 1 2 cos a
D )]
Frecuencia de tren o frecuencia de jaula (FT):
BF: v ?
[(
1
2
d
1 2 cos a
D )]
Frecuencia de giro de bola (BSF):
[ ( ( ) )]
2
1D d
BSF: v ? 11 cos a 2
2d D
Información necesaria: número de bolas (n), diámetro de bola (d), diámetro de
paso (D), ángulo de contacto (a) y velocidad de giro del eje (v).
Diagnóstico: una falla en rodamientos produce picos que no son múltiplos ente-
ros de la frecuencia de velocidad de giro 1x, es decir, son asincrónicos.
En la figura 24 se representa el espectro para un rodamiento con fisura, con una
frecuencia de giro del eje de 30 Hz y la representación en la escala RMS para la
velocidad y la escala 0-pk para la aceleración. Se evidencian en el espectro medido
sobre el punto 1 vertical las distintas frecuencias del rodamiento con excitación, lo
cual implica que existe un choque metal-metal de las bolas y la jaula, causadas por
la discontinuidad en la fisura, lo que genera un sobresalto a la frecuencia indicada.
Capítulo 3: Análisis de vibraciones en maquinaria rotativa 45
Figura 24. Representación del espectro de falla en rodamiento con fisura.
Punto: 1 HOR
1 80
60
0,5 40
20
0 0
0 200 400 600 800 1.000 0 200 400 600 800 1.000
Espectro ISO RMS (mm/s) vs RPS Espectro ISO RMS (mm/s) vs RPS
Aparecen las frecuencias de fallo, en los altos niveles del espectro
de aceleración. Visibles en todas las direcciones.
Punto: 1 VER
1,5 60
1 40
0,5 20
0 0
0 200 400 600 800 1.000 0 200 400 600 800 1.000
Espectro ISO RMS (mm/s) vs RPS Espectro ISO RMS (mm/s) vs RPS
Aparecen las frecuencias de fallo, en los altos niveles del espectro
de aceleración. Visibles en todas las direcciones.
Punto: 1 AXI
0,8 300
0,6 200
0,4 100
0,2
0 0
0 200 400 600 800 1.000 0 200 400 600 800 1.000
Espectro ISO RMS (mm/s) vs RPS Espectro ISO RMS (mm/s) vs RPS
Aparecen las frecuencias de fallo, en los altos niveles del espectro
de aceleración. Visibles en todas las direcciones.
Punto: 2 HOR
1 80
60
0,5 40
20
0 0
0 200 400 600 800 1.000 0 200 400 600 800 1.000
Espectro ISO RMS (mm/s) vs RPS Espectro ISO RMS (mm/s) vs RPS
Aparecen las frecuencias de fallo, en los altos niveles del espectro
de aceleración. Visibles en todas las direcciones.
Fuente: Autores, 2018.
46 Análisis de vibraciones en maquinaria
3.5.5. Caja de engranajes
»» Diagnóstico: hay tres frecuencias involucradas, la de velocidad de entra
da, la del engranaje y la de salida.
»» Frecuencia del engranaje: es igual al número de dientes por la velocidad
del eje.
»» Velocidad de salida: (velocidad de entrada por dientes de entrada) entre
dientes de salida.
»» El análisis de la onda en el tiempo es una herramienta poderosa cuando se
analizan mediciones de cajas de engranajes. Cuando un diente hace con-
tacto con otro hay un pulso en la onda y los dientes se pueden contar en
la onda de tiempo.
3.5.6. Recomendaciones básicas para el analista de vibraciones
»» Verificar el estado de sujeción de la máquina.
»» Identificar las solturas entre superficies.
»» Establecer cuáles son los datos de funcionamiento correcto.
»» Observar la tendencia (gráficos en cascada).
»» Identificar las fallas relacionadas.
»» Tener en cuenta un análisis de la fase, ya que nos dice cómo se mueve un
punto de medición con respecto a otro punto de medición.
»» La pregunta que todo analista debe hacerse es: ¿cuáles tendencias cam-
biaron cuando empezó a manifestarse la falla?
C APÍTULO 4
CASOS DE ESTUDIO:
E X PA N S I Ó N D E F O U R I E R
Figura 25. Señal de vibración en el tiempo.
1
mm/s
21
22
0 100 200 300 400 500 600 700 800 900 1.000
Fuente: Autores, 2019.
Para el caso de funciones periódicas es posible realizar un análisis a partir de datos
recolectados por el transductor de vibraciones que permitan idealizar el compor-
tamiento de la señal analógica, mediante el proceso de expansión de Fourier, ya
que no se conoce la función que describe el movimiento vibratorio combinado del
equipo, por tanto, al considerar que, en su mayor rango los datos obtenidos desde
la máquina girando a velocidad constante son periódicos, se utilizará un proceso
numérico para determinar dicha función.
48 Análisis de vibraciones en maquinaria
Como se describió en el capítulo anterior, la expansión de Fourier es:
a0 `
x(t) 5 1 ∑ (an cos nvt 1bn sin nvt)
2 n51
Para lo cual, se considera el conjunto de datos muestra desde la señal indicada en
la figura 13 y que aumentada en período corresponde a una tasa de diez muestras,
donde n 5 10 para un tiempo t 5 0,01 s, (figura 25).
Tabla 5. Valores muestrales para la expansión de Fourier.
Amplitud (RMS) Tiempo (s)
21,1 0,01
1,9 0,02
20,7 0,03
0,8 0,04
2,4 0,05
20,6 0,06
22,2 0,07
0,7 0,08
0,5 0,09
0 0,1
Fuente: Autores, 2018.
Figura 26. Fracción de la señal de vibración correspondiente a un período de 0,1 s.
Amplitud
Tiempo x1022
1 2 3 4 5 6 7 8 9 10
21.5
Fuente: Autores, 2018.
Capítulo 4: Casos de estudio: expansión fourier 49
Al usar Matlab podemos definir un algoritmo que calcula los coeficientes de la ex-
pansión de Fourier a partir de datos de muestra con un período y una frecuencia
específicos, según la tabla 5.
Frecuencia W 5 n 3 2 3 pi/T
A continuación se describe un código programado en Matlab R2016b, para la ob
tención la expansión:
%{ Cálculo de los coeficientes de Fourier desde el análisis numérico de muestras %}
t 5 0.01:0.01:0.1; % de intervalos de tiempo
p 5 [21.1,1.9,20.7,0.8,2.4,20.6,22.2,0.7,0.5,0]; % de muestras por intervalo
T 5 0.1; % de período
N 5 10; % cantidad de muestras
F50
An 5 0; Ao 5 0; Bn 5 0; P 5 0
n 5 1:10
for j 5 1: n
for I 5 1: length(t)
ao 5 (2/N) 3 p(i)/2; % de cálculo de coeficiente a0
Ao 5 Ao 1 ao;
a 5 (2/N) 3 p(i) 3 cos((2 3 n 3 pi 3 t(i)/T)); % de cálculo de coeficiente an
An 5 a 1 An
b 5 (2/N) 3 p(i) 3 sin((2 3 n 3 pi 3 t(i)/T)); % de cálculo de coeficiente bn
Bn 5 b 1 Bn
PS 5 p(i) 1 P
end
50 Análisis de vibraciones en maquinaria
W 5 n 3 2 3 pi/T
% visualizar variables
Ao
An
Bn
end
% de cálculo de la ecuación de la expansión
for ii 5 1: length(t)
for jj 5 1: 3
PP 5 Ao 1 (An(jj) 3 cos((2 3 pi/T) 3 t(ii) 3 jj)) 1 (Bn(jj) 3 sin((2 3 pi/T) 3 t(ii)
3 jj))
F 5 F 1 PP
end
FF(ii) 5 PP
end
plot(t,FF,’--r’)
%fprintf(‘la expansión de la serie es: \n %2.5g %2.5g %2.5g %2.5g %2.5g %2.5g
%2.5g’,Ao,An,Bn).
Capítulo 4: Casos de estudio: expansión fourier 51
Figura 27. Representación de la onda en la expansión de Fourier para n 5 3.
1.5
0.5
20.5
21
21.5
0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 0.11
Fuente: Autores, 2018.
Figura 28. Representación de la onda en la expansión de Fourier para n 5 6.
0.8
0.6
0.4
0.2
20.2
20.4
20.6
0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 0.11
Fuente: Autores, 2018.
52 Análisis de vibraciones en maquinaria
Figura 29. Representación de la onda en la expansión de Fourier para n 5 10.
0.8
0.6
0.4
0.2
20.2
20.4
0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 0.11
Fuente: Autores, 2018.
4.1. Observaciones
En las figuras 27 y 28 se observa que la expansion está desproporcionada, debido a
que no se aplicó ningún algoritmo de ajuste de datos para obtener una representa
ción más precisa de la forma real de la señal, como se ve en la figura 25. Conse
cuentemente, cuando se iguala la cantidad de armónicos calculados al número de
muestras, esto es cuando n 5 10, la expansión empieza a tomar la forma sinusoidal
y la representación de la señal original se vuelve pobre. De esta manera, este está
limitado por la resolución que sea posible obtener y de qué tan representativa es la
cantidad de las muestras. Además, para que el rango de amplitudes se corresponda
se debe ajustar un algoritmo en el proceso para los datos obtenidos por la expan-
sión, el cual está en escala cero-pico desde a la escala RMS que está configurada en
el equipo desde el que se hicieron las mediciones. Este proceso es bastante complejo
e implica un tratamiento de señales digitales, lo cual sale del propósito de este do-
cumento.
Capítulo 4: Casos de estudio: expansión fourier 53
4.2. Expansión de Fourier para describir el espectro de un
motor desbalanceado
Esta ecuación describe el movimiento de vibración libre sin amortiguamiento y
con una fuerza de excitación:
Mx 1 kx 5 F0 sen wt
Donde:
M 5 66,096 kg (masa total del sistema 2 conjunto motor)
m0 5 0,01198 kg (masa del desbalanceo)
e 5 0,074 m (radio del desbalanceo)
w 5 1.800 RPM 5 188,495 rad/s (velocidad de excitación)
fn 5 20 Hz (frecuencia natural)
wn 5 2pfn 5 2p(20 Hz) 5 125,663 rad/s
F0 5 m0ew2
Ecuación de amplitud máxima:
m0ew 2
Xmáx 5
k 2 Mw 2
k 5 Mwn2
Al remplazar:
k 5 1043734,46 N/m
Xmáx 5 20,00002414 m 5 2002414 mm
Ecuación total del movimiento:
x(t) 5 C1 cos(wnt) 1 C2 sen(wnt) 1 Xmáx cos(wt)
Se supone que la vibración parte desde el reposo en condiciones iniciales:
m0ew 2 x
C1 5 x0 2 ; C2 5 0
k 2 Mw 2
wn
t 5 0 → x0 5 0; x0 5 0
54 Análisis de vibraciones en maquinaria
La ecuación total se reduce a:
x(t) 5 Xmáx[cos(wt) 2 cos(wnt)] (1)
Al derivar 2 para x y x, se tiene:
x(t) 5 Xmáx[wn sen(wnt) 2 w sen(wt)] (2)
x(t) 5 Xmáx[wn2 cos(wnt) 2 w 2 cos(wt)] (3)
Al remplazar los valores en 1, 2 y 3 se tienen las ecuaciones del movimiento de la
señal en el tiempo:
x(t) 5 20,02414[cos(188.495 t) 2 cos(125.663 t)]
x(t) 5 23,0335 sen(125.663 t) 1 4,5502 sen(188.495 t)
x(t) 5 2381,1993 cos(125.663 t) 1 857,7081 cos(188.495 t)
Al obtener estas ecuaciones, se toman puntos en la señal de vibración. Para la am-
plitud de velocidad se promedian los datos con N 5 10 muestras por período (en
la tabla 6 se muestran los valores que se obtuvieron para cada estación de tiempo):
1 w
t5 →f5 5 30 Hz → t 5 0,033 s
f 2p
t 0,033 s
ti 5 5 5 0,0033 s
N 10
Tabla 6. Valores de la amplitud de la velocidad para N 5 10 muestras.
n ti(s) x i (mm/s)
1 0 0
2 0,0033 0,027443
3 0,0066 0,054882
4 0,0099 0,082311
5 0,0132 0,109726
6 0,0165 0,137123
7 0,0198 0,164496
8 0,0231 0,191841
9 0,0297 0,246429
10 0,033 0,273663
Fuente: Autores, 2018.
Capítulo 4: Casos de estudio: expansión fourier 55
Con los datos obtenidos se procede a analizar la expansión de la serie de Fourier:
a0 `
x(t) 5
2
1 ∑ (a cos nwt 1 b sin nwt )
n n
n51
Los datos de la tabla 7 se calculan a partir de los parámetros en la tabla 6, tenemos:
Tabla 7. Valores de la amplitud para cada muestra.
2npt i 2npt i
xi cos → xi sin
t t
n51 n52 n53
0 0 0 0 0 0
0,02220183 0,016130623 0,00848028 0,02609987 –0,00848047 0,02609981
0,016959317 0,052195934 –0,04440066 0,03225857 –0,04440019 –0,03225922
–0,025435843 0,078282301 –0,06659057 –0,04838178 0,06659164 –0,04838031
–0,088770578 0,064494803 0,03390843 –0,10435523 0,03390536 0,10435623
–0,137123 –1,00736E-06 0,137123 2,0147E-06 –0,137123 –3,0221E-06
–0,133079207 –0,096689496 0,0508293 0,15644589 0,0508362 –0,15644365
–0,059280253 –0,182452243 –0,15520495 0,11275812 0,15519915 0,1127661
0,199367164 –0,144844696 0,07615695 –0,23436589 –0,07614145 –0,23437093
0,273663 4,02088E-06 0,273663 8,0418E-06 0,273663 1,2063E-05
Fuente: Autores, 2018.
56
n54 n55 n56 n57 n58 n59 n 5 10
0 0 0 0 0 0 0 0 0 0 0 0 0 0
–0,02220195 0,01613046 –0,027443 –2,01608E-07 –0,02220171 –0,01613079 –0,00848008 –0,02609993 0,00848066 –0,02609974 0,02220207 –0,0161303 0,027443 4,03215E-07
0,01696008 –0,05219568 0,054882 8,06371E-07 0,01695855 0,05219618 –0,04440113 0,03225792 –0,04439971 –0,03225987 0,01696085 –0,05219544 0,054882 1,61274E-06
0,02543412 0,07828286 –0,082311 –1,81407E-06 0,02543757 –0,07828174 0,0665895 0,04838325 –0,0665927 0,04837884 –0,02543239 –0,07828342 0,082311 3,62814E-06
–0,08876868 –0,06449741 0,109726 3,22437E-06 –0,08877247 0,06449219 0,03391149 –0,10435423 0,03390229 0,10435722 –0,08876679 –0,06450002 0,109726 6,44874E-06
0,137123 4,0294E-06 –0,137123 –5,03681E-06 0,137123 6,0442E-06 –0,137123 –7,0515E-06 0,137123 8,0589E-06 –0,137123 –9,0663E-06 0,137123 1,00736E-05
–0,13308347 0,09668363 0,164496 7,25073E-06 –0,13307495 –0,09669536 0,05082241 0,15644813 0,05084309 –0,15644141 –0,13308773 0,09667776 0,164496 1,45015E-05
Análisis de vibraciones en maquinaria
0,05928964 –0,18244919 –0,191841 –9,8654E-06 0,05927087 0,18245529 0,15521075 –0,11275014 –0,15519335 –0,11277408 –0,05929902 0,18244614 0,191841 1,97308E-05
–0,19935759 –0,14485788 –0,246429 –1,62933E-05 –0,19937674 0,14483151 –0,07617244 0,23436086 0,07612596 0,23437596 0,19934801 0,14487106 0,246429 3,25866E-05
0,273663 1,6084E-05 0,273663 2,01044E-05 0,273663 2,4125E-05 0,273663 2,8146E-05 0,273663 3,2167E-05 0,273663 3,6188E-05 0,273663 4,02088E-05
Capítulo 4: Casos de estudio: expansión fourier 57
Se toman los datos de la tabla 7 y se suma, respectivamente, cada valor de n para
obtener los coeficientes an y bn, como sigue en la tabla 8:
Tabla 8. Valores de los coeficientes.
2 N 2npt i 2 N 2npt i
an 5 ∑
N i51
xi cos
t
bn 5 ∑
N i51
xi sen
t n n3w
0,013700486 20,042575952 1 188,495
0,06279295 20,01190608 2 376,99
0,06281005 20,04564459 3 565,485
0,01381163 20,05057662 4 753,98
20,016476 23,65065E-07 5 942,475
0,01380542 0,05057949 6 1.130,97
0,0628041 0,04565339 7 1.319,465
0,06279045 0,01191543 8 1.507,96
0,013693 0,04258258 9 1.696,455
0,2575828 2,58388E-05 10 1.884,95
Fuente: Autores, 2018.
Los términos de la serie de Fourier que corresponde a las funciones an cos(nwt)
y bn sen(nwt) se llaman armónicos de orden n de la función periódica x t. Estos
armónicos se trazan como líneas en un diagrama de amplitud contra la frecuencia
(nw), llamada espectro de frecuencia. Al graficar se obtiene (ver figuras 30 y 31):
Figura 30. Gráfico armónicos an del dominio de la frecuencia.
armónicos an
0,3
0,2
0,1
0
20,1 56 2 8 4 8 6 2 8 6 6
9 5 9 11 866 822 ,477 733 689 644 96,4 ,955
4 9
8, 76, 65, 4 9
3, 942 30, 9 , 4 , 9 16 884
18 3 5 75 19 507
11 13 1 1
armónicos an
Fuente: Autores, 2018.
58 Análisis de vibraciones en maquinaria
Figura 31. Gráfico armónicos bn del dominio de la frecuencia.
armónicos bn
0,06
0,04
0,02
0
20,02
20,04
6
20,06 4 8 36 92 48 ,46 955
56 12 68 22 77 73 68 64 9 6 ,
95 91 6 8 ,4 9 4
,48 53,9 942 30, 319, 507
,9 16 884
8,4 6,9 5
11
1
18 37 56 7 1 1
armónicos bn
Fuente: Autores, 2018.
A N E XO
Equipo utilizado para la medición y el análisis de vibración en máquinas rotativas,
en el laboratorio de la Universidad Francisco de Paula Santander Ocaña.
Imagen 4. Equipo de mediciones para prueba de vibración.
Fuente: Autores, 2018.
Imagen 5. Simulación de falla por desbalanceo dinámico.
Fuente: Autores, 2018.
60 Análisis de vibraciones en maquinaria
Imagen 6. Simulación de falla por desbalanceo estático.
Fuente: Autores, 2018.
BIBLIOGRAFÍA
Adamkowski, A., Henke, A. & Lewandowski, M. (2016). Resonance of torsional
vibrations of centrifugal pump shafts due to cavitation erosion of pump impe-
llers. Engineering Failure Analysis, 70, pp. 56-72.
Beer, F., Johnston, E., Mazurek, D. & Eisenberg, E. (2010). Mecánica vectorial para
ingenieros. México: Mc Graw Hill.
Crawford, J. (1977). Ondas. Barcelona: Editorial Reverté.
International Organization for Standardization. (1974). ISO [Link] Mecha
nical vibration of machines with operating speeds from 10 to 200 rev/s -- Basis
for specifying evaluation standards. Recuperado de [Link]
dard/[Link]
International Organization for Standardization. (2009). ISO 10816-3:2009(en):
Mechanical vibration – Evaluation of machine vibration by measurements on
non-rotating parts – Part 3: Industrial machines with nominal power above 15 kW
and nominal speeds between 120 r/min and 15.000 r/min when measured in situ.
Recuperado de [Link]
International Organization for Standardization. (2014). ISO 10816-[Link] Me
chanical vibration -- Evaluation of machine vibration by measurements on
non-rotating parts -- Part 8: Reciprocating compressor systems. Recuperado de
[Link]
Metrix. (s.f.). Vibration Severity Chart. Recuperado de [Link]
nanc [Link]/fileSendAction/fcType/0/fcOid/399590942963006681/
filePointer/399590942964827187/fodoid/399590942964827185/vibration_
severity_chart_and_other_guides.pdf
Robichaud, J. (2009). Reference Standards for Vibration Monitoring and Analysis.
Canadá: Bretech Engineering Ltd.
Roca Vila, R. y León, J. (1981). Vibraciones mecánicas. Ciudad de México: Limusa.
Serway, R. y Jewett, J. (2008). Física para ciencias e ingeniería. México D.F.: Cengage
Learning Editores.
Taylor, J. (1994). The Vibration Analysis Handbook. Estados Unidos: Vibration Con
sultants.
White, G. (2010). Introducción al análisis de vibraciones. Estados Unidos: Azima
DLI.
Este libro fue compuesto en caracteres Minion a 11
puntos, impreso sobre papel Bond de 75 gramos
y encuadernado con el método hot melt, en
diciembre de 2018, en Bogotá, Colombia.
Teoría de las vibraciones mecánicas
aplicadas en el diagnóstico de
maquinaria rotativa.
Introducción en el manejo de
equipos utilizados en la toma y
medición de vibraciones.
Formas de realizar la colección de
datos explicando los tipos de
espectros asociados a cada falla.
La teoría de las vibraciones, los distintos tipos de
vibraciones mecánicas con modelos matemáticos y
su aplicación en el análisis de maquinaria rotativa, Edwin Espinel Blanco
son el tema central de este libro. El texto presenta Ingeniero Mecánico, Especialista en Práctica
cada uno de los tipos de vibraciones, describe sus Docente Universitaria de la U. Francisco de Paula
Santander Ocaña, Especialista en Gerencia de
principios físicos y matemáticos, mientras señala los Mantenimiento Industrial y en Mantenimiento
límites definidos por normas ISO sobre los valores de Industrial de la EAFIT, Magister en Ingeniería de la
vibraciones aceptables en máquinas rotativas. Universidad EAFIT. Actualmente cursa su Doctora-
do en Ingeniería en la U. Pontificia Bolivariana.
Docente asistente vinculado al grupo de investi-
El libro está divido en cuatro capítulos, el primero, gación en Tecnología y Desarrollo en Ingeniería
aborda la concepción teórica sobre los distintos GITYD de la U. Francisco de Paula Santander
tipos de vibraciones, sus modelos matemáticos y sus Ocaña. Ha publicado diversos artículos en revistas
límites de vibración acorde a la normativa interna- científicas y ha sido ponente en diversos eventos
académicos a nivel nacional e internacional.
cional; el segundo, presenta la metodología para la
Andrés Camilo Melgarejo Ramos
toma de datos al momento de medir vibraciones en Ingeniero Mecánico egresado de la U. Francisco
máquina rotativa; el tercero, realiza una descripción de Paula Santander Ocaña. Co-Investigador
sobre el análisis de vibraciones y los distintos tipos adscrito al grupo de investigación en Tecnología
y Desarrollo ingeniería GITYD del Departamento
de espectros que corresponden a las fallas que se
de Ingeniería Mecánica de la misma institución.
presentan en la maquinaria; el cuarto, muestra Ha trabajo en proyectos de investigación relacio-
algunos casos de estudio sobre la aplicación de nados con el diseño y mantenimiento de sistemas
Fourier en la medición de vibraciones. mecánicos y el área de nuevos materiales.
Hugo Alberto Ramírez Galvis
Ingeniero Mecánico egresado de la U. Francisco
Dirigido a estudiantes, docentes e investigadores
de Paula Santander Ocaña. Co-Investigador
que deseen una aproximación al estudio de las adscrito al grupo de investigación en Tecnología
vibraciones mecánicas y a su aplicación en el estudio y Desarrollo en ingeniería GITYD del departamen-
y análisis de falla en maquinaria rotativa. Estudiantes to de ingeniería mecánica de la misma institución.
Ha trabajo en proyectos de investigación relacio-
de ingeniería mecánica, ingeniería de mantenimien-
nados con el diseño y mantenimiento de sistemas
to, o profesionales que trabajen en mantenimiento mecánicos y el área de nuevos materiales.
basado en condición y mantenimiento predictivo de
maquinaria.
ISBN 978-958-771-690-0
9 789587 716900