16.49 포트-해밀턴 시스템의 개요

1. 개요

포트-해밀턴 시스템(Port-Hamiltonian system, PHS)은 해밀턴 역학을 에너지 흐름과 외부 환경과의 상호작용 관점에서 일반화한 모델링 패러다임이다. 이 정식화는 1990년대 초반 마쉬케(Bernhard Maschke), 판 데르 샤프트(Arjan van der Schaft) 등에 의하여 체계화되었으며, 이후 수동성 기반 제어, 다물리(multiphysics) 시스템 모델링, 분산 매개변수 시스템 분석 등으로 광범위하게 확장되어 왔다.

전통적인 해밀턴 정식화가 보존 시스템을 위상 공간 위의 정준 방정식으로 기술하는 데 초점을 두는 반면, 포트-해밀턴 시스템은 외부와 에너지를 교환하고 내부 산일(dissipation)을 포함하는 개방형(open) 시스템을 통일적으로 다룬다. 이러한 특성은 포트-해밀턴 시스템을 로봇공학, 메카트로닉스, 전기-기계 결합 시스템, 그리고 인간-로봇 상호작용 시스템의 모델링과 제어에 매우 적합한 도구로 만든다.

본 절에서는 포트-해밀턴 시스템의 기본 정식화, 디랙(Dirac) 구조의 역할, 입력-출력 포트의 의미, 수동성 특성, 그리고 로봇공학적 응용을 학술적으로 기술한다.

2. 명시적 포트-해밀턴 시스템

가장 단순한 형태의 명시적 포트-해밀턴 시스템은 다음과 같이 표현된다.

$\dot{x} = \left[ J(x) - R(x) \right] \frac{\partial H}{\partial x}(x) + g(x) u$

$y = g(x)^T \frac{\partial H}{\partial x}(x)$

여기서 각 기호의 의미는 다음과 같다. $x \in \mathbb{R}^n$ 은 상태 벡터, $H(x): \mathbb{R}^n \to \mathbb{R}$ 은 해밀터니안 또는 에너지 함수, $J(x)$ 는 반대칭 인터커넥션 행렬로서 $J^T = -J$ , $R(x)$ 는 양의 준정부호 댐핑 행렬로서 $R = R^T \geq 0$ , $g(x)$ 는 입력 행렬, $u$ 는 입력 벡터, $y$ 는 출력 벡터이다.

이 정식화의 핵심적 특성은 입력과 출력의 짝 $(u, y)$ 가 외부와의 에너지 교환을 직접 표현한다는 점이다. 곱 $u^T y$ 는 시스템에 공급되는 일률(power)을 의미하며, 이러한 짝을 포트(port)라고 부른다. 포트의 개념은 시스템과 환경 사이의 인터페이스를 물리적 의미가 명확한 형태로 정식화한다.

3. 에너지 균형과 수동성

포트-해밀턴 시스템의 가장 중요한 성질 중 하나는 에너지 균형 방정식이다. 해밀터니안의 시간 미분을 계산하면 다음을 얻는다.

$\frac{dH}{dt} = \frac{\partial H}{\partial x}^T \dot{x} = \frac{\partial H}{\partial x}^T \left[ (J - R) \frac{\partial H}{\partial x} + g u \right]$

$J$ 가 반대칭이므로 $\frac{\partial H}{\partial x}^T J \frac{\partial H}{\partial x} = 0$ 이며, 따라서

$\frac{dH}{dt} = -\frac{\partial H}{\partial x}^T R \frac{\partial H}{\partial x} + \frac{\partial H}{\partial x}^T g u = -\frac{\partial H}{\partial x}^T R \frac{\partial H}{\partial x} + y^T u$

이 등식은 다음과 같이 해석된다. 시스템의 에너지 변화율은 외부로부터 공급되는 일률 $y^T u$ 에서 내부에서 산일되는 일률 $\frac{\partial H}{\partial x}^T R \frac{\partial H}{\partial x} \geq 0$ 을 뺀 것과 같다. 이는 열역학 제1법칙의 시스템 이론적 표현이며, 포트-해밀턴 시스템이 본질적으로 수동(passive)임을 보장한다.

해밀터니안이 유계 아래(bounded below)인 경우, 적절한 상수를 더하여 $H \geq 0$ 로 가정할 수 있으며, 이때 $H$ 는 시스템의 저장 함수(storage function)로서 작용한다. 수동성 부등식 $\dot{H} \leq u^T y$ 가 성립하므로 시스템은 입출력 쌍 $(u, y)$ 에 대하여 수동이다.

4. 디랙 구조와 일반 포트-해밀턴 시스템

명시적 포트-해밀턴 시스템은 디랙(Dirac) 구조라는 더 추상적이고 일반적인 개념의 특수한 경우이다. 디랙 구조는 노력(effort)과 흐름(flow)이라는 두 종류의 변수 사이의 기하학적 관계를 정의하며, 에너지 보존을 본질적인 성질로 내포한다.

흐름 공간 $\mathcal{F}$ 와 노력 공간 $\mathcal{E} = \mathcal{F}^*$ 가 주어질 때, 두 공간 위의 자연스러운 쌍대 짝 $\langle e | f \rangle$ 이 정의되며, 이는 일률을 표현한다. 디랙 구조 $\mathcal{D} \subset \mathcal{F} \times \mathcal{E}$ 는 다음 조건을 만족하는 부분 공간으로 정의된다.

$\mathcal{D} = \mathcal{D}^{\perp}, \quad \forall (f, e) \in \mathcal{D}: \langle e | f \rangle = 0$

이 조건은 디랙 구조 위에서 전체 일률의 합이 영임을 의미하며, 이는 에너지 보존의 추상적 표현이다.

일반 포트-해밀턴 시스템은 상태 다양체 $\mathcal{X}$ , 에너지 함수 $H: \mathcal{X} \to \mathbb{R}$ , 그리고 디랙 구조 $\mathcal{D}$ 로 구성된다. 시스템의 동역학은 다음 관계로 결정된다.

$\left( -\dot{x}, \frac{\partial H}{\partial x}, u, y \right) \in \mathcal{D}$

이 추상적 정식화는 명시적 형태가 다루기 어려운 제약 시스템, 비홀로노믹 시스템, 그리고 분산 매개변수 시스템에 적용 가능하다.

5. 인터커넥션 구조

포트-해밀턴 시스템의 가장 큰 강점은 인터커넥션(interconnection)에 대한 닫힘성(closedness)이다. 두 개의 포트-해밀턴 시스템이 그들의 포트를 통하여 결합되면, 결과로 나오는 시스템 또한 포트-해밀턴 시스템이다. 이러한 성질은 복잡한 다물리 시스템을 모듈화된 부품으로 나누어 각각을 포트-해밀턴 시스템으로 모델링한 후, 인터커넥션을 통하여 전체 시스템을 구성하는 체계적 모델링 방법을 가능하게 한다.

전형적인 인터커넥션은 다음과 같은 형태를 띤다. 두 시스템의 출력이 서로의 입력으로 연결되며, 이때 일률의 합은 보존된다. 이러한 인터커넥션은 디랙 구조의 합성으로 표현되며, 결과 시스템의 디랙 구조도 다시 디랙 구조의 성질을 만족한다.

6. 다물리 시스템과의 자연스러운 부합

포트-해밀턴 시스템은 서로 다른 물리 영역(기계, 전기, 유체, 열)에 속한 시스템을 통일된 언어로 기술할 수 있다. 각 영역에서 노력과 흐름의 짝은 다음과 같이 자연스럽게 식별된다.

기계 병진 시스템: 노력은 힘, 흐름은 속도
기계 회전 시스템: 노력은 토크, 흐름은 각속도
전기 시스템: 노력은 전압, 흐름은 전류
유체 시스템: 노력은 압력, 흐름은 부피 유량
열 시스템: 노력은 온도, 흐름은 엔트로피 유량

이러한 통일성은 본드 그래프(bond graph) 모델링과 자연스럽게 부합하며, 메카트로닉스 시스템과 액추에이터-센서 결합 시스템의 모델링에 매우 효과적이다.

7. 단순 예시: 질량-스프링-댐퍼 시스템

질량 $m$ , 스프링 상수 $k$ , 댐핑 계수 $b$ 를 가진 일차원 질량-스프링-댐퍼 시스템을 고려한다. 상태 변수를 $x_1 = q$ (변위), $x_2 = p = m\dot{q}$ (운동량)으로 선택하면, 해밀터니안은 다음과 같다.

$H(x_1, x_2) = \frac{1}{2 m} x_2^2 + \frac{1}{2} k x_1^2$

포트-해밀턴 형태는 다음과 같이 표현된다.

$\begin{bmatrix} \dot{x}_1 \\ \dot{x}_2 \end{bmatrix} = \left( \begin{bmatrix} 0 & 1 \\ -1 & 0 \end{bmatrix} - \begin{bmatrix} 0 & 0 \\ 0 & b \end{bmatrix} \right) \begin{bmatrix} \partial H / \partial x_1 \\ \partial H / \partial x_2 \end{bmatrix} + \begin{bmatrix} 0 \\ 1 \end{bmatrix} u$

$y = \begin{bmatrix} 0 & 1 \end{bmatrix} \begin{bmatrix} \partial H / \partial x_1 \\ \partial H / \partial x_2 \end{bmatrix} = \frac{x_2}{m} = \dot{q}$

여기서 입력 $u$ 는 외부 힘, 출력 $y$ 는 속도이며, 곱 $u y$ 는 외부에서 공급되는 일률에 정확히 대응한다.

8. 매니퓰레이터의 포트-해밀턴 정식화

$n$ 자유도 매니퓰레이터의 포트-해밀턴 형태는 다음과 같이 표현된다. 상태 $x = (q, p)$ , 해밀터니안 $H(q, p) = (1/2) p^T M^{-1}(q) p + V(q)$ 로 정의하면,

$\begin{bmatrix} \dot{q} \\ \dot{p} \end{bmatrix} = \begin{bmatrix} 0 & I \\ -I & 0 \end{bmatrix} \begin{bmatrix} \partial H / \partial q \\ \partial H / \partial p \end{bmatrix} + \begin{bmatrix} 0 \\ I \end{bmatrix} \tau$

$y = \begin{bmatrix} 0 & I \end{bmatrix} \begin{bmatrix} \partial H / \partial q \\ \partial H / \partial p \end{bmatrix} = M^{-1}(q) p = \dot{q}$

이 형태에서 입력은 관절 토크 $\tau$ 이고 출력은 관절 속도 $\dot{q}$ 이며, $\tau^T \dot{q}$ 는 액추에이터가 시스템에 공급하는 일률이다. 이러한 입력-출력 짝은 매니퓰레이터의 수동성 분석과 임피던스 제어 설계의 기초가 된다.

9. 분산 매개변수 시스템으로의 확장

포트-해밀턴 정식화는 유한 차원 시스템에 국한되지 않으며, 편미분 방정식으로 기술되는 분산 매개변수 시스템(distributed parameter system)으로 확장된다. 예를 들어, 유연 보(beam)와 유연 매니퓰레이터, 케이블, 유체-구조 상호작용 시스템 등이 무한 차원 포트-해밀턴 시스템으로 모델링될 수 있다.

분산 매개변수 포트-해밀턴 시스템에서는 슈테켈베르크-디랙(Stokes-Dirac) 구조라는 일반화된 디랙 구조가 사용되며, 경계에서의 에너지 흐름이 자연스럽게 포트로 표현된다. 이러한 정식화는 유연 로봇과 연속체 로봇(continuum robot)의 분석에 적용된다.

10. 제어 설계에서의 활용

포트-해밀턴 정식화는 제어 설계에서 다음과 같은 장점을 제공한다. 첫째, 수동성에 기반한 안정성 보장이 자연스럽게 도출되며, 모델 불확실성과 외란에 대한 견고성이 본질적으로 확보된다. 둘째, 인터커넥션에 대한 닫힘성을 활용하여 제어기 자체를 가상의 포트-해밀턴 시스템으로 설계하고 플랜트와 결합하는 모듈화된 제어 설계가 가능하다. 셋째, 에너지 형성과 댐핑 할당을 통하여 폐루프 시스템의 거동을 직관적으로 결정할 수 있다.

특히 IDA-PBC(Interconnection and Damping Assignment Passivity-Based Control)는 포트-해밀턴 시스템의 구조를 활용한 대표적 제어 설계 방법이며, 부족구동 시스템의 안정화에 효과적이다.

11. 로봇공학에서의 응용

포트-해밀턴 시스템은 로봇공학에서 다음과 같은 응용을 가진다. 첫째, 인간-로봇 상호작용에서 양방향 에너지 교환과 안전성을 보장하는 제어 설계에 활용된다. 둘째, 햅틱 인터페이스의 안정성 분석과 제어에서 통신 지연으로 인한 불안정성을 수동성 관점에서 처리한다. 셋째, 협동 로봇과 다중 로봇 시스템의 분산 제어 설계에 사용된다. 넷째, 유연 로봇과 연속체 로봇의 무한 차원 모델링과 제어에 적용된다. 다섯째, 전기-기계 결합 시스템, 특히 액추에이터의 전기 동역학과 기계 동역학을 통합 모델링하는 데 사용된다.

12. 본 절의 의의

본 절에서 다룬 포트-해밀턴 시스템은 해밀턴 역학을 외부 환경과 상호작용하는 개방 시스템으로 일반화한 강력한 모델링 패러다임이다. 디랙 구조를 통한 추상화, 인터커넥션에 대한 닫힘성, 수동성의 본질적 보장, 그리고 다물리 시스템에 대한 통일된 언어는 포트-해밀턴 시스템을 현대 로봇공학과 메카트로닉스의 핵심 도구로 만든다.

이 정식화는 후속 절에서 다루어질 로봇 시스템 설계 응용의 이론적 기초를 형성하며, 동시에 에너지 기반 제어와 수동성 기반 제어 방법론의 통합된 수학적 틀을 제공한다.

13. 학습 권장사항

본 절의 내용을 충분히 이해하기 위하여 해밀턴 역학의 정준 형식주의, 수동성과 디시패티비티 이론, 그리고 비선형 시스템의 안정성 분석에 대한 선행 학습이 권장된다. 본 절의 내용을 심화 학습하기 위해서는 디랙 구조의 미분 기하학, 본드 그래프 모델링, 분산 매개변수 시스템의 함수 해석학적 처리, 그리고 IDA-PBC 제어 설계에 대한 학습이 권장된다.

14. 참고 문헌

van der Schaft, A. J., and Jeltsema, D. (2014). “Port-Hamiltonian Systems Theory: An Introductory Overview.” Foundations and Trends in Systems and Control, 1(2-3), 173–378.
van der Schaft, A. J. (2017). L2-Gain and Passivity Techniques in Nonlinear Control (3rd ed.). Springer.
Maschke, B. M., and van der Schaft, A. J. (1992). “Port-Controlled Hamiltonian Systems: Modelling Origins and Systemtheoretic Properties.” IFAC Proceedings Volumes, 25(13), 359–365.
Duindam, V., Macchelli, A., Stramigioli, S., and Bruyninckx, H. (Eds.). (2009). Modeling and Control of Complex Physical Systems: The Port-Hamiltonian Approach. Springer.
Ortega, R., van der Schaft, A. J., Maschke, B., and Escobar, G. (2002). “Interconnection and Damping Assignment Passivity-Based Control of Port-Controlled Hamiltonian Systems.” Automatica, 38(4), 585–596.
Stramigioli, S. (2015). “Energy-Aware Robotics.” In Mathematical Control Theory I (pp. 37–50). Springer.
Macchelli, A., Melchiorri, C., and Stramigioli, S. (2009). “Port-Based Modeling and Simulation of Mechanical Systems with Rigid and Flexible Links.” IEEE Transactions on Robotics, 25(5), 1016–1029.

version: 1.0