任麗梅，徐 偉，李戰(zhàn)國(guó)

(西北工業(yè)大學(xué)理學(xué)院，陜西西安 710072)

在地震工程，橋梁工程等領(lǐng)域中動(dòng)力學(xué)系統(tǒng)的結(jié)構(gòu)可靠性一直是研究的熱點(diǎn)，而受到諸如地震、風(fēng)及海浪等隨機(jī)激勵(lì)的動(dòng)力學(xué)系統(tǒng)的首次穿越失效更是研究的焦點(diǎn)之一。首次穿越失效概率(簡(jiǎn)稱:首穿失效概率)是結(jié)構(gòu)安全性的重要指標(biāo)，是指在某個(gè)時(shí)間段內(nèi)結(jié)構(gòu)響應(yīng)首次超過(guò)指定安全域邊界的概率。盡管這個(gè)問(wèn)題已經(jīng)得到了高度重視，但仍是結(jié)構(gòu)安全性分析最具挑戰(zhàn)性的問(wèn)題，因?yàn)橹钡浆F(xiàn)在，即使是最簡(jiǎn)單的線性振子場(chǎng)合也無(wú)可利用的精確解。

自Rice［1-2］的開(kāi)拓性研究以來(lái)，研究者雖然提出了基于Fokker-Planck方程的路徑積分法、胞映射法、有限元等數(shù)值解法［3-5］，但這些方法的復(fù)雜度隨狀態(tài)空間維數(shù)的增加呈指數(shù)增加趨勢(shì)，因此普遍認(rèn)為等效線性化與基于方差縮減技巧的蒙特卡羅模擬方法是解決此類問(wèn)題的選擇［6］，研究者據(jù)此先后提出重要抽樣法，受控蒙特卡羅法［7］，子集模擬法［8］，區(qū)域分解法［9］等，在方差縮減技巧中基于Gisranov定理的重要抽樣法［10-15］已經(jīng)被廣泛應(yīng)用。

基于Gisranov定理的重要抽樣法可以描述為一種控制樣本路徑的方法:樣本被集中在樣本空間中最重要的部分，而不是均勻分布在樣本空間中。然而達(dá)到這樣的目的選擇怎樣的控制函數(shù)才是最關(guān)鍵的。在許多的研究中均是利用設(shè)計(jì)點(diǎn)來(lái)構(gòu)造控制函數(shù)的，但僅在穩(wěn)態(tài)高斯白噪聲激勵(lì)的線性系統(tǒng)場(chǎng)合，設(shè)計(jì)點(diǎn)具有解析表達(dá)式，在非線性系統(tǒng)場(chǎng)合，設(shè)計(jì)點(diǎn)的得到通常采用近似方法［16-18］?？墒窃O(shè)計(jì)點(diǎn)的準(zhǔn)確性直接影響到控制函數(shù)的控制功能，因此本文利用線性系統(tǒng)具有解析的設(shè)計(jì)點(diǎn)這一特性，首先根據(jù)Rice公式，得到與非線性系統(tǒng)方程具有相同平均上穿率的等效線性化系統(tǒng)方程，利用此等效線性化方程得到設(shè)計(jì)點(diǎn)的解析表達(dá)式，構(gòu)造控制函數(shù)。然后將控制函數(shù)運(yùn)用到非線性系統(tǒng)方程中。利用重要樣本法估計(jì)非線性系統(tǒng)的首穿失效概率。最后文章以Duffing振子為例，模擬結(jié)果顯示了方法的正確性與有效性。

1 重要樣本法

設(shè)結(jié)構(gòu)動(dòng)力學(xué)系統(tǒng)響應(yīng)是一p維向量X(t)=［X1(t)，X2(t)，…，Xp(t)］T且滿足 Ito^隨機(jī)微分方程:

Xi(0)=xi，其中W(t)=［W1(t)，W2(t)，…，Wq(t)］T是具有概率測(cè)度p的標(biāo)準(zhǔn)Wiener過(guò)程。m(t，X)，σ(t，X)是漂移、擴(kuò)散系數(shù)，分別滿足李普希茲條件和增長(zhǎng)性條件［19］。定義系統(tǒng)的安全域是D∈Rp。

首穿失效概率為

其中I［g(X)］是示性函數(shù)，

則在確定的初始條件下，基于概率測(cè)度 ~P的首穿失效概率的估計(jì)為

2 等效線性化

考慮零均值，具有單邊譜密度G0的穩(wěn)態(tài)高斯白噪聲激勵(lì)的Duffing振子:

(x(t)(t))的穩(wěn)態(tài)概率密度函數(shù)為

根據(jù)Poisson假設(shè)，利用Rica公式［19］，非線性方程的平均上穿率為:

即:

等效線性方程設(shè)為:

平均上穿率為:

其中:

則解方程:

可得關(guān)于安全域邊界b的等效線性化參數(shù)ωe值。當(dāng)非線性系統(tǒng)的平均上穿率無(wú)法得到解析表達(dá)式時(shí)，可以利用數(shù)值方法。蒙特卡羅方法也可以用來(lái)估計(jì)平均上穿率，在穩(wěn)態(tài)過(guò)程場(chǎng)合平均上穿率被定義為:

其中:μ(b)表示在時(shí)間間隔T內(nèi)，穩(wěn)態(tài)過(guò)程越出安全域的平均次數(shù)。

3 控制函數(shù)的構(gòu)造

在tk時(shí)刻，極限狀態(tài)方程為:

可靠性指標(biāo):

可得設(shè)計(jì)點(diǎn)為:

圖1 等效線性化參數(shù)與安全域邊界Fig.1 The frequencyvs threshold depending on different values of the non-linearity parameter ε.ε =0.1(*)ε=0.5(○)，ε=1(△)，ε=2(◇)

4 結(jié)果與討論

振子方程如式(6)所示，可令X1(t)=X(t)，X2(t)=(t)，其Duffing振子的It o^隨機(jī)微分方程為:

圖2 控制函數(shù)圖。ε=0.1虛線，ε=2實(shí)線Fig.2 Control function for an up-crossing of level b(t)=2.449 40 at time T=40 s.ε =0.1(dotted line)，ε=2(solid line)

圖3 加控制(實(shí)線)與不加控制(虛線)的系統(tǒng)(ε=0.1)響應(yīng)比較圖:Fig.3 The response of the non-linear system(ε =0.1)comparison chart:Add control(solid line)and not controlled(dotted line)

圖4 失效概率圖，b(t)=2.449 40，ε=0.01(* 線)，ε=2(○線)，原始蒙特卡羅法(實(shí)線)Fig.4 The first failure probability simulation comparison chart importance sampling procedure(ε =0.01(*)，ε=2(○))，crude Monte Carlo(solid line)

在穩(wěn)態(tài)高斯白噪聲激勵(lì)的Duffing振子中，令β=0.05，ω0=1.0(rad/s)，γ =0.3，T=40 s，Δt=0.05，n=800，安全域邊界取常數(shù)邊界b(t)=kσ0。圖1是當(dāng)非線性參數(shù)分別為 ε =0.1，ε =0.5，ε =1，ε =2 時(shí)，不同安全域邊界下等效線性化參數(shù)ωe與ω0的比值的示意圖。圖2是當(dāng)ε=0.1(虛線)，ε=2(實(shí)線)的控制函數(shù)圖。圖3是當(dāng)ε=0.1時(shí)Duffing振子加控制(實(shí)線)與不加控制(虛線)的系統(tǒng)相應(yīng)比較圖。從圖3中可看出在20 s之前加控制與不加控制的系統(tǒng)響應(yīng)大致相同，20 s之后控制函數(shù)力量強(qiáng)勁，起到了支配的作用，與振子達(dá)到共振頻率，使得振子在26 s時(shí)首次越過(guò)安全域邊界b(t)=2σ0=2.449 40。圖4是在安全域邊界為b(t)=2σ0時(shí)，利用重要樣本法 ε=0.01(* 線)，ε=2(○線)(N=103)時(shí)的首穿失效概率與原始蒙特卡羅(實(shí)線，N=106)法的比較圖，顯示了兩者很好的吻合，說(shuō)明本文所提方法的正確性與有效性。

5 結(jié)論

本文基于Gisranov定理，提出了非線性結(jié)構(gòu)動(dòng)力系統(tǒng)首穿失效概率的重要抽樣法。本文利用線性系統(tǒng)具有解析的設(shè)計(jì)點(diǎn)這一特性，根據(jù)Rice公式，得到非線性系統(tǒng)的等效線性化方程，等效線性化原理采用平均上穿率而非mean square原理，研究者A.Naess也在研究首次穿越失效問(wèn)題中采用了此線性化原理，他也認(rèn)為此線性化原理對(duì)研究首穿問(wèn)題是更加合理的。利用等效線性化方程得到解析的設(shè)計(jì)點(diǎn)，構(gòu)造控制函數(shù)。然后將控制函數(shù)運(yùn)用到非線性系統(tǒng)方程中，非線性系統(tǒng)不同非線性化參數(shù)ε下進(jìn)行了模擬，結(jié)果說(shuō)明了方法的正確性與有效性，且本文方法所需樣本量及計(jì)算時(shí)間遠(yuǎn)少于蒙特卡羅模擬。

［1］ Rice O C.Mathematical analysis of random noise［J］.Bell Syst Tech J，1944(23):282-332.

［2］ Rice O C.Mathematical analysis of random noise［J］.Bell Syst Tech J，1945(24):46-156.

［3］ Spencer B F，Bergman L A.On the numerical Solution of the Fokker-Planck equation for nonlinear stochastic systems［J］.Nonlinear Dynamic，1993(4):357-372.

［4］Li W，Xu W.Stochastic optimal control of first passage failure for coupled Duffing-Van derPol system under Gaussian white noise excitations ［J］.Chaos.Solitons＆ Fractal，2005，25(5):1221-1228.

［5］Li W，Xu W.First passage problem for strong nonlinear stochastic dynamical system［J］.Chaos.Solitons ＆ Fractal，2006，28(2):414-421.

［6］Proppe C，Pradlwarter H J，Schueller G I.Equivalent linearization and Monte Carlo simulation in Stochastic dynamics［J］.Probab Eng Mech，2003，18(3):1-15.

［7］ Pradlwarter H J，Schueller G I.Assessment of low probability events of dynamical systems by controlled Monte Carlo simulation［J］.Probab Eng Mech，1999，14(3):213-227.

［8］ Zuev K M，Katafygiotis L S.The Horseracing Simulation algorithm for evaluation of small failure probabilities［J］.Probab Eng Mech ，2011，26(2):157-164.

［9］Katafygiotis L，Cheung S H.Domain decomposrtion method for calculating the failure probability of linear dynamic systems subjected toGaussian stochasticloads［J］.Journalof Engineering Mechanics，2006，20:475-486.

［10］ Olsen A I，Naess A.An importance sampling procedure for estimating failure probabilities of Non-linear dynamic systems subjected to random noise［J］.J Non-linear Mech，2007(42):848-863.

［11］ Juu O，Hiroaki T.Importance sampling for stochastic systems under stationary noise having a specified power spectrum［J］.Probab.Eng.Mech.，2009，24(4):537-544.

［12］ Au S K，Beck J L.First excursion probability for linear systems by very efficient importance Sampling［J］.Probab.Eng.Mech.，2001，16(3):193-207.

［13］ Macke M，Bucher C.Importance sampling for randomly excited dynamical systems［J］.J.Sound Vibration.2003，268:269-290.

［14］ Olsen A I，Naess A.An importance sampling procedure for estimating failure probabilities of dynamic systems［C］//ASCE，Joint Speciality Conference on Probabilistic Mechanics and Structural Reliability，2004.

［15］ Pradlwarter H J，Schueller G I.Excursion probabilities of non-linear systems［J］.Int.J.Non Linear Mechanics，2004，39:1447-1452.

［16］ Der Kiureghian A.The geometry of random vibrations and solutions by FORM and SORM［J］.Probab.Engng.Mech.，2000，15(2):81-90.

［17］ Liu P L，Der Kiureghian A.Optimization algorithms for structural reliability［J］.Structural Safety，1991，9(3):161-177.

［18］ Zhang Y，Der Kiureghian A.Two improved algorithms for reliability analysis［C］//Rackwitz et al，editor.Proceedings of the 6thIfIP WG 7.5 conference on optimization of structural systems，1994:297-304.

［19］ Soong T T，Grigoriu M.Random vibration of mechancial and structural systems［M］.Prentice-Hall.Inc.Englewood Cliffs，NJ，1997.