■武金輝（蘭州大學土木工程與力學學院甘肅蘭州730000）

蘭州市某人工填土高邊坡穩(wěn)定性分析

■武金輝
（蘭州大學土木工程與力學學院甘肅蘭州730000）

人工填土在外力作用下使得土體原結(jié)構(gòu)遭受破壞，強度降低，工程地質(zhì)條件差。以蘭州市某填土不穩(wěn)定斜坡為例，根據(jù)邊坡結(jié)構(gòu)特點、變形特征和室內(nèi)土工實驗結(jié)果選取合理的巖土計算參數(shù)，利用FlAC3D有限差分法建立準三維模型模擬原始坡體變形特征并求解安全系數(shù)與傳統(tǒng)的極限平衡法進行比較。結(jié)果表明：兩種方法求解的結(jié)果發(fā)現(xiàn)坡體發(fā)生破壞時的危險滑動面的位置極為吻合，所求解的安全系數(shù)也非常相近，為坡體的穩(wěn)定性分析和后期的工程治理提供了很好的科學依據(jù)。

工程地質(zhì)人工填土邊坡穩(wěn)定性數(shù)值模擬

甘肅省為全國地質(zhì)災害多發(fā)區(qū)，蘭州市作為省會城市，近年來地質(zhì)災害頻頻發(fā)生，特別是人類工程活動對地質(zhì)災害的引發(fā)和加劇作用越來越明顯，經(jīng)統(tǒng)計，1949年以來，蘭州市共發(fā)生地質(zhì)災害238次，其中特大型9次，已造成773人死亡，累計直接經(jīng)濟損失達7.56億元，省、市政府高度重視地質(zhì)災害的防治工作。

文中以蘭州市九州開發(fā)區(qū)一典型填方邊坡為研究對象，采用極限平衡法和有限差差分法對比研究，定量評價其穩(wěn)定性及計算了一定安全系數(shù)下的下滑力為邊坡的治理提供了可靠地依據(jù)同時也為同類邊坡的勘察和設計提供了參考依據(jù)。

1　填土高邊坡概況

1.1邊坡基本特征

該邊坡位于蘭州市九州開發(fā)區(qū)春暉路，坡頂為訓練場地跑道及后勤中心，下部為雙信實業(yè)場地。斜坡坡體上部為人工填土，填土成分為以粉土為主并夾有砂巖、砂質(zhì)泥巖塊體，下部出露古近系泥巖、砂巖、泥質(zhì)砂巖，斜坡坡體呈弧形延伸，坡頂較為平坦起伏度小，坡向79°-130°，坡度45°—60°，坡面較順直，延伸長120m，坡體主要為人工填土，填土厚度大，填土高度最大達到35.5m，橫向厚度最大達40m，其主要軟弱結(jié)構(gòu)面為近水平基巖與人工填土，勘探結(jié)果未見地下水。

1.2現(xiàn)今變形跡象

現(xiàn)場調(diào)查發(fā)現(xiàn)坡腳（填土與基巖接觸面）處存在不連續(xù)土體隆起及土體的剝墜落現(xiàn)象，坡頂變形特征主要表現(xiàn)為坡頂?shù)孛娈a(chǎn)生不均勻沉降引發(fā)錯落現(xiàn)象，主要分布于食堂東側(cè)及東北側(cè)。其中食堂東側(cè)不均勻沉降緊挨食堂東側(cè)墻壁，產(chǎn)生的錯落錯距一般在2-3cm，最大可達5-6cm，造成食堂東側(cè)混凝土硬化地面沿接縫發(fā)生錯斷現(xiàn)象，近南北向延伸；食堂東北側(cè)不均勻沉降產(chǎn)生的錯落現(xiàn)象較嚴重，錯距一般在50-60cm左右，形成多級小型錯斷平臺，坡體頂部前緣可見掉塊現(xiàn)象。坡頂上部同時還發(fā)育有多條裂縫，裂縫張度最大達0.1m，最長延伸長度10—15m，距離坡體邊緣約1.5—2.5m處，部分已填埋，坡頂邊緣到場地8m范圍內(nèi)出現(xiàn)約8—10cm的下沉，致使威脅特警餐廳散水管開裂下沉，并有發(fā)展的趨勢。

1.3滑坡變形影響因素

其發(fā)生原因主要有三個方面：

（1）是由于邊坡為人工填土層，厚度大，而人工填土本身受到擾動，原狀結(jié)構(gòu)遭受破壞，力學性質(zhì)差，強度低，易在外力作用下發(fā)生失穩(wěn)破壞。

（2）是坡腳過度開挖后，坡度加大，坡體的水平支撐力降低或減小，坡體不同部位處于不同的應力釋放狀態(tài)：坡頂處于水平拉伸狀態(tài)，中部處于水平卸載狀態(tài)而處于拉伸和剪切狀態(tài)，坡底開挖部分則處于垂直向上的應力卸載狀態(tài)，臨空面部位應力相對集中，變形破壞為自身應力釋放。

（3）三是坡體變形破壞受自身軟弱結(jié)構(gòu)面的控制，坡體頂部覆蓋的人工素填土也是發(fā)生災害的主要位置。

2　基于FLAC3D有限差分法的數(shù)值模擬

2.1FLAC3D簡介

FLAC3D采用的顯式拉格朗日算法和混合離散分區(qū)技術(shù)，能夠非常準確地模擬材料的塑性破壞和流動，由于無需形成剛度矩陣，因此基于較小內(nèi)存空間就能夠求解較大范圍的三維問題。它適用于絕大多數(shù)的工程力學問題，尤其適用于材料的彈塑性、大變形分析和施工過程中巖土工程的數(shù)值模擬，因而在國內(nèi)外得到廣泛認可和應用。

2.2模型的建立和約束條件

為了消除邊界條件的影響，將上述邊坡按原始邊坡形態(tài)建立模型后按照張魯渝、鄭穎人等人（2003）的建議，即右邊界至坡腳等于1.5H（邊坡高度），坡頂至左邊界的距離等于2.5H，坡頂?shù)降撞窟吔绲扔?H的方式進行擴大，擴大后的模型如下，坡體高64m，底座長147m，向滑坡后部延伸80m；為了提高計算結(jié)果的可靠性適當加密坡體網(wǎng)格的密度，特別是填土層網(wǎng)格劃分較密，計算模型在ansys里通過指定單元體邊長來劃分網(wǎng)格，每個模型均按相同的單元體邊長劃分。

3　傳統(tǒng)極限平衡法驗算

為了增加坡體分析的可靠性，利用GEO5巖土計算軟件19版軟件內(nèi)的土質(zhì)邊坡穩(wěn)定性分析模塊對上述坡體的穩(wěn)定性進行驗算。根據(jù)《滑坡防治工程勘察規(guī)范》（DZ/T0218—2006）：穩(wěn)定系數(shù)：Fs<1. 00為不穩(wěn)定，1.00≤Fs<1.05欠穩(wěn)定,1.05≤Fs<1.15基本穩(wěn)定，F(xiàn)s≥1.15穩(wěn)定，將軟件安全系數(shù)設置為SFs=1.15并選擇驗算方法為畢曉普法（Bishop），采用圓弧滑動面法進行邊坡穩(wěn)定性分析，并利用GEO5中的自動搜索圓弧面的功能尋找邊坡的最危險滑動面并求解安全系數(shù)。

4　結(jié)論

（1）本文通過選取典型坡段，F(xiàn)LAC3D有限差分法計算的坡體的安全系數(shù)和傳統(tǒng)的極限平衡法對比可知：最終解算的安全系數(shù)和潛在滑動面的位置幾乎相同，可見基于FLAC3D有限差分法計算的坡體穩(wěn)定性是可靠和安全的。

（2）基于FLAC3D的強度折減法不需要假定滑移面的位置和形狀,能夠直觀地模擬邊坡漸進破壞過程及其滑移面的形狀、位置、范圍,同時能夠清楚的了解剪應力強弱帶的遞進分布特點,塑性區(qū)的分布和范圍及坡體位移情況，對于分析坡體穩(wěn)定性都有較為優(yōu)越的地方。

（3）坡體穩(wěn)定性的分析和后期的治理依賴于可靠地前期坡體的穩(wěn)定性分析，不同手段的相互結(jié)合和彌補對于后期的工程設計有利而無害，對于該不穩(wěn)定斜坡可以依據(jù)FLAC3D的強度折減法解算的安全系數(shù)和剪應力分布特點以及極限平衡法求解的下滑力，抗滑力的數(shù)值而采取分段，針對性治理的方案,同時可以利用兩種手段求解坡體在不同工況條件下的穩(wěn)定性和相應變形特點。

[1]丁祖全，黎志恒.蘭州市地質(zhì)災害與防治[M].蘭州：甘肅科學技術(shù)出版社，2009.

[2]馬金輝，年雁云，蔡迪花.蘭州地區(qū)滑坡因素及其與區(qū)域構(gòu)造的關系[J].自然災害學報，2006,15（3）：14～17.

[3]陳瑾.蘭州市區(qū)滑坡災害現(xiàn)狀及防治措施[J].甘肅科學學報，1999,11（2）58～62.

[4]彭文斌.FLAC3D實用教程.-北京：機械工業(yè)出版社，2007.8.

P5[文獻碼]B

1000-405X（2016）-1-393-1