何 強 王雪峰

(濟南市水利工程服務(wù)中心,山東 濟南 271199)

混凝土重力壩由于相對安全可靠、耐久性好、對地形和地質(zhì)條件適應(yīng)性好等特點得到了廣泛的應(yīng)用。但與此同時,重力壩壩體體積大,施工期溫度應(yīng)力和收縮應(yīng)力較大,易產(chǎn)生溫度裂紋[1]。在國內(nèi)外的混凝土大壩中普遍存在著裂縫,由于裂縫而失事或影響運行的實例很多[2-4]。

裂紋的產(chǎn)生和拓展不僅破壞混凝土壩的完整性,而且影響結(jié)構(gòu)穩(wěn)定性和安全性,威脅大壩的正常運行?；炷翂瘟芽p的形成原因十分復(fù)雜,現(xiàn)階段主要采用物理模型試驗[5]、工程經(jīng)驗[6]和數(shù)值仿真等[7-8]方法研究裂紋的成因。物理模型試驗在一定程度上能反映重力壩的破壞規(guī)律,但由于在實際工程中壩體的裂紋開度僅有數(shù)厘米,對比到模型中僅有毫米級別,較難觀測到裂紋的發(fā)展過程,難以滿足相似律要求。并且由于模型試驗實驗周期長,耗資較大,適用性不強,因此采用數(shù)值計算方法是研究裂紋成因的較優(yōu)方案。

在現(xiàn)階段的研究中,學(xué)者主要集中于分析在施工期溫度對于混凝土壩裂紋產(chǎn)生和拓展的影響[9-10]。針對運行期溫度對于混凝土壩壩體初始裂紋影響的研究不多。但是最近的研究結(jié)果表明,運行期氣溫日變化和寒潮等對運行期壩體溫度場和應(yīng)力場也具有較大的影響。

現(xiàn)階段我國很多已建的混凝土重力壩,由于所處的氣候條件不同,出現(xiàn)了各種不同程度的裂紋。本文針對廣東省楓樹壩水電站溢流壩段反弧段,在運行期出現(xiàn)貫穿性裂紋的現(xiàn)象進(jìn)行分析。研究主要采用數(shù)值模擬的方式,通過建立溢流壩段有限元模型,依據(jù)實測數(shù)據(jù)和工程經(jīng)驗確定模擬參數(shù),模擬大壩在運行期日平均氣溫下的溫度場和應(yīng)力場變化。通過模擬結(jié)果分析裂紋成因,從而為工程實踐提供參考。

1 算法理論

本文采用有限元對早期溢流壩反弧段的裂縫成因進(jìn)行分析。監(jiān)測數(shù)據(jù)表明大壩并未出現(xiàn)明顯的位移和沉降,推測為是由溫度變化導(dǎo)致的反弧段變電洞頂部裂縫。模擬時首先計算大壩內(nèi)部的溫度場分布,再依據(jù)溫度場計算結(jié)果求解應(yīng)力場。

1.1 溫度場計算

對于計算域R內(nèi)的任意一點,溫度場T(x,y,z,τ)必須滿足如下的熱傳導(dǎo)方程:

?2T=0

(1)

式中:T為混凝土溫度,℃。

在邊界上滿足第一類C或者第三類邊界條件C′:

(2)

式中:β為材料表面散熱系數(shù),kJ/(m2·h·℃);λ為導(dǎo)熱系數(shù),kg/(m·d·℃);Ta為環(huán)境溫度,℃;n為邊界的法向。

通過變分原理,得到溫度場的變分方程,并采用八節(jié)點六面體單元將計算域離散。利用經(jīng)典的伽遼金法,并采用高斯積分計算離散域的積分,得到溫度場計算公式:

(3)

式中:[H]為熱傳導(dǎo)矩陣;[R]為熱傳導(dǎo)補充矩陣;{Tn}和{Tn+1}為節(jié)點n和n+1時刻的溫度;{Fn+1}為節(jié)點的溫度荷載矩陣;Δτn為計算的時間步長。依據(jù)方程,{Tn+1}可以通過{Tn}得到。[H]和[R]是計算中的不變量,{Fn+1}是計算中的已知量。

1.2 應(yīng)力場計算

(4)

其中彈性應(yīng)變增量的表達(dá)式為

(5)

式中:[D]為彈性矩陣。

溫度應(yīng)變增量由非穩(wěn)定溫度場的計算結(jié)果推求得到,其表達(dá)式如下:

(6)

式中:α為材料的線性膨脹系數(shù),10-6/℃;ΔTn為Δτn時間段內(nèi)的溫度變化量,℃。

應(yīng)力場的計算區(qū)域和網(wǎng)格剖分與溫度場相同,依據(jù)虛功原理建立整體平衡方程:

(7)

得到節(jié)點的位移增量后,節(jié)點的應(yīng)力增量由式(8)計算:

(8)

累加每一步的增量結(jié)果,得到任意時刻的位移(δn)和應(yīng)力(σn)。

2 計算模型與參數(shù)

2.1 計算模型

本文研究對象為廣東省楓樹壩水電站,其在運行過程中溢流壩面反弧段區(qū)域出現(xiàn)裂紋,主要研究其裂縫成因。該溢流壩段高79.4m,上下游長度為86.8m,壩段寬度為17.0m。圖1為按照溢流壩段剖分出的有限元模型,總單元數(shù)量為67809個,總節(jié)點數(shù)為76984個。為保證計算的準(zhǔn)確性,考慮地基對于計算結(jié)果的影響,地基在大壩上下游各擴展100m,在大壩左右兩側(cè)各拓展54m,深度方向拓展150m,總的有限元計算模型見圖2。由于溢流壩段反弧段,即變電洞頂部出現(xiàn)裂紋,為突出研究重點,將此處的區(qū)域的網(wǎng)格細(xì)分,見圖3。

圖1 有限元計算模型(僅溢流壩段)

圖2 有限元計算模型(整體)

圖3 變電洞附近區(qū)域有限元模型

2.2 計算參數(shù)

根據(jù)大壩的設(shè)計資料并結(jié)合相關(guān)工程經(jīng)驗,計算時大壩溢流壩段混凝土以及地基的熱力學(xué)參數(shù),取值分別見表1和表2?？紤]到大壩已運行40多年,按偏于工程安全考慮,計算時不考慮混凝土徐變對計算結(jié)果的影響。

表1 壩體混凝土基本參數(shù)

表2 基巖基本參數(shù)

壩址處的月平均氣溫見表3,計算時多年平均日氣溫變化擬合為式(9),實測多年平均氣溫及其擬合值見圖4。

表3 月平均氣溫

(9)

式中:t為月份。

變電洞內(nèi)部由于在大壩內(nèi)部,并不與外界直接接觸,因此存在溫室效應(yīng),內(nèi)部溫度變幅較小。依據(jù)工程經(jīng)驗給出變電洞內(nèi)部溫度擬合公式:

(10)

式中:t為月份。

根據(jù)水庫水位觀測記錄,壩前最高水位為165.53m,低于水庫正常蓄水位166.0m;壩前最低水位為125.83m,低于水庫死水位128.0m。各月的平均水位見表4。

表4 月平均水位

計算溫度場時需要考慮庫水溫度,由于水庫無水溫監(jiān)測設(shè)備,采用式(11)～式(15)計算任意深度的任意時刻水溫[11]。

(11)

Tm(y)=c+(Ts-c)e-0.04y

(12)

c=(Tb-Ts)/(1-g),g=e-0.04H=(11-23.5)/

(1-e-0.04H)

(13)

A(y)=A0e-0.018y=8.57e-0.018y

(14)

ε=2.15-1.30e-0.085y

(15)

式中:T(y,τ)為水深y處在時間為τ時的溫度,℃;y為水深,m;τ為月份;τ0為氣溫最高的時間,計算中取為7.1;Tm(y)為任意深度的年平均水溫;Tb為庫底年平均水溫,廣東屬于華南地區(qū),取值為11℃;Ts為表面年平均水溫,其數(shù)值為年平均氣溫(20.5℃)加上溫度增量(2～4℃),計算時取值為23.5℃;H為水庫深度;A(y)為水溫年變幅,其中表面水溫年變幅A0與氣溫年變幅相同,取為8.57;ε為水溫相位差。

依據(jù)風(fēng)速觀測記錄,并參考文獻(xiàn)[11],得到壩體與外界直接接觸的表面,放熱系數(shù)為1095.6kJ/(m2·d·℃)。變電洞內(nèi)部為相對封閉的空間,假定表面風(fēng)速為0,放熱系數(shù)為400kJ/(m2·d·℃)。由于大壩內(nèi)部空間(包括廊道以及寬縫等)完全封閉,按照絕熱邊界進(jìn)行處理。

計算準(zhǔn)穩(wěn)定溫度場時,邊界條件為:上游壩基和上游壩面水面以下部分,按照第一類邊界條件處理,上游水位取多年平均水位149.24m,按照式(11)～式(15)計算各個節(jié)點的水溫;下游壩基和下游壩面水面以下部分,按照第一類邊界條件處理,受電站尾水?dāng)_動,假定下游水溫均勻分布,水溫取電站引水管進(jìn)口處對應(yīng)的水溫,水位按照正常尾水位計算;其余壩基面,按照絕熱邊界處理;壩體與外界大氣接觸的表面,按照第三類邊界條件處理,放熱系數(shù)為1095.6kJ/(m2·d·℃),溫度按照式(9)計算;變電洞內(nèi)部表面,按照第三類邊界條件處理,放熱系數(shù)為400.0kJ/(m2·d·℃),溫度按照式(10)計算;大壩內(nèi)部空間包括廊道以及寬縫等,由于完全封閉,按照絕熱邊界進(jìn)行處理。

在進(jìn)行應(yīng)力場計算時,假定地基四周以及底面為連桿垂直約束,考慮到重力壩屬于平面應(yīng)變問題,因此壩段側(cè)面仍采用連桿垂直約束。荷載施加如下:上游壩基和上游壩面水面以下部分,受到垂直于面的水壓力作用,水位與溫度邊界條件取值保持一致。下游壩基和下游壩面水面以下部分,受到垂直于面的水壓力作用,水位與溫度邊界條件取值保持一致。結(jié)構(gòu)自重按照體荷載進(jìn)行施加,豎直向下?？傆嬎銜r長為45年,即從大壩開始運行至如今。

3 裂縫成因分析

3.1 常規(guī)運行期

計算溢流壩段在常規(guī)運行期內(nèi),在多年平均日氣溫影響下的準(zhǔn)穩(wěn)定溫度場和相應(yīng)的應(yīng)力場。圖5～圖8分別為4月、7月、10月和1月大壩中心剖面溫度場和應(yīng)力場云圖。需要說明的是本文應(yīng)力云圖中的應(yīng)力均為第一主應(yīng)力。

圖5 常規(guī)運行期4月

圖7 常規(guī)運行期10月

圖8 常規(guī)運行期1月

從溫度云圖可以看出,溢流壩段上游尤其在水位以下區(qū)域的混凝土主要受到庫水溫度的影響,溫度較低且溫度變化小,不易受到外界氣溫變化的影響。溢流壩面以及變電洞附近區(qū)域的混凝土,易受氣溫變化的影響,夏季表面溫度高,冬季則溫度低,溫度年變幅較大。從應(yīng)力云圖可以看出,在4月時溢流壩段整體應(yīng)力較小,均低于0.6MPa,結(jié)構(gòu)安全。然而處于夏季7月時,變電洞頂部偏左側(cè)的內(nèi)部區(qū)域,拉應(yīng)力較大,達(dá)到了2.4MPa,超過混凝土的抗拉強度。大壩安全鑒定的結(jié)果表明該部分區(qū)域的混凝土出現(xiàn)了深層裂縫,與計算結(jié)果相近。處于秋季10月時,觀察到變電洞頂部、右側(cè)薄墻和右下方墻角以及溢流壩反弧段附近區(qū)域,拉應(yīng)力較大,達(dá)到了1.4MPa,但未超過抗拉強度。處于冬季1月時,溢流壩反弧段拉應(yīng)力較大,達(dá)到了1.6MPa,結(jié)構(gòu)比較危險,易產(chǎn)生裂縫。

3.2 成因分析

為進(jìn)一步分析兩條主要裂縫產(chǎn)生的原因,選取兩個特征點(見圖9),提取其溫度和應(yīng)力歷時曲線,見圖10。從圖10中可以看出位于變電洞頂部的1號特征點,溫度與應(yīng)力變化趨勢相同,即溫度升高,拉應(yīng)力增大。由于位于變電洞內(nèi),因此其溫度變幅較小,但是應(yīng)力變幅較大,拉應(yīng)力峰值甚至達(dá)到了2.64MPa,超過允許抗拉強度。圖11為實際拍攝得到的溢流壩段變電洞頂部出現(xiàn)的裂紋,可以發(fā)現(xiàn)模擬與觀測得到的裂紋所在位置相同。

圖9 特征點位置示意圖

位于反弧段附近的2號特征點,溫度與氣溫變化趨勢相反,即溫度降低,拉應(yīng)力減少。由于與空氣直接接觸,因此溫度變幅較大。同時在冬季溫度降低時,由于溢流面與空氣直接接觸導(dǎo)致混凝土因溫度降低而收縮,也產(chǎn)生了較大的拉應(yīng)力,拉應(yīng)力峰值達(dá)到了1.62MPa,雖未達(dá)到允許抗拉強度1.9MPa,但仍較為危險。

下面主要討論變電洞頂部拉應(yīng)力產(chǎn)生的原因。圖12為變電洞附近夏季7月的溫度和應(yīng)力場云圖,可見夏季溢流面附近溫度高,變電洞內(nèi)部由于處于室內(nèi)溫度相對較低,兩者溫差為4～5℃。變電洞頂部偏左側(cè)區(qū)域以及右下方角點處的拉應(yīng)力較大,結(jié)構(gòu)極易產(chǎn)生裂縫;而其余區(qū)域的拉應(yīng)力相對較小,不超過1.2MPa,相對安全。

圖12 應(yīng)力最大時刻變電洞附近模擬結(jié)果

提取應(yīng)力最大時刻變電洞附近變形和位移(見圖13),可以發(fā)現(xiàn)在夏季時,由于溫度的作用,變電洞附近尤其是頂部與右側(cè)薄壁,變形較大。變電洞頂部上側(cè)區(qū)域混凝土受熱膨脹,擠壓下游薄壁墻體向右變形,使得變電洞右下側(cè)轉(zhuǎn)折點處應(yīng)力較大。夏季下游薄壁墻體膨脹,頂托變電洞頂部區(qū)域向上變形,變電洞頂部偏左側(cè)下部區(qū)域?qū)⑹芾?同時,變電洞頂部上側(cè)混凝土受熱膨脹,將拉動頂部下側(cè)區(qū)域混凝土向右位移,因此變電洞頂部偏左側(cè)下部區(qū)域也將產(chǎn)生較大的拉應(yīng)力。

圖13 應(yīng)力最大時刻變電洞附近位移和變形(放大800倍)

3.3 監(jiān)測結(jié)果

為監(jiān)測主要的裂縫情況,于2017年1月在變電洞頂部水平裂縫處安裝4支測縫計(J01～J04)。測縫計監(jiān)測結(jié)果與氣溫水溫關(guān)系,見圖14。由于測縫計安裝于變電洞頂部水平裂縫處,不與大氣直接接觸,相當(dāng)于室內(nèi)環(huán)境,因此隨溫度變化并不明顯。其中J01和J03監(jiān)測點7月開合度最大,1月開合度較小,這與工況1的計算結(jié)果相吻合,即7月變電洞拱頂偏左測區(qū)域的混凝土將會受拉。J02監(jiān)測點裂縫開合度的變化并不明顯。J04監(jiān)測點變化趨勢與J03相近,但幅度較小。

圖14 變電洞頂部裂縫開合度監(jiān)測結(jié)果與氣溫關(guān)系

4 結(jié) 語

本文針對廣東省楓樹壩水電站混凝土重力壩溢流壩段變電洞頂部和溢流面反弧段出現(xiàn)裂縫的問題,采用數(shù)值模擬的方式探究裂縫成因。通過建立有限元,依據(jù)大壩設(shè)計資料和工程經(jīng)驗選取合理的計算參數(shù),模擬了大壩在常規(guī)運行期的溫度場和應(yīng)力場分布。

研究結(jié)果表明夏季高溫和強結(jié)構(gòu)約束是造成變電洞頂部偏左側(cè)區(qū)域產(chǎn)生深層裂縫的主要原因。由于夏季溫度升高,導(dǎo)致壩體內(nèi)部與表面存在較大的溫度差,且由此導(dǎo)致的溫度變形受到了較大的結(jié)構(gòu)約束,從而導(dǎo)致溢流壩反弧段即變電洞左側(cè)頂部區(qū)域拉應(yīng)力較大,達(dá)到了2.64MPa,超過混凝土允許抗拉強度(1.5MPa),產(chǎn)生了貫穿性的裂紋。

模擬結(jié)果與實際觀測到的裂縫出現(xiàn)位置相對應(yīng)。計算所得到應(yīng)力場隨溫度的變化規(guī)律與測縫計測量得到的裂縫開合度監(jiān)測結(jié)果相符合。驗證了計算結(jié)果的正確性和推斷原因的合理性。

當(dāng)然由于條件限制本文在模擬計算時,并未考慮混凝土澆筑后由于水化反應(yīng)導(dǎo)致的溫度應(yīng)力,這將在下一階段的研究中進(jìn)行完善。并考慮合理的補強方案,防止結(jié)構(gòu)破壞。