Cercetare de Ex - Institutul de Mecanica Solidelor
Cercetare de Ex - Institutul de Mecanica Solidelor
Cercetare de Ex - Institutul de Mecanica Solidelor
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
Testul <strong>de</strong> in<strong>de</strong>ntare uniaxială constă în apăsarea verticală a unui in<strong>de</strong>ntor rigid pe<br />
semispaţiul materialului ce se doreşte a fi caracterizat (fig.3a). Nanoin<strong>de</strong>terul înregistrează atât<br />
forţa <strong>de</strong> apăsare P cât şi adâncimea <strong>de</strong> pătrun<strong>de</strong>re hvarf a capului sferic al in<strong>de</strong>nterului în<br />
semispaţiul materialului. Capul sferic al in<strong>de</strong>nterului nu este complet rigid, fiind alcătuit dintr-un<br />
material hipoelastic având modulul Young E varf = 1016 GPa şi coeficientul Poisson ν varf = 0.07 .<br />
a)<br />
Fig. 3: a) Schema testului <strong>de</strong> in<strong>de</strong>ntare uniaxială; b) Mo<strong>de</strong>lare axisimetrică cu elemente finite a testului <strong>de</strong> in<strong>de</strong>ntare.<br />
Cu ajutorul unui program <strong>de</strong> element finit pentru <strong>de</strong>formaţii elasto-plastice finite, s-a reuşit<br />
simularea in<strong>de</strong>ntării unui semispaţiu izotrop alcătuit dintr-un material inelastic cu un in<strong>de</strong>nter<br />
având cap sferic. Acest program <strong>de</strong> element finit este intitulat SPPRc (SPrangPRozessioun fir<br />
contact-Mo<strong>de</strong>ller), programul fiind realizat <strong>de</strong> colaboratorul nostru Dr.-ing. Gaston Rauchs, <strong>de</strong> la<br />
Laboratoire <strong>de</strong> Technologies Industrielles, Centre <strong>de</strong> Recherche Public Henri Tudor, Luxemburg.<br />
Scris în limbaj Fortran, programul SPPRc mo<strong>de</strong>lează prin metoda elementelor finite contactul<br />
dintre capul sferic al in<strong>de</strong>nterului şi materialul inelastic. Mo<strong>de</strong>larea axisimetrică cu elemente<br />
finite a capului sferic al in<strong>de</strong>nterului şi a semispaţiului ce constituie materialul este prezentată în<br />
fig.3b, fiind folosite 112 elemente finite patrulatere pentru a mo<strong>de</strong>la jumătate din semispaţiul<br />
materialului in<strong>de</strong>ntat şi 47 elemente finite patrulatere pentru a mo<strong>de</strong>la un sfert din capul sferic al<br />
in<strong>de</strong>nterului. În privinţa contactului dintre capul sferic al in<strong>de</strong>nterului şi materialul in<strong>de</strong>ntat,<br />
mo<strong>de</strong>larea acestuia se face împiedicând întrepătrun<strong>de</strong>rea geometrică dintre cele două corpuri în<br />
contact, fapt ce se realizează prin aplicarea unor tracţiuni pe suprafaţa <strong>de</strong> contact. Aceste tracţiuni<br />
introduse <strong>de</strong> mo<strong>de</strong>larea contactului se calculează pe baza distanţei locale dintre nodurile<br />
semispaţiului materialului şi proiecţia acestor noduri pe suprafaţa in<strong>de</strong>nterului (Rauchs 2006).<br />
În programul <strong>de</strong> element finit SPPRc, ecuaţiile constitutive ale materialului inelastic sunt<br />
formulate incremental, sub forma unor legături diferenţiale între <strong>de</strong>formaţii şi tensiuni (Rauchs<br />
2006). Comportarea plastică a materialului in<strong>de</strong>ntat este mo<strong>de</strong>lată folosind teoria fluxului J2<br />
izotropic, curgerea plastică fiind <strong>de</strong>terminată <strong>de</strong> funcţia <strong>de</strong> curgere f :<br />
⎧ f < 0 pentru comportare elastica,<br />
y ⎪<br />
f = P :( σ −α) −K⎨f=<br />
0 pentru comportare plastica,<br />
⎪<br />
⎩ f > 0 exclus,<br />
y un<strong>de</strong> K este limita <strong>de</strong> curgere, α este aşa-numitul “back-stress”, o variabilă internă <strong>de</strong> răspuns a<br />
S<br />
materialului, σ este tensorul Cauchy al tensiunilor şi P= I −1I⊗I este operatorul <strong>de</strong>viatoric <strong>de</strong><br />
3<br />
proiecţie, I S fiind tensorul unitate simetric <strong>de</strong> ordinul 4. Variaţia limitei <strong>de</strong> curgere y K este<br />
<strong>de</strong>scrisă consi<strong>de</strong>rând ecruisarea izotropică:<br />
11