1、六庫水電站RCC縱向圍堰堰基應(yīng)力應(yīng)變及穩(wěn)定性分析             摘要采用二維和三維非線性有限元分析方法，研究了六庫水電站RCC縱向圍堰堰基應(yīng)力應(yīng)變狀態(tài)和穩(wěn)定性。在此基礎(chǔ)上對防滲墻各截面的剪力、軸力和彎矩進行計算研究，提出了配筋方案。計算表明，堰體位移大，堰體與防滲墻接頭處附近存在應(yīng)力集中現(xiàn)象，防滲墻的防滲效果顯著，分載作用明顯，右防滲墻承受的彎矩更大些，應(yīng)加大其配筋率。 關(guān)鍵詞有限元 圍堰 應(yīng)力 應(yīng)變 防滲墻 配筋率 Stress and Stability Ana

2、lysis of Longitudinal RCC Cofferdam Foundation of  Liuku Hydropower Station Abstract   The stress state and stability of longitudinal RCC cofferdam foundation of Liuku hydropower station are studied using 2D and 3D non-linear finite element method. Then the shear 、axial force and

3、 bending moment of governing sections of anti-seepage walls are studied. The result shows that the movement of cofferdam is large and the stress concentrates in joint of cofferdam and anti-seepage walls. Besides, leakage-prevention and load-carrying capacity of anti-seepage walls can meet engineerin

4、g requirement and the steel reinforcement rate of right anti-seepage wall should be increased due to its larger bending moment.Key words   finite element method，cofferdam，stress，deformation，anti-seepage walls, steel reinforcement rate 一 工程概況六庫水電站壩址區(qū)位于云南省怒江州首府六庫城上游4.5km處的怒江“U”型峽谷中，江水流向由北東轉(zhuǎn)為

5、南西流出壩址。六庫水電站為河床式電站。樞紐工程總體呈“一”字型布置，主體建筑物由混凝土重力壩、溢流壩、泄洪沖砂閘、主、副廠房等部分組成。縱向碾壓混凝土圍堰為永久建筑物，工程導(dǎo)流結(jié)束后不再拆除。縱向碾壓混凝土圍堰左側(cè)設(shè)有4孔7.5m×10.0m沖砂孔，底板高程798.0m，建基面高程794.5m。為了滿足加固邊坡、防止堰基掏刷及基礎(chǔ)防滲要求，在堰基部位設(shè)置兩道鋼筋混凝土防滲墻（連續(xù)墻）?？v向圍堰堰基座于沖積層上，沖積層厚度為15m20m，基巖為擠壓破碎巖體，厚25m30m，以散體結(jié)構(gòu)為主。沖積層滲透系數(shù)K=27.58 m/d5.69m/d,屬強中等透水性。沖積層以下為灰質(zhì)粉晶白云巖。該

6、地段分布的斷層有F4、F5、F2-1，它們和圍堰軸線成大角度斜交，呈陡傾角穿過堰基。由于特殊的地質(zhì)條件和巖層分布，六庫水電站縱向碾壓混凝土圍堰基礎(chǔ)穩(wěn)定是本工程的關(guān)鍵技術(shù)問題之一。二 有限元分析原理2.1巖體強度與本構(gòu)模型   六庫水電站縱向圍堰堰基主要由軟弱的砂卵礫石層和下伏散體結(jié)構(gòu)組成，其特點是變形和強度力學(xué)指標(biāo)低，易于發(fā)生塑性流動和屈服。為此，本次有限元分析采用低抗拉的巖體彈塑性本構(gòu)模型開展研究。   按低抗拉彈塑性模型分析，壩基巖體材料開裂條件用宏觀強度描述:表1 巖體物理力學(xué)參數(shù)表序號巖體代號變模E（GPa）泊桑比容重（t/m3）抗剪斷強度滲透系

7、數(shù)fC （MPa）K   10-7（m/s）1砂卵礫石QAL-20.0480.352.10.5250.014757.02散體結(jié)構(gòu)T2H0.0400.2802.700.3750.0578.73斷層0.0400.3502.100.3500.0350.014圍堰砼25.00.1672.451.100.9010-65鋼筋砼防滲墻26.00.1672.451.100.900.016鋼筋砼防滲隔板26.00.1672.451.100.900.017沖砂閘砼25.00.1672.451.100.9010-6825.00.1672.451.100.9010-69墊層砼25.00.1672.

8、451.100.9010-6     (i=1,2,3)式中ii表征應(yīng)力張量三個主應(yīng)力，分析中可能呈單向、雙向及三向開裂情況，由程序自行校核并進行剛度修正。巖體是否進入塑性狀態(tài)，按Druker-Prager準(zhǔn)則判別      （2）式中I1和J2分別為應(yīng)力張量的第一不變量和偏斜應(yīng)力張量的第二不變量，k是與巖體材料摩擦系數(shù)tg和凝聚力C有關(guān)的常數(shù)，由下式計算       （3）彈塑性矩陣為      

9、;  （4）        式中2.2軟弱結(jié)構(gòu)面非線性分析模型按層面法向抗拉材料分析，剪切滑移按Mohr-Coulomb條件校核         （5）    式中Cj和tgj分別為軟弱夾層抗剪強度參數(shù)。對于破碎帶寬度較大的斷層，按不抗拉彈塑性材料分析，是否進入塑性狀態(tài)的判別條件仍采用Drucker-prager準(zhǔn)則，只是材料摩擦系數(shù)和凝聚力改用斷層的相應(yīng)值，本構(gòu)矩陣仍沿用（4）式。 2

10、.3縱向圍堰穩(wěn)定安全系數(shù)計算方法在巖體穩(wěn)定性分析評價中，常常采用以下三種穩(wěn)定安全系數(shù)評價方法（1）點強度儲備安全系數(shù)（2）基礎(chǔ)整體抗滑安全系數(shù)：包括超載安全系數(shù)、強度儲備安全系數(shù)以及兼顧超載和強度儲備的綜合安全系數(shù)。（3）基礎(chǔ)抗滑安全系數(shù)的條分法，例如瑞典圓弧法。本次有限元分析以強度儲備法為穩(wěn)定安全性的主要評價方法，并輔以瑞典圓弧法。強度儲備法計算的基本方法是假定荷載不變，通過逐級下浮巖體強度參數(shù)，分析堰基變形破壞演變發(fā)展過程與超載倍數(shù)的關(guān)系，尋求堰基整體滑移時相應(yīng)的巖體強度參數(shù)下浮倍數(shù)Kc，即作為堰基整體抗滑穩(wěn)定的強度儲備安全系數(shù)。三  計算模型及計算參數(shù)3.1計算模型 

11、  本次二維計算主要針對六庫水電站縱向圍堰的四個典型剖面來進行。如圖1所示選定坐標(biāo)系統(tǒng)，X軸平行于縱向圍堰軸線，由下游指向上游；Y軸由左岸指向右岸；Z軸鉛直向上。圖2給出了剖面B的計算范圍，Z軸從720m取至地表，原點距縱向圍堰軸線水平距離111.749m，水平截取計算范圍170m。剖面B共計剖分節(jié)點541個，單元491個。沖積層、基巖、混凝土采用實體元模擬。在防滲墻與周圍基巖間布置了一層10cm厚夾層單元，以模擬剛度較大的防滲墻與相對軟弱的外周基巖間的相對錯動。 圖1 縱向圍堰有限元計算范圍 三維網(wǎng)格Y軸（橫河向）由左岸指向右岸，z軸 鉛直向上，計算范

12、圍與二維網(wǎng)格一致。X軸由下游指向上游。如圖1所示，共計切剖21個剖面，取順河向長220m。計算模型考慮了5道防滲隔板。共計剖分節(jié)點10431個，單元9126個。圖3繪出了有限元三維網(wǎng)格立體圖。 3.2巖體及結(jié)構(gòu)面物理力學(xué)參數(shù) 圖2  剖面B材料分區(qū)圖各類巖層物理力學(xué)參數(shù)計算選用值如表1所示。各類巖體抗拉強度Rt可由Mohr-Coulomb準(zhǔn)則導(dǎo)出   其中：C為凝聚力，為內(nèi)摩擦角    圖3縱向圍堰有限元三維網(wǎng)格立體圖        

13、四  研究方案 在本次六庫水電站縱向圍堰堰基穩(wěn)定分析中，其基本內(nèi)容包括以下幾個部分（1）根據(jù)質(zhì)量守恒及達西(Darcy)定律，對二維四個典型剖面分別進行汛期和枯期不同工況下的二維滲流場有限元計算以及在汛期工況下的三維滲流場有限元計算。（2）采用二維、三維線形有限元按彈性本構(gòu)模型，分別計算堰體自重、堰體自重汛期水頭、堰體自重枯期水頭工況下四個典型二維剖面和三維情況下的縱向碾壓混凝土圍堰變形及穩(wěn)定性。（3）采用二維非線形有限元按彈塑性本構(gòu)模型，分別計算和研究堰體自重汛期水頭、堰體自重枯期水頭工況下四個典型二維剖面分別在逐級下浮巖體強度參數(shù)Kc下的堰基變形破壞演變發(fā)展過程及整體穩(wěn)

14、定安全度。（4）采用三維非線形有限元按彈塑性本構(gòu)模型，研究堰體自重汛期水頭工況，巖體強度下降Kc=1.0，1.1，1.2，1.3，1.4，1.8六種方案下的堰基變形破壞演變發(fā)展過程及整體穩(wěn)定安全度。（5）分別計算四個典型二維剖面在自重、汛期、枯期水載作用下防滲墻及墊層所承受的軸力N，彎矩M和剪力Q，對防滲墻及墊層進行配筋設(shè)計。五、計算結(jié)果與分析5.1滲流分析   對比二維和三維滲流分析可知，平面及三維滲流場分析結(jié)果所得各剖面水頭分布和滲壓分布規(guī)律相近。均表現(xiàn)為：水頭在兩道防滲墻處迅速折減；而在砂卵沖積層和基巖則衰減緩慢。枯期滲壓均小于汛期滲壓。由滲壓分布還可發(fā)現(xiàn)，兩道防滲墻

15、之間的滲壓基本無折減（圖4）。 圖4 剖面B防滲墻汛期滲壓分布（MPa）5. 2變位分析堰體擋水后，在水推力和滲壓作用下，堰體和基礎(chǔ)發(fā)生變位。由于堰基巖變模很低（E00.040.048Gpa），致使堰體位移量值較大。以剖面B為例，汛期剖面B堰頂節(jié)點359向內(nèi)側(cè)變位-31.97cm，堰頂沉降-22.86cm。墊層底部左端節(jié)點213向內(nèi)側(cè)位移-16.86cm，沉降-29.68cm；墊層底部右端節(jié)點478向內(nèi)側(cè)位移-16.86cm，沉降-15.58cm(如圖5所示)。 圖5  剖面B汛期水平向位移uy(cm) 5.3應(yīng)力分析堰體及堰基大主壓應(yīng)力1和小主壓應(yīng)力3

16、隨高程降低而增大，在圍堰墊層與兩道防滲墻的接頭處附近存在應(yīng)力集中現(xiàn)象。剖面B汛期右側(cè)堰踵和左側(cè)堰趾拉應(yīng)力可達-0.35MPa；墊層中部拉應(yīng)力可達-0.4MPa。由于兩道防滲墻與圍巖相比剛度大，縱向圍堰的重量和荷載向防滲墻分載，墊層與兩道防滲墻的接頭處附近存在剪應(yīng)力集中現(xiàn)象。如圖6所示，剖面B汛期右側(cè)堰踵剪應(yīng)力約0.8MPa；左側(cè)堰趾剪應(yīng)力約1.3MPa。 圖6 剖面B汛期防滲墻剪應(yīng)力圖（MPa）         參 考 文 獻1   裴建良，張建海，陸民安等.百色RCC重力壩右岸擋水壩段滑移模式及滑移路徑研究J. 成都四川大學(xué)學(xué)報(工