Empuje Falcon 9 Merlin 1d calculado a través de cada momento de vuelo

Actualmente estoy tratando de recrear el lanzamiento de Falcon 9 en Unity. Para los cálculos necesito el empuje de la primera etapa (solo lo hago para una etapa por ahora).

¿Cómo calcularlo en función de la altitud? De valores solo tengo los de SpaceX (a nivel del mar y en el vacío), pero necesito valores en cada momento del vuelo. ¿Alguna sugerencia?

Estoy haciendo una simulación simple en 2D, donde el cohete simplemente se lanza directamente hacia arriba.

valores spacex

Consulte aquí space.stackexchange.com/a/43845/6944 y los enlaces que contiene.
Pero quiero decir, ¿cómo obtener valores en cualquier momento del vuelo? Para los cálculos, necesito una columna para los cambios de empuje, durante todo el vuelo. ¿Alguna sugerencia sobre eso?
Tenga en cuenta que los motores no funcionan a toda velocidad durante todo el vuelo. Se reducen alrededor de max-q, y también se pueden reducir cerca del final de la quema para reducir la carga g.

Respuestas (2)

Voy a expresar esto en unidades inglesas porque eso es lo que tengo claro y, por lo tanto, es menos probable que cometa un error vergonzoso.

Usemos la ecuación de empuje para motores que no respiran aire, copiada de aquí De la ecuación de empuje general hacia Tsiolkovsky, ¿cómo explicar la caída de estos términos en el camino?

      F = metro ˙ mi V mi + ( pag mi pag 0 ) A mi

El primer término a la derecha de = es el impulso de impulso . Es constante para una configuración de aceleración dada para lo que estamos haciendo aquí.

Das:

  • F v = 1,849,500 lbf
  • F s = 1,710,000 lbf

Podemos resolver el área del plano de salida completando lo que sabemos en las dos ecuaciones y restando una de la otra. Eso nos da un área del plano de salida (para los nueve motores) de 66 pie 2 .

Una respuesta (sin fuente) a esta pregunta La temperatura y la presión de los gases de escape del cohete dan como resultado una presión del avión de salida de Merlín de 0,7 atm.

Ahora podemos calcular el término de impulso de impulso metro ˙ mi V mi ser 1,751,703 lbf.

Comprobación rápida de la cordura: el impulso de empuje que calcule debe ser mayor que el empuje del nivel del mar y menor que el empuje del vacío.

Con eso y una mesa de ambiente práctico y elegante, tenemos todo lo que necesitamos.

Redondeé todo para que los números no coincidan exactamente. Solo repase lo que hice, usando las unidades de su elección. Luego, en su programa, use el empuje de impulso, la presión del plano de salida y el área del plano de salida que calculó en la ecuación de empuje, e introduzca una presión ambiental que obtuvo con el método que eligió.

Altitud Presión ambiental Empuje
0 pies 2117 lbf pie 2 1,709,793 lbf
50K pies 241 lbf pie 2 1,833,609 lbf
100K pies 23 lbf pie 2 1,847,997 lbf
infinito y más allá 0 lbf pie 2 1,849,515 lbf

Gráfico de empuje (lbf) vs altitud (ft)

ingrese la descripción de la imagen aquí

Y por cierto, todavía tengo una pregunta aquí. ¿Cómo obtuvo el área del plano de salida para todos los motores igual a 66 pies ^ 2? Usando la ecuación F(vacío) - F(nivel del mar) = p0*Ae obtuve Ae igual a 6.119 m^2, que es mucho menor que el tuyo Ae. Además, si observamos los números, podemos ver que falcon 9 El área de la sección transversal es +- 42 m^2 - mucho MÁS que el mío Ae, pero menor que el tuyo. ¿Me estoy perdiendo algo en la ecuación?
@mad.redhead recuerda que usé unidades inglesas. 66 pies ^ 2 está bastante cerca de 6 m ^ 2.
pero ¿por qué es esa gran diferencia entre Ae y el área de la sección transversal?
@ mad.redhead Este artículo de wikipedia le da al diámetro de salida de Merlín 3 pies. Eso es un poco más de 7 pies ^ 2, por lo que 9 de ellos tienen alrededor de 63 pies ^ 2. ¡Bastante cerca de mis 66! en.wikipedia.org/wiki/SpaceX_Merlin
Oooh, no eran metros. De acuerdo... estoy realmente confundido con otras personas y conmigo mismo también. Necesito estar más atento (

Con respecto a la dependencia de un empuje de una altitud, resolviendo una tarea similar, uso la fórmula barométrica para obtener la presión ambiental a una altitud dada. Puede obtener una diferencia entre el empuje conocido (o impulso específico, si lo desea) en el vacío ( Th1 ) y al nivel del mar ( Th0 ). Esta diferencia dTh = Th1 - Th0 se multiplicará por un coeficiente K dependiente de la altitud , de forma que a nivel del mar este coeficiente será igual a 1 y cero en el vacío. Luego, el producto del coeficiente y la diferencia de empuje se resta del empuje (o impulso específico) en el vacío: Thrust_at_given_altitude = Th1 - dTh * K .

Necesita obtener este coeficiente K . Usemos la fórmula barométrica para esto: K = presión_calculada_por_fórmula_barométrica / presión_a_nivel_del_mar .

Utilizo ambas fórmulas superiores para la presión en el artículo de Wikpedia , según la altitud de referencia, y la tabla para la altitud de referencia, por lo que la dependencia de la altitud de empuje tiene el patrón de escalón.

No estoy seguro de que sea una forma óptima, da resultados plausibles en mi caso (utilizo el programa GMAT para el modelado de ascenso de cohetes). Pero todavía estoy pensando en la corrección metódica de esta manera.

Otros enlaces útiles: Densidad del aire , Atmósfera estándar internacional .

Es un poco diferente de lo que esperaba calcular, ¡pero también lo intentaré de esta manera! Gracias por la respuesta, ¡ustedes realmente me están ayudando!