6.37 주축 정리와 관성 텐서 대각화

1. 관성 텐서의 정의

강체의 회전 운동을 기술하기 위해서는 질량 분포의 회전 관성을 정량화하는 텐서가 필요하다. 강체의 어떤 기준점 O에 대한 관성 텐서(inertia tensor)는 다음과 같이 정의된다.

\mathbf{I}_O = \int_{\mathcal{B}} \left( \|\mathbf{r}\|^2 \mathbb{I}_3 - \mathbf{r}\mathbf{r}^\top \right) dm

여기서 \mathbf{r} \in \mathbb{R}^3는 기준점 O로부터 질량 요소까지의 위치 벡터, \mathbb{I}_3는 3 \times 3 단위 행렬, \mathcal{B}는 강체의 점유 영역이다. 성분으로 표기하면 다음과 같다.

\mathbf{I}_O = \begin{pmatrix} I_{xx} & -I_{xy} & -I_{xz} \\ -I_{xy} & I_{yy} & -I_{yz} \\ -I_{xz} & -I_{yz} & I_{zz} \end{pmatrix}

대각 성분은 다음과 같이 정의된다.

I_{xx} = \int (y^2 + z^2) dm, \quad I_{yy} = \int (x^2 + z^2) dm, \quad I_{zz} = \int (x^2 + y^2) dm

비대각 성분은 관성곱(product of inertia)이라 하며 다음과 같이 정의된다.

I_{xy} = \int xy \, dm, \quad I_{yz} = \int yz \, dm, \quad I_{xz} = \int xz \, dm

2. 관성 텐서의 기본 성질

관성 텐서는 다음의 핵심 성질들을 가진다.

• 대칭성: 정의로부터 \mathbf{I}_O^\top = \mathbf{I}_O가 즉시 따른다.
• 양정치성: 임의의 비영벡터 \boldsymbol{\omega}에 대하여 회전 운동 에너지 T = \frac{1}{2} \boldsymbol{\omega}^\top \mathbf{I}_O \boldsymbol{\omega} > 0이므로 \mathbf{I}_O는 양정치 행렬이다. 단, 강체가 한 직선에 분포한 이상적 경우에는 양반정치이다.
• 양의 대각 성분: 모든 대각 성분은 양의 값이다.
• 삼각 부등식: 임의의 두 대각 성분의 합은 나머지 대각 성분보다 항상 크거나 같다.

이러한 성질들로 인하여 관성 텐서는 실대칭 양정치 행렬에 관한 모든 결과를 직접 적용할 수 있는 대상이 된다.

3. 주축 정리

정리 (주축 정리, principal axis theorem). 임의의 강체의 관성 텐서 \mathbf{I}_O에 대하여 다음을 만족하는 직교 좌표계가 존재한다.

1. 그 좌표계에서 관성 텐서는 대각 행렬로 표현된다.
2. 대각 성분 I_1, I_2, I_3는 모두 양수이며, 이들을 주관성 모멘트(principal moments of inertia)라 한다.
3. 그 좌표계의 좌표축들을 주축(principal axes)이라 한다.

이 정리는 실대칭 행렬의 스펙트럼 정리의 직접적인 결과이다. 관성 텐서가 실대칭이므로 직교 행렬 R이 존재하여

R^\top \mathbf{I}_O R = \mathrm{diag}(I_1, I_2, I_3)

이 성립한다. R의 열 벡터 \mathbf{e}_1, \mathbf{e}_2, \mathbf{e}_3가 주축의 방향을 결정한다.

4. 주축의 계산 절차

주어진 관성 텐서 \mathbf{I}_O로부터 주축과 주관성 모멘트를 계산하는 절차는 다음과 같다.

1. 특성 방정식 \det(\mathbf{I}_O - \lambda \mathbb{I}_3) = 0을 구성한다.
2. 특성 방정식의 세 근 \lambda_1, \lambda_2, \lambda_3를 구한다. 이들이 주관성 모멘트이다.
3. 각 고유값 \lambda_i에 대하여 동차 연립 방정식 (\mathbf{I}_O - \lambda_i \mathbb{I}_3)\mathbf{e}_i = \mathbf{0}을 풀어 고유벡터를 구한다.
4. 고유벡터들을 정규화하고, 필요한 경우 직교화를 통하여 정규 직교 기저를 구성한다.
5. 정규 직교 기저를 열로 가지는 직교 행렬 R을 구성한다. 이 행렬은 원래 좌표계에서 주축 좌표계로의 회전 변환에 대응한다.

5. 주축의 기하학적 의미

주축은 강체의 질량 분포가 가지는 대칭성과 밀접한 관계를 가진다. 강체가 어떤 평면에 대하여 거울 대칭을 가지면 그 평면의 법선 방향은 주축이 된다. 또한 강체가 어떤 축에 대하여 회전 대칭을 가지면 그 축은 주축이며, 이에 수직인 임의의 방향도 주축이 된다.

주관성 모멘트의 크기 순서에 따라 강체를 다음과 같이 분류한다.

• I_1 = I_2 = I_3: 구형 대칭(spherically symmetric) 강체. 임의의 축이 주축이다.
• I_1 = I_2 \neq I_3: 대칭 팽이(symmetric top). 비축 두 개는 임의로 선택 가능하다.
• I_1 \neq I_2 \neq I_3: 비대칭 팽이(asymmetric top). 세 주축이 모두 유일하게 결정된다.

6. 평행축 정리와 회전축 정리

주축 분석은 종종 다른 기준점에 대한 관성 텐서로부터 출발하여야 한다. 이때 다음의 두 정리가 활용된다.

평행축 정리 (parallel axis theorem). 강체의 질량 중심 G에 대한 관성 텐서를 \mathbf{I}_G라 하고, 평행한 다른 기준점 O가 질량 중심으로부터 \mathbf{d}만큼 떨어져 있다고 하자. 그러면

\mathbf{I}_O = \mathbf{I}_G + m\left( \|\mathbf{d}\|^2 \mathbb{I}_3 - \mathbf{d}\mathbf{d}^\top \right)

이 성립한다. 여기서 m은 강체의 총 질량이다.

회전축 변환 (rotation axis transformation). 좌표계 A에서의 관성 텐서가 \mathbf{I}^A이고, 좌표계 B가 A에 대하여 회전 행렬 R_B^A로 표현되는 자세를 가진다면, B 좌표계에서의 관성 텐서는 다음과 같이 변환된다.

\mathbf{I}^B = (R_B^A)^\top \mathbf{I}^A R_B^A

이는 텐서의 기본 변환 규칙이며, 좌표계 변환이 행렬의 닮음 변환과 일치한다.

7. 관성 타원체

주축과 주관성 모멘트의 시각화는 관성 타원체(inertia ellipsoid)를 통하여 이루어진다. 다음과 같이 정의된 곡면

\mathbf{x}^\top \mathbf{I}_O \mathbf{x} = 1

은 \mathbb{R}^3에서 타원체이며, 이를 관성 타원체라 한다. 주축 좌표계에서 이 타원체의 방정식은

I_1 x_1^2 + I_2 x_2^2 + I_3 x_3^2 = 1

이 되며, 타원체의 주축의 길이는 1/\sqrt{I_1}, 1/\sqrt{I_2}, 1/\sqrt{I_3}로 결정된다. 임의의 축 \hat{\mathbf{n}}에 대한 관성 모멘트는 그 축이 관성 타원체와 만나는 점까지의 거리의 제곱의 역수와 같다. 이 기하학적 표현은 강체의 회전 관성 분포를 직관적으로 이해하게 한다.

8. 오일러 방정식의 단순화

주축 좌표계에서의 회전 운동 방정식은 가장 단순한 형태인 오일러 방정식(Euler’s equations)으로 표현된다.

\begin{aligned} 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 \end{aligned}

여기서 \boldsymbol{\omega} = (\omega_1, \omega_2, \omega_3)^\top는 주축 좌표계에서 표현된 각속도 벡터이고, \boldsymbol{\tau}는 외부 토크이다. 임의의 좌표계에서는 비대각 관성곱 항으로 인하여 운동 방정식이 매우 복잡해지지만, 주축 좌표계에서는 위와 같이 간결하게 표현된다.

9. 운동 에너지와 각운동량

주축 좌표계에서 강체의 회전 운동 에너지와 각운동량은 다음과 같이 단순화된다.

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

\mathbf{L} = (I_1 \omega_1, I_2 \omega_2, I_3 \omega_3)^\top

이로부터 각속도 벡터와 각운동량 벡터가 일반적으로 평행하지 않으며, 오직 각속도가 주축 방향일 때에만 평행함을 알 수 있다. 이 사실은 비대칭 강체의 자유 회전이 보이는 복잡한 진동 운동의 원인이 된다.

10. 중간 축 정리

비대칭 팽이의 자유 회전에 관한 흥미로운 결과로 중간 축 정리(intermediate axis theorem) 또는 테니스 라켓 정리(tennis racket theorem)가 있다. 주관성 모멘트가 I_1 < I_2 < I_3로 정렬된 경우, 가장 작은 관성 모멘트의 주축(\mathbf{e}_1)과 가장 큰 관성 모멘트의 주축(\mathbf{e}_3) 주위의 회전은 안정인 반면, 중간 관성 모멘트의 주축(\mathbf{e}_2) 주위의 회전은 불안정하다. 이는 오일러 방정식의 평형점 부근 선형 안정성 해석을 통하여 증명될 수 있으며, 우주에서 자유 회전하는 강체의 자세 거동을 설명한다.

11. 로봇공학에서의 응용

11.1 로봇 링크의 동역학 모델링

다관절 로봇의 각 링크에 대한 관성 텐서는 동역학 방정식의 핵심 매개 변수이다. 일반적으로 컴퓨터 보조 설계 데이터로부터 추출된 관성 텐서는 임의의 좌표계에서 표현되어 있으므로, 주축 분석을 통하여 각 링크의 주축과 주관성 모멘트를 추출하면 동역학 계산의 효율성과 수치적 안정성이 향상된다. 또한 주축 정보는 링크의 균형 설계와 진동 특성 평가에 활용된다.

11.2 매니퓰레이터의 효율적 동역학 알고리듬

뉴턴-오일러 재귀 알고리듬과 합성 강체 알고리듬에서 각 링크의 관성 텐서를 주축 좌표계로 표현하면 행렬 곱셈의 부동 소수점 연산 횟수가 크게 감소한다. 특히 대각 형태의 관성 텐서는 영의 곱셈을 회피하게 하며, 실시간 동역학 계산의 효율성을 향상시킨다.

11.3 무인 항공기와 위성의 자세 제어

무인 항공기와 위성의 자세 제어 알고리듬은 일반적으로 주축 좌표계에서 설계된다. 주축 좌표계에서 표현된 오일러 방정식이 가지는 단순한 구조는 제어기 설계와 안정성 해석을 크게 단순화한다. 또한 중간 축 정리는 위성의 자유 회전 자세 안정성 평가와 자세 안정화 전략 수립에 활용된다.

11.4 두 발 로봇의 균형 제어

두 발 보행 로봇의 전체 시스템에 대한 합성 관성 텐서는 자세에 따라 변화하며, 각 자세에서의 주축 분석은 균형 유지를 위한 영점 모멘트 점 계산과 무게 중심 추적 제어에 필요한 정보를 제공한다. 또한 충격 시 발생하는 회전 운동량의 분배 전략 수립에도 활용된다.

11.5 충돌 동역학과 충격 해석

로봇과 환경의 충돌 시 발생하는 충격은 강체의 관성 텐서에 의하여 결정되는 충격 응답 행렬로 모델링된다. 이 행렬의 고유 구조는 충돌 후 회전 거동을 결정하며, 안전한 협동 로봇 설계와 충돌 회피 제어에 활용된다.

11.6 점 구름으로부터의 관성 추정

물체의 외형이 점 구름으로 표현되는 경우, 점 구름의 공분산 행렬은 균일 밀도 가정 아래에서 관성 텐서와 단순한 관계를 가진다. 공분산 행렬의 주성분 분석을 통하여 물체의 주축과 길이 방향, 폭 방향, 두께 방향을 추출할 수 있으며, 이는 물체 인식, 자세 추정, 그리고 로봇 파지 계획에 활용된다.

11.7 모바일 로봇의 동특성 분석

차량형 모바일 로봇의 차체 관성 텐서에 대한 주축 분석은 차량의 횡요(roll), 종요(pitch), 편요(yaw) 운동의 결합 특성을 결정한다. 주축이 차체 좌표축과 일치하지 않는 경우 운동들 사이에 결합이 발생하며, 이는 안정성 마진과 조향 응답 특성에 영향을 미친다.

참고문헌

• Goldstein, H., Poole, C., & Safko, J. (2014). Classical Mechanics (3rd ed.). Pearson.
• Greenwood, D. T. (2003). Advanced Dynamics. Cambridge University Press.
• Featherstone, R. (2008). Rigid Body Dynamics Algorithms. Springer.
• Spong, M. W., Hutchinson, S., & Vidyasagar, M. (2020). Robot Modeling and Control (2nd ed.). Wiley.
• Lynch, K. M., & Park, F. C. (2017). Modern Robotics: Mechanics, Planning, and Control. Cambridge University Press.
• Wie, B. (2008). Space Vehicle Dynamics and Control (2nd ed.). AIAA Education Series.

Version: 1.0