6.80 스크류 운동과 강체 변위의 표현

1. 도입

강체의 임의의 변위는 회전과 병진의 결합으로 표현되며, 일반적으로는 두 성분이 서로 독립적인 자유도로 다루어진다. 그러나 이 두 성분은 단일한 기하학적 운동, 즉 스크류 운동(screw motion)으로 통합 표현될 수 있다는 사실이 19세기 이래로 알려져 있다. 이 결과는 챠슬(Chasles)의 정리와 모치(Mozzi)의 정리로 알려진 고전적 명제이며, 임의의 강체 변위는 어떤 직선(스크류 축)을 중심으로 한 회전과 그 직선 방향으로의 병진의 결합으로 일의적으로 표현될 수 있음을 진술한다. 이러한 표현은 강체 운동을 단일 축에 대한 운동으로 환원함으로써 운동학과 동역학의 기술을 본질적으로 단순화하며, 매니퓰레이터 운동학의 지수 사상 표현, 트위스트 좌표, 리 군 위의 운동학적 형식화의 수학적 기반을 이룬다. 본 절에서는 스크류 운동의 정의, 스크류 축과 피치, 스크류 변위의 행렬 표현, 트위스트와의 관계, 그리고 로봇 공학에서의 활용을 학술적으로 기술한다.

2. 스크류 운동의 정의

2.1 스크류 운동의 기하학적 정의

스크류 운동은 어떤 직선(스크류 축) \ell 을 중심으로 한 회전과 동일한 직선 방향으로의 병진을 동시에 적용하는 강체 운동이다. 회전각이 \theta 이고 회전축 방향으로의 병진 거리가 d 일 때, 두 양의 비

h = \frac{d}{\theta}

를 스크류의 피치(pitch)라 한다. 피치는 회전 단위 각도당 진행하는 병진 거리의 척도이며, 나사못이 한 바퀴 돌 때 진행하는 거리에 비유할 수 있다.

2.2 챠슬-모치 정리

임의의 강체 변환 T \in SE(3) 는 다음 두 가지 경우 중 하나로 표현된다.

첫째, 회전 부분이 항등이 아닌 경우. 이 경우 유일한 스크류 축 \ell, 회전각 \theta, 그리고 피치 h 가 존재하여 T\ell 을 중심으로 한 각도 \theta 의 회전과 \ell 방향으로 거리 h \theta 의 병진의 결합으로 표현된다.

둘째, 회전 부분이 항등인 경우. 이 경우 강체 변환은 순수 병진이며, 피치를 무한대로 정의하고 스크류 축을 병진 방향으로 잡으면 된다.

이 정리는 모든 강체 변위가 본질적으로 단일 스크류 운동임을 보여주며, 6차원 자유도를 가지는 강체 변위가 (스크류 축의 4 자유도) + (회전각의 1 자유도) + (피치의 1 자유도) = 6 자유도로 분해됨을 보장한다.

3. 스크류 축과 스크류 매개변수

3.1 스크류 축의 기술

3차원 공간에서 직선은 한 점과 단위 방향 벡터로 기술된다. 스크류 축 \ell 의 단위 방향 벡터를 \hat{\mathbf{s}} \in \mathbb{R}^{3} 라 하고, 축 위의 한 점을 \mathbf{q} \in \mathbb{R}^{3} 라 하면, 스크류 축은 다음과 같이 매개변수화된다.

\ell = \{ \mathbf{q} + \lambda \hat{\mathbf{s}} \mid \lambda \in \mathbb{R} \}

이 표현은 스크류 축의 위치와 방향을 모두 담고 있으며, 동일한 직선이 무한히 많은 (\mathbf{q}, \hat{\mathbf{s}}) 쌍으로 표현될 수 있다는 모호성은 \mathbf{q} 를 축 위의 가장 가까운 점(원점에서의 수직 사영)으로 잡음으로써 제거할 수 있다.

3.2 스크류 매개변수

하나의 스크류 운동은 다음 네 가지 정보로 완전히 기술된다.

첫째, 스크류 축의 단위 방향 벡터 \hat{\mathbf{s}}. 둘째, 스크류 축 위의 한 점 \mathbf{q}. 셋째, 회전각 \theta. 넷째, 피치 h. 이 네 정보의 자유도를 모두 합하면 2 + 2 + 1 + 1 = 6 이며, 강체 변환의 자유도와 일치한다. 여기서 단위 방향 벡터의 자유도가 2인 이유는 단위 길이 제약 때문이고, 점의 자유도가 2인 이유는 축 위에서의 위치 조정이 결과에 영향을 주지 않기 때문이다.

4. 스크류 변위의 행렬 표현

4.1 회전과 병진의 결합 형태

스크류 운동에 대응하는 강체 변환을 명시적으로 구성하면 다음과 같다. 스크류 축이 원점을 지나가는 경우, 즉 \mathbf{q} = \mathbf{0} 인 단순한 경우에는

T = \begin{bmatrix} R(\hat{\mathbf{s}}, \theta) & h \theta \, \hat{\mathbf{s}} \\ \mathbf{0}^{T} & 1 \end{bmatrix}

으로 표현된다. 여기서 R(\hat{\mathbf{s}}, \theta) 는 단위 벡터 \hat{\mathbf{s}} 를 중심으로 각도 \theta 만큼 회전하는 회전 행렬이며, 로드리게스 공식으로 닫힌 형태로 주어진다.

4.2 일반적인 스크류 축의 경우

스크류 축이 원점을 지나가지 않는 일반적인 경우에는, 스크류 축 위의 점 \mathbf{q} 를 회전 중심으로 사용하여 다음과 같이 변환을 구성한다.

T = \begin{bmatrix} R & (I - R) \mathbf{q} + h \theta \, \hat{\mathbf{s}} \\ \mathbf{0}^{T} & 1 \end{bmatrix}

여기서 R = R(\hat{\mathbf{s}}, \theta) 이다. 병진 부분의 첫 번째 항 (I - R) \mathbf{q} 는 회전 중심이 원점이 아닌 데에서 비롯되는 보정 항이며, 두 번째 항 h \theta \hat{\mathbf{s}} 는 스크류 축 방향으로의 진행을 나타낸다.

4.3 강체 변환으로부터 스크류 매개변수의 추출

거꾸로 임의의 강체 변환 T = (R, \mathbf{t}) 가 주어졌을 때, 그에 대응하는 스크류 매개변수를 추출하는 절차도 닫힌 형태로 주어진다. 회전 부분 R 로부터 회전축 \hat{\mathbf{s}} 와 회전각 \theta 를 추출한 뒤, 병진 부분 \mathbf{t} 를 스크류 축 방향 성분과 수직 성분으로 분해하여 피치 h 와 축 위의 점 \mathbf{q} 를 결정한다. 이 절차는 강체 변환을 단일 스크류 운동의 형태로 정규화하는 과정이며, 운동의 본질적 기하학적 의미를 드러낸다.

5. 트위스트 표현과 지수 사상

5.1 트위스트의 정의

스크류 운동의 무한소 형태는 트위스트(twist)로 표현된다. 6차원 트위스트 벡터는 다음과 같이 정의된다.

\xi = \begin{bmatrix} \mathbf{v} \\ \boldsymbol{\omega} \end{bmatrix} \in \mathbb{R}^{6}

여기서 \boldsymbol{\omega} 는 각속도 부분이고 \mathbf{v} 는 선속도 부분이다. 트위스트는 리 대수 \mathfrak{se}(3) 의 원소를 6차원 벡터로 식별한 것이며, 강체의 순간 운동을 기술한다.

5.2 단위 트위스트와 스크류 축

회전 성분이 영이 아닌 경우, 트위스트는 단위 각속도를 갖도록 정규화될 수 있다. 단위 트위스트는 스크류 축의 매개변수와 다음과 같은 관계를 가진다.

\hat{\xi} = \begin{bmatrix} \mathbf{v} \\ \hat{\mathbf{s}} \end{bmatrix}, \qquad \mathbf{v} = -\hat{\mathbf{s}} \times \mathbf{q} + h \, \hat{\mathbf{s}}

여기서 첫 번째 항 -\hat{\mathbf{s}} \times \mathbf{q} 는 스크류 축이 원점에서 떨어진 데에서 비롯되는 모멘트 항이며, 두 번째 항 h \hat{\mathbf{s}} 는 피치에 의한 축 방향 성분이다.

5.3 지수 사상을 통한 강체 변환 생성

트위스트로부터 강체 변환은 행렬 지수 함수를 통하여 생성된다.

T = \exp(\theta \, \hat{\xi}^{\wedge}) = \exp\!\left( \theta \begin{bmatrix} [\hat{\mathbf{s}}]_{\times} & \mathbf{v} \\ \mathbf{0}^{T} & 0 \end{bmatrix} \right)

이 지수의 닫힌 형태는 회전 부분에서는 로드리게스 공식, 병진 부분에서는 그 일반화 형태로 주어진다.

\exp(\theta \hat{\xi}^{\wedge}) = \begin{bmatrix} e^{\theta [\hat{\mathbf{s}}]_{\times}} & G(\theta) \mathbf{v} \\ \mathbf{0}^{T} & 1 \end{bmatrix}

G(\theta) = I \theta + (1 - \cos \theta) [\hat{\mathbf{s}}]_{\times} + (\theta - \sin \theta) [\hat{\mathbf{s}}]_{\times}^{2}

이 공식은 단위 트위스트 \hat{\xi} 와 스칼라 매개변수 \theta 만으로 임의의 강체 변환을 생성할 수 있음을 보여주며, 스크류 운동의 수학적 본질을 명시적으로 드러낸다.

5.4 역사상과 강체 변환의 로그

지수 사상의 역사상은 강체 변환을 트위스트로 분해한다. 즉 임의의 T \in SE(3) 에 대하여 T = \exp(\xi^{\wedge}) 를 만족하는 \xi \in \mathfrak{se}(3) 를 닫힌 형태로 구할 수 있으며, 이 절차가 곧 챠슬-모치 정리의 대수적 구현이다.

6. 플뤼커 좌표와의 관계

6.1 직선의 플뤼커 좌표

스크류 축은 3차원 공간의 직선이며, 직선은 6차원 플뤼커 좌표(Plücker coordinates)로 표현될 수 있다. 직선 \ell 위의 한 점 \mathbf{q} 와 단위 방향 \hat{\mathbf{s}} 가 주어졌을 때, 플뤼커 좌표는 다음과 같이 정의된다.

(\hat{\mathbf{s}}, \, \mathbf{q} \times \hat{\mathbf{s}})

이 6차원 표현은 직선의 위치와 방향을 동시에 담고 있으며, 한 점 선택의 자유에 대해 불변이다. 트위스트의 단위 형태 \hat{\xi} 는 본질적으로 플뤼커 좌표의 자연스러운 일반화이며, 피치 항을 추가하여 6차원 표현 안에 회전과 병진의 결합 정보를 담는다.

6.2 스크류 이론의 통일적 관점

이러한 관점에서 트위스트와 스크류는 본질적으로 같은 객체의 두 가지 표현이며, 둘 다 강체의 무한소 운동을 단일 6차원 객체로 통합한다. 대응되는 쌍대 객체로는 렌치(wrench)가 있으며, 렌치는 힘과 모멘트를 결합한 강체 정역학의 6차원 표현이다. 트위스트와 렌치의 쌍대성은 강체 운동학과 정역학의 통일적 형식화의 기반을 이룬다.

7. 로봇 공학에서의 활용

7.1 매니퓰레이터의 지수 곱 형식

매니퓰레이터의 운동학을 표현하는 현대적 접근법인 지수 곱 형식(product of exponentials formulation)은 스크류 운동의 행렬 표현을 직접적으로 활용한다. 각 관절의 운동을 단위 트위스트와 관절 변수로 표현하고, 전체 매니퓰레이터의 운동학은 다음과 같은 곱으로 주어진다.

T(\boldsymbol{\theta}) = e^{\theta_{1} \hat{\xi}_{1}^{\wedge}} e^{\theta_{2} \hat{\xi}_{2}^{\wedge}} \cdots e^{\theta_{n} \hat{\xi}_{n}^{\wedge}} \, M

여기서 M 은 모든 관절 변수가 영일 때의 말단 자세이고, \hat{\xi}_{i}i 번째 관절의 단위 트위스트이다. 이 표현은 좌표계 부착 없이 운동학을 기술할 수 있게 해주며, 운동학의 기하학적 직관을 명료하게 드러낸다.

7.2 자코비안의 트위스트 표현

매니퓰레이터의 자코비안 행렬은 각 관절 트위스트를 열로 가지는 형태로 자연스럽게 구성되며, 이를 공간 자코비안 또는 물체 자코비안이라 한다. 트위스트 기반의 자코비안 표현은 운동학적 해석을 단순화하고, 특이점 분석과 조작성 평가를 자연스럽게 가능케 한다.

7.3 무인기와 자율 주행 차량의 자세 갱신

강체의 자세 시간 적분은 본질적으로 트위스트 적분이며, 스크류 운동의 행렬 표현을 활용한 리 군 위의 적분기는 직교성과 강체 변환의 군 구조를 자동으로 보존한다. 이는 장기간의 자세 추정과 제어에서 누적 오차를 억제하는 핵심 도구이다.

7.4 강체 보간과 경로 계획

두 강체 자세 사이의 경로를 생성할 때, 시작 자세와 목표 자세 사이의 상대 변환을 단일 스크류 운동으로 분해하고, 그 매개변수 \theta 를 시간에 따라 보간하면 자연스러운 강체 보간이 얻어진다. 이는 위치와 자세를 분리하지 않고 통합적으로 보간하는 우아한 방법이며, 작업 경로 계획과 시각화에서 표준적으로 사용된다.

7.5 강체 동역학의 기하학적 형식화

강체의 동역학을 트위스트와 렌치로 기술하는 형식화는 라그랑주 방정식을 리 군 위에서 자연스럽게 표현하게 해주며, 결과적으로 매니퓰레이터 동역학의 재귀적 알고리듬과 변분 적분기 등의 고급 기법의 기반을 이룬다.

7.6 시각-관성 항법과 슬램

시각-관성 항법과 슬램에서 강체의 움직임을 트위스트로 표현하고, 측정 갱신을 강체 변환의 곱으로 기술하는 절차는 추정의 일관성과 효율성을 동시에 확보한다. 트위스트 좌표 위에서 정의된 공분산은 자세의 회전 부분과 병진 부분의 결합 불확실성을 자연스럽게 표현한다.

8. 결론

스크류 운동은 강체의 임의의 변위를 단일 직선에 대한 회전과 그 직선 방향의 병진의 결합으로 통합 표현하는 고전적이면서도 강력한 개념이다. 챠슬-모치 정리에 의해 모든 강체 변환은 본질적으로 스크류 운동이며, 트위스트와 지수 사상을 통하여 리 군 SE(3) 의 모든 원소가 단일 스크류 매개변수로 환원될 수 있음이 보장된다. 이러한 통합적 표현은 매니퓰레이터의 지수 곱 운동학, 자코비안의 트위스트 표현, 리 군 위의 자세 적분, 강체 보간, 동역학의 기하학적 형식화, 시각-관성 항법 등 로봇 공학의 광범위한 영역에서 핵심적인 역할을 수행한다. 스크류 운동의 수학적 본질에 대한 이해는 강체 운동을 단순한 회전과 병진의 결합이 아닌, 단일한 기하학적 운동으로 바라보는 통찰을 제공하며, 이는 현대 로봇 운동학의 가장 중요한 토대 중 하나이다.


참고문헌

  • Murray, R. M., Li, Z., & Sastry, S. S. (1994). A Mathematical Introduction to Robotic Manipulation. CRC Press.
  • Lynch, K. M., & Park, F. C. (2017). Modern Robotics: Mechanics, Planning, and Control. Cambridge University Press.
  • Selig, J. M. (2005). Geometric Fundamentals of Robotics (2nd ed.). Springer.
  • Ball, R. S. (1900). A Treatise on the Theory of Screws. Cambridge University Press.
  • Chasles, M. (1830). “Note sur les propriétés générales du système de deux corps semblables entr’eux.” Bulletin des Sciences Mathématiques, 14, 321–326.
  • Mozzi, G. (1763). Discorso matematico sopra il rotamento momentaneo dei corpi. Stamperia di Donato Campo.
  • Spong, M. W., Hutchinson, S., & Vidyasagar, M. (2020). Robot Modeling and Control (2nd ed.). Wiley.
  • Siciliano, B., Sciavicco, L., Villani, L., & Oriolo, G. (2009). Robotics: Modelling, Planning and Control. Springer.
  • Bottema, O., & Roth, B. (1979). Theoretical Kinematics. North-Holland.
  • Hunt, K. H. (1978). Kinematic Geometry of Mechanisms. Oxford University Press.

Version: 1.0