Escuela de Ingeniería Civil Páá giná 1

INTRODUCCIÓN

Las Ecuaciones Diferenciales tienen una importancia fundamental en la Matemáticas

para la ingeniería debido a que muchos problemas se representan a través de leyes y
relaciones físicas matemáticamente por este tipo de ecuaciones. Es interés de este trabajo la
deducción de las Ecuaciones Diferenciales a partir de situaciones físicas que se presentan en
determinados problemas de carácter físico.

A esta transición del problema, al Modelo Matemático correspondiente se llama

Modelado. Este método tiene una gran importancia práctica para el ingeniero y se ilustra por
medio de ejemplos típicos. En estos ejemplos se ilustraran los pasos del modelado, es decir,
hacia un planteamiento matemático y su solución, y la interpretación física del resultado.

Se dedicará en este espacio la modelación de problemas que conduce a Ecuaciones

Diferenciales de segundo orden y esto lo justifica desde el punto de vista teórico y práctico
pues se verán más fáciles si uno se concentra primero en tales ecuaciones, pues de esta
manera los estudiantes familiarizado con los conceptos de segundo orden, resultaría más
fácil los conceptos, métodos y resultados hacia las de orden superior.

Escuela de Ingeniería Civil Páá giná 2

 ÍNDICE

INTRODUCCION 2

OBJETIVOS 4

- Aplicaciones de Ecuaciones Lineales a la Carrera 5

- Ley de Enfriamiento de Newton 11

- Crecimiento de Poblaciones 14

- Problema de Mezclas 17

Aplicaciones de Ecuaciones Diferenciales de Orden Superior 30

- Problemas Aplicados a los Resortes 43

Escuela de Ingeniería Civil Páá giná 3

 OBJETIVOS

 Mediante la primera y segunda ley de Hooke determinar ecuaciones diferenciales.

 Solucionar respectivamente los ejercicios usando las fórmulas adecuadas de la leyes de

Hooke y ecuación diferencial para hallar la función x(t).

 Hallar mediante la ecuación principal los valores y.

 Analizar y dar solución a los problemas planteados con referencia a las Ecuaciones
Diferenciales, pertenecientes a este tema de exposición.

Escuela de Ingeniería Civil Páá giná 4

 APLICACIONES DE LAS ECUACIONES LINEALES A LA
CARRERA
En esta sección vamos a considerar determinados fenómenos reales relacionados con diversos
ámbitos como la física, ingeniería, meteorología, sociología, etc., cuya modelización da lugar a una
ecuación diferencial de primer orden.

Para resolver las ecuaciones obtenidas aplicaremos los métodos estudiados.

TRAYECTORIA ORTOGONALES:

Dada una familia de curvas, consideramos el problema de encontrar otra familia de curvas que la
interceptan ortogonalmente en cada punto. Este problema aparece en el estudio de electricidad y
magnetismo y en la elaboración de cartas meteorológicas.
Consideramos una familia de curvas:

F (x, y) = C,

Siendo C un Parámetro. Derivando implícitamente esta ecuación obtenemos la ecuación diferencial

asociada a esta familia.

dF dF
dx+ dy = 0
dx dy
A partir de esta ecuación diferencial podemos obtener la pendiente de cada curva:

dF
dy dx
M= =-
dx dF
dy

Esto significa que las curvas ortogonales a la familia dada satisfacen la ecuación diferencial:

dF dF
dx- dy = 0
dY dx

Atendiendo al razonamiento anterior, dada una familia de curvas solución de la ecuación

diferencial:

M (x, y) dx + N (x, y) dy = 0
Planteamos la ecuación diferencial:
N (x, y) dx − M (x, y) dy = 0

Cuyas curvas solución son ortogonales a cada curva de la familia

Escuela de Ingeniería Civil Páá giná 5
 M (x, y) dx + N (x, y) dy = 0
Se dice que cada familia es una familia de curvas ortogonales a la otra.

Método de la resolución:

La familia de trayectorias ortogonales a una familia dada por la ecuación es la familia F (x, y) = C,
es la familia solución de la ecuación diferencial C.3.
Ejemplo:
Dada la familia de curvas x 2+ y 2 = C, veamos que la familia de trayectorias ortogonales a esta,
es la familia de rectas que pasan por el origen.

2 2
Calculamos la ecuación diferencial asociada a la familia x +y =C
2xdx + 2ydy = 0
Entonces la familia de trayectorias ortogonales satisface la ecuación:
2ydx − 2xdy = 0
Resolvemos la ecuación de variables separables.
X dy = y dx
1_dy = 1_dx
y x
Integrando:
ln |y| = ln |x| + C1

Por tanto:

|y| = |x| + e C 1 , y

Es decir:
y= C x, C = 0

Como y = 0 también es solución, la añadimos a la familia quitando la condición C=0

Por tanto, la familia de trayectorias ortogonales es:
y = C x, ∀ C ∈ R
Que presenta la familia de rectas que pasan por el origen:

Escuela de Ingeniería Civil Páá giná 6

 Familia de curvas Ortogonales

Hallemos la familia de trayectorias ortogonales a la familia de curvas definida implícitamente por:

4y + x2 + 1 − C e2y = 0.

Para hallar la ecuación diferencial asociada a esta familia, la reescribimos de la forma:

(4y + x2 + 1) e−2y = C
Y diferenciamos:

−2 y
2x e dx + (( 4 e−2 Y + X 2 +1 ) e−2Y (-2) dy)=0

Simplificando se tiene que la ecuación diferencial asociada a la familia dada es:

(−1 + 4y + x2 ) dx + xdy = 0 (2.35)

Resolvamos esta ecuación:

Sea M (x, y) = −1 + 4y + x2 y sea N (x, y) = x.

Podemos ver que no es una ecuación exacta ya que:

dM dN
=4= =1
dy dx

Pero admite un factor integrante de la forma u(x) puesto que la expresión.

dM dN
− 3
dY dx =
X
N

Solo depende de X. Por tanto:

3 dx
U(x) = ∫ 3

e x = e = X3
ln x

Multiplicando (2.35) por este factor integrante, obtenemos una ecuación exacta:

dF
=- x 3+ 4 x 2 y + x 5
dx

Escuela de Ingeniería Civil Páá giná 7

 df
= X4
dy

Comenzamos Integrando la segunda ecuación de este sistema respecto de y:

F(x,y) = ∫ X 4 dy = X 4 dy +ℵ ( I )

De donde despejamos:

ℵ ´ ( I ) =X 5− X 3

Integrando respecto de x:

I6 I4
(¿ X 5− X 3)dx= −
4 4
ℵ ( I )=∫ ¿

la solución de la ecuación diferencial exacta es:

X6 X4
I 4Y + − =c
6 4

Y la familia de curvas ortogonales viene dada por:

C 1 X2
y= + −
X4 4 6

PROBLEMAS DE APLICACIONES
EJERCICIO N° 01

Una columna perpendicularmente a una viga empotrada, donde la columna tiene una carga
5x
expresada con la siguiente función: 2y (ln +3seny) dy y la carga de la viga esta dada por:
x +3
3Y2
dX
X ( X +3 )

Determinar el coeficiente de elasticidad “c” en función de “x” y “y”

2
3Y
VIGA= dx
X ( X +3)
Escuela de Ingeniería Civil Páá giná 8
 5x
COLUMNA=2y (ln +3seny) dy
x +3

2
3y 5x
x ( x +3)
dx +2 y ln(x+3
+3 seny dy=0 )

6y
My (x, y)=
x ( x +3)

X +3 5 ( x+ 3 )−5 x 2y 15 6y
Nx(X, Y)= 2y ( x ¿= x =
5x ( x+3)
2
( x+ 3) 5 x x( x+3)

Luego la ecuación exacta es y:

2
F(x,y)= ∫ x (3x+3)
y 1 1
dx=3 y 2∫ ( −
1
3 x x+3 ) dx= y 2 ln
x
x+ 3
+ g ( y)

Como Fy=N

x 5x
Fy(x,y) = 2y ln +g´(y)= 2y( ln +3 seny ¿
x +3 x +3

Luego g´ (y) = 2 y ln5+6 yseny de donde

G (y)= y 2 ln 5+ 6∫ ysenydy

G (y)= y 2 ln 5+ 6 (−ycosy+ seny ) +C

Así, la solución general de la ecuación es:

5x
y 2 ln −6 ycosy+6 seny=C
x+ 3

EJERCICIO N° 02

Escuela de Ingeniería Civil Páá giná 9

 Un puente es diseñado para soportar una carga máxima: (2 xy dy ) . si el peso de la estructura es
( 2 y 2+3 x ) dx . Cual será la carga maxima “k”que puede soportar la estructura.

2
- Peso del puente :(2 y +3 x)dx
- Carga máxima: (2xy)dy

2
(2 y +3 x )dx +2xy dy= 0

M(x,y) = 2 ( y 2 )+ 3 x My ( x , y )=4 y

N(x,y)=2xy Nx(x,y)= 2y

My ≠ Nx ( No es exacta)

La ecuación no es exacta, por lo que buscamos un factor integrante:

dx
Mx−Nx 2 x 1 ∫
= = , solo depende de la x . Luego , e x =x es factor integrante
N 2 xy x

Ahora la ecuación (2x y 2+ 3 x 2 ¿ dx+ 2 x 2 ydy=0 es exacta

F(x,y)= ∫ 2 x 2 ydy=x 2 y 2+ h ( x ) y como Fx=M :

Fx(x,y)=2X y 2+ h ´ ( x ) =2 x y 2+ 3 x 2

Luego(x)= x 3+ K

Así, la ecuación la ecuación es

x 2 y 2 + x 3=K

EJERCICIO N° 03

UNA COMPUERTA FORMADA POR UNA PRESION: y ( 6 y 2−x−1 ¿ dx y su peso

2 xdy calcular el esfuerzo total “y”de toda la estructura:

y ( 6 y 2−x−1 ¿ dx +2 xdy =0

Escuela de Ingeniería Civil Páá giná 10

 −x +1 −3 3
La ecuación se puede escribir como y´ y= y , la que corresponde a una ecuación
2x x
de Bernoulli cuya solución es:

X +1 dx
x+1 dx e∫ 1 1 x
−2
−∫ x , es decir, 2 = x (C+6 e )
y =e x
(C +6 ∫ ) y x e
x

x ex
Por lo tanto, La solución de la ecuación es y= ± √
C+6 e x

LEY DE ENFRIAMIENTO DE NEWTON

La ley de Newton de enfriamiento menciona que la velocidad con que se enfría un objeto es
proporcional a la diferencia entre su temperatura y la temperatura del medio donde se encuentra.
dT
Supongamos que T representa la temperatura de un objeto en cual instante “t”, entonces es
dt
la velocidad con que se enfría el objeto, y está dado por:

dT
=k (T −T m)
dt

Donde k es una constante de proporcionalidad y T m es la temperatura del ambiente.

EJERCICIO N° 01

La temperatura de una muestra de suelo recién sacada del horno es de 110°C. Un minuto después se
ha enfriado a 100°C en un laboratorio que está a 20°C ¿Cuánto debe ser el período que debe
transcurrir antes de que la muestra de suelo alcance una temperatura de 60°C?

t o :T 1=110 ° C

Escuela de Ingeniería Civil Páá giná 11

 t 1 :T 2=100° C

T m = 20 °C

T x =60 ° C

dT
=k (T −T m)
dt

dT
=k (T −20)
dt

dT
∫ T −20 =∫ kdt
T
ln ⁡( T −20)¿ =kt

ln ( T −20 ) −ln ( 110−20 )=kt

De lo cual hallamos el valor de la constante: ln 90

ln ( T −20 ) −ln 90=kt

ln ( T −20 ) =kt+ ln ⁡( 90)

Reemplazamos que en el instante 1 la temperatura será de 100°C.

ln ( 80 ) =kt+ ln ⁡(90)

8
k =ln
9

Reemplazar el k en la ecuación anterior para una temperatura de 60°C

8
ln ( T −20 ) =ln ⁡( )t +ln ⁡( 90)
9
t
8
T −20=130 ()
9
t
8
60=130 ()
9
+20

Escuela de Ingeniería Civil Páá giná 12

 4
ln ⁡( )
13
t=
8
ln ⁡( )
9

t=10 min .
EJERCICIO N° 02

Una varilla de acero corrugado a una temperatura de 100°F se pone en un cuarto a una temperatura
constante 0°F.
Después de 20 minutos la temperatura de la barra es de 50°F.

a. ¿Cuánto tiempo tardara la barra para llegar a una temperatura de 25°F?

b. ¿Cuál será la temperatura de la barra después de 10 minutos?

a. Tiempo en 25°F

t o :T 1=100 ° F

T m = 0 °C

T x =25 ° F

dT
=k (T −T m)
dt

dT
=kT
dt

dT
∫ T ∫
= kdt

ln ⁡( T )¿T =kt

ln ( T )−ln ( 100 )=kt

El valor de la constante es ln (100) y reemplazamos el T para el instante t= 20min.

ln ( 50 )=k ( 20 ) +ln ⁡(100)

k ( 20 )=ln ( 12 )
k =−0.03465

Reemplazamos en la ecuación original:

Escuela de Ingeniería Civil Páá giná 13

 ln ( 50 )=−0.03465 ( 20 ) + ln ⁡( 100)

0.03465 t=4.605−3.219

t=40 min .

b. Temperatura de la barra después de 10 min

ln ( T )=kT +ln ( 100 )

ln ( T )=lnkT +ln ( 100 )

ln [ ln ( T ) ]=lnkT (100)
kt
T =100 e

T =100 e−0.3465

T =71° F

CRECIMIENTO DE POBLACIONES

Para estudiar el crecimiento de poblaciones se pueden seguir diferentes modelos. Por ejemplo, el
modelo maltusiano o exponencial se basa en que la tasa de crecimiento es proporcional a la
población; por tanto, esta descrito por una ecuación diferencial de la forma:

Un modelo de poblaciones más realista es el modelo logístico donde se supone que existe una tasa
de mortalidad debida a factores externos, llegándose a la ecuación diferencial de la forma:

EJERCICIO N° 01:

En 1790 Estados Unidos tenía una población de 3.93 millones de personas y en 1800 de 5.31
millones. Usando el modelo exponencial, estimemos la población de EEUU en función del tiempo.
Escuela de Ingeniería Civil Páá giná 14
 Sea x (t) la población de EEUU en un instante t (años). La solución de la ecuación diferencial
separable

Para el instante inicial t = 0, correspondiente al año 1790, se tiene que

X (0) = 3,93 millones. Para el instante t = 10, correspondiente al año 1800,

X (10) = 5,31 millones. Estos datos nos permiten calcular las constantes C y k:

Por tanto, la expresión que nos da la población de EEUU en función del tiempo es:

EJERCICIO N° 02:

Se sabe que la población de cierto país aumenta de una forma proporcional al número de habitantes
actuales. Si después de dos años la población se ha duplicado y después de tres años la población es
de 20.000 habitantes, hallar el número de habitantes que había inicialmente en el país.

X (t) = población en el instante t

dx
=kx ( t )
dt

dx
=kdt
x (t )

ln ( x (t ) ) =kt+ C

x ( t )=C ekt

x ( 0 )=x o

Escuela de Ingeniería Civil Páá giná 15

 x ( 0 )=x o ; x o =C

x ( 2 ) =2 x o
2k
2 x o =C e

2 k=ln|2|

k =0.346574

x ( 3 )=20.000 20.000=C e3 k
3 x 0.346574
20.000=x o e

: x o=7.071 habitantes

EJERCICIO N° 03:

Un cultivo tiene una cantidad inicial No de bacterias, cuando t = 1h. La cantidad medida de
bacterias es 3/2 No. Si la razón de reproducción es proporcional a la cantidad de bacterias
presentes, calcule el tiempo necesario para triplicar la cantidad inicial de los microorganismos.

dN
Se resuelve la E.D: =kN ( 1)sujeta a N (o)=No
dt
3
A continuación se define la condición empírica N (1 )= No , para hallar K .
2
dN
Tomemos (1) =KN
dt
dN
Pasando términos −k =0 que es una ecuaciónlineal con q(x )=0.
dt
Separando variables:

dN dN
=kN =kdt
dt N

dN
∫ N ∫
= kdt +C

Escuela de Ingeniería Civil Páá giná 16

 Dónde: N ( t ) :C e kt
Cuando:t=0 N ( o )=C e k 0=C
Pero N ( o)=No=C

Y por lo consiguiente: N ( t )=No . ekt

3
Cuandot=1 → No=N o . ek
2
3 3
=e k entonces :k =ln =0.4055
2 2

0.4055 t
Así: N ( t ) =N o .e

Para establecer el momento en que se triplica la cantidad de bacterias, despejamos t de:

3 No=No. e 0.4055t
0.4055 t=ln 3

ln ⁡(3)
t= =2.71hrs
0.4055

PROBLEMAS DE MEZCLAS

Al estudiar las aplicaciones de las ecuaciones diferenciales de primer orden vimos como modelizar,
mediante una ecuación diferencial, la velocidad de cambio de una sustancia disuelta en un líquido
contenido en un tanque, en el cual entraba un fluido con una cierta concentración de dicha sustancia
y donde la mezcla fluida hacia fuera del tanque.

dx
Recordemos que si x (t) es la cantidad de sustancia presente en el tanque en el instante t y es
dt
la rapidez con que x cambia respecto al tiempo, la ecuación diferencial que modeliza este problema
viene dada por:
dx
=Ve−Vs
dt

Dónde:

Escuela de Ingeniería Civil Páá giná 17

 [ ][ ]
velocidad de entrada concentraciòn al
cantidad
Ve ( t
=) del fluido
vol
( )
t
×
entrar
cantidad
(
vol )

[ )][ ]
velocidad de salida concentraciòn al
cantidad
Vs ( t
=) (
del fluido
vol
t
×
salir(cantidad
vol )

Siendo la concentración de salida, la cantidad de sustancia x(t) dividida por el volumen total en el
tanque en dicho instante t.

Ahora, vamos a considerar varios depósitos interconectados entre sí, de modo que se obtiene un
sistema de ecuaciones diferenciales lineales.

EJERCICIO N° 01

En dicha obra se cuenta con un cilindro que contiene 450 lt de líquido, que será empleado para el
curado de un concreto fabricado con cemento portland tipo 5 y además no cuenta con armadura. En
el cilindro se disuelven 30 gr de sal: Una salmuera que contiene 3 gr/lt se bombea al cilindro con
una intensidad de 6 lt/min, la solución adecuadamente mezclada se bombea hacia fuera con una
intensidad de 8 lt/min. Encuentre el número de gramos de sal y la concentración de sal, que hay en
el cilindro en un instante cualquiera.

Solución:

El volumen inicial de líquido en el cilindro es V 0 = 450 lt y la cantidad inicial de sal en el

tanque es x0 = 30 gr.

La salmuera que se bombea al tanque tiene una concentración C1 = 3 gr/lt y se bombea a una
razón Q1 = 6 lt/min. La solución, debidamente agitada y homogeneizada, se extrae del tanque a
razón Q2 = 8 lt/min.

La ecuación diferencial asociada a los

problemas de mezcla es

dx Q2
+ x =Q1 C (1)
dt v 0 +( Q1−Q 2)

Sustituyendo los datos en la ecuación (1)

Escuela de Ingeniería Civil Páá giná 18

 dx 8
+ x=18
dt 450+(6−8)t

Simplificando

dx 4
+ x=18
dt 225−t

dx
Despejando:
dt

dx 4
+ x (2)
dt 225−t

dx dx
Ya que la diferencial de la cantidad x de sal es dx = dt, sustituyendo dada por la
dt dt
ecuación (2):

4
dx= 18−( 225−t
x dt )

Reordenando la ecuación

4
dx + xdt=18 dt (3)
225−t

La ecuación (3) es una ecuación diferencial lineal de la forma x’ + F(t) x = G(t),

4
Donde F(t) = x , G ( t ) =18.
225−t

Para resolver la ecuación (3) debe determinarse un factor integrante μ = ℮∫ F (t )dt

4
μ = ℮∫ F (t )dt =℮∫ 225−t dt =℮−4 ln|225−t|= (225−t )−4

Multiplicando la ecuación (3) por el factor integrante

( 225−t )−4 dx + 4 ( 225−t )−5 x dt =18 ( 225−t )−4 dt (4)

Escuela de Ingeniería Civil Páá giná 19
 Puesto que

( 225−t )−4 dx + 4 ( 225−t )−5 x dt=d [ ( 225−t )−4 x ]

Sustituyendo en la ecuación (4)

d [ ( 225−t )−4 x ]=18 ( 225−t )−4 dt

Integrando

∫ d [ ( 225−t )−4 x ]=18 ∫ ( 225−t )−4 dt(5)

Ambas integrales son inmediatas

∫ d [ ( 225−t )−4 x ]= (225−t )−4 x+ k 1

( 225−t )−3
∫ ( 225−t )−4 dt= 3
+k 2

Sustituyendo los resultados de las integrales en la ecuación (5)

( 225−t )−4 x=6+ ( 225−t )−3 +k (6)

Para determinar el valor de la constante k de integración se utiliza la condición inicial x(0) = 30,
esto es, t0 = 0 min y x0 = 30 gr se sustituyen en la ecuación (6)

( 225 )−4 30=6+ ( 225 )−3 +k

Despejando k:
−4 −3
k =( 225 ) 30−6+ ( 225 )
−4
k =( 225 ) ( 30−6(225) )

k =( 225 )−4 (30−1350)

−4
k =−( 225 ) 1320

Este valor obtenido para k se sustituye en la ecuación (6)

( 225−t )−4 x=6+ ( 225−t )−3−( 225 )−4 1320

Multiplicando por ( 225−t )−4

Escuela de Ingeniería Civil Páá giná 20

 4
225−t
x(t)= 6 ( 225−t )−1320 ( 225
(7))
La ecuación (7) representa la ley de variación de la cantidad de sal en el cilindro en cualquier
instante t.

Para determinar la concentración de sal en el tanque en un instante t cualquiera, se debe

recordar que la concentración C(t) es el cociente entre la cantidad de sal y el volumen de
líquido en el tanque, en un instante t cualquiera, es decir

x (t )
C(t) = (8)
v (t)

Donde

V(t) = V0 + ( Q1 – Q2 ) t = 450 – 2t = 2 ( 225 – t ) (9)

Sustituyendo las ecuaciones (7) y (9) en la ecuación (8)

4
225−t
C(t) =
6 ( 225−t )−1320 (225 )
=3−690
( 225−t )−3
2(225−t) ( 225 )4

( 225−t )−3
C(t) =3−690 4
(10)
( 225 )

La ecuación (10) representa la ley de variación de la concentración de sal en el cilindro en

cualquier instante t.

EJERCICIO N° 02

Tenemos un par de tanques completamente interconectados A y B, que inicialmente tiene 50

galones de agua donde en el tanque A inicialmente se disolvieron 25lb y en tanque B tenemos
agua pura de dónde;

X1(0) = 25

X2(0) = 0

Plantear un sistema de ecuaciones que nos permita encontrar cuantas son las cantidades de sal
que hay en el tanque A y B.

Escuela de Ingeniería Civil Páá giná 21

 Solución:

gal X 2 gal X 1
X1`(t) = ( gal
3
min
×0
lb
+1 ×
gal min 50
− 4 )(×
min 50 )
X2 2 X1
X1`(t) = −
50 25

X2`(t) = ( 4
gal X 1 lb
×
min 50 gal ) ((
− 1
gal X 2 lb
×
min 50 gal
+3
gal X 2 lb
) ×
min 50 gal )
2 X1 2 X2
X2`(t) = −
25 50

Ordenando tenemos:
+2 X 1 X2
D X1 – =0
25 50
+2 X 2 2 X 1
D X2 −
50 25
2 X
⟹ D+( 25 )
X 1− 2 =0
50
−2 X 1 2 X2
D+ =0
25 25
Resolvemos por la regla Cramer; para el tanque A = X1:

Escuela de Ingeniería Civil Páá giná 22

 | |
2 −1
D+
25
−2
50 X = 0
D+
2 1 0 ||
25 25

, y por la regla de Surrus se obtiene;

| || |
2 −1 −1
D+ 0
25 50 X = 50
−2 2 1 2
D+ 0 D+
25 25 25

[( )D+
2 2 1
25
− X =0
625 1 ]
[ D 2+
4D 4
+ −
1
X =0
25 625 625 1 ]
[ D 2+
4D 3
+
]
X =0
25 625 1

[ 625 D2+100 D+3 ] X 1 =0

Resolvemos por la regla Cramer; para el tanque B = X2:

| || |
2 −1 2
D+ D+ 0
25 50 X = 25
2
−2 2 −2
D+ 0
25 25 25

[( )D+
2 2 1
25
−
]
X =0
625 2

[ D 2+
4D 3
+
]
X =0
25 625 2

[ 625 D2+100 D+3 ] X 2 =0

Tenemos;

[ 625 D2+100 D+3 ] X 1 =0

[ 625 D2+100 D+3 ] X 2 =0
2
para X 1=625 m +100 m+ 3

Escuela de Ingeniería Civil Páá giná 23

 X 1= (25 m )2 +4 ( 25 m) +3=0

(25+3) (25m+1) = 0

25m = -3 v 25m = -1

Entonces:
−3
m1=
25
1
m2=
25
Dónde:
m1 t m2 t
X 1=C 1 ℮ + C2 ℮
−3 −1
t t
X 1=C 1 ℮ 25 +C 2 ℮ 25
−3 −1
t t
X 2=C 3 ℮ 25 +C 4 ℮ 25

Como las ecuaciones diferenciales X 1 y X 2 , son las mismas ecuaciones auxiliares y se tiene;

X2 2 X1
X1`(t) = −
50 25
−3 −1 −3 −1 −3 −1
−3 t 1 t 1 t 1 t 2 t 2 t
C ℮ 25 − C ℮ 25 = C 3 ℮ 25 + C 4 ℮ 25 − C 1 ℮ 25 − C 2 ℮ 25
25 1 25 2 50 50 25 25

Agrupamos e igualamos a los del mismo exponente:

−3 1
C 1= C 3
25 50

C3 =50 ( −125 )C 1

C3 =−2 C 1

1 1
C2 = C4
25 50

C 4=50 ( 251 )C 2

C 4=2 C2

Escuela de Ingeniería Civil Páá giná 24

 Reemplazamos en las ecuaciones diferenciales y así reducirlas a dos constantes para que
cumpa con la ecuación diferencial original.
−3 −1
t t
X 1=C 1 ℮ 25 +C 2 ℮ 25
−3 −1
t t
X 2=−2 C1 ℮ 25 +2C 2 ℮ 25

De la condición inicial se tiene;

X1(0) = 25

X2(0) = 0

Entonces:

Cuando; X1(0) = 25
−3 −1
(0) (0)
25= C ℮ 25 +C ℮ 25
1 2

25= C1 ℮(0) +C 2 ℮(0)

25= C1 +C 2

Cuando; X2(0) = 0
−3 −1
(0) (0)
0= −2 C ℮ 25 +2 C ℮ 25
1 2

25= −2 C 1+2 C 2

Dónde:

C1 +C 2=25

−2 C 1+2 C 2=2

por 2 :

2C 1 +2C 2=50

−2 C 1+2 C 2=2

⇒ 4 C2=50

25
C2 =
2

Escuela de Ingeniería Civil Páá giná 25

 Reemplazando:

25
C1 =25−
2

25
C1 =
2

Obtenemos:
−3 −1
25 t 25 t
X 1= ℮ 25 + ℮ 25
2 2
−3 −1
25 25 t 25 25 t
X 2= ℮ + ℮
2 2

X 1=25 lb

X 2=25 lb

En un t=0 tenemos que en el tanque A y B entra la misma cantidad de sal = 25lb.

EJERCICIO N° 03

Consideremos dos tanques interconectados, conteniendo 1000 litros de agua cada uno de
ellos. El líquido fluye del tanque A hacia el tanque B a razón de 20 l/min y de B hacia A a razón
de 10l/min. Además, una solución de salmuera con una concentración de 2 kg/l de sal fluye
hacia el tanque A a razón de 20 l/min, manteniéndose bien agitado el líquido contenido en el
interior de cada tanque. La solución diluida fluye hacia el exterior del sistema, desde el tanque
A a razón de 10 l/min y desde el tanque B también a razón de 10 l/min. Si inicialmente el
tanque B solo contiene agua y el tanque A contiene 40 kg de sal, calculemos la concentración
de sal en el tanque B al cabo de 10 min.

Escuela de Ingeniería Civil Páá giná 26

 Solución:

Sea x1(t) la cantidad (kg) de sal en el tanque A en un instante t y sea x2(t) la cantidad (kg) de sal
en el tanque B en un instante t.
Como sabemos, la velocidad de cambio de la sustancia en un tanque para un tiempo t, debe
ser igual a la velocidad a la que dicha sustancia entra en el tanque menos la velocidad a la que
lo abandona, es decir,

Ve ( cantidad
t )−Vs (
cantidad
t )
Dónde:

[ ][ ]
velocidad de entrada concentraciòn al
cantidad
Ve ( t
=) del fluido
vol
( )
t
×
entrar
cantidad
(
vol )

[ )][ ]
velocidad de salida concentraciòn al
cantidad
Vs ( t
=) (
del fluido
vol
t
×
salir(cantidad
vol )

Siendo la concentración de salida, la cantidad de sustancia x(t) dividida por el volumen total en
el tanque en dicho instante t.

En este caso, tendremos:

X (t) kg X (t) kg
X1`(t) = ( 20
l
min
kg
×2 +10
l
l
× 2
min 1000 l
− 20+10
l
min
× 1 ) ((
1000 l ) )
Escuela de Ingeniería Civil Páá giná 27
 X 2 (t) 3 X 1 (t)
X1`(t) = 40 + −
100 100

X (t) kg X (t) kg
X2` (t) = ( 20
l
× 1
min 1000 l ) ((
− 10+10
l
min
× 2
1000 l ) )
2 X 1 (t ) 2 X 2 (t)
X2`(t) = −
100 100

Ordenamos las ecuaciones diferenciales:

−3 X 1 (t ) X 2 (t)
X1`(t) = + +40
100 100

X 1 (t) X 2 (t)
X2`(t) = −
50 50

3 X 1 (t ) X 2 (t)
DX1 ¿− + +40
100 100
X (t) X (t)
DX2 ¿ 1 − 2
50 50

X2
( D+ 1003 ) X − 100 =40
1

−X 1 1
50 (
+ D+ X =0
50 2 )
Obtenemos, por tanto, por la regla de Cramer el sistema de ecuaciones diferenciales:

| ||
3 −1
D+
−2
100 100 = x (t) + 40
D+
2 y (t ) 0 || |
100 100

Escuela de Ingeniería Civil Páá giná 28

 |( D 2+
5D
+
3
+ )1
100 5000 5000 |
|D + 100
5D
2
+
10000|
4

Por fórmula general obtenemos las raíces:

m1 ¿−0.01
m2 = −0.04

Y a continuación los vértices asociados:

| | | | ||
−3
−(−0.01) 0,01
H-0,01 = ( V V
1, 2 )= 100 × V1 = 0
0
−(−0.01) V 2
−1
0.02
50

| | | | ||
¿ −0,02 0,01 × V 1 = 0
0.02−0,01 V 2 0

¿−0.02V 1 +0,01 V 2=0

V 2=2 V 1

| | | |||
−3
−(−0.04)0,01
H-0,04 = ( V 1 , V 2 )=
100 × V1 = 0
0
−(−0.04) V 2
−1
0.02
50

| | | | ||
¿ 0,01 0,01 × V 1 = 0
0.02 0,02 V 2 0

¿ 0.01V 1+0,01 V 2=0

V 1+ V 2=0

Podemos tomar ⃗ ()
u1= 1 y ⃗
2 ( )
u2= 1 como vectores propos asociados a m1 y
−1
m2,
respectivamente. Por lo tanto la solución de la parte homogénea es:

X(t) = C1 e
m t 1
1

2 ()
+C2 e m t 1
2

−1( )
X(t) = C1 e
−0.01t 1
2 () ( )
+C2 e−0.04 t 1
−1

Para hallar una solución particular podemos aplicar el método de los coeficientes
indeterminados. Puesto que la parte no homogénea es constante, es decir, es un polinomio de
grado 0 y además 0 no es valor propio, suponemos que la solución tiene la forma:
Escuela de Ingeniería Civil Páá giná 29
 x 1(t) =⃗a =
⃗
() a1
a2

Derivando y sustituyendo en el sistema, tenemos:

−3 1

()
0
0
=
( )( ) | |
100 100 a1 + 40
2 −1 a2
100 50
0

Haciendo las correspondientes operaciones tenemos:

−3 1

( )( ) |
100 100 a 1 = −40
2 −1 a 2
100 50
0 |

De donde se obtiene: a1 = a2 = 2000.

Por tanto, la solución general del sistema es:

X(t) = ( x(t)
y (t ))=C e
1 ( )
1 +C e
−0.01 t
2 (2
−1 ) 2000
1 + 2000
−0.04 t

Sustituyendo las condiciones iníciales x(0) = 40 e y(0) = 0 obtenemos C1 y C2 :

(400)=C (12 )+C (−11 )+2000

1 2
2000

C
( )( )
¿ 1 1 + 1 + 2000
2 −1 C 2 2000

(12−11 )(CC )=(−1960

1

2
2000 )

De donde se obtiene: C1 = −1320 y C2 = −640.

Por tanto, la solución de este problema de valor inicial es:

Escuela de Ingeniería Civil Páá giná 30

 X(t) = ( x(t)
y (t )) (
−0.01t
=−1320 e
2 )
1 −640 e
( −1 ) 2000
1 + 2000
−0.04 t

Por lo tanto, en un tiempo de 10 minutos se tendrá:

X(t) = ( x(10)
y (10) )=−1320 e ( 1)−640 e ( 1 ) + 2000
−0.1
2
−0.4
−1 2000

La cantidad de sal en el tanque B será:

Y(10) = −1320 e−0.1 2 − 640 e−0.04 t (−1) + 2000 = 237,004 kg

y la concentración de sal en B al cabo de 10 minutos será:

237 . 004 kg
= 0,2370 kg/l.
1000 l

APLICACIONES DE ECUACIONES DIFERENCIALES EN INGENIERÍA

CIVIL

Las ecuaciones diferenciales son muy interesantes en cuanto a la posibilidad que presentan para
indagar sobre variedad de problemas de las ciencias físicas, biológicas y sociales. A partir de la
formulación matemática de distintas situaciones se describen procesos reales aproximados.

Dentro de los diversos campos de acción de la ingeniería civil, una de las múltiples aplicaciones
de ecuaciones diferenciales está relacionada con el estudio de las flexiones, un ejemplo es:

FLEXION DE UNA VIGA EN VOLADIZO PARA PEQUEÑAS FLEXIONES:

Escuela de Ingeniería Civil Páá giná 31

 Una viga o una barra delgada son sólidos
homogéneos e isótropos cuya longitud
es grande comparada con las dimensiones
de su sección trasversal.

Este ejemplo, nos permite practicar con procedimientos numéricos aplicados al

 Cálculo de la raíz de una ecuación.

 Integral definida.

Supongamos que

 La barra tiene una longitud L mucho mayor que las dimensiones de su sección trasversal, y
que la deformación debida a su propio peso es despreciable.
 Que la sección de la barra no cambia cuando se dobla. Cuando el espesor de la barra es
pequeño comparado con el radio de curvatura, la sección trasversal cambia muy poco.
En estas condiciones es aplicable la ecuación de Euler-Bernoulli que relaciona el momento
flector M de la fuerza aplicada y el radio de curvatura ρ de la barra deformada

COLUMNA DE EULER

Es interesante mencionar un tema muy aplicativo a columnas la cual se resuelve según Euler de la
siguiente forma:

Escuela de Ingeniería Civil Páá giná 32

 Consideremos una columna esbelta, de longitud L y sección transversal constante, cuyos extremos
están restringidos a mantenerse sobre una misma línea vertical pero que puede girar libremente.

La columna se encuentra sometida a una carga vertical P. Si esta carga sobrepasa cierto valor, la
columna resistirá sin experimentar deformaciones transversales. Cuando sobrepasa determinada
carga la columna se abomba y se producen deformaciones laterales.

El problema es una aplicación de la ecuación diferencial, pues en la sección SS de la columna se

d2 y
produce un momento dado por –PI el cual es igual a EI , esto es:
d x2

2
d y
M =EI
d x2

O también:

x´´ + k² x = 0

PROBLEMA N°1 DE APLICACIÓN

Escuela de Ingeniería Civil Páá giná 33

 Determinar la carga crítica para una barra delgada articulada en los extremos, cargada con
una fuerza de compresión axial en cada extremo. La línea de acción de las fuerzas pasa por el
centro de gravedad de la sección de la barra.

La carga crítica se define como la fuerza axial suficiente para mantener a la barra en una forma
ligeramente deformada. Bajo la acción de la carga P, la barra tiene la forma flexada representada en
la figura.

Para que se produzca la flexión lateral es necesario, individualmente, que un extremo de la barra
pueda moverse axialmente respecto al otro. La ecuación diferencial de la curva deformada es:

d2 y
M =EI → Ecuación 1
d x2

Aquí, el momento flector en el punto A de coordenadas (x,y) no es más que el momento de la fuerza
P aplicada en el extremo izquierdo de la barra, respecto a un eje por el punto A perpendicular al
plano.

El momento flector es:

M= -Py

Reemplazamos en 1

d2 y
−Py=EI → Ecuación 2
d x2

Si hacemos que:

P 2
=K → Ecuación 3
EI

Esta ecuación se transforma en:

d2 y 2
2
+k y=0 → Ecuación 4
dx

Resolviendo la ecuación:

Escuela de Ingeniería Civil Páá giná 34

 '' 2
y +k y=0 → Ecuación Homogénea de Coeficientes Constantes

2 2
Ecuación Auxiliar: m + k =0

2 2
m =−k

CASO III

p= 0

q = ki

Soluciones Particulares

px px
y 1=e cos ⁡( qx ) y 2=e sin ⁡( qx )

y 1=e 0 x cos ⁡( kx ) y 2=e 0 x sin ⁡(kx)

y 1=cos ⁡( kx) y 2=sin ⁡( kx )

Solución General

( qx )
C1 cos ( qx ) +C 2 sin ¿
px
y=e ¿

C1 cos ( kx ) +C 2 sin ( kx )
]
y=e0 x ¿

y=C 1 cos ( kx ) +C2 sin ( kx ) → Ecuacion 5

Determinamos C1 y C 2 . En el extremo izquierdo de la barra, y= 0 cuando x= 0, sustituyendo los

valores en ecuación 5 se tiene:

C1 =0

En el extremo derecho de la barra, y= 0 cuando x= L, sustituyendo los valores en la ecuación 5 se

tiene:

0 = C sinkL

Escuela de Ingeniería Civil Páá giná 35

 Evidentemente, C= 0 o sinkL= 0, pero si C= 0y es nulo en todos los puntos y tenemos solamente el
caso trivial de una barra recta, que es la configuración anterior a producirse el pandeo. Como esta
solución no es de nuestro interés, tomaremos:

sinkL = 0 → Ecuación 6

Para que sea cierto debemos tener:

kL= nπ radianes (n=1,2,3,….) → Ecuación 7

P
Sustituyendo =K 2 en la ecuación 7 obtenemos:
EI

√ P
EI
nπ

2 2
n π EI
P= → Ecuación 8
L2

Indudablemente, el menor valor de esta carga P corresponde a n= 1. Entonces tenemos el primer

modo del pandeo, en que la carga crítica está dad por:

π 2 EI
Pcr = → Ecuación 9
L2

Es la llamada carga de pandeo de Euler para una columna con extremos articulados. La forma
flexada correspondiente a esta carga es:

y=C sin
√ P
EI
x → Ecuación 10

Sustituyendo en esta ecuación el valor de (9), obtenemos:

πx
y=C sin → Ecuación 11
L

Por tanto, la deformación es un sinusoide. A causa de la aproximación adoptada en la deducción de

la ecuación 1, no es posible obtener la amplitud del pandeo, representada por C en la ecuación 11.

Como se puede verse en la ecuación 9, el pandeo de la barra se produce respecto al eje de la sección
para el cual I adopta un valor mínimo

Escuela de Ingeniería Civil Páá giná 36

 PROBLEMA N°2 DE APLICACIÓN

(METODO DE LA DOBLE INTEGRACION)

Una viga ligera de sección uniforme y longitud l esta soportada en los extremos y colocada
horizontalmente. Los extremos están empotrados horizontalmente y la viga se carga de manera que
la intensidad de carga aumenta uniformemente desde cero en un extremo hasta w en el otro. Hállese
una expresión para la flecha de la viga en cualquier punto de su longitud.

La carga a la distancia x de un extremo es wx / l y se tiene la ecuación diferencial Y

d 4 y wx
EI = = flecha a distancia x del extremo izquierdo
dx 4 l

Y
X

Integrando se obtiene

d 3 y wx 2
EI 3 = +A,
dx 2l

d 2 y wx3
EI = + Ax+ B ,
dx 2 6 l

4
dy wx 1
EI = + A x 2+ Bx+C ,
dx 24 l 2

Escuela de Ingeniería Civil Páá giná 37

 5
wx 1 1
EI y= + A x 3 + B x 2+ Cx+ D ,
120 l 6 2

Las condiciones iníciales que debe satisfacer la solución son

Y =0 cuando x=0,

y=0 cuando x=l ,

dy
=0 cuando x=0,
dx

dy
=0 cuando x=l ,
dx

En que las dos últimas condiciones expresan el hecho de que la viga es horizontal en los extremos.
Las condiciones (i) y (iii) hacen C=D=0, y (ii) y (iv) dan

1 1
w l2 + A l 3 +Bl=0,
24 2

1 1 1
w l 4 + A l 3 + Bl 2 =0,
120 6 2

Entonces, A=−3 wl /20, B=wl 2 /30 y se obtiene la solución

120 lEIy=w x2 ( x−l ) ( x2 +lx−2l 2 ) .

PROBLEMA N°3 DE APLICACIÓN

Utilizando las EDO de orden superior en el siguiente caso típico de carga de viga de concreto y
sección uniforme rectangular, cargada uniformemente con 2q ton/ml.

X =0
Y =0
Escuela de Ingeniería Civil Páá giná 38
 a) Determinar la ecuación analítica de la curva elástica.

b) La flecha máxima.

SOLUCION

a) Ecuación de curva elástica

a.1) Resolución estática de la viga.0

a.2) momento flector en sección 1-1. (DCL).

a.3) Ecuación diferencial de la elástica de una viga. El producto el que se llama rigidez a la
flexión, es normal. Constante a lo largo de la viga.

d2 y M
=
d x2 El

a.4) EDO de la elástica para la viga del presente caso, reemplazando el valor de M.
2
d y
El 2 = qLx – qx2
dx

a.5) Resolviendo la EDO de orden superior anterior, integrando sucesivamente dos veces

Integrando 2° vez
Escuela de Ingeniería Civil Páá giná 39
 a.6) Cálculo de constantes A, B.

Utilizando condiciones /restricciones de, extremo de viga. (La flecha o deformación en

apoyos flecha Y= 0 en apoyo izquierdo)

Flecha y=0 en apoyo derecha de viga

También se calcula el valor de A

Con la condición x=L/2 , y´=0; es decir que la tangente trigonométrica. Del ángulo θ que
hace la tangente geométrica a la elástica en su centro o cundo x=L/2 , es cero, o que y
´=dy/dx=0 en el punto de la mayor curva tiene de la elasticidad.

Se obtiene

a.7) Reemplazando valor de constantes en (∝) se obtiene la ecuación analítica de la curva

elástica.

Escuela de Ingeniería Civil Páá giná 40

 Asumiendo que Y es + hacia abajo:

qx
EI y = (L3 - 2Lx2 + x3)
12
………β

b) Flecha máxima (δ)

b.1) Valor máxima de la función β , ecuación de la curva elástica.

Razonamiento una función Iy=f(x) tiene un valor máximo /mínimo para y´=f(x)=0

Simplificando:

5 q L4
Ymax = ………..
192 EI
5 q L4
δ=
192 EI

Escuela de Ingeniería Civil Páá giná 41

Determinar la ecuación de la flecha “y” de la línea central (elástica)

R= qL
2
qL
M=
2

x=0

y=0
2 2
∑ M 1−1=Rx− qL2 − qx2 =0

x=0

Y’=0

Escuela de Ingeniería Civil Páá giná 42

 Ley de equilibrio: sumatoria de momentos de las cargas a la derecha respecto a la sección 1-1, es
Según “la ecuación diferencial de la elástica de una viga es:

d2 y
EI 2 =M
dx

2 2
qL qx
Se ha deducido que: M =qLx− −
2 2

d2 y qL2 qx 2
Se obtiene que: EI =qLx− −
dx 2 2 2

La EDO resultante es una EDO de orden superior de la forma:

dn y
n
=f (x)
dx

La solución de EDO de orden superior se obtiene por integración sucesiva

dy qLx2 qL2 x qx 3
EI ( )
dx
=
2
−
2
−
6
+A

3 2 2 4
qLx qL x qx
EIy= − − + Ax+ B … ..(∝)
6 4 24

Calculando constantes A, B en base a condiciones de extremo de viga en voladizo.

x=0

y=0

x=0

Y’=0

A=0

Escuela de Ingeniería Civil Páá giná 43

 B=0

Reemplazando Constantes En (∝)

qLx 3 qL2 x 2 qx 4
EIy= − − … …(∝)
6 4 24

La EDO de la flecha y:

3 2 2 4
qLx qL x qx
EIy= − − … …(∝)
6 4 24

Expresión de la flecha máxima ( δ )

En la viga en cantiliever/voladizo la mayor deformación de la elástica o la mayor flecha se presenta

en el extremo derecho; es decir matemáticamente cuando x es igual a L.

qL4 qL4 qL4

EIy= − −
6 4 24

4
qL
EIy=¿ -
8 EI

−q L 4
Luego δ=
8 EI

Escuela de Ingeniería Civil Páá giná 44

 PROBLEMAS APLICADOS A LOS RESORTES

1.1 Sistemas de resorte y masa:

 Movimiento Libre No Amortiguado:

Movimiento libre no amortiguado - Ley de Hooke Supongamos que, como en la figura (b), una
masa m1 está unida a un resorte flexible colgado de un soporte rígido. Cuando se reemplaza
m1 con una masa distinta m2 , el estiramiento, elongación o alargamiento del resorte
cambiará.

 Según la “Ley de Hooke”

Escuela de Ingeniería Civil Páá giná 45

 El resorte mismo ejerce una fuerza de restitución F, opuesta a la dirección del alargamiento y
proporcional a la cantidad de alargamiento S.
En pocas palabras decimos que, F = KX, donde “k” es una constante de proporcionalidad llamada
constante del resorte.

Aunque las masas con distintos pesos, estiran un resorte en cantidades distintas, este está
caracterizado esencialmente por su número “k”.

Por ejemplo:

Si una masa que pesa 10 libras, estira ½ pie en un resorte, utilizamos la ley de Hooke F = KX y
entonces reemplazamos: 10 = k (1/2), esto implica que k = 20 lb/pie.
Entonces, necesariamente, una masa cuyo peso sea de 8 libras estirará el resorte 2/5 de pie.

 Segunda Ley de Newton

Después de unir una masa “m” a un resorte, ésta lo estira a una longitud “s” y llega a una posición de
equilibrio, en la que su peso W, está equilibrado por la fuerza de restauración ks.

También debemos recordar que el peso se define por W = mg, donde la masa se expresa en slugs,
kilogramos o gramos y g = 32 ft/ s 2 , 9.8 m/ s 2 , 980 cm/ s 2 , respectivamente como se aprecia en
la figura (b) la condición de equilibrio es mg = ks ó mg–ks= 0. Si la masa se desplaza una distancia x
respecto de su posición de equilibrio, la fuerza de restitución del resorte es k(x + s).

Suponiendo que no hay fuerzas de retardo que actúen sobre el sistema y que la masa se mueve libre de otras
fuerzas externas (movimiento libre), entonces podemos igualar la segunda ley de Newton con la fuerza neta,
o resultante, de la fuerza de restitución y el peso:

d2 x
m =−k ( s+ x ) +mg=−kx +mg−ks … (1)
dt 2

d2 x
m =−kx …(1)
dt 2

El signo negativo de la ecuación (1) indica que la fuerza de restitución del resorte actúa en la dirección
opuesta del movimiento. Además, podemos adoptar la convención que los desplazamientos medidos abajo de
la posición de equilibrio son positivos como observamos.

Escuela de Ingeniería Civil Páá giná 46

 Ecuación diferencial del movimiento libre no amortiguado, Si dividimos la ecuación (1) por la masa “m”,
d2 x k
obtendremos la ecuación diferencial de segundo orden ( )
+
d t2 m
x=0 , 0 sea

2
d x
m 2
=−kx …(1)
dt
Obtendremos:

d2 x k
+( )
d t2 m
x=0

d2 x
2
+ ω2 x=0 … .(2)
dt

k
Donde ω2 = . Se dice que la ecuación (2) describe el movimiento simple o movimiento libre
m
no amortiguado, dos condiciones iniciales obvias relacionadas con la ecuación (2) son: x ( 0 )=α
y x ' ( 0)=β , el desplazamiento inicial y la velocidad inicial de la masa, respectivamente.

Solución y ecuación del movimiento, Para resolver la ecuación (2) observemos que las soluciones
de la ecuación auxiliar m2+ ω2=0 son los números complejos m1=ω i , m 2=−ω i
Así la solución general de (2) es:

x ( t )=C 1 cosωt + C2 senωt =0 …(3)

Escuela de Ingeniería Civil Páá giná 47

1.2 Sistemas de resorte y masa:

 Movimiento Amortiguado Libre:

El concepto del movimiento armónico libre no es realista porque el movimiento que describe la
ecuación (1) supone que no hay fuerzas de retardo que actúan sobre la masa en movimiento.

A menos que la masa esté colgada en un vacío perfecto, cuando menos habrá una fuerza de
resistencia debida al medio que rodea al objeto. Como observamos en la siguiente Imágenes, la
masa podría estar suspendida en un medio viscoso o conectado a un dispositivo amortiguador.

Ecuación diferencial del movimiento amortiguado libre: en mecánica se considera que las fuerzas de
amortiguamiento que actúan sobre un cuerpo son proporcionales a alguna potencia de la velocidad
instantánea. En particular, supondremos en el resto de la descripción que esta fuerza está expresada
por un múltiplo constante de dx/dt.

Cuando no hay otras fuerzas externas aplicadas al sistema, se sigue por la segunda ley de Newton:

d2 x dx
m 2
=−kx− β … .(1)
dt dt

Donde β es una constante de amortiguamiento positiva y el signo negativo es consecuencia del

hecho de que la fuerza amortiguadora que actúa en dirección opuesta a la del movimiento.

Al dividir la ecuación (1) por la masa “m”, la ecuación diferencial del movimiento amortiguado
libre es:

Escuela de Ingeniería Civil Páá giná 48

 d 2 x β dx k
+ ( ) ( )
+
dt 2 m dt m
x=0, es decir:

Donde:

2 λ= ( mβ ) ,ω = mk
2

d2 x dx
2
+2 λ +ω 2 x=0 … .(2)
dt dt

El símbolo 2 λ solo se usa por comodidad algebraica, porque así la ecuación auxiliar queda
2 2
m + 2 λ m+ w =0 y las raíces correspondientes son:

m1=−λ+ √ λ2 −ω2 m2=−λ− √ λ 2−ω 2

Ahora podemos distinguir tres casos posibles que dependen del signo algebraico de λ2−ω2 .
Puesto que cada solución contiene al factor de amortiguamiento e−λ t , λ> 0 , los
desplazamientos de la masa se vuelven insignificantes cuando el tiempo es grande.

CASO I: λ2−w2 >0 .

Aquí, se dice que el sistema está sobre amortiguado porque el coeficiente de amortiguamiento,
β , es grande comparado con la constante de resorte, k.

m1 t m2 t
La solución correspondiente de (2) es X ( t )=c 1 e + c2 e , o bien

−λ t 2 2 2 2

X ( t )=e (c 1 e √ λ −ω t + c 2 e−√ λ −ω t ) … .3

Esta ecuación representa un movimiento

suave y no oscilatorio. La figura muestra dos
gráficas posibles de x(t)

CASO II: λ2−w2=0

Se dice que el sistema está críticamente

amortiguado puesto que cualquier pequeña
disminución de la fuerza de
amortiguamiento originaría un movimiento oscilatorio. La solución general de la ecuación (ll) es

Escuela de Ingeniería Civil Páá giná 49

 x ( t )=c 1 e m t +c 2 t em t , es decir,
1 1

c 1+ c 2 t
x(t) )….(4)
¿ e−λ t ¿

CASO III: λ2−w2 <0

Se dice que el sistema está sub amortiguado porque el coeficiente de amortiguamiento es pequeño
en comparación con la constante del resorte. Ahora las raíces m1 y m2 son complejas:

m 1=−λ+ √ w 2−λ2 i , m 2=−λ− √ w2−λ 2 i

Entonces, la solución general de la ecuación (ll) es

x (t)=e−λt ( c 1 cos √ w 2−λ 2 t+ c 2 tsen √ w2− λ2 t ) … .(5)

Como se aprecia en la figura, el movimiento que

describe (9) es oscilatorio pero, a causa del coeficiente
e−λ t , las amplitudes de vibración tienden a cero
cuando t → ∞

Problema (01)

Una masa de que pesa 2 lb, hace que un resorte se estire 6 in, Cuando t=0, la masa se suelta desde un punto a

4 ft
8 in abajo de la posición de equilibrio con una velocidad inicial, hacia arriba de . Deduzca la
3 s

ecuación del movimiento Libre y Determinar la posición de la masa después de 2 segundos de soltarla masa

Solución

Escuela de Ingeniería Civil Páá giná 50

 6 in

2lb
8 in

2lb

Como se trata de un movimiento libre amortiguado, obedece a la siguiente ecuación diferencial:

2
d x k
+ x=0
d t2 m

Las medidas expresadas en pulgadas debemos convertirlas ó pasarlas a unidades pie, en el ejercicio

convertimos las unidades de in a” ft “.

12∈¿
1 ft 1
= ft
¿ 2
6∈¿ 6∈¿

Escuela de Ingeniería Civil Páá giná 51

 12∈¿
1 ft 2
= ft
¿ 3
8∈¿ 8∈¿

Otro Dato del Problema donde nos plantea, que tiene una velocidad inicial, el cual este se suelta hacia arriba

por lo tanto ira (-):

−4
x´ ( 0) = ft
3

La posición en el eje 0:

x ( 0 )=8 ft

Debemos convertir las unidades de peso, que están en libras, en unidades de masa partimos de:

W =m. g

2lb =m(32 ft /s 2)

1
m= slug
16

Aplicaremos También la Ley Hooke:

F=k (x )

1
2=k ( )
2

lb
k =4
ft

Esto implica que la constante del resorte es K=4 lb/ft ….Por lo tanto esta ecuación se transforma en

(aplicando la segunda ley de Newton):

Escuela de Ingeniería Civil Páá giná 52

 ∂2 x
m =−kx
d t2

1 ∂2 x
( ) 2 =−4 x
16 d t

2
∂ x
2
=−64 x
dt
2
∂ x
2
+ 64 x=0
dt

Como observamos ya tenemos armada la ecuación diferencial, damos paso a solucionar la ecuación

diferencial de la siguiente manera:

∂2 x
2
+ 64 x=0
dt

r +64=0

dx
Si decimos que r= , tenemos:
dt
2
r +64=0
2
r =−64

r=± √ −64

r=± 8 i

EL SISTEMA FUNDAMENTAL DE SOLUCIONES:

S . F . S={ e 0 t cos 8 t ; e 0 t sen 8 t }

LA SOLUCION GENERAL DE LA ECUACION DIFERENCIAL OBTENIDA:

x ( t )=C 1 cos 8 t+C 2 sen 8 t

Además: x ´ ( t )=−8 C1 sen 8t +8 C 2 cos 8 t

Escuela de Ingeniería Civil Páá giná 53

 * Pero por dato inicial del problema tenemos que:

2
x ( t )=C 1 cos 8 t+C 2 sen 8 t , pero x ( 0 )= ft
3

2
=C1 cos ⁡( 0)+C2 sen (0)
3

2
=C1
3

* Además:

4 −4 −4
x ´ ( t )=−8 C1 sen 8t +8 C 2 cos 8 t , pero x ´ (0)= ft x´ (0)= ft x´ ( 0) = ft
3 3 3

−4
=−8 C1 sen (0)+8 C2 cos ⁡( 0)
3

−4
=8 C 2
3

−1
=C 2
6

Y la ecuación diferencial queda de la siguiente manera:

2 1
x ( t )= cos 8 t− sen 8t
3 6

Finalmente la posición de la masa después de 2 segundos de soltar la masa es:

2 1
x ( 2 ) = cos 16− sen 16
3 6

x ( 2 ) =0.6 ft

REPUESTA: La posición de la masa después de 2 segundos de soltarla es de 0.6 ft

Problema (02)

Escuela de Ingeniería Civil Páá giná 54

 Una fuerza de 400 N estira un resorte 2 cm. Una masa de 50 kg se sujeta al extremo del resorte y se la suelta

desde la posición de equilibrio con una velocidad dirigida hacia arriba de 10 m/s. Halle la ecuación del

movimiento.

Solución:

Aplicando la ley de Hooke, 400=k (2) implican que la constante del resorte es k =200 , por tanto, la

ecuación:

F=400 Nx=2 mm=50 kg v=10 m/s

F=k (x )

400=k (2)

400
=k
2
k =200

Ecuación del movimiento

2
∂ x
m =−kx
d t2
2
∂ x
(50) =−200 x
d t2

r +4=

r
(¿¿ 2+4 )=0
¿
2
r =−4
r=± 2 i

2
ω =4
ω=2
dx
|
❑ =−10
dt t=0

Escuela de Ingeniería Civil Páá giná 55

 EL SISTEMA FUNDAMENTAL DE SOLUCIONES:

S . F . S={ e 0 t cos 2 t ; e 0 t sen 2 t }

LA SOLUCION GENERAL DE LA ECUACION DIFERENCIAL OBTENIDA:

x ( t )=C 1 cos 2t +C2 sen 2 t

Pero w= x(t)=2
x ( t )=c 1 cos 2 t +c 2 sen 2 t

2=c1 cos 2(0)+c 2 sen 2(0)

2=c1 ( 1 )+ 0

c 1=2

x ´ ( 0 )=[−10=−2C 1 sen 2 t+2 c 2 cos 2t ]

[−10=−2(2)sen 2 ( 0 ) +2 c2 cos 2 ( 0 ) ]

[−10=0+2 c 2 (1)]

x ´ ( 0 )=[−10=2 c 2 ]

c 2=−5
Reemplazamos:

x ( t )=c 1 cos 2 t +c 2 sen 2 t

x ( t )=2 cos 2 t−5 sen 2 t

Escuela de Ingeniería Civil Páá giná 56