33.18 매니퓰러빌리티 지수에 의한 특이점 근접도 평가

매니퓰러빌리티 지수(manipulability index)는 자코비안의 특이값 곱으로 정의되는 스칼라 지표로, 매니퓰레이터의 현재 관절 구성에서 특이점까지의 상대적 거리를 연속 함수로 정량화한다. Yoshikawa가 1985년에 도입한 이 지표는 특이점 근접도 평가, 매니퓰러빌리티 최적화, 경로 계획 비용 함수 구성, 실시간 모니터링의 표준 수단으로 자리 잡았다. 본 절에서는 매니퓰러빌리티 지수의 정의, 기하학적·수치적 해석, 특이점 근접도 평가 기법, 경사도 기반 최적화, 확장 변형, 실무 응용을 해라체로 기술한다.

1. 매니퓰러빌리티 지수의 정의

1.1 Yoshikawa 지수

Yoshikawa가 1985년 논문 “Manipulability of robotic mechanisms“에서 제안한 매니퓰러빌리티 지수는 다음과 같이 정의된다.

$w(\vec{q}) = \sqrt{\det\bigl( \mathbf{J}(\vec{q}) \mathbf{J}^\top(\vec{q}) \bigr)}$

$m \le n$ 인 자코비안에 대하여 정의되며, 여유 자유도 매니퓰레이터의 표준 지표이다.

33.18.1.2 특이값 곱 표현

SVD를 통하면 매니퓰러빌리티 지수는 특이값의 곱과 같다.

$w(\vec{q}) = \prod_{i=1}^{m} \sigma_i(\mathbf{J}(\vec{q}))$

이는 특이값 분포 전체가 지표에 반영됨을 보여준다.

1.2 정사각 자코비안의 경우

$m = n$ 인 정사각 자코비안의 경우, 매니퓰러빌리티 지수는 행렬식의 절댓값과 같다.

$w(\vec{q}) = \lvert \det \mathbf{J}(\vec{q}) \rvert$

33.18.1.4 특이점에서의 영점

특이점에서 하나 이상의 특이값이 영이 되므로, 매니퓰러빌리티 지수는 영이 된다.

$\vec{q}_s \in \mathcal{S} \Longleftrightarrow w(\vec{q}_s) = 0$

2. 기하학적 해석

2.1 타원체 부피

매니퓰러빌리티 지수는 속도 매니퓰러빌리티 타원체의 부피에 비례한다.

$\text{부피} = \frac{\pi^{m/2}}{\Gamma(m/2 + 1)} \cdot w(\vec{q})$

여기서 $\Gamma$ 는 감마 함수이다. 특이점에서 타원체가 저차원으로 퇴화하여 부피가 영이 된다.

33.18.2.2 운동 능력 평균

매니퓰러빌리티 지수는 모든 방향의 운동 능력의 기하 평균에 해당한다. 한 방향이라도 능력이 크게 저하되면 지수가 급격히 감소한다.

33.18.2.3 스크류 이론적 해석

관절 트위스트들의 외적(exterior product) 노름과 관련된다. 트위스트들이 선형 종속이 되면 외적이 영이 되어 매니퓰러빌리티가 영이 된다.

33.18.2.4 등방성 조건과의 비교

매니퓰러빌리티가 최대가 되는 구성과 조건수가 1(등방성)인 구성은 일반적으로 일치하지 않는다. 두 지표는 서로 다른 최적화 기준을 제공한다.

33.18.3 특이점 근접도 평가

33.18.3.1 연속적 지표

매니퓰러빌리티 지수는 관절 변수의 연속 함수이며, 특이점에서 영에 이른다. 따라서 특이점 근접도의 연속적 정량화에 적합하다.

33.18.3.2 임계치 설정

실무에서는 매니퓰레이터의 일반 구성에서의 매니퓰러빌리티 $w_0$ 를 기준으로 상대적 임계치를 설정한다.

$w(\vec{q}) < \alpha \cdot w_0 \Longrightarrow \text{특이점 근접}$

일반적으로 $\alpha = 0.1$ ~ $0.3$ 의 값이 사용된다.

2.2 근접도의 선형성

특이점으로부터의 관절 공간 거리 $d = \lVert \vec{q} - \vec{q}_s \rVert$ 에 대하여, 단순 특이점 근방에서

$w(\vec{q}) \approx c \cdot d$

의 선형 관계가 성립한다(단, 다른 특이값들이 유계로 유지되는 경우).

33.18.3.4 고차 특이점에서의 거동

복합 특이점 근방에서는 매니퓰러빌리티 지수가 $d^{r_d}$ 에 비례하여 수렴한다. 결손 차수가 높을수록 접근이 빠르게 감지된다.

33.18.4 다른 특이 지표와의 비교

33.18.4.1 최소 특이값 지표

$\sigma_{\min}$ 은 가장 약한 방향의 능력을 직접 반영하며, 매니퓰러빌리티 지수보다 민감하게 특이점에 반응한다. 수치 안정성 관점에서는 $\sigma_{\min}$ 이 선호된다.

33.18.4.2 조건수 지표

조건수 $\kappa = \sigma_{\max}/\sigma_{\min}$ 은 상대적 이방성을 측정하며, 매니퓰러빌리티 지수와 다른 정보를 제공한다. 스케일에 무관한 절댓값 지표이다.

33.18.4.3 프로베니우스 조건 지수

Angeles와 Lopez-Cajun의 프로베니우스 조건 지수는 모든 특이값을 고려하는 평균 기반 지표로, 매니퓰러빌리티와 조건수의 중간적 성격을 가진다.

33.18.4.4 지표 선택 기준

특이점까지의 거리 평가: 매니퓰러빌리티 지수.
수치 안정성 평가: 최소 특이값.
이방성 평가: 조건수.
종합 성능 평가: 프로베니우스 조건 지수.

응용 맥락에 따라 적절한 지표를 선택하거나 병용한다.

33.18.5 경사도 기반 최적화

33.18.5.1 매니퓰러빌리티 경사도

매니퓰러빌리티 지수의 관절 변수에 대한 경사도는 여유 자유도 해소의 중요한 구성 요소이다.

$\nabla w(\vec{q}) = \frac{\partial w}{\partial \vec{q}}$

2.3 해석적 계산

자코비안 행렬식의 편미분을 통해 해석적으로 계산할 수 있다. $\det(\mathbf{J}\mathbf{J}^\top)$ 의 편미분은 행렬 미분 공식을 이용한다.

$\frac{\partial w^2}{\partial q_i} = w^2 \cdot \operatorname{tr}\!\left[ (\mathbf{J}\mathbf{J}^\top)^{-1} \frac{\partial (\mathbf{J}\mathbf{J}^\top)}{\partial q_i} \right]$

33.18.5.3 수치적 계산

유한 차분 근사

$\frac{\partial w}{\partial q_i} \approx \frac{w(\vec{q} + \Delta q \hat{\vec{e}}_i) - w(\vec{q} - \Delta q \hat{\vec{e}}_i)}{2 \Delta q}$

은 구현이 단순하지만 수치 노이즈에 민감하다.

2.4 자동 미분

최근의 로봇 동역학 라이브러리는 자동 미분을 통해 매니퓰러빌리티 경사도를 정확하게 계산한다. Pinocchio, Drake, CasADi 등이 지원한다.

3. 여유 자유도 해소에의 활용

3.1 영공간 자체 운동

Yoshikawa의 여유 자유도 해소에서 자체 운동 벡터는 매니퓰러빌리티 경사도 방향으로 설정된다.

$\dot{\vec{q}}_0 = k \nabla w(\vec{q}), \quad k > 0$

이 자체 운동은 엔드 이펙터 운동을 유지하면서 매니퓰러빌리티를 국소적으로 최대화한다.

33.18.6.2 Liégeois 기법의 특수 경우

Liégeois의 구성 제어 기법에서 보조 목적 함수로 $-w(\vec{q})$ 를 사용하면, 매니퓰러빌리티 최대화 자체 운동이 자연스럽게 유도된다.

33.18.6.3 계층적 제어

작업 우선순위 해소에서 매니퓰러빌리티 최대화를 하위 작업으로 설정하면, 엔드 이펙터 추종이 우선되면서 특이점 회피가 보조적으로 수행된다.

33.18.6.4 수렴 특성

경사도 기반 최적화는 국소 최적해로 수렴한다. 전역 최적해 보장을 위해서는 추가적 탐색 기법이 필요하다.

33.18.7 확장 변형

33.18.7.1 위치 매니퓰러빌리티

위치 자코비안 $\mathbf{J}_v$ 에 기반한 위치 매니퓰러빌리티는 위치 방향 능력만을 평가한다.

$w_v(\vec{q}) = \sqrt{\det(\mathbf{J}_v \mathbf{J}_v^\top)}$

3.2 자세 매니퓰러빌리티

자세 자코비안 $\mathbf{J}_\omega$ 에 기반한 자세 매니퓰러빌리티는 자세 방향 능력을 평가한다.

3.3 가중 매니퓰러빌리티

작업 공간 가중 행렬 $\mathbf{W}$ 를 도입한 가중 매니퓰러빌리티는 방향별 중요도를 반영한다.

$w_W(\vec{q}) = \sqrt{\det(\mathbf{J} \mathbf{W} \mathbf{J}^\top)}$

33.18.7.4 동역학 일관 매니퓰러빌리티

관성 가중 매니퓰러빌리티 $w_M = \sqrt{\det(\mathbf{J} \mathbf{M}^{-1} \mathbf{J}^\top)}$ 는 동역학 관점에서의 능력을 평가한다. Khatib의 연산 공간 정식화와 연계된다.

33.18.7.5 방향별 매니퓰러빌리티

Chiu의 방향별 지표 $u(\hat{\vec{v}}) = 1 / \sqrt{\hat{\vec{v}}^\top (\mathbf{J}\mathbf{J}^\top)^{-1} \hat{\vec{v}}}$ 는 특정 방향 $\hat{\vec{v}}$ 에서의 능력을 평가한다. 작업 방향 정렬 분석에 유용하다.

33.18.8 실시간 모니터링

33.18.8.1 계산 비용

매니퓰러빌리티 지수의 계산은 $\mathbf{J} \mathbf{J}^\top$ 의 행렬식 평가를 요구하며, $\mathrm{O}(m^3)$ 의 복잡도이다. Cholesky 분해의 대각 원소 곱으로 효율적으로 계산 가능하다.

33.18.8.2 제어 주기 내 평가

6자유도 또는 7자유도 매니퓰레이터에서 매니퓰러빌리티 계산은 마이크로초 단위로 수행 가능하다. 1 kHz 이상의 제어 주기에서도 실시간 평가가 가능하다.

33.18.8.3 경사도 계산의 비용

경사도 계산은 각 관절 변수에 대한 편미분을 요구하므로 $\mathrm{O}(n \cdot m^3)$ 의 복잡도이다. 자동 미분 기법으로 효율 개선이 가능하다.

33.18.8.4 캐싱 활용

제어 주기 내에서 자코비안, 매니퓰러빌리티, 경사도 계산 결과를 공유하면 연산 비용이 절감된다.

33.18.9 경로 계획에의 활용

33.18.9.1 비용 함수 구성

경로 계획의 비용 함수에 매니퓰러빌리티 기반 페널티를 포함한다.

$C(\vec{q}) = C_{\mathrm{path}}(\vec{q}) + \beta / w(\vec{q})$

$w$ 가 작으면 비용이 커지므로 특이점 회피 경로가 선호된다.

3.4 RRT, PRM 기반 계획

샘플링 기반 경로 계획기에서 노드 간 연결 판정에 매니퓰러빌리티 임계치를 포함한다. 임계치 미만 구성은 연결에서 제외된다.

3.5 궤적 최적화

최적화 기반 궤적 계획은 매니퓰러빌리티 지수를 직접 목적 함수로 사용하거나 부등식 제약 $w(\vec{q}) \ge w_{\min}$ 으로 포함한다.

3.6 전역 매니퓰러빌리티

작업 공간 전반의 매니퓰러빌리티 분포는 작업 공간 해석의 지표이다. Gosselin과 Angeles의 전역 조건 지수가 관련된다.

4. 실무 응용

4.1 산업용 매니퓰레이터

산업용 컨트롤러는 매니퓰러빌리티 지수를 내장하여 특이점 근방 감지 및 속도 감속을 수행한다. ABB, KUKA, FANUC 등의 제품이 이를 표준 기능으로 제공한다.

4.2 협동 로봇

협동 로봇은 매니퓰러빌리티 기반 안전 속도 제한을 구현한다. 특이점 근방에서 자동 감속하여 인간 안전을 확보한다.

4.3 여유 자유도 매니퓰레이터

7자유도 협동 로봇은 매니퓰러빌리티 최대화 자체 운동을 지속적으로 수행하여 특이점 진입을 근본적으로 회피한다.

4.4 시뮬레이션과 교육

로봇 시뮬레이션 소프트웨어는 매니퓰러빌리티 타원체와 지수를 시각화한다. 이는 특이점 개념의 교육적 이해를 돕는 강력한 도구이다.

5. Salisbury–Craig의 공헌

5.1 운동학적 등방성 개념

Salisbury와 Craig가 1982년 논문 “Articulated hands: Force control and kinematic issues“에서 제시한 운동학적 등방성 개념은 매니퓰러빌리티 분석의 초기 기반을 제공하였다.

5.2 다지 그리퍼 해석

다지 로봇 손의 매니퓰러빌리티 분석은 이 논문의 주요 기여이며, 이후 매니퓰러빌리티 지수의 일반화 연구의 출발점이 되었다.

5.3 등방 구성의 설계

등방 구성(isotropic configuration)은 조건수가 1인 특수 구성으로, 매니퓰러빌리티 극대화와 관련된다. 설계 최적화의 이상적 기준으로 활용된다.

6. 본 절의 학술적 정리

본 절에서 다룬 매니퓰러빌리티 지수는 특이값의 곱으로 정의되는 스칼라 지표로, 특이점 근접도를 연속적으로 정량화하는 표준 수단이다. 속도 매니퓰러빌리티 타원체의 부피에 비례하며, 특이점에서 영으로 수렴한다. 최소 특이값, 조건수, 프로베니우스 조건 지수와의 비교에서 고유한 관점(부피 기반)을 제공하며, 실무에서는 상황에 따라 병용된다. 경사도 기반 최적화는 여유 자유도 매니퓰레이터의 자체 운동 설계와 경로 계획 비용 함수에 활용된다. 위치·자세 매니퓰러빌리티, 가중 매니퓰러빌리티, 동역학 일관 매니퓰러빌리티, 방향별 매니퓰러빌리티 등의 확장 변형이 응용 맥락에 따라 사용된다. 산업용 매니퓰레이터, 협동 로봇, 여유 자유도 매니퓰레이터, 시뮬레이션 도구에서 매니퓰러빌리티 지수는 특이점 근접도 평가의 표준 구성 요소로 기능한다.

7. 출처

Salisbury, J. K. and Craig, J. J., “Articulated hands: Force control and kinematic issues”, International Journal of Robotics Research, Vol. 1, No. 1, pp. 4–17, 1982.
Yoshikawa, T., “Manipulability of robotic mechanisms”, International Journal of Robotics Research, Vol. 4, No. 2, pp. 3–9, 1985.
Chiu, S. L., “Task compatibility of manipulator postures”, International Journal of Robotics Research, Vol. 7, No. 5, pp. 13–21, 1988.
Yoshikawa, T., Foundations of Robotics: Analysis and Control, MIT Press, 1990.
Gosselin, C. and Angeles, J., “A global performance index for the kinematic optimization of robotic manipulators”, Journal of Mechanical Design, Vol. 113, No. 3, pp. 220–226, 1991.
Angeles, J. and Lopez-Cajun, C. S., “Kinematic isotropy and the conditioning index of serial robotic manipulators”, International Journal of Robotics Research, Vol. 11, No. 6, pp. 560–571, 1992.
Siciliano, B., Sciavicco, L., Villani, L., and Oriolo, G., Robotics: Modelling, Planning and Control, Springer, 2009.
Spong, M. W., Hutchinson, S., and Vidyasagar, M., Robot Modeling and Control, 2nd edition, Wiley, 2020.

8. 버전

문서 버전: 2.0
작성일: 2026-04-21