ELEMENTE DE DINAMICA PĂMÂNTURILOR ELEMENTS OF SOILS DYNAMICS Lucrare sponsorizată de Consiliul Naţional al Cercetării Ştiinţifice din Învăţământul Superior în cadrul grantului nr. 33517/2002 – tema 1 DINU BRATOSIN ELEMENTE DE DINAMICA PĂMÂNTURILOR EDITURA ACADEMIEI ROMÂNE Bucureşti, 2002 Copyright © Editura Academiei Române, 2002. Toate drepturile asupra acestei ediţii sunt rezervate editurii. Adresa: EDITURA ACADEMIEI ROMÂNE Calea 13 Septembrie nr. 13, sector 5, Bucureşti, România, Tel: 4021–411 9008; 4021–410 3200 Fax: 4021–410 3983 e-mail: edacad@ear.ro web: www.ear.ro Redactor: PETRE MOCANU Tehnoredactor: MAGDALENA JINDICEANU Coperta: GIGI GAVRILĂ Bun de tipar: 5.10.2002 Format: 16/70×100 Coli de tipar: 12 CZ pentru biblioteci mari: 526 551 CZ pentru biblioteci mici: 52 ISBN: 973-27-0933-2 Imprimat în România Tipografia: STELIAN INVEST SRL Tel./Fax: +40–21–250 22 81 CUPRINS Prefaţă …………………………………………………………………….. 9 1. Noţiuni generale ………………………………………………….. 11 1.1. Caracteristici mecanice definitorii …………………………. 1.1.1. Neliniaritatea ………………………………………… 1.1.2. Efectele reologice ……………………………………. 1.1.3. Ireversibilitatea ………………………………………. 11 12 14 15 1.2. Geometria spaţiului invarianţilor …………………………... 1.2.1. Invarianţi ……………………………………………… 1.2.2. Geometria spaţiului invarianţilor ....………………….. 17 17 19 2. Sistem oscilant cu un grad de libertate dinamică ……………... 22 2.1. Sistem liniar ………………………………………………... 2.1.1. Vibraţia liberă ……………………………………….. 2.1.2. Regim staţionar armonic …………………………….. 22 23 27 2.2. Sistem neliniar ……………………………………………... 29 2.3. Soluţie numerică pentru oscilatorul neliniar ……………….. 2.3.1. Metoda Newmark ……………………………………. 2.3.2. Program FORTRAN ………………………………… 2.3.3. Aplicaţii ……………………………………………… 31 31 33 35 3. Teste de definire şi cuantificare ………………………………… 43 3.1. Introducere …………………………………………………. 43 3.2. Teste de conpresiune triaxială ……………………………… 44 3.3. Teste triaxiale de fluaj ……………………………………… 46 3.4. Modelarea testului din coloana rezonantă …………………. 3.4.1. Coloana rezonantă – descriere sumară ………………. 3.4.2.Torsiunea probei ……………………………………... 3.4.3. Determinarea amplitudinii rotaţiei şi deformaţiei torsionale ……………………………………………... 3.4.4. Determinarea modulului torsional …………………... 3.4.5. Estimarea amplitudinii momentului perturbator …….. 3.4.6. Determinarea raportului de amortizare torsională …... 47 47 49 51 52 54 54 4. Relaţii constitutive elastic neliniare …………………………….. 58 4.1. Forma generală a ecuaţiei constitutive ……………………... 58 4.2. Relaţii constitutive între invarianţi …………………………. 59 4.3. Funcţii-modul ………………………………………………. 62 5. Modelarea efectelor reologice în medii liniare …………………. 65 5.1. Fluaj şi relaxare …………………………………………….. 65 5.2. Răspunsul la vibraţii forţate ………………………………... 67 5.3. Legi de comportare liniare …………………………………. 69 5.4. Modele diferenţiale de ordinul întâi ………………………... 5.4.1 Modelul Kelvin-Voigt ……………………………….. 5.4.2 Modelul Maxwell ……………………………………. 5.4.3 Solidul vâscoelastic standard ………………………… 70 70 71 73 5.5. Modele diferenţiale generalizate .…………………………… 76 5.6. Modele vâscoelastice de tip integral ……………………….. 77 5.7. Răspunsul în regim armonic staţionar ……………………… 79 6. Model vâscoelastic neliniar ……………………………………… 83 6.1. Forma generală a ecuaţiei constitutive ……………………... 83 6.2. Funcţiile neliniare de relaxare ……………………………… 84 7. Modelarea răspunsului dinamic ………………………………… 89 7.1. Model dinamic vâscoelastic neliniar ……………………….. 89 7.2. Model dinamic pentru solicitări seismice ………………….. 97 7.3. Model Kelvin-Voigt neliniar …………………………..…... 100 7.4. Validarea modelului vâscoelastic neliniar …………………. 106 8. Model vâscoplastic ………………………………………………. 109 8.1. Forma generală a ecuaţiei constitutive ……………………... 109 8.2. Determinare prin încercări triaxiale ………………………... 109 8.3. Determinare prin încercări dinamice ……………………….. 114 8.4. Comparaţie între răspunsul lent şi cel dinamic …………….. 117 9. Echivalenţa liniară în dinamica materialelor neliniare ……….. 119 9.1. Introducere …………………………………………………. 119 9.2. Metode de liniarizare ………………………………………. 9.2.1. Minimizarea diferenţei dintre soluţii ………………… 9.2.2. Minimizarea diferenţelor dintre forţele elastice şi de amortizare ……………………………………………. 9.2.3. Liniarizarea globală …………………………………... 120 120 122 122 9.3. Liniarizarea globală în regim de rezonanţă ………………… 124 10. Ruperea dinamică la pământuri şi roci ………………………… 129 10.1. Introducere ……………………………………………….. 129 10.2. Ruperea în condiţii statice ………………………………... 10.2.1. Criterii de rupere ………………………………….. 10.2.2. Un nou criteriu tri-invariant ………………………. 130 130 135 10.3. Efectul proprietăţilor reologice asupra ruperii …………… 138 10.4. Determinarea ruperii prin oboseală la roci ……………….. 10.4.1. Corelaţia rezistenţă de rupere-densitate ………..… 10.4.2. Curba tip Wöhler …………………………………. 140 141 143 10.5. Efectul de degradare dinamică …………………………… 145 10.6. Ruperea în condiţii seismice ……………………………… 154 10.7. Lichefierea. Caz particular de rupere dinamică ……….….. 158 11. Model vâscoelastic neliniar pentru materiale antivibratorii ….. 161 11.1. Introducere ………………………………………………... 161 11.2. Modelarea răspunsului neliniar …………………………... 162 11.3. Utilizarea relaţiilor experimentale forţă-deplasare …...…... 11.3.1. Funcţiile elastice şi de amortizare ………………... 11.3.2. Detrminarea răspunsului neliniar ………………… 11.3.3. Liniarizarea ………………………………………. 167 169 171 172 12. Evaluarea efectelor structurale …………………………………. 174 12.1. Introducere ………………………………………………... 174 12.2. Funcţii de material ………………………………………... 175 12.3. Răspunsul structural ……………………………………… 178 Bibliografie ………………………………………………………. 181 Contents ………………………………………………………….. 187 PREFAŢĂ Pentru o ţară cu un potenţial seismic apreciabil, studiile de comportare di-namică a materialelor din amplasament sunt deosebit de necesare şi încă actuale. Reducerea riscului seismic este condiţionată şi de elucidarea unor numeroase as-pecte teoretice, dar cu un deosebit impact practic. Un astfel de aspect este şi de-terminarea răspunsului dinamic al sistemelor care includ materiale degradabile şi disipative cu mari neliniarităţi, cum sunt pământurile şi rocile din amplasamentele construcţiilor. Pentru rezolvarea optimă a acestei probleme sunt necesare clarificarea a două aspecte esenţiale - determinarea pentru aceste materiale cu compor-tare dinamică complexă a unor legi constitutive cât mai complete şi apropiate realităţii iar apoi, asigurarea unor posibilităţi de utilizare a acestor legi în calculul răspunsului dinamic structural. Realitatea impune în mod sigur ca legile constitutive să fie legi neliniare. Avem la această oră posibilitatea de a modela legi neliniare şi avem (parţial, din păcate) logistica necesară determinării parametrilor de material. Însă, legile constitutive neliniare introduc dificultăţi în calculul răspunsului structural, dificultăţi care pot fi totuşi depăşite apelând la calcu-lul numeric. Monografia de faţă este axată pe probleme de modelare a răspunsului dinamic al materialelor disipative şi degradabile cu mari neliniarităţi dar abordează şi posibilităţile de depăşire a dificultăţilor introduse de această modelare în calculul răspunsului structural. Însă, nu a fost în intenţia autoru-lui tratarea completă a acestor aspecte. Autorul ţine numai să aducă în discu-ţia celor interesaţi unele rezultate ale eforturilor sale din ultimii ani în încer-carea de a obţine unele modele constitutive şi metode, ca un compromis optim între completitudine şi aplicabilitate. Structura lucrării reflectă această intenţie. Sunt prezentate în principal „elemente” de dinamică axate pe rezultatele amintite mai sus dar cadrul general care include aceste elemente este redus la un minim necesar parcurgerii textului. Un prim capitol de „Noţiuni generale” pezintă câteva consideraţii privind principalele caracteristici mecanice care definesc comportarea dinamică a pămân-turilor. De asemenea, tot în acest capitol a fost inclus şi un „remember” al spaţiu-lui invarianţilor tensorului tensiunii deoarece pe parcursul lucrării toate legile constitutive sunt concepute ca relaţii între invarianţii tensorilor tensiunii şi cei ai tensorului deformaţiei iar ruperea, statică sau dinamică, este legată de spaţiul in-varianţilor tensorului tensiunii. Din vastul material al dinamicii sistemelor elastice a fost extras pentru capi-tolul al doilea oscilatorul cu un grad de libertate dinamică care ulterior a fost folosit pentru definirea unui oscilator neliniar aplicat în modelele dinamice prezentate. Tot aici este inclusă prezentarea unei metode de obţinere a soluţiei neliniare, metodă folosită, de asemenea, ulterior. În mecanică, în general, dar mai ales în mecanica pământurilor, experimen-tul are un rol predominant. Nu pot fi concepute, iar dacă se concep nu pot fi folo-site, legi constitutive fără o susţinere experimentală. Din acest motiv, capitolul al treilea conţine o succintă prezentare a unei minime aparaturi necesare definirii şi cuantificării relaţiilor constitutive expuse pe parcursul lucrării. Comportarea dinamică a unui material nu este complet distinctă faţă de răs-punsul la solicitări statice. Din acest motiv, în capitolele 4, 5 şi 6 sunt succint ex-puse modele elstic şi vâscoelastic, liniare şi neliniare, modele care stau la baza modelelor de comportare dinamică. De fapt, modelul vâscoelastic neliniar con-ceput pentru determinarea răspunsului dinamic, prezentat în capitolul 7, include drept cazuri particulare modelele statice din capitolele anterioare 4, 5 şi 6. A fost însă preferată expunerea de la simplu la complex. Modelul vâscoplastic de la capitolul 8 are rolul de a pune în discuţie efectele incuderii explicite în model a vitezei de încărcare şi oferă posibilitatea evidenţierii diferenţelor dintre efectele solicitărilor lente în raport cu cele rapide. Capitolul 9 prezintă o metodă de liniarizare echivalentă destinată unei aproximări acceptabile a soluţiei neliniare. Concepută împreună cu colegul meu de la Institutul de Mecanica Solidelor dr.Tudor Sireteanu, metoda s-a dovedit utilă în aplicaţii şi credem că potenţialul său aplicativ nu este epuizat. Problema epuizării capacităţii de rezistenţă mecanică la pământuri şi roci face obiectul capitolului 10, unde sunt prezentate particularităţile specifice aces-tor materiale atât în cazul solicitărilor statice, dar mai ales în cazul celor di-namice. Penultimul capitol se abate într-o oarecare măsură de la titlul general enun-ţat – nu pământurile fac obiectul modelului vâscoelastic neliniar prezentat ci mate-rialele folosite la izolări antivibratorii. Insă, comportarea dinamică a acestor ma-teriale a permis adaptarea modelelor dinamice ale pământurilor Ultimul capitol prezintă o metodă de evaluare a efectelor provocate de neliniaritatea caracteristicilor dinamice ale materialelor asupra răspunsului di-namic al amplasamentului. Autorul ţine să mulţumească Consilului Naţional al Cercetării Ştiinţifice din Învăţământul Superior - CNCSIS - care prin atribuirea grantului 33517/2002 – tema 1 a permis îndeplinirea condiţiei obligatorii de editare a prezentei lucrări. Autorul mulţumeşte, de asemenea, tuturor celor care l-au susţinut în perioada de realizare şi de redactare finală. AUTORUL 11 NOŢIUNI GENERALE 1.1. Caracteristici mecanice definitorii În prezent nu există ecuaţii constitutive care să descrie toate însuşirile reale ale comportării mecanice a corpurilor deformabile. Fiecare ecuaţie constitutivă reprezintă o idealizare a însuşirilor reale şi furnizează o descriere aproximativă a acestor însuşiri, acceptabilă numai între anumite limite. Un grad mai mare de generalitate a ecuaţiei constitutive permite o modelare mai apropiată de realitate dar poate crea dificultăţi, uneori insurmontabile, atât în determinarea constantelor de material, cât şi la rezolvarea problemelor la limită. Din acest motiv, ecuaţia constitutivă nu poate fi decât un compromis, restrângerea generalităţii fiind costul plătit aplicabilităţii. Includerea sau excluderea din modelul matematic al unora sau altora dintre proprietăţile reale ale materialelor poate avea un răsunet important sau, din contră, un efect neglijabil asupra răspunsului structurii în ansamblu. Această constatare cu valabibilitate generală capătă o pondere însemnată în cazul materialelor cu comportare atât de complexă cum sunt pământurile. Denumirea generică de „pământuri” acoperă o mare gamă de materiale apărute în decursul timpului ca rezultat al proceselor de degradare a rocilor. Clasificările din diverse tratate sau standarde cuprind zeci de grupe şi subgrupe. Insă, din punctul de vedere al comportării mecanice două mari grupe prezintă o importanţă deosebită – nisipurile şi argilele – materiale care deşi au multe proprietăţi mecanice comune au şi diferenţe de comportare care impun utilizarea unor modele diferite. În mod obişnuit nisipurile au proprietăţi reologice reduse şi pot fi acceptabil modelate (în anumite condiţii care vor fi precizate pe parcurs) cu un model elastic neliniar iar argilele care, în mod frecvent, prezintă modificări semnificative în timp pot fi modelate cu un model vâscoelastic neliniar. Insă, deoarece există o mare varietate de materiale cuprinse în şi între grupele „nisipuri” şi „argile” o delimitare netă a tipului de comportare nu poate fi făcută. In plus, un acelaşi material poate avea comportări diferite în funcţie de parametrii ce definesc starea fizică a materialului, cum sunt umiditatea, gradul de îndesare, porozitatea, etc. Din aceste motive, considerăm că numai analiza de laborator a comportării probelor prelevate din amplasament poate indica modelul matematic cel mai adecvat cazului concret analizat. Modelarea caracteristicilor mecanice ale pământurilor are o puternică influenţă asupra evaluării răspunsului structural, structura fiind în acest context atât pachetul de straturi din amplasament, cât şi construcţia propiu-zisă ce reazemă şi conlucrează cu materialele din amplasament. O enumerare succintă a principalelor caracteristici mecanice definitorii pentru construirea unei ecuaţii constitutive, într-o ordine pe care o estimăm descrescătoare ca pondere, poate fi următoarea : neliniaritatea (definită de dependenţa caracteristicilor mecanice de nivelul de deformaţie provocat de solicitarea exterioară); problemele date de proprietăţile reologice ale acestor materiale, care se manifestă prin fluaj şi relaxare în cazul solicitărilor exterioare statice şi prin disipare în timpul solicitărilor dinamice; ireversibilitatea deformaţiilor; frecvenţa sau viteza solicitărilor exterioare. Desigur, există numeroşi alţi factori cu răsunet în modelarea răspunsului static sau dinamic, cum sunt: natura şi proprietăţile fizice ale materialului (compoziţia mineralogică, granulometria, coeziunea, indicele porilor, densitate, umiditate etc.), starea de tensiuni “in situ”, perturbările provocate de operaţiile de prelevare şi prelucrare a probelor, tipul şi timpul de consolidare a probelor etc. Insă, din punctul de vedere al modelării aceşti factori nu definesc tipul relaţiei constitutive, ci intervin în valorea parametrilor care cuantifică ecuaţia constitutivă. 1.1.1. Neliniaritatea Caracteristica mecanică principală a răspunsului pământurilor şi rocilor la solicitări exterioare statice sau dinamice este, indiscutabil, dependenţa de nivelul de deformaţie provocat de aceste solicitări. La materialele care admit ipoteza de corp elastic liniar (de exemplu, metalele) proporţionalitatea dintre tensiuni şi deformaţii se conservă până în apropierea ruperii. Relaţiile constitutive sunt, în acst caz, relaţii liniare de forma : (1.1) unde modulii elastici: K – modulul volumic şi G – modulul de forfecare au valori constante pe parcursul procesului de deformare până în apropierea epuizării capacităţii de rezistenţă mecanică. Deşi se consideră că la nivele reduse de solicitare şi pământurile au această comportare, practic acest palier liniar este atât de redus încât este pe deplin neglijabil. Experienţa a demonstrat că valorile modulilor determinate pe probe de pământ diferă în raport cu amplitudinea solicitării aplicate probei şi că aceste diferenţe apar încă de la nivele extrem de mici. Prima modalitate de a ţine cont de această dependenţă de nivelul de solicitare, sau de nivelul de deformaţie impus, a fost realizată prin înlocuirea constantelor elastice prin funcţii de tensiune sau de deformaţie în cadrul modelelor numite cu "moduli variabili"[65], [106], [108] [119]. Un exemplu în acest sens poate fi scris: (1.2) Aceste modele sunt cunoscute şi sub denumirea de modele elastic neliniare, ele încadrându-se în domeniul elasticităţii neliniare. Însă, încadrarea în elasticitatea, fie ea şi neliniară, presupune admiterea existenţei unui potenţial elastic W sau W care să permită ca tensorii tensiunii şi deformaţiei să poată fi legaţi prin relaţii de forma [130], [132], [133]: sau (1.3) Această ipoteză conduce, însă, la acceptarea reversibilităţii deformaţiilor. Există materiale cu relaţii neliniare între tensiuni şi deformaţii şi care prezintă reversibilitate dar pământurile nu pot fi incluse în această categorie deoarece acumularea deformaţiilor ireversibile este un fenomen curent şi specific acestor materiale. Însă, în anumite condiţii, elasticitatea neliniară poate fi folosită şi la pământuri. Astfel, pentru o primă solicitare statică monoton crescătoare, modelele cu moduli variabili au o susţinere teoretică şi dau rezultate satisfăcătoare în ceea ce priveşte concordanţa model-experiment deoarece nu se pune problema descărcării în urma căreia apar evidente deformaţiile ireversibile. Cazul solicitărilor dinamice este, din acest punct de vedere, mult mai complicat. Solicitările dinamice sunt compuse, în principiu, dintr-o succesiune de încărcări - descărcări - reîncărcări, proces în care acumularea deformaţiilor este evidentă şi provoacă efecte evidente. În figurile 1.1 şi 1.2 este prezentat schematic răspunsul în tensiuni sau deformaţii al pământurilor solicitate ciclic în teste triaxiale cu tensiune controlată (fig.1.1) sau cu deformaţie controlată. (fig.1.2). Este evidentă în ambele cazuri prezenţa răspunsului histeretic, diminuarea rigidităţii reflectată de diminuarea valorilor modulului de forfecare G şi acumularea deformaţiilor ireversibile. Aceste caracteristici fac inacceptabilă utilizarea elasticităţii neliniare, dar cu acceptarea unor ipoteze, poate fi utilizată vâscoelasticitatea neliniară. În primul rând trebue admisă înlocuirea disipării histeretice cu o disipare de tip vâscos : ipoteza de echivalare reo - histeretică [29], [76]. Apoi, se acceptă că dependenţa de deformaţie, neliniaritatea, este modelată de o „curbă – schelet” definită ca locul geometric al vârfurilor buclelor de histereză. Astfel, modelul vâscoelastic neliniar este construit pe baza a două funcţii : o funcţie modul dinamic dată de funcţia modul a curbei - schelet şi o funcţie de amortizare cu valori legate de aria buclelor de histereză (fig.1.3). Fig. 1.1 Fig. 1.2 Fig. 1.3 1.1.2. Efectele reologice Pentru un material vâscoelasic, modificările în timp ale caracteristicilor sale mecanice reprezintă o trăsătură definitorie. Fenomenul de evoluţie în timp a deformaţiilor sub tensiuni constante este numit fluaj iar modificările tensiunilor sub deformaţii constante este cunoscut sub denumirea de relaxare. Ambele fenomene se întâlnesc şi la pământuri, concomitent cu dependenţa caracteristicilor mecanice de nivelul de deformaţie, ceea ce impune luarea în considerare a vâscoelasticităţii neliniare. Fenomene de fluaj şi relaxare sunt întâlnite şi la pământuri, în special la argile, ceea ce face ca aceste materiale să poată fi considerate vâscoelastice. Vom folosi şi în această lucrare aceeaşi denumire de material „vâscoelastic” cu rezervele datorate termenului de „elastic” menţionate anterior. Însă, neliniaritatea atât de prezentă în toate manifestările acestor materiale face ca pentru pământuri să nu fie admisibile decât modelele vâscoelastic neliniare. Modelele vâscoelastice s-au dovedit utile şi în cazul solicitărilor dinamice, şi anume la modelarea disipării de energie, fenomen prezent la pământuri. Desigur, amortizarea în pământuri nu are caracter vâscos ci este, mai curând, de natură histeretică. Însă, fenomenologic, amortizarea din aceste materiale poate fi cu succes modelată cu modele vâscoelastice. Aşa se face că utilizăm modele dinamice vâscoelastic neliniare inclusiv la pământurile la care nu prezintă fenomenele reologice „clasice” de fluaj sau relaxare (de exemplu, nisipurile). 1.1.3. Ireveresibilitatea Primele modelări ale ireversibilităţii au apărut în cadrul modelelor elasto-plastice [60], [61] prin adaptarea unor modele concepute pentru alte materiale, în special metalice. În esenţă, aceste modele presupun existenţa unor legi constitutive de tip incremental, cu una din formele : (1.4) unde şi dε sunt variaţiile incrementale ale tensorilor tensiunii σ şi deformaţiei ε, iar H şi M sunt tensorii de ordinul patru ale căror componente depind, în general, de starea de tensiune şi/sau de deformaţie. Admiţând izotropia materialului, numai două din componentele tensorilor H şi M rămân distincte. În cazul tensorului H componentele distincte sunt echivalente cu modulii elastici ( ; ; ) iar în cazul tensorului M componentele distincte sunt inversele acestor moduli. În construirea modelelor elastoplastice sunt folosite, în general, ipotezele [61]: • Deformaţiile sunt compuse aditiv din deformaţii incrementale, elastice şi plastice (ireversibile): (1.5) (ipoteză discutabilă la pământuri unde deformaţiile reversibile sunt extrem de reduse) • Variaţia deformaţiilor elastice poate fi calculată în funcţie de tensiuni: (1.6) caz în care componentele tensorului M sunt constante. • Deformaţiile plastice sunt ireversibile iar variaţia incrementală a acestora depinde de starea de tensiune σ şi de deformaţia plastică acumulată : (1.7) • Deformaţiile plastice apar când tensiunile ating o anumită limită, numită limită de curgere. În spaţiul componentelor tensorului tensiunii σ, tensiunile limită formează o suprafaţă, numită suprafaţă de curgere, definită de o funcţie F : (1.8) unde h este un parametru de ecruisare care guvernează modificările acestei suprafeţe pe parcursul proceselor de încărcare, introducându-se astfel, parţial, efectul istoriei încărcărilor. • Deformaţiile plastice derivă dintr-un potenţial plastic : (1.9) unde λ este un parametru. Această lege de curgere este numită asociată dacă sau neasociată dacă cele două suprafeţe au forme diferite. După modul în care evoluează suprafaţa de curgere, modelele elstoplastice pot fi izotrope (caz în care F se extinde sau se contractă izotrop de-a lungul axei tensiunilor sferice) sau anizotrope (când F se şi translatează). Dintre modelele elastoplastice cu ecruisare izotropă menţionăm modelul Roscoe [116], [117] seria de modele tip „cap” [119] şi modelul Lade [90], iar dintre modelele anizotrope pe cele ale lui Mroz [106] şi Prevost [115]. Aceste modele se deosebesc prin forma analitică pe care o atribue funcţiei de potenţial şi legilor de curgere. Modelele vâscoplastice care se obţin prin extinderea celor elastoplastice rezultă prin modificarea ec.(1.9) în: (1.10) unde viteza de deformare vâscoplastică ocupă locul deformaţiilor plastice incrementale d , iar parametrul λ este înlocuit cu o funcţie de forma: dacă (1.11) Modelele vâscoplastice astfel obţinute preiau, însă, şi deficienţele modelelor elastoplastice din care provin, printre care menţionăm: • Concepute iniţial pentru modelarea comportării uniaxiale a metalelor, modelele elasoplastice presupun existenţa unui anumit prag până la care deformaţiile sunt reversibile şi după care deformaţiile devin ireversibile. Modelele elastoplastice create pentru pământuri menţin în general acest cadru ipotetic, care în cazul acestor materiale apare ca neadecvat, deformaţiile remanente apărând încă de la primele solicitări. • Extinderea la solicitări triaxiale prin introducerea suprafeţelor de curgere rămâne la nivel de speculaţie teoretică greu verificabilă experimental. • Forma incrementală a ecuaţiilor constitutive elastoplastice este utilizabilă în programele de rezolvare „step by step” după o operaţie de integrare, ceea ce reprezintă o complicaţie faţă de forma în raport cu componentele totale ale tensorilor σ şi ε. Aceste motive fac să apară mai indicată calea aleasă în această lucrare, şi anume renunţarea la postulatul existenţei suprafeţelor de curgere, definirea şi utilizarea unor ecuaţii constitutive între invarianţii tensorilor tensiunii şi deformaţiei, ecuaţii definite şi cuantificate prin utilizarea unor teste-diagnostic atât pentru determinarea formei funcţiilor constituente cât şi a constantelor de material (cap.7). 1.2. Geometria spaţiului invarianţilor 1.2.1. Invarianţi Deoarece în relaţiile constitutive vor interveni invarianţii tensorilor tensiunii şi deformaţiei, în aceste paragraf vor fi expuse câteva definiţii, proprietăţi şi semnificaţii geometrice ale invarianţilor tensorului Cauchy al tensiunii T, noţiuni care pot fi extinse la tensorul deformaţiei E printr-o simplă schimbare de notaţii [130], [132], [133], [134]. Se numesc vectori proprii, vectorii nenuli , care au proprietatea: (1.12) unde coeficienţii sunt valorile proprii ale tensorului T. Se observă că vectorii proprii nu sunt rotiţi de T, deci sunt invarianţi în direcţie. Ei pot fi obţinuţi aducând ec.(1.12) în forma: (1.13) Sistemul de ecuaţii liniare (1.13) are soluţii nebanale numai dacă: (1.14) Dezvoltând ec.(1.14) şi ordonând după puterile lui se obţine ecuaţia caracteristică: (1.15) unde coeficienţii sunt: (1.16) iar între coeficienţii şi rădăcinile ec.(1.15) există relaţiile: (1.17) Tensorul T fiind simetric, cele trei valori proprii sunt reale, iar dacă sunt şi distincte din ec.(1.13) vor rezulta vectorii proprii care sunt în acest caz ortogonali între ei. Ca orice tensor simetric de ordinul doi, tensorul tensiunii T este complet definit de trei valori principale şi de trei vectori proprii . Dacă alegem un nou reper format din vectorii proprii , atunci T ia forma diagonală: (1.18) unde v este un tensor ortogonal ale cărui componente sunt proiecţiile vectorilor pe vectorii bazei iniţiale , iar componentele matricii diagonale σ sunt valorile proprii , numite curent componente principale ale tensorului T. Prin definiţie, o funcţie scalară izotropă a unui tensor de ordinul doi I este invariantă dacă : , unde Q este orice tensor ortogonal care caracterizează o rotaţie în urma căreia tensorul T devine . Se observă că: (1.19) deci coeficienţii (1.17) ai ecuaţiei caracteristice (1.15) sunt invarianţi, numiţi invarianţii principali ai tensorului tensiunii T. Notând deviatorul tensorului tensiunii cu S, componentele sale carteziene cu iar cele principale cu : ; (1.20) invarianţii principali ai deviatorului vor fi: (1.21) Tot invarianţi sunt şi funcţiile scalare izotrope - tensiunea medie: (1.22) (utilizată şi în definirea deviatorului) şi intensitatea: (1.23) 1.2.2. Geometria spaţiului invarianţilor Spaţiul componentelor principale ale tensorului tensiunii T poate fi conceput ca un spaţiu vectorial izomorf cu spaţiul euclidian tridimensional. Dacă se alege un reper cartezian drept a cărui bază o formează vectorii proprii v(k), k = 1, 2, 3, coordonatele unui punct din acest spaţiu vor fi tensiunile principale (fig. 1.4). Fig. 1.4 O stare de tensiune va fi reprezentată în acest spaţiu printr-un punct de coordonate , vectorul de poziţie al acestui punct reprezentând tensorul T în forma sa diagonală. Modulul acestui vector este deci, intensitatea tensorului tensiunii. Proiectând vectorul t pe bisectoarea principală (diagonala spaţiului:) obţinem o interpretare geometrică a părţilor – sferică şi deviatorică – a tensorului tensiunii. Astfel, proiecţia vectorului t pe diagonala spaţiului va reprezenta partea sferică a lui T, iar vectorul diferenţă s deviatorul lui T. De aici şi denumirea de plan deviatoric atribuită planului perpendicular pe bisectoarea principală . Pentru a calcula modulele vectorilor sferic şi deviator proiectăm pe diagonala spaţiului şi, apoi, pe planul deviatoric componentele principale . Deoarece bisectoarea principală este egal înclinată faţă de cele trei axe cu unghiul , proiecţia vectorului t pe această diagonală va fi: (1.24) Modulul vectorului deviator poate fi obţinut proiectând tensiunile principale pe un plan deviator . Se observă din fig.1.4 că toate proiecţiile au lungimile afectate cu iar direcţiile se conservă. Ataşând planului deviatoric un sistem local de coordonate yoz (fig.1.5) punctul de stare se va proiecta pe acest plan în punctul P de coordonate : Fig. 1.5 (1.25) Rezultă astfel: (1.26) deci modulul vectorului deviator este intensitatea deviatorului tensorului tensiunii. De asemenea, modulul vectorului deviator este legat şi de invariantul numit tensiuni octaedrale. Planul care include vectorii t şi s (deci şi bisectoare principală) îl vom numi planul intensităţilor deoarece punctele de stare din acest plan pot fi raportate la sistemul de referinţă format de intensităţile σ şi τ. Înclinarea vectorului t faţă de bisectoarea principală este dată de invariantul direcţional α: (1.27) iar înclinarea planului intensităţilor faţă de reperul cartezian este dată de invariantul direcţional β: (1.28) unde invariantul λ este: (1.29) Ec. (1.28) se obţine din ecuaţia caracteristică a deviatorului S : în care se substituie , unghiul β fiind unghiul pe care îl face proiecţia deviatorului (axa τ) cu cea mai apropiată proiecţie a unei axe principale pe planul deviator. Fig. 1.6 În cazul particular în care tensorul tensiunii are numai două valori proprii distincte, iar spaţiul componentelor principale se restrânge la un plan numit planul simetriei axiale (fig.1.6). În acest plan, invariantul direcţional β este nul. Planul intensităţilor include în acest caz şi una din axele principale ( în fig.1.6), iar intensităţile σ şi τ au expresiile: (1.30) relaţii care se pot obţine şi prin particularizarea ec.(1.24) şi (1.26). 22 SISTEM OSCILANT CU UN GRAD DE LIBERTATE DINAMICĂ În cele ce urmează vor fi expuse câteva noţiuni din dinamica sistemelor elastice, noţiuni necesare definirii elementelor de dinamica pământurilor utilizare pe parcursul lucrării. Ne vom axa în special pe comportarea oscilatorului cu un grad de libertate dinamică supus unui regim staţionar torsional, iar relaţiile necesare excitaţiei axiale vor fi deduse prin adaptările necesare aduse formei torsionale. 2.1. Sistem liniar Ecuaţia diferenţială a mişcării unui sistem cu un grad de libertate dinamică în regim staţionar de vibraţii torsionale armonice este [50]: (2.1) unde J0 este momentul de inerţie masic al vibratorului care imprimă mişcarea armonică, iar c, k sunt constantele modelului limiar Kelvin-Voigt: k = rigiditatea resortului şi c = vâscozitatea amortizorului. Ecuaţia (2.1) poate fi adusă la forma : (2.2) unde : (2.3) Soluţia acestei ecuaţii este compusă din două părţi : , prima, , fiind soluţia ecuaţiei diferenţiale omogene ataşată ec.(2.1) : iar a doua, , o soluţie particulară a ecuaţiei diferenţiale neomogene (2.1). Prima parte a soluţiei, θ1, modelează vibraţia liberă iar prin adăugarea celei de a doua părţi θ2 se include efectul momentului perturbator ce transformă mişcarea într-o vibraţie întreţinută. Datorită amortizării sistemului, după un scurt regim tranzitoriu, θ1 tinde la zero iar mişcarea intră într-un regim staţionar cu θ = θ2. 2.1.1. Vibraţia liberă În regim de vibraţie liberă ecuaţia de mişcare a oscilatorului (2.2) se reduce la ecuaţia omogenă: (2.4) care admite o soluţie de forma : (2.5) Prin introducerea acestei forme în ecuaţia omogenă.(2.4) rezultă ecuaţia caracteristică: , cu soluţiile : (2.6) Deci, soluţa ecuaţiei omogene este: (2.7) sau ţinând cont că : e, rezultă : (2.8) În lipsa amortizării () soluţile (2.6) se reduc la : (2.9) deci ω0 din ec.(2.2) este pulsaţia proprie sau naturală a sistemului neamortizat. Atunci când amortizarea este prezentă soluţia diferă în funcţie de semnul discriminantului , adică depinde de valoarea relativă a constantei de amortizare c în raport cu celelalte constante ale sistemului (k, J0) : • Dacă : , amortizarea este numită critică, discriminantul ∆ este nul, rădăcinile (2.6) sunt reale şi egale: iar soluţia ec.(2.4) este de forma : (2.10) deci, o mişcare neoscilatorie, descrescătoare exponenţial în timp (deoarece ). • Dacă : , cele două rădăcini (2.6) sunt reale şi pozitive, caz în care soluţia (2.5) pierde caracterul oscilatoriu şi devine o mişcare descrescătoare exponenţial în timp iar sistemul este denumit sistem supraamortizat ; • Dacă : , adică : , sistemul este subamortizat, rădăcinile (2.6) sunt complex conjugate : (2.11) iar soluţia ec.(2.4), are forma : (2.12) unde , iar ωd este pulsaţia proprie a sistemului cu amortizare : (2.13) Soluţia (2.12) este o vibraţie pseudoarmonică cu pulsaţia ωd şi pseudoperioada : (2.14) graficul soluţiei fiind cuprins între exponenţialele (fig.2.1): (2.15) Gradul de amortizare a sistemului în regim de vibraţie liberă amortizată poate fi definit de raportul a două amplitudini succesive (sau situate la un număr n de perioade): (2.16) Prin definiţie, decrementul logaritmic este logaritmul natural al acestui raport : (2.17) iar raportul : (2.18) este numit raport de amortizare sau fracţiune din amortizarea critică şi este o măsură a amortizării utilizată curent la evaluarea amortizării pământurilor. Fig.2.1 Între pulsaţia proprie a sistemului cu amortizare ωd şi pulsaţia proprie neamortizată ω0 există relaţia (obţinută din ec.(2.9), (2.13) şi relaţia de definiţie a amortizării critice ): (2.19) Se observă că, deoarece , pulsaţia proprie amortizată ωd este întotdeauna inferioară pulsaţiei proprii neamortizate ω0 şi că ωd tinde la zero pe măsură ce coeficientul de amortizare c tinde către valoarea sa critică cc (fig.2.2). De asemenea, se poate constata că pentru rapoarte de amortizare mari, de ordinul a 40%, întâlnite la pământuri din relaţia (2.19) rezultă că : , deci cele două pulsaţii proprii au valori destul de apropiate şi practic se poate folosi numai pulsaţia proprie neamortizată ω0. Fig.2.2 Când sistemul are o deplasare şi o viteză îniţială la timpul : (2.20) este posibilă evaluarea celor două constante, B1 şi B2, din ec.(2.8) şi determinarea soluţiei. Astfel, situându-ne în cazul amortizării subcritice, constantele au forma: (2.21) iar soluţia tranzitorie are expresia: (2.22) relaţie în care s-a folosit aproximaţia . 2.1.2. Regim staţionar armonic Admitem drept soluţie particulară a ecuaţiei diferenţiale (2.2) o soluţie de forma : (2.23) Substituind această soluţie în ec.(2.2) şi colectând coeficienţii funcţiilor trigonometrice se obţine un sistem liniar de două ecuaţii în cele două necunoscute B3 şi B4 : (2.24) din care rezultă : (2.25) Rearanjând soluţia în forma rezultă : (2.26) şi introducând pulsaţia proprie neamortizată : , coeficientul de amortizare critică , raportul de amortizare şi notând raportul dintre pulsaţia excitaţiei şi pulsaţia proprie neamortizată cu: , obţinem : (2.27) unde θst este este rotirea statică provocată de amplitudinea M0 : (2.28) F este „factorul” de amplificare al vibraţiei forţate, în fapt, o funcţie dependentă de frecvenţa relativă υ şi de factorul de amortizare D (fig.2.3): (2.29) iar defazajul capătă expresia: (2.30) Soluţia finală a ecuaţiei diferenţiale liniare (2.1) poate fi acum obţinută prin reunirea soluţiei ecuaţiei omogene (2.12) cu soluţia particulară (2.27) : . Fig. 2.3 2.2. Sistem neliniar Constantele elastice au în cazul corpurilor liniare expresii cunoscute. Astfel, constanta elastică de rigiditate k pentru o bară cilindrică în cazul mişcării torsionale, este: [Nm] (2.31) unde Ip este momentul de inerţie polar al secţiunii , h înălţimea barei, φ diametrul ei, iar G modulul torsional. Dacă bara este solicitată axial, expresia pentru k este: [N/m] (2.32) unde este aria secţuinii transversale iar E modulul longitudinal Young. De asemenea, coeficienţii de amortizare c pot fi exprimat în funcţie de raportul de amortizare D şi au expresiile: [Nms] – pentru mişcarea torsională [Ns/m] – pentru mişcarea longitudinală (2.33) relaţii în care inervine pulsaţia proprie ω0 , masa m sau o altă caracteristică a masei – momentul de inerţie masic J0. În cazul pământurilor comportarea neliniară a materialului poate fi modelată prin înlocuirea modululilor şi rapoartelor de amortizare prin funcţii-modul şi de amortizare, funcţii dependente de nivelul de deformaţie. Cum relaţiile au o formă asemănătoare, vom aborda în continuare cazul solicitărilor torsionale, caz în care funcţiile de material au ca argumente deformaţia torsională γ: şi sau, ţinând cont de relaţia de legătură dintre deformaţia γ şi rotaţia θ : , funcţiile dinamice de material pot fi definite şi în raport cu rotirea θ : şi . Prin urmare, putem admite că dependenţa de γ sau θ afectează atât caracteristica elastică k : (2.34) cât şi caracteristica de amortizare c : (2.35) În cazul sistemelor liniare pulsaţia proprie ω0 este definită în funcţie de caracteristica elastică k : Însă, în cazul unui sistem neliniar k este o funcţie şi nu o constantă. Din acest motiv, vom defini drept pulsaţie proprie a sistemului neliniar, pulsaţia corespunzătoare valorii iniţiale θ = 0 : (2.36) Cu aceste precizări privind comportarea neliniară a materialului, sistemul oscilant (2.1) ia forma : (2.37) Efectuând schimbarea de variabilă şi introducând o nouă funcţie [124]: (2.38) prin derivare în raport cu noua variabilă, se obţine : (2.39) astfel că ec.(2.37) ia forma adimensionalizată : (2.40) unde : (2.41) Soluţia staţionară a ecuaţiei adimensionalizate (2.40) se obţine pe aceiaşi cale. Admitem drept soluţie particulară a ecuaţiei diferenţiale (2.40) o soluţie de forma: (2.42) Substituind această soluţie în ec.(2.40) şi colectând coeficienţii funcţiilor trigonometrice se obţine un sistem liniar de două ecuaţii în cele două necunoscute şi : (2.43) din care rezultă : (2.44) Rearanjând soluţia în forma , unde: (2.45) se obţine: (2.46) unde F este factorul de amplificare in formă echivalentă relaţiei (2.29): (2.47) 2.3. Soluţie numerică pentru oscilatorul neliniar 2.3.1. Metoda Newmark O soluţie analitică pentru sistemul neliniar (2.37) sau (2.40) este dificil de obţinut în toate cazurile de interes practic, dar rezolvarea poate fi făcută numeric cu o precizie suficientă. Vom prezenta în continuare o posibilitate de a obţine o soluţie numerică cu metoda de integrare Newmark care foloseşte diferenţele finite [93], folosind în acest scop o formulare generică pentru ecuaţia de mişcare a unui sistem neliniar cu un grad de libertate dinamică: (2.48) Intervalul de timp se împarte în n subintervale egale . Admitem că deplasarea x, viteza şi acceleraţia sunt cunoscute la timput de debut al unui subinterval oarecare : . Metoda constă în a determina valorile acestor mărimi la sfârşitul intervalului pornind de la valorile cunoscute la timpul şi folosind aproximaţii ale evoluţiei acceleraţiilor pe parcursul subintervalului. Pornind de la primul subinterval unde valorile la capătul din stânga sunt date de condiţiile iniţiale şi parcurgând secvenţial toate subintervalele se obţin în final valorile căutate. Dacă pe un interval acceleraţia este constantă şi are o valoare estimată ca medie a valorilor acceleraţiilor de la capetele intervalului: , prin integrare se obţine viteza şi deplasarea la capătul final al intervalului: (2.49) Dacă se admite o variaţie liniară a acceleraţiei: , va rezulta: (2.50) Newmark a ales o variantă intermediară între (2.49) şi (2.50) prin introducerea unui parametru β: (2.51) Se observă că dacă se regăseşte prima variantă cu acceleraţia constantă pe interval iar dacă ec.(2.51) se reduce la ec.(2.50) corespunzătoare variaţiei liniare a acceleraţiei pe intervalul . Acum, din ec.(2.48) se obţine pentru : (2.52) unde: (2.53) iar şi sunt date de ec.(2.51). Dacă valoarea obţinută cu ec.(2.52) este diferită de valoarea estimată iteraţia se reia înlocuind valoarea estimată anterior cu valoarea obţinută cu ec.(2.52) până când diferenţa devine inferioară unei anumite erori acceptabile : (2.54) 2.3.2. Program FORTRAN Metoda este relativ simplă dar rezolvarea efectivă a cazurilor concrete necesită efectuarea computerizată a iteraţiilor. Vom prezenta în continuare un program FORTRAN bazat pe metoda Newmark, program adaptat pentru PC în sistem interactiv [39]. • Programul principal: c Rezolvarea ecuatiei neliniare a oscilatorului neliniar c (ecuatia adimensionalizata) c REAL*8 M0 COMMON/DATE/X(50000),DX(50000),D2X(50000),D2Z(50000) COMMON / TIMP / TT(50000),ACC(50000) c Date de intrare : c - M0 = amplitudinea excitatiei c - DT = pas de timp (tau=t*puls.pr.) c - NT = nr.pasi de timp (maximum 50.000) WRITE(*,101) 101 FORMAT(5X,'INTRODUCE:M0=AMPLITUDINEA XCITATIEI'/) READ(*,201) M0 201 FORMAT(F10.0) WRITE(*,202) M0 202 FORMAT(5X,E14.7) WRITE(*,102) 102 FORMAT(5X,'INTRODUCE FRECVENTA EXCITATIEI',I2/) READ(*,201)FRECV WRITE(*,202) FRECV WRITE(*,103) 103 FORMAT(5X,'INTRODUCE DT=PASUL DE TIMP',I2/) READ(*,201) DT WRITE(*,202) DT WRITE(*,104) 104 FORMAT(5X,'INTRODUCE NT=NR.MAX.DE PASI',I2/) READ(*,203) NT WRITE(*,204) NT 203 FORMAT(I10) 204 FORMAT(5X,I10) PI=3.14159265 OPEN(UNIT=10,FILE=' ',FORM='FORMATTED') c c Calculeaza in NT puncte: TT = timpii c X = deplasari c DX = viteze c D2X = acceleratii c si le introduce in fisierul x.dat c (x=denumire data de utilizator) c X(1)=0.0 DX(1)=0.0 D2X(1)=0.0 WRITE(10,'(4A14)') 'Tau','Depl','Vit','Acc' KOD = 1 DO 10 I=1,NT TT(I)=(I-1)*DT D2Z(I)=M0*SIN(FRECV*TT(I)) CALL NEWMARK(DT,D2Z(I),X(I),DX(I),D2X(I),KOD) ACC(I)=D2Z(I)+D2X(I) WRITE(10,'(4E14.7)') TT(I),X(I),DX(I),D2X(I) 10 CONTINUE CLOSE(UNIT=10) STOP END • Subrutina de integrare Newmark SUBROUTINE NEWMARK(DT,D2Z,X,DX,D2X,KOD) K=0 ERR=1.0E-6 EPS=1.E+14 IF(KOD.EQ.0) GO TO 2 KOD=0 X0=X DX0=0. D2X0=0. D2X1=0.1 GO TO 4 2 D2X1=D2X 4 CONTINUE DX=DX0+(D2X0+D2X1)*DT*0.5 X=X0+DX0*DT+(0.33*D2X0+0.17*D2X1)*DT*DT DD = C(X)*DX EE = RK(X)*X D2X=-D2Z-DD-EE IF(DX.GT.EPS) GO TO 10 CIF=ABS(D2X-D2X1) IF(ABS(D2X).GT.0.001) CIF=CIF/ABS(D2X) IF(CIF.LT.ERR) GO TO 13 K=K+1 IF(K.LE.19) GO TO 2 ERR=ERR*10. IF(CIF.LT.ERR) GO TO 13 K=0 GO TO 2 13 X0=X DX0=DX D2X0=D2X RETURN 10 WRITE(*,901) 901 FORMAT(' Neconvergent ') RETURN END • Subrutina de calcul a valorilor funcţiilor de material c ARGILA FUNCTION C(X) A=ABS(X) c C=0.415 C=0.308-0.269*EXP(-118.625*A) RETURN END FUNCTION RK(X) A=ABS(X) c RK=0.515 RK=0.25+0.75*EXP(-111.982*A) RETURN END 2.3.3. Aplicaţii Programul a fost aplicat în numeroase cazuri şi a dat rezultate corecte. În cele ce urmează vom prezenta trei aplicaţii. Prima constă în determinarea soluţiei unui oscilator liniar, a doua se va referi la determinarea soluţiei liniare a sistemului oscilant al coloanei rezonante pentru o treaptă de încărcare iar a treia aplicaţie va avea ca obiect determinarea soluţiei neliniare a aceluiaşi sistem oscilant al coloanei rezonante. În primul caz, cel al unui sistem liniar, verificarea este imediată deoarece există posibilitatea determinării analitice a soluţiei. Celelalte două aplicaţii au fost alese pentru verificare deoarece soluţia determinată cu programul de calcul poate fi confruntată cu datele experimentale. • Aplicaţia nr.1 Fie sistemul oscilant din fig.2.5, excitat de o forţă , cu următoarele caracteristici mecanice şi de încărcare: (2.55) Din aceste date pot fi obţinute alte caracteristici mecanice – pulsaţia proprie şi raportul de amortizare D: Fig. 2.5 (2.56) Ecuaţia de mişcare are în acst caz forma similară ec.(2.1): (2.57) ecuaţie cu aceiaşi forma adimensionalizată (2.40): (2.58) unde: (2.59) Soluţia analitică (staţionară) este (2.46): (2.60) unde, în acest caz la rezonantă, factorul de amplificare F din ec.(2.29) şi defazajul din ec.(2.30) au valoarile: (2.61) astfel că soluţia devine: (2.62) Revenind la variabila iniţială t: şi ţinând cont că: se obţine soluţia staţionară: (2.63) Programul de calcul trebuie să ruleze în acest caz cu datele iniţiale: k c F = F0 sin ωt m (2.64) Fig. 2.6 Fig. 2.7 Rezultatul calculului este introdus într-un fişier de date denumit de utilizator la solicitarea programului (xxxx.dat) şi poate fi transferat şi prelucrat într-o aplicaţie grafică. O astfel de prelucrate este dată în fig. 2.6. Se remarcă o bună concordanţă cu soluţia analitică staţionară din fig. 2.7. • Aplicaţia nr. 2 Fig. 2.8 Sistemul oscilant al coloanei rezonante Drnevich [142] format din ansamblul vibrator-probă este un sistem cu un grad de libertate dinamică. Pentru o treaptă de solicitare, de regulă în regim armonic staţionar, acest sistem poate fi considerat liniar. Considerăm un astfel de sistem (fig.2.8) supus unei excitaţii de forma . Vibratorul are un moment de inerţie masic iar proba (de argilă, în acest caz) are un modul de torsiune iniţial , iar pentru frecvenţa de rezonanţă , un modul torsional , un raport de amortizare , iar deformaţia impusă probei este . Din aceste date, preluate dintr-un experiment în coloana rezonantă, pot fi calculate următoarele caracteristici de material şi de încărcare: k c M = M0 sin ωt J0 (2.65) (2.65) Astfel, sistemul oscilant este defint şi se poate trece la calculul soluţiei în cele două variante – rezolvarea analitică şi rezolvarea cu ajutorul programului de calcul numeric. Soluţia analitică are în acest caz forma: (2.66) unde: (2.67) iar soluţia va fi: (2.68) deci, amplitudinea rezultă: (2.69) Apelând la programul de calcul, vom introduce datele iniţiale: (2.70) iar din prelucrarea datelor obţinute rezultă o amplitudine a rotaţiei staţionare de . Compararea celor două rezultate, date în figurile 2.9 şi 2.10, indică o funcţionare corectă a programului. În acest caz avem şi o verificare în plus, şi anume, valoarea amplitudinii deformaţiei se poate calcula direct din datele experimentale iar valoarea obţinută 0,035 rad este destul de apropiată de cele două valori obţinute anterior. Fig. 2.9 Fig. 2.10 • Aplicaţia nr.3 Vom folosi pentru această exemplificare tot sistemul oscilant al coloanei rezonante, sistem format din vibrator şi proba de pământ care, pe parcursul întregului proces de deformare, imprimă întergului sistem o comportare neliniară. Modelul analog va fi acelaşi ca la aplicaţia anterioară (fig.2.8), însă locul caracteristicilor liniare c şi k este preluat de funcţii de material determinate experimental, de forma şi : (2.71) În acest caz, pentru a verifica funcţionarea, progamul a fost rulat în varianta neliniară a subrutinei de tip FUNCTION cu funcţiile de material din relaţia (2.71). Pentru o frecvenţă relativă fixă ( , deci la rezonanţă) şi cu diverse valori ale excitaţiei adimensionalizate µ (aceleaşi valori ale excitaţiei folosite experimental), soliţiile neliniare staţionare sunt date în fig.2.11 iar comparaţia experiment – calcul este dată în fig.2.12. Se remarcă o apropiere satisfăcătoare. Fig. 2.11 Fig. 2.12 33 TESTE de DEFINIRE şi CUANTIFICARE 3.1. Introducere Toate tentativele de modelare sunt strict dependente de posibilităţile de determinare experimentală a funcţiilor sau parametrilor necesari cuantificării proprietăţilor mecanice ale materialului. Cu cât materialul este mai complex din punct de vedere al răspunsului mecanic, cu atât este mai dependent de experiment. Iar pămâturile şi rocile sunt un elocvent exemplu. Desigur, ar fi ideal dacă pentru fiecare material sau model în parte s-ar putea imagina şi utiliza experimentele cele mai indicate. Însă, realitatea ne obligă să utilizăm la maximum aparatura existentă şi numai în cazul în care aceasta se dovedeşte insuficientă şi inadecvată putem apela la conceperea unor noi aparaturi şi tehnologii de testare. Modelele care vor fi expuse în continuare, sunt astfel concepute încât să se bazeze pe o aparatură existentă, dar cu modificări ale tehnologiei de testare care să satisfacă necesităţile de cuantificare. Aparatele existente folosite sunt : triaxialul şi coloana rezonantă, iar posibilităţile oferite de aceste aparate larg utilizate în geotehnică şi mecanica pământurilor vor fi expuse, succint în continuare. Menţionăm că în determinarea proprietăţilor mecanice ale materialului aflat "in situ" prelevarea şi prelucrarea probelor prezintă o mare importanţă. Cu toate precauţiile care pot fi luate, prelevarea probei perturbă structura materialului şi produce anularea stării de tensiune pe care materialul a avut-o „in situ". Pentru a contracara, cel puţin parţial, aceste perturbări, înainte de a fi solicitată proba trebuie "consolidată" cel puţin 24 de ore, adică supusă unei stări constante de tensiune asemănătoare celei avute în configuraţia iniţială, stare care poate fi aproximată astfel : (3.1) unde este presiunea verticală produsă de sarcina geologică, iar coeficientul presiunii laterale în stare de repaus [97]. Decomprimarea şi recomprimarea probei modifică, într-o oarecare măsură, proprietăţile mecanice naturale ale materialului, acest fapt constituind un dezavantaj inevitabil al încercărilor de laborator. 3.2. Teste de compresiune triaxială În aparatul de compresiune triaxială o probă cilindrică (uzual, cu diametrul şi înălţimea ), este supusă unei presiuni hidrostatice , constante sau variabile, echivalentă tensiunilor sferice pe care materialul le-a suportat „in situ”, la care se adaugă o încărcare axială variabilă egală cu forţa axială verticală dată de piston, forţă repartizată uniform pe suprafaţa superioară a probei . Într-o experienţă triaxială uzuală forţa verticală este continuu crescută până la rupere impunând probei o anumită viteză de deformare şi se măsoară la anumite trepte de timp atât valoarea forţei axiale, cât şi unele din modificările geometrice ale probei. Dacă axa de simetrie a materialului coincide cu axa verticală a probei , starea de tensiune şi deformaţie a probei este axial-simetrică cu şi (fig. 3.1). Se observă că denumirea aparatului este impropie, dar fiind folosită în mod curent o vom menţine. Admiţând că starea de tensiune indusă în probă de solicitările exterioare este uniformă, tensiunile principale ale unei trepte de timp i pot fi obţinute direct din condiţiile de contur : (3.2) unde este valoarea presiunii din celulă la treapta i (care poate fi şi constantă la toate treptele), este valoarea forţei verticale transmisă de piston iar este secţiunea transversală medie a probei la treapta i: (3.3) relaţie în care s-a notat prin înălţimea iniţială a probei, prin volumul său iniţial iar prin şi variaţiile acestor mărimi pe parcursul treptelor de timp. De asemenea, admiţând că deplasările verticale ale fiecărei trepte au o variaţie liniară pe înălţimea probei iar cele radiale au variaţii liniare în sens radial (fig. 3.1), deformaţiile treptei rezultă în forma: (3.4) La aparatele obişnuite pot fi măsurate numai deplasările axiale . Deplasările radiale pot fi măsurate numai cu ajutorul unor dispozitive optice sau tensometrice care nu sunt în dotarea curentă. În plus, pe parcursul deformării probele nu rămân cilindrice datorită frecărilor care apar pe bazele probelor, astfel că diametrul probei variază pe înălţime. Din acest motiv este preferabilă determinarea unor deformaţii medii radiale cu ajutorul variaţiei de volum care poate fi măsurată. Astfel, ţinând cont că: , obţinem: (3.5) Numărul treptelor de timp poate fi astfel ales încât stările de tensiune şi de deformaţie ale unei trepte să poată fi considerate infinitezimale. Însă, o treaptă infinitezimală se dezvoltă în prezenţa tensiunilor şi deformaţiilor anterioare, care prin acumulare devin din ce în ce mai mari, asfel că principiul clasic al suprapunerii efectelor devine inaplicabil. Din acest motiv, este preferabil ca sumarea treptelor să urmeze un procedeu ameliorat [11], [41], procedeu bazat pe suprapunerea unor deformaţii infinitezimale peste o deformaţie finită [133]: (3.6) unde este tensorul deformaţiei neliniare de la sfârşitul treptei, este acelaşi tensor la începutul treptei, Ereprezintă tensorul deformaţiilor infinitezimale apărute pe parcursul treptei iar este gradientul deformaţiei datorate deplasării particulei materiale din configuraţia iniţială la configuraţia de la începutul treptei. Adaptând relaţia tensorilaă (3.5) la starea axial simetrică din triaxial obţinem: (3.7) După deteminarea componentelor principale ale tensorilor tensiunii şi deformaţiei putem obţine valorile de la fiecare treaptă ale invarianţilor tensorului tensiunii: (3.8) şi deformaţiei: (3.9) Fig. 2.1 În testele obişnuite, cu denumirea uzuală de “statice”, forţa axială este crescută continuu sau în trepte, cu o viteză constantă sau variabilă în timp, dar cu variaţie lentă. Viteza de deformare impusă astfel probei poate varia în limite relativ largi, ordinul de mărime fiind, de regulă, inclus în intervalul 10 . Datorită proprietăţilor neliniare şi reologice ale pământurilor, aspectul curbelor tensiune deformaţie este dependent de viteza în care este condus experimentul şi datorită acestui fapt, triaxialul poate fi folosit şi pentru obţinerea unor informaţii calitatative şi cantitative privind influenţa vitezei de deformaţie asupra răspunsului mecanic al acestor materiale. 3.3. Teste triaxiale de fluaj Un alt tip de încercări care pot fi efectuate în triaxial sunt încercările pentru determinarea caracteristicilor reologice ale pământurilor. După cum se ştie, la pământuri sunt prezente atât fenomenele de fluaj (creşterea deformaţiilor sub tensiune constantă), cât şi fenomenele de relaxare (scăderea tensiunilor sub deformaţie impusă constantă). Însă, datorită posibilităţilor tehnice pe care le oferă aparatul triaxial obişnuit, numai încercările de fluaj pot fi efectuate şi numai după înlocuirea sistemului de încărcare continuă cu un sisten de încărcare statică. În acest caz, încărcările sunt aplicate în trepte menţinute un anumit interval de timp, interval în care modificările geometrice ale probei sunt măsurate din timp în timp. Astfel, în fiecare interval de timp are loc o experienţă de fluaj. Modalitatea de utilizare a unor astfel de teste va fi prezentată în paragraful 6.2. 3.4. Modelarea testului din coloana rezonantă 3.4.1. Coloana rezonantă - descriere sumară Coloana rezonantă (resonant column) - tip Drnevich [142] - este un aparat destinat determinării experimentale a răspusului dinamic al pământurilor prin solicitarea unei probe cilindrice (column) cu vibraţii armonice staţionare longitudinale şi/sau torsionale în regim de rezonanţă (resonant). Proba utilizată în coloana rezonantă este asemănătoare în formă, dimensiuni şi mod de preparare cu proba triaxială. Probele sunt cilindrice, uzual cu două dimensiuni - una cu diametrul şi înălţimea cm şi alta cu şi cm. Baza inferioară a probei este fixată, iar la baza superioară se ataşează dispozitivul electromagnetic de vibrare. Proba împreună cu vibratorul, ansamblu ce constituie sistemul rezonant, sunt amplasate într-un tub cilindric transparent („celulă" asemănătoare celei triaxiale) în care se introduce o presiune hidrostatică cu ajutorul aerului, apei sau uleiului mineral. Mărimea acestei presiuni este asfel aleasă încât tensiunile sferice provocate probei să fie identice cu cele pe care materialul le-a suportat „in situ”. Dispozitivul de excitare torsională este format dintr-un sistem de bobine electromagnetice – magneţi permanenţi plasaţi în jurul bazei superioare iar vibratorul longitudinal constă dintr-o bobină montată coaxial cu proba. Semnalul de excitare este emis de un generator variabil de impulsuri sinusoidale, este apoi amplificat de un amplificator de putere şi introdus în bobinele vibratorului. Acest semnal de intrare este, totodată, conectat la axa orizontală a unui osciloscop şi măsurat cu un multimetru. Momentul de torsiune rezultat şi/sau forţa longitudinală sunt direct proporţionale cu intensitatea curentului excitant, raportul moment/intensitate şi/sau forţă/intensitate fiind constante ale aparatului. Mişcările, orizontală şi/sau verticală, sunt măsurate de accelerometre plasate în axa de rotaţie pentru mişcarea verticală şi la o distanţă de 3,175 cm de axa de rotaţie pentru oscilaţia rotaţională. Acceleraţiile pot fi convertite în deplasări prin împărţire cu ( în cazul mişcării verticale sau cu în cazul mişcării orizontale. Semnalul electric de ieşire dat de accelerometre este amplificat de un amplificator electronic (cu un factor de calibrare , unde n = 0,1 sau 2) şi apoi conectat la axa verticală a osciloscopului. Totodată, tensiunea acestui curent este măsurată cu multimetrul în (milivolţi „root-mean-square”) raportul fiind o constantă de calibrare. Pentru determinarea frecvenţei de rezonanţă este folosit ecranul osciloscopului X-Y, care are pe axa X semnalul de intrare şi pe axa Y pe cel de ieşire. Frecvenţa excitaţiei este mărită până când apare rezonanţa sistemului probă-generator, moment în care pe ecranul osciloscopului apare o elipsă Lissajous cu axele suprapuse pe X şi Y şi mărimea proporţională cu intensitatea semnalului de intrare. În acest moment, se înregistrează frecvenţa de rezonanţă [Hz], acceleraţia la rezonanţă A şi intensitatea curentului C. Apoi, se impune frecvenţa şi se înregistrează din nou acceleraţia A' şi intensitatea C'. Aceste valori sunt necesare determinării caracteristicilor dinamice ale materialului din probă : modulii elastici (torsional G şi/sau longitudinal E), rapoartele de amortizare (torsional Dt şi/sau longitudinal Dl) precum şi a amplitudinilor deformaţiilor (torsionale γ şi/sau longitudinale ε). Amplitudinea excitaţiei fiind suficient de redusă putem admite ipoteza [29], [47], [76] : • Pe parcursul vibraţiilor staţionare proba se comportă ca un corp vâscoelastic liniar echivalent cu aceeaşi masă, densitate, dimensiuni geometrice şi frecvenţă de rezonanţă. Această ipoteză, numită ipoteza reo-histeretică, stă la baza interpretării datelor experimentale. Ipoteza reo-histeretică, pe lângă echivalea fenomenologică dintre comportarea vâscoasă şi cea histeretică, implică faptul că încercarea din coloana rezonantă poate fi considerată o încercare nedistructivă, deci o aceeaşi probă poate fi folosită la evaluarea caracteristicilor dinamice corespunzătoare unor amplitudini tot reduse dar diferite. De asemenea, ipoteza permite ca, sub un regim staţionar al forţei perturbatoare, tensiunile şi deformaţiile din probă să poată fi considerate staţionare şi, în consecinţă, modulii şi rapoartele de amortizare determinate să rămână constante. Tot ca o consecinţă a ipotezei de mai sus, apare şi posibilitatea de a modela răspunsul materialului la un nivel constant al excitaţiei cu ajutorul unui model Kelvin-Voigt echivalent (fig. 3.2). Momentul torsional va provoca în probă unde torsionale iar forţa perturbatoare unde longitudinale. Deoarece în ambele cazuri ecuaţiile de mişcare au forme similare vom aborda cazul vibraţiilor torsionale. Fig. 3.2 3.4.2. Torsiunea probei Considerăm proba ca un cilindru circular vertical, încastrat la baza inferioară şi liber la cea superioară, solicitat de un moment de torsiune M al cărui vector este dirijat dupa axa verticală (fig. 3.3) Admiţând ipoteza secţiunilor plane a lui Bernoulli, ca efect al torsiunii, secţiunile transversale rămân plane, dar se rotesc (sau „lunecă”) una faţă de alta. Dacă deformaţiile sunt relativ mici, se poate scrie: (3.10) unde γ este lunecarea specifică, θ este rotirea secţiunii curente iar prin α s-a notat torsiunea specifică. Tensiunile tangenţiale au o variaţie liniară, sunt normale pe rază şi dirijate astfel pentru a echilibra momentul de torsiune M. Rezultă astfel : (3.11) Deoarece iar , unde este momentul polar de inerţie al secţiunii transversale, expresia momentului torsional devine : şi ţinând cont că , se obţine: (3.12) Produsul poartă numele de rigiditate la torsiune. Fig. 3.3 h = 8 cm = 3,57 cm R = 1,785 cm ra = 3,175 cm rm = 2R/3 Din (3.10) şi (3.12) se poate obţine prin integrare rotirea specifică θ : , deci : (3.13) Utilizând notaţiile din fig.2.3, avem : (3.14) rezultând astfel o relaţie de legătură între rotirea θ şi lunecarea specifică γ , relaţie care pentru proba uzuală (cu şi cm) este: (3.15) 3.4.3. Determinarea amplitudinii rotaţiei şi deformaţiei torsionale Deoarece proba este solicitată de un moment de torsiune cu variaţie armonică admitem că şi rotirea θ va fi de aceiaşi formă : , iar acceleraţia mişcării de rotaţie va fi de forma : . Notând cu traseul accelerometrului montat pe vibrator la o distanţă de ra = 3,175 cm de axa de rotaţie şi cu valoarea acceleraţiei furnizate de accelerometru, amplitudinea rotaţiei rezultă: (3.16) dacă valoarea acceleraţiei A este măsurată în rad/s2, sau : (3.17) dacă A este măsurat, ca deobicei, în mvrms . În relaţiile anterioare f este frecvenţa măsurată în Hz, ACF este factorul de amplificare al semnalului de ieşire, dat în pk-mv/pk-g, g fiind acceleraţia gravitaţională (g = 9,806 m/s2). Cu o amplificare obişnuită, ACF = 2500 pk-mv/pk-g, expresia amplitudinii rotaţiei se reduce la : [rad] (3.18) relaţie în care s-a renunţat la indicele "0" al amplitudinii θ0. Deformaţia torsională γ se poate obţine cu ajutorul ec.(3.15), rezultând : [ - ] sau [ %] (3.19) 3.4.4. Determinarea modulului torsional Considerăm ecuaţia propagării undelor torsionale prin proba considerată mediu vâscoelastic liniar [132]: (3.20) unde θ este unghiul de torsiune, este viteza de propagare: (3.21) G este modulul de torsiune iar ρ este densitatea materialului. Introducând în ec.(3.20) o soluţie de forma , rezultă: (3.22) unde raportul dintre pulsaţia naturală şi viteza de propagare este numit pulsaţie adimensională sau număr de undă. Tinând cont de condiţiile de capăt, la baza fixă x = 0 torsiunea fiind împiedicată iar la cea superioară x = h momentul de torsiune al probei trebuie să egaleze momentul de inerţie al vibratorului şi momentul perturbator: (3.23) din ec.(3.22) rezultă: (3.24) unde este momentul polar al secţiunii probei iar este raportul dintre momentul de inerţie masic al probei J: (3.25) şi cel al vibratorului : J0 = 2,8 10-3 N.m.s2 . De asemenea, în ec.(3.25) s-a notat cu m masa probei iar cu w greutatea sa. Frecvenţa excitaţiei poate fi variată până la apariţia rezonanţei, când din ec.(2.24) rezultă ecuaţia pulsaţiilor proprii : (3.26) Ecuaţia transcendentă (3.26) poate fi rezolvată prin dezvoltare în serie rezultând rădăcina în forma [22], [29]: (3.27) şi ţinând cont că iar, la rezonanţă, viteza de propagare se obţine în forma: (3.28) Din ec.(3.28) şi ec.(3.21) se obţine, apoi, modulul torsional : . În momentul în care greutatea probei este cunoscută, expresia modulului torsional G poate fi adusă la o formă mai simplă. Vom ilustra această modificare printr-un exemplu, în care admitem că urmează să introducem în aparat o probă de argilă cu o greutate de 158 g. În acest caz, momentul de inerţie masic poate fi obţinut : (3.29) deci şi raportul: (3.30) Radicalul din (3.28) are acum valoarea , iar viteza de propagare devine: (3.31) Densitatea probei ρ poate fi de asemenea calculată: [kg/m3 = Ns2/m4] (3.32) deci, modulul torsional G rezultă în forma: (3.33) relaţii modificate în mică măsură de variaţii ale greutăţii probelor. 3.4.5. Estimarea amplitudinii momentului perturbator Momentul care provoacă o stare de tensiune şi deformaţie torsională în probă este dat de: (3.34) relaţie rezultată din ec.(3.13) şi (3.14). Considerând proba de exemplificare rezultă: [m4] (3.35) [Nm] (3.36) 3.4.6. Determinarea raportului de amortizare torsională Calculul capacităţii de amortizare la un anumit nivel al amplitudinii excitaţiei se bazează pe ipoteza de echivalenţă reo-histeretică, ipoteză care permite utilizarea unui model echivalent Kelvin-Voigt în evaluarea energiei disipate în sistemul probă-vibrator. Ecuaţia diferenţială a mişcării modelului Kelvin-Voigt în regim staţionar de vibraţii torsionale este [50]: (3.37) unde : (3.38) β fiind factorul de amortizare, pulsaţia proprie a probei iar c, k constantele modelului Kelvin-Voigt (k = rigiditatea resortului şi c = vâscozitatea amortizorului). Cu ajutorul coloanei rezonante amortizarea se determină prin două metode: prin metoda vibraţiilor libere amortizate şi prin metoda vibraţiilor staţionare. În continuare aceste metode vor fi expuse succint, iar pentru detalii se poate apela la paragraful 2.1. Prima metodă constă în întreruperea bruscă a momentului perturbator ceea ce conduce la apariţia unor oscilaţii libere amortizate a sistemului probă-vibrator, oscilaţii descrise de ecuaţia omogenă ataşată ec.(3.37): (3.39) Soluţia acestei ecuaţii, pentru cazul amortizării subcritice al coloanei rezonante, este de forma : (3.40) deci, o vibraţie pseudoarmonică cu pseudopulsaţia şi pseudoperioada (fig.2.1). O măsură a intensităţii amortizării este decrementul logaritmic : (3.41) unde este raportul de amortizare (dintre amortizarea curentă c şi amortizarea critică ) sau fracţiunea din amortizarea critică, măsura amortizării utilizată în mod curent. Având înregistrată curba vibraţiei libere amortizate, se pot măsura două amplitudini succesive sau situate la n pseudoperioade şi ţinând cont că la amortizări mici , se obţine : (3.42) Cea de a doua metodă, numită metoda vibraţiilor staţionare sau metoda factorului de amplificare, utilizează factorul de amplificare (2.27): (3.43) care, la rezonanţă, când şi , devine , de unde rezultă: (3.44) Amplitudinea momentului perturbator M0 este proporţională cu intensitatea C a curentului ce trece prin bobinele vibratorului, raportul de proporţionalitate fiind o constantă de etalonare. De asemenea, amplitudinea mişcării la rezonanţă poate fi obţinută din acceleraţia mişcării la rezonanţă A afectată de o altă constantă de etalonare. Rezultă astfel : (3.45) În cazul coloanei rezonante Drnevich, factorul de calibrare kd poate fi exprimat în funcţie de valorile acceleraţiei şi a intensităţii curentului, valori culese la [142] : (3.46) Pentru proba uzuală valoarea expresiei din paranteza mare a ec.(3.46) este apropiată de 1, astfel că poate fi neglijată. Toate determinările expuse în acest paragraf presupun existenţa unei excitaţii cu o anumită amplitudine constantă şi a unei anumite presiuni constante în celula aparatului şi conduc la obţinerea unor valori ale modulilor dinamici şi factorilor de amortizare, valori corespunzătoare nivelului de deformaţie provocat probei de această amplitudine a excitaţiei. În cadrul ipotetic admis, încercarea din coloana rezonantă poate fi considerată nedistructivă, deci după terminarea unei încercări pe aceeaşi probă se poate efectua o nouă încercare cu o altă amplitudine a excitaţiei. Modificarea amplitudinii excitaţiei sau a presiunii din celulă conduce la obţinerea altor valori ale caracteristicilor dinamice, fapt care reflectă comportarea neliniară a materialului din probă. Dacă o singură încercare cu o anumită excitaţie constantă poate fi modelată de un model Kelvin-Voigt, diferenţa dintre rezultatele obţinute cu solicitări diferite nu mai poate fi explicată de acest model liniar. Apare deci, evidentă, necesitatea introducerii unui model neliniar. Cu titlu de exemplificare redăm în tabelul 3.1 datele experimentale obţinute în coloana rezonantă pe o probă de argilă, precum şi rezultatele prelucrării lor în vederea obţinerii valorilor γ, G şi D pentru 12 nivele de solicitare, valori calculate cu relaţiile: (3.47) Tabelul 2.1 Nivel γ G D [Hz] [mvrms] [mvrms] [mvrms] [mvrms] [%] [MPa] [%] 1 40,193 98,6 500 68 102 0,020 89,820 3,287 2 37,778 255 908 193 260 0,042 79,354 5,212 3 36,342 376 1130 303 390 0,056 73,433 6,463 4 31,746 812 1670 733 840 0,109 56,034 10,607 5 29,222 1160 2035 1010 1175 0,157 47,480 12,249 6 27,586 1360 2310 1170 1385 0,200 42,312 12,434 7 26,350 1560 2445 1345 1600 0,232 38,606 13,409 8 25,880 1700 2535 1370 1730 0,249 37,239 13,277 9 24,480 2030 2810 1600 2070 0,309 33,319 13,960 10 23,579 2380 3100 1800 2400 0,367 30,913 14,395 11 21,901 2775 3400 2100 2780 0,467 26,669 15,413 12 21,066 2985 3555 2180 2970 0,528 24,674 15,408 44 RELAŢII CONSTITUTIVE ELASTIC NELINIARE Modelele elastic neliniare ale pământurilor au fost primele modele care au luat în considerare comportarea neliniară a acestor materiale. Aceste modele au făcut obiectul a numeroase cercetări teoretice şi experimentale, dintre care menţionăm lucrările [68], [88], [106], [108]. Cu toată diversitatea punctelor de vedere conţinute, aceste lucrări au o trăsătură comună dată de utilizarea unor „moduli variabili”, autorii respectivi înlocuind constantele elastice , sau cu funcţii empirice care au ca argumente componente ale tensorilor tensiunii sau deformaţiei. În continuare va f expus un alt mod de abordare, bazat pe rolul de potenţial elastic al energiei de deformaţie ceea ce permite o analiză clară a dependenţei relaţiilor obţinute de condiţiile de obţinere şi, totodată, o analiză a posibilităţilor de extindere la alte tipuri de solicitare [10], [18]. 4.1. Forma generală a ecuaţiei constitutive Admitem că energia de deformaţie specifică a pământurilor care pot fi considerate corpuri elastice are rol de potenţial elastic, ceea ce permite utilizarea ecuaţiei constitutive neliniare în forma [132], [133]: (4.1) unde este cel de al doilea tensor Piola-Kirchhoff al tensiunii (tensor care raportează tensiunile din configuraţia curentă la unitatea de arie din configuraţia de referinţă existentă „in situ”), E este tensorul deformaţiei neliniare iar w este densitatea energiei de deformaţie specifică pe unitatea de volum pe care corpul a avut-o în configuraţia naturală K'. Prin şi au fost notate componentele tensorilor respectivi raportate la sistemul de coordonate ataşat configuraţiei naturale. Interesează, însă, forma acestei ecuaţii în cazul în care se alege drept configuraţie de referinţă configuraţia iniţială existentă „in situ” la momentul iniţial , configuraţie căreia îi ataşăm sistemul de coordonate . În această configuraţie corpul are o anumită stare de tensiune şi deformaţie şi . Admitem că energia de deformaţie poate fi dezvoltată în serie Taylor după puterile lui în jurul lui : (4.2) În această serie, componentele: (4.3) reprezintă constantele elastice de ordinele 1, 2 şi 3 ale materialului aflat în configuraţia iniţială, constante ce includ efectele proceselor de formare a materialului. Din ec.(4.1) şi (4.2) rezultă acum: (4.4) Cum diferenţele şi reprezintă tensiunile şi deformaţiile raportate la sistemul de coordonate al configuraţiei iniţiale, din ec.(4.4) obţinem forma căutată a ecuaţiei constitutive: (4.5) unde: (4.6) 4.2. Relaţii constitutive între invarianţi Considerăm pământurile care pot fi modelate cu un model elastic neliniar (nisipurile) drept materiale izotrope. Este, desigur, o ipoteză simplificatoare dar acceptabilă datorită faptului că lipsa de coeziune nu a permis fixarea unor direcţii preferenţiale. Energia de deformaţie a materialelor izotrope poate fi exprimată numai în funcţie de invarianţii principali ai tensorului deformaţiei. Aceeaşi dependenţă poate avea ca argumente şi alte trei combinaţii ale invarianţilor principali, combinaţii care sunt tot invariante. Astfel, dacă folosim invarianţii sfericului şi deviatorului tensorului deformaţiei: (4.7) unde este invariantul sferic al tensorului deformaţiei, iar sunt invarianţii deviatorului tensorului deformaţiei D, din ec.(4.2) va rezulta o formă polinomială a energiei de deformaţie : (4.8) unde constantele ; ; sunt noile constante elastice de ordinele 1, 2 şi 3, constante rezultate din combinarea constantelor elastice din ec.(4.3) Această formă a energiei de deformaţie poate fi utilizată ca potenţial pentru obţinerea invarianţilor analogi ai tensorului tensiunii: (4.9) unde este invariantul sferic al tensorului tensiunii T, iar sunt invarianţii deviatorului tensorului tensiunii S . Utilizând forma (4.8) a energiei de deformaţie drept potenţial se obţine, după câteva calcule simple, relaţia constitutivă între invarianţii tensorilor tensiunii şi deformaţiei: (4.10) formă aplicabilă solicitărilor ce conduc la deformaţii cu , caz în care este necesară relaţia constitutivă . Rapoartele dintre invarianţi din parantezele ec.(4.10) pot fi exprimate în raport cu două funcţii A şi B, funcţii proporţionale cu invarianţii direcţionali şi : (4.11) din ec.(4.10) rezultând relaţia căutată: (4.12) unde : (4.13) unde constantele sunt, combinaţii ale constantelor elastice din ec.(4.7) Printr-un procedeu similar se obţine şi relaţia : (4.14) unde: (4.15) Prezenţa invarianţilor direcţionali în relaţiile constitutive (4.12) şi (4.14) pune în evidenţă dependenţa de drumul în deformaţii sau tensiuni a relaţiilor constitutive tensiune-deformaţie ce se rezumă la dependenţa dintre invarianţii de acelaşi ordin ai tensorilor tensiunii şi deformaţiei sau . S-a observat experimental că aplicând unei probe de pământ o compresiune izotropă, în deformaţiile rezultate ponderea o deţin deformaţiile de volum, dar vor apare şi deformaţii de formă, iar aplicarea unor tensiuni tangenţiale provoacă în principal deformaţii de formă, dar apar şi deformaţii de volum [10], [18], [70], [90]. Faptul că în relaţiile (4.12) şi (4.14) invarianţii tensorului tensiunii depind de toţi cei trei invarianţi ai tensorului deformaţiei denotă că aceste relaţii pot modela atât dependenţa directă sau , cât şi efectele reciproce sau , mai reduse ca pondere dar prezente. Forma polinomială a relaţiilor (4.12) şi (4.14) reflectă şi ponderea diferită a invarianţilor tensorului deformaţiei. Astfel, invarianţii de acelaşi ordin sunt folosiţi drept „argumente” ale polinoamelor, iar ceilalţi invarianţi incluşi sub forma invarianţilor direcţionali afectează „coeficienţii” dezvoltărilor polinomiale. Cu alte cuvinte, relaţiile constitutive tensiune-deformaţie, de volum (4.12) şi de formă (4.14) au aceeaşi formă analitică, diferenţa dintre ele fiind numai de ordin cantitativ, diferenţe traduse prin valorile diferite atribuite "coeficienţilor" acestor polinoame. Aceste concluzii sunt confirmate de alura curbelor experimentale tensiune-deformaţie (fig.4.1) întâlnite în literatura de specialitate ([5], [51], [52], [57], [77], [83], [84], [102] etc.) sau obţinute din experimentări proprii [8], [15], [16]. 4.3. Funcţii-modul Din relaţiile constitutive (4.12) şi (4.14) pot fi acum definite funcţiile-modul, de volum: sau de formă: , în varianta „tangentă”: (4.16) sau în forma „secantă”: (4.17) denumiri derivate din semnificaţiile geometrice ale acestor funcţii. Menţionăm că semnul minus din faţa funcţiei este impus de curbura funcţiilor experimentale (fig. 4.1). Fig. 4.1 Pentru determinarea parametrilor relaţiilor constitutive (4.13) şi (4.15) sunt necesare încercări experimentale efectuate pe diferite drumuri în deformaţii sau tensiuni, încercări care să permită evaluarea atât a dependenţei dintre invarianţii de acelaşi ordin cât şi a influenţei invarianţilor direcţionali. Însă, în triaxial, aparatul cel mai „complet” existent în ţară, nu pot fi obţinute decât drumuri situate în planul simetriei axiale, plan pentru care co, deci, la ora actuală nu putem evalua experimental influenţa acestui invariant direcţional. Acest impediment restrânge aplicabilitatea modelului la rezolvarea problemelor la limită axial-simetrice şi, cu anumite modificări, la rezolvarea problemelor plane, arie de aplicabilitate suficient de extinsă pentru a prezenta interes practic. În planul simetriei axiale relaţiile constitutive (4.13) şi (4.15) au forma: (4.18) unde: (4.19) Desigur, noile constante diferă de constantele din (4.13) şi (4.15), însă au fost notate la fel pentru simplificarea notaţiilor. Se remarcă faptul că expresiile funcţiilor-modul din ec.(4.18) pot fi considerate forme polinomiale de „moduli variabili”. Deşi beneficiază de susţinere teoretică (legată de constantele elastice de ordine superioare), formele (4.18) pot provoca uneori dificultăţi în cuantificarea experimentală datorită inflexiunilor curbelor polinomiale de aproximare a distribuţiei datelor experimentale, inflexiuni care uneori nu au suport fizic. Din acest motiv, reţinând rolul de potenţial al energiei de deformaţie dar înlocuind dezvoltarea în serie Taylor a acesteia cu un alt tip de aproximare se pot obţine funcţii-modul cu expresii analitice „netede” pe întregul interval de definiţie. Mai mult, aceste expresii analitice pot fi astfel alese încât să asigure o aproximare optimă a distribuţiei datelor experimentale. Astfel, dacă se alege pentru energia de deformaţie a formei o expresie de forma: (4.20) funcţia-modul de forfecare va avea forma: (4.21) deci, o funcţie hiperbolică de tipul celor propuse în [77], [81], [88]. De asemenea, dacă expresia aleasă pentru este: (4.22) funcţia-modul de forfecare va lua o formă exponenţială: (4.23) formă cu o bună aplicabilitate practică. 55 MODELAREA EFECTELOR REOLOGICE ÎN MEDII LINIARE 5.1. Fluaj şi relaxare Fenomene de fluaj şi relaxare sunt întâlnite şi la pământuri, în special la argile, ceea ce face ca aceste materiale să poată fi considerate vâscoelastice. Vom folosi şi în această lucrare aceeaşi denumire de material „vâscoelastic" cu menţiunea că admitem, deocamdată, existenţa proprietăţilor elastice numai în cazul primei încărcări monoton crescătoare. Însă, neliniaritatea atât de prezentă în toate manifestările acestor materiale face ca pentru pământuri să nu fie admisibile decât modelele vâscoelastic neliniare. În acest paragraf vom aminti succint câteva definiţii şi noţiuni din mecanica corpurilor vâscoelastic liniare [56], [96], [111], [131], [135], elemente care vor fi utilizate ulterior în expunerea modelului vâscoelastic neliniar. Pentru un material vâscoelasic, modificările în timp ale caracteristicilor sale mecanice reprezintă o trăsătură definitorie. Fenomenul de evoluţie în timp a deformaţiilor sub tensiuni constante este numit fluaj iar modificările tensiunilor sub deformaţii constante este cunoscut sub denumirea de relaxare. Într-o experienţă de fluaj se măsoară evoluţia în timp a deformaţiilor apărute sub o tensiune constantă. La materialele vâscoelastice obişnuite (de exemplu. polimeri) asemenea experienţe sunt uniaxiale. La pământuri vom folosi teste care conduc la relaţii uniaxiale ca formă dar rolul componentelor tensorului deformaţiei este preluat de către invarianţii tensorului deformaţiei. Astfel, experienţele de fluaj la pământurile cu proprietăţi reologice vor consta în determinarea evoluţiilor în timp a invariantului sferic al tensorului deformaţiei: sau vom urmări evoluţia în timp al invariantului deviatorului tensorului deformaţiei: . Însă, pentru fixarea unor noţiuni şi definiţii, vom folosi o experienţă uniaxială oarecare în care o deformaţie uniaxială, fie Γ, se va modifica în timp sub acţiunuea unei tensiuni constante uniaxiale oarecare, fie , unde este funcţia lui Heviside (fig.5.1-a). În momentul aplicării tensiunii apare o deformaţie elastică instantanee (OA din fig.5.1), după care deformaţia creşte în timp (arcul ABE din fig.5.1-b) apropiindu-se asimptotic, în cazul unui fluaj amortizat, de o valoare OF denumită deformaţie întârziată. Dacă, după un anumit interval de timp , tensiunea este înlăturată, bara are o revenire instantanee marcată prin micşorarea deformaţiei (BC din fig.5.1-b), urmată de o revenire progresivă (CD din fig.5.1-b). În cazul materialelor cu caracteristici perfect elastice (fără deformaţii remanente) revenirea instantanee egalează deformaţia elastică instantanee iar revenirea progresivă tinde asimptotic la zero. Cum în realitate corpurile nu sunt perfect elastice, cele două segmente OA şi BC nu sunt egale iar deformaţia de revenire progresivă nu tinde la zero ci la o valoare diferită de zero numită deformaţie permanentă (sau remanentă) [135]. Fig. 5.1 Fig. 5.2 Experienţă de fluaj Experienţă de relaxare Funcţia care descrie evoluţia în timp a deformaţiei corespunzătoare tensiunii constante aplicate este numită funcţie de fluaj (fig.5.1-c): (5.1) funcţie cu un aspect grafic asemănător funcţiei . Valoarea iniţială a funcţiei de fluaj este numită fluaj instantaneu. Într-o experienţă de relaxare se impune o deformaţie constantă (fig.5.2-a) şi se măsoară evoluţia în timp a tensiunilor necesare pentru a menţine în timp deformaţia constantă. Se remarcă faptul că tensiunea necesară scade în timp tinzând la zero sau către o anumită valoare (arcul de curbă ABC din fig.5.2-b). Dacă la un moment dat deformaţia este suprimată, tensiunea din bară scade brusc (segmentul BC din fig.5.2-b) după care tinde asimptotic către o anumită valoare numită tensiune remanentă (arcul de curbă CD din fig.5.2-b). Funcţia care reflectă variaţia tensiunii necesare menţinerii unei deformaţii constante este numită funcţie de relaxare (fig.5.2-c): (5.2) Valoarea iniţială a acesti funcţii, , este numită modul instantaneu iar valoarea asimptotică este numită modul întârziat. 5.2. Răspunsul la vibraţii forţate Experienţele de fluaj şi relaxare pot furniza informaţii asupra comportării unui material vâscoelastic supus unor solicitări lente, numite solicitări statice. Datorită vitezelor lente, efectele inerţiale devin neglijabile în acest caz, ca şi cele disipative. Experienţele efectuate în regim dinamic, de regulă armonic, oferă informaţii suplimentare utile asupra răspunsului dinamic al materialelor vâscoelastice. Admitem că un material vâscoelastic este supus unei istorii a tensiunilor de forma: (5.3) Ca şi în cazul experienţelor de fluaj sau relaxate tesnsiunea T din ec.(5.3) este o tensiune uniaxială oarecare, deformaţiile corespondente fiind notate generic cu Γ. După un interval de timp, care este cu atât mai scurt cu cât proprietăţile disipative ale materialului sunt mai puternice, proba intră în regim staţionar, unde deformaţiile de răspuns au forma: (5.4) deci, deformaţia staţionară poate fi descompusă într-o parte „în fază” cu tensiunile, parte care caracterizează proprietăţile elastice ale materialului şi o a două parte, defazată, legată de capacităţile disipative ale materialului. Analiza vibraţiilor forţate este, de regulă, condusă în domeniul complex, caz în care relaţiile dintre tensiunea T şi deformaţia Γ poate fi scrisă: ; (5.5) Menţionăm că, în general, notaţia „Re” poate lipsi, însă chiar în lipsa acesteia relaţiile complexe de forma (5.5) se referă la partea reală a mărimilor implicate (T sau Γ). Raportul: (5.6) este numit modul complex, iar raportul invers: (5.7) este numit complianţă complexă. Modulele complexe (5.6) şi (5.7) pot fi scrise şi în forma: (5.8) unde indicile „re” semnifică partea reală a mărimii respective iar indicele „im” partea imaginară. este unghiul de întârziere a materialului cuprins în intervalul: . Valoarea corespunde unui maerial elastic iar valoarea: corespunde unui fluid newtonian vâscos. În aceste două cazuri extreme valoarea unghiului δ depinde de ω. Tangenta unghiului de întârziere , numită disipare interioară are expresia: (5.9) În cazul particular în care tensiunea este dată în forma (5.3) rezultă componenta în fază a deformaţiei sub forma: iar componenta defazată cu π/2 este . 5.3. Legi de comportare liniare O lege de comportare vâscoelastic liniară de forma: (5.10) este o lege de tip fluaj, iar o lege de comportare: (5.11) este o lege de tip relaxare. Legile de comportare vâscoelasic liniare, atât de tip fluaj cât şi de tip relaxare, pot fi diferenţiale sau integrale. O lege de comportare vâscoelastică diferenţială liniară are forma: (5.12) unde: (5.13) iar operatorii diferenţiali P şi Q sunt definiţi de: (5.14) cu Dacă , legea (5.12) este de ordinul n, dacă legea este de tip relaxare iar dacă legea este de tip fluaj. De asemenea, dacă există o constantă astfel ca: (5.15) atunci materialul este de tip elastic, deoarece legea de comportare poate fi obţinută prin derivare în raport cu timpul a unei legi de tip Hooke. Materialul pentru care este numit material vâscoelastic standard: (5.16) Cele mai simple modele ce aparţin acestei clase sunt modelele Kelvin-Voigt, pentru care în ec.(5.16) se introduce: : (5.17) şi Maxwell, cu în ec.(5.16): (5.18) Legile diferenţiale de comportare pot avea „ataşate” modele analoge compuse din resorturi care traduc proprietăţile elastice şi amortizoare ce simulează proprietăţile disipative. Astfel, modelul analog Kelvin-Voigt este format prin cuplarea în paralel a unui resort elastic de rigditate cu un amortizor de vâscozitate (fig.5.3) iar modelul analog Maxwell se obţine prin legarea în serie a unui resort de rigiditate cu un amortizor de vâscozitate (fig.5.4). Fig. 5.3 - Model analog Kelvin-Voigt Fig. 5.4 - Model analog Maxwell 5.4. Modele diferenţiale de ordinul întâi 5.4.1. Modelul Kelvin-Voigt Se observă că în cazul modelului Kelvin-Voigt deformaţia resortului şi amortizorului vor fi identice, în timp ce tensiunea se va compune aditiv din două părţi: una corespunzătoare tensiunii din resort proporţională cu deformaţia resortului ( ) şi alta corespunzătoare tensiunii din amortizor proporţională cu viteza de deformaţie (). Abordând cazul unei solicitări de torsiune pură (caz în care ), ecuaţia diferenţială a modelului Kelvin-Voigt poate fi scrisă: (5.19) Prin integrarea ec.(5.19), cu condiţia iniţială se obţine evoluţia deformaţiilor în forma: (5.20) care, pentru , se reduce la: (5.21) de unde funcţia de fluaj a modelului Kelvin-Voigt rezultă (fig.5.5): (5.22) unde este valoarea asimptotică a funcţiei de fluaj iar raportul este timpul de întârziere. Soluţia ec.(5.19) pentru , corespunzătoare unui experiment de relaxare, este: , de unde rezultă că modelul Kelvin-Voigt nu are funcţie de relaxare. Modelul Kelvin-Voigt este utilizat pentru o descriere în primă aproximaţie a comportării la fluaj a unui solid vâcoelastic fără fluaj instantaneu. Fig. 5.5 - Funcţia de fluaj a Fig. 5.6 - Funcţia de relxare a modelului Kelvin-Voigt modelului Maxwell 5.4.2. Modelul Maxwell Modelul Maxwell are ca model matematic ecuaţia diferenţială (5.18): (5.23) căreia îi corespunde modelul analogic constituit dintr-un resort montat în serie cu un amortizor din fig.5.4. Deci tensiunea va fi aceiaşi atât în resort cât şi în amortizor, iar deformaţia totală va fi formată din suma deformaţiilor resortului cu cele ale amortizorului. Abordând tot cazul solicitărilor de torsiune, unde şi notând cu deformaţia din resort şi cu deformaţia amortizorului, rezultă: (5.24) De aici se obţine ecuaţia diferenţială: (5.25) Pentru o tensiune constantă ecuaţia diferenţială (5.25) devine: (5.26) din care, prin integrare, deducem: (5.27) unde constanta de integrare C se detemină din condiţia iniţială, şi anume, pentru deformaţia γ trebue să fie egală cu deformaţia elastică . Rezultă în fina: (5.28) iar apoi funcţia de fluaj a modelului Maxwell: (5.29) deci, o dreaptă de pantă 1/ cu ordonata în origine . Comportarea la relaxare a modelului Maxwell se obţine introducând în ecuaţia diferenţială a modelului (5.25) deformaţia constantă . Se obţine: (5.30) cu soluţia : (5.31) unde constanta de integrare A se determină din condiţia iniţială, adică pentru momentul tensiunea este . Rezultă astfel, tensiunile de relaxare date de modelul Maxwell în forma exponenţialei ce tinde asimptotic la zero: (5.32) şi apoi funcţia de relaxare a modelului Maxwell (fig.5.6): (5.33) unde raportul este timpul de relaxare al modelului Maxwell. Modelul Maxwell este utilizat pentru o descriere în primă aproximaţie a unui fluid vâcoelastic. 5.4.3. Solidul vâscoelastic standard Solidul vâscoelastic standard are legea diferenţială (5.16), care mai poate fi scrisă şi în forma: (5.34) cuşi condiţiile iniţiale [56]: (5.35) Ecuaţiei diferenţiale (5.34) îi poate corespunde două modele analoge diferite, unul obţinut prin montarea în serie a unui resort cu un element Kelvin-Voigt (fig.5.7-a) iar al doilea prin montarea în paralel a unui resort cu un element Maxwell (fig.5.7-b). Existenţa a două modele analoge permite ca funcţiile de fluaj şi relaxare ale solidului vâscoelastic standard să poată fi obţinute direct, fară întegrarea ecuaţiilor diferenţiale. Modelul analog din fig. 5.7-a este format prin legarea în serie a unui resort cu un element Kelvin-Voigt. Deci, deformaţia totală a montajului va fi constituită din suma deformaţiilor elementelor componente. De asemenea, şi funcţia de fluaj va fi compusă din suma: (5.36) Fig. 5.7 Modelul anlog din fig.5.7-b este constituit prin legarea în paralel a unui resort cu un element Maxwell. Deci, tensiunea totală va fi suma tensiunilor parţiale din elementele componente, ca şi funcţia de relaxare: (5.37) Ecuaţiile diferenţiale ale acestor două modele analoge sunt: (5.38) pentru montajul din figura 5.7-a, şi: (5.39) pentru modelul analog din figura 5.7-b. Comparând ecuaţiile (5.34) şi (5.38) observăm că: (5.40) iar funcţia de fluaj (5.36) poate lua forma: (5.41) unde raportul q este denumit timp de întârziere. De asemenea, din ec.(5.34) şi (5.39) rezultă: (5.42) ceea ce conduce la următoarea formă a funcţiei de relaxare: (5.43) unde parametrul este timpul de relaxare. Tinând cont de: ; (5.44) rezultă că valorile iniţiale şi la limită ale funcţiilor de fluaj şi relaxare ale solidului vâscoelastic standard sunt inverse: iar cele două funcţii pot fi scrise în forma (fig.5.8 şi fig.5.9): (5.45) Acest model permite descrierea calitativă, într-o primă aproximaţie, a comportării unui solid vâscoelastic liniar, dar având numai un timp de întârziere şi numai un timp de relaxare nu permite utilizarea sa pentru o descriere cantitativă satifăcătoare. Fig. 5.8 - Funcţia de fluaj a Fig. 5.9 - Funcţia de relaxare a solidului standard solidului standard 5.5. Modele diferenţiale generalizate Pentru o descriere mai corectă a rezultatelor obţinute în experienţele de fluaj şi relaxare au fost utilizate modelele diferenţiale – Kelvin-Voigt şi Maxwell – generalizate. Modelul Kelvin-Voigt generalizat permite o descriere mai exactă a comportării la fluaj. Modelul său analog este constituit din n elemente Kelvin-Voigt (de vâscozităţi şi rigidităţi ) montate în serie cu un element Maxwell (de vâscozitate şi rigiditate ) care poate modela elasticitatea instantanee. Funcţia de fluaj a acestui model conţine n timpi de întârziere şi are expresia: (5.46) unde . Modelul posedă, deci, n timpi de întârzire . Modelul Maxwell generalizat permite o descriere adecvată a comportării la relaxare. Modelul său analog cuprinde n elemente Maxwell (de vâscozităţi şi rigidităţi ) grupate în paralel cu un resort (de rigiditate G0). Funcţia de relaxare a acestui model este: (5.47) cu n timpi de relaxare . O extindere imediată a acestor modele generalizate poate fi obţinută considerând că materialul posedă un spectru continuu de timpi de întârziere, funcţia de fluaj (5.46) luând forma: (5.48) sau dacă materialul posedă un spectru continuu de timpi de relaxare, funcţia de relaxare (5.47) devine: (5.49) scopul cercetărilor experimentale rămânând, în acest caz, determinarea funcţiilor sau pornind de la un spectru discret de timpi de întârziere sau relaxare. 5.6. Modele vâscoelastice de tip integral Amiţând existenţa unor funcţii tensoriale de fluaj şi relaxare se pot obţine legile de comportare vâscoelastic liniare integrale de tip fluaj: (5.50) sau de relaxare: (5.51) unde prin * s-a notat convoluţia Stieltjes a funcţiilor respective: (5.52) Evident, în virtutea simetriei tensorilor tensiunii şi deformaţiei (liniare) se poate admite că funcţiile tensoriale de fluaj J şi relaxare G satisfac relaţiile de simetrie: (5.53) iar dacă J şi G sunt funcţii de clasă H1, se pot înlocui convoluţiile Stieltjes (5.50) şi (5.51) prin convoluţii Riemann: (5.54) rezultând: (5.55) Dacă materialul este izotrop, se pot separa aditiv tensorii tensiunii şi deformaţiei în sferici şi deviatori: (5.56) iar tensorii de ordinul patru J şi G, care satisfac relaţiile de simetrie (5.53) pot fi puşi sub forma: (5.57) unde , sunt funcţii arbitrare de t. Introducând (5.56) şi (5.57) în ec.(5.55) şi separând sfericii de deviatori, se obţine după câteva transformări: (5.58) Deci, în cazul izotrop este posibilă separarea legilor de comportare vâscoelastic liniare în legi care guvernează modificările de volum (între sferici) şi legi care descriu schimbarea formei (între deviatori ). Această înpărţire justifică şi denumirile de funcţie de fluaj la dilatare pentru , funcţie de fluaj de forfecare pentru şi, analog, funcţie de relaxare a volumului pentru şi funcţie de relaxare a formei pentru . Dacă în locul tensorilor sferici şi deviatori se utilizează învarianţii sfericişi deviatori (: (5.59) atunci, relaţiile constitutive tensiune-deformaţie pentru un corp vâscoelastic izotrop (5.58), iau forma: (5.60) Se remarcă faptul că legile constitutive integrale (5.60) sunt relaţii între tensiunile prezente şi viteza a tensorului deformaţiei ε, viteză ponderată de funcţiile de material, şi anume de funcţiile-modul de relaxare şi . 5.6. Răspunsul în regim armonic staţionar Să analizăm acum comportarea unui solid vâscoelastic cu ecuaţia constitutivă (5.60)2 : (5.61) supus unui regim staţionar armonic: (5.62) unde este amplitudinea iar ω frecvenţa circulară (pulsaţia) excitaţiei. În cazul unui solid vâscoelastic, cu valoarea limită a relaxării diferită de zero, funcţia de relaxare poate fi descompusă în: (5.63) unde este valoarea la limită a funcţiei de relaxare: (5.64) iar funcţia diferenţă are limita nulă: (5.65) Un exemplu de astfel de funcţie de relaxare este dat în fig.5.8-b care reprezintă funcţia de relaxare a solidului vâscoelastic standard. Introducând relaţiile (5.62) şi (5.63) în ecuaţia constitutivă (5.61) rezultă: (5.66) iar după schimbarea de variabilă t: (5.67) Relaţia constitutivă poate fi pusă în forma alternativă: (5.68) unde modulul complex are expresia: (5.69) unde, părţile: reală şi imaginară , au forma: (5.70) Partea reală este denumită modul de acumulare (storage modulus) iar partea imaginară modul de atenuare (loos modulus). Este de remarcat comportarea ecuaţiei constitutive (5.61) în cazurile extreme, la frecvenţe foarte reduse sau foarte inalte. Deoarece: (5.71) rezultă că un solid vâscoelastic cu ecuaţia constitutivă (5.61) supus unor vibraţii foarte lente sau foarte rapide răspunde ca un solid elastic. Aspectul calitativ al variaţiilor modulilor dinamici în raport cu frecvenţa ω este dat în fig.5.10 [56], [131]. Fig. 5.10 Modulul complex poate lua şi expresia alternativă: (5.72) unde este valoarea absolută a modulului complex iar este unghiul de defazaj dintre tensiuni şi deformaţii: (5.73) Funcţiile de frecvenţă (5.70): (5.74) (rescrise aici pentru a facilita compararea) pot fi considerate ca transformări Fourier directe, iar prin transformare inversă rezultă funcţiile în domeniul timp: (5.75) Se observă că dacă funcţia de relaxare este cunoscută cu ajutorul transformărilor Fourier directe (5.74) pot fi obţinute părţile: realăşi imaginară ale modulului complex , iar dacă funcţiile dinamice sunt cunoscute, transformările Fourier inverse (5.75) pot furniza funcţia de relaxare. Însă, relaţiile (5.74) şi (5.75) sunt interconectate pe întregul lor domeniu de definiţie. Prin urmare obţinerea unor valori corecte pentru funcţia de relaxare ar necesita cunoaşterea modulului complex pe întreg spectrul de frecvenţe iar determinarea corectă a valorilor modulului complex necesită cunoaşterea funcţiei de relaxare pe întregul domeniu de timp. Cum experimental determinările sunt accesibile numai unor domenii restrânse este recomandabilă utilizarea unor teste specifice – de relaxare pentru funcţia de relaxare şi teste dinamice pentru modulul complex. 66 MODEL VÂSCOELASTIC NELINIAR 6.1. Forma generală a ecuaţiei constitutive Vectorul deplasărilor u, tensorii tensiunii T şi deformaţiei E sunt, în cazul materialelor vâscoelastic liniare, funcţii de poziţie x şi de timp t, funcţii ce definesc starea vâscoelastică a corpului. Pentru un timp dat şi fixat aceste funcţii vor defini o stare elastică a corpului. Această reducere a stărilor vâscoelastice la stări elastice - intuitivă, de fapt - este riguros definită şi demonstrată în [72], unde este dată şi condiţia necesară acestei reduceri, şi anume, condiţia ca pentru un timp t fixat, funcţiile de relaxare şi de fluaj să se reducă la constantele şi coeficienţii elastici corespunzători. Reducerea stărilor vâscoelastice la stări elastice este observabilă experimental şi în comportarea unor probe de argilă supuse în triaxial unei încercări de tip fluaj, izocronele şi/sau fiind curbe tensiune-deformaţie ce pot fi modelate cu un model elastic neliniar. Modelul prezentat în acest capitol se bazează pe reducerea stărilor vâscoelastice la stări elastice. Pornind de la forme generale ale ecuaţiei constitutive ale vâscoelasticităţii neliniare se ajunge la forme aplicabile la pământuri utilizând condiţia ca pentru un moment dat funcţiile neliniare de relaxare şi să se reducă la funcţiile-modul elastic neliniare şi din capitolul 4. În acest capitol va fi expus, de fapt, un procedeu de determinare a variaţiei în timp a acestor izocrone, ataşând astfel relaţiilor tensiune-deformaţie elastic neliniare o lege de variaţie în timp. În cazul materialelor vâscoelastic liniare, pentru un timp dat , funcţiile de relaxare şi se reduc la constantele elastice corespunzătoare K şi G iar relaţiile constitutive la relaţiile constitutive elastice corespunzătoare: (6.1) Admitem că această teoremă de reducere a stărilor vâscoelastice la cele elastice, demonstrată de Gurtin şi Sternberg pentru corpurile liniare [72], poate fi extinsă la materialele nelinare în sensul că pentru un timp dat , funcţiile de relaxare şi se reduc la funcţiile-modul corespunzătoare şi . Prin urmare, funcţiile de relaxare a unui material vâscoelastic neliniar pe lângă variabila timp trebuie să conţină drept argumente şi invarianţii tensorului deformaţiei: şi . În aceste condiţii, ecuaţiile constitutive vâscoelastic neliniare (5.60) iau forma: (6.2) 6.2. Funcţiile neliniare de relaxare Pentru ca legea vâscoelstic neliniară (6.2) să fie complet determinată trebuie ca funcţiile de material – funcţiile neliniare de relaxare şi să fie definite şi cuantificate prin teste experimentale de relaxare. În cazul pământurilor pot fi folosite în acest scop testele triaxiale. Însă, deoarece testele de relaxare sunt dificil de realizat în triaxial vom prezenta în continuare o metodă de determinare indirectă a funcţiilor de relaxare prin teste triaxiale de fluaj. Vom exemplifica metoda prin determinarea funcţiei de relaxare utilizând datele experimentale obţiute prin teste triaxiale pe probe de argilă [23]. Metoda este compusă din în două etape. În prima etapă se determină experimental funcţii de fluaj iar apoi, în etapa a doua, prin prelucrarea numerică a acestora se obţin funcţiile de relaxare căutate. Într-un experiment triaxial de fluaj o probă de pământ este solicitată de forţe cunoscute aplicate în n trepte de încărcare, fiecare treaptă i fiind menţinută constantă un anumit interval de timp . Deci, la fiecare treaptă proba va avea o stare de tensiune constantă iar deformaţiile corespunzătoare vor creşte, până la stabilizare sau rupere. Se realizează, astfel, la fiecare treaptă o experienţă de fluaj din datele căreia poate fi obţinută variaţia deformaţiilor treptei , adică funcţii de fluaj cu o formă analitică impusa de distribuţia datelor experimentale. De regulă, funcţia de fluaj a solidului vâscoelastic standard de formă exponenţială: (6.3) reprezintă o bună aproximaţie pentru fluajul unei trepte. In ec.(6.3) este valoarea invariantului deformaţiei treptei i la momentul iniţial al treptei, reprezintă valoarea asimptotică de stabilizare iar prin s-a notat timpul de întârziere al treptei (fig.6.1). Fig. 6.1 Acest procedeu este aplicabil în procesele de fluaj cu tendinţă de stabilizare, caz în care deformaţiile scad în timp (fluaj primar [56]). După un anumit număr de trepte, prin creşterea solicitărilor exterioare, viteza de deformaţie devine crescătoare iar materialul intră în fluaj accelerat către rupere (curgere). După determinarea funcţiilor de fluaj se poate trece la etapa următoare care constă în determinare unor izocrone folosind valorile τ şi γ culese din toate funcţiile de fluaj ale treptelor pentru un acelaşi timp t fixat. În acest scop, fiecare din intervalele de timp ale celor n trepte de timp se divizează în aceleaşi m trepte de timp, stabilindu-se la fiecare treaptă de încărcare i valori . Apoi, pentru un acelaşi timp tj, din toate cele n funcţii de fluaj ale celor n trepte de încărcare se pot obţine m relaţii constitutive elastic neliniare de forma (fig.6.2) : (6.4) Fig. 6.2 Cele m funcţii-modul ale acestor izocrone, de forma unor ecuaţii constitutive elastic neliniare (fig.6.3): (6.5) pot fi interpretate drept curbe pentru valori fixe ale variabilei timp pe suprafaţa generată de funcţia . Formele analitice ale relaţiilor constitutive (6.4) ca şi ale funcţiilor-modul (6.5) sunt impuse de distribuţia datelor experimentale. Pot fi utilizate în acest scop funcţii hiperbolice, exponenţiale, polinomiale, sigmoidale etc., aproximaţiile oferite nefiind mult diferite (un exemplu este dat în fig.6.4). Fig. 6.3 Fig. 6.3 Fig. 6.4 În final, din relaţiile (6.5) se obţin valori ale funcţiei-modul de relaxare: (6.6) din care, prin prelucrare statistică, se pot obţine expresii analitice pentru funcţia de relaxare căutată. O astfel de formă pentru funcţia de relaxare este: (6.7) unde, valoarea iniţialăeste constantă în raport cu γ şi t dar poate fi dependentă de alţi parametri fizici şi mecanici ai materialului, cum sunt: umiditatea, densitatea, indicile porilor, timpul de consolidare a probei sau raportul de supraconsolidare, presiunea de consolidare, etc. Pentru argila folosită ca exemplu a fost obţinută o funcţie de relaxare (fig.6.5), de forma ec.(6.7), unde funcţiile parţiale de deformaţie şi timp sunt: (6.8) Fig. 6.5 77 MODELAREA RĂSPUNSULUI DINAMIC 7.1. Model dinamic vâscoelastic neliniar Adimitem că un material vâscoelastic neliniar cu ecuaţia constitutivă (6.2)2 supus unei excitaţii armonice staţionare: (7.1) răspunde în deformaţii tot sub formă armonică staţionară: (7.2) În relaţiile de mai sus este amplitudinea excitaţiei torsionale, ω frecvenţa circulară (pulsaţia), amplitudinea răspunsului în deformaţii iar . Introducând în ecuaţia constitutivă (6.2)2 istoria deformaţiei (7.2) se obţine o ecuaţie constitutivă dinamică de forma[19], [22], [46]: (7.3) unde este funcţia-modul complex: (7.4) astfel denumită prin analogie cu dinamica liniară. Tot prin analogie vom denumi partea reală a funcţiei : , funcţie-modul de acumulare (storage modulus function), iar partea imaginară , funcţie-modul de atenuare (loss modulus function) : (7.5) iar este valoarea la limită: (7.6) Utilizând forma (6.7) a funcţiei neliniare de relaxare: (6.7) funcţiile de acumulare şi atenuare devin: (7.7) unde funcţia de deformaţie are aceiaşi formă ca la funcţia neliniară de relaxare (6.7) iar dependenţele de frecvenţă iau forma: (7.8) Este de remarcat comportarea ecuaţiei constitutive (7.3) în cazurile extreme, la frecvenţe foarte reduse sau foarte înalte. Deoarece: (7.9) rezultă că un solid vâscoelastic cu ecuaţia constitutivă neliniară (7.3) supus unor procese foarte lente sau foarte rapide răspunde ca un solid elastic neliniar. Primul caz, când , corespunde solicitărilor „statice” aplicate un timp suficient de îndelungat pentru ca fenomenele de relaxare şi fluaj să poată fi considerate „stinse”, deci este firesc ca în ec.(7.3) să intervină valorile stabilizate ale funcţiilor neliniare de relaxare din modelul vâscoelastic neliniar. În cazul în care ω→ datorită frecvenţei foarte mari a excitaţiei relaxarea este blocată iar în ec.(7.3) intervin valorile iniţiale ale funcţiilor neliniare de relaxare. Regăsim astfel, atât asemănarea cu dinamica liniară [56], [131], cât şi concordanţa cu observaţiile experimentale [6], [22], [95]. Funcţia-modul complex (7.4) poate lua forma alternativă: (7.10) unde: (7.11) Tot prin analogie cu dinamica liniară prima funcţie (7.11), numită funcţia-modul de torsiune dinamică, introduce în model rigiditatea iar cea de a doua funcţie (7.11), numită funcţia de amortizare torsională, caracterizează proprietăţile disipative ale materialului. Relaţiile inverse celor din ec.(7.11) sunt: (7.12) De regulă, pentru materialele cu amortizări reduse: . Datele experimentale au confirmat că această aproximare poate fi admisă şi la pământuri. Chiar în cazul argilei testate, pentru care în coloana rezonantă s-a obţinut o valoare maximă a raportului de amortizare de 0,22, valoarea maximă a funcţiei de sub radicalul din ec.(7.12) este de 1,10454, Deci, acceptarea aproximaţiei conduce la o eroare maximă de 9%, ceea ce în cazul pământurilor este acceptabilă. Utilizând expresiile (7.5) ale părţilor reală şi imaginară ale funcţiei-modul complex, funcţia-modul de torsiune dinamică şi funcţia de amortizare torsională devin: (7.13) Să examinăm acum aspectul calitativ al acestor funcţii. Dependenţa de frecvenţă a funcţiilor şi (fig.7.1) este în concordanţă cu comportarea calitativă a unui solid vâscoelastic şi cu observaţiile experimentale obţinute prin teste pe pământuri [56], [32]. De asemenea, dependenţa funcţiei-modul dinamic în raport cu nivelul de deformaţie este corectă – pentru o frecvenţă dată această funcţie dinamică se reduce la funcţia-modul elastic neliniară. Însă, funcţia de amortizare din (7.13) depinde numai de frecvenţă ceea ce implică o amortizare constantă în raport cu deformaţia, formă neadcvată modelării amortizării în pământuri. Datorită existenţei şi acumulărilor deformaţiilor ireversibile, disiparea de energie în pământuri este dependentă de amplitudinea deformaţiei, ceea ce impune modificarea formei analitice a funcţiei de amortizare. Acestă neconcordanţă calitativă se datorează formei (6.7) cu variabile separate admisă pentru funcţia neliniară de relaxare , formă utilă în prelucrarea datelor experimentale şi care nu a dat neconcordanţe cu distribuţia datelor experimentale. Însă, în relaţia (7.13) a condus la eliminarea funcţiei din raportul . Fig. 7.1 Din acest motiv, vom utiliza o altă formă a funcţiei de amortizare, obţinută ca o extindere în domeniul neliniar a raportului de amortizare din dinamica liniară, raport dependent atât de frecvenţă cât şi de nivelul de deformaţie: (7.14) dependenţe care au fost admise separabile, ca şi la funcţia-modul torsional. După cum se va observa mai departe, în paragraful 7.4, această formă este validată de experiment. Modelul vâscoelastic neliniar dinamic prezentat în acest paragraf a fost construit prin adaptarea controlată şi amendată experimental a modelului vâscoelastic neliniar bazat pe funcţia neliniară de relaxare obţinută prin teste triaxiale. Se poate spune că modelul dinamic este până acum tributar unor determinări experimentale în regim static şi nu dinamic cum ar fi de aşteptat. Desigur, între valorile pentru un anumit nivel γ ale funcţiilor de relaxare şi cele ale funcţiei-modul complex există o interdependenţă dată de transformările Fourier directe (5.73): (7.15) sau de transformările inverse (5.74): (7.16) Însă, aşa cum se observă din aceste relaţii de transformare, întreaga plajă de frecvenţe este pusă în corespondenţă cu întreaga plajă de timp ceea ce face extrem de dificilă determinarea experimentală. În plus, deşi modelul vâsoelastic neliniar obţinut prin teste experimentale statice s-a dovedit util pentru determinarea formei complexe a ecuaţiei constitutive dinamice, forma funcţiei de amortizare a trebuit modificată datorită inerentelor simplificări impuse de prelucrările datelor experimentale. Toate aceste motive fac recomandabilă menţinerea formei complexe a ecuaţiei constitutive, dar determinarea funcţiilor ce intervin în această formă - şi trebuie făcută nu prin utilizarea relaţiilor de transformare ci direct, prin teste dinamice în aparate destinate încercărilor dinamice cum sunt: coloana rezonantă [29], [76], triaxialul dinamic [79], [113], aparatul de torsiune dinamică [129] etc. În cele ce urmează, vom prezenta o astfel de determinare a funcţiilor dinamice prin utilizarea testelor în coloana rezonantă Drnevich [142] pe probe din aceiaşi argilă folosită şi la testele triaxiale de determinare a funcţiei neliniare de relaxare, determinare prezentată în capitolul 6. Prezentarea experimentului în coloana rezonantă a fost făcută în capitolul 3. Reamintim, că pentru un anumit nivel al excitaţiei torsionale , din datele furnizate de aparat se obţin valori ale frecvenţei ω, deformaţiei torsionale γ, ale modului torsional G şi ale raportului de amortizare D, valori corespunzătoare nivelului de deformaţie torsională γ . Modificând amplitudinea excitaţiei, se inpune probei un alt nivel de deformare γ căruia îi corespund alte valori G şi D. După n astfel de trepte de încărare se obţine un set de valori ale funcţiei-modul dinamic , şi un set de valori ale funcţiei de amortizare , din care prin prelucrare statistică se obţin funcţiile respective şi . Rezultatele sunt prezentate în fig.7.2 şi fig.7.3, unde datele furnizate de coloana rezonantă (la frecvenţe de peste 1 Hz) au fost completate pentru domeniul frecvenţelor joase prin încercări ciclice triaxiale. Fig. 7.2 Aplicând acelaşi procedeu statistic de determinare a dublei dependenţe (în raport cu deformaţia şi frecvenţa ω), procedeu prezentat în capitolul 5 cu ocazia determinării dublei dependenţe a funcţiei de relaxare (în raport cu deformaţia şi timpul t) şi admiţând separarea variabilelor s-a obţinut: = 0,001 Hz = 0,01 Hz = 0,1 Hz = 1 Hz = 10 Hz = 100 Hz (7.17) (7.18) Fig. 7.3 Aspectul calitativ al acestor funcţii dinamice este dat în următoarele diagrame spaţiale, figura 7.4 pentru funcţia-modul dinamic şi figura 7.5 pentru funcţia de amortizare Fig. 7.4 Fig. 7.5 7.2. Model dinamic pentru solicitări seismice Testele din coloana rezonantă indică faptul că ponderea majoră în variaţia valorilor funcţiilor dinamice şi o deţine variabila γ0, iar frecveţa ω are un aport neglijabil în plaja de frecvenţe de peste 1Hz, plajă de interes în ingineria seismică. Această comportare a fost semnalată în lucrări anterioare [22], [67], [74], [77], [81], [99], [123], [129], [136] şi poate fi observată şi în figurile 7.6 şi 7.7. Menţionăm că independenţa de ω a funcţiilor dinamice de material nu înseamnă independenţa de ω a răspunsului structural. Fig. 7.6 Deci, pentru frecvenţe dependenţa funcţiilor dinamice de ω poate fi neglijată, caz în care funcţia-modul complex (7.4) devine: (7.19) unde: (7.20) iar legea constitutivă (7.3) ia forma : (7.21) formă care poate fi privită ca o extindere în domeniul neliniar a modelului Kelvin-Voigt nevâscos [9]. Partea reală a ec.(7.21): (7.22) poate fi considerată o ecuaţie de tip elastic neliniar, ecuaţia curbei-schelet. De asemenea, funcţia de amortizare este legată de energia disipată de o buclă de histerezis de amplitudine cu relaţia cunoscută [80], [92]: (7.23) unde este energia potenţială acumulată până la amplitudinea iar este energia disipată de ciclul cu aplitudine . Fig. 7.7 Neglijând aportul frecvenţelor de peste 1 Hz, funcţia-modul din ec.(7.13) şi funcţia de amortizare din ec.(7.13) şi funcţia de amortizare din ec.(7.14) iau forma cunoscută: (7.24) unde indicele n semnifică forma normalizată a funcţiei în raport cu valoarea de normalizare notată cu indicele 0. De asemenea, funcţiile dinamice (7.17) şi (7.18) folosite ca exemplu de determinare a formei complete a ecuaţiei constitutive se reduc la: ; (7.25) Forma normalizată a funcţiilor dinamice are aspectul din fig.7.8. In această figură, ca şi în relaţiile şi figurile următoare, s-a renunţat la indicele "0", astfel că notaţia γ semnifică în cele ce urmază amplitudinea deformaţiei torsionale. Fig. 7.8 7.3. Model Kelvin-Voigt neliniar Modelul vâscoelastic neliniar cu caracteristici dinamice de material independente de frecvenţă din paragraful anterior 7.2, model bazat pe două funcţii de material – funcţia-modul dinamic care modelează rigiditatea materialului şi funcţia de amortizare legată de proprietăţile disipative ale materialului sugerează posibilitatea utilizării unui model vâscoelastic analog. Un asemenea model, construit pentru un solid vâscoelastic standard (deci, liniar) este şi modelul Kelvin-Voigt format din legarea în paralel a unui resort de rigiditate G cu un amortizor cu vâscozitatea η. Acum, paralela dintre modelul vâscoelastic neliniar şi modelul analog Kelvin-Voigt este evidentă. Un model analog neliniar poate fi realizat prin înlocuirea rigidităţii liniare G cu una neliniară bazată pe funcţia-modul dinamic şi a vâscozităţii liniare η cu o caracteristică de amortizare neliniară legată de funcţia de amortizare . Utilizând această analogie formală poate fi construit unui model analog Kelvin-Voigt neliniar cu caracteristici de rigiditate şi de amortizare dependente de nivelul de deformaţie sau de mărimea deplasării. Şi în acest caz vom porni de la date experimentale obţinute printr-un test în coloana rezonantă pe o probă de argilă, pentru care s-a obţinut experimental următoarele funcţii dinamice de material (fig.7.11 şi fig.7.12): Fig. 7.11 Fig. 7.12 (7.26) Asemănarea relaţiilor (7.26) şi a graficelor corespunzătoare 7.11 şi 7.12 cu funcţiile de relaxare şi fluaj a solidului vâscoelastic standard (relaţiile 5.45 şi figurile 5.8) este acum evidentă. Deci, o astfel de comportare a materialului sugerează că un model analog al solidului vâscoelastic standard (dar neliniar) este mai indicat decât un model Kelvin-Voigt neliniar. Însă, prin înlocuirea constantelor din modelul liniar Kelvin-Voigt cu funcţiile (7.26), care conţin constante libere cu rol de valori iniţiale, funcţionarea modelului Kelvin-Voigt neliniar devine identică cu a unui model standard neliniar iar diferenţele se menţin numai în alcătuirea grafică a modelelor analoge. Cum în acest context, diferenţa nu prezintă importanţă deoarece nu vom utiliza modelul în forma sa analogă, în cele ce urmează vom utiliza denumirea de model Kelvin-Voigt neliniar. Un astfel de model are şi un corespondent fizic în aparatul de încercări dinamice numit coloana rezonantă. Sistemul oscilant al coloanei rezonante este un sistem cu un grad de libertate dinamică format din probă şi vibratorul ataşat (fig.7.13) iar răspunsul dinamic structural depinde atât de caracteristicile dinamice ale materialului din probă cât şi de condiţiile concrete din aparat (geometria probei şi a vibratorului, tipul de solicitare dinamică, condiţiile limită şi de de capăt etc.). Această corespondenţă oferă avantajul unei verificări a modelului Kelvin-Voigt neliniar şi, implicit, a modelului vâsoelastic neliniar din care provine. Fig. 7.13 Curbele tensiune-deformaţie obţinute prin teste ciclice au forma unor bucle închise, ceea ce conduce la concluzia că amortizarea în pământuri are mai curând un caracter histeretic şi nu este de tip vâscos ca la modelul Kelvin-Voigt. Toate determinările din coloana rezonantă a capacităţii de amortizare corespunzătoare unui anumit nivel al excitaţiei sunt bazate pe ipoteza de echivalenţă dintre amortizarea histeretică a probei de pământ şi amortizarea vâscoasă a unei probe cu aceiaşi masă, densitate şi dimensiuni dar formată dintr-un material echivalent vâscos. În virtutea acestei ipoteze reo-histeretice modelul Kelvin-Voigt liniar este capabil de a descrie disiparea de energie din sistemul oscilant probă-vibrator pentru un anumit nivel al excitaţiei iar verificările au indicat o bună concordanţă model-experiment. O verificare suplimentară a ipotezei reo-histeretice poate fi realizată acum prin analiza capacităţii modelului Kelvin-Voigt neliniar de a modela buclele de histereză şi prin compararea funcţiilor de amortizare obţinute cu modelul Kelvin-Voigt neliniar şi cu modele histeretice directe. O serie de metode de determinare a „raportului de amortizare” D, (aşa cum este numit în dinamica liniară) pentru un anumit nivel de deformaţie sunt bazate pe relaţia cunoscută: (7.27) unde W este energia acumulată maximă iar ∆W este energia disipată pe un ciclu, reprezentată de aria inclusă în bucla de histerezis (fig.7.14). Unele metode folosesc buclele de histereză înregistrate experimental iar alte metode construiesc aceste bucle pornind de la curba-schelet şi utilizând regula lui Masing: ramura superiaoră şi cea inferioară a buclei de histerezis sunt obţinute din curba-schelet prin reducerea cu factorul doi în ambele direcţii [80]. Fig. 7.14 Pentru a testa capacitatea modelului Kelvin-Voigt de a modela bucle de histeresis, dar mai ales caracteristicile de amortizare a materialului putem apela la un „drum invers” – plecând de la funcţiile de material şi determinate experimental se calculează „forţa de revenire” şi se obţine cu ajutorul ei funcţiile care definesc bucla de histereză. Apoi apelând la relaţia de determinare a amortizării cu ajutorul ariei buclei de histereză (7.27) se determină valoarea funcţiei de amortizare corespunzătoare unei amplitudinii excitaţiei . Această valoare a funcţiei de amortizare este apoi comparată cu valoarea obţinută experimental în coloana rezonantă pentru acelaşi nivel al excitaţiei. Menţionăm că în coloana rezonantă valorile amortizării sunt determinate cu metode diferite, care nu utilizează relaţia (7.27), şi anume metode bazate fie pe înregistrarea oscilaţiei libere amortizate după „tăierea” excitaţiei, fie pe metoda factorului de amplificare (vezi, capitolul 3). Asfel, pentru modelul Kelvin-Voigt neliniar din fig.7.13 forţa de revenire este: (7.28) unde componenta elastică a forţei de revenire este (ecuaţia curbei-schelet), iar componenta de amortizare a forţei de revenire este . Pentru un anumit nivel al excitaţiei armonice , după stingerea oscilaţiei tranzitorii rotaţia de răspuns are forma : (7.29) Ţinând cont că: (7.30) prin eliminarea variabilei timp t din ecuaţiile (7.29) şi (7.30) rezultă: (7.31) iar forţa de revenire (7.28) rămâne cu o singură variabilă θ: (7.32) Funcţia (7.32) descrie într-un plan(bucla de histerezis corespunzătoare nivelului amplitudinii : cu semnul „+” va reprezenta ramura superioară a buclei iar cu semnul „ ” pe cea inferioară (fig.7.15). Pentru comparare, în fig.7.16 este dată bucla de histerezis construită cu ajutorul regulei lui Masing preluând din exemplul dat în fig.7.15 aceiaşi curbă-schelet şi aceiaşi amplitudine rad. Aşa cum rezultă din compararea celor două figuri (7.15 şi 7.16), aspectul geometric obţinut în cele două cazuri este diferit. Însă, valoarea raportului de amortizare D este dependent de aria buclei de histerezis şi nu de forma ei geometrică. Din acest punct de vedere, diferenţele persistă dar nu mai sunt atât de evidente. Calculând în cele două cazuri valori ale raportului de amortizare corespunzătoare unor aceleaşi amplitudini : , se obţin două funcţii de amortizare care nu sunt mult diferite între ele şi nici depărtate de curba determinată pe baza datelor experimentale obţinute direct în coloana rezonantă. (fig.7.17). Fig. 7.15 Fig.7.16 Fig.7.17 7.4. Validarea modelului vâscoelastic neliniar Validarea corectitudinii modelării comportării dinamice a materialului poate fi realizată dacă există posibilitatea comparării răspunsului dinamic structural calculat pe baza modelului cu răspunsul structural măsurat experimental. Cum asemenea experimente cu structuri la scară mare sunt dificil sau chiar imposibil de realizat preconizăm utilizarea în acest scop a unor experimente de laborator unde structura-test să fie constituită din sistemul oscilant al coloanei rezonante. Avantajul unei astfel de alegeri, pe lângă costul redus, constă în posibilitatea unui control riguros al răspunsului structural măsurat. Sistemul oscilant al coloanei rezonante este un sistem cu un grad de libertate dinamică format din probă şi vibratorul ataşat. Datorită materialului din probă, în toate cazurile, ecuaţia de mişcare a sistemului oscilant probă-vibrator va fi una neliniară. Pentru o excitaţie armonică staţionară de forma sistemul poate fi modelat cu un modelul Kelvin-Voigt neliniar din paragraful anterior (fig.7.13). În acest caz, ecuaţia de mişcare a sistemului neliniar poate fi scrisă în forma [4]: (7.33) unde este momentul de inerţie al vibratorului. Utilizând aceiaşi metodă de adimensionalizare prin schimbarea de variabilă rezultă ecuaţia de mişcare (vezi cap.2): (7.34) unde indicii superiori reprezintă derivata în raport cu noua variabilă de timp, iar: (7.35) Pentru o amplitudine normalizată µ dată şi o anumită pulsaţie relativă υ ecuaţia diferenţială neliniară (7.34) poate fi rezolvată numeric şi o soluţie de forma poate fi obţinută în n puncte. După eliminarea porţiunii iniţiale tranzitorii, reţinând numai soluţia staţionară se obţine amplitudinea rotaţiei (fig.7.18) : (7.36) Aceiaşi amplitudine poate fi direct obţinută din datele de ieşire ale coloanei rezonante: (7.37) unde A este valoarea acceleraţiei, ra este distanţa de la axa de rotaţie la axa accelerometrului (ra = 0,03175 m la coloana rezonantă Drnevici) iar pulsaţia semnalului de ieşire [142]. Compararea valorilor obţinute prin cele două căi distincte – direct din datele de ieşire a aparatului şi indirect prin calcul din soluţia sistemului oscilant – poate da o imagine a corectitudinii modelării. O astfel de comparaţie este dată în fig.7.19. Fig. 7.18 Fig.7.19 Date experimentale Fiting al datelor experimentale Rezolvări ale ecuaţiei neliniare Fiting al rezolvărilor neliniare 88 MODEL VÂSCOPLASTIC În capitolele anterioare au fost prezentate modelări ale principalelor caracteristici ale comportării mecanice, cum sunt : neliniaritatea, capacitatea disipativă, efectele reologice, toate inglobate, practic, într-un model vâscoelastic neliniar adaptabil atât solicitărilor statice cât şi celor dinamice. Capitolul de faţă are ca scop includerea în mod explicit a caracterului ireversibil al deformaţiilor şi a influenţei vitezei de deformare, ceea ce va conferi modelului un caracter vâscoplastic. Ne vom concentra atenţia asupra stabilirii unei ecuaţii constitutive vâscoplastice bazate pe utilizarea informaţiilor ce se pot obţine prin combinarea unor teste definitorii executate în triaxial şi coloana rezonantă. 8.1. Forma generală a ecuaţiei constitutive Adoptăm ca formă generală a ecuaţiei constitutive o formă directă între tensiuni şi deformaţii: (8.1) formă în care, influenţa vitezei de deformare a fost inclusă ca argument al funcţiei-modul torsional G. O asfel de formă a ecuaţiei constitutive vâscoplastice poate deveni operaţională cu condiţia ca să poată fi complet definită, calitativ şi cantitativ, prin teste experimentale. După cum se va vedea în cele ce urmează acestă condiţie poate fi satisfăcută utilizând atât teste lente în triaxial cât şi încercări dinamice în coloana rezonantă. 8.2. Determinare prin încercări triaxiale Apelând la toată gama de încercări triaxiale, cu viteze de deformare diferite, plecând de la cele mai lente cum sunt testele de fluaj şi continuând cu încercările monoton crescătoare pe probe din acelaşi material, pot fi obţinute curbe tensiune-deformaţie care diferă între ele numai datorită vitezei de deformare impusă experimental. Prelucrarea datelor experimentale oferite de aceste teste poate fi realizată în două etape, mai întâi determinând pentru fiecare viteză de deformare o ecuaţie tensiune-deformaţie cu aceiaşi formă analitică , iar apoi prin determinarea dependenţei în raport cu viteza de deformaţie a „parametrilor” funcţiilor analitice determinate în prima etapă. Poate fi astfel, evaluată calitativ şi cantitativ dubla dependenţă, atât în funcţie de nivelul de deformaţie atins cât şi în raport cu viteza cu care s-a ajuns la acest nivel. Vom exemplifica o astfel de determinare prin prelucrarea datelor experimentale obţinute pe probe de argilă în teste triaxiale, de fluaj şi monoton crescătoare. Astfel, în figurile 8.1 – 8.3 sunt date rezultatele a trei încercări monoton crescătoare cu viteze diferite. Impunând pistonului viteze de înaintare de 5, 1 şi 0,1 mm/min şi ţinând cont de dimensiunile geometrice ale probei vitezele de deformare impuse au fost de aproximativ 10-3, şi, respectiv 10-5 s-1. În fig.8.4 este dată curba fluajului stabilizat, adică curba obţinută de teste de fluaj cu tensiuni τ = ct. care au condus la stabilizarea în timp a deformaţiilor. Fig. 8.1 Fig. 8.2 Fig. 8.3 Fig. 8.4 Din aceste figuri se poate observa că o exponenţială de forma: (8.2) reprezintă o bună aproximare în toate cele 4 experimente, aliura calitativă este corespunzătoare iar ceea ce diferă sunt numai valorile parametrilor a şi b. Prelucrând aceste valori a şi b corespunzătoare vitezelor la care au fost obţinute se pot obţine funcţii , cu o formă impusă de distribuţia valorilor a şi b obţinute în pasul anterior. O astfel de prelucrare este dată în fig.8.5, în care au fost folosite valorile a şi b din fig.8.1 – 8.4. Determinarea calitativă şi cantitativă a funcţiilor şi conduce la explicitarea completă a ecuaţiei constitutive (8.1) în forma : (8.3) relaţie care pentru argila testată are aspectul din fig.8.6. De asemenea, din ec.(8.3) rezultă şi funcţia-modul torsional , al cărei aspect este dat în fig.8.7. Fig. 8.5 Fig. 8.5 Fig. 8.6 Fig. 8.7 8.3. Determinare prin încercări dinamice In coloana rezonantă o probă este solicitată în regim armonic: (8.4) ceea ce ne permite să admitem că şi deformaţia provocată are acelaşi aspect: (8.5) rezultând, astfel, că viteza de deformaţie impusă probei are forma: (8.6) Prin măsurarea vitezei de propagare a undelor prin probă experimentul furnizează, printre alte date utile, şi valoarea amplitudinii deformaţiei atinse de probă : (8.7) unde A este valoarea măsurată în mvrms a acceleraţiei mişcării iar fr este frecvenţa de rezonanţă a sistemului probă-vibrator. De aici rezultă că amplitudinea vitezei de deformaţie este: (8.8) Pentru o anumită valoare a amplitudinii momentului perturbator M0, în regim de rezonanţă, experimentul furnizează valorile corespunzătoare ale funcţiei-modul torsional G, deformaţiei γ şi vitezei de deformaţie . Repetând aceste determinări, pentru diferite valori M0, datorită cunoscutei comportări neliniare, se obţin alte valori G, γ şi , rezultând în final funcţii de forma (fig.8.8): (8.9) (fig.8.9): (8.10) sau, direct (fig.8.10), unde a fost utilizată ipoteza separării variabilelor: (8.11) Fig. 8.8 Fig. 8.9 Fig.8.10 Astfel, relaţia constitutivă (8.1) este complet determinată, aspectul ei din spaţiul fiind dat în fig.8.11. Fig.8.11 8.4 Comparaţie între răspunsul lent şi cel dinamic Comparând aspectul răspunsului în tensiuni determinat prin încercări lente (fig.8.6) cu răspunsul dinamic (fig.8.11) observăm diferenţe calitative în ceea ce priveşte variaţia în raport cu viteza de deformaţie. În cazul încercărilor lente tensiunile τ cresc în raport cu viteza de deformaţie în timp ce la solicitări dinamice tensiunile τ se micşorează pe măsură ce viteza de deformaţie este crescută. De asemenea, prin compararea funcţiilor-modul torsional obţinute prin încercări lente (fig.8.7) sau prin încercări dinamice în coloana rezonantă (fig.8.10) se observă că în timp ce în triaxial valorile G cresc cu la orice nivel de deformare γ, în coloana rezonantă valorile G scad în raport cu . O explicaţie plauzibilă a acestei comportări aparent bizare se obţine considerând fenomenul de degradare dinamică [25], [45] al cărui efect de diminuare a rezistenţelor mecanice pare că depăşeşte rigidizarea dată de creşterea vitezei de deformare. Experimental, în încercările triaxiale monoton crescătoare, se constată că modificarea vitezei de încărcare conduce la modificări în aspectul curbelor tensiune-deformaţie dar şi la obţinerea unor valori diferite ale rezistenţei de rupere [31]. Pe măsură ce viteza de încărcare creşte, procesele de deformare în timp datorate proprietăţilor reologice sunt stânjenite tot mai mult, iar rezistenţele de rupere obţinute au valori din ce în ce mai mari. În cazul încărcărilor aplicate în trepte, efectele reologice sunt şi mai evidente. Fiecare treaptă de încărcare menţinută un interval de timp dat constituie, de fapt, o încercare de fluaj în care deformaţiile evoluează până la stabilizare sau rupere. Este evident faptul că rezistenţele de rupere obţinue în acest tip de încercare vor avea valori mai reduse. Regăsim, astfel, o constatare experimentală anterioară – modulele dinamice au valori superioare celor statice. În cazul solicitărilor repetate sau ciclice viteza de încărcare prezintă aceiaşi importanţă. De regulă acest tip de solicitare presupune existenţa unor solicitări aplicate într-un mic interval de timp pe ciclu, ceea ce conduce la viteze de încărcare sau deformare mari. Ca şi în cazul solicitărilor monoton crescătoare, cu cât este mai mare viteza de încărcare cu atât fenomenele reologice sunt mai stânjenite, fapt ce se traduce prin creşterea rezistenţelor mecanice. Însă, tot experimental s-a constatat că succesiunea ciclurilor de încărcare-descărcare-reîncărcare conduce la un fenomen contrar de diminuare a rigidităţii, la o diminuare a rezistenţelor de rupere, fenomen cunoscut sub denumirea de „degradare dinamică". Toate aceste constatări experimentale conduc la concluiza că în cazul solicitărilor ciclice, cu viteze mari de deformare, degradarea dinamică a caracteristicilor mecanice ale materialului depăşeşte sporul de rigiditate adus de stânjenirea manifestărilor reologice prin creşterea vitezei de deformare. Diferenţe între răspunsul lent şi cel dinamic sugerează faptul că modelul vâscoplastic are mai mult o valoare teoretică şi mai puţin una practică. Pentru aplicaţii este mai sigură utilizarea variantelor modelului vâscoelastic neliniar, pentru solicitări statice sau dinamice, cu parametrii de material determinaţi prin teste statice sau, respectiv, dinamice. 99 ECHIVALENŢA LINIARĂ ÎN DINAMICA MATERIALELOR NELINIARE 9.1. Introducere Există probleme practice pentru a căror rezolvare nu este justificată generarea şi rezolvarea unor mari sisteme neliniare. De asemenea, există în mecanica structurilor metode eficiente de calcul bazate pe ipoteza de corp elastic liniar care sunt corecte pentru structurile propriu-zise dar neadecvate pentru materialele din amplasamentul acestor structuri. Amintim, în acest sens, metodele de interacţiune teren-structură care înlocuiesc răspunsul pământurilor pe o direcţie cu răspunsul dat de un model liniar tip Kelvin-Voigt. Un alt domeniu în care liniarizarea este utilă îl constituie metodele de control destinate protecţiei antiseismice, metode care operează în timp real şi unde nu pot fi folosite rezolvări ce necesită timp îndelungat. Toate aceste motive fac pe deplin justificată încercarea de a înlocui ecuaţia neliniară a unui oscilator cu un grad de libertate: (9.1) cu o ecuaţie liniar echivalentă de forma: (9.2) deci, de a înlocui funcţiile neliniare de material şi cu constantele respective c şi k, modificare care, însă, nu trebuie să afecteze sensibil soluţia sistemului. Avantajele aduse de liniarizare pentru rezolvarea problemelor structurale sunt evidente. Practic, dificultăţile sunt transferate determinării constantelor echivalente astfel ca soluţia liniarizată să fie cât mai apropiată de soluţia liniară. Posibilităţile de depăşire a acestor dificultăţi fac obiectul acestui capitol. Ecuaţia de mişcare a oscilatorului neliniar (9.1) poate fi scrisă în forma adimensionalizată (vezi capitolul 2): (9.3) iar pentru sistemul liniarizat forma adimensionalizată este: (9.4) După ce efectul condiţiilor iniţiale se anulează, soluţia stabilizată a sistemului liniar echivalent (9.4) poate fi scrisă: (9.5) unde: (9.6) Soluţia (9.5) se mai poate scrie şi sub forma: (9.7) unde: (9.8) 9.2. Metode de liniarizare Vom aborda calculul constantelor liniar echivalente, c şi k, prin trei metode relativ cunoscute [124] urmând ca pe baza acestora să prezentăm o nouă metodă. O primă metodă constă în minimizarea diferenţelor dintre soluţiile sistemelor (neliniar şi liniar echivalent), a doua constă în minimizarea diferenţelor dintre „forţele” elastice şi a celor de de amortizare ale celor două sisteme iar a treia metodă este o metodă de liniarizare globală [39]. 9.2.1. Minimizarea diferenţei dintre soluţii Prima metodă constă în calculul diferenţei dintre soluţia ecuaţiei neliniare (9.3) şi soluţia ecuaţiei linar echivalente (9.4). Mediind această diferenţă pe o perioadă obţinem o funcţie care depinde de cele două constante pe care dorim să le determinăm : (9.9) funcţie care în mod evident trebuie să fie minimă pentru a avea cea mai bună echivalenţă între cele două sisteme. Determinarea punctelor de minim ale funcţiei (9.9) se face rezolvând sistemul: (9.10) de unde, prin introducerea soluţiei stabilizate (9.7) şi gruparea convenabilă a termenilor, rezultă sistemul de două ecuaţii în cele două necunoscute c, k: (9.11) Menţionăm că la ambele ecuaţii ale sistemului (9.11) expresiile din membrul drept depind în mod indirect de constantele c şi k prin intermediul factorului de amplificare definit în ec.(9.8)1. De asemenea, se observă că ambele ecuaţii ale sistemului (9.11) depind, cum era de aşteptat, de amplitudinea µ şi frecvenţa excitaţiei ω. Rezolvarea numerică a sistemului (9.11), pentru o amplitudine şi frecvenţă adimensională , furnizează valorile corespunzătoare c şi k, valori dependente de µ şi υ. 9.2.2. Minimizarea diferenţelor dintre forţele elastice şi de amortizare O altă metodă de determinare a constantelor ecuaţiei liniar echivalente (9.4) impune ca mediile pe un ciclu ale diferenţelor dintre forţele elastice şi de amortizare din cele două sisteme, să fie nule. În acest caz se ajunge la sistemul: (9.12) în care ϕ (τ) este soluţia stabilizată a sistemului liniar echivalent. Introducând în sistemul (9.12) expresia (9.5) a soluţiei stabilizate, rezultă: (9.13) unde: (9.14) Şi sistemul (9.14) care determină constantele c şi k este neliniar, deci pentru obţinerea valorilor acestor constante este necesară utilizarea metodelor numerice de rezolvare. 9.2.3. Liniarizarea globală Această metodă corespunde unui criteriu de liniarizare ce utilizează numai o estimare a domeniului de variaţie a soluţiei ϕ. Notăm cu [ domeniul maxim de variaţie al soluţiei şi cu [] domeniul maxim de variaţie pentru . Ţinând seama de simetria caracteristicilor elastice şi de amortizare, putem determina coeficienţii liniar echivalenţi din rescrierea condiţiilor de liniarizare statistică sub formă globală. Astfel, din prima metodă de linarizare statistică, rescrierea în formă globală transformă relaţia de definiţie (9.9) în: (9.15) iar din cea de a doua metodă de linarizare statistică rezultă condiţiile globale: (9.16) Aplicând condiţiile (9.10) funcţiei (9.15) obţinem sistemul de ecuaţii: (9.17) unde: (9.18) Observăm că de data aceasta, termenii din dreapta ai ecuaţiilor sistemului (9.17) nu mai depind de constantele c şi k, deci sistemul (9.18) este liniar. Prin rezolvarea acestui sistem se obţine soluţia: (9.19) Acestă soluţie poate fi validată numai dacă şi , adică: (9.20) Dacă introducem soluţia stabilizată a sistemului liniar echivalent (9.7) în sistemul (9.17) se obţin constantele sistemului liniar echivalent într-o nouă formă: (9.21) Forma (9.21) a constantelor liniar echivalente are avantajul unui calcul imediat, evitând rezolvarea numerică a unor sisteme neliniare. Însă, necesită estimarea prealabilă a domeniului de variaţie a soluţiei, valoarea constantelor c şi k fiind, de fapt, dependentă de alegerea acestui domeniu. Dacă există soluţia neliniară se poate determina limta de integrare din integralele ec.(9.21), iar constantele liniarizate c şi k se determină direct. 9.3. Liniarizarea globală în regim de rezonanţă După cum se poate constata din analiza răspunsului sistemului neliniar, factorul de amplificare (9.8) este dependent de amplitudinea excitaţiei. Deci, este necesară determinarea unui sistem liniar echivalent pentru fiecare nivel de excitaţie. Ţinând cont că determinarea caracteristicilor de amortizare şi rigiditate C(ϕ) şi K(ϕ) a fost făcută în regim de rezonanţă, pentru diferite excitaţii µ, valorile şi pot fi, de asemenea estimate în regim de rezonanţă. Valoarea la rezonanţă a factorului de amplificare poate fi determinată pe calea obişnuită prin anularea derivatei: (9.22) reţinând soluţia cu semnificaţie fizică : la care corespunde următoarea valoare a factorului de amplificare: (9.23) Prin urmare: (9.24) Deoarece limita de integrare depinde c şi k din relaţiile (9.21) rezultă un sistem care poate fi rezolvat numeric, rezultând constantele liniar ehivalente căutate. Astfel, de exemplu, având funcţiile neliniare C(ϕ) şi K(ϕ) determinate experimental în forma exponenţială: (9.25) din ec.(9.21) cu limita de integrare din (9.24) se obţine sistemul: (9.26) care poate fi rezolvat numeric, pentru diferite valori µ date, obţinând pentru fiecare amplitudine a excitaţiei µ perechea de constante liniar echivalente c şi k. Această metodă permite ca prin calcularea soluţiei liniarizate să se poată realiza un control al liniarizării prin comparaţia cu soluţia neliniară . Pentru a exemplifica procedeul redăm în tabelul 9.2 constantele liniar echivalente calculate cu sistemul (9.26) cu parametrii funcţiilor (9.25) din tabelul 9.1. Tabelul 9.1 a1 a2 a3 b1 b2 b3 0,308 0,269 119 0,25 0,75 112 Tabelul 9.2 µ c k 0,041 0,995 0,063 0,914 0,108 0,760 0,141 0,649 0,246 0,334 0,277 0,273 0,301 0,251 0,305 0,250 Cu aceste constante liniar echivalente au fost determinate curbele de amplificare date în fig.9.1, care pot fi comparate cu cele din fig.9.2 determinate din soluţiile neliniare corespunzătoare. Fig. 9.1 Comparaţia dintre curbele de amplificare corespunzătoare aceleiaşi amplitudini ale excitaţiei µ relevă faptul că metoda de liniarizare conduce la rezultate satifăcătoare. Numai la valori µ mari apar inerente, dar rezonabile, diferenţe. În final, remarcăm că, în realitate, constantele liniar echivalente determinate sunt de fapt funcţii dependente de amplitudinea normalizată µ. Aparent, prin liniarizare, funcţiile neliniare de material şi ale sistemului neliniar au fost înlocuite cu alte funcţii neliniare şi iar efortul de a liniariza sistemul neliniar nu pare a fi justificat. Însă, trebuie să remarcăm că în timp ce argumentul ϕ al funcţiilor neliniare de material şi este şi necunoscută a problemei, argumentul „constantelor” de liniarizare – amplitudinea normalizată µ este o mărime dată de la început: (9.27) Prin urmare, cunoscând amplitudinea excitaţiei şi caracteristicile mecanice ale sistemului şi , µ este cunoscut, iar după determianrea constantelor liniar echivalente c şi k, amplitudinea de răspuns a sistemului liniarizat se obţine imediat: (9.28) Fig. 9.2 Pentru facilitatea aplicaţiilor practice pot fi determinate funcţii de material dependente de caracteristici ale excitaţiei, de forma CC ; sau, revenind la caracteristicile dimensionale, ; utilizând în acest scop datele rezultate din procesul de liniarizare. Astfel, utilizând valorile constantelor liniarizate din tabelul 9.2 se obţin funcţiile (fig. 9.3): (9.29) şi ţinând cont de caracteristicile mecanice şi de încărcare ale sistemului oscilant folosit la liniarizare () rezultă: (9.30) funcţii dependente de amplitudinea excitaţiei care pot fi considerate caracteristicile neliniare de material ale modelului Kelvin-Voigt neliniar. Fig. 9.3 1100 RUPEREA DINAMICĂ LA PĂMÂNTURI ŞI ROCI 10.1. Introducere Acest capitol dedicat ruperii la pământuri este structurat pe trei secţiuni principale. Mai întâi, va fi abordată problema determinării suprafeţei de rupere pentru pământuri în condiţiile care fac acceptabilă ipoteza comportării de corp elastic neliniar. Apoi, vor fi expuse câteva particularităţi ale determinării stărilor limită la pământurile cu manifestări reologice neneglijabile, apoi vom analiza efectele proprietăţilor reologice asupra rupeii, urmând ca, în final, să abordăm ruperea în condiţii de solicitare dinamică, mai întâi sub încărcări armonice apoi sub încărcări seismice. Nu pot exista stări de tensiune în tot cuprinsul spaţiului invarianţilor tensorului tensiunii (vezi paragraful 1.2). Suprafaţa care delimitează în interiorul său stările de tensiune fizic posibile este numită suprafaţa stărilor limită sau suprafaţa de rupere, stările de tensiune existente pe această suprafaţă provocând scoaterea din uz a materialului. Formularea condiţiilor mecanice care conduc la epuizarea capacităţii de rezistenţă a materialului sunt numite criterii de rezistenţă, acestea fiind de fapt, expresii analitice ale suprafeţei stărilor limită. Mai înainte de a ne referi la criteriile de rezistenţă utilizate la pământuri considerăm necesară precizarea sensului atribuit noţiunilor de „curgere” şi „rupere”. La pământuri epuizarea capacităţii de rezistenţă mecanică se face aparentă printr-o multitudine de forme de manifestare, în funcţie de proprietăţile fizice şi chimice ale diverselor materiale reunite sub denumirea generică de „pământuri”. Apar, ca efect al solicitărilor exterioare, atât distrugeri evidente prin decoeziune cât şi procese de deformare continuă sub tensiuni constante, caz în care coeziunea nu este complet anihilată dar materialul îşi pierde capacitatea de a prelua noi încărcări. Mai mult, se întâlnesc frecvent cazuri când procesul de deformare continuă este brusc întrerupt prin apariţia decoeziunii. Aspectul probelor triaxiale solicitate până la epuizarea capacităţii reflectă formele de cedare amintite. Astfel, când proba îşi pierde brusc coeziunea (naturală sau aparentă dată de presiunea hidrostatică din celulă) apar suprafeţe de separare distincte, în general sub forma unor plane înclinate cu un unghi faţă de bazele probei. În alte cazuri, deformaţiile cresc continuu iar proba, sub formă de „butoi”, îşi modifică continuu dimensiunile până când pistonul aparatului îşi epuizează cursa, sau până când, brusc, apare un plan de separare. Deteriorarea structurii materialului prin decoeziune bruscă este numită rupere casantă (fracture) [9]. Specifică rocilor, decoeziunea este întâlnită şi la unele categorii de pământuri cu densităţi şi rezistenţe mecanice mai ridicate. Procesul limită de deformare continuă sub tensiuni constante este cunoscut sub denumirea de curgere (flow) şi este un proces obişnuit la pământuri. Ambele procese care conduc la scoaterea din uz a materialului, fie prin epuizarea capacităţii de a prelua încărcări, fie prin depăşirea deformaţiilor admisibile, sunt reunite sub denumirea curentă de rupere (failure). Prin suprafaţa starilor limită sau suprafaţa de rupere vom înţelege, în cele ce urmeză, acea suprafaţă din spaţiul invarianţilor tensorului tensiunii constituită din locul geometric al stărilor de tensiune care provoacă ruperea, indiferent de procesul prin care materialul este scos din uz. Depistarea experimentală a limitelor până la care se pot solicita pământurile se face, deobicei, concomitent cu încercările necesare determinării parametrilor ecuaţiilor constitutive, când se utilizează drumuri în tensiuni diferite care pleacă sau nu din acelaşi punct de stare iniţial. O experienţă triaxială va furniza un punct al suprafeţei de rupere, şi anume, punctul în care drumul în tensiuni va întâlni suprafaţa de rupere. Valoarea invariantului τ din acest punct este numită rezistenţă de rupere. Prin urmare, pentru determinarea întregii suprafeţe de rupere sau numai a unor porţiuni sau curbe caracteristice ale ei sunt necesare, în principiu, mai multe încercări de rupere pe probe repetabile solicitate pe drumuri diferite. În cele ce urmeză vom raporta suprafaţa de rupere la "coordonatele" formate din invarianţii σ, τ şi c, forma generală a suprafeţei de rupere fiind, deci: (10.1) Urma suprafeţei de rupere în planul intensităţilor , planul în care , o vom numi secţiune meridiană, iar urma suprafeţei de rupere dintr-un plan deviator o vom numi secţiune transversală [30]. 10.2. Ruperea în condiţii statice 10.2.1. Criterii de rupere În decursul timpului, în funcţie de evoluţia cunoştinţelor despre comportarea pământurilor şi de aparatura disponibilă, au fost propuse sau utilizate o serie de criterii pentru ruperea pământurilor, multe din ele fiind „împrumutate”, cu sau fără modificări, de la alte materiale cu comportări mai simple şi mai cunoscute. Vom expune pe scurt, în continuare, câteva din criteriile de rupere folosite la pământuri, insistând îndeosebi asupra fondului experimental adus şi asupra cadrului ipotetic bazat pe acest fond experimental. Primele criterii folosite au fost cele bi-invariante bazate pe datele experimentale furnizate de clasicul „triaxial”, în fapt un aparat de compresiune biaxială. Au urmat, apoi, criterii tri-invariante în forma completă (10.1) propuse de cei care au putut dispune de un aparat capabil să impună drumuri cu . Fig. 10.1 Criteriul Mohr-Coulomb [9], [120] în forma sa iniţială admitea că ruperea apare atunci când tensiunile tangenţiale dintr-un plan ating valorile maxim posibile, valori ce depind de tensiunile normale pe planul de rupere. Extins sub formă invariantă , acest criteriu defineşte o curbă meridiană a cărei aproximaţie liniară este: (10.2) unde c este coeziunea iar φ unghiul de frecare internă, ambele fiind considerate constante de material (fig.10.1). Deci, secţiunea meridiană Mohr-Coulomb este liniară, iar secţiunile transversale prin piramida hexagonală neregulată definită de ec.(10.2) au forme de hexagoane neregulate cu laturi egale dar neparalele (fig.10.2). S-a observat că la presiuni hidrostatice mari rezistenţa la forfecare nu mai depinde de invariantul sferic σ, pământurile admiţând criteriul Huber-Misses-Hencky [9]: (10.3) Fig. 10.2 Evident, trecerea de la comportarea la presiuni hidrostatice mici la indiferenţa faţă de σ nu poate fi bruscă. Din acest motiv, s-au propus racordări continue între cele două suprafeţe (10.2) şi (10.3). Astfel, Nelson şi Baron [108] au utilizat un criteriu de rupere cu o secţiune meridiană de forma: (10.4) unde este o valoare-prag a invariantului σ ce corespunde nivelului presiunii hidrostatice care nu mai influenţează rezistenţa de rupere. Deşi adecvate fizic, ecuaţiile (10.4) sunt incomod de utilizat, motiv pentru care au fost propuse curbe de racordare cu ecuaţii unice. Astfel, DiMaggio şi Sandler [65] au folosit exponenţiala: (10.5) unde a, b şi c sunt constante cu semnificaţii ce reies din fig.10.1 şi care se determină experimental. Se observă că, pentru , , iar pentru , , deci exponenţiala (10.5) are pentru punct comun cu dreapta Mohr-Coulomb şi tinde asimptotic la dreapta Misses . Pentru determinarea experimentală a formei secţiunilor transversale sunt necesare încercări pe drumuri cu . Astfel de încercări au fost făcute de Haythornthwaite [78], Kirkpatrick [85], Wu, Loh şi Malvern [137] şi alţii, folosind în acest scop triaxiale astfel modificate încât să poată exercita, pe lângă compresiuni axial-simetrice, şi torsiuni. Ulterior, a fost pus la punct un aparat triaxial veritabil, aparat care poate solicita o probă cubică pe trei direcţii. Astfel de încercări au fost comunicate de Gudehuş [70], Goldscheider [69], Hambly [73], Ko şi Scott [86], Lade [89], [90], [91], Pearce [112] şi alţii. Rezistenţele maxime obţinute în aceste aparate s-au situat, în secţiuni transversale, între laturile Mohr-Coulomb şi cercul Misses (fig.10.2). Pornind de la aceste observaţii experimentale au fost propuse drept criterii de rupere ecuaţii care duc la secţiuni transversale curbe. Coleman [57] a propus un criteriu tri-invariant în forma: din care prin înlocuirea invarianţilor cu rezultă: (10.6) ecuaţie din care se obţine, pentru şi o curbă plasată între latura Mohr-Coulomb şi cercul Misses (fig.10.3). Criteriul sugerat de Gudehuş [71] conţine invarianţii direcţionali şi în forma: unde A şi B sunt constante care se determină prin încercări. Înlocuind invariantul se obţine: (10.7) deci, pentru , rezultă de asemenea o curbă . Fig. 10.3 Lade [91] a utilizat drept criteriu de rupere o relaţie simplă: unde k este o constantă determinată experimental. Trecând în invarianţii , criteriul Lade ia forma ceva mai complicată: (10.8) din care rezultă, pentru , o curbă care, deşi diferită ca expresie analitică de curbele Coleman sau Gudehuş, este asemănătoare grafic cu acestea. În teza sa de doctorat, Stutz [127] a propus drept criteriu de rupere o relaţie de forma: (10.9) unde m este o constantă experimentală subunitară. Curbele se depărtează de laturile hexagonului Mohr-Coulomb şi se apropie de cercul Misses pe măsură ce valorile parametrului m cresc. Criteriile Coleman, Gudehuş, Lade şi Stutz dau într-un plan deviator curbe asemănătoare între ele, cu puncte comune cu hexagonul Mohr-Coulomb în şi (compresiune şi extensie axial-simetrică). Însă, toate au secţiuni meridiane liniare, probabil datorită faptului că au fost obţinute prin încercări cu presiuni hidrostatice reduse, caz în care curba este apropiată de dreapta Mohr-Coulomb. 10.2.2. Un nou criteriu tri-invariant Din sumara examinare efectuată în subparagraful anterior asupra unor criterii de rupere utlizate la pământuri se pot reţine câteva constatări experimentale care permit adoptarea următoarelor ipoteze: • Secţiunea meridiană a suprafeţei de rupere este o curbă de racordare între dreapta Mohr-Coulomb şi dreapta Misses. • Secţiunea transversală este o curbă cu puncte comune cu hexagonul Mohr-Coulomb în şi . • Toate secţiunile transversale au aceiaşi formă, influenţa invariantului σ traducându-se numai prin modificarea mărimii secţiunii. • Admiţând că materialul are o comportare izotropă, secţiunea transversală are şase sectoare simetrice cu unghi la centru de 60. În consens cu aceste ipoteze, vom încerca obţinerea unui criteriu de rupere a cărei ecuaţie să fie formată prin combinarea ecuaţiei curbei meridiane cu ecuaţia curbelor transversale , criteriu care va trebui să satisfacă şi următoarele condiţii : • Forma suprafeţei de rupere să fie în concordanţă cu datele furnizate de aparatura adecvată investigării planului deviatoric ; • Constantele ce intervin în ecuaţia suprafeţei de rupere să poată fi determinate prin încercări în triaxialul obişnuit existent în ţară. Drept ecuaţie a curbei meridiane vom adopta exponenţiala DiMaggio-Sandler ec.(10.5) care satisface prima ipoteză iar constantele de material conţinute pot fi determinate prin încercări triaxiale obişnuite. Pentru a obţine forma secţiunii transversale, adică a dependenţei în punctele de rupere dintre invarianţii τ şi , vom apela la relaţia dintre aceşti invarianţi de până la rupere, adică la ecuaţia constitutivă elastic neliniară (4.19)2, o relaţie de forma : (10.10) unde: (10.11) (10.12) Scriind această ecuaţie pentru drumuri situate într-un plan deviator, drumuri cu , după o rearanjare se obţine: (10.13) unde prin au fost notate combinaţiile dintre constantele de material şi invariantul γ al tensorului deformaţiei. Deci, efectuând experienţe cu drumuri , se pot obţine valori ale combinaţiilor , iar relaţia (10.13) dintre invarianţii τ şi poate fi complet determinată. În lipsa aparaturii capabile să impună drumuri , ec.(10.13) poate fi înlocuită cu o alta derivată din ea, făcând apel la ipotezele admise şi condiţiile impuse. Se observă că, din ec.(10.13), pentru anumite valori ale invariantului direcţional c se regăsesc solicitările particulare: : compresiune triaxială : : forfecare pură : : extensie triaxială : de unde, rezultă: (10.14) Ca o consecinţă a ipotezei de coincidenţă cu hexagonul Mohr-Coulomb în punctele şi , raportul poate fi exprimat în funcţie de unghiul de frecare internă φ: (10.15) astfel încât ec.(10.14) devine: (10.16) În aceste condiţii, ec. (10.13) poate fi înlocuită cu: (10.17) relaţie ce poate fi complet determinată prin încercări triaxiale, iar funcţia de formă dă o secţiune meridiană curbă situată între curbele Gudehuş şi Lade (fig.10.3). Conform ipotezei care admite conservarea secţiunilor transversale , efectul invariantului σ este inclus în ecuaţia curbei meridiane . Deci, prin combinarea ecuaţiilor (10.5) şi (10.17) rezultă următorul criteriu tri-invariant: (10.18) Parametrii acestui criteriu : a, b, c, φ, au semnificatiile care reies din fig.10.1 şi pot fi obţinuţi prin încercări experimentale în triaxialul obişnuit. Atunci când încărcările exterioare provoacă tensiuni sferice reduse, curbura secţiunii meridiane este redusă, caz în care aproximaţia liniară a secţiunii meridiane dată de criteriului Mohr-Coulomb devine acceptabilă (tabelul 10.1) [17], ecuaţia suprafeţei de rupere luând forma: (10.19) Tabelul 10.1 Material σ τ Măsurat Liniar* Exponenţial** - [Mpa] [Mpa] [Mpa] [Mpa] Argilă 0 - 0,0879 0,0852 0,2245 0,1005 0,1040 0,1055 0,6149 0,1342 0,1288 0,1325 0,9792 0,1600 0,1531 0,1548 1,9325 0,1698 0,1766 0,1731 Nisip 0,4367 0,3756 0,3183 0,3621 1,1865 0,9431 0,8677 0,8217 1,8277 1,3595 1,3365 1,3444 2,4163 1,7024 1,7670 1,6922 Argilă: *) **) Nisip: *) **) 10.3. Efectul proprietăţilor reologice asupra ruperii Pământurile, cele argiloase în special, sunt cunoscute ca materiale cu proprietăţi reologice apreciabile ceea ce face ca procesele de deformare să continue şi după suprimarea cauzei exterioare care le-a provocat. Proprietăţile reologice îşi semnalează prezenţa şi în comportarea probei triaxiale atât în cazul încărcărilor în trepte cât şi atunci când proba este solicitată cu viteză de încărcare constantă şi continuă până la rupere. Experimental, s-a constat că modificarea vitezei de încărcare conduce la modificări în aspectul curbelor tensiune-deformaţie dar şi la obţinerea unor valori diferite ale rezistenţei de rupere [17], [31], [40]. Pe măsură ce viteza de încărcare creşte, procesele de deformare sunt stânjenite tot mai mult, iar rezistenţele de rupere obţinute au valori din ce în ce mai mari. Deci, acest tip de încercare conduce la o supraestimare a rezistenţei de rupere. În cazul încărcărilor aplicate în trepte, efectele reologice sunt şi mai evidente. Fiecare treaptă de încărcare menţinută un interval de timp dat constituie, de fapt, o încercare de fluaj în care deformaţiile evoluează până la stabilizare sau rupere. Este evident faptul că rezistenţele de rupere obţinue în acest tip de încercare vor avea valori mai reduse. În capitolul 4, în paragraful dedicat determinării experimentale a parametrilor ecuaţiei constitutive vâscoelastic neliniare, a fost expusă o metodă de determinare a constantelor din funcţiile neliniare de relaxare prin încărcări în trepte, fiecare treaptă fiind de fapt o încercare de fluaj. Acestă metodă poate folosi nu numai la cuantificarea comportării până la rupere ci oferă şi posibilitatea unei evaluări mai realiste a rezistenţei de rupere. Reamintim că după metoda amintită, din perechile de valori , ale celor n trepte de încărcare, colectate la aceleaşi trepte de timp , se obţin m izocrone în forma ecuaţiei constitutive elastic neliniare . Fig. 10.4 Se constată că izocronele diferă între ele în funcţie de timpul folosit şi conduce la estimarea unor valori de rupere diferite. Cu cât timpul este mai apropiat de momentul iniţial de aplicare a treptei de încărcare, cu atât valoarea rezistenţei de rupere este mai mare iar izocrona calculată cu valorile stabilizate ale deformaţiilor va da o valoare minimă a rezistenţei de rupere (fig.10.4). Aceste izocrone sunt asemănătoare curbelor tensiune-deformaţie obţinute din încercări cu viteze de încărcare diferite, dar constante, curbele obţinute cu viteze mari se apropie de izocrona iniţială, iar pe măsură ce viteza de încărcare scade, curbele se apropie de izocrona stabilizată. În concluzie, pentru a evita supraevaluarea rezistenţei de rupere datorată proprietăţilor reologice ale materialului, încercările experimentale trebuie efectuate fie cu o viteză de încărcare redusă, fie utilizând încărcările în trepte şi adoptând drept ultimă rezistenţă valoarea maximă a izocronei stabilizate. Pentru a ilustra diferenţele apreciabile de evaluare a rezistenţei de rupere în fig.10.4 sunt reunite câteva curbe tensiune-deformaţie obţinute cu viteze diferite, începând cu cele mai rapide, aşa cum sunt în testele din coloana rezonantă, continuând cu cele obţinute prin încercări triaxiale cu viteză de încărcare constantă şi terminând cu curba obţinută cu valorile stabilizate ale deformaţiilor dintr-o încercare triaxială în trepte. Toate aceste încercări au fost efectuate pe probe formate în laborator din aceeaşi argilă . 10.4. Determinarea ruperii prin oboseală la roci Deseori rocile sunt supuse unor solicitări dinamice datorate seismelor, exploziilor, traficului rutier şi feroviar, vibraţiilor fundaţiilor de maşini etc. În timp, astfel de solicitări pot provoca degradări mecanice ale materialului, degradări care se manifestă prin reducerea rezistenţelor mecanice, reducere care poate avansa până la rupere prin oboseală. Obosela materialelor a fost intens studiată în special la metalele din componenţa organelor de maşini. Există numeroase studii efectuate pe metale atât în scopul elucidării mecanismului de rupere prin oboselă cât şi în scopul cuantificării rezistenţelor de rupere dinamică [9]. Uzual, la materialele metalice, reducerea rezistenţelor mecanice datorită oboselii se determină experimental prin solicitarea dinamică a unei epruvete prelevate din materialul respectiv cu un număr de cicluri cu o amplitudine dată iar apoi, epruveta este solicitată până la rupere. Prin comparare cu rezistenţa de rupere determinată pe o altă epruvetă nesolicitată dinamic se poate evalua diminuarea rezistenţei de rupere datorată solicitării dinamice. Repetând procedeul cu o serie de epruvete supuse la numere diferite de cicluri se obţine o curbă , cunoscută sub denumirea de curba lui Wöhler, curbă care reflectă diminuarea rezistenţei de rupere în funcţie de amplitudinea şi numărul ciclurilor. Extinderea acestei tehnologii a fost încercată şi în cazul unor roci naturale (marmore, gresii) sau artificiale (ipsos, beton), încercări care au relevat diminuări importante ale rezistenţelor de rupere [63]. Însă, datorită neomogenităţilor pronunţate ale majorităţii geomaterialelor, existenţa unor „probe repetabile” este puţin probabilă, ceea ce introduce erori inacceptabile datorate comparării rezistenţei de rupere a unei probe solicitată dinamic cu rezistenţa de rupere statică determinată pe o altă probă. În acest paragraf prezentăm un procedeu de depăşire a acestor dificultăţi prin utilizarea unor rezistenţe de rupere determinate nedistructiv cu ajutorul corelaţiei cu densitatea materialului ceea ce permite o determinare mult mai precisă a curbei Wöhler [16]. 10.4.1. Corelaţia rezistenţă de rupere-densitate Pentru a evalua reducerea rezistenţei de rupere datorată exercitării unui anumit număr de cicluri fără a apela la comparaţia cu rezistenţa de rupere determinată pe o altă probă este necesară estimarea rezistenţei de rupere pe care proba respectivă a avut-o înaintea aplicării solicitării dinamice. În acest scop poate fi folosită relaţia de proporţionalitate care există între rezistenţa de rupere a unei probe şi densitarea sa ρ. Prin rezistenţă de rupere vom înţelege valoarea intensităţii tensiunilor tangenţiale τ care provoacă ruperea, valoare determinată în triaxial pe drumul în tensiuni corespunzător solicitării dinamice date. Deci, în spaţiul tensiunilor principale (, rezistenţa de rupere corespunde punctului în care drumul în tensiuni intersectează suprafaţa de rupere. În aparatul de compresiune axial-simetrică cunoscut sub denumirea de triaxial, unei probe cilindrice (uzual cu înălţimea şi diametrul ) i se impune o stare de tensiune axial-simetrică cu astfel că invariantul τ are expresia: (10.20) unde J2 este al doilea invariant al deviatorului tensorului tensiune: (10.21) care în cazul stării de tensiune axial-simetrică devine: (10.22) iar intensitatea tensiunilor tangenţiale se reduce la: (10.23) Drumurile în tensiuni care pot fi impuse în triaxial sunt situate în planul al spaţiului tensiunilor principale (fig.10.5) şi pot fi astfel alese încât să fie cât mai apropiate solicitării exterioare analizate. Fig. 10.5 Astfel, în cazul în care solicitarea dinamică este compusă din vibraţii verticale suprapuse peste o stare de tensiuni iniţiale – cazul solicitărilor date de fundaţiile de maşini sau de trafic – drumul în tensiuni indicat este drumul cu presiune laterală constantă (1 din fig.10.6), iar în cazul în care solicitarea exterioară constă în undele de forfecare date de seism sau explozii, mai indicat este ca în triaxial să se impună drumul cu presiune medie constantă (2 din fig.10.6). Ambele drumuri au ca punct de plecare punctul din spaţiul invarianţilor corespunzător stării iniţiale „in situ” (A din fig.10.6). Fig. 10.6 10.4.2. Curba tip Wöhler Solicitând până la rupere n probe cu densităţi diferite, prelevate din stratul şi amplasamentul analizat se poate evalua corelaţia dintre rezistenţa de rupere şi densitate. Ulterior, pentru determinarea unei curbe de tip Wöhler rezistenţa de rupere a unei probe, rezistenţă diminuată de solicitarea cu un anumit număr de cicluri, se poate compara cu rezistenţa de rupere evaluată înaintea solicitării în funcţie de densitatea iniţială. Vom exemplifica acest procedeu utilizând rezultatele experimentale obţinute pe probe de calcar [16]. Deoarece amplasamentul analizat era solicitat de vibraţiile generate de un turboagregat, atât încercările de rupere cât şi solicitările ciclice au fost efectuate pe drumul cu presiune laterală constantă, punctul de plecare fiind punctul corespunzător tensiunilor date de sarcina statică existentă „in situ” (). Prin ruperea a 7 probe de densităţi diferite s-a obţinut, mai întâi, corelaţia dintre rezistenţa de rupere statică, probabilă (fig.10.7): (10.24) Apoi, au fost solicitate cinci probe cu cicluri pulsatorii de amplitudini între 0,6 – 0,9 din rezistenţa de rupere statică probabilă a fiecărei probe determinând numărul de cicluri care au provocat ruperea. S-a obţinut astfel curba Wöhler normalizată în raport cu rezistenţa de rupere statică (fig.10.8): (10.25) Fig. 10.7 [kg s2 / cm4] Fig. 10.8 În cazul investigat, amplitudinea solicitării date raportată la rezistenţa de rupere a fost sub valoarea asimptotei dmin, ceea ce indică improbabilă apariţia ruperii prin oboseală. 10.5. Efectul de degradare dinamică În cazul încărcărilor monoton crescătoare epuizarea capacităţii de rezistenţă mecanică a materialului – ruperea în condiţii „statice” – apare atunci când tensiunile tangenţiale ating o valoare limită sau când deformaţiile scot materialul din exploatare. Atunci când solicitarea este repetată ciclic, chiar cu amplitudini mai mici decât amplitudinea care provoacă ruperea statică, capacitatea de rezistenţă se reduce cu fiecare ciclu, reducere ce poate avansa până la rupere. Valorile de rezistenţă dinamică sunt substanţial mai mici decât rezistenţele statice de rupere, solicitarea dinamică provocând o diminuare a rigidităţii materialului, o degradare a proprietăţilor de rezistenţă. În practica curentă rezistenţele de rupere ale diverselor materiale care pot fi considerate elastic liniare se determină prin încercări statice de laborator, completate în cazul solicitărilor dinamice cu încercări de oboseală. La pământuri şi roci, deşi fenomenul de degradare dinamică este semnalat în literatura de specialitate [55], [79], datorită proprietăţilor neliniare şi reologice ale acestor materiale, încercări de oboseală cu o metodă „clasică” preluată de la metale sunt deosebit de dificil de executat şi interpretat ceea ce conduce la evitarea lor şi la folosirea, în cel mai bun caz, a unor corecţii empirice. Dacă la roci dificultăţile date de neomogenitatea materialului au putut fi depăşite cu ajutorul relaţiei dintre rezistenţa de rupere şi densitate, la pământuri, materiale cu o comportare mecanică mult mai complexă, o astfel de metodă este ineficientă, fiind necesară o altă cale de abordare pentru determinarea rezistenţei de rupere şi a rezistenţei la oboseală. O astfel de metodă se bazează pe comportarea mecanică a materialului pe tot parcursul procesului de deformare, tratând ruperea ca un ultim stadiu al procesului de degradare a caracteristicilor mecanice, proces provocat de solicitarea dinamică [31]. În acest mod, ruperea este evaluată analizând comportarea dinamică pe tot parcursul solicitarii prin extindere şi nu abordând direct la ultimele resurse de rezistenţă. Metoda permite o aproximare mult ameliorată a ruperii dinamice şi se bazează pe datele experimentale obţinute în aparatele de laborator care imprimă probei solicitări torsionale (triaxialul dinamic, aparatul de torsiune dinamică). De asemenea, pot fi folosite şi încercările pulsatorii în aparatul triaxial obişnuit [25]. Răspunsul pământurilor la solicitări dinamice torsionale poate fi caracterizat de funcţia-modul de torsiune dinamică G. Prin definiţie, o valoare a acestei funcţii este dată de raportul dintre amplitudinile invarianţilor τ şi γ de la un moment dat t: (10.26) Această definiţie este în concordanţă cu sensul atribuit funcţiilor neliniare de relaxare a modelului vâscoelastic neliniar (vezi, cap.6) şi cu forma „secantă” a funcţiilor-modul elastic neliniare utilizată în cazul solicitărilor statice (vezi, cap.4). Deoarece amplitudinile invarianţilor τ şi γ definesc într-un plan (τ,γ) punctele de extrem ale curbei histeretice , funcţia G este legată de locul geometric al acestor extreme care, în planul (τ,γ) formeză o curbă numită curbă-schelet (backbone curve) [25], [45], [47], [79], [80] (fig.10.9). Fig. 10.9 La fel ca funcţia-modul de torsiune dinamică G, şi curba-schelet prezintă o formă neliniară, depinzând atât de nivelul de deformaţie γ cât şi de frecvenţa ω. Însă, s-a constatat experimental [17], [29], [67], [74], [75], [76], [77], [80] că pentru domeniul de frecvenţe ce prezintă interes în ingineria seismică , caracteristicile mecanice sunt puternic influenţate de numărul total al ciclurilor de încărcare n, numărul acestora pe unitatea de timp ω având o pondere redusă în cazul ciclurilor cu amplitudini mari. De asemenea, încercările pulsatorii triaxiale efectuate pe probe repetabile de argilă cu viteze de încărcare diferite dar cu aceeaşi amplitudine γ nu au condus la valori G diferite pentru acelaşi număr de cicluri [25]. Din acest motiv, vom utiliza în acest caz – cazul solicitărilor dinamice cu amplitudini mari - drept argumente ale funcţiei-modul de torsiune dinamică, amplitudinea γ şi numărul ciclurilor n, funcţie care va lua forma : . Influenţa numărului de cicluri asupra funcţiei-modul de torsiune dinamică poate fi abordată cu ajutorul oscilaţiilor staţionare în triaxial. Astfel, într-o experienţă triaxială ciclică cu tensiune controlată, menţinând amplitudinea τ constantă se observă o creştere a amplitudinilor deformaţiei γ (fig.10.10), iar în experienţele cu deformaţie controlată, menţinând constantă amplitudinea γ se constată o scădere a amplitudinilor tensiunii τ (fig.10.11). Fig. 10.10 Fig. 10.11 În ambele cazuri, valorile funcţiei-modul de torsiune dinamică scad iar influenţa numărului de cicluri se traduce printr-o diminuare a rigidităţii probei, printr-o degradare a proprietăţilor sale de rezistenţă mecanică. O măsură a degradării, de la primul ciclu până la ciclul n, este raportul: (10.27) dintre valoarea funcţiei-modul de torsiune dinamică calculată pentru ciclul n şi valoarea sa iniţială nedegradată. Cum valoarea acestui raport depinde atât de numărul ciclurilor cât şi amplitudinea lor şi reflectă degradarea mecanică a materialului, vom defini o funcţie de degradare: (10.28) unde G0 este valoarea iniţială, nedegradată a funcţiei-modul de torsiune dinamică. Determinarea acestei funcţii de degradare poate fi făcută experimental prin încercări triaxiale ciclice efectuate pe probe repetabile solicitate la oscilaţii staţionare dar cu amplitudini diferite. Astfel, într-un test cu „deformaţie controlată” o probă poate fi supusă unei oscilaţii pulsatorii cu o anumită amplitudine iar la fiecare ciclu n poate fi înregistrată forţa necesară menţinerii deformaţiei impuse, de unde se poate calcula valoarea degradării corespunzătoare numărului de cicluri n: (10.29) Prin prelucrarea statistică a perechilor de valori (),cicluri, pot fi găsite funcţii care să aproximeze optim distribuţia datelor experimentale. Repetând aceste determinări pe probe din acelaşi material, solicitate cu amplitudini γ diferite, se pot obţine un număr de funcţii , de unde se poate sintetiza o funcţie unică de degradare . În încercări cu „tensiune controlată” procedeul aste asemănător, rezultând funcţii de degradare de forma . Deoarece funcţia de degradare poate avea aparent forme diferite: (10.30) în încercări cu deformaţie controlată sau: (10.31) în încercări cu tensiune controlată, vom încerca obţinerea unei forme unice, utilizând valorile normalizate ale invarianţilor τ şi γ în raport cu valorile lor de rupere determinate în condiţii „statice”: sau (10.32) şi ţinând cont că: (10.33) rezultă că la rupere când : (10.34) deci, funcţia de degradare poate lua o formă unică , indiferent de tipul încercării folosite, iar argumentele ei devin amplitudinea normalizată r şi numărul de cicluri n. În continuare vom exemplifica metoda cu ajutorul datelor experimentale obţinute prin teste triaxiale pe probe de argilă [25], [31]. Astfel, prin încercarea a 8 probe repetabile de argilă cu 8 nivele diferite de deformaţii impuse s-au obţinut datele din tabelul 10.1. Prelucrarea statistică a datelor provenite din fiecare probă furnizează cele 8 funcţii de degradare (fig.10.12). Datele experimentale din tabelul 10.1 au permis o aproximare optimă a funcţiilor de degradare pentru fiecare palier de solicitare cu expresii analitice cu o formă comună: (10.35) dar cu constantele a şi b cu valori diferite, valori date în tabelul 10.2. Deoarece prin modificarea amplitudinii γ forma funcţiilor (10.35) se conservă, „parametrii” a şi b devin funcţii dependente de r (fig.10.13): (10.36) astfel încât funcţia de degradare devine: (10.37) funcţie cu aspectul din digrama spaţială 10.14. Tabelul 10.1 Test nr. Solicitare Valori măsurate γ r n τ G d = G/G0 (G0 = 109,5 MPa) [%] - - [kPa] [Mpa] 1 2 3 4 5 6 7 1 0,10 0,059 1 101,0 101,0 0,922 2 99,3 99,3 0,907 5 97,6 97,6 0,891 10 96,4 96,4 0,880 20 95,5 95,5 0,872 50 94,9 94,9 0,867 100 94,7 94,7 0,865 2 0,25 0,147 1 224,8 89,9 0,821 2 210,0 84,0 0,767 5 193,5 77,4 0,707 10 182,6 73,0 0,667 20 174,4 69,8 0,638 50 167,7 67,1 0,613 100 166,0 66,4 0,606 3 0,50 0,294 1 365,0 73,0 0,667 2 324,5 64,9 0,593 5 277,0 55,4 0,506 10 254,0 50,8 0,464 20 243,5 48,7 0,445 50 234,5 46,9 0,428 100 231,0 46,2 0,422 4 0,75 0,441 1 450,8 60,1 0,549 2 375,0 50,0 0,457 5 300,8 40,1 0,366 10 271,5 36,2 0,331 20 256,5 34,2 0,312 25 254,3 33,9 0,310 31 246,8 32,9 0,300 5 0,80 0,471 1 458,1 57,3 0,543 2 386,3 48,3 0,441 5 306,6 38,3 0,350 10 275,9 34,5 0,315 20 269,8 33,7 0,308 24 255,4 33,2 0,303 Tabelul 10.1 – continuare 1 2 3 4 5 6 7 6 0,90 0,529 1 479,9 53,3 0,487 2 358,3 42,8 0,391 5 316,3 35,1 0,321 7 307,5 34,2 0,312 10 292,7 32,5 0,297 7 1,00 0,588 1 514,4 51,1 0,467 2 403,0 41,7 0,381 4 361,4 36,1 0,330 6 325,2 34,5 0,315 8 1,25 0,735 1 522,9 41,8 0,382 2 436,6 34,9 0,319 3 406,3 32,5 0,297 Fig. 10.12 Tabelul 10.2 Nr. crt. γ [%] r a b n (la rupere) 1 0,10 0,059 0,581 0,963 - 2 0,25 0,147 0,794 0,979 - 3 0,50 0,294 1,249 1,008 - 4 0,75 0,441 1,743 0,981 38,7 5 0,80 0,471 1,761 0,928 27,4 6 0,90 0,529 1,887 0,872 9,5 7 1,00 0,588 2,000 0,814 5,5 8 1,25 0,731 2,090 0,633 2,7 1,000 1,0 Fig.10.13 Fig.10.14 Ruperea în condiţii dinamice poate fi acum definită de valorile minime ale funcţiei de degradare [31]: (10.38) valori, ce descriu în spaţiul o curbă obţinută prin intersecţia suprafeţei cu planul d(fig.10.14). În formulare explicită curba este o curbă tip Wöhler normalizată în raport cu rezistenţa de rupere statică [10]. În cazul folosit pentru exemplificare această curbă a fost determinată în forma (fig.10.15): (10.39) Fig.10.15 10.6. Ruperea în condiţii seismice Solicitările provocate de cutremure sunt solicitări aleatorii în amplitudine şi frecvenţă ce nu pot fi simulate în aparatura de laborator capabilă să impună probei numai vibraţii sau oscilaţii armonice. Din acest motiv, modelarea comportării dinamice a pământurilor se bazează pe determinarea în laborator a constantelor de material prin încercări cu frecvenţă şi amplitudine constantă. Desigur, ideal ar fi ca proba din aparatul de laborator să fie solicitată cu o forţă perturbatoare care să simuleze diagrama seismică, înregistrată sau sintetizată pentru amplasamentul analizat. Cum această cale nu este accesibilă, nu ne rămâne decât posibilitatea de a diminua discordanţa prin transformarea solicitării seismice într-o vibraţie armonică echivalentă cu care apoi, se solicită proba. În principiu, pentru a echivala o solicitare seismică este necesar, în primul rând, să se cunoască istoria tensiunilor şi/sau a deformaţiilor ce pot fi provocate în amplsament de un viitor cutremur. Datorită comportării neliniare a majorităţii materialelor din amplasament calculul starilor de tensiune şi de deformaţie provocate de cutremure rămâne o problemă deschisă. Însă, pentru procedeele de echivalare a mişcării seismice sunt acceptabile şi metode mai simplificate cum sunt cele bazate pe rezolvarea neliniarităţii cu metoda „elastic echivalentă” [58], [117]. Metodele de determinare a răspunsului în tensiuni a masivului de pământ sub acţiunea solicitărilor seismice sunt bazate pe ipoteza conform căreia efectele sunt datorate propagării verticale a undelor de forfecare care transmit o bună parte din energia seismică (~ 76%). Dintre aceste metode, cele care ţin cont de comportarea neliniară şi disipativă a pământurilor au dat rezultate în concordanţă cu observaţiile din teren [58], [122]. Însă, datorită dificultăţilor introduse de forma neliniară a ecuaţiei constiutive a pământurilor, aceste metode rezolvă neliniaritatea prin introducerea conceptului de comportare liniar echivalentă înlocuind materialul real cu un material vâscoelasic liniar echivalent care conduce la obţinerea aceloraşi tensiuni şi deformaţii. Caracteristicile mecanice ale acestui material echivalent sunt determinate printr-un procedeu iterativ, care va fi expus mai departe. În principiu, aceste metode implică următoarele etape [58] : • Detererminarea caracteristicilor mişcării seismice la nivelul rocii de bază : acceleraţia maximă, perioada predominantă, durata efectivă etc., utizându-se în acest scop accelerogramele unor cutremure puternice înregistrate anterior sau accelerograme generate sintetic. • Determinarea caracteristicilor dinamice ale fiecărui strat cuprins între roca de bază şi suprafaţa liberă a terenului, concretizate în determinarea completă a funcţiilor de material şi . • Calculul răspunsului seismic al masivului, concretizat şi în determinarea evoluţiei în timp a tensiunilor în punctele de interes, adică a tensogramei seismice. Au fost propuse mai multe metode de aproximare a vibraţiilor aleatorii prin vibraţii armonice echivalente [80], [113]. În esenţă, aceste metode constau în determinarea unei amplitudini armonice echivalente ( sau ) care, după un număr echivalent de cicluri , să provoace aceleaşi efecte ca şi solicitarea seismică. Un astfel de procedeu trebuie aplicat şi în cazul în care se urmăreşte determinarea rezistenţei de rupere dinamică. Vom prezenta în continuare o metodă de echivalare, utilizată iniţial de Seed şi Idris în calculul lichefierii [113], [123] şi adaptată pentru a putea fi folosită la determinarea degradării materialelor din amplasament, degradare provocată de solicitările seismice [31]. Această metodă presupune cunoscută diagrama seismică şi constă în următoarele etape : • Se determină din diagrama seismică tensiunile tangenţiale maxime provocate de seism; • Se determină o amplitudine echivalentă a tensiunilor tangenţiale armonice: (10.40) unde k este un factor empiric de conversie, cu valoarea 0,65 ; 0,75 sau 0,85 ; • Cunoscând rezistenţa de rupere statică , se calculeză o amplitudine normalizată echivalentă şi din ecuaţia curbei normalizate Wöhler (10.39), prin inversare, se obţine numărul de cicluri de amplitudine pe care le poate suporta materialul: (10.41) • Se stabilesc nivele de tensiune în diagrama seismică: (10.42) • Pentru fiecare nivel se calculează , iar apoi, apelând din nou la ecuaţia curbei normalizate Wöhler (10.39), din ec.(10.41) se obţine numărul echivalent de cicluri de amplitudine care pot provoca ruperea: (10.43) • Pentru fiecare nivel , se determină numărul de vârfuri ale diagramei seismice din acest nivel ; • Se transformă numărul de vârfuri din fiecare nivel într-un număr echivalent de cicluri de amplitudine : (10.44) • În final, numărul echivalent de cicluri rezultă prin sumarea numerelor echivalente de cicluri ale nivelelor: (10.45) Deci, având o tensogramă seismică se poate determina o oscilaţie armonică echivalentă compusă dintr-un număr de cicluri de amplitudine . Însă, materialul nu poate suporta decât un număr de cicluri de amplitudine . Prin urmare, dacă , rezultă că materialul poate suporta o încărcare cu o durată mai lungă decât cutremurul de calcul. În caz contrar, când sunt toate şansele ca un viitor cutremur să producă depăşirea rezistenţei de rupere a materialului. Fig. 10.16 Vom exemplifica metoda prezentată mai sus prin rezultatele obţinute pe probe de loess prelevate dintr-un amplasament [24], [101]. A fost utilizată tensorgama seismică din fig.10.16. În general, calculele au indicat o comportare bună a materialului solicitat de o oscilaţie armonică echivalentă formată din cicluri cu amplitudinea normalizată echivalentă . (tabelul 10.3) Însă, în cazul unui strat de loess submersat, din ecuaţia de rezistenţă dinamică (10.45) a rezultat că materialul poate suporta numai un număr maxim de cicluri cu amplitudinea . Deci, cum , apar depăşiri locale ale rezistenţei de rupere ce pot conduce la tasări suplimentare. Tabelul 10.4 Date din tensograma seismică Date din curba Wöhler ∆neq Nivel Ec.(10.45) 0,0485 0,681 0,5 4,458 8,907 0,0433 0,613 4,5 7,657 46,675 0,0385 0,545 6,5 16,106 32,052 0,0337 0,477 9,0 53,308 13,408 0,0289 0,409 14,0 1632,160 0,681 0,0241 0,341 32,5 - - = 102 10.7. Lichefierea. Caz particular de rupere dinamică Cazul cel mai spectaculos de depăşire a rezistenţei mecanice a pământurilor supuse acţiunii seismice este lichefierea. Multe din cutremurele majore au fost însoţite de distrugeri apreciabile datorate lichefierii [80], [113], Lichefierea apare deobicei în pământurile necoezive saturate, cu o anumită granulaţie şi densitate, situate în straturi cu adâncimi de peste 15 m[113]. Din cauza saturării şi a absenţei drenajului undele seismice de forfecare pot provoca, în funcţie de amplitudinea şi durata lor, o creştere a presiunii apei din pori peste presiunea medie efectivă preexistentă. Din acest moment, materialul îşi pierde capacitatea de a prelua tensiuni tangenţiale şi se comportă ca un fluid, de unde denumirea de „lichefiere” dată acestui proces. Au fost propuse multe modele pentru lichefiere. În marea lor majoritate, aceste modele reflectă mecanismul de producere menţionat mai sus constând în evaluarea creşterii presiunii apei din pori şi detectarea momentului în care egalează presiunea medie efectivă preexistentă. În acest scop, sunt evaluate tensiunile date de undele seismice de forfecare din stratul cu condiţii de lichefiere şi se echivalează tensograma obţinută cu o oscilaţie armonică echivalentă care se compară cu „curba de lichefiere” (amplitudinea excitaţiei versus numărul de cicluri care provoacă lichefierea), determinată experimental. Determinarea în laborator a curbei de lichefiere, de fapt o curbă de rezistenţă la lichefiere, se face prin încercări dinamice pe probe prelevate din amplasament. Probele sunt supuse unor tensiuni sau deformaţii armonice detectând numărul de cicluri la care presiunea apei din pori u egalează presiunea din celulă , moment în care proba lichefiază. Repetând aceste încercări cu diverse amplitudini τ (sau γ) se obţine curba de rezistenţă la lichefiere, utilizată de obicei în formă normalizată în raport cu tensiunea medie existentă „in situ”: (10.46) Curba cu ec.(10.46) poate fi privită ca o curbă Wöhler normalizată în raport cu . Ca şi în cazul materialelor nelichefiabile, lichefierea apare ca un ultim stadiu de degradare produsă de încărcările dinamice. Creşterea presiunii apei din pori apare treptat la fiecare ciclu iar capacitatea materialului de a prelua tensiuni tangenţiale se reduce treptat cu fiecare ciclu până la epuizarea totală numită lichefiere. Cu toate că degradarea şi ruperea materialelor nelichefiabile au un mecanism diferit de cel al lichefierii materialelor necoezive, se remarcă o similitudine fenomenologică în comportarea acestor materiale diferite dar supuse la aceleaşi încărcări dinamice. Din acest motiv, conceptul de degradare utilizat la modelarea ruperii materialelor nelichefiabile poate fi folosit ca atare şi în modelarea lichefierii, care apare ca un caz particular de rupere dinamică [31]. Astfel, transformând o curbă de tipul ec.(10.46) în una de forma , prin multiplicare cu , şi renormalizând în raport cu rezistenţa de rupere statică , se obţine o curbă de rezistenţă dinamică de forma ec.(10.39). Pentru exemplificare, în tabelul 10.5 este dată o astfel de transformare obţinută prin prelucrarea datelor prezentate în [83], [84]. După cum se observă, curbele obţinute prin transformare pot fi satisfăcător aproximate de exponenţiala (10.39) determinată pentru materialele nelichefiabile. Tabelul 10.5 Material n τ / σ0 r r ec.(10.39) Nisip de Fuji [83] = 0,10 MPa = 0,020 MPa 1 5 10 15 20 50 100 0,580 0,355 0,305 0,285 0,270 0,240 0,222 1,000 0,612 0,526 0,491 0,466 0,414 0,383 1,000 0,610 0,525 0,488 0,467 0,418 0,395 Nisip de Monterey [84] = 0,056 MPa = 0,020 MPa 1 4 9 16 63 0,357 0,237 0,206 0,183 0,148 1,000 0,664 0,576 0,511 0,415 1,000 0,679 0,573 0,507 0,416 1111 MODEL VÂSCOELASTIC NELINIAR PENTRU MATERIALE ANTIVIBRATORII 11.1. Introducere Materialele folosite în izolări antivibratorii au câteva caracteristici mecanice esenţiale, dintre care menţionăm, în primul rând, neliniaritatea şi capacitatea de disipare, motiv pentru care în modelarea răspunsului dinamic este necesar, în principiu, un model vâscoelastic neliniar (vezi, cap.7). Cum la pământuri, materiale cu mult mai dificile decât materialele antivibratorii, asemenea modele şi-au dovedit utilitatea vom prezenta în acest capitol adaptări ale modelelor vâscoelastic neliniare ale pământurilor, adaptări care să le facă aplicabile în evaluarea comportării materialelor antivibratorii. Deşi mult mai adecvat comportării reale, un model vâscoelastic neliniar conduce la relaţii constitutive ceva mai complicate. Neliniaritatea relaţiei constitutive tensiune-deformaţie-timp face ca ecuaţia de mişcare să fie puternic neliniară, iar datorită componentei vâscoase aceste ecuaţii devin ecuaţii diferenţiale neliniare. Rezolvarea numerică a acestor ecuaţii necesită liniarizarea, globală [58], [122] sau pe porţiuni [11], [33], [38], [41], [71], [141], procedeu laborios şi cu un control al erorilor dificil de realizat. În plus, în cazul materialelor izolatoare integrate într-un sistem de control al vibraţiilor timpul necesar rezolvării unor sisteme neliniare compromite metoda de control [54], [126]. Din acest motiv, considerăm pe deplin justificată existenţa în acest caz a unor metode de liniarizare controlată a comportării reale neliniare bazată pe utilizarea modelului vâscoelastic neliniar şi a posibilităţilor oferite de aparatura dinamică existentă . Controlul calitativ al liniarizării poate fi obţinut prin realizarea şi modelarea neliniară a experimentului dinamic, iar determinarea efectivă a parametrilor liniari ai modelului echivalent pot fi obţinuţi din condiţia ca soluţia liniarizată să fie apropiată celei neliniare. Şi în acest caz, metodele de liniarizare concepute pentru pământuri, metode prezentate anterior în capitolul 9, pot şi vor fi adaptate comportării materialelor antivibratorii. În principiu, liniarizarea echivalentã urmeză aceiaşi cale şi este bazată pe înlocuirea ecuaţiei diferenţiale neliniare a sistemului oscilant printr-o ecuaţie diferenţială liniară, cu condiţia ca determinarea coeficienţilor liniari echivalenţi sã fie astfel făcută încât soluţiile celor două ecuaţii să difere prin termeni de ordinul doi [39], [47], [124]. Problema determinării constantelor echivalente poate fi tratată ca o problemă de tip invers : cunoscând excitaţia (intrarea) ca şi soluţia sistemului neliniar (ieşirea) se determină caracteristicile sistemului. Avantajele utilizării la calculul ansamblului structură – izolaţie – fundaţie a unor constante în locul unor funcţii neliniare este evident, simplificând considerabil efortul de calcul prin înlocuirea ecuaţiilor neliniare cu unele liniare. De fapt, metoda de liniarizare propusă transferă dificultăţile de la utilizator la laboratorul care operează determinările şi prelucrările experimentale. Depăşirea dificultăţilor este astfel mult înlesnită, deoarece în laborator există posibilitatea verificării riguroase a modelării răspunsului probei prin compararea permanentă a datelor experimentale cu cele obţinute prin calcul, ceea ce pentru un sistem real, în mărime naturală, este imposibil de realizat. Structura acestui capitol este axată pe demonstrarea posibilităţii de realizare a unei astfel de liniarizări echivalente a comportării neliniare pe care o manifestă materialele utilizate ca izolatori antivibratorii. În primul paragraf se va prezenta o posibilitate de modelare a răspunsului neliniar la o solicitare dinamică, modelare asemănătoare celei utilizate la pământuri. Însă, la pământuri solicitarea dinamică predominantă este de torsiune în timp ce solicitările materialelor antivibratorii sunt, de regulă, date de forţe şi nu de momente. Apoi, va fi prezentată o metodă de liniarizare adaptată acestor materiale. Metoda de liniarizare este ilustrată cu determinări aplicative în care au fost folosite datele experimentale obţinute pe probe de cauciuc moale şi de elastomer, materiale curent folosite cu rol anivibrator. Menţionăm că metoda prezentată este relevantă pentru orice material izolator. Indiferent de natura concretă a materialului folosit, dacă răspunsul dinamic poate fi modelat cu un model vâscoelastic neliniar (şi marea majoritate a acestor materiale admit această modelare), atunci procedeul expus nu se modifică cu nimic. Sunt necesare numai date experimentale obţinute pe probe din materialul respectiv iar prelucrarea lor rămâne aceiaşi. 11.2. Modelarea răspunsului neliniar Considerăm un sistem oscilant cu un grad de libertate dinamică solicitat de o forţă unidirecţională: . Dacă proba conţine un material neliniar ecuaţia de mişcare a sistemului este [39]: (11.1) unde m este masa vibratorului, caracteristica de amortizare a materialului din probă, caracteristica de rigiditate a materialului din probă, amplitudinea excitaţiei iar ω pulsaţia ( , unde f este frecvenţa excitaţiei). Un sistem similar, dar liniar, are ecuaţia de mişcare de forma [48], [51]: (11.2) În dinamica liniară constanta elastică k este dată, în acest caz, de relaţia: [N/m] (11.3) unde E este modulul Young, A este aria secţiunii :, h înălţimea probei iar φ diametrul ei. Coeficientul de amortizare c exprimat în funcţie de raportul de amortizare D este : , şi cum , obţinem: [Ns/m] (11.4) Comportarea neliniară a materialului este reliefată de dependenţa modulului longitudinal E şi al raportului de amortizare D de nivelul de deformaţie ε, dependenţă definită prin introducerea funcţiei-modul şi a funcţiei de amortizare . De asemenea, ţinând cont de relaţia de legătură dintre deformaţia longitudinală ε şi deplasarea longitudinală x : funcţiile dinamice de material pot fi definite şi în raport cu deplasarea x : şi . Prin urmare putem admite că dependenţa de ε sau x afectează şi caracteristica elastică k: (11.5) şi caracteristica de amortizare c: (11.6) iar aceste expresii sunt folosite în ec.(11.1). În cazul sistemelor liniare pulsaţia proprie ω0 este definită în funcţie de caracteristica elastică k : , care, în cazul unui sistem neliniar este o funcţie şi nu o constantă. Din acest motiv, vom defini drept pulsaţie proprie a sistemului neliniar, pulsaţia corespunzătoare valorii iniţiale : (11.7) Efectuând schimbarea de variabilă şi introducând o nouă funcţie: (11.8) prin derivare în raport cu noua variabilă, se obţine: (11.9) astfel că ec.(11.1) devine: (11.10) unde: (11.11) În cele ce urmează vom exemplifica determinarea răspunsului neliniar utilizând rezultatele unui test în coloana rezonantă pe o probă de cauciuc moale în greutate de 147 g şi cu dimensiunile : şi . Prelucrarea datelor experimentale, a condus la obţinerea funcţiilor dinamice (fig. 11.1): (11.12) Ţinând cont de dimensiunile şi masa probei şi prin utilizarea funcţiilor dinamice de material şi din relaţiile 11.11 şi 11.12 se obţin funcţiile neliniare ale sistemului adimensionalizat din ec.(11.10) (fig.11.2): (11.13) Fig. 11.1 Fig. 11.1 Fig. 11.2 În funcţie de parametrii solicitării exterioare : amplitudinea normalizată µ şi pulsaţia relativă ν, ecuaţia neliniară (11.10) poate fi rezolvată numeric (vezi cap.3), obţinând o soluţie de forma dată în n valori discrete, iar prin prelucrarea acestora se pot determina curbele de amplificare (fig.11.3). Liniarizarea echivalentă urmează calea expusă în cap.9. În acest caz, utilizând constantele de material din tabelul 11.1 au fost obţinute constantele liniar echivalente din tabelul 11.2 iar curbele de amplificare liniarizate sunt date în fig.11.4. Şi în acest caz se poate constata că liniarizarea conduce la curbe de amplificare sufucient de apropiate de cele obţiute prin rezolvarea neliniară. Fig.11.3 Tabelul 11.1 a1 a2 a3 b1 b2 b3 0,32 0,28 105 0,25 0,75 105 Tabelul 11.2 µ c k 0,131 0,550 0,208 0,364 0,262 0,280 0,310 0,251 0,314 0,245 Fig. 11.4 11.3. Utilizarea relaţiilor experimentale forţă-deplasare Metoda de determinare a răspunsului neliniar şi de liniarizare expusă anterior presupune implicit posibilitatea de a realiza şi testa epruvete din materialul a cărei comportare dinamică se doreşte a fi cuantificată. Însă, în aplicaţiile practice nu întotdeauna operaţiile de prelevare-prelucrare a epruvetelor sunt accesibile. În plus, în cazul izolatorilor de mici dimensiuni care pot fi testaţi în aparatul de laborator sau în situ ca atare, confecţionarea de epruvete poate fi evitată iar determinările pot fi realizate pe baza datelor experimentale obţinute prin testarea întregului corp nu numai a unei părţi cum este epruveta. În cazul în care se testează întregul izolator, relaţiile experimentale obţinute sunt de tipul forţă-deplasare şi datorită geometriei mai complicate a corpului trecerea în relaţii de tip tensiune-deformaţie nu mai este atât de facilă ca în cazul epruvetei. De aici, necesitatea de a adapta metodele de cuantificare la acest tip de relaţii experimentale. O astfel de adaptare va fi prezentată în acest capitol, utilizând datele experimentale obţinute prin încercarea pe standul de încercări Schenck al Institutului de Mecanica Solidelor a unui izolator din elastomer. Elementul încercat face parte dintr-un distanţier-antivibrator pentru cablurile paralele ale liniilor electrice aeriene destinat atât menţinerii distanţei dintre două cabluri cât şi diminuării vibraţiilor induse de diverse cauze exterioare. Distanţierul este compus dintr-o tijă de metal care se prinde de fiecare cablu prin două bucşe din elastomer, elemente cu rol disipativ. Bucşele au o geometrie complexă, impusă de necesitatea unei prinderi ferme între tijă şi cable. Însă, dimensiunile relativ reduse ale ansamblului tijă-cable au permis montarea în aparatul de încercare şi încercarea întregului ansamblu la vibraţii axiale. Fig. 11.5 11.3.1. Funcţiile elastice şi de amortizare Considerând ansamblul tije-bucşe solicitat în regim staţionar de vibraţii axiale armonice ca un sistem vibrant neliniar cu un grad de libertate, ecuaţia de mişcare poate fi scrisă în forma (11.1): (11.14) unde m este masa celor 4 bucşe (160 g), este caracteristica neliniară de amortizare iar este caracteristica elastică, de asemenea neliniară. Datele experimentale obţinute prin solicitarea la vibraţii axiale, date care reprezintă răspunsul celor 4 bucşe de elastomer, sunt prezentate în fig.11.5 în forma unui răspuns histeretic forţă-deplasare, din care pot fi determinate caracteristicile elastice şi de amortizare. Fig. 11.6 Asfel, „curba-schelet” a buclelor de histerezis, curbă care este locul geometric al punctelor de extrem ale buclelor reprezintă forţa elastică, şi are în cazul de faţă o expresie analitică neliniară obţinută prin „best-fit” de tip sigmoidal (fig.11.6): (11.15) De aici, caracteristica elastică poate fi imediat obţinută în forma (fig.11.7): (11.16) Fig. 11.7 Caracteristica neliniară de amortizare poate fi obţinută prin intermediul formei neliniare a „raportului” de amortizare : (11.17) iar acest „raport”, variabil în cazul materialelor neliniare, se poate obţine din curbele de histerezis cu relaţia [80]: (11.18) În relaţiile de mai sus, (11.17) şi (11.18), au fost folosite notaţiile obişnuite: ω0 – pulsaţia proprie, definită prin ; – energia disipată, egală cu aria buclei de histerezis şi W – energia elastică maximă corespunzătoare nivelului de deformaţie atins. Utilizând curbele de histerezis din fig.11.5 s-au determinat, pentru fiecare buclă, cu relaţia (11.18) valori ale raportului de amortizare D iar prin prelucrarea acestor valori s-a obţinut funcţia de amortizare (fig.11.8): (11.19) Fig. 11.8 11.3.2. Determinarea răspunsului neliniar Prin aceiaşi schimbare de variabilă τ = ω0t (vezi capitolul 3) poate fi determinată forma adimensionalizată a ecuaţiei de mişcare (11.14): (11.20) unde, în acest caz, noile caracteristici neliniare pot fi determinate cu ajutorul funcţiilor neliniare şi obţinute anterior în forma din relaţiile (11.16) şi (11.19), astfel: (11.21) Ecuaţia neliniară (11.20) poate fi acum rezolvată numeric şi obţinută soluţia în n valori discrete. Folosind, ca în paragraful 11.1, factorul de amplificare maxim ca ilustrare a comportării sistemului dinamic, prin rezolvări numerice repetate cu diferite valori µ şi ν, se obţin funcţiile de amplificare din fig.11.9. Fig. 11.9 11.3.3. Liniarizarea Şi în acest caz liniarizarea caracteristicilor dinamice poate urma calea expusă în cap.9. Folosind aceiaşi metodă de liniarizare globală în regim de rezonanţă au fost obţinute constantele liniar echivalente din tabelul 11.3. Tabelul 11.3 F0 [N] µ C K 375 0,0005 0,162 0,648 750 0,001 0,182 0,616 1500 0,002 0,193 0,604 3000 0,004 0,200 0,601 6000 0,008 0,203 0,600 Funcţiile de amplificare liniarizate calculate cu aceşti coeficienţi liniar echivalenţi sunt date în fig.11.10 Fig. 11.10 1122 EVALUAREA EFECTELOR STRUCTURALE 12.1. Introducere Marea majoritate a materialelor din amplasamentele construcţiilor au caracteristici dinamice neliniare, mai mult sau mai puţin pronunţate. Numeroase teste, efectuate în special în coloana rezonantă [3], [6], [7], [8], au relevat modificări apreciabile ale valorilor iniţiale corespunzătoare unor deformaţii reduse, modificări care pot ajunge, în funcţie de materialul testat, până la 70 – 80%. Este de aşteptat ca asemenea modificări ale valorilor de rigiditate şi amortizare să aibă un efect apreciabil în răspunsul dinamic structural. În capitolele anterioare comportarea nelinară a fost modelată, în principal, prin considerarea acestor materiale drept materiale vâscoelastic neliniare. Modelul dinamic rezultat (capitolul 7) este construit pe baza a două funcţii dinamice, una care modelează caracteristicile de rezistenţă – funcţia-modul dinamic , iar a doua, funcţia de amortizare care include proprietăţile de disipare, ambele funcţii fiind dependente de nivelul de deformaţie γ. Aceste funcţii de material pot fi cuantificate experimental, în principal, prin teste efectuate în coloana rezonantă. Se constată experimental că evoluţia acestor funcţii de material în raport cu amplitudinea excitaţiei dinamice este diferită – funcţia-modul dinamic este descrescătoare (datorită fenomenului de degradare dinamică) iar funcţia de amortizare creşte pe măsură ce nivelul de deformaţie creşte. Aceaste evoluţii contradictorii ale funcţiilor de material au efecte contradictorii în răspunsul dinamic al sistemului oscilant în a cărui componenţă intră şi astfel de materiale. Este şi cazul masivelor din amplasamentul construcţiilor care sunt formate din materiale cu evidente caracteristici neliniare. În acest caz, răspunsul dinamic al întregului masiv este tributar atât neliniarităţilor materialelor constituente, cât şi efectelor contradictorii amintite mai sus, amplificarea sau disiparea structurală rezultând ca un compromis între aceste tendinţele divergente. De asemenea, amplificarea sau atenuarea într-un sistem oscilant neliniar prezintă caracteristici diferite faţă de sistemele liniare ceea ce face ca răspunsul neliniar să difere, atât calitativ cât şi cantitativ, de răspunsul obţinut prin neglijarea neliniarităţii. Pentru a testa modificările calitative şi cantitative pe care le aduce în răspunsul structural prezenţa şi sensul neliniarităţii caracteristicilor dinamice de material poate fi utilizat modelul Kelvin-Voigt neliniar bazat pe sistemul oscilant al coloanei rezonante (paragraful 7.3). Utilizând probe din materiale diferite prin compararea datelor de răspuns experimental cu soluţia ecuaţiei diferenţiale neliniare a sistemului oscilant se obţine o imagine pertinentă asupra modificărilor aduse răspunsului structural. 12.2. Funcţii de material Pentru a ilustra această metodă de evaluare a efectelor structurale vom folosi două materiale cu proprietăţi mecanice distanţate – nisip, material care face parte din grupa „pământuri” şi o rocă, calcarul. Prin testarea probelor prelevate din materialele respective au fost obţinute următoarele funcţii dinamice, pentru nisip (fig.12.1): (12.1) şi pentru calcar (fig.12.2): (12.2) După cum se observă din figurile 12.1 şi 12.2, nisipul are, în comparaţie cu calcarul, o funcţie-modul dinamic cu valori mai reduse dar o funcţie de amortizare cu valori superioare amortizărilor din calcar. Ca o consecinţă a acestor proprietăţi mecanice diferite şi caracteristicile neliniare ale modelului Kelvin-Voigt neliniar vor fi diferite. Astfel, pentru nisip caracteristicile neliniare, de rigiditate şi amortizare, au forma: (12.3) iar pentru calcar: (12.4) În relaţiile (12.3) şi (12.4) au fost folosite caracteristicile sistemului oscilant al coloanei rezonante (vezi cap.2): (12.5) O comparaţie a rigidităţilor este dată în fig.12.3, iar în fig.12.4 sunt ilustrate diferenţele dintre caracteristicile de amortizare. Fig. 12.1 Fig. 12.2 Fig. 12.3 Fig. 12.4 12.3. Răspunsul structural Prin rezolvarea repetată a ecuaţiei diferenţiale neliniare a sistemului oscilant probă-vibrator (2.40): pentru diferite valori ale amplitudinii normalizate µ şi pulsaţiei relative υ, se pot obţine soluţii de forma din care rezultă expresii ale „factorului de amplificare”, de fapt, funcţii de forma: (12.6) funcţii, care dau o imagine semnificativă a răspunsului dinamic. Astfel, în fig. 12.5 sunt date câteva asemenea funcţii de răspuns pentru nisip iar în fig. 12.6 pentru calcar. Fig. 12.5 Fig. 12.6 Linia punctelor de rezonanţă de rezonanţă Asemenea funcţii de amplificare structurală relevă câteva aspecte semnificative ale răspunsului dinamic structural deoarece au fost obţinute printr-un calcul neliniar, care a inclus ambele aspecte contradictorii ale comportării materialului – degradarea şi disiparea internă. La început trebuie remarcate câteva aspecte care ţin de neliniaritatea caracteristicilor dinamice ale materialului: • pentru amplitudini normalizate µ reduse (de exemplu, ) curba de răspuns are o formă apropiată de forma liniară cu un maximum de rezonanţă apropape de valoarea de rezonanţă a pulsaţiei relative (); • dacă amplitudinea excitaţiei µ creşte, amplitudinile de rezonanţă scad iar curba de amplificare ia o înclinare spre stânga, aspect caracteristic al materialelor neliniare cu caracteristică de rigiditate descrescătoare. În ceea ce priveşte comportarea diferită a celor două materiale analizate – nisip şi calcar – se observă câteva aspecte interesante care pot fi extinse la pământuri şi roci : • răspunsul dinamic structural este amplificat în toate cazurile dar valoarea amplificării diferă de la un material la altul ; în toate cazurile amplificarea este descrescătoare în raport cu amplitudinea excitaţiei ; • amplificarea dinamică este mult mai evidentă la materialele cu amortizări interne reduse, chiar dacă acestea au rezistenţe mecanice superioare; • în comparaţie cu rocile (calcar) pământurile (nisip), deşi materiale cu rezistenţe dinamice inferioare, au amplificări mai reduse datorită capacităţii superioare de disipare internă; • la toate materialele se constată o reducere a amplificărilor la rezonanţă pe măsură ce amplitudinea excitaţiei creşte; • la pământuri se constată o lărgire apreciabilă a plajei de frecvenţe care pot conduce la rezonanţă la diferite amplitudini ale excitaţiei; acest efect neliniar este prezent şi la roci însă plaja de valori υ este mult mai limitată. BIBLIOGRAFIE 1. ALLAM, M.M., SRIDHARAN, A., Efect of wetting and drying on shear stength, Journ. of the Geotechn. Eng.Div., 1981. 2. ANDERSEN, K.H., .JACOBUS, H.P., ROSENBRAND, W.F., Cyclic and static laboratory tests on dramen clay, Journ. of the Geotechn. Eng.Div., 1980. 3. ANDERSON, D.G., STOKOE, L.H., Shear modulus : a time-dependent soil property in dynamic geotechnical testing, Amer.Soc. of Testing and Materials, 1978. 4. ANNAKI, M., KENNETH, L.L., Equivalent uniform cycle concept for soil dynamics, Journ. of the Geotechn. Eng.Div., 1977. 5. AUBRY, D., KODAISSI, E., A viscoplastic constitutive equation for clay including a damage law, 5th Int.Conf. on Num.Meth. in Geomech., Nagoya, April, 1985. 6. BARDET, J.P., A viscoelastic model for the dynamic behavior of saturated poroelastic soils, Transaction of the ASME, 59, 128, March, 1992. 7. BELYTSCHKO, T., MULLEN, R., On dispersive properties of finite element solution in modern Problems in Elastic Wave Propagation, J.Willey, New York, 1977. 8. BESKOS, D.H., VGENOPOULOU, I., Wave propagation in satured, fissured, poroelastic rocks, Proc.of the European Conf. on Structural Dynamics, Bochum, June, 1990. 9. BILLINGTON, E.W., TATE, A., The physics of deformation and flow, McGraw-Hill, 1981. 10. BRATOSIN, D., Un model elastic neliniar pentru pământuri, St. şi Cerc. de Mec. Aplic., 4, 1977. 11. BRATOSIN, D., Un procedeu de aproximare a comportării unui corp elastic neliniar, St.şi Cerc. de Mec.Apl., 1, 1978. 12. BRATOSIN, D., Caracterizarea dinamică a pamânturilor prin încercări pulsatorii triaxiale, Vol.1 al Lucrărilor celei de a 5-a Conferinţe Naţionale de Geotehnică şi Fundaţii, Cluj-Napoca, 1983. 13. BRATOSIN, D., MĂRMUREANU, GH., COJOCARU, EMILIA, VASILE, I., BĂLAN, FL.ST.,Evaluarea în coloana rezonantă a caracteristicilor dinamice ale pământurilor, Construcţii, 11, 1983. 14. BRATOSIN, D., MĂRMUREANU, GH., VASILE, I., COJOCARU, EMILIA, BĂLAN, FL.ST.,Utilizarea coloanelor rezonante în caracterizarea dinamică a pământurilor, Vol.1 al Lucrărilor celei de a 5-a Conferinţe Naţionale de Geotehnică şi Fundaţii, Cluj-Napoca, 1983. 15. BRATOSIN, D., Modelarea comportării neliniare a probei din coloana rezonantă, Vol.1 al Lucrărilor celei de a 5-a Conferinţe Naţionale de Geotehnică şi Fundaţii, Cluj-Napoca, 1983. 16. BRATOSIN, D., O metodă de evaluare a oboselii rocilor, St.Cerc.Mec.Apl., 5, 1983. 17. BRATOSIN, D., Contribuţii la calculul stării de tensiune şi deformaţie în masive de pământ ţinând seama de caracteristicile neliniare şi reologice ale acestora,Teză de doctorat , Institutul Naţional de Fizică (ICEFIZ), l983. 18. BRATOSIN, D., Nonlinear relationships between stress and strain invariants of soils, Rev.Roum.Sci.Tech.-Mec.Appl., 30, 6, 1985. 19. BRATOSIN, D., On nonlinear behaviour of soils under dynamic loads, Proceedings of the EUROMEC 196 Symposium Rock and Soil Rheology, Bucharest, 1985. 20. BRATOSIN, D., VASILE, I., Vibration in nonlinear viscoelastic soils, Proceedings of the First Symposium Dynamics of Machine Foundations, Bucharest, 1985. 21. BRATOSIN, D., Mechanical model for dynamic behaviour of soils, Preprint ICEFIZ, EP - 30 – 1986. 22. BRATOSIN, D., Nonlinear behaviour of soils during harmonic vibration, Rev.Roum.Sci.Tech.-Mec.Appl.,31, 2, 1986. 23. BRATOSIN, D., Nonlinear viscoelastic model fors Soils,Rev.Roum.Sci.Tech.-Mec.Appl., 31, 1, 1986. 24. BRATOSIN, D, VASILE, I., BĂLAN, FL.ST., Rezistenţa dinamică a loessului, Vol.1 al Lucrărilor celei de a 6-a Conferinţe Naţionale de Geotehnică şi Fundaţii, Galaţi, 1987. 25. BRATOSIN, D., Pulsating torsional test of soils, Rev.Roum.Sci.Tech.-Mec.Appl., 32, 2, pp.229-240, 1987. 26. BRATOSIN, D., VASILE, I., COJOCARU, EMILIA, O lege constitutivă dinamică pentru loess, Vol.1 al Lucrărilor celei de a 6-a Conferinţe Naţionale de Geotehnică şi Fundaţii, Galaţi, 1987. 27. BRATOSIN, D., On failure strength of soils, Preprint ICEFIZ EP - 37 - 1988. 28. BRATOSIN, D., Rezistenţa limită a pămînturilor cu proprietăţi reologice, Simpozion ICEFIZ Progrese în Fizică, Constanţa, 1988. 29. BRATOSIN, D., Mechanics of resonant column test, Preprint ICEFIZ EP - 38 - l988. 30. BRATOSIN, D., On static failure of soils, Rev.Roum.Sci.Tech.-Mec.Appl., 33, 6, pp.605-618, l988. 31. BRATOSIN, D., On dynamic failure of soils, Rev.Roum.Sci.Tech.-Mec.Appl., 34, 1, 1989. 32. BRATOSIN, D., A dynamic constitutive equation for soils , Rev.Roum.Sci.Tech.-Mec.Appl., 38, 5, pp.515-519, 1993. 33. BRATOSIN, D., Dynamic finite elements for soils, Rev.Roum.Sci.Tech.-Mec.Appl., 39, 3, 1994. 34. BRATOSIN, D., On wave propagation through soils, Rev.Roum.Sci.Tech.-Mec.Appl., 39, 1, 1994. 35. BRATOSIN, D., Aspecte neliniare în mecanica pământurilor, Ed.Tehnică, 1996. 36. BRATOSIN, D., Nonlinear analysis of the dynamic response of soils, Proceedings of the Second Congress of Nonlinear Analysts , Atena – Grecia, 1996. 37. BRATOSIN, D., BĂLAN, F.S., The influence of soil fatigue on slope stability, Rev.Roum.Sci.Tech.-Mec.Appl., 42, 2, 1997. 38. BRATOSIN, D., Nonlinear analysis of the dynamic response of soils, Rev.Roum.Sci.Tech.-Mec.Appl., 42, 3-4, 1997. 39. BRATOSIN, D., SIRETEANU, T., BĂLAŞ, C., Equivalent linear model for nonlinear behaviour of soils, Rev.Roum.Sci.Tech.-Mec.Appl, 43, 4, 1998. 40. BRATOSIN, D., Strain rate effects on stress soil response, Rev.Roum.Sci.Tech.-Mec.Appl., 44, 4, 1999. 41. BRATOSIN, D., On the effects superposition in soils, Proceedings of the Romanian Academy, 1, 2000. 42. BRATOSIN, D., Elemente de mecanica structurilor, Ed.Fundaţiei România de Mâine, Bucureşti, 2001. 43. BRATOSIN, D., Structural damping in geological materials, Proceedings of the Romanian Academy, 1, 2001. 44. BRATOSIN, D., Evaluarea efectelor date de neliniarităţile de material asupra amplificării dinamice structurale, A doua Conferinţă Naţoinală de Inginerie Seismică, Bucureşti, 2001. 45. BRATOSIN, D., Soils fatigue modelling by dynamic degradation effect, Proceedings of the Romanian Academy, 3, 2001. 46. BRATOSIN, D., A dynamic constitutive law for soils, Proceedings of the Romanian Academy, 1-2, 2002. 47. BRATOSIN, D., SIRETEANU, T., A nonlinear Kelvin-Voigt model for soils, Proceedings of the Romanian Academy, 3, 2002. 48. BRATU, P., Vibraţiile sistemelor elastice, Ed.Tehnică, Bucureşti, 2000. 49. BUZDUGAN, GH., Izolarea antivibratorie a maşinilor, Ed.Academiei, 1980. 50. BUZDUGAN, GH., FETCU, L., RADEŞ, M., Vibraţiile sistemelor mecanice, Ed.Academiei, 1975. 51. BUZDUGAN, GH., FETCU, L., RADEŞ, M., Vibraţiile mecanice, Ed.Did.şi Ped., 1979. 52. BUZDUGAN, GH., MIHĂILESCU, ELENA, RADEŞ, M., Măsurarea vibraţiilor, Ed.Academiei, 1979. 53. CARNAHAN, B., LUTHER, H.A., WILKES, J.O., Applied numerical methods, J.Willey, New York, 1969. 54. CHASSIAKOS, A.G., MASRI, S.F., SMYTH, A.W., CAUGHEY, T.K., On-line identification hysteretic systems, Transactions of the ASME, 65, March, 1998. 55. CHEN, A.T.F., Application of modulus degradation model of clays, Journ. of the Geotechn. Eng.Div., Oct., 1982. 56. CHRISTENSEN, R.M., Theory of viscoelasticity, Academic Press, New York, 1971. 57. COLEMAN, J.D., Suction and yield and failure surface for soil in principle effective stress space, Geotechnique, 10, 1960. 58. CORNEA, I., ONCESCU, M., MĂRMUREANU, GH., BĂLAN, FL., Introducere în mecanica fenomenelor seismice şi inginerie seismică, Ed.Academiei, 1987. 59. CORR, D.T., STARR, M.J., VANDERBY JR., R., BEST, T.M., A nonlinear generalized fluid model for viscoelastic materials, Journal of Applied Mechanics, 68, September, 2001. 60. CRISTESCU, N., Dynamic plasticity, North Holland Publ., Amsterdam, 1967. 61. CRISTESCU, N., SULICIU, I., Vâscoplasticitate, Ed.Tehnică, 1977. 62. CRISTESCU, N., A viscoplastic equation for rocks, Preprint INCREST / ICEMAT, 49, 1979. 63. CRISTESCU, N., Mecanica rocilor. Modele matematice reologice, Ed.Ştiinţifică, Bucureşti, 1990. 64. DARVE, F., Contribution á la détermination de la loi rhéologique incrémentale des sols, These Univ.Grenoble, 1974. 65. DIMAGGIO, F.L., SANDLER, I.S., Material model for granular soils, Journ. of the Eng.Mech.Div, June, 1971. 66. DINCĂ, F., TEODOSIU, C., Vibraţii neliniare şi aleatoare, Ed.Academei, 1969. 67. DRNEVICH, V.P., HARDIN, B.O., SHIPPY, D.J., Modulus and damping of soils by resonant-column method, in Dinamic Geotechnical Testing, Amer.Soc. of Testing and Materials, 1978. 68. DUNCAN, J.M., CHANG C.Y., Nonlinear analysis of stress and strain in soils, Journ.of the Soil Mech.and Found.Div, 96, 1970. 69. GOLDSCHEIDER, M., GUDEHUS, G., Rectilinear extension of dry sand : testing apparatus and experimental results, 8-th Int.Conf. of Soil Mech. and Found.Eng., Moscova, 1973. 70. GUDEHUS, G., On constitutive laws for soils and rocks, 2-nd Iint.Conf. on Num.Meth. in Geomech., Blacksburg, Virginia, 1976. 71. GUDEHUS, G., (editor), Finite element in geomechanics, J.Willey, New York, 1979. 72. GURTIN, M.E., STERNBERG, E., On the linear theory of viscoelasticity, Arch.Rat.Mech.Anal., 4, 1962. 73. HAMBLY, E.C., A new true triaxial apparatus, Geotechnique, 19, 1969. 74. HARDIN B.O., BLACK W.L., Vibration Modulus of Normally Consolited Clays, Journ. of the Soil Mech. and Found.Div, 95, Nov., 1969. 75. HARDIN, B.O., DRNEVICH, V.P., Shear modulus and damping in Soils, Int.Rep.Univ.of Kentucky, 1970. 76. HARDIN, B.O., Suggested methods of tests for shear modulus and damping of soils by the resonant column, Amer.Soc. of Testing and Materials, 1970. 77. HARDIN, B.O., DRNEVICH V.P., Shear modulus and damping in soils : design equations and curves, Journ. of the Soil Mech. and Found.Div., July, 1972. 78. HAYTHORNTHWAITE, R.M., Mechanics of the triaxial test for soils, Journ. of the Soil Mech. and Found.Div., 86, 1960. 79. IDRISS, I.M., DOBRY, R., SINGH R.D., Nonlinear behavior of soft clays during cyclic loading, Journ. of the Geotechn. Eng.Div., Dec., 1978. 80. ISHIHARA, K., YOSHIDA, N., TSUJINO, S., Modelling of stress-strain relations of soils in cyclic loading, 5th Int.Conf. on Num.Meth. in Geomech., Nagoya, April, 1985. 81. ISHIHARA, K., Soil behaviour in earthquake geotechnics, Claredon Press, Oxford, 1996. 82. KAGAWA, T., KRAFT, L.M., Lateral pile response during earthquakes, Journ. of the Geotechn. Eng.Div., Dec., 1981. 83. KAGAWA, T., KRAFT, L.M., Modeling the liquefaction process, Journ. of the Geotechn. Eng.Div., Dec., 1981. 84. KATTIKAS, C.A., WYLIE, E.B., Sand liquefaction : inelastic effective stress model, Journ. of the Geotechn. Eng.Div., Ian., 1982. 85. KIRKPATRICK, W.M., The condition of failure for sands, 4-th Int.Congr. of Soil Mech.Found.Eng, 1957. 86. KO, H.Y., SCOTT, R.F., A new soil testing apparatus, Geotechnique, 17, 1967. 87. KO, H.Y., MASSON, R.M., Nonlinear characterisation and analysis of sand, Int.Conf. on Num.Meth. in Geomech., Blacksburg, Virginia, 1976. 88. KONDER, R.L., Hyperbolic stress-strain response : cohesive soils, Journ. of the Soil Mech. and Found.Div., 89, 1963. 89. LADE, P.V., DUNCAN, J.M., Cubical triaxial tests on cohesionless soils, Journ. of the Soil Mech. and Found.Div., 99, 1973. 90. LADE, P.V., DUNCAN, J.M., Stress-path dependent behavior of cohesionless soil, Journ. of the Geotechn. Eng.Div., 102, 1976. 91. LADE, P.V., MUSANTE, H.M., Failure condition in sand and remolded clay, 9-th Int.Conf. of Soil Mech. and Found.Eng., Tokio, 1977. 92. LAZAN, B.J., Damping of materials and membres in structural mechanics, Pergamon Press, New York, 1968. 93. LEVY, S., WILKINSON, J.P.D., The component element method in dynamics, McGraw-Hill, 1976. 94. LYSMER, J., UDAKA, T., TSAI, C.F., SEED B.S., FLUSH - A computer program for aproximate 3-D analysis of soil-structure interaction problems, Univ.of California, 1975. 95. MAKRIS, N., Time domain analysis of generalized viscoelastic models, Soil Dyn. and Earth.Eng., 14, 1995. 96. MALVERN, L.E., Introduction to the mechanics of a continuous medium, Prentince Hall, New Jersey, 1969. 97. MANOLIU, I., Încercarea pământurilor şi rocilor, în Încercarea materialelor, Ed.Tehnică, Bucureşti, 1982. 98. MARCUSON III, W.F., WHALS, H.E., Efects of time on damping ratio of clay, in Dynamic Geotechnical Testing, Amer.Soc. of Testing and Materials, 1978. 99. MARCUSON III, W.F., CURRO J.R., Field and laboratory determination of soil moduli, Journ. of the Geotechn. Eng.Div., Oct., 1981. 100. MĂRMUREANU, GH., BRATOSIN, D., BĂLAN, FL.ST., COJOCARU, EMILIA, The influence of shear modulus and damping ratio on amplification spectrum and response spectra, Rev.Roum.Sci.Tech. -Méc.Appl., 28, 2, l983. 101. MĂRMUREANU, GH., BRATOSIN, D., BĂLAN, FL., VASILE, I., Cercetări privind comportarea seismică a loessurilor din zona Brăila, Vol.1 al Lucrărilor Simpozionului Naţional Probleme geotehnice ale clădirilor rezistente la cutremur, Brăila, l989. 102. MĂRMUREANU, GH., BRATOSIN, D., MOLDOVEANU, CARMEN,The influence of the soil mechanical modelling on seismic earth massive response, Rev.Roum.Sci.Tech.-Mec.Appl., 37, 5, 1992. 103. MĂRMUREANU, GH., BRATOSIN, D, CIOFLAN, CARMEN, The dependence of Q with induced strain and frequences for surface layer soils from resonant columns, Proceedings of IASPEI'97, Salonic, Grecia, 1997. 104. MĂRMUREANU, GH., BRATOSIN, D., Cercetări privind comportarea seismică a masivelor de pământ şi stabilirea parametrilor antiseismici pentru protecţia clădirilor, Vol.1 Conferinţa Naţională de Inginerie Seismică, Bucureşti, sept.,1997. 105. MĂRMUREANU, GH., BRATOSIN, D., CIOFLAN, CARMEN,The dependence of Q with seismic-induced strains and frequencies for surface layers from resonant columns, Pure and Applied Geophysics, Birkhäuser Verlag, Basel, 157, 2000. 106. MROZ, Z., On stress-strain relations in soil mechanics, 1-st Baltic Conference on Soil Mechanics and Foundation Engineering, Gdansk, 1975. 107. MURAVSKII, G., FRYDMAN, S., Site response analysis using a non-linear hysteretic model, Soil Dyn. and Earth.Eng., 17, 4, June, 1998. 108. NELSON, I. , BARON, M.L., Aplication of variable moduli models to soil behavior, Int.Journ. of Solid Structures, 7, 1971. 109. NEWMARK, N.M., Failure hypothesis for soils, , Journ. of the Soil Mech. and Found.Div., 1960. 110. NOWACKI, W., Dinamica sistemelor elastice, Ed.Tehnică, 1969. 111. NOWACKI, W. , Théorie du fluage, Eyrolles, Paris, 1965. 112. PEARCE, J.A., A new triaxial apparatus, Proc.Roscoe Memorial, Cambridge, 1971. 113. PERLEA, V., PERLEA, M., Stabilitatea dinamică a pământurilor nisipoase, Ed.Tehnică, 1984. 114. PRESS, W.H., FLANNERY, B.P., TEUKOLSKY, S.A., VETTERLING, W.T., Numerical recipes. The art of scientific computing, Cambridge Univ.Press, 1987. 115. PREVOST, J.H., HÖEG, K., Plasticity model for undrained stress-strain behaviour, 9-th International Conference on Soil Mechanics and Foundation Engineering, Tokio, 1977. 116. ROSCOE, K.H., SCHOFIELD, A.N., WROTH, C.P., On the yielding of soils, Geotechnique, 1, 1958 117. ROSCOE, K.H., BURLAND, J.B., On the generalized stress-strain behaviour of wet clay, Engineering Plasticity, Cambridge Press, 1968. 118. QIAN, D., HANSEN, J.S., A tme domain substructure synthesis method for viscoelastic structures, Journ. of Appl.Mech., 62, 407, June, 1996. 119. SANDLER, I.S., DIMAGGIO, F.L., BALADI, G.Y., Generalized cap model for geological materials, Journal of the Geotechnical Engineering Division, July, 1976. 120. SCHAJER, G.S., Mohr-Coulomb failure criterion expressed in terms of stress invariants, Transactons o fthe ASME, 65, 6, December, 1998. 121. SCHANZ, M., CHENG, A.H.D., Dynamic analysis of a one-dimensional poroviscoelastic column, Transaction of the ASME, March, 2001. 122. SCHNABEL, P.B., LYSMER, J., SEED, H.B., SHAKE a computer program for earthquake response analysis of horizontally layered sites, Report No.EERC 71-12, Univ.of California, SUA, 1982. 123. SEED, H.B., IDRISS,I.M., Soil moduli and damping factors for dynamic response analyses, Int.Rep.Univ.of Clifornia, 1970. 124. SIRETEANU, T., Gündisch, O., Părăian, S., Vibraţiile aleatoare ale automobilelor, Ed.Tehnică, 1981. 125. SOÓS, E., TEODOSIU, C., Calculul tensorial cu aplicaţii în mecanica solidelor, Ed.Ştiinţifică şi Enciclopedică, 1983. 126. STEARNS, S.D., Digital signal analysis, Hayden Book Comp., New Jersey, 1975. 127. STUTZ, P., Contribution a l'étude de la loi rhéologique des milieux pulvérulents, Thése Univ.Grenoble, 1972. 128. TALAGANOV, K., Determination of liquefaction potential by cyclic strain approch, 13-th Regional Seminar on Earthquake Engineering, Istambul, 1987. 129. TALAGANOV, K, Stress-strain transformations and liquefaction of sands, Soil Dyn. and Earth.Eng., 15, 1996. 130. TEODORESCU, P.P., ILLE, V., Teoria elasticităţii şi introducere în mecanica solidelor deformabile, Ed.Dacia, 1976. 131. TEODOSIU, C., Thermodynamique, viscoélasticite et elastoviscoplasticité, Univ.de Paris, 1974. 132. TEODOSIU, C., Elastic Models of Cristal Defects, Springer-Verlag, Berlin, 1982. 133. TRUESDELL, C., NOLL W., The Nonlinear Field Theories of Mechanics, Handbuck der Physik III/3, Springer-Verlag, 1965. 134. VOINEA, R.P., VOICULESCU, D., SIMION, F.P., Introducere în mecanica solidului cu aplicaţii în inginerie, Ed.Academiei, 1989. 135. VOINEA, R.P., BRATOSIN D., Elemente de mecanica mediilor continue, Ed. EX PONTO, Constanţa, 2000. 136. WANG, G.X., KUWANO, J., Modeling of strain dependency of shear modulus and damping of clayey sand, Soil Dyn. and Earth.Eng., 18, 6, August, 1999. 137. WU, T.H., LOH, A.K., MALVERN L.E., Study of Failure Envelope of Soils, Journ. of the Soil Mech. and Found.Div., 89, 1963. 138. WU, M.M., BILLAH, K.Y.R., SHINOZUKA, M., Analitical study of the Duffing oscilator excited by colored noise using a systematic adiabatic expansion, Journ. of Appl.Mech., December, 1966. 139. ZEGHAL, M., ELGAMAL, A.W., ZENG, X., ARULMOLI, K., Mechanism of liquefaction response in sand-silt dynamic centrifuge tests, Soil Dyn. and Earth.Eng., 18, 1, January, 1999. 140. ZENG, X., ROSE, J.G., RICE, J.S., Stiffness and damping ratio of rubber-modified asphalt mixes: potemtial vibration attenuation for high-speed railway trackbeds, Journ.of Vibr. and Control, 7, 2001. 141. ZIENKIEWICZ, O.C., The finite element method in engineering science, McGraw-Hill, 1971. 142. *** Drnevich long-tor resonant column apparatus,Operating Manual, Soil Dynamics Instruments Inc, 1979. CONTENTS Preface .………………………………………………………... 9 1. General notions ………………………………………………. 11 1.1. Basic mechanical characteristics ……………………….. 1.1.1. Nonlinearity ……………………………………… 1.1.2. Rheological effects ………………………………. 1.1.3. Nonreversibility ………………………………….. 11 12 14 15 1.2. Geometry of the invariants space ………………………. 1.2.1. Some invariants of the second order tensor ……… 1.2.2. Space of invariants ……………………………….. 17 17 19 2. Single degree-of-freedom system ……………………………. 22 2.1. Linear system …………………………………………… 2.1.1. Free vibration …………………………………….. 2.1.2. Steady state conditions …………………………… 22 23 27 2.2. Nonlinear system ……………………………………….. 29 2.3. A numerical solution for nonlinear system …………….. 2.3.1. Newmark’s method ………………………………. 2.3.2. FORTRAN program ……………………………... 2.3.3. Applications ……………………………………… 31 31 33 35 3. Qualitative and quantitative evaluation tests ………………. 43 3.1. Introduction …………………………………………….. 43 3.2. Triaxial compression test ………………………………. 44 3.3. Triaxial creep test ………………………………………. 46 3.4. Resonant column test modelling ……………………….. 3.4.1. Resonant column – short description ……………. 3.4.2. Torsion of the specimen …………………………… 3.4.3. Determining torsional strain and rotation ………... 3.4.4. Determining torsional modulus …………………. 3.4.5. Torsional momentum estimate …………………… 3.4.6. Determining damping ratio ……………………… 47 47 49 51 52 54 54 4. Nonlinear elastic constitutive equations ……………………. 58 4.1. General form of the constitutive equations …………….. 58 4.2. Constitutive relationships between invariants ………….. 59 4.3. Modulus functions ……………………………………… 62 5. Modelling the rheological effects in linear media …………. 65 5.1. Creep and relaxation ……………………………………. 65 5.2. Forced vibration response ………………………………. 67 5.3. Linear constitutive laws ………………………………… 69 5.4. First order differential models ………………………….. 5.4.1 Kelvin–Voigt model ………………………………. 5.4.2 Maxwell model …………………………………… 5.4.3 Standard viscoelastic solid ………………………... 70 70 71 73 5.5. Generalized differential models ………………………… 76 5.6. Integral viscoelastic models …………………………….. 77 5.7. Steady state response …………………………………… 79 6. Nonlinear viscoelastic model ………………………………… 83 6.1. General form of the constitutive equation ……………… 83 6.2. Nonlinear relaxation functions ………………………….. 84 7. Dynamic response modelling ………………………………… 89 7.1. Nonlinear dynamic viscoelastic model …………………. 89 7.2. Dynamic model for seismic excitations ………………… 97 7.3. Nonlinear Kelvin-Voigt model ………………………….. 99 7.4. Validation of the nonlinear viscoelastic model …………. 106 8. Viscoplastic model ……………………………………………. 109 8.1. General form of the constitutive equation ……………… 109 8.2. Determining from triaxial tests …………………………. 109 8.3. Determining from dynamic tests ………………………... 114 8.4. Comparison between slow and dynamic response ……… 117 9. Linear equivalence in nonlinear dynamics …………………. 119 9.1. Introduction ……………………………………………... 119 9.2. Linearisation methods …………………………………... 9.2.1. Minimization of the differences between solutions 9.2.2. Minimization of the differences between damping and stiffness characteristics …………………….. 9.2.3. Global linearisation method …………………… 120 120 122 122 9.3. Global linearisation in the resonant conditions …………. 124 10. Dynamic failure of soils and rocks ………………………….. 129 10.1. Introduction …………………………………………… 129 10.2. Failure in static conditions ……………………………. 10.2.1. Failure criteria ………………………………… 10.2.2. A new tri-invariant criterion ………………….. 130 130 135 10.3. Effects of the rheological properties on failure ……….. 138 10.4. Determining the rocks fatigue ………………………… 10.4.1. Failure strength-density correlation …………… 10.4.2. Wöhler curve ………………………………….. 140 141 143 10.5. Dynamic degradation effect …………………………... 145 10.6. Failure in seismic conditions ………………………….. 154 10.7. Liquefaction – a particular failure case ……………….. 158 11. Nonlinear viscoelastic model for antivibratory materials …. 161 11.1. Introduction …………………………………………… 161 11.2. Nonlinear response modelling ………………………… 162 11.3. Use of the experimental force-displacement relations ... 11.3.1. Stiffness and damping functions ……………… 11.3.2. Determining the nonlinear response …………… 11.3.3. Liniarisation …………………………………... 167 169 171 172 12. Evaluation of the structural effects ………………………….. 174 12.1. Introductions …………………………………………... 174 12.2. Material functions ……………………………………... 175 12.3. Structural response ……………………………………. 178 References …………………………………………………….. 181 Contents ………………………………………………………. 187 PRINTED IN ROMANIA