20.04.2014 Views

VERIFICAREA MODELELOR DE CALCUL CU ELEMENTE FINITE

VERIFICAREA MODELELOR DE CALCUL CU ELEMENTE FINITE

VERIFICAREA MODELELOR DE CALCUL CU ELEMENTE FINITE

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

14.<br />

<strong>VERIFICAREA</strong> <strong>MO<strong>DE</strong>LELOR</strong> <strong>DE</strong> <strong>CAL<strong>CU</strong>L</strong><br />

<strong>CU</strong> <strong>ELEMENTE</strong> <strong>FINITE</strong><br />

Modelul de calcul şi rezultatele obŃinute cu ajutorul său trebuie supuse unor numeroase teste şi<br />

verificări. Scopul acestora este de a “valida” modelul, adică de a determina dacă acesta satisface<br />

exigenŃele impuse şi dacă rezultatele obŃinute cu ajutorul lui permit formularea unor răspunsuri<br />

neechivoce la întrebările clare puse de beneficiarul analizei cu elemente finite (FEA). Unele teste şi<br />

verificări sunt calitative şi globale, altele cantitative şi de detaliu.<br />

Dacă testele şi verificările duc la concluzii nefavorabile, modelul trebuie îmbunătăŃit şi procesul de<br />

verificare – îmbunătăŃire - verificare se continuă până când se obŃine un model satisfăcător, adică<br />

valid. În figura 8.1 este prezentată schema generală a procesului de verificare îmbunătăŃire a modelului<br />

de calcul cu elemente finite. În continuare se prezintă câteva metode şi procedee de verificare.<br />

Figura 8.1<br />

Verificările experimentale efectuate pe structura reală sunt cele mai concludente. Astfel de verificări<br />

sunt însă, de obicei, ulterioare calculului (după ce s-a proiectat şi executat structura) şi totdeauna sunt<br />

costisitoare. O situaŃie specială apare când se expertizează structuri vechi pentru care nu s-au efectuat<br />

analize cu elemente finite la proiectare. Se pot face verificări experimentale şi pe modele fizice reduse<br />

la scară ale structurii reale, în astfel de cazuri fiind necesară rezolvarea problemelor de modelare şi<br />

similitudine.<br />

Determinările experimentale permit verificări cantitative ale rezultatelor obŃinute prin calcul şi<br />

evaluarea preciziei acestora. Se pot verifica mărimi care provin din comportarea globală a structurii,<br />

cum sunt deplasările sau reacŃiunile în reazeme, sau mărimi cu caracter local, cum sunt tensiunile<br />

maxime.<br />

Efectuarea calculelor pe două sau mai multe modele şi compararea rezultatelor obŃinute.<br />

Modelele pot fi de acelaşi tip, adică elaborate pe baza aceleiaşi metode de calcul (de exemplu, metoda<br />

elemntelor finite - MEF) sau de tipuri diferite, adică elaborate pe baza unor metode de calcul diferite<br />

(de exemplu MEF şi metoda elementelor de frontieră sau o metodă analitică de calcul ).


De exemplu, pentru traversa de egalizare a unui excavator care se deplasează pe o cale de rulare s-au<br />

elaborat trei modele MEF: un model de bară (fig. 8.2,a), un model de tip stare plană de tensiuni<br />

(fig. 8.2,b) şi un model spaŃial (fg. 8.2,c).<br />

Figura 8.2<br />

Cele trei modele pot fi folosite în etape distincte ale proiectarii, respectiv pentru predimensionare,<br />

definitivarea formei constructive şi verificarea finală a structurii obŃinute. Unele dintre rezultatele<br />

obŃinute în cele trei variante de analiză, ca de exemplu deplasările maxime ( care sunt mărimi globale<br />

ale FEA ) trebuie să aibă valori apropiate. Valorile tensiunilor maxime, în special ale celor locale, pot<br />

să difere destul de mult de la un model la altul.<br />

Uneori este preferabil să se folosească în paralel modele elaborate pe baza unor metode de calcul<br />

diferite, ca de exemplu, modele cu elemente finite şi modele cu elemente de frontieră. Pentru trei dinŃi<br />

consecutivi ai roŃii dinŃate din figura 8.3,a s-a utilizat în paralel, în vederea verificării, un model cu<br />

elemente finite (fig. 8.3,b) şi un altul cu elemente de frontieră (fig. 8.3,c). În acest caz toate rezultatele<br />

obŃinute în urma celor două analize (deplasări, tensiuni, reacŃiuni, etc.) trebuie să aibă valori foarte<br />

apropiate.<br />

Figura 8.3<br />

Preprocesarea geometriei modelului MEF este cea mai utilizată şi cea mai eficientă metodă de<br />

verificare a geomeriei modelului, a corectitudinii definirii condiŃiilor de rezemare şi a aplicării<br />

sarcinilor. Se poate spune că este totdeauna obligatorie. Verificarea constă în citirea fişierului cu datele<br />

de intrare pentru programul MEF, preprocesarea informaŃiilor conŃinute în acest fişier şi trasarea<br />

unui desen al modelului structurii. Un astfel de exemplu se prezintă în figura 8.4, pentru modelul MEF<br />

al unei structuri industriale.<br />

Preprocesarea se face, de regulă, în programul MEF, care are module de elaborare a modelului,<br />

inclusiv de desenare a acestuia, în diverse condiŃii grafice. Foarte utilă este această verificare pe<br />

parcursul elaborării modelului, în diversele etape ale procesului. În situaŃii deosebite preprocesarea<br />

se poate face şi cu alte programe, preferate sau chiar elaborate de utilizator.


Figura 8.4 Figura 8.5<br />

Verificări ale condiŃiilor de simetrie. Pentru modele care prezintă proprietăŃi de simetrie sau<br />

antisimetrie geometrică şi mecanică, rezultatele obŃinute prin calcul - deplasări, tensiuni, reacŃiuni în<br />

reazeme, moduri proprii de flambaj (simetrice sau antisimetrice), moduri proprii de vibraŃii etc. -<br />

trebuie să aibă valori egale în punctele simetrice şi valori egale şi semn schimbat în punctele<br />

antisimetrice. Pentru suportul din figura 8.5, care are două plane de simetrie – XOY şi YOZ – atât<br />

pentru configuraŃia geometrică precum şi pentru reazeme şi sarcini, verificarea corectitudinii<br />

modelului MEF se face comparând valorile deplasărilor, tensiunilor şi reacŃiunilor în puncte şi<br />

elemente simetrice: aceste valori trebuie să fie egale, mărimile abaterilor putând fi un indiciu cantitativ<br />

al preciziei FEA, în ansamblu, adică a metodei, a programului şi a modelului.<br />

Verificări printr-un calcul simplu. Uneori este posibil să se verifice unele dintre rezultatele obŃinute<br />

cu un anumit model de calcul, considerând un caz de încărcare simplificat, de exemplu o sarcină<br />

concentratră, un moment etc. De exemplu, pentru grinzile longitudinale ale utilajului din figura 8.6,<br />

Figura 8.6 Figura 8.7<br />

se poate face un calcul la încovoiere cu relaŃia cunoscută din rezistenŃa materialelor, considerând<br />

grinda din figura 8.7 cu secŃiune constantă, încărcată la mijloc cu o sarcină uniform distribuită,<br />

rezemată la capete. În acest caz se are în vedere că rigiditatea la încovoiere în plan vertical a celor două<br />

grinzi longitudnale este foarte mare, adică efectul corpului recipientului poate fi neglijat în calculul<br />

menŃionat. Valorile săgeŃii maxime şi cea a tensiunii la mijlocul grinzii trebuie să fie apropiate de cele<br />

obŃinute cu FEA.


Discretizarea adaptivă. Acest procedeu nu este propriu-zis o cale de verificare a corectitudinii<br />

modelului FEA, dar poate oferi informaŃii consistente în această privinŃă. Pe de altă parte modelarea<br />

adaptivă se face automat, de către programul FEA, dacă sunt îndeplinite anumite condiŃii, ceea ce este<br />

foarte comod pentru utilizator. Se elaborează un model iniŃial MEF care se supune unui proces FEA şi<br />

se obŃin rezultatele corespunzătore. În programul MEF (dacă acesta are implementată procedura<br />

respectivă) se dau comenzile specifice analizei adaptive care constă în elaborarea, pentru modelul<br />

inŃial - printr-un proces iterativ - a unei discretizări mai fine (procedura h), utilizarea unor elemente<br />

finite de ordin superior, adică cu polinoame de interpolare de grad superior (procedura p), sau<br />

combinaŃii ale acestora (procedura h-p). Rezultatele obŃinute pentru noua variantă a modelului sunt<br />

mai precise decăt cele iniŃiale. Programul calculeză indicele de precizie al modelului şi când valoarea<br />

prescrisă a acestuia este atinsă, procesul iterativ de “rafinare” a reŃelei de discretizare se opreşte.<br />

Verificarea modelului constă în compararea rezultatelor obŃinute în cele două variante ale modelului şi<br />

anume rezultatele trebuie să fie suficient de apropiate.<br />

Figura 8.8 Figura 8.9<br />

Pentru exemplificare, se prezintă în figura 8.8 modelul iniŃial MEF al unei plăci dreptunghiulare<br />

plane, discretizată cu 55 de noduri şi 80 de elemente shell triunghiulare. După aplicarea unei proceduri<br />

de “rafinare“ a discretizării (procedura h), modelul are 449 de noduri şi 798 elemente, ca în figura 8.9.<br />

Se menŃionează faptul că o procedură de discretizare adaptivă, oricare ar fi ea, nu poate semnala<br />

eventuale greşeli sau neajunsuri ale modelului iniŃial, cum ar fi, de exemplu, configuraŃia geometrică<br />

de ansamblu, valorile dimensiunilor, alegerea tipului de element finit, impunerea condiŃiilor de<br />

reazeme, definirea sarcinilor, introducerea valorilor constantelor elastice şi fizice ale materialului etc.<br />

Verificarea greutăŃii structurii este o verificare globală, obligatorie. Trebuie verificate valorile<br />

reacŃiunilor din reazeme şi dacă acestea satisfac ecuaŃiile de echilibru scrise pentru întreaga structură.<br />

Dacă este posibil, este bine să se verifice şi poziŃia centrului de greutate al structurii.<br />

Verificări globale şi calitative ale modelului au în vedere configuraŃiile stărilor de tensiuni şi<br />

deplasări, semnele lor, ordinul de mărime şi chiar valorile rezultatelor obŃinute. Din practica<br />

inginerească şi din experienŃa altor analize se ştie unde sunt zonele cu tensiuni şi deplasări mari, care<br />

Figura 8.10


este configuraŃia structurii deformate şi între ce limite trebuie să se afle valorile mărimilor obŃinute<br />

prin FEA. Pentru exemplificare se prezintă în figura 8.10 proiecŃia în plan orizontal a structurii din<br />

figura 8.6 deformată, pentru un sistem de sarcini simetrice.<br />

Comentarii, observaŃii, concluzii:<br />

Prezentarea, mai sus, a unor modalităŃi de verificare a modelelor MEF, desigur că nu este<br />

exhaustivă. Fiecare utilizator poate să-şi imagineze şi alte tehnici şi metode de verificare.<br />

Se poate remarca faptul că în lucrarea de faŃă nu s-au făcut nici un fel de menŃiuni cantitative privind<br />

condiŃiile de precizie cu care se pot sau trebuie verificate modele MEF. Aceasta nu este o pierdere din<br />

vedere ci este rezultatul unei realităŃi şi anume că nu se poate stabili precizia unui model MEF în<br />

general, ci toate exigenŃele impuse modelului (inclusiv cele de precizie) depind de particularităŃile<br />

concrete, asociate problemei inginereşti care se rezolvă, ca, de exemplu, tipul structurii, scopul FEA,<br />

importanŃa structurii, gradul de pericol în cazul unei avarii, tipul solicitării, durata de exploatare etc.<br />

De asemenea, trebuie avut în vedere faptul că MEF este aproximativă, ceea ce înseamnă că nu se<br />

poate cere modelului mai mult decât poate oferi metoda, rezultatele obŃinute fiind determinate atât de<br />

performanŃele modelului cât şi de principiile, ipotezele şi procedurile matematice de calcul incluse în<br />

metoda şi în programul cu elemente finite.<br />

Toate verificările făcute modelelor cu elemente finite sunt validate în ultimă instanŃă de<br />

intuiŃia şi experienŃa utilizatorului şi atunci când este posibil, experimental.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!