CLASE DE ECUACIONES DIFERENCIALES N° 17

APLICACIONES DE LAS ECUACIONES LINEALES DE SEGUNDO

ORDEN.
Un oscilador masa-resorte amortiguado está formado por una masa m
unida a un resorte fijo en un extremo, como se muestra en la figura.
Diseñe una ecuación diferencial que gobierne el movimiento de este
oscilador, tomando en cuente las fuerzas que actúan sobre él debido a la
elasticidad del resorte, la fricción (amortiguamiento) y las posibles
influencias externas. Observe la figura (1)

Figura 1
La segunda ley de Newton: fuerza es igual a masa por aceleración (F =
ma); es sin duda la ecuación diferencial que aparece con más frecuencia
en la práctica. Es una ecuación diferencial de segundo orden, pues la
aceleración es la segunda derivada de la posición con respecto al tiempo
d2 y
( ''
y =
dt 2 ) .
Al aplicar la segunda ley de Newton a un oscilador masa-resorte, los
movimientos resultantes son experiencias cotidianas y podemos
aprovechar nuestra familiaridad con estas vibraciones para obtener una
descripción cualitativa de las soluciones de ecuaciones de segundo orden
más generales.
En la figura se muestra el oscilador masa-resorte, cuando el resorte no
está estirado y la masa inercial m está en reposo, el sistema está en
equilibrio; medimos la coordenada y de la masa mediante su
desplazamiento a partir de la posición de equilibrio. Al desplazar la masa
m con respecto al equilibrio, el resorte se estira o se comprime y ejerce
una fuerza que resiste el desplazamiento. Para la mayor parte de los
resortes, esta fuerza es directamente proporcional al desplazamiento y,
por lo está dada por la ley:
Fresorte =−ky (1)
donde la constante k es la rigidez y el signo negativo la naturaleza de la
oposición de la fuerza. Esta ley se conoce como la Ley HOOKE. Esta ley
es válida solo para desplazamientos suficientemente pequeños.
Todos los sistemas mecánicos experimentan la fuerza de fricción, por lo
general, para el movimiento de vibración, ésta fuerza se modela mediante
un término proporcional a la velocidad:

F fricción =−b
dy
dt
=−b y
'

(2)

donde b≥0 es el coeficiente de amortiguamiento y el signo negativo

indica la naturaleza de oposición de la fuerza.
Las otras fuerzas que actúan sobre el oscilador se consideran por lo
general como externas al sistema.
En consecuencia, la ley de Newton permite escribir la ecuación
diferencial para el oscilador masa-resorte:

m y ' ' =−ky −b y ' + F externa ( t ) o

'' '
m y + b y + ky= F externa (t )

(3)
¿Qué apariencia tienen los movimientos del sistema masa-resorte? La
experiencia cotidiana con amortiguadores gastados, cuando no hay
fricción (b = 0) o fuerzas externas, los movimientos (idealizados) son
vibraciones perpetuas. Estas vibraciones se parecen a las funciones
sinusoidales, donde su amplitud depende del desplazamiento y las
velocidades iniciales. La frecuencia de las oscilaciones aumenta para
resortes más rígidos y desaparece para masas mayores.
Ejemplo. Verificar que el sinusoide con amortiguamiento exponencial
dado por y (t )=e−3 t cos 4 t es una solución de la ecuación (3) si
F externa ( t )= 0 , m = 1, k = 25 y b = 6.
Solución.
La ecuación (3) tomara la forma:

y ' ' +6 y ' +25 y=0

Ahora, sustituyendo y (t )=e−3 t cos 4 t en la ecuación, tenemos:

( e−3t cos 4 t ) ″+6 ( e−3t cos4 t )′ +25 ( e−3 t cos4t )=

−7e−3t cos4 t+24 e−2 t sen4 t+6 ( 7 .3−3 t cos 4 t−4 e−2t sen 4 t ) +25 e−3t cos4 t=0
Si un sistema masa-resorte es controlado por una fuerza externa
sinusoidal con frecuencia ω , nuestra experiencia indica que aunque la
respuesta inicial del sistema un tanto errática, con el tiempo se
“sincronizará” con el controlador y oscilará con la misma frecuencia.

Figura (2)
Los ejemplos comunes de sistemas que vibran sincronizados con sus
controladores son las bocinas de sistemas de sonido, ciclistas que
utilizan su bicicleta sobre vías de ferrocarril, circuitos electrónicos
amplificadores y las mareas del océano (controladas por la influencia
gravitacional de la luna). Sin embargo, la historia tiene más de lo que se
ha revelado hasta ahora. Los sistemas pueden ser muy sensibles a la
frecuencia particular ω con la que se les controla. Así, las notas
musicales afinadas con precisión pueden hacer vibrar un cristal fino, las
vibraciones inducidas por el viento con frecuencia correcta (o incorrecta,
según se vea) pueden derribar un puente y una llave de agua goteando
puede causar dolores de cabeza fuera de lo común. Estas respuestas
“resonantes” (para las que las soluciones tienen amplitudes máximas)
pueden ser algo destructivas, y los ingenieros en estructuras deben tener
mucho cuidado para garantizar que sus productos no resonarán con
alguna de las vibraciones que podrían ocurrir en el ambiente de
operación. Por otro lado, los ingenieros de radio realmente quieren que
sus receptores resuenen en forma selectiva con el canal de transmisión
deseado.
La ecuación:
'' '
a y +b y +cy =f ( t )

Entonces, como motivación de los procedimientos matemáticos y la teoría

necesaria para resolver esta ecuación, con frecuencia se compara con el
paradigma masa-resorte:
[ inercia ] × y ' ' + [ amortiguamiento ] × y ' + [ regidez ] × y =Fexterna

SISTEMAS DE MASAS RESORTE-MASA.

Movimiento libre no amortiguado.

La ecuación diferencial que describe este movimiento es,
[ inercia ] × y ' ' + [ amortiguamiento ] × y ' + [ regidez ] × y =Fexterna
'' '
m y + b y +ky =F externa

En este caso, nos centramos en b=0 ; F externa =0 , de modo, que la

ecuación diferencial del movimiento libre no amortiguado es,

m y ' ' + ky=0 (1)

Esta ecuación diferencial es llamada movimiento no amortiguado; es
decir movimiento sin amortiguaciones.
Al dividir la ecuación diferencial (1) entre m , se obtiene,

'' k
y + y=0
m

y '' +ω 2 y=0 (2)

k
Donde
ω=
√ m . La ecuación auxiliar asociada a (2) es,
2 2
r + ω =0

tiene raíces complejas r=±ω i . Por tanto, una solución general de la

ecuación es,
y (t )=C 1 cos ωt +C 2 sen ωt (3)
Podemos expresar la solución en su forma alternativa:

y (t )= Asen ( ωt+φ ) , (4)

con A≥0 , haciendo C1 = Asen φ y C 2= A cos φ ; es decir,

Asen (ωt+φ )= A cosωtsen φ+ Asen ωt cosφ

=C1 cos ωt+C2 senωt

Al despejar A y φ en términos de C1 y C2 se tiene.

C1
A= √C 21 +C 22 y tan φ=
C2

el cuadrante donde se encuentra φ queda determinado por los signos

de C1 y C2 . Esto se debe a que el signo de sen φ tiene el mismo signo
 C1
que C1 ( senφ=
A ) y cosφ tiene el mismo signo que C2

C2
( cos φ=
A )
Por ejemplo si C1 >0 y C2 <0 , entonces el ángulo fase φ estará
en el segundo cuadrante.
De la ecuación (4) podemos deducir que el movimiento de una masa en
un sistema libre sin amortiguación no es más que una onda sinusoidal, o
lo que se llama un movimiento armónico simple. La constante A es la
amplitud del movimiento y φ se llama ángulo fase. El movimiento es
2π ω k
periódico con periodo ω , frecuencia natural
El periodo se mide en unidades de tiempo y la frecuencia natural en
2 π , donde
ω=
√ m .

periodos (o ciclos) por unidad de tiempo. La constante ω es llamada

frecuencia angular para la función seno en la ecuación (4) y se mide en
radianes por unidad de tiempo.
La gráfica del movimiento no amortiguado, donde se observa el periodo
la frecuencia y la amplitud del movimiento.

Figura 3
Observe que la amplitud y el ángulo fase dependen de las constantes
C1 y C2 , que a su vez quedan determinadas por la posición y velocidad
iniciales de la masa. Sin embargo, el periodo y la frecuencia dependen
sólo de k y m y no de sus condiciones iniciales.
Ejemplo 1. Se fija una masa de 1/8 kg a un resorte con rigidez k = 16
N/m, La masa se desplaza 1/2 m a la derecha del punto de equilibrio y
luego se le imparte una velocidad hacia afuera (hacia la derecha) de
√2 m/seg. Despreciando cualesquiera fuerzas de amortiguamiento o
externas que puedan estar presentes, determinar la ecuación de
movimiento de la masa, junto con su amplitud, periodo y frecuencia
natural. ¿Después de cuánto tiempo pasa la masa por el punto de
equilibrio?
Solución. Como se trata de un movimiento no amortiguado, la solución de
la ecuación diferencial es de la forma (3); de modo que la frecuencia
angular es,

k
ω=
√ m ,
Usando los datos del problema tenemos,

k 16
ω=
√ m
=
√
1
8
=8 √ 2

rad/seg.
Al sustituir este valor de ω en la ecuación (3) tenemos:

y(t )=C 1 cosωt +C 2 senωt =C1 cos(8 √ 2t )+C 2 sen(8 √ 2t ) (5)

Ahora usando las condiciones iniciales:
'
1
y (0 )= m y { y (0 )= √2 m/esg ¿
2 ,

para encontrar los valores de las constantes C1 y C 2 en la ecuación

(5);
1
C1 cos (8 √ 2(0))+C2 sen(8 √ 2(0 ))=
y (0 )= 2
1
C1 cos(0 )+C2 sen(0 )=
2
1
C1 =
2
y ' (0)=−8 √2 C 1 sen(0)+8 √ 2 C2 cos(0 )=√ 2
1
C2=
8
Por tanto, la ecuación del movimiento es,
1 1
y (t )= cos(8 √2 t )+ sen(8 √ 2t )
2 8
La forma alternativa de la solución es:

17
y (t )= √ sen ( 8 √2 t +φ )
8
Siendo,

1 1 √ 17
A= √C 21 +C 22 =
1
√ + =
4 64 8

C1 2
tan φ= = =4 ⇒ φ=Arc tan ( 4 )≈1. 326
C2 1
8

Como C1 y C 2ambas son positivas, entonces φ está en el primer

cuadrante. Así, tenemos que la amplitud es A=
√ 17 m, y el ángulo fase φ
8
2 π 2 π √2 π
P= = =
es aproximadamente 1.326 radianes. El periodo es ω 8√2 8
1 8
=
seg, y la frecuencia natural es P √2 π ciclos por segundo.
Finalmente, para determinar el momento en que la masa pasará por el
posición de equilibrio; y=0 , debemos resolver la ecuación
trigonométrica:

17
y (t )= √ sen ( 8 √2 t +φ ) =0
8 ,
en términos de t.

√17 sen ( 8 √2 t+φ )=0 ⇔ 8 √2 t+φ=nπ o t= nπ−φ ≈ nπ−1 .326

8 8 √2 8 √2
Donde n es entero. Al hacer n=1 determinamos el primer instante en que
la masa pasa por la posición de equilibrio:
nπ−φ π −1. 326
t= = ≈0 . 16 segundos
8 √2 8 √2 .

Ejemplo 2.
Una masa que pesa 2 libras alarga un resorte 6 pulgadas. En t = 0 se
libera la masa desde un punto que está 8 pulgadas debajo de la posición
4
de equilibrio con una velocidad ascendente de 3 pies/seg. Determine
la ecuación del movimiento.
Solución. Debido a que se está usando el sistema de unidades de
ingeniería, la medición dada en términos de pulgadas se debe convertir a
pies.
 1 2
6 pulgadas = pie; 8 pulgadas= pie
2 3 . Además, se debe convertir las
W
m=
unidades de peso dadas en libras a unidades de masa. De g se
2 1 1
obtiene que
=
32 16
m=
slug
. También de la ley de Hooke,
2=k
2 ()
significa que la constante de resorte es k =4 lib/pie. Por tanto,
reemplazando en la ecuación (3) se tiene,
1 ''
x =−4 x
16
'
¿
'' 2 4
x +64 x=0 ; x (0 )= y { x (0)=−
3 3
Resolviendo se obtiene:
x (t )=C1 cos8 t+C 2 sen 8 t

2 1
C1 = y C 2 =−
Aplicando las condiciones iniciales se obtiene que 3 6 .

Por consiguiente, la ecuación del movimiento es:

2 1
x(t )= cos8 t− sen8 t
3 6 .
En su forma alternativa, la ecuación será:

17
x(t )= √ sen ( 8 t +1. 816 )
6
Pues,

4 1 √ 17
A= √ C 21 +C 22 =
2
√ + =
9 36 6
pie

C1 3
tan φ= = =−4 ⇒φ= Arc tan (−4 )≈−1 . 326
C2 1
−
6

Observemos que Arc tan (−4 )≈−1 .326 se localiza en el cuarto cuadrante,
pero el ángulo fase φ de acuerdo a los valores de las constantes
C1 >0 y C 2 <0 , se localiza en el segundo cuadrante, por tanto, este
ángulo debe tomarse en el cuadrante que le corresponde,
φ=π +(−1. 326 )=1. 8116 radianes.