15.19 코리올리 행렬의 유도
1. 개요
매니퓰레이터의 동역학 방정식에서 코리올리 행렬 \mathbf{C}(\mathbf{q}, \dot{\mathbf{q}})는 코리올리력과 원심력을 함께 표현하는 핵심 항이다. 본 절에서는 라그랑주 방정식으로부터 코리올리 행렬의 유도와 그 형식을 자세히 다룬다.
2. 라그랑주 방정식의 분해
2.1 출발점
라그랑주 방정식은 다음과 같다.
\frac{d}{dt}\frac{\partial L}{\partial \dot{q}_i} - \frac{\partial L}{\partial q_i} = \tau_i
여기서 L = T - U이고 \tau_i는 일반화 힘이다.
2.2 운동 에너지 항의 처리
운동 에너지가 T = \frac{1}{2}\dot{\mathbf{q}}^T\mathbf{M}(\mathbf{q})\dot{\mathbf{q}}인 경우 다음의 미분이 필요하다.
\frac{\partial T}{\partial \dot{q}_i} = \sum_{j}M_{ij}(\mathbf{q})\dot{q}_j
이 양의 시간 미분은
\frac{d}{dt}\frac{\partial T}{\partial \dot{q}_i} = \sum_{j}M_{ij}(\mathbf{q})\ddot{q}_j + \sum_{j,k}\frac{\partial M_{ij}}{\partial q_k}\dot{q}_k\dot{q}_j
2.3 좌표에 대한 편미분
\frac{\partial T}{\partial q_i} = \frac{1}{2}\sum_{j,k}\frac{\partial M_{jk}}{\partial q_i}\dot{q}_j\dot{q}_k
3. 라그랑주 방정식의 정리
3.1 결합
위의 결과를 라그랑주 방정식에 대입하면
\sum_{j}M_{ij}\ddot{q}_j + \sum_{j,k}\frac{\partial M_{ij}}{\partial q_k}\dot{q}_k\dot{q}_j - \frac{1}{2}\sum_{j,k}\frac{\partial M_{jk}}{\partial q_i}\dot{q}_j\dot{q}_k + \frac{\partial U}{\partial q_i} = \tau_i
3.2 항의 분리
- 첫 번째 항: 관성 토크
- 두 번째와 세 번째 항: 코리올리력 및 원심력 토크
- 네 번째 항: 중력 토크
4. 코리올리 행렬의 정의
4.1 크리스토펠 기호
코리올리 행렬의 성분은 크리스토펠 기호의 1종(Christoffel symbols of the first kind)으로 정의된다.
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)
4.2 코리올리 행렬
코리올리 행렬은 크리스토펠 기호로부터 다음과 같이 정의된다.
C_{ij}(\mathbf{q}, \dot{\mathbf{q}}) = \sum_{k=1}^{n}c_{ijk}(\mathbf{q})\dot{q}_k
4.3 동역학 방정식의 형식
이 정의로 라그랑주 방정식은 다음과 같이 단순화된다.
\sum_{j}M_{ij}(\mathbf{q})\ddot{q}_j + \sum_{j}C_{ij}(\mathbf{q}, \dot{\mathbf{q}})\dot{q}_j + \frac{\partial U}{\partial q_i} = \tau_i
또는 행렬 형식으로
\mathbf{M}(\mathbf{q})\ddot{\mathbf{q}} + \mathbf{C}(\mathbf{q}, \dot{\mathbf{q}})\dot{\mathbf{q}} + \mathbf{g}(\mathbf{q}) = \boldsymbol{\tau}
여기서 \mathbf{g}(\mathbf{q}) = \nabla_{\mathbf{q}}U는 중력 항이다.
5. 크리스토펠 기호의 성질
5.1 대칭성
크리스토펠 기호는 두 번째와 세 번째 인덱스에 대해 대칭이다.
c_{ijk} = c_{ikj}
이는 정의의 대칭성에서 직접 확인된다.
5.2 결합 양의 분해
코리올리력 항 \sum_{j}C_{ij}\dot{q}_j를 풀어 쓰면
\sum_{j}C_{ij}\dot{q}_j = \sum_{j,k}c_{ijk}\dot{q}_j\dot{q}_k
여기서 j \neq k인 항이 코리올리력 항이고 j = k인 항이 원심력 항이다.
6. 코리올리 행렬의 다양한 형식
6.1 표준 형식
위에서 정의한 형식은 표준적이며, 안정성 분석에 적합한 반대칭성을 가진다.
6.2 다른 형식
코리올리 행렬을 다른 형식으로 정의할 수 있다. 다음의 두 형식이 자주 사용된다.
C_{ij}^{(1)} = \sum_{k}\frac{\partial M_{ij}}{\partial q_k}\dot{q}_k - \frac{1}{2}\sum_{k}\frac{\partial M_{jk}}{\partial q_i}\dot{q}_k
이는 표준 형식과 동등하다.
6.3 비유일성
코리올리 행렬은 유일하지 않다. 같은 \mathbf{C}\dot{\mathbf{q}} 항을 산출하는 여러 형식의 행렬이 가능하다.
7. 반대칭 성질
7.1 핵심 항등식
표준 형식의 코리올리 행렬은 다음의 반대칭 성질을 만족한다.
\dot{\mathbf{M}}(\mathbf{q}) - 2\mathbf{C}(\mathbf{q}, \dot{\mathbf{q}}) = -(\dot{\mathbf{M}}(\mathbf{q}) - 2\mathbf{C}(\mathbf{q}, \dot{\mathbf{q}}))^T
7.2 증명
\dot{\mathbf{M}}의 성분은 다음과 같다.
\dot{M}_{ij} = \sum_{k}\frac{\partial M_{ij}}{\partial q_k}\dot{q}_k
2C_{ij}를 빼면
\dot{M}_{ij} - 2C_{ij} = \sum_{k}\left(\frac{\partial M_{ij}}{\partial q_k} - \frac{\partial M_{ik}}{\partial q_j} - \frac{\partial M_{ij}}{\partial q_k} + \frac{\partial M_{jk}}{\partial q_i}\right)\dot{q}_k
= \sum_{k}\left(-\frac{\partial M_{ik}}{\partial q_j} + \frac{\partial M_{jk}}{\partial q_i}\right)\dot{q}_k
이 표현은 i, j 교환에 대해 부호가 변한다. 따라서 행렬 \dot{\mathbf{M}} - 2\mathbf{C}는 반대칭이다.
7.3 응용
이 반대칭 성질은 매니퓰레이터의 안정성 증명에서 핵심적으로 사용된다. 리야푸노프 함수의 시간 미분에서 코리올리 항이 사라지게 한다.
8. 응용 예시: 2자유도 평면 매니퓰레이터
8.1 관성 행렬
이전 절에서 유도한 관성 행렬을 다시 활용한다.
M_{12} = m_2(l_{c2}^2 + l_1 l_{c2}\cos q_2) + I_2
8.2 편미분
\frac{\partial M_{12}}{\partial q_2} = -m_2 l_1 l_{c2}\sin q_2
8.3 코리올리 행렬의 성분
C_{11} = c_{112}\dot{q}_2 = \frac{1}{2}\left(\frac{\partial M_{11}}{\partial q_2}\right)\dot{q}_2 = -m_2 l_1 l_{c2}\sin q_2 \dot{q}_2
다른 성분도 유사하게 계산된다.
9. 본 절의 의의
본 절은 라그랑주 방정식으로부터 코리올리 행렬의 유도를 자세히 다루었다. 이는 매니퓰레이터의 동역학 방정식의 표준 형식을 도출하고, 안정성 분석에 활용되는 중요한 결과이다.
10. 참고 문헌
- Spong, M. W., Hutchinson, S., & Vidyasagar, M. (2020). Robot Modeling and Control (2nd ed.). Wiley.
- Murray, R. M., Li, Z., & Sastry, S. S. (1994). A Mathematical Introduction to Robotic Manipulation. CRC Press.
- Slotine, J.-J. E., & Li, W. (1991). Applied Nonlinear Control. Prentice Hall.
- Siciliano, B., Sciavicco, L., Villani, L., & Oriolo, G. (2010). Robotics: Modelling, Planning and Control. Springer.
version: 1.0