多級載荷累積損傷下結(jié)構(gòu)的動態(tài)可靠性分析

孫劍萍 湯兆平 羅意平

1.華東交通大學(xué)交通運(yùn)輸與物流學(xué)院，南昌，3300132.中南大學(xué)交通運(yùn)輸工程學(xué)院，長沙,410075

0 引言

結(jié)構(gòu)的可靠性往往與時間密切相關(guān)，主要涉及強(qiáng)度退化、累積損傷等因素，近年來國內(nèi)外學(xué)者紛紛從概率、模糊及非概率的角度揭示其演變過程。王新剛等[1]研究了變幅隨機(jī)載荷和強(qiáng)度退化下結(jié)構(gòu)的動態(tài)可靠性問題。秦大同等[2]建立隨機(jī)風(fēng)引起的系統(tǒng)外部載荷激勵作用下的風(fēng)力發(fā)電機(jī)傳動系的動態(tài)可靠性模型。方永鋒等[3]運(yùn)用Turkstra方法分析了4種隨機(jī)應(yīng)力組合作用，以及結(jié)構(gòu)強(qiáng)度不退化和退化下的動態(tài)可靠性。徐芳等[4]對外部隨機(jī)風(fēng)載、內(nèi)部齒輪時變嚙合剛度、軸承時變剛度及綜合傳遞誤差等耦合進(jìn)行研究，并得出失效相關(guān)性和強(qiáng)度退化對可靠度的影響規(guī)律。HAJIALIZADEH等[5]研究了載荷的空間以及剩余強(qiáng)度變化下橋梁的動態(tài)可靠性評估。王磊等[6]將模糊變量進(jìn)行等概率當(dāng)量轉(zhuǎn)換，建立了模糊動態(tài)可靠性分析模型。高明君等[7]提出了一種模糊動態(tài)可靠性的建模方法。概率或模糊動態(tài)可靠性模型的求解主要基于Monte Carlo方法、靜態(tài)轉(zhuǎn)換和跨越率方法。張宇貽等[8]運(yùn)用Monte Carlo方法計算了平穩(wěn)二項(xiàng)動載和隨機(jī)恒載聯(lián)合作用下的時變可靠度。HU等[9]對時變可靠度問題進(jìn)行時不變轉(zhuǎn)換和求解，提高了計算效率。MOURELATOS等[10]結(jié)合隨機(jī)振動理論、總概率定理及使用向上跨越和接合點(diǎn)向上跨越率積分方程進(jìn)行振動梁的動態(tài)可靠性研究。

概率或模糊方法均需依賴大量的原始數(shù)據(jù)。工程實(shí)際中往往難以獲得參數(shù)的精確值，但其界限較易確定。近年來，區(qū)間方法成為處理有界不確定問題的有效途徑。姜潮等[11]針對概率與區(qū)間混合不確定問題，構(gòu)造一種含區(qū)間變量的動態(tài)可靠性分析方法。魏宗平等[12]提出了隨載荷作用次數(shù)而變化的區(qū)間動態(tài)可靠性指標(biāo)計算方法。方永鋒等[13]研究了階梯型、等幅交變型和任意時變型三種載荷作用下結(jié)構(gòu)的區(qū)間動態(tài)可靠性。張德權(quán)等[14]基于PHI2提出了一種求解算法，可獲得設(shè)計周期內(nèi)的可靠性區(qū)間。

考慮到可靠性隨載荷作用次數(shù)的影響，人們進(jìn)行了進(jìn)一步研究。方永鋒等[15]對多次區(qū)間載荷進(jìn)行當(dāng)量標(biāo)準(zhǔn)正態(tài)化處理，將區(qū)間隨機(jī)可靠性問題轉(zhuǎn)化為隨機(jī)可靠性問題求解。王正等[16]運(yùn)用概率微分方程和全概率公式建立了隨機(jī)載荷反復(fù)作用下的動態(tài)可靠性模型，為準(zhǔn)確評估可靠性提供了方法，但計算的是最大載荷下的可靠度。以上研究均未考慮不同級別的載荷作用次序而導(dǎo)致的非線性疲勞損傷對可靠性影響的差異，因此預(yù)測結(jié)果存在較大誤差[17]。

本文在現(xiàn)有動態(tài)可靠性理論的基礎(chǔ)上，研究不同級別載荷的不同作用次序而導(dǎo)致不同疲勞累積損傷的等效轉(zhuǎn)化問題，為多級載荷累積損傷下結(jié)構(gòu)的動態(tài)可靠性分析提供模型和方法。

1 剩余強(qiáng)度的描述及其退化模型

當(dāng)強(qiáng)度指標(biāo)(如疲勞強(qiáng)度)在載荷作用的整個過程中不斷變化時，系統(tǒng)的可靠性建模應(yīng)當(dāng)考慮強(qiáng)度退化的影響。以疲勞失效問題為例，強(qiáng)度指標(biāo)隨著載荷作用次數(shù)的增加而逐漸減小，即剩余強(qiáng)度不斷減小。通常，強(qiáng)度退化與載荷幅值、作用次數(shù)以及加載順序等因素有關(guān)。

材料的強(qiáng)度在實(shí)際工作中有一個隨時間或載荷作用次數(shù)的增加而下降的非線性過程，即強(qiáng)度退化過程。剩余強(qiáng)度不僅與結(jié)構(gòu)服役的時刻t有關(guān)，還與加載的應(yīng)力水平有關(guān)，即

r(t)=f(t,sp,R) 0≤t≤T

(1)

式中,T為服役周期(周期壽命)；r(t)為結(jié)構(gòu)在t時刻的剩余強(qiáng)度，r(0)表示t=0時刻結(jié)構(gòu)的強(qiáng)度即初始強(qiáng)度，r(T)表示結(jié)構(gòu)服役周期T時刻的強(qiáng)度；sp為疲勞應(yīng)力峰值；R為載荷的應(yīng)力比。

載荷應(yīng)力比R即為工作應(yīng)力的最小值smin與最大值smax之比，其表達(dá)式如下：

(2)

載荷應(yīng)力比R與服役周期T直接相關(guān)。當(dāng)實(shí)際載荷的應(yīng)力比變化時，根據(jù)Goodman疲勞極限圖將實(shí)際載荷轉(zhuǎn)化為實(shí)驗(yàn)載荷，再用實(shí)驗(yàn)s-N曲線進(jìn)行疲勞壽命評估。載荷應(yīng)力比R對材料的強(qiáng)度退化系數(shù)γ(γ>0)也產(chǎn)生直接影響，姚衛(wèi)星[18]對相同材料在不同載荷應(yīng)力比R下強(qiáng)度退化系數(shù)γ的變化情況進(jìn)行了研究。顯然，載荷應(yīng)力比R對剩余強(qiáng)度的影響可納入服役周期T及材料強(qiáng)度退化系數(shù)γ兩個參數(shù)中。

根據(jù)熱力學(xué)的觀點(diǎn)，損傷是不可逆的，對疲勞過程中強(qiáng)度的退化模型來說，剩余強(qiáng)度是單調(diào)遞減的函數(shù)。SCHAFF等[19]提出在常幅載荷作用下，剩余強(qiáng)度按冪指數(shù)規(guī)律變化，其關(guān)系表達(dá)式如下：

(3)

對于不同材料，其剩余強(qiáng)度退化規(guī)律不同，強(qiáng)度退化系數(shù)γ的取值就存在差異。當(dāng)γ取不同值時，剩余強(qiáng)度曲線也不同，如圖1所示。由圖1可以看出，當(dāng)γ=1時，式(3)退化成線性規(guī)律形式；當(dāng)γ>1時，r(t)呈現(xiàn)開始階段退化較慢，后期迅速下降的特點(diǎn)；當(dāng)γ→∞時，r(t)退化成折線式規(guī)律形式。對于金屬材料，從疲勞機(jī)理上看，在疲勞加載初期，元件在疲勞載荷作用下產(chǎn)生的缺陷如位錯、滑移、空洞等對元件強(qiáng)度的影響較小，剩余強(qiáng)度衰減得較慢，但在后期其缺陷如裂紋將使元件強(qiáng)度迅速減小，從而產(chǎn)生破壞現(xiàn)象,故金屬材料的強(qiáng)度退化系數(shù)γ>1。

圖1 剩余強(qiáng)度退化規(guī)律Fig.1 The law of residual strength degradation

2 多級載荷下疲勞損傷累積規(guī)律與剩余強(qiáng)度

累積損傷規(guī)律是疲勞研究中的關(guān)鍵問題之一。線性疲勞累積損傷理論認(rèn)為，對于任一應(yīng)力水平，不論在壽命的前期或后期，每次循環(huán)的疲勞損傷量應(yīng)是相等的。但事實(shí)上，疲勞損傷不是線性的，要想較好地反映結(jié)構(gòu)件的真實(shí)可靠性，需合理考慮非線性疲勞損傷累積規(guī)律。

有限壽命設(shè)計法主要基于疲勞累積損傷理論?？梢宰C明，式(3)必定滿足下列兩個邊界條件：

r(0)=σb

(4)

r(T)=sp

(5)

式中，σb為材料的靜拉伸強(qiáng)度。

(6)

可以驗(yàn)證：

D(0)=0

(7)

(8)

可見上述定義符合損傷的一般意義,能夠反映損傷演化的非線性特性。

圖2中，OA和OBC分別為在應(yīng)力峰值sp1和sp2作用下，材料發(fā)生的損傷量隨作用時間的變化曲線，其對應(yīng)的服役周期分別為T1和T2。對于金屬材料，當(dāng)sp1>sp2時，T1

(a)高-低載荷順序 (b)低-高載荷順序圖2 兩級載荷順序的疲勞損傷演化規(guī)律Fig.2 Fatigue damage evolution law of two levelloading sequence

由圖2可以看出，材料在sp1下作用時間t1和在sp2下作用時間T2-t2的疲勞損傷量相等。根據(jù)式(8)可得

(9)

可簡化為

(10)

(1)高-低順序加載時，有sp1>sp2，γ1<γ2,故

(11)

結(jié)合式(10)，則

(12)

故

(13)

(14)

假設(shè)結(jié)構(gòu)承受兩級載荷，在一級應(yīng)力sp1作用時間t1后，根據(jù)疲勞損傷量等效等值壽命原則，可換算成二級應(yīng)力sp2的等效作用時間t21，推導(dǎo)如下：

因?yàn)?/p>

所以

(15)

顯然，結(jié)構(gòu)受一級應(yīng)力和二級應(yīng)力作用之后，產(chǎn)生的損傷量等效于原結(jié)構(gòu)在二級應(yīng)力下連續(xù)作用時間t21+t2產(chǎn)生的總損傷量，即

(16)

則結(jié)構(gòu)受一級應(yīng)力和二級應(yīng)力作用之后，轉(zhuǎn)化為二級應(yīng)力下的剩余壽命為

(17)

在二級應(yīng)力下作用時間t2后，根據(jù)疲勞損傷量等效等值壽命原則，可換算成三級應(yīng)力的等效作用時間t32，推導(dǎo)如下：

因?yàn)?/p>

所以

(18)

同理，可推算出第i級應(yīng)力spi的等效作用時間ti(i-1)為

(19)

i=1,2,…，n

式中，ti-1為第i-1級應(yīng)力sp(i-1)的等效作用時間;t(i-1)(i-2)為第i-1級應(yīng)力的等效作用時間;Ti-1為結(jié)構(gòu)在應(yīng)力sp(i-1)作用下的生命周期。

根據(jù)式(3)，可計算出第i級應(yīng)力的等效作用時間下的剩余強(qiáng)度為

(20)

(21)

(22)

3 基于區(qū)間分析的動態(tài)非概率可靠性指標(biāo)

在實(shí)際服役過程中，受結(jié)構(gòu)自身的材料性能、使用時間、荷載效應(yīng)以及所處環(huán)境等因素影響，強(qiáng)度r和應(yīng)力s都是一個動態(tài)的過程，隨時間與結(jié)構(gòu)參數(shù)變化，且強(qiáng)度和應(yīng)力為區(qū)間變量，同時因強(qiáng)度和應(yīng)力還隨著疲勞載荷作用的次數(shù)或時間而變化，故強(qiáng)度和應(yīng)力的區(qū)間范圍也在不斷變化。強(qiáng)度r和應(yīng)力s可表示為t與x的函數(shù)：

r=r(x,t)∈rI(x,t)=[rl(x,t),ru(x,t)]

(23)

s=s(x,t)∈sI(x,t)=[sl(x,t),su(x,t)]

(24)

式中,ru(x,t)、rl(x,t)、su(x,t)、sl(x,t)分別表示結(jié)構(gòu)時變強(qiáng)度r(x,t)和時變應(yīng)力s(x,t)的上下界。

圖3所示為考慮結(jié)構(gòu)抗力與荷載效應(yīng)的時變性，用區(qū)間模型來描述不確定性變量的結(jié)構(gòu)時變非概率可靠性分析的輸入和輸出過程。

圖3 時變可靠性分析的輸入和輸出Fig.3 Input and output of time-dependent reliability analysis

當(dāng)考慮結(jié)構(gòu)荷載效應(yīng)與抗力的時變性時，基于應(yīng)力-強(qiáng)度干涉理論，以應(yīng)力極限狀態(tài)表示的結(jié)構(gòu)動態(tài)功能函數(shù)為

M=f(x,t)=g(r(x,t),s(x,t))=r(x,t)-s(x,t)

(25)

由區(qū)間非概率可靠性定義式，可推導(dǎo)得出在t時刻結(jié)構(gòu)處于安全狀態(tài)即M>0的非概率可靠性指標(biāo)為

(26)

(27)

(28)

(29)

(30)

式中，rc(x,t)、rr(x,t)分別為時變強(qiáng)度區(qū)間變量的中心和半徑；sc(x,t)、sr(x,t)分別為時變應(yīng)力區(qū)間變量的中心和半徑。

當(dāng)僅考慮強(qiáng)度隨作用時間衰減時，時變強(qiáng)度r(x,t)可以用r(t)代替,式(26)可簡化為

(31)

(32)

結(jié)構(gòu)剩余強(qiáng)度的中心和半徑分別為

(33)

則非概率可靠性指標(biāo)η為

(34)

M=sp-s(x,t)

(35)

在t時刻結(jié)構(gòu)處于安全狀態(tài)即M>0的非概率可靠性指標(biāo)η為可表示為

(36)

根據(jù)非概率可靠性指標(biāo)的概念[21]，η的值越大，結(jié)構(gòu)安全程度越高。當(dāng)η>1時，結(jié)構(gòu)安全可靠；當(dāng)η<-1時，結(jié)構(gòu)必然失效；當(dāng)-1<η<1時，結(jié)構(gòu)可能安全，也可能失效；當(dāng)η為1和-1時，結(jié)構(gòu)分別處于必定安全可靠和必定失效的臨界狀態(tài)。

4 多級載荷作用下結(jié)構(gòu)動態(tài)非概率可靠性

預(yù)測模型

基于應(yīng)力-強(qiáng)度干涉理論，以應(yīng)力極限狀態(tài)表示的結(jié)構(gòu)在ti時刻的動態(tài)功能函數(shù)為

M=g(r(x,t),s(x,t))=r(x,t)-s(t)=
r(ti(i-1)+ti)-spi

(37)

(38)

(39)

則非概率可靠性為

(40)

5 算例

算例研究材料為45鋼拉伸樣件在三級載荷作用下的疲勞拉伸動態(tài)可靠性。為便于對本文提出的多級載荷作用下結(jié)構(gòu)動態(tài)非概率可靠性預(yù)測模型與方法進(jìn)行對比分析與評估，選用文獻(xiàn)[18,22]中部分?jǐn)?shù)據(jù)，對結(jié)構(gòu)在一、二、三級載荷下的服役周期T1、T2、T3的區(qū)間寬度進(jìn)行設(shè)置，根據(jù)“常見通用設(shè)備經(jīng)濟(jì)壽命參考年限表”中經(jīng)濟(jì)壽命上下界關(guān)系[23]，區(qū)間上界值約為下界值的125%。具體如下：r(0)=[791.92,875.28]MPa，先在第一級載荷應(yīng)力區(qū)間[301.15,332.85]MPa下作用t1=640次循環(huán)，T1=[29 025,35 722]，γ1=4.092。再在第二級載荷應(yīng)力區(qū)間[285,315]MPa下作用t2=1 000次循環(huán)，T2=[64 116,80 812]，γ2=4.121。最后在第三級載荷應(yīng)力區(qū)間[264.1,291.9]MPa下作用t3=6 000,18 000,30 000,42 000,54 000,66 000次循環(huán)，直至達(dá)到安全臨界狀態(tài)，T3=[145 971,185 607]，γ3=4.186。

5.1 第一級載荷作用

根據(jù)式(6)，在第一級載荷應(yīng)力區(qū)間[301.15,332.85]MPa下作用t1=640次循環(huán)，產(chǎn)生損傷量：

(41)

則在第一級載荷作用時間t1=640次循環(huán)后，剩余強(qiáng)度為

(42)

由式(27)～式(30)，計算可得rc(t1)=833.599 9、sc(x,t)=317，rr(t1)=41.680 0、sr(x,t)=15.85，根據(jù)式(26)，非概率可靠性指標(biāo)為

η=8.979 7

(43)

5.2 第二級載荷作用

根據(jù)式(15)，將一級應(yīng)力sp1作用時間t1等效換算成二級應(yīng)力sp2的等效作用時間t21：

t21=[1 111.8,1 899.7]

(44)

本文計算的區(qū)間結(jié)果t21=[1 111.8,1 899.7]，該區(qū)間上下界的幾何平均數(shù)為1 453.3。而文獻(xiàn)[22]模型計算的t21=1 454次。顯然1 454∈[1 111.8,1 899.7]，且1 454與1 453.3非常接近，證明了本文區(qū)間模型的正確性。從區(qū)間寬度來看，參考統(tǒng)計學(xué)區(qū)間估計的計算，若對于樣本均值為1 454、標(biāo)準(zhǔn)差為100的樣本，要求置信度達(dá)到99.87%，則對應(yīng)的區(qū)間估計值為[1 154,1 754]，與本文運(yùn)算的區(qū)間[1 111.8,1 899.7]接近程度在可接受的范圍之內(nèi)，表明本文運(yùn)算的區(qū)間精確度可接受，在工程應(yīng)用中有切實(shí)的指導(dǎo)價值。

顯然，結(jié)構(gòu)受一級應(yīng)力和二級應(yīng)力作用之后，產(chǎn)生的損傷量等效于原結(jié)構(gòu)在二級應(yīng)力sp2下連續(xù)作用時間t21+t2產(chǎn)生的損傷量，即總時間為

(45)

(46)

即

(47)

結(jié)構(gòu)受一級應(yīng)力和二級應(yīng)力作用之后，產(chǎn)生的總損傷量為

(48)

剩余強(qiáng)度為

(49)

(50)

5.3 第三級載荷作用

第三級載荷[264.1,291.9]MPa，T3=[145 971,185 607]，γ3=4.186。根據(jù)式(18)，可得三級應(yīng)力等效作用時間t32：

(51)

本文計算的區(qū)間結(jié)果t32=[3 833.8,916 3.8]，該區(qū)間上下界的幾何平均數(shù)為5 927.24。而文獻(xiàn)[22]模型計算的t21=5 878次。顯然5 878∈[3 833.8,9 163.8]，且5 878與5 927.24非常接近，證明了本文區(qū)間模型的正確性。從區(qū)間寬度來看，參考統(tǒng)計學(xué)區(qū)間估計的計算，若對于樣本均值為5 878、標(biāo)準(zhǔn)差為570的樣本，要求置信度達(dá)到99.98%，則對應(yīng)的區(qū)間估計值為[3 826,7 930]，與本文運(yùn)算的區(qū)間[3 833.8,9 163.8]接近程度在可接受的范圍之內(nèi)，表明本文運(yùn)算的區(qū)間精確度可以接受，在工程應(yīng)用中有切實(shí)的指導(dǎo)價值。

當(dāng)t3=6 000時

(52)

(53)

(54)

結(jié)構(gòu)受一、二級應(yīng)力和三級應(yīng)力作用之后，產(chǎn)生的總損傷量為

(55)

剩余強(qiáng)度為

(56)

(57)

同理，計算出t3=18 000，30 000，42 000，54 000，66 000次循環(huán)直至達(dá)到安全臨界狀態(tài)時，結(jié)構(gòu)受一、二、三級應(yīng)力作用之后的非概率可靠性指標(biāo)。結(jié)構(gòu)受多級載荷應(yīng)力作用下指標(biāo)的變化情況見表1。圖4所示為本文算法與文獻(xiàn)[22]算法的剩余強(qiáng)度值對比。

表1 結(jié)構(gòu)受多級載荷作用下的非概率可靠性指標(biāo)Tab.1 Non-probability reliability index of structure under multi-level loads

圖4 本文算法與文獻(xiàn)[22]算法計算的剩余強(qiáng)度對比Fig.4 Compare residual strength of algorithm in this paper with the literature [22]

由表1和圖4可以看出，本文算法與文獻(xiàn)[22]算法的結(jié)果接近，文獻(xiàn)[22]算法的結(jié)果均處于本文運(yùn)算的區(qū)間內(nèi)，證明了在進(jìn)行多級載荷累積損傷下結(jié)構(gòu)的剩余強(qiáng)度計算中，用區(qū)間數(shù)表達(dá)變量，經(jīng)區(qū)間運(yùn)算以獲得區(qū)間范圍，包含了計算結(jié)果所有可能的取值，這正是區(qū)間數(shù)學(xué)的特點(diǎn)之一。同時還可發(fā)現(xiàn)本文算法所得剩余強(qiáng)度均值略大于文獻(xiàn)[22]算法的結(jié)果，表明文獻(xiàn)[22]算法相對本文區(qū)間算法更偏于保守。

圖5 直接運(yùn)算與等效轉(zhuǎn)化計算壽命比累積和對比Fig.5 Compare accumulation of the life proportion of direct operation with equivalent transformation

圖6 本文計算的非概率動態(tài)可靠性指標(biāo)與文獻(xiàn)[22]計算的可靠度Fig.6 Compare non-probability dynamic reliability index in this paper with degree of reliability in the literature [22]

圖6所示為本文計算的非概率動態(tài)可靠性指標(biāo)與文獻(xiàn)[22]計算的可靠度對比。由表1和圖6可知，本文計算的非概率動態(tài)可靠性指標(biāo)與文獻(xiàn)[22]計算的可靠度，在第三級載荷作用下，均隨著作用時間的推移呈減小趨勢。由于兩者的度量尺度不同(文獻(xiàn)[22]是基于剩余強(qiáng)度滿足正態(tài)分布的假設(shè)，根據(jù)一次二階矩法計算出結(jié)構(gòu)元件的疲勞可靠度)，無法直接比較,但仍然可以對兩個結(jié)果進(jìn)行對比分析，了解不同度量尺度下對結(jié)構(gòu)可靠度的判斷以及疲勞可靠度隨時間推移變化趨勢上的差異。本文計算的非概率動態(tài)可靠性指標(biāo)反映出，在第三級載荷作用時間達(dá)到129 972次循環(huán)時，非概率可靠性指標(biāo)η=1，疲勞試樣件處于安全臨界狀態(tài)；當(dāng)載荷作用時間小于129 972次循環(huán)時，結(jié)構(gòu)處于絕對安全狀態(tài)。由文獻(xiàn)[22]計算的可靠度結(jié)果可知，當(dāng)?shù)谌壿d荷作用時間小于90 000次循環(huán)時，結(jié)構(gòu)的可靠度達(dá)到1.00，結(jié)構(gòu)處于絕對安全狀態(tài)；而當(dāng)載荷作用時間大于90 000次循環(huán)時，結(jié)構(gòu)的可靠度持續(xù)減小，在載荷作用時間達(dá)到129 972次循環(huán)時，可靠度僅為0.85，這也印證了文獻(xiàn)[22]的算法對結(jié)構(gòu)的可靠度判斷偏于保守。此外，本文及文獻(xiàn)[22]的結(jié)果均認(rèn)為，當(dāng)?shù)谌壿d荷作用時間小于90 000次循環(huán)時，結(jié)構(gòu)處于絕對安全狀態(tài)，而當(dāng)載荷作用時間大于90 000次循環(huán)時，疲勞可靠度快速減小。兩者在對結(jié)構(gòu)可靠度的判斷以及疲勞可靠度隨時間推移變化趨勢上具有一致性，說明用本文建立的多級載荷累積損傷下結(jié)構(gòu)的動態(tài)可靠性模型與方法是行之有效的。

6 結(jié)論

(1)不同級別的載荷不同加載順序下對結(jié)構(gòu)的剩余疲勞強(qiáng)度將產(chǎn)生影響。金屬材料在高-低加載時存在“過載效應(yīng)”，最初施加的高載荷導(dǎo)致了材料的疲勞強(qiáng)度減小，加劇了后續(xù)施加低載荷造成的損傷，壽命比累積和將小于1；在低-高順序加載下，由于低應(yīng)力的鍛煉作用，產(chǎn)生低載強(qiáng)化效果，壽命比累積和將大于1?？筛鶕?jù)損傷等效等值壽命原則，在不同次序施加于結(jié)構(gòu)的不同級別應(yīng)力之間，對作用時間進(jìn)行等效轉(zhuǎn)化，采用高-低加載順序時，等效轉(zhuǎn)化運(yùn)算對剩余壽命的預(yù)測將更偏于保守，低-高順序加載則相反。

(2)利用區(qū)間理論，建立了動態(tài)非概率可靠性預(yù)測模型，解決了工程實(shí)際中樣本數(shù)據(jù)不足條件下多級載荷累積損傷作用的結(jié)構(gòu)動態(tài)可靠性預(yù)測問題，對概率與模糊的動態(tài)可靠性模型與方法進(jìn)行了有益的補(bǔ)充。

(3)在結(jié)構(gòu)的剩余強(qiáng)度的預(yù)測中，本文采用區(qū)間數(shù)表達(dá)變量，經(jīng)區(qū)間運(yùn)算以獲得的剩余強(qiáng)度區(qū)間范圍，包含了文獻(xiàn)[22]的計算結(jié)果；同時本文運(yùn)算得到的剩余強(qiáng)度均值略高于文獻(xiàn)[22]的計算結(jié)果。相比文獻(xiàn)[22]，本文的結(jié)果更偏于樂觀。

(4)由于度量尺度不同，無法直接比較本文計算的非概率動態(tài)可靠性指標(biāo)與文獻(xiàn)[22]計算的可靠度，但本文算法認(rèn)為，算例中，當(dāng)載荷作用時間小于129 972次循環(huán)時，結(jié)構(gòu)處于絕對安全狀態(tài)；而文獻(xiàn)[22]算法認(rèn)為在載荷作用時間達(dá)到129 972次循環(huán)時，可靠度僅為0.85。本文算法對結(jié)構(gòu)的可靠性判斷偏于樂觀。

參考文獻(xiàn)：

[1] 王新剛, 張義民, 王寶艷. 機(jī)械零部件的動態(tài)可靠性分析[J]. 兵工學(xué)報,2009,30(11):1510-1514.

WANG Xingang, ZHANG Yimin, WANG Baoyan. Dynamic Reliability Analysis of Mechanical Components[J]. Acta Armamentaril,2009,30(11):1510-1514.

[2] 秦大同,周志剛,楊軍,等. 隨機(jī)風(fēng)載作用下風(fēng)力發(fā)電機(jī)齒輪傳動系統(tǒng)動態(tài)可靠性分析[J].機(jī)械工程學(xué)報,2012,48(3):1-8.

QIN Datong, ZHOU Zhigang, YANG Jun, et al. Time-dependent Reliability Analysis of Gear Transmission System of Wind Turbine under Stochastic Wind Load[J]. Journal of Mechanical Engineering, 2012,48(3):1-8.

[3] 方永鋒, 陳建軍, 馬洪波. 共同隨機(jī)載荷下結(jié)構(gòu)動態(tài)可靠性分析[J]. 應(yīng)用力學(xué)學(xué)報,2016,33(2):345-351.

FANG Yongfeng, CHEN Jianjun, MA Hongbo. Analysis of Structural Dynamic Reliability under Together Stochastic Loads[J]. Chinese Journal of Applied Mechanics,2016,33(2):345-351.

[4] 徐芳, 周志剛. 隨機(jī)風(fēng)作用下風(fēng)力發(fā)電機(jī)齒輪傳動系統(tǒng)動載荷計算及統(tǒng)計分析[J]. 中國機(jī)械工程,2016,27(3):290-295.

XU Fang, ZHOU Zhigang. Dynamic Load Calculation and Statistical Analysis of Gear Transmission System of Wind Turbine under Stochastic Wind Loading[J].China Mechanical Engineering,2016,27(3):290-295.

[5] HAJIALIZADEH D, STEWART M G, ENRIGHT B, et al. Spatial Time-dependent Reliability Analysis of Reinforced Concrete Slab Bridges Subject to Realistic Traffic Loading[J]. Structure & Infrastructure Engineering,2015,12(9):1086385.

[6] 王磊, 李林龍, 馬亞飛,等. 模糊檢測數(shù)據(jù)下服役鋼筋混凝土橋梁的動態(tài)可靠度評估[J]. 防災(zāi)減災(zāi)工程學(xué)報,2014,34(4):529-536.

WANG Lei, LI Linlong, MA Yafei, et al. Dynamic Reliability Assessment of Existing Reinforced Concrete Bridges Based on Fuzzy Inspection Data[J]. Journal of Disaster Prevention & Mitigation Engineering,2014,34(4):529-536.

[7] 高明君, 張國義, 高家一,等. 基于強(qiáng)度退化的機(jī)構(gòu)模糊動態(tài)可靠性分析方法[J]. 強(qiáng)度與環(huán)境,2015(1):54-62.

GAO Mingjun, ZHANG Guoyi, GAO Jiayi, et al. Fuzzy Time-variant Reliability Analysis Method for Mechanisms [J]. Structure & Environment Engineering,2015(1):54-62.

[8] 張宇貽, 秦權(quán). 鋼筋混凝土橋梁構(gòu)件的時變可靠度分析[J]. 清華大學(xué)學(xué)報(自然科學(xué)版),2001,41(12):65-67.

ZHANG Yuyi, QIN Quan. Time-dependent Reliability Analysis of Deteriorating RC Bridge Girders[J]. Journal of Tsinghua University(Natural Science),2001,41(12):65-67.

[9] HU Z, DU X. A Sampling Approach to Extreme Value Distribution for Time-dependent Reliability Analysis[J]. Journal of Mechanical Design,2013,135(7):689-698.

[10] MOURELATOS Z P, MAJCHER M, GEROULAS V. Time-dependent Reliability Analysis of Vibratory Systems with Random Parameters[J]. Journal of Vibration & Acoustics,2016,138(3):1-9.

[11] 姜潮,黃新萍,韓旭,等. 含區(qū)間不確定性的結(jié)構(gòu)時變可靠度分析方法[J]. 機(jī)械工程學(xué)報,2013,49(10):186-193.

JIANG Chao, HUANG Xinping, HAN Xu, et al. Time-dependent Structural Reliability Analysis Method with Interval Uncertainty[J]. Journal of Mechanical Engineering,2013,49(10):186-193.

[12] 魏宗平,李天恩. 基于強(qiáng)度退化分析的結(jié)構(gòu)非概率動態(tài)可靠性模型[J]. 機(jī)械科學(xué)與技術(shù),2011,30(8):1397-1401.

WEI Zongping, LI Tianen. Non-probabilistic Time-dependent Reliability Model of a Structure Based on Strength Degradation Analysis [J]. Journal of Mechanical Science and Technology,2011,30(8):1397-1401.

[13] 方永鋒,陳建軍,閻彬,等. 結(jié)構(gòu)動態(tài)非概率可靠性預(yù)測模型[J]. 西安電子科技大學(xué)學(xué)報,2012,39(6):170-175.

FANG Yongfeng, CHEN Jianjun, YAN Bin, et al. Model for Prediction of Structural Dynamic Non-probabilistic Reliability[J]. Journal of Xidian University,2012,39(6):170-175.

[14] 張德權(quán), 韓旭, 姜潮,等. 時變可靠性的區(qū)間PHI2分析方法[J]. 中國科學(xué):物理學(xué)力學(xué)天文學(xué),2015,45(5):054601-13.

ZHANG Dequan, HAN Xu, JIANG Chao, et al. The Interval PHI2 Analysis Method for Time-dependent Reliability[J]. Scientia Sinica: Physica, Mechanica & Astronomica,2015,45(5):054601-13.

[15] 方永鋒,陳建軍,閻彬,等. 多次區(qū)間載荷下結(jié)構(gòu)混合動態(tài)可靠性的計算[J]. 華南理工大學(xué)學(xué)報(自然科學(xué)版),2013,41(11):91-95.

FANG Yongfeng, CHEN Jianjun, YAN Bin, et al. Calculation of Structural Hybrid Dynamic Reliability under Multiple Interval Loads [J]. Journal of South China University of Technology (Natural Science Edition),2013,41(11):91-95.

[16] 王正,謝里陽,李兵,等. 共因失效系統(tǒng)動態(tài)可靠性模型[J]. 中國機(jī)械工程,2008,19(1):5-9.

WANG Zheng, XIE Liyang, LI Bing, et al. Time-dependent Reliability Model of System with Common Cause Failure[J]. China Mechanical Engineering,2008,19(1):5-9.

[17] 蘇志霄,劉宏昭,王建平,等. 基于剩余強(qiáng)度退化規(guī)律的疲勞損傷非線性演化模型[J].機(jī)械強(qiáng)度,2000,22(3): 238-240.

SU Zhixiao, LIU Hongzhao, WANG Jianping , et al. Non-linear Model of Fatigue Damage Based on the Residual Strength Degradation Law [J]. Journal of Mechanical Strength,2000,22(3):238-240.

[18] 姚衛(wèi)星. 結(jié)構(gòu)疲勞壽命分析 [M].北京:國防工業(yè)出版社,2003.

YAO Weixing. The Statistical Analysis of Structure Fatigue Life [M].Beijing: National Defence Industry Press,2003.

[19] 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):158-181.

[20] 儲軍,鄭松林. 譜載荷下加載次序效應(yīng)定量評價的初步研究[J]. 中國機(jī)械工程,2015,26(3):365-368.

CHU Jun, ZHENG Songlin. Preliminary Study of Quantitative Evaluation for Spectrum Load Sequence Effect[J]. China Mechanical Engineering,2015,26(3):365-368.

[21] 郭書祥,呂震宙,馮元生. 基于區(qū)間分析的結(jié)構(gòu)非概率可靠性模型[J]. 計算力學(xué)學(xué)報,2001,18(1):56-60.

GUO Shuxiang, LYU Zhenzhou, FENG Yuan-sheng. A Non-probabilistic Model of Structural Reliability Based on Interval Analysis[J]. Chinese Journal of Computational Mechanics,2001,18(1):56-60.

[22] 李少宏,陳建軍,方永鋒.多級載荷作用下剩余強(qiáng)度模型的結(jié)構(gòu)元件疲勞可靠性分析[J].機(jī)械科學(xué)與技術(shù),2013,32(6):791-795.

LI Shaohong, CHEN Jianjun, FANG Yongfeng. Analysis of the Fatigue Reliability for Structural Element under Multi-level Loads Based on the Residual Strength Model[J]. Mechanical Science & Technology for Aerospace Engineering,2013,32(6):791-795.

[23] 編委會.最新資產(chǎn)評估常用數(shù)據(jù)與參數(shù)手冊2010～2011年[M].北京:經(jīng)濟(jì)科學(xué)出版社, 2010.

Editorial Board. The Latest Commonly Used Data and Parameters Manual of Asset Evaluation for 2010-2011 [M]. Beijing: Economic Science Press,2010.