張海燕，俞曉陽

（甘肅建筑職業(yè)技術(shù)學(xué)院 環(huán)境與市政工程系，甘肅 蘭州 730050）

0 引言

設(shè)備的故障檢測(cè)和定期維修是機(jī)械設(shè)備長期可靠、穩(wěn)定運(yùn)行必不可少的一個(gè)環(huán)節(jié)，與此相比，對(duì)設(shè)備故障的提前預(yù)測(cè)可以將對(duì)設(shè)備的維護(hù)和保養(yǎng)從被動(dòng)變?yōu)橹鲃?dòng)，從而提高設(shè)備的可靠性。目前，故障預(yù)測(cè)和健康管理技術(shù)已經(jīng)成為現(xiàn)代工業(yè)設(shè)備維修和保養(yǎng)以外提高可靠性的有效方式。預(yù)測(cè)和健康管理可以提高設(shè)備可靠性的前景是廣闊的，因此越來越受到人們的歡迎和使用。故障預(yù)測(cè)是健康管理的主要目的之一，根據(jù)從傳感器傳回的信息可以對(duì)設(shè)備受環(huán)境影響的累計(jì)效應(yīng)來預(yù)測(cè)其壽命的剩余部分。故障預(yù)測(cè)與健康管理通常是使用各類傳感器來采集系統(tǒng)中的時(shí)事數(shù)據(jù)，在計(jì)算機(jī)中通過一些算法對(duì)數(shù)據(jù)進(jìn)行處理以便對(duì)系統(tǒng)的健康狀態(tài)進(jìn)行預(yù)測(cè)和評(píng)估，在系統(tǒng)發(fā)生故障前能對(duì)該故障進(jìn)行可靠的預(yù)測(cè)可以提供給決策者一系列的維護(hù)和操作的建議，來避免嚴(yán)重故障帶來的損失。

從以往的經(jīng)驗(yàn)來看，機(jī)械設(shè)備的故障往往會(huì)造成很多大型機(jī)械系統(tǒng)的可靠性下降，帶來很多不可預(yù)測(cè)的后果。在一些高精度的數(shù)控機(jī)床、大型風(fēng)力發(fā)電設(shè)備、精密航空發(fā)動(dòng)機(jī)等領(lǐng)域，都對(duì)設(shè)備的精度、效能提出越來越高的要求，而這些需求都可以依靠預(yù)測(cè)故障來得到滿足。在規(guī)模龐大的機(jī)械設(shè)備系統(tǒng)中，機(jī)械元件長時(shí)間的使用累計(jì)了海量的數(shù)據(jù)，這些數(shù)據(jù)包含了多個(gè)傳感器對(duì)不同類型數(shù)據(jù)的采集信息，故障預(yù)測(cè)依托海量數(shù)據(jù)，可以為機(jī)械設(shè)備的可靠運(yùn)行帶來很大的提升。故障分析中處理的龐大數(shù)據(jù)已經(jīng)逐步成為分析和揭示機(jī)械設(shè)備故障原因及其本質(zhì)過程的重要資源。隨之而來的問題是如何對(duì)機(jī)械設(shè)備進(jìn)行精確的分析和預(yù)測(cè)，這由于機(jī)械設(shè)備本身的原理復(fù)雜性和其所屬環(huán)境及工況的不同都會(huì)成為機(jī)械設(shè)備故障預(yù)測(cè)及處理困難，如何來解決這些問題是現(xiàn)在故障分析領(lǐng)域遇到的一大難點(diǎn)。

在對(duì)保護(hù)軸承的研究過程中，最先提出的是將單層普通滾動(dòng)軸承作為保護(hù)軸承，之后也有學(xué)者提出在普通單層滾動(dòng)軸承外圈安裝阻尼器、雙層滾動(dòng)軸承等，其共同的特點(diǎn)是轉(zhuǎn)子與保護(hù)軸承之間存在固定的保護(hù)間隙。主動(dòng)磁軸承越來越多地應(yīng)用于壓縮機(jī)、渦輪機(jī)和發(fā)電機(jī)等工業(yè)機(jī)器中，因?yàn)樗鼈冎粫?huì)產(chǎn)生輕微的摩擦損失，無需潤滑，并且可以自適應(yīng)操作以優(yōu)化機(jī)器的可靠性和性能。盡管先進(jìn)的控制算法很大程度上提高了主動(dòng)磁軸承可靠性，但考慮到電源故障和軸承過載的故障情況，仍然需要保護(hù)軸承[1-3]。

許多研究成果提供了轉(zhuǎn)子動(dòng)力學(xué)系統(tǒng)在保護(hù)軸承上發(fā)生跌落事件后產(chǎn)生相應(yīng)的模擬結(jié)果。如Ishii和Kirk[4]，將保護(hù)軸承建模為線性彈簧、阻尼器和接觸摩擦，并試圖基于該模型優(yōu)化保護(hù)軸承性能。Sun 等[5-6]提出了一種具有熱增長的非線性滾珠軸承模型，為保護(hù)軸承提供了更精確的部件模型。

在有關(guān)保護(hù)軸承的壽命預(yù)測(cè)的文獻(xiàn)中，很少出現(xiàn)保護(hù)軸承失效的情況。一些標(biāo)準(zhǔn)，如API[7]規(guī)定了可接受的最小跌落事件數(shù)，但除了Sun 外，幾乎沒有關(guān)于保護(hù)軸承壽命的研究報(bào)告。本次研究使用Lundberg Palmgren 公式確定疲勞壽命，該公式僅適用于穩(wěn)定連續(xù)負(fù)載運(yùn)行。本次研究提出的雨流方法對(duì)隨機(jī)載荷有效，并且包括轉(zhuǎn)子和內(nèi)圈之間摩擦產(chǎn)生的剪切應(yīng)力的影響，該因素的加入，完善了文獻(xiàn)[8]對(duì)于保護(hù)軸承壽命的預(yù)測(cè)結(jié)果。

疲勞累計(jì)損傷的概念最早在1924 年就被提出，隨著人們對(duì)于工業(yè)機(jī)械設(shè)備的關(guān)注日益密切，研究人員提出了越來越多的疲勞累計(jì)損傷的理論，如Miner法則、修正Miner 法則、Corten Dolan 準(zhǔn)則和Manson雙線性準(zhǔn)則等[9]。Miner 法則由于其將疲勞損傷機(jī)理進(jìn)行了相當(dāng)?shù)暮喕?，所以這種方法在工程中被廣泛采納，但由于其將不同等級(jí)的載荷作用孤立起來，忽略了他們之間的相互作用，從而使估計(jì)的疲勞壽命和實(shí)際疲勞壽命之間有較大的誤差。疲勞累計(jì)損傷理論主要是以變幅疲勞載荷作用為研究對(duì)象，分析其在長時(shí)間作用的過程中由于疲勞得累計(jì)作用和破壞準(zhǔn)則下，累計(jì)損傷理論認(rèn)為每次循環(huán)載荷的作用會(huì)造成承載機(jī)構(gòu)的一定的損傷累計(jì)。現(xiàn)有的疲勞累計(jì)統(tǒng)計(jì)分析方法有：線性損傷累計(jì)理論、修正的線性損傷累積理論、非線性損傷累積理論以及概率累積損傷理論[10]。其中線性疲勞損傷累計(jì)損傷累計(jì)理論中計(jì)算疲勞破壞的方法一般是：通過統(tǒng)計(jì)循環(huán)載荷下每次累計(jì)的疲勞損傷線性相加值，并且認(rèn)為各應(yīng)力之間相對(duì)獨(dú)立沒有相互作用，一旦當(dāng)累計(jì)的疲勞損傷超過臨界值，則認(rèn)為元件會(huì)發(fā)生疲勞破壞。

本次研究的仿真模型包括主動(dòng)磁軸承、保護(hù)軸承和用有限元建模的水平轉(zhuǎn)子。接觸載荷、赫茲應(yīng)力、亞表面剪切應(yīng)力和滾動(dòng)軸承部件中的熱增長是使用非線性滾珠軸承模型計(jì)算的。應(yīng)用雨流計(jì)數(shù)法將考慮熱效應(yīng)影響的次表面剪切應(yīng)力的作用進(jìn)行累計(jì)，以預(yù)測(cè)保護(hù)軸承的疲勞壽命。次方法也適用于承受隨機(jī)載荷的任何滾動(dòng)元件軸承，例如風(fēng)力渦輪機(jī)傳動(dòng)系軸承。此外，還提供了參數(shù)研究，以確定摩擦系數(shù)、氣隙距離、轉(zhuǎn)子速度和靜態(tài)側(cè)載荷對(duì)保護(hù)軸承壽命的影響。

1 仿真模型

1.1 非線性滾珠軸承模型

考慮熱增長的非線性滾珠軸承模型排除了傾斜變形。軸承組件響應(yīng)于外力F 在圖1 所示的x、y和z方向上偏轉(zhuǎn)。內(nèi)圈、滾珠和外圈表示如圖1 所示。

圖1 非線性球軸承模型

內(nèi)圈、滾珠和外圈的有限元模型如下：

其中，Qi，e是接觸載荷，α是接觸角。內(nèi)圈、滾珠和外圈的模型如圖2。外圈插入軸承箱，軸承箱由沿著軸向被約束，剛度Ks為阻尼為Gs的支座支撐。此外，考慮熱膨脹計(jì)算接觸力。

圖2 轉(zhuǎn)子跌落仿真模型

1.2 轉(zhuǎn)子跌落模型

圖2 顯示了轉(zhuǎn)子下落和深溝球軸承、保護(hù)軸承模型。參考坐標(biāo)系固定在機(jī)械框架上。轉(zhuǎn)子和軸承內(nèi)圈的幾何中心分別為Qb和Qr。（x，y）是Qr的坐標(biāo)，（xb，yb）是固定參照系中Qb的坐標(biāo)。旋轉(zhuǎn)轉(zhuǎn)子和內(nèi)圈之間接觸點(diǎn)處的法向力為：

式中，δ為轉(zhuǎn)子和內(nèi)圈之間的距離。接觸系數(shù)Kc取決于材料類型和接觸幾何形式。cr是保護(hù)軸承徑向間隙，α值在0.08 和0.2 之間，p取決于接觸類型。

2 累計(jì)疲勞損傷的方法

本次研究中對(duì)于保護(hù)軸承的壽命預(yù)測(cè)以力學(xué)分析和數(shù)理統(tǒng)計(jì)相結(jié)合的方法為手段?，F(xiàn)有的關(guān)于隨機(jī)載荷的分析方法及理論有很多，雨流計(jì)數(shù)法、范圍對(duì)法、單參數(shù)計(jì)數(shù)法等。其中雨流計(jì)數(shù)法是在將動(dòng)強(qiáng)度和靜強(qiáng)度兩個(gè)變量統(tǒng)計(jì)到模型中進(jìn)行分析的一種以雙參數(shù)法為基礎(chǔ)的廣泛應(yīng)用的載荷分析方法。這種方法由于具有和疲勞載荷本身固有特性相適應(yīng)的特征，因此在預(yù)測(cè)疲勞載荷造成的損傷中應(yīng)用廣泛。

現(xiàn)有的對(duì)于設(shè)備壽命預(yù)測(cè)的方法大多數(shù)方法有兩個(gè)方向，一是依托于對(duì)設(shè)備檢測(cè)的各類數(shù)據(jù)，二是根據(jù)設(shè)備的物理模型進(jìn)行數(shù)學(xué)分析。其中：依托于數(shù)據(jù)的預(yù)測(cè)方法通常是通過大數(shù)據(jù)處理算法對(duì)狀態(tài)數(shù)據(jù)進(jìn)行處理，從而估算設(shè)備的疲勞壽命。根據(jù)設(shè)備的物理模型分析的方法則是通過建立對(duì)象的數(shù)學(xué)模型，對(duì)采集的數(shù)據(jù)進(jìn)行分析及計(jì)算，從而估算元件的疲勞壽命。本次研究的雨流計(jì)數(shù)法，是由英國工程師M.Matsuishi 和T.Endo 提出的一種將實(shí)際測(cè)量的設(shè)備所受載荷等效為載荷的循環(huán)作用來進(jìn)行元件的疲勞壽命的估計(jì)以及疲勞實(shí)驗(yàn)載荷譜編制的數(shù)學(xué)方法。具體過程為：將一段作用于設(shè)備的應(yīng)力曲線作為分析對(duì)象，將應(yīng)力作為縱坐標(biāo)，時(shí)間作為橫坐標(biāo)，旋轉(zhuǎn)90°，數(shù)據(jù)記錄點(diǎn)和屋檐上的雨點(diǎn)非常相似，而載荷在曲線上的作用過程就像從屋檐內(nèi)順流向下的雨滴一樣，直至最底層。

2.1 雨流循環(huán)計(jì)數(shù)方法

當(dāng)轉(zhuǎn)子跌落事件發(fā)生時(shí)，保護(hù)軸承受到可變應(yīng)力振幅的影響。雨水循環(huán)計(jì)數(shù)方法用于預(yù)測(cè)保護(hù)軸承的疲勞壽命。使用循環(huán)范圍直方圖和Miner 規(guī)則確定累積損傷D 和失效循環(huán)次數(shù)N。如果出現(xiàn)以下情況，預(yù)計(jì)會(huì)發(fā)生故障：

其中ni是統(tǒng)計(jì)的循環(huán)數(shù)Ni是在一定應(yīng)力振幅τi下失效的循環(huán)次數(shù)，在本研究中，公式（5）中D的臨界累積損傷值選擇為1。疲勞壽命L表示為；

2.2 S-N 曲線

N.Raje 和F.Sadeghi 表明，滾動(dòng)疲勞與扭轉(zhuǎn)疲勞相似。他們應(yīng)用扭轉(zhuǎn)疲勞的S-N 曲線來計(jì)算軸承的疲勞壽命，包括熱效應(yīng)的S-N 曲線表示為：

其中B是扭轉(zhuǎn)S-N 曲線的斜率。軸承鋼AISI-52100 的這些參數(shù)列于表1 中。

表1 軸承鋼AISI-52100 參數(shù)

2.3 作用在保護(hù)軸承上的剪切應(yīng)力

Palmgren 和Lundberg 表明，剪切應(yīng)力的振幅與赫茲應(yīng)力和橢圓率有關(guān)。在大多數(shù)滾動(dòng)軸承應(yīng)用中，表面剪切應(yīng)力與法向應(yīng)力相比可以忽略不計(jì)。然而，對(duì)于預(yù)測(cè)滾動(dòng)軸承的疲勞壽命，表面剪切應(yīng)力是不可忽視的，是決定滾動(dòng)軸承耐久性的最重要因素。

3 仿真結(jié)果

該系統(tǒng)由兩個(gè)主動(dòng)磁軸承、兩個(gè)保護(hù)軸承和水平轉(zhuǎn)子組成。圖3 顯示了轉(zhuǎn)子和有限元模型的幾何信息示意圖。本研究中使用的軸承為6016 止動(dòng)軸承。支撐結(jié)構(gòu)的剛度和阻尼分別為100000 N/m 和5000 Ns/m。為了研究保護(hù)軸承設(shè)計(jì)參數(shù)和運(yùn)行條件對(duì)保護(hù)軸承疲勞壽命的影響，考慮了摩擦系數(shù)、永磁側(cè)載荷、氣隙和轉(zhuǎn)子速度。圖4 顯示了轉(zhuǎn)子跌落時(shí)的軌道圖、剪切應(yīng)力時(shí)間曲線、雨流量直方圖和溫度曲線圖。

圖3 仿真模型

圖4 仿真結(jié)果

4 結(jié)論

在本研究中，使用考慮熱增長的非線性滾珠軸承模型和有限元建模的柔性轉(zhuǎn)子，研究了保護(hù)軸承和轉(zhuǎn)子的動(dòng)力學(xué)。此外，還計(jì)算了軸承部件中的接觸載荷、赫茲應(yīng)力、亞剪切應(yīng)力、表面剪切應(yīng)力和熱增長。通過將雨落計(jì)數(shù)方法應(yīng)用于剪切應(yīng)力時(shí)程，將保護(hù)軸承的疲勞壽命預(yù)測(cè)為轉(zhuǎn)子跌落次數(shù)。