多場耦合下大體積混凝土初次蓄水的溫度應(yīng)力問題研究摘要：在大體積混凝土壩初次蓄水時(shí)，溫度較低的庫水必然會對壩體溫度場產(chǎn)生較明顯的影響，從而影響壩體的變形，甚至產(chǎn)生溫度裂縫.為分析大體積混凝土初次蓄水的溫度應(yīng)力，本文將混凝土類多孔介質(zhì)視為連續(xù)介質(zhì)，綜合運(yùn)用水力學(xué)、熱學(xué)和固體力學(xué)等基本理論，根據(jù)動量守恒、質(zhì)量守恒和能量守恒方程建立了以位移、孔隙水壓力、孔隙氣壓力、溫度和孔隙率為未知量的多場耦合數(shù)學(xué)模型，在此基礎(chǔ)上編制了有限元計(jì)算程序，并對大體積碾壓混凝土塊的滲流場、溫度場和應(yīng)力場進(jìn)行了耦合分析，結(jié)果表明，考慮耦合后塊體溫降幅度及溫度大主應(yīng)力均較不考慮耦合條件下大.關(guān)鍵詞：多場耦合；數(shù)學(xué)模型；大體積混凝土；溫度應(yīng)力中圖分類號：TU43文獻(xiàn)標(biāo)識碼：A文章編號：1674-2974（2016）05-0030-09Abstract：Duringtheinitialimpoundmentofmassconcretedam，thereservoirwaterwithlowtemperaturewouldbecriticaltothetemperaturefieldofdambody，whichaffectsthedeformationofdambodyandevenresultsintemperaturecracks.Therefore，inordertoinvestigatethethermalstressdistributionofmassconcreteduringtheinitialimpoundment，theconcretewasassumedasacontinuousporousmediainthispaper.Inconsiderationofthebasictheoriesofsolidmechanics，hydraulics，andthermodynamics，themulti-fieldcouplingequationsofunsaturatedporousmediathatincludemomentumconservation，massconservation，andenergyconservationwereprovidedasthefunctionofdisplacements，poreliquidpressure，poregaspressure，temperature，andporosity.Thefiniteelementanalysisprogramwasthendeveloped.Amassroller-compactedconcreteblockwasconsideredforcouplinganalysesontheseepagefield，temperaturefield，andstressfield.Theanalysisresultsshowthatthetemperaturereductionandprinciplethermalstressoftheconcreteblockconsideringthemulti-fieldcouplingprocessaregreaterthanthosewithoutthecouplingeffect.Keywords：multi-fieldcoupling；mathematicalmodel；massconcrete；thermalstress在水利工程中，尤其是對大體積混凝土而言，較大的溫度拉應(yīng)力會引起混凝土的溫度裂縫，因此溫度應(yīng)力是大體積混凝土需要額外關(guān)注的.在水庫初次蓄水時(shí)，溫度較低的庫水必然會對壩體溫度場產(chǎn)生較明顯的影響，從而影響混凝土壩的變形，甚至產(chǎn)生溫度裂縫.如國際上，在Revslstoke[[1]，Dworshak[[2]和Russel[[3[]等重力壩的上游表面曾經(jīng)產(chǎn)生過嚴(yán)重的劈頭裂縫，深入壩內(nèi)幾十米，有的甚至將整個壩段一分為二，產(chǎn)生嚴(yán)重漏水.目前一般都認(rèn)為是在施工過程中壩體的上游側(cè)如果出現(xiàn)了表面裂縫，水庫蓄水之后，經(jīng)過一段時(shí)間，表面裂縫突然大范圍地?cái)U(kuò)展，成為劈頭裂縫，尤其是通倉澆筑的混凝土重力壩更容易出現(xiàn)這種劈頭縫，分析其原因主要是由于這類壩沒有布置縱縫，不進(jìn)行二期冷卻，在水庫初次蓄水時(shí)，壩體溫度還很高，與外界低溫的庫水之間形成較大的內(nèi)外溫差，容易使表面擴(kuò)展成為劈頭縫[[4].由此可見，初次蓄水對混凝土重力壩的溫度場和應(yīng)力場的影響是需要加以研究的，并為采取有效、合理地滲控措施提供科學(xué)依據(jù).事實(shí)上，如果將混凝土壩體也視為多孔介質(zhì)，則蓄水時(shí)滲流的產(chǎn)生也是必然的，若考慮庫水入滲，壩體非穩(wěn)定溫度場的變化將更加復(fù)雜.考慮滲流因素的存在，對大體積混凝土壩體溫控防裂及應(yīng)力狀態(tài)研究都有直接意義.國內(nèi)已經(jīng)有學(xué)者用多場耦合的方法來求解大體積混凝土的滲流場溫度場耦合作用[[5-7[]、滲流場溫度應(yīng)力耦合作用[[8[]等.本文在總結(jié)國內(nèi)外學(xué)者對多孔介質(zhì)多場耦合的機(jī)理、數(shù)學(xué)模型建立及求解方法的研究成果基礎(chǔ)上，將混凝土視為多孔介質(zhì)，綜合運(yùn)用固體力學(xué)、水力學(xué)、熱學(xué)等基本理論，結(jié)合多孔介質(zhì)的熱本構(gòu)關(guān)系以及孔隙流體的熱運(yùn)動規(guī)律，建立了混凝土多相、多場全耦合數(shù)學(xué)模型方程組，并編制了有限元求解程序，對某一大體積混凝土模型初次蓄水后一段時(shí)間內(nèi)的滲流場、溫度場、應(yīng)力場進(jìn)行耦合分析，以初步闡釋大體積混凝土初次蓄水時(shí)滲流場對壩體溫度場、應(yīng)力場的影響過程.1混凝土類多孔介質(zhì)概念模型及熱本構(gòu)關(guān)系模型1.1混凝土概念模型國內(nèi)外有些學(xué)者提出了關(guān)于多孔介質(zhì)的多場耦合模型[[9-13[].事實(shí)上，非飽和狀態(tài)下的混凝土材料也可視為連續(xù)性多孔介質(zhì)，其特征單元體主要由固相、液相和氣相三相構(gòu)成，其中固相介質(zhì)主要為經(jīng)硬化后形成具有堆聚結(jié)構(gòu)的復(fù)合物.特征單元體的體積V可以用式（1）表示：1.2混凝土類多孔介質(zhì)熱本構(gòu)關(guān)系模型1.2.1固相介質(zhì)的熱本構(gòu)關(guān)系當(dāng)溫度發(fā)生變化時(shí)，混凝土類多孔介質(zhì)將由于受熱發(fā)生膨脹、受冷而發(fā)生收縮，從而產(chǎn)生應(yīng)變.如果多孔介質(zhì)熱變形受到邊界的約束，就會產(chǎn)生應(yīng)力.混凝土的熱本構(gòu)關(guān)系表達(dá)方式比較多[[14-18[]，為了利用這些成熟的本構(gòu)關(guān)系，同時(shí)考慮溫度變化對多孔介質(zhì)的影響，可以給出以增量形式表述的熱本構(gòu)關(guān)系如式（3）所示：1.2.4溫度場的熱傳導(dǎo)定律在研究溫度場的問題時(shí)，廣義的Fourier定律[[21]作為基本定律，它是指在導(dǎo)熱過程中，單位時(shí)間內(nèi)通過給定截面所傳遞的熱量，正比例垂直于該截面方向上的溫度變化率，而熱量傳導(dǎo)的方向與溫度升高的方向相反，用數(shù)學(xué)表達(dá)式表示為：2混凝土類多孔介質(zhì)多場耦合的控制方程混凝土類多孔介質(zhì)非飽和狀態(tài)下多場耦合的控制方程將主要基于上述非飽和多孔介質(zhì)的概念模型及其相應(yīng)的熱本構(gòu)關(guān)系[[22]：2.1變形場控制方程對于一個非飽和多孔介質(zhì)的單元體來說，其準(zhǔn)靜態(tài)下的應(yīng)力平衡方程如式（7）：3.2有限元求解程序編制本文基于上述有限元離散成果利用FORTRAN語言研發(fā)了多孔介質(zhì)多場耦合求解程序THM-3D，該程序采用模塊式開發(fā)，每塊均具有獨(dú)立、明確的功能含義，以充分滿足實(shí)際工程中問題復(fù)雜的要求，主要算法流程圖如圖1所示，限于篇幅，該程序的驗(yàn)證將在其他文章給出.4工程算例假設(shè)有一大體積碾壓混凝土塊，如圖2所示，該混凝土塊高94.0m，上游側(cè)假定有蓄水至77.0m，下游側(cè)無水.該碾壓混凝土塊為澆筑式施工，在4月1日起開始澆筑，混凝土入倉溫度在氣溫的基礎(chǔ)上加3℃，分層澆筑至頂部后30d拆模.壩體在430d內(nèi)澆筑完成，之后將在10d內(nèi)分3段快速蓄水至正常蓄水位：在第3d蓄水至25.0m；第6d蓄水至50.0m；在第10d蓄至77.0m.計(jì)算周期為蓄水至2a，該混凝土塊體由3種級配的混凝土澆筑而成，從塊體上游至下游依次為變態(tài)混凝土、二級配混凝土和三級配混凝土.考慮到大體積混凝土主要關(guān)注的是溫度應(yīng)力的變化，因此應(yīng)力場將主要考慮溫度應(yīng)力.4.2計(jì)算結(jié)果分析4.2.1混凝土澆筑完成混凝土塊澆筑完成時(shí)的溫度和應(yīng)力計(jì)算結(jié)果如下圖3（a），（b）所示，其中溫度單位為℃，應(yīng)力單位為MPa.可以看出，在澆筑完成時(shí)，混凝土塊體溫度場的高溫區(qū)位于中下部和頂部位置，其中中下部的最高溫度達(dá)到44℃，主要是由于本算例中沒有考慮溫控設(shè)施，該處的澆筑溫度較高，且混凝土仍處于升溫階段，水泥水化熱產(chǎn)生的熱量無法散出；同時(shí)，塊體中均為拉應(yīng)力，但整體水平不高.4.2.2初次蓄水一段時(shí)間混凝土塊澆筑初次蓄水30d，2a的孔隙水壓力、溫度和應(yīng)力計(jì)算結(jié)果如圖4（a）～（c）和5（a）～（c）所示.從滲流場的計(jì)算結(jié)果來看，隨著蓄水時(shí)間的延長和庫水溫、氣溫的變化，庫水在塊體內(nèi)逐漸向下游側(cè)滲流，蓄水2a后庫水仍然沒有在塊體的下表面逸出，說明蓄水后混凝土塊體內(nèi)的滲流場在較長的時(shí)間內(nèi)不會達(dá)到穩(wěn)定狀態(tài)，也不會有水從下游面逸出.從溫度場的計(jì)算結(jié)果來看，在水庫蓄水后，隨著蓄水時(shí)間的延長，塊體上游側(cè)的溫度等值線逐漸向內(nèi)部延伸，這體現(xiàn)出了庫水對塊體溫度的影響過程.同時(shí)，塊體溫度場在庫水作用下，內(nèi)部溫度整體水平下降，并且上游側(cè)降溫比下游側(cè)明顯較快；塊體上游庫水以上部位在庫水和外界氣溫共同影響下溫度下降較快，此處的最高溫度由澆筑完成時(shí)的38℃降至蓄水2a后的18℃.從耦合條件下塊體溫度應(yīng)力分布的變化過程可以看出，隨著蓄水時(shí)間的延長，塊體上游側(cè)表面處的溫度由于受較低的庫水溫控制，與內(nèi)部相鄰部位的混凝土之間存在著溫差，蓄水后不同時(shí)期上表面出現(xiàn)了不同程度的拉應(yīng)力，局部拉應(yīng)力較大，可能對該部位的溫控防裂不利.4.2.3特征點(diǎn)溫度變化過程特征點(diǎn)1～4在耦合與非耦合條件下的溫度計(jì)算結(jié)果對比如圖6～9所示.通過對比可以看出，特征點(diǎn)1在非耦合條件下得到的溫度較耦合條件高，最大差值為1.5℃；特征點(diǎn)2在非耦合條件下的得到的溫度較耦合條件高，最大差值為1.8℃；特征點(diǎn)3在非耦合條件下得到的溫度較耦合條件高，最大差值為3.3℃，；特征點(diǎn)4在非耦合條件與耦合條件下得到的溫度計(jì)算結(jié)果吻合得較好，最大差值約0.56℃.可見，非耦合的方法對于上游面附近庫水浸沒區(qū)域以及其他庫水未浸沒區(qū)域的溫度計(jì)算結(jié)果較好，對于被庫水浸沒的塊體內(nèi)部區(qū)域的溫度計(jì)算結(jié)果則明顯高于耦合計(jì)算得到的結(jié)果.5結(jié)論與建議混凝土類多孔介質(zhì)多場耦合作用的數(shù)學(xué)模型，研究的主要是某一物理場方程中因變量或源匯項(xiàng)受其他物理場作用其變化的數(shù)學(xué)描述，也包括本構(gòu)規(guī)律的影響在控制方程中的反映，因此該類數(shù)學(xué)模型較為復(fù)雜，必定包含多個控制方程：1）本文根據(jù)連續(xù)介質(zhì)方法給出了混凝土類多孔介質(zhì)的概念模型，在給出基本假設(shè)的基礎(chǔ)上，依據(jù)熱本構(gòu)關(guān)系模型、流體運(yùn)動的廣義Darcy定律以及溫度場熱傳導(dǎo)的廣義Fourier建立了以位移、孔隙水壓力、孔隙氣壓力、溫度、孔隙率為未知量的非飽和多孔介質(zhì)多相多場全耦合研究的數(shù)學(xué)模型，并給出了相應(yīng)的定解條件，獲得了有限元格式的多場耦合求解方程組，并編制了有限元求解程序.2）對一大體積碾壓混凝土塊的滲流場、溫度場和應(yīng)力場進(jìn)行了耦合分析.結(jié)果表明，考慮耦合后塊體入滲區(qū)域內(nèi)特征點(diǎn)的溫降幅度較不考慮耦合條件下大，最大溫降差達(dá)到3.3℃；而且考慮耦合效應(yīng)得到的特征點(diǎn)溫度大主應(yīng)力也較非耦合大，最大差值達(dá)到0.32MPa.3）對于蓄水初期的實(shí)際大體積混凝土工程而言，影響其真實(shí)工作狀態(tài)的因素是極其復(fù)雜的，其中，作用在其上游側(cè)的低溫庫水則是眾多影響因素中最直接、最重要的因素.因此為了更為準(zhǔn)確地計(jì)算分析蓄水初期壩體溫度場和溫度應(yīng)力分布，應(yīng)該考慮滲流場的影響，本文采用耦合分析的方法為以后大體積混凝土壩蓄水后的溫控防裂研究提供了一種新的思路.參考文獻(xiàn)[1]BRUNNERWJ，WUKH.Crackingofrevelstokeconcretegravitydammassconcrete[R].London，UR：15thInternationalCongressonLargeDams，Vol.Ⅱ，1985：201-206.[2]HONGTONDL.Measuresbeingtakenforpreventionofcracksinmassconcreteatdworshakandlibbydam[R].Montreal，Canada：10thInternationalCongressonLargeDams，1979：1062-1066.[3]ZHUBF.Allowabletemperaturedifference，coolingcapacity，inspectionandtreatmentofcracks，andadministrationoftemperaturecontrol[M].Beijing：TsinghuaUniversityPress，2014：439-467.[4]朱伯芳.重力壩的劈頭裂縫[J].水力發(fā)電學(xué)報(bào)，1997（59）：85-93.[5]柴軍瑞.混凝土壩滲流場與穩(wěn)定溫度場耦合分析的數(shù)學(xué)模型[J].水力發(fā)電學(xué)報(bào)，2009，19（1）：27-35.[6]陳建余，朱岳明，張建斌.考慮滲流場影響的混凝土壩溫度場分析[J].河海大學(xué)學(xué)報(bào)：自然科學(xué)版，2003，3（2）：119-123.[7]崔皓東，朱岳明.蓄水初期的壩體非穩(wěn)定滲流場與溫度場耦合的理論模型及數(shù)值模擬[J].水利學(xué)報(bào)，2009，15（2）：238-243.[8]李守義，陳堯隆，王長江.碾壓混凝土壩滲流對溫度應(yīng)力的影響[J].西安理工大學(xué)學(xué)報(bào)，1996，12（5）：41-46.[9]HICKSTW，PINERJ，WILLIS-RICHARDSJ，etal.Ahydro-thermo-mechanicalnumericalmodelforHDRgeothermalreservoir[J].InternationalJournalofRockMechanicsandSciences，1996，33（5）：499-511.[10]CHAPULIOTS，GOURDINC，PAYENT，etal.Hydro-thermal-mechanicalanalysisofthermalfatigueinamixingtee[J].NuclearEngineeringandDesign，2005，235（2）：575-596.[11]MILLARDA，DURINM，STIETELA，etal.Discreteandcontinuumapproachestosimulatethethermo-hydro-mechanicalcouplingsinalarge，fracturedrockmass[J].InternationalJournalofRockMechanicsandSciences，1995，32（5）：409-434.[12]NGUYENTS，SELVADURAIAPS，ARMANDG.ModellingtheFEBEXTHMexperimentusingastatesurfaceapproach[J].InternationalJournalofRockMechanics&MiningSciences，2005，42（5）：639-651.[13]RUTQVISTJ，BRGESSONL，CHIJIMATSUM，etal.Thermohydromechanicsofpartiallysaturatedgeologicalmedia：governingequationsandformulationoffourfiniteelementmodels[J].InternationalJournalofRockMechanics&MiningSciences，2001，38（6）：105-127.[14]劉亞晨.核廢料貯存圍巖介質(zhì)THM耦合過程的數(shù)值模擬[J].地質(zhì)災(zāi)害與環(huán)境保護(hù)，2006，17（2）：78-82.[15]POLTRONIERIF，PICCOLROAZA，BIGONID，etal.Asimpleandrobustelastoplasticconstitutivemodelforconcrete[J].EngineeringStructures，2014，60（5）：81-84.[16]ZHOUYW，WUYF.GeneralmodelforconstitutiverelationshipsofconcreteandItscompositestructures[J].CompositeStructures，2012，94（3）：580-592.[17]THOMASG，ALAINM，JEANMACF，etal.Amultiaxialconstitutivemodelforconcreteinthefiresituation：theoreticalformulation[J].InternationalJournalofSolidsandStructures，2013，50（6）：3659-3673.[18]ZHAIY，DENGZC.Studyoncompressivemechanicalcapabilitiesofconcreteafterhightemperatureexposureandthermo-damageconstitutivemodel[J].ConstructionandBuildingMaterials，2014，68（3）：777-782.[19]張勇，薛禹群，謝春紅.高溫差條件下達(dá)西定律的理論推導(dǎo)[J].水科學(xué)進(jìn)展，1999，10（4）：362-367.[20]CHENYF，ZHOUCB，JINGLR.ModelingcoupledTHMprocessesofgeologicalporousmediawithmultiphaseflow：theoryandvalidatio