16.47 해밀턴 역학의 로봇 동역학 모델링 응용

1. 개요

해밀턴 역학은 로봇 동역학 모델링에서 라그랑주 정식화와 함께 가장 널리 사용되는 분석 역학적 도구이다. 라그랑주 정식화가 일반화 좌표와 일반화 속도를 기본 변수로 채택하는 반면, 해밀턴 정식화는 일반화 좌표와 일반화 운동량을 정준 변수의 짝으로 채택하여 위상 공간 위에서 운동을 기술한다. 이러한 차이는 단순한 변수 선택의 문제를 넘어, 동역학 분석, 수치 적분, 제어 설계, 그리고 시스템의 보존 법칙 활용 측면에서 본질적인 차이를 가져온다.

본 절에서는 해밀턴 역학이 로봇 시스템의 동역학 모델링에 적용될 때 어떠한 형태로 정식화되는지, 그리고 이 정식화가 매니퓰레이터, 모바일 로봇, 다족 보행 로봇, 부유체 로봇 등 다양한 로봇 클래스의 분석에 어떠한 학술적 가치를 제공하는지를 기술한다.

2. 매니퓰레이터의 해밀턴 정식화

$n$ 자유도를 가진 강체 매니퓰레이터의 라그랑지안은 다음과 같이 표현된다.

$L(q, \dot{q}) = \frac{1}{2} \dot{q}^T M(q) \dot{q} - V(q)$

여기서 $q \in \mathbb{R}^n$ 은 관절 좌표 벡터, $M(q)$ 는 대칭 양의 정부호 관성 행렬(inertia matrix), $V(q)$ 는 중력 위치 에너지이다. 일반화 운동량은 다음과 같이 정의된다.

$p = \frac{\partial L}{\partial \dot{q}} = M(q) \dot{q}$

관성 행렬이 양의 정부호이므로 이 관계는 가역적이며, 일반화 속도를 운동량으로 다음과 같이 표현할 수 있다.

$\dot{q} = M^{-1}(q) p$

해밀터니안은 르장드르 변환을 통하여 다음과 같이 도출된다.

$H(q, p) = \frac{1}{2} p^T M^{-1}(q) p + V(q)$

이 함수는 시스템의 총 기계 에너지(운동 에너지와 위치 에너지의 합)와 일치하며, 보존 시스템에서는 운동 상수가 된다.

3. 정준 운동 방정식

해밀턴의 정준 방정식을 매니퓰레이터에 적용하면 다음을 얻는다.

$\dot{q} = \frac{\partial H}{\partial p} = M^{-1}(q) p$

$\dot{p} = -\frac{\partial H}{\partial q} = -\frac{1}{2} \frac{\partial}{\partial q} \left( p^T M^{-1}(q) p \right) - \frac{\partial V}{\partial q}$

여기서 두 번째 방정식의 첫 번째 항은 코리올리 및 원심력 효과를 운동량 좌표로 표현한 것이며, 라그랑주 정식화에서의 코리올리 행렬 $C(q, \dot{q})$ 와 동등한 정보를 담고 있다. 외부 일반화 힘 $\tau$ 가 작용하는 경우 두 번째 방정식은 다음과 같이 확장된다.

$\dot{p} = -\frac{\partial H}{\partial q} + \tau$

이는 매니퓰레이터의 해밀턴 정식화에서의 입력 동역학을 기술하는 표준 형태이다.

4. 라그랑주 정식화와의 비교

라그랑주 정식화의 표준 형태는 다음과 같다.

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

해밀턴 정식화는 동일한 물리를 일차 미분 방정식의 짝으로 기술한다. 두 정식화의 주요 차이는 다음과 같다. 첫째, 해밀턴 정식화는 $2n$ 차원 위상 공간 위의 일차 시스템이므로 수치 적분과 상태 공간 분석에 직접적으로 적합하다. 둘째, 일반화 운동량 $p$ 는 외부 힘이 가해지지 않을 때 회전 대칭 또는 병진 대칭 좌표에 대하여 자연스럽게 보존되는 양이며, 보존 법칙의 활용이 직관적이다. 셋째, 푸아송 괄호 형식주의를 통하여 운동 상수와 대칭의 분석이 체계화된다.

다만 매니퓰레이터의 직접적인 제어 설계에서는 토크가 가속도와 직접 연결되는 라그랑주 형태가 더 익숙하기 때문에, 실용적 모델링에서는 두 정식화가 보완적으로 사용된다.

5. 비홀로노믹 시스템과의 호환성

차륜형 모바일 로봇과 같은 비홀로노믹(nonholonomic) 시스템은 속도에 대한 제약을 포함하므로 표준 해밀턴 정식화를 직접 적용하기 어렵다. 이 경우 다음 두 가지 접근이 가능하다. 첫째, 라그랑주 승수를 도입하여 제약 조건을 명시적으로 다루는 디랙 제약 시스템 이론을 적용한다. 둘째, 제약 조건과 양립하는 준 좌표(quasi-coordinate)를 도입하여 축약된 위상 공간 위에서 정식화한다.

후자의 접근은 비홀로노믹 환원(nonholonomic reduction) 이론의 핵심 주제이며, 해밀턴 형식주의를 일반 시스템으로 확장하는 데 중요한 역할을 수행한다.

6. 부동 베이스 시스템과 SE(3) 구조

수중 로봇, 우주 로봇, 비행체와 같은 부동 베이스(floating base) 시스템의 구성 공간은 회전과 병진을 모두 포함하는 특수 유클리드 군 $SE(3)$ 의 구조를 지닌다. 이 경우 위상 공간은 $SE(3)$ 의 코탄젠트 다발 $T^*SE(3)$ 이며, 자연스러운 정준 구조가 정의된다.

부동 베이스 시스템의 해밀터니안은 일반적으로 다음과 같은 형태를 띤다.

$H = \frac{1}{2} \xi^T \mathbb{I}^{-1} \xi + V$

여기서 $\xi \in \mathfrak{se}(3)^*$ 는 공간 운동량 또는 물체 운동량, $\mathbb{I}$ 는 일반화 관성 텐서이다. 이 정식화는 리(Lie) 군 위의 해밀턴 역학으로 일반화되며, 키르히호프(Kirchhoff) 방정식과 오일러-푸앵카레(Euler-Poincaré) 방정식의 통일된 도출을 가능하게 한다.

7. 다물체 시스템과 환원

다관절 매니퓰레이터, 보행 로봇, 우주 로봇 매니퓰레이터와 같은 다물체 시스템의 해밀턴 정식화는 시스템의 대칭성을 활용한 환원(reduction)을 통하여 효율적으로 다루어진다. 시스템이 군 $G$ 의 작용 하에서 대칭적이라면, 운동량 사상(momentum map) $J: T^*Q \to \mathfrak{g}^*$ 의 보존을 이용하여 위상 공간을 축약할 수 있다.

마던-와인스타인(Marsden-Weinstein) 환원 정리는 이러한 환원 절차를 일반화된 형태로 제공하며, 환원된 위상 공간 위에 자연스러운 심플렉틱 구조가 유도됨을 보장한다. 우주 로봇의 경우 전체 시스템의 선운동량과 각운동량이 보존되므로, 이 환원을 통하여 베이스의 자세와 매니퓰레이터의 결합 동역학을 효율적으로 분석할 수 있다.

8. 접촉 동역학과 충돌

다족 보행 로봇과 점프 로봇은 발과 지면 사이의 접촉이 동역학의 핵심을 차지한다. 해밀턴 정식화에서 접촉은 단위 법선 방향의 운동량 변화로 기술되며, 충돌 순간에는 해밀터니안이 비탄성 충돌의 경우 감소하고 탄성 충돌의 경우 보존된다. 이러한 분석은 충돌 충격의 운동량 점프(jump) 조건과 접촉 후 상태의 결정에 활용된다.

또한 푸아송 괄호 형식주의를 사용하여 접촉 제약과 양립하는 운동 방정식을 체계적으로 도출할 수 있으며, 이는 측정 미분 방정식(measure differential equation)과 보완성 문제(complementarity problem)의 해밀턴적 해석을 제공한다.

9. 운동량 보존의 활용

해밀턴 정식화의 가장 큰 장점 중 하나는 운동량 보존 법칙의 명시적 활용이다. 매니퓰레이터의 특정 관절이 순환 좌표(cyclic coordinate)인 경우, 즉 해밀터니안이 그 좌표에 의존하지 않는 경우, 대응하는 일반화 운동량은 보존된다. 이러한 보존 법칙은 다음과 같이 활용된다.

첫째, 시스템의 자유도를 효과적으로 축소하여 분석을 단순화한다. 둘째, 수치 적분에서 보존되어야 할 양을 명시적으로 추적할 수 있어 적분 정확도의 검증 기준을 제공한다. 셋째, 우주 로봇과 같이 외부 힘이 거의 존재하지 않는 시스템에서 베이스의 운동을 매니퓰레이터의 운동으로부터 직접 예측할 수 있다.

10. 에너지 일관성과 수동성

해밀턴 정식화는 시스템의 에너지 흐름을 명시적으로 추적할 수 있는 자연스러운 틀을 제공한다. 외부 입력 $\tau$ 와 일반화 속도 $\dot{q}$ 의 곱은 시스템에 공급되는 일률(power)이며, 다음 관계가 성립한다.

$\frac{dH}{dt} = \tau^T \dot{q} = \tau^T M^{-1}(q) p$

이 관계는 시스템의 수동성(passivity) 분석의 기초가 되며, 에너지 기반 제어 설계의 출발점이 된다. 입력-출력 쌍 $(\tau, \dot{q})$ 를 통하여 시스템은 에너지 저장 요소로 모델링되며, 이러한 관점은 후속 절에서 다루어질 에너지 기반 제어와 포트-해밀턴 정식화로 자연스럽게 확장된다.

11. 단순 예시: 2자유도 평면 매니퓰레이터

두 개의 회전 관절을 가진 평면 매니퓰레이터를 고려한다. 관절 좌표를 $q_1, q_2$ , 링크 길이를 $\ell_1, \ell_2$ , 링크 질량을 $m_1, m_2$ 로 표기하면, 관성 행렬은 다음과 같은 형태를 띤다.

$M(q) = \begin{bmatrix} a + 2b\cos q_2 & c + b\cos q_2 \\ c + b\cos q_2 & c \end{bmatrix}$

여기서 $a, b, c$ 는 질량과 길이의 조합으로 결정되는 상수이다. 일반화 운동량은 $p = M(q) \dot{q}$ 로 주어지며, 해밀터니안은 $H = (1/2) p^T M^{-1}(q) p + V(q)$ 로 표현된다. 정준 방정식을 적분하면 위상 공간 $(q_1, q_2, p_1, p_2)$ 위에서 시스템의 진화를 직접 추적할 수 있다.

12. 본 절의 의의

본 절에서 다룬 해밀턴 역학의 로봇 동역학 모델링 응용은 분석 역학의 추상적 형식주의가 실제 로봇 시스템의 모델링에 어떻게 구체화되는지를 보여 준다. 정준 변수의 짝과 위상 공간 형식주의는 수치 적분의 효율성, 보존 법칙의 명시적 활용, 다물체 시스템의 환원, 그리고 에너지 기반 분석의 통일된 틀을 제공한다.

특히 부동 베이스 시스템, 비홀로노믹 시스템, 접촉 동역학, 그리고 대칭성을 가진 다물체 시스템의 분석에서 해밀턴 정식화는 라그랑주 정식화로는 쉽게 접근하기 어려운 통찰을 제공한다. 이러한 모델링 도구는 후속 절에서 다루어질 에너지 기반 제어와 포트-해밀턴 시스템의 이론적 기초를 형성한다.

13. 학습 권장사항

본 절의 내용을 충분히 이해하기 위하여 라그랑주 매니퓰레이터 동역학, 르장드르 변환, 그리고 표준 해밀턴 역학의 정준 방정식에 대한 선행 학습이 권장된다. 본 절의 내용을 심화 학습하기 위해서는 리 군 위의 해밀턴 역학, 운동량 사상과 환원 정리, 비홀로노믹 환원, 그리고 디랙 제약 시스템에 대한 학습이 권장된다.

14. 참고 문헌

Murray, R. M., Li, Z., and Sastry, S. S. (1994). A Mathematical Introduction to Robotic Manipulation. CRC Press.
Spong, M. W., Hutchinson, S., and Vidyasagar, M. (2020). Robot Modeling and Control (2nd ed.). Wiley.
Marsden, J. E., and Ratiu, T. S. (1999). Introduction to Mechanics and Symmetry (2nd ed.). Springer.
Bullo, F., and Lewis, A. D. (2005). Geometric Control of Mechanical Systems. Springer.
Featherstone, R. (2008). Rigid Body Dynamics Algorithms. Springer.
Bloch, A. M. (2003). Nonholonomic Mechanics and Control. Springer.
Holm, D. D., Schmah, T., and Stoica, C. (2009). Geometric Mechanics and Symmetry: From Finite to Infinite Dimensions. Oxford University Press.

version: 1.0