Articulo de referencia

Mecánica orbital

Un satélite que orbita la Tierra tiene una velocidad tangencial y una aceleración hacia el interior . La mecánica orbital , la astrodinámica o la dinámica espacial consisten en ...

Un satélite que orbita la Tierra tiene una velocidad tangencial y una aceleración hacia el interior .

La mecánica orbital , la astrodinámica o la dinámica espacial consisten en la aplicación de la balística y la mecánica celeste a cohetes , satélites y otras naves espaciales . El movimiento de estos objetos se calcula generalmente a partir de las leyes del movimiento y de la gravitación universal formuladas por Isaac Newton . La astrodinámica es una disciplina fundamental en el diseño y control de misiones espaciales .

La mecánica celeste trata de forma más amplia la dinámica orbital de los sistemas bajo la influencia de la gravedad , incluyendo tanto naves espaciales como cuerpos astronómicos naturales como sistemas estelares , planetas , lunas y cometas . La mecánica orbital se centra en las trayectorias de las naves espaciales , incluyendo maniobras orbitales , cambios de plano orbital y transferencias interplanetarias, y es utilizada por los planificadores de misiones para predecir los resultados de las maniobras de propulsión .

La relatividad general es una teoría más exacta que las leyes de Newton para calcular órbitas, y a veces es necesario utilizarla para obtener mayor precisión o en situaciones de alta gravedad (por ejemplo, órbitas cerca del Sol).

Historia

Hasta el auge de los viajes espaciales en el siglo XX, apenas existía distinción entre mecánica orbital y celeste. En la época del Sputnik , el campo se denominaba dinámica espacial . [ 1 ] Por lo tanto , las técnicas fundamentales, como las utilizadas para resolver el problema kepleriano (determinar la posición en función del tiempo), son las mismas en ambos campos. Además, la historia de ambos campos es prácticamente compartida.

Johannes Kepler fue el primero en modelar con éxito las órbitas planetarias con un alto grado de precisión, publicando sus leyes en 1609. Isaac Newton publicó leyes más generales del movimiento celeste en la primera edición de Philosophiæ Naturalis Principia Mathematica (1687), que proporcionó un método para encontrar la órbita de un cuerpo que sigue una trayectoria parabólica a partir de tres observaciones. [ 2 ] Este método fue utilizado por Edmund Halley para establecer las órbitas de varios cometas, incluido el que lleva su nombre . El método de aproximación sucesiva de Newton fue formalizado en un método analítico por Leonhard Euler en 1744, cuyo trabajo fue a su vez generalizado a órbitas elípticas e hiperbólicas por Johann Lambert entre 1761 y 1777.

Otro hito en la determinación de órbitas fue la contribución de Carl Friedrich Gauss al descubrimiento del planeta enano Ceres en 1801. El método de Gauss permitía, con tan solo tres observaciones (pares de ascensión recta y declinación ), hallar los seis elementos orbitales que describen completamente una órbita. Posteriormente, la teoría de la determinación de órbitas se desarrolló hasta el punto de aplicarse hoy en día en receptores GPS, así como en el seguimiento y la catalogación de planetas menores recién descubiertos . La determinación y predicción de órbitas modernas se utilizan para operar todo tipo de satélites y sondas espaciales, ya que es necesario conocer sus posiciones futuras con gran precisión.

La astrodinámica fue desarrollada por el astrónomo Samuel Herrick a partir de la década de 1930. Consultó con el científico espacial Robert Goddard , quien lo animó a continuar su trabajo en técnicas de navegación espacial, ya que Goddard creía que serían necesarias en el futuro. Las técnicas numéricas de la astrodinámica se combinaron con nuevas y potentes computadoras en la década de 1960, y los humanos estaban listos para viajar a la Luna y regresar.

Técnicas prácticas

Reglas generales

Las siguientes reglas generales son útiles para situaciones aproximadas por la mecánica clásica bajo las suposiciones estándar de la astrodinámica que se describen a continuación. El ejemplo específico que se analiza es el de un satélite orbitando un planeta, pero estas reglas también podrían aplicarse a otras situaciones, como las órbitas de cuerpos pequeños alrededor de una estrella como el Sol.

  • Leyes de Kepler sobre el movimiento planetario :
    • Las órbitas son elípticas , con el cuerpo más masivo en uno de los focos de la elipse. Un caso especial de esto es una órbita circular (un círculo es un caso especial de elipse) con el planeta en el centro.
    • Una línea trazada desde el planeta hasta el satélite barre áreas iguales en tiempos iguales, independientemente de la porción de la órbita que se mida.
    • El cuadrado del período orbital de un satélite es proporcional al cubo de su distancia media al planeta.
  • Sin aplicar ninguna fuerza (como encender un motor de cohete), el período y la forma de la órbita del satélite no cambiarán.
  • Un satélite en una órbita baja (o en la parte baja de una órbita elíptica) se mueve más rápidamente con respecto a la superficie del planeta que un satélite en una órbita más alta (o en la parte alta de una órbita elíptica), debido a la mayor atracción gravitatoria que se produce más cerca del planeta.
  • Si se aplica propulsión en un solo punto de la órbita del satélite, este regresará a ese mismo punto en cada órbita subsiguiente, aunque el resto de su trayectoria cambie. Por lo tanto, no es posible pasar de una órbita circular a otra con una sola aplicación breve de propulsión.
  • Si se aplica un empuje en dirección opuesta al movimiento del satélite, la órbita se transforma en una elíptica. El satélite descenderá hasta alcanzar el punto orbital más bajo (el periapsis ) a 180 grados del punto de disparo, para luego ascender de nuevo. El periodo de la órbita resultante será menor que el de la órbita circular original. Si se aplica un empuje en la dirección del movimiento del satélite, se crea una órbita elíptica cuyo punto más alto ( apoapsis ) se encuentra a 180 grados del punto de disparo. El periodo de la órbita resultante será mayor que el de la órbita circular original.

Las consecuencias de las leyes de la mecánica orbital a veces resultan contraintuitivas. Por ejemplo, si dos naves espaciales se encuentran en la misma órbita circular y desean acoplarse, la nave que va detrás no puede simplemente encender sus motores para acelerar hacia la que va delante. Esto alteraría la forma de su órbita, provocando que gane altitud y disminuya su velocidad con respecto a la nave que va delante, alejándose así del objetivo. El encuentro espacial previo al acoplamiento normalmente requiere múltiples encendidos de motores calculados con precisión en varios periodos orbitales, lo que puede tardar horas o incluso días en completarse.

En la medida en que no se cumplen las suposiciones estándar de la astrodinámica, las trayectorias reales diferirán de las calculadas. Por ejemplo, la simple resistencia atmosférica es otro factor que complica la situación para los objetos en órbita terrestre baja .

Estas reglas generales resultan claramente imprecisas al describir dos o más cuerpos de masa similar, como un sistema estelar binario (véase el problema de los n cuerpos ). La mecánica celeste utiliza reglas más generales aplicables a una mayor variedad de situaciones. Las leyes de Kepler sobre el movimiento planetario, que pueden derivarse matemáticamente de las leyes de Newton, se cumplen estrictamente solo al describir el movimiento de dos cuerpos gravitatorios en ausencia de fuerzas no gravitatorias; también describen trayectorias parabólicas e hiperbólicas. En las proximidades de objetos grandes como las estrellas, las diferencias entre la mecánica clásica y la relatividad general también cobran importancia.

Leyes de la astrodinámica

Las leyes fundamentales de la astrodinámica son la ley de gravitación universal de Newton y las leyes del movimiento de Newton , mientras que la herramienta matemática fundamental es el cálculo diferencial .

En un marco newtoniano, las leyes que rigen las órbitas y las trayectorias son, en principio , simétricas en el tiempo .

En astrodinámica, las suposiciones estándar incluyen la ausencia de interferencia de cuerpos externos, una masa despreciable para uno de los cuerpos y fuerzas adicionales despreciables (como las del viento solar , la resistencia atmosférica, etc.). Si bien se pueden realizar cálculos más precisos sin estas simplificaciones, estos son más complejos. La mayor precisión a menudo no supone una diferencia significativa en el cálculo como para que valga la pena.

Las leyes de Kepler sobre el movimiento planetario pueden derivarse de las leyes de Newton, suponiendo que el cuerpo en órbita está sujeto únicamente a la fuerza gravitatoria del atractor central. Cuando existe un empuje o fuerza propulsora, las leyes de Newton siguen siendo válidas, pero las leyes de Kepler quedan invalidadas. Al cesar el empuje, la órbita resultante será diferente, pero volverá a estar descrita por las leyes de Kepler, expuestas anteriormente. Las tres leyes son:

  1. La órbita de cada planeta es una elipse con el Sol en uno de sus focos .
  2. Una línea que une un planeta con el Sol barre áreas iguales durante intervalos de tiempo iguales.
  3. Los cuadrados de los períodos orbitales de los planetas son directamente proporcionales a los cubos del semieje mayor de las órbitas.

Velocidad de escape

La fórmula para la velocidad de escape se deriva de la siguiente manera. La energía específica (energía por unidad de masa ) de cualquier vehículo espacial se compone de dos componentes: la energía potencial específica y la energía cinética específica . La energía potencial específica asociada a un planeta de masa M viene dada por ϵpag=GRAMOMETROr,{\displaystyle \epsilon _{p}=-{\frac {GM}{r}},} donde G es la constante gravitacional y r es la distancia entre los dos cuerpos. La energía cinética específica de un objeto viene dada por ϵk=v22,{\displaystyle \epsilon _{k}={\frac {v^{2}}{2}},} donde v es su velocidad. Por lo tanto, la energía orbital específica total es ϵ=ϵk+ϵpag=v22GRAMOMETROr.{\displaystyle \epsilon =\epsilon _{k}+\epsilon _{p}={\frac {v^{2}}{2}}-{\frac {GM}{r}}.}

Dado que la energía se conserva ,ϵ{\displaystyle \epsilon }no puede depender de la distancia,r{\displaystyle r}, desde el centro del cuerpo central hasta el vehículo espacial en cuestión, es decir, v debe variar con r para mantener constante la energía orbital específica. Por lo tanto, el objeto puede alcanzar infinitor{\displaystyle r}solo si esta cantidad es no negativa, lo que implica v2GRAMOMETROr.{\displaystyle v\geq {\sqrt {\frac {2GM}{r}}}.}

La velocidad de escape desde la superficie terrestre es de aproximadamente 11  km/s, pero resulta insuficiente para enviar el cuerpo a una distancia infinita debido a la atracción gravitatoria del Sol. Para escapar del Sistema Solar desde una ubicación a una distancia del Sol igual a la distancia Sol-Tierra, pero no cercana a la Tierra, se requiere  una velocidad de alrededor de 42 km/s. Sin embargo, se tendrá en cuenta la velocidad orbital de la Tierra para las naves espaciales lanzadas desde la Tierra, si su aceleración adicional (debido al sistema de propulsión) las lleva en la misma dirección que la Tierra en su órbita.

Fórmulas para órbitas libres

Las órbitas son secciones cónicas , por lo que la fórmula para la distancia de un cuerpo para un ángulo dado corresponde a la fórmula para esa curva en coordenadas polares , que es r=pag1+miporqueθ,μ=GRAMO(metro1+metro2),pag=h2/μ,{\displaystyle {\begin{aligned}r&={\frac {p}{1+e\cos \theta }},\\\mu &=G(m_{1}+m_{2}),\\p&=h^{2}/\mu ,\end{aligned}}} dóndeμ{\displaystyle \mu }se denomina parámetro gravitacional ,metro1{\displaystyle m_{1}}ymetro2{\displaystyle m_{2}}son las masas de los objetos 1 y 2, yh{\displaystyle h}es el momento angular específico del objeto 2 con respecto al objeto 1. El parámetroθ{\displaystyle \theta }se conoce como la verdadera anomalía ,pag{\displaystyle p}es el semi-lato recto , mientras quemi{\displaystyle e}es la excentricidad orbital , todo obtenible a partir de las diversas formas de los seis elementos orbitales independientes .

órbitas circulares

Todas las órbitas acotadas donde predomina la gravedad de un cuerpo central son elípticas. Un caso especial es la órbita circular, que es una elipse de excentricidad cero. La fórmula para la velocidad de un cuerpo en una órbita circular a una distancia r del centro de gravedad de masa M se puede derivar de la siguiente manera:

La aceleración centrífuga coincide con la aceleración debida a la gravedad. Por lo tanto, v2r=GRAMOMETROr2.{\displaystyle {\frac {v^{2}}{r}}={\frac {GM}{r^{2}}}.} Por lo tanto, v=GRAMOMETROr,{\displaystyle v={\sqrt {\frac {GM}{r}}},} donde G = 6,6743 × 10−11 m3/(kg·s2) es laconstante gravitacional.

Para utilizar correctamente esta fórmula, las unidades deben ser consistentes; por ejemplo,METRO{\displaystyle M}debe estar en kilogramos, yr{\displaystyle r}Si la medida debe estar en metros, la respuesta estará en metros por segundo.

La cantidadGRAMOMETRO{\displaystyle GM}A menudo se le denomina parámetro gravitacional estándar , que tiene un valor diferente para cada planeta o luna del Sistema Solar .

Una vez que se conoce la velocidad orbital circular, la velocidad de escape se encuentra fácilmente multiplicando por2{\displaystyle {\sqrt {2}}}: v=2GRAMOMETROr=2GRAMOMETROr.{\displaystyle v={\sqrt {2}}{\sqrt {\frac {GM}{r}}}={\sqrt {\frac {2GM}{r}}}.}

Para escapar de la gravedad, la energía cinética debe al menos igualar la energía potencial negativa. Por lo tanto, 12metrov2=GRAMOMETROmetror,v=2GRAMOMETROr.{\displaystyle {\begin{aligned}&{\frac {1}{2}}mv^{2}={\frac {GMm}{r}},\\&v={\sqrt {\frac {2GM}{r}}}.\end{aligned}}}

Órbitas elípticas

Si0<mi<1{\displaystyle 0<e<1}, entonces el denominador de la ecuación de órbitas libres varía con la anomalía verdaderaθ{\displaystyle \theta }, pero permanece positivo, sin llegar a ser nunca cero. Por lo tanto, el vector de posición relativa permanece acotado, teniendo su magnitud más pequeña en el periapsis.rpag{\displaystyle r_{p}}, que viene dado por rpag=pag1+mi.{\displaystyle r_{p}={\frac {p}{1+e}}.}

El valor máximor{\displaystyle r}se alcanza cuandoθ=180{\displaystyle \theta =180^{\circ }}. Este punto se llama apoapsis, y su coordenada radial se denotara{\displaystyle r_{a}}, es ra=pag1mi.{\displaystyle r_{a}={\frac {p}{1-e}}.}

Dejar2a{\displaystyle 2a}sea ​​la distancia medida a lo largo de la línea del ábside desde el periapsisPAG{\displaystyle P}a la apoapsisA{\displaystyle A}: 2a=rpag+ra.{\displaystyle 2a=r_{p}+r_{a}.} Sustituyendo las ecuaciones anteriores, obtenemos a=pag1mi2,{\displaystyle a={\frac {p}{1-e^{2}}},} donde a es el semieje mayor de la elipse. Resolviendo parapag{\displaystyle p}y sustituyendo el resultado en la fórmula de la curva de sección cónica anterior, obtenemos r=a(1mi2)1+miporqueθ.{\displaystyle r={\frac {a(1-e^{2})}{1+e\cos \theta }}.}

Período orbital

Bajo supuestos estándar, el período orbital (T{\displaystyle T}) de un cuerpo que viaja a lo largo de una órbita elíptica se puede calcular como T=2πa3μ,{\displaystyle T=2\pi {\sqrt {\frac {a^{3}}{\mu }}},} dóndeμ{\displaystyle \mu }es el parámetro gravitacional estándar ya{\displaystyle a}es la longitud del semieje mayor . Conclusiones:

Velocidad

Bajo supuestos estándar, la velocidad orbital (v{\displaystyle v}) de un cuerpo que viaja a lo largo de una órbita elíptica se puede calcular a partir de la ecuación vis-viva como v=μ(2r1a),{\displaystyle v={\sqrt {\mu \left({\frac {2}{r}}-{\frac {1}{a}}\right)}},} dóndeμ{\displaystyle \mu }es el parámetro gravitacional estándar ,r{\displaystyle r}es la distancia entre los cuerpos en órbita,a{\displaystyle a}es la longitud del semieje mayor .

La ecuación de velocidad para una trayectoria hiperbólica es v=μ(2r+1|a|).{\displaystyle v={\sqrt {\mu \left({\frac {2}{r}}+{\frac {1}{|a|}}\right)}}.}

Energía

Bajo supuestos estándar, energía orbital específica (ϵ{\displaystyle \epsilon }) de la órbita elíptica es negativa, y la ecuación de conservación de la energía orbital (la ecuación vis viva ) para esta órbita puede tomar la forma v22μr=μ2a=ϵ<0,{\displaystyle {\frac {v^{2}}{2}}-{\frac {\mu }{r}}=-{\frac {\mu }{2a}}=\epsilon <0,} dóndev{\displaystyle v}es la velocidad del cuerpo en órbita,r{\displaystyle r}es la distancia del cuerpo en órbita desde el centro de masa del cuerpo central ,a{\displaystyle a}es el semieje mayor ,μ{\displaystyle \mu }es el parámetro gravitacional estándar . Conclusiones:

  • Para un semieje mayor dado, la energía orbital específica es independiente de la excentricidad.

Utilizando el teorema del virial, encontramos:

  • El promedio temporal de la energía potencial específica es igual a2ϵ{\displaystyle 2\epsilon },
  • el promedio temporal der1{\displaystyle r^{-1}}esa1{\displaystyle a^{-1}},
  • El promedio temporal de la energía cinética específica es igual aϵ{\displaystyle -\epsilon }.

órbitas parabólicas

Si la excentricidad es igual a 1, entonces la ecuación de la órbita se convierte en: r=h2μ11+porqueθ,{\displaystyle r={\frac {h^{2}}{\mu }}{\frac {1}{1+\cos \theta }},} dónder{\displaystyle r}es la distancia radial del cuerpo en órbita desde el centro de masa del cuerpo central ,h{\displaystyle h}es el momento angular específico del cuerpo en órbita ,θ{\displaystyle \theta }es la verdadera anomalía del cuerpo en órbita,μ{\displaystyle \mu }es el parámetro gravitacional estándar .

A medida que la anomalía verdadera θ se aproxima a 180°, el denominador se aproxima a cero, de modo que r tiende hacia el infinito. Por lo tanto, la energía de la trayectoria para la cual e  =  1 es cero, y viene dada por ϵ=v22μr=0,{\displaystyle \epsilon ={\frac {v^{2}}{2}}-{\frac {\mu }{r}}=0,} dóndev{\displaystyle v}es la velocidad del cuerpo en órbita.

En otras palabras, la velocidad en cualquier punto de una trayectoria parabólica es v=2μr.{\displaystyle v={\sqrt {\frac {2\mu }{r}}}.}

órbitas hiperbólicas

Simi>1{\displaystyle e>1}la fórmula de la órbita r=h2μ11+miporqueθ{\displaystyle r={\frac {h^{2}}{\mu }}{\frac {1}{1+e\cos \theta }}} describe la geometría de la órbita hiperbólica. El sistema consta de dos curvas simétricas. El cuerpo en órbita ocupa una de ellas; la otra es su imagen matemática vacía. Claramente, el denominador de la ecuación anterior tiende a cero cuandoporqueθ=1/mi{\displaystyle \cos \theta =-1/e}Denotamos este valor de anomalía verdadera como θ=porque1(1mi),{\displaystyle \theta _{\infty }=\cos ^{-1}\left(-{\frac {1}{e}}\right),} ya que la distancia radial se aproxima al infinito a medida que la anomalía verdadera se aproximaθ{\displaystyle \theta _{\infty }}, conocida como la verdadera anomalía de la asíntota . Obsérvese queθ{\displaystyle \theta _{\infty }}se encuentra entre 90° y 180°. A partir de la identidad trigonométricapecado2θ+porque2θ=1{\displaystyle \sin ^{2}\theta +\cos ^{2}\theta =1}resulta que pecadoθ=1mimi21.{\displaystyle \sin \theta _{\infty }={\frac {1}{e}}{\sqrt {e^{2}-1}}.}

Energía

Bajo supuestos estándar, energía orbital específica (ϵ{\displaystyle \epsilon }) de una trayectoria hiperbólica es mayor que cero, y la ecuación de conservación de la energía orbital para este tipo de trayectoria toma la forma ϵ=v22μr=μ2a,{\displaystyle \epsilon ={\frac {v^{2}}{2}}-{\frac {\mu }{r}}={\frac {\mu }{-2a}},} dóndev{\displaystyle v}es la velocidad orbital del cuerpo en órbita,r{\displaystyle r}es la distancia radial del cuerpo en órbita desde el cuerpo central ,a{\displaystyle a}es el semieje mayor negativo de la hipérbola de la órbita ,μ{\displaystyle \mu }es un parámetro gravitacional estándar .

velocidad de exceso hiperbólica

Bajo supuestos estándar, el cuerpo que viaja a lo largo de una trayectoria hiperbólica alcanzará una velocidad de...r={\displaystyle r=\infty }una velocidad orbital llamada velocidad hiperbólica de exceso (v{\displaystyle v_{\infty }}) que se puede calcular como v=μa,{\displaystyle v_{\infty }={\sqrt {\frac {\mu }{-a}}},} dóndeμ{\displaystyle \mu }es parámetro gravitacional estándar ,a{\displaystyle a}es el semieje mayor negativo de la hipérbola de la órbita .

La velocidad de exceso hiperbólica está relacionada con la energía orbital específica o energía característica mediante 2ϵ=do3=v2.{\displaystyle 2\epsilon =C_{3}=v_{\infty }^{2}.}

Cálculo de trayectorias

La ecuación de Kepler

Un método para calcular órbitas (utilizado principalmente desde el punto de vista histórico) consiste en usar la ecuación de Kepler : METRO=miϵpecadomi,{\displaystyle M=E-\epsilon \sin E,} donde M es la anomalía media , E es la anomalía excéntrica yϵ{\displaystyle \epsilon }es la excentricidad .

Con la fórmula de Kepler, encontrar el tiempo de vuelo para alcanzar un ángulo ( anomalía verdadera ) deθ{\displaystyle \theta }El periapsis se divide en dos etapas:

  1. Calcular la anomalía excéntricami{\displaystyle E}de verdadera anomalíaθ{\displaystyle \theta }.
  2. Calcula el tiempo de vuelot{\displaystyle t}de la anomalía excéntricami{\displaystyle E}.

Encontrar la anomalía excéntrica en un momento dado ( el problema inverso ) es más difícil. La ecuación de Kepler es trascendental enmi{\displaystyle E}, lo que significa que no se puede resolvermi{\displaystyle E}algebraicamente . La ecuación de Kepler se puede resolver parami{\displaystyle E}analíticamente por inversión.

Una solución de la ecuación de Kepler, válida para todos los valores reales deϵ{\displaystyle \epsilon }es mi={norte=1METROnorte3norte¡límiteθ0(dnorte1dθnorte1[(θθpecadoθ3)norte]),ϵ=1,norte=1METROnortenorte¡límiteθ0(dnorte1dθnorte1[(θθϵpecadoθ)norte]),ϵ1.{\displaystyle E={\begin{cases}\displaystyle \sum _{n=1}^{\infty }{\frac {M^{\frac {n}{3}}}{n!}}\lim _{\theta \to 0}\left({\frac {\mathrm {d} ^{n-1}}{\mathrm {d} \theta ^{n-1}}}\left[\left({\frac {\theta }{\sqrt[{3}]{\theta -\sin \theta }}}\right)^{n}\right]\right),&\epsilon =1,\\\displaystyle \sum _{n=1}^{\infty }{\frac {M^{n}}{n!}}\lim _{\theta \to 0}\left({\frac {\mathrm {d} ^{n-1}}{\mathrm {d} \theta ^{n-1}}}\left[\left({\frac {\theta }{\theta -\epsilon \sin \theta }}\right)^{n}\right]\right),&\epsilon \neq 1.\end{cases}}}

Evaluar esto produce mi={incógnita+160incógnita3+11400incógnita5+125200incógnita7+4317248000incógnita9+12137207200000incógnita11+15143912713500800000incógnita13incógnita=(6METRO)13,ϵ=111ϵMETROϵ(1ϵ)4METRO33¡+(9ϵ2+ϵ)(1ϵ)7METRO55¡(225ϵ3+54ϵ2+ϵ)(1ϵ)10METRO77¡+(11025ϵ4+4131ϵ3+243ϵ2+ϵ)(1ϵ)13METRO99¡,ϵ1.{\displaystyle E={\begin{cases}\displaystyle x+{\frac {1}{60}}x^{3}+{\frac {1}{1400}}x^{5}+{\frac {1}{25200}}x^{7}+{\frac {43}{17248000}}x^{9}+{\frac {1213}{7207200000}}x^{11}+{\frac {151439}{12713500800000}}x^{13}\cdots \mid x=(6M)^{\frac {1}{3}},&\epsilon =1\\[1ex]\displaystyle {\frac {1}{1-\epsilon }}M-{\frac {\epsilon }{(1-\epsilon )^{4}}}{\frac {M^{3}}{3!}}+{\frac {(9\epsilon ^{2}+\epsilon )}{(1-\epsilon )^{7}}}{\frac {M^{5}}{5!}}-{\frac {(225\epsilon ^{3}+54\epsilon ^{2}+\epsilon )}{(1-\epsilon )^{10}}}{\frac {M^{7}}{7!}}+{\frac {(11025\epsilon ^{4}+4131\epsilon ^{3}+243\epsilon ^{2}+\epsilon )}{(1-\epsilon )^{13}}}{\frac {M^{9}}{9!}}\cdots ,&\epsilon \neq 1.\end{cases}}}

Alternativamente, la ecuación de Kepler se puede resolver numéricamente. Primero, hay que adivinar un valor demi{\displaystyle E}y calcular el tiempo de vuelo; luego ajustarmi{\displaystyle E}Según sea necesario, se ajusta el tiempo de vuelo calculado para que se aproxime al valor deseado hasta alcanzar la precisión requerida. Generalmente, se utiliza el método de Newton para lograr una convergencia relativamente rápida.

La principal dificultad de este enfoque es que puede tardar prohibitivamente mucho tiempo en converger para las órbitas elípticas extremas. Para órbitas casi parabólicas, la excentricidadϵ{\displaystyle \epsilon }es casi 1, y sustituyendomi=1{\displaystyle e=1}en la fórmula para la anomalía media,mipecadomi{\displaystyle E-\sin E}Nos encontramos restando dos valores casi iguales, lo que perjudica la precisión . Para órbitas casi circulares, es difícil hallar el periapsis (y las órbitas verdaderamente circulares carecen por completo de periapsis). Además, la ecuación se derivó bajo el supuesto de una órbita elíptica, por lo que no es válida para órbitas parabólicas o hiperbólicas. Estas dificultades fueron las que llevaron al desarrollo de la formulación de variables universales , descrita a continuación.

Órbitas cónicas

Para procedimientos sencillos, como el cálculo delta-v para elipses de transferencia coplanares, los métodos tradicionales son bastante eficaces. Otros, como el de tiempo de vuelo, son mucho más complejos, especialmente para órbitas casi circulares e hiperbólicas.

La aproximación cónica parcheada

La órbita de transferencia de Hohmann por sí sola es una aproximación deficiente para trayectorias interplanetarias porque ignora la gravedad propia de los planetas. La gravedad planetaria domina el comportamiento de la nave espacial en las proximidades de un planeta y, en la mayoría de los casos, Hohmann sobreestima gravemente delta-v y produce prescripciones muy imprecisas para los tiempos de encendido. Una forma relativamente sencilla de obtener una aproximación de primer orden de delta-v se basa en la técnica de la "aproximación cónica parcheada". Se debe elegir el cuerpo gravitatorio dominante en cada región del espacio por la que pasará la trayectoria y modelar solo los efectos de ese cuerpo en esa región. Por ejemplo, en una trayectoria de la Tierra a Marte, se comenzaría considerando solo la gravedad de la Tierra hasta que la trayectoria alcance una distancia donde la gravedad de la Tierra ya no domine la del Sol. A la nave espacial se le daría una velocidad de escape para enviarla al espacio interplanetario. Luego, se consideraría solo la gravedad del Sol hasta que la trayectoria alcance las proximidades de Marte. Durante esta etapa, el modelo de órbita de transferencia es apropiado. Finalmente, solo se considera la gravedad de Marte durante la parte final de la trayectoria, donde la gravedad de Marte domina el comportamiento de la nave espacial. La nave se aproximaría a Marte en una órbita hiperbólica, y una última maniobra retrógrada la desaceleraría lo suficiente como para ser capturada por Marte. Friedrich Zander fue uno de los primeros en aplicar el enfoque de cónicas parcheadas con fines astrodinámicos, al proponer el uso de la gravedad de cuerpos intermedios para viajes interplanetarios, en lo que hoy se conoce como asistencia gravitatoria . [ 3 ]

El tamaño de los "vecindarios" (o esferas de influencia ) varía con el radio. rSOI=apag(metropagmetroS)2/5.{\displaystyle r_{\text{SOI}}=a_{\text{p}}\left({\frac {m_{\text{p}}}{m_{\text{S}}}}\right)^{2/5}.} dóndeapag{\displaystyle a_{\text{p}}}es el semieje mayor de la órbita del planeta en relación con el Sol ;metropag{\displaystyle m_{\text{p}}}ymetroS{\displaystyle m_{\text{S}}}son las masas del planeta y del Sol, respectivamente.

Esta simplificación es suficiente para calcular estimaciones aproximadas de los requisitos de combustible y del tiempo de vuelo, pero generalmente no es lo suficientemente precisa como para guiar una nave espacial a su destino. Para ello, se requieren métodos numéricos.

La formulación de la variable universal

Para superar las limitaciones computacionales de los métodos tradicionales para resolver el problema de dos cuerpos, se desarrolló la formulación de variables universales . Esta funciona igualmente bien para los casos circular, elíptico, parabólico e hiperbólico, y las ecuaciones diferenciales convergen satisfactoriamente al integrarse para cualquier órbita. Además, se generaliza bien a problemas que incorporan la teoría de perturbaciones.

Perturbaciones

La formulación de variables universales funciona bien con la técnica de variación de parámetros, excepto que ahora, en lugar de los seis elementos orbitales keplerianos, utilizamos un conjunto diferente de elementos orbitales: concretamente, los vectores de posición y velocidad iniciales del satélite.incógnita0{\displaystyle x_{0}}yv0{\displaystyle v_{0}}en una época determinadat=0{\displaystyle t=0}En una simulación de dos cuerpos, estos elementos son suficientes para calcular la posición y la velocidad del satélite en cualquier momento futuro, utilizando la formulación de variables universales. A la inversa, en cualquier instante de la órbita del satélite, podemos medir su posición y velocidad, y luego usar el método de variables universales para determinar cuál habría sido su posición y velocidad iniciales en ese momento. En un movimiento perfecto de dos cuerpos, estos elementos orbitales serían invariantes (al igual que los elementos keplerianos).

Sin embargo, las perturbaciones hacen que los elementos orbitales cambien con el tiempo. Por lo tanto, el elemento de posición se escribe comoincógnita0(t){\displaystyle x_{0}(t)}y el elemento de velocidad comov0(t){\displaystyle v_{0}(t)}, lo que indica que varían con el tiempo. La técnica para calcular el efecto de las perturbaciones se convierte en una de encontrar expresiones, ya sean exactas o aproximadas, para las funciones.incógnita0(t){\displaystyle x_{0}(t)}yv0(t){\displaystyle v_{0}(t)}.

A continuación se describen algunos efectos que hacen que las órbitas reales difieran de los modelos simples basados ​​en una Tierra esférica. La mayoría de ellos pueden ser tratados en escalas de tiempo cortas (quizás menos de unos pocos miles de órbitas) mediante la teoría de perturbaciones, ya que son pequeños en relación con los efectos de dos cuerpos correspondientes.

  • Los abultamientos ecuatoriales provocan la precesión del nodo y del perigeo.
  • Los armónicos teserales [ 4 ] del campo gravitatorio introducen perturbaciones adicionales.
  • Las perturbaciones gravitatorias lunares y solares alteran las órbitas.
  • La resistencia atmosférica reduce el semieje mayor a menos que se utilice empuje de compensación.

En escalas de tiempo muy largas (quizás millones de órbitas), incluso pequeñas perturbaciones pueden predominar y el comportamiento puede volverse caótico . Por otro lado, los astrónomos expertos pueden orquestar las diversas perturbaciones para ayudar en las tareas de mantenimiento de órbita, como el mantenimiento de la posición , el mantenimiento o ajuste de la trayectoria terrestre o la sincronización del perigeo para cubrir objetivos seleccionados a baja altitud.

Maniobra orbital

En los vuelos espaciales , una maniobra orbital consiste en el uso de sistemas de propulsión para cambiar la órbita de una nave espacial . Para las naves espaciales que se encuentran lejos de la Tierra —por ejemplo, las que orbitan alrededor del Sol—, una maniobra orbital se denomina maniobra en el espacio profundo (DSM, por sus siglas en inglés) .

Transferencia orbital

Las órbitas de transferencia suelen ser órbitas elípticas que permiten a las naves espaciales pasar de una órbita (generalmente circular) a otra. Normalmente, requieren un encendido al inicio, otro al final y, a veces, uno o más encendidos en el medio.

  • La órbita de transferencia de Hohmann requiere un delta-v mínimo .
  • Una transferencia bielíptica puede requerir menos energía que la transferencia de Hohmann, si la relación de órbitas es de 11,94 o mayor, [ 5 ] pero conlleva un mayor tiempo de viaje que la transferencia de Hohmann.
  • Las transferencias más rápidas pueden utilizar cualquier órbita que se cruce con las órbitas de origen y destino, a costa de un mayor delta-v.
  • Utilizando motores de bajo empuje (como la propulsión eléctrica ), si la órbita inicial es supersíncrona con la órbita circular final deseada, la órbita de transferencia óptima se logra mediante un empuje continuo en la dirección de la velocidad en el apogeo. Sin embargo, este método requiere mucho más tiempo debido al bajo empuje. [ 6 ]

En el caso de transferencia orbital entre órbitas no coplanares, el empuje de cambio de plano debe realizarse en el punto donde se intersecan los planos orbitales (el "nodo"). Dado que el objetivo es cambiar la dirección del vector velocidad en un ángulo igual al ángulo entre los planos, casi todo este empuje debe realizarse cuando la nave espacial se encuentra en el nodo cerca del apoapsis, cuando la magnitud del vector velocidad es mínima. Sin embargo, una pequeña fracción del cambio de inclinación orbital puede realizarse en el nodo cerca del periapsis, inclinando ligeramente el empuje de inyección de la órbita de transferencia en la dirección del cambio de inclinación deseado. Esto funciona porque el coseno de un ángulo pequeño es muy cercano a uno, lo que hace que el pequeño cambio de plano sea prácticamente "gratuito" a pesar de la alta velocidad de la nave espacial cerca del periapsis, ya que el efecto Oberth debido al aumento del empuje ligeramente inclinado supera el costo del empuje en el eje normal a la órbita.

Una transferencia de Hohmann desde una órbita circular baja a una órbita circular más alta.
Una transferencia bielíptica desde una órbita circular baja inicial (azul oscuro) a una órbita circular más alta (rojo).
Transferencia elíptica genérica de dos impulsos entre dos órbitas circulares
Una transferencia general de una órbita circular baja a una órbita circular más alta.
Secuencia óptima para transferir un satélite de una órbita supersíncrona a una órbita geosíncrona mediante propulsión eléctrica.

Asistencia gravitatoria y el efecto Oberth

En una maniobra de asistencia gravitatoria , una nave espacial pasa cerca de un planeta y se aleja en una dirección y velocidad diferentes. Esto resulta útil para acelerar o desacelerar la nave en lugar de transportar más combustible.

Esta maniobra puede aproximarse a una colisión elástica a grandes distancias, aunque el sobrevuelo no implica contacto físico. Según la tercera ley de Newton (acción y reacción), cualquier momento que gane una nave espacial debe ser perdido por el planeta, o viceversa. Sin embargo, dado que el planeta es mucho más masivo que la nave, el efecto sobre su órbita es insignificante.

El efecto Oberth puede emplearse, en particular durante una operación de asistencia gravitatoria. Este efecto consiste en que el uso de un sistema de propulsión funciona mejor a altas velocidades, por lo que los cambios de rumbo se realizan mejor cerca de un cuerpo gravitatorio; esto puede multiplicar el delta-v efectivo .

Red de transporte interplanetario y órbitas difusas

Las rutas orbitales se pueden identificar mediante interacciones gravitacionales no lineales entre planetas y lunas del Sistema Solar. Estas trayectorias forman parte de la Red de Transporte Interplanetario y siguen órbitas caóticas que requieren un mínimo de combustible más allá de alcanzar un punto de Lagrange , con correcciones de rumbo periódicas. Las órbitas se pueden trazar desde la órbita terrestre alta hasta Marte a través de puntos troyanos terrestres . Las limitaciones incluyen largos tiempos de tránsito, que a menudo duran varios años, y ventanas de lanzamiento poco frecuentes.

Sin embargo, se han utilizado en proyectos como Génesis . Esta nave espacial visitó el punto L1 Tierra-Sol y regresó utilizando muy poco combustible.

Véase también

Referencias

  1. Thomson, William T. (1961). Introducción a la dinámica espacial . Nueva York: Wiley.
  2. Bate, RR; Mueller, DD; White, JE (1971). Fundamentos de astrodinámica . Courier Corporation. pág. 5. ISBN  978-0-486-60061-1.
  3. Negri, Rodolfo Batista; Prado, Antônio Fernando Bertachini de Alme (agosto de 2020). "Una revisión histórica de la teoría de las asistencias gravitatorias en la era pre-espacial" . Revista de la Sociedad Brasileña de Ciencias Mecánicas e Ingeniería . 42 (8): 406. doi : 10.1007/s40430-020-02489-x . S2CID 220510617 . 
  4. Weisstein, Eric W. "Armónico teseral" . MathWorld . Consultado el 7 de octubre de 2019 .
  5. Vallado, David Anthony (2001). Fundamentos de astrodinámica y aplicaciones . Springer. pág. 317. ISBN  0-7923-6903-3.
  6. Spitzer, Arnon (1997). Trayectoria óptima de órbita de transferencia mediante propulsión eléctrica . USPTO.

Lecturas adicionales

  • Lynnane George. Introducción a la mecánica orbital .
  • Sellers, Jerry J.; Astore, William J.; Giffen, Robert B.; Larson, Wiley J. (2004). Kirkpatrick, Douglas H. (ed.). Comprender el espacio: Una introducción a la astronáutica (2.ª  ed.). McGraw Hill. pág.  228. ISBN 0-07-242468-0.
  • "Introducción al espacio de la Universidad del Aire, Capítulo 8 - Mecánica orbital" (PDF) . USAF. Archivado del original (PDF) el 14 de febrero de 2013. Consultado el 13 de octubre de 2007 .
  • Bate, RR; Mueller, DD; White, JE (1971). Fundamentos de astrodinámica . Dover Publications, Nueva York. ISBN 978-0-486-60061-1.
  • Vallado, DA (2001). Fundamentos de astrodinámica y aplicaciones (2.ª  ed.). Springer. ISBN 978-0-7923-6903-5.
  • Battin, RH (1999). Introducción a las matemáticas y métodos de la astrodinámica . American Institute of Aeronautics & Astrodynamics, Washington, DC ISBN 978-1-56347-342-5.
  • Chobotov, VA, ed. (2002). Mecánica orbital (3.ª  ed.). American Institute of Aeronautics & Astronautics, Washington, DC ISBN 978-1-56347-537-5.
  • Herrick, S. (1971). Astrodinámica: Determinación de órbitas, navegación espacial, mecánica celeste, volumen 1. Van Nostrand Reinhold, Londres. ISBN 978-0-442-03370-5.
  • Herrick, S. (1972). Astrodinámica: Corrección orbital, teoría de perturbaciones, integración, volumen 2. Van Nostrand Reinhold, Londres. ISBN 978-0-442-03371-2.
  • Kaplan, MH (1976). Dinámica y control de naves espaciales modernas . Wiley, Nueva York. ISBN 978-0-471-45703-9.
  • Tom Logsdon (1997). Mecánica orbital . Wiley-Interscience, Nueva York. ISBN 978-0-471-14636-0.
  • John E. Prussing y Bruce A. Conway (1993). Mecánica orbital . Oxford University Press, Nueva York. ISBN 978-0-19-507834-3.
  • MJ Sidi (2000). Dinámica y control de naves espaciales . Cambridge University Press, Nueva York. ISBN 978-0-521-78780-2.
  • WE Wiesel (1996). Dinámica del vuelo espacial (2.ª  ed.). McGraw-Hill, Nueva York. ISBN 978-0-07-070110-6.
  • JP Vinti (1998). Mecánica orbital y celeste . American Institute of Aeronautics & Ast, Reston, Virginia. ISBN 978-1-56347-256-5.
  • P. Gurfil (2006). Astrodinámica moderna . Butterworth-Heinemann. ISBN 978-0-12-373562-1.
  • MECÁNICA ORBITAL (Tecnología de cohetes y espacio)
  • Kit de herramientas de astrodinámica Java
  • Grafo de conocimiento de eventos y tráfico espacial basado en astrodinámica