13.8 회전축 정리와 관성 모멘트 변환

1. 회전축 정리의 개요

회전축 정리(rotation axis theorem) 또는 좌표계 회전에 의한 관성 텐서의 변환은 강체의 관성 텐서가 좌표계의 회전에 대해 어떻게 변환되는지를 기술하는 정리이다. 평행축 정리가 좌표계의 병진에 대한 변환을 다루는 반면, 회전축 정리는 좌표계의 회전에 대한 변환을 다룬다. 두 정리를 결합하면 강체의 관성 텐서를 임의의 좌표계로 변환할 수 있으며, 이는 매니퓰레이터의 동역학 분석, 다체 시스템의 모델링, 자세 의존 동역학 분석 등에서 핵심적인 도구이다.

2. 회전 좌표 변환

2.1 좌표계 회전

두 좌표계 A와 B가 같은 원점을 가지지만 회전 행렬 \mathbf{R}로 관련되어 있다고 하자. 즉, 좌표계 B의 좌표가 좌표계 A의 좌표로 다음과 같이 변환된다.

\mathbf{r}_A = \mathbf{R}\mathbf{r}_B

여기서 \mathbf{R}은 직교 회전 행렬이다.

2.2 변환의 의미

이 변환은 같은 물리적 점이 두 좌표계에서 다른 좌표를 가짐을 의미한다.

3. 관성 텐서의 회전 변환

3.1 정리의 진술

회전 변환 정리: 좌표계 B에서의 관성 텐서를 \mathbf{I}_B라 하면, 좌표계 A에서의 관성 텐서 \mathbf{I}_A는 다음과 같다.

\mathbf{I}_A = \mathbf{R}\mathbf{I}_B\mathbf{R}^T

여기서 \mathbf{R}은 좌표계 B에서 좌표계 A로의 회전 행렬이다.

3.2 의미

관성 텐서는 회전 변환에 대해 텐서로서 변환된다. 이는 일반적인 2차 텐서의 변환 법칙이다.

3.3 역변환

역변환은 다음과 같다.

\mathbf{I}_B = \mathbf{R}^T\mathbf{I}_A\mathbf{R}

직교 회전 행렬에서 \mathbf{R}^{-1} = \mathbf{R}^T이기 때문이다.

4. 회전 변환 정리의 증명

4.1 설정

좌표계 A에서의 관성 텐서는

\mathbf{I}_A = \int_B\rho(\mathbf{r}_A^T\mathbf{r}_A\mathbf{I}_3 - \mathbf{r}_A\mathbf{r}_A^T)\,dV

4.2 변환 적용

\mathbf{r}_A = \mathbf{R}\mathbf{r}_B를 대입하면

\mathbf{I}_A = \int_B\rho((\mathbf{R}\mathbf{r}_B)^T(\mathbf{R}\mathbf{r}_B)\mathbf{I}_3 - (\mathbf{R}\mathbf{r}_B)(\mathbf{R}\mathbf{r}_B)^T)\,dV

4.3 단순화

직교 회전 행렬에서 \mathbf{R}^T\mathbf{R} = \mathbf{I}_3이므로

(\mathbf{R}\mathbf{r}_B)^T(\mathbf{R}\mathbf{r}_B) = \mathbf{r}_B^T\mathbf{R}^T\mathbf{R}\mathbf{r}_B = \mathbf{r}_B^T\mathbf{r}_B

또한

(\mathbf{R}\mathbf{r}_B)(\mathbf{R}\mathbf{r}_B)^T = \mathbf{R}\mathbf{r}_B\mathbf{r}_B^T\mathbf{R}^T

4.4 결과

따라서

\mathbf{I}_A = \int_B\rho(\mathbf{r}_B^T\mathbf{r}_B\mathbf{I}_3 - \mathbf{R}\mathbf{r}_B\mathbf{r}_B^T\mathbf{R}^T)\,dV

= \mathbf{R}\!\left[\int_B\rho(\mathbf{r}_B^T\mathbf{r}_B\mathbf{I}_3 - \mathbf{r}_B\mathbf{r}_B^T)\,dV\right]\mathbf{R}^T

= \mathbf{R}\mathbf{I}_B\mathbf{R}^T

여기서 \mathbf{R}\mathbf{I}_3\mathbf{R}^T = \mathbf{R}\mathbf{R}^T = \mathbf{I}_3의 사실을 사용하였다.

5. 평행축 정리와 회전 변환의 결합

5.1 일반적 변환

임의의 두 좌표계 사이의 관성 텐서 변환은 평행축 정리와 회전 변환의 결합으로 수행된다. 두 좌표계가 위치와 방향이 모두 다른 경우 다음의 단계로 변환한다.

1. 한 좌표계에서 그 좌표계의 원점에 대한 관성 텐서를 알고 있다.
2. 평행축 정리로 질량 중심에 대한 관성 텐서를 계산한다.
3. 회전 변환으로 다른 좌표계의 방향으로 변환한다.
4. 평행축 정리로 다른 좌표계의 원점에 대한 관성 텐서를 계산한다.

5.2 응용

매니퓰레이터의 동역학 분석에서 각 링크의 관성 텐서를 다양한 기준점과 방향에서 표현해야 한다. 이러한 변환이 평행축 정리와 회전 변환으로 수행된다.

6. 임의의 회전축에 대한 관성 모멘트

6.1 관성 텐서로부터의 계산

임의의 회전축 \hat{\mathbf{n}} (단위 벡터)에 대한 관성 모멘트는 다음과 같이 계산된다.

I_{\hat{\mathbf{n}}} = \hat{\mathbf{n}}^T\mathbf{I}\hat{\mathbf{n}}

여기서 \mathbf{I}는 회전축의 한 점을 원점으로 하는 좌표계에서의 관성 텐서이다.

6.2 명시적 형태

회전축의 단위 벡터를 \hat{\mathbf{n}} = (n_x, n_y, n_z)라 하면

I_{\hat{\mathbf{n}}} = I_{xx}n_x^2 + I_{yy}n_y^2 + I_{zz}n_z^2 + 2I_{xy}n_xn_y + 2I_{xz}n_xn_z + 2I_{yz}n_yn_z

6.3 의미

이 식은 관성 텐서가 임의의 회전축에 대한 관성 모멘트의 정보를 모두 담고 있음을 보여준다. 6개의 독립 성분으로 무한히 많은 회전축에 대한 관성 모멘트가 결정된다.

7. 회전 변환의 활용

7.1 좌표계 선택의 자유

관성 텐서의 회전 변환은 강체의 분석에 가장 편리한 좌표계를 선택할 수 있게 한다. 일반적으로 주축을 따른 좌표계가 선호된다(다음 절 참조).

7.2 운동 방정식의 단순화

적절한 좌표계 선택은 운동 방정식을 크게 단순화한다. 관성 텐서가 대각이면 분석이 쉽다.

7.3 시간에 따른 변화

본체 좌표계는 강체와 함께 회전하므로 강체의 관성 텐서는 본체 좌표계에서 시간에 무관하다. 그러나 세계 좌표계에서는 시간에 따라 변한다.

\mathbf{I}_{\text{world}}(t) = \mathbf{R}(t)\mathbf{I}_{\text{body}}\mathbf{R}^T(t)

여기서 \mathbf{R}(t)는 본체에서 세계로의 회전 행렬이다.

8. 매니퓰레이터에서의 응용

8.1 링크 좌표계

매니퓰레이터의 각 링크에는 본체 좌표계가 부착된다. 이 좌표계는 일반적으로 관절 위치와 방향에 의해 결정된다.

8.2 관성 텐서의 표현

각 링크의 관성 텐서는 본체 좌표계에서 시간에 무관하다. 동역학 분석에서 이를 다른 좌표계로 변환한다.

8.3 자세 의존 동역학

매니퓰레이터의 자세가 변하면 각 링크의 자세도 변한다. 세계 좌표계에서의 관성 텐서는 자세에 의존하며, 회전 변환으로 계산된다.

9. 동역학에의 영향

9.1 운동 방정식

매니퓰레이터의 운동 방정식에서 관성 행렬 \mathbf{M}(\mathbf{q})는 모든 링크의 관성 텐서가 적절히 변환되어 결합된 결과이다. 자세 \mathbf{q}에 의존한다.

9.2 관성 행렬의 계산

매니퓰레이터의 관성 행렬을 계산할 때 각 링크의 관성 텐서를 회전 변환으로 변환한 후 결합한다.

9.3 효율성

효율적인 알고리즘(예: Featherstone의 articulated body algorithm)은 이러한 변환을 효율적으로 수행한다.

10. 응용 예시: 회전된 직육면체

직육면체가 본체 좌표계에서 다음의 관성 텐서를 가진다.

\mathbf{I}_{\text{body}} = \frac{M}{12}\begin{bmatrix}b^2 + c^2 & 0 & 0 \\ 0 & a^2 + c^2 & 0 \\ 0 & 0 & a^2 + b^2\end{bmatrix}

본체가 z축 주위로 \theta만큼 회전하면 회전 행렬은

\mathbf{R}_z(\theta) = \begin{bmatrix}\cos\theta & -\sin\theta & 0 \\ \sin\theta & \cos\theta & 0 \\ 0 & 0 & 1\end{bmatrix}

세계 좌표계에서의 관성 텐서는 \mathbf{R}_z\mathbf{I}_{\text{body}}\mathbf{R}_z^T로 계산된다.

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

매니퓰레이터의 두 번째 링크의 관성 텐서가 본체 좌표계에서 알려져 있다. 첫 번째 관절이 회전함에 따라 두 번째 링크의 자세가 변한다. 세계 좌표계에서의 관성 텐서는 회전 변환으로 시간에 따라 변화한다.

12. 응용 예시: 우주선의 자세

우주선의 자세가 변할 때 본체의 관성 텐서는 본체 좌표계에서 일정하다. 그러나 관성 좌표계에서는 회전 변환으로 변환된다. 이는 자세 동역학 분석에 중요하다.

13. 응용 예시: 인간형 로봇

인간형 로봇의 각 부위가 다양한 자세를 가질 때, 전체 시스템의 관성 텐서는 각 부위의 관성 텐서를 회전 변환과 평행축 정리로 결합하여 계산된다.

14. 회전 변환의 수치 계산

14.1 행렬 곱

회전 변환은 단순한 행렬 곱이다.

def rotate_inertia(I_body, R):
    return R @ I_body @ R.T

이는 효율적으로 계산된다.

14.2 자코비안

매니퓰레이터의 동역학 분석에서 관성 텐서의 자세에 대한 자코비안이 필요할 수 있다. 회전 변환의 미분이 사용된다.

15. 본 절의 의의

본 절은 회전축 정리와 관성 텐서의 회전 변환을 다루었다. 평행축 정리와 함께 이 정리는 강체의 관성 텐서를 임의의 좌표계로 변환하는 완전한 도구를 제공한다. 매니퓰레이터, 우주선, 인간형 로봇 등 다양한 응용에서 이러한 변환이 일상적으로 사용된다.

16. 학습 권장사항

• 회전 변환 정리의 형태를 이해한다.
• 증명을 따라가 본다.
• 임의의 회전축에 대한 관성 모멘트를 계산해 본다.
• 평행축 정리와의 결합을 익힌다.
• 매니퓰레이터의 동역학에 적용해 본다.

17. 참고 문헌

• 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