Uso de los sistemas de coordenadas en la propagación en órbita

Para la propagación de la órbita en el sistema de coordenadas ECI, considerando los efectos de los armónicos zonales y las perturbaciones Sol/Luna, se formó el siguiente procedimiento:

  1. Convertir el vector de estado a ECEF
  2. Calcular la aceleración debida a los armónicos
  3. Convertir la aceleración a ECI
  4. Calcular la aceleración del Sol/Luna e integrar

Según tengo entendido , la conversión a ECEF se realiza porque la orientación de la Tierra afecta la aceleración por armónicos. La integración en ECI se realiza para evitar las fuerzas ficticias.

Preguntas

  1. ¿Es obligatorio convertir a ECEF? Me encantaría obtener una explicación detallada sobre cómo la orientación afecta la aceleración.
  2. ¿Cómo convertir la aceleración a un sistema de coordenadas diferente? Puedo convertir las coordenadas pero no el cambio de velocidad.

  • ECI: inercial centrada en la tierra (marco no giratorio)
  • ECEF: centrado en la Tierra, fijo en la Tierra (marco giratorio)
¡Gran pregunta! Acabo de agregar el acrónimo/la jerga para beneficiar a los lectores que no lo hayan visto antes.

Respuestas (1)

Tengamos en cuenta A γ un vector de magnitud | γ | = γ definido en el cuadro A . Dejar A y B ser dos marcos distintos, | A γ | = | B γ | . Además, puede ser útil recordar la definición de marco inercial: es un marco cuya aceleración es nula.

Finalmente, recordemos el teorema del transporte, según "Mecánica analítica de los sistemas espaciales" de Schaub y Junkins (esta es una cita exacta):

Dejar norte y B ser dos marcos con un vector de velocidad angular relativa de ω B / norte , y deja r ser un vector genérico; entonces la derivada de r en el norte frame se puede relacionar con la derivada de r en el B marco como:

norte d d t r = B d d t r + ω B / norte × r

Una cita posterior de gran interés es la siguiente:

... encontraremos que los vectores se diferencian típicamente con respecto a un marco inercial llamado norte .

Como se explica en la sección 2.4 de "GPS" de G. Xu e Y. Xu, 2016:

El movimiento de los satélites sigue la mecánica newtoniana, [y estos] solo son válidos y se expresan en un sistema de coordenadas inercial.

En el caso de la astrodinámica, las palabras "marco" y "sistema de coordenadas" casi siempre son intercambiables.

Ahora podemos responder a sus preguntas:

  1. ¿Es obligatorio convertir a ECEF? Me encantaría obtener una explicación detallada sobre cómo la orientación afecta la aceleración. La propagación del movimiento de la nave espacial debido a las fuerzas gravitatorias externas debe calcularse en un marco inercial. Si el movimiento de esa nave espacial se debe principalmente a la gravedad de la Tierra, entonces el movimiento celeste debe calcularse en el marco ECI. Sin embargo, es obligatorio convertir la posición y la velocidad del satélite al marco ECEF para calcular cómo los armónicos esféricos de la Tierra afectan el movimiento del satélite. Para una definición detallada del marco ECEF, recomiendo el Capítulo 2 de "GPS" de G. Xu y Y. Xu, 2016 (creo que este capítulo está disponible gratuitamente en el sitio web de Springer). Para una mayor explicación de por qué esto es importante, tenga en cuenta que si un marco es inercial pero no el otro, entonces la magnitud de un vector dado en un momento dado puede ser diferente, pero las componentes de dicho vector no lo serán. Por ejemplo, imagina dos marcos A y B que, en el momento t k están orientados exactamente de tal manera que la conversión de A a B corresponde a una rotación de + π / 2 sobre el eje Z. Además, deja A α = [ 1 ,   0 ,   0 ] T . Entonces, en el B marco, tenemos B α = [ 0 ,   1 ,   0 ] T . Por lo tanto, si una aceleración en el A marco conduce a un cambio de la X componente de 1 en un momento posterior de t norte , A α t norte = [ 0 ,   0 ,   0 ] T . Pero sería un error aplicar ese cambio desde el marco. A directamente a B α sin transformar ese vector en el A marco primero. De hecho, si lo hiciera, encontraría que su (incorrecto) actualizado B α t norte = [ 1 ,   1 ,   0 ] T . Eso es evidentemente incorrecto ya que | A α t norte | | B α t norte | .

  2. ¿Cómo convertir la aceleración a un sistema de coordenadas diferente? Puedo convertir las coordenadas pero no el cambio de velocidad. En realidad, no convierte la aceleración a un marco diferente. Convierta la posición y la velocidad al marco ECEF, calcule la aceleración debida a los armónicos en ECEF y luego reconstruya la posición y la velocidad actualizadas en el marco ECEF. Finalmente, vuelve a convertir su estado ECEF actualizado en ECI. Como prueba, así es como se hace en GMAT 2016a . (Estoy enlazando a Github porque el código fuente es significativamente más fácil de navegar allí).

Espero que esto ayude.

¡Gran respuesta! Pregunta: La posición y la velocidad actualizadas se reconstruyen a través de la integración. ¿Quieres decir que tengo que integrarme en ECEF?
@Leeloo, gracias, me alegro de que haya ayudado. No, no se integre en ECEF, integre en ECI. Pero, en cada paso de la integración, haga una copia de su estado y páselo a dos funciones diferentes. El primero calcula la gravedad newtoniana y el segundo convierte a ECEF, calcula los efectos de los armónicos, modifica el estado que se proporcionó, vuelve a convertir a ECI y devuelve ese nuevo estado de ECI. Luego, de vuelta en su integrador, sume los componentes de la dinámica newtoniana y la dinámica armónica.
@Leeloo, no, debes calcular el efecto de los armónicos con el estado anterior a los efectos newtonianos. Piénselo de esta manera: si usa el ECI posterior a la gravedad para calcular los armónicos, entonces está calculando los armónicos del estado futuro.
No, hice lo que dijiste, envié el estado inicial a 2 funciones diferentes. El punto es que la función de armónicos devuelve el nuevo estado ECI (no aceleración). no se como resumirlo con la aceleracion newtoniana
Ah ok ya veo. Una vez que tenga la aceleración newtoniana, puede agregarla al estado a partir del cual calculó la aceleración: X t k + 1 = X t k + X ˙ t k . Haz lo mismo con los armónicos (usando las ecuaciones de Pines u otro método). Y luego sume el estado actualizado de la gravedad newtoniana y el de los armónicos (asegúrese de que tanto la posición como la velocidad estén en ECI).
¡Excelente! Gracias. Comparé los resultados con GMAT y, curiosamente, para Xy Zcoordenadas obtuve 1mun error, pero para Y- 7m. El Yerror de coordenadas siempre es mucho mayor que para otros..
@Leeloo ¡esa es una muy buena precisión! Actualmente estoy trabajando para obtener esa buena precisión. ¿Cómo estás haciendo la transformación de coordenadas de ECI a ECEF? Si no me equivoco, GMAT usa un vector de rotación general por día que tiene en cuenta el movimiento polar, la nutación y la precesión, pero no coincide con el último marco de la IAU.
Configuré el año 2100 como inicial e integrado durante 10 días. Estoy escribiendo sobre Julia. Para el cálculo de la matriz rotacional utilicé SatelliteToolbox.jl, tiene en cuenta la nutación, la precesión y el movimiento polar. Usé IAU más antiguo
Interesante. Pensé que la matriz GMAT solo subía hasta 2100, pero debo estar equivocado.