17.25 접촉 동역학의 기본 원리

1. 개요

접촉 동역학은 로봇 시스템이 외부 환경, 작업 대상, 또는 다른 로봇과 물리적으로 접촉할 때 발생하는 힘과 운동의 상호 작용을 다루는 분야이다. 자유 운동 상태에서 로봇은 관절 공간 동역학 방정식에 의해 완전히 기술되지만, 접촉이 발생하는 순간부터는 접촉면에서의 제약, 접촉력, 마찰, 충돌 효과가 추가되어 동역학 구조 자체가 변화한다. 접촉 동역학은 힘 제어, 조립, 연삭, 표면 가공, 악수 등의 협동 작업뿐 아니라 이동 로봇의 보행, 다지 손의 파지, 그리고 인간-로봇 협업의 안전성에 이르기까지 광범위하게 적용된다. 본 절은 접촉 동역학을 지배하는 기본 개념과 수학적 구조를 체계적으로 서술한다.

2. 접촉 상태의 분류

접촉 상태는 접촉면의 상대 운동 특성에 따라 분류된다. 첫째, 점 접촉(point contact)은 두 물체가 단일 점에서 맞닿는 경우이며, 둘째, 선 접촉(line contact)과 면 접촉(surface contact)은 접촉 영역이 선이나 면을 이루는 경우이다. 접촉의 시간 특성에 따라 정적 접촉, 미끄럼 접촉(slipping), 굴림 접촉(rolling)으로 구분되며, 각각은 상대 속도 구속과 마찰 조건에서 상이한 수학적 취급이 요구된다. 또한 접촉 개시 순간에는 충돌(impact) 현상이 발생하여 이산적 속도 변화가 일어나며, 이후 지속 접촉으로 전이되는 과정이 단계적으로 기술되어야 한다.

3. 제약 방정식과 접촉 조건

접촉면이 매끄럽다고 가정할 때, 두 물체 사이의 접촉 조건은 관통 금지 제약(non-penetration constraint)으로 기술된다. 접촉점 사이의 간극 함수를 $\phi(q) \ge 0$ 로 정의하면 다음 조건이 성립한다.

$\phi(q) \ge 0, \qquad \lambda \ge 0, \qquad \phi(q)\,\lambda = 0$

여기서 $\lambda$ 는 접촉 법선 방향의 정규력 크기이다. 이 조건은 Signorini 조건으로 알려져 있으며, 상보성(complementarity) 제약의 형태를 가진다. 즉, 접촉이 발생했을 때( $\phi = 0$ )에만 접촉력이 작용할 수 있으며( $\lambda \ge 0$ ), 비접촉 상태( $\phi > 0$ )에서는 반드시 $\lambda = 0$ 이다.

접촉 야코비안과 제약의 동역학적 통합

접촉점에서 로봇 링크 부분의 속도와 관절 속도의 관계는 접촉 자코비안 $J_c(q)$ 를 통해 연결된다.

$v_c = J_c(q)\,\dot{q}$

이 자코비안은 법선 방향 성분 $J_n$ 과 접선 방향 성분 $J_t$ 로 분해된다. 접촉점의 법선 속도는 $\dot{\phi}(q) = J_n(q)\dot{q}$ 로 주어지며, 지속 접촉이 유지되려면 가속도 수준의 일관성 조건 $\ddot{\phi}(q) = J_n \ddot{q} + \dot{J}_n \dot{q} = 0$ 이 성립해야 한다. 이를 관절 공간 동역학에 포함하면 다음과 같은 라그랑주 승수 형태가 얻어진다.

$M(q)\ddot{q} + C(q,\dot{q})\dot{q} + g(q) = \tau + J_c(q)^\top \lambda_c$

$J_c(q)\ddot{q} + \dot{J}_c(q,\dot{q})\dot{q} = 0$

여기서 $\lambda_c$ 는 접촉 렌치(contact wrench) 벡터이다. 두 방정식을 결합하여 풀면 접촉력과 관절 가속도를 동시에 얻을 수 있다.

4. 마찰 조건과 마찰 원뿔

접선 방향의 상대 운동이 허용되지만 쿨롱 마찰이 존재할 경우, 접촉력은 마찰 원뿔(friction cone) 내부에 위치해야 한다.

$\|\lambda_t\| \le \mu_c\,\lambda_n$

여기서 $\lambda_n, \lambda_t$ 는 각각 법선 및 접선 성분이며, $\mu_c$ 는 마찰 계수이다. 미끄럼이 발생할 때 접선 마찰력은 다음 조건을 만족한다.

$\lambda_t = -\mu_c\,\lambda_n\,\frac{v_t}{\|v_t\|}$

미끄럼과 고착(stick) 상태는 동일한 방정식에서 상보성 조건으로 전환되며, 전체 접촉 조건은 선형 상보성 문제(linear complementarity problem, LCP)로 정식화된다. 정밀한 해석에서는 마찰 원뿔의 다면체 근사가 사용되어 LCP의 크기가 조정된다.

5. 접촉 모델의 두 접근

접촉 동역학의 수학적 처리는 크게 두 가지 접근으로 나뉜다.

5.1 강체 접촉 모델

강체 접촉 모델은 Signorini 조건을 엄격히 따르며, 접촉면을 이상적으로 딱딱한 제약으로 간주한다. 이 접근은 구속 조건을 해석적으로 부과하므로 이론적 정확성이 높지만, 접촉의 개시와 해제 시점에서 이산적 사건이 발생하고 LCP 해법이 요구되어 수치적 계산이 복잡해진다. Stewart-Trinkle 방식과 Anitescu-Potra 방식이 대표적 알고리즘이다.

5.2 규정된 준수 접촉 모델

규정된 준수(compliance) 접촉 모델은 접촉면이 탄성 변형을 일으킨다고 가정하여 관통 깊이 $d$ 에 따라 접촉력이 결정되도록 한다. 대표적 형태는 Kelvin-Voigt 모델과 Hunt-Crossley 모델이며, 후자는 에너지 손실을 자연스럽게 반영한다.

$\lambda_n = k\,d^n + b\,d^n\,\dot{d}$

여기서 $k, b, n$ 은 접촉 강성, 감쇠, 비선형 지수이다. 이 접근은 연속적이고 상미분 방정식 기반 시뮬레이션이 가능하지만, 강성 $k$ 가 큰 경우 강성 시스템이 되어 적분기의 안정성이 요구된다.

충돌 해석

접촉이 개시되는 순간에는 유한한 시간 내에 이산적 속도 변화가 발생하며, 이는 충돌(impact) 현상으로 기술된다. 충돌 전후의 속도는 운동량 보존 원리와 반발 계수(coefficient of restitution) $e$ 를 통해 연결된다.

$M(q)\bigl(\dot{q}^+ - \dot{q}^-\bigr) = J_n^\top \hat{\lambda}_n$

$J_n \dot{q}^+ = -e\,J_n \dot{q}^-$

여기서 $\hat{\lambda}_n$ 은 충격 임펄스이다. 충돌 모델은 에너지 감소 특성, 마찰 효과, 접촉점의 기하 조건을 포함하도록 확장된다. 다중 접촉점이 동시에 형성되는 경우 Moreau의 스위핑 프로세스 이론이 사용된다.

접촉 안정성과 일관성

접촉이 유지되려면 접촉력이 물리적으로 가능한 영역 내에 있어야 하며, 관통 금지 조건, 마찰 원뿔 조건, 양의 정규력 조건이 동시에 만족되어야 한다. 이러한 조건이 유지되지 않으면 접촉이 해제되거나 미끄럼으로 전환된다. 제어 관점에서는 접촉 안정성을 보장하기 위해 제어기가 이러한 조건을 명시적으로 인식하고 안전한 궤적을 계획해야 한다.

접촉 동역학의 관절 공간 통합

접촉이 존재하는 상황에서 관절 공간 동역학 방정식은 다음과 같이 확장된다.

$M(q)\ddot{q} + C(q,\dot{q})\dot{q} + g(q) + \tau_f = \tau + J_c(q)^\top \lambda_c$

여기서 $J_c^\top \lambda_c$ 는 외부 접촉이 관절 공간으로 투영된 등가 토크이다. 이 형식은 자유 운동과 접촉 운동을 통일적으로 다루며, 접촉이 없을 때 $\lambda_c = 0$ 으로 자연스럽게 환원된다. 제어기 설계에서는 접촉 상태의 발생 여부에 따라 제어 법칙을 전환하거나, 연속적 힘 제어 법칙을 통해 두 상태를 매끄럽게 연결하는 방식이 채택된다.

6. 본 절의 의의

본 절은 로봇이 환경과 상호 작용할 때 필수적으로 요구되는 접촉 동역학의 기본 원리를 체계적으로 정리한다. 접촉 조건, 제약의 동역학적 통합, 마찰 원뿔, 충돌, 강체 및 준수 모델은 모두 이후 절에서 다룰 접촉력 모델링, 힘 제어, 임피던스 및 어드미턴스 제어, 다체 접촉 해석, 충돌 시뮬레이션의 수학적 기초를 제공한다. 접촉 동역학의 정확한 이해는 안전한 인간-로봇 협업과 정밀한 상호 작용 작업의 구현에 필수적이다.

7. 학습 권장사항

독자는 단일 구가 평면과 수직으로 충돌하는 간단한 모델에 대해 강체 접촉 모델과 Hunt-Crossley 준수 모델을 각각 구현하고, 두 방식의 시뮬레이션 결과를 비교해 보는 것을 권장한다. 또한 평면 2자유도 매니퓰레이터가 벽에 접촉하는 상황에서 접촉 자코비안을 계산하고 라그랑주 승수 방법을 통해 접촉력을 풀어 보는 연습은 접촉 조건의 수학적 통합을 이해하는 데 도움이 된다. 마찰 원뿔의 다면체 근사와 LCP 해법의 관계를 이론적으로 살펴보는 것도 유익하다.

8. 참고 문헌

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.
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.
Hunt, K. H., & Crossley, F. R. E. (1975). Coefficient of restitution interpreted as damping in vibroimpact. Journal of Applied Mechanics, 42(2), 440–445.
Brogliato, B. (1999). Nonsmooth Mechanics: Models, Dynamics and Control. Springer.
Moreau, J. J. (1988). Unilateral contact and dry friction in finite freedom dynamics. In Nonsmooth Mechanics and Applications, 1–82. Springer.
Featherstone, R. (2008). Rigid Body Dynamics Algorithms. Springer.
Murray, R. M., Li, Z., & Sastry, S. S. (1994). A Mathematical Introduction to Robotic Manipulation. CRC Press.

version: 1.0