UNIVERZITET U BEOGRADU GRAĐEVINSKI FAKULTET Selimir V. Lelović KONSTITUTIVNE JEDNAČINE ZA PESAK I NJIHOVA PRIMENA U NUMERIČKOJ ANALIZI PONAŠANJA TRAKASTOG TEMELJA doktorska disertacija Beograd, 2012 UNIVERSITY OF BELGRADE FACULTY OF CIVIL ENGINEERING Selimir V. Lelović CONSTITUTIVE EQUATIONS FOR SAND AND THEIR APPLICATION IN NUMERICAL ANALYSIS OF STRIP FOUNDATION BEHAVIOR Doctoral Dissertation Belgrade, 2012 Mentor dr Miloš Lazović, vanredni profesor, Univerzitet u Beogradu, Građevinski fakultet Članovi komisije: dr Miloš Lazović, vanredni profesor, Univerzitet u Beogradu, Građevinski fakultet dr Mirjana Vukićević, docent, Univerzitet u Beogradu, Građevinski fakultet dr Slobodan Ćorić, redovni profesor, Univerzitet u Beogradu, Rudarsko geološki fakultet Datum odbrane: _____________________________ Posvećeno uspomeni na profesore Stevana Stevanovića i Savu Vukelića SADRŽAJ REZIME SUMMARY 1. UVOD 1 1.1 Uvodna razmatranja 1 1.2 Struktura rada 1 2. KONSTITUTIVNE RELACIJE 3 2.1. Uvod 3 2.2. Geometrija deformacije 4 2.2.1 Vektori deformacije i tenzori deformacije 5 2.2.2 Tenzori relativne deformacije 7 2.2.3 Prikaz tenzora deformacije preko gradijenata pomeranja 8 2.2.4 Glavne deformacije. Invarijante deformacije 10 2.3 Tenzor napona 11 2.3.1 Glavni naponi. Invarijante napona 12 2.3.2 Konvencija znaka 16 2.4. Veza između napona i deformacija 17 2.4.1 Elastični modeli 17 2.5. Plastičnost. Uslov tečenja 19 2.5.1 Osnovni pojmovi u matematičkoj teoriji plastičnosti 19 2.6 Elasto-plastični modeli za tlo 27 2.6.1 Mohr-Coulomb-ov uslov tečenja 27 2.6.2 Drucker-Prager-ov uslov tečenja 30 2.6.3 Kapa model 31 2.6.4 HISS model 33 3. INKREMENTALNA FORMULACIJA OSNOVNIH JEDNAČINA STATIČKE ANALIZE 35 3.1. Uvod 35 3.2 Formulisanje rešenja nelinearnih jednačina 37 3.2.1 Totalna Lagrange-ova formulacija 37 3.2.2 Korigovana Lagrange-ova formulacija 40 3.2.3 Linearizacija jednačina ravnoteže 40 3.3 Formulacija rešenja metodom konačnih elemenata 41 3.4 Numerička integracija konstitutivnih jednačina 46 3.4.1 Algoritam povratnog preslikavanja 47 4. NELINEARNI MOHR-COULOMB-OV MODEL 54 4.1 Uvod 54 4.2 Osnovne jednačine Mohr-Coulomb-ovog modela 55 4.3 Konstitutivne jednačine NMC modela 59 4.4 Verifikacija predloženog NMC modela 62 4.4.1 Opit direktnog smicanja - Zbijeni pesak 62 4.4.2 Londonska glina 63 4.4.3 Drobljeni antracit 64 4.4.4 Standardna triaksijalna kompresija (CTC) 65 5. PRORAČUN TRAKASTIH TEMELJA NA DEFORMABILNOJ PODLOZI 67 5.1 Uvod 67 5.2 Nosivost tla 67 5.2.1 Teorija nosivosti Terzaghi-a 68 5.3 Ilustracija primene MKE na konkretnim primerima 71 5.3.1 Kruta temeljna traka na pesku 72 5.3.2 Kruta temeljna traka na glini 81 6. KRITERIJUM STABILNOSTI ELASTO-PLASTIČNOG MATERIJALA 87 6.1 Uvod 87 6.2 Konstitutivne jednačine elasto-plastičnog materijala 88 6.3 Odredivanje kritične tačke 90 6.4 Ilustracija predloženog postupka na konkretnom primeru 93 6.4.1 Test standardne triaksijalne kompresije (CTC 95 6.4.2 Test redukovane triaksijalna ekstenzija (RTE 97 6.4.3 Test standardne triaksijalne ekstenzije (CTE) 99 6.4.4 Test redukovane triaksijalne kompresije (RTC) 101 6.4.5 Test triaksijalne kompresije (TC) 103 6.4.6 Test triaksijalne ekstenzije (TE) 105 6.4.7 Test triaksijalnog prostog smicanja (SS) 107 7. ZAKLJUČNA RAZMATRANJA 109 LITERATURA 111 REZIME Rešavanjem graničnih problema fundiranja i mehanike tla u građevinarstvu razvijen je veliki broj konstitutivnih modela i svaki je uspešan za određenu vrstu tla i za specifičnu putanju opterećenja. Zbog nehomogenosti i različitosti prirode tla, još uvek ne postoji generalni model koji uspešno opisuje ponašanje različitih vrsta tla pri proizvoljnim putanjama napona. Drugim rečima, još uvek ne postoji odgovarajući spoj matematičkih formulacija deformacije tla opisanih u mehanici kontinuuma i mehanici individualnih čestica. U ovom radu su korišćeni osnovni principi nelinearne mehanike kontinuuma i matematičke teorije plastičnosti. Temeljna konstrukcija i tlo su diskretizovani konačnim elementima a konstitutivne jednačine su primenjene za rešavanje granične nosivosti trakastog temelja oslonjenog na deformabilni sloj peska. Tlo je opisano kao homogeno, izotropno i elastoplastično telo koje se monotono opterećuje u dreniranim uslovima. Nelinearni Mohr-Coulombov (NMC) konstitutivni model je razvijen u cilju formulacije problema granične nosivosti i deformacije tla usled zadatog opterećenja. U radu je NMC model formulisan na osnovu četiri polazne pretpostavke: 1. Ugao smičuće otpornosti nelinearne anvelope loma ima hiperboličku zavisnost od normalnih napona (Maksimović, 1989.). 2. Ugao dilatancije se može odrediti kao odnos između ekvivalentne srednje zapreminske i ekvivalentne smičuće plastične deformacije. 3. Ugao dilatancije u linearnom modelu smičuće otpornosti jednak je maksimalnoj dilatanciji nelinearnog modela. 4. Bazni ugao nelinearne smičuće otpornosti se računa na osnovu aproksimacije da su tangentna smičuća otpornost za nelinearni i smičuća otpornost linearnog konstitutivnog modela jednake u određenoj vrednosti normalnog napona. Na osnovu NMC modela napravljen je računarski program da bi se analizirali različiti aspekti ponašanja trakastog temelja na deformabilnom sloju peska. Za verifikaciju predloženog modela korišćeni su rezultati laboratorijskih ispitivanja koji su publikovani u stručnoj literaturi. Doktorska disertacija predstavlja nastavak istraživačkog rada prikazanog u magistarskom radu: "Primena Mohr-Coulomb-ovog modela za rešavanje nekih problema u geotehnici". Ključne reči: Tlo, pesak, ugao smičuće otpornosti, dilatancija, elastoplastičnost, Mohr-Coulomb, nelinearna anvelopa loma, metod konačnih elemenata, povratno preslikavanje, granična nosivost NAUČNA OBLAST: GRAĐEVINARSTVO UŽA NAUČNA OBLAST: GRAĐEVINSKA GEOTEHNIKA UDK SUMMARY By solving boundary value problems in foundations and soil mechanics, a great number of constitutive models is developed in civil engineering each being successful for a certain kind of soil and for a specific path of loading. Due to inhomogeneous nature of different soils, there is no one general model to describe behavior of all kinds of soil and under all possible paths of loading. In other words, there is no accepted way to connect mathematical formulations of soil deformation by continuum mechanics and particle mechanics. In presented work, basic principles of nonlinear continuum mechanics and mathematics of theory of plasticity were used. Base construction and soil were netted with finite elements while constitutive formulations were applied to strip foundation on deformable sand. Soil is described as homogeneous, isotropic, and elastoplastic body which is monotonically loaded in drained conditions. Nonlinear Mohr-Coulomb (NMC) constitutive model is developed to describe soil deformation under applied load. In the PhD dissertation, NMC model is formulated based on four assumptions: 1. Angle of shearing resistance for nonlinear failure envelope has hyperbolic dependence on normal stresses (Maksimović, 1989). 2. Dilatation angle may be looked as the ratio between equivalent mean volumetric and equivalent shear plastic deformation. 3. In linear model, dilatation angle of shearing resistance is equal to the maximum dilatation angle in nonlinear model. 4. Basic angle in nonlinear model is calculated based on approximation that tangents of shear resistance for linear and nonlinear models were equal at one specific point of normal stress. Based on described NMC model, a computer program is developed to look at different behavior of strip foundation on deformable layer of sand. To verify new model, comparison was done with results from laboratory testing published in literature. PhD thesis is a continuation of investigation presented in MS thesis: "Application of Mohr- Coulomb model for solving certain problems in geo-technique". Key words: Soil, sand, angle of shear resistance, dilatation, elastoplastic, Mohr-Coulomb, nonlinear failure envelope, finite element method, return mapping, boundary loading CIVIL ENGINEERING GEOTECHNICAL ENGINEERING UDC 1 1. UVOD 1.1 Uvodna razmatranja U radu su analizirane konstitutivne jednačine za pesak sa posebnim osvrtom na njihovu primenu u numeričkoj analizi ponašanja trakastog temelja. Analiza je zasnovana na sledećim pretpostavkama: 1. mehaničko ponašanje krupnozrnog tla opisuje se korišćenjem modela homogenog, izotropnog i elasto-plastičnog tela; 2. razmatra se monotono opterećenje u dreniranim uslovima; i 3. temeljna konstrukcija i tlo se mogu na odgovarajući način diskretizovati konačnim elementima i numerički analizirati ponašanje trakastog temelja oslonjenog na deformabilni sloj peska. Primena nelinearne teorije tamo gde zato ima opravdanje i tamo gde postoji razumevanje posledica ovakvog proračuna, je nesumnjivo na mestu. Najčešći matematički okvir je teorija plastičnosti, gde se razvoj obavlja na dva načina - sa jedne strane pokušava se definisati generalna granica tečenja u prostoru glavnih napona za razne materijale i razne režime opterećenja, dok se sa druge strane pokušava matematički korektno i efikasno odrediti priraštaj plastičnih deformacija i pripadajuće disipacije energije. Do sada je predložen veliki broj konstitutivnih modela za tlo, ali još uvek ne postoji takav model koji uspešno opisuje ponašanje različitih vrsta tla pri proizvoljnim putanjama napona. Korišćenjem osnovnih principa nelinearne mehanike kontinuuma i matematičke teorije plastičnosti u radu je razvijen nelinearni Mohr-Coulombov (NMC) konstitutivni model u cilju opisivanja ponašanja tla. Postoji niz već prihvaćenih koraka prilikom testiranja konstitutivnog modela: matematička formulacija modela, nalaženje parametara koji definišu model (određuju se iz konvencionalnih opita), provera modela u laboratoriji za različite putanje napona i upoređivanje rezultata numeričke analize sa rezultatima merenja na terenu. Granični problemi u matematičkoj teoriji plastičnosti metodološki se slično rešavaju kao i odgovarajući problemi linearne teorije elastičnosti. Osnovne jednačine MKE za probleme plastičnosti mogu da se izvedu na isti način kao i za probleme elastičnosti. Jedina razlika, kvalitativne prirode, je u tome što se u matrici krutosti konačnog elementa umesto elastične konstitutivne matrice veze pojavljuje tzv. trenutna elasto–plastična konstitutivna matrica. Pošto elementi ove matrice posredno, kao funkcije napona, zavise od pomeranja čvorova, sistem jednačina MKE u problemima matematičke teorije plastičnosti je nelinearan. 1.2 Struktura rada Rad sadržajno obuhvata sedam poglavlja i jedan prilog. Posle uvodnih razmatranja i napomena u drugom poglavlju su obrađene konstitutivne relacije. U ovom poglavlju su date 2 teorijske postavke vezane geometriju deformacije i tenzor napona, kao i osnove matematičke teorije plastičnosti. U nastavku je dat prikaz nekih materijalnih modela koji se primenjuju u geotehnici. U trećem poglavlju su date teorijske osnove inkrementalne formulacije osnovnih jednačina statičke analize. U nelinearnoj analizi konstrukcija, umesto generalisanih pomeranja, za osnovne parametre u čvorovima usvajaju se inkrementi pomeranja tako da se uobičajenom primenom MKE dobija inkrementalni oblik jednačina ravnoteže. Pretpostavljeno je da su inkrementi dovoljno mali, tako da se može izvršiti linearizacija problema za svaki inkrement opterećenja. Ako je statičko-deformacijsko stanje na početku inkrementa poznato, onda se iz uslova ravnoteže može odrediti stanje na kraju inkrementa. Polazeći od nedeformisane konfiguracije sukcesivno se određuju naponsko-deformacijska stanja na krajevima svih inkremenata. U zavisnosti od izbora referentne konfiguracije u odnosu na koju se prati inkrementalni postupak postoje dve formulacije: totalna Lagrange-ova i korigovana Lagrange-ova formulacija. U drugom delu ovog poglavlja prikazane su osnovne teorijske postavke algoritma povratnog preslikavanja za intergraciju elasto-plastičnih konstitutivnih jednačina. Četvrto poglavlje predstavlja centralnu temu rada. U okviru ovog poglavlja formulisan je NMC materijalni model u cilju opisivanja ponašanja tla. Definisani su parametri modela kao i procedura za njihovo određivanje, koristeći rezultate standardnih laboratorijskih opita. Rezultati sprovedenih proračuna su upoređeni sa merenim vrednostima tla u laboratorijskim opitima. U petom poglavlju je izvršena verifikacija razvijenog NMC modela. Poglavlje je posvećeno praktičnim aspektima, sa prikazom dobijenih rezultata primenom matematičkog modela i numeričkog postupka u analizi napona i deformacija u tlu ispod trakastog temelja. U šestom poglavlju je određeno kritično opterećenje HISS modela uz korišćenje inkrementalno iterativnog algoritma za različite putanje napona u standardnom triaksijalnom opitu. Predloženim postupakom je pokazano da je određivanje kritične tačke u analizi materijalne stabilnosti elasto-plastičnog materijala moguće svesti na određivanje onog stanja pri kome dolazi do promene znaka osnovnih kvantiteta metrike površi tečenja. Sedmo poglavlje se odnosi na zaključne napomene. Ukratko je prikazana analiza relevantnih rezultata sa posebnim osvrtom na prednosti i mogućnosti razvijenog modela. Na kraju, dat je i spisak literature koja se odnosi na razmatranu problematiku u ovom radu. 3 2. KONSTITUTIVNE RELACIJE 2.1. Uvod Pri matematičkom modeliranju bilo koje fizičke pojave, uvek se kao osnovni problem nameće idealizacija tj. pronalaženje računskog modela koji na neki način odražava suštinu pojave, njen mehanizam i sve bitne uticaje, a pri tome je dovoljno jednostavan da omogućava rutinsku upotrebu u praktičnim proračunima. Osnovni principi (balansa mase, balansa količine kretanja i balansa momenta količine kretanja) važe za sve materijale, nezavisno od njihove strukture. Kako je skup nepoznatih veličina veći od broja jednačina koje slede iz ovih zakona, potrebno je uvesti dodatne jednačine da bi jednoznačno bilo određeno ponašanje materijala pri dejstvu spoljnih efekata, koji se matematički definišu kao granični i početni uslovi. U cilju uzimanja u obzir strukture materijala koja karakteriše ponašanje materijala, neophodno je odrediti dodatne jednačine koje nazivamo konstitivnim. Osobine materijala se uzimaju u obzir preko odgovarajućih konstitutivnih jednačina, za svaki materijal sa konstitutivnim promenjivim, ograničenim na domen njihove definisanosti koji je određen odgovarajućim fizičkim osobinama materijala. Teorija konstitutivnih jednačina zasnovana je na odgovarajućim fizičkim zakonima. S obzirom da su materijali za koje se određuju konstitutivne jednačine predstavljeni svojim matematičkim modelima, konstitutivne jednačine definišu idealne materijale. Za razliku od gotovo svih drugih materijala koji se koriste u realnim konstrukcijama, osobine tla su promenjive, tako da se uglavnom moraju meriti za svaki slučaj posebno. Ponašanje u oblasti radnih napona je vrlo kompleksno: plastične deformacije se mogu teško zanemariti, u većini slučajeva su nehomogenost i anizotropija jako naglašene, a osim toga na karakter pojave bitno utiče niz propratnih pojava, naročito prisustvo vode. Idealizovati ponašanje samog materijala u ovakvim slučajevima je izuzetno teško. Kako je nemoguće potpuno obuhvatiti složenu fizičku pojavu, ona se pojednostavljuje, shematizuje uvođenjem određenih ograničenja i pretpostavki. Na taj način izvedeni obrazac – konstitutivna jednačina, po svojoj suštini, ne predstavlja ništa drugo nego simbolički izraz tih pretpostavki. Stoga je primena opravdana samo utoliko, ukoliko su pretpostavke ostvarene. Ova činjenica se ne sme nikad gubiti iz vida. 4 2.2. Geometrija deformacije Radi potpunosti celokupnog izlaganja, u ovom odeljku se detaljnije definišu mehaničke veličine: tenzor deformacije, u vezi sa deformacijom, tenzor brzine deformacije i tenzor napona [ ]94 . Posmatramo model materijala u vidu neprekidne materijalne sredine koji zauzima određeni deo zapremine prostora ograničen konturnom površi. Takvu materijalnu sredinu nazivamo telo. Telo zamišljamo sastavljenim od beskonačno malih zapreminskih elemenata, raspoređenih neprekidno, pri čemi svakom zapreminskom elementu, koji čini okolinu posmatrane tačke tela, pripisujemo elementarnu količinu mase tela. Položaj tako definisanih materijalnih tačaka u referentnoj konfiguraciji obeležavamo koordinatama tačaka tela sa KX , (K 1, 2,3)= , sa KGJG odgovarajuće bazne vektore, sa V zapreminu tela i sa A konturu tela u referentnoj konfiguraciji. U skladu sa tim sa dV obeležavamo elementarnu zapreminu i sa dA elementarnu površinu u okolini proizvoljne tačke posmatranog tela. Koordinate KX nazivamo materijalne koordinate, jer su vezane za određene materijalne tačke tela i služe za njihovu identifikaciju, jer sve koordinantne linije u svakom trenutku vremena prolaze kroz jedne te iste tačke posmatranog tela. Sistem koordinata je vezan za telo i deformiše se zajedno sa njim. Važno svojstvo datog sistema krivolinijskih koordinata je to što se u njima koordinate proizvoljne tačke u nekom trenutku brojčano izražavaju istim koordinatama kao u nedeformisanom stanju. Potrebno je odrediti pomeranja tačaka tela, kojima je zadat prvobitni oblik, uslovi oslanjanja i opterećenja. Pri tome treba odrediti i oblik onih delova površine koji ograničavaju telo, a kojima pomeranja nisu unapred zadata. Konturni uslovi zadaju se na granicama tela kojima oblik zavisi od osnovnih veličina. Prema tome, nameće se, kao najpogodnija matematička forma – krivolinijske koordinate ix jer će u njima jednačine konture posmatranog tela imati isti oblik kao i pre deformacije. Zamislimo u posmatranom telu neki skup tačaka koji leži na jednoj liniji (materijalno vlakno) ili na nekoj površi (materijalni sloj). Pri datoj deformaciji vlakno (ili sloj) neprekidno menja svoj oblik i položaj u prostoru. Za određivanje položaja vlakna (ili sloja), kao geometrijskih tačaka, služi nam "prostorni " sistem koordinata, u tom smislu što se određene koordinate ne vezuju za određene materijalne tačke, već za određene tačke prostora. Prostorni koordinantni sistem može biti sastavljen od istih Dekartovih ili krivolinskih koordinata kao i materijalni sistem, tako da u tom slučaju materijalne koordinate predstavljaju "početne vrednosti" prostornih koordinata, vezane za referentnu konfiguraciju. Prostorni 5 koordinantni sistem obeležavamo oznakom ix , (i 1, 2,3)= , sa ig G odgovarajuće bazne vektore, sa v - zapreminu tela i sa a - konturu tela u trenutnoj konfiguraciji. U skladu sa tim sa dv obeležavamo elementarnu zapreminu i sa da - element površine tela. U toku procesa deformacije prostorne koordinate određene materijalne tačke tela predstavljaju funkciju vremena, tako da za određenu tačku tela, KX , prostorne koordinate predstavljaju njene konačne jednačine kretanja: ( ) ( )i i Kx x X , t , i, K 1, 2,3= = (2.1) gde je t vreme. Jasno je da oblik konačnih jednačina kretanja materijalne tačke zavise od početnih uslova. Relacija inverzna jednačinama (2.1) pokazuje, koja je to tačka tela za određini trenutak vremena "t" zauzima određeni položaj u prostoru. Jednačine (2.1) predstavljaju deformaciju kontunuuma. Deformacija može da se posmatra sa stanovišta materijalnih koordinata, tako i sa stanovišta prostornih koordinata, odakle sledi da moramo da uvedemo postojanje jednoznačne inverzne relacije oblika: ( ) ( )K K iX X x , t , i, K 1, 2,3= = (2.2) Oba izraza za deformaciju moraju da zadovoljavaju aksiom neprekidnosti, tj. transformacije (2.1) i (2.2) su u celoj oblasti posmatranja jednoznačne i dopuštaju neprekidne izvode do svakog potrebnog reda. Pri opisivanju procesa deformacije u deformabilnoj sredini primenjuju se dva pristupa. Pristup kada se za nezavisne promenljive usvajaju koordinate materijalnih tačaka tela označava se kao materijalna deskripcija i po tom pristupu sve veličine se definišu u određenoj tački tela za određeni vremenski trenutak. Za razliku od ovog pristupa, prostorna deskripcija za nezavisno promenljive ima prostorne koordinate. Ovaj drugi pristup opisivanja procesa deformacije određuje sve veličine u određenom položaju tačke u prostoru u trenutku vremena. 2.2.1 Vektori deformacije i tenzori deformacije Ako se usvoji materijalna deskripcija (2.1), tada je: i i K ;Kdx x dX= (2.3) gde se parcijalni izvodi i;Kx nazivaju materijalni gradijenti deformacije. Pri prostornoj deskripciji je: K K i ;idX X dx= (2.4) 6 Na taj način u materijalnoj deskripciji možemo napisati za elementarni vektor položaja: i K ;K idr x g dX= JJG JJG (2.5) pri čemu veličine, i K ;K iC x g= JG G (2.6) predstavljaju materijalne vektore deformacije, pomoću kojih je, K Kdr C dX= JJG G (2.7) Vektori deformacije KC JG zamenjuju bazne vektore pri prikazivanju elementarnog vektora deformisane konfiguracije u odnosu na materijalni sistem koordinata. U prostornoj deskripciji je: K Ki ;ic X G= JG JG (2.8) sa kojima prikazujemo i idR c dx= G G (2.9) Gradijenti deformacije i;Kx se obično označava: i i ;K Kx F≡ (2.10) pri čemu je F tenzor gradijenta deformacije. Pomoću tenzora deformacije izražava se kvadrat elementarnog rastojanja u okolini posmatrane tačke deformbilnog tela. Kvadrat rastojanja , 2 2(pp ) (dr)= JJJG JJG izračunavamo koristeći relacije (2.7) u obliku: 7 ( ) ( ) ( )2 K L K LK L KLdr C dX C dX C dX dX= ⋅ =G GG (2.11) gde tenzor, ( ) ( )i j i jKL ;K i ;L j ij ;K ;LC x g x g g x x= ⋅ =G G (2.12) predstavlja materijalni tenzor deformacije ili Grinov tenzor deformacije. Međutim, posmatrajući element rastojanja u referentnoj konfiguraciji i izražavajući kvadrat tog rastojanja polazeći od izraza (2.9), nalazimo ( ) ( )( )2 i j i ji j ijdR c dx c dx c dx dx= =G G G (2.13) gde tenzor, ( ) ( )K L K Lij ;i K ; j L KL ;i ; jc X G X G G X X= ⋅ =G G (2.14) predstavlja prostorni tenzor deformacije ili Košijev tenzor deformacije. Polazeći od izraza (2.11) i (2.13) jasno je da su oba tenzora deformacije (2.12) i (2.14) simetrični i pozitivno definitni. 2.2.2 Tenzori relativne deformacije Tenzori relativne deformacije u materijalnoj ili prostornoj deskripciji izražavaju razliku kvadrata dve beskonačno bliske tačke tela (ili tačaka u beskonačno bliskom prostornom položaju) u deformisanoj i u referentnoj konfiguraciji. Prema tome nalazimo: ( )2 2 K LKL KL(dr) (dR) C G dX dX− = −GG (2.15) gde je: KL KL KL2E C G= − (2.16) Tenzor KLE je materijalni ili Lagranžov tenzor relativne deformacije. Kada se KLC izrazi preko (2.12), sledi: ( )i jKL ij ;K ;L KL1E g x x G2= − (2.17) Izražavajući razliku (2.15) u prostornoj deskripciji, nalazimo: ( )2 2 i jij ij(dr) (dR) g c dx dx− = −GG (2.18) gde je: ij ij ij2e g c= − (2.19) 8 Tenzor ije je prostorni (Ojlerov) tenzor relativne deformacije. Kada se ijc izrazi preko (2.14), sledi: ( )K Lij ij KL ;i ; je 1 2 g G x x= − (2.20) Kako se izrazi (2.15) i (2.18) odnose na istu veličinu, to je moguće uspostaviti vezu između materijalnog i prostornog tenzora relativne deformacije. 2.2.3 Prikaz tenzora deformacije preko gradijenata pomeranja Na slici 2.2 je prikazan beskonačno blizak par tačaka P i P′odnosno p i p′ (položaji tačaka P i P′ u toku deformacije). Tačka P vrši pomeranje, ( )Ku u X , t=G G (2.21) dok je pomeranje tačke P′ : ( ) ( )K K Ku u X dX , t u X , t du= + = +G G G G (2.22) Između vektora dR G , duG i drG postoji veza: dr dR du= +GG G (2.23) Vektor pomeranja u materijalnoj deskripciji predstavljamo koordinatama Ku a u prostornom sistemu koordinatama iu , tako da je, ( )K Ku u X, t G=G G , odnosno ( )i iu u x, t g=G G (2.24) pri čemu je prema izrazima za operator paralelnog pomeranja i i K Ku g u= , odnosno K K iiu g u= (2.25) 9 Priraštaj pomeranja može se prikazati u dve varijante, na osnovu obrasca (2.24): K M ,M Kdu u G dX= GG , odnosno i j, j idu u g dx= GG (2.26) kao priraštaj pomeranja između dve tačke tela odnosno dve tačke prostora. U materijalnoj deskripciji vektorska jednačina (2.23) može biti predstavljena preko koordinata KdX , vektora dR G i izraza (2.5) za dr JJG : K M K K K ,K M KG dX u G dX C dX+ = G G G (2.27.1) i odavde ( )M M MK K ,K M M K ,KC G u G G u= + = δ +G G G G (2.27.2) gde je M,Ku kovarijantni parcijalni izvod vektora pomeranja. Na osnovu (2.27) i (2.16) dobija se: ( ) ( ) ( )M M N N MKL K K ,K N L ,L KL K,L L,K ,K M,L1 1E G u G u G u u u u2 2 = δ + ⋅ δ + − = + +  G G (2.28) gde zarezi označavaju kovarijantne parcijalne izvode po materijalnim koordinatam. Linearan deo izraza (2.28) predstavlja tenzor infinitezimalne deformacije: ( )KL K,L L,K K,L1E u u u2= + = (2.29) Na taj način je tačan izraz za tenzor deformacije dat u obliku: M KL KL ,K M,L 1E E u u 2 = + (2.30) Na sličan, u odnosu na prostorni koordinantni sistem nalazimo za Ojlerov tenzor deformacije sledeći izraz preko gradijenta pomeranja: ( )mij i, j j,i ,i m, j1e u u u u2= + − (2.31) Linearan deo izraza (2.31) predstavlja tenzor infinitezimalne deformacije, izražen preko prostornih koordinata: ( )ij i, j j,i i, j1e u u u2= + = (2.32) Tenzor konačne (relativne) deformacije je dat izrazom: m ij ij ,i m, j 1e e u u 2 = − (2.33) 10 2.2.4 Glavne deformacije. Invarijante deformacije Po definiciji, glavne deformacije su sopstvene vrednosti tenzora deformacije. Ose koje odgovaraju ovim deformacijama zovu se glavne ose, a ravni upravne na te ose su glavne ravni. Smičuće deformacije u glavnim ravnima su jednaki nuli. Komponente glavnih infinitezimalnih deformacija označićemo sa 1e , 2e i 3e uz pretpostavku da je 1e ; 2e ; 3e . Komponente glavnih infinitezimalnih deformacija su koreni jednačine: 3 2 e e ee I e II e III 0− + − = (2.34) Ovde smo sa eI , eII , eIII označili invarijante tenzora deformacije, definisane sa: e 1 2 3 ii iiI e e e e e tr(e)= + + = ≡ =∑ (2.35) ( )e 1 2 2 3 3 1 ii jj ij ji1II e e e e e e e e e e2= + + = − (2.36) ( )e kl 1 2 3III det e e e e= = (2.37) Pri deformaciji tela dolazi do promene i njegove zapremine. Odnos te promene i prvobitne (nedeformisane) zapremine: V 1 2 3 ee V V e e e I= ∆ = + + = (2.38) zove se kubna dilatacija. Tenzor deformacije možemo da razložimo na dva dela: kl 0 kl kle e e′= δ + (2.39) gde je: 0 e kke I 3 1 3e= = sferni deo tenzora deformacije, klδ je Kronekerov (Cronecker) δ − symbol, kle′ devijator tenzora deformacije. Devijatorski deo odgovara klizanju za koje je kubna dilatacija jednaka nuli. Dakle, prva invarijanta devijatora tenzora deformacije ( kle′ ) je identički jednaka nuli. Druge dve invarijante devijatora deformacija mogu da se izraze preko glavnih deformacija u obliku: ( ) ( ) ( ) ( ) ( ) ( )2 2 2 2 2 2de 1 2 2 3 3 1 1 2 31 1II e e e e e e e e e6 2   ′ ′ ′= − + − + − = + +    (2.40) ( )de kl 1 2 3III det e e e e′ ′ ′ ′= = ⋅ ⋅ (2.41) 11 Druga invarijante devijatora deformacija obuhvata uticaj smičuće deformacije, dok treća invarijante devijatora deformacija predstavlja pravac smičuće deformacije. 2.3 Tenzor napona Naponsko stanje u tački tela je u potpunosti određeno sa poznavanjem vektora napona za tri nekomplanarne ravni, koje prolaze kroz tu tačku. Stanje napona je određeno poznavanjem devet podataka (tri projekcije tri vektora), koji predstavljaju koordinate tenzora napona. Ova činjenica znači da se sa poznavanjem ovih devet podataka može odrediti vektor napona ( )( )nt za bilo koju presečnu ravan, sa normalom nG , kroz tu tačku, tj.: ( ) ij jn it t n= ⋅G (2.42) gde je: ijt tenzor napona, jn kosinusi pravca normale n G , ( )n it G vektor napona za presečnu ravan sa normalom nG . Napomenućemo da ćemo u daljim izlaganjima smatrati, osim kada naglasimo suprotno, da je tenzor napona simetričan, što sledi iz balansa momenta količine kretanja, u odsustvu naponskih spregova (tri skalarne jednačine). Prethodno definisan tenzor napona je Košijev (Cauchy) tenzor napona meren po jedinici površine deformisanog tela a referisan u sistemu baznih trijedara posle deformacije. Kako je početna konfiguracija tela poznata a samim tim poznate su i njegove granične površi, često je znatno pogodnije granične uslove izraziti u odnosu na ovaj referentni položaj. Iz tih razloga potrebno je definisati naponske veličine koje će biti merene po jedinici nedeformisane površi. Kada ovaj stav primenimo na Košijev napon dolazimo do pojma Piola-Kirhofovog (Piola-Kirchhoff) tenzora napona prve i druge vrste. Prvi Piola-Kirhofov tenzor napona može se definisati preko Košijevog napona na sledeći način: kK kl K kl 1 kL l ;l ;LT Jt X , t J T x −= = (2.43) Ovaj tenzor napona predstavlja napon u tački x , meren po jedinici površine u X . Prvi Piola- Kirhofov tenzor napona je vektor ( )T X, t, N koji je paralelan sa Košijevim vektorom napona ( )t x, t, n , ali je mera sile po jedinici nedeformisane površine. Uslov prostorne paralelnosti se može matematički izraziti na sledeći način: TdA tda= (2.44) 12 Ovako definisan tenzor napona je nesimetričan pa samim tim i nepogodan za praktičnu primenu. Nedostatak se eliminiše formulacijom simetričnog Piola-Kirhofov tenzor napona druge vrste merenog u odnosu na nedeformisanu konfiguraciju, odnosno KL K mL K L lm ;m ;l mS X T JX X t= = (2.45) Kako je naprezanje karakteristika samo trenutne konfiguracije, odnosno deformisanog stanja, i imajući na umu da su vektori napona na stranama elementarnog tetraedra za slučaj oba Piolina napona mereni u odnosu na normale nedeformisane površi, oni samim tim nemaju jasno fizičko značenje i, u praktičnim računskim analizama, moraju se izračunavati Cauchy- evi naponi. 2.3.1 Glavni naponi. Invarijante napona Služeći se analogijom možemo definisati i izvesti nekoliko pojmova koji su već dokazani u prethodnim razmatranjima o deformacijama. Tako se može dokazati da u svakoj tački tela posle deformacije postoje tri međusobno upravne ravni, za koje su svi smičući naponi jednaki nuli, dok normalni naponi poprimaju ekstremne vrednosti. Komponente glavnih napona (koje ćemo označavati sa 1σ , 2σ , 3σ i uz pretpostavku da je 1σ ; 2σ ; 3σ ) predstavljaju korene sekularne jednačine kao sopstvene veličine: ij ij 0σ −σδ = (2.46) ili u razvijenom obliku: 3 2 1 2 3I I I 0σ − σ + σ− = (2.47) Ovde smo sa 1I , 2I , 3I označili tri nezavisne invarijante tenzora napona, definisane sa: 1 1 2 3 ii iiI tr= σ +σ + σ = σ ≡ σ = σ∑ (2.48) ( )2 1 2 2 3 3 1 ii jj ij ji1I 2= σ σ +σ σ +σ σ = σ σ −σ σ (2.49) ( )3 ij 1 2 3I det= σ = σ σ σ (2.50) Tenzor napona možemo da razložimo na dva dela: ij ij ijp Sσ = δ + (2.51) 13 gde je: p sferni deo tenzora napona - koji se može izraziti kao 1p 1 3J= , klδ je Kronekerov δ − symbol i ijS devijatorski deo tenzora napona ili devijatorski tenzor napona. Devijatorski deo tenzora napona određuje odstupanje, devijaciju, od sfernog stanja napona. Ovaj deo napona izaziva promenu oblika, na pr. sfera prelazi u elipsoid, kocka u paralelopiped. Karakteristika ovog tenzora je da je njegova prva invarijanta uvek jednaka nuli. Glavni pravci devijatorskog tenzora napona, analogno tenzoru napona, određuju iz sistema jednačina: ( )ij ij jS S n 0− δ = (2.52) Na isti način se može pokazati da je karakteristična jednačina devijatorskog tenzora napona oblika: 3 2d 3d 2d 2dS I S I 0, I I− ⋅ − = = − (2.53) Ovde sa 2dI i 3dI označili drugu i treću invarijantu devijatorskog tenzora napona, koje mogu da se izraze preko glavnih devijatora napona u obliku: ( ) ( ) ( )2 2 22d 1 2 3 ij ij1 1I S S S S S2 2 = + + =  (2.54) ( )3d ij 1 2 3I det S S S S= = ⋅ ⋅ (2.55) Karakteristična jednačina devijatorskog tenzora napona može se rešiti smenim 2dS 2 I 3 sin= θ (2.56) Kubna jednačina se tada svodi na oblik: ( ) ( )3 2 32d 3d2 I 3 3sin 4sin Iθ− θ = (2.57) Koristeći trigonometrijsku identičnost: 3sin 3 3sin 4sinθ = θ− θ izraz u zagradi je jednak sin 3θ , pa se prema tome dobija 3 2 3d 2dsin 3 3 3 2 I I −θ = − ⋅ ⋅ (2.58) Uobičajeno je da se kriterijum tečenja geometrijski interpretira kao površina u naponskom prostoru i da se označi kao površ tečenja. Pošto tenzor napona σ ima šest različitih komponenata, generalno posmatrano ovakav slikoviti prikaz bi zahtevao šestodimenzionalni prostor. Praktičniji pristup je korišćenje trodimenzionalnog prostora u 14 kojem su referentne koordinate glavni naponi. Naponski prostor se definiše kao prostor čije su koordinate vrednosti glavnih napona (u svim tačkama tela), tj. svakoj tački ovog prostora odgovora stanje napona u nekoj tački tela. Ovakav prostor je prikazan na slici 2.3. Hidrostatička osa (prostorna dijagonala) je linija koja pravi jednake uglove sa koordinatnim osama (u našem slučaju sa glavnim pravcima napona). Jedinični vektor n , duž hidrostatičke ose, ima sledeće komponente: { }n 1 3 ,1 3 ,1 3= (2.59) Bilo koja ravan normalna na hidrostatičku osu zove se devijatorska ili okatedarska ravan. Devijatorska ravan koja prolazi kroz koordinantni početak zove se π -ravan. Naponsko stanje može biti prikazano tačkom P sa koordinatama 1 2,σ σ i 3σ i vektorom OP . Vektor položaja neke tačke, može da se razloži na dve komponente: OP ON NP= + (2.60) Gde je vektor ON duž prostorne dijagonale a vektor NP je u devijatorskoj ravni koja prolazi kroz tačke N i P . { }1 2 3ON OP n 1 3 3 p= ⋅ = ⋅ σ +σ +σ = ⋅ (2.61) gde je 1p 1 3 I= ⋅ srednji napon. Devijatorska ravan koja prolazi kroz tačku N je prikazana na slici 2.4. 15 Neka θ označava ugao meren od ose čistog klizanja { }2 1 31 2σ = ⋅ σ +σ , kao što je prikazano na slici 2.4, onda je: ( )1 1 2dI 3 2 3 I sin 2 3σ = + ⋅ ⋅ θ + π (2.62) 2 1 2dI 3 2 3 I sinσ = + ⋅ ⋅ θ (2.63) ( )3 1 2dI 3 2 3 I sin 2 3σ = + ⋅ ⋅ θ − π (2.64) gde je 6 6−π ≤ θ ≤ π i 1 2 3σ ≥ σ ≥ σ . Oktaedarska ravan je ravan koja gradi jednake uglove sa koordinantnim osama. 16 Napon smicanja u ovoj ravni naziva se oktaedarski napon i dat je izrazom ( ) ( ) ( )2 2 20 1 2 2 3 3 1 2d1 3 2 3 Iτ = σ −σ + σ −σ + σ −σ = ⋅ (2.65) 2.3.2 Konvencija znaka Za normalne komponente tenzora napona, s obzirom na njihovo dejstvo, koriste se i posebni nazivi i to: naponi pritiska, ako normalna komponenta deluje u suprotnom smeru od spoljne normale na površinski element, i suprotno njima naponi zatezanja. U mehanici tla naponi pritiska su dominantni, pa je usvojena konvencija računanja napona pritisaka kao pozitivnih i napona zatezanja kao negativnih normalnih napona. Za tenzor napona se koriste razna obeležavanja. Za Dekartov sistem u inženjerskoj praksi najčešći vid označavanja je: xσ , yσ , zσ , xyτ , zxτ i zyτ za simetričan tenzor napona. Na slici 2.6 je prikazana interpretacija pozitivnih komponenti tenzora napona. 17 2.4. Veza između napona i deformacija Konstatovali smo da su pri rešavanju praktičnih problema potrebne i dopunske jednačine, tj. veze između napona i deformacija. Proučavanje deformacija bilo je čisto geometrijsko: ispitujući promene oblika i zapremine tela nismo ulazili u njihovu zavisnost od spoljnih sila. Za uspostavljanje veza između deformacija i napona (a na taj način i između deformacija i spoljnih sila, jer je zavisnost između napona i spoljnih sila određena) treba da definišemo predmet našeg razmatranja, tj. realno telo, znači treba mu pripisati određene osobine. Te osobine moraju biti utvrđene neposrednim eksperimentima sa realnim telima. U tom slučaju će i zaključci koji proističu iz takve definicije moći da se primenjuju na realna tela. Zaključci i jednačine koje ćemo izvesti u daljim izlaganjima, mogu se primenjivati za realna tela ako se eksperimentalno utvrdi da imaju tu osobinu i samo u granicama u kojima je imaju. Za razliku od obrazaca matematike, koji operišu sa apstraktnim pojmovima, jednačine matematičke teorije plastićnosti postavljaju veze između stvarnih veličina: spoljnih sila, dimenzija tela i osobina materijala. Obično se usvajaju tri osnovna modela materijala: elastičan, plastičan i viskozni i odgovarajuće njihove kombinacije. U daljem tekstu će se spominjati samo elastični i plastični modeli, s obzirom da su dovoljni za odgovarajuće opisivanje problematike koja je predmet ovog rada. 2.4.1 Elastični modeli Engleski fizičar Robert Hooke je objavio 1678. god., kao rezultat svojih ogleda, da za mnoga realna tela, u izvesnim granicama, važi zakon proporcionalnosti između napona i deformacija. Hooke-ov zakon i iz njega izvedeni zaključci važe bez ograničenja samo za idealno elastično telo, koje u prirodi ne postoji. Međutim, njihova primena na realna tela, npr. tlo, je opravdana u oblastima u kojima naponi i deformacije ne prelaze određene vrednosti, koje se za pojedine materijale utvrđuju eksperimentom. Koeficijenti elastičnosti, uopšte, zavise od položaja tačke u telu i od izabranih pravaca osa kroz tu tačku. Pri izvođenju konstitutivne jednačine linearno elastičnog modela najčešće se polazi od generalisanog Hooke-ovog zakona koji se može prikazati sledećim jednačinama: ij ijkl klCσ = ⋅ε (2.66) ili inverzno 18 ij ijkl klDε = ⋅σ (2.67) gde su: ijσ komponente tenzora napona, klε komponente tenzora deformacija, ijklC komponente tenzora elastičnosti koje predstavljaju module krutosti materijala, ijklD komponente tenzora elastičnosti koje predstavljaju module fleksibilnosti materijala. Sva dalja razmatranja, u primenama, biće ograničena pretpostavkom da je telo homogeno, što znači da su koeficijenti elastičnosti isti u svim tačkama. Može se dokazati da između tih 36 koeficijenata postoji 15 veza i na taj način smanjiti broj nezavisnih koeficijenata do 21. Ograničavamo se još i pretpostavkom da je telo izotropno-elastično, što znači da ima iste elastične osobine za sve pravce kroz datu tačku. Ta pretpostavka daje niz veza među koeficijentima koje, između ostalog, obuhvataju i 15 ranije spomenutih veza. Za elastične modele napon u svakoj čestici tela jednoznačno je određen deformacijom u posmatranom trenutku tj. nezavisan je od kretanja pre dolaska u tekuću konfiguraciju. Za analizu Hooke-ovog zakona pogodno je i tenzor infinitezimalne deformacije i tenzor napona razložiti na sferi i devijatorski deo u obliku, kako je već prikazano u izrazima (2.39) i (2.51): kk kk3Kσ = ⋅ε (2.68) ij ijS 2G e= ⋅ (2.69) gde je: K - moduo kompresije i G - moduo smicanja Elastični izotropni model može biti predstavljen i preko Young-ovog modula i Poisson-ovog koeficijenta. Ove konstante se mogu dovesti u vezu sa K i G . 9KG 3K 2GE , 3K G 6K 2G −= ν =+ + (2.70) Može se dobiti i veza između konstati K i G s jedne strane i konstati E i ν sa druge strane, koristeći jednačinu (2.67), u sledećem obliku: ( ) ( ) E EK , G 3 1 2 2 1 = =− ν + ν (2.71) Iz predhodne jednačine može se dobiti tenzor izotropnih elastičnih konstanti u obliku: ( )ij ij kk ij ij kk ij ij kk ij kk ij ijkl klS 2Ge 3K 2G K Cσ = +σ δ = + ε δ = ε − ε δ + ε δ = ε (2.72) gde je : ( )ijkl ik jk ik il ij kl2C G K G3 = δ + δ + δ + δ + − δ δ   (2.73) 19 pri čemu je ( )ik jk ik il kl ij ik jk kl ij kk2 ,δ + δ + δ + δ ε = ε δ δ ε = δ ε (2.74) Gornja jednačina je simetrična po indeksima i j− , k l− ; i predstavlja najprostiji oblik veze između napona i deformacija u izotropnom elastičnom telu. Prećutno se pretpostavlja da je sferni deo tenzora napona proporcionalan sfernom delu tenzora deformacije, a isto tako i devijatorski deo tenzora napona odgovarajućem delu tenzora deformacije. Moduo kompresije predstavlja količnik iz srednjeg normalnog napona i kubne dilatacije, drugim rečima daje meru kojom se izotropni materijal suprostavlja promeni zapremine. Moduo klizanja se javlja kao koeficijent proporcionalnosti između odgovarajućih devijatorskih komponenta napona i deformacija. Govoriti o opravdanosti ili neopravdanosti primene ovog modela za proračun naponsko deformacijskog stanja u tlu nije suštinsko pitanje, već je to određivanje granica promene naponskog stanja kada je ova primena prihvatljiva [ ]141 . 2.5. Plastičnost. Uslov tečenja 2.5.1 Osnovni pojmovi u matematičkoj teoriji plastičnosti Kako se većina trodimenzionalnih teorija plastičnosti zasniva na uopštavanju jednodimenzionalne teorije, pa se pri razmatranju stanja napona i deformacija u plastično deformisanim delovima konstrukcije, polazi od idealizovanog dijagrama zavisnosti jednoosnog naprezanja i deformacija. Tipična kriva napon-deformacija dobijena u testu standardne triaksijalne kompresije je prikazana na slici 2.7.b. Opit triaksijalne kompresije se sprovodi na uzorku kod koga je visina jednaka dvostrukom prečniku, i kod koga je kontrolisan inkrementalni iznos aksijalnog pomeranja. Osim toga poznati su i granični uslovi po bočnim naponima, tokom trajanja čitavog opita, dok su nepoznate veličine devijator napona, kao i veličine bočnih pomeranja. Analizirajući pomenuti eksperiment može se uočiti sledeće ponašanje materijala uzorka: elastičan opseg – koji predstavlja linearnu vezu između napona i deformacije, pri čemu se po prestanku dejstva (zadato pomeranje ili zadato opterećenje) uzorak vraća u prvobitno stanje; plastičan opseg – slučaj gde postoje deformacije koje se ne vraćaju sa rasterećenjem uzorka. 20 Za materijal izložen spoljašnjim dejstvima (mehaničkim) kažemo da je trajno - plastično deformisan ako i posle prestanka spoljašnjeg dejstva ostaje promenjen oblik ili zapremina materijala u odnosu na neutralno stanje. Kako su komponente deformacije mehaničke veličine kojima opisujemo stanje, to znači da je u nekoj tački došlo do plastične deformacije ukoliko je, po prestanku dejstva uzroka deformisanja, bar jedna od komponenata deformacije različita od nule. Trajni deo deformacije zove plastična deformacija, a za odgovarajući uzorak kažemo da se u ovom opsegu ponaša plastično. Napon tečenja je definisan kao najveći napon koji se može postići pre nego što dođe do početne plastične deformacije. Ovakva definicija granice elastičnosti je više teoretski opis stanja jer se za realne materijale tačna vrednost napona tečenja veoma teško određuje eksperimentalnim putem. Ako su eksperimentalna posmatranja dovoljno precizna, onda se postavlja pitanje da li linearni region ili tačna vrednost napona tečenja uopšte i postoje za realni material (a naročito za peskove). Međutim, koncept tečenja ima centralnu ulogu u matematičkoj teoriji plastičnosti i standardno je usvojeno da se material idealizovano ponaša prvo u linearno-elastičnom regionu posle koga počinje plastično tečenje. Test jednoaksijalne kompresije se često koristi da se definišu neki aspekti ponašanja materijala, međutim taj test ne daje zadovoljavajuće odgovore za slučajeve višeosnog naponskog stanja u kome se najčešće nalazi temeljno tlo. Postavljaju se sledeća praktična pitanja: Šta se dešava ako na element određenog materijala deluju nekoliko operećenja u različitim pravcima? Koja kombinacija ovih opterećenja je merodavna i dovodi do pojave plastičnog tečenja? Kombinacija opterećenja 21 koja izaziva pojavu tečenja u materijalu naziva se kriterijum tečenja. U opštem slučaju, kriterijum tečenja trebalo bi da zavisi od dejstva primenjenog stanja napona i da odražava efekat prethodnih opterećenja na mikrostrukturu tog elementa. Čestice granularnog skeleta tla se pod uticajem smičućih napona pomeraju jedna u odnosu na drugu. Ako su ovi naponi smicanja veći od otpornosti tla dolazi do klizanja između čestica, što uslovljava tečenje materijala. Svesni činjenice da se ni u kom slučaju ne može izbeći idealizacija, usvajamo Hooke-ovu pretpostavku samo do izvesne granice i to do granice tečenja. Pod početkom tečenja podrazumeva se stanje pri kome je granica tečenja upravo dostignuta na jednom ili istovremeno na više mesta u tlu. Proširivanje procesa tečenja kroz konačno prostorno područje nastaje sa daljim povećanjem opterećenja, pri čemu prema usvojenoj idealizaciji deformacije postaju plastične i dobijaju veliki priraštaj. Ako je temelj preopterećen, ili ako tlo ima lošije osobine nego što je pretpostavljeno, tada se Hooke-ov zakon ne može primeniti u većim područjima jer je tlo isprožimano ravnima klizanja. Od izuzetnog značaja je proučavanje pojava koje predhode lomu tla, kao i praćenje tipa deformacije koje se javljaju u samoj fazi loma tla. Pored konstitutivnih veza kojima se opisuje ponašanje materijala u domenu elastičnosti, odnosno do pojave plastičnih deformacija, za opisivanje elasto-plastičnog ponašanja materijala neophodno je poznavati: uslov tečenja, zakon tečenja i zakon ojačanja. Zavisnost ( )σ = σ ε predstavlja uslov koji mora zadovoljiti napon pritiska ili zatezanja da bi došlo do pojave, odnosno povećanja plastične deformacije; ovaj uslov je poznat kao uslov tečenja. On određuje stanje napona pri kome počinje plastično tečenje i definiše površinu u naponskom prostoru koja odvaja oblast elastičnog ponašanja od oblasti u kojoj naponi pored elastičnih izazivaju i plastične deformacije. Do sada je bilo reči o samo jednoosnom naprezanju pa se može konstatovati da se u toj oblasti relativno jednostavno može utvrditi uslov tečenja materijala. Međutim u opštem slučaju, troosnog naprezanja problem je veoma složen. Do sada su vršena mnogobrojna eksperimentalna i teorijska istraživanja, ali ona nisu rezultirala stvaranjem opšte teorije uslova tečenja, koji zapravo predstavlja jedno od centralnih pitanja teorije plastičnosti. Ako sa q obeležimo unutrašnje promenjive vezane za strukturu materijala na makroskopskom nivou, koje se opisuju tenzorom drugog reda pijε  i skalarom k , a Košijev napon σ (spoljašnja promenljiva), onda se kriterijum tečenja može izraziti matematički sa skalarnom funkcijom tečenja f : 22 f ( ,q) 0σ = (2.75) gde je: p ijq q  = ε  (2.76) Veličina q zavisi od komponenata plastične deformacije pijε i načina deformisanja materijala (istorije deformacije) u posmatranoj tački. Pre dalje analize funkcije tečenja f , objasnimo ukratko pojam istorije deformacije materijala u posmatranoj tački. Naime, eksperimentima se može utvrditi da će povećanje komponenti plastične deformacije  pijε , pri određenim vrednostima pijε nastupiti u opštem slučaju pri različitim veličinama  ijσ za dva različita načina deformisanja materijala u posmatranoj tački. Ovo znači da funkcija q ne zavisi samo od vrednosti  pijε već se mora uvesti i parametar koji uzima u obzir način deformisanja. Prema Hill-u [ ]52 kao parametar se može koristiti rad plastičniog deformisanja po jedinici zapremine tela u posmatranoj tački ili ekvivalentna plastična deformacija. Neke od osobina funkcije ijf ( ,q)σ najbolje se mogu utvrditi geometrijskim predstavljanjem ove funkcije u naponskom prostoru. Kako svakoj plastičnoj deformaciji predhodi elastična, funkcija f predstavlja krivu tečenja u naponskom prostoru koja okružuje koordinatni početak ( ij 0σ = ) i elastičnu oblast, kako je prikazano na slici 2.7.a. Na istoj slici je sa 1 označena kriva inicijalne plastične deformacije, koja odgovara tački 1 u dijagramu ( ),σ ε pri jednoosnom naprezanju. Funkcija tečenja f se menja u opštem slučaju po obliku i veličini, usled ojačanja materijala. U literaturi se funkcija f zove i funkcija opterećenja jer se njome mogu utvrditi i uslovi povećanja plastične deformacije (povećanje "opterećenja") u posmatranoj tački materijala, odnosno uslovi pri kojima do toga neće doći. Ako posmatramo promenu napona ijσ u nekoj tački, menjaće se funkcija f , tako da se diferenciranjem po vremenu jednačine (2.75) dobija sledeći izraz: p ij ijp ij ij f f ff k k ∂ ∂ ∂= ⋅σ + ⋅ε + ⋅∂σ ∂ε ∂   (2.77) 23 Ukoliko je u nekom trenutku f 0 ≺ što znači da imamo "rasterećenje", na slici 2.7.a to bi odgovaralo putanji 2-3 u naponskom prostoru, odnosno liniji 2-3 u dijagramu ( ),σ ε , pa nema povećanja plastične deformacije i ojačanja, na osnovu toga sledi da je: ij ij f 0∂ ⋅σ∂σ  ≺ (2.78) kriterijum "rasterećenja". "Neutralno" opterećenje u nekoj tački postoji u slučaju kada pri promeni naponskog stanja ne dolazi do povećanja plastične deformacije ( )pij 0, k 0ε = = i pri tom je f 0= , odnosno: ij ij f 0∂ ⋅σ =∂σ  (2.79) "Neutralno" opterećenje odgovara putanji 2-4 čiji su elementarni delovi upravni na gradijent ijf∂ ∂σ površi f . Pri jednoosnom naprezanju „neutralno“ opterećenje se ne može definisati. Konačno ukoliko je ispunjen uslov f 0 ; , odnosno ij ij f 0∂ ⋅σ∂σ  ; (2.80) imamo povećanje "opterećenja" u posmatranoj tački (putanje 1-2 na slikama 2.7.a i 2.7.b). Zbog ojačanja materijala, kako je već predhodno rečeno , menja se u opštem slučaju i veličina i oblik površi f 0= . Može se pokazati, na osnovu fizičkog uslova da se pri plastičnoj deformaciji javlja gubitak mehaničke energije, da je površ f 0= konveksna, posmatrano iz vrha pozitivne normale-usmerene u smeru povećanja f odnosno u smeru gradijenta ijf∂ ∂σ . Za shvatanje pojma plastičnog materijala potrebno je imati na umu da se, za razliku od elastičnog tela, rad spoljašnjih sila utrošen na deformaciju ne može vratiti. Drugim rečima, rad koji se troši pri ciklusu opterećenje-rasterećenje ne može biti negativan (Drucker-ov postulat). Ovaj postulat uzima se često kao polazna tačka kod izvođenja veza između napona i brzina deformacija. Energija koja se troši pri plastičnoj deformaciji svedena na jedinicu zapremine i vremena da je izrazom p ij ijdR d 0= σ ⋅ ε ≥ (2.81) 24 pri čemu znak jednakosti važi za idealno plastičan materijal. Disipacija (rasipanje) energije data ovim izrazom može se shvatiti kao fizički zakon plastične deformacije. Može se pokazati da je tada, u skladu sa ovim, ispunjen i uslov tečenja ij ij ff 0∂= ⋅σ =∂σ   (2.82) Gornja jednačina ima sledeću geometrijsku interpretaciju: ako u naponskom prostoru shvatimo komponente ijσ kao komponente nekog vektora u devetodimenzijalnom prostoru, vektor ijf∂ ∂σ i vektor brzine ijσ napona su međusobno ortogonalni vektori. Isto tako, i uslov tečenja ima geometrijsku interpretaciju kao hiperpovršina. Vektor ijσ mora biti tangencijalan na tu površinu jer je moguća samo ona promena u naponima koja zadovoljava uslov tečenja - drugim rečima, priraštaj vektora napona mora ležati u tangencijalnoj ravni. Ako je to tako onda je vektor ijf∂ ∂σ normalan na tu površinu. Dodajmo sada uslovu (2.81) još i postulat da priraštaj napona ne vrši nikakav rad na brzinama komponentalnih deformacija plastične deformacije. Ovaj postupak geometrijski interpretirano znači da je vektor pijε normalan na vektor ijσ i shodno ranijoj činjenici sledi da su vektori pijε i ijf∂ ∂σ koaksijalni: p ij ij f∂ε = λ ⋅ ∂σ (2.83) gde je λ skalarni koeficijent. Da je to pozitivan koeficijent možemo videti iz izraza (2.81) ij ij fR 0∂= λ ⋅ ⋅σ∂σ   ; (2.84) Kako je ijσ vektor koji se meri od početka koordinatnog sistema, a ijf∂ ∂σ ima pravac normale na površ koja je konveksna, mora da je : ij ij f 0∂ ⋅σ∂σ  ; (2.85) pa je λ pozitivno i zavisi od naponskog stanja, priraštaja napona, veličine plastične deformacije i istorije deformacije materijala u posmatranoj tački. Treba primetiti da za sve materijale (na primer tlo) uvedeni fizički postulati nisu sasvim u skladu sa ponašanjem tih materijala, pa se smatra da ovo izjednačavanje nije sasvim opravdano. U slučaju da je uslov tečenja definisan sa više funkcija 1f , 2f , ... nf tada je svaka od tih funkcija važeća za jednu određenu oblast naponskog prostora. Geometrijski 25 predstavljeno, pojedini delovi ove hiperpovršine se seku i na tim prodorima vektor nije jednoznačno određen. Tada se može prikazati kao linearna kombinacija vektora koji su normalni na pojedine delove površine: p k ij k ij f , k 1, 2,...n∂ε = λ ⋅ =∂σ (2.86) Pri tome je svaka od pojedinačnih površina tečenja konveksna. Ovde neće biti izlagano razmatranje singularnih tačaka, uslova tečenja sa prekidnim izvodima, a u vezi sa tim i zakona tečenja (videti npr. [ ]15 , [ ]25 , [ ]57 , [ ]108 ). Za uspostavljanje veze između priraštaja plastičnih deformacija i priraštaja napona koristi se pretpostavka prema kojoj je inkrement plastičnih deformacija direktno proporcionalan gradijentu napona. Na taj način, opšti zakon tečenja može da se izvede iz plastičnog potencijala g , koji je između ostalog skalarna funkcija napona. Tada se inkrementi plastičnih deformacija dobijaju pomoću gradijenata plastičnog potencijala. p ij ij g∂ε = λ ⋅ ∂σ (2.87) gde je dλ koeficijent proporcionalnosti. Jednačinom (2.87) definisan je zakon tečenja. Zakonom tečenja uspostavlja se veza između pravaca vektora plastičnih deformacija i funkcije plastičnog potencijala. Obično se uvodi princip normalnosti koji kaže da je inkrement plastične deformacije upravan na površ plastičnog potencijala, odnosno da leži duž pravca spoljašnje normale za ovu površ. U slučaju asocijativne plastičnosti funkcija plastičnog potencijala usvaja se u istom obliku kao funkcija tečenja, pa je inkrement plastične deformacije upravan na površ tečenja. Zakon ojačanja poredi veličine inkrementa plastičnih deformacija sa veličinama inkremenata napona, pošto se stanje napona nalazi na površini tečenja. Ovim zakonom definiše se širenje ili skupljanje površi tečenja u zavisnosti od ukupnih plastičnih deformacija. Do sada su vršena mnogobrojna eksperimentalna i teorijska istraživanja, ali ona nisu rezultirala stvaranjem opšte teorije uslova tečenja, koji zapravo predstavlja jedno od centralnih pitanja teorije plastičnosti. Interesantno je napomenuti da eksperimentalni rezultati pokazuju da se kod frikcionih materijala mogu primeniti i neasocijativni modeli, u kojima se kao plastični potencijal ne koristi funkcija tečenja, već druga funkcija napona. Na teškoće u primeni asocijativne plastičnosti u rešavanju problema mehanike tla ukazao je Drucker [ ]36 , [ ]38 . Na osnovu rezultata velikog broja sprovedenih numeričkih analiza pojedini autori (Zienkiewicz [ ]148 , Vermer [ ]163 ) zaključuju da primena asocijativne plastičnosti 26 u modeliranju mehaničkog ponašanja frikcionih materijala ne daje rezultate koji odgovaraju ponašanju tog materijala u fizičkim eksperimentima. Prema njihovom mišljenju glavni razlog za navedeno odstupanje je pre svega u značajnom povećanju zapreminskih plastičnih deformacija, usled kojih primena asocijativne plastičnosti ne daje realne rezultate pri rešavanju određenih graničnih problema. Ovakvo ponašanje tla prevashodno se nalazi u pojavi dilatancije. Pojava povećanja zapremine tla pri promeni smičućih napona, u mehanici tla se koristi izraz dilatancija, kojim se objašnjava uticaj uzglobljenja zrna na ukupnu smičuću otpornost. Dilatancija je skalarni pokazatelj strukturnih karakteristika određene vrste tla. Dilatanciju zbijenih peskova prvi je uočio Reynolds, 1885.god., koji je eksperimentalno pokazao da se zbijeni peskovi šire pri uticaju smičućih napona. Iz ovih eksperimenata nije proizišla nikakva praktična mera koja bi mogla da ukaže na vezu između dilatancije i smičuće otpornosti zrnastih materijala. Casagrande 1936.god., je pokazao zavisnost ugla smičuće otpornosti peskova u funkciji koeficijenta poroznosti i zapreminskih promena uzorka peska izloženog smicanju. Pokazao je da se zbijeni peskovi šire pri smicanju i pokazuju relativno velike vrednosti ugla smičuće otpornosti, dok se rastresiti peskovi tokom smicanja sabijaju pokazujući pri tome značajno manje vrednosti ugla smičuće otpornosti. Efekti dilatancije su značajni u nižim područjima normalnih napona i postepeno nestaju sa porastom normalnih napona. Pri veoma visokim naponima, ponašanje većine frikcionih materijala zavisi od trenja koje je praćeno drobljenjem zrna. Za dovoljno visok nivo normalnih napona, efekti dilatancije se potpuno gube i materijal se plastično deformiše pri konstantnoj zapremini. Radi definisanja kriterijuma koji bi jedinstveno kvantifikovao pojavu dilatancije u ovom radu se koristi parametar koji se, zbog praktične primene, najčešće koristi u literaturi i predstavlja tzv. ugao dilatancije 27 2.6 Elasto-plastični modeli za tlo 2.6.1 Mohr-Coulomb-ov uslov tečenja Francuski vojni inženjer Coulomb (1773) je postavio hipotezu kojom definiše veličinu otpornosti tla na smicanje, pomoću poznatog izraza: f tan cτ = σ φ+ (2.88) gde su: σ - normalni napon koji deluje u ravni napona fτ c - specifična kohezija tanφ - konstanta proporcionalnosti koja izražava linearnu zavisnost između normalnog napona σ i smičućeg napona fτ . Sama formulacija izraza je posledica eksperimenata u kojima su smicani uzorci drveta i u kojima je pokazano da smičuća otpornost veze dva drvena elementa linearno zavisi od veličine normalnog napona. Njegovim radovima započelo je tumačenje graničnih stanja u tlu. Uglavnom je razmatrao problem pritisaka na potporne konstrukcije, ali je analizirao i stabilnost bočnih strana iskopa u prirodno vlažnom koherentnom tlu, pri čemu je koristeći pretpostavku o ravnim površinama klizanja, odredio izraz za graničnu visinu vertikalnog zaseka. Mohr (1882) je formulisao uslov loma za materijale sa nejednakom otpornošću na pritisak i zatezanje, izražen preko glavnih napona. Ograničio je domen važnosti ovako formulisanog uslova loma, tako da se ova linearna anvelopa posmatra kao granična kriva za naponska stanja pri jedoosnom pritisku i zatezanju. Stiče se utisak da Mohr i nije imao nameru da uspostavi opšti uslov loma, već je pokušao da odredi međusobni odnos osnovnih mehaničkih karakteristika, i to pre svega između čvrstoće na smicanje i čvstoće na pritisak. Mohr-ova hipoteza ne uzima uticaj srednjeg glavnog napona na formiranje uslova plastičnosti, tako da na granično stanje utiču samo algebarski najveći i najmanji napon. Mohr je relaciju (2.88) izrazio preko glavnih napona u sledećem obliku: ( ) ( )1 3 1 3 tan 2d 0σ −σ − σ +σ ⋅ α − = (2.89) gde je: d c cos , tan sin= ⋅ φ α = φ (2.90) 28 Empirijski Mohr-Coulomb-ov linearni kriterijum tečenja zadovoljavajuće opisuje stanje napona pri tečenju za određeni interval promene napona za sve vrste tla u konvencionalnoj mehanici tla. Presek ovako definisane prave linije sa nultim nivoom normalnih napona daje odsečak koji predstavlja komponentu smičuće čvrstoće koja ne obuhvata trenje. Ovaj parametar ima fizičkog smisla samo onda kad su čestice tla međusobno povezane značajnim silama, bez obzira na njihovu pravu prirodu; tada tlo ima koheziju. Parametre smičuće otpornosti treba shvatiti kao empirijske veličine i oni zavise od pripreme samih uzoraka i vrste ispitivanja. Potrebno je primeniti ona ispitivanja koja najpribližnije simuliraju stanje napona i deformacija i njihove promene u posmatranom uzorku tla. Linearna idealizacija, kod koje su kohezija i ugao smičuće čvrstoće nezavisni od nivoa napona, može se smatrati dobrom za rastresite peskove i normalne konsolidovane gline. Kod materijala koji pokazuju kruto ponašanje, tečenje se dešava progresivno, sa različitim stepenom angažovanja otpornosti u pojedinim tačkama same površi tečenja. Eksperimentalno je utvrđeno da se gore navedena zakonitost može primeniti sa konstantnim vrednostima parametra otpornosti samo za određene promene normalnog napona. Odstupanje od linearne anvelope može biti izraženo u području malih napona, što je posebno važno kod ispitivanja pri kojima su naponi niski. Izraze (2.62) i (2.64) možemo napisati u sledećem obliku: 1 3 2dJ cos2 σ −σ = θ (2.91) 29 1 3 1 2d 1 1I J sin 2 3 3 σ +σ = − θ (2.92) Ako u jednačinu (2.89) unesemo izraze (2.91) i (2.92), jednačinu (2.90) možemo, posle elementarnih transformacija, napisati u sledećem obliku: 1 2d 2d 1 If J cos J sin sin sin c cos 0 33 = θ + θ φ− φ− ⋅ φ = (2.93.1) gde je Lode-ov ugao: 1 3d 3/ 2 2d I1 3 3sin 3 2 I −  θ = − ⋅ ⋅    (2.93.2) Ako se uzmu u obzir i druge dve kombinacije napona ( 2 1σ σ; i 3 1σ σ; ) mogu se dobiti još dve relacije oblika (2.93). U trodimenzionalnom prostoru glavnih napona površ tečenja, koja se dobija kada se uzmu u obzir sve tri predhodne kombinacije, predstavlja nepravilnu šestostranu piramidu duž hidrostatičke ose čija je projekcija na ravan nepravilni šestougaonik (slika 2.9). Ovakav oblik površi tečenja koji se u ovom slučaju javlja, a ne cilindrični, posledica je učešća sfernog dela tenzora napona u uslovu tečenja. Za vezano tlo, ako uticaj kohezije zamenimo uporednim normalnim naponom cσ , jednačinu (2.93.1) za slučaj hidrostatičkog naponskog stanja ( 1 2 3σ = σ = σ ), možemo napisati u sledećem obliku 1 c c ctgσ = σ = ⋅ φ . Geometrijskom interpretacijom ovog izraza uočava se položaj vrha piramide na prostornoj dijagonali ( 1 2 3σ = σ = σ ) a na odstojanju od koordinatnog početka koje je jednako uporednom normalnom naponu c c ctgσ = ⋅ φ . 30 Vrednosti maksimalnog i minimalnog poluprečnika krugova koji prolaze kroz uglove piramide u π ravni jednake su: max 2d 2 6 с cos, R 2J 6 3 sin π ⋅ ⋅ φθ = − = = − φ (2.94.1) min 2d 2 6 с cos, R 2J 6 3 sin π ⋅ ⋅ φθ = = = + φ (2.94.2) 2.6.2 Drucker-Prager-ov uslov tečenja Ovaj model je nastao aproksimacijom Mohr-Coulomb-ove površi tečenja sa glatkom površi i kao modifikacija dobro poznatog Von Mises-ovog materijalnog modela uključujući sfernu komponentu tenzora napona. Drucker i Prager (1952) predložili su uslov tečenja u obliku: 2d 1f I I k 0= −α − = (2.95) gde su: 1I - prva invarijanta tenzora napona, 2dI - druga invarijanta devijatora tenzora napona, , kα - materijalne konstante. Ukoliko je 0α = jednačina (2.95) prelazi u Von Mises-ov uslov tečenja. Konstante α i k mogu se odrediti u zavisnosti od kohezije c i ugla smičuće otpornosti φ predpostavljajući da istovremeno Drucker-Prager-ov i Mohr-Coulomb-ov kriterijum tečenja daju identično granično opterećenje za date eksperimentalne uslove. Koristeći jednačine (2.94.1) i (2.95) zavisnost se može predstaviti, u triaksijalnim uslovima za CTC test: ( ) 2sin 3 3 sin φα = − φ (2.96.1) ( ) 6ccos ck tan3 3 sin φ= = αφ− φ (2.96.2) 31 U uslovima ravanskog stanja deformacija ova veza se može prikazati u obliku: 2 tan 9 12 tan φα = + φ (2.97.1) 2 3ck 9 12 tan = + φ (2.97.2) Ograničenja u praktičnim primenama Drucker-Prager-ovog modela su detaljno prikazana u radu [ ]21 . 2.6.3 Kapa model DiMaggio i Sandler (1971) su formulisali elasto-plastičan model koji se sastoji od dve površi tečenja koje se seku. Jedna površ tečenja opisuje idealno elasto-plastično ponašanje, a druga površ tečenja - kapa obuhvata izotropno ojačanje materijala pri povećanju plastične deformacije. Familija površi tečenja formira "kape" [ ]21 , pa otuda i ime modelu. Modifikovanom Drucker-Pragerovom površi tečenja je fiksirana u prostoru glavnih napona i može se napisati u obliku 1I f 2d 1f I e I 0 −β= + γ −θ −α = (2.98) gde je: 1I - prva invarijanta tenzora napona, 2dI - druga invarijanta devijatora tenzora napona i , , ,α β γ θ - materijalni parametri koji se eksperimentalno određuju. Druga površ tečenja je kapa koja se širi u prostoru glavnih napona u funkciji zapreminskih plastičnih deformacija. Posmatramo slučaj kape prikazane sledećom jednačinom: ( )22 2 2c 2d 1f R I I C R B 0= + − − = (2.99) koja predstavlja elipsu sa poluosama X C− i B . Parametar R predstavlja odnos između manje i veće poluose elipse, dok su C i X centar odnosno teme elipse na apcisnoj osi. Ojačanje materijala je definisano preko promene zapreminske plastične deformacije, u obliku: p v1X ln 1 Z D W  ε= − − +   (2.100) gde su: 32 X - vrednost 1I u preseku kape i hidrostatičke ose, p vε - zapreminska plastična deformacija koja odgovara tekućem X i D , Z i W - materijalni parametri modela. Veličina D određuje brzinu promene zapreminske plastične deformacije, Z određuje početnu elipsu, a W predstavlja maksimalnu plastičnu zapreminsku deformaciju koja smanjuje poroznost materijala na početnu vrednost. U elastičnoj oblasti model obuhvata ponašanje nelinearno elastičnog tela. U elasto- plastičnoj oblasti prirodu deformacija određuje eliptična kapa. Ako se tekuće naponsko stanje nalazi na površi kape u slučaju asocijativne plastičnosti vektor priraštaja plastičnih deformacija je upravan na eliptičnu površ tečenja. Zbog oblika kape, vektor priraštaja zapreminske plastične deformacije ima uvek pozitivnu komponentu u pravcu prve invarijante, pri čemu se kapa pomera unapred, odnosno materijal ojačava. U slučaju plastičnog deformisanja na graničnoj površi loma priraštaji zapreminskih plastičnih deformacija su negativni, odnosno sa razvojem dilatancije javlja se kontrakcija kape. Na taj način je kretanje kape određeno povećanjem ili smanjenjem zapreminske plastične deformacije. Za domen koji obuhvata male napone se može uvesti i tzv. početna (inicijalna) kapa čiji je centar u koordinatnom početku. Za definisanje modela potrebno je odrediti deset materijalnih parametara: dva elastična ( E,ν ), četiri za graničnu površ tečenja ( , , ,α β γ θ ) i na kraju još četiri konstante za eliptičnu kapu ( D , Z , W i R ). 33 Potrebno je napomenuti da se kod nelinearne elastičnosti broj elastičnih parametara povećava srazmerno usvojenom stepenu nelinearnosti. Svi materijalni parametri se mogu odrediti na osnovu standardnih laboratorijskih testova i numeričkih procedura [ ]31 . 2.6.4 HISS model Desai [ ]32 , [ ]33 je predložio hijerarhijski koncept definisanja materijalnih modela (HISS) u oblasti plastičnosti. Opšti oblik funkcije tečenja predstavljen je polinomom koji sadrži invarijante napona i devijatora napona. Funkcija tečenja menja oblik i veličinu u toku deformisanja, čime se obuhvata ojačanje ili omekšanje materijala. U ovom radu razmatra se funkcija tečenja za izotropno ojačanje asocijativne plastičnosti: 2 a 2d b sf p I F F 0 −= − ⋅ = (2.101.1) gde je: ( ) ( )n 2b 1 a 1 aF I p I p= −α ⋅ + γ ⋅ (2.101.2) ( )ms rF 1 S= −β⋅ (2.101.3) 3 2 r 3d 2dS 3 3 2 I I −= ⋅ ⋅ (2.101.4) α funkcija ojačanja, n , γ , β i m materijalne konstante i ap atmosferski pritisak. bF je osnovna funkcija koja prikazuje oblik funkcije tečenja u 1 2dI I− prostoru, sF funkcija koja reprezentuje oblik u oktaedarskoj ravni. Oblik funkcije ojačanja se definiše relacijom: 1a ηα = ξ (2.102) gde su: 1a i η materialne konstante kojima se opisuje ojačanje materijala i ξ je ekvivalentna (efektivna) plastična deformacija. Efektivna plastična deformacija definisana je relacijom: 34 ( )1/ 2p pij ijd dξ = ε ε∫ (2.103) Slika 2.11 Funkcija tečenja HISS (C.S.Desai, 1986): (a) u prostoru 1 2dI I− ; (b) u triaksijalnoj ravni; (c) u devijatorskoj ravni HISS uslov tečenja uspešno opisuje ponašanje frikcionih materijala za različite putanje napona. Prednost ovog modela, između ostalog, je što za razliku od većine kompleksnih elasto-plastičnih modela nema posebno definisanu graničnu površ loma. 35 3. INKREMENTALNA FORMULACIJA OSNOVNIH JEDNAČINA STATIČKE ANALIZE 3.1. Uvod Kada je potrebno formirati matematički model posmatranog mehaničkog sistema u mehanici tla ne mogu se izbeći nelinearni problemi, što podrazumeva korišćenje formulacija nelinearne mehanike kontinuuma. Nelinearna mehanika kontinuuma nije dostigla stepen razvoja u kome se većina problema mogu svesti na diferencijalne jednačine za koje postoji razrađene metode integracije, za razliku od teorije elastičnosti koja je dostigla izvestan stepen matematičke potpunosti [ ]89 . Odeljci 3.1, 3.2 i 3.3 ovog poglavlja urađeni su prema mongrafiji[ ]8 . Na slici 3.1 prikazane su tri konfiguracije proizvoljnog tela u toku deformacije: početna 0Ω , tekuća nΩ i naknadna konfiguracija n 1+ Ω . Predpostavljamo da je rešenje za kinematičke i statičke veličine za sve konfiguracije 0Ω , 1Ω , 2Ω , ..., nΩ poznato i da je potrebno odrediti rešenje za konfiguraciju n 1+ Ω . Koordinate proizvoljne tačke P tela u konfiguraciji 0Ω su 0 01 2x , x i 0 3x u nΩ su n n 1 2x , x i n 3x i u n 1+ Ω su n 1 n 11 2x , x+ + i n 1 3x+ . Na sličan način se obeležavaju i pomeranja, odnosno: ( )n 0 ni i in 1 0 n 1 i i i x x u , i 1, 2,3 x x u+ + = +  == +  (3.1) 36 Nepoznate inkremente pomeranja od nΩ do konfiguracije n 1+ Ω mogu se obeležiti sa ( )n 1 ni i iu u u , i 1,2,3+∆ = − = (3.2) U toku deformacije tela, njegova površina, naponi i deformacije menjaju se kontinualno. Površinu tela ćemo obeležiti sa 0 A , n A i n 1A+ . Pošto je konfiguracija n 1+ Ω nepoznata, sve spoljašnje sile, napone i deformacije izrazićemo u odnosu na neku poznatu ravnotežnu konfiguraciju. Komponente površinskih sila po jedinici površine u konfiguraciji n 1+ Ω , ali merene u konfiguraciji nΩ , su ( )n 1n ip , i 1, 2,3+ = . Slično je i kod napona. Komponente Piola-Kirchoff-ovog tenzora napona druge vrste koje odgovaraju konfiguraciji n 1+ Ω , ali merene u konfiguraciji nΩ , obeležićemo sa n 1n ijS+ . Ako se radi o veličini koja se posmatra u istoj konfiguraciji u kojoj se meri tada se levi donji indeks može izostaviti. To je slučaj kod Cauchy-jevog tenzora napona, tj. n 1 n 1n + +σ = σ . Komponente Cauchy-jevog infinitezimalnog tenzora deformacije koje odgovaraju konfiguraciji n 1+ Ω obeležićemo sa n 1 ije + . Komponente Green-Lagrange-ovog tenzora deformacije koje odgovaraju konfiguraciji n 1+ Ω u odnosu na konfiguraciju nΩ obeležićemo sa n 1n+ ε . Referentne konfiguracije u odnosu na koje ćemo određivati veličine spoljnih sila Piola- Kirchoff-ovog tenzora napona druge vrste i Green-Lagrange-ovog tenzora deformacije će biti 0Ω i nΩ . U formulaciji osnovnih jednačina ravnoteže pojaviće se i izvod pomeranja i koordinata. Prema usvojenoj notaciji zarezom ćemo označiti izvod po koordinati koja sledi iza nje, a levi donji indeks će označavati konfiguraciju u kojoj je posmatrana koordinata merena. U Lagrange-ovoj inkrementalnoj formulaciji ravnotežu tela u konfiguraciji (n+1) izrazićemo korišćenjem principa virtuelnih pomeranja. Osnovna pretpostavka od koje se polazi pri definisanju ovoga principa je da on važi za konstrukciju u ravnoteži. Iskaz da je algebarski zbir radova svih spoljašnjih i unutrašnjih sila pri bilo kojim virtuelnim pomeranjima tačaka tela jednak je nuli formuliše princip virtuelnih pomeranja. Pod pojmom virtuelna ili moguća pomeranja podrazumevamo beskonačno mala proizvoljna pomeranja koja su neprekidne i diferencijabilne funkcije koordinata tačaka tela i jednaka su nuli u svim onim tačkama u kojima su konturni uslovi zadati po pomeranjima. Princip virtuelnih pomeranja, odnosno slaba forma jednačine momentnog balansa, zadovoljava diferencijalnu jednačinu problema (uslove ravnoteže tela), granične uslove koji predstavljaju poznata 37 pomeranja i zadate sile na granici razmatranog domena. Princip virtuelnih pomeranja za konfiguraciju n 1+ Ω glasi: n 1 n 1 u sR R + +δ = δ (3.3) gde je: n 1 uR + δ - virtuelni rad unutrašnjih sila (rad stvarnih napona na virtuelnim deformacijama), a n 1 sR + δ - virtuelni rad spoljašnjih sila (rad zapreminskih i površinskih sila na virtuelnim pomeranjima), n 1 n 1 n 1 n 1 n 1 u ij ij V R e dV + + + + +δ = σ δ∫ (3.4.1) n 1 n 1 n 1 n 1 n 1 n 1 n 1 n 1 n 1 s n 1 s n 1 i i n 1 i i V A R f u dV p u dA + + + + + + + + + + + +δ = ρ δ + δ∫ ∫ (3.4.2) ( )n 1n 1 jn 1 n 1 n 1iij n 1 i, j n 1 j,in 1 n 1 j i u1 u 1e u u 2 x x 2 ++ + + + + ++ +  ∂δ∂δδ = + = δ +  ∂ ∂  (3.4.3) Primenom jednačine (3.3) dolazi se do jednačina ravnoteže. Rešenje jednačine (3.3) nije moguće dobiti direktno pošto je konfiguracija n 1+ Ω nepoznata. Moguće je dobiti približno rešenje ako se sve promenljive izraze preko odgovarajućih veličina u nekoj predhodno izračunatoj ravnotežnoj konfiguraciji i izvrši linearizacija tako dobijene jednačine. Bilo koja već sračunata konfiguracija može se usvojiti kao referentna konfiguracija. U praktičnim primenama koriste se dve varijante Lagrange-ove inkrementalne formulacije: totalna Lagrange-ova formulacija i korigovana Lagrange-ova formulacija. 3.2 Formulisanje rešenja nelinearnih jednačina 3.2.1. Totalna Lagrange-ova formulacija U totalnoj Lagrange-ovoj formulaciji za referentnu konfiguraciju usvajamo početnu konfiguraciju 0Ω . Ako u (3.4.1) izrazimo n 1 uR+ δ preko Piola-Kirchoff-ovog tenzora napona druge vrste i Green-Lagrange-ovog tenzora deformacije, koristeći relaciju (36) ispisanu za konfiguraciju n 1+ Ω : 0 n 1 n 1 n 1 0 u 0 ij 0 ij V R S dV+ + +δ = δ ε∫ (3.5) Spoljne sile u izrazu (3.4.2) možemo napisati u obliku 38 n 1 n 1 n 1 0 n 1 0 n 1 i 0 if dV f dV + + + + +ρ = ρ (3.6) n 1 n 1 n 1 0 n 1 i 0 ip dF p dA + + + + = (3.7) gde je usvojena pretpostavka da veličine i pravci sila 0 n 10 if +ρ i n 10 ip+ ne zavise od deformacije tela u konfiguraciji n 1+ Ω . Jednačine ravnoteže tela u konfiguraciju n 1+ Ω u odnosu na ravnotežnu konfiguraciju 0Ω možemo napisati u obliku: 0 n 1 n 1 0 n 1 0 ij 0 ij s V S dV R+ + +⋅δ ε ⋅ = δ∫ (3.8) gde je na osnovu (3.4.2), (3.6) i (3.7) 0 0 n 1 0 n 1 n 1 0 n 1 n 1 a 0 s 0 i i 0 i i V A R f u dV p u dA+ + + + +δ = ρ δ + δ∫ ∫ (3.9) U jednačini (3.8) n 10 ijS + -Piola-Kirchoff-ov tenzor napona druge vrste u konfiguraciji n 1+ Ω , a merene u konfiguraciji 0Ω je 0 n 1 0 n 1 0 0 ij n 1 i,k kl n 1 j,ln 1S x x + + + ++ ρ= ⋅ σ ⋅ρ (3.10) a n 10 ij +δ ε varijaciona komponenta Green-Lagrange-ov tenzor deformacije u konfiguraciji n 1+ Ω , a merene u konfiguraciji 0Ω Tenzor deformacije koji je u saglasnosti sa Piola-Kirchoff-ovim tenzorom napona druge vrste je Green-Lagrange-ov tenzor deformacije. Varijacija Green-Lagrange-ov tenzor deformacije može se prikazati jednačinom: ( )n 1 n 1 n 1 n 1 n 10 ij 0 i, j 0 j,i 0 k,i 0 k, j1 u u u u2+ + + + +δ ε = δ + + ⋅ (3.11) Kao što je već rečeno, ako su nam poznate komponente Piola-Kirchoff-ovog tenzora napona druge vrste i pomeranja u konfiguraciji nΩ , možemo odgovarajuće veličine vezane za konfiguraciju n 1+ Ω odrediti korišćenjem sledeće inkrementalne dekompozicije: n 1 n 0 ij 0 ij ijS S S + = + ∆ (3.12.1) n 1 n 0 ij 0 ij ij + ε = ε + ∆ε (3.12.2) n 1 n i i iu u u + = + ∆ (3.12.3) 39 gde su: ijS∆ i ij∆ε inkrementi Piola-Kirchoff-ovog tenzora napona druge vrste i Green- Lagrange-ovog tenzora deformacije koji se odnose na početnu konfiguraciju 0Ω , a iu∆ je inkrement pomeranja. Polazeći od relacije kojom je izražen Green-Lagrange-ov tenzor deformacije preko pomeranja, uz vođenje računa o jednačini (3.12), dobija se: n 1 n ij 0 ij 0 ij +∆ε = ε − ε (3.13) ( ) ( ) ( )n 1 n n 1 n n 1 n 1 n nij 0 i, j 0 i, j 0 j,i 0 j,i 0 i,k 0 k, j 0 i,k 0 k, j1 u u u u u u u u2 + + + + ∆ε = − + − + ⋅ − ⋅ =  ( ) ( )n 1 n n 1 ni i j j0 0 j i 1 u u u u 2 x x + + ∂ ∂− + − +∂ ∂ ( ) ( )n n n n0 i,k k,i 0 j,k k, j 0 i,k 0 k, ju u u u u u + ∆ + + ∆ + ⋅  (3.14) Odnosno; ( )n nij 0 i, j 0 j,i 0 i,k 0 k, j 0 j,k 0 k,i 0 i,k 0 k, j1 u u u u u u u u2∆ε = ∆ + ∆ + ⋅ ∆ + ⋅ ∆ + ∆ ⋅ ∆ (3.15) Inkrement Green-Lagrange-ovog tenzora deformacije ij∆ε može se prikazati kao odgovarajući zbir linearnog i nelinearnog dela: ij ij ije∆ε = ∆ + ∆η (3.16) gde su: ( )n nij 0 i, j 0 j,i 0 i,k 0 k, j 0 j,k 0 k,i1e u u u u u u2∆ = ∆ + ∆ + ⋅ ∆ + ⋅ ∆ (3.17) ij 0 i,k 0 k, j 1 u u 2 ∆η = ∆ ⋅ ∆ (3.18) Iz jednačine (3.12.2) sledi: ( )n 1 n0 ij 0 ij ij ij+δ ε = δ ε + ∆ε = δ∆ε (3.19) pa korišćenjem gornje jednačine i jednačine (3.8) uslov ravnoteže možemo napisati u obliku: 0 0 0 n 1 n 0 ij ij s 0 ij ij V V S dV R S dV+∆ ⋅δ∆ε ⋅ = δ − ⋅δ∆ε ⋅∫ ∫ (3.20) veza između inkremenata napona i deformacija je oblika: n ij 0 ijkl klS D∆ = ∆ε (3.21) 40 pa smenom jednačine (3.21) u jednačini (3.20) dobijamo oblik: 0 0 0 n 0 n 0 n 1 n 0 0 ijkl kl ij 0 ij ij s 0 ij ij V V V D dV S dV R S dV+∆ε ⋅δ∆ε ⋅ + ⋅δ∆η ⋅ = δ − ⋅δ∆ε ⋅∫ ∫ ∫ (3.22) koji predstavlja nelinearnu jednačinu za određivanje inkrementa pomeranja iu∆ u totalnoj Lagrange-ovoj formulaciji. 3.2.2.Korigovana Lagrange-ova formulacija Kod korigovane Lagrange-ove formulacije za referntnu konfiguraciju usvajamo tekuću konfiguraciju tela nΩ . Analogno predhodnom razmatranju, jednačina (3.3) dobija oblik: n n 1 n 1 n n 1 n ij n ij s V S dV R+ + +⋅δ ε ⋅ = δ∫ (3.23) n 1 n ijS + -Piola-Kirchoff-ov tenzor napona druge vrste u konfig. n 1+ Ω , a merene u konfig. nΩ n 1 n ij + ε -Green-Lagrange-ov tenzor deformacije u konfiguraciji n 1+ Ω , a merene u konfig. nΩ n 1 sR + δ -rad spoljašnjih sila pri virtuelnim pomeranjima Pretpostavljajući da je opterećenje nezavisno od deformacije tela, slično prethodnim razmatranjima, dobija se nelinearna jednačina za određivanje inkrementalnog pomeranja iu∆ u korigovanoj Lagrange-ovoj formulaciji u sledećem obliku: n n n n n n n n 1 n n n ijkl n kl n ij ij n ij s ij ij V V V D dV dV R dV+∆ε ⋅δ ∆ε ⋅ + σ ⋅δ ∆η ⋅ = δ − σ ⋅δ∆ε ⋅∫ ∫ ∫ (3.24) 3.2.3 Linearizacija jednačina ravnoteže U inkrementalno iterativnim postupcima princip virtuelnih pomeranja je zadovoljen na kraju ili početku svakog inkrementa opterećenja što znači da je primenom ovih postupaka obezbeđena ravnoteža razmatranog mehaničkog sistema (tela) samo u diskretnim koracima nanošenja spoljnog opterećenja. Da bi se postigla zahtevana ravnoteža tela u pojedinim inkrementima razmatranog nelinearnog mehaničkog sistema neophodno je korišćenje odgovarajućih iterativnih postupaka. Poslužićemo se geometrijskom interpretacijom, tj. u prostoru nezavisnih parametara ravnotežna stanja sistema formiraju krivu, pa konvergencija rešenja direktno zavisi od usvojenog reda aproksimacije te krive. Postupkom linearizacije 41 uvodimo aproksimaciju razmatrane krive njenom tangentom. Linearizacijom principa virtuelnih pomeranja dolazi se do pojma neuravnoteženih sila koje deluju na posmatrani mehanički sistem kao i do pojma tangentne matrice krutosti sistema koja karakteriše odgovor sistema na neuravnoteženu pobudu. Jednačine ravnoteže (3.22) i (3.24) su nelinearne po inkrementima pomeranja. Linearizaciju možemo izvršiti usvajanjem sledeće aproksimacije: 0 ij 0 ije∆ε ≈ ∆ (3.25) ij 0 ijkl 0 klS D∆ = ⋅ ∆ε (3.26) ij 0 ijeδ∆ε ≈ δ ∆ (3.27) Ako u jednačini (3.22) unesemo odgovarajuće izraze možemo napisati sledeću jednačinu ravnoteže u totalnoj Lagrange-ovoj formulaciji: 0 0 0 n 0 n 0 n 1 n 0 0 ijkl 0 kl ij 0 ij 0 ij e 0 ij 0 ij V V V D e e dV S dV R S e dV+⋅ ∆ ⋅δ∆ ⋅ + ⋅δ ∆η ⋅ = δ − ⋅δ ∆ ⋅∫ ∫ ∫ (3.28) Prema predhodno izvršenoj linearizaciji: ij n ijkl klS D e∆ ≈ ⋅∆ (3.29) ij ijeδ∆ε ≈ δ∆ (3.30) Jednačina ravnoteže u korigovanoj Lagrange-ovoj formulaciji (3.24) dobija sledeći oblik: n n n n n n n n 0 n 1 n ijkl n kl n ij ij 0 ij ij n ij s V V V D e e dV dV e dV R+⋅ ∆ ⋅δ ∆ ⋅ + σ ⋅δ ∆η ⋅ + σ ⋅δ ∆ ⋅ = δ∫ ∫ ∫ (3.31) Jednačine (3.28) i (3.31) su linearne jednačine po inkrementalnim pomeranjima. Izraz (3.31) predstavlja linearizovani oblik priraštaja unutrašnjeg virtuelnog rada u odnosu na priraštaj prostornih pomeranja u datoj tački prostora. Na osnovu njega se određuje linearizovani priraštaj unutrašnje energije sistema u odnosu na komponente prostornog vektora pomeranja tačaka razmatranog sistema, a u odnosu na trenutnu konfiguraciju. Može se odrediti putanja kojom će se sistem kretati iz tekuće konfiguracije. 42 3.3 Formulacija rešenja metodom konačnih elemenata Za rešavanje ovih jednačina primenićemo metodu konačnih elemenata. Usvajanjem izoparametarske diskretizacije, do odgovarajućih jednačina MKE dolazimo tako što koristimo iste interpolacione funkcije za koordinate i pomeranja unutar jednog konačnog elementa, kao što se to radi i standardnom procedurom u linearnoj analizi, tj.: N N N 0 0 n n n 1 n 1 i k ik i k ik i k ik k 1 k 1 k 1 x h x ; x h x ; x h x+ + = = = = ⋅ = ⋅ = ⋅∑ ∑ ∑ (3.32) N N n n i k ik i k ik k 1 k 1 u h u ; u h u = = = ⋅ ∆ = ⋅∆∑ ∑ (3.33) gde n ikx je koordinata čvora k u pravcu i u konfiguraciji nΩ , n iku je pomeranje čvora k u pravcu i u konfiguraciji nΩ , iku∆ je inkrement pomeranje čvora k u pravcu i u konfiguraciji nΩ , kh je interpolaciona funkcija polja pomeranja unutar elementa konačne zapremine V i površi A koja se odnosi na čvor k, i N je ukupan broj čvorova u elementu. Interpolaciona funkcija kh se definiše u sistemu prirodnih koordinata ( )ir i 1, 2,3= koje variraju u intervalu od -1 do 1. Interpolaciona funkcija kh ima osobinu da je njena vrednost u sistemu prirodnih koordinata jednaka jedinici u čvoru k i jednaka nuli u svim ostalim čvorovima. Da bismo omogućili efikasniju transpoziciju predhodnih izraza u MKE, jednačinu (3.33) za razmatrani konačni element m možemo napisati u matričnom obliku: n n mu H u= ⋅ (3.34.1) mu H u∆ = ⋅∆ (3.34.2) Iz jednačine (3.17) za linearni deo inkrementa Green-Lagrange-ovog tenzora deformacije u totalnoj Lagrange-ovoj formulaciji dobija se: n 0 0 de G u∆ = ⋅ ∆ (3.35) gde n0G matrica gradijenata deformacije u konfiguraciji nΩ a 0 du∆ vektor gradijenata pomeranja. Koristeći izraz (3.34.2) vektor 0 du∆ možemo prikazati pomoću pomeranja čvorova konačnog elementa: 0 d 0 d mu H u∆ = ⋅∆ (3.36) 43 gde je 0 dH matrica izvoda interpolacionih funkcija po koordinatama ( )0 ix i 1,2,3= . Na osnovu izraza (3.35) i (3.36) linearni deo inkrementa Green-Lagrange-ovog tenzora deformacije je: n 0 me B u∆ = ⋅∆ (3.37) gde je: n n 0 0 0 dB G H= ⋅ (3.38) U matrici n0G sadržani su gradijenati pomeranja u konfiguraciji nΩ ( )n0 m,iu koji se mogu odrediti iz izraza (3.34.1) kao: n n 0 d 0 d mu H u= ⋅ (3.39) Za konačni element m može se napisati: ( )T T n T n0 ijkl kl ij 0 m 0 0 0 mD e e e D e u B D B u⋅∆ ⋅δ∆ = δ∆ ⋅ ⋅∆ = δ∆ ⋅ ⋅ ⋅ ⋅∆ (3.40) odnosno: ( ) 0 0 T n 0 ijkl kl ij m 0 L m V D e e dV u K u⋅∆ ⋅δ∆ ⋅ = δ∆ ⋅ ⋅∆∫ (3.41) gde je n0 LK linearna matrica krutosti za element m koja je data izrazom: ( ) 0 n n T n 0 0 L 0 0 0 V K B D B dV= ⋅ ⋅ ⋅∫ (3.42) gde je 0 D inkrementalna konstitutivna matrica u konfiguraciji nΩ definisana u odnosu na u konfiguraciji 0Ω . Analognim razmatranjem, imajući u vidu jednačinu (3.18), može se izraz za geometrijsku matricu krutosti elementa m napisati u obliku: 0 n T n 0 0 G 0 d 0 0 d V K H S H dV= ⋅ ⋅ ⋅∫ (3.43) Za razmatrani konačni element m može se izraziti i integral: 0 0 0 n 0 T n 0 T n T n 0 T n 0 ij ij 0 m 0 0 m 0 V V V S e dV e S dV u B S dV u F⋅δ∆ ⋅ = δ∆ ⋅ = δ∆ ⋅ ⋅ ⋅ = δ∆ ⋅∫ ∫ ∫ (3.44) gde je: n0 F vektor čvornih sila koje su ekvivalentne naponima u elementu m u konfig. nΩ , tj.: 44 0 n n T n 0 0 0 0 V F B S dV= ⋅ ⋅∫ (3.45) a n0S vektor Piola-Kirchoff-ovog tenzora napona druge vrste za element m u konfig. nΩ . Izraz n 1 sR + δ za virtuelni rad spoljašnjih sila (rad zapreminskih i površinskih sila na virtuelnim pomeranjima), 0 0 n 1 0 n 1 n 1 0 n 1 n 1 a 0 s 0 i i 0 i i V A R f u dV p u dA+ + + + +δ = ρ δ + δ∫ ∫ (3.46) Kako je prema (3.12.3) ( )n 1 ni i i iu u u u+δ = δ + ∆ = δ∆ dobija se: 0 0 n 1 0 n 1 0 n 1 a 0 s 0 i i 0 i i V A R f u dV p u dA+ + +δ = ρ⋅ ⋅δ∆ ⋅ + ⋅δ∆ ⋅∫ ∫ (3.47) Matrični oblik ove jednačine, za razmatrani konačni element m, je: 0 0 Tn 1 0 T n 1 0 a n 1 0 s m 0 m 0 V A R u f dV u p dA+ + +δ = ρ⋅δ∆ ⋅ ⋅ + δ∆ ⋅ ⋅∫ ∫ (3.48) ili korišćenjem izraza (3.34.2) n 1 T n 1 s mR u P + +δ = δ∆ ⋅ (3.49) gde je n 1P+ vektor spoljnih sila koji odgovaraju elementu m u konfiguraciji n 1+ Ω tj. 0 0 n 1 0 T n 1 0 aT n 1 0 0 0 V A P H f dV H p dA+ + += ρ⋅ ⋅ ⋅ + ⋅ ⋅∫ ∫ (3.50) Unošenjem jednačina (3.41), (3.45) i (3.50) u jednačinu ravnoteže (3.28) u totalnoj Lagrange- ovoj formulaciji, za element m dobija se: ( )n n n 1 n0 L 0 G m 0K K u P F++ ⋅∆ = − (3.51) Slično totalnoj Lagrange-ovoj formulaciji možemo formulisati i rešenje jednačina (3.31) u korigovanoj Lagrange-ovoj formulaciji: ( )n n n 1 nn L n G m nK K u P F++ ⋅∆ = − (3.52) gde su: ( ) n n n T n n n L n n n V K B D B dV= ⋅ ⋅ ⋅∫ (3.53) 45 n n T n n n G n d n n d V K H S H dV= ⋅ ⋅ ⋅∫ (3.54) n n n T n n n n n V F B S dV= ⋅ ⋅∫ (3.56) dok za n 1P+ vektor spoljnih sila koji odgovaraju elementu m u konfiguraciji n 1+ Ω važi izraz (3.50). Elementi matrica u izrazima (3.42) do (3.56) su funkcije prirodnih koordinata elemenata i da se odgovarajući integrali najčešće izračunavaju numeričkim putem koristeći se transformacijom Descartes-ovih koordinata u prirodne koordinate. Integral prikazan izrazom (3.42) može se napisati u obliku: ( ) ( ) 0 0 n n T n 0 0 L 0 0 0 1 2 3 1 2 3 V V K B D B dV r , r , r dr dr dr= ⋅ ⋅ ⋅ = χ ⋅ ⋅ ⋅∫ ∫ (3.57) gde je: n T n 0 0 0 0B D B det( J)χ = ⋅ ⋅ ⋅ (3.58) a 0 J matrica transformacije, ili Jakobijan, kojim je izražena veza 0 j0 0 0 i i i x J , J r x r  ∂   ∂ ∂= =     ∂ ∂ ∂       (3.59) Kao što je predhodno rečeno, problemi nelinearnosti (geometrijske i materijalne) mogu se posmatrati nezavisno. Kao i u linearnoj teoriji elastičnosti možemo pretpostaviti da su pomeranja i dilatacije male veličine tako da jednačinu ravnoteže možemo postaviti u odnosu na nedeformisanu konfiguraciju ukoliko su bitni samo efekti materijalne nelinearnosti. Jednačine (3.28) i (3.31) prelaze u jednu jednačinu, oblika: n 1 n ijkl kl ij s ij ij V V D e e dV R e dV+⋅∆ ⋅δ∆ ⋅ = δ − σ ⋅δ∆ ⋅∫ ∫ (3.60) Rešavanjem jednačine (3.28) sračunavamo priraštaj pomeranja na osnovu kojeg sračunavamo: pomeranja, dilatacije i napone u konfiguraciji (n+1) koje ćemo označiti indeksom (i). Pošto smo izvršili linearizaciju jednačine ravnoteže, za pomeranja, deformacije i napone uzimamo približne vrednosti, pa je potrebno videti kolika je "neuravnoteženost opterećenja" između spoljašnjeg virtuelnog rada i virtuelnog rada sa sračunatim statičkim i kinematičkim veličinama u konfiguraciji (n+1). 46 n 1 n 1 (i) n 1 (i) e ij ij V P R e dV+ + +∆ = δ − σ ⋅δ∆ ⋅∫ (3.61) Da bi dalje redukovali neuravnoteženo opterećenje, potrebno je izvršiti još iteracija u okviru istog inkrementa opterećenja, dok ne dobijemo željenu tačnost po nekom unapred usvojenom kriterijumu konvergencije. Jednačina (3.60) prelazi u iterativnu jednačinu: n (m) (m) n 1 n (m 1) (m 1) ijkl kl ij e ij ij V V D e e dV R e dV+ − −⋅∆ ⋅δ∆ ⋅ = δ − σ ⋅δ∆ ⋅∫ ∫ (3.62) gde je priraštaj pomeranja odreden jednačinom: n 1 (m) n 1 (m 1) (m) i i iu u u + + −= + ∆ (3.63) Dobijeni sistem algebarskih jednačina rešavamo uobičajenom procedurom MKE. Kao što je već naglašeno, inkrementalno iterativni postupci rešavanja zadovoljavaju uslove ravnoteže samo u diskretnim tačkama. Kada konstitutivne relacije propisuju inkrementalnu vezu između napona i deformacija neizbežno je da postupci rešavanja koji se baziraju na inkrementalno iterativnim postupcima dovode do određene greške u rešenju. Izvor greške se nalazi u netačnosti integracije zakona tečenja duž putanje nezavisno promenljivih parametara, koji određuje iznos plastične deformacije tokom inkrementa a koja postoji u bilo kom vremenski diskretnom postupku rešavanja. U vremenski diskretizovanim postupcima uvek se implicitno usvaja da je priraštaj nezavisnih parametara sistema tokom trajanja inkrementa linearan. Kada je ravnoteža sila zadovoljena na početku i na kraju inkrementa i u slučaju kada priraštaj plastične deformacije nije linearan, rešenje u posmatranom inkrementu neće biti identično onom koje bi se dobilo kada bi se ovaj inkrement podelio na više sabinkremenata. 3.4 Numerička integracija konstitutivnih jednačina Na osnovu predhodnih razmatranja može se formulisati standardni Newton-Raphson- ov postupak za rešavanje nelinearnih problema [ ]27 , [ ]53 , [ ]62 , [ ]139 , [ ]161 , [ ]164 u kome se spoljašnje opterećenje ili zadato pomeranje nanosi u inkrementima. Do ravnoteže na kraju svakog od inkremenata dolazi se iterativnim postupkom, linearizacijom oko stanja koje odgovara predhodno uravnoteženom inkrementu spoljašnjeg opterećenja. U svakoj od iteracija pomeranja čvornih tačaka razmatranog tela izazivaju neuravnotežene čvorne sile, ( )in 1R+ , koje se određuje na osnovu izraza: 47 ( ) 0 ngausnelem nelem(i)T (i)Tin 1 n 1 n 1 n 1 0 n 1 n 1 n 1 0 ij ij ielem 1 ielem 1 igaus 1V R P B dV P B dV+ + + + + + + = = = = − σ = − σ∑ ∑ ∑∫ (3.64) pri čemu indeks (n+1) označava tekući increment, indeks (i) tekuću iteraciju u okviru inkrementa, a nelem i ngaus označavaju ukupan broj elemenata i Gausovih tačaka integracije po elementu. Vektor n 1P+ označava vektor spoljašnjeg opterećenja na kraju (n+1) inkrementa. Priraštaj čvornih pomeranja u okviru svake iteracije posmatranog inkrementa se određuje na osnovu rešenja sistema algebarskih jednačina ( ) ( )i in 1 n 1 tK u R + +⋅ = (3.65) pri čemu ( )in 1 tK + predstavlja tangentnu matricu krutosti sistema koja se dobija kao zbir matrica krutosti pojedinačnih elemenata. Tačnost rešenja sistema (3.64) u okviru unapred definisanog kriterijuma tačnosti, zavisi od postupka integracije konstitutivnih relacija u tzv. vremenskom koraku. Brzina konvergencije zavisi od načina računanja konstitutivne matrice. Kvadratna brzina konvergencije postiže se određivanjem tangentne konstitutivne matrice. 3.4.1 Algoritam povratnog preslikavanja [ ]99 Problem sračunavanja napona na osnovu poznatih inkremenata deformacija svodi se na integraciju konstitutivnih jednačina elasto-plastičnog problema. Istraživanja u ovoj oblasti su vrlo intezivna, jer od načina rešavanja integracije konstitutivnih jednačina zavisi tačnost rešenja razmatranog graničnog problema [ ]8 , [ ]28 , [ ]31 , [ ]42 , [ ]61 , [ ]100 , [ ]161 . Dalja razmatranja, u pogledu integracije konstitutivnih jednačina, prikazaćemo sa aspekta algoritma povratnog preslikavanja u plastičnosti. Za analizu elastičnih deformacija, nezavisno od toga da li se one posmatraju u sklopu infinitezimalnih ili konačnih deformacija, u računskoj mehanici najčešće se usvaja da je elastični tenzor konstantan i izotropan u svakoj ravnotežnoj konfiguraciji. Kod reševanja problema plastičnog deformisanja mehaničkog sistema kod koga su elastične deformacije relativno male u odnosu na odgovarajuće plastične deformacije, sasvim je svejedno na koji će način biti formulisan elastičan odgovor, pri čemu podrazumevamo korišćenje konstantnog tangencijalnog tenzora elastičnosti ili hiperelastičnog potencijala. 48 Analiza plastičnog ponašanja koja će biti sprovedena u okviru ovih razmatranja obuhvata matematički model materijala koga karakteriše povećanje elastičnog opsega u naponskom prostoru sa porastom plastične deformacije. Na taj način uvodi se familija površi tečenja kao funkcija napona i unutrašnjih promenljivih. Kao što je navedeno u predhodnom poglavlju, da bi se odredilo plastično deformisanje posmatranog tela potrebno je poznavati: zakon plastičnog tečenja, evolucionu jednačinu unutrašnjih promenljivih sistema i kriterijum opterećenja. Konstitutivne jednačine elasto-plastičnog materijala za vremenski nezavisnu plastičnost mogu se prikazati u sledećem obliku: e p ij ij ijε = ε + ε (3.66.1) ( )eij ij ,qσ = σ ε (3.66.2) ( )pij ijr , qε = ∆λ σ (3.66.3) ( )q h ,qα α= ∆λ ⋅ σ (3.66.4) Zahteva se da naponi pri bilo kom iznosu plastičnih deformacija leže na tačno utvrđenoj površi f 0= , poznatoj kao površ plastičnog tečenja. Priraštaj efektivne plastične deformacije odredićemo sa predhodnim ograničenjem i ograničenjima koja predstavljaju kriterijum opterećenja i rasterećenja. Ovaj kriterijum u Kuhn–Tucker–ovoj formi [ ]99 sadrži sledeća ograničenja: ( )f ,q 0σ ≤ (3.67) 0∆λ ≥ (3.68) f 0∆λ = (3.69) podrazumevajući njihovo simultano zadovoljenje u bilo kom trenutku vremena. Nejednakost (3.67) karakteriše elastični domen ponašanja materijala i usvaja se da je konveksnog oblika. Ova nejednakost može biti posmatrana i kao uslov da naponske veličine u odsustvu plastičnih deformacija ostanu unutar elastičnog domena. U toku procesa opterećenja uslovi (3.68) i (3.69) moraju biti simultano zadovoljeni. 49 Za f 0≺ iz uslova (3.69) sledi da je 0∆λ = što odgovara elastičnom ponašanju materijala. Plastično ponašanje materijala je određeno vrednošću 0∆λ ; odakle se iz uslova (3.69) zahteva zadovoljenje plastičnog kriterijuma f 0= . Obzirom na aditivnu strukturu elasto-plastičnih konstitutivnih jednačina, algoritam povratnog preslikavanja se javlja kao vrlo prikladan postupak za njihovu integraciju. Ovaj algoritam integracije napona u koraku zadovoljava uslov plastične konsistentnosti u svakom inkrementu i pokazuje dobre karakteristike u pogledu stabilnosti (pouzdanosti) postupka integracije, efikasnosti i tačnosti dobijenih rešenja. Njegova primena je uglavnom bila ograničena na jednostavne plastične modele sa konstantnim elastičnim modulom i linearnim zakonom ojačanja. Simo i Ortiz [ ]130 su usvojili sledeće pretpostavke: • elastičan odgovor materijala zasnovan je na generalnom hiperelastičnom modelu sa nekonstantnim modulom; • inkrementalne elastične konstitutivne relacije zamenjene su totalnim elastičnim relacijama čime je izbegnuta potreba za njihovom integracijom; • Plastični model materijala predstavljen je neasocijativnim zakonom tečenja i proizvoljnim plastičnim kriterijumom. U okviru ovog rada algoritam povratnog preslikavanja za integraciju napona je prikazan u materijalnoj deskripciji uvođenjem osnovnih naponskih i deformacijskih veličina modela merenih u odnosu na početnu konfiguraciju. Konstitutivne jednačine, saglasno njihovoj aditivnosti, mogu se razdvojiti na elastičan i plastičan deo. Ako pretpostavimo da je ukupni odgovor materijala elastičan onda se odgovarajuće jednačine možemo napisati u sledećem obliku: ( )tdpe =ε+ε=ε  (3.70.1) ( ) εε=σ  :eD (3.70.2) 0p =ε (3.70.3) 0q = (3.70.4) 50 Plastični deo odgovarajućih konstitutivnih jednačina, s obzirom na izraze (3.70.1 do 3.70.4), možemo izraziti u sledećem obliku: 0pe =ε+ε=ε  (3.71.1) ( )e pD :σ = − ε ε (3.71.2) ( )p r ,qε = ∆λ ⋅ σ (3.71.3) ( )q h ,q= ∆λ ⋅ σ (3.71.4) U elastičnim jednačinama predviđanja "zamrznut" je plastični odgovor materijala, odnosno plastične deformacije i unutrašnje promenljive tretirane su kao konstantne vrednosti, što inicira da je elastična deformacija ( ) ( ) pt tε = ε − ε poznata funkcija vremena. Proces plastične relaksacije definisan je jednačinama (3.71.1 do 3.71.4) i za njihovo reševanje potrebno je uvesti iterativni ciklus koji se nastavlja do zadovoljenja uslova plastičnosti ( )f ,q 0σ = u okviru unapred zadatog kriterijuma tačnosti. Polazeći od razdvojenih elasto-plastičnih konstitutivnih relacija (3.70) i (3.71) algoritam za integraciju napona se može postaviti u pogodnoj formi za konkretnu praktičnu primenu. Dobra predikcija plastične deformacije je važna radi ubrzanja konvergencije iterativnog rešenja. U nedostaku bolje informacije najčešće se polazi od nulte vrednosti, koja odgovara elastičnom rešenju. Jednačinama (3.70) definisano je elastično predviđanje ponašanja materijala i tako dobijeno stanje može se usvojiti kao početni uslov za proces plastične relaksacije određen jednačinama (3.71). Za jednostavne materijalne modele povratni put, definisan relaksacionim jednačinama, može biti određen analitički. Međutim, imajući u vidu činjenicu da je plastično ponašanje većine materijala definisano nelinearnim konstitutivnim jednačinama, analitička integracija relaksacionih jednačina je vrlo složena ili čak u opštem slučaju nije moguća pa smo prinuđeni na korišćenje nekog odgovarajućeg numeričkog postupka. Usvajanjem vremenski nezavisne plastičnosti, funkciju tečenja f možemo linearizovati u odnosu na zadnje sračunatu iteraciju na sledeći način: ( ) ( ) ( ) ( ) ( )n 1 (i) n 1 (i) n 1 (i) n 1 (i) n 1 (i) n 1 (i) n 1 (i) n 1 (i)f f , q , q : , q : q qq+ + + + + + + +∂φ ∂φ= σ + σ σ− σ + σ −∂σ ∂ (3.72) Diskretizacijom relaksacionih jednačina (3.71) konačnim priraštajima, odnosno: 51 n 1 (i 1) n 1 (i) n 1 (i) n 1 (i)D : r+ + + + +σ = σ −∆λ ⋅ (3.73.1) n 1 (i 1) n 1 (i) n 1 (i)q q h+ + + += + ∆λ ⋅ (3.73.2) i njihovom zamenom u linearizovanu jednačinu (3.72), zahtevajući pri tome zadovoljenje plastičnog uslova, određen je priraštaj efektivne plastične deformacije: ( ) ( ) ( ) ( ) ( ) ( ) in 1 i i i i in 1 n 1 n 1 n 1 n 1 f v : D : r : h + + + + + +∆λ = − ξ (3.74) pri čemu vrednosti ( ) ( )i in 1 n 1v , D ,+ + ( ) ( ) ( )i i in 1 n 1 n 1r , , h+ + +ξ odgovaraju stanju određenom naponima n 1 (i)+ σ i unutrašnjim promenljivim n 1 (i)q+ . U jednačini (3.74) je: v g= ∂ ∂σ parcijalni izvod funkcije plastičnog potencijala po naponima ( g je funkcija plastičnog potencijala) i f qξ = ∂ ∂ parcijalni izvod funkcije tečenja po unutrašnjim promenljivim. Numerički postupak, definisan jednačinama (3.73.1), (3.73.2) i (3.74) određuje povratni put napona preko diskretnih tačaka, dobijenih projekcijom predhodnog stanja (iteracije) na ravan f 0= . Iterativni ciklus nastavljamo sve dok trag ravni određene linearizacijom plastične funkcije f 0= ne postane tangenta na plastičnu površ u okviru unapred zadatog kriterijuma tačnosti. Geometrijska interpretacija predloženog algoritma prikazana je na slici 3.2. 52 Značajno je napomenuti da je u slučaju asocijativne idealne elasto-plastičnosti povratni put napona od elastičnog predviđanja do plastične površi određen linijom ("steepest descent"). Pravac ove putanje u svakom trenutku određen je lokalnom metrikom definisane tangentnim elastičnim modulom (slika 3.3). Konsistentnost ovog algoritma proističe iz metodologije razdvajanja elasto-plastičnog operatora koji obezbeđuje rezultujuće jednačine koje su u potpunosti saglasne sa konstitutivnim relacijama odgovarajuće konfiguracije. Njegova bezuslovna stabilnost je posledica odvojene bezuslovne stabilnosti elastičnog prediktora i plastičnih relaksacionih jednačina. Stabilnost elastičnog dela je očigledna, dok je bezuslovna stabilnost plastičnog dela omogućena time što su relaksacione jednačine disipativnog karaktera a odgovarajuće trajektorije obezbeđuju konveksnu plastičnu površ, odnosno plastični potencijal, odakle je definisan pravac plastičnog tečenja. Detaljni dokazi ovih tvrdnji prikazani su u radovima [ ]98 , [ ]99 , [ ]134 , [ ]135 , kao i više detalja može se naći u referentnom radu [ ]130 . Oblik i veličina konstitutivnog tenzora koji se primenjuje u numeričkoj integraciji jednačina zavisi od usvojenih konstitutivnih jednačina koji propisuju veze između priraštaja napona i pomeranja. Veličina elasto-plastičnog konstitutivnog tenzora može se odrediti na osnovu standardnih procedura u obliku: 53 n 1 n 1 n 1 n 1 n 1 e e e n 1 n 1 n 1 n 1 n 1 e ( D : r) ( D : v)D D v : D : r : h + + + + + + + + + + ⊗= − − ξ (3.75) gde je: n 1 eD + elastični konstitutivni tenzor, n 1r+ izvod funkcije tečenja po naponima, n 1v+ gradijent plastičnog potencijala po naponima, n 1+ ξ gradijent površi tečenja po unutrašnjim promenljivim, n 1h+ veličina koja propisuje funkciju ojačanja materijala. Relacijom (3.75) definisan je konstitutivni tangentni modul u najopštijem obliku, pretpostavljajući pri tome proizvoljnu plastičnu površ tečenja, proizvoljan pravac plastičnog tečenja kao i zakon ojačanja. Može se uočiti da je njegova simetričnost obezbeđena ako je u datom algoritmu integracije pretpostavljena asocijativna plastičnost sa konstantnim modulima ojačanja. Potrebno je napomenuti da se kontinualni (klasičan) tangencijalni operator [ ]134 razlikuje od takozvanog konsistentnog tangentnog operatora. Konsistentni tangentni operator direktno je zavistan od usvojenog algoritma za korekciju napona i on je usaglašen sa njim pa otuda i takav naziv. U prvom slučaju u okviru svakog koraka pravac plastičnog tečenja i moduli ojačanja su nam poznati i najčešće uzimaju konstantnu vrednost unutar celog inkrementa. Korekcija naponskih veličina implicitnim algoritmom koristi trenutne vrednosti pravca plastičnog tečenja, modula ojačanja i efektivne plastične deformacija, koje su nepoznate veličine jer i same zavise od trenutnog naponskog stanja, pa je samim tim jedan ovakav algoritam u svojoj osnovi iterativnog karaktera. Konsistentni tangentni operator ne utiče na tačnost konačnog rešenja ali omogućava kvadratnu konvergenciju Newton–Raphson – ovog postupka u rešavanju nelinearnih problema mehanike kontinuuma. Kod rešavanja nelinearnih jednačina inkrementalno iterativnim metodama, potrebno je definisati neki pogodan kriterijum za završetak iterativnog procesa. Na kraju svake iteracije potrebno je proveriti da li rešenje konvergira u dozvoljenim granicama tolerancije ili iterativni proces divergira. U proračunima najčešće se usvajaju kriterijumi konvergencije po pomeranjima, po neuravnoteženim silama i po inkrementalnoj unutrašnjoj energiji. 54 4. NELINEARNI MOHR-COULOMB-OV MODEL 4.1 Uvod Za razliku od gotovo svih drugih materijala koji se koriste u građevinskim konstrukcijama, kao što su beton ili opeka, osobine tla su promenljive, tako da se uglavnom moraju meriti za svaki slučaj posebno. Ponašanje tla u oblasti radnih napona je promenljivo sa izraženom anizotropijom na koju utiče niz propratnih pojava, naročito prisustvo vode. Pronalaženje računskog modela koji opisuje ponašanje tla u ovakvim slučajevima je izuzetno teško. Mehanička svojstva realnog tla pokazuju kompleksnu funkcionalnu zavisnost napona i deformacija. Veoma je težak i skup zadatak da se kontrolišu mnoge promenljive koje su vezane za tlo ili za opremu testiranja, da bi se dobila reprezentativna kriva napon deformacija za realni geomaterijal korišćenjem standardnih eksperimentalnih metoda. Statistička raspodela eksperimentalnih rezultata može biti primerena i primenjena u opisivanju kompleksnih funkcionalnih zavisnosti napona i deformacija. U pojedinim slučajevima za neki specifični test, moguće je dobiti veoma blizak matematički opis kompleksnog ponašanja realnog materijala. U domenu niskih i srednjih napona se odvijaju mikromehanički procesi: preraspodela čestica, klizanje, značajna redukcija u poroznosti, itd. U ovom regionu parametri materijala kao što su sadržaj vlage, početna gustina uzorka, su veoma značajni u opisivanju konstitutivne veze za gline [ ]74 i peskove [ ]81 . Određivanje mehaničkih karakteristika tla se sastoji od uzimanja uzoraka, njihovog opterećivanja po unapred definisanim putanjama napona i merenja nastalih pomeranja. Eksperimentalni podaci različitih opita prikazani u literaturi za neke vrste peskova, glina [ ]12 , šljunka [ ]56 i drugih geomaterijala [ ]18 koriste se za razvijanje matematičkih modela koji mogu opisati nelinerano ponašanje materijala. Pored jednostavnih naponskih putanja koje se karakterišu pravom linijom u sistemu glavnih napona, postoje i kompleksna opterećenja kod kojih se napon može menjati po krivoj liniji [ ]59 . Međutim, taj isti matematički model nije moguće primeniti u opisu ponašanja tog materijala pod drugim eksperimentalnim uslovima [ ]103 . Konstitutivne jednačine se razvijaju na osnovu empirijskih posmatarnja odnosa napona i deformacija. Formiranje konstitutivnog modela obuhvata sledeće korake: matematička formulacija, definisanje materijalnih parametara, način određivanja materijalnih 55 parametara iz laboratorijskih opita, poređenje rezultata modela sa rezultatima merenja i verifikacija modela u primeni rešavanja praktičnih problema. U radu je razvijen nelinearni Mohr-Coulombov konstitutivni model (u daljem tekstu NMC) u cilju opisivanja ponašanja tla. Model je formulisan na osnovu četiri polazne pretpostavke: • Ugao smičuće otpornosti nelinearne anvelope loma ( )ϕ ima hiperboličku zavisnost od normalnih napona; • Ugao dilatancije ( )δ se može odrediti kao odnos između ekvivalentne srednje zapreminske i ekvivalentne smičuće plastične deformacije; • Ugao dilatancije u linearnom modelu smičuće otpornosti jednak je maksimalnoj dilatanciji nelinearnog modela ∆ϕ ; • Bazni ugao nelinearne smičuće otpornosti ( )Bϕ se računa na osnovu aproksimacije da su tangentna smičuća otpornost za nelinearni i smičuća otpornost linearnog konstitutivnog modela jednake u određenoj vrednosti normalnog napona ( )Np . Prikazani hiperbolički konstitutivni model je generalisan za trodimenzionalni prostor napona i deformacija. Prvo je izložen kratak opis osnova matematičke formulacije hiperboličkog konstitutivnog modela. Zatim su prikazani rezultati numeričkog modeliranja i provera modela kroz poređenje sa referentnim eksperimentalnim rezultatima u opitima direktnog smicanja i triaksijalne kompresije. 4.2 Osnovne jednačine Mohr-Coulomb-ovog modela Radi preglednosti izlaganja u ovom odeljku, ponovićemo jednačinu (2.93.1), za LMC kriterijum tečenja: 1 2d Isin sinf I cos sin c cos 0 33  θ ⋅ φ= ⋅ θ + − ⋅ φ − ⋅ φ =   (4.1) gde je: θ Lode-ov ugao i može se napisati u obliku 1 3/ 23d 2d1 3 3sin I I3 2 − − θ = − ⋅ ⋅ ⋅    , 1I je prva invarijanta tenzora napona, 56 2dI i 3dI druga i treća invarijanta devijatorskog tenzora napona, c je kohezija i φ je ugao smičuće otpornosti. Ugao smičuće otpornosti i kohezija određuju svojstvo materijala. LMC model zahteva poznavanje pet materijalnih konstanti. Tri parametra: kohezija c , ugao smičuće otpornosti φ i ugao dilatancije δ  opisuju plastično deformisanje modela dok druga dva parametra E i ν , opisuju elastično ponašanje materijala. U referentnom radu [ ]109 može se naći diskusija o parametrima materijala, različitim modelima kao i odgovarajući primeri analize. Ako su primenjeni asocijativni uslovi matematičke teorije plastičnosti, onda se broj parametara redukuje na četiri jer je u tom slučaju δ = φ . Odgovarajući izvod LMC funkcije tečenja po naponu, može se izračunati korišćenjem pravila izvoda složene funkcije, prema sledećem izrazu: 2d1 ij 1 ij 2d ij ij IIf f f f I I ∂∂∂ ∂ ∂ ∂ ∂θ= ⋅ + ⋅ + ⋅∂σ ∂ ∂σ ∂ ∂σ ∂θ ∂σ (4.2 Prema tome, diferenciranjem funkcije tečenja (4.1), dobijaju se sledeće relacije: 1 12 2d1 ij ij 3d 2d ij ij3/ 2 23 3ij 2d 3d 2d2 r ( , ) If 1 r ( , ) 3 3sin S I I S t 3 2 22 I 1 ( I I ) − − − θ φ ⋅∂ θ φ  = − φ⋅δ + ⋅ + ⋅ ⋅ ⋅ ⋅ − ∂σ  − ⋅ ⋅ (4.3) gde je 1 sin sinr ( , ) cos 3 θ⋅ φθ φ = θ+ (4.4) 2 cos sinr ( , ) sin 3 θ⋅ φθ φ = θ+ (4.5) 2 ij ik kj 2d ij3t S S I= − δ (4.6) Da bi se kvantifikovao fenomen povećanja zapremine sa porastom smičuće deformacije, Taylor [ ]144 je uveo pojam ugao dilatancije δ . Kao što je već naglašeno, u ovom radu uvedena je pretpostavka da je ugao dilatancije definisan kao odnos između priraštaja ekvivalentne srednje zapreminske i priraštaja ekvivalentne smičuće plastične deformacije LMC. Ovo je osnovna pretpostavka u formulisanju NMC modela. Ekvivalentna plastična deformacija pε je definisana na osnovu priraštaja plastične deformacije pijdε : ( )1 2p p pij ijd dε = ε ⋅ ε∫ (4.7) 57 p vε i pdε su zapreminska i devijatorska komponenta ekvivalentne plastične deformacije p ,ε definisane sledećim izrazom: ( )1 2p p pv kk kk1 d d3ε = ε ⋅ ε∫ (4.8) ( )1 2p p pd ij ijde deε = ⋅∫ (4.9) gde je pkkdε inkrement zapreminske plastične deformacije i pijde je increment devijatorske plastične deformacije. Odnos između ekvivalentne plastične deformacije pε i odgovarajuće zapreminske i devijatorske komponente pvε i pdε dat je izrazom p p 2 p 2 v d( ) ( )ε = ε + ε (4.10) U jednačinu (4.10) uvodimo parametar - ugao dilatancije δ ( ) ( )2 2p p p p pd v d d1 3 sin 1ε = ε ε ε + = ε ⋅ δ + (4.11) Iz jednačine (4.11) parameter δ se može izraziti ( )21 p pv d1sin 13−  δ = ⋅ ε ε −   (4.12) Prema jednačinama (4.8) i (4.9), koristeći izraze (4.3), (4.4) i (4.5) odgovarajući priraštaji ekvivalentnih plastičnih deformacija su oblika: ( )1 2p p pv kk kk1 1d d d d sin3 3ε = ε ⋅ ε = λ ⋅ ⋅ φ (4.13) ( ) 1 2 1 2 2 p p p d ij ij 1 2 (sin )d de de d sin cos sin 2 63  φε = ⋅ = λ ⋅ + ⋅ θ⋅ θ ⋅ φ +   (4.14) Kao što je poznato, kod LMC odgovarajući priraštaji plastičnih deformacija su konstantni, pa se samim tim za triaksijalno stanje a na osnovu jednačine (4.11), odnos ukupne ekvivalntne srednje zapreminske plastične deformacije i ukupne ekvivalentne devijatorske plastične deformacije može se zameniti sa odgovarajućim priraštajima. Koristeći jednačine (4.12), (4.13), (4.14) za vrednost Lodeovog ugla 6θ = −π ugao dilatancije CTCδ za triaksijalno stanje kompresije može se izraziti u funkciji ugla smičuće otpornosti: CTC 1 2 2 sinsin 3 3 3 sin (sin ) −  φ δ = ⋅ − ⋅ φ + φ  (4.15.1) 58 Ako u jednačine (4.12), (4.13), (4.14) unesemo vrednost za Lodeov ugao 6θ = π za stanje redukovane triaksijalne ekstenzije ugao dilatancije RTEδ može se izraziti u funkciji ugla smičuće otpornosti: RTE 1 2 2 sinsin 3 3 3 sin (sin ) −  φ δ = ⋅ + ⋅ φ + φ  (4.15.2) Na osnovu predhodno navedenih jednačina i za vrednost Lodeovog ugla 0θ = ugao dilatancije SSδ za triaksijalni stanje prostog smicanja može se izraziti u obliku: SS 1 2 2 sinsin 3 3 (sin ) −  φ δ = ⋅ + φ  (4.15.3) Ovako definisan ugao dilatancije sa poznatom vrednošću za ugao smičuće otpornosti LMC poslužio nam je da definisemo materijalne parametre NMC bazni ugao smičuće otpornosti Bϕ , maksimalnu ugaonu razliku, ∆ϕ i odgovarajući normalni napon Np . Za triaksijalne opite: standardne kompresije (CTC), redukovane ekstenzije (RTE), i prostog smicanja (SS) tabelarno je prikazan ugao dilatancije ( )δ izračunat korišćenjem jednačina (4.15) na osnovu ugla smičuće otpornosti φ LMC za korak od pet stepeni. [ ]0φ CTCδ RTEδ SSδ 10 2.962 2.940 2.695 15 4.624 3.568 3.995 20 6.402 4.543 5.240 25 8.287 5.421 6.416 30 10.263 6.209 7.513 35 12.311 6.911 8.522 40 14.399 7.534 9.440 45 16.492 8.082 10.263 Tabela 4.1 Ugao dilatancije u funkciji ugla smičuće otpornosti LMC Na slici 4.1 je grafički prikazan ugao dilatancije u funkciji od ugla smičuće otpornosti LMC modela, za korak od jednog stepena, za prethodno navedene triaksijalne opite. 59 Slika 4.1 Ugao dilatancije u funkciji ugla smičuće otpornosti LMC 4.3 Konstitutivne jednačine NMC modela Široko je korišćena hiperbolička funkcija za simulaciju naponsko deformacionih krivih u analizi konačnih elemenata, koju su formulisali Duncan i Cheng [ ]40 , koristeći Kondner-ov [ ]63 dijagram napon –deformacija dobijen pri troosnom testu. Maksimović [ ]80 predlaže izraz za nelinearnu anvelopu u obliku hiperboličke funkcije: B 1 N1 I 3P ∆ϕϕ = ϕ + + (4.16) gde je: Bϕ bazni ugao; ∆ϕ maksimalna dilatancija; a NP normalni napon pri kome maksimalna dilatancija dostiže polovinu svoje vrednosti. 60 Funkcija tečenja NMC modela je direktno izvedena iz LMC kriterijuma tečenja zamenom ugla smičuće otpornosti Maksimovićevom [ ]75 , [ ]76 funkcionalnom vezom (4.16), u jednačini (4.1) tako da se dobija sledeći izraz: 1 2d sin sin If I cos sin c cos 0 33 θ⋅ ϕ = ⋅ θ+ − ⋅ ϕ− ⋅ ϕ =   (4.17) Shodno jednačini (4.2) odgovarajući parcijalni izvodi nelinearne funkcije tečenja (4.17) po naponima su: 2d1 ij ij 1 1 1 12 2d1 ij 3d 2d ij ij 3/ 2 23 32d 3d 2d2 IIf 1 sin sin cos c sin 3 3 3 I I r ( ) Ir ( ) 3 3S I I S t 2 22 I 1 ( I I ) − − −   ∂ ∂ϕ ∂ϕ= − ϕ− − ⋅ θ ⋅ ϕ ⋅ − ⋅ ϕ ⋅ δ +   ∂σ ∂ ∂    θ ⋅θ  + ⋅ + ⋅ ⋅ ⋅ ⋅ −  − ⋅ ⋅ (4.18.1) gde je ( )21 N 1 NI 3p 1 I 3P ∂ϕ ∆ϕ= −∂ + (4.18.2) 1 sin sinr ( , ) cos 3 θ⋅ ϕθ ϕ = θ+ (4.18.3) 2 cos sinr ( , ) sin 3 θ ⋅ ϕθ ϕ = θ + (4.18.4) i ijt određen izrazom (4.6). Za ravno stanje deformacije NMC uslov tečenja, izražen relacijom (4.17), možemo prikazati u sledećem obliku: 1 2d 1 If I tan d 0 3 = − ⋅ β − = (4.19.1) gde je: 1d c cos , tan sin= ⋅ ϕ β = ϕ (4.19.2) Geometrijska interpretacija parametra NMC modela za ravno stanje deformacije prikazana je na Slici 4.2. Prikazan bazni ugao Bβ je konstantan u odnosu na normalni napon, ugao smičuće otpornosti β teži vrednosti Bβ kada normalni napon teži beskonačnosti, za sve vrednosti normalnog napona bazni ugao Bβ je manji ugao od ugla smičuće otpornosti α LMC predstavljenog jednačinama (2.89) i (2.90). 61 Slika 4.2 NMC u prostoru 1 2dI I− Problemu integracije konstitutivnih relacija poklanja se posebna pažnja, jer od načina njegovog rešavanja zavisi uspešnost modeliranja ponašanja tla u plastičnom domenu deformisanja [ ]62 . 62 4.4 Verifikacija predloženog NMC modela 4.4.1 – Opit direktnog smicanja - Zbijeni pesak Opit direktnog smicanja počinje tako što se prvo izvrši konsolidacija uzorka opterećenjem hidrostatičke kompresije. Nakon toga, uzorak se podvrgava direktnom smicanju tako što se promena dimenzije uzorka vrši i meri u odnosu na primenjenu tangencijalnu silu. Dobija se eksperimentalni dijagram ili tabela, tangencijalni napon u odnosu na deformaciju za jednu vrednost inicijalnog hidrostatičkog napona. Da bi se dobila standardna devijacija usled nehomogenosti materijala i promenljivih vezanih za eksperimentalni mehanizam, potrebno je ponoviti taj isti eksperiment. Nakon toga, povećava se vrednost hidrostatičke kompresije i procedura se ponavlja. Da bi se dobile tri eksperimentalne vrednosti tangencijalnog napona u zavisnosti od normalnog napona potrebno je najmanje šest uzoraka materijala. U literaturi je prava retkost naći rezultate ponovljenih eksperimenata. Zbog toga je važno imati mogućnost da se na gore opisani način dođe do parametara nelinearnog hiperboličkog modela. Tabela 4.2 prikazuje rezultate merenja koje su Maksimović i Santrač [ ]79 naveli za opit direktnog smicanja materijal na uzorku zbijenog peska. Tabela 4.2 Rezultati merenja zbijenog peska Rezultati proračuna parametara NMC modela se dobijaju postavljanjem tri linearne jednačine sa tri nepoznate: bazni ugao smičuće otpornosti Bϕ , maksimalnu ugaonu razliku∆ϕ i normalni napon Np . Rešavanjem sistema linearnim jednačina dobijaju se: 0 B 31.4ϕ = , 014.2∆ϕ = , 2Np 559.3kN / m= Ovo su rezultati proračuna na osnovu eksperimentalno merenih podataka. S druge strane, primenimo prikazani model na sledeći način: prvo fitujemo eksperimentalne podatke sa LMC i dobijemo vrednost ugla smičuće otpornosti 035φ = . Iz Tabele 4.1 za odgovarajući ugao smičuće čvrstoće od 35° se dobija 012.3∆ϕ = za CTC opit. Za vrednost normalnog napona N2pσ = nalazi se vrednost baznog ugla od 0B 30.9ϕ = . [ ]kPaσ 250.0 500.0 1000.0 [ ]kPaτ 219.0 403.5 740.0 63 Ako uporedimo rezultate dobijene proračunom na osnovu eksperimentalnih podataka i na osnovu modela vidimo da se bazni ugao razlikuje za samo 0.5° a ∆ϕ se razlikuje za 1.9°. Ovi rezultati ukazuju da predloženi model dobro opisuje ponašanje zbijenog peska. Ovaj primer ukazuje na mogućnost da se iz LMC modela i na osnovu Tabele 4.1 dođe do materijalnih konstanti NMC modela u i slučaju kada ne posedujemo dovoljan broj eksperimentalnih podataka. 4.4.2 Londonska glina Atkinson i Farrar [ ]3 su prikazali rezultate triaksijalnog ispitivanja uzoraka zbijene visokoplastične gline. U intervalu efektivnih napona od 150 do 300 kPa, anvelopa je prilično linearna i definisana sa parametrima 2c 25kN / m= i 016φ = . Na osnovu prikazanog modela iz Tabele 4.1 za ugao smičuće otpornosti 016φ = dilatancija je izračunata 05∆ϕ = , a bazni ugao 0B 14.34ϕ = . Usvojena je računska vrednost za materijalnu konstantu 2 NP 28, 20 kN / m= na osnovu podataka prikazanih u knjizi [ ]80 (slika 5.40 na strani 201). Rezultati su normalizovani u odnosu na parametar Np i prikazani na slici 4.3.1. 64 4.4.3 Drobljeni antracit U knjizi [ ]80 su prikazani rezultati merenja drobljenog uglja antracita (slika 5.43.a na strani 205). Ugao smičuće otpornosti LMC je 037φ = . Dilatancija NMC modela jednaka je 013.1∆ϕ = (određena na osnovu slike 4.2 za CTC test). Odgovarajući bazni ugao NMC modela je 0B 32.6ϕ = . Rezultati dobijeni za nelinerane parametre a na osnovu eksperimentalnih podataka su bazni ugao 0B 32.2ϕ = , dilatancija 016.5∆ϕ = i normalni napon 2Np 104kN / m= . Rezultati su normalizovani u odnosu na parametar Np i prikazani na slici 4.3.2. Razlike između parametara dobijenih eksperimentalnim putem i na osnovu modela su sledeće: u baznom uglu je 0.4° a u parametru 03.4∆ϕ = . 65 4.4.4 Standardna triaksijalna kompresija (CTC-test) U testu standardne triaksijalne kompresije, uzorak materijala pripremljen u obliku cilindra je podvrgnut inicijalnom pritisku fluida 2 3 0σ = σ = σ (takozvani pritisak komore), a onda je vršeno povećanje normalnog napona 1σ u aksijalnom pravcu. Test se obično vrši na nekoliko različitih nivoa početnih napona 0σ . Usled nehomogenosti materijala i promenljivih vezanih za aparaturu, poželjno je ponoviti sva merenja kako bi se izračunale srednje vrednosti i standardna devijacija. U ovom primeru navodimo rezultate modeliranja i eksperimentalnog ispitivanja "rockfill" materijala (zaglinjena drobina) za Roadford branu [ ]104 . Parametri materijala navedeni u pomenutoj referenci prikazani su na slici 4.4.2. Materijal je podvrgnut CTC opitu na sledećim vrednostima početnih napona: 0σ = 225kN / m , 296kN / m i 2266kN / m . Parametri NMC modela su određeni na bazi LMC za ugao smičuće otpornosti 036φ = i Tabele 4.1: bazni ugao 0B 31.8ϕ = , 012.7∆ϕ = i 2Np 125kN / m= . U okviru ovog primera urađena je numerička analiza korišćenjem NMC modela zasnovanog na neasocijativnom uslovu tečenju. U funkciju plastičnog potencijala uključen je ugao dilatancije δ umesto ugla smičuće otpornosti ϕ . Na slici 4.4.1 su prikazane krive zavisnosti oktaedarskog napona octτ od aksijalne deformacije aε . Eksperimentalni rezultati zaglinjene drobine su prikazani simbolima a odgovarajući rezultati numeričke analize su prikazani linijama. Poređenje napona pri kojima octτ dostiže maksimum za tri vrednosti početnih napona: 0σ = 225kN / m , 296kN / m i 2266kN / m , pokazuje dobro slaganje rezultata proračuna dobijenih primenom NMC modela sa odgovarajućim eksperimentalnim krivama. Na slici 4.4.2 su prikazane krive zavisnosti zapreminske deformacije od aksijalne deformacije. Kao što se može videti sa dijagrama, u početku je nagib krive negativan jer dolazi do smanjenja zapremine materijala pod uticajem smičućih naprezanja. Negativni nagib je konstantan do određene vrednosti kada nagib postaje jednak nuli. U NMC modelu ta tačka odgovara naponu na kojem octτ dostiže maksimum (početak plastičnog deformisanja u modelu). Nakon minimuma, sa povećanjem priraštaja aksijalne deformacije dolazi do povećanja zapremine materijala. Na osnovu vrednosti prikazanih na slici 4.4.2 može se konstatovati zadovoljavajuće slaganje rezultata numeričkog modeliranja CTC opita sa odgovarajućim eksperimentalnim rezultatima. 66 67 5. PRORAČUN TRAKASTIH TEMELJA NA DEFORMABILNOJ PODLOZI 5.1 Uvod Temelj objekta je deo konstrukcije čija je namena da prenese i rasporedi opterećenje od objekta na tlo. Dimenzije i oblik temelja zavise od opterećenja koje se preko njega prenosi, konstrukcije objekta i potrebe da se opterećenje što ravnomernije prenese na tlo, koje je kao materijal manje otporno od materijala objekta. Dimenzije temelja zavise od osobina tla, vrste materijala temelja i veličine opterećenja koju temelj treba da prenese na tlo. Temelji mogu biti kruti i deformabilni. U grupu krutih temelja spadaju temelji čije se deformacije pri njihovom proračunu mogu zanemariti. Zatezanja usled savijanja temelja su neznatna. Otpornost materijala na zatezanje je u tom slučaju dovoljna za obezbeđenje potrebne sigurnosti temelja. Izvode se od nearmiranog betona. U grupu deformabilnih temelja spadaju temelji čije se deformacije pri njihovom proračunu ne mogu zanemariti. U njima se usled savijanja, javljaju naponi zatezanja o kojima se mora voditi računa pri dimenzionisanju deformabilnih temelja. Izvode se od armiranog betona. Trakasti temelji se po pravilu predviđaju ispod zidova. Ako je opterećenje po dužini zida konstantno, uslovi u svim poprečnim presecima su isti, pa su dimenzije ovakvih temelja na celoj njegovoj dužini iste. Samo ako se menja opterećenje po dužini zida menjaće se i širina temeljne trake. Temeljno tlo se može nalaziti u uslovima ravanskog stanja deformacija ili u uslovima prostornog stanja stanja napona i deformacija. Ako se podloga nalazi u uslovima ravanskog stanja deformacija, ma koji element određene širine izdvojili iz temeljne trake i analizirali u poprečnom pravcu nalaziće se u svim presecima u uslovima u kojima se nalazi i svaki drugi element te konstrukcije. Zato je za analizu dovoljno usvojiti i razmatrati samo jedan takav element. Ako je dužina temeljne konstrukcije veća od njene trostruke širine možemo smatrati da će se takva konstrukcija, na delu udaljenom za širinu temeljne konstrukcije od njenih krajeva nalaziti u uslovima ravanskog stanja deformacija. Elementi na krajevima nalaziće se u drugačijim uslovima. U svim drugim slučajevima imaćemo prostorno stanje napona i deformacija. 5.2 Nosivost tla Maksimalni pritisak temelja na tlo pri kome dolazi do loma zove se granična nosivost tla ili granični pritisak na tlo. Radi ocene bezbednosti objekta, s obzirom na stabilnost tla ispod temelja, potrebno je da se odredi granični pritisak na tlo. Teorijske metode određivanja granične nosivosti tla bazirane su na pojednostavljenim pretpostavkama osobina tla i 68 raspodele pritiska tla. Bez obzira na ova pojednostavljenja, poređenje teorijskih rezultata granične nosivosti tla i rezultata dobijenih na osnovu eksperimentalnih ispitivanja pokazuju relativno dobra slaganja. Greške izmedju teorijskih i eksperimentalnih ispitivanja nosivosti tla nisu mnogo veće od grešaka koje se javljaju pri određivanju stabilnosti konstrukcija od drugih građevinskih materijala. 5.2.1 Teorija nosivosti Terzaghi-a Pretpostavku da se tlo ponaša kao idealno plastično telo Terzaghi je prihvatio i primenio uz to da je otpornost određena Coulomb-ovim uslovom loma. Terzaghijeva metoda za određivanje graničnog pritiska je tako postala osnova za proračine nosivosti tla metodom graničnog stanja plastične ravnoteže. U klasičnoj mehanici tla trebalo je naći jednostavne metode za rešavanje nelinearnih graničnih problema, uz korišćenje osnovne pretpostavke da se tlo ponaša kao idealno plastično telo. Ovu pretpostavku prvi je uveo Prantl za metale, a Terzaghi je prihvatio i primenio na tlo čija je otpornost određena Coulomb-ovim uslovom loma. Metodu za određivanje graničnog pritiska na tlo Terzaghi je objavio 1943. god. i od tada ona je postala osnova za proračine nosivosti tla metodom graničnog stanja plastične ravnoteže. 69 Radi određivanja faktora nosivosti Terzaghi je uveo sledeće pretpostavke: Temelj je na kontaktu sa tlom hrapav, što znači da nema horizontalnih pomeranja tla u kontaktnoj površini i da se klin, formiran u tlu ispod temelja, može tretirati kao sastavni deo temelja; Uglovi između bočnih strana klina i horizontale ψ jednaki su uglu smičuće otpornosti tla φ ; Površine klizanja prostiru se do nivoa kontaktne površine, a uticaj tla iznad kontaktne površine se uzima da je ekvivalentan ravnomerno podeljenom pritisku od težine tla iznad kontaktne površine (slika 5.1). Opterećenje od trakastog temelja širine B prenosi se na tlo kao ravnomerno podeljeni pritisak, što znači da se uticaj krutosti temelja na raspodelu pritisaka u kontaktnoj površini zanemaruje. Ispod temelja se nalazi homogeno tlo čija je jedinična težina γ . Pri graničnom pritisku u tlu ispod temelja formira se klin. Opterećenje od temelja se prenosi preko ovog klina koji potiskuje tlo u stranu, formirajući sa obe strane oblasti u kojima su ispunjeni uslovi granične ravnoteže. Svaka od ovih oblasti sastoji se iz dva dela: zone radijalnih površi klizanja neposredno pored klina i zone ravnih površina klizanja koje se nalaze iza radijalnih. Pomeranju klina, pri graničnom pritisku na tlo, suprostavljaju se otpori na bočnim stranama klina. Otpori se sastoje od otpora usled trenja i kohezije koja deluje duž ravn AB klina pri njegovom sleganju. Rezultanta otpora pri pasivnom stanju granične ravnoteže zaklapa ugao φ sa normalom na bočnu ravan AB klina. Uglovi između bočnih strana klina i horizontale jednaki su uglu ψ . Izraz za granični pritisak na tlo, koji je predložio prof. K.Terzaghi, je oblika: f c qq c N q N 1 2 B Nγ= ⋅ + ⋅ + ⋅ γ ⋅ ⋅ (5.1) gde su: c qN , N i Nγ faktori nosivosti koji zavise od ugla smičuće otpornosti i pretpostavljenog oblika oblasti sa graničnim naponskim stanjem, fq D= γ ⋅ pritisak od težine tla iznad ravni koja se nalazi u nivou kontaktne površi temelja. Faktori nosivosti su dobijeni iz izraza: ( ) ( ) 3 2 tan q eN 2 cos 4 2 π −φ ⋅ φ = ⋅ π + φ (5.2) ( )c qN N 1 ctg= − ⋅ φ (5.3) Za faktor nosivosti Nγ ne postoji rešenje u zatvorenom obliku jer on zavisi od geometrije pretpostavljenih oblasti u kojima su ispunjeni uslovi granične ravnoteže tla. Brinch Hansen je predložio rešenje za faktor nosivosti Nγ oblika: ( )qN 1,8 N 1 tgγ = ⋅ − ⋅ φ (5.4) 70 Granični pritisak na tlo, kao što se vidi iz izraza (5.1), zavisi od ugla smičuće otpornosti tla φ , zapreminske težine tla γ , širine temelja B , kohezije c i pritiska od težine tla q iznad ravni koja se nalazi u nivou kontaktne površine. Na prvi član u izrazu za nosivost utiče, pored ugla smičuće otpornosti, i veličina kohezije. Ako je ugao smičuće otpornosti mali nosivost tla zavisi, uglavnom, od ovog člana. Na drugi član u izrazu (5.1) utiče, pored ugla smičuće otpornosti, i težina tla iznad kontaktne površine. Udeo ovog člana može biti veoma mali ako tlo ima mali ugao smičuće otpornosti. Treći član u izrazu (5.1) je proporcionalan jediničnoj težini tla koje se nalazi ispod temelja u oblastima sa graničnim naponskim stanjem. Terzaghi-evo rešenje pokazuje da u slučaju peska postoji nosivost i za temelje na samoj površini a što je saglasno sa fizičkom pojavom. Određivanjem nosivosti tla ispod plitkih temelja metodom graničnog stanja plastične ravnoteže bavili su se mnogi autori, i svaki autor je u svoju metodu uveo neku od različitih pretpostavki o obliku klizne ravni i osobinama tla, da bi pojednostavio rešenje za granični pritisak na tlo. Iz tih razloga sračunate veličine graničnog opterećenja se međusobno razlikuju. Postoji više metoda za proračun nosivosti temelja prema stanju granične ravnoteže. Navešćemo samo neke autore tih metoda: Brinch-Hansen, Meyerhof, Berezantsev, Vesić, Sokolovski, Buisman, Schwedler, Skempton i drugi. U cilju provere ponašanja trakastih temelja prema stanju granične ravnoteže urađen je veliki broj studija nosivosti tla metodom karakteristika. Ovom metodom rešavanje problema definisanih diferencijalnim jednačinama prvog i drugog reda svodi se na rešavanje obične diferencijalne jednačine jednim od numeričkih postupaka. Sokolovski (1960) je razmatrao čisto plastični problem tražeći rešenja po naponima, dok je određivanje deformacija i brzina deformisanja zanemario. Metodom karakteristika rešeni su zadovoljavajuće mnogi za inženjersku praksu značajni granični zadaci. Međutim, potrebno je istaći da ta rešenja podrazumevaju idealnu plastičnost i da se proračun izvodi za neki unapred određeni – pretpostavljeni mehanizam - potencijalnu kliznu ravan. Ovo pretpostavlja da tlo u zoni loma ima nepromenjenu gustinu pre i posle loma, a da maksimalnu vrednost ugao smičuće otpornosti dostiže istovremeno u svim tačkama. U realnim situacijama način i oblik (mehanizam) loma se znatno razlikuje od onog koji predviđaju konvencionalne metode proračuna. Nejednake specifične deformacije u zoni 71 smicanja dovode do mobilizacije smičuće otpornosti i loma duž linija loma – pravaca u različitim momentima. Sa porastom opterećenja, raste i neravnomernost specifičnih smičućih deformacija. Prema tome, lom nije do koga dolazi duž neke unapred utvrđene klizne površi – mehanizma već proces koji se progresivno razvija i širi u pravcima najmanjih otpora. 5.3 Ilustracija primene MKE na konkretnim primerima Opisane teorijske postavke, kao i numeričke procedure koje su zasnovane na njima (poglavlja 2 i 3), primenjene su u ovom odeljku za ilustraciju proračuna trakastih temelja na deformabilnoj podlozi metodom konačnih elemenata. Odabrana su dva konkretna primera proračuna trakastih temelja na deformabilnoj podlozi koji ilustruju osnovne karakteristike prikazane metodologije. Tlo je modelirano kao elasto-plastični kontinuum korišćenjem NMC modela uz uvođenje sledećih pretpostavki: 1. Raspodela pritisaka u kontaktnoj površini temeljne trake zavisi od njene krutosti; 2. Usvaja se da je elastični tenzor konstantan i izotropan u svakoj ravnotežnoj konfiguraciji za analizu elastičnih deformacija; 3. U svakoj tački kontaktne površine vertikalno pomeranje tačaka ose nosača jednako sleganju podloge u odgovarajućoj tački; 4. Kontaktna površina temelja je na površini tla. Uticaj dubine fundiranja se zamenjuje opterećenjem jednakim težini tla iznad kontaktne površi temelja, ravnomerno raspoređenim po celoj površini tla. Uticaj tla iznad kontaktne površine temelja zamenićemo sa ravnomerno podeljenim opterećenjem fDγ , gde je γ jedinična težina tla a fD dubina fundiranja; 5. Vertikalni napon od sopstvene težine tla u horizontalnoj ravni koja se nalazi na dubini z jednak je: ( )z fD zσ = γ ⋅ + (5.5) gde je z dubina posmatrane tačke ispod kontaktne površine temelja. Horizontalni napon od sopstvene težine tla, na dubini z ispod kontaktne površine temelja jednak je: y 0 zkσ = ⋅σ (5.6) gde je 0k koeficijent bočnog pritiska tla. 72 Primer 5.3.1 Kruta temeljna traka na pesku Kao prvi primer proračuna trakastog temelja na deformabilnoj podlozi MKE odabran je primer prikazan u radu: J. T. Christian, A. J. Hagmann and W.A. Marr: "Incremental Plasiticity Analysis of Frictional Soils", International Journal for Numerical and Analytical Methods in Geomechanics, Vol. 1, 343-375, (1977). Primer je odabran i zbog mogućnosti upoređenja rezultata sa rezultatima prikazanim u referentnom radu. Zbog činjenice da je ovo važna referenca navešću kratki rezime rada. Materijalni parametri prikazani u radu mogu se direktno koristiti kao ulazni parametri NMC modela, uključujući i određivanje odgovarajućih parametra smičuće otpornosti: bazni ugao 0B 26.6ϕ = , ugao dilatancije 010.3∆ϕ = i normalni napon 2Np 500kN / m= na osnovu ugla smičuće otpornosti 030φ = LMC (slika 4.2 - CTC test). Autori sa MIT-a su za proračun potporne konstrukcije i nosivosti trakastog temelja koristili Mohr-Coulomb-ov kriterijum loma primenom asocijativne i neasocijativne teorije plastičnosti. Utvrdili su veoma dobro slaganje između pritiska i otpora tla na osnovu rezultata dobijenih metodom graničnog stanja plastične ravnoteže i prikazanog numeričkog rešenja za proračun potporne konstrukcije. Međutim, kod proračuna nosivosti trakastog temelja na pesku konstatovali su odstupanje numeričkog rešenja u odnosu na rešenje nosivosti određeno metodom granične ravnoteže. Prikazanom analizom obuhvatili su uticaj različitih faktora na proračun nosivosti trakastog temelja: broj konačnih elemenata u diskretizovanom modelu, veličinu inkrementa opterećenja (pomeranja) i različite polazne pretpostavke u numeričkom postupku za integraciju konstitutivnih relacija. Kao što je poznato, diskretizovanje, izbor oblika i broj konačnih elemenata zavise od prirode problema i zahtevane tačnosti rešenja. U razmatranom graničnom zadatku nosivosti i sleganja trakastog temelja zadato opterećenje i obuhvaćena geometrija zapremene tla ispod opterećenog temelja su simetrični tako da je mreža konačnih elemenata formirana za odgovarajući simetričan problem. Posmatrali su ravansku deformaciju u dreniranim uslovima. Zamenu poluprostora modelom određenih dimenzija izvršili su diskretizacijom domena ispod temelja mrežom konačnih elemenata do dubine 2B i širine 2,5B, gde je B širina temeljne trake. Po bočnim konturama, sprečena su horizontalna a po osnovi i vertikalna i horizontalna pomeranja (fiksni oslonac). U radu je usvojen 2-D četvorougaoni izoparametarski konačni element. Uticaj broja konačnih elemenata na rešenje graničnog problema, analizirali su diskretizacijom na 30 i 96 elemenata. Zbog neznatnih razlika dobijenih na osnovu prikazanih rezultata proračuna autori su izveli zaključak da usvojeni broj konačnih elemenata ne utiče bitno na rešenje razmatranog graničnog problema. 73 Što se tiče veličine inkrementa opterećenja, zadavali su inkrementalno pomeranje u iznosima od 0,005, 0,01 i 0,02 m. Broj iteracija u okviru inkrementu je bio najviše do dve iteracije što proizilazi iz vrlo tačnog i efikasnog postupka integracije napona. Autori su testirali tri numeričke sheme integracije napona. Prvi postupak je obuhvatao "korekciju" vektora napona na površ tečenja uz održavanje konstantnog srednjeg napona, drugi korekciju napona u pravcu normale na anvelopu loma i treći je obuhvatao "korekciju" napona uz konstantan vertikalni napon. Radi diskretizacije domena ispod temeljne trake, u okviru ovog brojnog primera, usvojena je mreža konačnih elemenata dubine 2B=12,0 m i širine 5B=30,0 m, gde je B=6,0 m širina temeljne trake. Mreža konačnih elemenata, prikazana na slici 5.2, se sastoji od 1000 2- D 9-čvornih izoparametarskih elementa. Po bočnim konturama, sprečena su horizontalna a po osnovi i vertikalna i horizontalna pomeranja. Grafička interpretacija dobijenih rezultata numeričke analize urađena je na osnovu računskih uticaja u centru konačnih elemenata. Slika 5.3 predstavlja krivu nosivosti tla koja prikazuje zavisnost između sleganja centrično opterećenog krutog temelja i sile. Razvoj zona plastičnosti je prikazan na slikama 5.4.1 do 5.4.12 za inkrementalni priraštaj sleganja od 0,05 m. Prva zona plastifikacije se pojavljuje pored ivice temeljne trake pri neznatnom nivou sleganja. Analizom formiranja zona plastifikacije može se zaključiti da širenje, oblik i zone plastifikacije prate povećanje sleganja temeljne trake. Slika 5.5 prikazuje dijagrame sleganja slobodne površine a slika 5.6 dijagrame horizontalnih pomeranja duž ivice temelja, za različite nivoe mobilizacije nosivosti peska. Na osnovu upoređenja numeričkog rešenja za graničnu silu NMCfQ 9862 kN= i Terzaghi-jevog rešenja metodom granične ravnoteže TerzaghifQ 10296 kN= , može se konstatovati dobro slaganje. Potrebno je napomenuti da su numeričke vrednosti faktora nosivosti za određivanje granične sile po metodi graničnih stanja plastične ravnoteže određene prema knjizi Stevanović, S.: Fundiranje građevinskih objekata knjiga I, Izgradnja, Beograd, 1999. god., strana 91, tabela II.3. [ ]141 74 75 76 77 78 79 80 81 Primer 5.3.2 Kruta temeljna traka na glini U okviru ovog primera za potrebe primene matematičkog modela i numeričkog postupka u analizi napona i deformacija razmatra se trakasti temelj na glini. Usvojene su računske vrednosti ugla smičuće otpornosti 020φ = i kohezije 2c 15,0kN / m= prema podacima prikazanim u knjizi [ ]87 , strana 67, slika II.56. Na osnovu ugla smičuće otpornosti 020φ = LMC određeni su odgovarajući materijalni parametri NMC: bazni ugao 0B 17.9ϕ = , dilatancija 06.4∆ϕ = i normalni napon 2Np 100 kN / m= . Ostale materijalne konstante, potrebne za definisanje NMC modela, date su na slici 5.7. Diskretizacija domena temeljnog tla izvršena je usvojanjem mreže konačnih elemenata do dubine 2B=6,0 m i širine 4B=12,0 m, gde je B=3,0 m širina temeljne trake. Mreža konačnih elemenata, prikazana na slici 5.7, sastavljena je od 1152 2-D 9-čvorna izoparametarska elementa. Po bočnim konturama, sprečena su horizontalna a po osnovi i vertikalna i horizontalna pomeranja. Slika 5.8 predstavlja krivu nosivosti tla koja prikazuje zavisnost između sleganja centrično opterećene krute temeljne trake i sile. Razvoj zona plastičnosti je prikazan na slikama 5.9.1 do 5.9.8 za inkrementalni priraštaj sleganja od 0,035 m. Slika 5.10 prikazuje dijagrame sleganja slobodne površine a slika 5.11 dijagrame horizontalnih pomeranja duž ivice temelja, za različite nivoe mobilizacije nosivosti glinovitog tla. Na osnovu upoređenja numeričkog rešenja za graničnu silu NMCfQ 1752 kN= i Terzaghi-jevog rešenja metodom granične ravnoteže TerzaghifQ 1674 kN= , može se konstatovati veoma dobro slaganje. 82 83 84 85 86 87 6. KRITERIJUM STABILNOSTI ELASTO-PLASTIČNOG MATERIJALA 6.1 Uvod Problem stabilnosti predstavlja matematički zadatak koji razmatra uslove pod kojima postoji ili ne postoji njegovo rešenje tj., svodi se na problem određivanja granične vrednosti. Analiza stabilnosti geomehaničkih materijala za tačno utvrđenu geometriju, granične uslove i karakter opterećenja zahteva određivanje intervala parametara opterećenja pri kojem naponsko deformacijske veličine imaju konačne vrednosti. U ovom poglavlju je određeno kritično opterećenje posmatranog geomehaničkog materijala uz korišćenje inkrementalno iterativnog postupka za različite putanje napona. Kriva opterećenje deformacija sadrži prevojnu tačku iza koje neposredno dolazi do promene znaka krivine u odnosu na početni nedeformisani položaj. Određivanje kritične tačke moguće je svesti na određivanje onog stanja pri kome dolazi do promene znaka osnovnih kvantiteta metrike površi tečenja. Sama reč stabilnost ima niz različitih značenja, zavisno od problema, stanja ili procesa na koji se odnosi. Velika pažnja se posvećuje fenomenu lokalizacije deformacije u granularnim materijalima kako sa eksperimentalne tako i sa analitičke tačke gledišta. Klasični rezultati za elasto-plastični kontinuum izraženi preko kritične vrednosti modula ojačanja i odgovarajuća orijentacija "shear banda" prikazani su u radovima Hill [ ]52 , Rudnicki i Rice [ ]121 i Runesson [ ]122 . U radovima [ ]114 , [ ]115 definisani su statički i kinematički uslovi za formiranje shear band u granularnom materijalu. Prikazan je matematički precizan kriterijum za određivanje položaja kritične tačke izveden iz generalnog razmatranja svojstvenih vrednosti tangentne matrice krutosti. Ovako postavljen kriterijum za odredivanje kritične tačke je vrlo pogodan sa stanovišta praktičnih primena, mada je potrebno naglasiti da ovako formulisani uslov ima karakter potrebnog ali ne i dovoljnog uslova. Ova metodologija je primenjena od strane velikog broja drugih autora na granularne materijale [ ]15 , [ ]149 sa tom razlikom što su koristili Mohr-Coulomb-ov kriterijum kao funkciju tečenja. Od velikog praktičnog značaja je da se odredi da li će posmatrana plastična deformacija u trenutku kada je dostignuto ravnotežno stanje povećati svoju dužinu za konačan (mali) inkrement i prestati sa rastom (stabilan režim deformisanja) ili će nastaviti sa nekontrolisanim rastom dok ne dođe do velikih priraštaja posmatranih deformacija (nestabilan režim deformisanja). Prelaz fizičkog (mehaničkog) sistema iz stabilnog u nestabilan režim deformisanja često se naziva gubitak stabilnosti. Granica ovog prelaza poznata je kao kritično 88 stanje materijala a njemu odgovarajuće opterećenje je kritično opterećenje. Kritično opterećenje jeste granica grananja formi ravnoteže, kada se posle gubitka stabilnosti prvobitne forme ravnoteže pojavi nova forma kod koje malim priraštajima opterećenja odgovaraju veliki priraštaji odgovarajućih deformacija. 6.2 Konstitutivne jednačine elasto-plastičnog materijala Obično se u teoriji plastičnosti zanemaruju uticaji brzine deformacije i temperaturnog polja. Tako dobijena teorija zove se izotermička teorija plastičnosti. Ograničićemo naša razmatranja na izotermički process (T cons. 0)= > . U odeljku 2.5 su prikazani osnovni pojmovi matematičke teorije plastičnosti i odeljku 3.4 data je numerička integracija konstitutivnih jednačina, a čije delove i osnovne veze, u cilju preglednosti i dalje analize, ponovo navodimo. Pomeranje tačke tela označićemo sa u x( ), gde su x pravougle Dekartove koordinate posmatrane tačke. Usvajajući pretpostavku o malim pomeranjima, tenzor deformacije se može napisati u obliku: ε ij i j j iu u= +12 ( ), , (6.1) gde zarez označava parcijalni izvod po odgovarajućoj koordinati. Konstitutivne jednačine, matematičke plastičnosti izvedene su pod sledećim pretpostavkama: material posmatranog tela je izotropan, veza između napona i elastične deformacije je određena Hukovim zakonom i definisanju veze između priraštaja plastične deformacije i napona u plastičnoj oblasti. Na osnovu poslednje pretpostavke može se napisati sledeća relacija: p P( ,q), 0ε = λ σ λ > (6.2) Trenutni process plastične deformacije zavisi od istorije plastične deformacije, koja se najčešće uzima u obzir preko efektivne plastične deformacije koja predstavlja dužinu puta u prostoru plastičnih deformacija. Bolja mera, posebno za opisivanje višeosnih deformacija, umesto skalara bila bi neki tenzor [ ]89 . Definišu se odgovarajuće veze za iznalaženje unutrašnjih (strukturnih) koordinata sa odgovarajućim koeficijentima koji karakterišu fizičke osobine posmatranog materijala. Brzina priraštaja q može se predstaviti nelinearnom funkcijom h koja je funkcija tenzora napona σ i odgovarajućih unutrašnjih koordinata q , q h( ,q)= λ σ (6.3) Za infinitezimalne deformacije možemo izraziti tenzor deformacije kao zbir tenzora elastičnih i plastičnih deformacija: 89 e pε = ε + ε (6.4) gde indeks e označava elastičan deo, i indeks p označava plastičnu promenu u polju deformacija kada se unutrašnja koordinata q promeni u q q+  . Neka Q označava izvod funkcije tečenja f f q= ( , )σ po tenzoru napona i R odgovarajući izvod po unutrašnjoj koordinati q : f fQ( ,q) , R( ,q) q ∂ ∂σ = σ =∂σ ∂ (6.5) Ako je P normala na funkciju tečenja i predstavlja pravac tečenja u pravcu odgovarajućeg napona σ , onda uslov P Q= ograničava analizu na pridruženi zakon tečenja, a ako je P Q≠ onda je zakon tečenja neasocijativan. Ako funkciju tečenja f razvijemo u Tajlorov red, onda se Tajlorova formula za funkciju ( )f ,qσ može napisati u obliku: n 1 n 1 n 1 n 1 n 1 n 1 n 1 n 1 n 1 n 1 n 1 n 1 n 1 n 1 n 1 n 1 n 1 2 n 1 n 1 n 1 n 1 n 1 f ff f ( ,q ) ( ,q ) : ( ) ( ,q ) (q q ) q 1 f ( ,q ) : ( ) : ( ) 2 f2 ( ,q ) : ( ) (q q ) q f ( ,q ) (q q ) O( ,q ) q q + + + + + + + + + + + + + + + + + + + + + + ∂ ∂= σ + σ σ−σ + σ ⋅ − +∂σ ∂ ∂+ σ σ−σ σ−σ +∂σ∂σ ∂+ σ σ−σ ⋅ − +∂σ∂ ∂ + σ ⋅ − + σ∂ ∂ (6.6) Za dovoljno male vrednosti λ , u okviru ovog razmatranja, smatramo da su kvadratni članovi dovoljni dok članovi višeg reda predstavljaju ostatak n 1 n 1O( ,q )+ +σ koji ćemo zanemariti. Uvodimo sledeće oznake: 2 2 2f f fT( ,q) , U( ,q) ,V( ,q) q q q ∂ ∂ ∂σ = σ = σ =∂σ∂σ ∂σ∂ ∂ ∂ (6.7) gde je: T simetričan tenzor četvrtog reda, U simetričan tenzor drugog i V skalar. Jednačinu (6.6), s obzirom na uvedene oznake (6.7), možemo napisati u obliku: [n 1 n 1 n 1 n 1 n 1 n 1 n 1 n 1 2 n 1 n 1 n 1 n 1 n 1 1f f Q : ( ) R (q q ) T : ( ) : ( ) 2 2U : ( ) (q q ) V (q q ) + + + + + + + + + + + + + = + σ − σ + ⋅ − + σ − σ σ −σ + + σ − σ ⋅ − + ⋅ −  (6.8) Pri prelasku sa jedne na njoj infinitezimalno blisku površ tečenja, iz jednačine (6.8) uz korišćenje izraza (6.2) i (6.3) možemo napisati sledeću kvadratnu jednačinu po priraštaju efektivne plastične deformacije λ : 90 [ ] 21 2 3f a a a 0λ = λ − λ + = (6.9) gde su koeficijenti kvadratne jednačine 1 2a ,a i 3a određeni prema sledećim izrazima: 2 2 1 n 1 n 1 n 1 n 1 n 1 n 1 n 1 n 12a b T : (C : Q ) : (C : Q ) 2b U : C : Q h V h+ + + + + + + += ⋅ + ⋅ ⋅ − ⋅ (6.10.1) a b P C Q R hn n n n2 1 1 1 1= ⋅ −+ + + +: : (6.10.2) 3 n 1a f += (6.10.3) gde je: b R h P C Q R h n n n n n n = − − + + + + + + 1 1 1 1 1 1: : (6.10.4) Postoje dva pozitivna korena diskretnog uslova konsistencije predstavljenog jednačinom (6.9), koja ćemo označiti 1λ i 2λ : 1 2 1(a d ) / 2aλ = − (6.11.1) odnosno 2 2 1(a d) / 2aλ = + (6.11.2) gde je koeficijent d jednak: 2 2 1 3d (a ) 4a a= − ⋅ (6.11.3) Koeficijenti 1 2a ,a i 3a imaju isto značenje kao i u prikazanim izrazima (6.10.1-6.10.4). 6.3 Odredivanje kritične tačke Izraz (6.10.1) može geometrijski da se interpretira kao druga osnovna kvadratna forma površi tečenja f . T U V, , , izrazi određeni relacijama (6.7), su koeficijenti druge osnovne kvadratne forme površi tečenja. Izraz za normalnu krivinu posmatrane površi, s obzirom na jednačine ((6.10.1-6.10.4), možemo napisati u sledećem obliku: ( ) 2 2 n 1 n 1 n 1 n 1 n 1 n 1 n 1 n 1 2 n 1 n 1 n 1 n 1 b T : (C : Q ) : (C : Q ) 2b U : C : Q h V h1 b P : C : Q R h + + + + + + + + + + + + ⋅ + ⋅ ⋅ − ⋅=ρ ⋅ − (6.12) odakle je jasno da znak normalne krivine zavisi samo od znaka brojioca izraza na desnoj strani gornje jednačine. U proučavanju znaka normalne krivine, prema tome, mogu nastupiti tri slučaja: 21. U T V 0.δ = + ⋅ = (6.13.2) 91 Imenilac izraza na desnoj strani uvek je pozitivna veličina, jer predstavlja kvadrat. Tada je brojilac potpuni kvadrat i ne menja znak, tj. krivina opet ostaje stalnog znaka u svim pravcima, ali postoji jedan pravac u kome je krivina jednaka nuli (Slika 6.1). Ovakve tačke površine zovu se paraboličke, a za površinu se kaže da u takvoj tački ima paraboličnu krivinu. 22. U T V 0.δ = + ⋅ < (6.13.1) U tom slučaju brojilac ne menja znak i krivina normalnih preseka u svim pravcima ima isti znak, odn., glavne normale normalnih preseka u svim pravcima kroz datu tačku imaju smer na istu stranu. Takve tačke površine zovu se eliptičke tačke. Za površinu se kaže da u takvoj tački ima eliptičnu krivinu (Slika 6.2). U takvim tačkama površine ne postoji nijedan pravac kome bi odgovarao normalni presek krivine nula. Znak krivine je stalno pozitivan ili stalno negativan prema znaku koeficijenta T . 92 23. U T V 0.δ = + ⋅ ; (6.13.3) Brojilac menja znak, tj. u posmatranoj tački postoje glavne normale suprotnih smerova za normalne preseke raznih pravaca. Normalna krivina u istoj tački može u tom slučaju za razne pravce biti pozitivna i negativna, a postoje dva pravca u kojima je krivina nula (Slika 6.3). Takva tačka površi zovu se hiperbolička, a površina u toj tački ima hiperboličku krivinu. Sva ova razmatranja odnose se uvek na najbližu okolinu posmatrane tačke. Važno je napomenuti da se na istoj površini, u opštem slučaju, mogu nalaziti krivine raznih vrsta. Klasifikacija tačaka jedne površi, odnosno utvrđivanje oblika površi u okolini posmatrane tačke, jasno ističe značaj i druge osnovne kvadratne forme površi. Istaknimo još jedan konkretan geometrijski smisao druge osnovne kvadratne forme površi. Ona određuje meru odstupanja površi u okolini posmatrane tačke od njene tangentne ravni u toj tački. Kompletan nelinearan odgovor jednog elasto-plastičnog sistema formuliše njegovu ravnotežnu putanju, određujući pri tome i odgovarajuće singularne vrednosti. Zadatak stabilnosti je definisan preko oblasti dopustivih deformacija posmatranog elasto- plastičnog materijala. Sračunavanje diskretnih tačaka ravnotežne krive zajedno sa odgovarajućim vrednostima δ , definisane gornjim izrazima, nailazimo i na tačku u kojoj dolazi do promene znaka ove veličine. Sa numeričkog aspekta mogućnost poklapanja neke konfiguracije sa samom tačkom kritičnog stanja u okviru bilo kakvog algoritma gotovo je isključena. Pri vrednosti funkcije 0δ = izraz (6.11) daje λ →∞ , pa relacija 12a 0> postaje i potreban i dovoljan uslov za stabilnost ravnotežnog položaja posmatranog sistema za vreme opterećenja. Na ovaj način posmatranjem ravnotežne putanje sistema, odnosno analizom 93 njegovog ponašanja u okolini kritične tačke direktno se formulišu uslovi za egzistenciju granične tačke. 6.4 Ilustracija predloženog postupka na konkretnom primeru U ovom brojnom primeru razmatrane su različite naponske putanje u triaksijalnom opitu, grafički predstavljene na slici 6.4., za slučaj primene konstitutivnih jednačina HISS modela, koji je opisan u odeljku 2.6.4. Izvode funkcije tečenja, prikazane jednačinom 2.101, po odgovarajućim promenljivim možemo odrediti prema sledećim izrazima: ( ) ( ) ( ) 1 1 a a m 1 m 1 m m n 1I I1 a ijp p ij 2 13 a 2d ij ij2 f p Fs n 2 p mFs FbSrI S mFs FbA − − −− − − ∂  = α − γ δ +  ∂σ + − + (6.14.1) ( )1 a n I( 1) 1 p f a Fs− η+ ∂ = − ηξ∂ξ (6.14.2) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( )( )}{ ( ) 1 a m 1 1 1m a a m 1 1 1 m m m n 2I2 a ij klp ij kl n 1I I1 13 a 2d ij kl ij kl ij kl ij klp p 2 1 13 2d 2d ij kl ij kl ij kl4 f p Fs n n 1 2 p mFs n 2 SrI S S A A mFs FbI SrI 5 3 m 1 Fs Sr S S 2 Fs Sr 1 S A A S m m 1 Fs − − −− −− − − −− − ∂  = − α − γ δ δ +  ∂σ ∂σ    + α − γ ⋅ δ + δ − δ + δ +     + − − + − + −  − − ( )( ) ( ) m 2 m 1 m m 3m 1 m 2 2 13 1 ij kl a 2d ik jl ij kl2 3 3 3 2 2d ik jl ik jl ij kl ij kl2 3 FbA A p mFs FbSrI mFs Fb I S S S S − − − − − − + − δ δ − δ δ +  + β δ + δ − δ + δ  (6.14.3) ( ) ( ) ( )m 11 1 m a a 2 n 1 nI I( 1) 1 13 1 a ij 2d ij ijp p 2 ij f a p n Fs mFs SrI S A −−− η+ − −∂  = − ηξ δ + −  ∂ξ∂σ (6.14.4) ( )1 a n2 I( 2) 1 p2 f Fsa ( 1) − η+∂ = η η+ ξ∂ξ (6.14.5) gde je uvedena sledeća oznaka: 3 2 ij 2d ik kj 2d ij 3 3 2A I S S I 2 3 −  = β⋅ − δ   (6.14.6) 94 Slika 6.4 Putanje napona u triaksijalnoj ravni Pretpostavićemo da nema drugih napona izuzev normalnih, tako da su normalni naponi u ovom slučaju i glavni naponi. Ako sa xxσ označimo aksijalni napon a sa yyσ i zzσ bočne napone onda ćemo odgovarajuće veze izraziti prema ovim oznakama. Rešenje je dobijeno tako što je materijal prvo opterećen hidrostatičkim pritiskom od [ ]0 90 kPaσ = , pri atmosferskom pritisku [ ]ap 101 kPa= . U ovom slučaju odgovarajući izrazi (6.13 do 6.15) se svode na: ( )22 n n n 1 1 1  ηδ = γ + − − η+  (6.15) Asocijativni model je definisan sa sedam materijalnih konstanti, s tim da se može usvojiti konstantna vrednost za materijalni parametar m 0.5= − za većinu frikcionih materijala. Materijalni konstante HISS modela koje se odnose na jednu vrstu peska (Leighton Buzzard sand [ ]32 ) date su u tabeli 6.1. E ν β γ n a1 η 103770.30 0.29 0.442 0.089 3 0.00018 0.85 Tabela 6.1 Materijalni parametri HISS modela Više detalja o materijalnim parametrima i procedurama za njihovo određivanje može se naći u gore navedenoj Referenci. 95 6.4.1 Test standardne triaksijalne kompresije (CTC) Za test standardne triaksijalne kompresije komponente tenzora napona i komponente devijatora tenzora napona mogu se matrično predstaviti u obliku: 0 xx CTC 0 0 0 0 0 0 0 0 σ +σ  σ ⇒ σ  σ  (6.16.1) CTC xx 2 0 0 S 0 1 0 3 0 0 1  σ  ⇒ −  −  (6.16.2) Ako u jednačinu (6.13) unesemo odgovarajuće izraze (6.14) za naponsko stanje definisano izrazom (6.16), jednačinu (6.13) možemo napisati u obliku: ( ) ( )( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ){ ( ) ( ) ( ) ( ) } 2 14 2 2 a 1 0 xx 3 a2 24 a 1 0 xx 0 xx 3 a a2 0 xx a 0 xx a 2 n 12 p a 32 p1 n n 2p a 1 2 13 3 p 3 p1 n 21 3 3 p n 21 32 3 p 1 2k n n n 1 2 4k n n 1 2 2k 2n n 1 4 − η+− − η+− −η ξ σ +σ −β −η η+ ξ −βσ +σ σ +σ −β −−β σ +σ −−β σ +σ δ = − +  + + − α − γ −    − − + − α − γ +    + + − α − γ   (6.17.1) u kojoj su uvedene sledeće oznake: k c c c c c c c c c c = + − − + + + −+ + − + − − +1 2 1 3 1 2 9 2 1 2 2 3 1 5 4 9 1 2 1 2 3 2 2 1 1 3 2 2 2 1 2 1 c h b gn s b g b g b g ν ν ν ν ν( ) (6.17.2) ( ) ( )0 0 xx xx a n 23 3 1 1 p 1 c 1 n 2 −σ σ +σ σ −β  = + α − γ   (6.17.3) ( ) ( ) ( ) ( ) ( ) ( )21 a 0 xx 0 0 xx 31 a xx axx 2n 2 n 2a p 3 3 3 3 2 2 p p 3 11E c 1 n 2η+ −η σ +σ σ σ +σ σ −β−βξ σ  = + α − γ +   (6.17.4) Posle elementarnih transformacija jednačinu (6.17.1) možemo napisati u sledećem obliku: ( ) n 2 2 2 0 xx a 2 13 1 kn n n 1 2 0 1 p 1 2k 3 −  −β  σ +ση + δ = + − α − γ + =    η+ −     (6.18) 96 Rešenje karakteristične jednačine stabilnosti (6.18) je prikazano na slici 6.5.1. Najnepovoljnije ravni u pogledu kritičnog naponskog stanja biće one kod kojih ugao između totalnog napona i normale na ravan ima najveću vrednost. Posmatramo ravan čija normala zaklapa ugao θ sa većim glavnim naponom. Sa slike 6.5.2 je: cr CTC cr CTC xx cr CTC 0 xx 1 arcsin 4 2 2 π σθ = + σ + σ (6.19) 97 6.4.2 Test redukovane triaksijalna ekstenzija (RTE) Komponente tenzora napona i komponente devijatora tenzora napona za RTE test mogu se matrično predstaviti u obliku: 0 xx RTE 0 0 0 0 0 0 0 0 σ −σ  σ ⇒ σ  σ  (6.20.1) RTE xx 2 0 0 S 0 1 0 3 0 0 1 − σ  ⇒     (6.20.2) Ako u jednačinu (6.13) unesemo odgovarajuće izraze (6.14) za naponsko stanje određeno izrazom (6.20), jednačinu (6.13) možemo napisati u obliku: ( ) ( )( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ){ ( ) ( ) ( ) ( ) } 2 14 2 2 a 1 0 xx 3 a2 24 a 1 0 xx 0 xx 3 a a2 0 xx 0 xx a a 2 n 12 p a 32 p1 n n 2p a 1 2 13 3 p 3 p1 n 2 n 21 13 32 3 p 3 p 1 2k n n n 1 2 4k n n 1 2 2k 2n n 1 4 − η+− − η+− −η ξ σ −σ +β −η η+ ξ +βσ −σ σ −σ +β − −+β +βσ −σ σ −σ δ = − +  + + − α − γ −      − − + − α − γ + + − α − γ       (6.21.1) gde su: k c c c c c c c c c c = − − − − + + −− + + + − − +1 2 1 3 1 2 9 2 1 2 2 3 1 5 4 9 1 2 1 2 3 2 2 1 1 3 2 2 2 1 2 1 c h b gn s b g b g b g ν ν ν ν ν( ) (6.21.2) ( ) ( )0 0 xx xx a n 23 3 1 1 p 1 c 1 n 2 −σ σ −σ σ +β  = − α − γ   (6.21.3) ( ) ( ) ( ) ( ) ( ) ( )21 a 0 xx 0 0 xx 31 a xx axx 2n 2 n 2a p 3 3 3 3 22 p p 3 11Ec 1 n 2η+ −η σ −σ σ σ −σσ +β+βξ σ  = − α − γ +   (6.21.4) Posle elementarnih transformacija jednačinu (6.21.1) možemo napisati u sledećem obliku: ( ) n 2 2 2 0 xx a 2 13 1 kn n n 1 2 0 1 p 1 2k 3 −  +β  σ −ση + δ = + − α − γ + =    η+ −     (6.22) 98 Rešenje karakteristične jednačine stabilnosti (6.22) je prikazano na slici 6.6.1. Sa slike 6.6.2 vrednost kritičnog ugla θ je: cr RTE cr RTE xx cr RTE 0 xx 1 arcsin 4 2 2 π σθ = + σ − σ (6.23) 99 6.4.3 Test standardne triaksijalne ekstenzije (CTE) Komponente tenzora napona i komponente devijatora tenzora napona za CTE test mogu se matrično predstaviti u obliku: 0 CTE 0 yy 0 yy 0 0 0 0 0 0  σ σ ⇒ σ +σ  σ +σ  (6.24.1) yyCTE 2 0 0 S 0 1 0 3 0 0 1 − σ  ⇒     (6.24.2) Ako u jednačinu (6.13) unesemo odgovarajuće izraze (6.24) za naponsko stanje definisano izrazom (6.16), jednačinu (6.13) možemo napisati u obliku: ( ) ( )( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) 2 14 2 2 0 yya 1 3 a2 24 0 yy 0 yya 1 3 a a2 0 yy 0 yy a a 2 n 12 3 2p a2 p1 n n 23 2 3 2p a 1 2 1 p 3 p1 n 2 n 23 2 3 21 12 3 p 3 p 1 2k n n n 1 2 4k n n 1 2 2k 2n n 1 4 − η+− − η+− −σ + ση ξ +β −σ + σ σ + ση η+ ξ +β +β − −σ + σ σ + σ+β +β δ = − +  + + − α − γ −      − − + − α − γ + + − α − γ        (6.24.1) gde su: ( ) ( ) 2 5 42 1 1 23 9(1 ) 2 1 2 1 1 23 9 2c c 1 c k 2 c c c − ν +ν − + + −ν =  − − −ν  (6.24.2) ( ) ( )0 yy0yy a n 23 23 11 p 1c 2 n 2−σ + σσσ +β = + α − γ   (6.24.3) ( ) ( ) ( ) ( ) ( ) ( )2 0 yy 0 yy1 a 0 31 a yy ayy 2n 2 n 23 2 3 2a p 3 3 22 p p 3 11Ec 2 n 2η+ −σ + σ σ + ση σσ +β+βξ σ  = + α − γ +   (6.24.4) Posle elementarnih transformacija jednačinu (6.24.1) možemo napisati u sledećem obliku: ( ) n 2 2 0 yy2 a 3 2 2 11 kn n n 1 2 0 1 p 1 2k 3 −σ + σ  +β η + δ = + − α − γ + =    η+ −     (6.25) 100 Rešenje karakteristične jednačine stabilnosti (6.25) je prikazano na slici 6.7.1. Sa slike 6.7.2 vrednost kritičnog ugla θ je: cr CTE yycr CTE cr CTE 0 yy 1 arcsin 4 2 2 σπθ = + σ + σ (6.26) 101 6.4.4 Test redukovane triaksijalne kompresije (RTC) Za test redukovane triaksijalne kompresije komponente tenzora napona i komponente devijatora tenzora napona mogu se matrično predstaviti u obliku: 0 RTC 0 yy 0 yy 0 0 0 0 0 0  σ σ ⇒ σ −σ  σ −σ  (6.27.1) yyRTC 2 0 0 S 0 1 0 3 0 0 1  σ  ⇒ −  −  (6.27.2) Ako u jednačinu (6.13) unesemo odgovarajuće izraze (6.27) za naponsko stanje definisano izrazom (6.16), jednačinu (6.13) možemo napisati u obliku: ( ) ( )( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) 2 14 2 2 0 yya 1 3 a2 24 0 yy 0 yya 1 3 a a2 0 yy 0 yy a a 2 n 12 3 2p a2 p1 n n 23 2 3 2p a 1 2 1 p 3 p1 n 2 n 23 2 3 21 12 3 p 3 p 1 2k n n n 1 2 4k n n 1 2 2k 2n n 1 4 − η+− − η+− −σ − ση ξ −β −σ − σ σ − ση η+ ξ −β −β − −σ − σ σ − σ−β −β δ = − +  + + − α − γ −      − − + − α − γ + + − α − γ        (6.28.1) gde su: ( ) ( ) 2 5 42 1 1 23 9(1 ) 2 1 2 1 1 23 9 2c c 1 c k 2 c c c − ν +ν − + + −ν =  − − −ν  (6.28.2) ( ) ( )0 yy0yy a n 23 23 11 p 1c 2 n 2−σ − σσσ −β = − α − γ   (6.28.3) ( ) ( ) ( ) ( ) ( ) ( )2 0 yy 0 yy1 a 0 31 a yy ayy 2n 2 n 23 3 2a p 3 3 22 p p 3 11Ec 2 n 2η+ −σ +σ σ − ση σσ −β−βξ σ  = − α − γ +   (6.28.4) Posle elementarnih transformacija jednačinu (6.28.1) možemo napisati u sledećem obliku: δ ηη α σ σ γ β= + + − L NM O QP −F HG I KJ − + + − F HG I KJ − = − n n n p k k yy a n 2 0 2 2 1 1 3 2 2 1 1 2 2 1 3 0b g (6.29) 102 Rešenje karakteristične jednačine stabilnosti (6.29) je prikazano na slici 6.8.1. Sa slike 6.8.2 vrednost kritičnog ugla θ je: cr RTC yycr RTC cr RTC 0 yy 1 arcsin 4 2 2 σπθ = + σ − σ (6.30) 103 6.4.5 Test triaksijalne kompresije (TC) Za test triaksijalne kompresije komponente tenzora napona i komponente devijatora tenzora napona mogu se matrično predstaviti u obliku: 0 xx TC 1 0 xx2 1 0 xx2 0 0 0 0 0 0 σ +σ  σ ⇒ σ − σ  σ − σ  (6.31.1) TC xx 2 0 0 S 0 1 0 2 0 0 1  σ  ⇒ −  −  (6.31.2) Ako u jednačinu (6.13) unesemo odgovarajuće izraze (6.14) za naponsko stanje definisano izrazom (6.31), jednačinu (6.13) možemo napisati u obliku: ( ) ( )( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ){ ( ) ( ) ( ) ( ) } 2 14 2 2 a 1 0 3 a2 24 a 1 0 0 3 a a2 0 0 a a 2 n 12 p a 32 p1 n n 2p a 1 2 13 3 p 3 p1 n 2 n 21 13 32 3 p 3 p 1 2k n n n 1 2 4k n n 1 2 2k 2n n 1 4 − η+− − η+− −η ξ σ −β −η η+ ξ −βσ σ −β − −−β −βσ σ δ = − +  + + − α − γ −      − − + − α − γ + + − α − γ       (6.32.1) gde su: ( ) ( ) ( ) ( ) 23 5 4 1 2 1 1 22 4(1 ) 2 1 1 2 1 1 22 (2c 1) 1 c 2c c 1 c k 3(c 1) 2 1 c 2c 2c 1 c − ν +ν  − − + ν + + + −ν   =   + + + ν − + + −ν    (6.32.2) ( )0 0 xx a n 23 3 1 1 p 1 c n 2 −σ σ σ −β  = α − γ   (6.32.3) ( ) ( ) ( ) ( ) ( ) ( )21 a 0 0 0 31 a xx axx 2n 2 n 2a p 3 3 3 3 32 p p 4 11Ec n 2η+ −η σ σ σσ −β−βξ σ  = α − γ +   (6.32.4) Posle elementarnih transformacija jednačinu (6.32.1) možemo napisati u sledećem obliku: ( ) n 2 2 2 0 a 2 13 1 kn n n 1 2 0 1 p 1 2k 3 −  −β  ση + δ = + − α − γ + =    η+ −     (6.33) 104 Rešenje karakteristične jednačine stabilnosti (6.33) je prikazano na slici 6.9.1. Sa slike 6.9.2 vrednost kritičnog ugla θ je: cr TC cr TC xx cr TC 0 xx 1 3arcsin 4 2 4 π σθ = + σ + σ (6.34) 105 6.4.6 Test triaksijalne ekstenzije (TE) Za test triaksijalne ekstenzije komponente tenzora napona i komponente devijatora tenzora napona mogu se matrično predstaviti u obliku: 0 xx TE 1 0 xx2 1 0 xx2 0 0 0 0 0 0 σ −σ  σ ⇒ σ + σ  σ + σ  (6.35.1) TE xx 2 0 0 S 0 1 0 2 0 0 1 − σ  ⇒     (6.35.2) Ako u jednačinu (6.13) unesemo odgovarajuće izraze (6.14) za naponsko stanje definisano izrazom (6.35), jednačinu (6.13) možemo napisati u obliku: ( ) ( )( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ){ ( ) ( ) ( ) ( ) } 2 14 2 2 a 1 0 3 a2 24 a 1 0 0 3 a a2 0 0 a a 2 n 12 p a 32 p1 n n 2p a 1 2 13 3 p 3 p1 n 2 n 21 13 32 3 p 3 p 1 2k n n n 1 2 4k n n 1 2 2k 2n n 1 4 − η+− − η+− −η ξ σ +β −η η+ ξ +βσ σ +β − −+β +βσ σ δ = − +  + + − α − γ −      − − + − α − γ + + − α − γ       (6.36.1)) gde su: ( ) ( ) ( ) ( ) 23 5 4 1 2 1 1 22 4(1 ) 2 1 1 2 1 1 22 (2c 1) 1 c 2c c 1 c k 3(c 1) 2 1 c 2c 2c 1 c − ν +ν  + + + ν − + + −ν   =   − − + ν + + + −ν    (6.36.2) ( )0 0 xx a n 23 3 1 1 p 1 c n 2 −σ σ σ +β  = α − γ   (6.36.3) ( ) ( ) ( ) ( ) ( ) ( )21 a 0 0 0 31 a xx axx 2n 2 n 2a p 3 3 3 3 32 p p 4 11Ec n 2η+ −η σ σ σσ +β+βξ σ  = α − γ +   (6.36.4) Posle elementarnih transformacija jednačinu (6.36.1) možemo napisati u sledećem obliku: ( ) n 2 2 2 0 a 2 13 1 kn n n 1 2 0 1 p 1 2k 3 −  +β  ση + δ = + − α − γ + =    η+ −     (6.37) 106 Rešenje karakteristične jednačine stabilnosti (6.37) je prikazano na slici 6.10.1. Sa slike 6.10.2 vrednost kritičnog ugla θ je: cr TE cr TE xx cr TE 0 xx 1 3arcsin 4 2 4 π σθ = + σ − σ (6.38 107 6.4.7 Test triaksijalnog prostog smicanja (SS) Za test prostog smicanja komponente tenzora napona i komponente devijatora tenzora napona mogu se matrično predstaviti u obliku: 0 xx SS 0 0 xx 0 0 0 0 0 0 σ +σ  σ ⇒ σ  σ −σ  (6.39.1) SS xx 1 0 0 S 0 0 0 0 0 1   ⇒ ∆σ   −  (6.39.2) Ako u jednačinu (6.13) unesemo odgovarajuće izraze (6.14) za naponsko stanje definisano izrazom (6.39), jednačinu (6.13) možemo napisati u obliku: ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ){ ( ) ( ) ( ) ( ) ( ) ( ) } 0 a 0 0 a a 0 a 0 a 2 n 12 2 2 1 32 4 2 1 a 1 p n n 22 3 34 2 a 1 p 3 p n 231 1 1 3 p n 21 32 2 1 3 p 1 k k n p a p a 1 n n 1 2 2k 2k 2kk n n 1 2 k k 2 2n n 1 4 0 −− η+ σ− −− η+ σ σ− −σ −+β σ δ = + − η ξ +  + η η+ ξ + − α − γ −    − − + − + − α − γ +    + + + − α − γ =   (6.40.1) gde su: ( ) ( ) 2 2 1 1 2 1 1 21 1 32 21 1 1 1 2 1 1 21 1 c c 2 c c c 2 ck k 2c c 1 c 2c c 1 c ν ν +ν +ν +ν +ν + + + ν − − + ν= +− + + −ν − + + −ν (6.40.2) ( ) ( ) 2 1 1 1 21 1 3 2 1 1 1 21 2c 2c 1 c k k 2c c 1 c +ν +ν + + + − ν= − − + + −ν (6.40.3) ( ) ( ) ( ) 2 2 2 22 2 1 1 1 2 1 2 21 1 3 2 2 2 1 1 1 2 2 1 21 2c (3 c ) 3(1 )c c 2c (c ) 1 2 (c ) k 6c (1 c ) (5 )c c (c ) 2 1 2 (1 c ) 1 3 c +ν +ν +ν − + + ν − + + −ν − ν +=  + + + ν + + − ν + + − ν  (6.40.4) n 2 0 0 1 xx xx 3 3c n 2 −  σ σ = α − γ σ σ    (6.40.5) ( ) ( ) ( )21 a 0 01 xx axx n 22 n 2a p 3 30 2 pE xx 3c 3 n 2 1η+ −η σ σ σξ σ  σ  = α − γ +    σ  (6.40.6) 108 Posle elementarnih transformacija jednačinu (6.40.1) možemo napisati u sledećem obliku: ( ) ( ) n 2 2 0 1 1 2 a 1 3 2 k k kkn n n 1 2 1 3 0 1 p 3 1 k k −     ση − −δ = + − α − γ + − =   η+ + −       (6.41) Rešenje karakteristične jednačine stabilnosti (6.41) je prikazano na slici 6.11.1. Sa slike 6.11.2 vrednost kritičnog ugla θ je: cr SS cr SS xx 0 1 arcsin 4 2 π σθ = + σ (6.42) 109 7. ZAKLJUČNA RAZMATRANJA U radu je formulisan nelinearni Mohr-Coulombov konstitutivni model u cilju opisivanja ponašanja peska. U NMC modelu ugao smičuće otpornosti tla ϕ nije konstantan u odnosu na normalni napon, već ima hiperboličku funkcionalnu zavisnost. Za opisivanje elasto-plastičnog ponašanja tla, prema ovako definisanoj hiperboličkoj anvelopi loma, potrebno je odrediti šest parametara. Dva parametra opisuju elastične osobine tla: moduo elastičnosti tla E i Poisson-ov koeficijent tla ν , a četiri parametra definišu smičuću otpornost tla: kohezija c , bazni ugao Bϕ , maksimalna dilatancija ∆ϕ i normalni napon Np pri kome maksimalna dilatancija dostiže polovinu svoje vrednosti. NMC model ima mali broj parametara i relativno jednostavnu proceduru za njihovo određivanje. Na osnovu analize rezultata sprovedenih proračuna može se zaključiti da model pokazuje dobro slaganje sa referentnim rezultatima laboratorijskih ispitivanja tla koji su publikovani u stručnoj literaturi. Konstitutivni NMC model je iskorišćen za formulisanje odgovarajućeg numeričkog modela, koji je zasnovan na linearizaciji principa virtuelnog rada u materijalnoj deskripciji. Celokupno izlaganje i rešavanje graničnog problema, koji predstavlja matematički model analizirane pojave, usmereno je ka primeni u MKE, mada predloženi NMC model i postupak integracije konstitutivnih jednačina mogu biti korišćeni i u drugim pristupima nelinearne analize. Osnovne prednosti predloženog numeričkog modela nad metodama granične ravnoteže proračuna nosivosti temelja su: • nisu potrebne pretpostavke za određeno stanje (geometriju) u toku procesa plastičnog deformisanja; • konstitutivni NMC model uključuje i uticaj srednjeg glavnog napona; • omogućava utvrđivanje zona progresivnog loma; • homogenost je ograničena samo na Gausove integracione tačke; • opterećenje se može dodavati postepeno tokom rešavanja nelinearnog problema. Na osnovu analize rezultata sprovedenih proračuna, preko prikazanih dijagrama granične nosivosti i razvoja zona plastičnosti može se konstatovati da postoji kvalitativna i kvantitativna saglasnost sa odgovarajućim teorijskim rešenjima klasične mehanike tla. Mada je model prvenstveno razvijen za krupnozrne materijale, može se na isti način primeniti na sitnozrna tla u dreniranim uslovima kako je to i prikazano u primeru 5.3.2. 110 U šestom poglavlju je određeno kritično opterećenje HISS modela uz korišćenje inkrementalno iterativnog za rešavanje konstitutivnih jednačina. Predloženim postupakom je pokazano da je određivanje kritičnog opterećenja u analizi materijalne stabilnosti elasto- plastičnog modela moguće svesti na određivanje onog stanja pri kome dolazi do promene znaka osnovnih kvantiteta metrike površi tečenja. Upoređivanjem dobijenih numeričkih rezultata sa datim eksperimentalnim rezultatima standardnih triaksijalnih opita (CTC, TC, TE) [ ]32 jasno se uviđa dobra usaglašenost uspostavljenog kriterijuma sa prikazanim ponašanjem peska u triaksijalnom opitu. Na kraju, napomenimo da se predloženi postupak može dalje generalisati i na druge za inženjersku praksu interesantne materijalne modele, kao što su npr. modeli koji obuhvataju ciklično opterećenje materijala u nedreniranim uslovima, modeli anizotropne plastičnosti tla. 111 LITERATURA 1. Argyris, J.H., and Kelsey,S., "Energy Theorems and Structural Analyses", Butterworth, London. (1960) 2. Atkinson J.H. and Bransby P.L., "The mechanics of soils, An introduction to critical state soil mechanics", McGraw-Hill, London (1978) 3. Atkinson J.H. and Farrar D.M., "Stress path tests to measure soil strength parameters for shallow slips", Proceedings of the 11th Int. Conf. SMFE, San Francisco, Vol.2, 983-986, (1985) 4. Atkinson J.H., "The mechanics of soils and foundations", McGraw-Hill Int.,(1993) 5. Besseling J. F., "Non-linear analysis of structures by the finite element method as a supplement to linear analysis", Compo Meth. Appl. Mech. Eng., 3, 173 - 194, (1974) 6. Barton N., "The shear strength of rock and rock joints", International Journal of Rock Mechanics and Mining Sciences, Vol.13, 255-279 (1977) 7. Barton N. and Kjaernsli B., "Shear strength of rockfill", Journal of geotechnical engineering, ASCE, Vol.107, No.GT7, 873-891 (1981) 8. Bathe K. J., "Finite Element Procedures in Engineering Analysis", Prientice-Hall, Inc., Englewood Cliffs, New Jersey, (1982) 9. Besseling J. F., "Non-linear analysis of structures by the finite element method as a supplement to linear analysis", Compo Meth. Appl. Mech. Eng., 3, 173 - 194, (1974) 10. Bishop A.W. and Henkel D.J., "The measurements of soil properties in the triaxial test", Edward Arnold, London, (1962) 11. Bishop A.W, Webb D.L. and Lewin P.I., "Undisturbed samples of London clay from the Ashford Common shaft: Strength-effective stress relationship", Geotechnique, 15, No.1, 1-31, (1965) 12. Bićanić N. P., "Exat evaluation of contact stress state in computational elasto-plasticity", Eng. Comput., 6, 67 - 73, (1989) 13. Bolton M.D., "The strength and dilatancy of sands", Geotechnique, 36, No.1, 1986, 65-78, (1986) 14. Bolton M.D. and Lau C.K., "Scale effects in the bearing capacity of granular soils", XII ICSMFE, Vol.2, Rio de Janeiro, 895-898, (1989) 15. Borst, de R., "Integration of plasticity equations for singular yield functions", Computers and Structures, Vol. 26, pp. 823 - 829, (1987 ) 16. Borst, de R., "Bifurcations in Finite Elements Models with a Non-associated Flow Law", International Journal for Numerical and Analytical Methods in Geomechanics, 12(1), 99-117, (1988) 17. D. Božinović, Geotehničke osobine terena u područiju Beograda koje izgrađuju koherentni gornjemiocenski sedimenti, Izgradnja, Beograd, (1979) 18. Charles J.A. & Watts K.S., "The influence of confining pressure on the shear strength of compacted rockfill"; Geotechnique, 30 (4): 353-367, (1980) 19. Chen, W.F., "Limit analysis and soil plasticity", Elsevier, New York, (1975) 20. Christian J. T., A. J. Hagmann A. J. and Marr W.A: "Incremental Plasiticity Analysis of Frictional Soils", International Journal for Numerical and Analytical Methods in Geomechanics, Vol. 1 343-375, (1977) 21. Christian J. T., and Desai C.S., "Constitutive Laws for Geologic Media", in Num. Meth. in Geotech. Engin., Christian J. T., and Desai C.S. (Eds.), McGraw-Hill Book Company, New York, Chap. 2, (1977) 22. Clough R.W., "The finite element in plane stress analysis", Proceedings 2nd Conf. on Electronic Computation, Pittsburg, (1960) 23. Coulomb C.A., "Essai sur une application des regles des maximis et minimis a quelues problemes de statique relatifs a l'architecture", Mem. Acad. Roy. des Sciences, Paris, 3, 343-382, (1776) 24. Cook R. D., "Concepts and Applications of Finite Element Analysis", J.Wiley and Sons, (1974) 112 25. Crisfield M. A., "Plasticity computations using the Mohr-Coulomb yield criterion", Eng. Comput. 16. 4., December (1987) 26. Crisfield M. A., "A fast incremental/iterative solution procedure that handles snap through", Computers and Structures, Vol. 13, pp. 55 - 62, (1981 ) 27. Crisfield M. A., "Finite elements and solution procedures for structural analysis", Pineridge Press, Swansea, U.K. (1986) 28. Crisfield M. A., "A faster modified Newton-Raphson iteration", Comp. Meth. Appl. Mech. Eng., 20, 267 - 278, (1979) 29. Damjanić F. 8. Owen O. R. , "Implicit time integration of elasto-vicoplastic solids subject to the Mohr- Coulomb yield criterion", Int. J.Numer. Methods Eng., Vol. 18, 1873 - 1881, (1982) 30. Desai C. S. & Abel J. F., "Introduction to the Finite Element Method", Van Nostrand, 1972. 31. Desai, C.S., and Siriwardine, H.J., "Constitutive laws for engineering materials", Prentice Hall, Englewood- Cliffs, New Jersey (1984) 32. Desai, C.S., "Notes for Advanced School-Numerical Methods in Geomechanics Including Modelling", Udine, Italy, July 10-14, (1989) 33. Desai, C.S, Somasundaram, S. and G. Frantziskonis, G., "A hierarchical approach for constitutive modeling of geological materials", International Journal for Numerical and Analytical Methods in Geomechanics, 10(3), 225-257, (1987) 34. Dinkler, D., and Kroplin, B., "Global constraints applied to static instability analysis", Third Int. Conf. Numerical Methods for Non-linear Problems, 1, 41-60, Dubrovnik (1986) 35. DiMaggio, F.L., and Sandler, I.S.,"Material Models for Granular Soils", J.Eng.Mech. Div., ASCE ,97, 935- 950, (1971) 36. Drucker, D.C. , "Some Implications of Work Hardening and Ideal Plasticity", Quarterly of Applied Math., 7(4), 411-418, (1950) 37. Drucker, D.C. and Prager, W.,"Soil Mechanics and Plastic Analaysis of Limit Design", Quart.Appl. Mathematics,10,157-165(1952) 38. Drucker, D.C., "Stress-Strain Relations in the plastic Range of Metals-Experiments and Basic Concepts", Rheology: Theory and Applications,(Ed. F.R. Eirich), 1, 97-119, Academic Press, New York,(1956) 39. Drucker, D.C. ,"A Definition of Stable Inelastic Material", ASME Journal of Applied Mechanics, 26, 106- 112(1959) 40. Duncan J.M. and Chang C.Y., "Nonlinear analysis of stress and strain in soils", Journal of Soil Mech. and Found. Engineering, ASCE, Vol. 96, No. SM5, 1629-1653, (1970) 41. Duncan J.M. and Buchigani A.L., "An engineering manual for settlement studies", Geot. Eng. Report, Dept. of Civil Eng., Univ. of California at Berkeley, (1976) 42. Dunica Š., Kolundžija B., "Nelinearna analiza konstrukcija", Građevinski fakultet Univerziteta u Beogradu, Naučna knjiga, Beograd, (1986) 43. Fritz P., Numerical solution of rheological problems in rock", Proceedings Int. Symp. Numerical Models in Geomechanics, Zurich, Sept., (1982) 44. Fung, Y.C., "Foundations of Solid Mechanics", Prentice Hall, Englewood-Cliffs, New Jersey(1965) 45. Green A.E. and Atkins J.E., "Large Elastic Deformations and Nonlinear Continuum Mechanics", Clarendon Press, Oxford (1960) 46. Goodman R.E., "Introduction to Rock Mechanics", John Wiley & Sons, New York (1980) 47. Gurtin E.M., "An Introduction to Continuum Mechanics", Academic Press, New York (1981) 48. Haar M.E., "Foundations of theoretical soil mechanics", McGraw-Hill, New York, (1966) 49. Hajdin N., "Uvod u teoriju elastičnosti", Predavanja na III stepenu studija na konstruktivnom odseku Građevinskog fakulteta u Beogradu, Beograd, (1965) 113 50. Hajdin N., "Teorija površinskih nosača", Građevinski fakultet Univerziteta u Beogradu, Zavod za grafičku tehniku Tehnološko metalurškog fakulteta Beograd, (1984) 51. Hansen J.B., "A general formula for bearing capacity", Buletin No.11, Danish Technical Institute, Copenhagen, (1961). 52. Hill K., "The mathematical theory of plasticity", Oxford University Press, Oxford. (1950) 53. Hinton E. & Owen D. R. J., "Finite Element Programming", Academic Press, London, (1977) 54. Hinton E. & Owen D. R. J., "An Introduction to Finite Element Computations", Pineridge Press, Swansea (1979) 55. Hoeg K., Christian J. T. and Whitman R.V., "Settlement of Strip Load on Elastic-Plastic Soil", Journal of the Soil Mechanics and Foundation Division, ASCE, Vol. 94, No. SM2, March., 431-445, (1968) 56. Hoek E. & Brown T.E., "Empirical strength criterion for rock masses"; Jnter. J. Geotech. Eng. Div., ASCE, 106 (9): 1013, (1980) 57. Hofstetter, G., Simo J.C., Taylor R.L., "A modified cap model: closest point solution algorithms", Comput. & Structure, vol 48, 203-214, (1993) 58. Hughes, T.J.R., "Numerical implementation of constitutive models: Rate indipendent deviatoric plasticity", Workshop on Theor. Foundat. for large-scale Comput. of Nonlinear Material Behaviour, North-Western University, Evanston, IL (1983) 59. Качанов Л.М., "Основы теории пластичности", М.Наука, 420 с. (1969) 60. Kojić M., "Primenjena teorija plastičnosti", Biro za građevinarstvo, Beograd (1979) 61. Kojić, M., "A general concept of implicit integration of constitutive relations for inelastic material deformation", (in serbian), Cen. for scien. research SANU, Kragujevac (1993) 62. Kojić, M., Computational Procedures in Inelastic Analysis of Solids and Structures, Center for Scientific Research of Serbian Acadamy of Science and Arts and University of Kragujevac, Kragujevac, (1997) 63. Kondner R.L., "Hyperbolic stress-strain response: cohesive soils", Journal of SMFE, ASCE, Vol.89, No.SM1, February, 115-143, (1963) 64. Kondner R.L. and Zelasko J.S., "A hyperbolic stress-strain formulation for sands", 2nd. Pan -American Conf., SMFE, Brasil, (1963) 65. Krieg, R.D. and Krieg, D.B., "Accuracies of numerical solution methods for the elastic-perfectly plastic model", ASME J. Pressure Vessel Tech., 99, 510-515 (1977) 66. Lade P., Kim M., "Single hardening constitutive model for soil, rock and concrete", Int. Jour. Solids Structures, Vol 32 (Nu 14), 1963-1978, (1995) 67. Lade P.V., "Elastoplastic Stress - Strain Theory for Cohesionless Soils with Curved Yield Surface", Int. Jour. Solids Struct., ASCE, Vol 13, 1019-1035, (1977) 68. Lade P.V., Duncan J.M, "Elastoplastic Stress - Strain Theory for Cohesionless Soil", Jour. Geotech. Engng Div., ASCE, Vol.101, 1037-1053, (1975) 69. Lade, P.V., "Effects of Voids and Volume Changes on the Behaviour of Frictional Materials", Inter Jurn. for Numerical and Analytical Math. in Geomechanics, Vol 12 (351-370), (1988) 70. Lambe T.W. and Whitman,R.V., "Soil Mechanics", John Wiley & Sons, New York, (1969) 71. Lazović M., "Prllog nelinearnoj analizi aksijalno opterećenog šipa", Doktorska disertacija, Građevinski Fakultet Univerziteta u Beogradu, (1988) 72. Lazović M. i Vukićević M., "Konstitutivni modeli tla. I deo, Elastični modeli", G.F., Beograd, (1997). 73. Lelović S., "Primena Mohr-Coulomb-ovog modela za rešavanje nekih problema u geotehnici", Magistarski rad, Građevinski Fakultet Univerzitet u Beogradu, 1990 74. Lovel C.W. and J. M. Johnson J.M., "Shearing behavior of compacted clay after saturation" in the Laboratory Shear Strength of Soil, ASTM Special Technical Publication 740, Ed. R. N. Young and F. C. Townsend, Chicago IL, p. 277., (25 June 1980 ) 114 75. Maksimović M., "Nonlinear failure envelope for soils", Journal of Geotechnical Engineering, American Society of Civil Eng.,Vol.115 No.4 April, 581-586, (1989) 76. Maksimović M., "On the residual shearing strength of clays", Geotechnique 39, No.2, 347-351, (1989) 77. Maksimović M., "New description of the shear strength for rock joints", Rock Mechanics and Rock Engineering, 25(4), Springer Verlag, 275-284, (1992) 78. Maksimović M., "The shear strength components of a rough rock joint", International Journal of Rock Mechanics and Mining Sciences, Pergamon, Vol.33, No.8, 769-783, (1996) 79. Maksimović M. i Santrač, P., "Zbirka zadataka iz osnova mehanike tla", Peto izdanje, (2001) 80. Maksimović M., "Mehanika Tla" AGM Books 4th Ed. (in Serbian) p. 195., (2008) 81. Marachi N.D., Duncan J.M., Chan C.K. and Seed H.B., “Plane strain testing of sand” in the Laboratory Shear Strength of Soil, ASTM Special Technical Publication 740, Ed. R. N. Young and F. C. Townsend, Chicago IL, p. 294., (25 June 1980 ) 82. Marcal P. V. & King. 1. P., "Elasto-plastic analysis of twodimensional stress sistem by the finite element method", Int. J. Mech., Sci., 9, 143 - 155, 1967. 83. Marques J. M. M. C., "Stress computation in elastoplasticity", Eng. Comput., 1, 42 - 51, (1989) 84. Meyerhof G,G., "Penetration tests and bearing capacity of cohesionlee soils", Journal of the Soil Mechanics and Foundation Division, ASCE, Vol. 82, No. SM1, Jan., 1-19, (1956) 85. Meyerhof G,G., "Some recent research on the bearing capacity of foundations", Canadian Geotechnical Journal, No.1, 16-26, (1963) 86. Meyerhof G.G., "Shallow foundations", Journal Soil Mechanics and Foundation Eng, ASCE, Vol 91, No.SM2, 21-31, (1965) 87. Milović D., "Problemi fundiranja na lesnom tlu", Vojvođanska akademija nauka, Novi Sad, (1987) 88. Mitchell J.K., "Fundamentals of soil behaviour", John Wiley & Sons, (1976) 89. Mićunović M., "Primenjena mehanika kontinuuma", Naučna knjiga, Beograd (1980) 90. Mohr O., "Abhandlungen aus dem Gebiete der Technischen Mechanik", 2nd Ed., W.Ernst. Berlin, (1914) 91. Mrozz Z., Norris, V.A, Zienkiewicz O.C., "Application of on anistropic hardening model in the analysis of elasto - plastic deformation of soils", Geotechnique 29, No 1, 1-34, (1979) 92. Nayak G. C. & Zlenkiewicz, O. C., "Elasto-plastlc Stress Analysis", International Jurnal for Numerical Methods in Engineering, 5, 113 - 135, (1972) 93. Nayak G. C. & Zienkiewlcz, O. C. "Convenient forms of stress invariants for plasticity", Proceedings Am. Soc. Civ. Eng ., Proc ASCE, Vol. 98, 949 954, (1972) 94. Nearlović-Veljković N., "Pisana predavanja iz predmeta Mehanika kontinuuma", Poslediplomska nastava na konstruktivnom odseku Građevinskog fakulteta u Beogradu, (1986) 95. Nearlović-Veljković N., Plavšić M., "Teorija elastičnosti", Naučna knjiga, (1980) 96. Nonveiller E., "Mehanika tla i temeljenje", III izdanje, Školska knjiga, Zagreb, (1990) 97. Novožilov V.V., "Teorija uprugosti", GIZ Sudprom, Leningrad, (1958) 98. Ortiz, M., and Popov, E.P.,"Accuracy and Stability of Integration Algorithms for Elastoplastic Constitutive Relations", International Jurnal for Numerical Methods in Engineering, 21, 1561-1576 (1985) 99. Ortiz, M., and Simo, J.C., "An Analysis of a New Class of Integration Algorithms for Elastoplastic Constitutive Relations", International Jurnal for Numerical Methods in Engineering, 23, 353-366 (1986) 100. Owen D. R. J. & Hinton E., "Finite Elements in Plasticity", Pineridge Press Ltd ., U.K, (1980) 101. Paul B., "Modification of the Coulomb-Mohr theory of fractive", Journal of Applied Mechanics, 28, 256 - 268, (1961) 115 102. Pilkey W. ,Saczalski K. and Schaeffer H. (eds). "Structural Mechanics Computer Programs", University Press of Virginia, Charlottesville, (1974) 103. Potts D.M. and Zdravković L., "Finite Elements Analysis in Geotechnical Engineering - Theory", Thomas Telford Publishing, London, (1999) 104. Potts D.M. and Zdravković L., "Finite Elements Analysis in Geotechnical Engineering - Application", Thomas Telford Publishing, pp. 175-177, London, (1999). 105. Prandtl L., "Uber dei Harte plastischer Korper", Nachr. d. Ges. d. Wiss., Gottingen, (1920) 106. Prevost, J.H., "Plasticity Theory for Soil Stress-Strain Behaviour", Journal Engineering Mechanics Division, ASCE, 104, 1117-1194 (1978) 107. Prevost J., Popescu R, "Constitutive relations for soil materials", El. Jour. Geot. Enginer., 96-09, (1996) 108. Prevost J.H. "Two-surface vs. multi surface plasticity theories", International Journal for Numerical and Analytical Methods in Geomechanics, (6), 323-338, (1982) 109. Prevost J.H., "Modeling the Behaviour of Geomaterials", in Geotechnical Modeling and Applications, Ed. Sayed B.S.M., Gulf Publishing Co., 1987, p. 8-76, (1987) 110. Radenković D. i Ivković M., "O primeni teorije plastičnosti na određivanje nosivosti tla", Sedmi skup odeljenja tehničkih nauka, Beograd, 141-151, (1957) 111. Radenković G., "Konačne elastoplastične deformacije", Građev. fakultet Univerziteta u Beogradu, (1996) 112. Rankine W.J.M., "On the stability of loose earth", Phil.Trans.Roy.Soc.,London, 147, 1, 9-27, (1857) 113. Reissner H., "Zum Erddruckproblem", Proceedings 1 st Int. Congress for Applied Mechanics, eds. Biezeno and Burgers, Delft, 295-311, (1924) 114. Rice, J.R., "The Localization of Plastic Deformation", Proceedings 14th Int. Cong. Theoretical and Appl. Mech., Delft (Ed. W.T. Koiter), 1, 207-220, Delft, North-Holland, Amsterdam (1976). 115. Rice, J.R., "Plasticity in Soil Mechanics", Proceedings of the Symposium on the Role of Plasticity in Soil Mechanics., Ed. by A.C. Palmer, Depart. of Engineering, University of Cambridge, 263-280, (1973) 116. Roscoe K.H., Schofield A,N, & Wroth, C.P., "On the yielding of soils", Geotechnique, 8,1, 22-52, (1958) 117. Roscoe K.H. and Burland, J. B., "On the generalised stress-strain behaviour of an idealised wet clay", In Heineman & Leckie (eds.) Engineering plasticity, Cambridge University Press, 535-609, (1968) 118. Rowe P.W., "The stress-dilatancy relation for static equilibrium of an assembly of particles in contact", Proceedings Roy. Soc. A.269, (1962) 119. Rowe P.W., (1969). "The relation betwee the shear strength of sands in triaxial compression, plane strain and direct shear", Geotechnique 19, No.1, 75-86, (1969) 120. Rowe P.W., (1971). "Theoretical meaning and observed values of deformation parameters for soil, Stress-strain behavior of soils", Roscoe memorial symposium, 143-194, (1971) 121. Rudnicki, J.W., and Rice, J.R.,"Conditions for the Localization of Deformation in Pressure-Sensitive Dilatant Materials", J.Mech.Phys.Solids, 23, 371-394(1975) 122. Runesson, K., Ottosen, N.S., and Peric, D., "Discontinuous bifurcations of elastic-plastic solutions at plane stress and plane strain", Int. J. Plast., 7(1/2), 99-101, (1991) 123. Sandler, I.S., and Rubin, D., "An Algorithm and a modular subroutine for the cap model", Int. J. Numer.Anal.Meth.Geomech., 3, 173-186, (1979) 124. Savić, Lj.,"Contribution to the numerical investigation of some static problems of nonlinear structural behaviour",(in serbian), Doct.Thesis, Civil Engineering Faculty, Universite Belgrade (1990) 125. Schofield A.N. and Wroth, C.P., "Critical state soil mechanics", McGraw-Hill Book Co., London, (1968) 126. Schreyer, H.L., Kulak, R.F. and Kramer, M.M., "Accurate numerical solution for elasto-plastic model", Pressure Vessel Tech., ASME ,101, 226-234 (1979) 127. Scott R.F., "Principles of soil mechanics", Addison-Wesley Publ.Co., (1963) 116 128. Segerlind L. J., "Applied Finite Element Analysis", John Wiley Sons (1976) 129. Sekulović M., "Metod konačnih elemenata", Gradevinska knjiga, Beograd, (1984) 130. Simo, J.C., and Ortiz, M., "A Unified Approach to Finite Deformation Elastoplastic Analysis Based on the Use of Hyperelastic Constitutive Equations", Comp.Meth.Appl.Mech.Eng., 49, 221-245(1984) 131. Simo, J.C., and Govindjee, S., "Exact closed-form solution of the return mapping algorithm in plane stress elasto-viscoplasticity", Eng.Comput., 5, 254-258 (1988) 132. Simo, J.C., Ju, J.W., Pister, K.S., and Taylor, R.L., "Assessment of Cap Model: Consistent Return Algorithms and Rate-Dependent Extension", J.Eng.Mech., ASCE, 114,191-218 (1988) 133. Simo, J.C.,Kennedy J.G., and Govindjee, S., "Non-smooth multisurface plasticity and viscoplasticity. Loading/unloading conditions and numer. algorithms", Int. J.Numer. Methods Eng., 26, 2161-2186 (1988) 134. Simo J.C., Taylor R.L., "Consistent tangent operators for rate independent elasto-plasticity", Comp. Meth. Appl. Mech. Engng, Vol 48, 101-118, (1985) 135. Simo J.C., Taylor R.L., "A return mapping algorithm for plane stress elastoplasicity", International Jurnal for Numerical Methods in Engineering, Vol 22, 649-670, (1986) 136. Smith I.M., "Programing the Finite Element Method with aplications to geomechanics", John Wiley & Sons International, (1982) 137. Skempton A.W., "The bearing capacity of clays", Proceedings Building Research Congress, Vol.1¸(1951) 138. Sokolovski V.V., "Statics of granular media", Pergamon Press, Oxford, (1965) 139. Spasojević A., "Analiza stacionarnog širenja sfernih i cilindričnih šupljina u krupnozrnom tlu", Doktorska disertacija, Građevinski Fakultet Univerziteta u Beogradu, (1998) 140. Spenser, A.J.M., "Theory of Invariants", Continuum Physics, (Ed. A.C. Ernigen), 1, 239-353, Academic Press, New York (1972) 141. Stevanović S., "Fundiranje građevinskih objekata", Knjiga I, Građevinski Fakultet u Beogradu, Izgradnja, Beograd, (1999). 142. Stojadinović R., "Mehanika tla", Naučna knjiga, Beograd, (1986) 143. Stojanović R., "Uvod u Nelinearnu Mehaniku Kontinuuma", Zavod za izdav. udžbenika SR Srbije, (1965) 144. Taylor D.W., "Fundamentals of soil mechanics", John Wiley and Sons, New York, (1948) 145. Terzaghi K., "Theoretical soil mechanics", John Wiley and Sons, New York, (1943). 146. Terzaghi K. and Peck R.B., Soil mechanics in engineering practice", 2nd ed. ,John Wiley and Sons, New York, (1967) 147. Valanis, K.C., "Banding and Stability in Plastic Materials", Acata Mechanica, 79, 113-141, (1988) 148. Vermeer P.A., "A double hardening model for sand", Geotechnique 28 (4), 413-433, (1978) Vardoulakis I., "Shear band Inclination and Shear Modulus of Sand in Biaxial Tests," Int. J. Numer. Anal. Meth. Geomech.,4,2,pp103-121,(1980) 149. Vesić A.S., "Bearing capacity of deep foundation in sand", National Academy of Sciences, National Research Council, Highway Research Record, 39, 112-153, (1963) 150. Vesić A.S., "Ultimate loads and settlements of deep foundations in sand", Soil Mechanics Series No.5, Duke University, (1984) 151. Vesić A.S. and Clough G,W., "Behaviour of granular materials under high stresses", Proceedings ASCE, 94, No.SM3, 661-688, (1968) 152. Vesić A.S., "Analysis of ultimate loads of shallow foundations", Journal of Soil Mech and Found. Div. ASCE, 99, No.SM1, 45-73, (1973) 153. Vesić A.S., "Bearing capacity of shallow foundations", Chapter 3 in Winterkorn,H.F and Fang, H.Y. ed., Van Nostrand Reinhold, New York, 121-147, (1984) 154. Vujičić Č., "Fundiranje 1", Građevinski fakultet Beograd, Naučna knjiga, Beograd, (1984) 117 155. Vujić S., Berković M., Kuzmanović D., Milanović P., Sedmak A. i Mićić M., "Primena metode konačnih elemenata kod geostatičkih proračuna u rudarstvu", Rudarsko–Geološki fakultet u Beogradu, (1990) 156. Vukelić S., "Analiza mehaničkih stanja očvrslog betona", Doktorska disertacija, Građevinski Fakultet u Beogradu, (1981) 157. Vukićević M., "Primena elastoplastičnih modela za tlo u proračunu fleksibilnih potpornih konstrukcija", Doktorska disertacija, Građevinski Fakultet u Beogradu, (2006) 158. Wilkins, M.L., "Calculation of elastic-plastic flow," Methods of Computational Phisics, (Eds. B. Alder et al.), 3, Academic Press, (1964) 159. Wood D.M., "Soil behaviour and critical state soil mechanics", Cambridge University Press., (1990) 160. Zienkiewicz O.C. and Naylor D.J., "The adaptation of critical state soil mechanics theory for use in finite elements", Proceedings of the Roscoe Memorial Symposium, Stress Strain Behaviour of Soils, Cambridge, 537-547, (1971) 161. Zienkiewicz O.C., "The finite element method in engineering mechanics", 2nd Ed, John Wiley & Sons International, (1972) 162. Zienkiewicz O. C. & Nayak G. C., "Note on the 'alpha' -constant stiffness method for analysis of non-linear problems", International Jurnal for Numerical Methods in Engineering, 4, 579 - 582, (1972) 163. Zienkiewicz O. C., Humpheson C. & Lewis K. W., "Associated and non-associated visco-plasticity in soil mechanics", Geotechnique, 25, 671 - 685, (1975) 164. Zienkiewicz O. C., Valliappan S. & King I. P. "Elasto-plastic solutions of engineering problems. Initial- stress, finite element approach", Internat. Jurnal for Numerical Methods in Engineering, 1, 75 – l00, (1969) 165. Čabarkapa Ž., "Koeficijent Nγ u teoriji nosivosti plitkih temelja na pesku", Magistarski rad, Građevinski Fakultet Univerziteta u Beogradu, Beograd, (1996) 166. Ćorić S., "Geostatički proračuni", Rudarsko – Geološki fakultet Univerziteta u Beogradu, Beograd, (2006) Biografija Selimir Lelović rođen je 2.10.1961. god. u Kraljevu. Osnovnu školu i Gimnaziju (prirodno- matematičkog smera) završio je u Kraljevu. Po završenoj službi u vojsci upisao se, na Građevinski fakultet Univerziteta u Beogradu, školske 1981/82. god. Diplomirao je na Građevinskom fakultetu u Beogradu 15. septembra 1986. god., na Odseku za Konstrukcije (smer za betonske konstrukcije). Tokom redovnih studija postigao je srednju ocenu 8.19 (osam i 19/100) a diplomski zadatak iz predmeta Fundiranje odbranio je sa ocenom 10 (deset). Od 26.10.1986. god. do 2.02.1987. god. radio je u RO Izgradnja, OOUR Visokogradnja, SOUR Energoprojekt. Od 2.02.1987. god. do 15.01.1989. god. u RO Energodata, SOUR Energoprojekt. Na Građevinski fakultet Univerziteta u Beogradu prelazi 1989. god. u zvanju asistenta pripravnika na predmetu Fundiranje. Juna 1991. god. izabran je za asistenta na predmetu Fundiranje. Novembra 1986. god. upisao je poslediplomske studije na Građevinskom fakultetu u Beogradu, smer za Konstrukcije, Odsek za Tehničku mehaniku i teoriju konstrukcija. Ispite je položio, kao i sve obaveze predviđene programom za poslediplomske studije izvršio je 27.06.1989. god.. Poslediplomske studije završio je sa srednjom ocenom 9,70 (devet i 70/100), a magistarski rad sa temom "Primena Mohr- Coulomb-ovog modela za rešavanje nekih geotehničkih problema" odbranio je 4.02.1991. god. Položio je stručni ispit. U okviru stručne delatnosti bavi se projektima fundiranja inženjerskih objekata, projektovanjem zaštite temeljnih jama, sanacijama klizišta, projektima dinamički opterećenih temelja, različitim problemima geotehnike i numeričkim metodama u geotehici. Član je Društva za fundiranje i mehaniku tla Srbije, Međunarodnog društva za mehaniku tla i geotehniku (ISSMGE – International Society for Soil Mechanics and Geotehnical Engineering) i Društva građevinskih konstruktera Srbije. Član je republičke revizione komisije od 2000. godine. Obavlja funkciju sekretara Katedre za Građevinsku geotehniku na Građevinskom fakultetu.