VERIFICAREA MODELELOR DE CALCUL CU ELEMENTE FINITE
VERIFICAREA MODELELOR DE CALCUL CU ELEMENTE FINITE
VERIFICAREA MODELELOR DE CALCUL CU ELEMENTE FINITE
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.