7.68 해밀턴의 원리와 최소 작용 원리

1. 해밀턴의 원리

1.1 정식화

해밀턴의 원리(Hamilton’s principle)는 역학의 변분적 정식화로, 다음과 같이 기술된다.

해밀턴의 원리. 홀로노믹 시스템의 실제 운동 궤적 \mathbf{q}^*(t)는 동일한 경계 조건 \mathbf{q}(t_0) = \mathbf{q}_0, \mathbf{q}(t_f) = \mathbf{q}_f를 만족하는 모든 허용 궤적 중에서 작용 적분의 정상값을 달성하는 궤적이다.

\delta S = \delta\int_{t_0}^{t_f}\mathcal{L}(\mathbf{q}, \dot{\mathbf{q}}, t)\,dt = 0

이 원리는 뉴턴의 운동 법칙과 동등하나, 스칼라 함수(라그랑지안)만으로 시스템의 동역학을 완전히 기술한다는 점에서 벡터 역학과 구별된다.

1.2 보존계에 대한 적용

보존력만이 작용하는 시스템에서 라그랑지안은 \mathcal{L} = T - U이다. 해밀턴의 원리를 적용하면 오일러-라그랑주 방정식이 유도된다.

\frac{d}{dt}\frac{\partial\mathcal{L}}{\partial\dot{q}_i} - \frac{\partial\mathcal{L}}{\partial q_i} = 0, \quad i = 1, 2, \ldots, n

1.3 비보존력의 포함

비보존력(관절 토크, 마찰 등)이 존재하는 경우, 해밀턴의 원리는 확장된 형태를 갖는다.

\delta\int_{t_0}^{t_f}\mathcal{L}\,dt + \int_{t_0}^{t_f}\sum_{i=1}^{n} Q_i^{nc}\delta q_i\,dt = 0

여기서 Q_i^{nc}는 비보존 일반화 힘이다. 이로부터 다음을 얻는다.

\frac{d}{dt}\frac{\partial\mathcal{L}}{\partial\dot{q}_i} - \frac{\partial\mathcal{L}}{\partial q_i} = Q_i^{nc}

2. 최소 작용 원리

2.1 모페르튀의 원리

모페르튀의 최소 작용 원리(Maupertuis’ principle of least action)는 해밀턴의 원리와 관련되나 구별되는 변분 원리이다. 에너지가 보존되는 시스템에서, 실제 궤적은 축약 작용(abbreviated action)

S_0 = \int_A^B \mathbf{p}\cdot d\mathbf{q} = \int_A^B \sum_{i=1}^{n} p_i\,dq_i

을 최소화(또는 정상화)하는 형상 공간의 경로이다. 여기서 p_i = \partial\mathcal{L}/\partial\dot{q}_i는 일반화 운동량이며, 적분은 형상 공간의 경로를 따라 수행된다.

2.2 해밀턴의 원리와의 차이

해밀턴의 원리는 시간이 고정된 경로 사이에서의 변분이고, 최소 작용 원리는 에너지가 고정된 경로 사이에서의 변분이다.

3. 정상 조건과 극소 조건

3.1 정상 작용의 의미

해밀턴의 원리에서 작용이 반드시 최소가 되는 것은 아니다. 작용의 정상값은 극소, 극대, 또는 안장점일 수 있다. 그러나 충분히 짧은 시간 구간에서는 작용이 실제로 극소가 되며, 이 경우 “최소 작용 원리“라는 명칭이 정당화된다.

구체적으로, 야코비 조건이 만족되면(구간 내에 켤레점이 존재하지 않으면) 작용은 극소이다. 대부분의 로봇 제어 시나리오에서는 시간 구간이 충분히 짧아 극소 조건이 만족된다.

3.2 르장드르 조건과의 연결

운동 에너지 T = \frac{1}{2}\dot{\mathbf{q}}^T M(\mathbf{q})\dot{\mathbf{q}}에서 관성 행렬 M(\mathbf{q})가 양정치이므로, 르장드르 조건 \partial^2\mathcal{L}/\partial\dot{q}_i\partial\dot{q}_j > 0이 자동으로 만족된다. 따라서 역학 시스템의 라그랑지안은 강한 르장드르 조건을 만족하며, 이는 극소를 위한 필요 조건의 하나이다.

4. 해밀턴의 정준 방정식

4.1 르장드르 변환

라그랑지안 정식화에서 해밀턴 정식화로의 전환은 르장드르 변환(Legendre transformation)에 의한다. 일반화 운동량을 정의한다.

p_i = \frac{\partial\mathcal{L}}{\partial\dot{q}_i}

해밀토니안은 다음과 같이 정의된다.

H(\mathbf{q}, \mathbf{p}, t) = \sum_{i=1}^{n} p_i\dot{q}_i - \mathcal{L}(\mathbf{q}, \dot{\mathbf{q}}, t)

여기서 우변의 \dot{q}_i는 p_i = \partial\mathcal{L}/\partial\dot{q}_i의 역관계로부터 \mathbf{q}와 \mathbf{p}의 함수로 표현된다.

4.2 정준 방정식

해밀토니안으로부터 다음의 해밀턴 정준 방정식(Hamilton’s canonical equations)이 유도된다.

\dot{q}_i = \frac{\partial H}{\partial p_i}, \quad \dot{p}_i = -\frac{\partial H}{\partial q_i}, \quad i = 1, \ldots, n

이 2n개의 1차 상미분 방정식은 오일러-라그랑주 방정식의 n개의 2차 방정식과 수학적으로 동등하다. 정준 방정식은 위상 공간(phase space) (\mathbf{q}, \mathbf{p})에서의 시스템 거동을 기술하며, 심플렉틱 기하학과 최적 제어 이론의 기초를 형성한다.

4.3 해밀토니안의 물리적 의미

역학 시스템에서 운동 에너지가 일반화 속도의 동차 이차 형식이면, 해밀토니안은 시스템의 총 에너지와 같다.

H = T + U = E

\partial H/\partial t = 0이면 dH/dt = 0이므로 에너지가 보존된다.

5. 로봇 동역학에서의 의의

5.1 에너지 기반 정식화의 장점

해밀턴의 원리에 기반한 라그랑주/해밀턴 정식화는 다음의 장점을 갖는다.

1. 좌표 독립성: 어떤 일반화 좌표를 선택하더라도 동일한 운동 방정식을 산출한다.
2. 구속 조건의 자동 처리: 일반화 좌표로 표현하면 홀로노믹 구속이 자동으로 제거된다.
3. 스칼라 함수 기반: 에너지는 스칼라이므로 좌표 변환이 벡터 역학보다 단순하다.
4. 보존 법칙의 체계적 도출: 대칭성과 보존 법칙의 관계(뇌터 정리)를 명시적으로 활용할 수 있다.

5.2 최적 제어와의 연결

해밀토니안 역학의 구조는 최적 제어 이론에서 폰트랴긴의 최대 원리와 직접 대응한다. 최적 제어 해밀토니안은 역학 해밀토니안의 일반화이며, 공상태 변수(costate variable)는 일반화 운동량에 대응한다. 이 수학적 유사성은 로봇 동역학과 궤적 최적화를 통일된 틀에서 다룰 수 있게 한다.

6. 요약

해밀턴의 원리는 역학 시스템의 실제 운동이 작용 적분의 정상값을 달성한다는 변분 원리이다. 이로부터 오일러-라그랑주 방정식이 유도되며, 르장드르 변환을 통해 해밀턴의 정준 방정식으로 변환된다. 관성 행렬의 양정치성은 르장드르 조건을 자동으로 보장한다. 에너지 기반 정식화는 좌표 독립성, 구속 처리의 용이성, 보존 법칙의 체계적 도출 등의 장점을 가지며, 로봇 동역학의 유도와 최적 제어 이론의 수학적 기반을 동시에 제공한다.

참고 문헌

• Goldstein, H., Poole, C., & Safko, J. (2002). Classical Mechanics (3rd ed.). Addison-Wesley.
• Arnold, V. I. (1989). Mathematical Methods of Classical Mechanics (2nd ed.). Springer.
• Lanczos, C. (1986). The Variational Principles of Mechanics (4th ed.). Dover Publications.
• Siciliano, B., Sciavicco, L., Villani, L., & Oriolo, G. (2009). Robotics: Modelling, Planning and Control. Springer.
• Liberzon, D. (2012). Calculus of Variations and Optimal Control Theory: A Concise Introduction. Princeton University Press.

version: 0.1.0