一種與蒙特卡羅法及JC法相結(jié)合的響應(yīng)面法

丘晉文

目前,結(jié)構(gòu)可靠性分析中大多數(shù)方法如數(shù)值積分法和一次二階矩法及其改進方法等針對功能函數(shù)都是明確表達的,而實際工程中,由于結(jié)構(gòu)本身構(gòu)造復(fù)雜,作用形式多樣,往往不能給出功能函數(shù)的明確表達式,若直接應(yīng)用上述方法進行結(jié)構(gòu)可靠度分析就會遇見困難。響應(yīng)面方法就是用來處理這種結(jié)構(gòu)功能函數(shù)不能明確表達的一種有效方法。運用響應(yīng)面法進行結(jié)構(gòu)可靠度分析時,計算可靠指標(biāo)有幾種方法:1)結(jié)合JC法。2)結(jié)合幾何法。3)結(jié)合蒙特卡羅法直接進行數(shù)值模擬。本文簡要介紹響應(yīng)面法、蒙特卡羅法和JC法,然后提出了一種與蒙特卡羅法和JC法相結(jié)合的響應(yīng)面技術(shù),并結(jié)合數(shù)值算例說明具體應(yīng)用。

1 響應(yīng)面法

響應(yīng)面法(Response Surface Method)是處理結(jié)構(gòu)功能函數(shù)不能明確表達的一種有效方法。該方法采用有限的試驗通過回歸擬合解析表達式ˉZ=ˉg(X)代替真實功能函數(shù)曲面 Z=g(X),通過響應(yīng)面法,可將功能函數(shù)近似表示為隨機變量的顯式,再結(jié)合其他計算可靠度的方法進行結(jié)構(gòu)可靠度計算。對n個隨機變量X1,…,Xn的情況,通常采用不含交叉項的二次多項式形式:

其中,a,bi,ci均為表達式的待定系數(shù)。從響應(yīng)面函數(shù)表達式可以看出,如果考慮n個隨機變量,則有m=2n+1個待定系數(shù)。要得到待定系數(shù),需要取一定的狀態(tài)進行試驗或者有限元分析。擬合的響應(yīng)面是否為真實曲面的較好近似,目前比較好的方法是在驗算點附近展開的響應(yīng)面比較接近真實的曲面。

2 蒙特卡羅法

蒙特卡羅法(Monte-Carlo Method),又稱統(tǒng)計試驗法或隨機抽樣技巧方法。此方法是一種數(shù)值模擬方法,其過程如下:

設(shè)結(jié)構(gòu)功能函數(shù)為Z=g(X1,…,Xn),Xi(i=1,…,n)為任意分布的隨機變量,對Xi(i=1,…,n)進行 N次隨機抽樣,得 N組(j=1,…,N),將第 j組(j=1,…,N)的(i=1,…,n)的值代入結(jié)構(gòu)功能函數(shù)式,得到 N個Zj值(j=1,…,N),設(shè)在N個中存在0,則結(jié)構(gòu)構(gòu)件的失效概率為

3 JC法

JC法是解決結(jié)構(gòu)構(gòu)件功能函數(shù)由獨立任意分布基本隨機變量構(gòu)成時結(jié)構(gòu)可靠指標(biāo)計算問題。JC法的基本思路是對非正態(tài)基本隨機變量作當(dāng)量正態(tài)化處理,將其轉(zhuǎn)換為等效正態(tài)隨機變量,然后用改進的一次二階矩法求結(jié)構(gòu)可靠指標(biāo),其計算步驟如下:6)以新的,重復(fù)步驟 2)～5),直到前后兩次算出的β值之差小于容許誤差。

其中,φ(),Φ()分別為標(biāo)準(zhǔn)正態(tài)分布密度函數(shù)及分布函數(shù);fXi(),FXi()分別為非正態(tài)變量Xi的分布密度函數(shù)及分布函數(shù);μi′,σi′分別為等效正態(tài)隨機變量Xi′的均值和均方差。若結(jié)構(gòu)功能函數(shù)由相關(guān)基本隨機變量構(gòu)成,要將相關(guān)變量轉(zhuǎn)換為獨立變量,然后再利用JC法計算結(jié)構(gòu)可靠指標(biāo)。

4 與蒙特卡羅法和JC法相結(jié)合的響應(yīng)面法

1)對 n個隨機變量X1,…,Xn的情況,通常采用不含交叉項的二次多項式形式:

其中,a,bi,ci分別為表達式的待定和修正系數(shù)。假定迭代點),初次計算一般取均值點。

2)利用試驗(復(fù)雜結(jié)構(gòu)的可靠度分析一般采用數(shù)值模擬試驗,如有限元計算)計算功能函數(shù)Z=g(,…,,…)以及Z=g(,…±fσi,…,)得到(2n+1)個點估計值,其中系數(shù) f在第一輪估計中取 3(根據(jù)工程中的 3σ原則),在以后的迭代計算中取1,σi為隨機變量Xi的均方差。

3)由于表達式只有(2n+1)個待定系數(shù),利用2)求得(2n+1)個函數(shù)值解出待定系數(shù)a,bi,ci,從而確定結(jié)構(gòu)的功能函數(shù)。

4)利用蒙特卡羅法求解可靠度指標(biāo)β(k),并利用此可靠度指標(biāo)β(k)采用JC法計算出此極限狀態(tài)方程下的驗算點X*(k),其中上標(biāo)k表示第k次計算。

判斷收斂條件:｜β(k)-β(k+1)｜＜ε(ε為收斂精度),滿足則停止迭代,否則用插值法得到:

這種與蒙特卡羅法和JC法相結(jié)合的響應(yīng)面法,具有響應(yīng)面法可解決功能函數(shù)不能明確表達的優(yōu)點和蒙特卡羅法模擬相對精確的優(yōu)點以及JC法計算驗算點的優(yōu)點,具有方法上的正確性和可行性、較高的計算效率和精度。下面結(jié)合算例說明其應(yīng)用。

5 算例

圖1為門式平面框架,各單元的彈性模量為 E=2.0×106kN/m2,截面慣性矩與截面積的關(guān)系為 Ii=(i=1,2)。以節(jié)點3的水平位移作為需控制的結(jié)構(gòu)最大變形,可建立極限狀態(tài)方程:

其中,[u]為允許的最大水平位移,此處可取[u]=0.01 m;u3與隨機變量A1,A2以及P的關(guān)系不能用明確表達式表達。

通過自編程序用表1中的幾種方法進行計算,應(yīng)用本文的方法共迭代8次,其中直接使用筆者編制的確定性平面框架結(jié)構(gòu)分析程序共進行了56次運算,計算結(jié)果及與其他方法對比見表1。

同時筆者研究計算表明,可靠指標(biāo)越大,失效概率越小,本文方法就越接近精確值。

表1 計算結(jié)果

6 結(jié)語

對于工程中常見的功能函數(shù)不能明確表達的情況,采用響應(yīng)面方法結(jié)合蒙特卡羅法和JC法,本文提出了一種迭代格式的響應(yīng)面法,通過數(shù)值計算表明本文方法正確性和可行性,有一定的廣泛性和較好的實用性,也可為研究人員和工程師提供一定的啟示。

[1] 趙國藩.結(jié)構(gòu)可靠度理論[M].北京:中國建筑工業(yè)出版社,2000.

[2] 張建仁.結(jié)構(gòu)可靠度理論及其在橋梁工程中的應(yīng)用[M].北京:人民交通出版社,2003.

[3] 朱殿芳,陳健康,郭志學(xué).結(jié)構(gòu)可靠度分析方法綜述[J].中國農(nóng)村水利水電,2002(8):47-49.

[4] 陳鐵冰,孫 策.響應(yīng)面法以及其在桁架橋可靠度分析中的應(yīng)用[J].國外橋梁,2000(8):66-70.

[5] 吳 斌,歐麗萍,張紀(jì)剛,等.結(jié)構(gòu)動力可靠度的重要抽樣法[J].計算力學(xué)學(xué)報,2001,18(4):478-482.

[6] 佟曉利,趙國藩.一種與結(jié)構(gòu)可靠度分析幾何法相結(jié)合的響應(yīng)面法[J].土木工程學(xué)報,1997,30(4):51-57.

[7] 陳鐵冰,王書慶,石志源.計入結(jié)構(gòu)幾何非線性影響時斜拉橋可靠度分析[J].同濟大學(xué)學(xué)報,2000,28(4):407-412.

[8] L.Faravelli.A response surface approach for reliability analysis,J.Engrg.Mech.,ASCE,1989,115(12):23-24.

[9] C.G..Bucher,U.Bourgund.A fast and efficient response surface approach for structural reliability problems,Struct.Safety,1990,7(1):55-56.