METODA ELEMENTELOR FINITE L 7 Catedra de Constructii

METODA ELEMENTELOR FINITE – L 7 Catedra de Constructii Hidrotehnice CONDENSARE STATICA Procedeul de eliminare a unui nod interior al discretizarii prin reorganizarea sistemului global de ecuatii. Pentru exemplul de fata, are ca rezultat inlocuirea elementelor triunghilare 1 … 4 cu patrulaterul, prin eliminarea nodului interior 5. Partitionarea sistemului global de ecuatii: 3 y, v 2 Rezolvarea in i : 3 4 Substituirea: 2 5 1 1 4 x, u 1

METODA ELEMENTELOR FINITE – L 7 Catedra de Constructii Hidrotehnice - sistemul de ecuatii condensat - matricea de rigiditate condensata kc - vectorul incarcarilor condensate rc Nota: in mod obisnuit nu se asociaza conditii de margine (incarcari concentrate sau GDL cu valori impuse) nodurilor interioare ce urmeaza a fi condensate, astfel incat rc = re 2

METODA ELEMENTELOR FINITE – L 7 Catedra de Constructii Hidrotehnice SUBSTRUCTURARE - MACROELEMENTE Substructurarea este procedeul de formare a macroelementelor. Un macroelement este un element finit complex rezultat din asocierea unui numar de elemente finite obisnuite, prin condensare statica. Macroelementul va fi stocat in baza de date cu matrici caracteristice proprii, din care sunt eliminate gradele de libertate interioara. Macroelementele pot fi creeate pentru acele regiuni ale modelului care nu prezinta comportare neliniara, astfel incat matricile caracteristice nu necesita o reevaluare la fiecare iteratie. 3

METODA ELEMENTELOR FINITE – L 7 e – GDL exterioare Catedra de Constructii Hidrotehnice detaliu i – GDL interioare Macroelement - reducerea timpului de procesare pentru structuri (domenii) care contin paternuiri geometrice repetate, in special in cazul analizelor neliniare (se aplica regiunilor care raman in stadiul de comportare liniara); - rezolvarea modelelor mari cu resurse de calcul limitate 4

METODA ELEMENTELOR FINITE – L 7 Catedra de Constructii Hidrotehnice Scopul este de a mentine numai GDL exterioare pentru faza de procesare (solutie) prin eliminarea GDL interioare. Minimizarea functionalei EE la nivelul macroelementului Prin substituirea GDL interioara Procedeul de asamblare a matricii de rigiditate globale si a vectorului incarcarii iau in considerare numai variabilele asociate conturului macroelementului (asociate nodurilor de pe conturul acestuia). 5

METODA ELEMENTELOR FINITE – L 7 Catedra de Constructii Hidrotehnice Energia potentiala totala a sistemului se separa dupa cum urmeaza - n – numarul elementelor finite obisnuite din discretizare - m – numarul de macroelemente di discretizare Sistemul de ecuatii redus: Prin rezolvarea acestuia se obtin valorile GDL exterioare e. Revenind la primul set de ecuatii obtinut prin partitionare Reasambland E = e + i 6

METODA ELEMENTELOR FINITE – L 7 Catedra de Constructii Hidrotehnice Aplicare practica: o analiza prin substructurarea modelului implica 3 pasi distincti: - definirea (definition step); - solutia (solution step); - detalierea (expansion step). Definirea: elementele obisnuite sunt condensate in cadrul macroelementului, prin alegerea unor GDL master care au rolul de a defini interfata intre macroelement si zonele adiacente. Singura restrictie se refera la comportarea liniara a acestora. Alegerea GLD master: - trebuie sa permita conectarea tuturor elementelor finite obisnuite sau a altor macroelemente pe conturul macroelementului; - trebuie sa corespunda tuturor locatiilor (pozitiilor) conditiilor de margine (restrictiilor) si ale incarcarilor concentrate; - daca macroelementul este utilizat intr-o analiza dinamica, ele caracterizeaza comportarea dinamica a modelului (mase concentrate in aceste locatii). 7

METODA ELEMENTELOR FINITE – L 7 Catedra de Constructii Hidrotehnice Solutia: se efectueaza utilizand modelul asamblat din elemente obisnuite si macroelemente. Rezultatul consta in valorile GLD corespunzatoare nodurilor elementelor obisnuite si GLD master pentru macroelemente. Deoarece nu sunt disponibile rezultate privind GDL interioare, solutia este denumita solutia redusa. Detalierea: se calculeaza GDL interioare. Pasul se aplica independent tuturor macroelementelor pentru care se doresc rezultate interioare. 8

METODA ELEMENTELOR FINITE – L 7 Catedra de Constructii Hidrotehnice Substructurarea de jos in sus (bottom-up) (potrivita pentru modele mari): Fiecare macroelement este creeat separat intr-o etapa de generare individuala, iar macroelementele sunt asamblate in faza de solutie. Substructurarea de sus in jos (top-down) (potrivita pentru modele mici): - se construieste intregul model din elemente finite obisnuite, incluzand toate zonele (cu si fara viitoare macroelemente); - se genereaza macroelemente prin selectarea convenabila a elementelor finite obisnuite din discretizare si asocierea lor; - se executa solutia; - se executa detalierea separat pentru fiecare macroelement care przinta interes. 9

METODA ELEMENTELOR FINITE – L 7 Catedra de Constructii Hidrotehnice ELEMENTE FINITE PENTRU PLACI Placile sunt solide tridimensionale adesea utilizate in ingineria civila (diafragme, plansee, rezervoare, acoperisuri, etc). Placi plane (plates) Placi curbe (shells) Principala caracteristica: grosimea redusa t in raport cu suprafata. Abordari posibile: - utilizand elemente de clasa C 1, bazate pe teoria placilor subtiri (vezi curs T. E. ); in acest caz se defineste numai planul median al placii si sunt indeplinite o serie de ipoteze simplificatoare; - considerand placile drept cazuri speciale de solide 3 D, utilizand elemente patratice de clasa C 0. 10

METODA ELEMENTELOR FINITE – L 7 Catedra de Constructii Hidrotehnice Ipoteze de baza pentru placi plane si curbe subtiri (elemente de clasa C 1) : - eforturile unitare normale pe planul median ( z) sunt suficient de mici incat sa fie neglijate; - normala n la planul median ramane dreapta si normala pe suprafata placii in timpul si dupa deformarea acesteia; - deformatiile specifice si eforturile unitare au o distributie liniara pe grosimea placii; - deformatiile (sagetile) sunt mici in comparatie cu dimensiunile placii. In teoria placilor subtiri, deplasarile in planul placii (u, v) sunt neglijate in mod obisnuit, forma deformata fiind definita complet de deplasarea normala w (x, y). Pentru placile supuse la incovoiere pura, eforturile din planul placii sunt si ele neglijate (efect de membrana suprimat). 11

METODA ELEMENTELOR FINITE – L 7 m m xy Catedra de Constructii Hidrotehnice xy m x m y m xy m y x xy - mx si my momente incovoietoare pe unitatea de lungime - mxy moment de rasucire pe unitatea de lungime - q incarcarea distribuita pe unitatea de suprafata In formularea generala, deformatiile specifice sunt impartite in deormatiile specifice din planul placii, datorate incovoierii ( x, y, xy) si deformatiile specifice datorate forfecarii transversale ( xz, yz). 12

METODA ELEMENTELOR FINITE – L 7 Catedra de Constructii Hidrotehnice Relatii de baza: 13

METODA ELEMENTELOR FINITE – L 7 Catedra de Constructii Hidrotehnice Prin introducerea conceptului de deformatie specifica generalizata Relatia efort unitar – deformatie specifica devine: Tinand seama de distributia liniara a eforturilor si deformatiilor pe grosimea placii, rezulta: 14

METODA ELEMENTELOR FINITE – L 7 Catedra de Constructii Hidrotehnice Pentru o placa omogena si izotropa, de grosime constanta, conditia de echilibru se reduce la ecuatia bi-armonica: Principiul variational se poate exprima prin minimizarea energiei potentiale totale asociate efectului de incovoiere: In cazul placilor subtiri (placi Kirchoff), forfecare transversala se neglijeaza, = 0, rezulta: Functionala contine derivate de ordinul II ale functiei necunoscute w(x, y), fiind astfel necesara asigurarea clasei de continuitate C 1 pentru functiile de aproximare. 15

METODA ELEMENTELOR FINITE – L 7 Catedra de Constructii Hidrotehnice Un element de placa preia incarcarea prin ambele efecte – incovoiere si membrana (efectul de membrana introduce eforturi in planul placii datorate deplasarilor u si v). Energia potentiala totala este: - primul set de termeni exprima energia de deformatie datorata eforturilor si deformatiilor in planul placii; - al doilea set de termeni exprima energia de deformatie datorata deformatiilor si eforturilor generalizate (curbura si moment incovoietor unitar). 16

METODA ELEMENTELOR FINITE – L 7 Catedra de Constructii Hidrotehnice GDL sunt u si v pentru efectul de membrana si w, x, y pentru efectul de incovoiere. Solutia in elemente finite conduce la urmatoarea expresie a energiei potentiale totale: Expresia pune in evidenta faptul ca cele doua efecte sunt independente (efectul de incovoiere nu afecteaza efectul de membrana si invers). 17

METODA ELEMENTELOR FINITE – L 7 Catedra de Constructii Hidrotehnice Elemente de clasa C 0 – solide 3 D Datorita formei geometrice curbe, cel mai potrivit este un element finit isoparametric patratic - grosime t relativ redusa fata de celelalte dimensiuni ale planului median: sistemul de ecuatii global prost conditionat; - numar mare de GDL (20 3 = 60 GDL): creste memoria de stocare necesara si timpul de calcul; - nodul median plasat pe grosimea elementului violeaza ipoteza normalei drepte la suprafata mediana a elementului. 18

METODA ELEMENTELOR FINITE – L 7 Catedra de Constructii Hidrotehnice ELEMENTE DE CONTACT Din varii motive, solutii aflate initial in contact devin separate prin rosturi sau prezinta lunecari relative in lungul suprafetelor de contact. Un model detaliat trebuie uneori sa simuleze aceste efecte, pentru a pune in evidenta comportarea reala a structurii. Rost permanent Nivel coronament nivel fundatie plot i plot j plot k Starea rostului (inchisa sau deschisa, cu lunecare sau gripata) depinde de nivelul apei in amonte si de distributia campului termic, fiind variabila in timp. 19

METODA ELEMENTELOR FINITE – L 7 Catedra de Constructii Hidrotehnice La simularea unei excavatii sprijinite, structura de sprijin sufera deformatii in functie de nivelul excavatiei si de pozitia spraiturilor desprindere Nivel teren H 1 Uneori, contactul intre masiv si structura dispare pe anumite portiuni. forta sprait Strat coeziv Strat necoeziv H 2 Nivel excavatie In aceste zone, tensiuile intre masiv si structura trebuie sa dispara. 20

METODA ELEMENTELOR FINITE – L 7 Catedra de Constructii Hidrotehnice CONTACT NOD - NOD Permite mentinerea sau anularea contactului fizic intre doua suprafete, sau lunecarea lor relativa. z y x Suporta compresiune pe directie normala la suprafata de contact si forfecare in lungul directiei tangentiale. Orientarea suprafetei de contact (sau lunecare) este data de pozitia initiala a nodurilor care definesc elementul de contact (normala la directia i - j). 21

METODA ELEMENTELOR FINITE – L 7 Catedra de Constructii Hidrotehnice In conformitate cu modelul de frecare Coulomb, doua suprafate in contact pot suporta eforturi de forfecare relativa pana la un anumit nivel. Odata ce eforturile de forfecare echivalente depasesc aceasta limita, cele doua suprafete vor manifesta lunecare relativa Comportarea elementului depinde de starea lui curenta: - inchis si gripat; - inchis cu lunecare; - deschis. Rostul initial, definit in mod obisnuit de pozitia nodurilor i si j, se schimba odata cu depaasarile nodale calculate. 22

METODA ELEMENTELOR FINITE – L 7 Catedra de Constructii Hidrotehnice CONTACT NOD - SUPRAFATA Atunci cand un contact precis nod- nod nu se poate realiza. Contactul apare de fiecare data cand nodul de contact k penetreaza suprafata tinta i - j 23