14.23 코리올리력 항의 물리적 의미와 계산

1. 개요

매니퓰레이터의 동역학 방정식에서 코리올리력 항은 회전하는 좌표계에서 발생하는 가상의 힘에 해당하며, 두 관절의 속도의 곱에 비례하는 항이다. 이 항은 매니퓰레이터의 빠른 운동에서 중요한 역할을 하며, 정확한 동역학 모형에서 반드시 고려되어야 한다. 본 절에서는 코리올리력 항의 물리적 의미와 계산 방법을 다룬다.

2. 코리올리력의 물리적 배경

2.1 회전 좌표계에서의 가상력

회전 좌표계에서 운동하는 입자는 관성 좌표계의 관측자가 보면 단순한 운동이지만, 회전 좌표계의 관측자가 보면 가상의 힘이 작용하는 것처럼 보인다. 이러한 가상력 중 하나가 코리올리력이다.

2.2 코리올리력의 표현

회전 각속도 $\boldsymbol{\omega}$ 로 회전하는 좌표계에서 속도 $\mathbf{v}_{\text{rel}}$ 로 운동하는 질량 $m$ 의 입자에 작용하는 코리올리력은 다음과 같다.

$\mathbf{F}_{\text{Coriolis}} = -2m\boldsymbol{\omega} \times \mathbf{v}_{\text{rel}}$

이는 회전 운동과 병진 운동의 결합에서 발생한다.

3. 매니퓰레이터 동역학에서의 코리올리력 항

3.1 동역학 방정식의 구조

매니퓰레이터의 동역학 방정식은 다음의 형식이다.

$\mathbf{M}(\mathbf{q})\ddot{\mathbf{q}} + \mathbf{C}(\mathbf{q}, \dot{\mathbf{q}})\dot{\mathbf{q}} + \mathbf{g}(\mathbf{q}) = \boldsymbol{\tau}$

여기서 $\mathbf{C}(\mathbf{q}, \dot{\mathbf{q}})\dot{\mathbf{q}}$ 가 코리올리력과 원심력을 합한 항이다.

3.2 코리올리력 항의 분리

$\mathbf{C}(\mathbf{q}, \dot{\mathbf{q}})\dot{\mathbf{q}}$ 항을 풀어서 쓰면 다음과 같다.

$[\mathbf{C}(\mathbf{q}, \dot{\mathbf{q}})\dot{\mathbf{q}}]_i = \sum_{j=1}^{n}\sum_{k=1}^{n} c_{ijk}(\mathbf{q})\dot{q}_j\dot{q}_k$

이 중 $j \neq k$ 의 교차항이 코리올리력 항에 해당하고, $j = k$ 의 대각항이 원심력 항에 해당한다.

4. 크리스토펠 기호와 계수

4.1 정의

$c_{ijk}$ 는 다음과 같이 정의된다.

$c_{ijk}(\mathbf{q}) = \frac{1}{2}\left(\frac{\partial M_{ij}}{\partial q_k} + \frac{\partial M_{ik}}{\partial q_j} - \frac{\partial M_{jk}}{\partial q_i}\right)$

이는 미분 기하학의 크리스토펠 기호의 1종(Christoffel symbols of the first kind)와 같은 형태이다.

4.2 대칭성

$c_{ijk}$ 는 $j$ 와 $k$ 에 대해 대칭이다.

$c_{ijk} = c_{ikj}$

4.3 코리올리력 항의 명시적 표현

코리올리력 항은 다음과 같이 표현된다.

$[\mathbf{C}_{\text{Coriolis}}]_i = 2\sum_{j<k} c_{ijk}(\mathbf{q})\dot{q}_j\dot{q}_k$

여기서 인수 2는 $j$ 와 $k$ 의 순서를 고려한 것이다.

5. 코리올리력 항의 계산

5.1 관성 행렬의 미분

코리올리력 항의 계산은 관성 행렬 $\mathbf{M}(\mathbf{q})$ 의 관절 위치에 대한 편미분을 필요로 한다. 이는 해석적 또는 수치적으로 수행될 수 있다.

5.2 뉴턴-오일러 알고리즘에 의한 간접 계산

뉴턴-오일러 재귀 알고리즘에서 $\ddot{\mathbf{q}} = 0$ , 중력 무시( $\mathbf{g}_0 = 0$ )로 설정하면 출력되는 토크가 코리올리력과 원심력을 합한 항이다.

$\mathbf{C}(\mathbf{q}, \dot{\mathbf{q}})\dot{\mathbf{q}} = \text{NE}(\mathbf{q}, \dot{\mathbf{q}}, 0, 0)$

이 방법은 효율적이며 매니퓰레이터의 동역학 모의에서 일반적으로 사용된다.

6. 코리올리력 항의 성질

6.1 속도 의존성

코리올리력 항은 관절 속도의 곱에 비례하므로 속도가 증가하면 빠르게 커진다. 빠른 운동에서 그 효과가 두드러진다.

6.2 반대칭성

$\mathbf{C}$ 행렬은 다음과 같은 성질을 만족하도록 선택할 수 있다.

$\dot{\mathbf{M}}(\mathbf{q}) - 2\mathbf{C}(\mathbf{q}, \dot{\mathbf{q}})$

이 행렬은 반대칭이다. 이 성질은 매니퓰레이터의 안정성 증명에서 자주 사용된다.

6.3 에너지의 기여

코리올리력 항과 원심력 항은 일을 하지 않으며, 매니퓰레이터의 운동 에너지 보존에 영향을 주지 않는 형태로 동역학 방정식에 들어간다.

7. 응용

7.1 정밀 추종 제어

빠른 운동의 정밀 추종에서는 코리올리력 항을 보상하는 피드포워드 항이 필요하다. 이는 계산 토크 제어의 핵심이다.

7.2 동적 시뮬레이션

매니퓰레이터의 정확한 동적 시뮬레이션을 위해 코리올리력 항이 정확히 계산되어야 한다.

7.3 안정성 분석

매니퓰레이터의 제어 안정성 증명에서 코리올리력 항의 반대칭성이 핵심적으로 사용된다.

8. 응용 예시: 2자유도 평면 매니퓰레이터

2자유도 평면 매니퓰레이터에서 코리올리력 항은 다음과 같이 나타난다.

$[\mathbf{C}\dot{\mathbf{q}}]_1 = -m_2 l_1 l_{c2}\sin q_2(2\dot{q}_1\dot{q}_2 + \dot{q}_2^2)$

$[\mathbf{C}\dot{\mathbf{q}}]_2 = m_2 l_1 l_{c2}\sin q_2 \dot{q}_1^2$

여기서 $\dot{q}_1\dot{q}_2$ 항이 코리올리력 항이고, $\dot{q}_1^2, \dot{q}_2^2$ 항이 원심력 항이다.

9. 본 절의 의의

본 절은 코리올리력 항의 물리적 의미와 계산 방법을 다루었다. 코리올리력 항은 매니퓰레이터의 빠른 운동의 정확한 모형화에 필수적이며, 정밀 제어와 안정성 분석에서 핵심적이다.

10. 참고 문헌

Spong, M. W., Hutchinson, S., & Vidyasagar, M. (2020). Robot Modeling and Control (2nd ed.). Wiley.
Siciliano, B., Sciavicco, L., Villani, L., & Oriolo, G. (2010). Robotics: Modelling, Planning and Control. Springer.
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