Clase 8 - Ejercicio

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 7

Ejercicio Clase 8_ Vibración forzada amor guada con carga armónica

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 ).

m ≔ 450 tonne Masa


w ≔ m ⋅ g = 4412.993 kN Peso
tonne ⋅ g kN Rigidez
k ≔ 11 ―――= 10787.315 ――
cm m
ξ = 0.03

‾‾‾
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

Se consideran 2 condiciones de carga:


1) Frecuencia de la carga igual a 0.5 veces la frecuencia resonante.
2) Frecuencia de la carga igual a 1.5 veces la frecuencia resonante.

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 .

Resolviendo para la condición de carga 1 con condiciones iniciales iguales a 0:


rad
ω ≔ 0.5 ⋅ ωn = 2.448 ―― Frecuencia de la carga
s
p0 ≔ 8 tonne ⋅ g = 78.453 kN Amplitud de la carga
⎛ ω ⎞2
1 - ⎜―⎟
p0 ⎝ ωn ⎠
C ≔ ―⋅ ――――――――― 2
= 9.681 mm
k ⎛ 2⎞
2
⎛ω⎞ ⎛ ω⎞
⎜1 - ⎜―⎟ ⎟ + ⎜2 ⋅ ξ ⋅ ―⎟
⎝ ⎝ ωn ⎠ ⎠ ⎝ ωn ⎠
ω
-2 ⋅ ξ ⋅ ―
p0 ωn
D ≔ ―⋅ ――――――――― 2
= -0.387 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.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

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)) = 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

Para determinar si el movimiento está en fase o en desfase, se evalúa la razón de sintonía:


ω
― = 0.5
ωn Como la razón de sintonía es menor que 1 y ϕ es cercano a 0 , la
⎛ -D ⎞ respuesta de desplazamiento del sistema en estado permanente está
ϕ ≔ atan ⎜―― ⎟ = 0.04 rad
⎝ C ⎠ prác camente en fase total con la carga
1.2
12.5

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

Para determinar si el movimiento está en fase o en desfase, se evalúa la razón de sintonía:


ω
―=2 Como la razón de sintonía es mayor que 1 y ϕ es menor que 0, la
ωn
⎛ D⎞
respuesta de desplazamiento en estado permanente está desfasada
ϕ ≔ atan ⎜-―⎟ = -0.04 rad con la carga.
⎝ C⎠
1.2
7
6
5
4
3
u ((t)) ((mm)) 0.419
2
1
0
upermanente ((t)) ((mm)) -1 0 1 2 3 4 5 6 7 8 9 10
-2
-3
utransitoria ((t)) ((mm)) -4
-5
-6

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

A ≔ 0.897 mm B ≔ -4.753 mm Si condiciones iniciales son


C ≔ 9.613 mm D ≔ -0.897 mm nulas

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

Con esas condiciones iniciales, despejando 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 ―――――→ 51.2 ⋅ mm

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

Velocidad luego de 5 ciclos:

d ⎛ -ξ ⋅ ωn ⋅ t mm
―― ⎝e ⋅ ⎛⎝A ⋅ cos ⎛⎝ωD ⋅ t⎞⎠ + B sin ⎛⎝ωD ⋅ t⎞⎠⎞⎠ + C ⋅ sin ((ω ⋅ t)) + D cos ((ω ⋅ t))⎞⎠ = 3.448 ――
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))⎞⎠ = 973.965 ――
2 ⎝
dt s2

Situación con carga 2:


Condiciones finales Fase 1.
rad
ω ≔ 2 ⋅ ωn = 9.792 ―― Frecuencia de la carga
s
p0 ≔ 8 tonne ⋅ g = 78.453 kN Amplitud de la carga

A ≔ 0.224 mm B ≔ 4.831 mm Si condiciones iniciales son


nulas
C ≔ -2.403 mm D ≔ -0.224 mm

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

Con esas condiciones iniciales, despejando 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 ―――――→ 58.3 ⋅ mm

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

Velocidad luego de 5 ciclos:

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

También podría gustarte