Aplicaciones de Las Ecuaciones Diferenciales
Aplicaciones de Las Ecuaciones Diferenciales
Aplicaciones de Las Ecuaciones Diferenciales
11
5. M ovim iento Arm ónico Simple
Supóngase que un cuerpo de masa m está sujeto al extremo de un resorte flexible de
Resorte libre
uerpoe n movimiento
d 2x ,, x
m—2~ = m g — k\s + x)
= mg — ks —kx
kx,,
El
a lasigno
signo negativo indica que la ffuer
del movimiento. uerza
za de restitución del resorte actú a en direcció
direcciónn opue sta
Podemos escribir la ecuación (5.1) en la forma
2
d x k
— ¿ + -x = 0,
dt m
o bien
g ^-0 5.2)
2
donde u = k/m.
La ecuación (5.2) es la e c u a c ió n d if e rree n c ia l d e l m o v im ie n t o a r m ó n ic o s im p le
o movimiento vibratorio no amortiguado.
Hay dos condiciones iniciales asociadas con (5.2), a saber
f
x(0) = xo, x (0) = vo,
r
2
+ J2 = 0,
T\ = O>Z, T2 = — COÍ.
En consecuencia, la solución general de la ecuación diferencial (5.2) es
x(t) = c\ eeos
os or
oréé + C señaré, (5.3)
4 = (Jfe
de donde
Además m = W/g = 4/32 = 1/8 slug.
a) Luego, de (5.2), la ecuación diferencial que describe el movimiento es
cP x 8
o bien
64x = 0 (5.4)
dt 2
sujeta a las condiciones iniciales
z o) = o, Ao = \-
2
b) La ecuación auxiliar de (5.4) es r + 64
64 = 0, cuyas raíces son r = ±8z. En consecuencia
la solución general de (5.4) viene dada por
f
La condición incial x(0) = 0 implica que c\ = 0, mientras que x (0 ) = 1/3 conduce a
8c 2 = 1/3. De modo que c\ = 0, c 2 = 1/24 y la solución requerida es
x{t) — — s en 8 í.
Z4
Figura 5.2: Solución del ejemplo 1
c) La posición, velocidad y aceleración del peso 2 segundos después, están dadas, respec-
tivamente por
x {2 ) = — sen 16 = -0.011996,
i (2 ) = - eo
eoss 16 = -0.319 22,
ó
x {2) = - ^ sen 16 = 0.76774,
S o lu c ió n . Como antes
x(t) = C\ eos 8í + C sen 8í,
pero ahora, las condiciones iniciales son
x 0) = \, x'(0) = i
La condición x(0) = 1/4 exige de inmediato que c\ = 1/4, en tanto que usando x (0) =
1/3 se obtiene c^ — 1/24. Así que, la solución es
A = v/c? + 4 (5.7)
sen</>=~-
sen</>=~- y co s( /> = -j . (5.
(5.8)
Jx /x
T o ma n d o e n c u e n ta la s e x p re s io n e s (5 . 8 ), ju n to c o n e l h e c h o d e q u e e l ra n g o d e la
función f(x) = ar ct an z es el interv alo ((—
— —, — ) , concluimos que el valor de < > puede
calcularse simplemente como sigue
ar ct an ^ si C2 > 0,
t = { 5.9)
arctan ^ + TT si c<i < 0.
E J E M P L O 3 . R eescriba la sol
solución
ución ((5.5
5.5)) ddel
el eje
ejemplo
mplo 2 en la for
forma
ma alte rnativ a (5.6
(5.6)) y
determine el primer valor de t para el cual la masa pasa por la posición de equilibrio en
dirección hacia abajo.
Solución. En el ejemplo 2
A) V24/ 24
Además c 2 = 1/24 > 0, por lo cual de (5.9), se sigue que
1/4
f) = arctan - — - = arctan 6 = 1.4056 rad.
Por consiguiente
Los instantes en los cuales el cuerpo pasa por la posición de equilibrio, están deter-
minados por la condición
x(t) = 0,
es decir
sen (8í + 1.4056) = 0.
Las soluciones de esta ecuación vienen dadas por
8t +1.405 6 = 7ITT,
con
co n n un número entero. Despejando t y recordando que representa una cantidad positiva
(el tiempo), obtenemos la sucesión de valores
b) ¿En qué instante pasa el cuerpo por la posición de equilibrio en dirección hacia
arriba por tercera vez?
Solución Primero debemos observar que es necesario convertir a ft las longitudes ex-
presadas en pulgadas, usando la equivalencia
1 ft = 12 pulgadas.
9 ¡i 16
1
x(0) = -, x (0) = - 3 . 5.11)
i
x(t) = - c o s 4 > / 3 í - (5.12)
b) Escribimos la solución 5.12) en la for
forma
ma altern ativa. Tenemos que
1 3 __ 1
16 + 16 ~ 2
y como
c omo C2 < 0
/ 1 \ 7T 5
ó = arctan = +TT= h TI* = -TT.
V v3/ 6 6
Luego
= -sen
2
La gráfica de la ecuación del movimiento se muestra en la figura 5.4.
Los instantes en los cuales el cuerpo pasa por la posición de equilibrio vienen dados
p o r las soluciones de la ecuación x t) 0, es decir
1 / 5 \
-sen 4\/3í+-;7r) - 0 .
2 V 6 /
De aquí obtenemos la sucesión de valores de t
717T — | T
— 1, 2, o , . . .
í 5 = 1.8894 segundos.
c) Ahora, debemos hallar los valores de t para los cuales x t) — 1/4, esto es
se
senn
(i) _ 2 _^ + 2mr
P ^ = ^ « = 0 ,1
,1 ,2
,2 ,3
,3 ,
Obsérvese que en los tiempos tty tt y el cuerpo se mueve hacia aba jo d e la posición de
equilibrio, mientras que en los tiempos fflff l lo hace hacia arriba.
EJER CICIO S 5 1
2. U na fuer
fuerza
za de 10 N estira un resorte 0.12
0.1255 m. Después, al extrem o libre de ese
resorte se fija una masa de 5 kg.
b) Escriba la ecua
ecuación
ción de
dell movimiento en su forma
forma alternativ a.
c) ¿Cuántas oscilaciones completas realiza el cuerpo durante un intervalo de 8TT
segundos?
3. Cuando se sujeta una masa de 100 kg al extremo de un gran resorte, éste se estira
0.98
0.98 m. Se quita esta m asa y se reemplaza por un a de 40 kg, la cual se suelta desde
un punto que está 0.6 m debajo de la posición de equilibrio, con una velocidad
dirigida hacia arriba de 4 m/s.
a) Determine la ecuación del movimiento.
8. Int erp rete com o un movim iento armónico simple los siguientes proble ma s de valores
iniciales.
ld 2 x
a ) - ^ 2 + x = 0; z(0) - 2, z'(0) = - 4
1 d? x
b ) —^ + 25x = 0; s ( 0 ) = - 0 . 1 , x 0) = 3
+ 2Ar + u2 = 0,
2
r
xt
x {t)
{t ) - e~ (5.15)
x t) = e
Aí
(ci eo
eoss VUJ - X 2 t + c 2 se
senn V eo 2 - A £),
2
con C\ T
T¿¿ 0 y c2 ^ 0, puede escribirse como
donde
A= Jc\ 4 (5.19)
, c2
sen 0 = — y eos
eo s p — — .
Es decir
arctan si c 2 > 0,
5.20)
arctan ff-- + si c 2 < 0.
de donde
2
X t ¡) = 727T,
TÍ T —
n = 0,1,2,
(v^2- A2í + 0) = ±1 .
_ (2fc + 1)TT/
1) TT/2
2 - 4>
con k en N.
^ ¿
+ 2 0 ^ + 64x
64x = 0. (5.2
(5.21)
1)
dt dt
Las condiciones iniciales son
x(0) = ^, x'(0) = 0.
2
La ecuación auxiliar de (5.21) es r + 20r + 64 = 0 y sus raíces son r\ = — 4, r2 = —16,
—16,
de modo que
ci+c2 = i, (5.2
(5.22)
2)
4 1
Cl = 2=
9' ° 9
0.2 0.4 0.6 0.8 1 1.2
Por consiguiente
*<**> - 9 e
*< 9 e~ W - (5 24)
'
La gráfica de la solución (5.24) se da en la figura 5.8. Como se observa no ocurren
oscilaciones ya que el peso tiene tanto amortiguamiento que sólo retorna gradualmente a
la posi
posición
ción de equilibrio
equilibrio sin
sin pasar por esta. Se tr at a de un movimiento sobre am ortiguado .
La gráfica de la solución (5.26) se muestra con línea continua en la figura 5.9, donde
con línea interrumpida aparece también la solución del ejemplo 1, a fin de hacer una
1/ 3
0. 3
0 2
0 2 0 4 0 6 0
EJEMPLO 3. Si en el ejemplo 2 se
se cambia n las condiciones iniciales
iniciales a x(0) = 0 y
z'(0) = 1/3, determine x(t).
Solución. La solución general sigue siendo
8í t
x(t) = cie + 2cte ~
Sin embargo, las condiciones iniciales nos conducen a los valores de c\ = 0 y c2 = 1/3.
En consecuencia
x(t) = i í e 8 í . (5.27)
Se tiene que
A*)
A* ) = ^ - 8íe-8 t ) = ie- 8 t (l - 8í),
de dónde se observa que x(t) alcanza un desplazamiento máximo en t = 1/8 segundos
igual a max
mxa x = x(l/8) = 0.01533 ft. Ver figura 5.10.
EJEMPLO 4. Tomando en cuenta que (5= 1 rep ita el ejemplo 1. Esc riba la solución
en la forma alternativa (5.18). Determine los instantes en los que el cuerpo pasa por la
posición de equilibrio y realice la gráfica de la ecuación del movimiento.
Solución. La ecuación diferencial del movimiento es
4 2dx dx
3 2 dt
dt * * ~ di '
25
esto es
x x
t + 8 _ t+ 6 4 l = 0.
Ya as raíces
que llas carac terísticas son
sonr ^ = —4 ± 4\/3i, su solució
soluciónn general está da da porr
po
4í
x(t) = e (
1
3'
-4ci = 0.
x
(t ) = ie - 4 í (3cos4\/3í
(t) + V3 sen4V3t). (5.28)
1/3 7T
= arctan
V3/9 3
Luego, escribimos la
iH^-S*
solución (5.28)
2
en la forma
4í
Finalmente, la gráfica de x(t) es tan ge nte a las curvas exponenciales x\(t) —2>/3e /9
y 2(t) ~ —2\/3e~4 í /9 en los valores de t dados por la condición
se
senn
es decir
n<EN,
o bien
Gn +lO 11 \ X ,¿ J <> . . .
24^
La gráfica de (5.29) se muestra en la figura 5.11.
xx (t)
(t)
E J E M P L O 5. Después
Después de que un cuerpo que pesa 10 Ib se sujet sujetaa a un resorte de 5 ft
de largo, el resor te m ide 7 ft.
ft. Se qu ita el cuerpo d e 10 Ib y se
se le
le reempla za p or un o de 8
Ib . El sistema completo se coloca en un medio que ofrece una resistencia numéricamente
Ib.
igual a la velocidad instantánea.
2,
ó = arc tan 1 = —.
4
De modo que
-0.3
a) Sobreamortiguado.
b) Críticamente amortiguado.
c) Subamortiguado.
5. Una masa de 1/2 slug estira un resorte 4 ft y el medio que rodea al sistema masa-
resorte ofrece una resistencia al movimiento numéricamente igual a 4.5 veces la
velocidad ins tan tán ea. El peso se suelta 6 pulg adas aba jo de la posición de equilibrio
equilibrio
con una velocidad dirigida hacia arriba de VQ ft/s. ¿ Cómo debe ser vo para que la
masa pase por la posición de equilibrio ?
6. Dar una posible interpretación física del siguiente problema de valores iniciales
2
4 d dx 1
°
x(0)
x(0) = - 1 , 2/
2/((0) - - 1
7. Indique si las siguientes gráficas corresponden a un movimiento amortiguado, en
cuyo caso clasifíquelo.
b)
^ § (5-32)
o bien
^ + 2 A ^ + a ; 2 x = F (í
(í),
), (5.
(5.33)
33)
donde 2A = (3/m,uj2 = k /m y F(t) = f(t)/m.
Para resolver la ecuación no homogénea (5.33) podemos emplear el método de los
coeficientes indeterminados o el de variación de parámetros, según sea más conveniente.
EJEMPLO 1. Un resorte
de 1/2 slug. Se aplica vertical
una fuerza con constante
externa dada por de
f(t)6=lb/ft tiene
40 sen 2í, suspendida una emasa
í > 0. Supóngase
Supóngas que
actú a una fuerza
fuerza am ortigua dora n uméricam ente igual a dos
dos veces
veces la velocidad in stan táne a
y que inicialmen te el cuerp o está en reposo enen su posición de equilibrio. De term ine la
posición del cuerpo en cualquier tiempo t >0.
Solución. Con los valores de k =6 lb/ft, m = 1/2 slug y (3= 2, la ecuación diferencial
diferencial
de movimiento resultante es
^ + 4 ^ + l2 x = 80sen
80 sen 2í. (5.3
(5.34)
4)
¿
dt dt
La solución complementaria de (5.34) es
2 /
x c(t(t)) = e'\c x cos2\Í2t + c2 sen2v 2í).
Usando el método de los
los coe
coefi fici
cient
enteses ind etermina dos proponem os un a solución
solución p articula r
de (5.34) de la forma
(t ) = Aeos 2í + Bsen 2í.
x p(t)
En tal caso
x' p (t(t)) = -
x'¿(t) = -4Acos2t-4Bsen2t.
-4Acos2t-4Bsen2t.
Sustituyendo
Sustituyend o en 5.34), se sigue que
8 A + 8B) eos
eos 2í + 8B - 8A) sen 2í = 80 sen 2t.
8.4 + 8 5 = 0,
-8 ,4 + 8 5 = 80
80,,
x t) = 5e
5e~
~
x c t) = 5e~2í co
olución estacionaría
olución completa
5.3.1 Resonancia
Estudiaremos la ecuación (5.33) en el caso especial en que F(t) = Fo sen a i, i > 0, donde
F o y a son constantes positivas. La ecuación diferencial básica es
d x dx
-Tp + 2A
2A—— + u2x = F o sen a i , (5.35)
(5.35)
donde 2A = (3/m y u2 = k/m.
Supondremos que (3 es suficient
suficientemente
emente pequeño de mod o que el amortigu am iento es
menor que el crítico. En o tras pala bras, consideraremos que A < u. Luego, la solución
complementaria de (5.35) tiene la forma
(t)) = e xt(c 1 eos
x c (t /u2 - A2 í),
eo s yjuj2 - X H + c2 sen yy/u í) ,
con
co n C \ y C2
C2 con stan tes arb itra ria s, qu e depe nde n d e las condiciones
condiciones iniciales, o equivalen-
temente
(t)) = A e~xt se
x c (t senn {y/u2 - X H + 0 ),
donde A — ye2 + c2 , sen<f> = c\/A y eos
eos4> =
Ahora determinaremos
coeficientes indeterminados.unaSea
solución particular de (5.35), utilizando el método de los
Entonces
= —aB se
x' p (t
(t)) senn at + aC eos
eos a i,
t) = - a 2 B eos at - a2 C sen ai.
Zp(t)
Zp(
(- a 2 B
(-a + 2XaC + u 2B ) eos [{ u2 - a 2 ) C - 2aAJ5]
eos a i + [{u 2aAJ5] sen at = Fo sen
sen ai .