Sin ánimo de incomodarle señor John Carlin, deja usted al periódico "El País" y a los españoles en mal lugar.
Dice en su artículo que en el despegue el cohete Saturno V tenía un consumo de 3785 litros de combustible por segundo.
http://elpais.com/diario/2009/07/19/domingo/1247975558_850215.html
Soy ingeniero y le voy a demostrar rápidamente que ese número está EQUIVOCADO.
Voy a asumir el caso más favorable para su afirmación, que cuando habla usted de combustible SOLAMENTE se refiere al RP-1 y no a la mezcla RP-1 con el oxidante LOX. Si usted se refiriese a la mezcla el error sería muchísimo peor.
3785 litros de RP-1 son 3785*0.806 = 3050.71 kg
2.27 kg de LOX se mezclan con 1 kg de RP-1
3050.71 kg de RP-1 se mezclan con 3050.71*2.27 = 6925.1117 kg de LOX
En total tenemos 3050.71+6925.1117 = 9975.8217 kg/s
Teniendo en cuenta el impulso específico 265.4 segundos de los 5 motores F-1 de la primera etapa, ya podemos calcular el empuje obtenido en el momento del lanzamiento
Empuje = (9975.8217 kg/s)*(265.4*9.81 m/s) = 25972790.0068 N =~ 25.97E6 N
Un rápido cálculo del peso del Saturno V en el lanzamiento en la misión Apolo 11 resulta
Peso = 2941.12E3*9.81 = 28852387.2 N =~ 28.85E6 N
Como se puede comprobar EL PESO ES MAYOR QUE EL EMPUJE.
Por tanto con sus 3785 litros de combustible por segundo en el lanzamiento, cuando la cuenta atrás llegase a cero, el Saturno V se quedaría clavado sin elevarse ni un solo milímetro. Estaría así SIN MOVERSE desperdiciando combustible hasta que el empuje comenzase a ser mayor que el peso. Eso ocurriría exactamente después de (28852387.2-25972790.0068)/(9.81*9975.8217) = 29.43 segundos, casi medio minuto. Después de todo ese tiempo empezaría a elevarse EXTREMADAMENTE DESPACIO. Pero voy a dejar los cálculos aquí.
Ojalá este mensaje le sirva de ayuda.
Reciba usted un cordial saludo desde Madrid.
Sunday, August 26, 2012
Tuesday, August 14, 2012
In the Apollo 11 mission the Saturn V didn't have sufficient fuel
I have demonstrated [not in the strict mathematical sense] with the previous message (multistage.c) that the Saturn V in the Apollo 11 mission didn't have sufficient fuel to reach the escape velocity towards the Moon. Hereunder I comment briefly as it has been demonstrated.
Suppose that the Moon in its orbit is at the minimum distance from the Earth, that is the case most favorable or of the least fuel consumption. Suppose also that the Moon doesn't move relative to the Earth. Although this doesn't have physical sense it has mathematical sense, because it is a case more favorable or of less fuel consumption, as the Moon exerts its maximum attraction over the Saturn V constantly and this facilitates the ascent of the Saturn V. Suppose that the Saturn V travels in a straight line from the Earth towards the Moon. This is the trajectory most favorable or of the least fuel consumption to reach the escape velocity towards the Moon. Remember that the gravitational force is conservative. Besides, although the acceleration suffered by the astronauts in the Apollo 11 mission was variable below 5*g, most of the time below 3*g, suppose that the Saturn V moves in such a manner that the astronauts suffer an acceleration much higher constant 5*g, that sustained during more than 30 seconds surely is deadly. As the ascent acceleration is higher during all the time that the engines are turned on, the fuel consumption necessary to reach the escape velocity towards the Moon is less.
It results that IN THIS CASE MUCH MORE FAVORABLE (Moon at minimum distance, Moon without movement, Saturn V travels in a straight line, Saturn V travels with acceleration much higher constant 5*g over the astronauts) THE SATURN V DOESN'T HAVE SUFFICIENT FUEL TO REACH THE MOON. This can be checked running the program multistage.c of the previous message. As immediate consequence in any case less favorable than this the Saturn V doesn't have sufficient fuel to reach the Moon. THEREFORE, IN THE APOLLO 11 MISSION THE SATURN V DIDN'T HAVE SUFFICIENT FUEL TO REACH THE MOON.
The absence of fuel is produced in the stage 3 (S-IVB stage), when the J-2 engine exhausts its 107572 kg (237155 pounds) of fuel before to reach the escape velocity towards the Moon.
[In Spanish]
He demostrado [no en sentido matemático estricto] con el mensaje anterior (multistage.c) que el Saturn V en la misión Apollo 11 no tenía combustible suficiente para alcanzar la velocidad de escape hacia la Luna. A continuación comento brevemente como se ha demostrado.
Supón que la Luna en su órbita está a la distancia mínima de la Tierra, que es el caso más favorable o de menor gasto de combustible. Supón también que la Luna no se mueve respecto a la Tierra. Aunque esto último no tiene sentido físico sí tiene sentido matemático, porque es un caso más favorable o de menor gasto de combustible, ya que la Luna ejerce su atracción máxima sobre el Saturn V constantemente y esto facilita el ascenso del Saturn V. Supón que el Saturn V viaja en línea recta desde la Tierra hacia la Luna. Esta es la trayectoria más favorable o de menor gasto de combustible para alcanzar la velocidad de escape hacia la Luna. Recuerda que la fuerza gravitacional es conservativa. Además, aunque la aceleración sufrida por los astronautas en la misión Apollo 11 era variable por debajo de 5*g, la mayor parte del tiempo por debajo de 3*g, supón que el Saturn V se mueve de tal manera que los astronautas sufren una aceleración mucho mayor constante 5*g, que sostenida durante más de 30 segundos seguramente es mortal. Ya que es mayor la aceleración de ascenso durante todo el tiempo que los motores están encendidos, el consumo de combustible necesario para alcanzar la velocidad de escape hacia la Luna es menor.
Resulta que EN ESTE CASO MUCHO MÁS FAVORABLE (Luna a distancia mínima, Luna sin movimiento, Saturn V viaja en línea recta, Saturn V viaja con aceleración mucho mayor contante 5*g sobre los astronautas) EL SATURN V NO TIENE SUFICIENTE COMBUSTIBLE PARA ALCANZAR LA LUNA. Esto se puede comprobar ejecutando el programa multistage.c del mensaje anterior. Como consecuencia inmediata en cualquier caso menos favorable que este el Saturn V no tiene suficiente combustible para alcanzar la Luna. POR TANTO, EN LA MISIÓN APOLLO 11 EL SATURN V NO TENÍA SUFICIENTE COMBUSTIBLE PARA ALCANZAR LA LUNA.
La falta de combustible se produce en la etapa 3 (etapa S-IVB), cuando el motor J-2 agota sus 107572 kg (237155 pounds) de combustible antes de alcanzar la velocidad de escape hacia la Luna.
Suppose that the Moon in its orbit is at the minimum distance from the Earth, that is the case most favorable or of the least fuel consumption. Suppose also that the Moon doesn't move relative to the Earth. Although this doesn't have physical sense it has mathematical sense, because it is a case more favorable or of less fuel consumption, as the Moon exerts its maximum attraction over the Saturn V constantly and this facilitates the ascent of the Saturn V. Suppose that the Saturn V travels in a straight line from the Earth towards the Moon. This is the trajectory most favorable or of the least fuel consumption to reach the escape velocity towards the Moon. Remember that the gravitational force is conservative. Besides, although the acceleration suffered by the astronauts in the Apollo 11 mission was variable below 5*g, most of the time below 3*g, suppose that the Saturn V moves in such a manner that the astronauts suffer an acceleration much higher constant 5*g, that sustained during more than 30 seconds surely is deadly. As the ascent acceleration is higher during all the time that the engines are turned on, the fuel consumption necessary to reach the escape velocity towards the Moon is less.
It results that IN THIS CASE MUCH MORE FAVORABLE (Moon at minimum distance, Moon without movement, Saturn V travels in a straight line, Saturn V travels with acceleration much higher constant 5*g over the astronauts) THE SATURN V DOESN'T HAVE SUFFICIENT FUEL TO REACH THE MOON. This can be checked running the program multistage.c of the previous message. As immediate consequence in any case less favorable than this the Saturn V doesn't have sufficient fuel to reach the Moon. THEREFORE, IN THE APOLLO 11 MISSION THE SATURN V DIDN'T HAVE SUFFICIENT FUEL TO REACH THE MOON.
The absence of fuel is produced in the stage 3 (S-IVB stage), when the J-2 engine exhausts its 107572 kg (237155 pounds) of fuel before to reach the escape velocity towards the Moon.
[In Spanish]
He demostrado [no en sentido matemático estricto] con el mensaje anterior (multistage.c) que el Saturn V en la misión Apollo 11 no tenía combustible suficiente para alcanzar la velocidad de escape hacia la Luna. A continuación comento brevemente como se ha demostrado.
Supón que la Luna en su órbita está a la distancia mínima de la Tierra, que es el caso más favorable o de menor gasto de combustible. Supón también que la Luna no se mueve respecto a la Tierra. Aunque esto último no tiene sentido físico sí tiene sentido matemático, porque es un caso más favorable o de menor gasto de combustible, ya que la Luna ejerce su atracción máxima sobre el Saturn V constantemente y esto facilita el ascenso del Saturn V. Supón que el Saturn V viaja en línea recta desde la Tierra hacia la Luna. Esta es la trayectoria más favorable o de menor gasto de combustible para alcanzar la velocidad de escape hacia la Luna. Recuerda que la fuerza gravitacional es conservativa. Además, aunque la aceleración sufrida por los astronautas en la misión Apollo 11 era variable por debajo de 5*g, la mayor parte del tiempo por debajo de 3*g, supón que el Saturn V se mueve de tal manera que los astronautas sufren una aceleración mucho mayor constante 5*g, que sostenida durante más de 30 segundos seguramente es mortal. Ya que es mayor la aceleración de ascenso durante todo el tiempo que los motores están encendidos, el consumo de combustible necesario para alcanzar la velocidad de escape hacia la Luna es menor.
Resulta que EN ESTE CASO MUCHO MÁS FAVORABLE (Luna a distancia mínima, Luna sin movimiento, Saturn V viaja en línea recta, Saturn V viaja con aceleración mucho mayor contante 5*g sobre los astronautas) EL SATURN V NO TIENE SUFICIENTE COMBUSTIBLE PARA ALCANZAR LA LUNA. Esto se puede comprobar ejecutando el programa multistage.c del mensaje anterior. Como consecuencia inmediata en cualquier caso menos favorable que este el Saturn V no tiene suficiente combustible para alcanzar la Luna. POR TANTO, EN LA MISIÓN APOLLO 11 EL SATURN V NO TENÍA SUFICIENTE COMBUSTIBLE PARA ALCANZAR LA LUNA.
La falta de combustible se produce en la etapa 3 (etapa S-IVB), cuando el motor J-2 agota sus 107572 kg (237155 pounds) de combustible antes de alcanzar la velocidad de escape hacia la Luna.
Thursday, August 09, 2012
multistage.c
#define AA 49.05 // aceleración soportada por los astronautas 5*g (letal en pocos segundos)
#define SPECIFIC_IMPULSE_F1 265.4
#define SPECIFIC_IMPULSE_J2 426.0
#define SPECIFIC_IMPULSE_SPACECRAFT 452.0 // inventado por el momento (se ha buscado el valor más alto posible)
#define EV_F1 SPECIFIC_IMPULSE_F1*9.81
#define EV_J2 SPECIFIC_IMPULSE_J2*9.81
#define EV_XX SPECIFIC_IMPULSE_SPACECRAFT*9.81
#define THRUST_F1 6.7725E6
#define THRUST_J2 1033.5E3
#define MC 2941.12E3 // masa del cohete Saturn V (misión Apollo 11) 2941.12 toneladas
#define S_IC_STAGE_EXPENDABLES 2149.72E3
#define FIRST_STAGE_INERT_WEIGHT 130975.0+5200.44
#define S_II_STAGE_EXPENDABLES 444752.0
#define SECOND_STAGE_INERT_WEIGHT 36251.1+3665.03
#define S_IVB_STAGE_EXPENDABLES 107572.0
#define THIRD_STAGE_INERT_WEIGHT 11339.8+1952.72
#define SPACECRAFT_INERT_WEIGHT 20538.7
#define DT 0.1 // incremento de tiempo por iteración de cálculo 0.1 segundos
#define MT 5.97E24 // masa de la Tierra
#define RT 6.379E6 // radio de la Tierra
#define ML 7.35E22 // masa de la Luna
#define RL 1.739E6 // radio de la Luna
#define DTL 356400E3 // distancia Tierra-Luna (entre los centros de masas)
#define GG 6.67259E-11 // constante de gravitación universal
#include <math.h>
#include <stdio.h>
#include <stdlib.h>
int main(int argc, char **argv) {
double mm; // masa restante del cohete
double mg; // masa de gas expulsada en iteración
double mc; // masa de combustible consumido en la etapa
double mi; // masa inerte en la etapa
double ee; // empuje del cohete
double ff; // fuerza sobre el cohete
double aa; // aceleración del cohete
double vvi; // velocidad inicial del cohete
double vvf; // velocidad final del cohete
double rri; // posición inicial del cohete
double rrf; // posición final del cohete
double tt; // tiempo
double rr2; // posición donde se anulan las fuerzas
double vescape; // velocidad de escape
double EV; // "exhaust velocity" de los gases
int etapa;
rr2 = (DTL*(MT-sqrt(MT*ML)))/(MT-ML);
tt = 0.0;
mm = MC;
mc = 0.0;
vvi = 0.0;
rri = RT;
etapa = 1;
EV = EV_F1;
mi = FIRST_STAGE_INERT_WEIGHT + SECOND_STAGE_INERT_WEIGHT + THIRD_STAGE_INERT_WEIGHT + SPACECRAFT_INERT_WEIGHT;
do {
mg = (DT*mm*AA)/(EV+DT*AA);
mm = mm - mg;
mc = mc + mg;
ee = mg*EV/DT;
ff = ee + GG*(mm*ML)/pow(DTL-rri,2.0) - GG*(mm*MT)/pow(rri,2.0);
aa = ff/mm;
vvf = aa*DT + vvi;
rrf = aa*pow(DT,2.0)/2.0 + vvi*DT + rri;
vescape = sqrt(2.0*GG*(MT*(1.0/rrf - 1.0/rr2) + ML*(1.0/(DTL-rrf) - 1.0/(DTL-rr2)))); // MAL
tt = tt + DT;
printf("Etapa : %d\n", etapa);
printf("Tiempo : %g [s]\n", tt);
printf("Velocidad de salida de los gases respecto al cohete : %g [m/s]\n", EV);
printf("Masa de gas expulsada por segundo : %g [kg/s]\n", mg/DT);
printf("Masa restante del cohete : %g [kg]\n", mm);
printf("Masa restante del cohete (porcentaje) : %g [%%]\n", mm/MC*100.0);
printf("Masa inerte del cohete en esta etapa: %g [kg]\n", mi);
if (etapa == 1) {
printf("Empuje del cohete : %g [N] (%g motores F-1)\n", ee, ee/THRUST_F1);
} else if ((etapa == 2) || (etapa == 3)) {
printf("Empuje del cohete : %g [N] (%g motores J-2)\n", ee, ee/THRUST_J2);
} else {
printf("Empuje del cohete : %g [N]\n", ee);
}
printf("Fuerza total sobre el cohete : %g [N]\n", ff);
printf("Aceleración sufrida por los astronautas : %g [m/(s^2)]\n", AA);
printf("Aceleración del cohete : %g [m/(s^2)]\n", aa);
printf("Velocidad del cohete : %g [m/s]\n", vvf);
printf("Velocidad del cohete (km/h) : %g [km/h]\n", vvf/1000.0*3600.0);
printf("Velocidad de escape : %g [m/s]\n", vescape);
printf("Velocidad de escape (km/h) : %g [km/h]\n", vescape/1000.0*3600.0);
printf("Posición del cohete : %g [m]\n", rrf);
printf("Altura del cohete sobre la superficie de la Tierra : %g [m]\n", rrf-RT);
printf("Altura del cohete sobre la superficie de la Tierra (km) : %g [km]\n", (rrf-RT)/1000.0);
printf("Distancia del cohete a la superficie de la Luna : %g [m]\n", DTL-rrf-RL);
printf("Distancia del cohete a la superficie de la Luna (km) : %g [km]\n\n", (DTL-rrf-RL)/1000.0);
switch (etapa) {
case 1:
if (mc >= S_IC_STAGE_EXPENDABLES) {
printf("AGOTADO COMBUSTIBLE DE LA PRIMERA ETAPA: %g [kg]\n", S_IC_STAGE_EXPENDABLES);
printf("APAGADOS LOS 5 MOTORES F-1 DE LA PRIMERA ETAPA\n");
mm = mm - FIRST_STAGE_INERT_WEIGHT;
mi = SECOND_STAGE_INERT_WEIGHT + THIRD_STAGE_INERT_WEIGHT + SPACECRAFT_INERT_WEIGHT;
printf("LIBERADO PESO INERTE DE LA PRIMERA ETAPA: %g [kg]\n", FIRST_STAGE_INERT_WEIGHT);
printf("Masa restante del cohete : %g [kg]\n", mm);
printf("Masa restante del cohete (porcentaje) : %g [%%]\n", mm/MC*100.0);
EV = EV_J2;
printf("ENCENDIDOS LOS 5 MOTORES J-2 DE LA SEGUNDA ETAPA\n\n");
mc = 0.0;
etapa = 2;
}
break;
case 2:
if (mc >= S_II_STAGE_EXPENDABLES) {
printf("AGOTADO COMBUSTIBLE DE LA SEGUNDA ETAPA: %g [kg]\n", S_II_STAGE_EXPENDABLES);
printf("APAGADOS LOS 5 MOTORES J-2 DE LA SEGUNDA ETAPA\n");
mm = mm - SECOND_STAGE_INERT_WEIGHT;
mi = THIRD_STAGE_INERT_WEIGHT + SPACECRAFT_INERT_WEIGHT;
printf("LIBERADO PESO INERTE DE LA SEGUNDA ETAPA: %g [kg]\n", SECOND_STAGE_INERT_WEIGHT);
printf("Masa restante del cohete : %g [kg]\n", mm);
printf("Masa restante del cohete (porcentaje) : %g [%%]\n", mm/MC*100.0);
printf("ENCENDIDO EL MOTOR J-2 DE LA TERCERA ETAPA\n\n");
mc = 0.0;
etapa = 3;
}
break;
case 3:
if (mc >= S_IVB_STAGE_EXPENDABLES) {
printf("AGOTADO COMBUSTIBLE DE LA TERCERA ETAPA: %g [kg]\n", S_IVB_STAGE_EXPENDABLES);
return 0;
printf("APAGADO MOTOR J-2 DE LA TERCERA ETAPA\n");
mm = mm - THIRD_STAGE_INERT_WEIGHT;
mi = SPACECRAFT_INERT_WEIGHT;
printf("LIBERADO PESO INERTE DE LA TERCERA ETAPA: %g [kg]\n", THIRD_STAGE_INERT_WEIGHT);
printf("Masa restante del cohete : %g [kg]\n", mm);
printf("Masa restante del cohete (porcentaje) : %g [%%]\n", mm/MC*100.0);
EV = EV_XX;
printf("ENCENDIDO EL MOTOR DE LA CUARTA ETAPA (ASTRONAVE)\n");
printf("¡¡¡USANDO COMBUSTIBLE NECESARIO PARA EL VIAJE DE RETORNO!!!\n\n");
mc = 0.0;
etapa = 4;
}
}
vvi = vvf;
rri = rrf;
} while ((vvf < vescape) && (mm > mi));
if (vvf >= vescape) {
printf("ALCANZADA VELOCIDAD DE ESCAPE\n\n");
} else if (mm <= mi) {
printf("COMBUSTIBLE ACABADO ANTES DE ALCANZAR LA VELOCIDAD DE ESCAPE\n\n");
}
return 0;
}
Etapa : 3
Tiempo : 235.852 [s]
Velocidad de salida de los gases respecto al cohete : 4179.06 [m/s]
Masa de gas expulsada por segundo : 947.176 [kg/s]
Masa restante del cohete : 80699.4 [kg]
Masa restante del cohete (porcentaje) : 2.74383 [%]
Masa inerte del cohete en esta etapa: 33831.2 [kg]
Empuje del cohete : 3.95831e+006 [N] (3.83 motores J-2)
Fuerza total sobre el cohete : 3.38454e+006 [N]
Aceleración sufrida por los astronautas : 49.05 [m/(s^2)]
Aceleración del cohete : 41.94 [m/(s^2)]
Velocidad del cohete : 9489.73 [m/s]
Velocidad del cohete (km/h) : 34163 [km/h]
Velocidad de escape : 10183.7 [m/s] MAL
Velocidad de escape (km/h) : 36661.3 [km/h]
Posición del cohete : 7.48513e+006 [m]
Altura del cohete sobre la superficie de la Tierra : 1.10613e+006 [m]
Altura del cohete sobre la superficie de la Tierra (km) : 1106.13 [km]
Distancia del cohete a la superficie de la Luna : 3.47176e+008 [m]
Distancia del cohete a la superficie de la Luna (km) : 347176 [km]
AGOTADO COMBUSTIBLE DE LA TERCERA ETAPA: 107572 [kg]
LA EXPRESIÓN DE LA VELOCIDAD DE ESCAPE ESTÁ MAL.
TAMPOCO SE PUEDE DESPRECIAR EL EFECTO DEL SOL.
NOTA: He usado DT con valor 0.001 segundos en lugar de 0.1 segundos.
NOTA: He actualizado el impulso específico de los motores F-1 y J-2 con unos valores un poco más altos que he encontrado en un documento de la NASA. Pongo a continuación los enlaces a los parámetros del motor F-1 y del motor J-2 extraídos de ese documento.
#define SPECIFIC_IMPULSE_F1 265.4
#define SPECIFIC_IMPULSE_J2 426.0
#define SPECIFIC_IMPULSE_SPACECRAFT 452.0 // inventado por el momento (se ha buscado el valor más alto posible)
#define EV_F1 SPECIFIC_IMPULSE_F1*9.81
#define EV_J2 SPECIFIC_IMPULSE_J2*9.81
#define EV_XX SPECIFIC_IMPULSE_SPACECRAFT*9.81
#define THRUST_F1 6.7725E6
#define THRUST_J2 1033.5E3
#define MC 2941.12E3 // masa del cohete Saturn V (misión Apollo 11) 2941.12 toneladas
#define S_IC_STAGE_EXPENDABLES 2149.72E3
#define FIRST_STAGE_INERT_WEIGHT 130975.0+5200.44
#define S_II_STAGE_EXPENDABLES 444752.0
#define SECOND_STAGE_INERT_WEIGHT 36251.1+3665.03
#define S_IVB_STAGE_EXPENDABLES 107572.0
#define THIRD_STAGE_INERT_WEIGHT 11339.8+1952.72
#define SPACECRAFT_INERT_WEIGHT 20538.7
#define DT 0.1 // incremento de tiempo por iteración de cálculo 0.1 segundos
#define MT 5.97E24 // masa de la Tierra
#define RT 6.379E6 // radio de la Tierra
#define ML 7.35E22 // masa de la Luna
#define RL 1.739E6 // radio de la Luna
#define DTL 356400E3 // distancia Tierra-Luna (entre los centros de masas)
#define GG 6.67259E-11 // constante de gravitación universal
#include <math.h>
#include <stdio.h>
#include <stdlib.h>
int main(int argc, char **argv) {
double mm; // masa restante del cohete
double mg; // masa de gas expulsada en iteración
double mc; // masa de combustible consumido en la etapa
double mi; // masa inerte en la etapa
double ee; // empuje del cohete
double ff; // fuerza sobre el cohete
double aa; // aceleración del cohete
double vvi; // velocidad inicial del cohete
double vvf; // velocidad final del cohete
double rri; // posición inicial del cohete
double rrf; // posición final del cohete
double tt; // tiempo
double rr2; // posición donde se anulan las fuerzas
double vescape; // velocidad de escape
double EV; // "exhaust velocity" de los gases
int etapa;
rr2 = (DTL*(MT-sqrt(MT*ML)))/(MT-ML);
tt = 0.0;
mm = MC;
mc = 0.0;
vvi = 0.0;
rri = RT;
etapa = 1;
EV = EV_F1;
mi = FIRST_STAGE_INERT_WEIGHT + SECOND_STAGE_INERT_WEIGHT + THIRD_STAGE_INERT_WEIGHT + SPACECRAFT_INERT_WEIGHT;
do {
mg = (DT*mm*AA)/(EV+DT*AA);
mm = mm - mg;
mc = mc + mg;
ee = mg*EV/DT;
ff = ee + GG*(mm*ML)/pow(DTL-rri,2.0) - GG*(mm*MT)/pow(rri,2.0);
aa = ff/mm;
vvf = aa*DT + vvi;
rrf = aa*pow(DT,2.0)/2.0 + vvi*DT + rri;
vescape = sqrt(2.0*GG*(MT*(1.0/rrf - 1.0/rr2) + ML*(1.0/(DTL-rrf) - 1.0/(DTL-rr2)))); // MAL
tt = tt + DT;
printf("Etapa : %d\n", etapa);
printf("Tiempo : %g [s]\n", tt);
printf("Velocidad de salida de los gases respecto al cohete : %g [m/s]\n", EV);
printf("Masa de gas expulsada por segundo : %g [kg/s]\n", mg/DT);
printf("Masa restante del cohete : %g [kg]\n", mm);
printf("Masa restante del cohete (porcentaje) : %g [%%]\n", mm/MC*100.0);
printf("Masa inerte del cohete en esta etapa: %g [kg]\n", mi);
if (etapa == 1) {
printf("Empuje del cohete : %g [N] (%g motores F-1)\n", ee, ee/THRUST_F1);
} else if ((etapa == 2) || (etapa == 3)) {
printf("Empuje del cohete : %g [N] (%g motores J-2)\n", ee, ee/THRUST_J2);
} else {
printf("Empuje del cohete : %g [N]\n", ee);
}
printf("Fuerza total sobre el cohete : %g [N]\n", ff);
printf("Aceleración sufrida por los astronautas : %g [m/(s^2)]\n", AA);
printf("Aceleración del cohete : %g [m/(s^2)]\n", aa);
printf("Velocidad del cohete : %g [m/s]\n", vvf);
printf("Velocidad del cohete (km/h) : %g [km/h]\n", vvf/1000.0*3600.0);
printf("Velocidad de escape : %g [m/s]\n", vescape);
printf("Velocidad de escape (km/h) : %g [km/h]\n", vescape/1000.0*3600.0);
printf("Posición del cohete : %g [m]\n", rrf);
printf("Altura del cohete sobre la superficie de la Tierra : %g [m]\n", rrf-RT);
printf("Altura del cohete sobre la superficie de la Tierra (km) : %g [km]\n", (rrf-RT)/1000.0);
printf("Distancia del cohete a la superficie de la Luna : %g [m]\n", DTL-rrf-RL);
printf("Distancia del cohete a la superficie de la Luna (km) : %g [km]\n\n", (DTL-rrf-RL)/1000.0);
switch (etapa) {
case 1:
if (mc >= S_IC_STAGE_EXPENDABLES) {
printf("AGOTADO COMBUSTIBLE DE LA PRIMERA ETAPA: %g [kg]\n", S_IC_STAGE_EXPENDABLES);
printf("APAGADOS LOS 5 MOTORES F-1 DE LA PRIMERA ETAPA\n");
mm = mm - FIRST_STAGE_INERT_WEIGHT;
mi = SECOND_STAGE_INERT_WEIGHT + THIRD_STAGE_INERT_WEIGHT + SPACECRAFT_INERT_WEIGHT;
printf("LIBERADO PESO INERTE DE LA PRIMERA ETAPA: %g [kg]\n", FIRST_STAGE_INERT_WEIGHT);
printf("Masa restante del cohete : %g [kg]\n", mm);
printf("Masa restante del cohete (porcentaje) : %g [%%]\n", mm/MC*100.0);
EV = EV_J2;
printf("ENCENDIDOS LOS 5 MOTORES J-2 DE LA SEGUNDA ETAPA\n\n");
mc = 0.0;
etapa = 2;
}
break;
case 2:
if (mc >= S_II_STAGE_EXPENDABLES) {
printf("AGOTADO COMBUSTIBLE DE LA SEGUNDA ETAPA: %g [kg]\n", S_II_STAGE_EXPENDABLES);
printf("APAGADOS LOS 5 MOTORES J-2 DE LA SEGUNDA ETAPA\n");
mm = mm - SECOND_STAGE_INERT_WEIGHT;
mi = THIRD_STAGE_INERT_WEIGHT + SPACECRAFT_INERT_WEIGHT;
printf("LIBERADO PESO INERTE DE LA SEGUNDA ETAPA: %g [kg]\n", SECOND_STAGE_INERT_WEIGHT);
printf("Masa restante del cohete : %g [kg]\n", mm);
printf("Masa restante del cohete (porcentaje) : %g [%%]\n", mm/MC*100.0);
printf("ENCENDIDO EL MOTOR J-2 DE LA TERCERA ETAPA\n\n");
mc = 0.0;
etapa = 3;
}
break;
case 3:
if (mc >= S_IVB_STAGE_EXPENDABLES) {
printf("AGOTADO COMBUSTIBLE DE LA TERCERA ETAPA: %g [kg]\n", S_IVB_STAGE_EXPENDABLES);
return 0;
printf("APAGADO MOTOR J-2 DE LA TERCERA ETAPA\n");
mm = mm - THIRD_STAGE_INERT_WEIGHT;
mi = SPACECRAFT_INERT_WEIGHT;
printf("LIBERADO PESO INERTE DE LA TERCERA ETAPA: %g [kg]\n", THIRD_STAGE_INERT_WEIGHT);
printf("Masa restante del cohete : %g [kg]\n", mm);
printf("Masa restante del cohete (porcentaje) : %g [%%]\n", mm/MC*100.0);
EV = EV_XX;
printf("ENCENDIDO EL MOTOR DE LA CUARTA ETAPA (ASTRONAVE)\n");
printf("¡¡¡USANDO COMBUSTIBLE NECESARIO PARA EL VIAJE DE RETORNO!!!\n\n");
mc = 0.0;
etapa = 4;
}
}
vvi = vvf;
rri = rrf;
} while ((vvf < vescape) && (mm > mi));
if (vvf >= vescape) {
printf("ALCANZADA VELOCIDAD DE ESCAPE\n\n");
} else if (mm <= mi) {
printf("COMBUSTIBLE ACABADO ANTES DE ALCANZAR LA VELOCIDAD DE ESCAPE\n\n");
}
return 0;
}
Etapa : 3
Tiempo : 235.852 [s]
Velocidad de salida de los gases respecto al cohete : 4179.06 [m/s]
Masa de gas expulsada por segundo : 947.176 [kg/s]
Masa restante del cohete : 80699.4 [kg]
Masa restante del cohete (porcentaje) : 2.74383 [%]
Masa inerte del cohete en esta etapa: 33831.2 [kg]
Empuje del cohete : 3.95831e+006 [N] (3.83 motores J-2)
Fuerza total sobre el cohete : 3.38454e+006 [N]
Aceleración sufrida por los astronautas : 49.05 [m/(s^2)]
Aceleración del cohete : 41.94 [m/(s^2)]
Velocidad del cohete : 9489.73 [m/s]
Velocidad del cohete (km/h) : 34163 [km/h]
Velocidad de escape : 10183.7 [m/s] MAL
Velocidad de escape (km/h) : 36661.3 [km/h]
Posición del cohete : 7.48513e+006 [m]
Altura del cohete sobre la superficie de la Tierra : 1.10613e+006 [m]
Altura del cohete sobre la superficie de la Tierra (km) : 1106.13 [km]
Distancia del cohete a la superficie de la Luna : 3.47176e+008 [m]
Distancia del cohete a la superficie de la Luna (km) : 347176 [km]
AGOTADO COMBUSTIBLE DE LA TERCERA ETAPA: 107572 [kg]
LA EXPRESIÓN DE LA VELOCIDAD DE ESCAPE ESTÁ MAL.
TAMPOCO SE PUEDE DESPRECIAR EL EFECTO DEL SOL.
NOTA: He usado DT con valor 0.001 segundos en lugar de 0.1 segundos.
NOTA: He actualizado el impulso específico de los motores F-1 y J-2 con unos valores un poco más altos que he encontrado en un documento de la NASA. Pongo a continuación los enlaces a los parámetros del motor F-1 y del motor J-2 extraídos de ese documento.
Thursday, August 02, 2012
La altura de los saltos en la Luna
Las personas en la Luna pesan aproximadamente un 16.57 % de lo que pesan en la Tierra. Esto se debe a que la Luna tiene menos masa que la Tierra.
(ML*RT^2)/(MT*RL^2)*100 = 16.56605307 %
Un astronauta de 80 kg en la Luna tendría la impresión de que tiene una masa de 13.25 kg, y obviamente seguiría siendo igual de grande e igual de fuerte.
Un salto vertical de hT metros de altura realizado por una persona en la Tierra en la Luna tendría una altura hL aproximada 6 veces mayor. La fórmula exacta se calcula teniendo en cuenta que la variación de la energía potencial del salto en la Tierra es igual que la variación de la energía potencial del salto en la Luna. La ecuación que se obtiene es la siguiente:
-GG*mm*MT/(RT+hT)+GG*mm*MT/RT = -GG*mm*ML/(RL+hL)+GG*mm*ML/RL
Despejando hL en función de hT se obtiene lo siguiente:
hL = MT*RL^2*hT/(-MT*RL*hT+ML*RT^2+ML*RT*hT)
hL = -(47000.*(-1.28839000*10^8+1.540375613*10^11*hT)) / (-1.199343378*10^15+3.975163007*10^9*hT)
Así, un salto vertical de 0.4 metros de altura en la Tierra en la Luna tendría 2.41 metros.
Un salto vertical de 0.2 metros de altura en la Tierra en la Luna tendría 1.2 metros.
Un salto vertical de 0.1 metros de altura en la Tierra en la Luna tendría 0.6 metros.
Un salto vertical de 0.05 metros de altura en la Tierra en la Luna tendría 0.3 metros.
Viendo los videos de la misión Apollo 11 con los astronautas en la Luna se observa que los astronautas saltan muy poco. Casi nunca saltan más de 10 o 20 centímetros. ¿Por qué?
Según Nasa:
"24. A spacesuit weighs approximately 280 pounds on the ground -- without the astronaut in it. In the microgravity environment of space, a spacesuit weighs nothing."
Fuente: http://www.nasa.gov/audience/foreducators/spacesuits/facts/facts-index.html
280 pounds son unos 127 kg (127.006 kg).