嚴(yán)　飛，鄧　琴

(1.云南省交通規(guī)劃設(shè)計研究院， 云南 昆明 650041；2.陸地交通氣象災(zāi)害防治技術(shù)國家工程實驗室， 云南 昆明 650041；3.中國科學(xué)院武漢巖土力學(xué)研究所 巖土力學(xué)與工程國家重點實驗室， 湖北 武漢 430071)

?

華麗高速公路華坪岸邊坡穩(wěn)定性的DDA模擬研究

嚴(yán)飛1,2，鄧琴3

(1.云南省交通規(guī)劃設(shè)計研究院， 云南 昆明 650041；2.陸地交通氣象災(zāi)害防治技術(shù)國家工程實驗室， 云南 昆明 650041；3.中國科學(xué)院武漢巖土力學(xué)研究所 巖土力學(xué)與工程國家重點實驗室， 湖北 武漢 430071)

為了對華麗金安金沙江大橋華坪岸的邊坡穩(wěn)定性進行評價，并能夠直觀了解邊坡變形及可能的破壞過程，采用非連續(xù)變形分析方法(DDA)結(jié)合強度折減法對華坪岸z1剖面進行模擬。在計算模型中利用真實節(jié)理和虛擬節(jié)理的概念充分考慮了已有的斷層、軟弱結(jié)構(gòu)面的表征。研究結(jié)果表明：在自然狀態(tài)下華坪岸邊坡是穩(wěn)定的；在強度參數(shù)降低的情形下，華坪岸邊坡可能的破壞形式為牽引式破壞，即沿著t3凝灰?guī)r夾層和坡體前緣卸荷帶陡傾結(jié)構(gòu)面滑動。

非連續(xù)變形分析(DDA)；邊坡穩(wěn)定性；強度折減法；邊坡破壞過程

邊坡是人工邊坡和自然斜坡的統(tǒng)稱。在水利水電工程、港口工程，橋梁工程、房屋建筑及深基坑工程都會遇到邊坡穩(wěn)定問題，分析邊坡穩(wěn)定性的方法主要由極限平衡法、有限元法、離散元法、非連續(xù)變形分析方法(DDA)[1]。其中以極限平衡法和有限元法最為常用，不過極限平衡法假定巖土體為剛體，不考慮土體變形和應(yīng)力的關(guān)系，且無法分析邊坡的破壞過程[2]；有限元法則可以分析巖土體的變形及應(yīng)力[3]，且通過節(jié)理單元可處理少量的不連續(xù)面，但對于含大量的不連續(xù)面的巖土體問題則顯得無能為力。離散元方法[4-5]克服了有限元在處理大量非連續(xù)面所存在的困難，同時可以自動識別出新的接觸，可處理大位移、大轉(zhuǎn)動的問題。正是由于離散元的這些優(yōu)點，其在巖土工程領(lǐng)域獲得了廣泛的應(yīng)用，例如董金海等[6]用UDEC對三峽庫區(qū)軟硬互層高切坡崩塌進行了研究取得了較好的效果；李世海等[7]基于三維離散元剛性塊體模型研究三峽永久船閘開挖穩(wěn)定性問題；雷遠(yuǎn)見等[8]將離散元和強度折減法結(jié)合研究了邊坡穩(wěn)定性問題；胡亞東等[9]利用3DEC程序?qū)λ娬具吰麻_挖后的位移變形特征進行了分析。不過該方法存在的問題是，需要采用很小的時間步長，計算耗時巨大。非連續(xù)變形分析方法[10-12]，便于對塊體的大位移、大轉(zhuǎn)動等進行模擬，可以自動識別新的接觸，因此便于處理各種接觸問題。除此以外，DDA基于最小勢能原理，力學(xué)上非常嚴(yán)格，且其用開閉迭代進行接觸處理，接觸力計算較離散元更為嚴(yán)格準(zhǔn)確。因此該方法很適合用來分析含有軟弱夾層或者結(jié)構(gòu)面控制的邊坡穩(wěn)定狀況。本文就是采用DDA方法來對華坪岸的邊坡穩(wěn)定性進行分析。

1　華坪岸邊坡工程地質(zhì)條件

擬建金安金沙江大橋位于金沙江中游干流河段的五郎河支流至美河沖溝之間的7 km長的河道上。金安金沙江大橋主墩和錨墩分別置于金沙江兩岸斜坡之上，如圖1所示，是5個備選橋位的三維示意圖。所處位置總體呈近SN向發(fā)育，河道較為順直局部略微彎曲。河谷呈“V”型，為縱向單斜谷，兩岸地形基本對稱。因受地層產(chǎn)狀和順河向陡傾角結(jié)構(gòu)面影響，兩岸岸坡多形成階梯狀地形。華坪岸為順向坡，總體坡度20°～30°，變化較小。岸坡腰部斜向的發(fā)育一系列連續(xù)的30 m～60 m高的陡坎，陡坎的坡度可達45°以上。出露地層主要為二疊系上統(tǒng)玄武巖組上段(P2β3)，出露基巖為致密玄武巖、杏仁狀玄武巖、火山角礫熔巖夾等多層凝灰?guī)r軟弱夾層，且受后期構(gòu)造錯動影響巖體破碎，凝灰?guī)r層間卸荷帶局部存在泥化夾層。

華坪岸坡地層整體產(chǎn)狀近SN，傾向西，傾角12°～30°。同時發(fā)育有與河道近于平行和正交的兩組陡傾角節(jié)理。另外，有傾角較緩且有局部泥化的t3凝灰?guī)r軟弱夾層。出露的基巖為致密玄武巖、杏仁狀玄武巖、火山角礫熔巖夾等多層凝灰?guī)r軟弱夾層，且受后期構(gòu)造錯動影響巖體破碎，凝灰?guī)r層間局部存在泥化夾層，邊坡穩(wěn)定性可能主要受t3凝灰?guī)r軟弱夾層的控制。

華坪岸的崩塌堆積體主要集中發(fā)育于岸坡腰部的陡崖之下，其沿著向坡外緩傾的凝灰?guī)r層面形成臺階狀堆積，目前處于穩(wěn)定狀態(tài)。

上岸坡巖體以物理風(fēng)化作用為主。強風(fēng)化帶厚度一般小于15 m，局部可達20 m～25 m；中風(fēng)化帶垂直深度為35 m～65 m。岸坡巖體卸荷現(xiàn)象明顯，且淺表部位的巖體沿著陡傾節(jié)理張開，順坡向分布的凝灰?guī)r層或流面發(fā)生剪切破壞。正是由于岸坡巖體的結(jié)構(gòu)松弛，加速了巖體的風(fēng)化和岸坡的局部變形崩塌發(fā)生。

圖1工程區(qū)三維展示圖

2　DDA基本原理及邊坡模型的建立

2.1DDA基本原理

非連續(xù)變形分析方法(DDA)是石根華博士[12]提出的一種專門處理含大量非連續(xù)面的巖體問題的數(shù)值分析方法。其以塊體的形心位移(u,v)、塊體繞形心的轉(zhuǎn)動(r0)和塊體應(yīng)變(εx，εy，εxy)為基本未知量，基于最小勢能原理，形成系統(tǒng)剛度方程并求解。

2.1.1塊體內(nèi)任意一點位移的表征

塊體i內(nèi)任意一點(x,y)的位移U=(u,v)T，用完全一階近似的位移函數(shù)來表征：

U=TDi

(1)

其中：Di為塊體i的位移矢量；T為位移變換矩陣：

Di=(x,y,r,εx,εy,εxy)

2.1.2系統(tǒng)剛度方程的建立

系統(tǒng)的總勢能包括塊體的應(yīng)變能，初始應(yīng)力的勢能，外荷載作用下的勢能，慣性力勢能和黏聚力勢能。由最小勢能原理，當(dāng)勢能取最小值時系統(tǒng)達到平衡。

塊體系統(tǒng)的總勢能可寫成：

(2)

KD=F

(3)

若系統(tǒng)由n個塊體組成，則方程(3)可進一步寫為：

(4)

式中，Kij(i,j=1,2,…n)為6×6的子矩陣，當(dāng)i=j時，Kij只與塊體i的材料屬性有關(guān)，否則，Kij與塊體i和j之間的相互作用情況有關(guān)；Di代表塊體i的位移矢量；Fi代表作用于塊體i上的荷載矢量[13]。結(jié)合邊界條件和塊體系統(tǒng)的運動學(xué)條件，即可對上述方程求解，獲得各個塊體的位移、變形及應(yīng)力狀態(tài)。

2.1.3塊體系統(tǒng)的運動學(xué)條件

塊體之間相互接觸時，塊體之間應(yīng)該滿足無拉伸和侵入很小的條件。當(dāng)塊體相互接觸且有侵入發(fā)生時，需要采用較大剛度的彈簧施加于塊體之間來限制塊體的侵入。初始時，根據(jù)塊體之間的相對位置找出所有可能的接觸對，并根據(jù)接觸類型施加接觸彈簧。塊體之間有角-角接觸、角-邊接觸兩種基本接觸類型，而邊-邊接觸可轉(zhuǎn)化為兩個角-邊接觸[14]。

對角-角接觸，只需要施加一個法向彈簧，而角-邊接觸則需施加一個法向彈簧和一個切向彈簧，并在總剛度矩陣中的相應(yīng)位置加入接觸子矩陣，求解式(3)。如果塊體間的接觸力為拉力，則撤銷法向彈簧和切向彈簧，并修改剛度矩陣，重新求解式(3)。這個過程稱為開-閉迭代，基于開-閉迭代不斷調(diào)整接觸狀態(tài)，直到不再有接觸力為拉力的情形發(fā)生，便進入下一個時間步的計算，此時的接觸狀態(tài)為下一個時間步中接觸的初始狀態(tài)[13]。

2.2邊坡模型的建立

選取華坪岸的z1剖面為例(如圖1所示的z1方案)，來分析華坪岸邊坡的穩(wěn)定性。以已有的節(jié)理、軟弱夾層、斷層為界面將整個模型劃分為一系列的三角形塊體，其中真實的節(jié)理、軟弱夾層、斷層處的接觸參數(shù)(真實節(jié)理)采用該處實際的材料參數(shù)，而其余部分的接觸參數(shù)(虛擬節(jié)理[11])則采用高強度的參數(shù)，以便于表征完整巖石。整個計算模型的網(wǎng)格如圖2所示。

圖2華坪岸z1剖面塊體劃分

研究該邊坡在自然狀態(tài)的穩(wěn)定性及強度參數(shù)折減后的邊坡破壞過程。為了便于了解邊坡在不同折減系數(shù)下邊坡的位移情況，我們在邊坡的不同位置設(shè)置了3個監(jiān)測點，如圖2所示。

本文采用類似有限元中的強度折減法求穩(wěn)定系數(shù)，即將節(jié)理面的強度參數(shù)c(黏聚力)和φ(內(nèi)摩擦角)都除以折減系數(shù)K，得到一組新的強度參數(shù)值ck和φk，即：

(5)

然后用該組參數(shù)進行試算，直到邊坡達到失穩(wěn)臨界狀態(tài)，這時的折減系數(shù)K即為邊坡的穩(wěn)定系數(shù)。

2.3結(jié)構(gòu)面參數(shù)

根據(jù)現(xiàn)場勘查及設(shè)計單位提供的勘察報告可知,該邊坡可能會沿著t3凝灰?guī)r夾層滑動。涉及到塊體的滑動，因此結(jié)構(gòu)面參數(shù)的選取就顯得尤為關(guān)鍵。本次計算的參數(shù)根據(jù)基本巖石力學(xué)試驗是及反分析綜合確定，t3凝灰?guī)r夾層的參數(shù)為c=24°，φ=50kPa，其余節(jié)理面的參數(shù)取為c=35°，φ=150kPa，重度取為29kN/m3。

3　模擬結(jié)果及分析

為了得到邊坡在自然狀態(tài)的穩(wěn)定系數(shù)，采用不斷增大的折減系數(shù)進行試算，看邊坡的最終狀態(tài)是否穩(wěn)定。當(dāng)折減系數(shù)K=1.4時，計算得到的邊坡狀態(tài)、位移矢量如圖3所示。從圖3中可以明顯地看出，主要是t3凝灰?guī)r夾層上方兩個塊體發(fā)生明顯的滑動，最頂部的部分基本保持不動。由此，折減系數(shù)K=1.4時，邊坡沿著t3凝灰?guī)r夾層及坡體前緣卸荷帶陡傾結(jié)構(gòu)面出現(xiàn)明顯的滑動，前緣部分先滑動，后緣部分失去支撐，也開始滑動，華坪岸邊坡可能的破壞形式為牽引式破壞。這說明，華坪岸在自然狀態(tài)的穩(wěn)定系數(shù)約為1.4。

圖4展示的是折減系數(shù)K=1.4時圖2中的3個監(jiān)測點的位移-時間步曲線。從圖4中可知，監(jiān)測點1的位移基本保持不變，說明監(jiān)測點所在的塊體是穩(wěn)定的；監(jiān)測點2、監(jiān)測點3所在的塊體恰位于t3凝灰?guī)r夾層的上方，其位移不斷增大，沒有收斂的趨勢。說明在該折減系數(shù)下，t3凝灰夾層上方的兩個塊體失穩(wěn)，向下不斷滑動，這說明在該折減系數(shù)下，華坪岸邊坡失去穩(wěn)定。

圖3　折減系數(shù)K=1.4時華坪岸的位置狀態(tài)

圖4折減系數(shù)K=1.4時，監(jiān)測點的位移時間步曲線

圖5展示的折減系數(shù)K=1.3時3個監(jiān)測點的位移時間步曲線。由圖5可知，隨著時間步的增大，3個監(jiān)測點的位移均趨于穩(wěn)定，并未出現(xiàn)圖4所示的位移逐漸逐漸增大情況。

因此，在折減系數(shù)K=1.3時，邊坡保持穩(wěn)定狀態(tài)。

圖5折減系數(shù)K=1.3時，監(jiān)測點的位移時間步曲線

綜合上述分析比較，可認(rèn)為華坪岸的邊坡穩(wěn)定系數(shù)是1.4。

4　結(jié)　論

綜合上面的分析可以得出如下結(jié)論：

(1)DDA中結(jié)合真實節(jié)理和虛擬節(jié)理的概念可充分考慮了已有的斷層、軟弱結(jié)構(gòu)面的表征。結(jié)合強度折減法，DDA除可評價邊坡的穩(wěn)定性外，還可對邊坡變形、破壞的全過程進行預(yù)測。

(2) 對華坪岸的模擬研究結(jié)果表明，在自然狀態(tài)，華坪岸邊坡是穩(wěn)定的，在強度參數(shù)降低的情形下，華坪岸邊坡可能的破壞形式為牽引式破壞，即沿著t3凝灰?guī)r夾層和坡體前緣卸荷帶陡傾結(jié)構(gòu)面滑動。

致謝：感謝云交科2014(A)01—《金沙江中游庫區(qū)高速公路大跨徑橋址高邊坡災(zāi)害調(diào)控技術(shù)及應(yīng)用示范研究》項目的資助及課題組全體人員的配合、幫助和支持。

[1]方建瑞,朱合華,蔡永昌.邊坡穩(wěn)定性研究方法與進展[J].地下空間與工程學(xué)報,2007,3(2):343-349.

[2]FangJR,ZhuHH,CaiYC.Advancementandmethodsoftheslopestabilityresearch[J].ChineseJournalofUndergroundSpaceandEngineering, 2007,3(2):343-349.

[3]陳祖煜.土質(zhì)邊坡穩(wěn)定分析:原理方法程序[M]．北京:中國水利水電出版社,2003．

[4]鄭宏,李春光,李焯芬,等.求解穩(wěn)定系數(shù)的有限元法[J].巖土工程學(xué)報,2002,24(5):626-628.

[5]CundallPA.Acomputermodelforsimulatingprogressivelargescalemovementinblockrocksystem[C]//SymposiumofInternationalSocietyofRockMechanics.Nancy:SymposiumISRM, 1971:129-136．

[6]ItascaConsultingGroupInc.UDECusermanualversion4.0[R]. [S.l]:ItascaConsultingGroupInc., 2004.

[7]董金玉,楊繼紅,伍法權(quán),等.三峽庫區(qū)軟硬互層近水平地層高切坡崩塌研究[J].巖土力學(xué),2010,31(1):151-157.

[8]DongJY,YangJH,WuFQ,etal.Researchoncollapseofhighcuttingslopewithhorizontalsoft-hardalternantstratainThreeGorgesreservoirarea[J].RockandSoilMechanics, 2010,31(1):151-157.

[9]李世海,高波,燕琳.三峽永久船閘高邊坡開挖三維離散元數(shù)值模擬[J].巖土力學(xué)，2002，23(3)：272-277.

[10]雷遠(yuǎn)見,王水林.基于離散元的強度折減法分析巖質(zhì)邊坡穩(wěn)定性[J].巖土力學(xué),2006,27(10):1693-1698.

[11]胡亞東,肖燃,王劍梅.某水電站右壩肩傾倒變形邊坡三維離散元數(shù)值分析[J].地質(zhì)災(zāi)害與環(huán)境保護,2014,25(4):51-55.

[12]ShiGH.Discontinuousdeformationanalysis-anewnumericalmodelforthestaticsanddynamicsofblocksystem[D].Berkeley：UniversityofCalifornia-Berkeley, 1988.

[13]張秀麗．?dāng)嗬m(xù)節(jié)理巖體破壞過程的數(shù)值分析方法研究[D].武漢：中國科學(xué)院武漢巖土力學(xué)研究所，2007．

[14]焦玉勇,張秀麗,劉泉聲，等.用非連續(xù)變形分析方法模擬巖石裂紋擴展[J].巖石力學(xué)與工程學(xué)報,2007,26(4):682-691.

DDA Simulation of Huaping Shore Slope Stability of Huali Expressway

YAN Fei1,2， DENG Qin3

(1.EngineeringSurveyBranch,BroadvisionEngineeringConsultants,Kunming,Yunnan650041,China；2.NationalEngineeringLaboratoryforSurfaceTransportationMeterologicalHazardPreventionTechnology,Kunming,Yunnan650041,China； 3.StateKeyLaboratoryofGeomechanicsandGeotechnicalEngineering,InstituteofRockandSoilMechanics，ChineseAcademyofSciences，Wuhan,Hubei430071,China)

In order to evaluate the slope stability of Huaping shore of the Jinsha River Bridge, and intuitively understand the slope deformation and possible damage process, we use the discontinuous deformation analysis (DDA) method combined with strength reduction to simulate the slope of Huaping shore. Based on the engineering geological characteristics of Huaping shore, a calculation model of z1 profile was developed. In addition the concept of real joint and virtual joints joint was used, the faults or weak structure plane may take into consideration directly, the fracture of the rock itself also can be model naturally. The results show that: in the natural state, Huaping shore slope is stable. When the strength parameters reduced, Huaping shore slope is possible failure due to traction damage, the slope may slip along t3 tuff sandwich and steep joint, after that the small block behind the big block also began to slip because of loss support.

discontinuous deformation analysis (DDA); slope stability; strength reduction; slope failure

10.3969/j.issn.1672-1144.2016.04.022

2016-03-27

2016-04-26

云南省交通運輸廳科技項目(云交科2014(A)01)

嚴(yán)飛(1979—)，男，云南昆明人，高級工程師，主要從事公路工程方面的工作。E-mail：yfmail621@163.com

O39；TU457

A

1672—1144(2016)04—0107—04