A Splitting Architecture for Exact Reduced Coulomb Friction

Resumo Técnico: Uma Arquitetura de Divisão para o Atrito de Coulomb Reduzido Exato

Definição do Problema
O artigo aborda o desafio numérico de simular a dinâmica de contato friccional para corpos rígidos utilizando a lei de atrito de Coulomb reduzida exata. Embora a lei exata (formulada por De Saxcé e Feng [DSF98]) capture com precisão comportamentos macroscópicos como empilhamento estável, suporte de arcos e transições limpas de stick-slip, ela apresenta dificuldades computacionais significativas. Especificamente, a lei introduz uma "velocidade aumentada" que acopla a velocidade de deslizamento tangencial na direção normal. Esse acoplamento torna a lei de atrito não associada, o que significa que o problema não pode ser formulado como um único problema de minimização de energia ou de otimização convexa.

Abordagens existentes normalmente lidam com isso através de:

1. Suavização (Smoothing): Substitui a lei não suave por uma força regularizada, introduzindo um compromisso entre rigidez e precisão (stiffness-accuracy trade-off).
2. Relaxação Convexa: Descarta o termo de acoplamento não associado para resolver um Problema de Complementaridade de Cone (CCP). Isso é eficiente, mas produz artefatos visíveis, como o afastamento de contatos na fronteira do cone durante o deslizamento sustentado.
3. Newton Não Suave Monolítico: Resolve o sistema não suave completo diretamente. Embora exato, esses métodos são difíceis de globalizar, paralelizar ou modularizar.

Os autores argumentam que a lei reduzida exata possui uma divisão estrutural natural que tem sido subutilizada em solvers de computação gráfica: uma resposta linear restrita por cone acoplada a um termo escalar não associado.

Metodologia
Os autores propõem uma arquitetura de divisão (splitting architecture) que trata a lei de Coulomb reduzida exata decompondo-a em dois operadores distintos:

1. Operador A (Resposta Linear Restrita por Cone): Um Programa de Otimização de Segundo Ordem (SOCP) fortemente convexo que representa a resposta de contato linear restrita pelo cone de atrito.
2. Operador B (Acoplamento Escalar): O termo de acoplamento de De Saxcé e Feng, que é barato de avaliar por contato, mas globalmente não monotônico.

O solver emprega um esquema de iteração externa Forward-Backward-Forward (FBF) do tipo Tseng:

• Loop Externo: Avalia explicitamente o acoplamento não associado $B(\lambda)$ na iteração atual.
• Resolução Interna (Inner Solve): Resolve um QP (Programação Quadrática) de cone fortemente convexo (Operador A) com o termo de acoplamento congelado. Este problema interno é bem condicionado, independentemente do posto da matriz de Delassus, devido à regularização proximal.
• Correção: Aplica uma correção explícita leve para contabilizar a mudança no termo de acoplamento entre a iteração atual e a intermediária.

Principais Contribuições

1. Arquitetura de Solver de Divisão: O artigo apresenta um solver que adota a decomposição cone-QP/acoplamento-escalar identificada por Acary et al. [ACLM11], mas substitui a iteração externa de Picard pelo esquema FBF do tipo Tseng. Essa estrutura torna a resolução do cone interna modular, permitindo que diferentes solvers de SOCP convexos (ex: Gauss-Seidel de bloco, ADMM, método de pontos interiores) sejam substituídos sem alterar o loop externo.
2. Regra de Passo Adaptativo Protegida: Como o acoplamento de De Saxcé é não monotônico, os teoremas de convergência padrão para o método de Tseng não se aplicam diretamente. Os autores introduzem um mecanismo de controle de passo combinando uma estimativa baseada em Lipschitz com um teste de aceitação de variação local e uma política de crescimento assimétrica ajustada à fronteira de stick-slip.
3. Implementação Livre de Matriz e Paralela por Contato: Implementado sobre o NVIDIA Warp e o motor de física Newton, o solver evita a formação de matrizes densas. Ele utiliza operadores livres de matriz para a matriz de Delassus e paraleliza a avaliação e a correção de acoplamento sobre os contatos.

Resultados
O método foi avaliado em testes de referência de corpo rígido projetados para estressar suporte estático, transmissão friccional e transições de modo:

• Cubo em Inclinação e Bola com Backspin: O solver reproduz a transição analítica exata de deslizamento para aderência (sliding to sticking) e lida corretamente com a geração de força normal proveniente do backspin tangencial. Em contraste, o solver de MuJoCo (com relaxação convexa) exibiu separação normal espúria (deriva) e falhou em manter o contato.
• Mesa Giratória e Caixa de Painlevé: O método captura corretamente a transição entre captura e ejeção sob carga centrípeta e a transição de deslizamento para tombamento, correspondendo ao comportamento do solver de alta performance Kamino.
• Casa de Cartas e Arco de Alvenaria: Em cenários estáticos altamente acoplados, o solver de divisão manteve a integridade estrutural e o suporte friccional local. O MuJoCo falhou em manter a configuração de pé (afastando-se), enquanto o método proposto preservou o equilíbrio. Notavelmente, o solver equilibrou com sucesso um arco de 101 pedras onde o Kamino falhou sob a mesma configuração.
• Desempenho: Em testes de referência ricos em contatos, o método foi competitivo com o Kamino e uma ordem de magnitude mais rápido que o MuJoCo. A resolução do cone interna escala linearmente com a contagem de contatos e a adjacência de corpos.

Significância e Alegações
O artigo afirma que esta arquitetura de divisão reproduz com sucesso a lei de Coulomb reduzida exata sem introduzir parâmetros de suavização ou relaxações convexas. Ao separar o acoplamento não associado da geometria do cone, o método alcança:

• Precisão Física: Elimina os artefatos de "elevação normal" comuns em solvers de relaxação convexa, garantindo que pilhas e arcos estáveis permaneçam estáveis.
• Modularidade: O solver interno é desacoplado da iteração externa, permitindo que solvers específicos de hardware ou de tamanho de problema sejam usados eficientemente.
• Convergência: Embora a natureza não monotônica do acoplamento impeça uma garantia formal de taxa de convergência global, o método demonstra convergência empírica até o limite de precisão simples (single-precision floor) em todos os testes realizados.

Os autores reconhecem que o tratamento explícito do acoplamento não associado pode perder contratividade em sistemas globalmente acoplados (como grandes arcos), exigindo mais iterações externas. No entanto, eles afirmam que o método permanece robusto e preciso onde outros solvers exatos têm dificuldade, oferecendo um caminho viável para simulação de atrito de alta fidelidade em computação gráfica.