Guia de Laboratorio 2 Vibracion Amortiguada
Guia de Laboratorio 2 Vibracion Amortiguada
Guia de Laboratorio 2 Vibracion Amortiguada
x(t)
❖ OBJETIVOS
General: Identificar la vibración resultante de un sistema mecánico amortiguado
cuando se modifica la constante de amortiguación, utilizando el programa Matlab
simulink u otro programa.
Específicos:
1- Establecer la expresión matemática para el análisis, que represente el
modelo físico del sistema.
2- Identificar los operadores y su representación del programa de simulink o
de otro programa
3- Aplicar las variables de entrada al sistema las variables de salida.
4- Verificar las señales de entrada y de salida de un sistema vibrante.
pág. 1
❖ MARCO TEÓRICO
Definición.
El proceso por el cual la vibración disminuye continuamente de amplitud se le
llama amortiguación. En la amortiguación, la energía del sistema se disipa en
calor, sonido o cualquier combinación de ellas.
pág. 2
LABORATORIO 2 VIBRACIÓN AMORTIGUADA
𝑋 = 𝐷 𝑒 𝑠𝑡 (2)
Donde D y S son constantes. Las derivadas son:
Ẋ = 𝐷𝑆 𝑒 𝑠𝑡
Ẍ = 𝐷𝑆 2 𝑒 𝑠𝑡
Sustituyendo en la ecuación 1:
𝑚 𝐷𝑆 2 𝑒 𝑠𝑡 + 𝐶𝐷𝑆 𝑒 𝑠𝑡 + 𝐾𝐷 𝑒 𝑠𝑡 = 0
𝑚 𝑆 2 + 𝐶𝑆 + 𝐾 = 0
𝐶 𝐾
𝑆 2 + ( )𝑆 + =0
𝑚 𝑚
Determinando las raíces de la ecuación de segundo grado:
𝐶 𝐶 𝐾 1/2
𝑆1,2 = − 2𝑚 ± [(2𝑚) − ] (3)
𝑚
Por lo tanto, con dos valores de S la solución general está dada por la ecuación
2 se modifica:
𝑋 = 𝑎 𝑒 𝑆1 𝑡 + 𝑏 𝑒 𝑆2 𝑡 (4)
Tipos de amortiguamiento.
Para facilitar el manejo de la función y tomando como parámetro la ec. 3, se
utilizará el termino factor de amortiguamiento ζ.
𝐶
𝐶
𝜁 = 2𝑚 =
𝐾
( 𝑚)1/2 2𝑚𝜔𝑛
La ecuación 3 quedará:
𝑆1,2 = − 𝜁 𝜔𝑛 ± 𝜔𝑛 ( 𝜁 2 − 1 )1/2 (5)
a) Amortiguamiento crítico ζ = 1.
𝐶
𝜁=1=
2 𝑚 𝜔𝑛
𝐶𝑐 = 2 𝑚 𝜔𝑛 = ( 4 𝐾 𝑚)1/2
Cc es el amortiguamiento crítico. Hace que el valor ζ sea 1. Los valores de S1 y
S2 serán menores que cero.
𝑆1 = 𝑆2 = − 𝜔𝑛
La ecuación 4 quedará:
𝑋 = (𝑎 + 𝑏 𝑡) 𝑒 − 𝜔𝑛 𝑡 (6)
pág. 3
LABORATORIO 2 VIBRACIÓN AMORTIGUADA
Fig. 2. Amortiguamiento crítico.
b) Sobre amortiguado ζ > 1.
𝐶𝑅 > (4 𝑚 𝐾)1/2 = 2 𝑚 𝜔𝑛
𝐶𝑅 > 𝐶𝑐
Fig. 3. Sobreamortiguado.
c) Sub amortiguado ζ < 1.
𝐶 < 2 𝑚 𝜔𝑛 = (4 𝑚 𝐾)1/2
𝐶𝑆𝑅 < 𝐶𝑐
𝜔𝑑 = 𝜔𝑛 (1 − 𝜁 2 )1/2
2𝜋 1
𝑇𝑑 = =
𝜔𝑑 𝑓𝑑
ωd es la frecuencia de la amortiguación.
Td es el período de la amortiguación.
fd es la frecuencia de la amortiguación.
pág. 4
LABORATORIO 2 VIBRACIÓN AMORTIGUADA
ωd será siempre menor que ωn porque ζ oscila entre mayor que 0 y menor que
1.
La ecuación 4 quedará:
𝑋 = 𝑎 𝑒 𝑆1 𝑡 + 𝑏 𝑒 𝑆2 𝑡
Al utilizar Euler, sustituyendo y desarrollando, la ecuación resulta así:
Todas las variables se determinan por las condiciones iniciales. Se concluye que
el desplazamiento tiende a cero cuándos t tiende al infinito.
❖ METODOLOGIA DE SOLUCION
pág. 5
LABORATORIO 2 VIBRACIÓN AMORTIGUADA
6- Aplicar la Transformada de Laplace a la ecuación diferencial resultante
del sistema.
7- Elaborar el diagrama de bloques en Simulink.
8- Verificar la señal resultante para observar el tipo de vibración que
experimenta el cuerpo según sea el caso de amortiguamiento.
❖ EJEMPLO
Enunciado
Un bloque de 50 N con las condiciones iniciales en t = 0 se encuentra en X (0)
= 10 mm por encima de su punto de equilibro y se suelta hacia abajo con una
velocidad Ẋ (0) = 350 mm/s. Determine la ecuación diferencial que rige el
movimiento, determine la frecuencia natural, determine las tres ecuaciones de
desplazamiento para cada caso de amortiguamiento, además realice el diagrama
de bloques con su respectiva simulación de los tres casos mediante un diagrama
de operadores en Simulink.
Desarrollo
pág. 6
LABORATORIO 2 VIBRACIÓN AMORTIGUADA
1000 N/m (x (t)
)
1333 N/m (x (t))
+
x (t)
C 𝑥̇ (𝑡)
a) Amortiguamiento crítico ζ = 1.
𝐶𝑐
𝜁=1=
2 𝑚 𝜔𝑛
𝐶𝑐 = 2 𝑚 𝜔𝑛
𝐶𝑐 = 2 (5.1 𝑘𝑔)(21.38 𝑟𝑎𝑑⁄𝑠) = 218.07
𝑥(𝑡) = (𝑎 + 𝑏 𝑡) 𝑒 − 𝜔𝑛 𝑡
Sustituyendo valores.
pág. 7
LABORATORIO 2 VIBRACIÓN AMORTIGUADA
𝑥(𝑡) = (𝑎 + 𝑏 𝑡) 𝑒 − 21.38 𝑡
𝑥(0 𝑠) = 10 𝑚𝑚 = 0.01 𝑚
(𝑎 + 𝑏(0 𝑠))𝑒 − 21.38 (0 𝑠) = 0.01
𝑎 = 0.01
𝑚𝑚
𝑥̇ (0 𝑠) = 350
= 0.35 𝑚/𝑠
𝑠
[𝑏 − (21.38)(0.01 + 𝑏(0))]𝑒 − 21.38 (0) = 0.35 𝑚/𝑠
𝑏 = 0.5638
𝐶𝑅 = 2 𝐶𝑐 = 2(218.07) = 436.15
𝑥(𝑡) = 𝑎 𝑒 𝑆1𝑡 + 𝑏 𝑒 𝑆2 𝑡
𝑥(𝑡) = 𝑎 𝑒 −5.73 𝑡 + 𝑏 𝑒 −79.79 𝑡
pág. 8
LABORATORIO 2 VIBRACIÓN AMORTIGUADA
Derivando la ecuación.
𝑎 + 𝑏 = 0.01 𝑚 (𝑒𝑐𝑢𝑎𝑐𝑖𝑜𝑛 1)
𝑎 = 0.0155 ∧ 𝑏 = −0.0055
c) Subamortiguado ζ < 1.
𝐶𝑆𝑅 < 𝐶𝑐
𝑥̇ (𝑡) = −𝜁𝜔𝑛 𝑒 −𝜁𝜔𝑛 𝑡 [𝑎 cos 𝜔𝑑 𝑡 + 𝑏 sen 𝜔𝑑 𝑡] + 𝑒 −𝜁𝜔𝑛𝑡 [−𝑎𝜔𝑑 sen 𝜔𝑑 𝑡 + 𝑏𝜔𝑑 cos 𝜔𝑑 𝑡]
pág. 9
LABORATORIO 2 VIBRACIÓN AMORTIGUADA
𝑥(𝑡) = 𝑒 − 𝜁𝜔𝑛𝑡 (𝐴 𝐶𝑜𝑠 𝜔𝑑 𝑡 + 𝐵 𝑆𝑒𝑛 𝜔𝑑 𝑡)
𝑥(𝑡) = 𝑒 − (0.5)(21.38)𝑡 (𝐴 𝐶𝑜𝑠 18.51𝑡 + 𝐵 𝑆𝑒𝑛 18.51𝑡)
𝑥(𝑡) = 𝑒 −10.69 𝑡 (𝐴 𝐶𝑜𝑠 18.51𝑡 + 𝐵 𝑆𝑒𝑛 18.51𝑡) (Ecuación 1)
𝑥̇ (𝑡) = −𝜁𝜔𝑛 𝑒 −𝜁𝜔𝑛 𝑡 [𝐴 cos 𝜔𝑑 𝑡 + 𝐵 sen 𝜔𝑑 𝑡] + 𝑒 −𝜁𝜔𝑛𝑡 [−𝐴𝜔𝑑 sen 𝜔𝑑 𝑡 + 𝐵𝜔𝑑 cos 𝜔𝑑 𝑡]
𝑥̇ (𝑡) = −(0.5)(21.38)𝑒 −(0.5)(21.38)𝑡 [𝐴 cos 18.51𝑡 + 𝐵 sen 18.51𝑡] + 𝑒 −(0.5)(21.38)𝑡 [−𝐴(18.51) sen 18.51𝑡 + 𝐵(18.51) cos 18.51𝑡]
𝑥̇ (𝑡) = −10.69𝑒 −10.69𝑡 [𝐴 cos 18.51𝑡 + 𝐵 sen 18.51𝑡] + 𝑒 −10.69𝑡 [−18.51𝐴 sen 18.51𝑡 + 18.51𝐵 cos 18.51𝑡] (Ecuación 2)
𝑥(0 𝑠) = 0.1 𝑚
−10.69 (0) (𝐴
𝑒 𝐶𝑜𝑠 [18.51(0 𝑠)] + 𝐵 𝑆𝑒𝑛 [18.51(0 𝑠)]) = 0.1
𝑒 −0 (𝐴) = 0.1
𝐴 = 0.1
𝑥̇ (0 𝑠) = 0.35 𝑚/𝑠
−10.69𝑒 −10.69 (0) [𝐴 cos 18.51(0) + 𝐵 sen 18.51(0)] + 𝑒 −10.69(0) [−18.51𝐴 sen 18.51(0) + 18.51𝐵 cos 18.51(0)] = 0.35
pág. 10
LABORATORIO 2 VIBRACIÓN AMORTIGUADA
𝑠 2 𝑋(𝑠) + 457.45 𝑋(𝑠) + 0.1960 𝐶𝑠𝑋(𝑠) = 0
Como se puede observar solo existe una señal en este sistema. Lo que significa
que dicha señal es de entrada y salida, que existe un retorno hacia el inicio.
Despejando a nuestra conveniencia.
1
𝑋(𝑠) = ( 2 ) [−457.45 𝑋(𝑠) − 0.1960 𝐶𝑠𝑋(𝑠) ]
𝑠
𝟎. 𝟏𝟗𝟔𝟎𝑪 𝒔
-
X(s)
𝟏
𝒔𝟐
-
𝟒𝟓𝟕. 𝟒𝟓
pág. 11
LABORATORIO 2 VIBRACIÓN AMORTIGUADA
Diagrama de bloques original Diagrama de bloque equivalente
𝟎. 𝟏𝟗𝟔𝟎𝑪 𝒔
-
𝟏 𝟏 X(s)
𝒔 𝒔
-
𝟒𝟓𝟕. 𝟒𝟓
∫ 𝑉(𝑡) = 𝑥(𝑡)
Aplicando la transformada de Laplace.
𝑉(𝑠)
= 𝑋(𝑠)
𝑠
Eso quiere decir que la señal de salida del primer bloque integrador en serie es
la velocidad:
pág. 12
LABORATORIO 2 VIBRACIÓN AMORTIGUADA
𝟎. 𝟏𝟗𝟔𝟎𝑪
-
𝟏 𝟏 X(s)
𝒔 𝒔
-
V(s)
𝟒𝟓𝟕. 𝟒𝟓
pág. 13
LABORATORIO 2 VIBRACIÓN AMORTIGUADA
Como se puede notar se ingresan los
valores de la posición inicial y
velocidad inicial
pág. 14
LABORATORIO 2 VIBRACIÓN AMORTIGUADA
1. Se da dos clics izquierdos sobre el bloque para abrir el cuadro de
parámetros.
2. Se introduce el valor inicial que queremos en el cuadro “Initial condition”.
pág. 15
LABORATORIO 2 VIBRACIÓN AMORTIGUADA
2) Sobre amortiguado ζ > 1, C = 436.15.
pág. 16
LABORATORIO 2 VIBRACIÓN AMORTIGUADA
3) Subamortiguado ζ < 1, C = 109.04.
pág. 17
LABORATORIO 2 VIBRACIÓN AMORTIGUADA
❖ TALLER DE SIMULACIÓN 2 VIBRACIÓN AMORTIGUADA
pág. 18
LABORATORIO 2 VIBRACIÓN AMORTIGUADA