18.18 다점 접촉의 역학 해석

1. 개요

다점 접촉은 하나의 접촉체가 여러 개의 접촉점이나 접촉 영역을 통해 다른 접촉체와 상호작용하는 상황을 의미한다. 로봇의 여러 파지 핑거가 물체와 동시에 접촉하는 경우, 이동 로봇의 여러 바퀴가 지면과 접촉하는 경우, 보행 로봇의 여러 발이 접촉하는 경우가 대표적 예이다. 다점 접촉의 역학 해석은 각 접촉점에서의 힘과 모멘트를 전체 강체 동역학과 일관되게 결정해야 하며, 이는 일반적으로 여러 개의 제약 조건과 불확정성이 결합된 복잡한 문제가 된다. 본 절은 다점 접촉의 수학적 정식화, 정역학적 조건, 불확정성과 과결정성, 상보성 정식화, 수치 해법, 로봇 공학적 응용을 학술적으로 정리한다.

2. 수학적 정식화

다점 접촉에서 n개의 접촉점이 존재하는 경우, 각 접촉점 i에서의 법선력 f_{n,i}와 접선력 \mathbf{f}_{t,i}는 체의 강체 운동에 기여한다. 각 접촉점에서의 힘과 위치 \mathbf{r}_i가 결합되어 전체 반력 \mathbf{F}_c와 모멘트 \mathbf{M}_c가 다음과 같이 구성된다.

\mathbf{F}_c = \sum_{i=1}^n \bigl(f_{n,i}\,\mathbf{n}_i + \mathbf{f}_{t,i}\bigr)

\mathbf{M}_c = \sum_{i=1}^n \mathbf{r}_i \times \bigl(f_{n,i}\,\mathbf{n}_i + \mathbf{f}_{t,i}\bigr)

여기서 \mathbf{n}_i는 접촉점 i에서의 법선 벡터이다. 체의 병진과 회전 방정식은 이러한 반력과 외력의 합에 의해 결정되며, 여기서 각 접촉점의 법선력은 비음수이고 접선력은 마찰 원추 내에 있어야 한다는 제약이 추가된다.

3. 정역학적 조건과 해의 존재

정적 평형 상태의 다체가 여러 점에서 접촉하는 경우, 외력과 중력을 반력으로 균형시키는 접촉력의 집합이 존재하여야 한다. 이러한 접촉력 집합은 법선력의 비음수 조건, 접선력의 마찰 원추 조건, 전체 힘과 모멘트의 균형 조건을 만족해야 한다. 이 조건 집합은 선형 계획 문제(LP)나 2차 계획 문제(QP)로 정식화될 수 있으며, 해의 존재는 접촉 구성과 외력의 기하학적 관계에 의해 결정된다. 해가 존재하는 경우에도 일반적으로 여러 개의 가능한 접촉력 분포가 존재하므로, 해의 불확정성을 다루는 절차가 필요하다.

4. 과결정성과 불확정성

다점 접촉에서 접촉점의 수가 체의 자유도에 비해 많은 경우, 접촉력 분포가 유일하게 결정되지 않는 정적 과결정성이 발생한다. 예를 들어 네 개의 다리로 지탱된 강체의 경우, 세 점만으로도 평면에서의 균형이 유지되므로, 네 번째 접촉점의 반력은 고유하게 결정되지 않는다. 이러한 불확정성은 재료의 탄성 변형을 고려한 탄성 접촉 해석을 통해 해결될 수 있으며, 실제 접촉 강성의 분포가 힘의 분배를 결정한다. 강체 가정에 기반한 해석에서는 추가적인 기준(예: 반력의 합의 최소화, 특정 함수의 최적화)이 도입되어야 한다.

5. 상보성 정식화

비평활 역학의 관점에서 다점 접촉은 상보성 조건으로 엄밀히 정식화된다. 각 접촉점 i에서 법선 방향의 상대 속도 v_{n,i}와 법선력 f_{n,i}는 다음의 상보 관계를 만족한다.

0 \le f_{n,i} \perp v_{n,i} \ge 0

이는 접촉이 유지되는 경우(v_{n,i} = 0) 법선력이 양일 수 있고, 분리되는 경우(v_{n,i} > 0) 법선력이 영이 되는 조건을 의미한다. 접선 방향에서도 유사한 상보 관계가 마찰과 결합되어 정식화된다. 이러한 상보성 체계는 선형 상보 문제(LCP)나 비선형 상보 문제(NCP)로 수치적으로 해결되며, Moreau와 Stewart-Trinkle의 연구가 그 이론적 기반을 제공하였다.

마찰을 포함한 경우

접촉에 마찰이 존재하는 경우, 각 접촉점에서 접선력은 마찰 원추 \|\mathbf{f}_{t,i}\| \le \mu\,f_{n,i}의 제약을 받는다. 이 원추는 일반적으로 비선형 제약이므로, 수치 해법을 위해 다면체 근사(polyhedral approximation)가 자주 사용된다. 원추를 m개의 선형 경계면을 가진 피라미드로 근사하면, 전체 문제는 LCP로 환원되며, Stewart와 Trinkle의 단계별 접근에서 구현된다. 더 정확한 해석을 위해서는 원추를 직접 처리하는 NCP 정식화가 사용된다. Anitescu와 Potra는 이러한 정식화의 수학적 성질을 체계적으로 분석하였다.

수치 해법

다점 접촉의 수치 해는 상보성 문제의 해법에 의존한다. Lemke의 피벗 알고리즘은 LCP의 대표적 해법이며, 이론적으로 유한 단계에서 해를 찾는다. 반복 기반 해법으로는 Projected Gauss-Seidel, Projected Jacobi, Nonlinear Gauss-Seidel 등이 있으며, 대규모 접촉 문제에 적합하다. 다체 동역학 엔진인 Bullet, ODE, MuJoCo, DART, Drake는 이러한 해법을 내부적으로 구현하며, 각 엔진의 선택은 정확성, 안정성, 계산 속도의 균형에 영향을 미친다. Todorov의 convex cone projection 기반 접근은 MuJoCo의 핵심 알고리즘이다.

탄성 접촉과의 결합

강체 기반 해석의 불확정성을 피하기 위해, 각 접촉점에 탄성 응답을 부여하는 접근이 사용된다. 이 경우 각 접촉점의 법선력과 접선력은 침투 깊이와 접선 변위의 함수로 결정되며, 반력 분포가 자연스럽게 고유해진다. 벌칙 기반 구현에서는 접촉 강성이 큰 값으로 설정되어 강체 근사에 근접하며, 동시에 해의 유일성을 유지한다. 이러한 접근은 시뮬레이션 구현에서 널리 채택된다.

로봇 공학적 응용

다점 접촉의 역학 해석은 로봇 공학의 여러 실무에 적용된다. 보행 로봇의 발이 여러 점에서 지면과 접촉하는 경우, 지면 반력 분포의 해석은 균형 유지와 제어의 기초가 된다. 여러 핑거를 가진 그리퍼가 물체를 파지하는 경우, 각 핑거의 접촉력 분포는 파지력 제어의 대상이 된다. 이동 로봇의 네 바퀴와 지면의 접촉에서는 타이어 반력의 분포가 주행 안정성과 조향 응답에 직접 영향을 미친다. 조립 작업에서 여러 접촉점의 상호작용은 접촉 상태의 판정과 힘 제어의 핵심 주제이다.

접촉 상태의 변화

다점 접촉의 해석에서 중요한 도전은 접촉 상태가 시간에 따라 변화할 수 있다는 점이다. 접촉점이 새로 생기거나, 기존 접촉점이 분리되거나, 고착에서 미끄럼으로 또는 그 반대로 전이하는 사건들이 시뮬레이션 중에 발생한다. 각 상태 변화는 상보성 조건의 활성 집합을 재구성해야 하므로, 효율적 상태 관리 알고리즘이 요구된다. 시뮬레이션 엔진은 이러한 상태 전이를 각 시간 단계에서 판정하고 해결하는 구조를 가진다.

접촉 힘 분배의 최적화

과결정 상황에서 접촉 힘의 분배를 결정하기 위한 기준으로 여러 최적화 목적 함수가 제안되었다. 반력의 크기의 합의 최소화, 개별 접촉력의 균형, 마찰 소비의 최소화, 에너지 소산의 최소화 등이 사용되며, 각 기준은 서로 다른 해를 제공한다. 로봇 공학에서는 파지 안정성의 최대화, 바퀴 슬립의 최소화, 지면 반력의 균형 유지와 같은 응용 중심의 목적 함수가 주로 선택된다. 이러한 최적화는 실시간 제어와 오프라인 분석 모두에 활용된다.

본 절의 의의

본 절은 다점 접촉의 역학 해석을 수학적, 수치적, 실무적 관점에서 체계적으로 정리한다. 다점 접촉은 로봇의 보행, 파지, 주행, 조립과 같은 주요 작업의 기본 해석 단위이며, 상보성 정식화와 수치 해법은 이러한 해석의 공통 기반을 제공한다. 본 절의 내용은 이후 절들에서 다룰 파지, 구름 접촉, 접촉 안정성, 시뮬레이터 구현의 실무적 기반을 형성한다.

학습 권장사항

독자는 단순한 네 바퀴 이동 로봇 또는 네 다리 보행 로봇의 정적 반력 분포를 LP로 수식화하고, 중력과 외력의 변화에 따른 반력의 변화를 관찰해 볼 것을 권장한다. 상보성 조건을 포함한 동적 시뮬레이션을 구현하여 접촉 상태의 전이를 재현하는 것은 수치적 이해를 심화시킨다. 다체 시뮬레이션 엔진의 접촉 해석 결과를 단순 예제의 해석 결과와 비교하면 실무적 감각을 얻을 수 있다.

참고 문헌

• Stewart, D. E., & Trinkle, J. C. (1996). An implicit time-stepping scheme for rigid body dynamics with inelastic collisions and Coulomb friction. International Journal for Numerical Methods in Engineering, 39(15), 2673–2691.
• Moreau, J. J. (1988). Unilateral contact and dry friction in finite freedom dynamics. In Nonsmooth Mechanics and Applications, Springer, 1–82.
• Anitescu, M., & Potra, F. A. (1997). Formulating dynamic multi-rigid-body contact problems with friction as solvable linear complementarity problems. Nonlinear Dynamics, 14(3), 231–247.
• Todorov, E. (2014). Convex and analytically-invertible dynamics with contacts and constraints: theory and implementation in MuJoCo. In Proceedings of the IEEE International Conference on Robotics and Automation, 6054–6061.
• Cottle, R. W., Pang, J.-S., & Stone, R. E. (2009). The Linear Complementarity Problem. SIAM.
• Brogliato, B. (2016). Nonsmooth Mechanics: Models, Dynamics and Control (3rd ed.). Springer.
• Featherstone, R. (2008). Rigid Body Dynamics Algorithms. Springer.

version: 1.0