13.12 강체의 회전 운동 방정식

1. 회전 운동 방정식의 개요

강체의 회전 운동 방정식(rotational equations of motion)은 외부에서 가해지는 토크가 강체의 회전 운동을 어떻게 변화시키는지를 기술하는 미분 방정식이다. 병진 운동에서의 뉴턴 제2법칙( $\mathbf{F} = m\mathbf{a}$ )의 회전 형태이며, 직선 운동에 비해 본질적으로 더 복잡한 구조를 가진다. 이 복잡성은 관성 텐서의 텐서 성질과 회전 동역학의 비선형성에서 비롯된다. 본 절에서는 강체의 회전 운동 방정식의 유도와 다양한 형태를 다룬다.

2. 각운동량의 정의

2.1 점 질량의 각운동량

기준점에 대한 점 질량의 각운동량은 다음과 같이 정의된다.

$\mathbf{L} = \mathbf{r} \times \mathbf{p} = m\mathbf{r} \times \mathbf{v}$

여기서

$\mathbf{r}$ : 기준점에서 점 질량까지의 위치 벡터
$\mathbf{p} = m\mathbf{v}$ : 점 질량의 운동량

2.2 강체의 각운동량

강체의 각운동량은 모든 점 질량의 각운동량의 합이다.

$\mathbf{L} = \int_B\rho(\mathbf{r})\mathbf{r} \times \mathbf{v}(\mathbf{r})\,dV$

2.3 관성 텐서를 사용한 표현

강체가 한 점 주위로 회전할 때, 그 점에 대한 각운동량은 관성 텐서와 각속도의 곱으로 표현된다.

$\mathbf{L} = \mathbf{I}\boldsymbol{\omega}$

여기서 $\mathbf{I}$ 는 그 점에 대한 관성 텐서이다.

3. 토크의 정의

3.1 점 질량에 작용하는 토크

기준점에 대한 점 질량에 작용하는 토크는 다음과 같이 정의된다.

$\boldsymbol{\tau} = \mathbf{r} \times \mathbf{F}$

여기서 $\mathbf{F}$ 는 점 질량에 작용하는 힘이다.

3.2 강체에 작용하는 토크

강체에 작용하는 외부 토크의 합은 다음과 같다.

$\boldsymbol{\tau}_{\text{ext}} = \sum_i\mathbf{r}_i \times \mathbf{F}_i^{\text{ext}}$

여기서 $\mathbf{r}_i$ 는 기준점에서 $i$ 번째 점 질량까지의 위치이고, $\mathbf{F}_i^{\text{ext}}$ 는 그 점에 작용하는 외력이다.

4. 각운동량과 토크의 관계

4.1 점 질량

점 질량에 대해 다음의 관계가 성립한다.

$\frac{d\mathbf{L}}{dt} = \boldsymbol{\tau}$

이는 토크가 각운동량의 시간 변화율과 같음을 의미한다. 이는 회전 운동의 뉴턴 제2법칙이다.

4.2 강체

강체에 대해서도 같은 관계가 성립한다.

$\frac{d\mathbf{L}}{dt} = \boldsymbol{\tau}_{\text{ext}}$

여기서 모든 양은 같은 기준점에 대해 측정된다.

4.3 증명 개요

내력 쌍이 토크의 합에서 상쇄됨을 보임으로써 증명된다. 내력은 두 점 질량을 잇는 직선을 따라 작용하므로 그 토크의 합이 0이다.

5. 운동 방정식의 일반 형태

5.1 관성 좌표계

관성 좌표계(고정된 좌표계)에서 회전 운동 방정식은 다음과 같다.

$\boldsymbol{\tau} = \frac{d\mathbf{L}}{dt} = \frac{d}{dt}(\mathbf{I}\boldsymbol{\omega})$

여기서 관성 텐서 $\mathbf{I}$ 는 일반적으로 시간에 의존한다(강체가 회전하면서 텐서의 표현이 변하기 때문).

5.2 미분의 전개

미분을 전개하면

$\boldsymbol{\tau} = \dot{\mathbf{I}}\boldsymbol{\omega} + \mathbf{I}\dot{\boldsymbol{\omega}}$

여기서 $\dot{\mathbf{I}}$ 는 시간 변화율이다. 이 항이 회전 동역학의 비선형성을 도입한다.

5.3 본체 좌표계

본체에 고정된 좌표계에서는 관성 텐서가 시간에 무관하다.

$\frac{d\mathbf{I}_{\text{body}}}{dt} = \mathbf{0}$

그러나 본체 좌표계는 비관성 좌표계이므로 추가적인 항이 등장한다.

6. 본체 좌표계에서의 운동 방정식

6.1 유도

본체 좌표계에서 관성 텐서가 일정하지만, 본체 좌표계의 회전으로 인해 추가 항이 등장한다. 결과는 다음과 같다.

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

여기서 $\boldsymbol{\omega}$ 는 본체 좌표계에서의 각속도이고, $\boldsymbol{\tau}_{\text{body}}$ 는 본체 좌표계에서의 토크이다.

6.2 자이로스코프 항

$\boldsymbol{\omega} \times (\mathbf{I}_{\text{body}}\boldsymbol{\omega})$ 항은 자이로스코프 효과를 나타낸다. 이는 각운동량이 일정하지 않을 때 추가로 필요한 토크이다.

6.3 비선형성

이 항의 존재가 회전 동역학을 비선형으로 만든다. 직선 운동의 단순한 $\mathbf{F} = m\mathbf{a}$ 에 비해 회전 운동은 본질적으로 더 복잡하다.

7. 오일러 회전 운동 방정식

7.1 주축 좌표계

본체 좌표계를 주축 좌표계로 선택하면 관성 텐서가 대각이 되어 운동 방정식이 단순화된다.

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

7.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$

여기서 $\omega_i$ 와 $\tau_i$ 는 각각 $i$ 번째 주축 방향의 각속도와 토크 성분이다.

7.3 자유 강체

외부 토크가 없는 자유 강체( $\boldsymbol{\tau} = \mathbf{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$

이는 자유 강체 회전의 비선형 미분 방정식이다.

8. 회전 운동 방정식의 해석

8.1 직선 운동과의 비교

직선 운동	회전 운동
$\mathbf{F} = m\mathbf{a}$	$\boldsymbol{\tau} = \mathbf{I}\boldsymbol{\alpha} + \boldsymbol{\omega} \times \mathbf{I}\boldsymbol{\omega}$
선형	비선형
단순	복잡
결합 항 없음	자이로스코프 결합 항 있음

8.2 비선형성의 원인

회전 동역학의 비선형성은 다음의 요인에서 비롯된다.

관성 텐서가 텐서이므로 각속도와 각운동량이 일반적으로 평행하지 않다.
본체 좌표계가 회전하므로 추가적인 코리올리 효과가 등장한다.
회전군 $SO(3)$ 이 비가환이므로 회전의 합성이 단순한 합산이 아니다.

9. 응용 예시: 자유 회전 강체

외부 토크 없이 자유롭게 회전하는 강체를 고려하자. 오일러 방정식이 적용된다. 운동의 분석은 다음의 보존량을 사용한다.

각운동량 $\mathbf{L}$ (방향과 크기 모두)
운동 에너지 $T = \frac{1}{2}\boldsymbol{\omega}^T\mathbf{I}\boldsymbol{\omega}$

9.1 안정 회전

가장 큰 또는 가장 작은 주관성 모멘트의 주축 주위의 회전은 안정적이다. 중간 주관성 모멘트의 축 주위의 회전은 불안정하다(테니스 라켓 정리).

10. 응용 예시: 자이로스코프

자이로스코프는 회전하는 강체의 흥미로운 동역학을 보여준다. 빠르게 회전하는 자이로스코프는 외부 토크에 대해 흥미로운 응답을 보인다.

10.1 진동(Precession)

회전하는 자이로스코프에 외부 토크가 작용하면 회전 축이 진동(precess)한다. 이 운동은 다음과 같이 기술된다.

$\boldsymbol{\omega}_{\text{prec}} = \frac{\boldsymbol{\tau}}{L}$

여기서 $L$ 은 자이로스코프의 각운동량의 크기이다.

10.2 응용

자이로스코프는 항법 시스템, 자세 안정화 등에 사용된다.

11. 응용 예시: 회전하는 우주선

우주선의 자세 동역학은 오일러 방정식으로 기술된다. 외부 토크(태양 광압, 자기 토크 등)와 내부 토크(운동량 휠, 추력기)가 자세를 변화시킨다.

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

매니퓰레이터의 각 링크는 강체로 모델링되며, 회전 운동 방정식이 적용된다. 관절 토크와 다른 링크로부터의 반력이 외부 토크이다.

13. 응용 예시: 떨어지는 고양이

떨어지는 고양이가 자세를 바꾸어 항상 발로 착지하는 현상은 강체 회전 동역학의 흥미로운 예이다. 각운동량 보존 하에서 자세 변화가 가능함을 보여준다.

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

드론의 자세 동역학은 회전 운동 방정식으로 기술된다. 모터의 회전이 토크를 생성하여 자세를 변화시킨다.

15. 회전 운동 에너지

15.1 정의

강체의 회전 운동 에너지는 다음과 같다.

$T_{\text{rot}} = \frac{1}{2}\boldsymbol{\omega}^T\mathbf{I}\boldsymbol{\omega}$

이는 직선 운동 에너지 $T = \frac{1}{2}m\mathbf{v}^T\mathbf{v}$ 의 회전 형태이다.

15.2 보존

외부 토크가 없는 경우(또는 토크가 회전과 수직인 경우), 회전 운동 에너지가 보존된다.

15.3 일

토크가 회전 운동에 한 일은 다음과 같다.

$W = \int\boldsymbol{\tau} \cdot \boldsymbol{\omega}\,dt$

이는 토크와 각속도의 내적의 시간 적분이다.

16. 각운동량 보존

16.1 보존 법칙

외부 토크의 합이 0이면 강체의 각운동량이 보존된다.

$\frac{d\mathbf{L}}{dt} = \mathbf{0} \implies \mathbf{L} = \text{const}$

16.2 응용

떨어지는 고양이의 자세 변화
빙판 위 회전하는 피겨 스케이터
자유 강체의 회전 운동
우주선의 자세 분석

17. 임펄스와 충격 토크

17.1 충격 토크

매우 짧은 시간에 작용하는 큰 토크를 충격 토크라 한다. 그 시간 적분이 각 임펄스(angular impulse)이다.

$\mathbf{J}_{\text{ang}} = \int\boldsymbol{\tau}\,dt$

17.2 각운동량 변화

각 임펄스는 각운동량의 변화와 같다.

$\mathbf{J}_{\text{ang}} = \Delta\mathbf{L}$

17.3 응용

충돌 분석에서 사용된다. 충격 시간 동안의 다른 토크는 무시 가능하다고 가정한다.

18. 회전 운동 방정식의 수치 적분

18.1 도전

회전 운동 방정식은 비선형이므로 일반적으로 수치적으로 풀어야 한다. 단순한 오일러 적분은 부정확할 수 있다.

18.2 방법

다양한 수치 적분 방법이 사용된다.

룽게-쿠타 방법 (RK4 등)
심플렉틱 적분기
변분 적분기

특히 자세는 다양체 위에서 표현되므로(SO(3) 또는 단위 쿼터니언), 다양체 구조를 고려한 적분이 필요하다.

18.3 응용

물리 시뮬레이터(Gazebo, MuJoCo, PyBullet 등)는 효율적인 적분 알고리즘을 사용한다.

19. 매니퓰레이터의 회전 운동

19.1 다체 시스템

매니퓰레이터는 여러 강체가 관절로 연결된 다체 시스템이다. 각 링크에 대한 회전 운동 방정식이 결합되어 전체 동역학을 형성한다.

19.2 관성 행렬

다체 시스템의 관성 효과는 관성 행렬 $\mathbf{M}(\mathbf{q})$ 로 표현된다. 이는 단일 강체의 관성 텐서의 다체 일반화이다.

19.3 라그랑주와 뉴턴-오일러

다체 시스템의 동역학은 라그랑주 방법이나 뉴턴-오일러 방법으로 유도된다. 두 방법 모두 본 절에서 다룬 단일 강체의 운동 방정식에 기반한다.

20. 시뮬레이션

20.1 강체 시뮬레이션

강체 시뮬레이션 엔진은 회전 운동 방정식을 수치적으로 풀어 강체의 운동을 시뮬레이션한다. 정확하고 안정적인 알고리즘이 필요하다.

20.2 응용

게임의 물리 엔진
로봇 시뮬레이터
영화의 시각 효과
공학 분석

21. 학습 권장사항

각운동량과 토크의 관계를 이해한다.
본체 좌표계와 관성 좌표계의 차이를 인식한다.
오일러 회전 운동 방정식을 학습한다.
회전 동역학의 비선형성을 인식한다.
자유 강체 회전의 분석을 연습한다.

22. 본 절의 의의

본 절은 강체의 회전 운동 방정식을 다루었다. 이는 강체 동역학의 핵심이며, 후속 절에서 다룰 다양한 응용의 토대이다. 회전 운동 방정식의 비선형성은 강체 동역학을 직선 운동보다 더 복잡하게 만들지만, 동시에 흥미로운 현상(자이로스코프 효과, 떨림 등)을 발생시킨다.

23. 참고 문헌

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.
Murray, R. M., Li, Z., & Sastry, S. S. (1994). A Mathematical Introduction to Robotic Manipulation. CRC Press.

version: 1.0