궤도 역학

물리 계산물리 마지막 수정: 2026-08-13 04:26:41

1. 개요[편집]

궤도 역학 Orbital mechanics / Astrodynamics
기본 문제	이체 문제 (two-body problem)
지배식	r̈ = −μ r / r³
보존량	비각운동량 h, 비에너지 ε, 이심률 벡터 e
해의 형태	원뿔곡선 r = p / (1 + e cos ν)
지구 중력상수	μ_⊕ = 398600.4418 km³/s²
정립	케플러(1609) → 뉴턴(1687) → 라그랑주·가우스

우주에서는 브레이크를 밟으면 앞으로 나가고 액셀을 밟으면 뒤로 처진다. 직관을 통째로 버리고 방정식만 믿어야 하는 몇 안 되는 분야.

궤도 역학(orbital mechanics, 천체역학의 공학 판이라 astrodynamics라고도 한다)은 중력장 안에서 움직이는 물체, 특히 인공위성과 우주선의 운동을 예측·설계·제어하는 학문이다. 출발점은 딱 하나 — 질량이 무시할 만한 물체가 중심 천체의 역제곱 중력만 받는 이체 문제이며, 이것은 닫힌 해가 존재하는 몇 안 되는 역학 문제다.

\ddot{\mathbf r} = -\frac{\mu}{r^3}\mathbf r, \qquad \mu = G(M+m)

궤도 역학의 실전은 이 깔끔한 해를 기준 궤도로 삼고, 현실이 그것과 어긋나는 만큼을 섭동으로 얹는 방식으로 굴러간다. 이 문서는 이체 해와 그 위의 임무 설계(전이·영향권·섭동·발사창)를 다룬다. 중력원이 셋 이상이 되어 닫힌 해가 사라지는 순간부터는 N체 문제의 영역이고, 그 수치적분을 장기간 안정적으로 굴리는 기법은 심플렉틱 적분기·변분 적분기에 있다. 추진제 예산을 질량비로 환산하는 이야기는 치올콥스키 로켓 방정식에 있으니 여기서는 $\Delta v$ 가 얼마나 필요한지만 계산한다.

2. 이체 문제의 적분[편집]

역제곱 중력장은 보존량을 세 벌 내놓는다. 먼저 $\mathbf r \times$ 를 취하면 비각운동량이 나온다.

\mathbf h = \mathbf r \times \mathbf v = \text{const}

$\mathbf h$ 가 상수라는 것은 운동이 $\mathbf h$ 에 수직인 한 평면 안에 갇힌다는 뜻이고, $|\mathbf h|$ 가 상수라는 것이 케플러 제2법칙(면적 속도 일정)이다. 다음으로 $\mathbf v\cdot$ 를 취해 적분하면 비궤도에너지가 나온다.

\varepsilon = \frac{v^2}{2} - \frac{\mu}{r} = -\frac{\mu}{2a}

마지막으로 역제곱 법칙에만 있는 특별한 보존량 — 이심률 벡터(라플라스-룽게-렌츠 벡터)다.

\mathbf e = \frac{\mathbf v \times \mathbf h}{\mu} - \frac{\mathbf r}{r}

이 벡터는 근점 방향을 가리키고 크기가 이심률과 같다. 이게 보존되기 때문에 궤도가 닫힌 타원으로 그려진다. 역제곱에서 조금이라도 벗어나면 $\mathbf e$ 가 천천히 돌기 시작하고, 그 회전이 곧 근점 이동이다.¹ 세 보존량을 조합하면 궤도 방정식이 떨어진다.

r = \frac{p}{1 + e\cos\nu}, \qquad p = \frac{h^2}{\mu}

$e < 1$ 이면 타원, $e = 1$ 이면 포물선, $e > 1$ 이면 쌍곡선. 그리고 에너지식과 궤도 방정식을 합치면 궤도 역학에서 가장 자주 두드리는 한 줄, 비스-비바(vis-viva) 방정식이 나온다.

v^2 = \mu\left(\frac{2}{r} - \frac{1}{a}\right)

“어느 지점에서 얼마나 빠른가”를 반장축 하나로 답해 준다. 전이 궤도 설계의 90%는 이 식을 두 번 쓰고 빼는 것으로 끝난다.

3. 궤도 6요소[편집]

위치·속도 벡터 6개 성분 대신, 물리적 의미가 붙은 6개 스칼라로 궤도를 기술하는 것이 관례다.

요소	이름	무엇을 정하나
$a$	반장축	궤도의 크기·주기·에너지
$e$	이심률	궤도의 납작함
$i$	경사각	궤도면이 기준면에서 기운 정도
$\Omega$	승교점 적경	궤도면이 기준면과 만나는 방향
$\omega$	근점 인수	궤도면 안에서 타원이 놓인 방향
$\nu$ 또는 $M$	진근점각 / 평균근점각	지금 어디쯤 있나

앞의 다섯은 이체 문제에서 상수이고 여섯 번째만 시간에 따라 변한다. 이 분리가 섭동 이론의 골격이 된다 — 섭동이 있으면 앞의 다섯이 “천천히 변하는 상수”가 되고, 그 변화율을 쓰는 것이 가우스/라그랑주 행성방정식이다.

실무의 함정은 특이점이다. $e \to 0$ 이면 근점이 없어 $\omega$ 가 정의되지 않고, $i \to 0$ 이면 승교점이 없어 $\Omega$ 가 정의되지 않는다. 원궤도에 가까운 위성일수록 요소가 수치적으로 요동치므로, 이럴 땐 특이점이 없는 춘분점 요소(equinoctial elements) 같은 조합으로 갈아탄다.

4. 케플러 방정식과 뉴턴 반복[편집]

“시각 $t$ 에 위성이 어디 있나”는 궤도 역학에서 가장 자주 계산하는 질문인데, 놀랍게도 닫힌 해가 없다. 평균 운동 $n = \sqrt{\mu/a^3}$ 과 평균근점각 $M = n(t - t_p)$ 를 놓으면 이심근점각 $E$ 가 만족하는 초월방정식이 나온다.

M = E - e\sin E

이것이 케플러 방정식이다. $E$ 에 대해 대수적으로 못 푼다는 것이 400년 전부터 알려져 있고, 그래서 인류는 이 한 줄을 풀려고 수백 가지 반복법을 발명했다. 표준은 뉴턴법 반복이다.

E_{k+1} = E_k - \frac{E_k - e\sin E_k - M}{1 - e\cos E_k}

수렴은 보통 3~4회면 배정밀도 한계까지 간다. 다만 두 군데서 사고가 난다.

** $e \to 1$ 이고 $M \to 0$ **인 구간. 분모 $1 - e\cos E$ 가 0에 가까워지며 뉴턴 스텝이 폭주한다. 초기값을 $E_0 = M$ 으로 잡으면 흔히 발산하므로, $E_0 = M + e\sin M$ 같은 개선된 시작값이나 댄비(Danby)식 3차 시작값을 쓴다.
쌍곡선 궤도는 $M_h = e\sinh H - H$ 로 형태가 바뀌고, $\sinh$ 가 지수적으로 커져 초기값이 조금만 나빠도 오버플로가 난다.

$E$ 를 얻으면 나머지는 대수다. $r = a(1 - e\cos E)$ 이고

\tan\frac{\nu}{2} = \sqrt{\frac{1+e}{1-e}}\,\tan\frac{E}{2}

로 진근점각을 뽑는다. 이심률 구간별로 방정식을 갈아 끼우는 게 싫으면 보편 변수(universal variable) 정식화를 쓴다. 타원·포물선·쌍곡선을 스텀프프 함수 $C(z), S(z)$ 하나로 통합해 버리는 방법이고, 범용 궤도 전파 코드는 거의 다 이걸 쓴다.

5. 궤도 전이[편집]

5.1. 호만 전이[편집]

같은 평면 위 반지름 $r_1$ 원궤도에서 $r_2$ 원궤도로 옮기는 가장 유명한 방법. 두 원에 동시에 접하는 타원( $a_t = (r_1+r_2)/2$ )을 타고, 근점과 원점에서 한 번씩 총 두 번 점화한다.

\Delta v_1 = \sqrt{\frac{\mu}{r_1}}\left(\sqrt{\frac{2r_2}{r_1+r_2}} - 1\right), \qquad \Delta v_2 = \sqrt{\frac{\mu}{r_2}}\left(1 - \sqrt{\frac{2r_1}{r_1+r_2}}\right)

전이 시간은 반주기 $t = \pi\sqrt{a_t^3/\mu}$ 다. 고도 200 km 저궤도( $r_1 = 6578$ km)에서 정지궤도( $r_2 = 42164$ km)로 가면 $\Delta v_1 \approx 2.46$ km/s, $\Delta v_2 \approx 1.48$ km/s, 합 약 3.93 km/s이고 걸리는 시간은 약 5.3시간이다. 여기에 경사각을 28.5°에서 0°로 바꾸는 비용이 더 붙으므로, 실제 운용에서는 원점 점화에 평면 변경을 합쳐서 벡터 합으로 절약한다.

5.2. 바이엘립틱 전이[편집]

호만이 항상 최적은 아니다. 아주 먼 중간 지점 $r_b \gg r_2$ 까지 한 번 나갔다가 되돌아오는 바이엘립틱 전이는 점화가 세 번이지만, 목표 궤도가 충분히 멀면 총 $\Delta v$ 가 더 작다. 경계는 고전적으로 알려져 있다.

$r_2/r_1 < 11.94$ — 호만이 언제나 이긴다.
$r_2/r_1 > 15.58$ — 바이엘립틱이 (중간점을 충분히 크게 잡으면) 언제나 이긴다.
그 사이 — 중간점 $r_b$ 를 얼마로 잡느냐에 달렸다.

직관은 오베르트 효과의 역이다. 멀리 나가면 속도가 느려져서 방향을 바꾸는 데 드는 $\Delta v$ 가 싸진다. 대신 전이 시간이 며칠~몇 달로 늘어나므로, 시간이 공짜인 임무에서만 쓸 수 있다.²

5.3. 램버트 문제[편집]

호만은 “출발지와 도착지”만 주어졌을 때의 답이다. 현실의 임무 설계는 날짜까지 정해져 있다 — 위치 $\mathbf r_1$ 에서 $\mathbf r_2$ 로 정확히 $\Delta t$ 만에 가는 궤도를 찾아라. 이것이 램버트 문제이고, 램버트 정리에 따르면 전이 시간은 $a$ , 현의 길이 $c$ , $r_1+r_2$ 세 개에만 의존한다. 라그랑주 정식화나 보편 변수 기반의 반복(배틴, 구딩, 이조 알고리즘 등)으로 풀며, 궤도 역학 라이브러리의 심장부에 해당한다.

램버트 솔버를 출발일·도착일 격자 위에서 수천 번 돌려 소요 에너지 $C_3 = v_\infty^2$ 의 등고선을 그린 것이 그 유명한 포크찹 플롯이고, 발사창은 그 그림에서 골짜기를 찾는 작업이다. 지구-화성 회합주기는 약 780일이라 화성 창은 대략 26개월마다 열리며, 호만급 전이 시간은 편도 8~9개월이다.

6. 원뿔곡선 이어붙이기와 영향권[편집]

행성간 비행은 태양·지구·화성이 동시에 당기는 명백한 다체 문제지만, 실무는 이것을 원뿔곡선 이어붙이기(patched conic)로 쪼갠다. 공간을 천체별 영향권(SOI)으로 나누고, 각 구간에서는 그 천체 하나만 있는 이체 문제로 푼 다음 경계에서 위치·속도를 이어 붙이는 것이다. 영향권 반지름은 라플라스의 추정으로

r_{SOI} \approx a_p \left(\frac{m_p}{M_\odot}\right)^{2/5}

지구는 약 92만 km, 화성은 약 58만 km, 달은 지구 기준 약 6.6만 km다. 지구 SOI 안에서는 쌍곡선 탈출 궤도를 풀고, 나오면 태양 중심 타원으로 갈아타고, 화성 SOI에 들어가면 다시 쌍곡선 접근 궤도로 갈아탄다. 근사임이 명백한데도 오차가 임무 설계 초기 단계에서 감당할 만한 수준이라, 지금도 1차 설계는 전부 이걸로 한다. 정밀 궤적은 이 해를 초기 추정값으로 넣고 전체 중력 모델로 수치적분해 다듬는다.

7. 섭동 — J2가 하는 일[편집]

지구는 자전 때문에 적도가 부푼 회전타원체이고, 그 결과 중력 퍼텐셜에 $J_2 = 1.08263\times10^{-3}$ 항이 붙는다. 다른 모든 섭동(고차 중력항, 대기 항력, 3체, 태양복사압)보다 세 자릿수쯤 크기 때문에, 지구 궤도 위성 설계에서는 사실상 $J_2$ 만 따로 취급한다. 장주기 성분을 평균해서 남는 세속 변화율은 이렇다.

\dot\Omega = -\frac{3}{2}\,n J_2 \left(\frac{R_\oplus}{p}\right)^2\cos i, \qquad \dot\omega = \frac{3}{4}\,n J_2 \left(\frac{R_\oplus}{p}\right)^2\left(5\cos^2 i - 1\right)

이 두 줄이 실무에서 통째로 도구가 된다.

태양동기 궤도. $\dot\Omega$ 가 지구 공전과 같은 하루 약 0.9856°가 되도록 경사각을 맞추면, 궤도면이 태양을 늘 같은 각도로 보게 된다. $\cos i < 0$ 이 필요하므로 역행 궤도여야 하고, 고도 700~800 km에서 $i \approx 98^\circ$ 가 나온다. 지구관측 위성이 죄다 이 궤도를 쓰는 이유 — 매번 같은 지방시에 지나가야 영상 비교가 되기 때문이다.
임계 경사각. $5\cos^2 i - 1 = 0$ , 즉 $i = 63.435^\circ$ 에서는 근점이 안 돈다. 근점을 남반구에 고정해 두고 원점 쪽에서 고위도를 오래 비추는 몰니야 궤도가 정확히 이 각도를 쓴다.

수치적으로 섭동을 다루는 방식은 세 갈래다. 전체 가속도를 그냥 다 넣고 적분하는 코웰법, 이체 해와의 차이만 적분해 오차 누적을 줄이는 엔케법, 궤도 요소의 변화율을 적분하는 매개변수 변화법. 장기 전파에서는 에너지 표류가 문제라 심플렉틱 적분기를 쓰고, 관측을 섞어 상태를 추정하는 궤도 결정에서는 칼만 필터와 최소자승법이 주력이다.

8. 오베르트 효과와 평면 변경[편집]

임무 설계자가 $\Delta v$ 를 아끼는 방식은 대개 두 가지 감각에서 나온다.

오베르트 효과. 같은 $\Delta v$ 를 써도 빠를 때 쓰는 편이 에너지 이득이 크다. 비에너지가 $\varepsilon = v^2/2 - \mu/r$ 이므로 속도 방향으로 $\Delta v$ 를 더하면 에너지 증가는

\Delta\varepsilon = v\,\Delta v + \frac{(\Delta v)^2}{2}

즉 현재 속도 $v$ 에 비례하는 항이 붙는다. 그래서 탈출 점화는 궤도에서 가장 빠른 지점, 곧 근점에서 몰아서 하는 것이 언제나 유리하다. 얕은 궤도에서 여러 번 나눠 태우는 것보다 근점 한 번에 몰아 태우는 쪽이 이기는 이유이며, 반대로 저추력 추진이 근본적으로 손해를 보는 지점이기도 하다.

평면 변경은 비싸다. 속도 $v$ 인 원궤도에서 경사각만 $\Delta i$ 만큼 바꾸려면

\Delta v = 2v\sin\frac{\Delta i}{2}

가 필요하다. 저궤도에서 $v \approx 7.8$ km/s이므로 경사각 30°를 바꾸는 데만 4 km/s가 든다 — 저궤도에서 정지궤도까지 가는 비용보다 크다. 그래서 평면 변경은 (a) 속도가 느린 원점에서, (b) 고도 변경 점화와 벡터 합으로 합쳐서 하는 것이 원칙이고, 애초에 발사장 위도가 목표 경사각을 사실상 결정해 버린다.³ 적도 근처 발사장이 프리미엄을 받는 이유가 여기 있다.

9. 상대 운동과 랑데부[편집]

두 물체가 거의 같은 궤도에 있을 때는 절대 궤도 요소보다 상대 위치로 기술하는 편이 편하다. 기준 원궤도(각속도 $n$ )에 붙은 회전좌표계에서 선형화하면 클로히시-윌트셔 방정식이 나온다.

\ddot x - 3n^2 x - 2n\dot y = 0, \qquad \ddot y + 2n\dot x = 0, \qquad \ddot z + n^2 z = 0

( $x$ 는 반경 방향, $y$ 는 진행 방향, $z$ 는 궤도면 법선.) 여기서 우주 랑데부의 악명 높은 반직관이 그대로 튀어나온다. 진행 방향으로 가속하면 $\dot x$ 가 생겨 궤도가 커지고, 주기가 길어져서 결국 뒤로 처진다. 앞차를 따라잡으려면 감속해야 한다는 그 이야기다. 반대로 면외 성분 $z$ 는 다른 축과 결합하지 않고 단순 조화진동만 하므로, 궤도면이 다른 물체를 만나려면 진동을 억지로 죽여야 한다 — 그게 위에서 본 비싼 평면 변경이다.

CW 방정식은 선형·원궤도 근사라 수백 km를 넘어가면 어긋나고, 이심률이 있으면 시간 가변 계수를 갖는 초-훈-슈트라입 계열로 바꿔야 한다. 실제 접근 유도는 이 해석해를 초기 추정으로 삼고 최적 제어와 칼만 필터 기반 항법을 얹어 굴린다.

10. 이체 밖으로[편집]

제한 삼체 문제. 두 무거운 천체가 원궤도로 돌고 세 번째 물체의 질량이 무시될 만한 경우. 닫힌 해는 없지만 회전좌표계에서 야코비 적분이라는 보존량이 남고, 그 등위면이 힐 영역을 그린다. 평형점 5개(L1~L5)와 그 주위의 헤일로 궤도는 제한 삼체 문제의 결과물이며, 여기서부터는 궤도가 카오스 이론의 언어로 기술된다.
저추력 궤적. 점화를 순간으로 근사할 수 없는 전기추진은 위 공식들이 통째로 무의미해진다. 수천 바퀴를 나선으로 도는 궤적은 최적 제어 문제로 풀어야 한다.
발사 구간. 지상에서 궤도 속도까지 올라오는 구간은 중력장 문제가 아니라 대기·추력 문제이므로 중력 선회가 담당한다.

11. 관련 문서[편집]

수성 근일점의 이상 이동 43″/세기가 뉴턴 중력으로 설명이 안 된 것이 일반상대성이론의 초기 증거가 되었다. 즉 이심률 벡터가 도는 속도를 정밀하게 재는 것만으로 중력 이론을 검증할 수 있었다는 뜻인데, 궤도 역학이 물리학의 시험대였던 시절 이야기다. ↩
그리고 시간이 공짜인 임무는 사실상 없다. 바이엘립틱이 이론적으로 유명한 것에 비해 실제 사용 사례가 드문 이유다. 다만 저추력 나선 궤적은 사실상 무한 다중 임펄스라서 이 논리의 극한에 해당한다. ↩
발사장 위도 $\phi$ 보다 작은 경사각으로는 직접 못 올린다. 궤도면이 지구 중심을 지나야 하므로 발사 지점이 궤도면 위에 있어야 하기 때문이다. 케네디 우주센터가 28.5°, 바이코누르가 45.6°, 기아나 우주센터가 5.2°. 정지궤도 위성 사업에서 기아나가 유독 인기 있는 데는 낭만이 아니라 산수가 있다. ↩