李建舉 李玉瑩

(1.中國航空港建設(shè)第九工程總隊，成都 611430； 2. 長安大學(xué) 公路學(xué)院，西安 710064)

邊坡穩(wěn)定分析中圓弧與折線滑動面對比分析研究

李建舉1李玉瑩2

(1.中國航空港建設(shè)第九工程總隊，成都 611430； 2. 長安大學(xué) 公路學(xué)院，西安 710064)

以分析圓弧和折線兩種滑動面以及簡化Bishop法和通用條分法兩種方法對邊坡穩(wěn)定分析結(jié)果的影響為目標(biāo)，基于極限平衡理論，采用邊坡穩(wěn)定分析軟件Zslope，以均質(zhì)邊坡、加筋邊坡以及成層土非均質(zhì)邊坡為例，計算中采用圓弧和折線兩種滑動面，以及簡化Bishop法和通用條分法兩種極限平衡方法，分析比較其安全系數(shù)和臨界滑動面的差異．研究結(jié)果表明，當(dāng)滑動面為圓弧滑動面時，簡化Bishop法與通用條分法計算得到的安全系數(shù)和臨界滑動面位置基本相同；與圓弧滑動面相比，折線滑動面得到的安全系數(shù)偏小，對于加筋邊坡以及非均質(zhì)邊坡，折線與圓弧滑動面搜索到的安全系數(shù)差別可能會較大，圓弧滑動面容易陷入局部極值，而折線滑動面則更有利于搜索到全局極值．研究成果對于邊坡工程設(shè)計具有一定的參考價值．

邊坡穩(wěn)定； 圓弧滑動面； 折線滑動面； 安全系數(shù)； 通用條分法

邊坡穩(wěn)定分析是巖土工程領(lǐng)域研究的熱點問題，其分析方法也一直是學(xué)者們研究的重點之一．自從瑞典人彼德森最早提出極限平衡條分法的概念之后[1]，許多學(xué)者對于極限平衡理論及其在邊坡穩(wěn)定分析中的應(yīng)用展開研究，豐富了這一領(lǐng)域的研究成果．基于極限平衡理論的條分法仍然是目前為止工程中常用的邊坡穩(wěn)定分析方法，為廣大工程師所熟悉．

在瑞典條分法之后的研究過程中出現(xiàn)了一系列數(shù)值分析算法，主要包括適用于非圓弧滑動面的通用條分法、更適用于實際工程的傾斜條分法、能夠準(zhǔn)確進行理論計算的三維極限平衡法等[2]．盧坤林等[3]將空間內(nèi)土條劃分后條柱間的作用力看作滑動面上的正應(yīng)力，根據(jù)整個滑動土體的靜力平衡條件建立方程，提出了一種適用于一般三維滑動面形式的極限平衡法；白桃等[4]為了考察空間內(nèi)土體參數(shù)變異時極限平衡法的適應(yīng)性，對極限平衡法4種不同算法下得到的安全系數(shù)進行了比較，并提出采用安全系數(shù)和破壞概率2個指標(biāo)能夠定量化地表述邊坡穩(wěn)定性．王旭等[5]考慮到土體參數(shù)變異性和隨機性對邊坡穩(wěn)定性的影響，提出將可靠度方法得到的可靠度和破壞概率等指標(biāo)與確定性方法相結(jié)合，可以更好地評價邊坡的穩(wěn)定性；鄧東平等[6]基于摩爾-庫倫強度理論和滑動面法向力的計算，根據(jù)空間內(nèi)土體的靜力平衡條件，提出對于三維滑動土體可以采用不完全滿足所有平衡條件的準(zhǔn)嚴(yán)格法求解安全系數(shù)．代凌輝等[7]采用有限元軟件計算分析了抗滑樁的樁長和樁徑等參數(shù)對于抗滑樁加固后的邊坡的安全系數(shù)和臨界滑動面的影響．但是以往的研究主要集中在極限平衡分析方法對于計算結(jié)果的影響上，滑動面的形狀對于邊坡安全系數(shù)以及臨界滑動面的影響仍無系統(tǒng)研究，對于均質(zhì)土坡研究較多，對于加筋邊坡以及非均質(zhì)邊坡的研究較少，本文基于邊坡穩(wěn)定分析軟件ZSlope，分別設(shè)置圓弧滑動面和折線滑動面為初始滑動面，采用極限平衡法中的簡化Bishop法和通用條分法兩種算法計算，通過均質(zhì)土邊坡、加筋邊坡和非均質(zhì)土邊坡算例的計算，比較不同算法得到的安全系數(shù)和臨界滑動面，歸納總結(jié)相應(yīng)規(guī)律，為邊坡穩(wěn)定分析方法的選取提供參考．

1 基本理論

極限平衡法是將滑動土體劃分為若干個垂直土條后，分析所劃分土條的力和力矩平衡條件，根據(jù)靜力平衡方程進行分析和計算．該方法假設(shè)滑動面發(fā)生滑動造成土體破壞，滑動土體服從破壞條件，且土體內(nèi)滑動面已知[8]．不同的算法采用不同的土條間相互作用關(guān)系，也適用于不同的滑動面形狀．根據(jù)條間力假定的不同，極限平衡法可分為兩大類：一類是只滿足力與力矩平衡的部分靜力平衡條件的方法，例如瑞典條分法、簡化Bishop法、簡化Janbu法，另一類是嚴(yán)格滿足力與力矩平衡的所有靜力平衡條件的方法，如Spencer法、Morgenstern-Price法、通用條分法等．

1.1 簡化Bishop法

最初畢肖普提出將土體的抗剪強度指標(biāo)，即粘聚力c和內(nèi)摩擦角φ的正切值tanφ分別降低為c/F和tanc/F，然后根據(jù)靜力平衡條件建立方程，其中F即為穩(wěn)定安全系數(shù)；簡化Bishop法則是以圓弧滑動面為基礎(chǔ)，假設(shè)只考慮土條間的水平作用力，忽略切向剪應(yīng)力的影響，根據(jù)滑動面的靜力平衡條件求解安全系數(shù)[9]．

1.2 通用條分法

通用條分法滿足條間力、力與力矩平衡的所有靜力平衡條件，并將其合力的大小和方向作為未知數(shù)，通過多次定義條間力的傾角函數(shù)，經(jīng)數(shù)值迭代求解后，判斷是否滿足靜力平衡條件進而判定其安全系數(shù)，該方法既適用于圓弧滑動面，也適用于折線滑動面[9-11]．

1.3 ZSlope軟件簡介

ZSlope是一款基于極限平衡理論的邊坡穩(wěn)定分析軟件[12]．該軟件系統(tǒng)不僅提供了規(guī)范中所有安全系數(shù)求解算法，如瑞典法、簡化Bishop法、簡化Janbu法、不平衡推力法、通用條分法GLE、Morgenstern-Price法、Spencer法等．軟件提供了5種臨界滑面搜索算法，其中，圓弧滑動面搜索方法為窮舉法[12]，折線形滑動面搜索方法采用常值轉(zhuǎn)角法、隨機轉(zhuǎn)角法、節(jié)點移動法和變常值轉(zhuǎn)角法[13-15]．本文對于折線滑動面的搜索方法是改進的Monte-Carlo搜索算法，即常值轉(zhuǎn)角法，該方法把隨機轉(zhuǎn)角變?yōu)榱司哂懈蟮乃阉鞣秶透叻€(wěn)定性的常值轉(zhuǎn)角，使之受初始滑面的影響更小，因此具有更強的適應(yīng)性[15]．

2 算例分析

本文以均質(zhì)邊坡、加筋邊坡以及成層土非均質(zhì)邊坡為例，對比分析不同形狀的初始滑動面和不同的算法對安全系數(shù)及臨界滑動面的影響．

2.1 均質(zhì)邊坡

2.1.1 基本計算模型及參數(shù)

首先分析一均質(zhì)邊坡，坡比為1∶0.25，邊坡高度12 m，地基厚度6 m，其基本計算模型尺寸如圖1所示，土體參數(shù)見表1．

圖1 算例1計算模型(單位：m)

算例容重γ/kN·m-3飽和容重γsat/(kN·m-3)內(nèi)摩擦角φ/°粘聚力c/kPa1、218．020．020．035．0

以均質(zhì)邊坡為例的條分法劃分如圖2～3所示，土條劃分?jǐn)?shù)n為25條．

圖2 圓弧滑動面條分法劃分簡圖(單位：m)

圖3 折線滑動面條分法劃分簡圖(單位：m)

2.1.2 結(jié)果分析

經(jīng)分析計算后，得到的安全系數(shù)及其相應(yīng)臨界滑動面(見表2和圖4)．

表2 算例1安全系數(shù)計算結(jié)果

圖4 算例1臨界滑動面計算結(jié)果(單位：m)

由上述計算結(jié)果可知，初始滑動面設(shè)置為圓弧滑動面時，簡化Bishop法與通用條分法計算得到的安全系數(shù)和臨界滑動面位置基本相同；與圓弧滑動面相比，折線滑動面得到的安全系數(shù)偏小，原因在于折線滑動面沒有規(guī)定具體的滑動面形狀，在搜索過程中可以找到更為多樣的潛在滑動面，這與文獻[16]的研究結(jié)果是相同的．

為了驗證上述結(jié)果的合理性和可行性，進而給出下列兩種算例進行分析．

2.2 加筋邊坡

2.2.1 基本計算模型及參數(shù)

算例2為一土工格柵加筋邊坡，即在算例1的基礎(chǔ)上增加拉筋帶，筋材長度為12 m，基本計算模型尺寸如圖5所示，土體參數(shù)與算例1相同(見表1)．筋材的抗拉屈服強度為80 kN/m，計算分析中將筋材簡化為力增加到土條的力與力矩平衡分析中，計算方法詳見文獻[17]．

圖5 算例2計算模型(單位：m)

2.2.2 結(jié)果分析

經(jīng)分析計算后，得到的安全系數(shù)及其相應(yīng)臨界滑動面(見表3和圖6)．

表3 算例2安全系數(shù)計算結(jié)果

圖6 算例2臨界滑動面計算結(jié)果(單位：m)

由算例2的分析結(jié)果可知，邊坡加筋后的安全系數(shù)總體提高，同時也滿足在圓弧滑動面下兩種算法安全系數(shù)和滑動面十分接近；折線滑動面的安全系數(shù)比圓弧滑動面小，滑動面有所差別但總體趨勢一致．

2.3 成層土非均質(zhì)邊坡

2.3.1 基本計算模型及參數(shù)

算例3為延安機場遷建工程的非均質(zhì)邊坡，遷建項目場區(qū)處于延安市的東南方，地處寶塔區(qū)柳林鎮(zhèn)，離市中心約15 km．整體而言，該區(qū)西北向偏低，東南向偏高，場區(qū)的山體大多被沖溝沖刷為形狀各異的侵蝕和剝蝕黃土峁與黃土梁，起伏不平，深切綿延．工程現(xiàn)場的山梁和山峁寬從250 m到450 m不等，梁頂和峁頂?shù)牡貏菹鄬Ρ容^平坦，溝谷寬從75 m至300 m之間變化，谷坡陡峭嶙峋，溝谷底部與梁峁頂部相差大約147 m．依據(jù)勘測資料和鉆探結(jié)果，地層從上到下順次是植物層、坡洪積粉質(zhì)黏土、第四系上更新統(tǒng)風(fēng)積黃土、殘積古土壤、中更新統(tǒng)風(fēng)積黃土和殘積古土壤．具體土層劃分及基本計算模型尺寸如圖7所示，土體參數(shù)見表4．

圖7 算例3計算模型(單位：m)

算例地層容重/(kN·m-3)內(nèi)摩擦角φ/°粘聚力c/kPa3①填土16．430．050．0②黃土14．620．030．0③古土壤16．120．028．0④黃土14．821．028．0⑤古土壤16．820．030．0⑥黃土16．820．035．0⑦黃土19．020．035．0⑧黃土19．030．070．0⑨填方邊坡17．630．050．0

2.3.2 結(jié)果分析

經(jīng)分析計算后，得到的安全系數(shù)及其相應(yīng)臨界滑動面(見表5和圖8)．

表5 算例3安全系數(shù)計算結(jié)果

圖8 算例3臨界滑動面計算結(jié)果(單位：m)

算例3依舊保留上述結(jié)果分析的規(guī)律，但折線滑動面與圓弧滑動面計算得到的安全系數(shù)相差較大，其主要原因是，由于折線滑動面和圓弧臨界滑動面的搜索方法不同，在成層土的情況下，折線和圓弧的臨界滑動面差異較大，因此其安全系數(shù)差別較大．說明當(dāng)土體介質(zhì)較為復(fù)雜時，圓弧滑動面在搜索過程中容易陷入局部極值，而折線滑動面容易搜索到全局極值．

3 結(jié) 論

1)通過算例對比分析，當(dāng)滑動面為圓弧滑動面時，簡化Bishop法與通用條分法計算得到的安全系數(shù)和臨界滑動面位置基本相同．

2)不同初始滑動面形式下得到的臨界滑動面有一定的差異，但總體趨勢相接近；對于幾何形狀復(fù)雜的非均質(zhì)邊坡而言，只有初始滑動面充分接近全局臨界滑動面時，搜索過程中才不易陷入局部極值．

3)由于折線滑動面沒有規(guī)定具體的滑動面形狀，因此可以搜索到更為多樣和合理的滑動面，故其安全系數(shù)比圓弧滑動面方法計算的小，對于均質(zhì)土坡，折線滑動面與圓弧滑動面的安全系數(shù)和臨界滑動面相差不大，但對于加筋邊坡與非均質(zhì)邊坡，兩者有一定的差別，此時應(yīng)用折線滑動面更容易搜索到全局極值．

[1] 李領(lǐng)先．基于極限平衡法的邊坡穩(wěn)定性計算[J]．科技信息，2009(12)：614-615．

[2] 蘇利軍．基于極限平衡法和有限元法的土質(zhì)邊坡穩(wěn)定分析研究[D]．成都：西華大學(xué)，2012．

[3] 盧坤林，朱大勇，甘文寧，等．一種邊坡穩(wěn)定性分析的三維極限平衡法及應(yīng)用[J]．巖土工程學(xué)報，2013，35(12)：2276-2282．

[4] 白 桃，黃曉明．考慮土體參數(shù)空間變異性的邊坡極限平衡法求解[J]．西南交通大學(xué)學(xué)報，2014，49(4)：662-667．

[5] 王 旭，劉東升，宋強輝，等．基于極限平衡法的邊坡穩(wěn)定性可靠度分析[J]．地下空間與工程學(xué)報，2016，12(3)：839-844．

[6] 鄧東平，李 亮．基于滑動面應(yīng)力假設(shè)下的三維邊坡穩(wěn)定性極限平衡法研究[J]．巖土力學(xué)，2017，38(1)：189-196．

[7] 代凌輝，張 胤，李 建．抗滑樁在加固邊坡穩(wěn)定性中的應(yīng)用及影響分析[J]．水利水電技術(shù)，2016，47(12)：52-55，65．

[8] 劉金龍，陳陸望，王吉利．邊坡穩(wěn)定性分析方法簡述[J]．水電能源科學(xué)，2008(1)：133-137．

[9] 崔 蕾．簡化畢肖普法在不同安全系數(shù)體系下應(yīng)用的對比分析[J]．赤峰學(xué)院學(xué)報：自然科學(xué)版，2008，24(5)：105-108．

[10] Fredlund D G, Krahn J． Comparison of Slope Stability Methods of Analysis[J]． Canadian geotechnical journal， 1977, 14(3)：429-439．

[11] Chugh A K． Variable Interslice Force Inclination in Slope Stability Analysis[J]． Soils and foundations, 2008, 26(1)： 115-121．

[12] 張魯渝．一個用于邊坡穩(wěn)定分析的通用條分法[J]．巖石力學(xué)與工程學(xué)報，2005，24(3)：496-501．

[13] 張魯渝．邊坡穩(wěn)定分析軟件ZSlope的開發(fā)[J]．巖石力學(xué)與工程學(xué)報，2004，23(16)：2830-2835．

[14] 張魯渝，張建民．基于MonteCarlo技術(shù)的臨界滑面搜索算法的實現(xiàn)及改進[J]．巖土工程學(xué)報，2006,28(7)：857-862．

[15] 張魯渝，張建民．基于滑面段旋轉(zhuǎn)的Monte Carlo搜索技術(shù)及改進(I)：轉(zhuǎn)角大小隨機[J]．巖土力學(xué)，2006,27(11)：1902-1908．

[16] 張魯渝，張建民．基于滑面段旋轉(zhuǎn)的Monte Carlo搜索技術(shù)及改進(Ⅱ)：常值轉(zhuǎn)角[J]．巖土力學(xué)，2006,27(12)：2197-2202．

[17] 肖 杰，吳 宏，張錦濤．折線型滑裂面邊坡的穩(wěn)定性問題研究[J]．工程與試驗，2015，55(3)：26-29．

[18] Song F, Chen R Y, Ma L Q, et al． A New Method for the Stability Analysis of Geosynthetics-reinforced Slopes[J]． Journal of mountain sciences, 2016, 13(11)：2069-2078．

ComparisonbetweenArcandPolylineFailureSurfacesinSlopeStabilityAnalysis

Li Jianju1Li Yuying2

(1. The 9th General Team of Airport Construction of China, Chengdu 611430, China; 2. School of Highway Engineering, Chang'an University, Xi'an 710064, China)

In this paper, the effects of two kinds of sliding surfaces of arc and polyline and simplified Bishop method and generalized limit equilibrium method (GLE) for slope stability analysis are studied.Taking the homogeneous slope, the reinforced slope and the stratified soil heterogeneous slope as examples, the circular and polyline sliding surfaces are calculated by two algorithms to analyze the difference between the safety factor and the critical slip surface. The results show that when the sliding surface is an arc sliding surface, the simplified Bishop method and the generalized limit equilibrium method have the same safety factor and critical slip surface. Compared with the arc sliding surface, the safety factor obtained by the polyline sliding surface is smaller. It is shown that the arc is easy to fall into local extremum for complex slopes; while the polyline sliding surface is easy to search for global extremum. The study results can provide theoretical references for the design of engineering slopes.

slope stability; arc sliding surface; polyline sliding surface; factor of safety; general limit equilibrium method

10.13393/j.cnki.issn.1672-948X.2017.05.004

2017-08-08

住房與城鄉(xiāng)建設(shè)部研究開發(fā)項目(2015-K2-008)

李建舉(1962-)，男，高級工程師，碩士研究生，研究方向為機場工程及管理．E-mail：Lijj1962@outlook.com

TU43

A

1672-948X(2017)05-0017-05

[責(zé)任編輯周文凱]

土木水電論壇欄目征文

《三峽大學(xué)學(xué)報(自然科學(xué)版)》是科技部中國科技論文統(tǒng)計源期刊和俄羅斯《文摘雜志》、美國《化學(xué)文摘》、英國《INSPEC》收錄期刊，并已加入《中國學(xué)術(shù)期刊(光盤版)》和“中國期刊網(wǎng)”，在水利電力行業(yè)有一定的影響力和知名度。本刊自2004年開設(shè)了“水電論壇”欄目(現(xiàn)更名為“土木水電論壇”)，主要刊載反映國內(nèi)外土木水電建設(shè)中的最新科研成果和工程應(yīng)用的論文。歡迎土木水電工程科研和建設(shè)單位的廣大作者踴躍投稿，我們將以最大的熱情和最快的速度為您服務(wù)。對您所投的“土木水電”類論文，屬省、部級以上基金資助的，我刊將免收版面費，且稿酬從優(yōu)！