17.11 동역학 운동 방정식의 일반 형식

1. 개요

매니퓰레이터와 다체 로봇 시스템의 운동은 시간에 따라 변하는 관절 변수의 궤적으로 표현되며, 이 궤적이 만족해야 하는 미분 방정식 체계를 동역학 운동 방정식(equations of motion)이라 한다. 정역학 해석이 한 순간의 평형 상태만을 다루는 반면, 동역학은 관성력과 속도 의존 힘을 명시적으로 포함하여 시간에 따른 상태 진화를 기술한다. 로봇 제어, 궤적 계획, 시뮬레이션, 매개변수 식별 등 대부분의 고급 로봇공학 과제는 이 운동 방정식을 출발점으로 삼는다.

본 절에서는 동역학 운동 방정식의 표준적인 일반 형식과 그 구성 요소, 각 항의 물리적 의미, 정식화 방법의 등가성, 그리고 방정식이 지니는 구조적 성질을 학술적으로 기술한다. 구체적 계산 기법과 각 항의 성질에 대한 심층적 논의는 후속 절들에서 개별적으로 다루어진다.

2. 표준 형식의 유도 배경

강체 다체 시스템의 운동 방정식은 뉴턴-오일러 방법, 라그랑주 방법, 해밀턴 방법, 가상 일의 원리 등 여러 정식화 방법을 통하여 유도될 수 있다. 각 방법은 출발점과 중간 단계의 형식이 다르지만, 동일한 물리적 시스템에 대해서는 결과적으로 같은 동역학 방정식을 제공한다. 이는 물리 법칙이 표현 방식과 무관하게 일관된다는 사실의 반영이다.

본 절에서는 $n$ 자유도 강체 매니퓰레이터를 대상으로 하여, 관절 변수 $\mathbf{q} \in \mathbb{R}^n$ 을 일반화 좌표로 채택한 표준 형식을 제시한다. 관절은 회전 관절 또는 병진 관절이며, 각 자유도에 하나의 액추에이터가 부착되어 일반화 힘 $\tau_i$ 를 공급한다고 가정한다. 링크는 강체로 모형화되며, 유연성의 처리는 본 장의 후속 절들에서 별도로 다룬다.

3. 동역학 방정식의 일반 형식

$n$ 자유도 강체 매니퓰레이터의 동역학 운동 방정식은 다음의 표준 형식으로 표현된다.

$\mathbf{M}(\mathbf{q}) \ddot{\mathbf{q}} + \mathbf{C}(\mathbf{q}, \dot{\mathbf{q}}) \dot{\mathbf{q}} + \mathbf{g}(\mathbf{q}) + \mathbf{f}(\dot{\mathbf{q}}) = \boldsymbol{\tau} + \mathbf{J}^T(\mathbf{q}) \mathbf{F}_e$

방정식 각 항의 의미는 다음과 같다. $\mathbf{M}(\mathbf{q}) \in \mathbb{R}^{n \times n}$ 은 관성 행렬(inertia matrix) 또는 질량 행렬(mass matrix)로, 관절 가속도와 그에 필요한 일반화 힘 사이의 선형 관계를 정의한다. $\mathbf{C}(\mathbf{q}, \dot{\mathbf{q}}) \in \mathbb{R}^{n \times n}$ 은 코리올리-원심력 행렬(Coriolis and centrifugal matrix)로, 속도 의존 관성 효과를 포함한다. $\mathbf{g}(\mathbf{q}) \in \mathbb{R}^n$ 은 중력 토크 벡터이며, 각 자세에서 중력이 관절에 유발하는 일반화 힘이다. $\mathbf{f}(\dot{\mathbf{q}}) \in \mathbb{R}^n$ 은 마찰 토크 벡터로, 점성 마찰, 쿨롱 마찰, 그리고 기타 속도 의존 저항을 포함한다. $\boldsymbol{\tau} \in \mathbb{R}^n$ 은 관절 액추에이터가 공급하는 일반화 힘이며, $\mathbf{F}_e \in \mathbb{R}^6$ 은 말단에서 작용하는 외부 렌치, $\mathbf{J}(\mathbf{q}) \in \mathbb{R}^{6 \times n}$ 은 말단 자코비안이다.

이 형식은 정역학 방정식의 자연스러운 확장이며, 정적 극한( $\dot{\mathbf{q}} = \mathbf{0}$ , $\ddot{\mathbf{q}} = \mathbf{0}$ )에서 앞 절들의 정적 평형 방정식으로 환원된다.

관성 행렬 항의 물리적 의미

관성 행렬 $\mathbf{M}(\mathbf{q}) \ddot{\mathbf{q}}$ 는 시스템의 운동 에너지와 관련되는 항으로, 매니퓰레이터가 가속도를 얻기 위하여 요구하는 일반화 힘을 표현한다. 운동 에너지 $K$ 는 다음과 같이 주어진다.

$K(\mathbf{q}, \dot{\mathbf{q}}) = \frac{1}{2} \dot{\mathbf{q}}^T \mathbf{M}(\mathbf{q}) \dot{\mathbf{q}}$

이 이차 형식은 관성 행렬이 양의 정부호라는 사실을 함축한다. 관성 행렬의 원소 $M_{ij}(\mathbf{q})$ 는 일반적으로 관절 변수의 함수이며, 매니퓰레이터의 자세가 변함에 따라 운동 에너지의 기여가 달라진다. 대각 원소 $M_{ii}$ 는 관절 $i$ 의 효과적 관성에 해당하고, 비대각 원소 $M_{ij}$ ( $i \neq j$ )는 관절 $i$ 와 $j$ 사이의 관성적 결합을 나타낸다. 이러한 결합은 한 관절의 가속이 다른 관절에 일반화 힘을 유발함을 의미한다.

4. 코리올리 및 원심력 항의 물리적 의미

$\mathbf{C}(\mathbf{q}, \dot{\mathbf{q}}) \dot{\mathbf{q}}$ 항은 관성 행렬이 자세에 의존한다는 사실로부터 발생하는 속도 의존 관성 효과를 표현한다. 이 항은 두 유형의 물리적 현상을 포함한다. 첫째, 원심력(centrifugal force) 성분은 단일 관절의 회전이 다른 관절 또는 같은 관절에 유발하는 원심 효과이며, 속도 제곱에 비례한다. 둘째, 코리올리 힘(Coriolis force) 성분은 두 관절이 동시에 회전할 때 발생하는 결합 효과이며, 두 관절 속도의 곱에 비례한다.

원심력과 코리올리 힘은 관성 행렬의 요소에 대한 편미분으로 명시적으로 표현되며, 크리스토펠 기호(Christoffel symbol)를 이용한 표현이 표준이다.

$C_{ij}(\mathbf{q}, \dot{\mathbf{q}}) = \sum_{k=1}^{n} c_{ijk}(\mathbf{q}) \dot{q}_k, \qquad c_{ijk} = \frac{1}{2} \left( \frac{\partial M_{ij}}{\partial q_k} + \frac{\partial M_{ik}}{\partial q_j} - \frac{\partial M_{jk}}{\partial q_i} \right)$

이 표현은 코리올리-원심력 행렬이 유일하지 않음을 시사한다. 왜냐하면 서로 다른 $\mathbf{C}$ 행렬이 동일한 $\mathbf{C}(\mathbf{q}, \dot{\mathbf{q}}) \dot{\mathbf{q}}$ 벡터를 재현할 수 있기 때문이다. 크리스토펠 표현은 추가적 성질을 보장하는 대표적 선택이며, 그 성질은 후속 절에서 상세히 다루어진다.

중력 토크 항

$\mathbf{g}(\mathbf{q})$ 항은 매니퓰레이터와 부착물의 중력이 각 관절에 유발하는 일반화 힘이다. 이 항은 오직 자세만의 함수이며, 라그랑주 역학의 위치 에너지 $U(\mathbf{q})$ 로부터 다음과 같이 유도된다.

$\mathbf{g}(\mathbf{q}) = \frac{\partial U(\mathbf{q})}{\partial \mathbf{q}}$

중력 토크의 구체적 계산 방법과 성질은 본 장의 별도 절에서 상세히 다룬다.

5. 마찰 토크 항

$\mathbf{f}(\dot{\mathbf{q}})$ 항은 관절 내부의 마찰과 감쇠 효과를 반영한다. 가장 단순한 모형은 점성 마찰과 쿨롱 마찰의 합으로 표현된다.

$\mathbf{f}(\dot{\mathbf{q}}) = \mathbf{F}_v \dot{\mathbf{q}} + \mathbf{F}_c \, \mathrm{sign}(\dot{\mathbf{q}})$

여기서 $\mathbf{F}_v$ 와 $\mathbf{F}_c$ 는 각각 점성 마찰 계수와 쿨롱 마찰 계수의 대각 행렬이다. 실제 관절 마찰은 더 복잡한 특성을 가지며, Stribeck 효과, 정적 마찰, 히스테리시스 등을 포함할 수 있다. 마찰 모형의 상세는 본 장의 별도 절에서 다루어진다.

외부 렌치 항

$\mathbf{J}^T(\mathbf{q}) \mathbf{F}_e$ 항은 말단 효과기 또는 다른 링크에 작용하는 외부 렌치가 관절에 미치는 일반화 힘이며, 앞 절에서 유도한 정적 힘 관계가 동역학 방정식에 삽입된 형태이다. 여러 접촉점이 있는 경우 이 항은 각 접촉 기여의 합으로 확장된다.

$\sum_{k=1}^{N} \mathbf{J}_k^T(\mathbf{q}) \mathbf{F}_{e, k}$

이 항은 로봇과 환경의 상호작용을 동역학 방정식 내부에서 명시적으로 다룰 수 있게 한다.

6. 정식화 방법의 등가성

같은 시스템에 대한 동역학 방정식은 여러 방법으로 유도될 수 있으며, 다음의 대표적 접근이 있다. 첫째, 라그랑주 방법은 운동 에너지와 위치 에너지로부터 라그랑지안 $L = K - U$ 를 구성하고, 오일러-라그랑주 방정식을 적용하여 운동 방정식을 유도한다. 이 방법은 수식적으로 간결하며, 관성 행렬과 중력 항의 물리적 의미를 명료하게 드러낸다.

둘째, 뉴턴-오일러 방법은 각 링크에 대한 선형 운동량과 각운동량의 변화율에 뉴턴의 제2법칙과 오일러의 회전 방정식을 직접 적용하여 재귀적으로 운동 방정식을 도출한다. 이 방법은 수치 계산에 효율적이며, 관절 반력과 내부 하중을 부수적으로 제공한다.

셋째, 해밀턴 방법은 일반화 좌표와 그 켤레 운동량을 상태 변수로 하여 해밀턴 방정식을 도출한다. 이 방법은 시스템의 심플렉틱 구조를 명시적으로 드러내며, 에너지 보존 성질과 수치 적분에서 유리한 성질을 가진다.

넷째, 가상 일의 원리와 달랑베르 원리를 결합한 방법은 관성력을 외력처럼 다루어 가상 일의 원리를 동역학에 확장한다. 이 방법은 구속 조건이 있는 시스템의 해석에 특히 유용하다.

이들 방법은 서로 다른 계산 절차를 거치지만, 결과로 얻어지는 방정식은 동일한 표준 형식을 가지며 동일한 해를 제공한다.

7. 운동 방정식의 구조적 성질

동역학 방정식의 일반 형식은 다음과 같은 중요한 구조적 성질을 지닌다. 이러한 성질은 제어기 설계와 안정성 해석의 기초를 이룬다.

첫째, 관성 행렬 $\mathbf{M}(\mathbf{q})$ 는 대칭이고 양의 정부호이다. 이는 운동 에너지가 항상 양인 스칼라 형식이라는 물리적 사실의 반영이다. 양의 정부호성은 또한 관성 행렬의 역행렬이 항상 존재함을 보장하여, 가속도 해를 유일하게 결정할 수 있게 한다.

둘째, 특정한 코리올리 행렬 선택 아래에서 $\dot{\mathbf{M}}(\mathbf{q}) - 2 \mathbf{C}(\mathbf{q}, \dot{\mathbf{q}})$ 는 반대칭 행렬이다. 즉 임의의 벡터 $\mathbf{v}$ 에 대하여 다음이 성립한다.

$\mathbf{v}^T \left( \dot{\mathbf{M}}(\mathbf{q}) - 2 \mathbf{C}(\mathbf{q}, \dot{\mathbf{q}}) \right) \mathbf{v} = 0$

이 성질은 리아푸노프 함수 기반 안정성 증명에서 핵심적으로 활용된다.

셋째, 에너지 보존 또는 에너지 흐름의 형식이 성립한다. 외력이 없고 마찰이 없는 경우, 시스템의 총 에너지 $E = K + U$ 는 액추에이터가 공급하는 일률에 의하여 시간적으로 변화하며, 이는 시스템의 수동성(passivity)과 관련된다.

넷째, 운동 방정식은 매개변수에 대하여 선형(linear in parameters)이다. 즉 관성 행렬, 코리올리 행렬, 중력 벡터는 링크의 질량, 질량 중심, 관성 모멘트 등의 동역학 매개변수의 선형 결합으로 표현될 수 있다. 이 성질은 매개변수 식별과 적응 제어에서 결정적으로 활용된다.

상태 공간 형식으로의 변환

제어 이론에서 흔히 사용되는 상태 공간 형식으로 표현하기 위하여 상태 벡터를 $\mathbf{x} = [\mathbf{q}^T, \dot{\mathbf{q}}^T]^T$ 로 정의하고, 운동 방정식을 다음과 같이 재정리한다.

$\dot{\mathbf{x}} = \begin{bmatrix} \dot{\mathbf{q}} \\ \mathbf{M}^{-1}(\mathbf{q}) \left[ \boldsymbol{\tau} + \mathbf{J}^T(\mathbf{q}) \mathbf{F}_e - \mathbf{C}(\mathbf{q}, \dot{\mathbf{q}}) \dot{\mathbf{q}} - \mathbf{g}(\mathbf{q}) - \mathbf{f}(\dot{\mathbf{q}}) \right] \end{bmatrix}$

이 형식은 1차 미분 방정식 시스템의 형태로 수치 적분과 수치적 시뮬레이션에 직접적으로 사용된다. 또한 피드백 선형화, 비선형 제어, 그리고 관측기 설계의 출발점이 된다.

8. 구속 조건의 통합

매니퓰레이터가 환경과 접촉하는 등 구속 조건이 추가된 경우, 라그랑주 승수(Lagrange multiplier)를 도입하여 다음의 형식으로 확장된다.

$\mathbf{M}(\mathbf{q}) \ddot{\mathbf{q}} + \mathbf{C}(\mathbf{q}, \dot{\mathbf{q}}) \dot{\mathbf{q}} + \mathbf{g}(\mathbf{q}) = \boldsymbol{\tau} + \mathbf{A}^T(\mathbf{q}) \boldsymbol{\lambda}$

여기서 $\mathbf{A}(\mathbf{q}) \dot{\mathbf{q}} = \mathbf{0}$ 과 같은 구속 방정식이 함께 주어지며, $\boldsymbol{\lambda}$ 는 구속 힘을 나타내는 라그랑주 승수이다. 이러한 형식은 폐쇄 체인 메커니즘, 고정된 접촉, 그리고 롤링 없이 굴러가는 바퀴 등을 포함하는 다양한 시스템에 적용된다. 구속된 다체 시스템의 구체적 처리는 본 장의 후속 절에서 다루어진다.

운동 방정식의 역학적 해석

운동 방정식을 역으로 활용하면, 원하는 궤적 $\mathbf{q}_d(t)$ 와 그 속도, 가속도가 주어졌을 때 해당 궤적을 실현하기 위한 관절 토크를 직접 계산할 수 있다.

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

이러한 역동역학(inverse dynamics) 계산은 궤적 전향 제어, 계산 토크 제어, 시간 최적 궤적 생성 등에 사용된다. 이와 반대로 순동역학(forward dynamics)은 주어진 토크와 상태에서 가속도를 계산하여 시간 적분을 통해 상태 진화를 시뮬레이션하는 방향이다.

9. 본 절의 의의

본 절에서 제시한 동역학 운동 방정식의 일반 형식은 본 장의 후속 내용과 로봇공학 전반의 동역학 관련 주제를 체계적으로 전개하기 위한 기본 틀이다. 각 항의 상세한 계산, 행렬의 구조적 성질, 그리고 이에 기반한 제어기 설계는 후속 절들에서 각각 심도 있게 다루어질 것이며, 본 절은 그 시작점으로서 전체 방정식이 어떻게 구성되고 어떠한 의미를 가지는지를 정리한다.

또한 이 표준 형식은 로봇 시뮬레이션 소프트웨어, 궤적 계획 알고리즘, 그리고 제어 시스템 설계 도구의 공통 수학적 언어로 기능한다. 이를 정확히 이해하는 것은 이론과 실무 모두에서 필수적이다.

10. 학습 권장사항

본 절의 내용을 충분히 이해하기 위하여 강체 역학, 뉴턴-오일러 방정식, 라그랑주 역학의 기초, 그리고 선형 대수와 미분 방정식에 대한 선행 학습이 권장된다. 본 절의 내용을 심화 학습하기 위해서는 해석 역학, 리 군에 기반한 다체 동역학 정식화, 그리고 수치 적분 이론에 대한 추가 학습이 권장된다.

11. 참고 문헌

Spong, M. W., Hutchinson, S., and Vidyasagar, M. (2020). Robot Modeling and Control (2nd ed.). Wiley.
Siciliano, B., Sciavicco, L., Villani, L., and Oriolo, G. (2009). Robotics: Modelling, Planning and Control. Springer.
Craig, J. J. (2018). Introduction to Robotics: Mechanics and Control (4th ed.). Pearson.
Murray, R. M., Li, Z., and Sastry, S. S. (1994). A Mathematical Introduction to Robotic Manipulation. CRC Press.
Featherstone, R. (2008). Rigid Body Dynamics Algorithms. Springer.
Lynch, K. M., and Park, F. C. (2017). Modern Robotics: Mechanics, Planning, and Control. Cambridge University Press.
Khalil, W., and Dombre, E. (2002). Modeling, Identification and Control of Robots. Hermes Penton.

version: 1.0