鮑俊文，胡欣宇，邢明源，趙少飛

(華北科技學(xué)院 建筑工程學(xué)院，北京 東燕郊 065201)

0 引言

自20世紀(jì)中期巖土工程領(lǐng)域開展可靠度研究以來，巖土工程可靠度研究成果在工程實(shí)踐中產(chǎn)生了重大影響。目前，蒙特卡羅(Monte Carlo)方法是可靠度分析中的主要手段之一，還是檢驗(yàn)其它可靠度計(jì)算方法精度的唯一途徑。由于巖土工程存在很多不確定性因素，該方法也大量應(yīng)用到巖土工程可靠度分析。文獻(xiàn)[1]中詳細(xì)的介紹了該方法的原理。楊明等[2]應(yīng)用蒙特卡羅法，分析了均質(zhì)土坡的穩(wěn)定可靠度。吳振君等[3]應(yīng)用該方法對(duì)邊坡工程將進(jìn)行了可靠度分析，討論了拉丁超立方抽樣(LHS)方法和蒙特卡羅傳統(tǒng)抽樣計(jì)算可靠度指標(biāo)的各種方法，得到LHS法較傳統(tǒng)抽樣方法效率上有很大改善。針對(duì)蒙特卡羅傳統(tǒng)抽樣方法計(jì)算量大的問題，伍國軍等[4]基于拉丁超立方抽樣方法編制了有限元可靠度程序。由于蒙特卡羅方法的結(jié)果精度主要取決于抽出樣本的質(zhì)量,符蓉[5]對(duì)比了拉丁超立方抽樣方法與傳統(tǒng)抽樣方法，得出拉丁超立方抽樣方法計(jì)算效率明顯優(yōu)于傳統(tǒng)抽樣方法。曹偉等[6]應(yīng)用LHS方法進(jìn)行了邊坡可靠度的有限元分析，認(rèn)為LHS方法相較于傳統(tǒng)抽樣方法更加高效，且計(jì)算結(jié)果較為精確。Nguyen和Chowdhury[8]提出一種應(yīng)用于巖土工程的風(fēng)險(xiǎn)分析方法，得到了良好的結(jié)果。

為了對(duì)比蒙特卡羅方法中傳統(tǒng)抽樣方法和拉丁超立方抽樣結(jié)果的差別，采用Matlab軟件進(jìn)行抽樣編程，分析了抽樣方法對(duì)條形基礎(chǔ)極限承載力可靠度結(jié)果的影響。

1 蒙特卡羅方法的原理

蒙特卡羅方法需要首先產(chǎn)生隨機(jī)變量樣本，然后計(jì)算得到函數(shù)的樣本，進(jìn)而統(tǒng)計(jì)失效區(qū)域的樣本數(shù)來估計(jì)失效概率，該方法的優(yōu)點(diǎn)是概念清晰，易于使用。這種方法的主要缺點(diǎn)是用于準(zhǔn)確估計(jì)失效概率的樣本數(shù)應(yīng)足夠大。這必然導(dǎo)致過大的計(jì)算量[1]。

可以寫成

(1)

(2)

其中，xi是sx(x)的一個(gè)樣本，N為樣本總數(shù)。

蒙特卡羅法實(shí)施即為上述計(jì)算積分過程。對(duì)于工程失效概率Pf

(3)

(4)

因?yàn)镮(x)服從二項(xiàng)分布且參數(shù)為Pf，其方差為NPf(1-Pf)。Pf是I(x)的均值，Pf的方差可以估計(jì)為

(5)

(6)

2 傳統(tǒng)舍選法和拉丁超立方法原理

蒙特卡羅傳統(tǒng)抽樣方法有逆變換法和舍選法，由于逆變換法有一定的條件限制，適用于隨機(jī)變量X具有累積分布函數(shù)表達(dá)式情況。而舍選法僅需要隨機(jī)變量的概率密度函數(shù)表達(dá)式，本文選用舍選法與拉丁超立方抽樣方法進(jìn)行對(duì)比。

2.1 舍選法

假如hX(x)和sY(y)為兩個(gè)概率密度函數(shù)。先從sY(y)抽取樣本，然后按下面步驟對(duì)生成的樣本進(jìn)行選擇，使得選出的樣本分從hX(x)的分布。具體抽樣步驟如下：

10.秋季仔豬容易發(fā)生豬瘟、霉菌中毒、圓環(huán)病毒病、副傷寒、豬流感、喘氣病、藍(lán)耳病及球蟲病，故建議在此期間針對(duì)性的不間斷投服“金泰妙（45%泰妙菌素可溶性粉）”、 “金奇氟(20%替米考星)”、伊愛爾（0.6%伊維菌素預(yù)混劑）、“氟力佳（20%氟苯尼考可溶性粉）”，連用5～7 d。

(1) 在sY(y)中抽取一個(gè)樣本為y*。

(2) 在均勻分布[0, 1]區(qū)間隨機(jī)抽一個(gè)樣本u，若u

(3) 再重復(fù)(1)步驟就可以得到需要的樣本數(shù)量。

2.2 拉丁超立方法

在給定抽樣函數(shù)情況下, 拉丁超立方抽樣方法更能代表抽樣函數(shù)。對(duì)于n個(gè)相互獨(dú)立的隨機(jī)變量X1，X，…，Xn，抽樣的具體實(shí)施步驟如下:

(1) 每一個(gè)隨機(jī)變量Xk被分成M個(gè)區(qū)域，使得xk1

(2) 假定U是服從[0,1]區(qū)間上的均勻分布，從中抽取一個(gè)樣本u，在每個(gè)區(qū)域中抽取Xk在第m個(gè)區(qū)間段抽取的樣本為

(7)

在整個(gè)區(qū)間，共抽取隨機(jī)變量Xk的M個(gè)隨機(jī)變量。

3 抽樣方法的對(duì)比

3.1 條形基礎(chǔ)極限承載力計(jì)算

巖土工程中，條形基礎(chǔ)的極限承載力常用Hansen公式計(jì)算，極限承載力qu為

qu=0.5γBNγ+cNc+γmdNq

(8)

圖1 承受均布荷載的淺基礎(chǔ)

以某一條形基礎(chǔ)為例[1]，基礎(chǔ)寬度為2m，埋深為0.5m,土的重度均取20kN/m3，黏聚力和內(nèi)摩擦角為隨機(jī)變量，服從正態(tài)二元分布，黏聚力的均值和標(biāo)準(zhǔn)差分別為20kPa和5kPa，內(nèi)摩擦角的均值和標(biāo)準(zhǔn)差分別為30°和6°，兩者的相關(guān)系數(shù)為-0.5。根據(jù)Hansen公式計(jì)算的極限承載力qu的均值為1148.5kPa，基礎(chǔ)承受的荷載q為1000kN/m。

3.2 舍選法計(jì)算

用傳統(tǒng)舍選法，隨機(jī)抽樣數(shù)量分別為3,000、5,000、50,000、200,000和500,000，計(jì)算該淺基礎(chǔ)相應(yīng)的失效概率和變異系數(shù)，結(jié)果如表1所示。

表1 傳統(tǒng)舍選法結(jié)果

3.3 拉丁超立方法計(jì)算

對(duì)于上述算例，假設(shè)粘聚力和內(nèi)摩擦角相互獨(dú)立，在拉丁超立方抽樣中M取5，為了對(duì)比兩種抽樣方法，抽樣數(shù)量一樣，分別為3,000、5,000、50,000、200,000和500,000，失效概率和變異系數(shù)如表2所示。

表2 拉丁抽樣法結(jié)果

3.4 結(jié)果對(duì)比分析

兩種方法的結(jié)果如圖2所示。由圖2可知，當(dāng)抽樣個(gè)數(shù)達(dá)到50,000時(shí)，拉丁超立方抽樣方法所計(jì)算的失效概率已經(jīng)趨于穩(wěn)定，而使用傳統(tǒng)舍選法抽樣需要達(dá)到200,000時(shí)，失效概率才趨于穩(wěn)定。采用蒙特卡羅進(jìn)行可靠度分析時(shí)，要得到失效概率穩(wěn)定值，傳統(tǒng)舍選法抽樣數(shù)量為拉丁超立方抽樣數(shù)量的4倍，表明使用拉丁抽樣方法可以大大減少抽樣數(shù)量。

兩類抽樣方法所得到失效概率的變異系數(shù)差別顯著，如圖3所示。傳統(tǒng)舍選法變異系數(shù)隨著抽樣數(shù)量的增加，減小緩慢，而拉丁超立方抽樣法得到變異系數(shù)減小明顯。隨著抽樣數(shù)量的增加，兩類抽樣方法所得變異系數(shù)的比值由1.0，增加到2.0以上。對(duì)于達(dá)到同樣的變異系數(shù)，需要的樣本數(shù)量差別懸殊，例如當(dāng)抽樣數(shù)量為500,000時(shí)，變異系數(shù)為8.0%，而拉丁超立方法的樣本數(shù)量50,000時(shí)，達(dá)到了相同的變異系數(shù)。這表明，達(dá)到相同的變異系數(shù)，拉丁超立方抽樣數(shù)量可以減小到樣本數(shù)量的1/10。

圖2 失效概率隨抽樣數(shù)量的變化

圖3 變異系數(shù)隨樣本數(shù)量的變化

除了達(dá)到穩(wěn)定失效概率時(shí)，兩種抽樣數(shù)量差別懸殊外，在配置為Intel酷睿雙核i3-6100臺(tái)式機(jī)上,計(jì)算時(shí)間有很大的差異。以抽樣數(shù)量為5,000 為例，計(jì)算過程表明，傳統(tǒng)舍選法所用的時(shí)間約為65.17 s，而拉丁超立方抽樣方法計(jì)算所用的時(shí)間為2 s，這表明在相同抽樣數(shù)量的情況下，拉丁超立方抽樣方法計(jì)算所用的時(shí)間要遠(yuǎn)小于傳統(tǒng)舍選法。

傳統(tǒng)舍選法計(jì)算的失效概率為0.35%，而拉丁超立方抽樣方法為0.32%，兩種方法的結(jié)果差異，是由于舍選法中土的粘聚力和內(nèi)摩擦角被認(rèn)為是二元正態(tài)分布的隨機(jī)變量，而拉丁超立方抽樣法中，假設(shè)粘聚力和內(nèi)摩擦角相互獨(dú)立，這是導(dǎo)致兩者結(jié)果略有差異的原因。

4 結(jié)論

(1) 采用蒙特卡羅法計(jì)算工程的失效概率時(shí)，拉丁超立方法得到失效概率所需要抽樣數(shù)量顯著小于傳統(tǒng)舍選法所需的抽樣數(shù)量，需要的樣本數(shù)量約為舍選法數(shù)量的25%。

(2) 拉丁超立方法通過增加樣本數(shù)量能夠更有效減小失效概率的變異系數(shù)，相同抽樣數(shù)量時(shí)其失效概率變異系數(shù)小于舍選法所得結(jié)果的一半。