Dynamic brittle fracture using Lip-field approach in an explicit dynamics context

Resumo Técnico: Fratura Frágil Dinâmica Usando a Abordagem de Campo Lip em um Contexto de Dinâmica Explícita

Definição do Problema
Este artigo aborda a simulação numérica de fratura frágil dinâmica em materiais sujeitos a carregamento de alta taxa de deformação, como impactos ou explosões. Enquanto a mecânica da fratura clássica (por exemplo, o modelo de Griffith) apresenta dificuldades com a iniciação e ramificação de trincas, abordagens modernas como Modelos de Zona Coesiva (CZM) e modelos de dano difuso foram desenvolvidas. No entanto, modelos de dano difuso locais sofrem de sensibilidade patológica à malha, exigindo regularização para introduzir uma escala de comprimento e garantir resultados independentes da malha. As técnicas de regularização existentes, como os métodos de Campo de Fase (Phase-Field) e de Nível de Conjunto Espesso (Thick Level Set - TLS), frequentemente envolvem funcionais de energia não locais que podem ser computacionalmente dispendiosos. O desafio específico abordado aqui é a extensão da abordagem de campo Lip — um método variacional que introduz escalas de comprimento via restrições de Lipschitz no campo de dano — de suas aplicações anteriores em regime quase-estático e dinâmico unidimensional para a fratura dinâmica bidimensional dentro de um arcabouço de integração temporal explícita.

Metodologia
Os autores propõem uma formulação variacional para a mecânica do dano dinâmica baseada na minimização de um potencial incremental. O núcleo da metodologia envolve:

1. Formulação Variacional: O problema mecânico é definido pela minimização de um potencial incremental $F$ que inclui energia cinética, energia livre e potencial de dissipação. As variáveis de estado são o deslocamento ($u$) e uma variável escalar de dano ($d$). O modelo constitutivo utiliza materiais de padrão generalizado com funções específicas de degradação de energia ($g(d)$) e amolecimento ($h(d)$).
2. Regularização por Campo Lip: Ao contrário dos modelos de Campo de Fase que modificam o funcional de energia com um termo de gradiente, a abordagem de campo Lip mantém o funcional de energia local (não regularizado). Em vez disso, impõe restrições de Lipschitz ao campo de dano, exigendo que o gradiente de dano seja limitado por um valor crítico ($1/l$). Isso restringe o campo de dano a um espaço Lipschitz específico $\mathcal{L}$, introduzindo efetivamente uma escala de comprimento $l$ sem alterar a densidade de energia local.
3. Esquema Estagiado Explícito: O artigo emprega um esquema de integração temporal de Newmark explícito (método de diferença central). Um esquema estagiado é utilizado onde os campos de deslocamento e dano são resolvidos sequencialmente. Para um dado passo de tempo, o deslocamento é atualizado primeiro, seguido pela atualização do dano.
4. Solucionador de Dano Eficiente: Uma contribuição algorítmica fundamental é o uso de estimativa de limites (derivada de trabalhos anteriores) para otimizar o cálculo do dano. O processo envolve:
   • Resolver um problema de minimização puramente local para encontrar um campo de dano local ($d_{loc}$).
   • Projetar $d_{loc}$ em limites superiores e inferiores baseados na restrição de Lipschitz.
   • Identificar regiões onde esses limites diferem.
   • Resolver o problema de minimização não local (com restrição de Lipschitz) apenas na pequena porção do domínio onde os limites não são iguais. Em regiões onde os limites coincidem, a solução local é aceita. Isso reduz significativamente o custo computacional em comparação com a resolução do problema não local em todo o domínio.
5. Discretização: A discretização espacial utiliza elementos triangulares lineares para o campo de deslocamento. Uma "malha-lip" separada (construída a partir dos centroides dos elementos) é usada para armazenar e restringir o campo de dano, garantindo que o número de restrições de Lipschitz seja minimizado.

Principais Contribuições

• Extensão para Dinâmica 2D: O artigo estende com sucesso a abordagem de campo Lip, validada anteriormente em fragmentação dinâmica e casos quase-estáticos 1D, para problemas de fratura dinâmica bidimensional.
• Integração com Dinâmica Explícita: Ele descreve uma implementação específica da abordagem de campo Lip utilizando um esquema de integração temporal explícita, que é computacionalmente eficiente para problemas de impacto de alta velocidade.
• Eficiência Computacional via Limites: Os autores demonstram que a utilização da estimativa de limites permite que a minimização não local dispendiosa seja restrita a pequenos subdomínios, tornando a abordagem de campo Lip computacionalmente competitiva para simulações dinâmicas.
• Distinção Variacional: O trabalho esclarece a distinção entre as abordagens de campo Lip e de Campo de Fase: o campo Lip mantém um funcional de energia local com restrições não locais, enquanto o Campo de Fase utiliza um funcional de energia não local.

Resultados de Simulação
O modelo proposto foi validado através de exemplos numéricos bidimensionais:

1. Tração com Entalhe de Borda Única (SENT): O modelo simulou com sucesso a propagação e ramificação de trincas dinâmicas. Os resultados mostraram convergência de malha em potencial e energia dissipada. O fenômeno de ramificação de trincas foi capturado sem critérios ad-hoc, e as velocidades da ponta da trinca permaneceram geralmente abaixo da velocidade da onda de Rayleigh, exceto por picos transientes durante eventos de ramificação. Os resultados foram comparáveis às simulações de Campo de Fase e TLS encontradas na literatura.
2. Influência da Intensidade de Carga: Ao variar a tração aplicada, o modelo demonstrou a capacidade de capturar múltiplas ramificações de trincas, com a ramificação ocorrendo mais próxima da ponta do entalhe sob cargas mais altas.
3. Teste de Kalthoff-Winkler: O modelo simulou a iniciação e propagação de trincas sob carregamento de cisalhamento dinâmico.
   • Em baixas velocidades de impacto ($33$ m/s), previu um ângulo de propagação de trinca Modo-I de $\approx 62,4^\circ$, consistente com observações experimentais de $\approx 70^\circ$.
   • Em altas velocidades de impacto ($100$ m/s), capturou ramificações sucessivas de trincas e múltiplas trincas.
   • Limitação: Os autores observam que, embora experimentos mostrem uma transição para falha dúctil (fratura por cisalhamento em $\approx -10^\circ$) em altas velocidades, o atual modelo elástico-dano não captura esse comportamento dúctil, pois carece de plasticidade.

Significância e Alegações
O artigo afirma que a abordagem de campo Lip oferece um arcabouço variacional robusto para a fratura dinâmica que preserva a relevância matemática do problema ao mesmo tempo em que introduz as escalas de comprimento necessárias. Sua principal significância reside na sua eficiência computacional alcançada através do solucionador baseado em limites, que evita a resolução de equações não locais globalmente. Os autores afirmam que a abordagem lida efetivamente com fenômenos complexos de fratura, como iniciação, ramificação e coalescência de trincas, sem critérios ad-hoc.

O artigo conclui modestamente que, embora a abordagem de campo Lip seja eficaz para fratura frágil e fragmentação, ela é limitada ao comportamento elástico-dano. Reconhece que simular a falha dúctil em altas taxas de deformação exigiria uma extensão para modelos de dano elástico-plástico, o que está fora do escopo do presente trabalho. O método é apresentado como uma alternativa viável aos métodos de Campo de Fase e TLS, particularmente onde a eficiência computacional em dinâmica explícita é crítica.