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

2014-07-24 05:01:04姜英碩姜祺周佳新汪東華

軸承 2014年9期

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

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

1 前言

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

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

(1)

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

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

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

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

表1 集材用拖拉機(jī)軸承載荷頻譜

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

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

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

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

(2)

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

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

(3)

(4)

{Zi}NSP)ε]，

(5)

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