李 珂 姜夢林

(1．西安公路研究院，陜西 西安 710054; 2．機(jī)械工業(yè)勘察設(shè)計研究院，陜西 西安 710043)

安全系數(shù)是由確定的方法得到的一個定值，它沒有考慮設(shè)計參數(shù)可能發(fā)生的內(nèi)在的變異性，所以其大小并不能完全代表一個邊坡的安全程度。因此，邊坡工程的可靠度分析被引入巖土工程領(lǐng)域，經(jīng)過多年的發(fā)展，邊坡工程可靠性研究已經(jīng)有了許多研究成果，并逐步用于設(shè)計實踐［1］。目前常用的可靠度分析計算方法有蒙特卡羅法(Monte Carlo法)、一次二階矩法(FOSM法)、點(diǎn)估計法。蒙特卡羅法適用于隨機(jī)變量的概率密度分布形式已知或符合假定的情況，在目前可靠度計算中，是一種相對精確的方法。隨著計算機(jī)技術(shù)的不斷發(fā)展，蒙特卡羅法在工程中廣為流行。點(diǎn)估計法中應(yīng)用最廣的為Rosenblueth法，當(dāng)各狀態(tài)變量的概率分布為未知時，利用其均值和方差，便可以求得狀態(tài)函數(shù)(安全系數(shù)或安全儲備)的一階矩(均值)、二階矩(方差)及三、四階矩，從而求得邊坡的可靠指標(biāo)，因而應(yīng)用廣泛。蒙特卡羅模擬法是理論比較成熟、精度比較高的一種數(shù)值計算方法，由于其模型和程序簡單，且不受狀態(tài)變量分布類型以及變量間關(guān)系的限制，所以在多個狀態(tài)變量概率分布已知時，對于陜縣1號點(diǎn)，臨猗縣北景鄉(xiāng)2號點(diǎn)、臨猗縣北景鄉(xiāng)3號點(diǎn)、萬榮縣4號點(diǎn)，借助GeoStudio軟件可以較快實現(xiàn)計算結(jié)果。而在狀態(tài)變量概率分布數(shù)據(jù)缺乏時，為了保證可靠度的精度，優(yōu)先選用Rosenblueth法進(jìn)行計算。

1 蒙特卡羅法計算原理

一般情況下，我們把邊坡安全系數(shù)定義為一個狀態(tài)函數(shù):

其中，F(xiàn)為穩(wěn)定系數(shù);X1，X2，…，Xm為m個有一定分布規(guī)律、統(tǒng)計獨(dú)立的隨機(jī)變量?？梢允侨葜?、內(nèi)摩擦角、粘聚力等隨機(jī)變量。隨機(jī)地從樣本Xi中抽取n個同分布變量 x1，x2，…，xn，將每一個變量代入式(1)，便能得到一個隨機(jī)的安全系數(shù) Fi，經(jīng)過N次重復(fù)計算后達(dá)到預(yù)期的精度要求，可得到N個安全系數(shù)樣本觀測值F1，F(xiàn)2，…，F(xiàn)n，其計算結(jié)果相互獨(dú)立。因邊坡失穩(wěn)的極限狀態(tài)為F=1，所以設(shè)在N次抽樣計算過程中F≤1的次數(shù)為N，則

李萍等［2］對黃土地區(qū)公路路塹高邊坡進(jìn)行了大量研究，本文參考其調(diào)查結(jié)果，選取4個邊坡進(jìn)行可靠度分析計算。本次所研究的4個邊坡工點(diǎn)位于黃土地區(qū)的三門峽—河津亞區(qū)，其上覆土層主要為馬蘭黃土Q3，下層屬于離石黃土Q2，其中Q3黃土厚度在10 m左右，根據(jù)已有數(shù)據(jù)建立地質(zhì)模型(見表1)。

表1 黃土邊坡地質(zhì)模型

根據(jù)實測數(shù)據(jù)，結(jié)合三門峽—河津亞區(qū)的黃土特性(見表2)［2］，建立概率模型。用GeoStudio軟件中的Slope計算模塊，采用蒙特卡羅法對每一個工點(diǎn)進(jìn)行10 000次的抽樣計算，得到每一處邊坡的可靠度和失效概率(見表3)，以及安全系數(shù)的概率分布圖(見圖1)。

表2 三門峽—河津亞區(qū)黃土物理特性

表3 蒙特卡羅法對邊坡的計算結(jié)果

圖1 四個邊坡的安全系數(shù)概率分布圖

2 Rosenblueth法

Rosenblueth法是由Rosenblueth于1975年提出的一種矩估計的近似方法，其原理簡單，可快速地對邊坡可靠度進(jìn)行分析。它要求在某幾個點(diǎn)上估計功能函數(shù)的值，根據(jù)這些數(shù)據(jù)即可通過簡單的計算公式確定可靠度指標(biāo)［3］。

在缺少狀態(tài)變量x1，x2，…，xn的分布函數(shù)的情況下，不需要考慮它們變化形態(tài)，僅需在(xmin，xmax)區(qū)間內(nèi)對稱地選擇2個取值點(diǎn)，一般取均值μx的正負(fù)一個標(biāo)準(zhǔn)差σx，即:

對于n個狀態(tài)變量，可得到2n個取值點(diǎn)，全部取值點(diǎn)組合后可以得到2n個。在2n個狀態(tài)變量組合下，可根據(jù)狀態(tài)方程得到2n個狀態(tài)函數(shù)Z，即2n個安全系數(shù)。若n個狀態(tài)變量相互獨(dú)立，各個組合的出現(xiàn)概率相等，則Z的均值為:

在黃土邊坡可靠度分析中，粘聚力c和內(nèi)摩擦角φ是影響邊坡穩(wěn)定性的重要指標(biāo)，利用式(8)，式(9)對邊坡進(jìn)行可靠度計算。

當(dāng)黃土邊坡只有一種土層時，粘聚力c和內(nèi)摩擦角φ有4種組合，進(jìn)而可以得到4個不同的穩(wěn)定系數(shù)，如式(10):

此次所選用的4個工點(diǎn)為涇陽地區(qū)黃土邊坡，屬于三門峽—河津亞區(qū)［2］黃土，上部覆蓋有大約10 m厚的Q3馬蘭黃土，下部為Q2離石黃土。根據(jù)表2可知，上下土層的粘聚力c、內(nèi)摩擦角φ和標(biāo)準(zhǔn)差，可以得到16種組合(見表4)，其計算結(jié)果見表5。

表4 上下土層組合后得到的穩(wěn)定系數(shù)F

表5 Rosenblueth法對邊坡的計算結(jié)果

3 計算比較

Rosenblueth法是基于隨機(jī)變量在有限的對稱點(diǎn)處取值的點(diǎn)估計，所以計算精度沒有蒙特卡羅法方法高［4］。以蒙特卡羅法的計算結(jié)果作為參照值，將Rosenblueth法的結(jié)果代入式(11)計算其相對誤差。

其中，X為相對誤差;YM為蒙特卡羅法計算值;YR為Rosenblueth法計算值。其計算結(jié)果見表6。由表6可以看出，較蒙特卡羅法，Rosenblueth法計算得到的穩(wěn)定系數(shù)偏高，但相對誤差均不超過1%，可靠度相對誤差約為5%，標(biāo)準(zhǔn)差相對誤差約為2%，與蒙特卡羅法結(jié)果非常接近，能夠滿足工程精度要求，但是失效概率的相對誤差達(dá)到了25%，且普遍偏低。

表6 相對誤差計算值

4 結(jié)語

1)Rosenblueth法計算次數(shù)少、效率高，不需要了解各種狀態(tài)變量的概率分布，只要利用它的均值和方差，就可以求出狀態(tài)函數(shù)的一階、二階矩。在邊坡工程需要應(yīng)用數(shù)值計算的情況下可以方便地進(jìn)行可靠度分析，但其可靠性需要得到進(jìn)一步驗證。2)蒙特卡羅法原理簡單，計算量大，但借助Geoslope軟件可以快速實現(xiàn)，計算結(jié)果精度高，誤差與樣本容量有關(guān)，所以在條件允許的情況下盡可能地進(jìn)行多次計算，對于給定的誤差限，可求出所需模擬的次數(shù)。3)考慮受現(xiàn)場條件的限制，在無相關(guān)軟件進(jìn)行現(xiàn)場自動分析而計算樣本不是很多的情況下，可以優(yōu)先考慮選用Rosenblueth法進(jìn)行初步計算，蒙特卡羅模擬法可作為一種校驗方法，對其計算結(jié)果進(jìn)行校驗和對比。4)在一般情況下，土體邊坡受力情況較為復(fù)雜，尤其在土質(zhì)為非均質(zhì)土質(zhì)邊坡中，選取大量的狀態(tài)變量發(fā)現(xiàn)其變化差值很大，由此得出的極限狀態(tài)方程大多是非線性的。

蒙特卡羅模擬法不受分析條件的限制，可以忽略其極限狀態(tài)方程的線性關(guān)系，是否屬正態(tài)分布，通過大量的模擬結(jié)果對此類土質(zhì)邊坡進(jìn)行精確的可靠度分析。Rosenblueth法則不宜應(yīng)用于此類非均質(zhì)邊坡中，由于其選取的狀態(tài)變量較少，土體的復(fù)雜分布會導(dǎo)致分析結(jié)果與實際差異較大，不能對土質(zhì)邊坡進(jìn)行準(zhǔn)確地判斷。

［1］祝玉學(xué)．邊坡可靠性分析［M］．北京:冶金工業(yè)出版社，1993．

［2］李 萍，王秉綱，李同錄，等．陜西地區(qū)黃土路塹高邊坡可靠度研究［J］．中國公路學(xué)報，2009，22(6):18-25．

［3］王劉洋．Rosenblueth法在土坡穩(wěn)定可靠度分析中的應(yīng)用［J］．巖土工程技術(shù)，2005，19(2):72-74．

［4］唐亞明．基于可靠度的黃土斜坡穩(wěn)定性分析［J］．地質(zhì)通報，2008，27(8):1217-1222．