În mecanica cerească , problema Lambert privește determinarea unei orbite pornind de la doi vectori de poziție și timpul de zbor. A fost rezolvat de matematicianul elvețian Johann Heinrich Lambert . Are aplicații importante în domeniul întâlnirilor și manevrelor orbitale .
Să presupunem că un corp, sub influența unei forțe gravitaționale, este observat pe măsură ce se deplasează de la un punct P 1 la un punct P 2 , de-a lungul unei traiectorii conice date, într-un timp T. În această ipoteză, teorema lui Lambert afirmă că timpul de zbor depinde de unii parametri geometrici ai problemei. Mai precis, enunțul teoremei este următorul:
„Timpul necesar unui corp pentru a se deplasa între două puncte ale unei traiectorii conice este doar o funcție a sumei distanțelor celor două puncte de la originea forței, a distanței liniare dintre cele două puncte și axa semi-majoră a conicii. " |
Exprimată într-un alt mod, problema Lambert poate fi formulată prin impunerea unor condiții limită asupra ecuației diferențiale a problemei cu doi corpuri
- {\ displaystyle {\ ddot {\ bar {r}}} = - \ mu \ cdot {\ frac {\ hat {r}} {r ^ {2}}}} ,
pentru care orbita Kepleriană este soluția generală.
Formulare exactă
Formularea precisă a problemei Lambert este următoarea.
Sunt date două momente diferite în timp {\ displaystyle t_ {1}} Și {\ displaystyle t_ {2}} și doi vectori de poziție {\ displaystyle {\ bar {r}} _ {1} = r_ {1} {\ hat {r}} _ {1}} Și {\ displaystyle {\ bar {r}} _ {2} = r_ {2} {\ hat {r}} _ {2}} .
Găsește soluția{\ displaystyle {\ bar {r}} (t)} care satisface ecuația diferențială de mai sus deci
- {\ displaystyle {\ bar {r}} (t_ {1}) = {\ bar {r}} _ {1}}
- {\ displaystyle {\ bar {r}} (t_ {2}) = {\ bar {r}} _ {2}} .
Analiza geometrică preliminară
Figura 1:
{\ displaystyle F_ {1}} , originea forței gravitaționale;
{\ displaystyle P_ {1}} , punctul corespunzător vectorului
{\ displaystyle {\ bar {r}} _ {1}} ;
{\ displaystyle P_ {1}} , punctul corespunzător vectorului
{\ displaystyle {\ bar {r}} _ {2}} .
Cele trei puncte
{\ displaystyle F_ {1}} , originea forței gravitaționale,
{\ displaystyle P_ {1}} , punct corespunzător vectorului {\ displaystyle {\ bar {r}} _ {1}} ,
{\ displaystyle P_ {1}} , punct corespunzător vectorului {\ displaystyle {\ bar {r}} _ {2}} ,
formează un triunghi în planul definit de vectori {\ displaystyle {\ bar {r}} _ {1}} Și {\ displaystyle {\ bar {r}} _ {2}} și reprezentată în Figura 1. Distanța dintre puncte {\ displaystyle P_ {1}} Și {\ displaystyle P_ {2}} Și {\ displaystyle 2d} , distanța dintre {\ displaystyle P_ {1}} Și {\ displaystyle F_ {1}} Și{\ displaystyle r_ {1} = r_ {m} -A} și distanța dintre {\ displaystyle P_ {2}} Și {\ displaystyle F_ {1}} Și{\ displaystyle r_ {2} = r_ {m} + A} . Valoarea a {\ displaystyle A} este pozitiv sau negativ în funcție de care dintre puncte {\ displaystyle P_ {1}} Și {\ displaystyle P_ {2}} este mai departe de {\ displaystyle F_ {1}} . Problema geometrică care trebuie rezolvată este găsirea unei elipse care traversează punctele {\ displaystyle P_ {1}} Și {\ displaystyle P_ {2}} și să ai un foc înăuntru {\ displaystyle F_ {1}} .
Figura 2: hiperbola cu focare în puncte
{\ displaystyle P_ {1}} Și
{\ displaystyle P_ {2}} și trecând prin punct
{\ displaystyle F_ {1}} .
Punctele {\ displaystyle F_ {1}} , {\ displaystyle P_ {1}} Și {\ displaystyle P_ {2}} definesc o hiperbolă care traversează punctul {\ displaystyle F_ {1}} și are focare în puncte {\ displaystyle P_ {1}} Și {\ displaystyle P_ {2}} , așa cum se arată în figura 2. Punctul {\ displaystyle F_ {1}} poate fi fie pe ramura dreaptă, fie pe cea stângă a hiperbolei, în funcție de semnul {\ displaystyle A} . Axa semi-majoră a hiperbolei este {\ displaystyle \ left | A \ right |} și excentricitate {\ displaystyle E} Și {\ displaystyle {\ frac {d} {\ left | A \ right |}}} .
Ecuația hiperbolei , în coordonate carteziene, este
- {\ displaystyle {\ frac {x ^ {2}} {A ^ {2}}} - {\ frac {y ^ {2}} {B ^ {2}}} = 1 \ quad (1)}
cu
- {\ displaystyle B = \ left | A \ right | {\ sqrt {E ^ {2} -1}} = {\ sqrt {d ^ {2} -A ^ {2}}} \ quad (2)}
Pentru orice punct de pe aceeași ramură a hiperbolei căreia îi aparține {\ displaystyle F_ {1}} diferența dintre distanțele până la puncte {\ displaystyle P_ {1}} Și {\ displaystyle P_ {2}} Și
- {\ displaystyle r_ {2} -r_ {1} = 2A \ quad (3)}
Pentru fiecare punct {\ displaystyle F_ {2}} pe cealaltă ramură a hiperbolei relația corespunzătoare este
- {\ displaystyle s_ {1} -s_ {2} = 2A \ quad (4)}
acesta este
- {\ displaystyle r_ {1} + s_ {1} = r_ {2} + s_ {2} \ quad (5)}
Figura 3: elipsa cu focare în puncte
{\ displaystyle F_ {1}} Și
{\ displaystyle F_ {2}} și trecând prin puncte
{\ displaystyle P_ {1}} Și
{\ displaystyle P_ {2}} .
Aceasta înseamnă punctele {\ displaystyle P_ {1}} Și {\ displaystyle P_ {2}} ambii aparțin elipsei care are punctele ca focare {\ displaystyle F_ {1}} Și {\ displaystyle F_ {2}} și axa semi-majoră
- {\ displaystyle a = {\ frac {r_ {1} + s_ {1}} {2}} = {\ frac {r_ {2} + s_ {2}} {2}} \ quad (6)}
Elipsa rezultată, a ales un punct {\ displaystyle F_ {2}} arbitrar, este prezentat în Figura 3.
Soluția problemei Lambert pentru o orbită eliptică
În primul rând, luăm în considerare două cazuri separate, în funcție de faptul dacă polul orbital are direcție {\ displaystyle {\ bar {r}} _ {1} \ times {\ bar {r}} _ {2}} sau {\ displaystyle - {\ bar {r}} _ {1} \ times {\ bar {r}} _ {2}} . În primul caz unghiul de transfer {\ displaystyle \ alpha} pentru prima trecere prin {\ displaystyle {\ bar {r}} _ {2}} va fi în raza de acțiune {\ displaystyle 0 <\ alpha <180 ^ {\ circ}} , în timp ce în al doilea caz va fi inclus în gamă {\ displaystyle 180 ^ {\ circ} <\ alpha <360 ^ {\ circ}} . Prin urmare{\ displaystyle {\ bar {r}} (t)} va continua să treacă {\ displaystyle {\ bar {r}} _ {2}} cu fiecare revoluție orbitală.
În cazul în care {\ displaystyle {\ bar {r}} _ {1} \ times {\ bar {r}} _ {2} = 0} , cum ar fi în cazul {\ displaystyle {\ bar {r}} _ {1}} Și {\ displaystyle {\ bar {r}} _ {2}} au direcție opusă, toate planurile orbitale care conțin linia corespunzătoare sunt la fel de potrivite și unghiul de transfer {\ displaystyle \ alpha} pentru prima trecere prin {\ displaystyle {\ bar {r}} _ {2}} va fi de {\ displaystyle 180 ^ {\ circ}} .
Pentru fiecare {\ displaystyle \ alpha} așa să fie {\ displaystyle 0 <\ alpha <\ infty} , în triunghiul format de puncte {\ displaystyle P_ {1}} , {\ displaystyle P_ {2}} Și {\ displaystyle F_ {1}} , așa cum se arată în Figura 1, avem
- {\ displaystyle d = {\ frac {\ sqrt {{r_ {1}} ^ {2} + {r_ {2}} ^ {2} -2r_ {1} r_ {2} \ cos \ alpha}} {2 }} \ quad (7)}
iar axa semi-majoră (semnată) a hiperbolei este
- {\ displaystyle A = {\ frac {r_ {2} -r_ {1}} {2}} \ quad (8)}
Excentricitatea (cu semn) a hiperbolei este
- {\ displaystyle E = {\ frac {d} {A}} \ quad (9)}
iar axa semi-minoră este
- {\ displaystyle B = \ left | A \ right | {\ sqrt {E ^ {2} -1}} = {\ sqrt {d ^ {2} -A ^ {2}}} \ quad (10)}
Coordonatele punctului {\ displaystyle F_ {1}} relativ la sistemul de coordonate canonice pentru hiperbolă sunt
- {\ displaystyle x_ {0} = - {\ frac {r_ {m}} {E}} \ quad (11)}
(Rețineți că {\ displaystyle E} are semnul de {\ displaystyle r_ {2} -r_ {1}} )
- {\ displaystyle y_ {0} = B {\ sqrt {{\ left ({\ frac {x_ {0}} {A}} \ right)} ^ {2} -1}} \ quad (12)}
unde este
- {\ displaystyle r_ {m} = {\ frac {r_ {2} + r_ {1}} {2}} \ quad (13)}
Alegând coordonata y a punctului {\ displaystyle F_ {2}} pe cealaltă ramură a hiperbolei ca parametru liber, coordonata x a {\ displaystyle F_ {2}} Și
- {\ displaystyle x = A {\ sqrt {1 + {\ left ({\ frac {y} {B}} \ right)} ^ {2}}} \ quad (14)}
(Rețineți că {\ displaystyle A} are semnul {\ displaystyle r_ {2} -r_ {1}} )
Axa semi-majoră a elipsei care trece prin puncte {\ displaystyle P_ {1}} Și {\ displaystyle P_ {2}} și având focarele în puncte {\ displaystyle F_ {1}} Și {\ displaystyle F_ {2}} Și
- {\ displaystyle a = {\ frac {r_ {1} + s_ {1}} {2}} = {\ frac {r_ {2} + s_ {2}} {2}} = {\ frac {r_ {m } + Ex} {2}} \ quad (15)}
Distanța dintre focare este
- {\ displaystyle {\ sqrt {(x_ {0} -x) ^ {2} + (y_ {0} -y) ^ {2}}} \ quad (16)}
iar excentricitatea este în consecință
- {\ displaystyle e = {\ frac {\ sqrt {(x_ {0} -x) ^ {2} + (y_ {0} -y) ^ {2}}} {2a}} \ quad (17)}
Adevarata anomalie {\ displaystyle \ theta _ {1}} în sens {\ displaystyle P_ {1}} depinde de direcția mișcării, adică de semnul {\ displaystyle \ sin \ alpha} . În ambele cazuri avem asta
- {\ displaystyle \ cos \ theta _ {1} = - {\ frac {(x_ {0} + d) f_ {x} + y_ {0} f_ {y}} {r_ {1}}} \ quad (18 )}}
unde este
- {\ displaystyle f_ {x} = {\ frac {x_ {0} -x} {\ sqrt {(x_ {0} -x) ^ {2} + (y_ {0} -y) ^ {2}}} } \ quad (19)}
- {\ displaystyle f_ {y} = {\ frac {y_ {0} -y} {\ sqrt {(x_ {0} -x) ^ {2} + (y_ {0} -y) ^ {2}}} } \ quad (20)}
este vectorul unitar care reprezintă direcția de la {\ displaystyle F_ {2}} la {\ displaystyle F_ {1}} exprimată în coordonate canonice.
De sine {\ displaystyle \ sin \ alpha} atunci este pozitiv
- {\ displaystyle \ sin \ theta _ {1} = {\ frac {(x_ {0} + d) f_ {y} -y_ {0} f_ {x}} {r_ {1}}} \ quad (21) }
De sine {\ displaystyle \ sin \ alpha} în schimb este negativ
- {\ displaystyle \ sin \ theta _ {1} = - {\ frac {(x_ {0} + d) f_ {y} -y_ {0} f_ {x}} {r_ {1}}} \ quad (22 )}}
Fiind
- axa semi-majoră {\ displaystyle a} ,
- excentricitatea orbitală {\ displaystyle e} ,
- anomalia adevărată inițială {\ displaystyle \ theta}
funcțiile cunoscute ale parametrului y, timpul necesar pentru ca adevărata anomalie să crească cu o valoare {\ displaystyle \ alpha} este, de asemenea, o funcție cunoscută a lui y. De sine {\ displaystyle t_ {2} -t_ {1}} este în intervalul care poate fi obținut cu o orbită Kepleriană eliptică, valoarea y corespunzătoare poate fi obținută cu un algoritm iterativ .
În cazul particular pe care îl aveți {\ displaystyle r_ {1} = r_ {2}} (sau foarte asemănător), avem asta {\ displaystyle A = 0} iar hiperbola este redusă la o singură linie ortogonală la îmbinarea {\ displaystyle P_ {1}} cu {\ displaystyle P_ {2}} cu ecuație
- {\ displaystyle x = 0 \ quad (1 ')}
Ecuațiile (11) și (12) sunt apoi înlocuite cu
- {\ displaystyle x_ {0} = 0 \ quad (11 ')}
- {\ displaystyle y_ {0} = {\ sqrt {{r_ {m}} ^ {2} -d ^ {2}}} \ quad (12 ')}
(14) se înlocuiește cu
- {\ displaystyle x = 0 \ quad (14 ')}
și (15) cu
- {\ displaystyle a = {\ frac {r_ {m} + {\ sqrt {d ^ {2} + y ^ {2}}}} {2}} \ quad (15 ')}
Exemplu numeric
Să presupunem următoarele valori pentru o orbită kepleriană centrată pe Pământ:
- r 1 = 10 000 km
- r 2 = 10 000 km
- α = 100 °
Alegând o valoare de y = 30 000 km, se obține un timp de transfer de 3072 secunde, presupunând constanta gravitațională egală cu μ = 398 603 km 3 / s 2 . Parametrii orbitali corespunzători sunt:
- axa semi-majoră = 23 001 km;
- excentricitate = 0,566613;
- anomalie adevărată la momentul t 1 = -7,577 °;
- anomalie adevărată la momentul t 2 = 92,423 °.
Plasând în schimb:
- r 1 = 10 000 km
- r 2 = 10 000 km
- α = 260 °
obținem aceeași orbită eliptică, dar rotația are loc în direcția opusă și avem:
- anomalie adevărată la momentul t 1 = 7,577 °;
- anomalie adevărată la momentul t 2 = 267,577 ° = 360 ° - 92,423 °
și un timp de transfer de 31 645 secunde.
Componentele radiale și tangențiale ale vitezei sunt date de:
- {\ displaystyle V_ {r} = {\ sqrt {\ frac {\ mu} {p}}} \ cdot și \ cdot \ sin \ theta}
- {\ displaystyle V_ {t} = {\ sqrt {\ frac {\ mu} {p}}} \ cdot (1 + e \ cdot \ cos \ theta)}
Aplicații practice
În mod normal, acest algoritm este utilizat pentru a rezolva problema Lambert în proiectarea misiunilor interplanetare. De exemplu, o navă spațială care se deplasează de pe Pământ pe Marte poate fi considerată ca o primă aproximare într-o orbită Kepleriană eliptică eliocentrică de la poziția Pământului la lansare până la poziția lui Marte la sosire. Prin compararea vectorilor vitezei inițiale și finale pentru Pământ și Marte, se poate obține o estimare bună a energiei de lansare necesare și a manevrelor necesare pentru a intra pe orbita marțiană.
Elemente conexe
linkuri externe