32.28 매니퓰러빌리티 지수를 이용한 성능 평가

매니퓰러빌리티 지수는 매니퓰레이터의 국소 운동학적 능력을 단일 스칼라로 정량화하는 표준 지표이나, 이를 실제 성능 평가에 활용하기 위해서는 평가 영역의 선정, 정규화, 통계적 집계, 응용 맥락에 따른 가중치 선택 등의 추가적 고려가 요구된다. 본 절에서는 매니퓰러빌리티 지수를 이용한 국소 성능 평가, 전역 성능 평가, 작업 의존적 성능 평가, 설계 최적화에의 활용, 그리고 실험적·수치적 평가 절차를 학술적으로 체계화한다.

1. 국소 성능 평가

국소 성능 평가(local performance evaluation)는 매니퓰레이터의 특정 구성 \vec{q}에서의 운동학적 능력을 매니퓰러빌리티 지수 w(\vec{q}) = \sqrt{\det(\mathbf{J}(\vec{q}) \mathbf{J}^\top(\vec{q}))}로 정량화하는 절차를 의미한다.

국소 평가의 표준 절차는 다음과 같이 정리된다.

첫째, 자코비안의 계산: 주어진 구성 \vec{q}에서 기하학적 자코비안 \mathbf{J}(\vec{q})를 계산한다. 6자유도 이상의 공간 매니퓰레이터의 경우 \mathbf{J} \in \mathbb{R}^{6 \times n}이며, 평면 매니퓰레이터의 경우 \mathbf{J} \in \mathbb{R}^{3 \times n}이다.

둘째, 자코비안 곱 행렬의 구성: \mathbf{J} \mathbf{J}^\top \in \mathbb{R}^{m \times m}을 계산한다. 이 행렬은 완전 행 계수의 정규 구성에서 대칭 양정치이다.

셋째, 매니퓰러빌리티 지수의 계산: w(\vec{q}) = \sqrt{\det(\mathbf{J} \mathbf{J}^\top)}을 수치적으로 계산하거나, 특이값 분해를 사용하여 w = \prod_{i=1}^{m} \sigma_i로 계산한다.

넷째, 구성 간 비교: 서로 다른 구성 \vec{q}_1, \vec{q}_2에 대해 w(\vec{q}_1)과 w(\vec{q}_2)를 비교하여 국소 운동학적 능력의 우열을 판정한다.

국소 평가는 작업 실행 중 매 순간의 성능 모니터링, 실시간 특이점 근접도 측정, 자세 선택의 기준으로 활용된다.

2. 전역 성능 평가

전역 성능 평가(global performance evaluation)는 작업 공간 전체 또는 관심 영역에 걸친 매니퓰러빌리티의 통계적 집계를 의미하며, 매니퓰레이터 구조의 본질적 성능을 평가한다.

2.1 전역 매니퓰러빌리티 지수

Gosselin과 Angeles(1991)에 의해 제안된 전역 성능 지수(Global Performance Index, GPI)는 작업 공간 \mathcal{W}에 걸친 매니퓰러빌리티의 적분 평균으로 정의된다.

\eta_w = \frac{\int_{\mathcal{W}} w(\vec{q}) \, dV}{\int_{\mathcal{W}} dV}

여기서 dV는 작업 공간의 미분 요소이다. 이 지표는 매니퓰레이터 구조가 작업 공간 전역에서 평균적으로 얼마나 우수한 운동학적 능력을 제공하는지를 단일 수치로 표현한다.

32.28.2.2 최소 매니퓰러빌리티

전역 평균과 상호 보완적인 지표로서 작업 공간 내 최소 매니퓰러빌리티 w_{\min} = \inf_{\vec{q} \in \mathcal{Q}(\mathcal{W})} w(\vec{q})가 사용된다. 이 지표는 작업 공간 내에서 가장 성능이 저하되는 지점에서의 매니퓰러빌리티를 나타내며, 특이점 근접 영역의 존재 여부를 판정하는 기준이 된다.

32.28.2.3 매니퓰러빌리티 편차

Kim과 Khosla(1991)가 제안한 매니퓰러빌리티 편차는 작업 공간 전역에서 매니퓰러빌리티의 변동을 측정한다.

\sigma_w^2 = \frac{\int_{\mathcal{W}} (w(\vec{q}) - \eta_w)^2 \, dV}{\int_{\mathcal{W}} dV}

편차가 작으면 작업 공간 전역에서 매니퓰러빌리티가 균일하게 유지됨을 의미하며, 일관된 성능을 요구하는 응용에서 바람직한 성질이다.

2.2 수치적 적분 절차

전역 지표의 계산은 일반적으로 해석적 적분이 불가능하므로 수치적 적분을 요구한다. 표준 절차는 다음과 같다.

격자 샘플링: 작업 공간 \mathcal{W}를 균일 격자로 이산화하여 N개의 점 \{\vec{x}_k\}_{k=1}^{N}을 생성한다.

역기구학 해법: 각 점 \vec{x}_k에 대해 역기구학으로 구성 \vec{q}_k를 계산한다. 다중 해가 존재하는 경우 특정 분지 또는 매니퓰러빌리티가 최대인 구성을 선택한다.

매니퓰러빌리티 평가: 각 \vec{q}_k에서 w(\vec{q}_k)를 계산한다.

통계적 집계: 평균, 최소, 편차 등의 통계량을 계산한다.

\eta_w \approx \frac{1}{N} \sum_{k=1}^{N} w(\vec{q}_k), \quad w_{\min} \approx \min_{k} w(\vec{q}_k)

격자의 조밀도는 작업 공간의 기하학적 복잡도와 매니퓰러빌리티의 공간적 변동에 따라 결정되며, 수렴성 검증을 위해 격자를 점진적으로 세분화하는 절차가 병행된다.

32.28.3 작업 의존적 성능 평가

매니퓰러빌리티 지수는 모든 운동 방향에 대한 능력을 대칭적으로 취급하지만, 실제 작업은 특정 방향에서의 운동 능력을 우선적으로 요구하는 경우가 많다. 이러한 작업 의존적 성능 평가(task-dependent performance evaluation)를 위한 확장 지표가 제안되어 왔다.

32.28.3.1 작업 호환성 지표

Chiu(1988)의 작업 호환성 지표(task compatibility index)는 작업 요구 방향 \vec{u}에 대한 매니퓰레이터의 운동 능력을 정량화한다. 속도 호환성 지표는 다음과 같이 정의된다.

p_v(\vec{q}, \vec{u}) = \frac{1}{\sqrt{\vec{u}^\top (\mathbf{J} \mathbf{J}^\top)^{-1} \vec{u}}}

이 지표는 매니퓰러빌리티 타원체의 \vec{u} 방향 반지름과 일치하며, 해당 방향으로 생성 가능한 최대 선속도를 의미한다. 작업 요구 방향이 타원체의 장축과 정렬되면 p_v가 크고, 단축과 정렬되면 p_v가 작다.

2.3 가중 매니퓰러빌리티

Park과 Brockett(1994)이 제안한 가중 매니퓰러빌리티는 작업 공간의 방향별 중요도를 반영한 행렬 \mathbf{W}를 도입한다.

w_{\text{weighted}}(\vec{q}) = \sqrt{\det(\mathbf{J} \mathbf{W} \mathbf{J}^\top)}

여기서 \mathbf{W}는 관절 공간의 가중 행렬 또는 작업 공간의 변환 행렬로 해석된다. 이를 통해 작업 방향 또는 관절별 선호도를 반영한 성능 평가가 가능하다.

32.28.3.3 조작 능력의 방향 해석

작업이 특정 방향의 힘 또는 속도를 요구할 때, 매니퓰러빌리티 타원체의 형상(이방성)이 성능에 결정적 영향을 미친다. 이 경우 매니퓰러빌리티 지수 단독보다는 조건수 및 타원체의 주축 방향이 함께 고려된다. 타원체의 장축 방향이 작업 요구 방향과 정렬되고 단축 방향이 구속 방향과 정렬된 구성이 이상적인 작업 자세이다.

32.28.4 성능 평가의 정규화 문제

매니퓰러빌리티 지수의 비교 해석에는 정규화(normalization)가 요구된다. 매니퓰러빌리티 지수는 링크 길이, 관절 변수의 스케일, 작업 공간의 크기에 민감하게 의존하므로, 서로 다른 매니퓰레이터 간의 직접적인 수치 비교는 물리적 의미를 가지지 않는다.

32.28.4.1 링크 길이 정규화

Angeles(2002)의 제안에 따라 매니퓰레이터의 특성 길이 L(예: 총 링크 길이의 합 또는 최장 링크 길이)을 사용하여 매니퓰러빌리티 지수를 무차원화한다.

\bar{w}(\vec{q}) = \frac{w(\vec{q})}{L^m}

여기서 m은 작업 공간의 차원이다. 무차원화된 \bar{w}는 서로 다른 크기의 매니퓰레이터 간 비교에 사용될 수 있다.

2.4 최대 매니퓰러빌리티 기준 정규화

특정 매니퓰레이터의 작업 공간 전역에서의 최대 매니퓰러빌리티 w_{\max}로 나누어 상대적 매니퓰러빌리티를 정의한다.

\hat{w}(\vec{q}) = \frac{w(\vec{q})}{w_{\max}} \in [0, 1]

이 정규화는 개별 매니퓰레이터 내에서 구성별 상대 성능을 직관적으로 표현하며, 성능 저하 영역의 식별에 유용하다.

32.28.4.3 단위 일관성 정규화

선속도와 각속도가 혼재된 기하학적 자코비안에서의 매니퓰러빌리티는 단위가 일관되지 않는 문제가 있다. 이를 해결하기 위해 특성 길이 L_c를 도입하여 각속도 블록을 선속도 차원으로 변환한다.

\mathbf{J}_{\text{norm}} = \begin{bmatrix} \mathbf{J}_v \\ L_c \, \mathbf{J}_\omega \end{bmatrix}

이 정규화된 자코비안으로부터 계산된 매니퓰러빌리티 지수는 물리적 일관성을 갖는다.

3. 설계 최적화에의 활용

매니퓰러빌리티 지수는 매니퓰레이터의 설계 단계에서 구조 매개변수의 결정 기준으로 활용된다.

3.1 링크 길이 최적화

2자유도 평면 매니퓰레이터의 경우 매니퓰러빌리티 지수는 w = a_1 a_2 |\sin \theta_2|로 주어진다. 작업 공간 전체에 걸친 평균 매니퓰러빌리티를 최대화하는 링크 길이 비율은 해석적으로 유도될 수 있으며, 고전적 결과로서 a_1 = a_2에서 작업 공간의 평균 매니퓰러빌리티가 최적화된다는 결과가 알려져 있다.

일반적 직렬 매니퓰레이터에서 링크 길이 벡터 \vec{a} = (a_1, a_2, \ldots, a_n)에 대해 최적화 문제는 다음과 같이 설정된다.

\max_{\vec{a}} \; \eta_w(\vec{a}) \quad \text{subject to} \quad \sum_i a_i = L_{\text{total}}, \quad a_i > 0

이 문제는 일반적으로 수치적 최적화를 요구하며, 경사법 또는 진화 알고리즘이 적용된다.

32.28.5.2 관절 배치 최적화

DH 매개변수 중 링크 오프셋, 비틀림 각, 관절 축의 위치 등의 배치 매개변수도 매니퓰러빌리티 기반으로 최적화된다. 각 설계 변수에 대한 매니퓰러빌리티의 민감도 분석이 선행되며, 이를 바탕으로 설계 공간의 축소 또는 국소 탐색이 수행된다.

32.28.5.3 작업 배치 설계

매니퓰레이터의 구조가 고정된 경우, 작업 대상 또는 공작물의 위치 \vec{x}_{\text{work}}를 설계 변수로 하여 매니퓰러빌리티를 최적화한다.

\max_{\vec{x}_{\text{work}}} \; \int_{\mathcal{T}(\vec{x}_{\text{work}})} w(\vec{q}(\vec{x})) \, d\vec{x}

여기서 \mathcal{T}는 작업 대상의 위치에 의존하는 작업 경로이다. 이 최적화는 조립 셀 설계, 용접 공정 레이아웃, 의료 로봇 수술 환경 배치 등에 활용된다.

3.2 다목적 최적화

실제 설계에서는 매니퓰러빌리티 단독이 아닌, 작업 공간 크기, 강성, 속도 한계, 관성, 질량 등 다수의 설계 기준이 동시에 고려된다. 이 경우 매니퓰러빌리티는 다목적 최적화의 한 구성 요소로 포함되며, 파레토 최적해 또는 가중합 방법이 적용된다.

4. 작업 자세 선택에의 활용

여유 자유도 매니퓰레이터 또는 다중 역기구학 해가 존재하는 매니퓰레이터에서 매니퓰러빌리티 지수는 자세 선택의 기준으로 활용된다.

4.1 다중 역기구학 해의 선택

6자유도 매니퓰레이터의 일반 구성에서 동일 엔드 이펙터 자세에 대해 최대 8개의 역기구학 해가 존재할 수 있다. 이 중 매니퓰러빌리티가 최대인 해를 선택하면 특이점 근접도가 낮고 경로 상의 성능 저하가 최소화된다.

4.2 여유 자유도의 팔 각도 선택

7자유도 매니퓰레이터의 경우 팔 각도(arm angle)라는 1차원 여유 자유도가 존재한다. 팔 각도 \phi의 선택 기준으로 매니퓰러빌리티가 사용된다.

\phi^* = \arg\max_{\phi} \; w(\vec{q}(\vec{x}, \phi))

이를 통해 동일 엔드 이펙터 자세에서도 매니퓰러빌리티가 최대인 구성을 선택하여 운동 능력을 최대화한다.

32.28.6.3 실시간 자세 전환

작업 중 특이점에 접근하는 경우 매니퓰러빌리티 지수를 모니터링하여 임계값 이하로 감소하면 다른 역기구학 분지로 자세를 전환한다. 이 절차는 연속적 경로 추종을 보장하면서 성능 저하를 회피한다.

32.28.7 궤적 계획에의 활용

엔드 이펙터의 경로가 결정된 상태에서 관절 공간의 궤적을 매니퓰러빌리티 기반으로 최적화하는 문제가 제기된다.

32.28.7.1 경로 적분형 성능 지표

경로 \mathcal{P} 상의 평균 매니퓰러빌리티는 다음과 같이 정의된다.

\bar{w}_{\mathcal{P}} = \frac{1}{T} \int_0^T w(\vec{q}(t)) \, dt

여유 자유도가 있는 경우 이 지표를 최대화하는 관절 공간 궤적을 선택한다.

4.3 경로 최소 매니퓰러빌리티

경로 상의 최소 매니퓰러빌리티는 경로 전체의 최악 성능을 나타낸다.

w_{\min, \mathcal{P}} = \min_{t \in [0, T]} w(\vec{q}(t))

이 지표를 최대화하는 궤적은 특이점 회피를 보장한다.

32.28.7.3 제약 조건이 있는 최적 제어

매니퓰러빌리티 기반 궤적 최적화는 제약 조건이 있는 최적 제어 문제로 공식화된다.

\max_{\vec{q}(t)} \; \int_0^T w(\vec{q}(t)) \, dt \quad \text{subject to} \quad \mathbf{f}(\vec{q}(t)) = \vec{x}(t)

여기서 \mathbf{f}는 순기구학 사상이고 \vec{x}(t)는 규정된 엔드 이펙터 경로이다. 이 문제는 폰트랴긴 최대 원리 또는 직접 수치법으로 해결된다.

5. 실시간 성능 모니터링

매니퓰레이터의 실제 운용 중 매니퓰러빌리티 지수를 실시간 모니터링하여 성능 저하를 조기에 감지하고 대응하는 절차가 표준 관행으로 확립되어 있다.

5.1 특이점 경고 시스템

매니퓰러빌리티 지수가 미리 설정된 임계값 w_{\text{threshold}} 이하로 감소하면 경고 신호를 생성한다. 임계값은 일반적으로 해당 매니퓰레이터의 작업 공간 내 중간값 매니퓰러빌리티의 일정 비율로 설정된다.

5.2 적응적 속도 제한

매니퓰러빌리티 지수에 따라 엔드 이펙터 속도를 적응적으로 제한하는 방법이 사용된다. 매니퓰러빌리티가 감소하면 동일 엔드 이펙터 속도에 대응하는 관절 속도가 증가하므로 관절 속도 한계를 초과할 위험이 커진다. 다음과 같은 속도 스케일링이 사용된다.

\vec{v}_{\text{scaled}}(t) = \min\!\left(1, \frac{w(\vec{q}(t))}{w_{\text{threshold}}}\right) \cdot \vec{v}_{\text{desired}}(t)

32.28.8.3 경사 기반 회피 제어

매니퓰러빌리티의 경사 \nabla w를 이용한 영공간 기반 회피 제어는 특이점 근방에서 \nabla w 방향의 자기 운동을 실행하여 매니퓰러빌리티를 증가시킨다. 제어 법칙은 다음과 같다.

\dot{\vec{q}} = \mathbf{J}^+ \dot{\vec{x}} + (\mathbf{I} - \mathbf{J}^+ \mathbf{J}) \, k \, \nabla w

여기서 \mathbf{J}^+는 의사 역행렬이고 k > 0은 이득 계수이다. 두 번째 항은 영공간 투영이므로 엔드 이펙터 운동에 영향을 주지 않으면서 매니퓰러빌리티를 최대화한다.

6. 매니퓰레이터 간 비교 평가

서로 다른 매니퓰레이터 구조의 성능 비교에서 매니퓰러빌리티 기반 평가는 다음의 절차를 따른다.

첫째, 비교 대상 매니퓰레이터를 동일한 작업 공간 기준(예: 도달 가능 부피) 또는 동일한 총 링크 길이로 정규화한다.

둘째, 동일한 샘플링 격자 또는 통계적 샘플링 방법으로 작업 공간을 이산화한다.

셋째, 각 매니퓰레이터에 대해 매니퓰러빌리티의 평균, 최소, 표준 편차 등의 통계량을 계산한다.

넷째, 무차원화된 매니퓰러빌리티 \bar{w} 또는 상대적 매니퓰러빌리티 \hat{w}를 비교 척도로 사용한다.

다섯째, 성능 지표만으로 부족한 경우 조건수, 작업 공간 부피, 이방성 등의 보조 지표를 함께 보고한다.

이러한 체계적 비교는 새로운 매니퓰레이터 설계의 평가 또는 기존 설계의 선정에 활용된다.

7. 수치적 고려 사항

매니퓰러빌리티 기반 성능 평가의 실무적 구현에서는 다음의 수치적 고려 사항이 요구된다.

행렬식 계산의 안정성: \det(\mathbf{J} \mathbf{J}^\top)의 직접 계산은 차원이 큰 경우 수치적으로 불안정할 수 있다. 특이값 분해를 통한 w = \prod \sigma_i 계산이 보다 안정적이다.

특이점 근방의 수치 문제: 매니퓰러빌리티가 영에 접근하는 경우 부동 소수점 표현의 한계로 정확한 값의 계산이 어렵다. 로그 변환 \log w = \sum_i \log \sigma_i을 사용하거나 감쇠 매니퓰러빌리티를 정의한다.

역기구학의 다중 해 처리: 전역 평가 시 각 작업 공간 점에 대해 어느 역기구학 분지를 선택할지 명시적으로 정의해야 한다. 이는 최대 매니퓰러빌리티를 갖는 분지 또는 특정 분지(팔꿈치 위 등)를 선택하는 규칙으로 처리된다.

계산 비용: 대규모 작업 공간 격자의 모든 점에서 매니퓰러빌리티를 계산하는 것은 계산 비용이 크다. 적응적 샘플링, 병렬 계산, 대칭성 활용 등의 기법으로 효율성이 개선된다.

8. 한계와 주의 사항

매니퓰러빌리티 기반 성능 평가의 유용성에도 불구하고 다음의 한계와 주의 사항이 존재한다.

이방성 무시: 매니퓰러빌리티 지수는 타원체 부피만을 측정하므로, 크게 이방적인 타원체를 가진 구성이 높은 매니퓰러빌리티로 평가될 수 있다. 이방성의 반영은 조건수 또는 최소 특이값의 병용으로 해결된다.

작업 비의존성: 매니퓰러빌리티 지수 자체는 모든 방향의 운동 능력을 대칭적으로 취급하므로, 특정 방향의 운동을 우선시하는 작업에는 직접 적용이 부적절하다. 작업 호환성 지표 또는 가중 매니퓰러빌리티가 대안으로 사용된다.

동역학적 효과 미반영: 매니퓰러빌리티 지수는 기구학 기반 지표이므로 관성, 중력, 토크 한계 등의 동역학적 제약을 반영하지 않는다. 동역학적 성능을 요구하는 응용에서는 동역학적 매니퓰러빌리티가 함께 사용된다.

전역 평가의 분지 의존성: 전역 매니퓰러빌리티 평가는 선택된 역기구학 분지에 의존한다. 실제 매니퓰레이터는 작업 중 분지 전환이 제한되므로, 분지별 전역 평가가 보다 현실적이다.

정규화의 모호성: 링크 길이 정규화 또는 특성 길이 선택에 대한 보편적 기준이 존재하지 않으며, 비교 결과가 정규화 방법에 의존할 수 있다.

단위 혼재 문제: 선속도와 각속도의 단위 혼재는 기하학적 자코비안 기반 매니퓰러빌리티의 해석에 지속적인 어려움을 제기한다.

9. 학술적 의의

본 절에서 정리한 매니퓰러빌리티 지수 기반 성능 평가 방법론은 로봇 공학에서 다음의 학술적 의의를 가진다. 첫째, 국소 평가는 실시간 성능 모니터링과 자세 선택의 이론적 기초를 제공한다. 둘째, 전역 평가는 매니퓰레이터 구조의 본질적 성능 비교와 설계 최적화의 기준을 제공한다. 셋째, 작업 의존적 평가는 매니퓰러빌리티 지수를 실제 작업 맥락에 적합하도록 확장하는 방법론을 제공한다. 넷째, 정규화 기법은 서로 다른 매니퓰레이터 간 비교의 수학적 기반을 확립한다. 다섯째, 설계 최적화, 자세 선택, 궤적 계획, 실시간 모니터링의 응용은 매니퓰러빌리티 지수의 실무적 유용성을 입증한다. 여섯째, 매니퓰러빌리티 기반 평가의 한계에 대한 인식은 조건수, 최소 특이값, 동역학적 매니퓰러빌리티 등 보완적 지표의 필요성을 학술적으로 정당화한다. 일곱째, 본 절은 후속 절에서 다룰 자코비안의 조건수 및 등방성 조건이 왜 필요한지에 대한 동기를 제공한다.

10. 출처

• Yoshikawa, T., “Manipulability of robotic mechanisms”, International Journal of Robotics Research, Vol. 4, No. 2, pp. 3–9, 1985.
• 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.
• Kim, J. O. and Khosla, P. K., “Dexterity measures for design and control of manipulators”, Proceedings of the IEEE International Conference on Robotics and Automation, pp. 758–763, 1991.
• Chiu, S. L., “Task compatibility of manipulator postures”, International Journal of Robotics Research, Vol. 7, No. 5, pp. 13–21, 1988.
• Park, F. C. and Brockett, R. W., “Kinematic dexterity of robotic mechanisms”, International Journal of Robotics Research, Vol. 13, No. 1, pp. 1–15, 1994.
• Angeles, J., Fundamentals of Robotic Mechanical Systems: Theory, Methods, and Algorithms, 2nd edition, Springer, 2002.
• Angeles, J. and López-Cajún, 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.
• Spong, M. W., Hutchinson, S., and Vidyasagar, M., Robot Modeling and Control, 2nd edition, Wiley, 2020.
• Siciliano, B., Sciavicco, L., Villani, L., and Oriolo, G., Robotics: Modelling, Planning and Control, Springer, 2009.
• Lynch, K. M. and Park, F. C., Modern Robotics: Mechanics, Planning, and Control, Cambridge University Press, 2017.

11. 버전

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