삼체 문제의 수학적 모델링: 방정식에서 주기 궤도까지
삼체 문제는 뉴턴의 만유인력 법칙에 따라 세 개의 질량 있는 물체의 궤적을 계산하는 문제입니다. 초기 조건은 3차원 공간에서 물체들의 위치와 속도를 지정하며, 이는 18개의 비선형 2차 미분방정식 시스템으로 이어집니다.
전체 문제에서 각 물체의 좌표와 속도는 상호 거리에 의존하는 가속도와 연결됩니다. 정규화된 질량(최대 질량 대비 비율)이 계산을 단순화합니다. 다음은 1차 미분 형태의 시스템입니다:
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간소화된 변형: 평면 문제와 제한 문제
평면 문제(P3BP)는 2차원 공간으로 제한되어 방정식을 12개로 줄입니다. 여기서 물체들은 초기 조건으로 정의된 평면에 머무릅니다.
평면 문제의 경우:
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제한 문제(PR3BP)는 세 번째 물체의 질량이 작다는 가정을 합니다. 그 운동 방정식:
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원형 제한 문제(CR3BP)는 두 물체 시스템에서 작은 물체가 원형 궤도에 있는 우주선을 모델링합니다. 매개변수 μ = 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이 시스템들은 Mathematica, Python(SciPy), 또는 MATLAB에서 수치 방법으로 풀립니다.
해석적 해: 오일러와 라그랑주
알려진 유일한 해석적 해는 오일러의 직선 배치 해(1767)와 라그랑주의 삼각형 해(1772)입니다. 라그랑주 경우에는 물체들이 정삼각형 꼭짓점에서 시작하며 속도가 일치합니다. 궤적은 120° 회전된 타원이며, 물체들은 삼각형 구성을 유지합니다.
이 해들은 불안정합니다: 초기 조건의 작은 편차가 카오스를 초래합니다. 실제 구현에는 10^{-10} 수준의 정밀도가 필요합니다.
평면 문제에서의 주기 궤적
앙리 푸앵카레는 동일 질량에서 무한히 많은 주기 해의 존재를 증명했습니다. 컴퓨터 탐색으로 대칭성을 가진 패밀리가 발견되었습니다.
- Broucke-Enò-Hadjidimitriou family (1975): 8- 및 16-심볼 궤도. R4 예시 (주기 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
```
- Shuvakov-Dmitrashinovich family (2013): Dragonfly (주기 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
```
기타: Butterfly, Moth, Goggles, Yin-Yang.
- Shin family (2016): Oval, cat's face (주기 5.1).
재현하려면 적응형 스텝 크기의 8차 룬게-쿠타법을 사용하세요.
주요 포인트
- 전체 문제는 18개 방정식이 필요합니다; 간소화로 실용 시뮬레이션에서 4–6개로 줄입니다.
- 오일러/라그랑주 해석적 해는 불안정하며 이상 조건에만 적용됩니다.
- 주기 궤적은 수치 방법으로 발견되며, Broucke와 Shuvakov 패밀리에 수백 개 변형이 있습니다.
- 주기성을 위해 초기 조건을 높은 정밀도로 지정해야 합니다.
- CR3BP는 NASA/ESA 임무에서 중요합니다: 지구-달 및 태양-목성 시스템에서의 기동.
— Editorial Team
아직 댓글이 없습니다.