Problem trzech ciał: równania i okresowe orbity

Matematyczne modelowanie problemu trzech ciał: od równań do okresowych orbit

Problem trzech ciał wymaga obliczenia trajektorii trzech masywnych obiektów pod wpływem grawitacyjnego przyciągania zgodnie z prawem Newtona. Warunki początkowe określają pozycje i prędkości ciał w przestrzeni 3D, prowadząc do układu 18 nieliniowych równań różniczkowych drugiego rzędu.

Dla pełnego problemu współrzędne i prędkości każdego ciała są powiązane przyspieszeniami zależnymi od wzajemnych odległości. Znormalizowane masy (stosunek do maksymalnej) upraszczają obliczenia. Oto system w formie pierwszych pochodnych:

X1'  = VX1
Y1'  = VY1
Z1'  = VZ1
X2'  = VX2
Y2'  = VY2
Z2'  = VZ2
X3'  = VX3
Y3'  = VY3
Z3'  = VZ3
VX1' = -m2*(X1-X2)/((X1-X2)^2+(Y1-Y2)^2+(Z1-Z2)^2)^1.5-m3*(X1-X3)/((X1-X3)^2+(Y1-Y3)^2+(Z1-Z3)^2)^1.5
VY1' = -m2*(Y1-Y2)/((X1-X2)^2+(Y1-Y2)^2+(Z1-Z2)^2)^1.5-m3*(Y1-Y3)/((X1-X3)^2+(Y1-Y3)^2+(Z1-Z3)^2)^1.5
VZ1' = -m2*(Z1-Z2)/((X1-X2)^2+(Y1-Y2)^2+(Z1-Z2)^2)^1.5-m3*(Z1-Z3)/((X1-X3)^2+(Y1-Y3)^2+(Z1-Z3)^2)^1.5
VX2' = -m3*(X2-X3)/((X2-X3)^2+(Y2-Y3)^2+(Z2-Z3)^2)^1.5-m1*(X2-X1)/((X1-X2)^2+(Y1-Y2)^2+(Z1-Z2)^2)^1.5
VY2' = -m3*(Y2-Y3)/((X2-X3)^2+(Y2-Y3)^2+(Z2-Z3)^2)^1.5-m1*(Y2-Y1)/((X1-X2)^2+(Y1-Y2)^2+(Z1-Z2)^2)^1.5
VZ2' = -m3*(Z2-Z3)/((X2-X3)^2+(Y2-Y3)^2+(Z2-Z3)^2)^1.5-m1*(Z2-Z1)/((X1-X2)^2+(Y1-Y2)^2+(Z1-Z2)^2)^1.5
VX3' = -m1*(X3-X1)/((X1-X3)^2+(Y1-Y3)^2+(Z1-Z3)^2)^1.5-m2*(X3-X2)/((X3-X2)^2+(Y3-Y2)^2+(Z3-Z2)^2)^1.5
VY3' = -m1*(Y3-Y1)/((X1-X3)^2+(Y1-Y3)^2+(Z1-Z3)^2)^1.5-m2*(Y3-Y2)/((X3-X2)^2+(Y3-Y2)^2+(Z3-Z2)^2)^1.5
VZ3' = -m1*(Z3-Z1)/((X1-X3)^2+(Y1-Y3)^2+(Z1-Z3)^2)^1.5-m2*(Z3-Z2)/((X3-X2)^2+(Y3-Y2)^2+(Z3-Z2)^2)^1.5

Uproszczone warianty: płaski i ograniczony problem

Płaski problem (P3BP) ogranicza się do przestrzeni 2D, skracając równania do 12. Tutaj ciała pozostają w płaszczyźnie określonej warunkami początkowymi.

Dla płaskiego problemu:

X1'  = VX1
X2'  = VX2
X3'  = VX3
Y1'  = VY1
Y2'  = VY2
Y3'  = VY3
VX1' = m2*(X2-X1)/((X1-X2)^2+(Y1-Y2)^2)^1.5+m3*(X3-X1)/((X1-X3)^2+(Y1-Y3)^2)^1.5
VX2' = m3*(X3-X2)/((X2-X3)^2+(Y2-Y3)^2)^1.5+m1*(X1-X2)/((X2-X1)^2+(Y2-Y1)^2)^1.5
VX3' = m1*(X1-X3)/((X3-X1)^2+(Y3-Y1)^2)^1.5+m2*(X2-X3)/((X3-X2)^2+(Y3-Y2)^2)^1.5
VY1' = m2*(Y2-Y1)/((X1-X2)^2+(Y1-Y2)^2)^1.5+m3*(Y3-Y1)/((X1-X3)^2+(Y1-Y3)^2)^1.5
VY2' = m3*(Y3-Y2)/((X2-X3)^2+(Y2-Y3)^2)^1.5+m1*(Y1-Y2)/((X2-X1)^2+(Y2-Y1)^2)^1.5
VY3' = m1*(Y1-Y3)/((X3-X1)^2+(Y3-Y1)^2)^1.5+m3*(Y2-Y3)/((X3-X2)^2+(Y3-Y2)^2)^1.5

Ograniczony problem (PR3BP) zakłada małą masę trzeciego ciała. Równania dla jego ruchu:

X' = VX
Y' = VY
VX' = X+2*Y4-m1*(X+m2)/((X+m2)^2+Y^2)^1.5-m2*(X-m1)/((X-m1)^2+Y^2)^1.5
VY' = Y-2*Y3-m1*Y/((X+m2)^2+Y^2)^1.5-m2*Y/((X-m1)^2+Y^2)^1.5

Kołowy ograniczony problem (CR3BP) modeluje statek w systemie dwóch ciał z kołową orbitą mniejszego. Stosowany jest parametr μ = m2 / (m1 + m2):

X' = VX
Y' = VY
z' = VZ
VX' = 2*VY+X-(1-mu) * (X+mu)/((X+mu)^2+Y^2+Z^2)^1.5-mu * (X-1+mu)/((X-1+mu)^2+Y^2+Z^2)^1.5
VY' = -2*VX+Y-(1-mu) * Y/((X+mu)^2+Y^2+Z^2)^1.5-mu * Y/((X-1+mu)^2+Y^2+Z^2)^1.5
VZ' = -(1-mu) * Z/((X+mu)^2+Y^2+Z^2)^1.5-mu * Z/((X-1+mu)^2+Y^2+Z^2)^1.5

Te układy rozwiązuje się numerycznymi metodami w Mathematica, Python (SciPy) lub MATLAB.

Analityczne rozwiązania: Euler i Lagrange

Jedynymi znanymi analitycznymi rozwiązaniami są kolinearne Eulera (1767) i trójkątne Lagrange’a (1772). W przypadku Lagrange’a ciała wyruszają z wierzchołków równobocznego trójkąta z prędkościami zsynchronizowanymi. Trajektorie to elipsy obrócone o 120°; ciała zachowują trójkątną konfigurację.

Rozwiązania są niestabilne: małe odchylenia warunków początkowych prowadzą do chaosu. W praktyce wymagana jest dokładność rzędu 10^{-10}.

Okresowe orbity w płaskim problemie

Henri Poincaré udowodnił istnienie nieskończonej liczby okresowych rozwiązań dla mas równych. Poszukiwania komputerowe ujawniły rodziny o symetrii.

• Rodzina Bruka-Eno-Hadżidimitriou (1975): orbity 8- i 16-kształtne. Przykład R4 (okres 5.4):

```

T0 = 0, T1 = 5.4

X1(T0) = 0.8733047091, X2(T0) = -0.6254030288, X3(T0) = -0.2479016803

Y1(T0) = 0, Y2(T0) = 0, Y3(T0) = 0

VX1(T0) = 0, VX2(T0) = 0, VX3(T0) = 0

VY1(T0) = 1.010776444, VY2(T0) = -1.683353346, VY3(T0) = 0.6725769022

```

• Rodzina Szuwakowa-Dmitraszinowicza (2013): Dragonfly (okres 22):

```

T0 = 0, T1 = 22

X1(T0) = -1, X2(T0) = 1, X3(T0) = 0

Y1(T0) = 0, Y2(T0) = 0, Y3(T0) = 0

VX1(T0) = 0.08058, VX2(T0) = 0.08058, VX3(T0) = -0.16116

VY1(T0) = 0.58884, VY2(T0) = 0.58884, VY3(T0) = -1.17768

```

Inne: Butterfly, Moth, Goggles, Yin-Yang.

• Rodzina Šina (2016): Ówal, kocia mordka (okres 5.1).

Do odtworzenia użyj metody Rungego-Kutty 8. rzędu z adaptacyjnym krokiem.

Co ważne

• Pełny problem wymaga 18 równań; uproszczenia redukują liczbę do 4–6 w praktycznych symulacjach.
• Analityczne rozwiązania Eulera/Lagrange’a są niestabilne i stosowane tylko w warunkach idealnych.
• Okresowe orbity odkrywane są metodami numerycznymi; rodziny Bruka i Szuwakowa zawierają setki wariantów.
• Warunki początkowe określa się z wysoką dokładnością, aby zapewnić okresowość.
• CR3BP jest kluczowy dla misji NASA/ESA: manewry w systemach Ziemia–Księżyc, Słońce–Jowisz.

— Editorial Team