基于橢球模型的巖坡穩(wěn)定非概率可靠評價研究

(南京水利科學(xué)研究院，江蘇 南京 210098)

·建筑與土木工程·

基于橢球模型的巖坡穩(wěn)定非概率可靠評價研究

韓孝峰，陳海軍

(南京水利科學(xué)研究院，江蘇 南京 210098)

在考慮影響因素的不確定性條件下，快速有效地評價巖坡的穩(wěn)定狀態(tài)成為目前研究的熱點。傳統(tǒng)概率可靠方法評價要求數(shù)據(jù)充足且服從某種概率分布，工程中難以實測大量數(shù)據(jù)支持，從而阻礙了該方法的推廣應(yīng)用。本文針對以上不足，結(jié)合平面滑動型巖坡在只需確定變化范圍的前提下應(yīng)用橢球模型表達不確定影響因素的變異性質(zhì)，采用相應(yīng)的非概率可靠指標(biāo)評價了巖石邊坡的穩(wěn)定狀態(tài)，與蒙特卡洛模擬法計算結(jié)果的對比顯示兩種評價結(jié)果一致?；跈E球模型的非概率可靠評價方法所需數(shù)據(jù)量小，計算簡便快捷，在巖坡穩(wěn)定問題的評價中具有一定的適用性。

橢球模型；非概率可靠；巖石邊坡；穩(wěn)定性評價；蒙特卡洛模擬

隨著我國經(jīng)濟建設(shè)的發(fā)展，各大工程領(lǐng)域中都會涉及到巖石邊坡穩(wěn)定問題[1]。巖坡的安全穩(wěn)定對工程具有不小影響，快速準(zhǔn)確地評價巖坡的穩(wěn)定性很重要。實際工程中，巖坡穩(wěn)定影響因素的取值往往不是一個固定值。在分析巖坡穩(wěn)定問題的定性分析、定量分析和不確定可靠分析這3大類方法中[2]，不確定可靠分析方法能夠考慮影響因素的隨機變異特點[3]；但目前運用較多的概率不確定可靠性分析方法對巖坡穩(wěn)定進行分析評價時，需要明確各個影響參數(shù)的統(tǒng)計概率分布規(guī)律[4]。這在實際工況中是無法實現(xiàn)的，若采用小樣本的概率分布模型計算，則會出現(xiàn)較大偏差，影響評價結(jié)果。針對以上存在的問題，本文將橢球模型與可靠性分析原理相結(jié)合，建立一種在數(shù)據(jù)不足、樣本數(shù)目很小的情況下也能評價巖石邊坡穩(wěn)定性的方法，該評價方法有別于傳統(tǒng)概率分析方法，可劃入非概率可靠分析方法范疇中。由于該方法不用建立影響因素的概率分布，使計算量得到縮減，可十分快速地得到評價結(jié)果。結(jié)合經(jīng)典算例評價，并與蒙特卡洛模擬結(jié)果進行對比，顯示評價結(jié)果十分準(zhǔn)確。

1 不確定參數(shù)的橢球表達

橢球模型能夠準(zhǔn)確地表達出不確定參數(shù)的變化特點。在巖坡穩(wěn)定分析問題中，不確定參數(shù)指的就是影響邊坡穩(wěn)定的因素。橢球模型用來表達巖坡穩(wěn)定影響因素變化性質(zhì)時，可以看作是每個影響參數(shù)不同取值的組合。從幾何上看，橢球模型具有獨特的形狀和尺度，形狀反映了該模型對不確定參數(shù)的已知程度，尺度反映了不確定參數(shù)的偏離程度[5]。

假設(shè)存在n個不確定的參數(shù)，用X∈Rn表示不確定參數(shù)可能取值組成的向量。橢球模型用一個超橢球體來界定n個不確定參數(shù)，可用下式表達：

(1)

2 基于橢球模型的非概率可靠性分析方法的建立

2.1可靠性分析理論

結(jié)構(gòu)分析的核心思想是使結(jié)構(gòu)在使用期限內(nèi)能夠滿足安全性、適用性和耐久性等3個方面的要求[6]，結(jié)構(gòu)的可靠指的是滿足上述3項要求。結(jié)構(gòu)的可靠性程度一般用可靠度來定量描述。結(jié)構(gòu)的可靠性往往受到各種外部荷載、內(nèi)部介質(zhì)強度以及結(jié)構(gòu)自身幾何尺寸等因素的影響，這些因素往往是在一定范圍內(nèi)變化且其變化規(guī)律具有不確定性?？煽啃苑治鲋蟹Q這些影響結(jié)構(gòu)可靠性的因素為基本變量[7]，即對應(yīng)巖坡穩(wěn)定問題中提到的各個影響因素。設(shè)Xi,i=1,2,…,n為結(jié)構(gòu)的n個基本變量，則結(jié)構(gòu)安全狀態(tài)的功能函數(shù)可用基本變量的復(fù)合函數(shù)來表示：

(2)

上式是結(jié)構(gòu)基本變量的復(fù)合函數(shù)，這一函數(shù)表示了結(jié)構(gòu)可靠性的狀態(tài)，隨著基本變量組合的不同，其表達的可靠性狀態(tài)也不一樣。結(jié)構(gòu)系統(tǒng)是否可靠取決于結(jié)構(gòu)的可靠性狀態(tài)，結(jié)構(gòu)的可靠性狀態(tài)是用其功能函數(shù)的大小來表示的：

2.2基于橢球模型的非概率可靠指標(biāo)

將橢球模型表達的不確定參數(shù)應(yīng)用于可靠分析理論中可獲得如下評判方法。在標(biāo)準(zhǔn)化平面中，橢球模型表示的是一個平面橢球形，如圖1所示。

該坐標(biāo)系中，基本變量的邊界可用‖Q‖2=1來表示，即圖中的單位圓。基本變量的組合點都位于單位圓區(qū)域以內(nèi)，極限狀態(tài)曲線將平面劃分為安全域和失效域。

圖1 橢球模型下的非概率可靠性度量

(3)

根據(jù)以上定義，基于橢球模型的非概率可靠性指標(biāo)取值可以為任意實數(shù)。結(jié)合可靠性分析理論，當(dāng)η≥1時，結(jié)構(gòu)的基本變量取值都位于安全域內(nèi)，此時結(jié)構(gòu)完全可靠。當(dāng)-1<η<1時，基本變量的組合點有部分落在安全域內(nèi)，也有部分落在失效域內(nèi)，此時結(jié)構(gòu)存在可靠風(fēng)險；但在該區(qū)間內(nèi)，結(jié)構(gòu)的可靠性隨著η的增大而增大。當(dāng)η<-1時，基本變量的組合點全部落在失效域內(nèi)，此時的結(jié)構(gòu)完全不可靠，對于巖坡穩(wěn)定問題來說指的就是發(fā)生滑坡。

2.3巖坡穩(wěn)定問題基于橢球模型的非概率可靠分析步驟

對巖坡穩(wěn)定問題進行非概率可靠性分析的具體步驟如圖2所示。

圖2 巖坡穩(wěn)定問題基于橢球模型的非概率可靠性分析步驟圖

3 應(yīng)用算例

3.1計算模型的選取

本文選取經(jīng)典的二維巖石邊坡計算模型[8]，具有一定的代表性。該巖坡模型中存在一組軟弱結(jié)構(gòu)面，假設(shè)巖坡發(fā)生破壞時的滑動面與該結(jié)構(gòu)面重合，發(fā)生平面型滑動破壞，如圖3所示。該巖坡模型一旦發(fā)生破壞，上部滑體沿結(jié)構(gòu)面下滑，坡頂后緣受拉產(chǎn)生張拉裂縫，裂縫會充水。從圖中可以看出該計算模型存在著以下幾何特征：

(4)

圖3 巖石邊坡示意圖

以上一些參數(shù)取值如表1所示，邊坡寬度取單位寬度[8]。

表1 巖石邊坡參數(shù)取值表

本文在進行非概率可靠評價的基礎(chǔ)上，同時采用了蒙特卡洛模擬對該巖坡的穩(wěn)定性能進行計算，旨在對比2種方法的結(jié)果是否一致。

為了更清晰地對比2種方法的評價效果，計算過程中分別考慮了有水平地震力和無水平地震力2種情況下，該巖石邊坡的穩(wěn)定狀態(tài)。

3.2巖坡平面滑動問題的橢球模型非概率可靠性分析

3.2.1 巖石邊坡穩(wěn)定性不確定參數(shù)的橢球模型

該巖石邊坡的2個不確定參數(shù)選取為潛在滑動面的黏聚力及內(nèi)摩擦角，它們的取值均服從標(biāo)準(zhǔn)正態(tài)分布，如表2所示[8]。

表2 巖石邊坡不確定參數(shù)取值分布

值得一提的是，上述2個參數(shù)的標(biāo)準(zhǔn)正態(tài)分布是在進行蒙特卡洛模擬時才需用到，而采用橢球模型表達時只需知道參數(shù)的變化范圍。這2個基本變量可用橢球模型表示為

(5)

3.2.2 功能函數(shù)的確定

巖坡沿潛在滑動面的抗滑力大于下滑力時，巖坡能保持穩(wěn)定狀態(tài)。巖坡的功能函數(shù)可以表達為

Z=R-S。

(6)

其中：R表示巖石邊坡的抗滑力；S表示巖石邊坡的下滑力。

1)不考慮水平地震作用。

對圖3的巖石邊坡模型進行受力分析，下滑力S的表達式為

S=Wsinθp+Vcosθp。

(7)

抗滑力R的表達式為

(8)

整理可得該坡穩(wěn)定性能所對應(yīng)的功能函數(shù)為

Z=12 711tanφ+70c-11 962。

(9)

2)考慮水平地震作用。

地震作用簡化為施加在物體上的一個水平慣性力。設(shè)水平地震加速度的取值在0～0.4之間。

相應(yīng)功能函數(shù)可表達為

(10)

3.2.3 非概率可靠性指標(biāo)的計算

1)不考慮水平地震作用。

巖坡穩(wěn)定功能函數(shù)可標(biāo)準(zhǔn)化為

(11)

采用直接迭代法求得橢球模型非概率可靠性指標(biāo)為1.745 7，迭代曲線如圖4所示。該指標(biāo)值表明巖石邊坡穩(wěn)定狀態(tài)可靠。從曲線看出僅用了2步迭代就得到了接近精確解的值。

圖4無水平地震作用下巖坡穩(wěn)定非概率可靠指標(biāo)迭代歷史

2)考慮水平地震作用。

考慮地震作用產(chǎn)生的水平力時，巖坡穩(wěn)定可靠性對應(yīng)的功能函數(shù)可標(biāo)準(zhǔn)化為

323qαqtanφ+2 319。

(12)

采用迭代法可快速計算出非概率可靠性指標(biāo)為0.439 2，過程如表3所示。此時巖坡穩(wěn)定不可靠。

表3 水平地震作用下巖坡穩(wěn)定非概率可靠指標(biāo)迭代歷史

通過迭代計算出的可靠性指標(biāo)，不難看出地震作用對巖坡穩(wěn)定可靠性的影響極其明顯。在水平地震作用下，該巖坡的穩(wěn)定狀態(tài)不可靠，存在失穩(wěn)破壞的風(fēng)險。

3.3巖坡平面滑動問題的Monte Carlo模擬法

蒙特卡洛模擬法能達到的精度與樣本數(shù)量存在必然的聯(lián)系，是目前概率可靠分析中較精確的方法。采用該方法的前提是要掌握有關(guān)參數(shù)的大量實測值以確定該參數(shù)的概率分布，否則蒙特卡洛模擬無法進行。

3.3.1 不考慮水平地震作用

不考慮水平地震作用時，蒙特卡洛模擬中的隨機變量為軟弱結(jié)構(gòu)面的黏聚力和內(nèi)摩擦角，它們的均值和方差見表2。

圖5為巖坡安全系數(shù)的累積密度分布圖，圖中與橫坐標(biāo)安全系數(shù)為1時所對應(yīng)的縱坐標(biāo)值大約為0.02，表明該巖石邊坡的失效概率不到2%。

圖5 無地震作用下巖坡安全系數(shù)累積概率密度分布圖

3.3.2 考慮水平地震作用

考慮水平地震作用時，隨機變量增加一個水平地震加速度。圖6為巖坡安全系數(shù)的累積密度分布圖，圖中與橫坐標(biāo)安全系數(shù)為1時所對應(yīng)的縱坐標(biāo)值大約為0.31左右，表明該巖坡的失效概率為31%左右。

蒙特卡洛模擬計算的結(jié)果表明該巖石邊坡在不受水平地震作用時的穩(wěn)定狀態(tài)很好，失效概率不到2%，在受水平地震作用時其穩(wěn)定失效概率會增至31%左右。

圖6 地震作用下巖坡安全系數(shù)累積概率密度分布圖

3.4結(jié)果討論

通過對2種方法的計算結(jié)果進行對比可以看出：

1)對于本文選取的巖石邊坡計算模型，在無水平地震作用時，基于橢球模型的非概率可靠性分析方法結(jié)果為該巖石邊坡穩(wěn)定狀態(tài)可靠。通過蒙特卡洛模擬計算得出該巖石邊坡的失效概率不到2%，該巖坡幾乎不會破壞。在有水平地震作用下，基于橢球模型非概率可靠性分析結(jié)果表明巖石邊坡穩(wěn)定狀態(tài)不可靠。此時采用蒙特卡洛方法計算得出巖石邊坡的失效概率高達31%，該巖坡存在很大的破壞風(fēng)險。

2)考慮水平地震作用和未考慮水平地震作用這2種情況下，基于橢球模型的非概率可靠性分析結(jié)果和蒙特卡洛模擬法計算結(jié)果對巖坡穩(wěn)定狀態(tài)的評價保持一致。

3)蒙特卡洛模擬法計算時需要隨機變量的概率分布，且計算耗時較長，而基于橢球模型的非概率可靠性分析方法只需明確隨機變量的變化范圍且迭代計算過程較快。

4 結(jié)論

本文在總結(jié)國內(nèi)外有關(guān)巖坡穩(wěn)定分析、橢球模型理論及可靠性分析研究現(xiàn)狀[2-10]的基礎(chǔ)上，開展了將基于橢球模型的非概率可靠性分析方法應(yīng)用于巖坡穩(wěn)定問題評價的研究工作。結(jié)合經(jīng)典的巖石邊坡模型，建立不確定參數(shù)的橢球模型表達，采用相應(yīng)的非概率指標(biāo)計算方法，利用計算結(jié)果評價分析了巖石邊坡的穩(wěn)定狀態(tài)，并通過與現(xiàn)已成熟的蒙特卡洛模擬法結(jié)果進行對比，驗證了該評價方法在巖坡穩(wěn)定問題中的適用性。得到以下結(jié)論：

1)基于橢球模型的非概率可靠性分析方法評價巖坡穩(wěn)定問題時，僅需知道不確定參數(shù)的變化范圍，避免了傳統(tǒng)概率可靠性分析時必須確定參數(shù)數(shù)學(xué)概率分布的限制條件，使計算簡便快捷。

2)在對樣本數(shù)據(jù)量要求不高的前提下，本文方法計算出的可靠性指標(biāo)值與蒙特卡洛模擬法利用不確定參數(shù)的概率分布計算出的失效概率對評價平面滑動型巖坡穩(wěn)定狀態(tài)的意義保持一致，具備對巖石邊坡穩(wěn)定狀態(tài)的評價功能。

[1]張俊云.巖石邊坡植被護坡系統(tǒng)的水分平衡及控制[J].巖石力學(xué)與工程學(xué)報 ,2013(9):1729-1735.

[2]劉章軍,陳飛,周宜紅，等.巖質(zhì)路塹深邊坡穩(wěn)定性評價的模糊概率方法[J].巖土力學(xué),2008,29(增刊1):368-372.

[3]劉曉,唐輝明,熊承仁，等.邊坡動力可靠性分析方法的模式、問題與發(fā)展趨勢[J].巖土力學(xué),2013(5):1217-1234.

[4]Johari A,Javadi A A. Reliability Assessment of Infinite Slope Stability Using the Jointly Distributed Random Variables Method[J]. Scientia Iranica,2012,19(3):423-429.

[5]羅陽軍,亢戰(zhàn),Alex Li,等.基于凸模型的結(jié)構(gòu)非概率可靠性指標(biāo)及其求解方法研究[J].固體力學(xué)學(xué)報,2011,32(6):646-654.

[6]姜潮,劉麗新,龍湘云,等.一種概率-區(qū)間混合結(jié)構(gòu)可靠性的高效計算方法[J].計算力學(xué)學(xué)報 ,2013(5):605-609.

[7]趙凱.土質(zhì)邊坡穩(wěn)定性的可靠度研究[D].西安：西安理工大學(xué),2013.

[8]李典慶,周創(chuàng)兵,陳益峰，等.邊坡可靠度分析的隨機響應(yīng)面法及程序?qū)崿F(xiàn)[J].巖石力學(xué)與工程學(xué)報,2010,29(8):1513-1523.

[9]陳帆，劉立，陸?？?，等.巖質(zhì)高陡邊坡穩(wěn)定性分析[J].西華大學(xué)學(xué)報：自然科學(xué)版，2010，29(3)：54-55.

[10]王玉平，曾志強，潘樹林.邊坡穩(wěn)定性分析方法綜述[J].西華大學(xué)學(xué)報：自然科學(xué)版，2012，31(2)：101-105.

(編校：葉超)

Non-probabilisticReliabilityAnalysisBasedontheEllipsoidalModelofRockSlopeStability

HAN Xiao-feng， CHEN Hai-jun

(NanjingHydraulicResearchInstitute，Nanjing210098China)

For evaluating the slopes’ stability accurately and quickly, non-probabilistic reliability based on the ellipsoidal model is used to overcome the existing disadvantages. Firstly, uncertain parameters are considered via the ellipsoidal model of the rock slope. And then the corresponding methods are used to calculate the non-probabilistic reliability index to evaluate the slope’s stability. Finally, the Monte-Carlo simulation were carried out to evaluate the rock slope’s stability as a comparison. The result shows that the non-probabilistic reliability analysis is more efficient and applicable.

ellipsoidal model；non-probabilistic reliability；rock slope；stability analysis；Monte-Carlo simulation

2014-11-10

中央級公益性科研院所基本科研業(yè)務(wù)費重點基金項目(Y313011)；江蘇省科技支撐計劃項目(BE2011373)

韓孝峰(1989—)，男，工程師，碩士，主要研究方向為結(jié)構(gòu)體系可靠性分析、軟土地基處理等。

*通信作者：陳海軍(1973—)，男，教授級高級工程師，注冊咨詢工程師，博士，主要研究方向為工程地質(zhì)及巖土工程。E-mail：hjchen@nhri.cn

TB114.3;TP18

：A

：1673-159X(2015)06-0096-05

10.3969/j.issn.1673-159X.2015.06.020