非線性Mohr—Coulomb破壞準則下邊坡可靠度上限
2017-03-29 12:29:06賀志軍曹吉趙煉恒瞿召乾楊勝博
土木建筑與環(huán)境工程 2016年6期
賀志軍 曹吉 趙煉恒 瞿召乾 楊勝博
摘要：傳統(tǒng)邊坡可靠度分析往往在巖土參數(shù)服從線性Mohr-coulomb（簡稱線性M-c）破壞準則的假設條件下進行，并且常常采用極限平衡法或有限元法計算安全系數(shù)。然而，巖土介質(zhì)破壞準則具有一定的非線性。為能更加實際地描述巖土破壞機理和得到嚴格精確的解，基于非線性Mohr-coulomb（簡稱非線性M-c）破壞準則，結(jié)合極限分析上限法和蒙特卡洛法，進行邊坡可靠度上限分析。當非線性參數(shù)m=1時，與等效的線性M-c破壞準則進行對比計算，驗證了方法的可行性。同時，將初始粘聚力、內(nèi)摩擦角arctan（Co/σt）和非線性參數(shù)作為隨機變量且服從截斷正態(tài)分布，進行了參數(shù)變異性和敏感性影響分析。研究表明：非線性M-c破壞準則下，邊坡可靠度隨初始粘聚力、內(nèi)摩擦角arctan（Co/σt）和非線性參數(shù)變異性的增大而減小；邊坡可靠度隨初始粘聚力和內(nèi)摩擦角arctan（Co/σt）的增大而增大，隨非線性參數(shù)的增大而減小。
關鍵詞：邊坡；可靠度；破壞準則；極限分析；蒙特卡洛法
中圖分類號：TU 457
文獻標志碼：A
文章編號：1674-4764（2016）06-0001-09
邊坡穩(wěn)定性問題一直是巖土工程的一個重要研究內(nèi)容。目前，邊坡穩(wěn)定性分析主要有兩大體系：確定性體系與不確定性體系（可靠度體系）。確定性體系使用極限平衡法、數(shù)值模擬法或極限分析法等方法分析邊坡的穩(wěn)定性，求得邊坡最小安全系數(shù)，以此作為邊坡穩(wěn)定性評價指標。然而，邊坡是一個極其復雜的系統(tǒng)，巖土參數(shù)具有明顯的隨機性，采用確定性體系分析邊坡穩(wěn)定性不符實際。文獻[2-3]也指出：由于安全系數(shù)沒有考慮參數(shù)隨機性和離散型對結(jié)果的影響，導致實際工程中很多結(jié)構(gòu)在滿足安全系數(shù)的條件下依然出現(xiàn)了破壞現(xiàn)象。以概率論為基礎的可靠度體系可考慮邊坡系統(tǒng)內(nèi)部的隨機關系，可給出邊坡穩(wěn)定程度，因而可以彌補用單一安全系數(shù)分析邊坡穩(wěn)定性的局限性。
邊坡可靠度分析主要有兩大步驟：一是構(gòu)建計算邊坡穩(wěn)定性安全系數(shù)的模型（功能函數(shù)）；二是使用可靠度分析方法計算邊坡可靠度（失效概率或可靠度指標）。在構(gòu)建模型方面，目前應用最廣泛的是極限平衡法，該方法理論簡單、易于實施，但所作假設較多，根據(jù)塑性理論可知，該方法所獲解答不是嚴格的上下限解。另一應用廣泛的方法是數(shù)值方法，該方法將有限元技術(shù)應用到邊坡穩(wěn)定性分析中，可以考慮土體與其中結(jié)構(gòu)物的共同作用，但其所得到的極限荷載值仍不夠精確。相比上述兩種方法，極限分析法可以得到邊坡極限荷載的嚴格上限解，在此基礎上進行可靠度分析，可以得到嚴格邊坡可靠度上限值，這對于分析邊坡穩(wěn)定性可能具有重要意義。在可靠度分析方法方面，常用的分析方法有一次二階矩法、JC法、MonteCarlo法等。其中，Monte Carlo法被認為是一種相對精確的方法，根據(jù)大數(shù)定律，只要抽樣次數(shù)足夠大，其精度就能足夠高。目前，眾多學者應用可靠度分析理論對邊坡穩(wěn)定性進行研究均是在線性M-C破壞準則假設下進行的。而事實上，巖土介質(zhì)服從非線性破壞準則，線性破壞準則只是一種特例：將更為符合實際的略微彎曲巖土材料強度線簡化成直線形狀，該方法雖簡單易于分析，但無法準確表述巖土強度特性。因此，考慮巖土破壞準則為非線性情況下的邊坡可靠度研究顯得十分必要。1987年，Zhang等提出了冪函數(shù)非線性破壞準則，爾后，大量學者對基于該破壞準則下的邊坡穩(wěn)定性問題進行了深入研究，非線性M-C破壞準則得到快速發(fā)展且已較為成熟。但多年來，鮮見基于非線性M-C破壞準則下邊坡可靠度的深入研究。
基于以上考慮，在非線性Mohr-Coulomb破壞準則下，結(jié)合外切線技術(shù)和強度折減技術(shù)計算邊坡安全系數(shù)上限解，進一步運用蒙特卡洛法計算邊坡可靠度。通過與基于線性M-C破壞準則的計算結(jié)果進行對比，驗證了方法的可行性。并視初始粘聚力Co、內(nèi)摩擦角arctan（Co/σt）與非線性參數(shù)m為隨機參數(shù)且服從截斷正態(tài)分布，研究非線性M-C破壞準則下各參數(shù)變異性、敏感性對均質(zhì)各向同性邊坡可靠度的影響。
1.非線性M-C破壞準則及抗剪強度參數(shù)引入方法
2.基于強度折減技術(shù)的邊坡極限分析上限法
邊坡可靠度分析中，安全系數(shù)計算式表示了邊坡安全系數(shù)與土工參數(shù)之問的關聯(lián)關系，作為功能函數(shù)的構(gòu)建主體至關重要。極限分析上限法基于虛功率原理推導，根據(jù)外力做功和內(nèi)部耗能相等原理獲得目標函數(shù)并根據(jù)能量耗散最小化原理獲得極限荷載的最小值，進一步結(jié)合強度折減技術(shù)可獲得嚴格精確的安全系數(shù)上限解。
選取符合簡單邊坡的對數(shù)螺旋面破壞機構(gòu)作為破壞模式，以通過坡趾下的對數(shù)螺旋線旋轉(zhuǎn)間斷機構(gòu)為例進行分析，如圖2。剛性塊體ABCCA繞旋轉(zhuǎn)中心O相對對數(shù)螺旋面BC以下的靜止材料作剛體旋轉(zhuǎn)，BC'是速度問斷面。ABCCA區(qū)重力做功率為外荷載功率；問斷面BC上滑面耗損率為
在已知邊坡幾何尺寸及巖土材料參數(shù)條件下，邊坡安全系數(shù)Fa可由θCh、θCD、β、φt等4個未知參數(shù)確定，且安全系數(shù)計算公式是Fs的隱函數(shù)。將Fs作為目標函數(shù)，通過非線性數(shù)學規(guī)劃方法，利用Matlab軟件，可求得邊坡最小安全系數(shù)Fs。
3.基于非線性M-C破壞準則的邊坡可靠度上限計算
3.1巖土隨機參數(shù)的選定及其分布類型
線性M-C破壞準則下巖土參數(shù)分布類型已有大量研究成果。馬建全等認為不同環(huán)境下巖土參數(shù)具有不用分布類型，同時，比較了巖土參數(shù)服從正態(tài)分布、對數(shù)正態(tài)分布等不同分布類型下的邊坡可靠度大小；張繼周等從概率分布類型的產(chǎn)生背景、所描述對象的物理意義人手，研究各分布類型對可靠度分析的影響，得出正態(tài)分布和對數(shù)正態(tài)分布較為合理的結(jié)論；其他眾多邊坡可靠度研究亦在巖土參數(shù)服從正態(tài)分布下進行。因此，筆者也選取正態(tài)分布作為巖土隨機參數(shù)的分布類型。
線性M-C破壞準則下，巖土參數(shù)c與φ服從正態(tài)分布。由式（2）可知，當m=1時，非線性M-C破壞準則指數(shù)型表達式中的參數(shù)Co與無量綱參數(shù)Co/σt的反正切值arctan（Co/σt）服從正態(tài)分布。當m≠1時，假設參數(shù) Co與內(nèi)摩擦角arctan（Co/σt）服從正態(tài)分布，并進一步考慮非線性參數(shù)的變異性，假設m亦服從正態(tài)分布。由此，選取初始粘聚力c。內(nèi)摩擦角arctan（Co/σt）與非線性參數(shù)m作為隨機參數(shù)。
桂勇等在邊坡穩(wěn)定二元體系的建立中，采用同時考慮材料指標的統(tǒng)計分布和區(qū)間分布的方法，即截斷分布；Johari等在基于畢紹普法的地震邊坡穩(wěn)定概率模型建立中，采用截斷正態(tài)分布作為隨機變量的概率分布。截斷分布可以考慮巖土參數(shù)具體變化范圍，更符合工程實際。選取截斷的正態(tài)分布作為巖土隨機參數(shù)的分布類型，隨機參數(shù)概率密度函數(shù)為
4.對比分析
目前，鮮有非線性M-C破壞準則方面的邊坡可靠度研究。為驗證方法的可行性，基于線性M-C破壞準則將邊坡可靠度上限解與Johari等、蔡寧等及呂楊的計算結(jié)果進行對比分析。參數(shù)信息如表1所示，對比結(jié)果如表2所示。
由表2可知，線性M-C破壞準則下本文結(jié)算結(jié)果與Johari等、蔡寧等及呂楊的計算結(jié)果失效概率PF最大相差5.33%。最小相差0.90%，平均相差2.455%；可靠指標口最大相差1.77%，最小相差0.57%，平均相差1.105%，表明該方法可行。
5.算例分析
選取一典型邊坡作為研究對象，如圖3所示。研究依據(jù)文獻確定巖土材料參數(shù)Co和σt獲得參數(shù)c。與內(nèi)摩擦角arctan（Co/σt）的均值，分別為90kN/m°與20°（σt為247.3 kN/m。）。并假設坡角α=60。坡高H=25m，巖土容重y=20kN/m3，非線性參數(shù)均值um=1.6。隨機參數(shù)變異系數(shù)取值情況為：1）參數(shù)Co與內(nèi)摩擦角arctan（Co/σt）沿用Johari等的c與ψ的變異系數(shù)，分別為0.20與0.10。2）非線性參數(shù)m對于不同巖土材料有一定變化范圍，因此，非線性參數(shù)m取較小變異系數(shù)符合實際工程，取值為0.05。截斷的正態(tài)分布隨機參數(shù)如表3所示。
5.1參數(shù)變異性影響
為分析非線性M-C破壞準則下隨機參數(shù)（以下簡稱參數(shù)）的變異性對邊坡可靠度的影響，進行參數(shù)變異性影響分析。當參數(shù)變異系數(shù)的增加值δ=0、COV/3、2COV/3、3COV/3、4COV/3、5COV/3、6COV/3時，邊坡失效概率Pf的變化情況如圖4所示；邊坡可靠指標β的變化情況如圖5所示。對應（δ=0、δ=COV）的安全系數(shù)統(tǒng)計結(jié)果對比如圖6～7所示；對應（δ=0、δ=COV）的可靠度相對變化如表4所示。
由圖4～7與表4可知，初始粘聚力Co的變異性對邊坡可靠度具有較大影響；內(nèi)摩擦角arctan（Co/σt）與非線性參數(shù)m的變異性對邊坡可靠度的影響較小。隨著初始粘聚力變異系數(shù)的增大，邊可靠度減小得較為明顯，且這種變化幅度隨著變異系數(shù)的增大逐漸減小；隨著內(nèi)摩擦角arctan（Co/σt）與非線性參數(shù)變異系數(shù)的增大，邊坡可靠度減小且較為緩慢。
對于該算例，巖土材料參照文獻[12]選取，其初始粘聚力為90kN/m2初始內(nèi)摩擦角arctan（Co/σt）為20°。已有研究表明，結(jié)合上限理論與外切線技術(shù)后，所獲得的瞬時粘聚力Ct隨非線性參數(shù)的增大而增大，而瞬時內(nèi)摩擦角ψt隨非線性參數(shù)的增大而減小，可知該算例中粘聚力參數(shù)對邊坡穩(wěn)定性的影響占主導作用。同時，加之在計算過程中選取了較大的變異系數(shù)（20%），使得初始粘聚力的變異性對邊坡可靠度影響較其它參數(shù)大得多。
5.2敏感性分析
為研究在非線性M-C破壞準則下隨機參數(shù)（以下簡稱參數(shù)）的變化對邊坡可靠度的影響，進行參數(shù)敏感性分析。當參數(shù)均值的增加值δ=0、std/3、2std/3、3std/3、4std/3、5std/3、6std/3時，邊坡失效概率Pf的變化情況如圖8所示；邊坡可靠指標β的變化情況如圖9所示。對應（δ=0、δ=td）的安全系數(shù)統(tǒng)計結(jié)果對比如圖10～11所示；對應（δ=0、δ=std）的可靠度相對變化如表5所示。
由圖8～11和表5可知，初始粘聚力Co的均值對邊坡可靠度有較大影響；無量綱參數(shù)arctan（Co/σt）與非線性參數(shù)m的均值對邊坡可靠度的影響相對較小。隨著初始粘聚力和無量綱參數(shù)arctan（Co/σt）均值的增大，抗剪強度參數(shù)Ct和ψt均增大，邊坡可靠度增大；隨著非線性參數(shù)均值的增大，抗剪強度參數(shù)Ct和ψt發(fā)生變化，邊坡可靠度減小。
6.結(jié)論
基于非線性M-C破壞準則，結(jié)合外切線技術(shù)和強度折減技術(shù)，采用極限分析上限法求得邊坡安全系數(shù)上限解，進一步運用蒙特卡洛法計算邊坡的可靠度（失效概率Pf與可靠指標β）。在線性M-C破壞準則下，通過與已有算例對比計算，驗證了方法的可行性。同時，基于非線性破壞準則分析參數(shù)變異性和參數(shù)敏感性對邊坡可靠度的影響，得到如下結(jié)論：
1）邊坡可靠度隨初始粘聚力Co、內(nèi)摩擦角arctan（Co/σt）和非線性參數(shù)m變異系數(shù)的增大而減小。
2）邊坡可靠度隨初始粘聚力c。和內(nèi)摩擦角arctan（Co/σt）均值的增大而增大，隨非線性參數(shù)m均值的增大而減小。
工程實際中，巖土材料參數(shù)分布形態(tài)眾多，筆者研究的基于非線性M-C破壞準則的邊坡可靠度上限分析僅采用截斷正態(tài)分布這一種，其他分布形態(tài)的研究有待深入；巖土材料參數(shù)取值范圍和變異系數(shù)差異較大，研究的基于非線性M-C破壞準則的邊坡可靠度上限分析采用了已有文獻中參數(shù)分析的取值，符合工程實際的其他參數(shù)取值范圍和變異性的研究有待深入。
