許海嘯,王延壽
(甘肅省有色金屬地質(zhì)勘查局天水礦產(chǎn)勘查院,甘肅 天水 741020)
中國黃土分布面積之廣及其深度之深,堪稱世界之最,黃土分布面積約為6.4×105km2[1],最大深度約為400 m[2]。黃土邊坡分布極為廣泛,加之大量的公路、鐵路、建筑、機場工程等建在黃土地區(qū)[3-7],不可避免地開挖邊坡坡腳,進而涉及黃土邊坡穩(wěn)定性問題更為廣泛。另外,由于灌溉作用也引起了大量的黃土邊坡失穩(wěn)問題,如涇陽南塬[8-9]和蘭州黑方臺[10]地區(qū)經(jīng)常發(fā)生灌溉誘發(fā)型黃土滑坡。
黃土是分布在干旱、半干旱地區(qū)的一種多孔隙、結(jié)構(gòu)疏松的第四紀沉積物,其骨架顆粒接觸關(guān)系有架空、架空-鑲嵌和分散接觸3種,骨架顆粒連接形式有接觸、膠結(jié)和半膠結(jié)連接3種,形成的孔隙主要以架空和鑲嵌孔隙為主。這種具有特殊微結(jié)構(gòu)的土體導(dǎo)致其具有特殊的水理性質(zhì)和力學(xué)性質(zhì),如收縮性大、滲透性好、崩解性強、壓縮性高和抗剪性低[2]等。
由此可見,黃土邊坡很容易發(fā)生變形失穩(wěn),尤其是在地震、降水、開挖坡腳以及多種作用下極易發(fā)生滑坡[11-15]。因黃土邊坡失穩(wěn)而造成人民財產(chǎn)損失和威脅人身安全的事故常有發(fā)生,故對黃土邊坡穩(wěn)定性進行正確分析尤為重要??捎糜谶吰路€(wěn)定性分析的方法有多種[16],研究主要討論采用基于極限平衡條分法的瑞典法、簡化畢肖普法和簡化簡布法分析黃土邊坡穩(wěn)定性時的適用性,即以條分法的土條寬度、搜索時的圓心步長和半徑步長為變量時計算得到的穩(wěn)定性系數(shù)的可靠性。
極限平衡法是當(dāng)前國內(nèi)外應(yīng)用最廣泛的邊坡穩(wěn)定性分析方法。該方法是在已知滑移面上對邊坡進行靜力平衡計算,從而求出邊坡穩(wěn)定性系數(shù),即必須事先知道滑移面的位置和形狀。當(dāng)滑移面為一簡單平面時,靜力平衡可采用解析法計算,獲得解析解。然而黃土邊坡滑移面為一圓弧,無法獲得解析解,通常采用條分法求解,此時,坡體為一靜不定問題,通過對某些未知量作假定,使方程的數(shù)目與未知數(shù)數(shù)目相等,從而使問題成為靜定,求出邊坡穩(wěn)定性系數(shù),如瑞典法、簡化畢肖普法和簡化簡布法等。
由于黃土邊坡穩(wěn)定性分析問題是一個高次超靜定問題,因而各種條分法都需要對未知量的數(shù)值或分布作一些假定,以便求解未知數(shù)。
瑞典法在求條底反力時忽略了條間力的作用,相當(dāng)于假定條間力Pi和Pi-1的合力ΔPi的方向平行于條底,即
Θi=αi,
(1)
其中:Θi為條間力Pi和Pi-1的合力ΔPi的方向與水平方向的傾角;αi為條底傾角。
簡化畢肖普法和簡化簡布法假定條間力的方向為水平方向,相當(dāng)于假定條間剪力差為0或條間力Pi和Pi-1的合力ΔPi的方向為水平方向,即
Θi=0。
(2)
(1) 瑞典法 圓弧滑動面的靜力方程為
(3)
條底剪力Si方程為
(4)
將方程(4)代入方程(3),得穩(wěn)定性系數(shù)方程:
(5)
其中:Fs為滑動圓弧的穩(wěn)定性系數(shù);τf為滑裂面上的平均抗剪強度;li為條底i的長度;c′i、tanφ′i為條底i的有效抗剪強度;Wi為體力;Ks為地震影響系數(shù);Qi為坡面外力;Pi為條間力;Si為條底剪力;Ni為條底法向力。
(2) 簡化畢肖普法 簡化畢肖普法亦是假定滑裂面為圓柱面(剖面圖上是圓弧),同時還假定條間力的方向為水平方向,得到的穩(wěn)定性系數(shù)方程與瑞典法的相同。但是當(dāng)Θi= 0時,條底法向力Ni可簡化為
(Wi+Qicosθi-c′milisinαi+
Uitanφ′misinαi)。
(6)
條底法向力Ni方程與穩(wěn)定性系數(shù)Fs有關(guān),因此,簡化畢肖普法的穩(wěn)定性系數(shù)方程是一個隱式表達式,求解需要迭代。令
(7)
則方程(5)可以寫為
F=Φ(F)。
(8)
方程(8)是一個便于利用不動點迭代法進行迭代的形式。不動點迭代法的迭代公式為
F(k+1)=Φ(F(k))。
(9)
如果得到的序列{F(k)}滿足
(10)
則F(*)就是Φ的不動點。由此可求出簡化畢肖普法的穩(wěn)定性系數(shù)Fs。
(3) 簡化簡布法 按靜力平衡條件,∑Fz=0,得
Nicosθi=Wi+ΔHi-Tisinθi,
(11)
其中:Ti為條底切向力;ΔHi為切向力Hi、Hi-1的合力。
∑Fx=0,得
ΔPi=Ticosθi-Nisinθi。
(12)
將式(11)代入式(12)后整理得
(13)
根據(jù)極限平衡條件,考慮穩(wěn)定系數(shù)Fs,有
(14)
將式(11)代入式(14),得
(15)
將式(15)代入式(13),得
ΔHi)tanφi]-(Wi+ΔHi)tanθi。
(16)
ΔHi)tanφi]-∑(Wi+ΔHi)tanθi=0,
(17)
整理后得
(18)
某黃土邊坡剖面如圖1所示,坡體由第四系上更新統(tǒng)馬蘭黃土構(gòu)成。坡高為6 m,坡度62°,坡長為6.9 m,天然容重r=18.5 kN/m3,飽和容重rsr=19 kN/m3,粘聚力c=20 kPa,內(nèi)摩擦角φ=24°。
圖1 算例計算模型Fig.1 Example computational model
對上述邊坡模型分別采用瑞典法、簡化畢肖普法和簡化簡布法進行穩(wěn)定性系數(shù)計算,并且對每一種方法在土條寬度、圓心步長和半徑步長分別等于0.1,0.2,…,1.0時計算穩(wěn)定性系數(shù),進行對比分析。
穩(wěn)定性系數(shù)計算結(jié)果見表1。由表1可知,瑞典法計算得到的穩(wěn)定性系數(shù)隨土條寬度、圓心步長和半徑步長的增大而增大,但是變化較小,最大值與最小值相差0.023;簡化畢肖普法計算得到的穩(wěn)定性系數(shù)亦是隨土條寬度、圓心步長和半徑步長的增大而增大,在土條寬度、圓心步長和半徑步長等于0.1時,穩(wěn)定性系數(shù)最小,為1.356,等于0.3時突然變大,為1.571,之后變化相對平穩(wěn);簡化簡布法計算得到的穩(wěn)定性系數(shù)隨土條寬度、圓心步長和半徑步長的增大先增大后變小,在土條寬度、圓心步長和半徑步長等于0.1時,穩(wěn)定性系數(shù)最小,為1.475,等于0.3時突然變大,為1.608,之后又變小,但是變化相對平穩(wěn)。3種方法求得的穩(wěn)定性系數(shù)在土條寬度、圓心步長和半徑步長等于0.8時基本一致(見圖2),分別為1.601、1.600、1.599(見表1)。
表1 穩(wěn)定性系數(shù)
為驗證3種方法計算結(jié)果的準確性,利用邁達斯軟件中巖土和隧道分析系統(tǒng)(GTS,geotechnical and tunnel analysis system)進行對比,計算模型如圖3所示,節(jié)點共3 647個,單元共3 249個。
圖2 穩(wěn)定性系數(shù)變化Fig.2 Stability coefficient change graph
由基于有限元法的邁達斯GTS數(shù)值模擬軟件(強度折減法)求得的穩(wěn)定性系數(shù)為1.825。根據(jù)最大剪應(yīng)變云圖(見圖4)和有效塑性應(yīng)變云圖(見圖5)可知,最大剪應(yīng)變和有效塑性應(yīng)變只發(fā)生在坡腳處,并未形成貫通區(qū)域,故該邊坡尚未發(fā)生變形破壞。值得一提的是,強度折減法考慮了土體的非線性本構(gòu)關(guān)系,考慮了應(yīng)力與應(yīng)變的關(guān)系[17],故由其計算的結(jié)果稍大一些,而極限平衡條分法計算的結(jié)果略偏于安全。經(jīng)綜合判斷,該黃土邊坡的穩(wěn)定性系數(shù)為1.600,處于穩(wěn)定狀態(tài)。
對比分析結(jié)果表明,當(dāng)采用極限平衡條分法進行分析黃土邊坡穩(wěn)定性時,簡化畢肖普法和簡化簡布法在土條寬度、圓心步長和半徑步長≥0.3時,計算結(jié)果較為可靠,瑞典法受其影響較小,即當(dāng)簡化畢肖普法和簡化簡布法計算結(jié)果小于瑞典法時,由瑞典法計算的穩(wěn)定性系數(shù)較為可靠;當(dāng)3種方法計算的結(jié)果較接近時,應(yīng)進一步進行計算,直至基本一致時,此值為該邊坡的穩(wěn)定性系數(shù)。
圖3 計算模型和網(wǎng)格劃分Fig.3 Computational model and meshing
圖4 最大剪應(yīng)變云圖Fig.4 Maximum shearing strain cloud map
圖5 有效塑性應(yīng)變云圖Fig.5 Effective plasticity strain cloud map
(1) 瑞典法計算的穩(wěn)定性系數(shù)隨土條寬度、圓心步長和半徑步長的增大而增大,但是受其變化影響不大,即當(dāng)分析黃土邊坡穩(wěn)定性時,瑞典法幾乎不受土條寬度、圓心步長和半徑步長變化的影響,其計算結(jié)果可直接用之;
(2) 簡化畢肖普法和簡化簡布法計算的穩(wěn)定性系數(shù)隨土條寬度、圓心步長和半徑步長的變化而變化,當(dāng)取值<0.3時,計算結(jié)果不可用,取值≥0.3時,計算結(jié)果較為可靠;
(3) 當(dāng)簡化畢肖普法和簡化簡布法計算的穩(wěn)定性系數(shù)小于瑞典法者時,應(yīng)以瑞典法計算結(jié)果為準,當(dāng)3種方法計算結(jié)果基本一致時,此時的數(shù)據(jù)宜作為該黃土邊坡的穩(wěn)定性系數(shù)。