Deformation gradient averaging regularization for third medium contact

Resumo Técnico: Regularização por Média do Gradiente de Deformação para Contato de Terceiro Meio

1. Problema e Contexto

O método de contato de terceiro meio (Third Medium Contact - TMC) tem ganhado popularidade como alternativa aos métodos de contato tradicionais (como penalidade ou multiplicadores de Lagrange), especialmente em problemas onde a busca por superfícies de contato é complexa, como na otimização topológica.

• Conceito: O método preenche o espaço vazio entre corpos com um material fictício extremamente complacente. Este material não interfere na solução antes do contato, mas endurece drasticamente quando comprimido a volumes próximos de zero, impedindo a penetração.
• Desafio: Em simulações de grandes deformações (finite strain), o material de terceiro meio tende a deformar excessivamente na fase pré-contato, levando a instabilidades numéricas.
• Solução Existente: Modelos anteriores utilizam termos de regularização baseados em segundos gradientes de deslocamento (derivadas de segunda ordem). Embora eficazes, essas abordagens exigem:
  • Elementos finitos de segunda ordem (quadráticos); ou
  • Elementos de primeira ordem com graus de liberdade adicionais (DOFs).
  • Isso aumenta a complexidade de implementação e o custo computacional.

2. Metodologia Proposta

Os autores propõem uma nova abordagem de regularização que elimina a necessidade de calcular segundos gradientes de deslocamento explicitamente, permitindo o uso de elementos finitos de primeira ordem (lineares) sem graus de liberdade adicionais.

A formulação do modelo de terceiro meio é composta por três partes principais na densidade de energia de deformação ($W_{3M}$):

1. Termo de Contato (Volumétrico): Baseado na parte volumétrica da lei hiperelástica neo-Hookeana. Reage apenas à compressão volumétrica extrema ($J \to 0$), impedindo a penetração. A parte isocórica (distorsional) é omitida para melhorar a estabilidade.
2. Termo de Rigidez Constante (Pré-contato): Utiliza uma lei elástica linear de pequenas deformações adaptada para grandes deformações. Isso garante que o material tenha uma rigidez constante e previsível antes do contato, evitando distorções excessivas de elementos de baixa rigidez.
3. Termo de Regularização por Média (O Núcleo da Contribuição):
   • Em vez de penalizar o gradiente do gradiente de deformação ($\nabla F$), o método penaliza a desvio do gradiente de deformação ($F$) em relação à sua média no elemento.
   • A energia de regularização é dada por:
     $$W_{\bar{F}} = \frac{1}{2} \kappa_{\bar{F}} \| F - \bar{F} \|^2$$
     Onde:
     • $F$ é o gradiente de deformação calculado nos pontos de integração de Gauss.
     • $\bar{F}$ é o gradiente de deformação avaliado no centroide do elemento.
     • $\kappa_{\bar{F}}$ é um parâmetro de penalidade.
   • Mecanismo: Ao penalizar a diferença entre os valores de $F$ nos pontos de integração e o valor médio no centroide, o método força uma uniformidade do gradiente de deformação dentro do elemento. Isso suprime modos de deformação espúrios (como "hourglassing") sem exigir derivadas de segunda ordem.

3. Principais Contribuições

• Simplicidade de Implementação: A abordagem permite o uso de formulações padrão de elementos finitos de primeira ordem (bilineares), que são amplamente disponíveis em códigos de elementos finitos, sem a necessidade de adicionar graus de liberdade extras ou usar elementos quadráticos.
• Eliminação de Segundos Gradientes: Remove a necessidade de operadores de gradiente de segunda ordem, simplificando a formulação matemática e a implementação numérica.
• Robustez em Grandes Deformações: O termo de rigidez constante combinado com a regularização por média estabiliza o comportamento do material fictício na fase pré-contato.
• Conservação de Energia: Todo o modelo (termos de contato, rigidez e regularização) deriva de um potencial de energia, preservando a simetria do tensor tangente de rigidez global, o que é crucial para a convergência de solvers Newton-Raphson.

4. Resultados e Validação

O método foi implementado no software de elementos finitos de código aberto OOFEM e validado através de quatro benchmarks clássicos:

1. Benchmark em Formato de "C" (C-shape):
   • Estudo de convergência de malha comparando diferentes refinamentos.
   • Os resultados mostraram que, ao ajustar o parâmetro de penalidade $\kappa_{\bar{F}}$ com o tamanho da malha, o método produz resultados consistentes e comparáveis ao método tradicional de contato "nó-a-segmento" (Node-to-Segment), com erros de contato controlados.
2. Caixa Fechada (Closed Box):
   • Simulação de auto-contato de uma caixa sob carregamento mecânico.
   • O método reproduziu com alta precisão a curva força-deslocamento e os estados deformados, coincidindo quase perfeitamente com a solução de contato tradicional.
3. Flambagem Restrita e Auto-contato:
   • Problema de instabilidade envolvendo uma viga bi-apoiada comprimida por um pistão.
   • O modelo capturou corretamente a carga crítica de Euler, o endurecimento pós-flambagem ao tocar a parede, e o contato subsequente entre partes da própria viga, demonstrando robustez em cenários de bifurcação e contato complexo.
4. Forças Configuracionais:
   • Validação contra uma solução analítica fechada para a força configuracional em um sólido elástico em contato com um indentador rígido.
   • O método de terceiro meio recuperou com excelente precisão a força horizontal reativa e a energia liberada, validando a capacidade do modelo de lidar com grandes deformações e deslizamento.

5. Significado e Conclusão

O artigo apresenta uma evolução significativa na aplicação do método de terceiro meio. Ao substituir a complexa regularização de segundos gradientes por uma regularização baseada na média do gradiente de deformação, os autores criaram uma ferramenta:

• Universal: Aplicável a problemas 2D e 3D.
• Eficiente: Não aumenta o número de graus de liberdade do sistema.
• Fácil de Implementar: Pode ser integrada em códigos existentes de elementos finitos com modificações mínimas (apenas no cálculo da energia e das matrizes de rigidez elementares).

A técnica é particularmente valiosa para a comunidade de otimização topológica, onde a detecção de contato é frequentemente problemática, e também para simulações multiphísicas que envolvem o mapeamento de espaços "tipo ar" ou vazios. A robustez demonstrada nos benchmarks confirma que o método é uma alternativa viável e superior em termos de facilidade de uso aos modelos de regularização existentes.