基于降維算法的結(jié)構(gòu)混合可靠性分析方法

孟廣偉，馮昕宇，周立明，李鋒

基于降維算法的結(jié)構(gòu)混合可靠性分析方法

孟廣偉，馮昕宇，周立明，李鋒

(吉林大學(xué) 機(jī)械與航空航天工程學(xué)院，吉林 長(zhǎng)春，130025)

針對(duì)工程實(shí)際中同時(shí)帶有模糊變量與隨機(jī)變量的結(jié)構(gòu)混合不確定性問題，提出1種基于降維算法的混合不確定性變量的可靠性分析模型。首先，利用模糊數(shù)學(xué)中的?截集概念，將模糊變量轉(zhuǎn)變?yōu)樗浇丶孪鄳?yīng)的區(qū)間變量，再借助于降維算法，將含有個(gè)隨機(jī)變量的結(jié)構(gòu)功能函數(shù)展開為個(gè)一維隨機(jī)變量函數(shù)；將所得到的降維表達(dá)式進(jìn)行泰勒展開，得到結(jié)構(gòu)功能函數(shù)的上、下界表達(dá)式；運(yùn)用變量轉(zhuǎn)換方法將其中的隨機(jī)變量轉(zhuǎn)換為均值為0，方差為0.5的正態(tài)分布變量，并結(jié)合二項(xiàng)式展開定理、Gauss?Hermite積分方法與變量轉(zhuǎn)換方法計(jì)算出結(jié)構(gòu)功能函數(shù)上下界的統(tǒng)計(jì)矩；將所得的矩信息應(yīng)用到Edgeworth級(jí)數(shù)展開式中，計(jì)算得到對(duì)應(yīng)于?截集的失效概率區(qū)間，從而獲得失效概率的隸屬度函數(shù)。研究結(jié)果表明：本文提出的方法不僅計(jì)算精度高，而且計(jì)算量小。

結(jié)構(gòu)可靠性分析；混合可靠性；降維算法；模糊變量；Edgeworth級(jí)數(shù)

結(jié)構(gòu)可靠度分析是工程實(shí)踐中的重點(diǎn)研究?jī)?nèi)容之一[1?6]。帶有隨機(jī)變量的失效概率計(jì)算問題是概率可靠性設(shè)計(jì)中的重點(diǎn)研究?jī)?nèi)容之一，帶有模糊變量的失效概率計(jì)算問題同樣也是模糊可靠性設(shè)計(jì)中的重點(diǎn)問題。然而在工程實(shí)踐中，往往僅帶有單一類型變量的結(jié)構(gòu)比較少，混合類型變量共同存在的情況較多。目前，已有很多方法用于求解帶有隨機(jī)變量和模糊變量的結(jié)構(gòu)失效概率f，例如，CHAKRABORTY等[7?8]基于信息熵方法，將模糊隸屬函數(shù)轉(zhuǎn)換為等效概率密度函數(shù)進(jìn)行可靠性分析；SMITH等[9]基于尺度變換方法，將傳統(tǒng)概率可靠性分析理論應(yīng)用于混合可靠性分析；ADDURI等[10?12]基于概率理論利用一次二階矩法和一次二階矩法的改進(jìn)方法求解混合失效概率等。在利用上述方法計(jì)算時(shí)，信息熵法在理論上與模糊可靠度的結(jié)果存在一定偏差，且喪失了模糊性。另外，尺度變換法對(duì)模糊數(shù)類型的要求較為嚴(yán)格；一次二階矩法在非線性情況下，存在迭代收斂慢甚至不收斂的情況。因此，本文作者結(jié)合降維算法和Edgeworth級(jí)數(shù)展開法[13]，針對(duì)同時(shí)含有隨機(jī)變量和模糊變量的混合結(jié)構(gòu)，提出1種簡(jiǎn)便的計(jì)算失效概率方法：通過?截集概念將模糊變量轉(zhuǎn)化為水平截集下相應(yīng)的區(qū)間變量；利用降維算法，將功能函數(shù)分解為個(gè)一維隨機(jī)變量函數(shù)疊加的形式；借助于泰勒展開，得到功能函數(shù)的上、下界；再結(jié)合變量轉(zhuǎn)換法、Gauss-Hermite數(shù)值積分計(jì)算獲得功能函數(shù)上、下界的矩信息，并將其應(yīng)用到Edgeworth級(jí)數(shù)展開式中，最終計(jì)算得到結(jié)構(gòu)功能函數(shù)的失效概率隸屬度。

1 模糊變量λ?截集

對(duì)于同時(shí)帶有隨機(jī)變量和模糊變量且兩者相互獨(dú)立的結(jié)構(gòu)，其結(jié)構(gòu)功能函數(shù)可表示為

圖1 模糊變量與λ?截集關(guān)系

2 降維算法

h與0的上、下界分別為

式中：(x)為隨機(jī)變量x的概率密度函數(shù)。

3 變量轉(zhuǎn)換

針對(duì)式(19)，可運(yùn)用傳統(tǒng)的數(shù)值積分公式進(jìn)行計(jì)算。計(jì)算積分的方法有很多種，其中值得提出的有快速傅里葉變換法和蒙特卡羅方法(MCS)。前者的前提條件較為嚴(yán)格，要求存在解析解且要求功能函數(shù)相互獨(dú)立；后者是較為普遍的方法，雖然簡(jiǎn)單易行但計(jì)算量相當(dāng)大，耗時(shí)較長(zhǎng)，其工程應(yīng)用范圍受到了一定限制。因此，本文作者引入變量轉(zhuǎn)換的方法，將相互獨(dú)立的各隨機(jī)變量均轉(zhuǎn)換為正態(tài)分布變量(即均值為0，方差為0.5的正態(tài)分布變量)。為確保準(zhǔn)確性，假設(shè)Z(=1，2，…，)是相互獨(dú)立的，則變量轉(zhuǎn)換的表達(dá)式可寫為

4 Edgeworth級(jí)數(shù)

計(jì)算結(jié)構(gòu)失效概率的前提是具備功能函數(shù)的概率密度函數(shù)或者其聯(lián)合概率密度函數(shù)表達(dá)式，但是在工程實(shí)踐中，很難通過大量的數(shù)據(jù)得出精確的概率密度函數(shù)表達(dá)式。因而，可選用擬合功能函數(shù)的概率密度函數(shù)或累積分布函數(shù)的方式近似求解復(fù)雜結(jié)構(gòu)的失效概率。本文選用Edgeworth級(jí)數(shù)方法進(jìn)行擬合，將式(17)與式(18)計(jì)算得到的功能函數(shù)原點(diǎn)矩信息代入式(28)與式(29)，最終將所得中心矩信息作為Edgeworth級(jí)數(shù)展開式系數(shù)，擬合結(jié)構(gòu)功能函數(shù)相應(yīng)的累積分布函數(shù)表達(dá)式。

5 數(shù)值算例

十桿桁架結(jié)構(gòu)如圖2所示。桿件長(zhǎng)為=3.6 m；于節(jié)點(diǎn)2與節(jié)點(diǎn)4處施加豎直向下載荷(為外載荷，其服從均值為710 kN，變異系數(shù)為0.03的正態(tài)分布)；水平方向的桿①~④面積為1m2，豎直方向的桿⑤~桿⑥面積為2m2，斜向的桿⑦~桿⑩面積為3m2(1，2，3和均為模糊變量)；為材料的彈性模量，其服從均值為2.1×1011Pa，變異系數(shù)為0.01的正態(tài)分布。當(dāng)在節(jié)點(diǎn)2處的垂直位移max不大于允許位移allow(allow=0.004 2 m)時(shí)，結(jié)構(gòu)功能函數(shù)可表示為：

MCS方法在工程實(shí)踐中被廣泛認(rèn)可為精確解。表1所示為本文方法與MCS方法的十桿桁架失效概率對(duì)比。圖3所示為本文方法與MCS方法的失效概率?隸屬度的關(guān)系。由表1和圖3可知：本文方法與MCS方法計(jì)算結(jié)果非常接近。但在每個(gè)?截集處，本文方法的樣本點(diǎn)數(shù)僅為6個(gè)，遠(yuǎn)低于MCS方法的106個(gè)，體現(xiàn)了本文方法的簡(jiǎn)單性和有效性。

圖2 平面十桿桁架結(jié)構(gòu)

表1 2種方法十桿桁架失效概率對(duì)比

1—MCS；2—本文方法(即Edgeworth)。

6 結(jié)論

1) 提出1種針對(duì)同時(shí)帶有模糊變量和隨機(jī)變量不確定性問題的結(jié)構(gòu)失效概率計(jì)算方法。

2) 借助于?截集概念，將模糊變量轉(zhuǎn)變?yōu)樗浇丶孪鄳?yīng)的區(qū)間變量，將隨機(jī)模糊問題有效地轉(zhuǎn)變?yōu)殡S機(jī)區(qū)間問題。

3) 利用降維算法建立結(jié)構(gòu)功能函數(shù)的一維隨機(jī)變量降維表達(dá)式，無需借助于多重積分求解結(jié)構(gòu)功能函數(shù)的統(tǒng)計(jì)矩，無需求解結(jié)構(gòu)功能函數(shù)矩陣的逆。

4) 高效且穩(wěn)定地降低計(jì)算成本，有效避免了由于迭代而存在的收斂慢甚至不收斂情況的問題。

[1] 柳詩(shī)雨, 呂震宙, 員婉瑩, 等. 小失效概率情況下的全局可靠性靈敏度分析的高效方法[J]. 航空學(xué)報(bào), 2016, 37(9): 2766?2774. LIU Shiyu, Lü Zhenzhou, YUAN Wanying, et al. Efficient method for global reliability sensitivity analysis with small failure probability[J]. Acta Aeronautica Et Astronautica Sinica, 2016, 37(9): 2766?2774.

[2] 孫文彩, 楊自春, 李軍. 含隨機(jī)－區(qū)間耦合變量的結(jié)構(gòu)可靠性分析方法[J]. 華中科技大學(xué)學(xué)報(bào)(自然科學(xué)版), 2015, 43(4): 44?48. SUN Wencai, YANG Zichun, LI Jun. Reliability analysis for structures with random-interval coupling variables[J]. Huazhong University of Science & Technology(Natural Science Edition), 2015, 43(4): 44?48.

[3] 鄭燦赫, 孟廣偉, 李鋒, 等. 一種基于 SAPSO-DE 混合算法的結(jié)構(gòu)非概率可靠性優(yōu)化設(shè)計(jì)[J]. 中南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2015, 46(5): 1628?1634. JONG Chanhyok, MENG Guangwei, LI Feng, et al. Structural non-probabilistic reliability optimization design based on a SAPSO-DE hybrid algorithm[J]. Journal of Central South University (Science and Technology), 2015, 46(5): 1628?1634.

[4] MOJSILOVI? N, STEWART M G. Probability and structural reliability assessment of mortar joint thickness in load-bearing masonry walls[J]. Structural Safety, 2015, 52: 209?218.

[5] 韓志杰, 王璋奇. 基于區(qū)間有限元的吊梁非概率可靠性研究及敏感性分析[J]. 中南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2012, 43(5): 1746?1742. HAN Zhijie, WANG Zhangqi. Non-probabilistic reliability research and sensitivity analysis of hanging beam based on interval finite element method[J]. Journal of Central South University (Science and Technology), 2012, 43(5): 1746?1742.

[6] 李少宏, 陳建軍. 一種新的改進(jìn)響應(yīng)面法的結(jié)構(gòu)可靠性計(jì)算方法[J]. 中南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2013, 44(5): 1837?1841. LI Shaohong, CHEN Jianjun. A new improved response surface method for structural reliability computing[J]. Journal of Central South University (Science and Technology), 2013, 44(5): 1837?1841.

[7] CHAKRABORTY S, SAM P C. Probabilistic safety analysis of structures under hybrid uncertainty[J]. International Journal for Numerical Methods in Engineering, 2007, 70(4): 405?422.

[8] 譚巍, 李冬, 樊照遠(yuǎn), 等. 基于模糊信息熵的航空發(fā)動(dòng)機(jī)性能評(píng)估和可靠性分析[J]. 航空發(fā)動(dòng)機(jī), 2011, 37(5): 45?48. TAN Wei, LI Dong, FAN Zhaoyuan, et al. Aeroengine performance synthetic estimation and reliability analysis based on fuzzy information entropy[J]. Aeroengine, 2011, 37(5): 45?48.

[9] SMITH S A, KRISHNAMURTHY T, MASON B H. Optimized vertex method and hybrid reliability[C]//Americal Institute of Aerospace and Astronautics: AIAA-2002-1645. 43rd AIAA/ ASME/ASCE/AHS/ASC Structures: Structural Dynamics and Materials Conference. Denver, Colorado, USA: American Institute of Aerospace and Astronautics Inc, 2002: 22?25.

[10] ADDURI P, PENMETSA R, GRANDHI R. Membership function development for reliability analysis with mixed uncertain variables[C]//American Institute of Aerospace and Astronautics: AIAA-2005-2072. 46th AIAA/ASME/ASCE/ AHS/ASC Structures: Structural Dynamics and Materials Conference: Texas, USA: American Institute of Aerospace and Astronautics Inc, 2005: 18?21.

[11] ZHAO Yangang G, ONO T. Moment methods for structural reliability[J]. Structural Safety, 2001, 23(1): 47?75.

[12] 董玉革, 趙征權(quán). 模糊可靠性分析的一次二階矩法[J]. 合肥工業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版), 2005, 28(9): 980?984. DONG Yuge, ZHAO Zhengquan. Fuzzy reliability analysis based on the FOSM method[J]. Journal of Hefei University of Technology (Natural Science), 2005, 28(9): 980?984.

[13] KONG Chunyang, SUN Zhigang, NIU Xuming, et al. Moment methods for C/SiC woven composite components reliability and sensitivity analysis[J]. Science and Engineering of Composite Materials, 2014, 21(1): 121?128.

[14] LI Guijie, LU Zhenzhou, XU Jia. A fuzzy reliability approach for structures based on the probability perspective[J]. Structural Safety, 2015, 54: 10?18.

[15] DAI Hongzhe, ZHANG Hao, WANG Wei. A Multiwavelet neural network-based response surface method for structural reliability analysis[J]. Computer-aided Civil and Infrastructure Engineering, 2015, 30(2): 151?162.

[16] CHATTERJEE S, SINGH J B, ROY A. A structure-based software reliability allocation using fuzzy analytic hierarchy process[J]. International Journal of Systems Science, 2015, 46(3): 513?525.

[17] PERI?ARO G A, SANTOS S R, RIBEIRO A A, et al. HLRF-BFGS optimization algorithm for structural reliability[J]. Applied Mathematical Modelling, 2015, 39(7): 2025?2035.

[18] RAHMAN S, XU Heqin. A univariate dimension-reduction method for multi-dimensional integration in stochastic mechanics[J]. Probabilistic Engineering Mechanics, 2004, 19(4): 393?408.

[19] YOUN B D, XI Zhimin, WANG Pingfeng. Eigenvector dimension reduction (EDR) method for sensitivity-free probability analysis[J]. Structural and Multidisciplinary Optimization, 2008, 37(1): 13?28.

[20] LI Gang, ZHANG Kai. A combined reliability analysis approach with dimension reduction method and maximum entropy method[J] Structural and Multidisciplinary Optimization, 2011, 43(1): 121?143.

[21] ZHANG Xufang, PANDEY M D, ZHANG Yimin. A numerical method for structural uncertainty response computation[J]. Science China Technological Sciences, 2011, 54(12): 3347?3357.

(編輯 伍錦花)

Hybrid structural reliability analysis method based on dimension reduction algorithm

MENG Guangwei, FENG Xinyu, ZHOU Liming, LI Feng

(School of Mechanical and Aerospace Engineering, Jilin University, Changchun 130025, China)

A new hybrid reliability analysis method based on dimension reduction algorithm was proposed, aiming at solving the hybrid uncertain problems of structures with both fuzzy variables and random variables in engineering practice. Firstly, fuzzy variables were transformed into corresponding interval variables within the level set using the concept of?cut set in fuzzy mathematics. Taking advantage of dimension reduction method, structural performance functionwithrandom variables was expanded tofunctions, each with one random variable. This dimension reduction formula is then expanded by Taylor method to get the expression of the upper and lower limits of the performance function. By using the variable transformation method, the random variables could be changed to mutual independent normal variables, with mean value zero and the variance 0.5. Combining the binomial theorem, Gauss-Hermite integration method and variable transformation method, the statistical moments of the upper and lower limits of the structural performance function were computed. After substituting the moment information into the Edgeworth series expansion formula, the intervals of structural failure probability corresponding to the λ?cut set were obtained. Therefore, the membership function of the failure probability was calculated. The results show that the proposed method has high precision with low computational cost.

structural reliability; hybrid reliability; dimension reduction method; fuzzy variable; Edgeworth series

TB114.3

A

1672?7207(2018)08?1944?06

10.11817/j.issn.1672?7207.2018.08.015

2017?08?26；

2017?10?26

國(guó)家自然科學(xué)基金資助項(xiàng)目(51305157)；吉林省科技廳基金資助項(xiàng)目(20160520064JH)(Project(51305157) supported by the National Natural Science Foundation of China; Project(20160520064JH) supported by the Science and Technology Department Fund of Jilin Province)

李鋒，博士，副教授，從事結(jié)構(gòu)可靠性研究：E-mail：fengli@jlu.edu.cn