15.28 에너지 보존과 리아푸노프 안정성 해석
1. 개요
매니퓰레이터의 동역학에서 에너지 보존은 시스템의 본질적 성질이며, 리아푸노프 안정성 해석의 기반이다. 에너지 함수가 리아푸노프 함수의 후보로 사용되어 매니퓰레이터의 안정성을 증명하는 데 활용된다. 본 절에서는 에너지 보존과 리아푸노프 안정성 해석을 다룬다.
2. 매니퓰레이터의 에너지 보존
2.1 총 에너지
매니퓰레이터의 총 에너지는 운동 에너지와 위치 에너지의 합이다.
E(\mathbf{q}, \dot{\mathbf{q}}) = T(\mathbf{q}, \dot{\mathbf{q}}) + U(\mathbf{q}) = \frac{1}{2}\dot{\mathbf{q}}^T\mathbf{M}(\mathbf{q})\dot{\mathbf{q}} + U(\mathbf{q})
2.2 에너지 변화율
총 에너지의 시간 변화율은 다음과 같다.
\dot{E} = \dot{\mathbf{q}}^T\mathbf{M}\ddot{\mathbf{q}} + \frac{1}{2}\dot{\mathbf{q}}^T\dot{\mathbf{M}}\dot{\mathbf{q}} + \dot{\mathbf{q}}^T\frac{\partial U}{\partial \mathbf{q}}
2.3 동역학 방정식의 활용
매니퓰레이터의 동역학 방정식 \mathbf{M}\ddot{\mathbf{q}} = \boldsymbol{\tau} - \mathbf{C}\dot{\mathbf{q}} - \mathbf{g}를 대입하면
\dot{E} = \dot{\mathbf{q}}^T(\boldsymbol{\tau} - \mathbf{C}\dot{\mathbf{q}} - \mathbf{g}) + \frac{1}{2}\dot{\mathbf{q}}^T\dot{\mathbf{M}}\dot{\mathbf{q}} + \dot{\mathbf{q}}^T\mathbf{g}
= \dot{\mathbf{q}}^T\boldsymbol{\tau} + \frac{1}{2}\dot{\mathbf{q}}^T(\dot{\mathbf{M}} - 2\mathbf{C})\dot{\mathbf{q}}
2.4 반대칭성의 활용
\dot{\mathbf{M}} - 2\mathbf{C}가 반대칭이므로 두 번째 항이 0이다.
\dot{\mathbf{q}}^T(\dot{\mathbf{M}} - 2\mathbf{C})\dot{\mathbf{q}} = 0
따라서
\dot{E} = \dot{\mathbf{q}}^T\boldsymbol{\tau}
2.5 에너지 보존
외력이 없으면 (\boldsymbol{\tau} = 0) 총 에너지가 보존된다.
\dot{E} = 0 \Rightarrow E = \text{const}
3. 일률 정리
3.1 에너지 변화와 일
총 에너지의 변화율은 일반화 힘이 일반화 속도에 대해 한 일률과 같다.
\dot{E} = \dot{\mathbf{q}}^T\boldsymbol{\tau}
3.2 적분 형식
이를 시간에 대해 적분하면
\Delta E = \int_{t_1}^{t_2}\dot{\mathbf{q}}^T\boldsymbol{\tau}\, dt = W
여기서 W는 일반화 힘에 의한 일이다.
3.3 물리적 해석
일률 정리는 에너지의 변화가 외부에서 한 일과 같다는 물리적 사실을 표현한다.
4. 리아푸노프 함수
4.1 정의
리아푸노프 함수 V(\mathbf{x})는 다음의 성질을 만족하는 스칼라 함수이다.
- V(\mathbf{x}_e) = 0 (평형점에서 0)
- V(\mathbf{x}) > 0, \forall \mathbf{x} \neq \mathbf{x}_e (다른 곳에서 양수)
- \dot{V}(\mathbf{x}) \leq 0 (감소 또는 일정)
이러한 함수가 존재하면 평형점이 안정하다(또는 점근적으로 안정).
4.2 리아푸노프 정리
리아푸노프 정리는 안정성 분석의 기본 결과이다.
- \dot{V} \leq 0: 평형점이 안정
- \dot{V} < 0 (평형점 외에서): 평형점이 점근적으로 안정
5. 매니퓰레이터의 안정성 해석
5.1 정적 평형점
매니퓰레이터의 정적 평형점은 다음을 만족한다.
\dot{\mathbf{q}}_e = 0, \quad \mathbf{g}(\mathbf{q}_e) = \boldsymbol{\tau}
5.2 리아푸노프 함수의 후보
리아푸노프 함수의 후보로 다음의 형태가 자주 사용된다.
V = \frac{1}{2}\dot{\mathbf{q}}^T\mathbf{M}(\mathbf{q})\dot{\mathbf{q}} + \tilde{U}(\mathbf{q})
여기서 \tilde{U}는 평형점에서 0이 되는 위치 에너지의 변형이다.
5.3 PD + 중력 보상의 안정성
PD 제어 + 중력 보상의 매니퓰레이터를 고려한다.
\boldsymbol{\tau} = \mathbf{g}(\mathbf{q}) + \mathbf{K}_p(\mathbf{q}_d - \mathbf{q}) - \mathbf{K}_v\dot{\mathbf{q}}
리아푸노프 함수 후보
V = \frac{1}{2}\dot{\mathbf{q}}^T\mathbf{M}(\mathbf{q})\dot{\mathbf{q}} + \frac{1}{2}(\mathbf{q}_d - \mathbf{q})^T\mathbf{K}_p(\mathbf{q}_d - \mathbf{q})
미분하여 \dot{V} = -\dot{\mathbf{q}}^T\mathbf{K}_v\dot{\mathbf{q}} \leq 0임을 보일 수 있다. 이는 평형점이 안정적임을 보여준다.
5.4 점근적 안정성
추가 조건(예: 라살의 불변 원리)을 통해 점근적 안정성도 증명될 수 있다.
6. 라살의 불변 원리
6.1 진술
라살의 불변 원리는 \dot{V} \leq 0인 경우에 점근적 안정성을 증명하는 도구이다.
\Omega = \{\mathbf{x}: V(\mathbf{x}) \leq c\}
가 양의 불변 집합이라면 \Omega 안의 가장 큰 불변 집합 M이 존재한다. 모든 해는 M으로 수렴한다.
6.2 매니퓰레이터의 응용
PD + 중력 보상의 경우 \dot{V} = 0이려면 \dot{\mathbf{q}} = 0이어야 한다. 이는 시스템이 평형점에 수렴함을 의미한다.
7. 응용
7.1 제어 안정성 증명
매니퓰레이터의 다양한 제어 알고리즘의 안정성 증명에 에너지 함수와 리아푸노프 분석이 사용된다.
7.2 적응 제어
적응 제어의 안정성 증명에서 에너지 함수를 확장한 리아푸노프 함수가 사용된다.
7.3 임피던스 제어
임피던스 제어의 수동성 증명에 에너지 분석이 활용된다.
8. 본 절의 의의
본 절은 에너지 보존과 리아푸노프 안정성 해석을 다루었다. 에너지 함수는 매니퓰레이터의 안정성 분석의 핵심 도구이며, 다양한 제어 기법의 안정성 증명에 활용된다.
9. 참고 문헌
- Spong, M. W., Hutchinson, S., & Vidyasagar, M. (2020). Robot Modeling and Control (2nd ed.). Wiley.
- Khalil, H. K. (2002). Nonlinear Systems (3rd ed.). Prentice Hall.
- Slotine, J.-J. E., & Li, W. (1991). Applied Nonlinear Control. Prentice Hall.
- Takegaki, M., & Arimoto, S. (1981). A new feedback method for dynamic control of manipulators. ASME Journal of Dynamic Systems, Measurement, and Control, 103(2), 119-125.
version: 1.0