17.26 접촉력과 접촉 모멘트의 모델링

1. 개요

접촉 동역학의 기본 원리가 접촉 조건과 제약의 통합 방식을 다룬 반면, 실제 제어와 시뮬레이션에서는 접촉점에서 발생하는 힘과 모멘트를 정량적으로 기술하는 구체적 수학 모델이 요구된다. 접촉력과 접촉 모멘트는 접촉면의 기하, 재질, 변형 특성, 마찰 조건, 상대 운동의 형태에 의해 결정되며, 정확한 모델링은 힘 제어, 파지 안정성 해석, 조립 작업, 보행 로봇의 지면 반력 해석, 충돌 안전성 평가에 직접적인 영향을 미친다. 본 절은 접촉점에서 작용하는 렌치(wrench) 표현, 법선 및 접선 성분의 모델, 점·선·면 접촉에서 모멘트가 발생하는 조건, 분포 압력의 집중 표현, 접촉 렌치의 수치적 계산 방법을 체계적으로 제시한다.

2. 접촉 렌치의 정의

접촉점이 하나의 공간 점 \mathbf{p}_c로 이상화될 때, 접촉점에 작용하는 물리적 상호 작용은 힘 \mathbf{f}_c \in \mathbb{R}^3과 모멘트 \mathbf{n}_c \in \mathbb{R}^3의 결합인 접촉 렌치 \mathcal{F}_c = [\mathbf{f}_c^\top, \mathbf{n}_c^\top]^\top \in \mathbb{R}^6로 표현된다. 접촉이 이상적인 점 접촉인 경우 접촉점 주위의 모멘트는 영이 되지만, 접촉 영역이 유한한 면이라면 분포 압력과 분포 마찰에서 유도되는 모멘트 성분이 함께 나타난다. 이를 통해 접촉은 점 접촉, 마찰 점 접촉, 소프트 핑거 접촉, 면 접촉 등의 유형으로 분류되며, 각 유형은 렌치 공간에서 서로 다른 구속 구조를 가진다.

3. 법선 접촉력의 모델

접촉면에 수직한 법선 방향의 힘은 관통 금지 제약과 물체 변형의 관계로부터 유도된다. 가장 단순한 형태는 선형 탄성 모델이다.

\lambda_n = k\,d,\qquad d = \max(0, -\phi(q))

여기서 d는 관통 깊이, k는 접촉 강성이다. 감쇠 효과를 포함하는 Kelvin-Voigt 모델은 다음과 같다.

\lambda_n = k\,d + b\,\dot{d}

이 모델은 구현이 단순하지만 분리 순간에 비물리적 음의 힘을 허용하는 한계가 있다. 이를 보완하는 Hunt-Crossley 모델은 비선형 형태로 다음과 같이 정의된다.

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

여기서 \chi는 감쇠 계수, n은 접촉 기하에 따른 지수이다. 구 접촉 해석에서는 Hertz 이론에 의해 n = 3/2가 나타나며, 에너지 보존과 분리 조건이 자연스럽게 만족된다. 이러한 비선형 모델은 분포 압력을 집중 힘으로 근사한 결과로서, 접촉 동역학 시뮬레이터에서 널리 채택된다.

접선 마찰력의 모델

접선 방향의 힘은 접촉면 사이의 상대 운동과 마찰 조건에 의해 결정된다. 고착 조건에서는 접선 상대 속도가 영이므로 마찰력은 마찰 원뿔 내부의 임의 값을 취하며, 미끄럼 조건에서는 쿨롱 법칙에 의해 크기가 결정된다.

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

여기서 v_t는 접선 상대 속도이다. 쿨롱 모델은 속도가 영점에 접근하면서 불연속이 되므로 수치적 안정성을 확보하기 위해 다음과 같은 규화가 자주 사용된다.

\lambda_t = -\mu_c\,\lambda_n\,\tanh(\kappa \|v_t\|)\,\frac{v_t}{\|v_t\|}

보다 정교한 접근에서는 LuGre, Dahl, Bristle 등의 동적 모델이 도입되어, 미시적 변위 z를 상태로 포함하는 상미분 방정식 형태로 접선 마찰력이 계산된다. 이 방법은 스틱-슬립 전환을 연속적으로 포착하며, 정밀 접촉 시뮬레이션과 햅틱 렌더링에서 이점이 크다.

마찰 원뿔과 접촉 렌치 공간

마찰이 존재하는 접촉에서 실현 가능한 접촉력 집합은 3차원 원뿔로 표현된다.

\mathcal{C}_\mu = \{\lambda \in \mathbb{R}^3 \mid \lambda_n \ge 0,\; \sqrt{\lambda_{t1}^2 + \lambda_{t2}^2} \le \mu_c\,\lambda_n\}

수치 해석과 최적화에서는 이 원뿔을 다면체로 근사하여 선형 제약 조건으로 변환하며, k개의 생성 벡터를 통해 표현된다. 선형 근사는 LCP, QP, LP 기반 해법과 결합되어 다체 접촉 해석을 가능케 한다. 다수의 접촉점이 존재할 때 전체 접촉 렌치는 각 접촉점의 개별 렌치의 합이며, 합력 렌치 공간은 개별 원뿔의 Minkowski 합으로 기술된다. 이 구조는 파지 안정성 해석에서 결정적 역할을 한다.

4. 모멘트 성분과 소프트 핑거 접촉

이상적 점 접촉에서는 접촉점 주위의 모멘트가 영이지만, 고무와 같은 탄성 재질이나 유한한 접촉 영역을 가지는 경우 접촉 영역 전반에 걸친 분포 압력이 존재하여 모멘트가 발생한다. 특히 파지에서 자주 사용되는 소프트 핑거 접촉 모델은 법선 축 주위의 마찰 토크를 허용한다.

|n_c| \le \mu_t\,\lambda_n

여기서 \mu_t는 비틀림 마찰 계수이다. 이 모델은 이상적 점 접촉보다 현실적 파지 특성을 반영하며, 파지력 폐쇄(force-closure) 판정에서 추가적인 자유도를 제공한다. 면 접촉에서는 접촉 영역의 압력 분포를 통해 법선 모멘트, 롤링 저항, 굴림 마찰 등이 파생되며, 이는 보행 로봇의 발 접촉이나 연삭 공정의 공구 접촉 모델링에서 필수적이다.

분포 압력의 집중화

유한한 접촉 면적에 걸친 분포 압력은 해석과 시뮬레이션의 편의를 위해 등가 집중 힘과 집중 모멘트의 쌍으로 변환된다. 압력 분포 p(\mathbf{r})와 마찰 분포 \mathbf{t}(\mathbf{r})가 주어지면 등가 렌치는 다음과 같이 계산된다.

\mathbf{f}_c = \int_{\mathcal{A}}\bigl(p(\mathbf{r})\mathbf{n} + \mathbf{t}(\mathbf{r})\bigr)\,dA

\mathbf{n}_c = \int_{\mathcal{A}} \mathbf{r}\times\bigl(p(\mathbf{r})\mathbf{n} + \mathbf{t}(\mathbf{r})\bigr)\,dA

여기서 \mathcal{A}는 접촉 영역, \mathbf{n}은 접촉면 법선이다. 특히 법선력의 작용점은 압력 중심(center of pressure, CoP)으로 정의되며, 보행 로봇에서 발 접촉의 안정성 평가에 활용된다. CoP가 발 지지 영역의 경계를 벗어나면 기울어짐이 발생하므로, 이 지점은 보행 제어의 핵심 지표로 사용된다.

감쇠와 에너지 손실

접촉에서 에너지 손실은 재료의 히스테리시스, 마찰 소산, 그리고 파동의 방사로 인해 발생한다. Hunt-Crossley 모델은 재료 히스테리시스에 의한 손실을 자연스럽게 반영하는 특성을 가지며, 반발 계수와 감쇠 계수 사이의 관계가 다음과 같이 유도된다.

e \approx 1 - \frac{\chi}{k}\,v_0

여기서 v_0는 충돌 직전 법선 속도이다. 마찰에 의한 소산은 접선 상대 속도와 접선 힘의 내적으로 계산되며, 패시비티 기반 제어 설계에서 수동 소자로 취급된다. 정확한 에너지 손실 기술은 충돌 안정성과 햅틱 렌더링의 현실성에 직접적인 영향을 미친다.

5. 측정과 외란 관측

실제 시스템에서 접촉력과 접촉 모멘트는 엔드 이펙터에 장착된 6축 힘/토크 센서, 관절 토크 센서, 또는 외란 관측기를 통해 측정 또는 추정된다. 6축 힘/토크 센서는 직접 측정을 제공하지만 센서 드리프트, 온도 의존성, 자중 보상이 요구되며, 관절 토크 센서 기반 추정은 전체 매니퓰레이터 동역학 모델을 이용하여 \tau - \hat{\tau}_{\text{model}} \approx J^\top \mathcal{F}_c의 관계로부터 접촉 렌치를 역추정한다. 후자는 별도의 센서 없이도 접촉력을 추정할 수 있는 장점이 있으나, 모델 불확실성과 마찰에 의한 편향에 취약하다.

6. 제어기 설계와의 연결

접촉력 모델은 힘 제어기, 임피던스 제어기, 어드미턴스 제어기의 설계에 직접 반영된다. 예컨대 임피던스 제어에서 환경의 접촉 강성을 k_e라 하고 제어기가 원하는 임피던스를 k_d, b_d, m_d로 설정할 때, 접촉력은 엔드 이펙터의 편차와 환경 강성의 곱으로 모델링되며, 제어기는 이 관계를 통해 원하는 상호 작용 거동을 구현한다. 접촉 모델이 부정확하면 힘 오버슈트, 접촉 손실, 진동 불안정성이 나타나므로, 모델과 제어기의 일관성이 요구된다.

7. 본 절의 의의

본 절은 접촉점에서 실제로 발생하는 힘과 모멘트의 수학적 모델링을 체계적으로 제시하여, 이상화된 접촉 조건을 구체적 수치 해석과 제어 설계에 연결하는 다리를 제공한다. 법선과 접선 힘의 모델, 마찰 원뿔, 모멘트 성분, 분포 압력의 집중 표현, 에너지 손실과 측정 방법은 후속 절에서 다룰 힘 제어, 임피던스 제어, 다체 접촉 해석, 파지 안정성 평가의 수학적 기반이 된다.

8. 학습 권장사항

독자는 Hertz 접촉 이론을 이용하여 구-평면 접촉에서 Hunt-Crossley 모델을 유도하고, 반발 계수와 감쇠 계수의 관계를 수치적으로 검증해 볼 것을 권장한다. 또한 3점 지지 파지에서 마찰 원뿔의 다면체 근사를 사용하여 파지력 폐쇄 여부를 LP로 판정하는 실습은 접촉 렌치 공간의 구조를 구체적으로 이해하는 데 도움이 된다. 6축 힘/토크 센서의 자중 보상 절차를 구현하고 계측값에서 접촉 렌치를 추출해 보는 실험은 이론과 실제의 간격을 좁히는 좋은 기회가 된다.

9. 참고 문헌

• Johnson, K. L. (1985). Contact Mechanics. Cambridge University Press.
• Hunt, K. H., & Crossley, F. R. E. (1975). Coefficient of restitution interpreted as damping in vibroimpact. Journal of Applied Mechanics, 42(2), 440–445.
• Mason, M. T. (2001). Mechanics of Robotic Manipulation. MIT Press.
• Murray, R. M., Li, Z., & Sastry, S. S. (1994). A Mathematical Introduction to Robotic Manipulation. CRC Press.
• Bicchi, A., & Kumar, V. (2000). Robotic grasping and contact: A review. IEEE International Conference on Robotics and Automation.
• Flores, P., Machado, M., Silva, M. T., & Martins, J. M. (2011). On the continuous contact force models for soft materials in multibody dynamics. Multibody System Dynamics, 25(3), 357–375.
• Haddadin, S., De Luca, A., & Albu-Schäffer, A. (2017). Robot collisions: A survey on detection, isolation, and identification. IEEE Transactions on Robotics, 33(6), 1292–1312.

version: 1.0