INSTITUTO TECNOLÓGICO Y DE ESTUDIOS SUPERIORES DE MONTERREY

Escuela de Ingeniería y Ciencias

Campus Ciudad de México

PRIMER PROYECTO PARCIAL

​“Modelando la resistencia del aire”
Materia:
Ecuaciones Diferenciales

Semestre:
Febrero - Junio 2020

Profesor:
Maurice Oliva Leyva
Daniel Eduardo Malpica Cuevas A01656812
Marco Ortega Maldonado A01656512
Enrique Vázquez Bobadilla A01650630
Jesús Martín Torres Ramírez A01339712
1. ¿Qué tan alto? (Sin resistencia del aire) (20 puntos)
Suponga que una pequeña bala de cañón de masa m = 5 kg se dispara verticalmente hacia
arriba, con una velocidad inicial de Vo = 90 m/s. La respuesta a la pregunta “¿Qué tanto sube
la bala de cañón?”, depende de si se considera la resistencia del aire.

(a) Suponga que se desprecia la resistencia del aire. Si la dirección es positiva hacia arriba, la
altura ​s(t) de la bala respecto al suelo al tiempo ​t est´a dada por la ED ds2 / dt2 = −g, donde
g = 9.8 m/s 2 . Pero ds/dt = v(t), entonces de la ED anterior se tiene que dv/dt = −g.
Encuentre la velocidad ​v(t)​ de la bala al tiempo​ t​.

SOLUCIÓN:

ds 2
dt 2
=− g
ds
dt = v (t)
dv
dt =− g
ds 2 dv
dt 2
= dt =− g = a
m
g=-9.8 s
m
Vo=90 s

Dada las fórmulas de física que describen el movimiento de un cuerpo en tiro vertical hacia
arriba, despreciando obviamente la resistencia del aire tenemos:

Vy = -g*t + Vo

donde Vy es la velocidad vertical, g es la gravedad, t es el tiempo y Vo es la velocidad inicial

con la que se avienta el objeto, en este caso con la que de dispara la bala.

De esta manera, se obtiene que:

V(t)=-9.8t+90
(b) Utilice el resultado que se obtuvo en el inciso a) para determinar la altura s(t) de la bala
medida desde el nivel del suelo, y grafíquela. Determine la altura máxima que alcanza la bala.

SOLUCIÓN:

Para determinar la altura máxima, es necesario conocer el tiempo de vuelo de la bala lanzada
desde el suelo hasta que llega a su punto máximo, donde la velocidad se vuelve cero, por esa
razón igualamos la ecuación que obtuvimos en el inciso anterior a cero:
-9.8t+90=0
90
t= 9.8
t=9.1836 segundos

ds
Ahora, dado que dt = v (t) , entonces al integrar la ecuación se obtiene s(t):

ds
dt = v (t)
ds=v(t) dt

∫ ds = ∫ v(t) dt

s(t)= ∫− 9.8t + 90 dt

s(t) =− 4.9t2 + 90t +c

En t=0, la posición es también cero, puesto que la constante es 0. Esto resulta, finalmente en
la ecuación de posición:

s(t) =− 4.9t2 + 90t

Finalmente, sustituyendo el tiempo encontrado para el ascenso en la ecuación, obtenemos que

la altura máxima:

90 2 90
s(t) =− 4.9( 9.8 ) + 90( 9.8 )
Smax = 413.2653 metros.
La gráfica 1 muestra en azul, la trayectoria descrita por la bala desde su lanzamiento a nivel
de suelo hasta que vuelve a topar con él:

Gráfica 1. Desplazamiento de bala sin resistencia del aire.

2. ¿Qué tan alto? (Resistencia lineal del aire) (20 puntos)
Repita el problema 1, es decir, (a) encuentre la velocidad ​v(t) y​ (b) la altura ​s(t) de la bala,
pero esta vez suponga que la resistencia del aire es proporcional a la velocidad instantánea.
Esta es la razón por la que la altura máxima que alcanza la bala debe ser menor que la del
inciso b) del problema 1. Demuestre esto suponiendo que la constante de proporcionalidad es
k = 0.0025. Además, graque la función s(t) obtenida en este caso,
SOLUCIÓN:
g *m−k *v
m dv
dt
= − m
dv k *v
dt = −g− m
dv k
dt = m *v = −g
Utilizando Factor integrante

∫ mk dt k
u=e = emt

∫ d(ek/m *t * v)dt = ∫− g(ek/m *t )dt

k
k −g e m *t
e m *t v = k/m
+ C1
Usando Condición Inicial V(0)=V0
−g *m k
Vo= k
+ C 1 e− m (0)
g *m
V0+ k = C 1 = 19, 690
 m*g m*g − k *t
V (t) = − k
+ (V 0 + k
)ε m
m*g m*g
m ln( /V 0+ )
t= − k
k
k

5*9.8 9.8
−5ln( 0.0025 /90+5* 0.0025 )
t= 0.0025
= 9.1626 segundos
Integrar V(t) para obtener S(t)
m*g m − mk *t
s(t) = − k * t − C 1( k e ) + C2
5(9.8) 5 0.0025
0= − 0.0025
(0) − 19690( 0.0025 e− 5
(0)
) + C2
C 2 = 39, 380, 000
5(9.8) 5 0.0025
s(9.1626) = − 0.0025
(9.1626) − 19, 690( 0.0025 e− 5
(9.1626)
) + 39, 380, 000
s(9.1626) = 412 metros

Gráfica 2 (ecuación s(t)), obtenida con el software de Geogebra

3. ¿Qué tan alto? (Resistencia no lineal del aire) (20 puntos)

Análogamente al problema anterior, repita el problema 1, pero esta vez suponiendo que la
resistencia del aire es proporcional al cuadrado de la velocidad instantánea. Suponga que la
dirección positiva es hacia arriba y tome k = 0.0003.
SOLUCIÓN:
Tiempo para alcanzar altura máxima
dv k *v 2
dt
= (− g − m
) * dt

∫ g+dv = − ∫ dt
k*v 2
m
 √m
(tan−1 ( √vg**√km )) = −t+C
√*
g k

√ √
g *m kg
v (t) = k
tan[arctan( V o − m
t)]
√mg/k

t=
√ 5
9.8*0.0003
tan−1 (
√
90
5*9.8
0.0003
)

t = 9.0362 segundos

Altura máxima
m dv
dt = − m * g − k * v
2

∫ m*g+k
m*v
*v
= − ∫ dy
2

m
2k ln(m * g + k * v 2 ) = − y + C 1
−2*k*y
m * g + k * v 2 = C 2 e( m )

2*k*y
m*g
v 2 = e (− m )
− k

Usando t=0 ; v(0)=V0

m*g
V o2 = C 3 − k
m*g
C3 = V o + k
k *V o2 +m*g ( −2k*y ) m*g
V2 =( k
)e m − k

Usando, donde y es la altura en metros

m k *V o2 +m*g
y= 2*k
ln( m*g
)

y = 403.34 metros
Usando la distancia desde el punto máximo al suelo, con td como tiempo de descenso:
dv
m dy = − m * g + k * v2

∫ m*g−k
m*v
*v
= ∫ dy
2

m 2
− 2k ln(m * g − k v ) = y + h1
−2k*y
m * g − k * v 2 = h2 * e m
)

m*g −2k*y
v2 = k (1 − h 3 e( m
)
)
−2k*y
m*g
v2 = k
(1 −e m ) )
 −k*V 02 +mg
mln( mg
)
y= 2k

y = − 403.34 metros
El signo es negativo debido a que la bala se encuentra cayendo.

4. Todo lo que sube... (30 puntos)

(a) Es bien conocido que el modelo que desprecia la resistencia del aire, inciso (a) del
problema 1, predice que el tiempo que tarda la bala en alcanzar su altura máxima es el mismo
tiempo td que tarda la bala en llegar al suelo. Además, la magnitud de la velocidad de
impacto vi es igual a la velocidad inicial v0 de la bala. Compruebe ambos resultados.
(b) Utilizando el modelo lineal del problema 2 que considera la resistencia del aire, compare
el valor de ta con td y el valor de la magnitud de vi con v0. Aquí puede ser útil un programa
para determinar raíces de un sistema algebraico computacional (o una calculadora
graficadora).
(c) Por ´ultimo, compare el valor de ta con td y el valor de la magnitud de vi con v0, para el
caso del modelo no lineal del problema 3.
SOLUCIÓN:
(a)
Usando la ecuación del problema 1
V y = − gt + V 0
V (t) = − 9.8 t + 90
0 = − 9.8 t + 90
90
t= 9.8

t = 9.1836 segundos

s(t) = − 4.9 t2 + 90 t
s(t) = − 4.9 (9.1836)2 + 90(9.1836)
s(t) = 413.26 metros
Comprobación de la velocidad y tiempo, como la bala se dirige hacia abajo la gravedad
cambia de signo negativo a positivo y se presenta la siguiente ecuación:
v (t) = 9.8t + V 0
Como la bala se encuentra arriba se plantean las siguientes condiciones iniciales: V0=0,
V(t)=90 (debido a que V(i)= V0 del proceso de ascenso)
 90 = 9.8t + 0
90
t= 9.8

t = 9.1836 segundos
Con este resultado podemos confirmar que el t de ascenso es el mismo que el t de
descenso, de igual forma usaremos este t para obtener v(9.1836)
v (t) = 9.8(9.1836) + 0
v (9.1836) = 90 m/s = V 0 de ascenso
(b)

Gráfica 3 (ecuación s(t)), obtenida con el software de Geogebra

En esta gráfica podemos observar el comportamiento de la bala en un periodo de tiempo

0<=t<=18.339, esto nos quiere decir que el trayecto total de la bala hasta alcanzar la altura
máxima y después disminuir hasta llegar al suelo nuevamente es de 18.339 segundos. A
continuación realizaremos la comprobación, modelando la ecuación desde el punto máximo
hasta llegar al suelo.
dv
dt = − g + k v/m
dv
dt − mk v = − g

−∫ mk dt k
u=e = e− m *t
Multiplicando por el factor integrante e integrando obtenemos
k
v (t) = g mk + C 1e m *t
Usando V(0)=V0
 k
V 0 = g mk + C 1e m (0)
gm
V0= k + C1
gm
C1 = V 0 − k

C 1 = − 19600
gm gm k *t
V (t) = k
+ (V 0 − k
)e m
Despejando para t y usando la condición inicial V(t)=0
mg mg
mln(− /V 0− )
t= k
k
k

Sustituyendo para t con los siguientes valores, g=9.8, m=5, k=0.0025, V0=90, obtenemos
t descenso = 9.2048 segundos
Sumando el valor de t subiendo + t descenso, tenemos
t subiendo + t descenso = 9.1626 s + 9.2048 s = t total
t total = 18.3674 segundos
Este valor corresponde al obtenido en la gráfica 3
Comparación de las velocidades subiendo y en descenso
La siguiente gráfica (4) está descrita mediante la ecuación
0.0025
5*9.8 5*9.8
V (t) = − 0.0025
+ (90 + 0.0025
)e− 5 *t la cual representa a la bala cuando se encuentra
subiendo, por ello nos encontramos con que la velocidad parte de 90 m/s en el tiempo 0 (esto
de acuerdo a las condiciones iniciales impuestas por el problema), posterior a ello nos
encontramos con la siguiente intersección con el eje x la cual representa la altura máxima a la
que llega la bala al alcanzar​ t=9.1626 segundos.

Gráfica 4, representa la ecuación V(t) cuando la bala se encuentra subiendo

Por otra parte la siguiente gráfica (5) está descrita mediante la ecuación
0.0025
5*9.8 5*9.8
V (t) = 0.0025
+ (0 − 0.0025
)e 5 *t la cual representa a la bala cuando se encuentra bajando,
por ello nos encontramos con que la velocidad parte de 90 m/s en el tiempo 0 (esto de
acuerdo a las condiciones iniciales), posterior a ello nos encontramos con la siguiente
intersección con el eje x la cual representa la altura en la que la bala llega al suelo y alcanza
su altura inicial a la que llega la bala, este tiempo corresponde a ​t=9.2048 segundos.

Gráfica 5, representa la ecuación V(t) cuando la bala se encuentra bajando

c)
Como se utilizó en los incisos a y b, sabemos el tiempo necesario para llegar al punto más
alto con la primer ecuación fue de ​t = 9.1836 s.
Para obtener el tiempo que tardó en llegar al punto más alto con la ecuación del problema 3:
dv k *v 2
dt
= (− g − m
) * dt
Es necesario obtener esta fórmula en términos de la velocidad.

∫ g+dv = − ∫ dt
k*v 2
m

√m
(tan−1 ( √vg**√km )) = − t + C
√g*k

√ √
g *m kg
v (t) = k
tan[arctan( V o − m
t)]
√mg/k
Sabemos que en el punto más alto la velocidad es cero, por lo que al sustituir t, obtendremos
el tiempo que transcurrió para llegar al punto más alto.

t=
√ 5
9.8*0.0003
tan−1 (
√
90
5*9.8
0.0003
)

t = 9.0362 segundos
De esta manera nos damos cuenta que ciertamente tener la resistencia al cuadrado afectó la
altura máxima de la tercer ecuación. Podemos notar que aunque la diferencia solamente son
décimas de segundo, sí existe una diferencia.
Grafica 6. Tiempo hasta llegar a la altura máxima
Después deberemos obtener el tiempo en el que impactó con el suelo y así obtener la
velocidad final para poder comparar. El tiempo final en el que cayó parte de la ecuación de
velocidad 2 lo que nos da un tiempo muy similar al primero ya que la diferencia es mínima.

√ tanh −1 ( vo√k 2
mk
t = g )
√mg+kvo
td = 9.0389
td>ta
Al sumar los dos obtenemos un tiempo total de 18.075 s
Con la segunda fórmula de velocidad sustituimos el valor de td.
De esta manera obtenemos que vf = 87.1894
Por lo que podemos concluir que en con esta fórmula Vi no será igual a Vo.

5. Las gotas de lluvia siguen cayendo... (10 puntos)

La velocidad promedio de una gota de lluvia cayendo es de 11 km/h. En una búsqueda rápida
por la internet se encontró que el meteorólogo Jeff Haby ofreció información adicional de que
una gota de lluvia esférica en “promedio” tiene un radio de 0.1 cm y un volumen aproximado
de 0.0000000044 m3. Utilice estos datos, y si se necesita investigue más, y haga otras
suposiciones razonables para determinar si “la velocidad promedio de . . . 11 kilómetros por
hora es consistente con los modelos de resistencia del aire de los problemas anteriores.
SOLUCIÓN:
Usando el primer modelo
 m*g m*g − k *t
v (t) = k
+ (V 0 − k
)ε m
Obteniendo masa de la gota
m = V olumen * densidad
El volúmen es proporcionado en el problema y la densidad es de 997 kg/m^3, por ello:
m = 4.4x10−9 m3 * 997kg/m3
m = 4.38x10−6 kg
Para obtener k, usamos la siguiente fórmula:
k = 12 P * A * S
Donde P=1
A= área
A = 4Π(1x10−3 m)2 = 0.012566m2
Donde S=0.4
Por ello:
k = 2.5132x10−3
Encontramos que el tiempo estimado en que tarda una gota de agua desde la nube al
suelo es de 270 segundos, por ello sustituimos los valores antes mencionados en la
primera ecuación.
Usando V0=11km/hr=3.0556 m/s
2.5132x10−3
4.3868x10−6 *9.8 4.3868x10−6 *9.8 − 4.3868x10−6 (270)
v (270) = 2.5132x10−3
+ (3.0556 − 2.5132x10−3
)e
v (270) = 0.0171m/s = 0.061 km/hr

Usando el segundo modelo

√
g *k
√
m
tan−1 (V 0 k
v (t) =
√ g *k
tan(− m*g
− m * t))

Usando las constantes mencionadas arriba obtenemos que:

v (270) =
√ 4.3868x10−6
9.8*2.5132x10−3
tan(− tan(3.0556
√ 2.5132x10−3
4.3868x10−6 *9.8
−
√ 9.8*2.5132x10−3
4.3868 * 270))
v (270) = 16.11m/s = 57.9 km/hr
Comparamos nuestros resultados con los obtenidos en 1904, el físico Philipp Lenard, él
construyó un túnel de viento vertical para calcular la velocidad de las gotas al caer y encontró
que la velocidad de las gotas se encuentra entre (8-32) km/hr, por lo que podemos concluir
que el segundo modelo es el más acercado pero se tiene que mejorar para dar un resultado
más real.