方永鋒 陳建軍 閻彬 曹鴻鈞

(西安電子科技大學(xué) 電子裝備結(jié)構(gòu)設(shè)計教育部重點實驗室，陜西 西安 710071)

有時結(jié)構(gòu)在其服役期承受的載荷和其強度都具有不確定性.隨機分析和模糊分析是目前結(jié)構(gòu)工程中處理不確定性問題時應(yīng)用最為普遍的兩種方法［1-5］，但不論是隨機分析還是模糊分析都對已知信息有較強的依賴性.然而在實際工程中，往往可得到的數(shù)據(jù)貧乏，不足以確定不確定參數(shù)的概率分布函數(shù)或模糊隸屬度函數(shù)，且可靠度的預(yù)測結(jié)果對分布參數(shù)的準確性較為敏感，分布參數(shù)很小的誤差有時可能導(dǎo)致可靠度計算結(jié)果呈現(xiàn)較大的誤差［6］，這在很大程度上限制了隨機和模糊的方法在工程中的應(yīng)用.為此，文獻［7-8］認為，當掌握的不確定性參數(shù)信息較少時，則宜采用區(qū)間或凸集模型來描述問題的不確定性.文獻［9-10］提出了基于區(qū)間分析的結(jié)構(gòu)非概率可靠性模型，并預(yù)測了區(qū)間靜載荷作用一次或特定次數(shù)時結(jié)構(gòu)的可靠度.實際上在結(jié)構(gòu)服役期內(nèi)，隨機型和區(qū)間型不確定因素往往是共同存在的，對于同時含有隨機變量和區(qū)間變量的結(jié)構(gòu)可靠性問題，利用混合模型獲得結(jié)構(gòu)的可靠度，似比給出一確定的單值具有更高的可信度［11］.文獻［12］提出了結(jié)構(gòu)的隨機-區(qū)間問題的順序單回路法(SSLM)處理方法，但只考慮了靜力學(xué)問題.文獻［13］從理論上證明了區(qū)間變量轉(zhuǎn)化為隨機變量與直接計算隨機-區(qū)間變量對于結(jié)構(gòu)可靠度的計算結(jié)果是相同的.文獻［14-15］利用區(qū)間因子研究了隨機-區(qū)間結(jié)構(gòu)系統(tǒng)的靜力學(xué)可靠性問題，得到的結(jié)構(gòu)可靠度是一個常量.然而，由于振動、沖擊、疲勞、腐蝕、老化等外在因素以及一些不確定性內(nèi)在因素的共同影響，結(jié)構(gòu)的可靠性往往是與時間有關(guān)的函數(shù).文獻［16］給出了一種考慮載荷服從某種概率分布且多次作用下零件的可靠性預(yù)測模型，但該模型在實際應(yīng)用時，當作用荷載等幅的條件難以滿足時，導(dǎo)致其可靠性的預(yù)測結(jié)果誤差偏大且難以避免.文獻［17-19］對結(jié)構(gòu)的動力可靠性問題進行了分析，但僅限于作用載荷為平穩(wěn)隨機過程的情況.許多情況下，結(jié)構(gòu)在服役過程中受到的載荷既不是一次也不是僅一種，往往是多次與多種.

文中基于對結(jié)構(gòu)動態(tài)概率可靠性研究的基礎(chǔ)，針對實際問題中掌握作用載荷的統(tǒng)計信息量較少的情況，利用隨機-區(qū)間理論，對結(jié)構(gòu)承受的多次區(qū)間載荷隨時間變化的情況進行了分析.根據(jù)應(yīng)力-強度干涉理論，建立了當作用載荷隨時間變化且結(jié)構(gòu)強度隨時間退化時結(jié)構(gòu)的動態(tài)隨機-區(qū)間可靠性預(yù)測模型，并對模型中的區(qū)間參數(shù)進行隨機化處理，利用一次二階矩法獲得結(jié)構(gòu)的動態(tài)可靠性指標和可靠度.最后通過兩個算例說明了文中預(yù)測模型的合理性、實用性和精確性.

1 動態(tài)載荷與動態(tài)強度

1.1 動態(tài)載荷

作用于結(jié)構(gòu)上的載荷為任意的時變載荷，且載荷的幅值為區(qū)間變量，為了求得在此種情況下結(jié)構(gòu)在任意t 時刻的可靠度，文中對結(jié)構(gòu)承受的載荷作如下處理:

(4)令ΦsI(t)=max{Φ(t)}，ΦsI(t)對應(yīng)的區(qū)間記為sI(t)，將sI(t)稱為載荷作用多次時的最大區(qū)間應(yīng)力，它對應(yīng)的應(yīng)力即為最大區(qū)間應(yīng)力.

若結(jié)構(gòu)在最大區(qū)間載荷下不失效，則結(jié)構(gòu)在這多次區(qū)間載荷下也不會失效.故區(qū)間載荷作用多次時的結(jié)構(gòu)可靠度等價于這些載荷中最大區(qū)間載荷的可靠度.

記

sI(t)的均值與離差分別為

1.2 動態(tài)強度

考慮結(jié)構(gòu)的強度為隨機變量，在結(jié)構(gòu)服役期內(nèi)，由于振動、沖擊、疲勞、腐蝕、老化等外在因素以及一些不確定性內(nèi)在因素的共同影響，導(dǎo)致結(jié)構(gòu)的強度是隨時間退化的.文獻［20］認為可用weibull 分布來描述剩余強度，而且給出了結(jié)構(gòu)在t 時刻的剩余強度為

式中，r(0)是結(jié)構(gòu)的初始強度，s(t)是結(jié)構(gòu)的動態(tài)應(yīng)力，T 是結(jié)構(gòu)的生命周期，e 是材料的退化指數(shù).

根據(jù)式(4)和weibull 分布的性質(zhì)，可得到結(jié)構(gòu)在t 時刻的剩余強度的均值為

式中，μ(r(0))為結(jié)構(gòu)的初始強度的均值，μ(r(1))=μ(s(t))，為結(jié)構(gòu)服役期末的強度的均值.

結(jié)構(gòu)剩余強度的均方差為

式中，σ(r(0))為結(jié)構(gòu)的初始強度的均方差，σ(r(1))為結(jié)構(gòu)服役期末的強度的均方差.

2 結(jié)構(gòu)的區(qū)間-隨機動態(tài)可靠性分析

文獻［13］從理論上證明了如下結(jié)論:對于只含有區(qū)間變量與隨機變量的可靠性問題，把區(qū)間變量轉(zhuǎn)化為均勻分布的隨機變量，就把區(qū)間-隨機可靠性問題轉(zhuǎn)化為傳統(tǒng)的隨機可靠性問題，用一次二階矩法計算可知，兩者的可靠度結(jié)果是等價的.

基于上述結(jié)論，可將式(1)對于結(jié)構(gòu)的可靠度計算轉(zhuǎn)化為

由式(8)可得s(t)的均值和方差分別為

再根據(jù)應(yīng)力-強度干涉理論，則可得到在多次區(qū)間載荷作用下且結(jié)構(gòu)強度退化時的功能函數(shù)為

為了保證用一次二階矩法計算可靠度的精確性，將s(t)和r(t)在任意t 時刻當量標準正態(tài)化處理為(X(t)，Y(t))，如此可求得結(jié)構(gòu)的動態(tài)可靠性指標為

初始值由(10)得到，經(jīng)過反復(fù)迭代計算式(12)，可得任意t 時刻的結(jié)構(gòu)動態(tài)可靠性指標為與(X*(t)，Y*(t))相對應(yīng)的(r(t)，s(t))即為式(11)的最大失效點.

在迭代過程中，為了加速收斂，根據(jù)經(jīng)典響應(yīng)面法的迭代選點方法，則在第k+1 迭代過程中有

這里a 第一次取2.25，以后各次迭代取0.75［18］，

對于r(t)，若服從正態(tài)分布則直接計算，否則采用下式:

這里b 第1 次取3，以后各次迭代取1［19］，r(0)(t)=μ(r(t))

當滿足:

時迭代終止.式中，β(k+1)(t)、β(k)(t)分別為第k +1次和第k 次可靠性指標，ε 為迭代終止要求的精度.

綜上，給出多次區(qū)間載荷作用下且結(jié)構(gòu)隨機強度退化時結(jié)構(gòu)動態(tài)隨機-區(qū)間可靠度的算法如下:

(1)由2.1 求區(qū)間載荷的等價區(qū)間應(yīng)力;

(2)由式(7)將區(qū)間應(yīng)力轉(zhuǎn)化為均勻分布的隨機應(yīng)力;

(3)把上一步轉(zhuǎn)化的隨機應(yīng)力與隨機強度當量標準正態(tài)化;

(4)應(yīng)用一次二階矩法計算可靠性指標β(1)(t);

(5)用式(14)、(15)計算可靠性指標β(k+1)(t);

(6)計算式(16);

(7)若前一步達到了收斂精度，計算終止，否則繼續(xù)由步驟(5)進行迭代，直至滿足收斂精度，迭代終止，輸出β(*)(t)，并據(jù)β(*)(t)值由標準正態(tài)概率分布函數(shù)表獲得結(jié)構(gòu)的動態(tài)可靠度.

3 算例

例1 某汽車后軸在不考慮強度退化的情況下，其強度r(0)服從均值為160 MPa、均方差為15 MPa的正態(tài)分布.由2.1 可得到載荷的工作應(yīng)力s(t)的范圍為45～115 MPa，在載荷作用多次時，用文中方法得到的和利用SSLM 及蒙塔卡羅方法(MCM)模擬106次計算得到的可靠度結(jié)果見表1.

表1 例1 3 種方法計算結(jié)果比較Table 1 Comparison of computation results of example 1 by using three methods

從表1 可以看出，SSLM 和文中方法的計算結(jié)果誤差相當，且與MCM 的計算結(jié)果基本一致，但文中方法的迭代次數(shù)要比SSLM 幾乎少一半，表明該方法的斂速快、精度高.

例2 某車體底座，其強度r(t)隨時間退化，其初始強度r(0)服從均值為150 MPa、均方差為15 MPa的正態(tài)分布，T=10000h，e=4.092.由2.1 可得到載荷的工作應(yīng)力s(t)的范圍為45～115MPa，在載荷作用多次且考慮強度退化時，利用文中方法和MCM模擬計算106次在時間t=0，500，1500，5000，9000 h時的結(jié)構(gòu)動態(tài)可靠度，結(jié)果見表2.

表2 例2 不同時刻兩種方法計算結(jié)果比較Table 2 Comparison of computation results of example 2 by using two methods

從表2 中可以看出，文中提出的方法與MCM兩者的計算結(jié)果幾乎一致;在考慮結(jié)構(gòu)強度退化時，結(jié)構(gòu)在多次區(qū)間載荷作用下的可靠度是隨時間持續(xù)下降的.該結(jié)果可為結(jié)構(gòu)生命的全周期檢測與動態(tài)設(shè)計提供參考依據(jù).

例3 某懸臂管結(jié)構(gòu)如圖1 所示，受到外力F1(t)、F2(t)、P(t)和扭矩W(t)的共同作用.圖中，d 為懸臂管的外直徑，ι 為懸臂管的厚度，L1為懸臂管的長度，L2為F2的力臂長度.給定T=10000 h，e =4.092.該結(jié)構(gòu)中各個不確定參數(shù)的取值見表3，對隨機變量，則分布參數(shù)1 和2 分別表示均值與均方差，對區(qū)間變量，則分布參數(shù)1 和2 分別表示區(qū)間變量的下界和上界.若該結(jié)構(gòu)可靠，則施加于管端面的等效應(yīng)力smax(t)應(yīng)小于結(jié)構(gòu)的強度r(t)，對應(yīng)的極限狀態(tài)方程為

g(r(t)，s(t))=r(t)-smax(t).

圖1 某懸臂管Fig.1 A cantilever tube

表3 懸臂管不確定參數(shù)Table 3 Uncertain variables for the cantilever tube

其中，

式中，θ1、θ2分別是F1、F2與z 軸的夾角，D 為截面積，計算公式為

M(t)為彎矩，計算公式如下

本算例在t=0 時刻，也就是不考慮強度退化的時候，僅考慮一次載荷作用下，文獻［13］用SSLM 方法，迭代次數(shù)16 次，計算可靠度為0.971 6，文中方法迭代6 次即可達到該結(jié)果.在考慮強度退化時，多次載荷作用時文中方法與MCM 模擬106在不同時刻的計算結(jié)果見表4.

從表4 中可以看出，文中方法與MCM 在0、500、1500、5000、9000 h 時兩者的計算結(jié)果基本相同，但迭代次數(shù)明顯要少.同時也說明，在考慮強度退化，外部作用多次時懸臂管的可靠度是隨時間不斷下降的.

表4 例3 不同時刻兩種方法計算結(jié)果比較Table 4 Comparison of computation results of example 3 by using two methods

4 結(jié)論

文中針對多次區(qū)間載荷作用、且結(jié)構(gòu)隨機強度退化的情況構(gòu)建了結(jié)構(gòu)動態(tài)可靠性預(yù)測模型，經(jīng)過理論推理和算例驗證得到如下的結(jié)論:

(1)針對實際問題中掌握作用載荷的統(tǒng)計信息量較少的情況，利用隨機-區(qū)間理論，對結(jié)構(gòu)承受的多次區(qū)間載荷隨時間的變化情況進行了分析和等效處理;

(2)建立了當作用區(qū)間載荷隨時間變化且結(jié)構(gòu)隨機強度隨時間退化時結(jié)構(gòu)的動態(tài)隨機-區(qū)間混合可靠性預(yù)測模型，并對模型中的區(qū)間參數(shù)進行隨機化處理，利用一次二階矩法得到了結(jié)構(gòu)的動態(tài)可靠度.最后通過3 個算例說明了文中模型的合理性、易行性與實用性.

［1］Ditlevsen O，Madsen H O.Structural reliability methods［M］.Chichester:Wiley，1996:135-158.

［2］Madsen H O.First order vs.second order reliability analysis of series structures ［J］.Structures Safety，1985，2(3):207-214.

［3］Madsen H O，Krenk S，Lind N C.Methods of structural safety［M］.Englewood Cliffs:Prentice-Hall，1986:59-78.

［4］Ayyub B M，Lai K L.Structural reliability assessment with ambiguity and vagueness in failure［J］.Naval Engineering Journal，1992，104(3):21-35.

［5］Jiang Qimi，Chen Chun Hsien.A numerical algorithm of fuzzy reliability ［J］.Reliability Engineering ＆ System Safety，2003，80(3):299-307.

［6］Ellishkoff I.Essay on uncertainties in elastic and viscoelastic structures:form A M Freudenthal's criticisms to modern convex modeling ［J］.Computers ＆ Structures，1995，56(6):871-895.

［7］Ben-Haim Y.A non-probabilistic concept of reliability［J］.Structural Safety，1994，14(4):227-245.

［8］Ben-Haim Y，Elishakoff I.Discussion on:a non-probabilistic concept of reliability［J］.Structural Safety，1995，17(3):195-199.

［9］邱志平.非概率集合理論凸方法及其應(yīng)用［M］.北京:國防工業(yè)出版社，2005:23-37.

［10］郭書祥，呂震宙，馮元生.基于區(qū)間分析的結(jié)構(gòu)非概率可靠性模型［J］.計算力學(xué)學(xué)報，2001，18(1):56-60.Guo Shu-xiang，Lu Zhen-zhou，F(xiàn)eng Yuan-sheng.A nonprobabilistic model of structural reliability based on interval analysis ［J］.Chinese Journal of Computational Mechanics，2001，18(1):56-60.

［11］朱增青，陳建軍，梁震濤，等.隨機-區(qū)間型天線結(jié)構(gòu)有限元及可靠性分析［J］.高技術(shù)通訊，2008，16(6):624-629.Zhu Zeng-qing，Chen Jian-jun，Liang Zhen-tao，et al.Finite element and reliability analyses for antenna structures with the mixture of random and interval variables［J］.Chinese High Technology Letters，2008，16(6):624-629.

［12］Du Xiao-ping，Agus Sudjianto，Huang Bei-qing.Reliability-based design with the mixture of random and interval variables［J］.Journal of Mechanical Design，2005，127(12):1068-1077.

［13］Jiang C，Lu G Y，Han X，et al.A new reliability analysis method for uncertain structures with random and interval variables［J］.International Journal of Machine Material Design，2012，6(8):169-182.

［14］Gao W.Interval finite element analysis using interval factor method ［J］.Computational Mechanics，2007，39(4):709-717.

［15］Gao Wei.Natural frequency and mode shape analysis of structures with uncertainty［J］.Mechanical Systems and Signal Processing，2007，21(1):24-39.

［16］王正，謝里陽，李兵.考慮載荷作用次數(shù)的零部件可靠性模型［J］.機械強度，2008，33(1):68-71.Wang Zheng，Xie Li-yang.Li Bing.Reliability model of component considering tomes of random load acting［J］.Journal of Mechanical Strength，2008，30(1):68-71.

［17］胡太彬，陳建軍，高偉，等.平穩(wěn)隨機激勵下隨機桁架結(jié)構(gòu)動力可靠性分析［J］.力學(xué)學(xué)報，2004，36(2):241-246.Hu Tai-bin，Chen Jian-jun，Gao Wei，et al.Dynamic reliability analysis of stochastic truss structures under stationary random excitation［J］.Acta Mechanica Sinica，2004，36(2):241-246.

［18］趙俊釗，陳務(wù)軍，潘欽，等.梁柱結(jié)構(gòu)協(xié)同作用的車輻式索桿體系成形分析［J］.華南理工大學(xué)學(xué)報:自然科學(xué)版，2011，39(12):152-158.Zhao Jun-zhao，Chen Wu-jun，Pan Qin，et al.Analysis of formation of spoke-shaped cable-strut system with synergistic effect of column-beam structures ［J］.Journal of South China University of Technology:Natural Science Edition，2011，39(12):152-158.

［19］韓大建，宋雄彬，王海濤.帶彈性支撐桿張拉結(jié)構(gòu)整體找形分析的動力松弛法［J］.華南理工大學(xué)學(xué)報:自然科學(xué)版，2005，33(11):66-71.Han Da-jian，Song Xiong-bin，Wang Hai-tao.Dynamic relaxation for integer form-finding analysis of tension structures with elastic supporting bars ［J］.Journal of South China University of Technology:Natural Science Edition，2005，33(11):66-71.

［20］Schaff J R，Davidson B D.Life prediction methodology for composite structures，Part I-constant amplitude and two-stress level fatigue［J］.Journal of Composite Materials，1997，31(2):128-157.