[논문 리뷰] Rapidly computable viscous friction and no-slip rigid contact models

Evan Drumwright|arXiv (Cornell University)|2015. 04. 03.

Dynamics and Control of Mechanical Systems참고 문헌 19인용 수 3

한 줄 요약

이 논문은 쿠론 마찰 모델의 높은 복잡도를 피하기 위해 점성 마찰 또는 슬립이 없는 조건에서의 강체 접촉을 계산적으로 효율적으로 시뮬레이션하는 알고리즘을 제시한다. 문제를 양의 준정적 선형 보완 문제(LCP)로 재구성함으로써 O(m³ + m²n)의 예상 시간 복잡도를 달성하여 해법 실패를 줄이고 로봇의 보행 및 집게 작업에서 더 빠른 시뮬레이션을 가능하게 한다.

ABSTRACT

This article presents computationally efficient algorithms for modeling two special cases of rigid contact---contact with only viscous friction and contact without slip---that have particularly useful applications in robotic locomotion and grasping. Modeling rigid contact with Coulomb friction generally exhibits $O(n^3)$ expected time complexity in the number of contact points and $2^{O(n)}$ worst-case complexity. The special cases we consider exhibit $O(m^3 + m^2n)$ time complexity ($m$ is the number of independent coordinates in the multi rigid body system) in the expected case and polynomial complexity in the worst case; thus, asymptotic complexity is no longer driven by number of contact points (which is conceivably limitless) but instead is more dependent on the number of bodies in the system (which is often fixed). These special cases also require considerably fewer constrained nonlinear optimization variables thus yielding substantial improvements in running time. Finally, these special cases also afford one other advantage: the nonlinear optimization problems are numerically easier to solve.

연구 동기 및 목표

• 쿠론 마찰을 고려한 강체 접촉 시뮬레이션의 높은 계산 비용을 해결하되, 이는 O(n³)으로 증가하며 최악의 경우 2^O(n)에 이를 수 있다.
• 순수한 점성 마찰과 슬립이 없는 강체 접촉이라는 두 가지 특수한 접촉 케이스에 대해 계산적으로 효율적인 대안을 개발한다.
• 쿠론 마찰 모델을 사용하지 않을 경우 비선형 최적화 변수의 수를 줄이고 접촉 해법의 수치적 안정성을 향상시킨다.
• 복잡한 쿠론 마찰 모델을 대체하여 로봇의 보행 및 집게 작업에서 더 빠르고 신뢰할 수 있는 시뮬레이션을 가능하게 한다.

제안 방법

• 접촉 문제를 혼합 선형 보완 문제(MLCP)로 공식화한 후, 행렬 역행렬을 이용해 비제약 변수를 제거함으로써 양의 준정적 선형 보완 문제(LCP)로 간소화한다.
• 결과적으로 얻어진 LCP를 해결하기 위해 수정된 사영 원단일-이중(PPM) 알고리즘을 사용하며, 이는 표준 쿠론 마찰 솔버에서 사용하는 공액성-더블(LCP)보다 수치적으로 더 안정적이다.
• LCP 공식화를 두 가지 특수 케이스에 적용한다: 슬립이 없는 접촉(상대 접선 속도가 0)과 점성 마찰(마찰력이 상대 속도에 비례함).
• 시간 복잡도를 O(m³ + m²n)으로 유도하며, 여기서 m은 독립 좌표의 수이고 n은 접촉점의 수이다. 이는 복잡도를 접촉점의 수에서 시스템 내 물체의 수로 분리한다.
• 심플렉틱 오일러 및 룬게-쿠타 적분을 사용하여 Moby 시뮬레이터에 알고리즘을 구현하였으며, 엄격한 침투 방지를 위한 이벤트 기반 시뮬레이션을 수행한다.
• 효율적인 행렬 연산을 활용하여 PPM 솔버를 최적화하였으며, MEX 파일 컴파일을 피함에도 불구하고 상당한 속도 향상을 달성하였다.

실험 결과

연구 질문

• RQ1 슬립이 없는 상태 또는 점성 마찰을 고려한 강체 접촉 모델은 다항 시간 복잡도로 해결될 수 있는가? 쿠론 마찰의 지수적 최악의 경우 복잡도를 피할 수 있는가?
• RQ2 쿠론 마찰을 사용하지 않는 접촉 모델에서 비선형 최적화 변수의 수를 상당히 줄일 수 있는가?
• RQ3 공액성-더블 LCP를 양의 준정적 LCP로 대체하면 수치적 안정성과 솔버 수렴성이 향상되는가?
• RQ4 제안된 방법은 물리적 정확도를 희생시키지 않고도 로봇의 보행 및 집게 시뮬레이션에서 상당한 런타임 향상을 달성할 수 있는가?

주요 결과

• 제안된 알고리즘은 O(m³ + m²n)의 예상 시간 복잡도를 달성하여 계산 부담을 접촉점의 수에서 물체의 수로 이동시켰다.
• 슬립이 없는 접촉 모델의 시뮬레이션은 기준값(416.284초 대비 301.796초)보다 28% 빠르게 실행되었으며, 평균 접촉점 수는 6개, 표준편차는 8이었다.
• 점성 마찰 모델은 37%의 속도 향상을 기록했으며(2936.36초 대비 2139.73초), 그러나 ODE의 민감성으로 인한 급격한 운동으로 인해 더 긴 런타임이 관찰되었다.
• MEX 최적화 없이도 수정된 PPM 알고리즘이 quadprog보다 약 20% 빠른 성능을 보였으며, 평균 피벗 연산 수는 5.5개였다.
• 양의 준정적 LCP 공식화는 표준 쿠론 마찰 솔버보다 해법 실패를 줄이고 수치적 안정성을 향상시켰다.
• 이 방법은 침투가 없는 검증 가능한 시뮬레이션을 가능하게 하여 고정밀도 로봇 역학 시뮬레이션에 적합하며 성능을 희생시키지 않는다.