Viaje de la Tierra a Marte siguiendo una trayectoria en forma de espiral logarítmica.

La espiral logarítmica es una de las curvas notables o maravillosas junto a la catenaria, la cicloide, etc. La ecuación de la espiral logarítmica en coordenadas polares es

r=r0·exp(b·θ)

Donde r0 es el radio inicial, b es un parámetro y θ es el ángulo en radianes.

La principal característica de la espiral logarítmica, es que el radio vector r y la tangente a la espiral forman un ángulo π/2-γ que se mantiene constante.

El vector velocidad de la partícula (en coordenadas polares) que describe una espiral logarítmica es

v = dr dt r ^ +r dθ dt θ ^ =r(b r ^ + θ ^ )

Calculamos el ángulo π/2-γ entre los vectores v y r , empleando la definición de producto escalar de dos vectores.

cos( π 2 γ )=sinγ= v r v·r = b b 2 +1 ,b=tanγ

Cuando b→0, φ→π/2 y rr0, la espiral logarítmica se convierte en una circunferencia de radio r0.

Longitud de un arco de la espiral logarítmica

La longitud de un arco infinitesimal ds es el módulo del vector desplazamiento dr

ds= (dr) 2 + (r·dθ) 2

La longitud del arco de espiral logarítmica entre θ=0 y θ es

0 s ds = 0 θ (dr) 2 + (r·dθ) 2 = 0 θ 1+ b 2 r 0 e bθ dθ= 1+ b 2 b r 0 ( e bθ 1 )= 1+ b 2 b (r r 0 )

Ecuaciones del movimiento

En todos los casos analizados en este capítulo, a partir de la expresión de la fuerza sobre el cuerpo celeste, obtenemos la ecuación de la trayectoria. En esta sección, conocida la ecuación de la trayectoria (espiral logarítmica), obtenemos la fuerza tangente a la trayectoria que actúa sobre le cuerpo celeste, además de la atracción inversamente proporcional al cudrado de la distancia al centro de fuerzas.

Supondremos que el motor de la nave espacial de masa m proporciona una fuerza de empuje F que tiene la misma dirección que la velocidad de la nave (tangente a la trayectoria)

Las fuerzas que actúan sobre la nave espacial son:

En coordenadas polares, la ecuación del movimiento se escribe

m( d 2 r d t 2 r ( dθ dt ) 2 )=FsinγG Mm r 2 m( r d 2 θ d t 2 +2 dr dt dθ dt )=Fcosγ

Como la ecuación de la trayectoria es r=r0·exp(θ·tanγ), calculamos la derivada primera y segunda de r con respecto del tiempo t

dr dt =tanγ·r dθ dt d 2 r d t 2 =tanγ·r d 2 θ d t 2 + tan 2 γ·r ( dθ dt ) 2

Sustituimos en la primera ecuación diferencial

tanγ·r d 2 θ d t 2 + tan 2 γ·r ( dθ dt ) 2 r ( dθ dt ) 2 = F m sinγG M r 2 r d 2 θ d t 2 = F m cosγG M r 2 1 tanγ + r tanγ ( dθ dt ) 2 tanγ·r ( dθ dt ) 2

Sustituimos en la segunda ecuación diferencial, rd2θ/dt2 y dr/dt llegando a

r ( dθ dt ) 2 =G M r 2 cos 2 γ

Ecuaciones paramétricas de la trayectoria

Sabiendo que la ecuación de la trayectoria r=r0·exp(θ·tanγ), integramos esta ecuación diferencial para obtener el ángulo θ en función del tiempo t.

0 θ exp( 3 2 θ·tanγ )dθ = GM cosγ r 0 3/2 0 t dt exp( 3 2 θ·tanγ )=1+ 3 2 sinγ GM r 0 3 t

Obtenemos r en función del tiempo t

dr dt =tanγ·r dθ dt = GM r 1/2 sinγ r 0 r r dr = GM sinγ 0 t dt r 3/2 = r 0 3/2 + 3 2 GM sinγ·t

Las ecuaciones paramétricas (en función del tiempo t) de la trayectoria son

r= ( r 0 3/2 + 3 2 GM sinγ·t ) 2/3 θ= 2 3tanγ ln( 1+ 3 2 sinγ GM r 0 3 t )

Velocidad de la nave espacial

El módulo de la velocidad de la nave espacial en coordenadas polares es

v= ( dr dt ) 2 + ( r dθ dt ) 2 = GM sin 2 γ r + r 2 GM cos 2 γ r 3 = GM r

es igual a la velocidad de una nave espacial que describa una órbita circular de radio r.

Fuerza de empuje tangente a la trayectoria

La segunda ecuación del movimiento, nos permite calcular el módulo de la fuerza F de empuje de los motores de la nave.

( r d 2 θ d t 2 +2r ( dθ dt ) 2 tanγ )= F m cosγ

Conocemos la expresión de la derivada primera dθ/dt y calculamos la derivada segunda  d2θ/dt2 del ángulo θ respecto del tiempo t. Para ello, combinamos las ecuaciones

r 3 ( dθ dt ) 2 =GM cos 2 γ

r=r0·exp(θ·tanγ),

para obtener la expresión

exp(3θ·tanγ) ( dθ dt ) 2 = GM cos 2 γ r 0 3

Derivando con respecto del tiempo y teniendo en cuenta que, el segundo miembro es constante.

d 2 θ d t 2 = 3tanγ 2 ( dθ dt ) 2 = 3GMsinγcosγ 2 r 3

Ahora, se despeja de la segunda ecuación del movimiento, la fuerza de empuje

F m = 1 2 GMsinγ r 2

El empuje F necesario para que la nave espacial describa una trayectoria en forma de espiral logarítmica va disminuyendo de forma inversamente proporcional al cuadrado de la distancia r al Sol. Su dirección es la misma que la velocidad (tangente a la trayectoria).

Trabajo-energía

El trabajo de las fuerzas que actúan sobre el cuerpo celeste es igual a la variación de energía cinética

La resultante de las fuerzas que actúan sobre la partícula de masa m es

F m =( FsinγG M r 2 ) r ^ +( Fcosγ ) θ ^ =( 1 2 GM sin 2 γ r 2 G M r 2 ) r ^ +( 1 2 GMsinγcosγ r 2 ) θ ^

El vector desplazamiento (en coordendas polares)

dr =dr· r ^ +rdθ· θ ^

Efectuamos el producto escalar e integramos para calcular el trabajo (por unidad de masa)

W= A B F m dr = A B { ( 1 2 GM sin 2 γ r 2 G M r 2 )dr+ 1 2 GMsinγcosγ r 2 r·dθ } = GM( 1 1 2 sin 2 γ ) r 0 r dr r 2 + 1 2 GMsinγcosγ 0 θ dθ r 0 exp( θ·tanγ ) = GM( 1 1 2 sin 2 γ )( 1 r 1 r 0 ) 1 2 GMsinγcosγ 1 tanγ ( 1 r 1 r 0 )= GM( 1 r 1 r 0 )( 1 1 2 sin 2 γ 1 2 cos 2 γ )= 1 2 GM( 1 r 1 r 0 )

La variación de energía cinética (por unidad de masa) es

Δ E k = 1 2 GM r 1 2 GM r 0 = 1 2 GM( 1 r 1 r 0 )

Tiempo de viaje

El tiempo de viaje T a un planeta que dista r1 del Sol partiendo de r0 es

T= 2( r 1 3/2 r 0 3/2 ) 3 GM sinγ

Dado el tiempo de viaje T entre dos posiciones r1 y r0 calculamos γ, tal como veremos en el ejemplo más abajo

Las ecuaciones paramétricas (en función del tiempo t) de la trayectoria son

r= r 0 { 1+( ( r 1 r 0 ) 3/2 1 ) t T } 2/3 θ= 2 3tanγ ln{ 1+( ( r 1 r 0 ) 3/2 1 ) t T }

Cuando llega a la posicíon final r=r1

θ 1 = 1 tanγ ln( r r 0 )

Ejemplo

Si queremos que el viaje de la nave espacial entre la Tierra y Marte se realice en T=3 años=36 meses=1080 días. De la fórmula del tiempo T de viaje, obtenemos el valor del ángulo γ=1.8172º.

El desplazamiento angular de la nave espacial es θ1=13.2819 rad=761º

M=1.98e30; %masa del Sol
r0=1.496e11; %radio de la órbita de la Tierra (1 UA)
G=6.67e-11; %constante
r1=1.524*r0; %Marte (1.524 UA)
T=1080*24*60*60; %tiempo de viaje 3 años
gamma=asin(2*(r1^(3/2)-r0^(3/2))/(3*T*sqrt(G*M)));
r=@(t) (1+((r1/r0)^(3/2)-1)*t/T).^(2/3);
th=@(t) 2*log(1+((r1/r0)^(3/2)-1)*t/T)/(3*tan(gamma));
fplot(@(t) r(t).*cos(th(t)),@(t) r(t).*sin(th(t)),[0,T])
grid on
axis equal
hold off
xlabel('x (UA)')
ylabel('y (UA)')
title('Trayectoria Tierra-Marte')

El ángulo γ es positivo, el empuje de la nave espacial F tiene el mismo sentido que la velocidad, que disminuye a medida que se incrementa la distancia radial. La posición angular final θ1 es

>> gamma*180/pi
ans =    1.8172
>> th(T)
ans =   13.2819

También es posible el viaje a los planetas interiores como Venus, cambiamos la línea de código r1=0.723*r0;, obtenemos la siguiente trayectoria

El ángulo γ es negativo, el empuje de la nave espacial F se opone a la velocidad, que se incrementa a medida que disminuye la distancia radial

>> gamma*180/pi
ans =   -0.7940

Viaje de la Tierra a Marte

Aplicando la ecuación de la dinámica del movimiento circular uniforme calculamos la velocidad angular de Marte en su órbita circular alrededor del Sol.

G Mm r M 2 =m ω M 2 r M

ωM=1.055·10-7 rad/s=0.5225º/día

Cálculos similares se realizan para obtener la velocidad angular de la Tierra ωT en su órbita circular alrededor del Sol.

Si la posición inicial (en el instante t=0) del planeta Marte es φ0M, su posición cuando llega la nave espacial a las proximidades de su órbita circular, al cabo de un tiempo T es

φM0MMT

La nave espacial parte de la Tierra en el instante t=0, cuando su posición es φ0T

Para llegar a la órbita de Marte, la nave se desplaza un ángulo θ1. Su posición angular final será φ0T+ θ1

Para que la posición de la nave espacial y de Marte coincidan se tiene que cumplir la igualdad

φ0T+ θ10MMT

La diferencia entre las posiciones angulares de la Tierra y de Marte en el momento del lanzamiento de la nave espacial (t=0) deberá ser de

φ0T - φ0MMT- θ1
φ0T - φ0M =
0.5225·1080-761=-196.6º

En la figura, se muestra las posiciones de la Tierra y Marte en el momento de partida de la nave espacial para un viaje de 3 años.

Actividades

Se introduce

Se pulsa el botón titulado Nuevo, para establecer la posición inicial angular de los planetas, dados por dos números aleatorios comprendidos entre 0 y 360. A continuación, se pulsa para que comiencen a moverse describiendo órbitas circulares.

Se calcula el intervalo angular (φ0T - φ0M) entre las posiciones de la Tierra y de Marte en el momento de la partida de la nave espacial, tal como se ha hecho en el ejemplo.

Se pulsa el botón pausa ||, para parar el movimiento y el botón paso a paso >|, para examinar las posiciones angulares de los planetas y acercarnos a la posición de lanzamiento.

Se pulsa el botón titulado Lanza.

Si la nave no llega a Marte, se pulsa el botón titulado Lanza y se vuelve a intentar

En la parte derecha, se nos proporciona los valores del tiempo (en días), la posición angular y la velocidad de los tres cuerpos: la nave espacial, la Tierra y Marte.

Observamos que la nave espacial parte de las proximidades de la Tierra con su misma velocidad orbital y llega a las proximidades de Marte con su misma velocidad orbital. 

Referencias

Bacon R. H. Logarithmic spiral: An ideal trajectory for interplanetary vehicle with engines of low sustained thrust. Am. J. Phys. 27 (1959), pp. 164-165.