An Investigation of Stabilization Scaling in Finite-Strain Virtual Element Methods for Hyperelasticity

Este trabalho propõe uma nova estratégia de estabilização livre de submalha para o Método de Elemento Virtual em grandes deformações hiperelásticas, que desacopla os canais volumétrico e deviatório para eliminar a sensibilidade à tesselação interna e garantir robustez numérica no regime quase incompressível, superando as limitações das abordagens baseadas em energia de substituição.

O essencial

Imagine que você está tentando construir uma casa usando blocos de formas estranhas e irregulares (pentágonos, hexágonos, etc.), em vez dos quadrados e retângulos perfeitos que usamos normalmente. Essa é a ideia por trás do Método de Elemento Virtual (VEM): uma forma inteligente de simular como materiais (como borracha ou tecidos biológicos) se deformam sob pressão.

No entanto, há um problema. Quando usamos formas irregulares, o computador pode "alucinar" e criar modos de deformação que não deveriam existir (como se a parede estivesse tremendo sem motivo). Para corrigir isso, os engenheiros adicionam um "remendo" matemático chamado estabilização.

O artigo que você enviou discute como esse "remendo" está sendo feito hoje e propõe uma maneira muito melhor de fazê-lo, especialmente para materiais que são quase impossíveis de comprimir (como a água ou o corpo humano).

Aqui está a explicação simples, usando analogias:

1. O Problema: O "Remendo" Quebrado

Hoje, a maneira padrão de fazer esse remédio (estabilização) é como se você estivesse tentando consertar um buraco no chão de um quarto irregular.

• O jeito antigo: O engenheiro pega um pedaço de papel, corta em triângulos pequenos (sub-triângulos) dentro do quadrado irregular e cola o remendo nesses triângulos.
• O defeito: O remendo fica dependente de como você cortou o papel. Se você cortar os triângulos de um jeito ou de outro, a força do remendo muda. Pior ainda, para materiais que não podem ser espremidos (como uma bola de água), esse remendo antigo começa a usar "medidas de volume" para consertar problemas de "cisalhamento" (distorção).
• A consequência: É como tentar segurar uma esponja molhada apertando-a com uma força de "volume" quando você só deveria estar segurando a "forma". O material fica artificialmente duro e não se deforma como deveria. Isso é chamado de travamento (locking). O computador acha que o material é muito mais rígido do que realmente é.

2. A Solução Proposta: O "Remendo" Inteligente e Limpo

Os autores (Paulo Enabe e Rodrigo Provasi) propõem uma nova maneira de fazer o remendo, sem precisar cortar o papel em triângulos extras.

• Sem sub-triângulos: Em vez de olhar para dentro do bloco, eles olham apenas para as bordas e os vértices. É como verificar se a moldura da janela está reta, sem precisar entrar dentro da parede.
• Separação de Funções (O Segredo): Eles dividem o remendo em duas caixas separadas:
  1. Caixa de Cisalhamento (Forma): Controla apenas a distorção da forma. É escalada apenas pela "rigidez ao cisalhamento" do material.
  2. Caixa de Volume (Compressão): Controla apenas se o material está sendo espremido.
• A Mágica: No jeito antigo, se o material fosse quase incompressível, a "Caixa de Volume" vazava para dentro da "Caixa de Cisalhamento", endurecendo tudo. No novo método, elas são desacopladas. A caixa de cisalhamento nunca recebe "lixo" da caixa de volume.

3. A Analogia do Orquestra

Pense na simulação como uma orquestra tocando música.

• O jeito antigo: O maestro (o computador) está tentando ouvir os violinos (a física real), mas os tímpanos (o remendo de estabilização) estão tocando muito alto e de forma errada, cobrindo a música. Se o maestro tentar ajustar o volume dos tímpanos baseando-se no som dos violinos, ele acaba estragando a melodia inteira, especialmente quando a música fica muito lenta (regime quase incompressível).
• O novo método: O maestro coloca os tímpanos em uma sala separada com isolamento acústico perfeito. Os tímpanos tocam apenas o que precisam para manter o ritmo (estabilidade), sem interferir no som dos violinos. A música fica limpa, precisa e a orquestra não "trava".

4. Por que isso importa? (O Teste da Membrana de Cook)

Os autores testaram isso em um cenário clássico chamado "Membrana de Cook" (uma peça de borracha em forma de trapézio que é torcida).

• Com o método antigo: Quando a borracha é quase incompressível (como borracha de silicone), o método antigo diz que a peça é super dura e quase não se move. É como se a borracha virasse pedra.
• Com o novo método: A peça se deforma suavemente e corretamente, independentemente de quão estranhas sejam as formas dos blocos usados na simulação.

Resumo Final

Este artigo é sobre limpar a bagunça matemática.
Os autores criaram uma fórmula matemática que age como um filtro inteligente. Ela garante que o computador corrija apenas os erros que realmente existem, sem adicionar "peso" extra desnecessário que faria o material parecer mais duro do que é.

Isso permite que engenheiros e cientistas simulem materiais reais (como tecidos humanos, pneus ou borrachas) com muito mais precisão, usando malhas de formas irregulares, sem que o computador fique "travado" ou dê resultados errados quando o material é muito elástico ou incompressível. É uma melhoria fundamental para a precisão de softwares de engenharia do futuro.

Resumo técnico

Resumo Técnico: Investigação da Escala de Estabilização em Métodos de Elemento Virtual (VEM) de Deformação Finita para Hiperelasticidade

1. O Problema

O Método de Elemento Virtual (VEM) é uma estrutura conformada para resolver equações diferenciais parciais em malhas poligonais e poliedrais gerais. Em formulações de baixa ordem para hiperelasticidade de deformação finita, a contribuição consistente da energia é calculada através de projeções polinomiais. No entanto, os modos de deformação que residem no núcleo (kernel) do projetor — ou seja, os modos não capturados pela projeção polinomial — não contribuem para a energia consistente e devem ser controlados exclusivamente por um termo de estabilização.

A prática atual de estabilização em VEM hiperelástico enfrenta três problemas críticos:

1. Dependência de Sub-malha: Muitas formulações avaliam a estabilização integrando uma energia de "surroga" (substituta) sobre uma triangulação interna auxiliar (sub-malha). Isso introduz uma dependência arbitrária na tesselação interna, que não é intrínseca ao espaço VEM.
2. Contaminação Volumétrica em Modos Isochóricos: Em regimes quase incompressíveis (onde o coeficiente de Poisson $\nu \to 1/2$), as formulações clássicas frequentemente injetam proxies volumétricos (dependentes do módulo de compressibilidade $\lambda$) nas penalidades de tipo cisalhante. Isso resulta na atribuição de energia volumétrica excessiva a modos de deformação que deveriam ser puramente isochóricos (sem mudança de volume), causando um endurecimento artificial (locking) dos modos não resolvidos.
3. Complexidade de Linearização: A integração sobre sub-malhas e o uso de energias não lineares complexas complicam a linearização consistente de Newton, tornando o processo computacionalmente mais custoso e menos robusto.

O objetivo deste trabalho é desenvolver uma estratégia de estabilização que seja livre de sub-malha, desacople os canais deviatório e volumétrico, e escale de forma espectralmente equivalente à energia tangente de Newton atual no núcleo do projetor, mantendo robustez em malhas distorcidas e no limite incompressível.

2. Metodologia

Os autores propõem uma nova formulação de estabilização baseada nos seguintes pilares:

2.1. Estabilização Apenas no Núcleo (Kernel-Only)

Em vez de usar sub-malhas internas, a estabilização é construída diretamente a partir do resíduo do projetor nos graus de liberdade (nós) do elemento.

• Define-se o resíduo nodal $r_i = u_i - (\Pi^\nabla_E u)_i$.
• Como o projetor $\Pi^\nabla_E$ mapeia campos polinomiais para si mesmos, o resíduo é zero para campos polinomiais. Assim, a estabilização atua exclusivamente no espaço complementar (o kernel), preservando a consistência polinomial sem necessidade de integração interna.

2.2. Desacoplamento Deviatório/Volumétrico

A energia de estabilização é dividida em dois termos independentes:

• Termo Deviatório: Escalado estritamente pelo módulo de cisalhamento $\mu_E$. Incorpora pesos direcionais baseados na geometria do elemento (usando os autovetores do tensor de segundo momento geométrico) para lidar com anisotropia em polígonos distorcidos.
• Termo Volumétrico: Escalado por um módulo de bulk $\kappa_E$. É formulado como uma penalidade apenas nas componentes normais do resíduo nas bordas do elemento.
• Controle no Limite Incompressível: O termo volumétrico pode ser suprimido ou limitado (capado) quando $\nu \to 1/2$, evitando que proxies volumétricos inflacionem a rigidez deviatória.

2.3. Abordagem Espectral

Os autores estabelecem um quadro teórico onde a estabilidade é caracterizada por quocientes de Rayleigh generalizados e limites de autovalores.

• A condição de estabilidade exige que a forma bilinear de estabilização seja espectralmente equivalente à energia tangente de Newton restrita ao kernel.
• Isso permite diagnosticar se a estabilização atribui energias físicas corretas aos modos não resolvidos, evitando modos "hourglass" (energia zero) ou endurecimento artificial.

3. Principais Contribuições

1. Perspectiva Espectral para VEM de Deformação Finita: Reformulação do requisito de estabilidade como uma equivalência espectral entre a estabilização e a energia tangente no kernel, fornecendo uma métrica independente de base para avaliar a atribuição de rigidez.
2. Novo Termo de Estabilização Desacoplado e Livre de Sub-malha: Proposta de uma estabilização baseada apenas em resíduos nodais, que é exatamente linearizável (tangente consistente de Newton) e evita a necessidade de triangulações auxiliares.
3. Escalonamento Direcional Consciente de Modos: Incorporação de pesos de anisotropia baseados na geometria para redistribuir a rigidez em polígonos distorcidos, corrigindo a falha de escalonamento escalar clássico em elementos anisotrópicos.
4. Resultados Teóricos de Escalonamento Uniforme em $\nu$: Prova teórica de que o termo deviatório proposto é equivalente a $\mu_E |\cdot|_{1,E}^2$ no kernel, com constantes independentes do tamanho da malha ($h$) e do coeficiente de Poisson ($\nu$).
5. Diagnóstico e Validação Numérica: Introdução de testes de diagnóstico de nível de elemento que isolam a contaminação volumétrica, demonstrando que as formulações clássicas falham em regimes quase incompressíveis, enquanto a proposta permanece robusta.

4. Resultados

Os resultados foram validados através de diagnósticos de nível de elemento e benchmarks globais:

• Diagnóstico de Modo Isocórico de Elemento Único:

  • Ao testar um modo de deformação isocórico (sem mudança de volume) em um kernel, a estabilização clássica mostrou um aumento drástico na energia atribuída à medida que $\nu \to 1/2$, indicando que ela está escalando o modo com um parâmetro de cisalhamento efetivo inflado por proxies volumétricos.
  • A estabilização proposta manteve-se constante e escalada apenas pelo módulo de cisalhamento físico $\mu$, confirmando a ausência de contaminação volumétrica.

• Análise Modal Restrita ao Kernel:

  • A análise de autovalores confirmou que a estabilização proposta anula-se exatamente em modos polinomiais.
  • A separação deviatória/volumétrica foi clara: modos puramente deviatórios não foram afetados pela variação do parâmetro de bulk $\kappa_E$.
  • A condicionamento da matriz de estabilização mostrou-se robusto, com crescimento previsível em elementos retangulares distorcidos, controlável pelos pesos de anisotropia.

• Benchmark da Membrana de Cook (Regime Quase Incompressível):

  • Simulações com $\nu = 0.499$ em diversas famílias de malhas (quadriláteros regulares, distorcidos, altamente distorcidos e Voronoi) revelaram um comportamento distinto.
  • Estabilização Clássica: Exibiu sintomas severos de locking volumétrico, subestimando drasticamente o deslocamento da ponta em malhas grosseiras e convergindo lentamente para o valor de referência à medida que a malha era refinada.
  • Estabilização Proposta: Produziu uma resposta coerente e livre de locking em todas as famílias de malhas. O deslocamento da ponta convergiu rapidamente para um valor de referência estável (aprox. 8.5), independentemente da distorção geométrica ou do tipo de malha, demonstrando superioridade no regime quase incompressível.

5. Significância

Este trabalho resolve um problema fundamental na aplicação de VEM a materiais hiperelásticos quase incompressíveis. Ao demonstrar que a dependência de sub-malhas e a mistura de escalas volumétricas e deviatórias são fontes de instabilidade e locking, os autores propõem uma alternativa teoricamente fundamentada e computacionalmente eficiente.

A principal contribuição é a mudança de paradigma: tratar a estabilização não como uma energia de surroga complexa, mas como um mecanismo de atribuição de energia espectral no kernel do projetor. A nova formulação garante:

• Robustez Geométrica: Funciona bem em polígonos distorcidos e irregulares.
• Robustez Material: Mantém a precisão e a convergência no limite incompressível, eliminando o locking artificial.
• Eficiência Computacional: Elimina a necessidade de gerar e integrar sobre sub-malhas internas, simplificando a implementação e a linearização de Newton.

Os resultados sugerem que esta abordagem deve ser adotada como padrão para VEM de deformação finita em aplicações de engenharia que envolvem borrachas, polímeros e tecidos biológicos, onde a incompressibilidade e a distorção geométrica são comuns.