6.128 관성 텐서의 정의와 계산

1. 개요

관성 텐서(inertia tensor)는 강체의 질량 분포가 회전 운동에 미치는 영향을 기술하는 2차 대칭 텐서이다. 로봇공학에서 각 링크의 관성 텐서는 동역학 모델링의 핵심 매개변수이며, 뉴턴-오일러 방정식과 라그랑주 방정식 모두에서 필수적으로 사용된다. 본 절에서는 관성 텐서의 수학적 정의, 주요 성질, 그리고 다양한 형상에 대한 계산 방법을 다룬다.

2. 관성 텐서의 정의

2.1 질점계의 관성 텐서

원점 $O$ 에 대한 질점계의 관성 텐서는 다음과 같이 정의된다. $N$ 개의 질점 $m_k$ 가 위치 $\mathbf{r}_k = (x_k, y_k, z_k)^T$ 에 있을 때,

$I_{ij} = \sum_{k=1}^{N} m_k \left( \lVert \mathbf{r}_k \rVert^2 \delta_{ij} - r_{k,i} \, r_{k,j} \right)$

여기서 $\delta_{ij}$ 는 크로네커 델타이고, $r_{k,i}$ 는 $k$ 번째 질점 위치 벡터의 $i$ 번째 성분이다.

2.2 연속체의 관성 텐서

연속적인 질량 분포를 가진 강체에 대해, 관성 텐서는 다음의 적분으로 정의된다.

$I_{ij} = \int_V \rho(\mathbf{r}) \left( \lVert \mathbf{r} \rVert^2 \delta_{ij} - r_i \, r_j \right) dV$

여기서 $\rho(\mathbf{r})$ 는 질량 밀도 함수이고, $V$ 는 강체가 차지하는 체적이다.

2.3 행렬 표현

관성 텐서를 $3 \times 3$ 행렬로 나타내면 다음과 같다.

$\mathbf{I} = \begin{bmatrix} I_{xx} & I_{xy} & I_{xz} \\ I_{yx} & I_{yy} & I_{yz} \\ I_{zx} & I_{zy} & I_{zz} \end{bmatrix}$

대각 성분은 **관성 모멘트(moment of inertia)**이다.

$I_{xx} = \int_V \rho(y^2 + z^2) \, dV$

$I_{yy} = \int_V \rho(x^2 + z^2) \, dV$

$I_{zz} = \int_V \rho(x^2 + y^2) \, dV$

비대각 성분은 **관성 상승적(product of inertia)**이며, 음의 부호를 포함한다.

$I_{xy} = I_{yx} = -\int_V \rho \, x \, y \, dV$

$I_{xz} = I_{zx} = -\int_V \rho \, x \, z \, dV$

$I_{yz} = I_{zy} = -\int_V \rho \, y \, z \, dV$

3. 관성 텐서의 성질

3.1 대칭성

관성 텐서는 대칭이다.

$I_{ij} = I_{ji}$

따라서 $3 \times 3$ 관성 텐서는 6개의 독립 성분을 가진다.

3.2 양의 정부호성

관성 텐서는 양의 준정부호(positive semi-definite) 행렬이다. 질량이 0이 아닌 물체에 대해, 3개 이상의 질점이 비공선(non-collinear)이면 양의 정부호(positive definite)이다.

$\boldsymbol{\omega}^T \mathbf{I} \boldsymbol{\omega} \geq 0, \quad \forall \boldsymbol{\omega} \in \mathbb{R}^3$

이 이차 형식은 회전축 $\boldsymbol{\omega}$ 에 대한 관성 모멘트를 나타낸다.

$I_{\boldsymbol{\omega}} = \hat{\boldsymbol{\omega}}^T \mathbf{I} \hat{\boldsymbol{\omega}}$

여기서 $\hat{\boldsymbol{\omega}} = \boldsymbol{\omega} / \lVert \boldsymbol{\omega} \rVert$ 는 회전축 방향의 단위 벡터이다.

3.3 주관성 모멘트와 주축

관성 텐서는 대칭 행렬이므로 직교 대각화가 가능하다.

$\mathbf{I} = \mathbf{R} \, \text{diag}(I_1, I_2, I_3) \, \mathbf{R}^T$

여기서 $I_1, I_2, I_3$ 은 **주관성 모멘트(principal moments of inertia)**이고, $\mathbf{R}$ 의 열벡터는 **주축(principal axes)**이다. 주관성 모멘트는 관성 텐서의 고유값이며, 주축은 대응하는 고유벡터이다.

주관성 모멘트 사이에는 삼각 부등식이 성립한다.

$I_1 \leq I_2 + I_3, \quad I_2 \leq I_1 + I_3, \quad I_3 \leq I_1 + I_2$

4. 평행축 정리

4.1 슈타이너 정리

질량 중심 $C$ 에 대한 관성 텐서 $\mathbf{I}_C$ 를 알고 있을 때, 임의의 점 $O$ 에 대한 관성 텐서 $\mathbf{I}_O$ 는 다음과 같이 계산된다.

$I_{O,ij} = I_{C,ij} + m\left( \lVert \mathbf{d} \rVert^2 \delta_{ij} - d_i \, d_j \right)$

여기서 $\mathbf{d}$ 는 질량 중심 $C$ 에서 점 $O$ 까지의 벡터이고, $m$ 은 물체의 총 질량이다.

행렬 형태로 표현하면 다음과 같다.

$\mathbf{I}_O = \mathbf{I}_C + m\left( \mathbf{d}^T\mathbf{d} \, \mathbf{E}_3 - \mathbf{d}\mathbf{d}^T \right)$

여기서 $\mathbf{E}_3$ 은 $3 \times 3$ 단위 행렬이다. 반대칭 행렬을 이용하면 다음과 같이 쓸 수도 있다.

$\mathbf{I}_O = \mathbf{I}_C - m \, [\mathbf{d}]^2$

여기서 $[\mathbf{d}]$ 는 $\mathbf{d}$ 의 반대칭 행렬이다.

5. 좌표 변환

좌표계를 회전 행렬 $R$ 로 변환할 때, 관성 텐서는 다음과 같이 변환된다.

$\mathbf{I}' = R \, \mathbf{I} \, R^T$

지표 표기법으로는 다음과 같다.

$I'_{ij} = R_{ik} \, R_{jl} \, I_{kl}$

이는 $(0,2)$ -형 텐서의 변환 규칙과 일치한다.

6. 기본 형상의 관성 텐서 계산

6.1 균질 직육면체

변의 길이가 $a$ , $b$ , $c$ 이고 질량이 $m$ 인 균질 직육면체의 질량 중심에 대한 관성 텐서:

$\mathbf{I}_C = \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}$

6.2 균질 원기둥

반지름 $r$ , 높이 $h$ , 질량 $m$ 인 균질 원기둥의 질량 중심에 대한 관성 텐서(축이 $z$ 축 방향):

$\mathbf{I}_C = \begin{bmatrix} \frac{m}{12}(3r^2 + h^2) & 0 & 0 \\ 0 & \frac{m}{12}(3r^2 + h^2) & 0 \\ 0 & 0 & \frac{m}{2}r^2 \end{bmatrix}$

6.3 균질 구

반지름 $r$ , 질량 $m$ 인 균질 구의 질량 중심에 대한 관성 텐서:

$\mathbf{I}_C = \frac{2}{5}mr^2 \begin{bmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{bmatrix} = \frac{2}{5}mr^2 \, \mathbf{E}_3$

구는 등방성(isotropy)을 가지므로 관성 텐서가 단위 행렬의 스칼라 배이다.

6.4 균질 타원체

반축 길이가 $a$ , $b$ , $c$ 인 균질 타원체의 질량 중심에 대한 관성 텐서:

$\mathbf{I}_C = \frac{m}{5} \begin{bmatrix} b^2 + c^2 & 0 & 0 \\ 0 & a^2 + c^2 & 0 \\ 0 & 0 & a^2 + b^2 \end{bmatrix}$

7. 공간 관성 행렬

7.1 차원 공간 관성 행렬

로봇 링크의 공간 관성 행렬(spatial inertia matrix)은 병진 및 회전 운동을 통합하여 $6 \times 6$ 행렬로 표현된다.

$\mathcal{G} = \begin{bmatrix} \mathbf{I} & m[\mathbf{p}_c] \\ m[\mathbf{p}_c]^T & m\mathbf{E}_3 \end{bmatrix}$

여기서 $\mathbf{I}$ 는 원점에 대한 $3 \times 3$ 관성 텐서, $m$ 은 질량, $\mathbf{p}_c$ 는 원점에서 질량 중심까지의 벡터, $[\mathbf{p}_c]$ 는 그 반대칭 행렬이다.

7.2 운동 에너지와의 관계

공간 관성 행렬을 이용하면 강체의 운동 에너지를 다음과 같이 표현할 수 있다.

$K = \frac{1}{2} \mathcal{V}^T \mathcal{G} \, \mathcal{V}$

여기서 $\mathcal{V} = (\boldsymbol{\omega}, \mathbf{v})^T$ 는 물체의 트위스트이다.

8. 합성체의 관성 텐서

여러 부분체로 구성된 물체의 관성 텐서는 각 부분체의 관성 텐서를 더하여 구한다. 각 부분체의 관성 텐서를 공통 좌표계와 공통 기준점으로 변환한 후 합산한다.

$\mathbf{I}_{\text{total}} = \sum_{k=1}^{N} \mathbf{I}_k$

여기서 $\mathbf{I}_k$ 는 $k$ 번째 부분체의 공통 기준점에 대한 관성 텐서이다. 각 부분체의 질량 중심이 다를 경우 평행축 정리를 먼저 적용해야 한다.

9. 수치 계산 예제

9.1 예제: 로봇 링크의 관성 텐서

길이 $L = 0.5\,\text{m}$ , 반지름 $r = 0.02\,\text{m}$ , 질량 $m = 2\,\text{kg}$ 인 원기둥 형태의 로봇 링크를 고려한다. 질량 중심에 대한 관성 텐서는 다음과 같다.

$I_{xx} = I_{yy} = \frac{m}{12}(3r^2 + L^2) = \frac{2}{12}(3 \times 0.0004 + 0.25) = 0.0419 \, \text{kg} \cdot \text{m}^2$

$I_{zz} = \frac{m}{2}r^2 = \frac{2}{2} \times 0.0004 = 0.0004 \, \text{kg} \cdot \text{m}^2$

링크의 한쪽 끝(관절 위치)에 대한 관성 텐서는 평행축 정리를 적용하여 구한다. $\mathbf{d} = (0, 0, L/2)^T$ 이면,

$I'_{xx} = I'_{yy} = I_{xx} + m\left(\frac{L}{2}\right)^2 = 0.0419 + 2 \times 0.0625 = 0.1669 \, \text{kg} \cdot \text{m}^2$

$I'_{zz} = I_{zz} = 0.0004 \, \text{kg} \cdot \text{m}^2$

10. 참고 문헌

Goldstein, H., Poole, C., & Safko, J. (2002). Classical Mechanics. 3rd ed., Addison-Wesley.
Craig, J. J. (2005). Introduction to Robotics: Mechanics and Control. 3rd ed., Pearson Prentice Hall.
Featherstone, R. (2008). Rigid Body Dynamics Algorithms. Springer.
Lynch, K. M., & Park, F. C. (2017). Modern Robotics: Mechanics, Planning, and Control. Cambridge University Press.
Siciliano, B., Sciavicco, L., Villani, L., & Oriolo, G. (2009). Robotics: Modelling, Planning and Control. Springer.

v 0.1