程 丹，蘇 凱，石怡安

(武漢大學水資源與水電工程科學國家重點實驗室，湖北武漢 430072)

隨著高水頭電站以及抽水蓄能電站的大量興建，為降低工程造價，減小土石方的開挖，常采用聯(lián)合供水和分組供水的方案，這樣就需要采用洞徑較大的發(fā)電引水隧洞和分岔管。由于岔管位于壓力管道的末端并靠近廠房，內(nèi)水壓力較大，傳統(tǒng)的鋼板襯砌結(jié)構(gòu)在施工中存在洞內(nèi)運輸、拼接、焊接等方面難題，且用鋼量較大，為降低工程造價及施工難度，實際工程中產(chǎn)生了鋼筋混凝土襯砌結(jié)構(gòu)這一需求。鋼筋混凝土岔管的體型、荷載、邊界以及圍巖聯(lián)合承載的機理極其復(fù)雜，傳統(tǒng)的結(jié)構(gòu)力學法因其計算條件和計算方法不完善，無法準確描述結(jié)構(gòu)的變形和受力特征［1-2］。文獻［1-9］采用有限單元數(shù)值計算法對鋼筋混凝土岔管圍巖穩(wěn)定及岔管結(jié)構(gòu)進行分析，指出圍巖開挖后塑性區(qū)主要位于主洞、岔洞交叉處，由于隧洞體型連續(xù)性遭到破壞，塑性區(qū)開展深度較大［9］，是施工中洞室成形難度較大的部位，除了對這些部位及時噴錨支護外，同時建議對岔襠處尖角以及主管與支管轉(zhuǎn)折處進行修圓處理［3，8］。肖明等［10-11］對鋼筋混凝土岔管進行了優(yōu)化分析，指出岔角結(jié)構(gòu)修圓處理對岔角處的應(yīng)力集中有較大影響，能有效地緩解應(yīng)力集中現(xiàn)象，但以上研究均未對修圓處理前后的圍巖穩(wěn)定性和襯砌結(jié)構(gòu)安全進行系統(tǒng)的對比分析。本文應(yīng)用Ansys前處理模塊，建立了鋼筋混凝土岔管數(shù)值模型，并對鋼筋混凝土岔管襯砌銳角區(qū)進行修圓處理，采用有限單元法開展了圍巖穩(wěn)定性及襯砌結(jié)構(gòu)應(yīng)力變形的計算分析，系統(tǒng)地分析了銳角區(qū)不同半徑修圓處理對圍巖穩(wěn)定及岔管襯砌結(jié)構(gòu)受力的影響。

1 工程概況

某抽水蓄能電站裝機容量2400 MW，安裝6臺單機容量為400 MW的水泵水輪機組，引水系統(tǒng)采用一洞三機供水方式，主管直徑5.5 m，支管直徑3.0 m。引水岔管段圍巖大部分為Ⅱ類，局部為Ⅲ類。鋼筋混凝土岔管設(shè)計采用主、支管圓弧的最低點置于同一水平面上，岔管段襯砌厚度為1m，主管與支管段襯砌混凝土厚度為0.6 m，外表面銳角修圓半徑分別采用1m和2m兩種情況(外表面銳角修圓半徑是指岔襠銳角區(qū)襯砌外表面腰部圓弧半徑)，分別對應(yīng)修圓1和修圓2，引水岔管結(jié)構(gòu)圖見圖1。

圖1 岔管結(jié)構(gòu)示意圖 (單位:mm)

根據(jù)岔管結(jié)構(gòu)尺寸建立了三維有限元模型，模型的頂部、底部及兩側(cè)圍巖均取5倍以上的引水主管直徑，整體模型計算范圍和修圓后的引水岔管襯砌網(wǎng)格及局部放大圖見圖2。引水岔管處完整圍巖力學參數(shù)按Ⅱ2類圍巖選用，各材料參數(shù)見表1和表2。采用笛卡兒直角坐標系，其整體坐標系的x軸與主管管道的軸線一致，以指向下游為正;鉛直方向為z軸，以向上為正，y軸以右手法則確定。引水岔管坐標原點位于1號岔管分岔點上游6 m管道中心處。模型左右兩側(cè)、上下游側(cè)和底部施加法向位移約束，頂部為自由面。岔管結(jié)構(gòu)上的荷載包括襯砌自重、內(nèi)水壓力、外水壓力。

圖2 引水岔管網(wǎng)格 (單位:m)

表1 巖體材料物理力學參數(shù)

表2 襯砌材料力學參數(shù)

非線性分析中，巖體材料采用理想彈塑性本構(gòu)關(guān)系和 Drucker-Prager屈服準則［12］;岔管襯砌混凝土采用多線性彈性本構(gòu)關(guān)系，并取組合破壞準則作為開裂準則;鋼筋采用分布鋼筋模式。在混凝土開裂前其應(yīng)力應(yīng)變關(guān)系采用多線性彈性關(guān)系，通過破壞準則(采用組合破壞準則，三向受拉應(yīng)力狀態(tài)時取最大主應(yīng)力準則，三向受壓應(yīng)力狀態(tài)取William-Warnke五參數(shù)準則［13］，其他應(yīng)力狀態(tài)取兩準則的過渡形式)來判斷是否達到破壞曲面，當達到破壞曲面時，按拉壓不同破壞形式相應(yīng)改變應(yīng)力應(yīng)變關(guān)系;受拉開裂后，混凝土應(yīng)力應(yīng)變矩陣沿著破壞面和垂直于破壞面的方向建立，并設(shè)置相應(yīng)參數(shù)反映混凝土開裂后的應(yīng)力應(yīng)變關(guān)系。材料在受壓破壞后所有方向發(fā)生應(yīng)變軟化，單元完全喪失承載力。

2 施工期圍巖穩(wěn)定性分析

巖體鉛直方向的應(yīng)力符合自重應(yīng)力場規(guī)律，水平向應(yīng)力側(cè)壓力系數(shù)為1.2～1.3，施工開挖均按一次性開挖考慮。

a．施工開挖期間圍巖第三主應(yīng)力σ3等值線見圖3。從圖3可以看出，1號和2號岔管段圍巖的應(yīng)力均為壓應(yīng)力，最小壓應(yīng)力出現(xiàn)在1號岔管分岔處圍巖的頂部和底部，最大壓應(yīng)力均出現(xiàn)在靠近主管與岔管以及岔管與支管的交界處開挖面的頂部和底部。施工期岔管段開挖圍巖應(yīng)力與鉛直方向的位移見表3，由表3可以看出，3種方案的圍巖應(yīng)力數(shù)值比較接近。

圖3 施工開挖期圍巖σ3等值線 (單位:MPa)

b．圍巖開挖后形成臨空面，洞周巖體向洞內(nèi)變形，其中鉛直方向的位移最為明顯，最大位移出現(xiàn)在岔管底部。修圓后的岔管段圍巖鉛直方向位移有所增大，且隨著修圓半徑的增大，岔管段圍巖鉛直向位移逐漸增大，如1號岔管處圍巖的位移由5.55mm增加到7.38 mm，增大了33.0%。

表3 施工期岔管段開挖圍巖應(yīng)力與位移分析

c．圍巖開挖后，圍巖的塑性區(qū)集中在岔襠和1號岔管洞室頂部、底部，深度在1.0～1.5 m之間，見圖4。銳角區(qū)修圓后，圍巖的塑性集中程度減緩。由于分岔處跨度(最大開挖橫斷面的寬度)加大，圍巖塑性區(qū)范圍增加，且隨著修圓半徑的增大，塑性區(qū)范圍增大，如圖5所示(圖中橫坐標為底數(shù)取10的對數(shù)坐標)。

圖5 不同修圓條件下施工開挖期間塑性區(qū)范圍

3 運行期襯砌結(jié)構(gòu)分析

3．1 運行期襯砌結(jié)構(gòu)線性分析

運行期內(nèi)水壓力為3.85 MPa，在內(nèi)水壓力作用下，岔管襯砌主要承受拉應(yīng)力，岔襠出現(xiàn)了明顯的應(yīng)力集中現(xiàn)象。由于運行期2號岔管襯砌拉應(yīng)力較1號岔管小，受力規(guī)律與1號岔管相同，表4僅給出了1號岔管襯砌應(yīng)力與位移。由表4可知，修圓后岔管襯砌結(jié)構(gòu)最大拉應(yīng)力減小，岔襠和岔頂應(yīng)力變化趨向平緩，應(yīng)力集中程度有所緩解。相對修圓1，修圓2修圓半徑增大，岔管應(yīng)力集中程度惡化，故需要選擇適當?shù)男迗A半徑。襯砌主要向洞外變形，最大值均出現(xiàn)在岔管分岔處的正上方頂板處，整體合位移不大，襯砌合位移最大值隨著修圓半徑的增大而增大。

表4 運行期1號岔管襯砌結(jié)構(gòu)應(yīng)力與位移分析

3．2 運行期襯砌結(jié)構(gòu)非線性分析

襯砌混凝土材料為非線性材料，襯砌主管、支管和岔管環(huán)向鋼筋采用雙層φ32mm@200mm，軸向鋼筋采用雙層φ25 mm@250 mm，計算時采用線性加載方式，各內(nèi)水壓力作用下的內(nèi)外層環(huán)向鋼筋最大值見圖6。可以看出鋼筋應(yīng)力隨著內(nèi)水壓力的增加而增大，但鋼筋應(yīng)力的增大與水頭的增加不成比例。

圖6 運行期1號、2號岔管環(huán)向鋼筋應(yīng)力分析

對于未修圓體型，當襯砌所受內(nèi)水水頭為33 m時，內(nèi)外層環(huán)向鋼筋應(yīng)力均較小，襯砌未開裂;內(nèi)水水頭增加到55 m時，環(huán)向鋼筋應(yīng)力最大值為10.17MPa，但此時2號岔管襯砌內(nèi)表面銳角區(qū)已開裂，隨著水頭的增加，1號岔管襯砌銳角區(qū)開裂，開裂范圍向襯砌腰部和外表面延伸;當內(nèi)水水頭達到165 m時，襯砌內(nèi)層單元除分岔處頂部和底部未開裂，其他絕大部分已開裂，外層單元僅腰部開裂，其中1號岔管與主管交界處3 m內(nèi)襯砌均開裂，開裂深度為0.8～1.0 m，2號岔管腰部開裂深度為0.4～0.6 m。隨著水頭的繼續(xù)增加，開裂范圍進一步擴展，當水頭達到385 m時，僅岔管頂部和底部未開裂，其他已基本全裂。修圓后的開裂規(guī)律與未修圓時類似，而對于修圓1和修圓2體型，當水頭為55m時，兩岔管襯砌內(nèi)層單元銳角修圓區(qū)均開裂，水頭增加到385m時，未開裂部位集中在分岔處頂部和底部，開裂范圍與未修圓時相差不大，如圖7所示。

對于修圓1體型，總體上1號岔管和2號岔管內(nèi)層環(huán)向鋼筋應(yīng)力較未修圓時減小，1號岔管外層鋼筋應(yīng)力增大，2號岔管外層鋼筋應(yīng)力變化幅度很小，襯砌鋼筋應(yīng)力分布更加均勻。比較修圓2和修圓1兩種修圓體型的鋼筋應(yīng)力值，可以發(fā)現(xiàn)前者的內(nèi)層鋼筋應(yīng)力大于后者，而前者的外層鋼筋應(yīng)力小于后者，且前者鋼筋均勻性較后者差。修圓后鋼筋整體應(yīng)力有所減小，應(yīng)力分布更加均勻，對襯砌受力是有利的，但修圓半徑并非越大越好。

圖7 運行期岔管襯砌開裂情況(水頭385 m)

4 檢修期襯砌結(jié)構(gòu)線性分析

檢修期外水壓力為1.11 MPa，作用在襯砌的外表面，岔管襯砌主要承受壓應(yīng)力，局部出現(xiàn)拉應(yīng)力。未修圓時，拉應(yīng)力集中在1號岔管岔襠外表面和岔管襯砌頂部;修圓后拉應(yīng)力數(shù)值和集中范圍略有增大，主要集中在岔管頂部，如圖8和圖9所示。

圖8 檢修期岔管襯砌σ1等值線 (單位:MPa)

圖9 檢修期岔管襯砌σ3等值線 (單位:MPa)

由于檢修期2號岔管襯砌壓應(yīng)力較1號岔管小，受力規(guī)律與1號岔管相同，表5僅給出了1號岔管襯砌應(yīng)力與位移。由表5可知:修圓1相對未修圓體型，1號岔管岔襠、岔頂內(nèi)外表面應(yīng)力差減小，應(yīng)力集中程度也減小，而修圓2相對未修圓體型，1號岔襠、岔頂內(nèi)外層應(yīng)力差反而增大，即修圓半徑過大，應(yīng)力集中并未得到緩解。襯砌主要向洞內(nèi)變形，修圓后合位移增大，且隨著修圓半徑的增大，合位移增大。

表5 檢修期1號岔管襯砌應(yīng)力與位移分析

5 結(jié)語

由于岔管銳角區(qū)存在不連續(xù)面，因而圍巖穩(wěn)定性與襯砌結(jié)構(gòu)受力情況復(fù)雜。從水流條件看，未進行修圓處理時，銳角區(qū)岔襠襯砌可起到類似導(dǎo)流板的作用，對水流條件是有利的;對銳角區(qū)進行修圓處理后，對岔襠處的水流條件是不利的，但對改善襯砌受力是有利的。本文對岔管銳角區(qū)進行了不同程度的修圓，并從結(jié)構(gòu)受力方面進行了施工、運行、檢修期的對比分析，根據(jù)計算分析結(jié)果可以得出:①施工開挖期，岔襠處應(yīng)力集中較為明顯，塑性區(qū)集中在岔襠和岔管分岔處的頂部、底部。修圓后，塑性集中程度降低，由于岔角處跨度加大，開挖對圍巖擾動增大，圍巖塑性區(qū)范圍增加，變形略有增大。②運行期襯砌主要承受拉應(yīng)力，適當?shù)男迗A處理可使襯砌應(yīng)力分布更加均勻，襯砌開裂區(qū)范圍減小，鋼筋應(yīng)力有所降低，對襯砌受力是有利的。但襯砌合位移隨著修圓半徑的增大而增大，認為緩變的岔角結(jié)構(gòu)更合理，更有利于襯砌結(jié)構(gòu)的安全。③檢修期間襯砌主要承受壓應(yīng)力，岔襠和岔管頂部出現(xiàn)了局部拉應(yīng)力集中。適當?shù)男迗A處理后，拉應(yīng)力集中范圍和數(shù)值有所增大，但岔襠壓應(yīng)力集中現(xiàn)象有所改善，對整體的變形影響相對較小。因而在岔襠外圍圍巖穩(wěn)定的情況下盡量減少開挖，選擇適當?shù)男迗A半徑，有利于圍巖的穩(wěn)定性和襯砌結(jié)構(gòu)的受力。

［1］黃麗偉，鄧柏旺，梁幼．地下鋼筋混凝土岔管三維有限元分析［J］．華北水利水電學院學報，2007，28(2):40-42．(HUANG Liwei，DENG Bowang，LIANG You．3D finite element analysis of underground reinforced concrete bifurcation［J］．Journal of North China Institute of Water Conservancy and Hydroelectric Power，2007，28(2):40-42．(in Chinese))

［2］韓前龍．大型地下鋼筋混凝土岔管結(jié)構(gòu)研究［D］．武漢:武漢大學，2005．

［3］韓前龍，伍鶴皋，蘇凱，等．地下鋼筋混凝土岔管應(yīng)力分析［J］．武漢大學學報:工學版，2005，38(3):31-35．(HAN Qianlong，WU Hegao，SU Kai，et al．Stress analysis of underground reinforced concrete bifurcation［J］．Engineering Journal of Wuhan University:Engineering Edition，2005，38(3):31-35．(in Chinese))

［4］蘇凱，伍鶴皋，李青麒．抽水蓄能電站鋼筋混凝土分岔管有限元分析［J］．水電能源科學，2005，23(5):51-54．(SU Kai，WU Hegao，LI Qingqi．Finite element method analysis of reinforced concrete manifold of pumped storage power station ［J］．Water Resources and Power，2005，23(5):51-54．(in Chinese))

［5］鄧柏旺．地下鋼筋混凝土高壓岔管圍巖穩(wěn)定與襯砌結(jié)構(gòu)分析［D］．南京:河海大學，2007．

［6］蘇凱，伍鶴皋，英鵬濤，等．斷裂影響下的高壓鋼筋混凝土岔管結(jié)構(gòu)［J］．水利水電技術(shù)，2008，39(4):41-45．(SU Kai，WU Hegao，YING Pengtao，et al．Study on highpressure reinforced concrete manifold structure under impact of fault［J］．Water Resource and Hydropower Engineering，2008，39(4):41-45．(in Chinese))

［7］任旭華，徐海奔，束加慶．地下水工隧洞圍巖及鋼筋混凝土岔管結(jié)構(gòu)數(shù)值分析［J］．四川大學學報:工程科學版，2009，41(4):8-13．(REN Xuhua，XU Haiben，SHU Jiaqing．Numerical analysis on surrounding rock and reinforced concrete bifurcation structure of undergrounding tunnels［J］．Journal of Sichuan University:Engineering Science Edition，2009，41(4):8-13．(in Chinese))

［8］王明．地下鋼筋混凝土岔管三維有限元分析研究［D］．長沙:長沙理工大學，2009．

［9］張繼勛，任旭華．地下鋼筋混凝土岔管圍巖穩(wěn)定及結(jié)構(gòu)數(shù)值分析［J］．三峽大學學報:自然科學版，2011，33(5):2-6．(ZHANGJixun，RENXuhua．Numerical analysis of surrounding rock stability and lining structure of underground reinforced concrete bifurcation ［J］．JournalofChina Three Gorges University:Natural Science，2011，33(5):2-6．(in Chinese))

［10］肖明．大型地下鋼筋混凝土岔管結(jié)構(gòu)優(yōu)化分析［J］．水利學報，2001，32(12):8-13．(XIAO Ming．Optimization analysis of large-scale underground reinforced concrete branching conduit［J］．Journal of Hydraulic Engineering，2001，32(12):8-13．(in Chinese))

［11］龔玉鋒，丁旭柳．地下鋼筋混凝土岔管結(jié)構(gòu)優(yōu)化分析［J］．水電站設(shè)計，2002，18(1):26-30．(GONG Yufeng，DING Xuliu．Structure optimization analysis on underground reinforced concrete branch pipe ［J］．Design of Hydroelectric Power Station，2002，18(1):26-30．(in Chinese))

［12］徐干成，鄭穎人． 巖石工程中屈服準則應(yīng)用的研究［J］．巖土工程學報，1990，12 ( 2 ) : 93-99． ( XU Gancheng， ZHANG Yinren． Study on application of yield criterions in geotechnical engineering ［J ］． Chinese Journal of Geotechnical Engineering，1990，12 ( 2 ) : 93-99． ( in Chinese) )

［13］董哲仁．鋼筋混凝土非線性有限元法原理與應(yīng)用［M］．北京:中國鐵道出版社，1993．