(1.中國電建集團(tuán)北京勘測(cè)設(shè)計(jì)研究院有限公司，北京 100024 2.水利部建設(shè)管理與質(zhì)量安全中心，北京 100038)

傳統(tǒng)的邊坡問題依靠極限平衡法進(jìn)行穩(wěn)定性分析，但極限平衡法很難獲得邊坡的變形情況。有限元法不但滿足力的平衡條件，而且考慮了材料的應(yīng)力應(yīng)變關(guān)系，使得計(jì)算更加精確合理。有限元強(qiáng)度折減法是將巖土強(qiáng)度參數(shù)凝聚力c和摩擦角φ值除以某折減系數(shù)f得到一組新值代入有限元試算，邊坡失穩(wěn)時(shí)對(duì)應(yīng)的折減系數(shù)f為邊坡的穩(wěn)定安全系數(shù)。在有限元強(qiáng)度折減法中，主要對(duì)強(qiáng)度參數(shù)c、φ值進(jìn)行折減，但對(duì)于其他參數(shù)中的彈性模量E和泊松比μ的討論還不是很詳細(xì)。

本文利用有限元通用軟件ANSYS，基于Drucker-Prager屈服準(zhǔn)則(D-P準(zhǔn)則)，運(yùn)用有限元強(qiáng)度折減法對(duì)蘇洼龍水電站廠房后邊坡進(jìn)行邊坡穩(wěn)定計(jì)算，并重點(diǎn)分析了彈性模量和泊松比對(duì)邊坡應(yīng)力應(yīng)變及穩(wěn)定安全系數(shù)的影響。

1 基本原理

強(qiáng)度折減系數(shù)法的基本原理[1]：抗剪強(qiáng)度折減系數(shù)定義為在外荷載保持不變的條件下，邊坡內(nèi)巖土體所發(fā)揮的最大抗剪強(qiáng)度與外荷載在邊坡內(nèi)所產(chǎn)生的實(shí)際剪應(yīng)力之比。即保持巖土體的重力加速度為常數(shù)，將巖土體的強(qiáng)度參數(shù)c、φ值同時(shí)除以一個(gè)折減系數(shù)f，得到一組新的c′和φ′，然后作為新的材料參數(shù)進(jìn)行試算，當(dāng)邊坡處于臨界狀態(tài)時(shí)，即f再稍大一些，邊坡將發(fā)生破壞，此時(shí)對(duì)應(yīng)的f值被稱為邊坡的安全系數(shù)，其具體公式如下：

(1)

(2)

屈服條件是物體內(nèi)一點(diǎn)進(jìn)入屈服時(shí)，其應(yīng)力狀態(tài)所滿足的條件[2]。有限元分析方法所采用的本構(gòu)模型為理想彈塑性模型，巖土工程分析中常用的準(zhǔn)則有摩爾-庫倫準(zhǔn)則(M-C準(zhǔn)則)和D-P準(zhǔn)則，ANSYS軟件采用的是D-P準(zhǔn)則，即M-C不等角六邊形外接圓。D-P準(zhǔn)則可定義為：

(3)

式中，I1和J2是應(yīng)力張量的第一不變量和應(yīng)力偏張量的第二不變量，α、k是與巖土材料強(qiáng)度參數(shù)有關(guān)的常數(shù)。

針對(duì)邊坡有限元強(qiáng)度折減計(jì)算方法，目前常用的邊坡失穩(wěn)判據(jù)有三種[3]，?靜力平衡計(jì)算的不收斂性；?塑性區(qū)或者等效塑性應(yīng)變從坡腳到坡頂貫通狀態(tài)；?巖土體滑動(dòng)塊無限移動(dòng)，此時(shí)滑移面上位移發(fā)生突變且無限發(fā)展。本文在判斷邊坡失穩(wěn)時(shí)，考慮以上三種判據(jù)。

2 工程背景

2.1 工程概況

蘇洼龍水電站位于四川省甘孜藏族自治州和西藏自治區(qū)昌都藏族自治州的界河上，總庫容6.96億m3，多年平均流量929 m3/s，總裝機(jī)容量1200MW。水電站左岸引水發(fā)電系統(tǒng)包括引水系統(tǒng)、地面廠房系統(tǒng)和尾水建筑物。其中地面廠房系統(tǒng)位于一片緩坡地帶，為岸邊式地面廠房。

2.2 地質(zhì)條件

廠房工程區(qū)周圍均有覆蓋層發(fā)育，上下游側(cè)及坡頂部與冰水泥石流堆積的碎石土相接，坡下部與河床沖洪積砂卵礫石層相接，冰水泥石流堆積碎石土厚度約8～22m，上游側(cè)稍薄，頂部及下游側(cè)稍厚。各土石層對(duì)應(yīng)的地質(zhì)力學(xué)參數(shù)見表1。

表1 各土石層物理力學(xué)參數(shù)

選取廠房后邊坡某地質(zhì)剖面為分析對(duì)象，開挖邊坡高約40m，上一級(jí)邊坡高20m，坡比為1∶1，馬道寬3m，下一級(jí)邊坡高20m，坡比為1∶0.5，地質(zhì)剖面如圖1所示。

3 計(jì)算分析

3.1 計(jì)算模型

根據(jù)選用的廠房后邊坡地質(zhì)剖面，運(yùn)用有限元分析軟件，采用四節(jié)點(diǎn)等參單元，建立邊坡計(jì)算模型。開挖邊坡高40m，模型坡頂?shù)阶髠?cè)邊界為100m，坡腳到右側(cè)邊界為60m，左側(cè)豎向邊界高80m，考慮合理的網(wǎng)格劃分，開挖前單元總數(shù)12902，節(jié)點(diǎn)總數(shù)13036，開挖后單元總數(shù)12108，節(jié)點(diǎn)數(shù)10141，左右及底部邊界為鏈桿約束邊界。開挖后有限元網(wǎng)格模型如圖2所示。

圖1 邊坡地質(zhì)剖面

圖2 開挖后有限元網(wǎng)格模型

3.2 開挖邊坡穩(wěn)定計(jì)算

基于材料參數(shù)表1，對(duì)廠房后邊坡按平面應(yīng)變進(jìn)行應(yīng)力應(yīng)變分析計(jì)算，開挖后邊坡的合位移、最大主應(yīng)力、最小主應(yīng)力和等效塑性應(yīng)變?nèi)鐖D3～圖6所示。

圖3 開挖后邊坡合位移云圖(單位：m)

圖4 開挖后邊坡最大主應(yīng)力云圖 (單位：N)

圖5 開挖后邊坡最小主應(yīng)力云圖 (單位：N)

圖6 開挖后邊坡等效塑性應(yīng)變?cè)茍D

由計(jì)算結(jié)果可以看出，當(dāng)邊坡開挖完成以后，邊坡向開挖臨空區(qū)變形，上部的冰水堆積塊碎石土層區(qū)域變形最大，最大位移為0.067m；坡腳位置最大應(yīng)力為0.651MPa；開挖后坡頂后20m處出現(xiàn)局部塑性屈服區(qū)，馬道內(nèi)側(cè)7m處也出現(xiàn)小范圍塑性屈服區(qū)。

采用強(qiáng)度折減法，對(duì)開挖后的邊坡進(jìn)行應(yīng)力應(yīng)變計(jì)算，選取具有代表性的三個(gè)節(jié)點(diǎn)，統(tǒng)計(jì)其合位移變化規(guī)律，如圖7所示。塑性區(qū)貫通時(shí)的等效塑性應(yīng)變?nèi)鐖D8所示。有限元計(jì)算不收斂、節(jié)點(diǎn)特征位移突變和塑性區(qū)貫通作為失穩(wěn)判據(jù)時(shí)的各邊坡穩(wěn)定安全系數(shù)如表2所列。

圖7 節(jié)點(diǎn)合位移統(tǒng)計(jì)

圖8 塑性區(qū)貫通時(shí)的等效塑性應(yīng)變?cè)茍D (f=1.79)

表2 邊坡穩(wěn)定安全系數(shù)

由計(jì)算結(jié)果可以看出，有限元計(jì)算不收斂、特征位移突變和塑性區(qū)貫通作為判據(jù)時(shí)，各邊坡穩(wěn)定安全系數(shù)為1.86、1.85、1.79，大于規(guī)范要求值(A類Ⅰ級(jí)邊坡1.25～1.30)，故本計(jì)算邊坡剖面整體穩(wěn)定。

3.3 彈性模量和泊松比對(duì)邊坡計(jì)算的影響

為探討彈性模量E和泊松比μ對(duì)邊坡應(yīng)力應(yīng)變及穩(wěn)定安全系數(shù)的影響，針對(duì)冰積塊碎石土層的材料參數(shù)，基于單因素控制法則，考慮彈性模量和泊松比的合理取值范圍，首先進(jìn)行彈性模量分別為0.1E0和10E0時(shí)的邊坡穩(wěn)定計(jì)算，計(jì)算結(jié)果見表3和圖9～圖10，然后進(jìn)行泊松比為μ0-0.1和μ0+0.1時(shí)的邊坡穩(wěn)定計(jì)算，計(jì)算結(jié)果見表4和圖11～圖12。

圖9 等效塑性應(yīng)變?cè)茍D (E=5MPa, μ=3.50, f=1.80)

圖10 等效塑性應(yīng)變?cè)茍D (E=500MPa, μ=3.50,f=1.75)

表3 彈性模量對(duì)邊坡分析的影響

圖11 等效塑性應(yīng)變?cè)茍D (E=50MPa, μ=2.50, f=1.50)

圖12 等效塑性應(yīng)變?cè)茍D (E=50MPa, μ=4.50, f=1.79)

由表3可以看出，彈性模量相差10倍時(shí)對(duì)應(yīng)的坡腳最大應(yīng)力差值約為13%，因此彈性模量對(duì)應(yīng)力計(jì)算影響很?。挥捎趶椥阅Ａ勘揪褪亲冃文芰Φ闹笜?biāo)，故三種彈性模量對(duì)應(yīng)的最大位移相差較大，但彈性模量較大時(shí)，位移與應(yīng)力并不對(duì)應(yīng)成比例關(guān)系；此外，彈性模量的變化對(duì)邊坡穩(wěn)定安全系數(shù)影響很小，但彈性模量較大時(shí)，塑性區(qū)較易形成貫通，計(jì)算中彈性模量取值500MPa時(shí)，塑性區(qū)在1.75的折減系數(shù)下即達(dá)到貫通狀態(tài)，明顯小于其他相應(yīng)安全系數(shù)值。

表4 泊松比對(duì)邊坡分析的影響

由表4可以看出，泊松比相差0.1時(shí)，對(duì)應(yīng)的坡腳最大應(yīng)力差值約為14%和8%，因此彈性模量對(duì)應(yīng)力影響較??；三種泊松比對(duì)應(yīng)的最大位移相差亦不大，且位移值與泊松比成負(fù)相關(guān)；此外，泊松比較大時(shí)，泊松比的變化對(duì)邊坡穩(wěn)定安全系數(shù)影響很小，但是泊松比較小時(shí)，泊松比對(duì)邊坡穩(wěn)定安全系數(shù)影響較大，一方面，由圖11可以看出，塑性區(qū)達(dá)到貫通狀態(tài)時(shí)塑性區(qū)范圍較大，另一方面，相應(yīng)邊坡穩(wěn)定安全系數(shù)明顯較小，計(jì)算中泊松比取值0.25時(shí)，邊坡穩(wěn)定安全系數(shù)明顯小于其他兩組相應(yīng)系數(shù)值。

4 結(jié) 論

a.本計(jì)算邊坡整體穩(wěn)定，邊坡開挖后，上部冰水堆積塊碎石土層中存在小范圍的塑性屈服區(qū)，可采用部分邊坡支護(hù)措施或減緩邊坡坡度。

b.彈性模量的性質(zhì)決定了其對(duì)應(yīng)變的影響很大，但是對(duì)應(yīng)力和邊坡穩(wěn)定安全系數(shù)影響較小，且彈性模量較大時(shí)，塑性區(qū)較易形成貫通。

c.泊松比對(duì)應(yīng)力和應(yīng)變影響均較小。泊松比較大時(shí)，泊松比對(duì)邊坡穩(wěn)定安全系數(shù)影響很??；但泊松比較小時(shí)，相應(yīng)塑性區(qū)貫通時(shí)邊坡內(nèi)部的塑性區(qū)范圍較大，且邊坡穩(wěn)定安全系數(shù)較小。

d.本文計(jì)算分析了蘇洼龍水電站廠房后邊坡的應(yīng)力應(yīng)變及邊坡穩(wěn)定安全系數(shù)，并探討了彈性模量和泊松比對(duì)邊坡應(yīng)力應(yīng)變及穩(wěn)定安全系數(shù)的影響，得到了一些定性的結(jié)論，還需要進(jìn)一步探索本構(gòu)模型中各因素之間的相互影響，從理論機(jī)制上尋找更多的依據(jù)。