Formulario básico de mecánica orbital (generalidades)

Fuerza de atracción entre dos cuerpos MmG

Fg =
r2
Fuerza centrífuga MV 2
Fc =
r
Energía específica  V 2

= −
2 r

 =−
2a

Momento específico H = rV cos 

Relación excentricidad-Momento H2
e = 1−
a
Velocidades para diferentes órbitas 2 
V= −
r a

V= circular a = r
r
2 
V= − elíptica a  0
r a
2
V= parabolica a = 
r
2 
V= + hipérbolia a  0
r a

Formulario básico de mecánica orbital (órbitas elípticas)

Excentricidad e c
e=
a
ra − rp
e=
ra + rp
ra
e= −1
a
rp
e = 1−
a
r2 − r1
e=
r1 cos 1 − r2 cos  2
Ángulo de “Flight path”  e sin 
tan  =
1 + e cos
Velocidad media “n” 
n=
a3
Periodo “P” 2
P=
n
P = 2 a 3 
Radio “r”
r=
(
a 1− e2 )
1 + e cos
rp (1 + e )
r=
1 + e cos
Radio del apoapsis “ra” ra = a(1 + e )
ra = 2a − rp

ra = rp
(1 + e)
(1 − e)
Radio del periapsis “rp” rp = a(1 − e )
rp = 2a − ra

rp = ra
(1 − e )
(1 + e )
Semieje mayor “a”
a=
(ra + rp )
2
r
a=
2 − V 2 r
rp
a=
(1 − e)
ra
a=
(1 + e)
Tiempo desde el periapsis “t” t = (E − e sin E ) n
Anomalía excéntrica elíptica E e + cos
cos E =
1 + e cos
Anomalía verdadera  rp (1 + e ) 1
cos = −
re e
a (1 + e ) 1
2
cos = −
re e
Velocidad V 2 
V= −
r a
rpV p = raVa
Formulario básico de mecánica orbital (órbitas hiperbólicas)

Ángulo de asíntota  tan  = b a

tan  = bVHE
2

2brp
tan  =
b 2 − rp2
cos  = 1 e
Excentricidad e e= 1
cos 
e = 1 + (rp a )
e = 1+ b2 a2 ( )

Ángulo de “Flight path”  e sin 

tan  =
1 + e cos
Velocidad media “n” 
n=
a3
Radio “r”
r=
(
a e2 −1 )
1 + e cos
Radio del periapsis “rp” rp = b (e − 1) (e + 1)
rp = a(e − 1)
rp = c − a

rp = b tan  
 2
2  +  (e − 1)
rp =
V p2
2

   
rp = − +   + b2
V 2
V 2 
HE
 HE 
rp = − a + a 2 + b 2
Semieje mayor “a” a=b e2 −1
rp
a=
e −1

a= 2
VHE
(
a = b 2 − rp2 2rp )
rp
a=
rpV p2 − 2 
Semi ejemenor “b” b = rp (e + 1) (e − 1)
b = a e2 −1
2
b = rp 2
+1
rpVHE
Tiempo desde el periapsis “t” t = (e sinh F − F ) n
Anomalía excéntrica hiperbólica F e + cos
cosh F =
1 + e cos
(
F = ln cosh F + cosh 2 F − 1 )
sinh F = exp(F ) − exp(− F )
1
2
Anomalía verdadera 
cos =
(
a e2 −1 1
−
)
re e
Anomalía verdadera de la asíntota a  a = 180º  
1
cos a = −
e
V = VHE = 
Velocidad V
a

V =  2   + VHE
2

 r
V= (2 r ) +  a
C 3 = VHE
2

Formulario básico de mecánica orbital (Maniobras orbitales)

Concepto general Vcircular + V = Velíptica

Maniobra complanar (ecuación general) V 2 = Vi 2 + V f2 − 2ViV f cos 

Transferencia de Hohmann rpt = ri

rat = r f
V1 = V pt − Vi
V2 = Vat − V f
Cambio de plano orbital V = 2Vi sin ( 2 )
Vi = V f
Ángulo de cambio de plano orbital a cos  = cos ii cos i f + sin ii sin i f cos( )
Argumento de la latitud para el cambio de sin i f sin ( )
plano ALa sin ALa =
sin 
ALa = i +  i
Sitio de lanzamiento, campo visual de un satélite y groud track (traza de un satélite)

Máxima inclinación obtenible en función cos i = cos La  sin Az

de la latitud y azimut del sitio de
lanzamiento
Semiángulo de visión de un satélite Rs
sin  h =
h
Semiángulo de área cubierta por el satélite Rs
cos  h =
h
Distancia del satélite al horizonte Dh = r cos  h
Arco cubierto por el satélite S w = 2 h Rs
Latitud de un satélite La sin La = sin i sin ALa
Longitud de un satélite Lo Lo =  + lo
Longitud de un satélite lo medida desde el sin lo = tan La tan i
nodo ascendente
Ajuste de la longitud por rotación de la Re = 360t = 0.0041781t
tierra Re 86164
Ajuste de la longitud por regresión de los 
nodos
Ajuste global de la longitud Orbita directa Lo =  + lo −  − Re
hemisferio este
Ajuste global de la longitud Orbita directa Lo =  + lo +  + Re
hemisferio oeste
Ajuste global de la longitud Orbita Lo =  + lo +  − Re
retrograda hemisferio este
Ajuste global de la longitud Orbita Lo =  + lo −  + Re
retrograda hemisferio oeste

Perturbaciones orbitales

Variación del nodo  cos i

 SOL = −0,00154
ascendente en una órbita n
debido a la presencia del sol
Variación del nodo  cos i
 LUNA = −0,00338
ascendente en una órbita n
debido a la presencia de la
luna
Variación del argumento del 4 − 5  sen 2i
perigeo en una órbita debido  SOL = 0,00077
n
a la presencia del sol
Variación del argumento del 4 − 5  sen 2i
perigeo en una órbita debido  LUNA = 0,00169
n
a la presencia del sol
  = −1,5  n  J   RE   (cos i )
Variación del nodo 2


( )
ascendente en una órbita J2 2
 a  1− e2
2
debido a la no esfericidad de
la tierra
 −2,06474 1014 
1

(cos i )
a
7
( ) (
1− e2
2
)
Variación del argumento del
 J
2
(
 R  4 − 5  sen i
= 0,75  n  J 2   E  
2
)
( )
perigeo en una órbita debido
 a 
2
1− e2
2

a la no esfericidad de la tierra
 −1,03237 10 14

1

(4 − 5  sen i ) 2

( a ) (1 − e )
7 2 2

CD  A 2
 a   p  exp(− c )  I 0 + 2  e  I1 
Variación del semieje mayor
arev = −2 
en una orbita debido a la m
resistencia aerodinámica
Variación de la excentricidad CD  A 2  
 a   p  exp(− c )   I 0 +  (I 0 + I 2 )
e
en una orbita debido a la erev = −2 
m  2 
resistencia aerodinámica
Variación del semieje mayor CD  A 2
arev = −2  a 
en una orbita circular debido m
a la resistencia aerodinámica
Variación de la excentricidad erev = 0
en una orbita circular debido
a la resistencia aerodinámica
Variación del periodo en una CD  A 2
orbita circular debido a la − 6  2 
a 
Prev = m
resistencia aerodinámica V
Variación de la velocidad en C A
Vrev =   D  a   V
una orbita circular debido a m
la resistencia aerodinámica
Desaceleración producida 1 C A 2
aD = −    D V
por efecto de la resistencia 2 m
aerodinámica
Coeficiente balístico m
(CD  A)
Desaceleración producida A
aR  −4,5  10−8 
por efecto de la resistencia m
aerodinámica

Diseño de orbitas

Orbita Geosincrona P = 23h56 min 4.09s = 86164.09s

Fases de ejecución 1. Orbita Circular baja de aparcamiento.
2. Orbita de transferencia de Hohmann
3. Cambio de plano a plano ecuatorial.
4. Circularización de la orbita de Hohmann a la
orbita final.
Órbita heliosíncrona Permite obtener un ángulo constante entre el plano
orbital y la posición del sol, para lograrlo debo hacer
coincidir que mi regresión de nodos coincida a la
velociad de rotación de la tierra respecto al sol.
Órbita Molniya

Órbita Tundra Misma inclinación que la molniya, pero con un periodo

orbital de 23h y 56 min y 4.09s

Ventanas de lanzamiento

Orbita Geosincrona P = 23h56 min 4.09s = 86164.09s

Fases de ejecución 5. Orbita Circular baja de aparcamiento.
6. Orbita de transferencia de Hohmann
7. Cambio de plano a plano ecuatorial.
8. Circularización de la orbita de Hohmann a la
orbita final.
Órbita heliosíncrona Permite obtener un ángulo constante entre el plano
orbital y la posición del sol, para lograrlo debo hacer
coincidir que mi regresión de nodos coincida a la
velociad de rotación de la tierra respecto al sol.
Órbita Molniya

Orbitas interplanetarias

25
Radio de la  M planeta 
esfera de Rs  R 
influencia  M sol 
Periodo sinódico 2 1
S= =
tierra −  planeta 1 − 1 Pplaneta
Trayectoria de sin  sin b
sin it =
partida (cambio sin c
de plano orbital) a = ( + 180 ) − Ltierra
b = L planeta − ( + 180 )
cos c = cos a cos b + sin a sin b cos 
 = 180 − i planeta
Trayectoria de cos  = cos it cos( tierra +  vehículo )
partida
C 3 = VHE = Vtierra − sol + Vvehículo− sol − 2Vtierra − solVvehículo− sol cos 
2 2 2
(calculando VHE
y C3)
Trayectoria de sin  sin a
sin itp =
llegada (cambio sin c
de plano orbital)
Trayectoria de cos  a = cos itp cos( tierra −  vehículo )
llegada
V2 = Vtierra − sol + Vvehículo− sol − 2Vtierra − solVvehículo− sol cos  a
2 2
(calculando V∞)
Trayectorias V = 2V cos 
ayudadas por la
gravedad V = 2V
e
(cambio de V planeta − sol sin  a
velocidad) sin  =
Va
− salida − sol = Vvehículo −llegada − sol + V − Vvehículo −llegada − sol V cos(    )
2 2 2
Vvehículo

− sol + Vvehículo − salida − sol − V

2 2 2
V planeta
cos  d =
2V planeta − solVvehículo− salida − sol