滾動軸承壽命的統(tǒng)計模擬

姜英碩，姜祺，周佳新，汪東華

(1.沈陽建筑大學(xué)，沈陽 110168；2.清華大學(xué) 工業(yè)工程系，北京 100084；3.中交煤氣設(shè)計院有限公司， 沈陽 110026)

1 前言

拖拉機等工程運輸車輛上使用的滾動軸承(以下簡稱軸承)通常有幾套、幾十套，甚至多達(dá)上百套，構(gòu)成了軸承子系統(tǒng)。工作過程中任何一套軸承發(fā)生疲勞、磨損或永久變形等故障，均會導(dǎo)致軸承子系統(tǒng)乃至整臺主機發(fā)生故障。因此，預(yù)測軸承可靠性指標(biāo)，及時采取可能的措施提高軸承使用壽命是十分必要的。

目前軸承選型和壽命計算普遍采用ISO 281—2007和GB/T 6391—2010確定的滾動軸承額定壽命計算公式[1],經(jīng)變換公式形式可得軸承額定壽命計算式為

(1)

式中：C為軸承額定動載荷，N;P為軸承當(dāng)量動載荷，N；n為軸承工作轉(zhuǎn)速，r/min；ε為壽命指數(shù)。

用(1)式可求出軸承可靠壽命L10的值。但應(yīng)當(dāng)注意的是，在拖拉機等工程運輸車輛上，軸承載荷和轉(zhuǎn)速實際上都是變量，變化幅度較大，且具有隨機性。

2 軸承載荷和轉(zhuǎn)速的分布

現(xiàn)場試驗和工作實際表明，拖拉機等工程運輸車輛上的軸承載荷為隨機變量，載荷峰值服從正態(tài)分布[2]。根據(jù)集材用拖拉機模擬工況下的軸承“載荷-時間歷程”樣本記錄，用峰值計算法求得軸承載荷頻譜[3-4]，結(jié)果見表1。

表1 集材用拖拉機軸承載荷頻譜

表2 及的正態(tài)分布參數(shù)

3 軸承壽命的統(tǒng)計模擬

(1)用數(shù)學(xué)方法產(chǎn)生統(tǒng)計獨立的(0，1)均勻分布的隨機數(shù)序列{ri}N。

(2)用標(biāo)準(zhǔn)正態(tài)分布隨機數(shù)抽樣式，成對產(chǎn)生標(biāo)準(zhǔn)正態(tài)分布隨機數(shù)序列{Zi}N,其抽樣式為

(2)

式中：i=1，3，5，…，N。對于標(biāo)準(zhǔn)正態(tài)分布的分布參數(shù)，均值μ=0，均方差σ=1。

(3)利用線性變換，對軸承隨機載荷和隨機轉(zhuǎn)速進(jìn)行模擬，產(chǎn)生正態(tài)分布的隨機數(shù)序列{Pi}N，{ni}N，其模擬式為

(3)

(4)

{Zi}NSP)ε]，

(5)

軸承參數(shù)主要有：型號為N210(圓柱滾子軸承)，額定動載荷C=4.2×104N，依據(jù)GB/T 6391—2010取壽命指數(shù)ε=3.33。為保證模擬精度取樣本容量N=2.2×103次，軸承載荷和轉(zhuǎn)速的均值、均方差見表2。將各已知常數(shù)、隨機變量和已知參數(shù)代入(5)式，通過編程計算，可得軸承壽命的隨機數(shù)序列{Lhi}N。對壽命樣本進(jìn)行統(tǒng)計處理，并用χ2檢驗法對軸承壽命分布進(jìn)行擬合度檢驗,結(jié)果見表3。

(6)

(7)

表3 軸承壽命分布的χ2檢驗結(jié)果

(5)用可靠性特征量計算公式[8]求取軸承可靠度R(L)、故障率λ(L)、可靠壽命L(R)及中值壽命L50等指標(biāo)。

R(L)=1-Φ(Z)=

(8)

式中：Φ(Z)為標(biāo)準(zhǔn)正態(tài)分布函數(shù)，可查標(biāo)準(zhǔn)正態(tài)分布表。

(9)

(10)

式中：ZP為標(biāo)準(zhǔn)正態(tài)分布分位點，可查標(biāo)準(zhǔn)正態(tài)分布表。

(11)

圖1 軸承R(L)-Lh曲線

圖2 軸承λ(L)-Lh曲線

表4 滾動軸承可靠性特征量

4 軸承統(tǒng)計模擬結(jié)果及分析

5 結(jié)論

(2)如果軸承隨機載荷和隨機轉(zhuǎn)速服從正態(tài)分布，軸承壽命將服從對數(shù)正態(tài)分布。

(3)用統(tǒng)計模擬方法求取軸承壽命是可行的。模擬結(jié)果可作為預(yù)測軸承可靠度等指標(biāo)和備件數(shù)量的依據(jù)。