匡俊，陳芋如，朱仁民，曹建林，居俊，唐強*

（1.江蘇蘇州地質(zhì)工程勘察院，江蘇蘇州215011；2.江蘇地質(zhì)礦產(chǎn)局第四地質(zhì)大隊，江蘇蘇州215004；3.蘇州大學，江蘇蘇州215131；4.常熟理工學院，江蘇蘇州215500）

0 引言

近年來，隨著計算機技術(shù)的飛速發(fā)展，滑坡等地質(zhì)災害的演化規(guī)律及控制取得了一系列重要成果和突破。低成本高效率的仿真技術(shù)發(fā)展和應用對于預防和控制此類災害有顯著效果[1-4]。呂文斌等[5]以西寧市張家灣某滑坡為研究對象，結(jié)合極限平衡法以及Midas/GTS有限元分析方法，研究了滑坡形成機理和發(fā)展過程，定量評價了該滑坡3種不同工況條件下的穩(wěn)定狀態(tài)。李效萌等[6]采用折線滑動法計算滑坡穩(wěn)定性系數(shù)等方法，對滑坡在2種工況下的穩(wěn)定性進行了定量探討，在此基礎上分析了滑坡成因和可能的失穩(wěn)破壞模式。劉洪波等[7]基于有限差分方法建立滑坡三維模型，分析了邊坡剖面、剪切帶處土顆粒位移場分布后發(fā)現(xiàn)：當土體完全失穩(wěn)時剪切帶完全貫通，滑移體內(nèi)土體向下加速滑動?？苷鐫萚8]對秦嶺某滑坡研究分析后表明：飽和工況下坡體將失穩(wěn)，坡體失穩(wěn)時基巖面附近滑動面最先發(fā)生剪切破壞。

離散元技術(shù)因其能夠有效模擬滑坡機理而引起眾多學者的關(guān)注。王燾[9]針對滑坡堵江預判問題，分析了滑坡體方量與堵江程度相關(guān)關(guān)系。一些學者[10-17]對滑坡破壞方式及發(fā)展過程進行模擬研究，結(jié)果與實際情況基本吻合。

本文應用離散元技術(shù)模擬鳳凰山棄土場邊坡破壞過程，對關(guān)鍵部位的位移、應力進行分析，研究其變形破壞機理。

1 工程概況

蘇州市科技城鳳凰山棄土場位于高新區(qū)元方路和逢春路交叉口西側(cè)，原鳳凰山采石宕口。場區(qū)所在地的原始地貌屬蘇州西部低山丘陵，山體高度在10～130 m，自然坡度在18°～25°。經(jīng)采石開挖后宕口邊坡發(fā)生較大變化，邊坡多在14°～58°之間，局部呈小范圍直立狀，宕口地形西高東低，宕口填土邊坡分為南北兩部分。南側(cè)邊坡坡度較北坡緩，局部較陡，在其坡腳有水滲出，形成溪流，該坡局部存在明顯的雨水沖刷坍塌現(xiàn)象。北側(cè)填土邊坡整體為一坡到底，局部見明顯的蠕滑破裂縫隙，寬度在5～20 mm左右。南、北邊坡面未見巖石裸露，無滾石、崩塌等不良地質(zhì)作用?？碧狡陂g測得場地各勘探孔孔口高程在14.20～44.95 m之間，地勢起伏大。圖1為邊坡的地質(zhì)剖面圖。

圖1 邊坡地質(zhì)剖面圖Fig.1 Geological profile of the slope

由圖1可知，邊坡土層間有成分為粉質(zhì)黏土的軟弱夾層，粉質(zhì)黏土構(gòu)成邊坡失穩(wěn)變形的潛在滑動面。邊坡前緣開挖卸荷引起坡面臨空，應力場發(fā)生變化，形成的不穩(wěn)定塊體極有可能發(fā)生滑移破壞。當變形達到一定程度時，坡面將產(chǎn)生裂縫，發(fā)生失穩(wěn)破壞。由于松散體為均質(zhì)體，各個顆粒在各方向的抵抗能力相同，滑動沿著力矩相同的軌跡滑出，形成圓弧形滑動面。

圖2為鳳凰山現(xiàn)場勘探結(jié)果。由于碎石上部覆蓋大量植被及土壤，土體的黏聚力顯著增大，這延緩了邊坡失穩(wěn)破壞需要的時間，破壞形式表現(xiàn)為漸進式的蠕動變形。邊坡上出現(xiàn)的多處裂縫表明，該處邊坡正處于滑移過程之中。一旦受到雨水沖刷等外界作用，邊坡將迅速發(fā)生破壞，有可能造成安全隱患。為近一步研究邊坡的穩(wěn)定性，本文利用離散元軟件對邊坡的破壞過程展開數(shù)值分析。

圖2 現(xiàn)場踏勘資料Fig.2 Site investigation data

2 邊坡破壞過程數(shù)值模擬

2.1 模型建立

采用離散元技術(shù)模擬邊坡滑動，在建模過程中，ball-ball接觸與ball-facet接觸會貫穿始終，線性模型具有較好的模擬效果。根據(jù)杜永斌[18]的研究，線性模型接觸力通過成對出現(xiàn)并相接觸顆粒之間恒定法向和切向剛度的線性彈簧產(chǎn)生，假設法向剛度與切向剛度之比為1，兩種剛度的計算公式如下：

式中：Kn為兩個接觸顆粒間的法向接觸剛度；Ks為切向接觸剛度；A為顆粒間的接觸面積；E為楊氏模量；L為兩個顆粒的半徑之和。由于本文中的地下水處于滑移面下方，根據(jù)何倩等[16]和宋浩燃等[17]的研究，此時地下水對于邊坡滑移影響較小，因此在本次模擬中忽略其影響。本次模擬主要接觸參數(shù)具體取值如表1所示。

表1 碎石土的細觀參數(shù)Table 1 Mesoscopic parameters of gravelly soil

模擬具體步驟為：成樣，預壓，加自重以及切坡。首先采用預壓來釋放試樣中過大的內(nèi)應力，隨后放大重力場。在程序中經(jīng)過多次重力加速度的調(diào)整，最終在80倍的重力加速度下，顆粒能夠向下自然滑落。在加自重前進行了土體平衡，并定義了邊坡底部坐標。

顆粒流模型計算的下部邊界是有限的。顆粒流模型高度應大于地基的塑性區(qū)開展深度，即H＜Zmax，塑形區(qū)的最大開展深度為：

式中：p為地基的條形均布荷載；γ為土體容重；h為邊坡高度；H為模型高度；c為土體的黏聚力；φ為土體的內(nèi)摩擦角。

為確?；瑒蛹羟忻娉霈F(xiàn)在邊坡內(nèi)，邊坡模型最終滿足H≥1.5h的條件。根據(jù)式（3），建立模型長40 m，高20 m，顆粒最小半徑為0.05 m，最大半徑為0.1 m，顆?？倲?shù)為64 892個。最終得到滑坡數(shù)值計算模型。圖3為滑坡數(shù)值計算模型及監(jiān)測點。

圖3 滑坡數(shù)值計算模型及監(jiān)測點Fig.3 Numerical landslide model and monitoring points

2.2 臨界坡腳確定

鳳凰山棄土場的邊坡多在14°～58°之間，經(jīng)過程序多次模擬，當坡角為30°時，邊坡未發(fā)生滑動位移。當坡角為35°時，邊坡有大部分顆粒滑落。而在選擇32°坡角時，邊坡只有少量顆粒發(fā)生了較大位移，邊坡達到極限平衡狀態(tài)，即該邊坡發(fā)生滑動的臨界坡角為32°。圖4為臨界坡角的計算云圖。

圖4 臨界坡角的計算云圖Fig.4 Calculation clouds for critical slope angles

2.3 坡度最大處失穩(wěn)過程分析

2.3.1 位移矢量分析

經(jīng)實地踏勘，現(xiàn)場坡角最大值為58°。由于蠕動變形，邊坡坡頂已出現(xiàn)明顯的裂縫。通過數(shù)值模擬對該處邊坡的失穩(wěn)過程進行分析。

位移是邊坡發(fā)生破壞的直觀表現(xiàn)，由圖5可知，該邊坡的失穩(wěn)模式是蠕動性滑坡。計算初始，邊坡開始變形，靠近坡面的顆粒逐漸滑落，巖體蠕動變形。隨著時間的推移，顆?；瑒游恢猛聝?nèi)移動；當計算到10 s時，剪切帶基本貫通，形成穩(wěn)定斜坡，滑動面呈圓弧狀，坡趾圓也基本形成；當t=30 s時最終平衡，顆粒的最大位移達10.4 m。

圖5 邊坡破壞過程模擬Fig.5 Simulation of the slope failure process

2.3.2 應力變化分析

在斜面上設置3個監(jiān)測點，以觀察對比斜坡頂部、中部以及底部的應力變化。監(jiān)測點的具體位置以及名稱序號如圖3所示。監(jiān)測曲線圖如圖6所示，分別監(jiān)測斜坡x方向正應力（σx），剪應力（τxy），以及y方向的正應力（σy）。

圖6 應力變化曲線Fig.6 Curves of stress changes

由圖6可知，監(jiān)測點1處的顆粒在模擬過程中位于坡體的頂部，并不斷發(fā)生向下的位移，因此此時的應力曲線呈總體減小趨勢。

相反，監(jiān)測點2處的x，y方向正應力以及剪應力均呈上升趨勢，且y方向的正應力與x方向正應力相比較大。這是因為斜坡中部隨著變形量的增加，不斷受到碎石土向下的擠壓，顆粒不斷堆積，導致y向應力逐漸增加，峰值應力為115.24 kPa。

坡腳處的監(jiān)測點在滑動初始階段應力驟降，隨后趨于平緩，這是因為坡腳處在一開始的滑動過程中，顆粒間不斷調(diào)整位置，發(fā)生剪切破壞。隨著上覆土體的不斷滑下，監(jiān)測點1處土體被上覆土體覆蓋，坡體滑動對其影響減小，坡腳處達到應力平衡。剪應力在整個滑動過程中變化幅度不明顯，坡腳的剪應力最大值為70.84 kPa。

2.3.3 結(jié)果討論

根據(jù)瑞典條分法，邊坡的安全系數(shù)可以由公式（4）進行計算：

式中：NWi為孔隙水壓力；TDi為滲透壓力產(chǎn)生的平行滑動面分力；RDi為滲透壓力產(chǎn)生的垂直滑動面分力；Wi為土體自重；Ci為土體內(nèi)聚力；φi為土體內(nèi)摩擦角；Li為滑動面長度；αi為滑動面傾角；A為地震加速度；Kf為穩(wěn)定系數(shù)。

根據(jù)石崇[19]的研究，植被根系及黏土含量的增加會增大邊坡填土的黏聚力。公式（4）表明，黏聚力的增大會顯著提高土體的抗滑力，提升邊坡的穩(wěn)定性。在黏聚力的作用下，該邊坡緩慢發(fā)生蠕動變形。一旦受到雨水沖刷等外界作用，土體的黏聚力將迅速減小，邊坡容易發(fā)生失穩(wěn)破壞。

3 結(jié)語

本文通過離散元技術(shù)對蘇州市科技城鳳凰山棄土場邊坡穩(wěn)定性進行了模擬，得到以下結(jié)論。

1）邊坡達到32°的坡角時，達到極限平衡狀態(tài)。從現(xiàn)場勘探情況看，最大坡角處邊坡已經(jīng)出現(xiàn)明顯的裂縫，可以通過削方減載、擋土墻法、抗滑樁法、格構(gòu)錨桿法等來實施滑坡整治措施。

2）模擬現(xiàn)場最大坡角58°處的邊坡過程中，隨著時間的推移，顆?；瑒拥奈恢孟蛐逼聝?nèi)部移動，滑動剪切帶基本形成并貫通以后，顆粒發(fā)生塑形流動，形成穩(wěn)定的圓弧形滑坡面。

3）坡體中部顆粒的正應力由于不斷受到上覆顆粒的擠壓呈非線性增長，坡腳位置處的應力呈總體下降趨勢，隨著時間的推移最終達到平衡。