李炎隆，涂 幸，王海生，王瑞科

(1 西安理工大學(xué) 水利水電學(xué)院，陜西 西安 710048；2 陜西省水利廳，陜西 西安 710004；3 中國(guó)水電顧問(wèn)集團(tuán)西北勘測(cè)設(shè)計(jì)研究院，陜西 西安 710065)

混凝土面板承擔(dān)著向下游堆石體傳遞水壓力及大壩防滲的重任，對(duì)大壩的安全與工程的運(yùn)行起著決定性的作用[1]。從空間結(jié)構(gòu)看，面板的長(zhǎng)度較長(zhǎng)，寬度較小，而厚度相對(duì)來(lái)說(shuō)很小，是一塊長(zhǎng)條形的薄板，受環(huán)境溫度的影響較大，尤其是處于高寒地區(qū)的混凝土面板堆石壩[2]，在氣溫驟變條件下，面板內(nèi)將產(chǎn)生拉應(yīng)變，再加上壩體的沉降及其他不確定性因素，面板很容易出現(xiàn)裂縫。從國(guó)內(nèi)外修建的面板堆石壩來(lái)看，面板開(kāi)裂的情況比較普遍，個(gè)別的壩還相當(dāng)嚴(yán)重[3-6]。因此，考慮溫度荷載的作用，系統(tǒng)研究混凝土面板的應(yīng)力變形特性具有重要的意義。

國(guó)內(nèi)外關(guān)于面板堆石壩面板裂縫成因的研究結(jié)果表明：面板裂縫尤其是貫穿性裂縫的產(chǎn)生，主要是壩體的不均勻沉降及面板的溫度應(yīng)力共同作用的結(jié)果[7]。目前，關(guān)于混凝土面板堆石壩面板應(yīng)力變形特性的研究，主要集中在其于水荷載和自重作用下的應(yīng)力變形計(jì)算方面[8-10]，但這些研究多認(rèn)為面板尺寸相對(duì)較小，故將溫度應(yīng)力忽略不計(jì)，這顯然與混凝土面板堆石壩的實(shí)際應(yīng)力變形情況不符。面板的結(jié)構(gòu)特點(diǎn)決定了其溫度應(yīng)力取決于面板本身的溫度變化引起的應(yīng)力及墊層對(duì)面板的接觸約束2個(gè)方面的因素[11-14]。關(guān)于面板與墊層之間接觸約束的模擬問(wèn)題，傳統(tǒng)的面板應(yīng)力計(jì)算方法將面板與墊層之間的接觸面按彈性約束處理，然后計(jì)算面板的溫度應(yīng)力。而將面板與墊層之間的關(guān)系按“彈性約束”處理時(shí)，夸大了墊層對(duì)面板的約束作用，因此計(jì)算結(jié)果不夠準(zhǔn)確。

本研究擬基于混凝土面板與墊層之間接觸面的數(shù)值計(jì)算模型，考慮溫度荷載作用，模擬混凝土面板堆石壩的施工過(guò)程，并以公伯峽面板堆石壩為例進(jìn)行有限元仿真計(jì)算，以期為全面系統(tǒng)地研究堆石壩面板的應(yīng)力變形特性提供參考。

1 面板與墊層之間接觸面單元計(jì)算模型的基本原理

1.1 溫度場(chǎng)計(jì)算模型

如果將面板與墊層之間的接觸視作點(diǎn)-面接觸問(wèn)題，則接觸摩擦單元的熱傳導(dǎo)矩陣可表示為：

(1)

式中：[Kc]e為接觸摩擦單元的熱傳導(dǎo)矩陣；Rc為接觸摩擦單元的接觸熱阻，可由試驗(yàn)確定；{N}e為接觸摩擦單元形函數(shù)列陣；T為矩陣轉(zhuǎn)置符號(hào)。

節(jié)點(diǎn)熱荷載向量可表示為：

{Fc}e=qe{N}e。

(2)

式中：{Fc}e為節(jié)點(diǎn)熱荷載向量，qe為面板單元表面到墊層接觸節(jié)點(diǎn)的熱流量。

1.2 溫度應(yīng)力計(jì)算模型

考慮變溫作用時(shí)，接觸摩擦單元的等效單元?jiǎng)偠?約束方程可表示為：

(3)

(4)

(5)

接觸摩擦單元的等效單元?jiǎng)偠?約束矩陣Kc及等效荷載向量fcT分別為[15]：

(6)

(7)

式(6)與式(7)所表示的接觸摩擦單元的等效單元?jiǎng)偠?約束矩陣Kc及等效荷載向量fcT，可按標(biāo)準(zhǔn)的有限元集成規(guī)則迭加到整體溫度應(yīng)力計(jì)算的總剛度矩陣和總荷載向量中。

實(shí)際計(jì)算時(shí)，對(duì)于每個(gè)時(shí)間步ti，首先假定單元處于某種接觸狀態(tài)(固定、滑動(dòng)及自由)，由此確定相應(yīng)的約束荷載矢量a*，然后按式(3)進(jìn)行增量節(jié)點(diǎn)接觸應(yīng)力及接觸位移的試算，檢驗(yàn)是否與原假定狀態(tài)一致，若一致則計(jì)算結(jié)束，否則采用試算解為新的假定狀態(tài)，進(jìn)行新一輪迭代直至收斂。進(jìn)行接觸狀態(tài)判定時(shí)，法線方向的容許應(yīng)力[σ]取最大拉應(yīng)力，切線方向的容許應(yīng)力[τ]按Mohr-Coulomb準(zhǔn)則確定[16]。

2 工程實(shí)例分析

2.1 工程概況

公伯峽混凝土面板堆石壩最大壩高139 m，壩頂全長(zhǎng)429 m，壩頂寬10 m，上游壩坡坡比1∶1.4，下游局部壩坡比1∶1.5～1∶1.4，綜合壩坡比為1∶1.81。鋼筋混凝土面板頂端厚0.3 m，底部最大計(jì)算厚度0.76 m。面板按壩體應(yīng)力變形計(jì)算結(jié)果設(shè)置豎向縫，受拉區(qū)豎縫間距6 m，受壓區(qū)豎縫間距12 m，沿高程方向不設(shè)縫。面板混凝土標(biāo)號(hào)為C25，面板內(nèi)配置一層雙向鋼筋[17]。

2.2 計(jì)算模型

根據(jù)大壩橫剖面圖與沿壩軸線縱剖面圖，計(jì)算整體模型取面板F9所在的整個(gè)壩段，沿壩軸線方向?qū)挾葹?2 m。計(jì)算域壩基范圍為：在壩踵向上游和壩趾向下游各取約1倍壩高120 m，地基深度也取120 m。整體坐標(biāo)系的坐標(biāo)原點(diǎn)在壩踵面板F9中心位置。x軸正向?yàn)檠貕屋S線指向右岸，y軸正向指向下游，z軸正向?yàn)殂U直向上。

溫度場(chǎng)計(jì)算的邊界條件：壩基4個(gè)側(cè)面及底面、壩體x軸方向的2個(gè)側(cè)面按絕熱邊界條件施加；水位以下的上、下游面為固-水邊界，水位以下為固-氣邊界，固-水邊界按第一類邊界條件施加，固-氣邊界按第三類邊界條件施加。應(yīng)力場(chǎng)計(jì)算的邊界條件為：壩基底面及4個(gè)側(cè)面為全約束，除面板以外的壩體的x軸方向2個(gè)側(cè)面按x向簡(jiǎn)支處理，其余為自由邊界。壩體有限元計(jì)算模型如圖1所示。采用空間8節(jié)點(diǎn)六面體等參單元對(duì)計(jì)算模型進(jìn)行剖分，共剖分了25 040個(gè)單元，29 095個(gè)節(jié)點(diǎn)。

圖 1 公伯峽混凝土面板堆石壩壩體有限元計(jì)算模型

2.3 計(jì)算參數(shù)

壩址區(qū)氣溫、水溫、地溫、日溫差及寒潮等氣象資料均采用壩址下游的循化縣氣象站統(tǒng)計(jì)資料。其中，公伯峽元月中旬氣溫最低，7月中下旬氣溫最高，月平均氣溫在全年內(nèi)呈周期性變化。

以面板F9混凝土澆筑開(kāi)始時(shí)刻(2004-04-05)作為時(shí)間起點(diǎn)的氣溫?cái)M合函數(shù)表達(dá)式為：

(8)

式中：t為時(shí)間，h。

根據(jù)水庫(kù)蓄水計(jì)劃，仍以面板F9混凝土澆筑開(kāi)始時(shí)刻作為時(shí)間起點(diǎn)，擬合的水溫函數(shù)為：

(9)

式中：h為水深，m。

由于缺乏壩體材料及基巖熱力學(xué)參數(shù)的試驗(yàn)資料，因此結(jié)合公伯峽工程的具體情況，并參照類似工程的試驗(yàn)成果，熱力學(xué)參數(shù)選用結(jié)果見(jiàn)表1。

表 1 公伯峽混凝土面板堆石壩壩體材料及基巖的熱力學(xué)參數(shù)

混凝土材料及基巖用線彈性本構(gòu)模型模擬。其中，基巖的彈性模量為10 GPa，泊松比為0.25；混凝土絕熱溫升表達(dá)式為θ=42.72×(τ/(2.04+τ))，混凝土彈性模量E(τ)與齡期τ的關(guān)系式為：

E(τ)=25×(τ/(6.64+τ))。

(10)

式中:E(τ)為彈性模量(GPa)，泊松比為0.17。

墊層、過(guò)渡層以及主、次堆石區(qū)采用鄧肯-張(E-B)材料本構(gòu)模型，其力學(xué)參數(shù)如表2所示。

表 2 公伯峽混凝土面板堆石壩壩料的鄧肯-張模型(E-B)參數(shù)

2.4 計(jì)算工況

設(shè)計(jì)2個(gè)計(jì)算工況。工況1(不考慮溫度荷載)：上游為正常蓄水位2 005.0 m，下游無(wú)水；荷載為自重、水荷載，不考慮外界溫度變化及混凝土的水化熱對(duì)面板應(yīng)力的影響。工況2(考慮溫度荷載)：上游為正常蓄水位2 005.0 m，下游無(wú)水；荷載為自重、水荷載及由外界溫度變化和混凝土水化熱引起的溫度荷載。

3 計(jì)算結(jié)果及分析

3.1 不考慮溫度荷載作用時(shí)面板的應(yīng)力變形特性

對(duì)于混凝土面板堆石壩而言，混凝土面板的順坡向的應(yīng)力及撓度是反映其應(yīng)力變形特性的主要指標(biāo)。本研究計(jì)算的應(yīng)力以拉為正，以壓為負(fù)；面板撓度以指向壩內(nèi)為負(fù)。不考慮溫度荷載時(shí)面板撓度及順坡向應(yīng)力的分布分別如圖2、3所示。從圖2和圖3可以看出，面板撓度絕對(duì)值最大為0.24 m，基本位于面板的中心位置。正常蓄水位下，整塊面板在水壓力作用下成為一個(gè)凹曲面，在面板頂部及底部撓度的變化梯度較大，此時(shí)，面板的撓度主要受水壓力與堆石體變形的影響。面板的順坡向應(yīng)力基本為壓應(yīng)力，只有面板頂部與底部出現(xiàn)了極小的拉應(yīng)力，面板表面點(diǎn)順坡向壓應(yīng)力絕對(duì)值最大為9.05 MPa，底面點(diǎn)順坡向壓應(yīng)力絕對(duì)值最大為8.96 MPa，均出現(xiàn)在約1/2壩高處。面板表面點(diǎn)的順坡向壓應(yīng)力絕對(duì)值基本比底面點(diǎn)大，而在分布規(guī)律上，二者是一致的。

圖 2 不考慮溫度荷載時(shí)公伯峽混凝土面板堆石壩面板的撓度

圖 3 不考慮溫度荷載時(shí)公伯峽混凝土面板堆石壩面板順坡向應(yīng)力的分布

3.2 考慮溫度荷載作用時(shí)面板的應(yīng)力變形特性

混凝土面板的撓度主要受水壓力與堆石體變形的影響，溫度荷載引起的面板撓度微小，可以忽略不計(jì)。因此，本研究重點(diǎn)分析溫度荷載對(duì)面板順坡向應(yīng)力的作用。根據(jù)公伯峽水庫(kù)水位實(shí)測(cè)結(jié)果可知，在2005年2月下旬，水庫(kù)水位基本為正常蓄水位 2 005.0 m。圖4為正常蓄水位時(shí)對(duì)應(yīng)的面板順坡向溫度應(yīng)力沿高程的分布，圖5為考慮溫度荷載作用時(shí)面板順坡向應(yīng)力的分布。

圖 4 公伯峽混凝土面板堆石壩面板順坡向溫度應(yīng)力的分布

由圖4可知，面板的溫度應(yīng)力主要表現(xiàn)為拉應(yīng)力，拉應(yīng)力最大值為1.46 MPa，出現(xiàn)在大約1/2壩高處的面板表面。高程1 990.0 m以下面板的順坡向溫度應(yīng)力值均大于1.0 MPa；當(dāng)高程大于1 990.0 m時(shí)，接近壩頂?shù)拿姘屙樒孪驕囟葢?yīng)力呈現(xiàn)出迅速降低趨勢(shì)，這種變化趨勢(shì)主要受當(dāng)時(shí)氣溫、水溫的影響。同一高程處面板表面點(diǎn)、中心點(diǎn)及底面點(diǎn)的溫度應(yīng)力值略有差異，其中表面點(diǎn)的溫度應(yīng)力最大，底面點(diǎn)的溫度應(yīng)力最小。

圖 5 考慮溫度荷載時(shí)公伯峽混凝土堆石壩面板順坡向應(yīng)力的分布

從圖5可以看出，面板的順坡向應(yīng)力分布呈現(xiàn)出中部受壓、兩端受拉的變化趨勢(shì)。面板底面點(diǎn)與表面點(diǎn)的應(yīng)力分布基本一致，在同一高程處，面板底面點(diǎn)的應(yīng)力較表面點(diǎn)大。面板的順坡向壓應(yīng)力絕對(duì)值最大為7.9 MPa，出現(xiàn)在約1/2壩高處的面板底面，面板底部與頂部拉應(yīng)力的數(shù)值相差不大，約為1.3 MPa。與不考慮溫度荷載時(shí)相比，面板順坡向應(yīng)力分布規(guī)律變化不大，但在數(shù)值上有明顯增大，主要是受到累積溫度拉應(yīng)力的影響所致。

3.3 面板應(yīng)力變形計(jì)算結(jié)果與實(shí)測(cè)結(jié)果的對(duì)比分析

公伯峽面板堆石壩布設(shè)有三向應(yīng)變計(jì)和無(wú)應(yīng)力計(jì)組合用來(lái)觀測(cè)面板的應(yīng)力應(yīng)變，采用面板F9所在壩段的面板撓度、順坡向應(yīng)力實(shí)測(cè)數(shù)據(jù)(圖6和圖7)，與本研究計(jì)算結(jié)果進(jìn)行對(duì)比分析。

圖 6 不同高程下公伯峽混凝土面板堆石壩實(shí)測(cè)面板撓度的歷時(shí)變化

圖 7 不同高程下公伯峽混凝土面板堆石壩實(shí)測(cè)面板順坡向應(yīng)力的歷時(shí)變化

由圖6可知，公伯峽混凝土面板堆石壩實(shí)測(cè)撓度最大值為22.5 cm，本研究仿真計(jì)算的面板撓度最大值為24.0 cm，與實(shí)測(cè)結(jié)果基本一致。

從圖7可以看出，在運(yùn)行期，公伯峽混凝土面板堆石壩面板底端與頂端出現(xiàn)了拉應(yīng)力，底端平均拉應(yīng)力約2 MPa，面板頂部平均拉應(yīng)力約1 MPa；面板中部呈受壓狀態(tài)，壓應(yīng)力絕對(duì)值的最大值約為7.5 MPa，位于1 902.622 m高程。本研究計(jì)算的考慮溫度荷載時(shí)面板順坡向應(yīng)力也呈現(xiàn)出中部受壓、兩端受拉的變化趨勢(shì)，中部壓應(yīng)力絕對(duì)值的最大值為7.9 MPa，位于約1 925 m高程；底部最大拉應(yīng)力約為1.4 MPa，頂部最大拉應(yīng)力約為1.2 MPa。對(duì)比可知，計(jì)算結(jié)果與實(shí)測(cè)結(jié)果基本一致，可見(jiàn)在考慮溫度荷載時(shí)，混凝土面板的應(yīng)力變形數(shù)值計(jì)算結(jié)果將更為合理與準(zhǔn)確。

4 結(jié) 論

本研究基于接觸摩擦單元數(shù)值計(jì)算模型，考慮溫度荷載的作用，以公伯峽混凝土面板堆石壩為例，進(jìn)行了混凝土面板應(yīng)力變形有限元仿真計(jì)算，并將計(jì)算結(jié)果與實(shí)測(cè)結(jié)果進(jìn)行了對(duì)比分析，研究結(jié)果表明：1)在水壓力與溫度荷載共同作用下，面板撓度呈現(xiàn)出中部大、兩端小的分布規(guī)律；面板的順坡向應(yīng)力分布呈現(xiàn)出中部受壓、兩端受拉的變化趨勢(shì)；面板底面與表面的應(yīng)力分布一致，在同一高程處，面板底面點(diǎn)的應(yīng)力比表面點(diǎn)大。2)考慮溫度荷載作用時(shí)，混凝土面板的順坡向應(yīng)力分布規(guī)律與不考慮溫度荷載作用時(shí)基本一致，但拉應(yīng)力值明顯增大，并且與實(shí)測(cè)結(jié)果更為接近。3)本研究未考慮混凝土面板可能出現(xiàn)的脫空現(xiàn)象，關(guān)于面板脫空情況下的應(yīng)力變形特性尚有待進(jìn)一步研究。

[參考文獻(xiàn)]

[1] 蔣國(guó)澄,傅志安,鳳家驥.混凝土面板壩工程 [M].武漢:湖北科學(xué)技術(shù)出版社,1997.

Jiang G C,Fu Z A,Feng J J.Concrete faced rock-fill dam engineering [M].Wuhan:Hubei Science and Technology Press,1997.(in Chinese)

[2] Cooke J B.Development in high concrete face rock-fill dams [J].Hydropower & Dams,1997,4(4):69-73.

[3] 馬洪琪.300 m級(jí)面板堆石壩適應(yīng)性研究及對(duì)策研究 [J].中國(guó)工程科學(xué),2011,13(12):4-8.

Ma H Q.300 m grade concrete faced rock-fill dam adaptability and countermeasures [J].China Engineering Science,2011,13(12):4-8.(in Chinese)

[4] 酈能惠,楊澤燕.中國(guó)混凝土面板堆石壩的技術(shù)進(jìn)步 [J].巖土工程學(xué)報(bào),2012,34(8):1361-1368.

Li N H,Yang Z Y.Technical advances in concrete face rock-fill dams in China [J].Chinese Journal of Geotechnical Engineering,2012,34(8):1361-1368.(in Chinese)

[5] 徐澤平.混凝土面板堆石壩應(yīng)力變形特性研究 [M].鄭州:黃河水利出版社,2005.

Xu Z P.Study on characteristics of stress and deformation of concrete face rock-fill dam [M].Zhengzhou:The Yellow River Water Conservancy Press,2005.(in Chinese)

[6] 張國(guó)新,厲易生.堆石壩面板收縮性貫穿裂縫的理論分析及防裂措施 [J].水力發(fā)電學(xué)報(bào),2005,24(3):30-33.

Zhang G X,Li Y S.Faced rockfill dam:Theoretical analysis of through shrinkage cracks and prevention measures [J].Journal of Hydroelectric Engineering,2005,24(3):30-33.(in Chinese)

[7] 王瑞駿,王黨在,陳堯隆.混凝土面板堆石壩施工期面板溫度應(yīng)力仿真分析 [J].西北農(nóng)林科技大學(xué)學(xué)報(bào)：自然科學(xué)版,2004,32(10):123-126.

Wang R J,Wang D Z,Chen Y L.Simulation analysis on slabs thermal stress of CFRD in construction period [J].Journal of Northwest A&F University：Natural Science Edition，2004,32(10):123-126.(in Chinese)

[8] 王 勇,殷宗澤.面板堆石壩堆石流變對(duì)面板應(yīng)力變形的影響分析 [J].河海大學(xué)學(xué)報(bào),2000,28(6):5-10.

Wang Y,Yin Z Z.The rockfill creep on stress and deformation analysis of the effect of panel [J].Journal of Hohai University,2000,28(6):5-10.(in Chinese)

[9] Maranda das Neves E.Advances in rockfill structure [M].London:Kluwer Academic Publishers，1991．

[10] 李守義,馬成成,李炎隆,等.混凝土面板開(kāi)裂情況下堆石壩應(yīng)力場(chǎng)與滲流場(chǎng)耦合分析 [J].西北農(nóng)林科技大學(xué)學(xué)報(bào)：自然科學(xué)版,2012,40(3):206-212.

Li S Y,Ma C C,Li Y L,et al.Coupling analysis of stress field and seepage field for concrete face rock-fill dam with slabs cracking [J].Journal of Northwest A&F University：Natural Science Edition，2012,40(3):206-212.(in Chinese)

[11] Yalin A.Investigation of the cracking of CFRD face plates [J].Computers and Geotechnics,2011,38:905-916.

[12] Alemdar B,Murat E K,Suleyman A.The effect of concrete slab-rockfill interface behavior on the earthquake performance of a CFR dam [J].International Journal of Non-Linear Mechanics,2011,46:35-46.

[13] Alemdar B,Murat E K.Linear and nonlinear response of concrete slab on CFR dam during earthquake [J].Soil Dynamics and Earthquake Engineering,2010,30:990-1003.

[14] Zhang G,Zhang J M.Numerical modeling of soil-structure interface of a concrete-faced rockfill dam [J].Computers and Geotechnics,2009,36:762-772.

[15] 王瑞駿,李章浩,王黨在,等.面板與墊層之間接觸面的溫度應(yīng)力計(jì)算模型研究 [J].水力發(fā)電學(xué)報(bào),2006,25(3):58-61.

Wang R J,Li Z H,Wang D Z,et al.Research on computation model of thermal stresses of interface between slabs and cushion layer [J].Journal of Hydroelectric Engineering,2006,25(3):58-61.(in Chinese)

[16] 雷曉燕,G Swoboda,杜慶華.接觸摩擦單元的理論及其應(yīng)用 [J].巖土工程學(xué)報(bào),1994,16(3):23-32.

Lei X Y,Swoboda G,Du Q H.Theory and application of contact-friction interface element [J].Chinese Jounal of Geotechnical Engineering,1994,16(3):23-32.(in Chinese)

[17] 國(guó)家電力公司西北勘測(cè)設(shè)計(jì)研究院.黃河公伯峽水電站工程混凝土面板堆石壩設(shè)計(jì)說(shuō)明 [R].西安:國(guó)家電力公司西北勘測(cè)設(shè)計(jì)研究院,2001.

Northwest Institute of Survey and Design of Hydropower Engineering.Design notes of Yellow River Gongboxia hydropower project CFRD [R].Xi’an:Northwest Institute of Survey and Design of Hydropower Engineering,2001.(in Chinese)