Métodos de Balanceo
Métodos de Balanceo
ISBN: 978-607-95309-5-2 << pag. 915 >> Derechos Reservados © 2011, SOMIM
MEMORIAS DEL XVII CONGRESO INTERNACIONAL ANUAL DE LA SOMIM
21 al 23 DE SEPTIEMBRE, 2011 SAN LUIS POTOSÍ, MÉXICO
A4_63
2
casos, los puntos solución delimitan un área Figura No. 2 Función de ajuste
común formada por tres vértices.
Las causas de error que producen la no El máximo de la función de ajuste
intersección de los círculos en un punto común mostrada en la figura número 2 representa la suma
impiden la rápida convergencia de la vibración a vectorial de dos vectores con la misma dirección
valores seguros de operación. Tratándose de un (las fuerzas de desbalance inicial y de prueba)
método que utiliza muchas corridas de prueba, el mientras que el mínimo de la función representa la
tiempo utilizado para la realización del balanceo suma vectorial de dos vectores con direcciones
cuando éste tiene una lenta convergencia es mayor opuestas entre sí.
debido a dichos errores. De acuerdo a lo anterior, el ángulo de
fase de la lectura de vibración en condiciones
Aproximación numérica de los ángulos de fase iniciales es aproximadamente el mismo que el
perteneciente al valor máximo de la función.
El objetivo del ajuste numérico de
El procedimiento para la realización del
funciones es obtener, mediante la utilización de
balanceo dinámico en un plano, utilizando este
una cantidad finita de puntos o datos, expresiones
algoritmo de ajuste, requiere de la utilización del
analíticas que representen aproximaciones al
método de coeficientes de influencia.
comportamiento real de las variables de estudio.
El método de coeficientes de influencia
El modelado del comportamiento de la amplitud
para el balanceo dinámico de rotores rígidos en un
de vibración causada por desbalances se realiza
plano requiere de una lectura de vibración en
proponiendo la utilización de una combinación
condiciones iniciales y una lectura con un peso de
lineal de funciones no lineales de la siguiente
prueba conocido colocado en el plano a balancear
forma:
[3].
𝑓 𝑥 = 𝑎0 + 𝑎1 cos 𝜔𝑥 + 𝑏1 sen 𝜔𝑥 (10)
Tabla No. 1 Información para el balanceo en
un plano
donde: 𝑎0 , 𝑎1 y 𝑏1 son constantes por determinar y
𝜔 es la frecuencia circular de giro de la función. Magnitud Ángulo de fase
Condiciones
Utilizando el método de ajuste de X0 δ0
funciones por mínimos cuadrados es posible iniciales
Peso de
calcular los valores de las constantes en función Wpi βpi
de los datos obtenidos en las tres corridas de prueba
prueba [2]. Las lecturas de vibración Con peso de
Xi γi
seleccionadas para realizar el ajuste corresponden prueba
a pruebas realizadas con una separación espacial
de ciento veinte grados entre sí dentro del plano Utilizando la información anterior se
de pruebas, con el fin de abarcar la totalidad de los calcula el coeficiente de influencia y el peso de
posibles comportamientos en la amplitud de corrección de acuerdo con las siguientes fórmulas:
vibración.
𝑋𝑖 ∠𝛾𝑖 − 𝑋0 ∠𝛿0
𝐴= (11)
𝑊𝑝𝑖 ∠𝛽𝑝𝑖
ISBN: 978-607-95309-5-2 << pag. 916 >> Derechos Reservados © 2011, SOMIM
MEMORIAS DEL XVII CONGRESO INTERNACIONAL ANUAL DE LA SOMIM
21 al 23 DE SEPTIEMBRE, 2011 SAN LUIS POTOSÍ, MÉXICO
A4_63
−𝑋0 ∠𝛿0 valor obtenido de manera experimental (con
𝑊𝑐 = (12)
𝐴 errores) 𝑋2 ;
ISBN: 978-607-95309-5-2 << pag. 917 >> Derechos Reservados © 2011, SOMIM
MEMORIAS DEL XVII CONGRESO INTERNACIONAL ANUAL DE LA SOMIM
21 al 23 DE SEPTIEMBRE, 2011 SAN LUIS POTOSÍ, MÉXICO
A4_63
2 2
𝑊𝑝 𝑋𝑟
𝑘= 2
(21) −𝑋0 𝑊𝑝
𝑋0 𝑊𝑐 =
𝑋2 (28)
− 𝑋0
𝑝
donde k es un factor de proporcionalidad. La
ecuación (20) representa la relación de
proporcionalidad entre dos conjuntos de fasores: −𝑝𝑋0 𝑊𝑝
𝑊𝑐 = (29)
los fasores representativos de los pesos colocados 𝑋2 − 𝑝𝑋0
en el plano de balanceo y los fasores
representativos de la desviación de los valores El fasor p, que relaciona valores reales
experimentales respecto de los valores reales. Se con valores con errores, puede utilizarse de tal
define: manera que afecte al fasor 𝑋0 , transformándose
éste en un valor erróneo 𝑋0 . De manera análoga a
𝑊𝜖 𝑒 𝑖𝜃 𝑊 𝜖 = 𝑊𝑐 𝑒 𝑖𝜃𝑐 − 𝑊𝑝 𝑒 𝑖𝜃𝑝 (22) la ecuación (14), se define:
𝜖 2 𝑋0 = 𝑝 ∗ 𝑋0 (30)
𝑘 (23)
= 2
𝑊𝑐 2 + 𝑊𝑝 − 2 𝑊𝑐 𝑊𝑝 cos 𝜃𝑐 − 𝜃𝑝 Igualando (14) y (30) se obtiene la siguiente
relación:
𝜃𝜖 = −𝜃𝑊𝜖 (24) 𝑋0 𝑋2
𝑝= = (31)
𝑋0 𝑋2
𝜃𝓊 = 𝜃𝑐 − 𝜋 (25)
Reescribiendo:
El ángulo del fasor unitario es tal que, al
calcular el valor de p, éste represente la dirección 𝑋0 𝑋2 𝑒 𝑖(𝜃0 +𝜃2 ) = 𝑋0 𝑋2 𝑒 𝑖(𝜃0 +𝜃2 ) (32)
de referencia sobre la cual se presentan las
desviaciones producidas por posibles errores.
De acuerdo a (32) se pueden manejar de
Realizando las consideraciones anteriores es
manera independiente las magnitudes y los
posible obtener al fasor p, el cual se calcula de la
ángulos de fase para la realización de las
siguiente manera:
perturbaciones. Sean ε y δ pequeñas
perturbaciones realizadas a la amplitud y al ángulo
𝑝 𝑒 𝑖𝜃 = 1 𝑒 𝑖𝜃𝓊 − 𝜖 2 𝑒 𝑖𝜃𝜖 (26)
de fase respectivamente, se definen:
El peso de corrección corregido, 𝑋0 = 𝑋0 1 + 𝜀 (33)
considerando al fasor p se calcula con la siguiente
𝑋2 = 𝑋2 1 − 𝜀 (34)
fórmula:
𝜃0 = 𝜃0 1 + 𝛿 (35)
𝑊𝑐 𝜃2 = 𝜃2 1 − 𝛿 (36)
𝑊𝑐 =
𝜖 𝑊𝑐 − 𝑊𝑝 (27)
1+ Aplicando las perturbaciones a (32) se obtiene:
𝑝 𝑊𝑐
2𝛿 2
𝑖 (37)
El nuevo peso de corrección representa 𝑋0 𝑋2 = 𝑋0 𝑋2 1 − 𝜀 2 𝑒 1−𝛿 2
una aproximación válida del valor obtenido de
manera convencional; sin embargo, la dirección Despejando para los dos casos posibles:
del fasor unitario es referida sobre la dirección del
peso de corrección con errores, por tanto se realiza 2𝛿 2
𝑖 (38)
un procedimiento de perturbaciones sobre la 𝑋0 = 𝑋0 ∗ 𝑝 1 − 𝜀 2 𝑒 1−𝛿 2
amplitud y la fase del fasor p para mejorar la 𝑝
exactitud del peso de corrección. 𝑋2 = 𝑋2 ∗ 2𝛿 2 (39)
𝑖
Un peso de corrección sin errores, 1 − 𝜀 2 𝑒 1−𝛿 2
calculado con las ecuaciones (11) y (12), puede
expresarse de las siguientes maneras, al considerar Las ecuaciones anteriores permiten
la ecuación (14): aplicar las perturbaciones del fasor p
ISBN: 978-607-95309-5-2 << pag. 918 >> Derechos Reservados © 2011, SOMIM
MEMORIAS DEL XVII CONGRESO INTERNACIONAL ANUAL DE LA SOMIM
21 al 23 DE SEPTIEMBRE, 2011 SAN LUIS POTOSÍ, MÉXICO
A4_63
indistintamente a las lecturas en condiciones Tabla No. 2 Lecturas de vibración antes de las
iniciales o con peso de prueba. pruebas
Para discriminar cual de los dos valores Vibración residual inicial (mm/s)
perturbados de p producirá la convergencia del Apoyo 1 0.7366
peso de corrección calculado al valor real del Apoyo 2 1.016
mismo se propone calcular un valor de vibración
residual producida por la colocación de un peso de Utilizando el método analítico para el
corrección con errores 𝑊𝑐 (el peso calculado antes balanceo en un plano se procede a la realización
de la perturbación) y el nuevo peso de corrección del balanceo.
𝑊𝑐 obtenido con la ecuación (27) al utilizar el Reduciendo el error en los datos de la
fasor p obtenido después de las perturbaciones. tabla 3 con el método propuesto en la sección
Esto se muestra en la ecuación (40). anterior, se obtienen los pesos de corrección de la
tabla 4.
ISBN: 978-607-95309-5-2 << pag. 919 >> Derechos Reservados © 2011, SOMIM
MEMORIAS DEL XVII CONGRESO INTERNACIONAL ANUAL DE LA SOMIM
21 al 23 DE SEPTIEMBRE, 2011 SAN LUIS POTOSÍ, MÉXICO
A4_63
Reduciendo el error en los datos de la
Tabla No. 5 Balanceo utilizando el método de tabla 7 con el método propuesto se obtiene:
aproximaciones
Amplitud de vibración Tabla No. 8 Pesos de corrección ajustados
(mm/s) 𝑊𝑐
Condiciones iniciales Wp1= 250º gr. 𝑊𝑐1 = 25.455 ∠233.07°
X0=20.2438
(𝑊𝑑 = 20𝑔𝑟. ∠200°) Wp2= 25120º gr. 𝑊𝑐2 = 31.3 ∠227.97°
Wp1= 250º gr. X1 =8.8138 Wp3= 25240º gr. 𝑊𝑐3 = 29.73 ∠224.61°
Wp2= 25120º gr. X2 =33.5534
Wp3= 25240º gr. X3 =41.2242 Utilizando 𝑊𝑐3 , por ser este el peso
calculado con el peso de prueba que disminuye la
𝑊𝑐 = 17.7 ∠14.4° gr. Xr = 2.4892 vibración de acuerdo a la tabla 7, se obtiene una
vibración residual de 3.8862 mm/s.
El peso de corrección original reduce la
Reduciendo el error en los datos de la
vibración inicial en un 43%, contra un 88% de
tabla 5 con el método propuesto se obtiene:
reducción al utilizar el peso de corrección
ajustado.
Tabla No. 6 Pesos de corrección ajustados
𝑊𝑐
CONCLUSIONES
Wp1= 250º gr. 𝑊𝑐1 = 19.5798 ∠17.3317°
Wp2= 25120º gr. 𝑊𝑐2 = 18.1322 ∠15.1718° El análisis desarrollado en el presente
Wp3= 25240º gr. 𝑊𝑐3 = 18.08 ∠15.0793° trabajo permite una descripción cuantitativa de los
errores introducidos al peso de corrección debidos
principalmente a inexactitudes en las lecturas de
Utilizando 𝑊𝑐1 , se obtiene una vibración
amplitud de vibración. El método de ajuste
residual de 1.1684 mm/s, consiguiendo una desarrollado permite utilizar la información
reducción del 94% del valor inicial, contra un obtenida por algoritmos de balanceo sin fase,
87% de la aproximación inicial del peso de centrándose principalmente en el comportamiento
corrección. de los pesos de prueba y la amplitud de vibración
Los balanceos mostrados con residual para mejorar el peso de corrección que
anterioridad fueron realizados en condiciones disminuirá los niveles de vibración hacia cero.
Para tal efecto se vuelve indispensable una buena
controladas de desbalance. El siguiente
selección de la magnitud y la posición de los
experimento se realizó en condiciones
pesos de prueba a utilizar durante el balanceo.
desconocidas de desbalance utilizando el método
Los experimentos realizados permiten aseverar
por aproximaciones numéricas. que el procedimiento de mejora de datos es
satisfactorio, sin importar el algoritmo de
Tabla No. 7 Balanceo utilizando el método de balanceo utilizado, al mostrar valores de vibración
aproximaciones residual de hasta un 45% menos, al usar los pesos
Amplitud de vibración de corrección ajustados.
(mm/s)
Condiciones REFERENCIAS
X0=33.274
iniciales
[1]. García Reynoso, A.C.; “Método de balanceo
Wp1= 200º gr. X1 =48.26
dinámico de rotores rígidos sin la medición del
Wp2= 20120º gr. X2 =42.4688 ángulo de fase de la vibración”, XIII Congreso
Wp3= 20240º gr. X3 =15.0622 internacional anual de la SOMIM, Durango,
Dgo., septiembre 2007.
[2]. Nakamura, Shoichiro; “Métodos numéricos
𝑊𝑐 =
Xr = 18.9738 aplicados con software”, 1992, Pearson
13.85 ∠230.4° gr. Prentice-Hall.
ISBN: 978-607-95309-5-2 << pag. 920 >> Derechos Reservados © 2011, SOMIM
MEMORIAS DEL XVII CONGRESO INTERNACIONAL ANUAL DE LA SOMIM
21 al 23 DE SEPTIEMBRE, 2011 SAN LUIS POTOSÍ, MÉXICO
A4_63
[3]. Den Hartog, Jacob P., "Mechanical
Vibrations", 4th ed. Reprint, DOVER
PUBLICATIONS INC., 1985
[4]. García Reynoso, A.; “Análisis del error y
convergencia optima del método de
coeficientes de influencia para el balanceo
dinámico de rotores”, Memoria en extenso,
segundo coloquio del IPN, 1989.
[5]. E. Ladrón de Guevara D., A.C. García R., E.
Hernández M., B. Peña P., J. Camargo M., U.
Uscanga B., “Efecto de la captura automática
de datos de vibración en el balanceo dinámico
de rotores”, XIII Congreso Internacional
Anual de SOMIM, Durango, Dgo.,
Septiembre, 2007.
ISBN: 978-607-95309-5-2 << pag. 921 >> Derechos Reservados © 2011, SOMIM