metoda sklejanych krzywych stożkowychprędkość nadmiarowa v∞sfera wpływu
Astronautyka 7
Sfery wpływu i podróż międzyplanetarna
Sonda lecąca na Marsa jest przyciągana jednocześnie przez Ziemię, Słońce i Marsa. Dokładny ruch w polu trzech ciał nie ma rozwiązania w postaci wzorów, ale prawie zawsze jedno ciało przeważa. Metoda sklejanych stożkowych dzieli lot na trzy odcinki, każdy z jednym ciałem centralnym, i skleja je na granicach sfer wpływu. Tak liczy się okna startowe, budżety Δv i czasy lotu wszystkich misji planetarnych, zanim przyjdzie czas na dokładne całkowanie.
Czas
ok. 2 h teorii + ok. 2 h zadań
Zadania
14 zadań w treści + 5 do przemyślenia, z odpowiedziami
Źródła
H. D. Curtis, Orbital Mechanics for Engineering Students, rozdz. 8; C. A. Kluever, Space Flight Dynamics, rozdz. o trajektoriach międzyplanetarnych; U. Walter, Astronautics, rozdz. o lotach międzyplanetarnych
0 Prolog: Mariner 4 i pierwsze zdjęcia Marsa
W listopadzie 1964 roku z przylądka Canaveral (nazywanego wtedy Przylądkiem Kennedy’ego) wystartował Mariner 4. Po prawie ośmiu miesiącach lotu, w lipcu 1965 roku, minął Marsa w odległości ok. 10 000 km i przesłał na Ziemię 21 zdjęć. Pokazywały kratery i martwą powierzchnię, a pomiar przejścia sygnału radiowego przez atmosferę wykazał, że jest ona bardzo rzadka. Skończyły się marzenia o kanałach.
Start nie mógł się odbyć w dowolnym momencie. Mars i Ziemia muszą stać w odpowiednim ustawieniu, które powtarza się co ok. 26 miesięcy. Ta lekcja liczy, skąd biorą się te liczby: 26 miesięcy, 8 miesięcy lotu i ok. 3,6 km/s z niskiej orbity okołoziemskiej.
Pytanie na start. Ucieczka z orbity 300 km kosztuje 3,20 km/s, a sonda ma opuścić pole Ziemi z prędkością 2,95 km/s. Czy potrzeba 3,20+2,95=6,15 km/s? Odpowiedź w części 3.
1 Sfera wpływu
W pobliżu Ziemi sondę przyciąga głównie Ziemia, daleko głównie Słońce. Gdzie przebiega granica? Laplace zaproponował kryterium: liczymy ruch sondy względem Ziemi z Słońcem jako zaburzeniem i względem Słońca z Ziemią jako zaburzeniem. Granica leży tam, gdzie stosunek zaburzenia do przyspieszenia głównego jest w obu opisach taki sam. Wynik:
Promień sfery wpływu
rSOI=a(Mm)2/5,
a odległość planety od Słońca (lub księżyca od planety), m/M stosunek mas. Ziemia: 925 000 km (145 promieni Ziemi); Księżyc względem Ziemi: 66 000 km; Mars: 577 000 km; Jowisz: 48 mln km.
Sfera wpływu Ziemi jest olbrzymia w porównaniu z Ziemią, ale mała w porównaniu z odległością do Słońca: 0,6% jednostki astronomicznej. To dlatego metoda sklejanych stożkowych działa: w skali Układu Słonecznego sonda wychodzi z „punktu”, a w skali Ziemi pole Słońca jest drobnym zaburzeniem. Granica nie jest fizyczną powierzchnią, tylko wygodnym miejscem na zmianę opisu; wykładnik 2/5 wynika z porównania przyspieszeń pływowych.
Zadania · część 1
3 zadania · ok. 20 min
Policz promień sfery wpływu Ziemi. Masa Ziemi to 3,00⋅10−6 masy Słońca (z μ⊕/μ⊙), a=1,496⋅108 km.
Odpowiedź
(3,00⋅10−6)0,4=6,18⋅10−3, rSOI=9,25⋅105 km: ok. 2,4 odległości do Księżyca.
Sfera wpływu Księżyca ma promień 66 000 km, a Księżyc jest 384 400 km od Ziemi. Opisz trajektorię lotu na Księżyc w języku sklejanych stożkowych.
Odpowiedź
Najpierw elipsa wokół Ziemi z apogeum w okolicy orbity Księżyca. Po wejściu w sferę wpływu Księżyca (ok. 320 000 km od Ziemi) statek ma względem Księżyca prędkość nadmiarową i porusza się po hiperboli wokół niego; bez hamowania minie go i wyjdzie z jego sfery. Dla Księżyca przybliżenie jest gorsze niż dla planet, bo jego sfera nie jest mała względem odległości Ziemia–Księżyc.
Dlaczego sfera wpływu Jowisza (48 mln km) jest prawie 50 razy większa od ziemskiej, mimo że Jowisz jest tylko 5 razy dalej od Słońca?
Odpowiedź
Jowisz jest 318 razy masywniejszy od Ziemi: 3180,4=10,0; do tego czynnik odległości 5,2. Razem ok. 52 razy.
2 Metoda sklejanych stożkowych
Lot z Ziemi na Marsa dzielimy na trzy odcinki, w każdym zagadnienie dwóch ciał z lekcji 1:
Sklejane stożkowe
1. Wylot: hiperbola względem Ziemi, od orbity parkingowej do granicy sfery wpływu; na granicy prędkość względem Ziemi to v∞.
2. Lot heliocentryczny: elipsa wokół Słońca z prędkością początkową v⊕+v∞ (sfera wpływu traktowana jak punkt w miejscu Ziemi).
3. Przylot: hiperbola względem Marsa z prędkością nadmiarową v∞M=vsondy−vMars.
Dla orbit kołowych w jednej płaszczyźnie najtańszy odcinek heliocentryczny to elipsa Hohmanna (lekcja 5 Manewry orbitalne: transfer Hohmanna i zmiana płaszczyzny orbity) z peryhelium przy Ziemi i aphelium przy Marsie. Ziemia ma prędkość 29,78 km/s; sonda w peryhelium elipsy o a=1,262 AU ma 32,73 km/s, więc v∞=2,95 km/s w kierunku ruchu Ziemi. Przy Marsie sonda ma w aphelium 21,48 km/s, Mars 24,13 km/s: sonda przylatuje wolniej, a Mars „dogania” ją z v∞=2,65 km/s. Lot trwa pół okresu elipsy: 259 dni.
Symulacja 1 Transfer Hohmanna z Ziemi do innych planet (orbity kołowe w jednej płaszczyźnie). Puste kółka: położenia Ziemi i celu w chwili startu, pełne: w wybranej chwili lotu. Zwróć uwagę na kąt fazowy: gdzie musi być cel w chwili startu.
Zadania · część 2
3 zadania · ok. 25 min
Policz v∞ przy Ziemi dla transferu Hohmanna na Marsa (aM=1,524 AU, μ⊙=1,327⋅1011 km³/s²).
Ile trwa lot Hohmanna na Marsa i na Wenus (aW=0,723 AU)?
Odpowiedź
Mars: at=1,262 AU, pół okresu =21⋅1,2621,5 roku =0,709 roku =259 dni. Wenus: at=0,862 AU, 21⋅0,8621,5=0,40 roku =146 dni.
Dlaczego lot na Wenus wymaga hamowania względem Ziemi, a mimo to v∞ jest dodatnie (2,50 km/s)?
Odpowiedź
Sonda ma w aphelium elipsy przejściowej 27,29 km/s, mniej niż Ziemia (29,78 km/s): względem Ziemi leci do tyłu z prędkością 2,50 km/s. v∞ to długość wektora; jego kierunek jest przeciwny do ruchu Ziemi. Hiperbola wylotu jest więc obrócona: perygeum po stronie przeciwnej niż przy locie na Marsa.
3 Wylot z Ziemi
Sonda zaczyna na orbicie parkingowej o promieniu rp. W perygeum hiperboli wylotu prędkość wynika z energii: vp2=v∞2+2μ/rp. Koszt wylotu:
Δv=v∞2+rp2μ−rpμ.
Odpowiedź na pytanie z prologu: nie. Dla orbity 300 km i v∞=2,95 km/s: 8,67+119,38−7,73=3,59 km/s, a nie 6,15. Prędkości nie dodają się, bo dodają się energie (efekt Obertha z lekcji 5): impuls w głębi studni potencjału daje dużo więcej energii niż ten sam impuls daleko od Ziemi. Hiperbola ma mimośród e=1+rpv∞2/μ=1,145, a jej asymptota wyjściowa jest 150,8∘ od perygeum; silnik zapala się więc ok. 151∘ (licząc wzdłuż orbity) przed kierunkiem, w którym sonda ma odlecieć.
Symulacja 2 Hiperbola wylotu z orbity parkingowej (Ziemia w skali). Asymptota wskazuje kierunek v∞, czyli kierunek ruchu Ziemi wokół Słońca przy locie na planetę zewnętrzną. Zmień v∞ i porównaj koszt z naiwnym „ucieczka + v∞”.
Ile kosztowałby wylot na Marsa z orbity geostacjonarnej (bez kosztu dolotu na nią)? Dlaczego z niskiej orbity jest taniej?
Odpowiedź
r=42164 km: 8,67+18,91−3,075=2,18 km/s, mniej niż 3,59 km/s, ale dolot na GEO kosztuje 3,9 km/s. Z niskiej orbity: 3,59 km/s łącznie, bo cały impuls jest w głębi studni potencjału, gdzie prędkość jest największa (efekt Obertha).
4 Okna startowe i czas misji
Sonda ma spotkać Marsa po 259 dniach, po drugiej stronie Słońca (180∘ od miejsca startu). Mars pokonuje w tym czasie nM⋅259 d =135,7∘, więc w chwili startu musi być 180∘−135,7∘=44,3∘ przed Ziemią. Takie ustawienie powtarza się co okres synodyczny:
Tsyn=∣n1−n2∣2π,Ziemia–Mars: 780 dni (2,14 roku),Ziemia–Wenus: 584 dni.
Misja z powrotem musi czekać na planecie na okno powrotne. Po przylocie (Hohmann) Ziemia jest 75∘ przed Marsem, a do startu powrotnego musi być 75∘ za nim. Ziemia dogania Marsa o 0,46∘ na dobę, więc czekanie trwa 454 dni, a cała misja 259+454+259=972 dni, ok. 2,7 roku. Tak wygląda podstawowy scenariusz załogowej misji marsjańskiej.
Prawdziwe orbity planet są eliptyczne i lekko nachylone (Mars: e=0,093, i=1,85∘), więc okna różnią się między sobą: wymagane C3 zmienia się od ok. 8 do ok. 16 km²/s². Dokładne okna liczy się z położeń planet i problemu Lamberta (lekcja 9 Problem Lamberta i wykres porkchop).
Zadania · część 4
3 zadania · ok. 25 min
Policz okres synodyczny Ziemia–Mars (okres Marsa 687 dni).
Odpowiedź
1/Tsyn=1/365,26−1/687=1,282⋅10−3 /d, Tsyn=780 dni.
Jaki kąt fazowy jest potrzebny przy starcie na Wenus? Gdzie jest Wenus w chwili startu?
Odpowiedź
Wenus pokonuje w 146 dni 146⋅360∘/224,7=234∘. Ma być 180∘ od miejsca startu, więc przy starcie jest 180∘−234∘=−54∘: 54∘ za Ziemią.
Dlaczego okno na Jowisza wypada co 399 dni, prawie co rok?
Odpowiedź
Jowisz porusza się wolno (okres 11,86 roku), więc Ziemia dogania go prawie co rok: 1/Tsyn=1−1/11,86 /rok, Tsyn=1,09 roku. Im dalsza planeta, tym okres synodyczny bliższy rokowi.
5 Przylot i wejście na orbitę
Przy Marsie sonda ma v∞=2,65 km/s i porusza się po hiperboli. Odległość perycentrum wybiera się małą korektą toru daleko przed planetą. Żeby zostać na orbicie, trzeba w perycentrum zahamować:
Δvwej=v∞2+rp2μ−vorb(rp).
Wejście na orbitę kołową 300 km nad Marsem kosztuje 2,09 km/s, a na wydłużoną elipsę 300×33000 km tylko 0,91 km/s. Dlatego sondy najpierw wchodzą na wydłużoną orbitę, a potem obniżają apocentrum hamowaniem w atmosferze (ang. aerobraking): setki przelotów przez górne warstwy atmosfery zamiast setek kilogramów paliwa. Cały budżet z niskiej orbity okołoziemskiej na kołową orbitę 300 km wokół Marsa (Hohmann, wszystko silnikiem): 3,59+2,09=5,68 km/s.
Zadania · część 5
2 zadania · ok. 20 min
Policz Δv wejścia na orbitę kołową 300 km wokół Marsa (μM=42828 km³/s², RM=3396 km).
Dlaczego przy Jowiszu wejście na niską orbitę kołową kosztowałoby prawie 17 km/s, a sondy wchodzą na orbitę za ok. 0,5–0,6 km/s?
Odpowiedź
Przy rp=1,1RJ prędkość w perycentrum hiperboli to ok. 57 km/s, a kołowa ok. 40 km/s. Zatrzymanie na bardzo wydłużonej orbicie (perycentrum 1,1RJ, apocentrum 100RJ, okres 44 doby) kosztuje 0,59 km/s: wystarczy zmniejszyć energię z +v∞2/2 do lekko ujemnej, a efekt Obertha przy ogromnej prędkości w perycentrum robi resztę.
6 Metoda i przykłady rozwiązane
Metoda: wstępny projekt misji planetarnej
Odcinek heliocentryczny: elipsa (Hohmann albo z problemu Lamberta); v∞ na obu końcach.
Wylot:Δv=v∞2+2μ/rp−μ/rp z orbity parkingowej; kierunek asymptoty wyznacza miejsce zapłonu.
Przylot: hamowanie w perycentrum, najlepiej na wydłużoną orbitę; ewentualnie hamowanie atmosferą.
Terminy: kąt fazowy, okres synodyczny, czekanie na okno powrotne.
Sprawdzenie: sfera wpływu dużo mniejsza od odległości planety od Słońca; potem dokładne całkowanie.
Przykład 6.1 Budżet misji na Marsa
Sonda startuje z orbity 300 km wokół Ziemi i ma wejść na orbitę kołową 300 km wokół Marsa. Policz budżet i terminy (Hohmann).
Odcinek heliocentryczny: v∞=2,95 km/s (Ziemia), 2,65 km/s (Mars), 259 dni.
Wylot 3,59 km/s, wejście na orbitę 2,09 km/s: razem 5,68 km/s.
Start, gdy Mars jest 44∘ przed Ziemią; okna co 780 dni.
Sprawdzenie: naiwne dodawanie prędkości dałoby 6,15 km/s tylko na wylot; różnica to efekt Obertha ✓.
Przykład 6.2 Lot na Wenus
Policz podstawowe liczby lotu Hohmanna na Wenus i wejścia na orbitę 300 km.
v∞: 2,50 km/s (Ziemia, do tyłu względem ruchu Ziemi), 2,71 km/s (Wenus); 146 dni; Wenus 54∘ za Ziemią przy starcie; okna co 584 dni.
Wylot z 300 km: 3,48 km/s.
Wejście na orbitę kołową 300 km nad Wenus (μ=324859 km³/s²): 3,32 km/s.
Wynik: ok. 6,8 km/s; Wenus jest prawie tak „blisko” jak Mars na wylocie, ale droższa na przylocie (większa masa planety).
Sprawdzenie: większe μ zwiększa i prędkość w perycentrum hiperboli, i prędkość kołową; ich różnica rośnie ✓.
Przykład 6.3 Bezpośrednio na Jowisza
Ile kosztuje bezpośredni lot Hohmanna na Jowisza i ile trwa?
v∞=8,79 km/s przy Ziemi, C3=77 km²/s²; wylot z 300 km: 6,30 km/s.
Czas 2,73 roku; przy Jowiszu v∞=5,64 km/s.
Wejście na orbitę 1,1×100RJ: 0,59 km/s.
Wynik: ok. 6,9 km/s i prawie 3 lata. Takie C3 (i większe) osiągały lekkie sondy lecące bezpośrednio: Pioneery 10 i 11, Voyagery, New Horizons. Ciężkie sondy (Galileo, Cassini, Juno) przy takim C3 nie zmieściłyby się w możliwościach rakiet, dlatego leciały przez asysty grawitacyjne (lekcja 8 Asysta grawitacyjna).
Sprawdzenie:v∞ przy Ziemi rośnie z odległością celu (Mars 2,9, Jowisz 8,8, Saturn 10,3 km/s), a granica 2−1 prędkości Ziemi, ok. 12,3 km/s, to ucieczka z Układu Słonecznego ✓.
7 Błędne intuicje, których trzeba się oduczyć
„Na Marsa leci się, gdy jest najbliżej.”
Startuje się, gdy Mars jest 44∘ przed Ziemią, a dolatuje po drugiej stronie Słońca. Największe zbliżenie (opozycja) wypada w trakcie lotu.
„Do ucieczki i do v∞ potrzeba sumy prędkości.”
Dodają się energie, nie prędkości: Δv=v∞2+vu2−vk.
„Sfera wpływu to fizyczna granica.”
To umowne miejsce zmiany opisu; obie siły działają wszędzie, zmienia się tylko, która jest przeważająca.
„Lot na Wenus jest łatwy, bo jest bliżej Słońca.”
Trzeba stracić prędkość względem Słońca; do Merkurego trzeba jej stracić tyle, że przylot jest droższy niż lot na Jowisza.
„Misja na Marsa i z powrotem trwa dwa razy po 8 miesięcy.”
Plus ok. 15 miesięcy czekania na okno powrotne: razem ok. 2,7 roku.
8 Zadania do przemyślenia
Te zadania wymagają więcej myślenia niż liczenia.
Zadania do przemyślenia
5 zadań · ok. 1,5 godziny
Skąd wykładnik 2/5? Naszkicuj, dlaczego kryterium Laplace’a daje rSOI∝(m/M)2/5.
Rozwiązanie
Względem planety: przyspieszenie główne Gm/r2, zaburzenie od Słońca (pływowe) ∼GMr/a3, stosunek ∼(M/m)(r/a)3. Względem Słońca: główne GM/a2, zaburzenie od planety ∼Gm/r2, stosunek ∼(m/M)(a/r)2. Przyrównanie: (M/m)(r/a)3=(m/M)(a/r)2, czyli (r/a)5=(m/M)2 i r=a(m/M)2/5.
Szybciej niż Hohmann. Załoga nie powinna lecieć 259 dni. Co się dzieje z v∞ przy szybszym transferze na Marsa?
Rozwiązanie
Szybszy transfer to elipsa z aphelium dalej niż Mars, przecinająca orbitę Marsa przed aphelium: większe v∞ na wylocie i przy Marsie (i nieprostopadły przylot).
Koszt rośnie szybko: skrócenie lotu o kilka miesięcy podnosi v∞ o 1–3 km/s na obu końcach. Wymagana zmiana prędkości przy Marsie może wymusić hamowanie aerodynamiczne. Rachunek dla konkretnych dat robi się problemem Lamberta (lekcja 9).
Kierunek asymptoty. Sonda na Marsa ma odlecieć w kierunku ruchu Ziemi wokół Słońca. W jakiej porze dnia (czas lokalny na orbicie parkingowej) zapala się silnik?
Rozwiązanie
Kierunek ruchu Ziemi wokół Słońca jest prostopadły do kierunku do Słońca i wskazuje stronę porannego terminatora (ok. 6:00 lokalnego czasu słonecznego; dlatego nad ranem widać więcej meteorów). Tam ma wskazywać asymptota wylotu.
Perygeum hiperboli leży 151∘ wcześniej wzdłuż orbity, czyli ok. 10 h lokalnego czasu wcześniej: ok. 20:00, po nocnej stronie Ziemi. Dlatego zapłony w kierunku Marsa wypadają nad nocną półkulą.
Księżyc a sklejane stożkowe. Dlaczego przybliżenie sklejanych stożkowych jest dla lotu na Księżyc wyraźnie gorsze niż dla Marsa?
Rozwiązanie
Stosunek promienia sfery wpływu do odległości: Mars–Słońce 0,25%, Księżyc–Ziemia 17%. Statek spędza długi czas w obszarze, gdzie obie siły są porównywalne.
Do tego masa Księżyca to 1/81 Ziemi (Mars to 3⋅10−7 Słońca): zaburzenia są duże. Trajektorie księżycowe (np. powrót „swobodny” Apollo) liczy się od razu numerycznie w zagadnieniu trzech ciał.
Wenus w drodze na Marsa. Istnieją trajektorie z Ziemi na Marsa przez przelot obok Wenus. Co mogą dać i jaką mają wadę?
Rozwiązanie
Przelot obok Wenus zmienia kierunek v∞ względem Wenus za darmo (lekcja 8): może skrócić czekanie na okno, zmniejszyć v∞ przy przylocie albo umożliwić powrót bez długiego postoju na Marsie.
Wada: trzy planety muszą stać w odpowiednim ustawieniu (rzadkie okna), lot jest dłuższy, a sonda zbliża się do Słońca (termika).
9 Sprawdź się
Odpowiedz w myślach, zanim rozwiniesz odpowiedź.
1 Wzór na promień sfery wpływu?
rSOI=a(m/M)2/5; Ziemia ok. 925 000 km.
2 Jakie są trzy odcinki metody sklejanych stożkowych?
Hiperbola wylotu, elipsa heliocentryczna, hiperbola przylotu; sklejone przez v∞.
3 Ile kosztuje wylot na Marsa z orbity 300 km?
3,59 km/s (v∞=2,95 km/s, C3=8,7 km²/s²).
4 Co to jest okres synodyczny i ile wynosi dla Marsa?
Czas powtarzania się ustawienia: 2π/∣n1−n2∣; Mars 780 dni.
5 Jak tanio wejść na orbitę planety?
Hamować w perycentrum na wydłużoną orbitę, potem ewentualnie hamowanie atmosferyczne.
Jowisz: 8,79/5,64 km/s, 2,73 roku, co 399 d, wylot 6,30 km/s
Tsyn=2π/∣n1−n2∣
Do przeczytania
Obowiązkowo: H. D. Curtis, Orbital Mechanics for Engineering Students, rozdz. 8 (trajektorie międzyplanetarne, sfera wpływu, okres synodyczny, wylot i przylot); C. A. Kluever, Space Flight Dynamics, rozdz. o trajektoriach międzyplanetarnych.
Dla chętnych: R. R. Bate, D. D. Mueller, J. E. White, Fundamentals of Astrodynamics, rozdz. o trajektoriach międzyplanetarnych; materiały NASA o misjach Mariner i Mars Exploration Program.
Hiperbola przylotu ma jeszcze jedną własność: obraca wektor v∞. Względem planety długość v∞ się nie zmienia, ale względem Słońca prędkość sondy może wzrosnąć o kilka kilometrów na sekundę. Kolejna lekcja zamienia to w napęd: asystę grawitacyjną, dzięki której Voyagery doleciały do granic Układu Słonecznego.