Clase 8 - Ejercicio
Clase 8 - Ejercicio
Clase 8 - Ejercicio
Determinar el desplazamiento de un sistema con ξ ≔ 3% , k=11 (Ton/cm) y w=450 (Ton) para un empo de
1.2 s, si se aplica una carga armónica de 8 sin ((ωt)) (Ton) con una frecuencia igual a 0.5 veces y 1.5 veces la
frecuencia resonante. Determinar si el sistema vibra en fase o en desfase.
Se sabe además que la carga aplicada ene una duración de 20 s, después de los cuales se aplica otra carga
armónica con un 50% de la amplitud pero que hace que el sistema entre en resonancia. Determinar el
desplazamiento 5 ciclos después que se inicia la carga resonante.
Solución:
En primer lugar se deben determinar las propiedades del sistema a estudiar (en unidades compa bles ).
‾‾‾
k rad
ωn ≔ ―= 4.896 ―― Frecuencia natural
m s rad
ωD ≔ ωn ⋅ ‾‾‾‾‾
1 - ξ 2 = 4.894 ―― Frecuencia natural amor guada
s
2⋅π 2⋅π
Tn ≔ ―― = 1.283 s TD ≔ ―― = 1.284 s
ωn ωD
En ambos casos la carga depende de la "frecuencia resonante", la que se define como la frecuencia que hace
que el sistema entre en resonancia.
Se sabe que un sistema entrará en resonancia cuando cuando es excitado por una carga armónica con
frecuencia igual a la frecuencia natural, por lo tanto la frecuencia resonante es igual a ωn .
mm
Determinar A y B considerando condiciones iniciales iguales a 0 ( uinicial ≔ 0 mm , vinicial ≔ 0 ――) para t ≔ 0
s
solve , A , float , 3
A ≔ e -ξ ⋅ ωn ⋅ t ⋅ ⎛⎝A ⋅ cos ⎛⎝ωD ⋅ t⎞⎠ + B sin ⎛⎝ωD ⋅ t⎞⎠⎞⎠ + C ⋅ sin ((ω ⋅ t)) + D cos ((ω ⋅ t)) = uinicial ―――――→ 0.387 ⋅ mm
d ⎛ -ξ ⋅ ωn ⋅ t solve , B , float , 3
B ≔ ―― ⎝e ⋅ ⎛⎝A ⋅ cos ⎛⎝ωD ⋅ t⎞⎠ + B sin ⎛⎝ωD ⋅ t⎞⎠⎞⎠ + C ⋅ sin ((ω ⋅ t)) + D cos ((ω ⋅ t))⎞⎠ = vinicial ―――――→ -4.83 ⋅ mm
dt
u ((t)) ≔ e -ξ ⋅ ωn ⋅ t ⋅ ⎛⎝A ⋅ cos ⎛⎝ωD ⋅ t⎞⎠ + B sin ⎛⎝ωD ⋅ t⎞⎠⎞⎠ + C ⋅ sin ((ω ⋅ t)) + D cos ((ω ⋅ t))
Posición:
u ≔ e -ξ ⋅ ωn ⋅ t ⋅ ⎛⎝A ⋅ cos ⎛⎝ωD ⋅ t⎞⎠ + B sin ⎛⎝ωD ⋅ t⎞⎠⎞⎠ + C ⋅ sin ((ω ⋅ t)) + D cos ((ω ⋅ t)) = 4.253 mm
Velocidad:
d ⎛ -ξ ⋅ ωn ⋅ t mm
v ≔ ―― ⎝e ⋅ ⎛⎝A ⋅ cos ⎛⎝ωD ⋅ t⎞⎠ + B sin ⎛⎝ωD ⋅ t⎞⎠⎞⎠ + C ⋅ sin ((ω ⋅ t)) + D cos ((ω ⋅ t))⎞⎠ = -40.836 ――
dt s
Aceleración:
d 2 ⎛ -ξ ⋅ ωn ⋅ t mm
a ≔ ―― e ⋅ ⎛⎝A ⋅ cos ⎛⎝ωD ⋅ t⎞⎠ + B sin ⎛⎝ωD ⋅ t⎞⎠⎞⎠ + C ⋅ sin ((ω ⋅ t)) + D cos ((ω ⋅ t))⎞⎠ = -54.66 ――
2 ⎝
dt s2
10
7.5
4.659
5
u ((t)) ((mm)) 2.5
0
upermanente ((t)) ((mm)) -2.5
0 1 2 3 4 5 6 7 8 9 10
-5
utransitoria ((t)) ((mm)) -7.5
-10
-12.5
t ((s))
90
75
60
45
30
15
0
p ((t)) ((kN)) -15
0 1 2 3 4 5 6 7 8 9 10
-30
-45
-60
-75
-90
t ((s))
Como se observa en los gráficos, la respuesta total de desplazamiento (línea negra, corresponde a la suma de
u transitorio + u permanente) se mueve en el mismo sen do de la carga (gráfico con p(t)), lo que ra fica que
el movimiento está en fase.
Para carga 2, se ene:
rad
ω ≔ 2 ⋅ ωn = 9.792 ―― Frecuencia de la carga
s 2
⎛ω⎞
1 - ⎜―⎟
p0 ⎝ ωn ⎠
C ≔ ―⋅ ――――――――― 2
= -2.42 mm
k ⎛ 2⎞
2
⎛ω⎞ ⎛ ω⎞
⎜1 - ⎜―⎟ ⎟ + ⎜2 ⋅ ξ ⋅ ―⎟
⎝ ⎝ ωn ⎠ ⎠ ⎝ ωn ⎠
ω
-2 ⋅ ξ ⋅ ―
p0 ωn
D ≔ ―⋅ ――――――――― 2
= -0.097 mm
k ⎛ 2⎞
2
⎛ω⎞ ⎛ ω⎞
⎜1 - ⎜―⎟ ⎟ + ⎜2 ⋅ ξ ⋅ ―⎟
⎝ ⎝ ωn ⎠ ⎠ ⎝ ωn ⎠
mm
Determinar A y B considerando condiciones iniciales iguales a 0 ( uinicial ≔ 0 mm , vinicial ≔ 0 ――) para t ≔ 0
s
solve , A , float , 3
A ≔ e -ξ ⋅ ωn ⋅ t ⋅ ⎛⎝A ⋅ cos ⎛⎝ωD ⋅ t⎞⎠ + B sin ⎛⎝ωD ⋅ t⎞⎠⎞⎠ + C ⋅ sin ((ω ⋅ t)) + D cos ((ω ⋅ t)) = uinicial ―――――→ 0.0968 ⋅ mm
A = 0.097 mm
d ⎛ -ξ ⋅ ωn ⋅ t solve , B , float , 3
B ≔ ―― ⎝e ⋅ ⎛⎝A ⋅ cos ⎛⎝ωD ⋅ t⎞⎠ + B sin ⎛⎝ωD ⋅ t⎞⎠⎞⎠ + C ⋅ sin ((ω ⋅ t)) + D cos ((ω ⋅ t))⎞⎠ = vinicial ―――――→ 4.85 ⋅ mm
dt
B = 4.85 mm
Solución con condiciones iniciales nulas:
u ((t)) ≔ e -ξ ⋅ ωn ⋅ t ⋅ ⎛⎝A ⋅ cos ⎛⎝ωD ⋅ t⎞⎠ + B sin ⎛⎝ωD ⋅ t⎞⎠⎞⎠ + C ⋅ sin ((ω ⋅ t)) + D cos ((ω ⋅ t))
Evaluando para t ≔ 1.2 s :
Posición:
u ≔ e -ξ ⋅ ωn ⋅ t ⋅ ⎛⎝A ⋅ cos ⎛⎝ωD ⋅ t⎞⎠ + B sin ⎛⎝ωD ⋅ t⎞⎠⎞⎠ + C ⋅ sin ((ω ⋅ t)) + D cos ((ω ⋅ t)) = 0.148 mm
Velocidad:
d ⎛ -ξ ⋅ ωn ⋅ t mm
v ≔ ―― ⎝e ⋅ ⎛⎝A ⋅ cos ⎛⎝ωD ⋅ t⎞⎠ + B sin ⎛⎝ωD ⋅ t⎞⎠⎞⎠ + C ⋅ sin ((ω ⋅ t)) + D cos ((ω ⋅ t))⎞⎠ = 1.699 ――
dt s
Aceleración:
d 2 ⎛ -ξ ⋅ ωn ⋅ t mm
a ≔ ―― e ⋅ ⎛⎝A ⋅ cos ⎛⎝ωD ⋅ t⎞⎠ + B sin ⎛⎝ωD ⋅ t⎞⎠⎞⎠ + C ⋅ sin ((ω ⋅ t)) + D cos ((ω ⋅ t))⎞⎠ = -131.003 ――
2 ⎝
dt s2
t ((s))
90
75
60
45
30
15
0
p ((t)) ((kN)) -15
0 1 2 3 4 5 6 7 8 9 10
-30
-45
-60
-75
-90
t ((s))
Como se observa en los gráficos, la respuesta total de desplazamiento (línea negra, corresponde a la suma de
u transitorio + u permanente) se mueve en sen do opuesto a la carga (gráfico con p(t)), lo que ra fica que el
movimiento está en desfase.
Ahora, se sabe que la carga anterior se aplica por 20 segundos, y que después de ese instante comienza a
aplicarse una carga de la mitad de la amplitud pero con frecuencia resonante.
En este caso se debe separa el problema en 2 fases:
Fase 1 (de 0 a 20 segundos): vibración forzada con carga p ((t)) ≔ 8 tonne ⋅ g sin ((ω ⋅ t)) con condiciones
iniciales nulas.
Fase 2 (desde los 20 segundos): vibración forzada en resonancia con carga p ((t)) ≔ 4 tonne ⋅ g sin ⎛⎝ωn ⋅ t⎞⎠ con
condiciones iniciales iguales a las condiciones finales de la fase 1.
Se resolverá primero la carga 1.
Condiciones finales Fase 1.
rad
ω ≔ 0.5 ⋅ ωn = 2.448 ―― Frecuencia de la carga
s
p0 ≔ 8 tonne ⋅ g = 78.453 kN Amplitud de la carga
Evaluando para t ≔ 20 s :
Posición:
ufinal_Fase1 ≔ e -ξ ⋅ ωn ⋅ t ⋅ ⎛⎝A ⋅ cos ⎛⎝ωD ⋅ t⎞⎠ + B sin ⎛⎝ωD ⋅ t⎞⎠⎞⎠ + C ⋅ sin ((ω ⋅ t)) + D cos ((ω ⋅ t)) = -9.434 mm
Velocidad:
d ⎛ -ξ ⋅ ωn ⋅ t mm
vfinal_Fase1 ≔ ―― ⎝e ⋅ ⎛⎝A ⋅ cos ⎛⎝ωD ⋅ t⎞⎠ + B sin ⎛⎝ωD ⋅ t⎞⎠⎞⎠ + C ⋅ sin ((ω ⋅ t)) + D cos ((ω ⋅ t))⎞⎠ = 5.262 ――
dt s
Fase 2:
rad
ω ≔ ωn = 4.896 ―― Frecuencia de la carga
s
p0 ≔ 4 tonne ⋅ g = 39.227 kN Amplitud de la carga
uinicial ≔ ufinal_Fase1 = -9.434 mm
mm
vinicial ≔ vfinal_Fase1 = 5.262 ――
s
p0
ust_0 ≔ ―= 3.636 mm
k
ust_0 Solución para movimiento en
C ≔ 0 mm D ≔ -―― = -60.606 mm
2⋅ξ resonancia
A = 51.2 mm
d ⎛ -ξ ⋅ ωn ⋅ t solve , B , float , 3 1.54 ⋅ mm ⋅ rad + 1.08 ⋅ mm
B ≔ ―― e ⋅ ⎛⎝A ⋅ cos ⎛⎝ωD ⋅ t⎞⎠ + B sin ⎛⎝ωD ⋅ t⎞⎠⎞⎠ ↲⎞ = vinicial ―――――→ ―――――――――
dt ⎜+ C ⋅ sin (ω ⋅ t) + D cos (ω ⋅ t) ⎟ rad
⎝ ( ) ( ) ⎠
B = 2.62 mm
Para calcular el desplazamiento luego de 5 ciclos iniciada la fase 2, se debe calcular el empo que transcurre
durante 5 ciclos, el que es igual a t ≔ 5 ⋅ Tn = 6.417 s
NOTA: todas las fórmulas son derivadas par endo desde empo cero, por lo que, si bien la fase 2 inicia en el
segundo 20, para usar las ecuaciones se debe considerar un nuevo origen t=0 al comienzo de la fase 2.
Dado esto, se calcula para t = 6.417 s , aun cuando en realidad el empo total es t' ≔ t + 20 s = 26.417 s .
Posición luego de 5 ciclos:
u ≔ e -ξ ⋅ ωn ⋅ t ⋅ ⎛⎝A ⋅ cos ⎛⎝ωD ⋅ t⎞⎠ + B sin ⎛⎝ωD ⋅ t⎞⎠⎞⎠ + C ⋅ sin ((ω ⋅ t)) + D cos ((ω ⋅ t)) = -40.672 mm
d ⎛ -ξ ⋅ ωn ⋅ t mm
―― ⎝e ⋅ ⎛⎝A ⋅ cos ⎛⎝ωD ⋅ t⎞⎠ + B sin ⎛⎝ωD ⋅ t⎞⎠⎞⎠ + C ⋅ sin ((ω ⋅ t)) + D cos ((ω ⋅ t))⎞⎠ = 3.448 ――
dt s
d 2 ⎛ -ξ ⋅ ωn ⋅ t mm
―― e ⋅ ⎛⎝A ⋅ cos ⎛⎝ωD ⋅ t⎞⎠ + B sin ⎛⎝ωD ⋅ t⎞⎠⎞⎠ + C ⋅ sin ((ω ⋅ t)) + D cos ((ω ⋅ t))⎞⎠ = 973.965 ――
2 ⎝
dt s2
Evaluando para t ≔ 20 s :
Posición:
ufinal_Fase1 ≔ e -ξ ⋅ ωn ⋅ t ⋅ ⎛⎝A ⋅ cos ⎛⎝ωD ⋅ t⎞⎠ + B sin ⎛⎝ωD ⋅ t⎞⎠⎞⎠ + C ⋅ sin ((ω ⋅ t)) + D cos ((ω ⋅ t)) = -2.342 mm
Velocidad:
d ⎛ -ξ ⋅ ωn ⋅ t mm
vfinal_Fase1 ≔ ―― ⎝e ⋅ ⎛⎝A ⋅ cos ⎛⎝ωD ⋅ t⎞⎠ + B sin ⎛⎝ωD ⋅ t⎞⎠⎞⎠ + C ⋅ sin ((ω ⋅ t)) + D cos ((ω ⋅ t))⎞⎠ = -10.535 ――
dt s
Fase 2:
rad
ω ≔ ωn = 4.896 ―― Frecuencia de la carga
s
p0 ≔ 4 tonne ⋅ g = 39.227 kN Amplitud de la carga
uinicial ≔ ufinal_Fase1 = -2.342 mm
mm
vinicial ≔ vfinal_Fase1 = -10.535 ――
s
p0
ust_0 ≔ ―= 3.636 mm
k
ust_0 Solución para movimiento en
C ≔ 0 mm D ≔ -―― = -60.606 mm
2⋅ξ resonancia
A = 58.3 mm
d ⎛ -ξ ⋅ ωn ⋅ t solve , B , float , 3 1.75 ⋅ mm ⋅ rad - 2.15 ⋅ mm
B ≔ ―― e ⋅ ⎛⎝A ⋅ cos ⎛⎝ωD ⋅ t⎞⎠ + B sin ⎛⎝ωD ⋅ t⎞⎠⎞⎠ ↲⎞ = vinicial ―――――→ ―――――――――
dt ⎜+ C ⋅ sin (ω ⋅ t) + D cos (ω ⋅ t) ⎟ rad
⎝ ( ) ( ) ⎠
B = -0.4 mm
Para calcular el desplazamiento luego de 5 ciclos iniciada la fase 2, se debe calcular el empo que transcurre
durante 5 ciclos, el que es igual a t ≔ 5 ⋅ Tn = 6.417 s
Posición luego de 5 ciclos:
u ≔ e -ξ ⋅ ωn ⋅ t ⋅ ⎛⎝A ⋅ cos ⎛⎝ωD ⋅ t⎞⎠ + B sin ⎛⎝ωD ⋅ t⎞⎠⎞⎠ + C ⋅ sin ((ω ⋅ t)) + D cos ((ω ⋅ t)) = -37.889 mm
d ⎛ -ξ ⋅ ωn ⋅ t mm
―― ⎝e ⋅ ⎛⎝A ⋅ cos ⎛⎝ωD ⋅ t⎞⎠ + B sin ⎛⎝ωD ⋅ t⎞⎠⎞⎠ + C ⋅ sin ((ω ⋅ t)) + D cos ((ω ⋅ t))⎞⎠ = -2.527 ――
dt s
Aceleración luego de 5 ciclos:
d 2 ⎛ -ξ ⋅ ωn ⋅ t mm
―― ⎝ e ⋅ ⎛⎝A ⋅ cos ⎛⎝ωD ⋅ t⎞⎠ + B sin ⎛⎝ωD ⋅ t⎞⎠⎞⎠ + C ⋅ sin ((ω ⋅ t)) + D cos ((ω ⋅ t))⎞⎠ = 909.008 ――
d t2 s2