Analisis Informe 1234567
Analisis Informe 1234567
Analisis Informe 1234567
Laboratorio n°1:
“Vibraciones libres y forzadas en un
sistema de un grado de libertad”
El análisis de vibraciones es un área que todo ingeniero debería conocer, debido a que está
presente en casi todo sistema conocido, ya sea en estructuras o sistemas mecánicos como
motores o turbinas.
Debido a esto es que nace la necesidad de experimentar con estos fenómenos para así lograr un
entendimiento de modo que dichos conocimientos puedan ser usados a nuestro favor, o bien,
evitarlos en la medida de lo posible.
El resorte se deforma ∆ hasta llegar a un equilibrio estático igualándose la fuerza del resorte con el
peso de la masa 𝑘∆= 𝑊 = 𝑚𝑔
Σ𝐹 = 𝑚𝑥̈ = 𝑚𝑔 − 𝑘(∆ + 𝑥)
Como 𝑘∆= 𝑚𝑔, nos queda
̈ 𝑘𝑥 = 0
𝑚𝑥 + (1)
La frecuencia natural del sistema es:
𝑘
𝜔𝑛 2 =
𝑚
Escribiendo la ecuación (1) en términos de 𝜔𝑛 2 :
𝑥̈ + 𝜔𝑛 2 𝑥 =0 (2)
𝑥 = 𝐶1 𝑠𝑒𝑛 𝜔𝑛 𝑡 + 𝐶2 𝑐𝑜𝑠 𝜔𝑛 𝑡
Utilizando las condiciones iniciales 𝑥0 y 𝑥̇ 0 podemos encontrar las constantes 𝐶1 y 𝐶2 :
𝑥̇ 0
𝑥= 𝑠𝑒𝑛 𝜔𝑛 𝑡 + 𝑥0 𝑐𝑜𝑠 𝜔𝑛 𝑡 (2.1)
𝜔𝑛
Vibraciones libres con amortiguador viscoso
𝑥(𝑡) = 𝐶𝑒 𝑟𝑡
donde C y r son constantes indeterminadas. La inserción de esta función en la ecuación (3) nos lleva
a la ecuación característica
𝑚𝑟 2 + 𝑐𝑟 + 𝑘 = 0
Cuyas raíces son
−𝑐 ± √𝑐 2 − 4𝑚𝑘 −𝑐 𝑐 2 𝑘
𝑟1,2 = = ± √( ) − (4)
2𝑚 2𝑚 2𝑚 𝑚
𝑥1 (𝑡) = 𝐶1 𝑒 𝑟1 𝑡 𝑦 𝑥2 (𝑡) = 𝐶2 𝑒 𝑟2 𝑡
Por lo tanto, la solución general de la ecuación (3) es:
𝑥 (𝑡) = 𝐶1 𝑒 𝑟1 𝑡 + 𝐶2 𝑒 𝑟2 𝑡
Con 𝐶1 𝑦 𝐶2 constantes arbitrarias que se determinan con las condiciones iniciales
Amortiguamiento critico 𝑐𝐶 : se define como el valor de la constante de amortiguamiento c con la
cual el radical de la ecuación (4) se vuelve cero.
𝑐𝑐 2 𝑘 𝑘
( ) − =0 → 𝑐𝑐 = 2𝑚√ = 2𝑚𝜔𝑛 = 2√𝑘𝑚
2𝑚 𝑚 𝑚
Factor de amortiguamiento:
𝑐
𝜉= (5)
𝑐𝑐
𝑐
𝜉=
2𝑚𝜔𝑛
Escribiendo la ecuación (3) en términos de 𝜁 y 𝜔𝑛 nos queda:
𝑥̇ 0 + 𝜉𝜔𝑛 𝑥0
𝑥 = 𝑒 −𝜉𝜔𝑛 𝑡 (𝑥0 cos √1 − 𝜉 2 𝜔𝑛 𝑡 + 𝑠𝑒𝑛√1 − 𝜉 2 𝜔𝑛 𝑡)
𝜔𝑛 √1 − 𝜉 2
𝑥̇ 0 + (𝜉 + √𝜉 2 − 1)𝜁𝜔𝑛 𝑥0
𝐶1 =
2𝜔𝑛 √𝜉 2 − 1
−𝑥̇0 − (𝜉 − √𝜉 2 − 1)𝜁𝜔𝑛 𝑥0
𝐶2 =
2𝜔𝑛 √𝜉 2 − 1
Decremento logarítmico
El decremento logarítmico hace referencia al logaritmo natural de la relación de dos amplitudes
sucesivas de una vibración libre con amortiguamiento.
Fig 3. Decremento Logarítmico, con 𝑋𝑛 . 𝑋𝑛+1 dos amplitudes consecutivas y 𝑇𝑑 el periodo del
sistema con amortiguamiento.
Donde:
𝑋𝑛
𝛿 = 𝑙𝑛 = 𝜉𝜔𝑛 𝑇𝑑
𝑋𝑛+1
2𝜋
𝑇𝑑 =
𝜔𝑛 √1 − 𝜉 2
Sustituyendo 𝑇𝑑 en la relación:
2𝜋𝜉
𝛿= (7)
√1 − 𝜉 2
𝑥𝑝 (𝑡) = 𝑋0 𝑆𝑖𝑛(Ω𝑡 − 𝜙)
Con 𝑋0 y 𝜙 iguales a:
𝐹0 𝑐Ω
𝑋0 = 𝑘 ; 𝑇𝑎𝑛𝜙 = 𝑘
2 𝑀Ω2
𝑀Ω2 𝑐Ω 2 1−
√(1 − ) +( ) 𝑘
𝑘 𝑘
𝑘
Y sabiendo que √𝑀 = 𝜔𝑛 , se obtiene finalmente:
Ω
𝑋0 1 2𝜉 (𝜔 )
𝑛
= ; 𝜙 = 𝑡𝑎𝑛−1
𝐹0 2 2 2 Ω 2
1−( )
𝑘 √[1 − ( Ω ) ] + [2𝜉 ( Ω )] [ 𝜔𝑛 ]
𝜔 𝑛 𝜔 𝑛
Fig 5.Curva de amplitud y fase en una respuesta forzada.
1 𝑚𝑟 𝑦𝑥̇ 2
𝑑𝑇 = ( 𝑑𝑦) ( )
2 𝑙 𝑙
Energía cinética total del sistema
𝑙
1 1 𝑚𝑟 𝑦𝑥̇
𝑇= 𝑚𝑥̇ 2 + ∫ ( 𝑑𝑦) ( )
2 0 2 𝑙 𝑙
1 1 𝑚𝑟 2
𝑇= 𝑚𝑥̇ 2 + 𝑥̇
2 2 3
Por lo tanto, la masa equivalente del sistema es
𝑚𝑟
𝑚𝑒𝑞 = (𝑚 + )
3
Energía potencial total del sistema
1
𝑈= 𝑘𝑥 2 → 𝑘𝑒𝑞 = (𝑘)
2
3. Objetivos
Objetivo general:
Objetivo general: Analizar las vibraciones libres y forzadas de un sistema viga – resorte.
La viga se considera rígida y pivoteada en uno de sus extremos (sin restricción a la
rotación).
Objetivos específicos:
Comenzando el laboratorio se tomaron las medidas de las distancias desde el agarre hasta
diferentes puntos de interés a lo largo de la viga, estos puntos son el motor, amortiguador, sensor
de desplazamiento y el resorte.
Luego se procede a elegir 3 masas con las cuales se realiza la experiencia, estas se ubicaron a la
altura de el sensor. Para comenzar el movimiento, con la mano se sostiene la viga manteniéndola
en una posición inicial así al soltarla la viga comienza a vibrar. Es banco de ensayo (Img.1) está
conectado mediante un sistema versátil de adquisición de datos (VDAS, Img.2). se realizo tres
veces el procedimiento, primero no se utilizó el amortiguador, segundo se insertó el amortiguador
con agua, y por ultimo se cambio el liquido del amortiguador por silicona, esta vez se utilizo los
orificios del amortiguador abiertos ya que al ocuparlo cerrado la viga vuelve demasiado rápido al
equilibrio y dejaba muy pocos datos para analizar.
Img. 1 Banco de ensayo utilizado en el laboratorio Img. 2 sistema versátil de adquisición de datos
La segunda parte del laboratorio consiste en vibraciones forzadas, para esto se procede a
encender el motor el cual tiene una masa a 35mm del centro la cual provoca una excitación y así al
girar produce vibraciones en la viga, como el movimiento de la masa en el motor es rotatorio, esto
produce una fuerza senoidal. Para esta parte del laboratorio también se realizo 3 veces cambiando
el amortiguamiento de la misma forma anterior.
4. RESULTADOS
Asociación del sistema físico a un modelo de un grado de libertad:
1 2 1 𝑳𝒗 2 1 2 1 𝑴𝑘 2 1 2 1 2
𝑴𝒎 (𝑳𝒎 𝜃̇ ) + 𝑴𝒗 ( 𝜃̇ ) + 𝑴𝑐 (𝑳𝒄 𝜃̇ ) + (𝑳𝒌 𝜃̇ ) + 𝑴𝐴𝑘 (𝑳𝒌 𝜃̇ ) = 𝑱𝑬𝑸 𝜃̇
2 2 𝟐 2 2 3 2 2
1 𝑳𝒗 2 𝑴𝑘 2 2 1 2
⌊𝑴𝒎 𝑳𝒎 2 + 𝑴𝒗 + 𝑴 𝑐 𝑳𝒄 2 + 𝑳𝒌 + 𝑴𝐴𝑘 𝑳𝒌 2 ⌋ 𝜃̇ = 𝑱𝑬𝑸 𝜃̇
2 4 3 2
Luego por semejanza obtenemos que la inercia rotacional equivalente es:
𝑳𝒗 2 𝑴𝑘
𝑴 𝒎 𝑳𝒎 2 + 𝑴 𝒗 + 𝑴 𝑐 𝑳𝒄 2 + 𝑳𝒌 2 + 𝑴𝐴𝑘 𝑳𝒌 2 = 𝑱𝑬𝑸
4 3
Reemplazando los valores:
0.8152 0.388
4.500 ∗ 0.3802 + 1.650 ∗ + 0.530 ∗ 0.6502 + ∗ 0.7552 + 0.090 ∗ 0.7552 = 𝑱𝑬𝑸
4 3
𝑱𝑬𝑸 = 𝟏. 𝟐𝟕𝟑 𝒌𝒈 ∗ 𝒎𝟐
M1=0.2kg
M2=0.5kg
M3=1kg
Fig 7.
Como se puede observar, los 3 datos crean una recta casi perfecta entre sí. Los círculos
representan los 3 puntos tabulados en la tabla 2 cuando se grafican en la función descrita en la
ecuación (9).
Al realizar un ajuste lineal a los datos graficados, se obtiene la pendiente de la gráfica.
Fig 8.
Obteniendo una pendiente de la recta de 3.1113, es decir K=3.1113 N/mm. La recta ajustada no
pasa por el origen, la causa de esto se discutirá en el apartado discusión y análisis de resultados.
Sabiendo que 𝑥 = 𝑟𝑆𝑒𝑛𝜃, pero para valores pequeños de 𝜃, aproximamos 𝑆𝑒𝑛𝜃 = 𝜃, resultando
que 𝑥 = 𝑟𝜃, despejando para 𝜃, y sabiendo que 𝑟 = 𝐿𝑘 obtenemos que:
𝑥
𝜃=
𝐿𝑘
Reemplazando en (n):
𝑥
𝑘 𝑇𝐸𝑄 = 𝑘𝐸𝑄 𝑥𝐿𝑘
𝐿𝑘
Simplificando, se puede obtener la rigidez torsional equivalente:
𝑘 𝑇𝐸𝑄 = 𝑘𝐸𝑄 𝐿𝑘 2
𝒌𝑻𝑬𝑸 = 𝟏𝟕𝟕𝟑. 𝟓𝟐 𝑵𝒎
Teniendo los valores de 𝐽𝐸𝑄 y 𝑘 𝑇𝐸𝑄 podemos obtener el valor de 𝜔𝑛 con la fórmula:
𝑘 𝑇𝐸𝑄
𝜔𝑛 = √
𝐽𝐸𝑄
1773.52
𝜔𝑛 = √
1.273
Finalmente:
𝒓𝒂𝒅
𝝎𝒏 = 𝟑𝟕. 𝟑𝟐𝟓
𝒔
Con los valores de 𝐽𝐸𝑄 y 𝜔𝑛 se puede obtener el valor de la constante de amortiguamiento critico
𝑐𝑐 :
𝑐𝑐 = 2𝐽𝐸𝑄 𝜔𝑛
𝑐𝑐 = 2 ∗ 1.273 ∗ 37.325
𝒄𝒄 = 𝟗𝟓. 𝟎𝟑 𝑵𝒎𝒔
Fig 9.
Forma general:
1 𝑥𝑛 2𝜋𝜁 𝛿2
→ 𝛿= ln ( ) → 𝛿= → 𝜉= √
𝑁 𝑥𝑛+𝑁 √1 − 𝜉 2 (2𝜋)2 + 𝛿 2
Con los valores del factor de amortiguamiento 𝜉𝑎𝑔𝑢𝑎 y 𝜉𝑠𝑖𝑙𝑖𝑐𝑜𝑛𝑎 el amortiguamiento critico
𝑐𝑐 podemos calcular los valores de los amortiguamientos mediante la ecuación (5):
𝑐 = 𝑐𝑐 𝜉
Reemplazando valores:
𝜔𝑑 = 𝜔𝑛 √(1 − 𝜉 2 )
𝜔𝑑 𝑎𝑔𝑢𝑎 = 𝜔𝑛 √(1 − 𝜉𝑎𝑔𝑢𝑎 2 ) 𝜔𝑑 𝑠𝑖𝑙𝑖𝑐𝑜𝑛𝑎 = 𝜔𝑛 √(1 − 𝜉𝑠𝑖𝑙𝑖𝑐𝑜𝑛𝑎 2 )
𝜔𝑑 𝑎𝑔𝑢𝑎 = 37.325 √(1 − 0.03532 ) 𝜔𝑑 𝑠𝑖𝑙𝑖𝑐𝑜𝑛𝑎 = 37.325 √(1 − 0.05112 )
𝒓𝒂𝒅 𝒓𝒂𝒅
𝝎𝒅 𝒂𝒈𝒖𝒂 = 𝟑𝟕. 𝟑𝟎𝟏 𝝎𝒅 𝒔𝒊𝒍𝒊𝒄𝒐𝒏𝒂 = 𝟑𝟕. 𝟐𝟕𝟔
𝒔 𝒔
Calculados los valores de 𝐽𝐸𝑄 , 𝐾𝑇𝐸𝑄 , 𝑐𝑎𝑔𝑢𝑎 y 𝑐𝑠𝑖𝑙𝑖𝑐𝑜𝑛𝑎 podemos proceder a asociar el sistema físico a
un modelo de un grado de libertad,
𝑥̇ 0 + 𝜉𝑠𝑖𝑙𝑖𝑐𝑜𝑛𝑎 𝜔𝑛 𝑥0
𝑥(𝑡) = 𝑒 −𝜉𝑠𝑖𝑙𝑖𝑐𝑜𝑛𝑎 𝜔𝑛 𝑡 (𝑥0 Cos 𝜔𝑑𝑠𝑖𝑙𝑖𝑐𝑜𝑛𝑎 𝑡 + 𝑆𝑒𝑛𝜔𝑑𝑠𝑖𝑙𝑖𝑐𝑜𝑛𝑎 𝑡)
𝜔𝑑𝑠𝑖𝑙𝑖𝑐𝑜𝑛𝑎
Al graficar los datos obtenidos en el laboratorio para las vibraciones forzadas obtenemos:
Fig 10.
Fig 11.
Fig 12.
Sin embargo, al realizar los gráficos de Bodé para las vibraciones forzadas:
Fig 13.
5. Discusión, análisis de resultados y conclusiones
Fig 14.
Fig 15.
Fig 16.
La segunda parte del laboratorio se pudo apreciar que al encender el motor e ir subiendo
la frecuencia llegaba a un punto donde la frecuencia rondaba entre 6 y 7 Hz donde se
produce el fenómeno de “resonancia”. En este instante la amplitud del sistema creció
notoriamente (como se puede apreciar en las Fig 12. 13. 14), el motivo por que creció la
amplitud es porque la relación entre la frecuencia del motor y la frecuencia natural del
Ω
sistema es cercana a 1 ( ≈ 1), después de cierto tiempo la frecuencia fue interrumpid
𝜔𝑛
6. Bibliografía