32.17 스크류 이론 기반 자코비안 유도

스크류 이론(screw theory)은 강체의 순간 운동과 외력을 스크류 축(screw axis)과 피치(pitch)라는 통일된 기하학적 원소로 기술하는 수학 이론이며, 19세기 Mozzi와 Chasles의 기하학적 정리로부터 Ball의 체계화(1900)로 이어지는 고전적 기원을 가진다. 로봇 기구학에서 스크류 이론은 자코비안의 각 열을 관절 축에 대응하는 스크류의 특정 좌표 표현으로 해석하며, 특히 제품 지수 공식(product of exponentials formula)과 결합하여 데나빗-하르텐버그 매개변수 없이 순기구학과 속도 기구학을 통일적으로 기술하는 현대적 수학 틀을 제공한다. 본 절에서는 스크류, 트위스트, 피치의 학술적 정의, 제품 지수 공식, 공간·물체 자코비안의 스크류 기반 구성, 수반 변환과의 관계를 체계적으로 다룬다.

1. 스크류, 트위스트, 피치의 학술적 정의

스크류(screw)는 3차원 공간의 방향을 가진 직선(스크류 축)과 스칼라 양(피치 $h$ )의 쌍으로 구성된 기하학적 객체이다. 스크류는 순간 나선 운동(instantaneous screw motion)을 기술하며, 스크류 축 주위의 회전과 그 축을 따라가는 병진이 특정 비율로 결합된 합성 운동이다.

트위스트(twist)는 강체의 순간 운동을 기술하는 스크류의 6차원 좌표 표현이다. 선속도 $\vec{v} \in \mathbb{R}^3$ 와 각속도 $\vec{\omega} \in \mathbb{R}^3$ 를 결합하여 다음과 같이 표현된다.

$\mathcal{V} = \begin{bmatrix} \vec{\omega} \\ \vec{v} \end{bmatrix} \in \mathbb{R}^6$

(일부 교재에서는 상·하 블록의 순서를 바꾸어 $[\vec{v}^\top, \vec{\omega}^\top]^\top$ 로 표기하기도 한다. 본 절에서는 Lynch와 Park의 Modern Robotics 규약에 따라 각속도를 상단에 배치한다.)

피치 $h$ 는 스크류의 나선 운동에서 회전 한 바퀴당 축 방향 병진 거리의 비율이다. 트위스트로부터 피치는 다음과 같이 계산된다.

$h = \frac{\vec{\omega} \cdot \vec{v}}{\|\vec{\omega}\|^2} \quad (\vec{\omega} \neq \vec{0})$

피치의 특수한 값에 따라 스크류는 다음과 같이 분류된다.

영 피치( $h = 0$ ): 순수 회전 운동. 회전 관절의 관절 스크류에 대응.
유한 피치( $0 < |h| < \infty$ ): 나사 운동(helical motion). 나사 관절의 관절 스크류에 대응.
무한 피치( $h \to \infty$ , 또는 $\vec{\omega} = \vec{0}$ ): 순수 병진 운동. 직동 관절의 관절 스크류에 대응.

2. 트위스트의 행렬 표현과 리 대수

트위스트는 $4 \times 4$ 반대칭 행렬로도 표현된다. 트위스트 $\mathcal{V} = [\vec{\omega}^\top, \vec{v}^\top]^\top$ 에 대응하는 행렬 표현은 다음과 같다.

$[\mathcal{V}] = \begin{bmatrix} [\vec{\omega}]_\times & \vec{v} \\ \vec{0}^\top & 0 \end{bmatrix} \in \mathfrak{se}(3)$

여기서 $[\vec{\omega}]_\times$ 는 각속도 벡터의 반대칭 행렬이며, $[\mathcal{V}]$ 는 특수 유클리드 군 $SE(3)$ 의 리 대수 $\mathfrak{se}(3)$ 의 원소이다. 이 리 대수적 해석은 스크류 이론을 현대 미분 기하학의 리 군-리 대수 언어와 직접 연결한다.

동차 변환 행렬 $\mathbf{T}(t) \in SE(3)$ 의 시간 미분과 트위스트 사이에는 다음 관계가 성립한다.

$\dot{\mathbf{T}}(t) \, \mathbf{T}^{-1}(t) = [\mathcal{V}^s](t)$

여기서 $\mathcal{V}^s$ 는 공간 좌표계에서 표현된 공간 트위스트(spatial twist)이며, 물체 좌표계 표현의 물체 트위스트(body twist) $\mathcal{V}^b$ 와는 $\mathbf{T}^{-1} \dot{\mathbf{T}} = [\mathcal{V}^b]$ 의 관계로 연결된다.

3. 제품 지수 공식

제품 지수 공식(product of exponentials formula, PoE)은 직렬 매니퓰레이터의 순기구학을 관절 스크류와 관절 변수의 지수 사상 곱으로 표현하는 현대적 수학 틀이다. 기본 구성은 다음과 같다.

각 관절 $i$ 에 대해 기저 좌표계에서 표현된 관절 스크류(screw axis)를 $\mathcal{S}_i = [\vec{\omega}_i^\top, \vec{v}_i^\top]^\top \in \mathbb{R}^6$ 로 정의한다. 회전 관절의 경우 관절 축 방향 단위 벡터 $\hat{\omega}_i$ 와 관절 축 위의 한 점 $\vec{q}_{i,0}$ 을 이용하여

$\mathcal{S}_i = \begin{bmatrix} \hat{\omega}_i \\ -\hat{\omega}_i \times \vec{q}_{i,0} \end{bmatrix}$

로 주어지며, 직동 관절의 경우 관절 축 방향 단위 벡터 $\hat{v}_i$ 를 이용하여

$\mathcal{S}_i = \begin{bmatrix} \vec{0} \\ \hat{v}_i \end{bmatrix}$

로 주어진다. 이러한 관절 스크류는 영 구성(home configuration) $\vec{q} = \vec{0}$ 에서 고정된 값으로 평가된다.

엔드 이펙터의 영 구성에서의 동차 변환을 $\mathbf{M} \in SE(3)$ 로 표기하면, 순기구학은 다음의 제품 지수 공식으로 표현된다.

$\mathbf{T}(\vec{q}) = e^{[\mathcal{S}_1] q_1} \, e^{[\mathcal{S}_2] q_2} \cdots e^{[\mathcal{S}_n] q_n} \, \mathbf{M}$

여기서 $e^{[\mathcal{S}_i] q_i} \in SE(3)$ 은 리 대수 $\mathfrak{se}(3)$ 의 원소 $[\mathcal{S}_i] q_i$ 의 지수 사상이며, 닫힌 형식으로 계산 가능하다. 이 공식은 Brockett(1984)에 의해 로봇 공학에 도입되었으며, 데나빗-하르텐버그 매개변수의 모호성 없이 순기구학을 기술하는 장점을 가진다.

32.17.4 공간 자코비안의 스크류 기반 유도

공간 자코비안(space Jacobian) $\mathbf{J}_s(\vec{q}) \in \mathbb{R}^{6 \times n}$ 은 공간 좌표계에서 표현된 엔드 이펙터의 공간 트위스트 $\mathcal{V}^s$ 와 관절 속도를 연결한다.

$\mathcal{V}^s = \mathbf{J}_s(\vec{q}) \, \dot{\vec{q}}$

제품 지수 공식의 시간 미분을 통해 공간 자코비안의 각 열이 유도된다. $\dot{\mathbf{T}} \mathbf{T}^{-1}$ 의 계산 과정에서 각 관절 스크류가 이전 관절들의 지수 사상을 통해 수반 변환되는 구조가 나타나며, 이는 다음의 닫힌 형식으로 정리된다.

$\vec{J}_{s,i}(\vec{q}) = [\text{Ad}_{e^{[\mathcal{S}_1] q_1} \cdots e^{[\mathcal{S}_{i-1}] q_{i-1}}}] \, \mathcal{S}_i$

여기서 $[\text{Ad}_\mathbf{T}] \in \mathbb{R}^{6 \times 6}$ 은 동차 변환 $\mathbf{T} = (\mathbf{R}, \vec{p})$ 에 대응하는 수반 사상으로, 다음과 같이 정의된다.

$[\text{Ad}_\mathbf{T}] = \begin{bmatrix} \mathbf{R} & \mathbf{0} \\ [\vec{p}]_\times \mathbf{R} & \mathbf{R} \end{bmatrix}$

이 표현은 $\mathcal{V} = [\vec{\omega}^\top, \vec{v}^\top]^\top$ 규약에 해당하며, 각속도 블록을 상단에 배치할 때 성립한다. $\vec{J}_{s,1} = \mathcal{S}_1$ 이고, $i \geq 2$ 의 경우 관절 $i$ 이전의 관절 운동에 의해 변환된 관절 스크류가 공간 자코비안의 $i$ 번째 열이 된다.

4. 물체 자코비안의 스크류 기반 유도

물체 자코비안(body Jacobian) $\mathbf{J}_b(\vec{q}) \in \mathbb{R}^{6 \times n}$ 은 물체 좌표계(엔드 이펙터 좌표계)에서 표현된 물체 트위스트 $\mathcal{V}^b$ 와 관절 속도를 연결한다.

$\mathcal{V}^b = \mathbf{J}_b(\vec{q}) \, \dot{\vec{q}}$

제품 지수 공식의 $\mathbf{T}^{-1} \dot{\mathbf{T}}$ 계산 과정에서 각 관절 스크류가 이후 관절들의 지수 사상의 역으로 수반 변환되는 구조가 나타나며, 물체 자코비안의 열 벡터는 다음의 닫힌 형식으로 정리된다.

$\vec{J}_{b,i}(\vec{q}) = [\text{Ad}_{e^{-[\mathcal{B}_n] q_n} \cdots e^{-[\mathcal{B}_{i+1}] q_{i+1}}}] \, \mathcal{B}_i$

여기서 $\mathcal{B}_i$ 는 물체 좌표계에서 표현된 관절 스크류이며, 공간 좌표계의 관절 스크류와는 $\mathcal{B}_i = [\text{Ad}_{\mathbf{M}^{-1}}] \mathcal{S}_i$ 의 관계를 가진다. $\vec{J}_{b,n} = \mathcal{B}_n$ 이고, $i < n$ 의 경우 관절 $i$ 이후의 관절 운동의 역변환이 적용된다.

공간 자코비안과 물체 자코비안은 동차 변환 $\mathbf{T}(\vec{q})$ 의 수반 변환에 의해 직접 연결된다.

$\mathbf{J}_s(\vec{q}) = [\text{Ad}_{\mathbf{T}(\vec{q})}] \, \mathbf{J}_b(\vec{q}), \qquad \mathbf{J}_b(\vec{q}) = [\text{Ad}_{\mathbf{T}^{-1}(\vec{q})}] \, \mathbf{J}_s(\vec{q})$

32.17.6 수반 변환의 학술적 역할

수반 변환 $[\text{Ad}_\mathbf{T}]$ 은 스크류 이론 기반 자코비안 유도의 중심 연산이다. 수반 변환의 주요 성질은 다음과 같다.

리 군 준동형: $[\text{Ad}_{\mathbf{T}_1 \mathbf{T}_2}] = [\text{Ad}_{\mathbf{T}_1}] [\text{Ad}_{\mathbf{T}_2}]$ .

역변환: $[\text{Ad}_\mathbf{T}]^{-1} = [\text{Ad}_{\mathbf{T}^{-1}}]$ .

트위스트 좌표 변환: 한 좌표계에서 표현된 트위스트를 다른 좌표계에서 표현된 트위스트로 변환한다.

수반 변환은 자코비안의 좌표계 변환에서도 핵심적 역할을 한다. 좌표계 $\{A\}$ 에서 $\{B\}$ 로의 동차 변환을 ${}^B\mathbf{T}_A$ 라 할 때, 자코비안의 좌표계 변환은 다음과 같다.

${}^B\mathbf{J} = [\text{Ad}_{{}^B\mathbf{T}_A}] \, {}^A\mathbf{J}$

이는 선속도와 각속도 블록을 독립적으로 변환하는 단순 회전 행렬 적용과 다르며, 병진 오프셋에 기인하는 블록 간 결합이 반영된다.

5. 스크류 기반 자코비안의 수치적 계산

스크류 이론 기반 자코비안의 수치적 계산 절차는 다음과 같다.

단계 1: 관절 스크류 $\mathcal{S}_i$ 와 영 구성 엔드 이펙터 변환 $\mathbf{M}$ 을 매니퓰레이터의 기하학적 설계로부터 사전 결정한다.

단계 2: 관절 변수 $\vec{q}$ 가 주어지면 각 관절의 지수 사상 $e^{[\mathcal{S}_i] q_i}$ 을 Rodrigues 공식의 확장으로 계산한다.

단계 3: 공간 자코비안의 경우 누적 수반 변환 $[\text{Ad}_{e^{[\mathcal{S}_1] q_1} \cdots e^{[\mathcal{S}_{i-1}] q_{i-1}}}]$ 을 관절별로 재귀 계산하며, 각 관절 스크류에 적용하여 $\vec{J}_{s,i}$ 를 얻는다.

단계 4: 물체 자코비안의 경우 엔드 이펙터로부터 기저 방향으로 역재귀적 수반 변환을 계산하여 $\vec{J}_{b,i}$ 를 얻는다.

단계 5: 계산된 열 벡터들을 가로로 배열하여 자코비안 행렬을 조립한다.

이 알고리즘의 시간 복잡도는 $\mathcal{O}(n)$ 이며, 기하학적 자코비안의 닫힌 형식 재귀와 동일한 효율을 가진다. 현대 로봇 공학 라이브러리(Pinocchio, Drake 등)는 스크류 기반 계산을 표준적으로 채택한다.

6. 기하학적 자코비안과의 동등성

스크류 기반 자코비안과 32.6절에서 유도된 기하학적 자코비안은 동일한 수학적 객체의 서로 다른 좌표 표현이다. 구체적으로 $[\vec{\omega}^\top, \vec{v}^\top]^\top$ 규약의 공간 자코비안을 $[\vec{v}^\top, \vec{\omega}^\top]^\top$ 규약으로 성분 순서만 교환하면 기하학적 자코비안의 표현과 일치한다.

이러한 동등성은 회전 관절의 공간 자코비안 열 벡터가 $\vec{J}_{s,i} = [\hat{\omega}_i^\top, (-\hat{\omega}_i \times \vec{q}_i(\vec{q}))^\top]^\top$ 의 형태를 가지며, 관절 축 위의 점 $\vec{q}_i(\vec{q})$ 와 $\hat{\omega}_i$ 의 외적 관계가 기하학적 자코비안의 $\hat{z}_{i-1} \times (\vec{p}_n - \vec{p}_{i-1})$ 구조와 기하학적으로 일치함에서 확인된다. 두 표현은 수학적으로 동치이며, 어느 형태를 채택할지는 해석의 편의와 응용 맥락에 따른다.

7. 스크류 기반 유도의 학술적 장점

스크류 이론 기반 자코비안 유도는 다음의 학술적 장점을 제공한다.

데나빗-하르텐버그 매개변수로부터의 독립성: 제품 지수 공식은 관절 스크류와 영 구성 변환만으로 매니퓰레이터의 순기구학을 기술하므로, 데나빗-하르텐버그 규약의 비유일성과 모호성을 피할 수 있다.

기하학적 투명성: 각 자코비안 열이 관절 스크류의 수반 변환으로 명시적으로 표현되므로 기하학적 해석이 명료하다.

리 군-리 대수 구조와의 통합: $SE(3)$ 와 $\mathfrak{se}(3)$ 의 리 군-리 대수 구조에 자연스럽게 대응하므로, 현대 미분 기하학 기반 제어 이론과 직접 연결된다.

좌표계 변환의 체계성: 수반 변환을 통한 좌표계 변환이 일관된 수학적 틀로 수행된다.

확장성: 스크류 이론은 강체 동역학, 병렬 기구, 연성체 로봇 등 다양한 분야로 일관되게 확장된다.

병렬 기구와의 결합: 병렬 기구의 기구학 해석에서 스크류 이론은 구속 조건의 체계적 기술을 제공한다.

8. 학술적 의의

본 절에서 정립한 스크류 이론 기반 자코비안 유도는 로봇 공학에서 다음의 학술적 의의를 가진다. 첫째, 스크류·트위스트·피치의 통일된 기하학적 객체를 통해 회전 관절과 직동 관절의 관절 운동을 하나의 수학 틀로 기술한다. 둘째, 제품 지수 공식은 순기구학의 현대적 표현으로 데나빗-하르텐버그 매개변수의 모호성을 해소한다. 셋째, 수반 변환을 매개로 한 공간 자코비안과 물체 자코비안의 동등한 구성은 좌표계 선택의 체계적 취급을 가능하게 한다. 넷째, 리 군 $SE(3)$ 와 리 대수 $\mathfrak{se}(3)$ 의 언어로 자코비안을 기술함으로써 현대 미분 기하학 기반 기구학·제어 이론과의 자연스러운 통합을 제공한다. 다섯째, 스크류 기반 계산 알고리즘은 $\mathcal{O}(n)$ 복잡도로 기하학적 자코비안과 동등한 효율을 가지며, 현대 로봇 공학 소프트웨어의 표준적 구현 기반이 된다. 여섯째, 본 절의 내용은 후속 절 트위스트 표현과 자코비안의 관계, 공간·물체 자코비안 구분의 이론적 토대를 이룬다.

9. 출처

Ball, R. S., A Treatise on the Theory of Screws, Cambridge University Press, 1900.
Brockett, R. W., “Robotic manipulators and the product of exponentials formula”, Mathematical Theory of Networks and Systems, Springer, pp. 120–129, 1984.
Murray, R. M., Li, Z., and Sastry, S. S., A Mathematical Introduction to Robotic Manipulation, CRC Press, 1994.
Lynch, K. M. and Park, F. C., Modern Robotics: Mechanics, Planning, and Control, Cambridge University Press, 2017.
Selig, J. M., Geometric Fundamentals of Robotics, 2nd edition, Springer, 2005.
Park, F. C., “Computational aspects of the product-of-exponentials formula for robot kinematics”, IEEE Transactions on Automatic Control, Vol. 39, No. 3, pp. 643–647, 1994.
Featherstone, R., Rigid Body Dynamics Algorithms, Springer, 2008.
Duffy, J., “The fallacy of modern hybrid control theory that is based on ‘orthogonal complements’ of twist and wrench spaces”, Journal of Robotic Systems, Vol. 7, No. 2, pp. 139–144, 1990.

10. 버전

문서 버전: 1.1
작성일: 2026-04-19