¿Calcular el tiempo de vuelo entre dos anomalías da un resultado negativo?

Quiero calcular el tiempo de vuelo (TOF) entre dos anomalías verdaderas en una órbita con elementos orbitales conocidos. Basado en este artículo de wikipedia y la publicación Tiempo de vuelo entre dos anomalías en una trayectoria orbital conocida , terminé usando la siguiente fórmula.
Conociendo las verdaderas anomalías v 1 y v 2 (con v 2 > v 1 ) en algunos momentos (desconocidos) t 1 y t 2 ( t 2 > t 1 ), expreso el TOF t 2 t 1 de v 1 a v 2 por:

t 2 t 1 = a 3 m ( mi 2 mi 1 + mi ( pecado mi 1 pecado mi 2 ) )
Dónde mi 1 y mi 2 son las anomalías excéntricas correspondientes calculadas usando la ecuación:
mi = 2 arcán ( 1 mi 1 + mi broncearse v 2 )

Sin embargo, esta expresión en la práctica me da un TOF negativo incorrecto.

Esto se puede verificar con un ejemplo numérico simple. Desde a 3 / m es siempre positivo, podemos centrarnos sólo en el signo de la pieza que implica las anomalías excéntricas. Considerando :

mi = 0.11 (órbita ligeramente elíptica)
v 1 = 2.43   radical
v 2 = 3.78   radical

Calculamos mi 1 y mi 2 :

mi 1 2.35   radical
mi 2 2.43   radical

Que luego da:

( mi 2 mi 1 + mi ( pecado mi 1 pecado mi 2 ) ) = ( 2.43 2.35 + 0.11 × ( pecado ( 2.35 ) pecado ( 2.43 ) ) ) 4.63 < 0

Este problema surgirá siempre que una anomalía verdadera sea mayor que π debido a la π -periodicidad de broncearse . ¿Hay alguna forma de contrarrestar esto? ¿Existe una fórmula para TOF que funcione bien para todos los casos? Si es posible, uno que también funcione con una órbita hiperbólica (a través de cambios de signo y el uso de la anomalía hiperbólica).

Gracias de antemano.

@Uwe No entiendo tu comentario. Nunca mencioné dos o más cuerpos. Me estoy enfocando en una única órbita elíptica perfecta (estoy trabajando en una simulación). Una nave espacial se mueve en esa órbita. Quiero calcular el tiempo de vuelo de esa nave espacial que va desde un punto A en esa órbita especificada por su verdadera anomalía. v 1 , y un punto B de anomalía verdadera v 2 .
Como está trabajando con una órbita elíptica, continúe trabajando en el cálculo. Si el tiempo de vuelo resultante es negativo, agregue el período orbital.
@novotny Alternativamente agregue 2 π a la anomalía excéntrica final ( mi 2 ) cuando mi 2 < mi 1 .
Lo siento, confundí anomalías y mascons.

Respuestas (1)

Cuando la anomalía verdadera está entre 0 y π la fórmula para E está bien, pero arriba π tiene que ajustar el valor de Arctan para obtener E en el cuadrante correcto. Parece que solo agrega 2 π a un valor negativo E.

Sí, hay un sistema para no usar E, llamado variables universales. Ver Fundamentos de astrodinámica por Bate et al. , y Fundamentos de Astrodinámica y Aplicaciones de Vallado, ambos en el Capítulo 4.

Utilice asteriscos para títulos de libros y artículos. El uso de signos de dólar es para matemáticas, y aunque el modo matemático es maravilloso para matemáticas, es pésimo para poner palabras en cursiva. Te arreglaré esto.