Achieving Material Robustness via Symmetric Stress Finite Element Discretizations

Resumo Técnico: Alcançando Robustez Material via Discretizações de Elementos Finitos de Tensão Simétrica

Enunciação do Problema
Na mecânica do contínuo, formulações variacionais como a elasticidade de Hellinger–Reissner e as formulações velocidade–tensão para escoamento incompressível envolvem tensores de tensão simétricos. A simetria do tensor de tensão de Cauchy está fisicamente ligada à conservação do momento angular. Ao discretizar esses problemas usando elementos finitos conformes a $H(\text{div})$, os praticantes devem escolher entre impor a simetria fortemente (exigindo que os tensores discretos sejam simétricos ponto a ponto) ou fracamente (impondo a simetria via um multiplicador de Lagrange). Embora a simetria fraca seja frequentemente preferida na literatura devido à complexidade percebida na construção de elementos fortemente simétricos, este artigo investiga o impacto dessa escolha na precisão numérica, particularmente para materiais com leis constitutivas anisotrópicas.

Metodologia
Os autores desenvolvem um quadro teórico e uma série de benchmarks numéricos para comparar esquemas simétricos fortes e fracos.

1. Formulação Matemática: O estudo considera um problema estático linearizado geral para fluidos e sólidos, onde a tensão $\sigma$ relaciona-se com a deformação $\varepsilon(u)$ e um termo de campo externo $\tilde{F}$ (representando anisotropia). O problema é formulado como um sistema de ponto de sela.
2. Esquemas de Discretização:
   • Fortemente Simétricos: Os esquemas de Johnson-Mercier-Křížek (JMK) e de Hu-Zhang (HZk) são utilizados. Estes impõem $\sigma_h \in \Sigma^{\text{sim}}_h$ ponto a ponto.
   • Fracamente Simétricos: Os esquemas PEERS e de Arnold-Falk-Winther (AFWk) são utilizados. Estes introduzem um multiplicador de Lagrange $\omega_h$ (representando a parte antissimétrica do gradiente) para impor a simetria no sentido fraco.
3. Construção de Benchmarks: Os autores constroem "soluções fabricadas" onde a tensão exata é identicamente zero ($\sigma \equiv 0$), mas os campos de deslocamento e rotação são não triviais. Essas soluções são derivadas de leis constitutivas inspiradas em:
   • Sólidos isotrópicos lineares.
   • Sólidos transversalmente isotrópicos (materiais reforçados com fibras).
   • Fluidos polares (redes de polímeros de cristal líquido).
     A magnitude do termo de força anisotrópica é controlada por um parâmetro de escala $\delta$.
4. Análise Teórica: Uma teoria abstrata de espaços de Hilbert é desenvolvida para problemas genéricos de ponto de sela. Os autores definem discretizações "que preservam estrutura" e derivam estimativas de erro a priori para explicar o comportamento numérico observado.

Principais Contribuições

• Definição de Robustez Material: O artigo introduz o conceito de "robustez material", definido como a propriedade pela qual um esquema numérico produz uma aproximação de tensão zero para qualquer configuração exata de tensão zero, independentemente da lei constitutiva ou da magnitude da força anisotrópica.
• Identificação de Modos de Falha: Os autores demonstram que esquemas fracamente simétricos podem produzir aproximações de tensão arbitrariamente ruins (erros escalando com o parâmetro de força $\delta$) mesmo quando a tensão exata é zero. Isso ocorre especificamente quando o tensor de rotação exato $\omega = \text{skw}(\nabla u)$ não pertence ao espaço discreto do multiplicador de Lagrange $\Xi_h$.
• Teoria Unificadora: Uma teoria rigorosa é apresentada mostrando que esquemas fortemente simétricos que satisfazem propriedades específicas de inclusão de núcleo ($\ker B_h \subseteq \ker B$) são inerentemente robustos materialmente. Em contraste, esquemas fracamente simétricos carecem dessa invariância, a menos que o espaço de rotação discreto contenha o tensor de rotação exato.
• Estimativas de Erro: Novas estimativas de erro a priori são derivadas para discretizações parcialmente que preservam estrutura, mostrando explicitamente como o erro de tensão depende da qualidade da aproximação do tensor de rotação no caso fracamente simétrico.

Resultados
Experimentos numéricos em malhas não estruturadas confirmam as previsões teóricas:

• Caso Isotrópico: Para sólidos isotrópicos lineares onde a rotação exata é constante (e, portanto, facilmente representável), ambos os esquemas forte e fraco performam de forma similar, com erros dominados por arredondamento.
• Casos Anisotrópicos: Para sólidos transversalmente isotrópicos e fluidos polares, onde o tensor de rotação é complexo:
  • Esquemas Fortemente Simétricos (JMK, HZ3): Produzem consistentemente erros de tensão ao nível da precisão de máquina (estado de tensão zero), independentemente do parâmetro de escala $\delta$.
  • Esquemas Fracamente Simétricos (PEERS, AFWk): Produzem erros de tensão que escalam linearmente com $\delta$. Mesmo para esquemas de ordem muito alta (por exemplo, AFW3), se o tensor de rotação exato não estiver contido no espaço discreto $\Xi_h$, o erro de tensão permanece grande.
• Comportamento Transiente: Um exemplo de fluido polar transiente mostra que, embora ambos os esquemas concordem inicialmente, o esquema fracamente simétrico não robusto falha em relaxar para o estado estacionário de tensão zero correto, enquanto o esquema fortemente simétrico robusto o faz.

Significado e Alegações
O artigo alega que a escolha entre imposição de simetria forte e fraca não é meramente uma questão de complexidade de implementação, mas afeta fundamentalmente a fidelidade física da solução numérica. Os autores argumentam que impor a simetria fortemente é crucial para a robustez material, garantindo que esquemas numéricos não gerem tensões espúrias em configurações livres de tensão, particularmente em aplicações envolvendo invariantes rotacionais ou materiais anisotrópicos.

O trabalho defende a adoção mais ampla de elementos finitos fortemente simétricos (como JMK e Hu-Zhang), observando que esses elementos estão cada vez mais disponíveis em softwares de código aberto como Firedrake e NGSolve. Os autores concluem que, para problemas sensíveis à conservação do momento angular, discretizações fortemente simétricas fornecem uma garantia necessária de precisão que esquemas fracamente simétricos não podem oferecer sem condições restritivas sobre os espaços discretos.