6.36 고유값 분해의 로봇 동역학 응용

1. 로봇 동역학 방정식의 일반 형식

다관절 로봇의 운동 방정식은 라그랑주 역학에 기반하여 다음의 일반적인 형식으로 표현된다.

$M(\mathbf{q})\ddot{\mathbf{q}} + C(\mathbf{q}, \dot{\mathbf{q}})\dot{\mathbf{q}} + G(\mathbf{q}) = \boldsymbol{\tau}$

여기서 $\mathbf{q} \in \mathbb{R}^n$ 는 일반화 좌표, $M(\mathbf{q})$ 는 대칭 양정치 질량 행렬, $C(\mathbf{q}, \dot{\mathbf{q}})$ 는 코리올리스 및 원심력 행렬, $G(\mathbf{q})$ 는 중력 항, $\boldsymbol{\tau}$ 는 일반화된 입력 토크이다. 이 방정식은 본질적으로 비선형이지만, 정적 평형점 또는 명목 궤적 주변에서의 선형화를 통하여 선형 시간 불변 시스템으로 근사할 수 있다. 이렇게 얻은 선형 시스템의 해석에서 고유값 분해는 핵심적인 역할을 수행한다.

2. 질량 행렬과 강성 행렬의 동시 대각화

평형점 $\mathbf{q}_0$ 에서 운동 방정식을 선형화하면 다음의 형태를 얻는다.

$M_0 \ddot{\delta \mathbf{q}} + K_0 \delta \mathbf{q} = \mathbf{0}$

여기서 $M_0 = M(\mathbf{q}_0)$ 는 대칭 양정치 행렬이고, $K_0$ 는 강성 행렬로서 대칭 양반정치 행렬이다. 이 경우 일반화된 고유값 문제

$K_0 \mathbf{v} = \omega^2 M_0 \mathbf{v}$

를 풀면 고유값 $\omega_i^2$ 와 고유벡터 $\mathbf{v}_i$ 를 얻는다. 이 고유벡터들은 $M_0$ 가중 내적에 대하여 직교한다.

$\mathbf{v}_i^\top M_0 \mathbf{v}_j = \delta_{ij}, \qquad \mathbf{v}_i^\top K_0 \mathbf{v}_j = \omega_i^2 \delta_{ij}$

이 직교성은 두 행렬을 동시에 대각화하는 모달 변환 행렬 $V = (\mathbf{v}_1, \mathbf{v}_2, \ldots, \mathbf{v}_n)$ 의 존재를 보장한다.

3. 모달 좌표 변환

모달 변환 $\delta \mathbf{q} = V \boldsymbol{\eta}$ 를 적용하면 운동 방정식은 다음과 같이 변환된다.

$\ddot{\boldsymbol{\eta}} + \Omega^2 \boldsymbol{\eta} = V^\top \boldsymbol{\tau}_{\delta}, \qquad \Omega^2 = \mathrm{diag}(\omega_1^2, \omega_2^2, \ldots, \omega_n^2)$

이 형식에서 각 모드 $\eta_i$ 는 독립적인 단일 자유도 진동 방정식을 따르며, 모드 사이의 결합이 완전히 제거된다. 모달 좌표계에서의 분리는 다음과 같은 이점을 제공한다.

각 모드의 자연 진동수 $\omega_i$ 가 명시적으로 드러난다.
입력 토크가 각 모드에 미치는 영향을 독립적으로 분석할 수 있다.
감쇠가 비례 감쇠 형태인 경우 감쇠비를 모드별로 계산할 수 있다.
시간 응답이 모드 응답의 일차 결합으로 표현된다.

4. 고유 진동 모드와 노드선

각 고유벡터 $\mathbf{v}_i$ 는 해당 고유 진동수 $\omega_i$ 에서 로봇이 보이는 변형 형상을 나타낸다. 이를 모드 형상(mode shape)이라 한다. 모드 형상에서 진폭이 영이 되는 점이나 선을 노드(node)라 하며, 노드의 위치는 진동 측정 센서와 능동 제동 액추에이터의 배치 결정에 활용된다. 일반적으로 저차 모드는 전체적인 강체 운동에 가까운 형태를 보이는 반면, 고차 모드는 국부적인 변형을 강조하는 경향이 있다.

5. 안정성 해석에서의 고유값 분해

선형화된 상태 공간 모델

$\dot{\mathbf{x}} = A\mathbf{x} + B\mathbf{u}$

에서 시스템 행렬 $A$ 의 고유값은 시스템의 안정성을 결정한다. 모든 고유값이 음의 실수부를 가지면 시스템은 점근적 안정이며, 양의 실수부를 가지는 고유값이 존재하면 불안정이다. 순허수 고유값은 한계 안정 또는 진동 응답을 의미한다.

또한 고유값과 고유벡터로부터 모드 분해된 응답을 얻을 수 있다.

$\mathbf{x}(t) = \sum_{i=1}^{n} c_i e^{\lambda_i t} \boldsymbol{\phi}_i$

여기서 $\boldsymbol{\phi}_i$ 는 우고유벡터, $c_i$ 는 좌고유벡터를 통하여 결정되는 계수이다. 이 분해는 응답의 빠른 모드와 느린 모드를 분리하는 데 활용되며, 모델 차수 축약과 분리 제어 설계의 토대가 된다.

6. 모달 제어 설계

모달 좌표계에서 시스템이 분리되므로 각 모드에 대한 제어기를 독립적으로 설계할 수 있다. 모달 제어(modal control)의 절차는 다음과 같다.

시스템 행렬을 대각화하여 모달 좌표계를 구성한다.
제어 목표에 따라 각 모드에 대한 원하는 폐루프 극점을 설정한다.
각 모드의 폐루프 극점이 원하는 위치가 되도록 모달 이득을 결정한다.
모달 이득을 원래 좌표계로 역변환하여 실제 제어기를 구성한다.

이 접근법은 폴 배치(pole placement) 기법의 자연스러운 형태이며, 시스템의 물리적 의미를 보존하면서 제어기를 설계할 수 있게 한다.

7. 비례 감쇠와 감쇠 모달 해석

감쇠 행렬 $D$ 가 질량 행렬과 강성 행렬의 일차 결합으로 표현되는 경우, 즉 $D = \alpha M + \beta K$ ( $\alpha, \beta \geq 0$ )인 경우, 모달 변환을 통하여 감쇠 행렬도 동시에 대각화된다. 이를 비례 감쇠 또는 레일리 감쇠(Rayleigh damping)라 한다. 이 경우 모달 좌표계에서

$\ddot{\eta}_i + 2\zeta_i \omega_i \dot{\eta}_i + \omega_i^2 \eta_i = f_i(t)$

의 형태를 가지며, 모드별 감쇠비 $\zeta_i$ 를 명시적으로 정의할 수 있다.

8. 차수 축약

고차 시스템의 모든 모드를 제어기 설계에 포함시키는 것은 계산적으로 비효율적이며 수치적으로 불안정할 수 있다. 고유값 분해는 차수 축약(model order reduction)을 위한 자연스러운 도구를 제공한다. 일반적으로 다음의 기준으로 모드를 선별한다.

관심 주파수 대역 내의 모드만을 유지한다.
입력으로부터 잘 제어되고 출력에서 잘 관측되는 모드, 즉 가관측성과 가제어성이 높은 모드를 우선적으로 유지한다.
한켈 특이값(Hankel singular value)이 큰 모드를 보존하는 균형 절단(balanced truncation)을 적용한다.

축약된 모델은 실시간 제어 구현과 실시간 시뮬레이션에서 필수적이다.

9. 유연 매니퓰레이터에 대한 응용

링크가 유연한 매니퓰레이터의 동역학은 일반적으로 분포 매개 변수 시스템으로 모델링되지만, 가정 모드법(assumed mode method)이나 유한 요소법을 통하여 유한 차원 모델로 이산화된다. 이산화된 모델의 질량 및 강성 행렬에 대한 일반화된 고유값 문제를 풀면 유연 진동 모드와 그 진동수가 추출된다. 이를 활용하여 다음의 분석이 수행된다.

작업점 추종 정밀도를 결정하는 주요 진동 모드의 식별
액추에이터와 센서의 비병치(non-collocated) 배치에 따른 영점 분석
잔류 진동 억제를 위한 입력 정형(input shaping) 기법의 설계 매개 변수 결정

10. 다관절 로봇의 동특성 분석

직렬 매니퓰레이터의 질량 행렬은 자세에 따라 변화하므로, 작업 공간의 여러 자세에 대하여 고유값 분해를 수행함으로써 동특성의 변화를 정량적으로 분석할 수 있다. 질량 행렬의 고유값은 일반화 좌표 공간에서의 관성 분포를 나타내며, 가장 큰 고유값은 가장 큰 관성 방향을, 가장 작은 고유값은 가장 작은 관성 방향을 의미한다. 두 고유값의 비는 동적 조작성(dynamic manipulability)의 척도로 활용된다.

11. 비선형 동역학의 국소 해석

비선형 동역학 시스템의 평형점 부근에서의 거동은 야코비 선형화로부터 얻은 시스템 행렬의 고유값에 의하여 국소적으로 결정된다. 하트만-그로브만 정리(Hartman-Grobman theorem)에 의하여, 고유값의 실수부가 모두 영이 아닌 쌍곡 평형점에서는 비선형 시스템이 선형화 시스템과 위상적으로 동등하게 거동한다. 이 결과는 로봇의 평형 자세 안정성 분석과 보행 로봇의 한계 사이클 안정성 평가에서 광범위하게 활용된다.

12. 두 발 로봇의 영점 모멘트 점 동역학

두 발 보행 로봇의 선형 도립진자 모델과 같은 단순화된 동역학 모델에서, 시스템 행렬의 고유값과 고유벡터는 안정적 보행을 위한 영점 모멘트 점(zero moment point) 궤적의 설계에 활용된다. 양의 실수 고유값에 대응하는 모드는 발산하는 거동을 나타내며, 이 모드를 제어 입력으로 안정화하는 것이 균형 유지의 핵심이다.

13. 회전 강체의 자세 동역학

자유 회전 강체의 운동 방정식인 오일러 방정식은 비선형이지만, 평형 회전 상태 주변에서의 선형화를 통하여 자세 안정성을 해석할 수 있다. 관성 텐서의 고유값과 고유벡터, 즉 주관성 모멘트와 주축은 안정적 회전축과 불안정한 회전축을 결정한다. 일반적으로 가장 큰 주관성 모멘트와 가장 작은 주관성 모멘트에 대응하는 주축 주위의 회전은 안정인 반면, 중간 주관성 모멘트의 주축 주위의 회전은 불안정하다. 이를 중간 축 정리(intermediate axis theorem)라 한다.

14. 임피던스 제어와 고유 구조 해석

임피던스 제어에서 원하는 임피던스 행렬은 일반적으로 대칭 양정치 형태로 설계되며, 그 고유값과 고유벡터는 작업 공간에서 로봇이 환경에 대하여 보이는 강성과 감쇠 특성의 방향성을 결정한다. 고유 구조의 적절한 설계를 통하여 특정 방향으로는 부드럽게 반응하면서 다른 방향으로는 단단하게 반응하는 등의 비등방성 임피던스 거동을 구현할 수 있다.

참고문헌

Spong, M. W., Hutchinson, S., & Vidyasagar, M. (2020). Robot Modeling and Control (2nd ed.). Wiley.
Siciliano, B., Sciavicco, L., Villani, L., & Oriolo, G. (2010). Robotics: Modelling, Planning and Control. Springer.
Lynch, K. M., & Park, F. C. (2017). Modern Robotics: Mechanics, Planning, and Control. Cambridge University Press.
Featherstone, R. (2008). Rigid Body Dynamics Algorithms. Springer.
Goldstein, H., Poole, C., & Safko, J. (2014). Classical Mechanics (3rd ed.). Pearson.
Meirovitch, L. (2010). Methods of Analytical Dynamics. Dover Publications.
Khalil, H. K. (2014). Nonlinear Systems (3rd ed.). Pearson.

Version: 1.0