萬(wàn)海龍，齊一帆，張 萍，呂從聰，熊勃勃

（1.中國(guó)電建集團(tuán)華東勘測(cè)設(shè)計(jì)研究院有限公司,浙江 杭州 311122；2.三峽大學(xué)水利與環(huán)境學(xué)院,湖北 宜昌443002）

1 引言

降雨入滲是誘發(fā)邊坡失穩(wěn)的主要原因之一,是一種典型的非飽和流固耦合現(xiàn)象,且研究表明絕大多數(shù)的滑坡均發(fā)生在雨季［1-3］。究其原因在于,降雨入滲造成坡體內(nèi)含水量增加,導(dǎo)致非飽和區(qū)域的基質(zhì)吸力減小,土顆粒間的有效應(yīng)力降低［4］。由于降雨入滲對(duì)邊坡穩(wěn)定性的影響與降雨強(qiáng)度、歷時(shí)、雨量、雨型以及邊坡的地形地貌等諸多因素有關(guān)［5］,因而難以全面考慮降雨對(duì)邊坡穩(wěn)定的影響因素［2］。

目前,分析降雨條件下邊坡穩(wěn)定的方法主要包括：簡(jiǎn)化計(jì)算滲流的極限平衡法、有限元方法和極限分析法［6］。國(guó)內(nèi)外學(xué)者［7-9］已對(duì)降雨入滲條件下邊坡的非飽和滲流過(guò)程及其對(duì)邊坡穩(wěn)定性的影響進(jìn)行了較多研究,然而這些研究多是針對(duì)加固前邊坡,對(duì)錨桿加固后降雨邊坡的穩(wěn)定性研究涉及較少。年延凱等［10］利用考慮樁-土-邊坡相互作用的強(qiáng)度折減有限元程序,探討了抗滑樁-邊坡體系的計(jì)算模型尺度,并探討了抗滑樁參數(shù)對(duì)邊坡穩(wěn)定安全系數(shù)及臨界滑動(dòng)面的影響；唐曉松等［11］結(jié)合有限元強(qiáng)度折減法分析了抗滑樁間距對(duì)邊坡加固效果的影響；他們雖然研究了抗滑樁加固對(duì)邊坡穩(wěn)定性的影響,但卻未考慮降雨因素。CAI等［12］利用有限元強(qiáng)度折減法,研究了降雨條件下水平排水對(duì)抗滑樁加固邊坡穩(wěn)定性的影響。李寧等［13］提出了降雨作用下抗滑樁邊坡穩(wěn)定性分析的有限元強(qiáng)度折減法,并應(yīng)用該方法研究了降雨前后抗滑樁對(duì)邊坡穩(wěn)定性的影響；郭震山等［14］研究了降雨入滲對(duì)邊坡孔隙水壓力分布及邊坡穩(wěn)定性的影響,并對(duì)降雨強(qiáng)度、基質(zhì)吸力及樁頂約束形式對(duì)邊坡穩(wěn)定性的影響進(jìn)行了參數(shù)分析；次仁拉姆［15］基于有限元分析,對(duì)不同降雨工況下的孔隙水壓分布、錨索軸力及采用強(qiáng)度折減法計(jì)算的邊坡安全系數(shù)變化規(guī)律進(jìn)行了研究。雖然這些學(xué)者對(duì)降雨條件下的抗滑樁邊坡開(kāi)展了相應(yīng)研究,但未涉及錨桿形式、布置方式及其數(shù)量對(duì)加固后邊坡穩(wěn)定性影響的研究。

本文以會(huì)東縣松坪山塘邊坡為例,基于ABAQUS有限元軟件進(jìn)行降雨條件下滲流-應(yīng)力耦合分析,并采取有限元強(qiáng)度折減法求解其安全系數(shù),分析降雨入滲對(duì)山塘邊坡穩(wěn)定性的影響,并在此基礎(chǔ)上進(jìn)一步探討砂漿錨桿的長(zhǎng)度、位置、數(shù)量和縱向間距等因素對(duì)降雨邊坡穩(wěn)定性的影響。

2 研究理論基礎(chǔ)

2.1 巖土的非飽和滲流理論

實(shí)現(xiàn)非飽和滲流計(jì)算的關(guān)鍵在于確定非飽和區(qū)滲透系數(shù)和飽和度隨基質(zhì)吸力的變化關(guān)系,文中采用Cho和Lee經(jīng)驗(yàn)公式［7］,即非飽和區(qū)滲透系數(shù)隨基質(zhì)吸力的變化關(guān)系為：

式中：Kw為非飽和區(qū)材料的滲透系數(shù)；Kws為土體飽和時(shí)的滲透系數(shù)；ua為孔隙氣壓力；uw為孔隙水壓力；aw、bw和cw為材料常數(shù)。

飽和度隨基質(zhì)吸力的變化關(guān)系為：

式中：Sr為飽和度；Si為殘余飽和度；Sn為最大飽和度；as、bs、cs均為材料常數(shù)。

由于沒(méi)有上述參數(shù)的試驗(yàn)數(shù)據(jù),所以在后續(xù)數(shù)值分析過(guò)程中,上述參數(shù)均參考文獻(xiàn)［2］和［7］進(jìn)行取值。

2.2 非飽和土邊坡的強(qiáng)度折減法

由于極限平衡法和傳統(tǒng)強(qiáng)度折減法均不能考慮孔隙水壓力和基質(zhì)吸力的影響,因此難以直接用來(lái)求解考慮降雨和地下水位作用下的邊坡穩(wěn)定性問(wèn)題。Fredlund在飽和土的Mohr-Coulomb強(qiáng)度公式基礎(chǔ)上考慮基質(zhì)吸力對(duì)土體強(qiáng)度的影響,并提出了經(jīng)典的非飽和土抗剪強(qiáng)度計(jì)算公式：

式中：c＇為有效粘聚力；σ為正應(yīng)力；ua為孔隙氣壓力；uw為孔隙水壓力；ua-uw為基質(zhì)吸力；φ＇為有效內(nèi)摩擦角；tanφb為抗剪強(qiáng)度隨基質(zhì)吸力增加的速率。

從式（3）中可看出,基質(zhì)吸力對(duì)非飽和土抗剪強(qiáng)度有一定影響,通常將基質(zhì)吸力視為黏聚力的一部分,則粘聚力變?yōu)椋?/p>

因此,基于Fredlund理論的非飽和土強(qiáng)度折減法中粘聚力和內(nèi)摩擦角的折減公式分別變?yōu)椋?/p>

計(jì)算中假設(shè)不同的強(qiáng)度折減系數(shù),根據(jù)折減后的強(qiáng)度參數(shù)進(jìn)行有限元分析,并以坡頂點(diǎn)位移拐點(diǎn)作為邊坡失穩(wěn)的判據(jù)［16］。

3 降雨條件下松坪山塘邊坡穩(wěn)定性分析

3.1 工程概況

松坪山塘位于四川省涼山彝族自治州南端會(huì)東縣松坪鄉(xiāng)集鎮(zhèn)內(nèi),該地地理環(huán)境復(fù)雜多樣,垂直差異明顯。境內(nèi)地形呈中間低、四周高的山間盆地狀,海拔高度在640 m~3331.8 m之間。場(chǎng)地原始地面高程為1840 m~1880 m,開(kāi)挖后山塘底部高程為1947.8 m。經(jīng)地質(zhì)勘測(cè)知邊坡內(nèi)地層巖性不一,由上往下分別是含礫質(zhì)黏土、粉土和千枚巖,地質(zhì)情況較為復(fù)雜,且氣候變化大,降雨頻發(fā),在自然降雨條件下可能會(huì)引起坡內(nèi)土體的輕微滑動(dòng)。

選取典型剖面建立二維有限元模型,見(jiàn)圖1。模型的底部采用固定約束,兩側(cè)采用鉸支座約束,邊坡的頂面和坡面（水位以上）設(shè)置為降雨入滲邊界。有限元模型共劃分單元2503個(gè),節(jié)點(diǎn)7722個(gè)。各層材料采用的物理力學(xué)參數(shù)見(jiàn)表1［17-18］。

圖1 典型剖面二維有限元模型

表1 材料的物理力學(xué)指標(biāo)

3.2 考慮降雨對(duì)山塘邊坡內(nèi)滲流場(chǎng)的影響

得到不考慮降雨及降雨72 h之后兩種工況下邊坡內(nèi)孔壓等值線云圖,見(jiàn)圖2。

圖2 不同工況下邊坡內(nèi)滲流場(chǎng)分布

對(duì)比分析可以發(fā)現(xiàn),在降雨72 h之后,山塘邊坡的孔隙水壓力分布仍呈層狀分布,但邊坡表面的孔隙水壓力有所提高。究其原因在于,在邊坡的降雨入滲過(guò)程中,邊坡表層逐漸變濕,基質(zhì)吸力隨著降低,邊坡表層的孔隙水壓力在不斷增加并且形成一個(gè)暫態(tài)飽和區(qū),由初始的負(fù)孔隙水壓力逐漸上升為正孔隙水壓力。

3.3 考慮降雨對(duì)山塘邊坡穩(wěn)定的影響

降雨72 h之后邊坡的水平位移分布及折減系數(shù)隨山塘邊坡水平位移U1的變化關(guān)系見(jiàn)圖3~圖4。

圖3 邊坡失穩(wěn)時(shí)的水平位移分布

圖4 邊坡的安全系數(shù)

通過(guò)強(qiáng)度折減法的評(píng)價(jià)標(biāo)準(zhǔn)及特征部位（坡頂）的位移拐點(diǎn)可以判斷出,不考慮降雨時(shí)邊坡的安全系數(shù)為1.348,大于《建筑邊坡工程技術(shù)規(guī)范》規(guī)定安全系數(shù)允許值1.30,說(shuō)明降雨前邊坡處于穩(wěn)定狀態(tài)；而降雨72 h之后邊坡的安全系數(shù)降為1.194,小于規(guī)范建議的允許值。綜上可知,研究坡段受降雨影響劇烈,在降雨72 h之后不穩(wěn)定,并產(chǎn)生較大滑動(dòng),最大水平位移達(dá)11.03 cm。

4 邊坡加固措施及錨桿參數(shù)對(duì)降雨邊坡穩(wěn)定性的影響分析

4.1 邊坡加固措施

由邊坡加固前穩(wěn)定性分析結(jié)果可知,研究坡段在降雨工況下處于不穩(wěn)定狀態(tài),需要進(jìn)行邊坡加固。結(jié)合項(xiàng)目區(qū)地質(zhì)災(zāi)害發(fā)育特點(diǎn)、相似工程的治理經(jīng)驗(yàn),主要采用砂漿錨桿加固的方式,并輔以框格梁并碼放植生袋進(jìn)行坡面綠化,同時(shí)在邊坡周邊設(shè)置排水溝以減少降雨對(duì)邊坡的不利影響。沿坡面布置砂漿錨桿,錨孔直徑取為110 mm,角度下傾34°。

為進(jìn)一步探究錨桿參數(shù)對(duì)降雨邊坡穩(wěn)定性的影響,分別針對(duì)不同錨桿長(zhǎng)度、布置位置、以及錨桿數(shù)量開(kāi)展數(shù)值研究。

4.2 錨桿長(zhǎng)度對(duì)結(jié)果的影響

降雨72 h之后邊坡塑性區(qū)分布見(jiàn)圖5,從圖中可以看出,邊坡塑性區(qū)完全貫通,邊坡失穩(wěn),此時(shí)預(yù)設(shè)值錨桿處的塑性區(qū)深度為8.3 m。為探究降雨條件下錨桿長(zhǎng)度不同對(duì)邊坡安全系數(shù)的影響,在邊坡中間位置分別布置桿長(zhǎng)為6 m、9 m、12 m和15 m的錨桿,并與不設(shè)置錨桿情況進(jìn)行比較。

圖5 降雨邊坡塑性區(qū)分布

計(jì)算得到不同錨桿長(zhǎng)度下邊坡的安全系數(shù)隨位移變化分布見(jiàn)圖6。從圖中可以看出,隨著錨桿長(zhǎng)度的不斷增加,降雨邊坡的安全系數(shù)逐漸增大,說(shuō)明在一定范圍內(nèi),樁長(zhǎng)越長(zhǎng),錨桿的抗滑作用越明顯,邊坡安全系數(shù)越高。如果設(shè)置錨桿長(zhǎng)度小于塑性區(qū)深度8.3 m,比如設(shè)置6 m錨桿,則對(duì)提高邊坡穩(wěn)定性效果不理想,這是因?yàn)殄^桿長(zhǎng)度較短時(shí),滑動(dòng)面會(huì)越過(guò)錨桿底部,造成錨桿和滑動(dòng)土體一起下滑。

圖6 錨桿長(zhǎng)度對(duì)安全系數(shù)的影響

4.3 錨桿位置對(duì)結(jié)果的影響

分別將15 m的錨桿布置在圖7 中的邊坡中點(diǎn)I處、中間偏上3 m的II處和中間偏上6 m的III處,并與不布置錨桿情況進(jìn)行比較。計(jì)算得到錨桿位置不同對(duì)邊坡安全系數(shù)的影響見(jiàn)圖8。

圖7 降雨邊坡錨桿位置

圖8 錨桿位置對(duì)安全系數(shù)的影響

從圖中可以看出,錨桿分別設(shè)置在I、II、III處時(shí)邊坡的安全系數(shù)依次為1.258、1.244和1.236,說(shuō)明考慮降雨條件下,錨桿位置距離坡腳越遠(yuǎn),錨桿加固邊坡的安全系數(shù)越小,加固效果越不理想。分析其原因可能在于降雨引起坡腳處滲透力增大,進(jìn)而導(dǎo)致坡腳處率先發(fā)生破壞,錨桿距離坡腳較遠(yuǎn)起不到穩(wěn)固坡腳的作用。

4.4 錨桿數(shù)量對(duì)結(jié)果的影響

上述分析表明,無(wú)論是增加錨桿長(zhǎng)度或合理調(diào)整錨桿的位置均有利于提高邊坡的穩(wěn)定,但計(jì)算得到的安全系數(shù)均小于規(guī)范建議值1.30,仍不能滿足要求。為此,文中進(jìn)一步探討錨桿數(shù)量對(duì)邊坡穩(wěn)定性的影響,沿邊坡分別布置1根、2根和3根錨桿,并與不布置錨桿情況進(jìn)行比較。計(jì)算得到布置多根錨桿對(duì)降雨邊坡安全系數(shù)的影響,見(jiàn)圖9。

圖9 多根錨桿對(duì)安全系數(shù)的影響

從圖9（a）中可以看出,隨著錨桿數(shù)量的增加,降雨邊坡的安全系數(shù)逐漸增大,當(dāng)布置3根錨桿時(shí),安全系數(shù)提高至1.341,大于1.30,說(shuō)明布置3根及以上錨桿后邊坡處于穩(wěn)定狀態(tài)。從圖9（b）中可以看出,考慮降雨情況下,錨桿布置的過(guò)密,不利于雨水的排放。分析原因可能在于,樁后地下水位的抬升將增加滑坡體的飽和狀態(tài),使得抗滑力減小,進(jìn)而降低錨桿抗滑效果。

5 結(jié)論

以松坪山塘邊坡為研究對(duì)象,基于巖土非飽和滲流理論和非飽和土的強(qiáng)度折減法,并借助ABAQUS軟件研究了降雨對(duì)山塘邊坡穩(wěn)定性的影響,并在此基礎(chǔ)上進(jìn)一步對(duì)降雨條件下錨桿加固邊坡進(jìn)行研究,得到如下結(jié)論：

（1）與不考慮降雨情況相比,考慮降雨之后,山塘邊坡的孔隙水壓力分布仍呈層狀分布,但邊坡表面的孔隙水壓力有所提高；邊坡受降雨影響強(qiáng)烈,在降雨之后易邊坡失穩(wěn)。因此,需要對(duì)降雨邊坡進(jìn)行加固處理。

（2）在一定范圍內(nèi),錨桿長(zhǎng)度越長(zhǎng),降雨邊坡的安全系數(shù)越大,錨桿的抗滑作用越明顯。如果設(shè)置錨桿長(zhǎng)度過(guò)短,且小于塑性區(qū)最大深度,則對(duì)提高邊坡穩(wěn)定性效果不理想,所以在進(jìn)行錨桿加固時(shí),應(yīng)根據(jù)規(guī)范和實(shí)際情況適當(dāng)加大錨桿長(zhǎng)度。

（3）對(duì)于常降雨位置邊坡,可布置多根錨桿以提高邊坡的安全系數(shù)；同時(shí),應(yīng)合理布置多錨桿之間的間距,以便更好地發(fā)揮錨固效果。