周祥志,陳 昊

(四川省水利水電勘測(cè)設(shè)計(jì)研究院,成都,610072)

1 概述

俄公堡水電站工程在四川省涼山州木里縣境內(nèi),位于雅礱江一級(jí)支流木里河干流上,是木里河干流(上通壩～阿布地)梯級(jí)開(kāi)發(fā)的第四級(jí)電站。該電站為閘壩引水式開(kāi)發(fā),引水隧洞長(zhǎng)14853.9m,電站額定水頭77.5m,額定引用流量191.1m3/s。電站裝機(jī)3臺(tái),單機(jī)容量44MW,總裝機(jī)容量132MW。工程靜態(tài)總投資162550萬(wàn)元,單位千瓦投資12314元/kW,單位電能投資2.578元/kW?h,引水隧洞工程投資占總投資的38.7%。實(shí)施階段圍巖類別以Ⅳ、Ⅴ類為主,在保證工程安全的前提下,為控制投資優(yōu)化,引水隧洞永久襯砌結(jié)構(gòu)的配筋計(jì)算采取了常用的邊值法、公式法及西安理工大學(xué)的新型有限元法進(jìn)行對(duì)比分析,并提出在實(shí)施階段適合工程實(shí)際情況的計(jì)算成果。

2 “邊值法”與“彈性力學(xué)法”計(jì)算分析

“邊值法”與“彈性力學(xué)法”的計(jì)算結(jié)果十分相近。“彈性力學(xué)法”和“邊值法”相同之處是:都假定混凝土襯砌是一個(gè)完整的厚壁圓筒,可以承擔(dān)拉應(yīng)力?！皬椥粤W(xué)法”和“邊值法”的不同之處是:前者是建立在彈性力學(xué)基礎(chǔ)上的解析解;后者是建立在結(jié)構(gòu)力學(xué)基礎(chǔ)上的微分方程組,通過(guò)龍格-庫(kù)塔法逐次漸進(jìn)地計(jì)算出圍巖抗力、混凝土襯砌的軸力、彎矩和切力。因而,它們計(jì)算成果之間略有差異。

不管采用“邊值法”還是“彈性力學(xué)法”,混凝土襯砌內(nèi)、外側(cè)應(yīng)力,都大于混凝土拉應(yīng)力的允許值;內(nèi)、外側(cè)應(yīng)變也都大于混凝土極限拉伸值,襯砌必定開(kāi)裂?？傊?圓形有壓水工隧洞在大多情況下,一般都要開(kāi)裂。此外,隧洞的混凝土襯砌,在圍巖的強(qiáng)約束下,除了承擔(dān)均勻內(nèi)水壓力、自重、山巖壓力、滿洞水重等荷載外,還要承擔(dān)溫度應(yīng)力、干縮應(yīng)力、施工冷縫、開(kāi)挖面不良引起的應(yīng)力集中和養(yǎng)護(hù)不良等其他因素的影響。所以,以前有“逢洞必裂”一說(shuō)。在《水工隧洞設(shè)計(jì)規(guī)范》DL/T5195-2004的條文說(shuō)明中(127頁(yè)),也有“根據(jù)工程實(shí)際和大量試驗(yàn)資料表明,按抗裂設(shè)計(jì)由于混凝土抗拉強(qiáng)度設(shè)計(jì)值較低,即使襯砌設(shè)計(jì)的很厚,仍然控制不了裂縫出現(xiàn)”。

隧洞混凝土襯砌開(kāi)裂以后,厚壁圓筒已經(jīng)變成幾個(gè)弧形“瓦片”,其結(jié)構(gòu)與受力條件,當(dāng)然不能稱為“厚壁圓筒”,但在“邊值法”配筋計(jì)算中,仍按混凝土襯砌未開(kāi)裂的“厚壁圓筒”計(jì)算,其結(jié)果必然是配筋量很大(因?yàn)橐r砌承擔(dān)的荷載占內(nèi)水壓力的40%以上)。這既不符合混凝土襯砌開(kāi)裂后實(shí)際情況,又造成不必要的浪費(fèi),同時(shí)也與“新奧法”以圍巖為承載主體的思路不合。所以,如果隧洞襯砌是用“鋼”或者其他抗拉強(qiáng)度很高的物質(zhì)做成的,或者水頭很低,混凝土襯砌不會(huì)開(kāi)裂時(shí),用“邊值法”計(jì)算的結(jié)果肯定正確。但當(dāng)厚壁圓筒已經(jīng)開(kāi)裂,變成幾個(gè)弧形“瓦片”之后,再按“邊值法”計(jì)算就不合適了。

3 公式法計(jì)算分析

“公式法”的前提是:首先“承認(rèn)”混凝土襯砌已經(jīng)開(kāi)裂,混凝土襯砌不承擔(dān)任何拉力,在此前提下推導(dǎo)出計(jì)算圍巖變形、混凝土襯砌內(nèi)層鋼筋(或外層鋼筋)的徑向位移、應(yīng)變、應(yīng)力和鋼筋面積公式。

“公式法”的基本假定,較“彈性力學(xué)法”和“邊值法”更要符合水工隧洞鋼筋混凝土襯砌受力的實(shí)際狀況。但問(wèn)題在于:不管圍巖實(shí)際位移大小,一律按鋼筋的允許應(yīng)力[σs]代替鋼筋的實(shí)際應(yīng)力[σ],計(jì)算所需的鋼筋面積,有時(shí)會(huì)得到“負(fù)”值。其實(shí),在大部分情況下,鋼筋應(yīng)力是達(dá)不到鋼筋的允許應(yīng)力[σs]的。所以,有人研究了改進(jìn)方法。現(xiàn)有改進(jìn)方法可歸結(jié)為兩類,但其對(duì)鋼筋應(yīng)力的假設(shè)仍然存在問(wèn)題。第一類是當(dāng)按最小配筋率配筋計(jì)算出鋼筋應(yīng)力大于允許應(yīng)力時(shí),又回到了原來(lái)的公式法。第二類是先計(jì)算出混凝土應(yīng)變,根據(jù)混凝土應(yīng)變來(lái)計(jì)算鋼筋應(yīng)力,據(jù)此對(duì)襯砌配筋。該法忽略了裂縫處鋼筋對(duì)混凝土的錨栓作用,計(jì)算結(jié)果偏于危險(xiǎn)。

在“公式法”中,圍巖力學(xué)性能用彈性抗力來(lái)表征,考慮襯砌混凝土的開(kāi)裂,經(jīng)過(guò)推導(dǎo)得出計(jì)算公式,這對(duì)水工高壓隧洞襯砌混凝土的考慮比較符合實(shí)際。但對(duì)于水工高壓隧洞洞周圍巖,肯定存在一定范圍塑性區(qū)。眾所周知,當(dāng)圍巖進(jìn)入塑性以后,變形模量會(huì)大幅降低,圍巖分擔(dān)荷載也會(huì)隨之降低,使得襯砌要分擔(dān)更大的荷載。而在“公式法”中,圍巖力學(xué)性能用彈性抗力來(lái)表征,就相當(dāng)于把圍巖按完全彈性的材料來(lái)考慮,那么我們用彈性抗力來(lái)算得的結(jié)果,就不能反映真實(shí)的情況。綜上所述,“公式法”有以下兩點(diǎn)不足:(1)圍巖按彈性計(jì)算;(2)鋼筋應(yīng)力按變形平均值給出,不符合在襯砌裂開(kāi)處局部增大情況。

4 新型有限元法計(jì)算分析

4.1 數(shù)值分析模型的建立

許多的數(shù)值分析算法中,雖然圍巖按彈塑性考慮,但襯砌在大多時(shí)候仍采用梁?jiǎn)卧M。梁?jiǎn)卧且环N彈性單元,其剛度矩陣計(jì)算以橫梁?jiǎn)卧碚摓橐罁?jù),在此考慮所有的斷面參數(shù),如基準(zhǔn)面、剪切面和慣性力矩。其基本假定和“邊值法”的基本假定相似,其不合理性也就出現(xiàn)了。

在認(rèn)識(shí)到各方法的不足后,考慮到圍巖和襯砌等材料在水工高壓隧洞實(shí)際工程表現(xiàn)出性狀,大多情況下不是彈性體,而是彈塑性體,并且大多情況下混凝土襯砌會(huì)開(kāi)裂,即襯砌和圍巖是有條件的聯(lián)合作用。鋼筋能夠承受環(huán)向拉力,但不能承受彎矩。在襯砌開(kāi)裂處,由于混凝土已經(jīng)開(kāi)裂,不再承受拉應(yīng)力,襯砌所受拉力全部由鋼筋承擔(dān)。

考慮到以上問(wèn)題后,我院委托西安理工大學(xué)采用從奧地利引進(jìn)的大型巖土工程數(shù)值仿真分析軟件FINAL,以及西安理工大學(xué)自行研制開(kāi)發(fā)的大型巖土工程軟件包ROCKS進(jìn)行數(shù)值分析。圍巖和襯砌混凝土按彈塑性考慮,并考慮襯砌開(kāi)裂影響。圍巖和襯砌混凝土采用三角形6結(jié)點(diǎn)等參實(shí)體單元模擬;在裂縫處的接觸面模擬,采用FINAL有限元分析程序獨(dú)有的單元——COJO單元;鋼筋用FINAL有限元分析程序獨(dú)有的一種只能承受拉應(yīng)力,不能承受彎矩和壓應(yīng)力的梁?jiǎn)卧狟EAm3單元模擬。襯砌和鋼筋計(jì)算模型見(jiàn)圖1。

圖1 襯砌和鋼筋計(jì)算模型

4.2 邊界條件

在隧洞的埋深足夠大(可認(rèn)為隧洞埋深大于內(nèi)水水頭高度)時(shí),假定模型底側(cè)為不可動(dòng),即底側(cè)面的水平位移和豎直位移是完全約束的;模型的兩側(cè)為上下可動(dòng)但左右不可動(dòng),即兩側(cè)邊界的水平位移是約束的,豎直方向是自由的;模型頂部約束為上覆巖體壓力。

4.3 最大裂縫條數(shù)確定

4.3.1 根據(jù)鋼筋的許用應(yīng)變確定。Ⅱ級(jí)鋼筋的許用應(yīng)力[σs]=229.63MPa,相應(yīng)的應(yīng)變值為[εs]=11.4815×10-4。

4.3.2 根據(jù)鋼筋混凝土中兩條裂縫間的最小距離確定。

鋼筋混凝土中相鄰裂縫最小間距為:L=3c+0.1d/ρte

其中,配筋率ρte=As/Ate,當(dāng)ρte＜0.03時(shí),取ρte=0.03;Ate=2asb。

4.4 最大裂縫寬度確定

鋼筋混凝土襯砌限裂驗(yàn)算規(guī)定,最大裂縫寬度為0.25mm。

5 幾種計(jì)算方法配筋結(jié)果對(duì)比

表1 Ⅳ類圍巖區(qū)配筋結(jié)果對(duì)比

表2 Ⅴ類圍巖區(qū)配筋結(jié)果對(duì)比

以俄公堡電站Ⅳ類圍巖區(qū)配筋結(jié)果對(duì)比表(表1)、Ⅴ類圍巖區(qū)配筋結(jié)果對(duì)比表(表2)為例,經(jīng)過(guò)對(duì)比可以看出,按照“邊值法”和“彈性力學(xué)法”計(jì)算結(jié)果非常大。原因是襯砌按厚壁圓筒考慮,沒(méi)有考慮到襯砌開(kāi)裂?！肮椒捌涓倪M(jìn)方法”考慮到了襯砌開(kāi)裂,但沒(méi)考慮到裂縫處鋼筋應(yīng)力局部增大的情況,計(jì)算結(jié)果則非常小。有限元法克服了以上兩種方法的缺點(diǎn),使計(jì)算邊界條件與襯砌實(shí)際受力情況更吻合。實(shí)施階段采用有限元法配筋成果,局部斷層破碎帶進(jìn)行加強(qiáng),采用有限元法的計(jì)算成果對(duì)控制整個(gè)工程投資起到了關(guān)鍵性的作用。