Oscilaciones Forzadas
Oscilaciones Forzadas
Oscilaciones Forzadas
ma=-kx-λv+F0·cos(ωf t)
Expresamos la ecuación del movimiento en forma de ecuación diferencial que describe las
oscilaciones forzadas
dx F0 k λ
dt m m m
2 2
x1 = (C cos(ω t) + Dsin(ω t)) ⋅ exp(−γ t) ω = √ω − γ
0
x2=Acos(ωf t)+Bsin(ωf t)
Obtendremos los valores de A y B haciendo que cumpla la ecuación diferencial lineal completa
2 2
F (ω − ω ) 2γωf F
0 f
A = B =
2 2
2 2 2 2 2 2 2 2
m ((ω − ω ) + 4γ ω ) m ((ω − ω ) + 4γ ω )
0 f f 0 f f
El estado estacionario,
independiente de las condiciones
iniciales y es el que permanece
dx F ωf
v = = ω0 (−C sin(ω0 t) + D cos(ω0 t)) − sin(ωf t)
dt 2 2
m(ω −ω )
0 f
⎧ x0 = C +
F .
2 2
m(ω −ω )
t = 0 ⎨ 0 f
⎩
v0 = ω0 D
v0 F
x = x0 cos(ω0 t) + sin(ω0 t) + (cos(ωf t) − cos(ω0 t))
ω0 2 2
m(ω −ω )
0 f
Cociente F/m=1
Frecuencia angular propia, w0=100 rad/s y
Frecuencia de la fuerza oscilante, wf=120 rad/s
v0 F t sin(ω0 t)
x = x0 cos(ω0 t) + sin(ω0 t) +
ω0 m 2ω0
>> syms w0 wf t;
>> limit((cos(wf*t)-cos(w0*t))/(w0^2-wf^2),wf,w0)
ans =(t*sin(t*w0))/(2*w0)
O bien, resolvemos de nuevo la ecuación diferencial con ωf=ω0. Establecemos los siguientes
valores:
Rozamiento γ<ω0
Si las condiciones iniciales son t=0, x=x0, v=v0.
C = x0 − A
1
D = (v0 + γC − ωf B)
ω
Las condiciones iniciales más sencillas son x=0, y dx/dt=0 en el instante t=0. La partícula de
masa m parte del origen con velocidad inicial nula.
2 2 2 2
F (ω − ω ) F γ (ω + ω )
0 f 0 f
C = − D = −
2 2
2 2 2 2 2 2 2 2
m ((ω − ω ) + 4γ ω ) mω ((ω − ω ) + 4γ ω )
0 f f 0 f f
La posición de la partícula x que experimenta una oscilación forzada en función del tiempo t es
2 2
(ω − ω ) (cos(ωf t) − exp(−γ t) cos(ω t)) +
⎛ 0 f ⎞
F
x = ⎜ 2
ω +ω
2 ⎟
⎜ 0 f ⎟
2
2 2 2 2 2γωf (sin(ωf t) − exp(−γ t)sin(ω t))
m ((ω − ω ) + 4γ ω ) ⎝ 2ωωf ⎠
0 f f
Cuando ω0≠ωf
Establecemos los siguientes valores:
Cuando ω0=ωf
Establecemos los siguientes valores:
Actividades
Se introduce
La frecuencia propia del oscilador en ω0=100 rad/s, no se puede cambiar, pero se puede variar
la frecuencia de la fuerza oscilante ωf alrededor de dicha frecuencia.
x=-5, v=0
x=+5, v=0
x=0, v=+500
x=0, v=-500
Ejemplos de condiciones iniciales en las proximidades de la resonancia (ωf =90)
x=-5, v=0
x=+5, v=0
x=0, v=+500
x=0, v=-500
Sin rozamiento (γ =0)
Ejemplos de condiciones iniciales en la resonancia (ωf =100)
x=1.5, v=0
x=0, v=-50
x=0, v=150
Ejemplos de condiciones iniciales en las proximidades de la resonancia (ωf =90)
x=5, v=0
x=0, v=-10
Frec. fuerza: 100 Amplitud: 5
Posición: 0 Velocidad: 0
Amortiguamiento: 10 Nuevo ► ||
Solución numérica
Escribimos la ecuación diferencial de segundo orden que describe las oscilaciones forzadas en
forma de dos ecuaciones diferenciales de primer orden, para resolverlas utilizando la función
ode45 de MATLAB
2
d x dx 2 F0
+ 2γ + ω x = cos (ωf t)
dt2 dt 0 m
dx
⎧ = v
dt
⎨
⎩ dv 2 F0
= −2γv − ω x + cos (ωf t)
dt 0 m
En este ejemplo, ωf=100. Las condiciones iniciales son: x0=0, v0=-0.08. El estado transitorio
cambia, pero no lo hace la amplitud constante de la oscilación que se establece al cabo de un
cierto tiempo (teóricamente infinito)