3658x
000434
2024-01-26
Struttura

Elementi finiti in RFEM

Elementi 1D

Per gli elementi dell'asta, si presume che la sezione trasversale rimanga piana durante la deformazione. Gli "elementi 1D" sono utilizzati per rappresentare travi, travature reticolari, nervature, funi e giunti rigidi. Un elemento dell'asta 1D ha un totale di dodici gradi di libertà - sei all'inizio e sei alla fine dell'elemento, ciascuno per gli spostamenti (ux, uy, uz ) e rotazioni (φx, φy, φz ). Quando si calcolano i dati strutturali linearmente, trazione, compressione e torsione sono espresse come funzioni lineari dell'asse x dell'asta, indipendentemente dalla flessione e dal taglio. Queste funzioni sono approssimate da un polinomio del terzo ordine in x, inclusa l'influenza delle tensioni tangenziali dalle forze di taglio Vy e Vz. La matrice di rigidezza KL (12, 12) descrive il comportamento lineare degli elementi 1D. L'interazione reciproca tra forza assiale e flessione in caso di problemi geometricamente non lineari è espressa nella matrice di rigidezza KNL (12, 12). Per ulteriori informazioni, vedere [1] e [2].

Per i calcoli secondo l'analisi a grandi spostamenti, si consiglia di utilizzare un Infittimento della mesh della linea in modo che i risultati dell'asta possano essere determinati con la precisione corrispondente.

Elementi 2D

Per le superfici, vengono creati elementi 2D. Laddove non è possibile utilizzare elementi quadrangolari, il generatore di mesh inserisce elementi triangolari.

I gradi di libertà degli elementi quadrangolari e triangolari nei nodi d'angolo sono gli stessi degli elementi 1D: gradi di libertà di spostamento (ux, uy, uz ) e di rotazione (φx, φy, φz ). Ciò garantisce la compatibilità tra gli elementi 1D e 2D nei nodi. I parametri sono definiti nel sistema di coordinate locali planari degli elementi e sono convertiti nel sistema di coordinate globali quando si compone la matrice di rigidezza globale.

L'immagine Elementi quadrangolari per gusci mostra l'approccio per elementi gusci piani. Si basa sulla teoria di Mindlin/Reissner. Per garantire l'accoppiamento diretto con gli elementi dell'asta, viene scelto un approccio quadrato nel piano del guscio (ux, uy ). Eliminando i nodi intermedi, viene creato un elemento a quattro nodi con gradi di libertà rotazionali aggiuntivi nelle direzioni x, y e z. Pertanto, per gli elementi a diaframma è possibile un accoppiamento diretto con elementi a trave. Basato su un'interpolazione combinata di spostamenti trasversali, rotazioni della sezione trasversale e deformazioni tangenziali trasversali, i metodi di Bathe e Dvorkin elementi MITC presentati da efer [3] ( Mixed Interpolation of Tensorial Components ) sono stati utilizzati: MITC3+ per triangoli, MITC4 per quadrangoli.

Al momento, gli elementi asta sono considerati risolvendo direttamente l'equazione differenziale secondo l'analisi del secondo ordine. Gli effetti di ingobbamento non possono essere considerati quando si utilizza la torsione di Saint Venant.

I calcoli delle membrane si basano sui principi di Bergan [4] [5] [6]. Per gli elementi triangolari, ad esempio, le funzioni di base sono divise in tre deformazioni del corpo rigido, tre condizioni di deformazione costante e tre speciali gradienti lineari di tensione e deformazione. Il campo di deformazione all'interno di un elemento è quadratico e il campo di tensione è lineare. La matrice di rigidezza dell'elemento KL è successivamente trasformata in nove parametri collettivi dei tipi ux, uy, φz. I componenti di questa matrice sono inseriti insieme ai componenti per flessione e taglio nella matrice di rigidezza complessiva (18, 18). Questa matrice è il risultato del concetto di Lynn/Dhillon. In questo modo, vengono applicate le cosiddette piastre Mindlin (piastre con deformazione a taglio distinta) e analizzate secondo Timoshenko. In questo modo, RFEM trova la soluzione corretta sia per le piastre spesse che per quelle sottili (piastre di Navier).

Nel caso di problemi geometricamente non lineari, la suddetta divisione della condizione tensione-deformazione in uno stato planare e in flessione/taglio non è possibile. Le influenze reciproche di questi stati sono considerate nella matrice KNL. RFEM utilizza una forma piuttosto semplice ma efficace della matrice KNL che si basa sugli approcci di Zienkiewicz [7]. Viene applicata la componente quadrata ε2 del tensore di deformazione di Green/Lagrange ε = ε1 + ε2. Si assume una distribuzione lineare di uz (x, y) della condizione di tensione planare e distribuzioni lineari di ux (x, y) e uy (x, y) dell'interazione con la flessione. Questa ipotesi è possibile, perché l'effetto principale dell'interazione dipende dalla prima derivazione dell'equazione differenziale e perché l'influenza delle componenti di ordine superiore diminuisce rapidamente con la divisione in elementi più piccoli. La correttezza di questa procedura è stata dimostrata in diverse analisi numeriche.

Info

Lo spessore degli elementi per gli elementi shell deve essere significativamente inferiore alle dimensioni. In caso contrario, si consiglia di modellare gli oggetti come solidi. Inoltre, per quanto riguarda gli elementi a guscio, evitare di introdurre tensioni torsionali selettivamente: Il grado di libertà rotazionale intorno alla normale alla superficie reagisce in modo molto sensibile.

Elementi 3D

Gli elementi 3D sono utilizzati per i solidi. I seguenti tipi di elementi sono implementati in RFEM: tetraedro, pentaedro (prisma, piramide) ed esaedro. Per informazioni dettagliate sugli elementi applicati e sulle matrici, vedere Sevcik's 3D Elements with Rotational Degrees of Freedoms [8] per i dettagli. È possibile richiedere questa documentazione a Dlubal Software.

In generale, tutti i gradi di libertà rotazionali devono essere considerati critici per i solidi. Poiché la deformazione di un solido è determinata esclusivamente dai vettori di spostamento, la rotazione di un nodo della mesh (ad esempio, a causa della torsione introdotta singolarmente) non influenza la deformazione all'interno del solido.


Bibliografia
  1. Vladimír Kolář et al. Bemessung von zwei- und dreidimensionalen Strukturen mit FEM. Springer-Verlag, New York / Wien, 1975. Kapitel 1 (1D-Element) und 6 (Variationsprinzip)
  2. Vladimír Kolář und Ivan Němec. Finite Element Analysis of Structures. United Nations Development Program, Economic Com. for Europe, Workshop on CAD Techniques, Prague - Geneva, 1984.
  3. Eduardo N. Dvorkin und Klaus-Jürgen Bathe. A continuum mechanics based four-node shell element for general non-linear analysis. Engineering Computations, 1, 1984.
  4. P. G. P.rgan. P.nite Elements Based on Energy Orthogonal Functions. P.ternational Journal for Numerical Methods in Engineering, 15, 1980.
  5. P. G. P.rgan und M. P. P.gård. P.nite Elements With Increased Freedom in Choosing Shape Functions. P.ternational Journal for Numerical Methods in Engineering, 20, 1984.
  6. P. G. P.rgan und Carlos A. P.lippa. P.Triangular Membrane Element With Rotational Degrees of Freedom. P.mputer Methods in Applied Mechanics and Engineering, 50, 1985.
  7. Olgierd Cecil Zienkiewicz. The Finite Element Method in Engineering Science. Mc Graw-Hill, London 3. Auflage, 1979. Chapter 18 - 19 (Nonlinear Problems).
  8. I. Sevčík. I. Finite Elements with Rotational Degrees of Freedom. I.M Consulting s.r.o, Brno.
Sezione originaria