Un software para la propagación de la órbita, considerando la Luna, el Sol y los armónicos zonales.

Tengo que propagar la órbita de un satélite en 1 año, cumpliendo con el Sol, la Luna y los efectos armónicos zonales.

¿Hay preferiblemente un software libre, que podría ser utilizado? Me sugirieron usar PyEphem . ¿Realmente da resultados precisos, considerando la Luna, el Sol, los armónicos zonales y, probablemente, Júpiter y Venus?

Lo que realmente necesito es calcular el cambio en el argumento del perigeo de la órbita en 1 año.

Respuestas (2)

Recomiendo encarecidamente el GMAT (Herramienta de análisis de misión general) de NASA Goddard. Es bastante fácil de usar, tiene varios tutoriales y se ha utilizado en operaciones de naves espaciales.

Pero, ¿cómo considerar el Sol, la Luna, los armónicos zonales allí? ¿Es posible?
Además, aquí se describe el PyEphem, que puede considerar esas perturbaciones. ¿Qué piensas? rhodesmill.org/pyephem/quick.html#bodies-with-orbital-elements
Sí, GMAT le permite considerar todos los objetos celestes disponibles en cualquier archivo DE*.bsp, incluidas lunas específicas de Júpiter si es necesario. Por lo general, no necesita más que el Sol, la Luna, la presión de la radiación solar y armónicos zonales de 21 por 21, al menos para órbitas terrestres bajas (en varios de los análisis de alta fidelidad que hago en el trabajo, solo usamos ese tipo de modelo ).
PyEphem parece proporcionar solo efemérides planetarias, y supongo que necesitarías codificar tu propia propagación de órbita. GMAT hace todo eso y está completamente validado. No dudaría ni un segundo en usar este software de la NASA, especialmente porque es gratuito y no tiene restricciones ITAR.
Aquí he mostrado mi configuración del modelo de fuerza GMAT. ¿Podrías, por favor, verificarlo? space.stackexchange.com/q/28327/25120

Descargo de responsabilidad completo: soy el autor y principal desarrollador de poliastro.

En poliastro hay varias perturbaciones ya definidas, y entre ellas tienes los armónicos zonales, la presión solar y el efecto gravitatorio de la Luna. Puedes ver la lista completa de perturbaciones aquí:

http://docs.poliastro.space/en/latest/api.html#module-poliastro.twobody.perturbations

Y un ejemplo simple de cómo usarlos aquí:

http://docs.poliastro.space/en/latest/user_guide.html#studying-non-keplerian-orbits-perturbations

Este es un trabajo en progreso y es posible que encuentre algunos problemas en el camino: específicamente, que el rendimiento no es muy bueno y que extraer la historia temporal de los elementos keplerianos es mucho más difícil que hacerlo con los elementos cartesianos. ¡ Siéntete libre de unirte al chat y hacer cualquier tipo de pregunta!

Gracias. ¿Resuelve numéricamente (ecuaciones diferenciales)? ¿Cómo se explican los armónicos zonales? ¿Podría, por favor, mostrar fórmulas/enlaces?
Usamos el método de Cowell ( github.com/poliastro/poliastro/blob/master/src/poliastro/… ) y las fórmulas analíticas para los primeros armónicos ( github.com/poliastro/poliastro/blob/master/src/poliastro/… )