UNIVERSIDAD MAYA

ECUACIONES DIFERENCIALES

ECUACIONES DIFERENCIALES
Y
SU APLICACIÓN A LA
INGENIERÍA

ING. ELIZABETH MÉNDEZ

SÁNCHEZ

CARLOS GUZMAN MORENO

APLICACIONES DE ECUACIONES DIFERENCIALES LINEALES

Movimiento vibratorio de sistemas mecánicos

El sistema más simple disponible para estudiar el movimiento vibratorio

consiste en un resorte de peso despreciable [figura(a)] suspendido verticalmente
de un soporte fijo.
Suponga que un peso W se cuelga del resorte [figura (b)].

Cuando el peso está en reposo describimos su posición como la posición de

equilibrio. Si el peso se desplaza ejerciendo una fuerza vertical y hacia abajo y se
deja de ejercer cuando tiene una cierta distancia respeto de la posición de
equilibrio, estará bajo un movimiento vibratorio alrededor de la posición de
equilibrio [figura(c)].Queremos averiguar el movimiento que realiza el cuerpo en
su desplazamiento respecto de la posición de equilibrio. Para conseguir este
propósito, tendremos que conocer las fuerzas que actúan sobre el peso durante
su movimiento. Por la experiencia vemos que hay una fuerza que hace regresar
o restaurar un peso desplazado a su posición de equilibrio. Esta fuerza se llama
la fuerza restauradora.
La ley que gobierna esta fuerza es un caso especial de la ley generalizada de
Hooke. Nos referiremos a este caso especial como la ley de Hooke, la cual se
enuncia como sigue:
La fuerza ejercida por un resorte, tendiente a restaurar el peso W a la posición de
equilibrio, es proporcional a la distancia de W a la posición de equilibrio. (“la fuerza
es proporcional al alargamiento”).
Denotamos la magnitud de la fuerza restauradora por | f |, y sea x la posición de W
medida desde la posición de equilibrio. Se supone la dirección positiva hacia abajo,
de modo que x es positivo cuando W está por debajo de la posición de equilibrio y
negativo cuando W esté por encima de esta posición.

De acuerdo a la ley de Hooke,

fαx o f = k x donde k > 0

Es una constante de proporcionalidad, que depende de la dureza del resorte y se

llama constante del resorte o módulo de elasticidad del resorte. Para determinar
la dirección de la fuerza, observamos que cuando x > 0 la fuerza está dirigida
hacia arriba y por tanto es negativa. Cuando x < 0 la fuerza está dirigida hacia
abajo y es por tanto positiva. Esto se puede satisfacer, sólo si la fuerza está dada
tanto en magnitud como dirección por - kx, de modo que la ley de Hooke es:
f=-kx

Cuando el peso W se coloca en el resorte, se estira una distancia s como se ve en

la anterior figura en la posición (b), según la ley de Hooke, la tensión T en el resorte
es proporcional al elongamiento, y así T1 = ks; puesto que el resorte y el peso están
en equilibrio se tiene que
T1 = ks= W

Cuando el peso se baja más y se suelta, su posición en cualquier tiempo se muestra

en la figura (c). La tensión T, en el resorte en este tiempo es, de acuerdo a la ley
de Hooke, T2 = k(s + x)

Sigue que la fuerza neta en la dirección positiva está dada por:

W-T2=W-ks-kx=-kx

Por la ley de Newton, la ecuación diferencial que describe el movimiento es:

W d x2
2 = −kx
g dt

Así, la fuerza neta es simplemente la fuerza restauradora y no depende del peso

W.

3
El tiempo para un ciclo completo se llama periodo. Del gráfico vemos que el período
es

= T ; otra forma de ver que el período es ,sin el gráfico, es determinando cuándo

el peso está en un extremo de su trayectoria (esto es, ya sea el punto más alto, o
más bajo). Cuando tomamos el punto más bajo dado por x = A cos Et, esto ocurrirá
cuando cos Et = 1, esto es, Et = 0, 2 π, 4 π, 6π, . . o t = 0; π /4, 2 π /4, 3 π /4,
. . .la primera vez que x = max es cuando t = 0, la segunda cuando t = π /4, la
tercera cuando t = 2 π /4, etc. La diferencia entre tiempos sucesivos es π /4, la cual
es el período. El número de ciclos por segundos se llama frecuencia.

Período (T) = Número unidades de tiempo por ciclo = π /4

1 4
Frecuencia (f) = Número de ciclos por unidad de tiempo = =

1
En general si T es el período, la frecuencia f está dada por f =
T

Si ahora de la solución general x = A cos Et + B sen Et damos las condiciones

iniciales: dx

≠ 0 x ≠ 0 en t = 0 .
dt
Despejando las constantes A y B obtenemos una ecuación de la forma:

x = a cos Et + b sen Et llamando a E = ω tendremos x = a cos ω t + b sen ω t

esta ecuación la podemos poner de la forma:

a cos ωt + bsen ωt = a2 + b sen2 (ω φt + )

obtenida de a2 +b sen2(ω φt + =)a2 +b2 (sen tω φ ω φcos + cos tsen

)=
 2 2⎡ ⎛ b ⎞ ⎛ a ⎞⎤
= a +b ⎢(sen tω)⎜⎜ 2 2 ⎟⎟+(cosωt)⎜⎜ 2 2 ⎟⎟⎥ = acosω ωt +bsen t
⎢⎣ ⎝ a +b ⎠ ⎝ a +b ⎠⎥⎦

El ángulo ø a menudo se llama el ángulo de fase.

Si un movimiento se puede describir de la forma:

1ω 2π
x = Asen(ωt +φ) entonces f== T=
T 2π ω

Obteniendo la relación ω = 2πf

El resorte vibrante con amortiguamiento (movimiento

amortiguado).

Los resortes son modelos idealizados fuera de los casos reales, puesto que las
oscilaciones no disminuían, como uno esperaría por experiencia, sino por el
contrario, se mantenían constantes en el tiempo. En la práctica, fuerzas de fricción
y otras (tales como la resistencia del aire) actúan para decrecer la amplitud de las
oscilaciones y finalmente traer el sistema al reposo. Una manera para obtener una
mejor aproximación a la realidad, es suponiendo una fuerza amortiguadora. La ley
exacta para esta fuerza no se conoce, puesto que depende de muchas variables,
pero se ha encontrado mediante la experimentación, que para velocidades

5
pequeñas, la magnitud de la fuerza amortiguadora es aproximadamente
proporcional a la velocidad instantánea del peso en el resorte.

La magnitud por tanto está dada por:

dx
β dt donde β es la constante de proporcionalidad llamada la constante de
amortiguamiento.

La fuerza amortiguadora se opone al movimiento, de modo, que cuando el peso va

bajando la fuerza amortiguadora actúa hacia arriba, mientras que actúa hacia abajo
cuando el peso va subiendo. Asumiendo hacia abajo como la dirección positiva,
vemos que la fuerza amortiguadora debe ser negativa cuando dx/dt sea positiva, y
debe ser positiva cuando dx/dt sea negativa. Así, con β > 0, es claro que la fuerza
amortiguadora debe estar dada tanto en magnitud como en dirección por -β dx/dt
.

Cuando se tienen en cuenta las fuerzas restauradoras ya consideradas, según la

ley de Newton, la ecuación diferencial del movimiento es:

W d x2 dx W d x2 dx
2 = −β −kx o 2 +β + kx = 0 g
dt dt g dt dt

Ejemplo aclaratorio:

Un peso W estira un resorte x1 unidades de longitud. Si el peso se halla x2 unidades

de longitud por debajo de la posición de equilibrio y se suelta:

(a) Establecer la ecuación diferencial y condiciones asociadas que describan el

movimiento.

(b) Encontrar la posición del peso como una función del tiempo.
Formulación matemática:

Teniendo en cuenta la fuerza

Amortiguadora

La ecuación del movimiento será 2 +β + kx = 0

dx
dt dt

donde sacamos que (m2 +mF + E e) mt = 0 para la cual (m2 +mF + =E) 0 y tiene
raíces
− F ± F2 −4E
m = De donde la ecuación diferencial tiene la solución de la forma: 2
2 2
x = c e1 mt1 +c e2 −m t2 donde m1 = − F +F −4E m2 = − F − F −4E
2 2

De otra forma m =− ±a bi si F 2 −4E < 0 x = c e1 t(− +a bi) +c e2 t(− −a bi) = e−ta[c

e1 tbi +c e2 −tbi ] siendo la expresión x =[c e1 tbi +c e2 −tbi ] cuya solución general
debe ser como la ecuación x A sen bt B cos bt = + , son las dos ecuaciones
generales. Debe haber una correlación

A sen bt B cos bt + ≡ c e1 tbi +c e2 −tbi si consideramos que en t=0 e±tbi =1

tendremos

B = +c1 c2 diferenciando en ambos lados de la ecuación A sen bt B cos bt +

≡ c e1 +c e2 −tbi
tbi

Æ b(A cos bt B sen bt )− ≡ bi c e[ 1 bit −c e2 −bit ] para t=0 b(A ) ≡ bi c[ 1 −c2]

con lo que
A = −i c[ 1 c2] con los valores obtenidos antes de A y B , B = c1 +c2 , A = i c[ 1

−c2], tendremos , c1(cosbt+isen bt) +c2(cosbt−isen bt) ≡ c e1 bit −c e2 −bit con lo cual

ebit = cosbt +isen bt .

7
De la condición x = x2 en t = 0, encontramos A, así que sustituyendo en
dx
x A sen bt B cos bt = +
y derivando respecto de t e introduciendo la condición = 0 en dt
t=0 podemos encontrar el valor de B y con lo cual la x para las condiciones
dadas.

Por analogía de las ecuaciones:

A sen bt B cos bt +≡ c e1 tbi +c e2 −tbi

[
A sen bt B cos bt +≡ e−ta c e1 tbi +c e2 −tbi ]
[ ]
Tendremos que la segunda ecuación se convierte en x = e−ta c e1 tbi +c e2 −tbi donde

su ecuación solución será: x = e−ta[Acosbt + Bsenbt]

Determinando las constantes A y B sujetas a las condiciones citadas anteriormente,

encontramos la ecuación de la posición en función del tiempo, que determina una gráfica de
este tipo:

2 2
a cos ω t + bsen ω t = a + bsen t+ )
(ωφ
Si hacemos uso de la identidad
y de la ecuación x = e−ta[Acosbt + Bsenbt] se puede escribir:

x = e−ta c2 +d2[sen(ω φt + )]

El gráfico de la figura está entre los gráficos de líneas punteadas puesto que el
seno varía entre - 1 y + 1. El punto P representa un mínimo relativo, mientras que
el punto Q está sobre la curva punteada. La diferencia constante en tiempos entre
máximos (o mínimos) sucesivos se llama cuasi período, aunque se puede
referenciar éste como el período. El adjetivo cuasi se usa, puesto que los valores
funcionales no se repiten como lo harían si realmente hubiera periodicidad.

El movimiento descrito en este ejemplo se llama movimiento oscilatorio

amortiguado o movimiento amortiguado.

Se debería notar que x = e−ta c2 +d2[sen(ω φt + )] tiene la

forma x =α( )t sen [ (ωt +φ)]

El cuasi período está dado por: T =

El resorte con fuerzas externas

Tales fuerzas pueden ocurrir, por ejemplo, cuando el soporte que sostiene el
resorte se mueve arriba y abajo en una manera especificada tal como en
movimiento periódico, o cuando al peso se le da un pequeño empuje cada vez que
alcanza la posición más baja. Si denotamos la fuerza externa por F(t), la ecuación
diferencial para el movimiento es:

W d x2 dx W d x2 dx
2 = −β −kx + F t( ) o 2 +β + kx = F t( ) g dt

dt g dt dt 2

9
 dx dx dt dt
Donde la última ecuación puede escribirse: a 2 +b + =cx F t( )
(donde a = W/g, b= β, c = k), a menudo llamada la ecuación de vibraciones
forzadas.

Lineal

Un conjunto de funciones y1 (x), y2 (x),…, yn (x) denotadas como y1, y2,…, yn , se

dice que son linealmente dependientes en un intervalo si existe un conjunto de n
constantes, no todas cero, tales que en el intervalo α1 1y +α2y2 +...+αn yn ≡ 0 en
caso contrario se dice que el conjunto es linealmente independiente.

La dependencia o independencia se averigua mediante el Wronskiano.

Teorema: Si y1, y2 son linealmente dependientes en un intervalo, y si sus derivadas

y´1, y´2 existen en el intervalo, entonces el Wronskiano de y1, y2 dado por W y y( 1,

2) = dety y1, 2, ´ , ´y 1 y 2 = y y1 2´ −y y2 1´

es idénticamente cero (esto es, W= 0) en el intervalo.

Este teorema puede enunciarse en términos de independencia lineal como sigue:

Si el Wronskiano de y1, y2 no es idénticamente cero (esto es, W≠ 0) en un intervalo,

entonces y1, y2 son linealmente independientes en el intervalo.

si y1, y2 fueran linealmente dependientes en el intervalo su Wronskiano sería

idénticamente cero en el intervalo según el teorema. Esta contradicción muestra
que las funciones no son linealmente independientes, esto es, son linealmente
dependientes en el intervalo. Con lo cual, si W = y y1 ´2−y y2 ´1≠ 0 tenemos α1 = 0 α2
= 0.

Otros teoremas importantes nos dicen:

-Si x1,x2 son soluciones linealmente independientes de la ecuación (aD2 + +bD c x)
= 0, entonces x = c x1 1 +c x2 2 es una solución de la ecuación para cualesquiera
constantes c c1, 2 .

-Si x1,x2 son soluciones linealmente independientes de (aD2 + +bD c x) = 0 y xp es

una solución particular de (aD2 + +bD c x) = F t( ) entonces x = c x1 1 +c x2 2 + xp es

una solución de la ecuación diferencial (aD2 + +bD c x) = F t( ) para cualesquiera

constantes c1,c2 .

La resonancia mecánica

Cuando la frecuencia de una fuerza externa periódica aplicada a un sistema

mecánico está relacionada de una manera sencilla con la frecuencia natural del
sistema, puede ocurrir el fenómeno de la resonancia mecánica, la cual eleva las
oscilaciones tales que el sistema puede desplomarse o colapsarse. Una compañía
de soldados marchando en fila a través de un puente puede de esta manera hacer
que el puente colapse. En una manera análoga, puede ser posible que una nota
musical de una frecuencia característica propia estalle un cristal.
La resonancia mecánica puede también servir para propósitos útiles. Por ejemplo,
si un automóvil se quedara atascado en la nieve (o barro), puede, al “mecerlo”,
ser puesto a vibrar con su frecuencia natural. Entonces, al aplicarle una fuerza
con esta misma frecuencia, la resonancia mecánica resultante puede ser algunas
veces suficiente para liberarlo.
Un peso P estira un resorte x1 unidades de longitud. Si el peso se halla x2
unidades de longitud por debajo de la posición de equilibrio y se suelta aplicándole
una fuerza externa dada por F t( ) = AcosBt .
Describir el movimiento que resulta si se supone que inicialmente el peso está en
la posición de equilibrio (x = 0) y que su velocidad inicial es cero.

Puesto que inicialmente (t = 0) el peso está x2 por debajo de la posición de

equilibrio, tenemos x = x2 en t = 0 .
También, puesto que el peso se suelta (esto es, tiene velocidad cero) en t = 0,
dx
= 0 en t = 0 . dt
Solución:
Ejemplo aclaratorio:

11
Como F t( ) = AcosBt contiene una función de los términos cosBt, con lo que
podemos utilizar el método de los coeficientes indeterminados para determinar la
solución particular.
Ensayando como solución particular x p = c senBt1 + c2 cos
Bt de
P 2
( D +k x) = F t( ) g
xp = c senBt1 +c2 cosBt x´p= Bc1cosBt − Bc senBt2 x´´p=−B c senBt2 1
− B c2 2 cosBt
P
Sustituyendo en la ecuación x´´+ =Kx F t( ) calculamos el valor de las
constantes c c1 2, y g
por tanto la solución particular.
P d x2
La solución complementaria de 2 = −kx+ F t( ) se obtiene de la misma forma
que en el g dt
Apartado anterior 1.1 con lo que sería xc = AcosBt + BsenBt con lo que la solución
general será de la forma:
x = + =xc xp AcosBt + BsenBt +c senBt1 +c2 cosBt

Usando las condiciones iniciales, podemos observar una solución será cuando A
= B = 0 de donde podemos sacar que: x = c senBt1 +c2 cosBt

El gráfico de la ecuación x = c senBt1 +c2 cosBt está entre los gráficos de x = t y

x = - t como se muestra en la figura. Se puede ver en el gráfico que las oscilaciones
van creciendo sin límite. Naturalmente, el resorte está limitado a romperse dentro
de un corto tiempo.