TEMPERATURNA ANALIZA POSTOPNE GRADNJE MASIVNIH BETONSKIH KONSTRUKCIJ THERMAL ANALYSIS OF SUCCESSIVE CONSTRUCTION OF MASS CONCRETE STRUCTURES Anka Ilc, univ. dipl. inž. grad. Znanstveni članek Primorje, d. d., Vipavska cesta 3, 5270 Ajdovščina UDK: 624.012.4 Gregor Trtnik, univ. dipl. inž. grad. IGMAT, d. d., Polje 351 c, 1000 Ljubljana izr. prof. dr. Igor Planinc, univ. dipl. inž. grad. Univerza v Ljubljani, Fakulteta za gradbeništvo in geodezijo, Jamova 2, 1000 Ljubljana prof. dr. Goran Turk, univ. dipl. inž. grad. Univerza v Ljubljani, Fakulteta za gradbeništvo in geodezijo, Jamova 2, 1000 Ljubljana Povzetek l V članku je prikazan nov numerični postopek za določanje temperatur strjujočega se betona pri masivnih betonskih konstrukcijah. Model je zasnovan na osnovi metode končnih elementov in omogoča relativno preprosto modeliranje postopne gradnje. Dobro ujemanje med numeričnimi in eksperimentalnimi rezultati temperatur strjujočega se betona pri masivnem armiranobetonskem vodnjaku dokazuje veliko natančnost in s tem primernost predstavljenega numeričnega postopka za analizo postopne gradnje vseh vrst masivnih betonskih konstrukcij. Summary l A new numerical procedure for the prediction of temperature development in early age mass concrete structures is briefly presented. This model is implemented in the finite element code and enables modelling of successive conctruction. In order to analyze the adequacy of the model, temperature profile was both measured on site and calculated with the proposed numerical model. Calculated results are in good agreement with the experimentally determined temperature profile. Therefore, the presented numerical procedure can be effectively used to model the temperature development in mass concrete structures. 1«UVOD Sproščanje hidratacijske toplote med strjevanjem betona pri gradnji večine gradbenih konstrukcij ne povzroča večjih težav, nasprotno pa lahko pri masivnih betonskih konstrukcijah zaradi oteženega odvajanja sproščene hidratacijske toplote povzroča pregrevanje betona in posledično nezaželene poškodbe. Tehnološki postopki gradnje masivnih betonskih konstrukcij zato pogosto zahtevajo postopno gradnjo z manjšimi betonskimi bloki oziroma sloji, kar tovrstne poškodbe prepreči. Zaradi pomembnosti masivnih betonskih konstrukcij v gradbeništvu zasledimo v literaturi številne raziskave, ki se ukvarjajo s temperaturno analizo tovrstnih konstrukcij in s tem povezanimi pojavi. Raziskave lahko razdelimo v dve glavni skupini. V prvo uvrščamo raziskave o kemijskih in fizikalnih lastnostih betona med strjevanjem ([Hakarai, 2007], [Ishida, 2007], [Scherer, 2007], [Snelson, 2008]), drugo skupino pa predstavljajo numerični postopki za temperaturno analizo betona med strjevanjem oziroma med gradnjo. Z razvojem računalništva dobivajo te raziskave vse večji pomen. Glede na uporabljeno numerično metodo za reševanje osnovne enačbe prevajanja toplote v trdni snovi ločimo numerične metode, zasnovane na diferenčni metodi ([Turk, 1991], [Ballim, 2004]) oziroma na metodi končnih elementov ([Turk, 1991], [Tanabe, 1985], [Hue, 2000], [Faria, 2006]). Glavna pomanjkljivost teh raziskav je v tem, da razviti numerični postopki ne omogočajo določitve temperatur strjujočega se betona pri postopni gradnji masivnih betonskih konstrukcij, kar je ena od pomembnih značilnosti gradnje tovrstnih konstrukcij. V članku prikazujemo nov numerični postopek za določitev temperatur strjujočega se betona pri postopni gradnji masivnih betonskih konstrukcij. Novost prikazanega postopka, ki je zasnovan na metodi končnih elementov, je v postopnem vključevanju novih končnih elementov v globalno mrežo. S primerjavo med numeričnimi in eksperimentalnimi rezultati temperatur strjujočega se betona pri postopni gradnji masivnega armiranobetonskega vodnjaka smo analizirali natančnost in s tem ustreznost predstavljenega numeričnega postopka za analizo postopne gradnje masivnih betonskih konstrukcij. Dodatno smo s tem postopkom določili optimalni režim gradnje armiranobetonskega masivnega vodnjaka. Optimalne debeline betonskih slojev vodnjaka in časovni presledek med betoniranji le-teh smo določili z upoštevanjem zahtevanih temperaturnih omejitev skladno s prEN 13670:2007 oziroma drugimi dostopnimi eksperimentalnimi rezultati. 2'TEORETIČNE OSNOVE Časovno spreminjanje temperatur v strjujočem se betonu pri postopni gradnji masivne betonske konstrukcije določimo z rešitvijo znane Fourierove parcialne diferencialne enačbe prevajanja toplote v trdni snovi: PC> dt ~dx{ xdx)+dy V (1) sevajočih teles v okolici obravnavane konstrukcije, o Stefan - Boltzmanova konstanta, in e emisivnost sevajoče površine. TB v enačbi (3) je predpisana temperatura na površini betonskega elementa, T0 v enačbi (4) pa je začetna temperatura betona. Skladno z metodo končnih elementov lahko enačbo (1) zapišemo v diskretizirani obliki [Cook, 1989]: Kt + Ct, =f, (5) Ker v naravnem okolju adiabatnih pogojev običajno ni, moramo količino sproščene hid-ratacijske toplote betona ustrezno popraviti. S tem namenom vpeljemo zrelostno funkcijo fA s katero v analizi upoštevamo vpliv temperature in časa na sproščanje hidratacijske toplote. V svojem modelu smo uporabili Arrheniusovo zrelostno funkcijo, ki jo določa enačba: /A(r,eJ=exp R 1 Pcp) + T" Enačbo (1) rešimo z upoštevanjem naslednjih robnih in začetnih pogojev: , dr , dr dT k —n +k —n +k —n + Xdx x ydyy z dz 2 + h(T-TA)+ae{T*-TR4)=0, (2) kjer so K, C in f globalna prevodnostna matrika, globalna kapacitetna matrika in globalni vektor desnih strani, t je vektor vozliščnih temperatur, tt pa vektor časovnih odvodov vozliščnih temperatur. Sproščeno toploto Qa(t) v adiabatnih razmerah določimo na osnovi adiabatne krivulje Ta(t) z enačbo: (8) V enačbi (8) R predstavlja plinsko konstanto (R= 8,315 J/mol K), EA je aktivacijska energija, Ta pa začetna temperatura betona pri adiabatnem preizkusu. Količino sproščene hidratacijske toplote v betonski konstrukciji lahko tako izrazimo z enačbo: T{x,y,z,t)=TB{t), (3) T{x,y,z$) = T0(x,y,z), (4) kjer sta p in cp gostota in specifična toplota betona, kx, ky in kz pa so toplotne prevodnosti betona v smereh x, y in z. V enačbah (1-4) je T= T(x, y z, t) temperatura betona pri času t, Q pa je sproščena hidratacijska toplota betona, ki jo določimo na osnovi adi-abatnega preizkusa. V enačbi (2) so z nx, ny in nz označene komponente vektorja normale na robno ploskev, TA je temperatura okolice, h prestopni koeficient zraka, TR temperatura Qa(t)=Ta{t)pcp. (6) V literaturi zasledimo različne enačbe za določitev adiabatne krivulje ([Morabito, 2001 ], [Ammar, 1973], [Van Breugel, 1997]). V prikazanem numeričnem postopku za analizo postopne gradnje masivnih betonskih konstrukcij izberemo za adiabatno krivuljo krivuljo sigmoidne oblike [Ammar, 1973]: Ta(t) + (7) kjer so a, b, c in d materialni parametri. Določimo jih z metodo najmanjših kvadratov za vsako betonsko mešanico posebej. Q = fAQa- (9) Temperatura strjujočega se betona pri masivnih betonskih konstrukcijah je v veliki meri odvisna od postopne gradnje masivnih blokov. Zato smo razvili numerični model s pripadajočim računalniškim programom, s katerim s podajanjem začetnega in končnega časa aktivnosti posameznega končnega elementa in ustreznih robnih pogojev relativno preprosto modeliramo postopno gradnjo. Bistveni operaciji, ki jih mora računalniški program izvajati za modeliranje postopne gradnje, sta dodajanje in odvzemanje končnih elementov, vozlišč in robnih pogojev v globalni mreži končnih elementov. Ko se v nekem časovnem koraku zaradi postopne gradnje pojavi nov končni element, ga mora program dodati v globalno mrežo končnih elementov, kar shematsko prikazujemo na sliki 1. Tedaj program preveri, če novi element vključuje nova vozlišča. Če jih vključuje, se začetne temperature teh vozlišč dodajo v globalni sistem enačb kot nove pro-stostne stopnje. Novi končni element s podano začetno temperaturo vozlišč ne vpliva samo na nova vozlišča, ampak tudi na že obstoječa. V vozliščih, kjer je novi končni element v stiku z že obstoječimi, se temperatura izračuna kot obteženo povprečje, kjer je utež toplotna vdornost snovi b [Kladnik, 1983]. V pravokotni mreži dvodimenzionalnih končnih elementov se lahko stikajo največ štirje elementi. Temperaturo skupnega vozlišča T v tem primeru izračunamo z enačbo: J1 _ i=1_ "v bi=pfpi - (10) kjer je T temperatura obravnavanega vozlišča /-tega končnega elementa, b, X,, p in c pa so toplotna vdornost, toplotna prevodnost, gostota in specifična toplota /-tega elementa (/ = 1 ... nV). Z nV smo označili število elementov, ki imajo skupno obravnavano vozlišče. V večini primerov je nV enako 4, lahko pa je tudi 2 ali 3. Na sliki 1 so z zeleno barvo označena vozlišča, ki so aktivna ves čas analize, z rdečo pa tista, ki so postala aktivna po vključitvi novih elementov. Z modro barvo so označena vozlišča, ki so bila pred dodajanjem novih elementov na robu, ob vključitvi novih elementov pa moramo njihovo temperaturo določiti po enačbi (10). Ob vključitvi novih končnih elementov postanejo nekateri zunanji robovi meje med Slika 1* Shematski prikaz dodajanja novih končnih elementov v globalno mrežo končnih elementov elementi (modra črta), zato tam robnih pogojev ne moremo predpisati, nove robne pogoje pa dobimo na zunanjih robovih novih elementov (rdeča črta). Pomembna prednost prikazanega postopka za modeliranje postopne gradnje je tudi v tem, da poleg dodajanja novih končnih elementov v globalno mrežo omogoča tudi njihovo odstranjevanje. To bi na primer lahko uporabili za modeliranje toplotne izolacije ali pa opaža, ki ga seveda po nekem času odstranimo. 3 * PRIMERJAVA MED EKSPERIMENTALNIMI IN NUMERIČNIMI REZULTATI 3.1 Meritve Pri gradnji odseka hitre ceste Vipava-Razdrto preko pobočja Rebernice so bili kot temelji hitre ceste pri premoščanju plazovitega območja zgrajeni trije masivni armiranobetonski vodnjaki. Vodnjake s premerom 8 m in globino približno 40 m s približno pol metra debelim armiranobetonskim opažem (slika 2) so gradili po slojih. Za zagotavljanje ustrezne kvalitete vgrajenega betona je bilo potrebno med gradnjo skladno s tehnološkim elaboratom določiti najvišje temperature strjujočega se betonu. Tako se je med gradnjo vodnjakov z ustrezno eksperimentalno opremo in računalniško podporo merilo časovni potek temperatur strjujočega se betona v spodnjih štirih slojih enega izmed vodnjakov. Debeline betonskih slojev in razporeditev merskih točk prikazujemo na sliki 2. Poleg merjenja temperatur betona v 12 točkah je bila merjena tudi temperatura zraka, vedno pol metra nad trenutno najvišjim slojem. Betoniranje prvih štirih slojev je potekalo 23. 6., 3. 7., 15. 7. in 22. 7. 2008, temperature pa so bile z nekaj prekinitvami merjene med 23. 6. in 25. 7. 2008, torej v poletnem času. Zaradi omejenega števila merilnih kanalov uporabljenega merilnega inštrumenta se je merjenje v nekaterih merskih točkah spodnjih slojev opuščalo, kar je omogočilo meritve temperatur v zgornjih slojih. O .Z P6 m L si WM