13.14 오일러 회전 운동 방정식의 유도

1. 개요

오일러 회전 운동 방정식(Euler’s rotation equations)은 강체의 회전 동역학을 본체 좌표계의 주축에서 기술하는 비선형 미분 방정식이다. 1758년 레온하르트 오일러(Leonhard Euler)에 의해 제시되었으며, 강체 회전 운동의 가장 기본적이고 우아한 표현 중 하나이다. 본 절에서는 오일러 회전 운동 방정식의 체계적인 유도를 다룬다.

2. 유도의 출발점

2.1 각운동량과 토크의 관계

관성 좌표계에서 강체에 대한 회전 운동 방정식은 다음과 같다.

$\frac{d\mathbf{L}}{dt}\bigg|_{\text{inertial}} = \boldsymbol{\tau}$

여기서 $\mathbf{L}$ 은 각운동량, $\boldsymbol{\tau}$ 는 외부 토크이다.

2.2 본체 좌표계의 도입

본체 좌표계에서는 관성 텐서가 시간에 무관하므로 분석이 단순해진다. 그러나 본체 좌표계가 회전하므로 시간 미분에 추가 항이 필요하다.

3. 시간 미분의 변환

3.1 회전 좌표계의 미분

벡터 $\mathbf{A}$ 의 시간 미분을 두 좌표계에서 비교하면

$\frac{d\mathbf{A}}{dt}\bigg|_{\text{inertial}} = \frac{d\mathbf{A}}{dt}\bigg|_{\text{body}} + \boldsymbol{\omega} \times \mathbf{A}$

여기서 $\boldsymbol{\omega}$ 는 본체 좌표계의 각속도이다.

3.2 각운동량에 적용

각운동량 $\mathbf{L}$ 에 적용하면

$\frac{d\mathbf{L}}{dt}\bigg|_{\text{inertial}} = \frac{d\mathbf{L}}{dt}\bigg|_{\text{body}} + \boldsymbol{\omega} \times \mathbf{L}$

4. 운동 방정식의 변환

4.1 본체 좌표계에서

운동 방정식 $\dot{\mathbf{L}} = \boldsymbol{\tau}$ 에 위의 변환을 대입하면

$\frac{d\mathbf{L}}{dt}\bigg|_{\text{body}} + \boldsymbol{\omega} \times \mathbf{L} = \boldsymbol{\tau}$

4.2 관성 텐서의 도입

본체 좌표계에서 $\mathbf{L} = \mathbf{I}_{\text{body}}\boldsymbol{\omega}$ 이고 $\mathbf{I}_{\text{body}}$ 는 시간에 무관하므로

$\mathbf{I}_{\text{body}}\dot{\boldsymbol{\omega}} + \boldsymbol{\omega} \times (\mathbf{I}_{\text{body}}\boldsymbol{\omega}) = \boldsymbol{\tau}$

4.3 의미

이는 본체 좌표계에서의 회전 운동 방정식이다. 자이로스코프 항 $\boldsymbol{\omega} \times (\mathbf{I}\boldsymbol{\omega})$ 가 운동 방정식의 비선형성을 도입한다.

5. 주축 좌표계로의 단순화

5.1 주축의 선택

본체 좌표계를 주축 좌표계로 선택하면 관성 텐서가 대각 행렬이 된다.

$\mathbf{I}_{\text{principal}} = \begin{bmatrix}I_1 & 0 & 0 \\ 0 & I_2 & 0 \\ 0 & 0 & I_3\end{bmatrix}$

5.2 성분 형태

각속도 벡터 $\boldsymbol{\omega} = (\omega_1, \omega_2, \omega_3)$ 와 토크 $\boldsymbol{\tau} = (\tau_1, \tau_2, \tau_3)$ 를 사용하여 운동 방정식의 성분을 계산한다.

5.2.1 자이로스코프 항의 계산

$\mathbf{I}\boldsymbol{\omega} = (I_1\omega_1, I_2\omega_2, I_3\omega_3)$

$\boldsymbol{\omega} \times (\mathbf{I}\boldsymbol{\omega}) = \begin{vmatrix}\hat{\mathbf{e}}_1 & \hat{\mathbf{e}}_2 & \hat{\mathbf{e}}_3 \\ \omega_1 & \omega_2 & \omega_3 \\ I_1\omega_1 & I_2\omega_2 & I_3\omega_3\end{vmatrix}$

전개하면

$\boldsymbol{\omega} \times (\mathbf{I}\boldsymbol{\omega}) = ((I_3 - I_2)\omega_2\omega_3, (I_1 - I_3)\omega_3\omega_1, (I_2 - I_1)\omega_1\omega_2)$

5.2.2 운동 방정식의 성분

이를 사용하여 각 성분의 운동 방정식을 작성한다.

$I_1\dot\omega_1 + (I_3 - I_2)\omega_2\omega_3 = \tau_1$

$I_2\dot\omega_2 + (I_1 - I_3)\omega_3\omega_1 = \tau_2$

$I_3\dot\omega_3 + (I_2 - I_1)\omega_1\omega_2 = \tau_3$

이것이 오일러 회전 운동 방정식이다.

6. 오일러 방정식의 형태

6.1 일반 형태

$I_i\dot\omega_i + (I_k - I_j)\omega_j\omega_k = \tau_i, \quad (i, j, k) \in \{(1,2,3), (2,3,1), (3,1,2)\}$

6.2 자유 강체

외부 토크가 없는 경우( $\tau_i = 0$ ), 오일러 방정식은 다음과 같이 단순화된다.

$I_1\dot\omega_1 = (I_2 - I_3)\omega_2\omega_3$

$I_2\dot\omega_2 = (I_3 - I_1)\omega_3\omega_1$

$I_3\dot\omega_3 = (I_1 - I_2)\omega_1\omega_2$

이는 자유 강체 회전의 비선형 동역학을 결정한다.

7. 보존량

7.1 운동 에너지

오일러 방정식의 첫 번째 적분은 운동 에너지의 보존이다.

$T = \frac{1}{2}(I_1\omega_1^2 + I_2\omega_2^2 + I_3\omega_3^2) = \text{const}$

자유 강체에서 토크가 없으면 운동 에너지가 보존된다.

7.2 각운동량의 크기

각운동량의 크기도 보존된다.

$L^2 = I_1^2\omega_1^2 + I_2^2\omega_2^2 + I_3^2\omega_3^2 = \text{const}$

7.3 두 보존량과 운동의 분류

이 두 보존량이 오일러 방정식의 해를 분류하는 데 사용된다.

8. 자유 강체의 회전 운동 분석

8.1 안정 회전

가장 큰 또는 가장 작은 주관성 모멘트의 주축 주위의 회전은 안정적이다. 작은 교란이 가해져도 회전이 그 축 근처에 머문다.

8.2 불안정 회전

중간 주관성 모멘트의 주축 주위의 회전은 불안정하다. 작은 교란이 가해지면 회전 축이 다른 안정 축으로 전환된다. 이를 테니스 라켓 정리라 한다.

8.3 분석

오일러 방정식을 작은 교란에 대해 선형화하면 안정성을 분석할 수 있다.

9. 오일러 방정식의 해

9.1 적분 가능성

자유 강체의 오일러 방정식은 야코비 타원 함수를 사용하여 해석적으로 풀 수 있다. 이는 적분 가능 시스템의 한 예이다.

9.2 일반 해

해는 운동 에너지와 각운동량의 두 보존량으로 매개변수화된다. 운동의 형태가 이 두 양의 비율에 의존한다.

9.3 수치적 해

대부분의 응용에서 수치적 적분이 사용된다. 룽게-쿠타 방법 등이 적용된다.

10. 응용 예시: 균일 원기둥의 자유 회전

균일 원기둥은 두 주관성 모멘트가 같다( $I_1 = I_2 \neq I_3$ ). 이 경우 오일러 방정식은 다음과 같이 단순화된다.

$I_1\dot\omega_1 + (I_3 - I_1)\omega_2\omega_3 = 0$

$I_1\dot\omega_2 + (I_1 - I_3)\omega_3\omega_1 = 0$

$I_3\dot\omega_3 = 0$

세 번째 방정식에서 $\omega_3$ 가 일정함을 알 수 있다. 처음 두 방정식이 결합되어 $\omega_1$ 과 $\omega_2$ 가 주기 운동을 한다.

11. 응용 예시: 자유 강체의 자이로스코프 효과

자이로스코프(빠르게 회전하는 강체)는 외부 토크에 대해 흥미로운 응답을 보인다. 오일러 방정식이 이러한 응답을 정확히 기술한다.

12. 응용 예시: 위성의 자세 동역학

지구 궤도의 위성은 약한 토크(대기 항력, 중력 경사 등)의 영향을 받는다. 오일러 방정식이 자세 동역학의 분석에 사용된다.

13. 응용 예시: 매니퓰레이터의 링크

매니퓰레이터의 각 링크의 회전 동역학은 오일러 방정식의 형태로 기술된다. 다체 시스템의 복잡한 결합이 추가된다.

14. 응용 예시: 무인 항공기

드론의 자세 동역학은 오일러 방정식으로 기술된다. 모터의 토크와 공기력의 영향이 외부 토크로 작용한다.

15. 본 절의 의의

본 절은 오일러 회전 운동 방정식의 유도를 다루었다. 이는 강체 회전 동역학의 가장 기본적이고 우아한 표현이며, 후속 절과 다양한 응용의 토대이다. 매니퓰레이터, 우주선, 무인 항공기 등 모든 회전 시스템의 분석에서 오일러 방정식이 사용된다.

16. 학습 권장사항

회전 좌표계의 시간 미분을 이해한다.
본체 좌표계와 관성 좌표계의 차이를 인식한다.
오일러 방정식의 유도를 따라가 본다.
보존량과 운동의 관계를 학습한다.
단순한 강체에 대해 오일러 방정식을 풀어 본다.

17. 참고 문헌

Euler, L. (1758). “Du mouvement de rotation des corps solides autour d’un axe variable.” Mémoires de l’Académie des Sciences de Berlin.
Goldstein, H., Poole, C., & Safko, J. (2002). Classical Mechanics (3rd ed.). Addison-Wesley.
Symon, K. R. (1971). Mechanics (3rd ed.). Addison-Wesley.
Greenwood, D. T. (2003). Advanced Dynamics. Cambridge University Press.
Featherstone, R. (2008). Rigid Body Dynamics Algorithms. Springer.

version: 1.0