The semi-explicit nonsmooth Newmark time integrator for robust unilateral contact in dynamic fragmentation simulations

Este artigo introduz e valida um esquema de integração temporal Newmark semissimplícito (NSN) que lida de forma robusta com o contato unilateral em simulações de fragmentação dinâmica ao impor fortemente as restrições, alcançando assim uma estabilidade e precisão superiores aos métodos baseados em penalidade, ao mesmo tempo em que revela que a dissipação de contato pode paradoxalmente aumentar a contagem de fragmentos ao melhorar a localização do dano.

O essencial

Imagine que você está tentando simular o que acontece quando um objeto sólido, como um prato de cerâmica ou uma rocha, é atingido com tanta força que se estilhaça em milhares de pequenos pedaços. Isso não é apenas uma quebra simples; é uma explosão caótica onde os pedaços voam, colidem uns com os outros, batem nas paredes e se esmagam.

O artigo apresenta um novo programa de computador (um "integrador de tempo") projetado para simular esse caos sem que o computador trave ou apresente resultados sem sentido. Veja como ele funciona, dividido em conceitos simples:

1. O Problema: A Armadilha da "Mola"

Para simular a quebra, os cientistas geralmente usam um método onde fingem que o material é feito de pequenas molas. Quando o material quebra, as molas se rompem. Quando os pedaços colidem, eles usam "molas de penalidade" para empurrá-los para longe, para que não passem uns pelos outros.

A Analogia: Imagine tentar parar uma bola de boliche com um elástico.

• Se o elástico for muito frouxo (baixa rigidez), a bola passa direto por ele (não físico/irrealista).
• Se o elástico for super apertado (alta rigidez) para parar a bola perfeitamente, ele age como uma parede rígida. Mas, se você o tornar muito apertado, o computador terá que dar passos minúsculos, minúsculos para calcular o rebote, fazendo a simulação demorar uma eternidade.
• A Alegação do Artigo: O método antigo (usando essas molas apertadas) é instável. Ele faz o computador derivar, perder energia ou travar, especialmente quando há milhões de colisões acontecendo ao mesmo tempo.

2. A Solução: O "Guarda de Trânsito" (Nonsmooth Newmark)

Os autores criaram um novo método chamado Nonsmooth Newmark (NSN). Em vez de usar elásticos para afastar os pedaços, este método age como um guarda de trânsito rigoroso em um cruzamento movimentado.

A Analogia:

• O Volume (O Carro): O corpo principal do objeto move-se livre e suavemente. O computador prevê para onde o carro iria se não houvesse obstáculos. Esta parte é calculada de forma muito rápida (explicitamente).
• O Contato (O Cruzamento): Se o carro atingir uma parede ou outro carro, o "guarda de trânsito" intervém. Em vez de empurrar o carro de volta com uma mola, o guarda diz instantaneamente: "Pare! Você não pode ir para lá". Ele impõe uma regra rígida: Não passar através de outro objeto.
• A Magia: Este método trata a regra de "não passar" como uma lei dura da física, em vez de uma mola suave. Isso permite que o computador dê passos de tempo muito maiores porque não precisa se preocupar com o elástico ficando apertado demais.

3. A Abordagem de "Personalidade Dividida"

O artigo descreve este método como "semi-explícito". Pense nisso como uma dança de dois passos:

1. Passo A (A Previsão): O computador adivinha onde tudo estará no próximo momento, ignorando as colisões.
2. Passo B (A Correção): Se a previsão mostrar dois pedaços sobrepostos, o computador corrige instantaneamente a velocidade e a posição deles para corrigir a sobreposição, exatamente como uma bola de bilhar atingindo outra e mudando de direção instantaneamente.

Isso permite que a simulação seja rápida (como a previsão) mas estável e precisa (como a correção).

4. O Que Eles Descobriram (Os Experimentos)

Os autores testaram este novo método de "guarda de trânsito" contra os antigos métodos de "elástico" usando três cenários:

• A Bola Saltitante: Uma bola simples saltitando no chão. O novo método foi tão preciso quanto os melhores métodos existentes, mas lidou com os saltos sem perder energia ou apresentar instabilidades (jitter).
• A Barra de Impacto: Uma barra de metal atingindo uma parede. Os métodos antigos tiveram dificuldades com a velocidade do impacto, mas o novo método lidou com o "esmagamento" perfeitamente, mantendo os cálculos de energia corretos.
• A Barra Estilhaçada: Uma barra que já possui rachaduras e então atinge uma parede. Os métodos antigos exigiam passos de tempo tão minúsculos para manter a estabilidade que eram incrivelmente lentos. O novo método pôde dar passos enormes, rodando 27 vezes mais rápido e sendo mais preciso.

5. A Descoberta Surpreendente: Estilhaçamento Confinado

A parte mais interessante do artigo envolve um experimento "confinado". Imagine uma barra se estilhaçando dentro de uma caixa pequena em vez de em um espaço aberto.

• A Intuição Antiga: Você poderia pensar que, se os pedaços quicarem nas paredes e perderem energia (dissipação), haveria menos energia restante para quebrar o material, resultando em menos pedaços e pedaços maiores.
• O Achado do Artigo: O oposto aconteceu. Quando os pedaços quicaram nas paredes e perderam um pouco de energia (dissipação de contato), o material na verdade se quebrou em mais pedaços e menores.
• Por quê? Os autores explicam que o "quicar" atua como um filtro. Em um mundo perfeitamente elástico (com muito rebote), as ondas de tensão ricocheteiam descontroladamente, fazendo o material ficar "confuso" e desenvolver muitas rachaduras minúsculas e fracas que não se separam totalmente. Quando as paredes absorvem parte dessa energia, as ondas se acalmam. Isso permite que a tensão se concentre em pontos específicos, conduzindo as rachaduras por todo o caminho para criar fragmentos limpos e separados.

Resumo

O artigo apresenta uma nova ferramenta matemática que simula a quebra de objetos tratando as colisões como regras duras e instantâneas, em vez de molas suaves. Isso torna a simulação por computador:

1. Mais Estável: Não trava nem sofre desvios.
2. Mais Rápida: Pode dar passos de tempo maiores.
3. Mais Precisa: Prevê corretamente em quantos pedaços um objeto se quebrará.

Os autores concluem que esta ferramenta está pronta para ser usada em simulações 3D complexas, como entender como detritos espaciais se despedaçam ou como rochas se estilhaçam em avalanches, ao fornecer uma maneira robusta de lidar com a dança caótica de milhões de fragmentos colidindo.

Resumo técnico

Resumo Técnico: O Integrador de Tempo Newmark Semie-Explícito para Contato Unilateral Robusto em Simulações de Fragmentação Dinâmica

Definição do Problema
Simulações numéricas de fragmentação dinâmica envolvem um desafio duplo: capturar a rápida evolução topológica de um contínuo sólido fragmentando-se em fragmentos discretos e resolver os contatos densos e simultâneos que ocorrem entre esses fragmentos e as faces das trincas. Embora o Método de Elementos Finitos (FEM) com modelos de zona coesiva (CZM) extrínsecos trate a fratura de forma eficaz, os tratamentos de contato convencionais apresentam dificuldades significativas. Métodos baseados em penalidade, frequentemente usados pela sua compatibilidade com integração de tempo explícita, sofrem de um compromisso fundamental: a baixa rigidez permite interpenetração não física, enquanto a alta rigidez restringe o passo de tempo estável a níveis computacionalmente proibitivos. Além disso, a interação entre a fratura e o contato regularizado introduz não-suavidade (nonsmoothness), o que leva a uma deriva de energia numérica persistente e instabilidade, mascarando a física subjacente.

Metodologia
Os autores propõem e validam um novo esquema de integração de tempo semie-explícito: o método Newmark Não-Suave (NSN). Esta abordagem adapta o framework de Dinâmica de Contato Não-Suave (NSCD) para a fragmentação dinâmica, desacoplando o tratamento da dinâmica de volume e das interações de contato:

1. Formulação Semie-Explícita: O esquema divide a dinâmica em uma parte "suave" (deformação de volume, carregamento externo e fratura coesiva) e uma parte "não-suave" (reações de contato).
   • A dinâmica suave é integrada explicitamente usando um esquema Newmark-$\beta$ de segunda ordem ($\beta=0, \gamma=1/2$), mantendo alta precisão durante as fases de voo livre.
   • A dinâmica não-suave (contato unilateral) é tratada implicitamente ao nível da velocidade usando uma abordagem de Euler implícito de primeira ordem. Isso impõe fortemente as restrições de lacuna de Signorini sem complacência artificial.
2. Implementação Algorítmica: O problema de contato é formulado como um Problema de Complementaridade Linear (LCP), que é resolvido como um problema de Programação Quadrática (QP) convexa restrito ao conjunto de contato ativo. Isso permite soluções eficientes e paralelizáveis que escalam com o número de contatos ativos, em vez do tamanho total do sistema.
3. Lei de Tração-Separação (TSL) Modificada: Para abordar o desafio numérico de rigidez coesiva infinita conforme o dano se aproxima de zero em esquemas explícitos, os autores introduzem um limite de rigidez ($\tilde{k}$). Esta regularização evita que o passo de tempo se torne nulo, preservando a convexidade do problema de contato.
4. Validação e Aplicação: O framework é implementado na biblioteca FEM de código aberto Akantu. Ele é validado contra soluções analíticas e métodos não-suaves estabelecidos (Moreau-Jean e CD-Lagrange) usando benchmarks como a bola saltitante e a barra de impacto. É subsequentemente aplicado à fragmentação dinâmica 1D sob expansão livre e confinada.

Principais Contribuições e Resultados

• Precisão e Estabilidade: Testes de benchmark demonstram que o esquema NSN atinge precisão comparável a métodos não-suaves estabelecidos (Moreau-Jean e CD-Lagrange) e supera significativamente as abordagens baseadas em penalidade. No teste da barra de impacto, o NSN mantém precisão de segunda ordem para deslocamento e precisão superior de velocidade em comparação ao CD-Lagrange, particularmente em regimes elásticos.
• Eficiência Computacional: Embora o método NSN incorra em um custo computacional por passo de tempo mais elevado (aproximadamente 5 vezes maior que um esquema de penalidade puramente explícito devido à resolução do QP), sua estabilidade robusta permite passos de tempo significativamente maiores. Para benchmarks 1D, isso resulta em um tempo total de solução comparável ou melhor que abordagens puramente explícitas, com erros de conservação de energia reduzidos em ordens de magnitude (próximo à precisão de máquina).
• Insights Físicos na Fragmentação Confinada:
  • Mudança no Orçamento de Energia: Em cenários de expansão confinada (onde os fragmentos colidem com paredes rígidas), o orçamento de energia para a fratura desloca-se da energia cinética local dos fragmentos individuais para a energia cinética global de todo o sistema. Isso leva a uma redução drástica no tamanho médio dos fragmentos e a um aumento na dissipação de energia de fratura total em comparação com a expansão livre.
  • Papel Contraintuitivo da Dissipação de Contato: O estudo revela que a introdução de dissipação de contato (colisões inelásticas, $e < 1$) reduz a energia total dissipada pela fratura, mas aumenta a contagem final de fragmentos. Isso ocorre porque a dissipação de contato atua como um filtro passa-baixas, suprimindo oscilações de ondas de tensão de alta frequência. Este filtragem promove uma propagação de ondas de tensão mais limpa e uma melhor localização de dano, conduzindo as trincas à separação total em vez de permitir que o dano permaneça distribuído em zonas parcialmente danificadas ($d < 1$).

Significância
O artigo estabelece o esquema NSN como uma ferramenta robusta para gerar estatísticas de fragmentação de alta fidelidade em cenários de contato denso. Ao tratar rigorosamente o contato como não-suave enquanto retém a eficiência explícita para a dinâmica de volume, o método supera as instabilidades numéricas inerentes às formulações de penalidade. Os autores afirmam que este framework é essencial para simular interações complexas de múltiplos corpos, como as encontradas em cenários de detritos orbitais, onde interações frequentes entre fragmentos impulsionam a fragmentação secundária. O trabalho representa um passo crítico para estender a análise de fragmentação de alta fidelidade para dimensões superiores, aproveitando a capacidade do método de preservar a fidelidade física sem o custo proibitivo de resoluções implícitas globais.