王桂金

(原鋼鐵研究總院，北京 100083)

1 疲勞壽命的Weibull分布

服從Weibull分布的疲勞壽命在可靠性理論中占有重要地位[1-2]。與其相關(guān)的實(shí)際疲勞壽命數(shù)據(jù)應(yīng)該是一組滿足Weibull分布的隨機(jī)變量,并應(yīng)給出Weibull分布定義的統(tǒng)計(jì)量。帶量綱的統(tǒng)計(jì)量有均值即期望值μ, 方差σ2, 中值median和峰值mode。無量綱統(tǒng)計(jì)量有標(biāo)準(zhǔn)誤差σ/μ，斜度γ1和過盈峭度γ2。對于具有形狀參數(shù)κ和尺寸參數(shù)λ的二參數(shù)Weibull分布, 其計(jì)算公式為[1]

中值median=λ(ln2)1/κ,

式中：Γi=Γ(1+i/κ)，為Gamma函數(shù)；κ為形狀參數(shù)；λ為尺寸參數(shù)。

值得注意的是，無量綱統(tǒng)計(jì)量都是形狀參數(shù)κ的函數(shù), 而與尺寸參數(shù)λ無關(guān)。 因此只要Weibull分布的形狀參數(shù)相同,其標(biāo)準(zhǔn)誤差、斜度和過盈峭度亦相同。通常研究者從試驗(yàn)數(shù)據(jù)出發(fā)，用極大似然法[1]或最大線性不變法[3]求出κ和λ,并計(jì)算擬合的可靠度。 由于在滿足最小偏差的條件下，壽命分布形狀不一定正確， 因此有必要從試驗(yàn)數(shù)據(jù)算出和分布形狀有關(guān)的斜度γ1及過盈峭度γ2，并計(jì)算出對應(yīng)的κ(γ1)及κ(γ2)，再同擬合的κ值比較。為了避免實(shí)測壽命的系統(tǒng)誤差和人為誤差，可以首先隨機(jī)產(chǎn)生一定數(shù)量的Weibull壽命(例如100個(gè))，按極大似然法計(jì)算κ和λ，并和設(shè)定值比較。然后由隨機(jī)壽命數(shù)據(jù)按常規(guī)方法計(jì)算斜度和過盈峭度，進(jìn)而推算相應(yīng)的Weibull形狀參數(shù)κ(γ1)及κ(γ2) ，并和極大似然法的形狀參數(shù)κ相對照。最后還應(yīng)該對實(shí)測壽命數(shù)據(jù)做同樣的處理，以便對實(shí)際壽命數(shù)據(jù)是否為Weibull隨機(jī)分布做出判斷。

2 隨機(jī)壽命和極大似然法擬合[1]

Weibull分布的隨機(jī)壽命L為

L=λ(-lna)1/κ,

(1)

式中：a為在[0,1]區(qū)間內(nèi)均勻分布的隨機(jī)變量，為易于和鉻鋼軸承試驗(yàn)數(shù)據(jù)相比較，兩者皆取為1。從文獻(xiàn)[4]的附表中依次選取100個(gè)隨機(jī)變量ai,i=1, 2, …,100，按(1)式算出的隨機(jī)壽命Li見表1。根據(jù)下述極大似然法公式可求出這組數(shù)據(jù)的Weibull分布參數(shù)κ和λ

表1 Weibull隨機(jī)壽命Li(i=1，2，…，100，κ=1，λ=1)

(2)

(3)

式中：N為所取樣本大小，可根據(jù)需要選擇，本例為10≤N≤100；Li為第i個(gè)試樣的隨機(jī)壽命。首先經(jīng)對κ迭代計(jì)算使(3)式接近于零(小于|1×10-5|)得到κ(N)，然后將其代入(2)式得到λ(N)。另外，按統(tǒng)計(jì)公式計(jì)算隨機(jī)壽命分布的斜度γ1(N)和過盈峭度γ2(N)，并計(jì)算出相應(yīng)的形狀參數(shù)κ(γ1,N) 和κ(γ2,N)。

3 計(jì)算結(jié)果

3.1 隨機(jī)壽命分布的λ(N)，κ(N)，κ(γ1, N)和κ(γ2, N)

為檢查隨機(jī)壽命數(shù)據(jù)組的隨機(jī)特性，計(jì)算了N=10～100 的λ(N)，κ(N)，κ(γ1,N)和κ(γ2,N)，結(jié)果如圖1所示(橫坐標(biāo)為樣本大小N，下同)。當(dāng)采用全部隨機(jī)壽命數(shù)據(jù)，即N=100時(shí)，λ(100)=1.016 57，κ(100)=0.996 07，κ(γ1,100)和κ(γ2,100)分別為1.028 910和1.083 690，都很接近設(shè)定值1。結(jié)果表明這100個(gè)隨機(jī)壽命符合κ=1和λ=1的Weibull分布，不僅偏差最小，而且分布具有κ=1的不對稱性。極大似然法算出λ和κ在整個(gè)N范圍變動不大，只在N<20時(shí)變動幅度略大。如果只取部分隨機(jī)壽命計(jì)算斜度和過盈峭度及其形狀參數(shù)和尺寸參數(shù)，則在N≥48時(shí)仍相當(dāng)接近1。但是當(dāng)N<48時(shí)，斜度和過盈峭度的形狀參數(shù)顯著增大，并且波動幅度大。 查表1可知，第48號是數(shù)據(jù)組的最大值 5.632 390，應(yīng)該是它的出現(xiàn)終止了急劇的波動。因此，當(dāng)隨機(jī)壽命樣本少于48時(shí)，雖然極大似然法能給出接近1的形狀參數(shù)，但是分布不對稱性未能達(dá)到κ=1的Weibull分布理論值。直到第48號數(shù)據(jù)加入，3種形狀因子才得以一致。 其原因是極大似然法不涉及三、四次矩，對分布的形狀不很敏感。

圖1 隨機(jī)壽命試樣的極大似然法和由數(shù)據(jù)的斜度和過盈峭度得出的形狀參數(shù)

3.2 隨機(jī)壽命由小到大排列的λ(N)，κ(N)，κ(γ1,N)和κ(γ2, N)

通常實(shí)測壽命是由小到大按次序記錄，為便于比較，把隨機(jī)壽命由小到大排列，重新計(jì)算N=10～100的λ(N)，κ(N)及κ(γ1,N)，κ(γ2,N) ，結(jié)果如圖2所示，可以看出：

圖2 由小到大排列的隨機(jī)壽命的極大似然法和直接由斜度及過盈峭度得出的形狀參數(shù)

(1)與圖 1相同， 全樣本的隨機(jī)壽命(即N=100)，由極大似然法和由數(shù)據(jù)的斜度和過盈峭度等3種方法得出的λ值和各κ值均接近1 ，這是合理的；

(2)在10≤N≤100時(shí)，極大似然法得出的κ(N)由1逐步增大到1.671 50(N=21)，隨后稍微減小，再接著增大到2.822 30(N=10)；

(3)從N=100到N=10由斜度推算的κ(γ1,N) 增長較快并有幾個(gè)峰，例如在N=21達(dá)到峰值3.897 70, 在N=12處取極小值1.781 70，幾乎和極大似然法的κ值1.743 30重合。這是極大似然法的κ(N)和由斜度推算的κ(γ1,N)第2次重合；

(4)由數(shù)據(jù)組的過盈峭度換算的κ(γ2,N)值，在所討論的N范圍相當(dāng)離散。為清晰起見，在圖上只標(biāo)出κ(γ2,N)小于4.5 的值(其他圖也作類似處理)，不過仍然可以看到，在N=12,N=99和N=100處，κ(γ2,N)比較接近相應(yīng)的κ(N)；

(5)極大似然法的尺寸參數(shù)λ(N)從1.016 70(N=100)平穩(wěn)地減小到0.056 83(N=10)。

4 試驗(yàn)分析

4.1 1208K+H208軸承的疲勞壽命[5](37個(gè)試樣)

試驗(yàn)條件：載荷5 kN，轉(zhuǎn)速 3 000 r/min，實(shí)測疲勞壽命見表2，計(jì)算結(jié)果如圖3所示。圖中：λ(37)=1 387 h，κ(37)= 0.955 50，κ(γ1, 37)=1.101 04，κ(γ2, 37)=1.149 38。

表2 H208軸承實(shí)測疲勞壽命(37個(gè)試樣) h

圖3 1208K+H208軸承實(shí)測疲勞壽命按極大似然法和按斜度及過盈峭度得出的λ(N) ，к(N)，к(γ1, N)和к(γ2, N)

4.2 7208軸承的疲勞壽命[6](60個(gè)試樣)

試驗(yàn)條件：徑向載荷為8.281 kN, 軸向載荷為4.673 kN, 轉(zhuǎn)速 4 500 r/min，結(jié)果如圖4所示，λ(60)=267.12 h，κ(60)=1.024 43，κ(γ1, 60)=0.926 16，κ(γ2, 60)=0.931 83。

圖4 7208軸承實(shí)測疲勞壽命按極大似然法和按斜度及過盈峭度得出的λ(N) ，κ(N)，κ(γ1, N) 和κ(γ2, N)

4.3 6104軸承(50個(gè)試樣)和6307軸承(20個(gè)試樣)的疲勞壽命[7]

6104軸承試驗(yàn)條件：徑向載荷 3.2 kN, 轉(zhuǎn)速 2 800 r/min，結(jié)果如圖5所示，λ(50)=125.28×106r，κ(50)=1.027 600，κ(γ1, 50)=1.336 360，κ(γ2, 50)=1.595 200。

圖5 6104軸承實(shí)測疲勞壽命按極大似然法和按斜度及過盈峭度得出的λ(N) ，κ(N)，κ(γ1, N) 和κ(γ2, N)

6307軸承試驗(yàn)條件：徑向載荷 7.8 kN, 轉(zhuǎn)速 1 900 r/min，結(jié)果如圖6所示。λ(20)=175.2×106r，κ(20)=1.577 70，κ(γ1, 20)=2.175。但κ(γ2，20)在κ=1～10時(shí)無解，因?yàn)棣?(20) =-0.725 960，低于γ2最小值-0.289(κ=3.42)。

圖6 6307軸承實(shí)測疲勞壽命分布按極大似然法和直接由斜度及過盈峭度得出的λ(N) ，κ(N)，κ(γ1, N) 和κ(γ2, N)

4.4 結(jié)果分析

總的看來，1208K+H208和7208軸承的實(shí)測疲勞壽命按極大似然法和直接由斜度及過盈峭度得出的各形狀參數(shù)κ值在全樣本處幾乎重疊，與圖2的結(jié)果很相近，因此其實(shí)測疲勞壽命分布十分接近隨機(jī)分布。而6104和6307軸承的3種形狀參數(shù)κ都還未完成相互靠攏的過程，可能需要更大的樣本或者更嚴(yán)格控制樣品制作和試驗(yàn)條件作進(jìn)一步驗(yàn)證。因此認(rèn)為，極大似然法、斜度和過盈峭度的形狀參數(shù)κ相互重疊的出現(xiàn)，應(yīng)是滿足Weibull壽命隨機(jī)特性的一個(gè)條件，具體而言，一組符合Weibull分布的疲勞壽命數(shù)據(jù)應(yīng)有如下特征。

(1)當(dāng)全樣本的試樣數(shù)N足夠大，應(yīng)在一個(gè)或一個(gè)以上的N值處，實(shí)測疲勞壽命按極大似然法和由斜度及過盈峭度得出的κ比較接近甚至重疊。這樣，疲勞壽命分布和理想的Weibull分布既有最小的偏差， 而且分布的形狀(即斜度和峭度)也一致。 如果有多個(gè)重合點(diǎn)，建議采用最大N值的形狀參數(shù)κ和尺寸參數(shù)λ。

(2)隨著N值增大，由于長壽命數(shù)據(jù)進(jìn)入擬合計(jì)算，Weibull壽命分布的3種形狀參數(shù)κ都逐漸變小。而最大似然法的尺寸參數(shù)則幾乎呈線性增大。

(3)同一材料在同一試驗(yàn)條件下，如果有些試樣因故不能試驗(yàn)到失效而提前舍去，可能會得出偏大的形狀參數(shù)和偏小的尺寸參數(shù)，從而導(dǎo)致L10的偏差。

5 結(jié)束語

一般認(rèn)為，最大似然法在大樣本試驗(yàn)中比現(xiàn)有其他方法優(yōu)越[1]。本研究表明，如果能輔之以數(shù)據(jù)的斜度和過盈峭度，取得有關(guān)壽命分布形狀的信息，可以把Weibull分布參數(shù)評估得更合理。