Ecuaciones Diferenciales en Movimiento Armónico
Ecuaciones Diferenciales en Movimiento Armónico
LINEALES
as
extremo superior fijado a una superficie horizontal y el otro extremo libre al
atic
cual se le fija una carga de masa m, la fuerza de recuperación del resorte esta
dada por la Ley de Hooke:
tem
F = ks
Ma
donde
s: elongación del resorte.
de
tuto
k: constante elástica del resorte.
nsti
F
ive
s s
posición de F
m
Un
•
equilibrio (P.E.) 0
x
" #mg m
La x bajo la posición de equilibrio
se considera positiva; x = 0 es P.E. mg
x+
Figura 4.2
4.11. APLICAC. DE LA E.D. DE SEGUNDO ORDEN 141
d2 x
m = −F + mg = −k(x + s) + mg
dt2
⇒ −F + mg = −kx −ks + mg
| {z }
= 0: en reposo
as
d2 x d2 x
atic
m 2 = −kx ⇒ m 2 + kx = 0
dt dt
tem
d2 x k
⇒ 2 + x=0
Ma
dt m
llamemos
k
= ω2 ⇒
d2 x
+ ω2x = 0 de
tuto
m dt 2
| {z }
E.D. segundo orden con coeficientes ctes.
nsti
Ecuación caracterı́stica
a, I
√
m2 + ω 2 = 0 ⇒ m = −ω 2 = ±ωi
qui
ntio
Sol. general
x(t) = C1 cos ωt + C2 sen ωt
eA
Condiciones iniciales:
x(0) = α
dd
x′ (0) = β
Si β > 0 el cuerpo inicia su movimiento con velocidad hacia abajo.
Llamemos q
C12 + C22 = A
y si hacemos
C1 C2
sen φ = , cos φ =
as
A A
atic
o sea que
C1
tan φ =
tem
C2
La constante A se le llama la amplitud del Movimiento Armónico Simple
Ma
(M.A.S.). y el ángulo φ se le llama ángulo de Fase.
de
Como tuto
A(C1 cos ωt + C2 sen ωt) C1 C2
x(t) = =A cos ωt + sen ωt
A A A
nsti
2π
Periodo de Vibraciones Libres T se define ası́: T = ω
qui
es decir: f = T1 = 2π
ω
eA
tenemos:
d2 x dx
Un
m 2 = −k(x + s) + mg − β
dt dt
donde β constante de fricción o constante de amortiguamiento.
d2 x dx
m 2
+β + kx = 0
dt dt
Dividiendo por m:
d2 x dx
2
+ 2λ + ω 2 x = 0
dt dt
4.11. APLICAC. DE LA E.D. DE SEGUNDO ORDEN 143
con carga y
con movimiento
as
P.E.: x = 0 •
atic
x F
tem
m
lı́quido
Ma
x+ mg
Figura 4.3
de
tuto
nsti
β k
Donde, 2λ = m
> 0 y ω2 = m
a, I
Ecuación caracterı́stica:
qui
p2 + 2λp + ω 2 = 0
√ √
ntio
√ √
λ2 −ω 2 )t λ2 −ω 2 )t
x(t) = C1 ep1 t + C2 ep2 t = C1 e(−λ+ + C2 e(−λ−
ive
por lo tanto
lı́m x(t) = 0
Un
t→∞
as
atic
0 t
tem
Figura 4.4 Sobreamortiguado
Ma
x(t)
de
tuto
1
nsti
a, I
qui
ntio
0 t
eA
add
√ √
C1 e−λt cos ω 2 − λ2 t + C2 e−λt sen ω 2 − λ2 t
x(t) = A
A
4.11. APLICAC. DE LA E.D. DE SEGUNDO ORDEN 145
x(t)
Ae−λt
as
atic
tem
−Ae−λt
Ma
Figura 4.6 Subamortiguado
de
tuto
−λtC1 √ C2 √
= Ae cos ω 2 − λ2 t + sen ω 2 − λ2 t
nsti
A A
−λt
√
2 2
= Ae sen ( ω − λ t + φ)
a, I
sumo una vez la posición de equilibrio como se ve en las gráficas 4.4 y 4.5
dd
d2 x dx
m = −k(x + s) + mg − β + F (t)
dt2 dt
d2 x dx
m 2
= −kx − ks + mg − β + F (t)
dt dt
d2 x dx
m 2 +β + kx = F (t)
| dt dt{z }
E.D. segundo orden de coef. ctes, no homogénea
146 CAPÍTULO 4. TEORIA DE LAS E.D.O. LINEALES
as
s
atic
P.E.: x = 0 •
tem
x F
m
Ma
lı́quido
mg F (t) (fuerza externa)
de
x+ tuto
Figura 4.7
nsti
a, I
d2 x
+ ω 2 x = F0 cos γt,
eA
dt2
donde F0 = constante, x(0) = 0 y x′ (0) = 0.
dd
m2 + ω 2 = 0 ⇒ m = ±ωi
rsid
1 1
xp (t) = 2 2
F0 cos γt = F0 cos γt
D +ω −γ + ω 2
2
F0
= cos γt
ω − γ2
2
Solución general:
F0
x(t) = C1 cos ωt + C2 sen ωt + cos γt
ω2 − γ 2
4.11. APLICAC. DE LA E.D. DE SEGUNDO ORDEN 147
x(0) = 0, x′ (0) = 0
entonces
F0
C1 = − , C2 = 0
ω2 − γ2
as
luego
atic
F0 F0
x(t) = (− cos ωt + cos γt) = 2 (cos γt − cos ωt)
tem
−γ ω2
2 ω − γ2
2F0 ω−γ ω+γ
Ma
= sen t sen t
ω2 − γ 2 2 2
de
ω−γ
Periodo de sen se calcula ası́:
tuto
2
ω−γ 4π
T1 = 2π ⇒ T1 =
nsti
2 ω−γ
a, I
Periodo de sen ω+γ
2
se calcula ası́:
qui
ω+γ 4π
T2 = 2π ⇒ T2 =
ntio
2 ω+γ
x(t)
dd
T2 2F0
ω 2 −γ 2
sin( ω−γ
2
)t
a
rsid
t
ive
Un
T1 = 4π
ω−γ − ω22F−γ0 2 sin( ω−γ
2
)t
as
atic
xh (t) = C1 cos ωt + C2 sen ωt
tem
1 1
xp (t) = 2 2
F0 sen ωt = 2 F (I eiωt )
2 0 m
D + ω D + ω
Ma
1 iωt 1
= F0 (Im 2 2
e ) = F0 (Im eiωt 1 )
D +ω (D + ωi)2 + ω 2
de
1
= F0 (Im eiωt 2 1 )
tuto
D + 2iωD − ω 2 + ω 2
1
nsti
iωt
= F0 (Im e 1 )
(D + 2iω)D
a, I
iωt 1 iωt 1 D
= F0 (Im e t ) = F0 (Im e 1− t )
D + 2iω 2iω 2iω
qui
1 1
= F0 (Im (cos ωt + i sen ωt) t− )
ntio
2iω 2iω
F0 1
= (−t cos ωt + sen ωt)
eA
2ω 2ω
F0
descartamos el término (2ω) 2 sen ωt ya que esta en xh ; la solución general es
dd
F0
x(t) = xh + xp = C1 cos ωt + C2 sen ωt − 2ω t cos ωt. Con las condiciones
a
rsid
F0 sen ωt F0 t→∞
|x(t)| = − t cos ωt −→ ∞
2ω 2 2ω
NOTA: la E.D. del movimiento vibratorio de los resortes también se
aplica en:
a). Circuitos en serie, en este caso, la E.D. es
d2 q dq 1
L 2
+ R + q = E(t)
dt dt C
4.11. APLICAC. DE LA E.D. DE SEGUNDO ORDEN 149
x(t)
F0
2ω t
as
t
atic
F0
− 2ω t
tem
Ma
de
tuto
Figura 4.9 Resonancia
nsti
a, I
L
ive
Un
E(t) R
Figura 4.10
150 CAPÍTULO 4. TEORIA DE LAS E.D.O. LINEALES
as
atic
tem
Eje elástico
Ma
de
tuto
θ(t)
nsti
a, I
qui
Figura 4.11
ntio
donde
eA
• c es la constante de amortiguación
a
rsid
θ a
α
as
m
atic
tem
s
Ma
Figura 4.12
de
tuto
d2 s 2
pero s = aθ, luego dt2
= a ddt2θ y sustituyendo en 4.12 se obtiene
nsti
d2 θ
a = −g sen θ (4.13)
dt2
a, I
d2 θ g
ntio
+ sen θ = 0 (4.14)
dt2 a
eA
d2 θ g
+ θ=0 (4.15)
a
dt2 a
rsid
as
atic
tem
Ma
b
Wv
170
T
θ de
tuto
b +
nsti
W
a, I
+
y
qui
dy dθ d2 y d2 θ
a
v= =r y a= = r
rsid
dt dt dt2 dt2
Por la segunda ley de Newton:
ive
d2 y
Un
X
de fuerzas en la dirección del eje y = m (4.18)
dt2
luego
W dy d2 y
−T +W − =m 2 (4.19)
170 dt dt
X d2 θ
: de torques alrededor del eje g = Ig 2
dt
4.11. APLICAC. DE LA E.D. DE SEGUNDO ORDEN 153
W r 2 1 d2 y W r d2 y
Tr = = (4.20)
g 2 r dt2 2g dt2
W d2 y W dy d2 y W d2 y
− + W − = m =
2g dt2 170 dt dt2 g dt2
as
3 W d2 y W dy
atic
2
+ =W
2 g dt 170 dt
tem
luego
d2 y g dy 2g
Ma
2
+ =
dt 255 dt 3
′
las condiciones iniciales son y(0) = 0 y y (0) = 0.
de
tuto
La solución general es
255
nsti
g
y = C1 + C2 e− 255 t + 170(t − )
g
a, I
y la solución particular es
qui
g g
la velocidad es
eA
g
y ′ = −170e− 255 t + 170
dd
la velocidad lı́mite
lı́m y ′ = 170 = vl
a
rsid
t→∞
g
y ′ (20) −170e− 255 20 + 170 g 4
Un
as
(Rta.: x(t) = − 14 cos 4 6t)
atic
Ejercicio 3. Un extremo de una banda elástica está sujeto a un punto
tem
A. Una masa de 1 kg. atada al otro extremo, estira la banda verticalmente
hasta un punto B, de tal manera que la longitud AB es 16 cm. mayor que la
Ma
longitud natural de la banda. Si la masa se estira más, hasta una posición
de 8 cm. debajo de B y se suelta; cual será su velocidad (despreciando la
de
resistencia),
√
al pasar por B?
g
tuto
(Rta.: ± 5 mts./seg.)
nsti
(Rta.: 0 < m ≤ 2
eA
hasta que su cara superior está a nivel de la superficie del agua y luego se
suelta. Se encuentra que el periodo es 1 seg.. Despreciando la resistencia, ha-
llar la E.D. del movimiento del bloque y hallar el peso P del bloque?(Ayuda:
la fuerza hacia arriba ejercida sobre el bloque es igual al peso del agua des-
alojada por el bloque; densidad de peso del agua: 62,4lb/pie2 )
2
(Rta.: ddt2x + 62,4g
P
x = 0, 62,4
4π 2
g
)
cuyo peso es 12,48 libras, flota en el agua (densidad del agua 62,4 libras/pie3 ),
manteniendo su eje vertical. Si se hunde de tal manera que la superficie del
agua quede tangente al bloque y luego se suelta. Cúal será el perı́odo de
vibración y la ecuación de su movimiento? Desprecie la resistencia.
Ayuda: La fuerza hacia arriba ejercida sobre el bloque es igual al peso del
agua desplazada por el bloque. √
(Rta.: √2π seg., x = ( 1
− 1) cos 5πg t pie.)
as
5πg 5π
atic
Ejercicio 8. Un peso de 10 libras sujeto a un resorte lo estira en 2 pies. El
tem
peso se sujeta a un mecanismo de amortiguación que ofrece una resistencia
numérica igual a β veces la velocidad instantánea (β > 0). Determinar los
Ma
valores de la constante de amortiguación β de modo que el movimiento sea:
a sobreamortiguado,
b crı́iticamente amortiguado,
c subamortiguado.
de
5 5 5
(Rta.:
a β > 2,
b β = 2,
c 0 < β < 2) tuto
Ejercicio 9. Después que un peso de 5 libras se sujeta a un resorte de
nsti
to que está a 21 pie bajo la posición de equilibrio con una velocidad dirigida
ntio
hacia arriba. √
(Rta.: a) x = 12 e−2t cos 4t − 21 e−2t sen 4t; b) x = 22 e−2t sen (4t + 3π
4
); c)
dd
nπ 3
t = 4 − 16 π, n = 1, 3, 5, . . .)
a
rsid
sistema está sobre una mesa que opone una fuerza de roce numéricamente
igual a 32 veces la velocidad instantánea. Inicialmente el peso está desplazado
Un
as
Ejercicio 12. Una masa de una libra sujeta a un resorte cuya constante
atic
es 9 libras/pie. El medio ofrece una resistencia al movimiento numéricamente
igual a 6 veces la velocidad instantánea. La masa se suelta desde un punto
tem
que está 8 pulg. sobre la posición de equilibrio con una velocidad dirigida
hacia abajo de v0 pies/seg.. Determinar los valores de v0 de modo que poste-
Ma
riormente la masa pase por la posición de equilibrio.
(Rta.: v0 > 2 pies/seg.)
de
tuto
Ejercicio 13. Un peso de 16 libras estira un resorte en 83 pies. Inicialmen-
te el peso parte del reposo desde un punto que está a 2 pies bajo la posición
nsti
una fuerza exterior f (t) = 8 sen 4t se aplica al sistema. Hallar la ecuación del
movimiento, si el medio que rodea al sistema opone una fuerza de amorti-
guación igual a 8 veces la velocidad instantánea.
(Rta.: x = 14 e−4t + te−4t − 14 cos 4t)
Ejercicio 17. Hallar las E.D. del sistema de resortes acoplados con cons-
tantes k1 y k2 y masas m1 y m2 respectivamente, como se muestra en la
as
figura 4.14
atic
2 2
(Rta.: m1 ddt2x = −k1 x + k2 (y − x), m2 ddt2y = −k2 (y − x))
tem
con carga y con carga y
Ma
en equilibrio en movimiento
k1 de
tuto
k1
nsti
P.E. m1 •
x
a, I
k2 m1
qui
ntio
P.E. m2 •
k2
eA
y
dd
m2
x+ y+
a
rsid
Figura 4.14
ive
Un
Ejercicio 18. Hallar las E.D. del sistema de tres resortes acoplados con
constantes k1 , k2 , k3 y masas m1 y m2 respectivamente, como se muestra
en la figura 4.15
2 2
(Rta.: m1 ddt2x = −k1 x + k2 (y − x), m2 ddt2y = −k2 (y − x) − k3 y)
k1
k1
m1 •
as
P.E.
x
atic
m1
k2
tem
k2
m2 x+•
P.E.
Ma
y
m2
de
k3
k3 y+
tuto
nsti
Figura 4.15
a, I
qui
√
mueve de acuerdo a y = sen 3gt. Encontrar la E.D. del movimiento del
ntio
peso y resolverla. √ √
2 √
(Rta.: g1 ddt2x = −4(x − y), x = −2 3 sen 2 g t + 4 sen 3g t)
eA
dd
t = 0 los resortes están sin estirar y el de peso 96 tiene una velocidad de 600
pies/seg. alejándose del muro y el otro esta en reposo. Encontrar las E.D. del
ive
P.E. + P.E. +
x y
m1 m2
K1 K2
96 64
as
Figura 4.16
atic
tem
Ma
l
k
de
tuto
nsti
b
C
x
a, I
+
qui
Figura 4.17
ntio
eA
constante k. Hallar
q el periodoqy la frecuencia de vibración del cilindro.
a
3m 1 8 k
rsid
(Rta.: T = 2π 8 k
, f= 2π 3m
, donde m es la masa del cilindro.)
ive
Ejemplo 42. Hallar con el paquete Maple las raı́ces del polinomio ca-
racterı́stico de la E.D. y (5) + 5y (4) − 2y (3) − 10y ′′ + y ′ + 5y = 0
Solución: la ecuacióm caracterı́stica de esta E.D. es
T1 b
B
θ T2
as
b +
atic
tem
+ W
y
Ma
Figura 4.18
>eq := m^5+5*m^4-2*m^3-10*m^2+m+5=0; de
tuto
eq := m5 + 5 ∗ m4 − 2 ∗ m3 − 10 ∗ m2 + m + 5 = 0
nsti
a, I
>solve(eq,m);
qui
ntio
>sols := [solve(eq,m)];
dd
a
rsid
> restart;
> DE(tools):diff(y(x),x,x,x)-3*diff(y(x),x,x)+9*diff(y(x),x)+13*y(x)=0;