18.12 마찰의 수치적 모델링과 시뮬레이션

1. 개요

마찰의 수치적 모델링과 시뮬레이션은 앞서 논의된 여러 마찰 모델을 수치 해법으로 안정적이고 정확하게 구현하는 기법을 다룬다. 마찰 현상은 불연속성, 비평활성, 강성, 수동성 제약과 같은 수치적 도전 과제를 동반하며, 이를 적절히 처리하지 못하면 시뮬레이션 결과가 비물리적 진동, 에너지 불균형, 상태 드리프트를 나타낸다. 본 절은 이러한 도전 과제를 체계적으로 정리하고, 주요 수치 기법인 평활화 근사, 속도 불감대, 상보성 정식화, 암시적 적분, 이벤트 구동 기법을 비교하며, 로봇 시뮬레이션의 실무적 구현에 적용되는 지침을 제공한다.

2. 수치적 도전 과제

마찰 모델의 수치적 구현은 여러 도전 과제를 동반한다. 첫째, 부호 함수의 불연속성은 속도가 영을 지나는 순간에 마찰력이 급격히 변화시키며, 이는 표준 명시적 적분기에서 고주파 진동과 에너지 오차를 유발한다. 둘째, 상태 변수 기반 동적 모델의 강성은 강성 계수 $\sigma_0$ 이 큰 경우에 시간 상수가 매우 작아져, 수치 적분의 시간 단계를 크게 제한한다. 셋째, 고착 상태의 정확한 유지와 전이의 판정은 명시적 방법으로는 어려우며, 상태 판정 논리가 요구된다. 넷째, 수동성 조건의 수치적 위반은 에너지 드리프트와 비물리적 거동을 유발한다. 이러한 문제들을 종합적으로 고려한 수치 기법의 선택이 요구된다.

3. 평활화 근사

평활화 근사는 부호 함수의 불연속성을 매끄러운 함수로 대체하여 표준 적분기와 호환되도록 하는 단순하고 효과적인 기법이다. 자주 사용되는 근사는 $\tanh(v/\varepsilon)$ , $v/\sqrt{v^2 + \varepsilon^2}$ , $2\arctan(v/\varepsilon)/\pi$ 등이며, 매개변수 $\varepsilon$ 의 크기가 근사의 정확성과 수치 안정성을 결정한다. 작은 $\varepsilon$ 은 불연속성을 더 정확히 근사하지만 수치 강성을 유발하며, 큰 $\varepsilon$ 은 수치 안정성을 향상시키지만 저속 영역의 거동을 왜곡한다. 평활화 기법은 구현이 단순하고 미분 가능성이 보장되므로, 제어 설계와 민감도 해석에서 유리하다.

4. 속도 불감대 기법

속도 불감대 기법은 상대 속도의 크기가 특정 임계값 이하인 경우를 고착 상태로 간주하여, 이 구간에서 마찰력을 외부 접선력과 평형을 이루도록 설정하는 방법이다. 이 기법은 Karnopp의 고전적 연구에서 체계화되었으며, 평활화 근사보다 고착 상태를 더 정확히 포착할 수 있다는 장점이 있다. 구현은 조건부 스위칭 논리를 포함하므로 매끄러움은 감소하지만, 단순한 수치 적분기와 결합하여 실용적 결과를 제공한다. 불감대 크기의 선택은 시스템 특성과 시간 단계에 맞게 조정되어야 한다.

5. 상보성 정식화

상보성 정식화는 고착과 미끄럼의 이상적 구분을 수학적으로 엄밀하게 표현하는 접근이다. 이 방법은 각 시간 단계에서 상대 속도와 마찰력 사이의 상보성 조건을 만족하는 해를 찾도록 정식화되며, 선형 상보 문제(LCP)나 비선형 상보 문제(NCP)로 수치적으로 해결된다. Moreau의 소인 과정과 Stewart-Trinkle의 단계별 방법은 대표적 예이며, 충격을 포함한 비평활 역학의 일관된 수치 처리에 유용하다. 상보성 정식화는 정확성이 높지만, 각 시간 단계에서 상보 문제를 풀어야 하므로 계산 비용이 크고 구현이 복잡하다.

6. 암시적 적분 기법

상태 변수 기반 동적 마찰 모델의 강성을 다루기 위해 암시적 적분 기법이 자주 사용된다. 후진 오일러(backward Euler)와 같은 1차 암시적 방법은 단순하고 안정성이 뛰어나며, 2차 이상의 Runge-Kutta 암시적 방법은 정확성을 향상시킨다. 암시적 방법은 각 시간 단계에서 비선형 방정식을 풀어야 하므로 계산 비용이 증가하지만, 큰 시간 단계에서도 안정적으로 적분할 수 있다는 장점이 있다. 로봇 시뮬레이션에서 LuGre 모델의 $\sigma_0$ 이 매우 큰 경우, 암시적 방법의 사용이 특히 권장된다.

7. 이벤트 구동 기법

이벤트 구동 기법은 고착과 미끄럼의 전이를 정확한 이벤트로 취급하여, 전이 순간에 적분을 중단하고 상태를 재초기화한 후 새로운 상태에서 적분을 계속하는 방법이다. 이 접근은 전이의 물리적 의미를 명확히 포착하며, 에너지 손실과 상태 드리프트를 최소화한다. 그러나 이벤트 검출 알고리즘의 구현 복잡도가 크고, 이벤트가 빈번한 경우 계산 비용이 급증할 수 있다. 정밀한 해석이 요구되는 연구 목적의 시뮬레이션에서 주로 사용된다.

8. 에너지 보존과 수동성의 수치적 유지

시뮬레이션의 장시간 거동은 수치 적분이 에너지 보존과 수동성을 얼마나 정확히 유지하는지에 따라 결정된다. 명시적 방법은 에너지 오차가 시간에 따라 누적될 수 있으며, 특히 마찰과 같은 소산 요소가 있는 경우에 인위적 에너지 주입이 발생할 수 있다. 대칭적 적분기와 변분 적분기는 이러한 문제를 완화하며, 장시간 해석의 정확성을 향상시킨다. 마찰 모델의 수치적 구현에서는 수동성을 위반하지 않도록 시간 단계를 선택하고, 에너지 균형을 정기적으로 검증하는 것이 권장된다.

9. 실시간 시뮬레이션의 고려 사항

실시간 로봇 시뮬레이션에서는 계산 예산이 엄격히 제한되므로, 정밀도와 계산 비용의 균형이 중요하다. 평활화 기법과 명시적 적분의 결합은 계산 비용이 작고 구현이 단순하여 가장 널리 사용된다. 고품질이 요구되는 응용에서는 속도 불감대와 단순 암시적 방법의 결합이 선호된다. 상보성 정식화와 이벤트 구동 기법은 일반적으로 실시간 예산을 초과하므로, 오프라인 분석과 검증에서 사용된다. 로봇 시뮬레이터의 구현은 이러한 트레이드오프를 고려한 설계를 요구한다.

10. 매개변수와 시간 단계의 상호 작용

마찰 모델의 매개변수와 시간 단계의 선택은 상호 의존적이다. 상태 변수 기반 모델에서 시간 상수 $F_c/\sigma_0$ 이 시간 단계보다 작은 경우에는 수치 강성이 발생하며, 이 경우 시간 단계를 줄이거나 암시적 방법을 사용해야 한다. 평활화 매개변수 $\varepsilon$ 은 일반적으로 시간 단계와 호환되는 크기로 선택된다. 이러한 상호 작용에 대한 이해는 시뮬레이션 설정의 안정성을 확보하는 데 필수적이며, 과도한 시간 단계 축소 없이 정확한 결과를 얻는 실무적 감각을 제공한다.

11. 검증과 수치적 품질 관리

수치 마찰 시뮬레이션의 품질은 해석적 해가 존재하는 단순 예제와의 비교, 에너지 균형 검사, 공간적 대칭성 검사, 시간 단계 민감도 분석을 통해 검증된다. Dahl 모델의 해석적 해, 고착-미끄럼 시스템의 주기와 진폭, 주파수 응답의 이론적 예측 등이 검증의 기준으로 사용된다. 이러한 검증 절차는 구현의 정확성을 보장하며, 이후의 복잡한 시스템 시뮬레이션에 대한 신뢰를 제공한다.

12. 본 절의 의의

본 절은 마찰 현상의 수치적 모델링과 시뮬레이션에서 발생하는 도전 과제와 주요 기법을 체계적으로 정리한다. 마찰의 수치 구현은 선택된 마찰 모델과 수치 기법의 조합에 따라 결과가 크게 달라지며, 적절한 선택과 검증은 시뮬레이션 품질의 결정적 요소이다. 본 절의 내용은 이후 절들에서 다룰 접촉 역학의 수치 시뮬레이션과 시뮬레이터 구현에 직접 연결되는 실무적 기반을 제공한다.

13. 학습 권장사항

독자는 동일한 마찰 모델에 대해 평활화 근사, 속도 불감대, 암시적 적분을 각각 구현하고, 저속 고착-미끄럼 시뮬레이션에서 결과를 비교해 볼 것을 권장한다. 시간 단계를 변화시키면서 에너지 균형과 해의 수렴을 관찰하면, 수치 기법의 특성을 정량적으로 이해할 수 있다. 또한 Dahl 모델의 해석적 해를 기준으로 수치 오차를 측정하는 실습은 검증 절차의 감각을 길러준다.

14. 참고 문헌

Karnopp, D. (1985). Computer simulation of stick-slip friction in mechanical dynamic systems. Journal of Dynamic Systems, Measurement, and Control, 107(1), 100–103.
Moreau, J. J. (1988). Unilateral contact and dry friction in finite freedom dynamics. In Nonsmooth Mechanics and Applications, Springer, 1–82.
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.
Hairer, E., Lubich, C., & Wanner, G. (2006). Geometric Numerical Integration (2nd ed.). Springer.
Brogliato, B. (2016). Nonsmooth Mechanics: Models, Dynamics and Control (3rd ed.). Springer.
Acary, V., & Brogliato, B. (2008). Numerical Methods for Nonsmooth Dynamical Systems. Springer.
Olsson, H., Åström, K. J., Canudas de Wit, C., Gäfvert, M., & Lischinsky, P. (1998). Friction models and friction compensation. European Journal of Control, 4(3), 176–195.

version: 1.0