A posteriori analysis of a virtual element approach on polytopal meshes for the buckling eigenvalue problem

Título: Análise a posteriori de uma abordagem de Elemento Virtual em malhas poliedrais para o problema de autovalor de flambagem

1. Problema Investigado

O artigo aborda o problema de autovalor de flambagem de placas finas submetidas a carregamentos no plano. Este fenômeno de instabilidade é modelado matematicamente como um problema de autovalor para o deslocamento transversal da placa, governado por um operador de flexão de quarta ordem ($\Delta^2$) acoplado a termos que representam os efeitos das tensões no plano (que podem ser não lineares).

• Objetivo: Determinar a carga crítica (o menor autovalor) na qual a placa perde a estabilidade e descrever os modos de flambagem (autofunções).
• Desafio: A complexidade geométrica das estruturas modernas e a necessidade de lidar com tensões não lineares exigem métodos numéricos flexíveis. Além disso, como a solução exata é geralmente desconhecida, é crucial ter estimadores de erro confiáveis para guiar o refinamento da malha.

2. Metodologia

Os autores propõem uma abordagem baseada no Método de Elemento Virtual (VEM - Virtual Element Method) conformante de classe $C^1$, aplicado tanto em duas (2D) quanto em três dimensões (3D).

• Formulação Variacional: Utiliza-se a formulação variacional em $H^2$ para o problema de autovalor. Diferentemente dos Elementos Finitos Clássicos ($C^1$ FEM), que exigem graus de liberdade (DoFs) elevados e polinômios de alto grau para garantir a continuidade $C^1$, o VEM atinge isso com poucos DoFs (3 por vértice em 2D e 4 por vértice em 3D).
• Malhas: O método é implementado em malhas poliedrais/policonais gerais, permitindo o uso de elementos não convexos e facilitando o tratamento de nós pendentes (hanging nodes) em estratégias adaptativas.
• Estimador de Erro A Posteriori: O núcleo da contribuição é o desenvolvimento de um estimador de erro baseado em resíduos para o problema de flambagem. O estimador global ($\eta$) é composto por:
  • Resíduo de Volume ($\Xi$): Relacionado ao termo de carga crítica e tensões no plano.
  • Resíduo de Salto ($J$): Relacionado à descontinuidade das derivadas normais nas interfaces dos elementos.
  • Termo de Estabilização ($S$): Essencial no VEM para garantir a coercividade.
  • Oscilação de Dados ($\Lambda$): Mede o erro de aproximação polinomial das tensões no plano.
• Implementação: O método foi implementado na biblioteca de código aberto vem++, utilizando o solver SLEPc para resolver o problema de autovalor generalizado.

3. Principais Contribuições

1. Primeiro Estimador para VEM $C^1$ em Flambagem: Este trabalho apresenta o primeiro estimador de erro a posteriori confiável e eficiente para o método VEM conformante $C^1$ aplicado ao problema de flambagem em 2D e 3D.
2. Inclusão de Tensões Não Lineares: O estimador incorpora explicitamente os efeitos de tensões no plano não lineares, permitindo a previsão precisa de modos de flambagem sob carregamentos multiaxiais complexos.
3. Provas Teóricas Rigorosas:
   • Confiabilidade (Reliability): Foi provado que o estimador fornece um limite superior para o erro real, utilizando projeções polinomiais, contribuições de estabilização e estimativas de interpolação.
   • Eficiência (Efficiency): Foi provado que o estimador também fornece um limite inferior (até constantes multiplicativas e termos de ordem superior), utilizando argumentos de funções de bolha (bubble functions).
4. Código Aberto: Disponibilização da implementação no vem++, suportando malhas poliedrais gerais.

4. Resultados Numéricos

Os autores realizaram extensos experimentos numéricos para validar a teoria:

• Refinamento Uniforme: Testes em domínios quadrados (2D) e cúbicos (3D) com tensões isotrópicas mostraram que o método recupera a taxa de convergência ótima esperada ($O(h^2)$ para o erro do autovalor e do estimador), confirmando a confiabilidade e a eficiência do estimador.
• Refinamento Adaptativo:
  • Domínios com Singularidades: Foram testados domínios em "L" (2D) e cubo de Fichera (3D), além de domínios perfurados. O esquema adaptativo, guiado pelo estimador, conseguiu recuperar a taxa de convergência ótima evitando o refinamento desnecessário em regiões onde a solução é suave.
  • Aplicação de Engenharia (Pás de Turbina): Um caso de estudo realista simulou a flambagem de pás de turbina a jato com passagens de resfriamento interno (furos axiais).
    • Foram comparados dois casos de temperatura (ambiente e alta temperatura) e duas condições de contorno (apoiada e engastada).
    • Eficiência Computacional: O refinamento adaptativo alcançou soluções com a mesma precisão do refinamento uniforme, mas utilizando aproximadamente 94% menos graus de liberdade.
    • Física: Os resultados mostraram uma redução de ~34% na capacidade de carga da pá em altas temperaturas e confirmaram que condições de contorno engastadas (CP) suportam cargas críticas maiores que as apoiadas (SSP).

5. Significado e Impacto

Este trabalho é significativo por preencher uma lacuna na literatura numérica, fornecendo uma ferramenta robusta para a análise de estabilidade de estruturas complexas.

• Flexibilidade Geométrica: A capacidade de usar malhas poligonais/poliedrais permite modelar geometrias complexas (como perfurações internas em turbinas) sem a necessidade de gerar malhas de triângulos/tetraedros conformes complexas.
• Eficiência: A estratégia adaptativa baseada no novo estimador reduz drasticamente o custo computacional, tornando viável a análise de precisão de estruturas críticas em engenharia aeroespacial e automotiva.
• Generalidade: A abordagem lida com não-linearidades nas tensões no plano, o que é crucial para aplicações de engenharia real onde o estado de tensão prévio à flambagem não é uniforme.

Em resumo, o artigo estabelece uma base teórica sólida e prática para o uso do VEM $C^1$ em problemas de flambagem, oferecendo um estimador de erro confiável que permite simulações precisas e eficientes em geometrias complexas.