Método de Separación de Variables Aplicado a un Modelo de Caída Libre Retardada por la Resistencia del Aire (Enfoque Analítico Matemático y Numérico Computacional)

in stem-espanol •  6 years ago  (edited)

Saludos estimados seguidores de Steemit y comunidad de #stem-espanol en general, continuando con las aplicaciones de las matemáticas en la ciencia, el día de hoy abordaré uno de los modelos utilizados para describir el movimiento de caída libre retardado por la resistencia del aire, uno de mis artículos anteriores se enfocó en un problema similar referente a la caída de un paracaidista, sin embargo, el presente artículo aborda una situación diferente.

esquema.png
Movimiento de una bala en ascenso. Fuente: Elaboración propia.

Las matemáticas suelen ser utilizadas para estudiar diversos fenómenos del mundo real, por ejemplo, cuando se dispara cualquier tipo de proyectil, la leyes de la mecánica pueden utilizarse para describir este movimiento, muchas veces al aplicar dichas leyes aparecen ecuaciones diferenciales ordinarias que deben ser resueltas para obtener una solución que pueda interpretarse y ser interpretada en el mundo real.

De esta manera, las abstracciones matemáticas permiten modelar situaciones reales y describir su evolución en el tiempo, sin embargo, no basta solo con modelar los fenómenos, también se deben aplicar técnicas matemáticas para resolver los modelos obtenidos, en el caso del disparo de una bala hacia arriba el movimiento puede ser descrito teniendo en cuenta que se trata de un típico caso de caída libre retardada por la resistencia del aire, en este sentido, la física ha construido una serie de modelos para describir este movimiento, uno de los cuales se explica en el presente artículo.

El movimiento de una bala que es disparada hacia arriba puede ser modelado utilizando algunos postulados básicos de la mecánica teniendo en cuenta que se trata de un ejemplo de caída libre en el que interviene la resistencia del aire, diversas expresiones matemáticas pueden ser aplicadas para modelar este fenómeno.

Por ejemplo, la segunda ley de newton establece que la fuerza que actúa sobre la bala es siempre igual a

img1.png

Donde (m) es la masa de la bala y (a) su aceleración, definiendo la dirección positiva hacia arriba se establece que una fuerza positiva actúa acelerando la bala hacia arriba y una fuerza negativa actúa generando un movimiento hacia abajo.

En el caso del movimiento de la bala existen dos fuerzas que componen la fuerza total que actúa sobre la bala, la primera es la fuerza de gravedad y la segunda la resistencia del aire, matemáticamente esto se expresa como:

img2.png

Donde Ft es la fuerza total, Fg la fuerza debida a la gravedad y Fr la resistencia del aire. La fuerza generada por la gravedad puede ser modelada teniendo en cuenta que dicha fuerza arrastra a la bala hacia abajo es decir en dirección negativa por lo tanto es igual a

img3.png

En el caso de la resistencia del aire se utilizan diversas expresiones para modelarla, en un artículo anterior se modelo utilizando la expresión

img4.png

Donde v es la velocidad del objeto y k una constante de proporcionalidad, sin embargo, esta expresión no es la única utilizada para modelar la resistencia del aire, en el caso del disparo de proyectiles se suele utilizar otra fórmula para la resistencia del aire la cual se expresa como:

img5.png

Se observa que la resistencia del aire es negativa debido a que actúa hacia abajo al detener a la bala en su ascenso, combinando las ecuaciones 1,2 y 3 se obtiene

img6.png

Teniendo en cuenta la segunda ley de newton

img7.png

Recordando que la aceleración de un objeto en movimiento es equivalente a la razón de cambio de su velocidad en el tiempo (derivada), se puede sustituir la aceleración por:

img8.png

Obteniendo

img9.png

Ahora bien, esta expresión analiza el ascenso de la bala desde la perspectiva de la dinámica de las fuerzas involucradas, este análisis omite el cálculo de la distancia recorrida por la bala, si se quiere analizar la distancia se puede partir del hecho de que la velocidad de la bala equivale a la razón de cambio (derivada) de la distancia en el tiempo, es decir:

img10.png

Según la regla de la cadena

cad.png

Luego (5) puede reescribirse como

img13.png

Esta Ecuación Diferencial Ordinaria involucra una derivada que incluye a la distancia, el método utilizado para resolverla se conoce como separación de variables y consiste en despejar en cada lado de la ecuación cada variable (v con dv) y (x con dx) y luego integrar ambo lados.

img14.png

Dividiendo entre m

img15.png

Integrando ambas expresiones se observa que

img16.png

Realizando la sustitución

img17.png

y

img18.png

Se obtiene

img19.png

Cuando la bala es disparada se encuentra a cero metros de altura de su punto de partida, es decir, cuando v=v0 se cumple que x=0 por lo tanto

img20.png

Al sustituir la constante de integración en la solución general se obtiene

img21.png

Ejemplo

Una bala es disparada hacia arriba con una velocidad inicial vertical de 988 m/s, si la bala posee una masa de 0,00356 kg y una constante de resistencia del aire de k=0,0000073 kg/m determine la altura máxima que alcanza la bala.

Para encontrar la altura x se necesita conocer la velocidad de la bala en ese instante de tiempo debido a la naturaleza del movimiento en ese instante la velocidad de la bala es de 0 m/s ya que en ese punto comenzará a caer, (se omiten las unidades para simplificar los cálculos).

img22.png

Es decir, la bala alcanza una altura máxima de 1.298 metros antes de empezar su descenso, adicionalmente también se puede analizar la variación de la velocidad de la bala en el tiempo, para este fin se utilizará la ecuación número 5 obtenida con anterioridad

img23.png

Dividiendo ambos lados de la expresión entre m

img24.png

La solución de esta ecuación diferencial permitirá obtener el valor de la velocidad en cada instante de tiempo. En el problema planteado la masa de la bala es de 0,00356 kg y la constante de resistencia del aire de k=0,0000073 kg/m, por lo tanto, la ecuación diferencial se expresa como

img25.png

Aplicando factor común y separación de variables

img26.png

Integrando ambos miembros

img27.png

La integral del lado izquierdo puede ser resuelta utilizando la siguiente fórmula

img28.png

Para realizar la sustitución se debe observar que

img29.png

Luego al sustituir la ecuación número 11 en la número 10 se obtiene

img30.png

Reemplazando los valores de u y a

img31.png

Debido a que según los datos del problema en el instante de tiempo 0 la velocidad inicial de la bala es de 988 m/s se cumple que

img32.png

Al sustituir la constante de integración en la solución general de la ecuación diferencial ordinaria se obtiene

img33.png

img34.png
Fórmula para calcular la velocidad de la bala en función del tiempo de recorrido.

Esta función permite calcular la velocidad de la bala en cualquier instante de tiempo, si se desea determinar el tiempo que tarda la bala en ascender hasta el punto más elevado de su trayectoria se debe recordar que en dicho punto su velocidad es cero debido a que está a punto de empezar a caer , es decir:

img35.png

Por lo tanto, la bala tarda 10,58799025 segundos en llegar hasta su punto más elevado que como se demostró se encuentra ubicado a 1298,234140 metros de altura.


Solución Numérica (Método de Runge-Kutta de 4° orden)


En el ejemplo mostrado en el presente artículo se pudo resolver la ecuación diferencial ordinaria aplicando métodos provenientes del análisis matemático, sin embargo, existen casos en que dichos métodos no son suficientes, ante esta realidad se han desarrollado una serie de métodos numéricos que permiten obtener aproximaciones bastante fiables a las soluciones analíticas de las ecuaciones diferenciales ordinarias.

Con el fin de garantizar la precisión en la solución se suele utilizar un método desarrollado en los años 1900 por los matemáticos Carl David Tolmé Runge y Martin Wilhelm Kutta el cual se conoce como el método de Runge-Kutta de 4° orden, este método es uno de los más utilizados en la actualidad debido a su alto grado de precisión y por ser programado de forma elemental, para obtener la solución numérica de la ecuación diferencial ordinaria mediante el método de Runge-Kutta de 4° orden solo se requiere una condición inicial.

Las formulas utilizada por este método son las siguientes:

img36.png

Donde h es el tamaño de paso, un valor bastante pequeño.

En líneas generales se inicia con la condición inicial y(x0)=y0, se estiman las cuatro pendientes s1, s2, s3 y s4 y mediante el promedio ponderado de estas pendientes se aproxima el siguiente valor de yi a partir del valor anterior.

Este método puede ser aplicado al problema mostrado en el presente artículo, recordando que la ecuación número 9 y la velocidad inicial de 988 metros por segundo permiten construir el problema de valor inicial:

img37.png

Debido a que la bala asciende hasta que t=10,58799025 segundos el intervalo llega hasta dicho punto es decir podemos aplicar el método de Runge-Kutta de 4° orden hasta un valor de t=10,5; si se selecciona un tamaño de paso de 0,25 se puede aplicar el método de forma eficiente.

img38.png

En el presente artículo se utilizó el Sistema de Algebra Computacional WxMaxima el cual permite aplicar este tipo de métodos numéricos, a continuación se muestra el Script utilizado


/*************************************/
/* MÉTODO DE RUNGE-KUTTA de 4° ORDEN /
/
Licdo Ysmael González /
/
U.C. Cálculo Numérico /
/
UNEFM 2016 /
/*************************************/
f(t,v):=-9.8-0.002050562
v^2$ /Ecuación diferencial a resolver/
t0:0$ v0:988$ /Condición inicial t0 y v0/
tf:10.5$ delta:0.25$ /Valor final del intervalo tn y tamaño de paso/
fpprintprec:5$ /Número de cifras decimales/
sol(t):=69.13159633tan((10.58799025-t)/7.054244523)$ /Solución analítica de la EDO/
tn:t0$ n:0$ vn:v0$ iter:(tf-t0)/delta$
print("");
while n<=iter do
(
s1:f(tn,vn),
s2:f(tn+(delta/2),vn+(s1
delta/2)),
s3:f(tn+(delta/2),vn+(s2delta/2)),
s4:f(tn+delta,vn+(s3
delta)),
g:sol(tn),
ep:abs(g-vn)100/g,
print("It=",n," t=",float(tn)," V(t)=",float(vn)," V.E.=",float(g)," E.P.=",float(ep),"%"),
n:n+1,tn:tn+delta,vn:vn+(delta/6)
(s1+2s2+2s3+s4)
);

Este programa genera la siguiente salida:

Tiempov(t)Sol. Aprox.Error
0.0988.0988.00.000%
0.25654.1654.110.001654%
0.5487.98487.960.003%
0.75388.35388.330.003938%
1.0321.81321.790.003897%
1.25274.12274.110.003628%
1.5238.18238.170.003332%
1.75210.07210.060.003064%
2.0187.43187.420.002833%
2.25168.75168.750.002639%
2.5153.06153.050.002477%
2.75139.64139.640.002343%
3.0128.02128.020.00223%
3.25117.84117.830.002138%
3.5108.81108.810.002062%
3.75100.74100.740.002%
4.093.45893.4560.001952%
4.2586.84486.8430.001915%
4.580.79580.7930.001888%
4.7575.22775.2250.001872%
5.070.07270.0710.001865%
5.2565.27565.2740.001867%
5.560.78960.7880.001879%
5.7556.57556.5730.001901%
6.052.59752.5960.001933%
6.2548.82948.8280.001976%
6.545.24545.2440.002032%
6.7541.82441.8230.002102%
7.038.54638.5450.002189%
7.2535.39535.3940.002295%
7.532.35732.3560.002424%
7.7529.41829.4170.002583%
8.026.56626.5650.002779%
8.2523.79123.790.003025%
8.521.08221.0820.003336%
8.7518.43218.4310.003739%
9.015.83115.8310.004277%
9.2513.27313.2720.005026%
9.510.74810.7480.006131%
9.758.25188.25120.007909%
10.05.77635.77570.01122%
10.253.31553.31480.01946%
10.50.8630.86230.07464%

Tabla N° 1 Variación de la velocidad en el tiempo según la solución Analítica y Numérica. Fuente. Elaboración propia.

En la tabla n° 1 se observa como disminuye la velocidad en el tiempo, se evidencia que tanto la solución analítica obtenida al resolver la ecuación diferencial ordinaria como la solución numérica al aplicar el método de Runge-Kutta de 4° orden arrojan valores muy cercanos. Recordando que la función de velocidad es la siguiente:

38b.png

Se puede representar gráficamente como

d1.JPG
Gráfica N° 1 Velocidad de la bala(m/s) en función del tiempo(s) durante el ascenso. Fuente: Elaboración propia mediante la herramienta https://www.desmos.com/

La altura de la bala se puede calcular teniendo en cuenta la ecuación número 8

img39.png

Al sustituir los datos del problema planteado

img40.png

Al sustituir la ecuación 12 en la expresión anterior se obtiene

img41.png

Gráficamente la altura en función del tiempo se observa en la siguiente imagen

d2.JPG
Gráfica N° 2 Altura de la bala(m) en función del tiempo(s) durante el ascenso. Fuente: Elaboración propia mediante la herramienta https://www.desmos.com/

Conclusiones


  1. Existen diversos modelos que estudian la resistencia del aire, dependiendo del tipo de problema estudiado suele aplicarse cada uno de ellos.
  2. Las técnicas analíticas de resolución de Ecuaciones Diferenciales Ordinarias permiten obtener soluciones exactas pero no siempre pueden ser aplicadas.
  3. El Método de Runge-Kutta de 4° orden tiene una gran variedad de aplicaciones en la ciencia en general permitiendo obtener respuestas muy precisas a modelos matemáticos.

REFERENCIAS BIBLIOGRÁFICAS Y LECTURAS RECOMENDADAS


  1. Anton, Bivens y Davis (2010), Cálculo de una variable, Trascendentes tempranas. 2da edición Editorial Limusa Willey.

  2. González (2018) Modelado Matemático mediante Ecuaciones Diferenciales Ordinarias, Caso de Estudio: Caída Libre retardada por la Resistencia del Aire (Resolución mediante SCILAB) Disponible aquí.

  3. Resnick, Halliday y Krane (1993), Física. 3ra edición Compañía Editorial Continental México Volumen 1.

  4. Saenz (2009), Cálculo Integral con Funciones Trascendentes Tempranas para Ciencias e Ingeniería. 2da edición Editorial Hipotenusa.

  5. Sauer (2013), Análisis Numérico. 2da edición Editorial Pearson.


Imagen de dominio público. Fuente: @steemstem

Authors get paid when people like you upvote their post.
If you enjoyed what you read here, create your account today and start earning FREE STEEM!
Sort Order:  

una gran referencia matemática

Gracias por el apoyo @karimuddinn

Felicitaciones! @ydavgonzalez. Gran talento en esta comunidad.

Gracias @elvigia por el apoyo, me alegra que revises mi perfil con frecuencia, Saludos.

Muy buen articulo. Me hubiera servido mucho en mis tiempos universitarios.

Así es, gracias por el apoyo @marcusantoniu26.

Excelente trabajo @ydavgonzalez, enalteciendo como siempre la aplicación de las matemáticas a problemas reales, contextualizado en este caso al problema de caída libre retardada por la resistencia del aire; en definitiva un problema de mucho de interés, muy bien abordado, con ese toque didáctico. Disfruté mucho de su lectura, saludos y muchos éxitos!!!

Gracias @eliaschess333 he visto que en tus artículos también usas mucha matemática, lo cual me agrada bastante, te agradezco que leas mis publicaciones.

Saludos.

Interesante método.

Felicitaciones por tu artículo, como siempre abordas los temas con gran cuidado y con bases sólidas que justifican tus procedimientos. Éxitos!

Gracias por tu comentario y por tomar algo de tu tiempo para leer y analizar mi trabajo, así es, en la ciencia las matemáticas permiten abordar diversos fenómenos físicos.

Saludos @isgledysduarte.

Espero seguir contando con tu apoyo en mis futuros artículos.

Excelente, la aplicación que le encuentras a la matemática en un problema de física, me gusto mucho tu artículo @ydavgonzalez.

Saludos @veronicacoli, que bueno que hayas leído mi post, gracias por apoyarme.

Espero seguir contando con tu apoyo en mis futuros artículos.

Te felicito @ydavgonzalez muy interesante la aplicación de este modelado matemático a través de un software. Un trabajo impecable y que nos muestra la importancia de las matemáticas para la descripción y abordaje de temas en las distintas ramas de la ciencia.

¡Felicitaciones!

Gracias por el apoyo, como químico también has tenido oportunidad de utilizar la matemática al estudiar diversos fenómenos de la ciencia.

Saludos @joseleogon.

Espero seguir contando con tu apoyo en mis futuros artículos.

Bonito problema de mecánica clásica, uno de los libros que posee esa teoría es el Marion (Dinámica Clásica). Buen trabajo @ydavgonzalez

Gracias por el apoyo, cuando tenga algo de tiempo lo revisaré.

Saludos @germanmontero.

Espero seguir contando con tu apoyo en mis futuros artículos.

Impecable como siempre!, es bastante practico mostrar el resultado numerico de este tipo de problemas.

Así es, los Sistemas de Álgebra Computacional permiten desarrollar estos métodos mediante el uso del ordenador y obtener una solución numérica muy aproximada a la solución real. De esta forma se combinan la matemática y la programación en pro de obtener resultados a problemas que describen fenómenos físicos.

Saludos @vjap55.

Espero seguir contando con tu apoyo en mis futuros artículos.

Muy buena explicación donde aplicas tus conocimientos matemáticos, donde se emplea la fuerza y existe un movimiento caracterizado por un desplazamiento, así que también podemos ver la intervención de un tipo de energía en todo este proceso. Me fascinan los fenómenos físicos!

Gracias, claro también se podría realizar un análisis desde la perspectiva de la energía involucrada, saludos @azulear.

Excelente post @ydavgonzalez. Lo que más me agrada es que usas software libre :D

Así es, WxMaxima es un programa interesante. que permite abordar este tipo de problemas.

Saludos @eniolw.

Espero seguir contando con tu apoyo en mis futuros artículos.

Muy buen artículo @ydavgonzalez, excelente trabajo como siempre. Saludos!

Gracias por el apoyo estimada amiga, es bueno saber que hay personas que se interesan por el trabajo que con esmero realizamos.

Saludos @monserratt3.

Espero seguir contando con tu apoyo en mis futuros artículos.

Interesante aplicación de las matemáticas en un problema de física, buen artículo sobre todo la parte en la que muestras el resultado numérico obtenido mediante un programa. Éxito :).

Gracias por el apoyo estimada @cataplum, las herramientas de Software Libre son muy útiles en este tipo de problemas.

Saludos.

Si, el uso del programa te permite darle otro enfoque a la resolución del problema.

Así es, estimada @cataplum desde la perspectiva analítica y numérica.

Saludos :)

Excelente, saludos.

Que inteligente te admiro. Es difícil comprender del todo... Tu mente es brillante..

Hi @ydavgonzalez!

Your post was upvoted by utopian.io in cooperation with steemstem - supporting knowledge, innovation and technological advancement on the Steem Blockchain.

Contribute to Open Source with utopian.io

Learn how to contribute on our website and join the new open source economy.

Want to chat? Join the Utopian Community on Discord https://discord.gg/h52nFrV

Hola @ydavgonzalez, me gusta mucho que uses el computador para mostrar este tipo de problemas, sobre todo su resolución numérica.

Con mucho respeto a tu trabajo, te dejo algunas sugerencias,

  1. La forma de resuelves el problema de variables separables, no es la mejor forma de presentarlo, ya que el símbolo de la derivada (dy/dx) no puedes usarlo como si fuese en realidad una fracción que cambia de posición, ya que como símbolo, es eso, un solo símbolo que denota una operación, la manera correcta de hacerlo, debería ser mostrando que lo puedes intercambiar usando cambios de variables adecuados. Cuando se integra una igualdad, se integra con respecto a una misma variable ambos lados, no puedes integrar una parte de la igualdad con respecto a una variable y la otra con respecto a otra variable distinta. Te dejo un link a una publicación que realicé algunos días atrás, precisamente abordando este punto (Introducción a las Ecuaciones Diferenciales Ordinarias | Lección #2).
  2. Sería interesante resolver problemas como el que presentas, donde no se conoce la solución analítica, o donde es complicado encontrar una expresión analítica del problema planteado.

Igual te felicito por usar software libre para mostrar las aplicaciones, pocas personas han adoptado esta filosofía.

Saludos y exitos.

Loading...