A Splitting Architecture for Exact Reduced Coulomb Friction

Technische Samenvatting: Een Splitting-architectuur voor Exacte Gereduceerde Coulomb-wrijving

Probleemstelling
Het artikel behandelt de numerieke uitdaging van het simuleren van wrijvingscontactdynamica voor starre lichamen met behulp van de exacte gereduceerde Coulomb-wrijvingswet. Hoewel de exacte wet (geformuleerd door De Saxcé en Feng [DSF98]) macroscopische gedragingen zoals stabiele stapeling, boogondersteuning en zuivere stick-slip-overgangen accuraat weergeeft, brengt het aanzienlijke computationele moeilijkheden met zich mee. Specifiek introduceert de wet een "geaugmenteerde snelheid" die de tangentiële glijsnelheid koppelt aan de normale richting. Deze koppeling maakt de wrijvingswet niet-geassocieerd, wat betekent dat het probleem niet kan worden geformuleerd als een enkel energie-minimalisatieprobleem of een convex optimalisatieprobleem.

Bestaande benaderingen gaan hier doorga af aan door:

1. Smoothing (Gladstrijken): Het vervangen van de niet-gladde wet door een geregulariseerde kracht, wat een afruil tussen stijfheid en nauwkeurigheid introduceert.
2. Convex Relaxatie: Het weglaten van de niet-geassocieerde koppelingsterm om een Kegel Complementariteits Probleem (CCP) op te lossen. Dit is efficiënt, maar produceert zichtbare artefacten, zoals contacten die loslaten bij de kegelgrens tijdens aanhoudende glijbewegingen.
3. Monolithische Nonsmooth Newton: Het direct oplossen van het volledige niet-gladde systeem. Hoewel deze methoden exact zijn, zijn ze moeilijk te globaliseren, te parallelliseren of te modulariseren.

De auteurs stellen dat de exacte gereduceerde wet een natuurlijke structurele splitsing bezit die onderbenut is in computer graphics solvers: een kegel-beperkte lineaire respons gekoppeld aan een scalaire, niet-geassocieerde term.

Methodologie
De auteurs stellen een splitting-architectuur voor die de exacte gereduceerde Coulomb-wet behandelt door deze te ontleden in twee afzonderlijke operatoren:

1. Operator A (Kegel-beperkte Lineaire Respons): Een sterk convexe Second-Order Cone Program (SOCP) die de lineaire contactrespons representeert, beperkt door de wrijvingskegel.
2. Operator B (Scalaire Koppeling): De niet-geassocieerde De Saxcé en Feng koppelingsterm, die per contact goedkoop te evalueren is maar globaal niet-monotoon is.

De solver maakt gebruik van een Tseng-stijl Forward-Backward-Forward (FBF) outer iteratie-schema:

• Outer Loop: Evalueert expliciet de niet-geassocieerde koppeling $B(\lambda)$ bij de huidige iteratie.
• Inner Solve: Lost een sterk convex kegel-beperkt QP (Operator A) op waarbij de koppelingsterm bevroren is. Dit innerlijke probleem is goed geconditioneerd, ongeacht de rang van de Delassus-operator, dankzij proximale regularisatie.
• Correctie: Past een lichte expliciete correctie toe om rekening te houden met de verandering in de koppelingsterm tussen de huidige en de tussenliggende iteraties.

Belangrijkste Bijdragen

1. Splitting Solver Architectuur: Het artikel presenteert een solver die de kegel-QP/scalaire-koppeling decompositie adopteert zoals geïdentificeerd door Acary et al. [ACLM11], maar vervangt hun Picard outer iteratie door een Tseng-stijl FBF-schema. Deze structuur maakt de innerlijke kegel-solve modulair, waardoor verschillende convexe SOCP-solvers (bijv. block Gauss-Seidel, ADMM, interior-point) kunnen worden gesubstitueerd zonder de outer loop te wijzigen.
2. Safeguarded Adaptive Step-Size Rule: Omdat de De Saxcé-koppeling niet-monotoon is, zijn standaard convergentie-theorema's voor de methode van Tseng niet direct van toepassing. De auteurs introduceren een stapgrootte-controlemechanisme dat een Lipschitz-gebaseerde schatting combineert met een lokale variatie-acceptatietest en een asymmetrisch groeibeleid, afgestemd op de stick-slip grens.
3. Matrix-vrije, Contact-parallele Implementatie: Geïmplementeerd bovenop NVIDIA Warp en de Newton physics engine, vermijdt de solver het vormen van dichte matrices. Er wordt gebruikgemaakt van matrix-vrije operatoren voor de Delassus-matrix en de evaluatie van koppeling en correctie wordt over contacten geparallelliseerd.

Resultaten
De methode is geëvalueerd op rigid-body benchmarks die ontworpen zijn om statische ondersteuning, wrijvingsoverdracht en modus-transities te testen:

• Kubus op een helling & Bal met backspin: De solver reproduceert de exacte analytische overgang van glijden naar plakken (sticking) en handelt de generatie van normale krachten door tangentiële backspin correct af. In tegen plaats vertoonde de convex-gerelaxeerde MuJoCo-solver spuurse normale separatie (drifting) en slaagde er niet in contact te behouden.
• Draaitafel & Painlevé Box: De methode vangt de overgang tussen capture en ejectie onder centripetale belasting correct op, evenals de slide-to-tumble transitie, waarbij de resultaten overeenkomen met het gedrag van de industriële Kamino-solver.
• Kaartenhuis & Metselwerkboog: In sterk gekoppelde statische scenario's behield de splitting-solver de structurele integriteit en lokale wrijvingsondersteuning. MuJoCo slaagde er niet in de staande configuratie te handhaven (het uiteen drijven), terwijl de voorgestelde methode het evenwicht bewaarde. Opvallend genoeg slaagde de solver erin een boog van 101 stenen te balanceren waar Kamino onder dezelfde setup faalde.
• Prestaties: Op contact-rijke benchmarks was de methode concurrerend met Kamino en een orde van grootte sneller dan MuJoCo. De innerlijke kegel-solve schaalt lineair met het aantal contacten en de nabijheid van lichamen (body adjacency).

Betekenis en Claims
Het artikel claimt dat deze splitting-architectuur succesvol de exacte gereduceerde Coulomb-wet reproduceert zonder de introductie van smoothing-parameters of convex relaxaties. Door de niet-geassocieerde koppeling te scheiden van de kegelgeometrie, bereikt de methode:

• Fysieke Nauwkeurigheid: Het elimineert de "normal-lifting" artefacten die veel voorkomen bij convex-gerelaxeerde solvers, waardoor wordt gegarandeerd dat stabiele stapels en bogen stabiel blijven.
• Modulariteit: De innerlijke solver is ontkoppeld van de outer iteratie, waardoor hardware-specifieke of probleemgrootte-specifieke solvers efficiënt kunnen worden gebruikt.
• Convergentie: Hoewel de niet-monotone aard van de koppeling een formele globale convergentiegarantie verhindert, demonstreert de methode empirische convergentie naar de single-precision floor op alle geteste benchmarks.

De auteurs erkennen dat de expliciete behandeling van de niet-geassocieerde koppeling contractiviteit kan verliezen in globaal gekoppelde systemen (zoals grote bogen), wat meer outer iteraties vereist. Ze stellen echter dat de methode robuust en accuraat blijft waar andere exacte solvers moeite hebben, wat een levensvatbaar pad biedt voor hoogwaardige wrijvingssimulatie in computer graphics.