Download - Centro "E.Piaggio"
Download
Report
Transcript Download - Centro "E.Piaggio"
Elementi Finiti:
Analisi Strutturale
[email protected] + Analisi agli elemen4 fini4 • Il FEM è un metodo numerico (pertanto approssimato) che perme>e la risoluzione di equazioni differenziali alle derivate parziali. • Il metodo degli elemen4 fini4 consiste nella discre'zzazione di un assegnato dominio in elemen% fra loro connessi in un numero finito di pun4, ver4ci degli elemen4 chiama4 nodi, in corrispondenza dei quali sono valutate le componen4 della funzione incognita. • Il valore della funzione all'interno del singolo elemento è o>enuto sulla base dei valori dei parametri nodali a>raverso l'uso di opportune funzioni di forma. • La scelta di tali funzioni, come pure del 4po di mesh con cui discre4zzare il dominio è di importanza cruciale per una corre>a convergenza della soluzione. + Nota • Esistono altri strade che possono portare alla formulazione della “matrice fondamentale” – Metodi variazionali (principio dei lavori virtuali) (vedi dispensa) – Formulazione dire>a (vedi dispensa) – Minimizzazione di un funzionale (energia potenziale totale) + Elemento asta ELEMENTO ASTA/1
• Travature re4colari piane e spaziali • solo sforzo normale ELEMENTO ASTA/3 - TRALIC
• 2 nodi • 2 o 3 g.d.l /nodo • carichi applicabili Travature reticolari piane e spaziali
• solo sforzo normale
solo nei nodi • 2 nodi
• Car. geometriche: A• 2 o 3 g.d.l /nodo
• carichi applicabili solo nei nodi
• Car. geometriche: A
+ Elemento trave • Equazione della linea elas4ca • 2 nodi • 3 gdl/nodo • Carichi concentra4 e distribui4 • Cara>eris4che geometriche (sezione, momento d'inerzia, ...) + Lastra (plane stress) • Stati piani di tensione: • sono caratterizzati
dall’avere una delle
componenti principali di
tensione identicamente
nulla • si verificano tipicamente
in corpi piani, di spessore
piccolo rispetto alle altre
dimensioni caratteristiche
del problema, caricati nel
loro piano medio. Possibilità di inserire lo spessore del corpo
+ Modelli di omogenizzazione • Modello di Voigt A1
E1 F
E2 • Modello di Reuss A2
F
l
A1
ν1 =
A1 + A2
E1 l1
E2 l2
F
F
E = E1ν1 + E2 (1− ν1 )
A
E1E2
E=
E1 (1− f1 ) + E2 f1
l1
f1 =
l
l
+ Esercizio 1 • Valutare il modulo elas4co complessivo dei seguen4 corpi della precedente diaposi4va con il modello anali4co e con quello ad elemen4 fini4 (u4lizzare l'analisi plane stress). + Nota esercizio 1 • I modelli di Reuss e Voigt non prendono in considerazione carichi di 4po trasversale. Per introdurre questo conce>o è necessario porre il modulo di Poisson pari a 0, in modo tale che deformazioni normale provochino deformazioni (e quindi carichi) trasfersali. • Il carico da imporre nel modello di modello di Voigt (o di isoderformazione) è quello di uno spostamento in direzione normale in modo da avere una isodeformazione su entrambi blocchi. + Esempio soluzione esercizio 1 (1/4) Modello di Voigt
A1
E1 F
E2 A2
l
F
E = E1ν1 + E2 (1− ν1 )
A1
ν1 =
A1 + A2
L = 0.1 m Spessore = 0.1 m A1 = 0.03*0.1 m2 = 0.003 m2 A2 = 0.07*0.1 m2 = 0.007 m2 ν1 = 0.3 ν2 = 0.7 Se E1= 10 E2= 100 GPa E = 37 GPa Per avere una deformazione del -‐10% lungo la direzione y devo applicare una forza pari a = F = (E * ε) * A = 73 GPa * (-‐0.1) * 0.01 m2 = = -‐3.7 * 107 N + Esempio soluzione esercizio 1 (2/4) + Esempio soluzione esercizio 1 (3/4) Modello di Reuss
F
E1 E2 A
l1
l2
F
E1E2
E=
E1 (1− f1 ) + E2 f1
l1
f1 =
l
l
L = 0.1 m Spessore = 0.1 m l1 = 0.03 m l2 = 0.07 m f1 = 0.3 f2 = 0.7 A = 0.1 * 0.1 m2 = 0.01 m2 Se E1= 10 E2= 100 GPa E = 13.7GPa Applicando una forza pressione in direzione y di 1 kPa o>engo uno spostamento totale di Δy = 0.1m * (-‐1 kPa / 13.7 GPa) = 7.3*10-‐9 m + Esempio soluzione esercizio 1 (4/4) Spostamento dell’intera stru>ura valutato lungo la direzione y + Considerazioni di simmetria (1/5) • L’uso di considerazioni di simmetria consente di ridurre le dimensioni del modello. I più comuni 4pi di simmetria sono: – Simmetria speculare o di riflessione – Simmetria polare o di rotazione + Considerazioni di simmetria (2/5) Sfruttando la simmetria è possibile includere nel modello so
parte della struttura, sostituendo la parte mancante con oppo
vincoli posti sul piano di divisione.
• Sfru>ando la simmetria è possibile includere nel modello solo una parte della stru>ura, sos4tuendo la parte mancante con opportuni vincoli pos4 sul piano di divisione + •
Considerazioni di simmetria (3/5) I carichi non devono necessariamente essere simmetrici, dato che
condizione di carico qualsiasi può essere scissa in una
I carichi non devono una
componente simmetrica ed in una antisimmetrica.
necessariamente F
essere simmetrici, dato che una condizione di carico qualsiasi può essere scissa in una componente simmetrica ed in una an4simmetrica. F/2
F/2
F/2
F/2
+ Considerazioni di simmetria (4/5) • Simmetria di riflessione • La stru>ura viene tagliata Simmetria
di riflessione
La struttura
in corrispondenza del viene tagliata in corrispondenza del piano di simmetria
piano di simmetria Piano di simmetria
Y
X
Z
VINCOLI SUI NODI
Carichi
simmetrici
Uz=0
ROTx=0
ROTy=0
Carichi
antisimm.
Uy=0
Ux=0
ROTz=0
+ Considerazioni di simmetria (5/5) Corpi assial-‐simmetrici l Geometria assial-‐simmetrica (rotazione di una sezione a>orno ad un asse fisso) l Carichi a simmetria cilindrica • Fissato un sistema di riferimento cilindrico “r, θ, z”, per simmetria lo stato di tensione/deformazione risulta indipendente da θ e le componen4 di spostamento in direzione circonferenziale (θ) risultano nulle: il problema può di conseguenza essere studiato come piano. l
+ Esercizio 2 P • Lastra intagliata in trazione – Schema4zzare la lastra di figura sfru>ando i piani di simmetria – Misure in mm – Spessore: 5 mm – Modulo Elas4co 109 Pa – P = 3000 Pa 20 5 60 + Link u4li • h>p://www.uniroma2.it/didalca/
Calc_Aut_Sis_Mec/deposito/08-‐Elemen4-‐
Formulazione-‐Generale_V1.pdf • h>p://www.aero.polimi.it/~ls075775/
bacheca/032FormulazioneFEM.pdf