Djinious
Analisi strutturaleProduzione e robotica

FEA strutturale e modale

Lo strumento di simulazione più usato nell’ingegneria — discretizza una struttura in elementi, risolve per spostamento e sollecitazione sotto carico, ed estrae le sue frequenze naturali e forme modali da un problema agli autovalori. Sollecitazione di una trave a sbalzo, convergenza della mesh, analisi modale, concentrazione di sollecitazione, instabilità e una staffa reale.

DjiniousLab
Il campo di sollecitazione di von Mises di una trave a sbalzo caricata in punta, massimo alla radice incastrata e decrescente verso l’estremità libera
deformazione della trave a sbalzo contro la teoria delle travi
0.58%deformazione della trave a sbalzo contro la teoria delle travi
frequenze modali contro soluzione analitica
0.0001%frequenze modali contro soluzione analitica
fattore di concentrazione di sollecitazione a un foro
3.3fattore di concentrazione di sollecitazione a un foro
carico di instabilità contro Eulero
0.00%carico di instabilità contro Eulero
notebook di progettazione
10notebook di progettazione
requisiti PASS
5 / 5requisiti PASS

Si romperà, e come vibra?

L’analisi a elementi finiti risponde alle due domande a cui ogni struttura deve rispondere: quanto si deforma e dove si rompe sotto carico, e come vibra. Discretizzi il pezzo in migliaia di piccoli elementi, assembli una matrice di rigidezza gigante, e risolvi — per spostamento e sollecitazione, o per le frequenze naturali e le forme modali che escono da un problema agli autovalori. È lo strumento di simulazione più diffuso nell’ingegneria (ANSYS Mechanical, Nastran), e questo programma lo costruisce dall’elemento in su — validato, a ogni passo, contro le formule che puoi verificare a mano.

Dove il materiale lavora più duramente.

Carica una trave a sbalzo alla sua estremità e la sollecitazione non è uniforme — si concentra alla radice incastrata, sulle superfici superiore e inferiore dove la flessione è maggiore, e cade a zero all’asse neutro e all’estremità libera. Il risolutore a elementi finiti calcola l’intero campo di sollecitazione di von Mises che mostra esattamente questo, e la deformazione in punta che predice cade entro lo 0,58% della formula della trave di Eulero-Bernoulli. Il campo ti dice dove aggiungere materiale; l’accordo ti dice di fidartene.

DjiniousLab
Il campo di sollecitazione di von Mises di una trave a sbalzo caricata in punta, concentrato alla radice incastrata
Il protagonista: il campo di sollecitazione di von Mises di una trave a sbalzo caricata in punta. La sollecitazione si concentra sulle fibre superiore e inferiore della radice incastrata — dove una trave si rompe — e decade verso l’estremità libera. La deformazione predetta corrisponde alla teoria delle travi di Eulero-Bernoulli entro lo 0,58%. È questo lo scopo della FEA: vedere dove un pezzo lavora più duramente.

Ogni struttura ha note a cui vuole risuonare.

Oltre alla resistenza statica, una struttura ha frequenze naturali — le velocità a cui vuole vibrare — e una forma modale per ciascuna. Trovale risolvendo un problema agli autovalori generalizzato sulle matrici di rigidezza e massa, e ne escono le risonanze da tenere lontane dal traffico di un ponte o dalla vibrazione di un motore. Su una trave a sbalzo le frequenze calcolate corrispondono ai valori analitici fino alla quarta cifra decimale, e le forme modali sono le classiche forme di flessione da manuale.

DjiniousLab
Le prime quattro forme modali di flessione di una trave a sbalzo, con un numero crescente di nodi
Analisi modale: le prime quattro forme modali di vibrazione di una trave a sbalzo, ciascuna con un nodo in più della precedente, da un problema agli autovalori sulle matrici di rigidezza e massa. Le frequenze naturali corrispondono alla teoria analitica delle travi entro lo 0,0001%. Conoscere queste risonanze è il modo per evitare che una struttura si scuota fino a disgregarsi.

Fori, angoli e colonne snelle.

I pezzi reali si rompono nei loro dettagli. Un foro in una piastra caricata triplica la sollecitazione locale — il fattore di concentrazione di sollecitazione di circa 3 che predice la soluzione analitica di Kirsch, e che il risolutore recupera al restringersi del foro. Una colonna snella non snerva, si instabilizza, al carico di Eulero che il risolutore riproduce alla cifra. E l’angolo rientrante di una staffa a L è un innalzatore di sollecitazione che ne fissa il fattore di sicurezza. Questi sono i dettagli che decidono se un pezzo sopravvive.

DjiniousLab
Una piastra caricata con un foro, che mostra una concentrazione di sollecitazione al bordo del foro circa tre volte la sollecitazione a campo lontano
Concentrazione di sollecitazione: un foro in una piastra caricata alza la sollecitazione locale a circa 3 volte quella a campo lontano — il fattore di concentrazione di sollecitazione che ogni progettista conosce, recuperato qui (Kt ≈ 3,3, che converge al valore analitico di 3 al restringersi del foro). Il motivo per cui un foro, e non il materiale in massa, è di solito dove un pezzo si incrina.
DjiniousLab
Un traliccio a ponte Warren con le aste colorate per trazione e compressione sotto un carico a metà campata
Un traliccio, l’altro classico: un ponte Warren con ogni asta colorata a seconda che sia in trazione o compressione sotto un carico a metà campata. Lo stesso FEM, applicato a un telaio di membri snelli — il pane quotidiano dell’ingegneria strutturale.

Da un campo di sollecitazione a un fattore di sicurezza.

Mettilo insieme su un componente reale — una staffa a L incastrata lungo la parte superiore con un carico sul suo braccio. Il risolutore trova la sollecitazione di picco all’angolo rientrante, e dividendo per essa la resistenza di snervamento del materiale si ottiene il fattore di sicurezza, il singolo numero su cui ruota una revisione di progetto. Esplora il carico e leggi l’ammissibile: quanto puoi spingere il pezzo prima che il margine di sicurezza si esaurisca. Questo è ciò che la FEA fa davvero di mestiere — trasformare una geometria e un carico in un go/no-go.

DjiniousLab
Il campo di sollecitazione di von Mises di una staffa a L con un picco all’angolo rientrante
Un pezzo reale: il campo di von Mises di una staffa a L, con la sollecitazione di picco all’angolo rientrante — 78,7 MPa sotto carico, un fattore di sicurezza di 3,2 contro lo snervamento dell’acciaio. Esplora il carico e ottieni l’ammissibile. Il passaggio da geometria-più-carico a go/no-go che la FEA esiste per fornire.

Ogni numero ridedotto in fase di collaudo finale.

Il notebook V&V ricostruisce il solver da zero e ridefinisce ogni requisito rispetto alla teoria in forma chiusa, stampando un quadro PASS/FAIL.

Risultato

  • Deflessione della mensola vs teoria delle travi: 0,58%
  • Frequenze modali vs soluzione analitica: 0,0001%
  • Convergenza della mesh: monotona
  • Fattore di concentrazione delle tensioni: 3,3 → 3
  • Carico di instabilità vs Euler: 0,00%

Requisito

  • Deflessione della mensola vs teoria delle travi: R-01 < 3%
  • Frequenze modali vs soluzione analitica: R-02 < 0,1%
  • Convergenza della mesh: R-03
  • Fattore di concentrazione delle tensioni: R-04
  • Carico di instabilità vs Euler: R-05 < 1%

Un solver vero, a fedeltà didattica.

Il deliverable è composto da dieci notebook e dal dossier — un campo strutturale discretizzato con un problema agli autovalori, non una rete concentrata acausale, quindi non c’è alcun blocco personalizzato né canvas. Utilizza elementi Q4 plane-stress ed elementi trave di Eulero-Bernoulli, deformazione elastica lineare a piccole deformazioni, un solver diretto denso su mesh di dimensioni modeste e recupero delle tensioni nel punto centrale — nessun contatto, plasticità o grandi deformazioni. Ciò che il programma dimostra è la capacità che mancava di più al portfolio: la meccanica strutturale agli elementi finiti, sia le tensioni statiche sia le vibrazioni modali, validate rispetto alla teoria in forma chiusa a ogni passaggio — ed è lo stesso solver FEM che il flagship di ottimizzazione topologica utilizza nel proprio loop, la controparte per il campo strutturale dei solver di campo CFD ed elettromagnetico.