Un cohete de empuje constante

Consideremos un cohete de masa inicial m que lleva una velocidad v respecto a un Sistema de Referencia Inercial (por ejemplo, la Tierra).

En el instante t+Δt, una masa Δμ se expulsa con una velocidad constante –u relativa al cohete, como consecuencia la masa restante (m-Δμ) del cohete se incrementa en v+Δv.

En el instante t el cohete de masa m lleva una velocidad v. Su momento lineal es

p(t)=mv

En el instante t+Δt

El momento lineal en este instante es

p(t+Δt)=(m-Δμ)(v+Δv)+ Δμ(–u+ v)

Si el cohete está en el espacio exterior, el sistema formado por el cohete y el combustible que expulsa es aislado, el momento lineal p permanece constante. La  ecuación del movimiento del cohete es

Δp= p(t+Δt)- p(t)=m·Δv- Δμ-ΔΔμ=0

En el límite cuando Δt→0

m dv dt u dμ dt =0

La masa M del sistema formado por el cohete m y el combustible expulsado μ es constante M=μ+m, por lo que dμ+dm=0. La masa del cohete disminuye en dm y aumenta la masa del combustible expulsado en la misma cantidad.

m dv dt +u dm dt =0

Despejando dv

dv=u dm m v 0 v dv =u m 0 m dm m

cuya integración entre los instantes 0 y t conduce a la siguiente expresión

v= v 0 +u·ln m 0 m v= v 0 +u·ln m 0 m 0 Dt

Suponemos que la cantidad de combustible quemado en la unidad de tiempo, D, es constante, D=-dm/dt. La masa m del cohete en el instante t valdrá m=m0-D·t. Donde m0 es la suma de la carga útil más el combustible inicial, y D·t es el combustible quemado al cabo de un cierto tiempo t.

Interpretamos la ecuación del movimiento, afirmando que el cohete se comporta como una partícula de masa variable m que se mueve bajo la acción de una fuerza de empuje constante u·D.

m dv dt =u·D

Para hallar el desplazamiento x del cohete en el tiempo t, se integra la expresión de la velocidad

x x 0 = 0 t v·dt

para lo cual, es útil conocer la integral

lnz·dz =z·lnzz

>> syms x;
>> int('log(x)',x)
ans =x*(log(x) - 1)

resultando la expresión

x= x 0 + v 0 t+utln m 0 + u D [ ( m 0 Dt )ln( m 0 Dt )+Dt m 0 ln m 0 ]

Cuando se agota el combustible, el cohete se mueve con velocidad constate.

Momento lineal

Un cohete con una masa inicial m0 empieza a expulsar los gases con velocidad u relativa al cohete y por consiguiente, comienza a moverse en el espacio exterior. Al cabo de un cierto tiempo, alcanza una velocidad v, expulsando los gases con una velocidad u relativa al cohete o con una velocidad v-u relativa al observador terrestre.

Si el cohete parte del reposo v0=0, la velocidad de los gases respecto del observador terrestre v-u no es constante. Al principio tiene sentido contrario al movimiento del cohete v-u<0, pero al cabo de un cierto tiempo t la velocidad de los gases v-u>0 es del mismo sentido que la del cohete.

u·ln m 0 m 0 Dt u>0 m 0 m 0 Dt >e

t>0.632 m 0 D

En el instante t=0.632m0/D los gases expulsados están en reposo respecto del observador terrestre.

El momento lineal del cohete en el instante t, es

mv=( m 0 Dt)·u·ln m 0 m 0 Dt

El momento lineal de los gases expulsados entre el instante t=0, y el instante t, es

0 t dm(vu) =D 0 t ( uln m 0 m 0 Dt u )dt =u( m 0 Dt)ln m 0 m 0 Dt

Como el cohete y los gases forman un sistema aislado el momento lineal del cohete es igual y de sentido contrario al de los gases expulsados desde el instante inicial t=0, al instante t. Como la velocidad de los gases respecto del observador terrestre no es constante es necesario calcular una integral para obtener el momento total de los gases expulsados hasta el instante t y comprobar de este modo, que se cumple el principio de conservación del momento lineal, principio en el que nos hemos basado por otra parte, para obtener la ecuación del movimiento del cohete.

Energías

La energía cinética del cohete en el instante t es

E c = 1 2 m v 2 = 1 2 ( m 0 Dt ) v 2

La energía cinética de los gases expulsados desde el instante t=0, al instante t es

E g = 0 t 1 2 dm ( vu ) 2 = 1 2 Dt· u 2 1 2 ( m 0 Dt ) ( uln m 0 m 0 Dt ) 2 E g = 1 2 Dt· u 2 1 2 ( m 0 Dt ) v 2

Para llegar a esta expresión, es necesario conocer el resultado de las siguientes las integrales:

lnz·dz =z·lnzz ln 2 z·dz =z· ln 2 z2zlnz+2z

>> syms x;
>> int('log(x)^2',x)
ans =x*(log(x)^2 - 2*log(x) + 2)

La energía cinética total del sistema aislado formado por el cohete y los gases en el instante t es

E c + E g = 1 2 Dt· u 2

El rendimiento del cohete en el instante t es el cociente entre la energía cinética del cohete y la energía cinética total (cohete más los gases expulsados).

Ejemplo

Masa total del cohete=carga útil+combustible+masa del recipiente

m0=800+9000+0.05·9000=10250 kg

  1. Tiempo que tarda en agotarse el combustible

    Como hay 9000 kg de combustible que se queman a razón de 1000 kg/s, el combustible se agota en 9 s.

  2. Velocidad máxima alcanzada por el cohete

  3. v=uln m 0 m 0 Dt =2000·ln( 10250 102509000 )=4208m/s

  4. Distancia recorrida hasta que se agota el combustible

  5. x=utln m 0 + u D [ ( m 0 Dt )ln( m 0 Dt )+Dt m 0 ln m 0 ] x=2000·9·ln10250+ 2000 1000 [ 1250ln1250+900010250ln10250 ]=12740m

  6. Energías

  7. La energía total proporcionada por el combustible es

    E i = 1 2 Dt· u 2 = 1 2 9000· 2000 2 J

    La energía cinética del cohete cuando se ha agotado el combustible es

    E k = 1 2 ( 102509000 )· 4208 2 J

    El rendimiento es el cociente entre Ek/Ei=61.5%

u=2000; %velocidad de escape de los gases (respecto al cohete)
combustible=9000; %combustible
carga=800;  %carga útil
m0=carga+1.05*combustible; %masa total
D=1000; %kg de combustible quemado por segundo
tf=combustible/D; %tiempo hasta que se agota el combustible

figure
t=0:0.05:tf;
v=u*log(m0./(m0-D*t));
plot(t,v)
grid on
xlabel('t(s)')
ylabel('v(m/s)')
title('Velocidad')

figure
x=u*t*log(m0)+u*((m0-D*t).*log(m0-D*t)+D*t-m0*log(m0))/D;
plot(t,x)
grid on
xlabel('t(s)')
ylabel('x(m)')
title('Posición')
fprintf('Velocidad final %4.1f, posición final %5.1f\n',v(end),x(end));
Velocidad final 4208.3, posición final 12739.7

Velocidad final 4208.3, posición final 12739.7

Actividades

Se introduce.

Se pulsa el botón titulado Nuevo.

El cohete parte con velocidad inicial cero v0=0 desde el origen x0=0.

La masa inicial m0 es la suma de la carga útil, más el combustible y más la masa del recipiente cilíndrico que será proporcional a la masa del combustible que contiene

masa inicial m0  =carga útil +(1+r) ·combustible.

donde r es del orden del 5% ó 0.05

El tiempo tmáx que tarda en agotarse el combustible es igual al cociente entre la masa de combustible c y la cantidad D que se quema por segundo

tmáx=c/D

Cuando se agota el combustible c, el cohete describe un movimiento rectilíneo y uniforme, ya que no actúan fuerzas sobre el mismo.

v=uln m 0 m 0 c

En la cola del cohete se dibuja una flecha que indica la fuerza de empuje. Como la velocidad de los gases es constante (en el Sistema de Referencia en el cohete) el empuje es constante.

Comprobamos que el cohete alcanza el mismo valor de la velocidad máxima, independientemente de la cantidad D de combustible quemado en la unidad de tiempo. Mantener constantes la cantidad de de combustible c, la carga útil y variar la cantidad de combustible quemado por segundo. Anotar en una tabla la velocidades finales vmáx, una vez agotado todo el combustible, el tiempo empleado en alcanzar la velocidad máxima tmáx, y el desplazamiento del cohete x hasta este instante.

Se muestra el balance energético del cohete. Un círculo muestra la energía inicial del combustible en el cohete y cómo esta energía se va transformando en energía cinética de los gases expulsados y en energía cinética del cohete. Al agotarse todo el combustible, una parte de la energía inicial del combustible se ha transformado en energía cinética del cohete.