MPM을 위한 정밀 마찰 접촉 해결 방법

초록

본 논문은 암시적 물질점법(MPM) 시뮬레이션에서 접촉 검출과 마찰 접촉 해석을 일관된 파이프라인으로 통합한다. 입자 중심의 기하학 프리미티브를 이용해 서브셀 수준에서 접촉점을 정확히 찾고, 접촉 충동을 비선형 보완 문제(NCP) 형태로 정의한 뒤 ADMM 알고리즘으로 해결한다. 기존의 그리드 기반 검출·벌점 방식과 달리 전역적인 접촉 행렬을 재사용해 효율성을 높이고, 재료 모델·보간·전달 방식에 독립적인 일반성을 제공한다. 7가지 실험을 통해 탄성·소성, 복합 형상, 다양한 마찰 조건에서 높은 정확도와 안정성을 입증한다.

상세 분석

이 연구는 MPM의 핵심 약점인 마찰 접촉 처리 문제를 두 단계로 명확히 구분하고, 각각에 최적화된 기법을 적용한다. 첫 번째 단계인 충돌 검출에서는 입자마다 변형된 테트라헤드론을 프리미티브로 사용하고, GJK 알고리즘 변형을 통해 입자‑입자 혹은 입자‑메시 간의 근접성을 서브셀 정밀도로 판단한다. 이때 얻어지는 접촉점, 법선, 접촉 프레임은 입자 변형에 직접 연동되므로, 격자 노드 기반 검출에서 발생하는 확산 현상을 원천적으로 제거한다. 두 번째 단계인 접촉 해석에서는 접촉 충동 λ를 변수로 하는 비선형 보완 문제(NCP)를 수립한다. 여기서 Signorini 조건과 쿠롱 마찰 법칙, 최대 소산 원리를 결합해 2차원 원뿔(K_μ)과 그 쌍대 원뿔(K*_μ)으로 표현한다. 핵심은 MPM의 암시적 시간 적분 과정에서 이미 구성된 어드미턴스 행렬 A⁻¹을 재활용해 Delassus 연산자 G = H A⁻¹ Hᵀ 를 효율적으로 얻는 점이다. 이렇게 얻은 G와 자유 접촉 속도 g를 이용해 ADMM 반복을 수행한다. ADMM은 선형 단계에서 기존의 전역 행렬을 활용해 빠르게 속도 보정을 계산하고, 비선형 단계에서는 접촉 프레임별로 마찰 원뿔에 투영하는 근접 연산만 수행한다는 장점이 있다. 알고리즘은 원시적인 프리미티브와 전역 행렬을 공유함으로써 메모리 사용량을 최소화하고, 각 타임스텝마다 행렬을 한 번만 팩터링하면 된다는 점에서 기존의 순차적 혹은 지역적 벌점 방식보다 월등히 효율적이다. 또한, 입자‑그리드‑접촉 프레임 간의 전이 과정을 바리센트리 가중치를 이용해 4개의 정점 프록시로 분산시킴으로써, 접촉이 발생한 셀 내부에서도 활성 노드에만 힘이 전달되도록 보장한다. 실험 섹션에서는 선형, St. Venant‑Kirchhoff, 코로테이셔널, Neo‑Hookean 등 다양한 재료 모델과 APIC, MLS, CK‑MPM 등 여러 전송·보간 스킴을 조합해도 동일한 정확도와 수렴성을 유지함을 보여준다. 특히, 다물체 시뮬레이션에서 각 물체마다 독립적인 그리드를 유지하면서도 전역 접촉 행렬을 하나로 결합하는 설계는 확장성을 크게 향상시킨다. 전체적으로 이 논문은 MPM에 특화된 마찰 접촉 해결 체계를 제시함으로써, 로봇 제어·학습 등 실시간 혹은 차분 가능한 시뮬레이션에 필요한 정확도와 안정성을 동시에 달성한다는 점에서 큰 의의를 가진다.