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

Titolo

Analisi a posteriori di un approccio agli elementi virtuali su mesh poliedriche per il problema agli autovalori di instabilità (buckling).

1. Il Problema

Il lavoro si concentra sul problema agli autovalori di instabilità (buckling) di piastre sottili soggette a carichi nel piano. Questo fenomeno fisico è modellato matematicamente come un problema agli autovalori per lo spostamento trasversale della piastra.

• Equazione governante: L'equazione coinvolge un operatore di flessione del quarto ordine ($\Delta^2$) combinato con termini aggiuntivi che tengono conto degli effetti degli sforzi nel piano (che possono essere non lineari).
• Significato fisico: A differenza dei problemi vibrazionali, l'autovalore qui rappresenta un parametro di carico critico. Il buckling si verifica quando il carico raggiunge valori per cui il problema omogeneo ammette soluzioni non banali. L'autovalore più piccolo determina il primo carico di instabilità.
• Sfida computazionale: La risoluzione di questi problemi richiede discretizzazioni che rispettino la continuità $C^1$ (necessaria per gli operatori del quarto ordine). I metodi agli elementi finiti classici (FEM) $C^1$ richiedono solitamente gradi di libertà (DoF) elevati e polinomi di alto grado. Inoltre, la gestione di geometrie complesse e fenomeni localizzati (come singolarità negli angoli rientranti o bordi interni) richiede strategie di adattività efficienti.

2. Metodologia

Gli autori propongono un approccio basato sul Metodo degli Elementi Virtuali (VEM) conforme di ordine $C^1$, applicato sia in due (2D) che in tre dimensioni (3D).

• Formulazione Variazionale: Viene utilizzata una formulazione variazionale nello spazio $H^2$. Il metodo VEM permette di utilizzare mesh poligonali (2D) e poliedriche (3D) generali, inclusi elementi non convessi e nodi "hanging" (sospesi), semplificando notevolmente l'implementazione di strategie di raffinamento adattivo.
• Discretizzazione:
  • Viene impiegata la versione a ordine più basso del VEM conforme ($k=2$), che richiede solo 3 DoF per vertice in 2D e 4 in 3D, riducendo significativamente il costo computazionale rispetto ai FEM classici.
  • Vengono definiti proiettori polinomiali ($\Pi^0, \Pi^\nabla, \Pi^{\nabla^2}$) e termini di stabilizzazione per rendere computabili le forme bilineari discrete.
• Stima dell'Errore A Posteriori:
  • Viene sviluppato un estimatore di errore basato sul residuo per il problema agli autovalori.
  • L'estimatore è composto da quattro termini principali:
    1. Residuo di volume ($\Xi$): Legato all'equazione differenziale interna.
    2. Residuo di salto ($J$): Legato alla discontinuità delle derivate normali attraverso le facce degli elementi.
    3. Oscillazione dei dati ($\Lambda$): Misura l'errore di approssimazione polinomiale degli sforzi nel piano.
    4. Termino di stabilizzazione ($S$): Derivato dal termine di stabilizzazione del VEM.
  • Proprietà Teoriche: Viene dimostrata rigorosamente la affidabilità (upper bound dell'errore reale) e l'efficienza (lower bound, a meno di termini di ordine superiore e oscillazione) dell'estimatore. La dimostrazione utilizza proiezioni polinomiali, funzioni di bolla e stime di interpolazione.
• Implementazione: Il metodo è stato implementato nella libreria open-source vem++, utilizzando il solver SLEPc per il problema agli autovalori generalizzato.

3. Contributi Chiave

Questa ricerca rappresenta il primo lavoro che:

1. Sviluppa un estimatore di errore a posteriori affidabile ed efficiente basato sul residuo per il metodo VEM conforme $C^1$ applicato al problema di buckling in 2D e 3D.
2. Include esplicitamente gli effetti degli sforzi nel piano non lineari nell'estimatore, permettendo una previsione accurata delle modalità di buckling in piastre sottoposte a carichi multi-assiali complessi.
3. Fornisce un'implementazione open-source completa che supporta mesh poliedriche generali, facilitando l'analisi di strutture con geometrie complesse.

4. Risultati Numerici

Gli esperimenti numerici confermano la teoria e dimostrano l'efficacia del metodo:

• Raffinamento Uniforme (Esempio 1):
  • Su domini semplici (quadrato, cubo) con mesh Voronoi, perturbate e quadrate, il metodo recupera il tasso di convergenza ottimale $O(h^2)$ per l'errore dell'autovalore e dell'estimatore globale.
  • L'indice di effettività (rapporto tra errore stimato e errore reale) rimane limitato e costante, confermando l'affidabilità dell'estimatore.
• Raffinamento Adattivo (Esempio 2):
  • Su domini con singolarità geometriche (domini a L, cubo di Fichera, domini perforati), l'algoritmo adattivo (SOLVE-ESTIMATE-MARK-REFINE) riesce a recuperare l'ottimalità della convergenza concentrandosi sulle regioni critiche (angoli rientranti, bordi dei fori).
  • Vengono testati diversi valori di tensori degli sforzi nel piano ($\kappa$), dimostrando la robustezza del metodo anche in presenza di carichi non uniformi.
• Applicazione Ingegneristica (Esempio 3):
  • Studio del comportamento di buckling di una pala di turbina con passaggi di raffreddamento interno (fori assiali).
  • Vengono analizzati due casi di temperatura (ambiente e alta temperatura) e due condizioni al contorno (piastra semplicemente appoggiata SSP e incassata CP).
  • Risultati significativi:
    • La strategia adattiva richiede circa il 94% in meno di gradi di libertà rispetto al raffinamento uniforme per raggiungere la stessa accuratezza.
    • Viene quantificato l'impatto della temperatura: la capacità di carico della pala diminuisce di circa il 34% passando dalla temperatura ambiente a quella elevata (a causa della riduzione del modulo di Young).
    • Le condizioni al contorno incassate (CP) forniscono fattori di carico critico superiori rispetto alle condizioni semplicemente appoggiate (SSP), in linea con il comportamento fisico atteso.

5. Significato e Impatto

Questo lavoro è significativo per diversi motivi:

• Avanzamento Teorico: Colma una lacuna nella letteratura fornendo la prima analisi a posteriori completa per il problema di buckling con VEM $C^1$ in 3D, includendo effetti di sforzo non lineari.
• Efficienza Computazionale: Dimostra che l'uso di mesh poliedriche e metodi adattivi può ridurre drasticamente il costo computazionale (riduzione dei DoF) senza sacrificare l'accuratezza, rendendo fattibili simulazioni di strutture complesse che sarebbero proibitive con FEM tradizionali.
• Applicabilità Industriale: L'applicazione alle pale delle turbine mostra come il metodo possa essere utilizzato per l'analisi predittiva di sicurezza in settori critici come l'aerospaziale, permettendo di valutare l'impatto di condizioni operative estreme (temperature elevate, geometrie interne complesse) sulla stabilità strutturale.
• Open Source: La disponibilità del codice nella libreria vem++ favorisce la riproducibilità e l'adozione di queste tecniche avanzate nella comunità scientifica e ingegneristica.