18.16 접촉 강성 및 감쇠 모델링

1. 개요

접촉 강성과 감쇠의 모델링은 로봇과 환경의 상호작용을 동역학적으로 기술하기 위한 핵심 요소이다. 접촉 강성은 접촉력이 상대 변위에 어떻게 반응하는지를 정량화하며, 접촉 감쇠는 접촉 과정에서 발생하는 에너지 소산을 기술한다. 이 두 요소는 접촉 진동, 충격 응답, 힘 제어 폐회로의 안정성, 시뮬레이터의 수치적 거동에 결정적 영향을 미친다. 본 절은 접촉 강성의 정의와 선형 및 비선형 모델, 감쇠의 물리적 기원과 모델, 스프링-감쇠 근사의 구조, 고차 비선형 모델, 매개변수 식별, 로봇 공학적 의의를 체계적으로 정리한다.

2. 접촉 강성의 정의

접촉 강성은 접촉력의 상대 접근 거리에 대한 변화율로 정의된다.

k_c = \frac{dF_n}{d\delta}

여기서 F_n은 법선 접촉력, \delta는 접근 거리이다. 탄성 접촉의 경우 강성은 재료의 탄성 계수, 접촉 기하학, 접근 거리에 의존하며, 헤르츠 구 접촉에서는 k_c \propto \delta^{1/2}의 비선형 관계를 가진다. 접선 방향의 강성은 접선력과 접선 변위의 관계로 정의되며, 법선 강성과는 다른 값을 가진다. 이러한 강성은 국소 접촉 영역의 탄성 응답을 특징짓는 핵심 매개변수이다.

선형 스프링 모델

가장 단순한 접촉 강성 모델은 선형 스프링이다. 이 모델에서 접촉력은 침투 깊이에 비례하며, F_n = k_c\,\delta로 표현된다. 선형 모델은 구현과 해석이 단순하고 수치 적분에서 안정적이지만, 접촉 기하학의 영향을 반영하지 못하며 하중 범위의 변화에 둔감하다는 단점이 있다. 또한 하중 제거 시 침투 깊이가 영이 되는 순간에 접촉력이 갑자기 영으로 떨어지므로, 이탈 순간의 연속성이 부자연스럽다. 그럼에도 불구하고 계산 효율과 단순성 덕분에 많은 실시간 로봇 시뮬레이션에서 기본 선택으로 사용된다.

비선형 헤르츠 기반 모델

헤르츠 이론을 기반으로 한 비선형 접촉 모델은 F_n = k\,\delta^{3/2}의 형태로 표현되며, 구-평면이나 구-구 접촉의 실제 탄성 응답을 정확히 기술한다. 이 형태는 강성이 접근 거리에 따라 증가하는 자연스러운 거동을 제공하며, 접촉의 이탈 순간에서 힘이 매끄럽게 영으로 수렴하는 장점이 있다. 비선형 강성 모델은 수치 적분에서 다소 높은 계산 비용을 요구하지만, 고품질 시뮬레이션에서 물리적 현실성이 우수하다. 지수 3/2는 구 접촉에 해당하며, 다른 접촉 기하에서는 기하학에 따라 다른 지수가 적용된다.

접촉 감쇠의 물리적 기원

접촉 감쇠는 접촉 과정에서 에너지가 소산되는 현상이며, 여러 물리적 기원을 가진다. 첫째, 재료의 내부 감쇠는 변형 에너지의 일부가 재료 내부의 마이크로 구조 운동으로 소산되는 과정을 기술한다. 둘째, 점탄성 재료는 시간 의존적 변형 응답을 보이며, 이에 따른 이력 루프가 감쇠로 나타난다. 셋째, 소성 변형은 비가역적 에너지 소산을 유발하며, 충격 과정에서 큰 감쇠 기여를 한다. 넷째, 미세 마찰과 표면의 미시적 미끄럼은 기계적 에너지를 열로 변환한다. 이러한 기원들은 접촉 감쇠 모델의 선택에서 물리적 근거를 제공한다.

선형 스프링-감쇠기 모델

가장 단순한 접촉 감쇠 모델은 선형 점성 감쇠기이며, 접촉력은 침투 깊이에 비례하는 탄성 성분과 침투 속도에 비례하는 감쇠 성분의 합으로 표현된다.

F_n = k_c\,\delta + c_c\,\dot{\delta}

이 모델은 단순하지만 몇 가지 비물리적 특성을 가진다. 접촉 개시 순간에 속도가 영이 아니면 감쇠력에 의해 접촉력이 불연속적으로 시작되며, 이탈 순간에는 탄성력이 영이지만 감쇠력이 음의 방향으로 물체를 끌어당기는 비물리적 거동이 발생할 수 있다. 이러한 단점에도 불구하고, 단순성 덕분에 널리 사용된다.

3. Kelvin-Voigt와 Hunt-Crossley 모델

Hunt와 Crossley는 1975년 Coefficient of restitution interpreted as damping in vibroimpact에서 선형 스프링-감쇠기 모델의 비물리적 특성을 보완하는 비선형 감쇠 모델을 제안하였다. 이 모델에서 접촉력은 다음과 같이 표현된다.

F_n = k\,\delta^{n} + \lambda\,\delta^{n}\,\dot{\delta}

여기서 감쇠항이 침투 깊이의 함수로 변조되어, 접촉 개시와 이탈 순간에서 힘이 매끄럽게 영이 된다. 이 모델은 반발 계수와 감쇠 매개변수 사이의 이론적 관계를 제공하며, 실험 데이터와의 적합성이 우수하다. Hunt-Crossley 모델은 로봇 시뮬레이션과 다체 동역학 해석에서 접촉 감쇠의 표준적 선택 중 하나이다.

고차 비선형 모델

더 고차의 비선형 모델은 점탄성, 히스테리시스, 의존 시간 응답을 포함하는 여러 형태를 가진다. Lankarani와 Nikravesh는 에너지 손실을 반발 계수의 함수로 표현한 모델을 제시하였으며, 다체 충돌 해석에서 널리 사용된다. 점탄성 접촉 모델은 Maxwell, Kelvin-Voigt, Zener와 같은 일반적 점탄성 단위의 결합으로 구성되며, 소프트 접촉과 유연 표면 해석에 적합하다. 이들 모델은 복잡한 이력 거동을 정확히 재현하는 대신 매개변수의 수와 식별 난이도가 증가한다.

매개변수의 식별

접촉 강성과 감쇠의 매개변수는 충격 응답 실험, 정적 압입 실험, 주파수 응답 실험을 통해 식별된다. 충격 실험에서 얻어진 속도 이력과 반발 계수는 감쇠 매개변수의 정보를 제공하며, 정적 압입 실험은 강성 매개변수를 직접 제공한다. 주파수 응답 실험은 접촉 공진과 손실 각도를 측정하여 점탄성 매개변수를 식별한다. 로봇 공학에서는 힘 센서와 위치 센서를 결합한 측정 장치가 이러한 실험의 기본 장비로 사용된다.

접촉 강성과 힘 제어의 관계

힘 제어 매니퓰레이터의 폐회로 안정성은 환경의 접촉 강성에 크게 의존한다. 높은 접촉 강성은 폐회로의 대역폭을 제한하며, 불안정성의 원인이 될 수 있다. 임피던스 제어와 어드미턴스 제어는 접촉 강성에 대응하여 설계되며, 접촉 강성의 추정치는 제어기 이득 선택에 직접 사용된다. 이러한 맥락에서 접촉 강성의 정확한 추정과 온라인 적응은 힘 제어 성능의 핵심 요소이다.

시뮬레이션에서의 구현

로봇 시뮬레이터는 접촉 강성과 감쇠를 벌칙 함수(penalty function)로 구현하는 것이 일반적이다. 접촉 제약을 엄격히 만족시키는 대신, 침투를 허용하면서 벌칙력을 부과하는 이 방법은 구현이 단순하고 계산 효율이 우수하다. 벌칙 계수는 접촉 강성으로 해석될 수 있으며, 침투 속도에 비례하는 감쇠가 추가된다. 벌칙 기반 구현의 안정성은 시간 단계와 강성 계수의 관계에 민감하며, 이는 수치 기법의 선택에 영향을 미친다.

로봇 공학적 응용

접촉 강성과 감쇠 모델은 로봇 공학의 여러 실무에 적용된다. 힘 제어 매니퓰레이터의 안정성 해석, 이동 로봇 바퀴와 지면 사이의 지지 응답, 파지 작업의 미세 힘 조절, 촉각 센서의 기계적 응답 해석이 대표적 예이다. 또한 시뮬레이터에서 접촉 강성의 비현실적 선택은 비물리적 튕김 현상, 수치적 불안정, 에너지 폭발을 유발할 수 있으므로, 신중한 매개변수 조정이 요구된다. 실무에서는 검증된 기본 설정과 사용자 조정의 조합이 일반적이다.

본 절의 의의

본 절은 접촉 강성과 감쇠의 모델링을 체계적으로 정리하고, 로봇 공학에서의 실무적 의의를 제시한다. 접촉 강성은 힘 응답의 크기를, 감쇠는 에너지 소산과 안정화를 결정하는 핵심 요소이며, 두 요소의 적절한 선택은 시뮬레이션의 품질과 제어의 성능을 결정한다. 본 절의 내용은 이후 절들에서 다룰 접촉 기하학, 충돌 해석, 소프트 접촉, 수치 시뮬레이션의 기반이 된다.

학습 권장사항

독자는 동일한 충격 실험에 대해 선형 스프링-감쇠기, Hertz 비선형 강성, Hunt-Crossley 모델을 각각 구현하고, 접촉 힘의 시간 이력과 반발 계수를 비교해 볼 것을 권장한다. 각 모델의 접촉 개시와 이탈 순간의 연속성을 관찰하면, 모델 선택의 물리적 근거를 구체적으로 이해할 수 있다. 또한 힘 제어 시뮬레이션에서 환경 강성을 변화시키면서 폐회로의 안정성을 평가하는 실습은 제어 설계의 실무 감각을 제공한다.

참고 문헌

• Hunt, K. H., & Crossley, F. R. E. (1975). Coefficient of restitution interpreted as damping in vibroimpact. Journal of Applied Mechanics, 42(2), 440–445.
• Lankarani, H. M., & Nikravesh, P. E. (1990). A contact force model with hysteresis damping for impact analysis of multibody systems. Journal of Mechanical Design, 112(3), 369–376.
• Johnson, K. L. (1985). Contact Mechanics. Cambridge University Press.
• 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.
• Marhefka, D. W., & Orin, D. E. (1999). A compliant contact model with nonlinear damping for simulation of robotic systems. IEEE Transactions on Systems, Man, and Cybernetics, Part A, 29(6), 566–572.
• Machado, M., Moreira, P., Flores, P., & Lankarani, H. M. (2012). Compliant contact force models in multibody dynamics: evolution of the Hertz contact theory. Mechanism and Machine Theory, 53, 99–121.
• Gilardi, G., & Sharf, I. (2002). Literature survey of contact dynamics modelling. Mechanism and Machine Theory, 37(10), 1213–1239.

version: 1.0