Ecuación para un cuerpo que cae incluyendo la velocidad terminal

Estoy creando una aplicación que multiplica el tiempo que tarda en caer una piedra y luego calcula la distancia que cayó.

Me di cuenta de que el simple F ( t ) = 1 2 gramo t 2 se estaba volviendo cada vez más inexacto a medida que aumentaba la distancia, por lo que tengo curiosidad por saber si existe una fórmula estándar que explique la velocidad terminal.

Lo que he encontrado hasta ahora

La página de wikipedia sobre la velocidad terminal enumera la fórmula para la velocidad terminal ( V t ) como:

V t = 2 metro gramo ρ A C d

Estoy aproximando que pesa el guijarro metro = 5 gramo con un área proyectada A = 2 C metro 2 . C d = 0.47 para esferas y ρ = 1.204 en 20 C , por lo que nos da:

V t = 2 0.005 9.81 1.204 0.0004 0.47 20.82 metro / s

Parece bastante honesto hasta ahora. Pero ahora necesito combinar F ( t ) = 1 2 gramo t 2 con V t = 20.82 metro / s entonces F ( t ) "crece asintóticamente" a V t . No sé cómo hacer eso, pero jugando con el programa de gráficos de mi computadora obtuve esto:

F ( t ) = ( v t 1 5 ) X 2 + 1 X 2 + 5 1 5

gráfica de la gravedad ideal frente a la fórmula que inventé

( 1 2 gramo t 2 es verde, V t es discontinua, mi fórmula inventada es azul)

Mi fórmula aproximada parece... ¿cerca? Podría tomar algunas medidas y validar esta nueva fórmula experimentalmente, pero no puedo imaginar que sea la primera persona que necesite aproximar la distancia dada en el tiempo para la caída de un objeto.

Además: este gráfico dejó bastante claro que después de ~ 2 segundos de la caída del guijarro 1 2 gramo t 2 (la línea verde) comienza a volverse extremadamente imprecisa.

tl; dr : ¿cuál es la fórmula para la velocidad de un objeto que cae con respecto al tiempo dado? gramo , ρ , A , y C d ?

Has probado a modificar gramo tal que representa la fuerza de arrastre? Es decir, ¿por qué no intentarlo? gramo gramo F d ( v ( t ) ) / metro ?
@honeste_vivere Oh, interesante, parece que debería funcionar. Lo intentaré ahora - F d ( ) parece que podría ser difícil de conseguir.
en.wikipedia.org/wiki/Trajectory_of_a_projectile , "Descripción analítica simple del movimiento de un proyectil en un medio con fuerza de arrastre cuadrática", PS Chudinov, VA Eltyshev, Yu. A. Barykin, "Movimiento de proyectil con resistencia del aire cuadrático en la velocidad" GW Parker...
Entonces gramo ahora es una función del tiempo, obtienes una solución más precisa si expresas v = v 0 + gramo d t
tal vez esto pueda ser útil keisan.casio.com/exec/system/1231475371

Respuestas (3)

Primero, aclaremos el origen de las diferentes expresiones de velocidad terminal y velocidad en función del tiempo para un cuerpo que cae.

Se espera que la fuerza de fricción sea una función creciente de la velocidad del cuerpo y, como resultado, existe una velocidad para la cual esta fuerza equilibra exactamente la gravedad. metro gramo . Ahora bien, para calcular esta velocidad terminal, o para recuperar t v ( t ) , es necesario saber más sobre la expresión de la fuerza de rozamiento. Los modelos más comunes de esta fuerza son:

  • La ley cuadrática, F = 1 2 A C d ρ v 2 . Esas son las mismas variables que en tu pregunta. Desde C d es aproximadamente constante en el régimen de alta velocidad, la v 2 factor representa aproximadamente toda la dependencia de la velocidad, y la fuerza es proporcional a v 2 . También podemos escribir F = k s q tu a r mi v 2 . Pero entonces, esto se aplica a movimientos lo suficientemente rápidos, aquellos para los cuales el número de Reynolds R mi es largo. Como es proporcional a la inversa de la viscosidad del medio, y como el aire de la atmósfera es poco viscoso/denso, R mi es efectivamente grande para la mayoría de los experimentos de caída en el aire y esta es la expresión correcta para F
  • la ley lineal, F = k yo i norte mi a r v . Esta ley se aplica al régimen de baja velocidad y es principalmente relevante para fluidos de alta viscosidad.

Si bien la primera expresión produce la misma expresión que empleó para la velocidad terminal (la raíz cuadrada se origina a partir de la potencia cuadrada en v en F ( v ) ), no es consistente con una dependencia temporal exponencial. La expresión correcta es, según lo dado por wikipedia :

v ( t ) = 2 metro gramo ρ A C d bronceado ( t gramo ρ C d A 2 metro )
La evolución temporal exponencial de la velocidad se obtendría utilizando la expresión lineal de la fuerza (entonces la ecuación de movimiento es una ecuación diferencial lineal de primer orden en t v ( t ) ).

Puedes integrar la expresión correcta para v ( t ) recuperar la distancia de caída con el tiempo (será de la forma d ( t ) en ( aporrear ( t / τ ) ) )

Con respecto a su comparación en http://i.stack.imgur.com/j7Uhh.png , no tiene sentido, ya que está comparando una velocidad (que tiende a ser constante) con una distancia ( 1 2 gramo t 2 ).

Así que creo que tengo una respuesta, pero no estoy completamente convencido de que sea la correcta.

Un amigo (¡gracias Ricky!) me señaló varios lugares en línea donde la gente había resuelto este problema:

Cada gráfica de v ( t ) contra t se ve más o menos así:

gráfico típico de v(t) frente a t: crecimiento lineal de 0 a alrededor de 1,5 segundos, luego crecimiento asintótico a 4 segundos y plano a 8 m/s después de 4 segundos

y sigue la fórmula:

v ( t ) = v t ( 1 mi t τ )

dónde τ = v t gramo y v t es la constante de velocidad terminal calculada en mi pregunta original.

Nota al margen: la posición se obtiene integrando esa fórmula, lo que te da:

y ( t ) = v t ( t + τ ( mi t τ 1 ) )

Eso se ve bien hasta ahora, y es probablemente lo más cerca que voy a estar sin entrar en matemáticas locas, pero la parte que no me sienta bien es cuando esta función de velocidad se grafica con el ideal 1 2 gramo t 2 función:

gráfico de esta función de velocidad que es más rápida que la función ideal durante los primeros 1,5 segundos

La función ideal es más rápida que la función ideal durante los primeros ~1,5 segundos. Eso no tiene un sentido intuitivo para mí: ¿ cómo podría la adición de la fuerza de arrastre hacer que la velocidad inicialmente se acelere más rápido?

De todos modos, haré un trabajo de campo tomando datos experimentales de dejar caer guijarros desde lugares altos y bajos y crearé una función híbrida que se aproxime a la realidad lo más cerca que pueda. Puedo terminar haciendo una función por partes que usa la función ideal durante los primeros 1.5 segundos, luego cambia a esta nueva función asintótica.

Comenzaría midiendo el guijarro y pesándolo para determinar su tamaño. Luego preguntaría "¿Simplemente lo estoy dejando caer?" como Galileo o estoy tratando de hacer una inferencia basada en lanzarlo y luego estimar los efectos al aterrizar.

Si es lo primero, entonces debería obtener un número, ya que no importa el tamaño de la roca, todas caen a la misma velocidad.