基于Samcef Rotor和Cobra的某滾動軸承-轉(zhuǎn)子系統(tǒng)固有頻率分析

洛陽軸承研究所有限公司 于曉凱 徐俊 謝鵬飛

河南科技大學機電工程學院 屈馳飛

基于Samcef Rotor和Cobra的某滾動軸承-轉(zhuǎn)子系統(tǒng)固有頻率分析

洛陽軸承研究所有限公司 于曉凱 徐俊 謝鵬飛

河南科技大學機電工程學院 屈馳飛

建立轉(zhuǎn)子系統(tǒng)Samcef Rotor有限元模型，通過軸承分析軟件Cobra分析得到軸承5×5剛度矩陣，將軸承剛度矩陣與轉(zhuǎn)子動力學軟件Samcef Rotor耦合分析，可精確得到軸承-轉(zhuǎn)子系統(tǒng)固有頻率，通過試驗驗證了分析方法的有效性和正確性。

一、前言

隨著電力、航空航天、石油化工及機械制造業(yè)的飛速發(fā)展，各種旋轉(zhuǎn)機械向高速方向發(fā)展，因此對轉(zhuǎn)子-軸承系統(tǒng)的動力學特性提出了更高的要求。隨之對軸承的設計和使用提出了越來越高的要求，因此，研究滾動軸承-轉(zhuǎn)子系統(tǒng)的動態(tài)特性具有重要的實用價值和意義。

滾動軸承支承的轉(zhuǎn)子系統(tǒng)的動力學特性取決于支承轉(zhuǎn)子的滾動軸承的動力學特性，尤其是軸承的徑向剛度和軸向剛度。在高速旋轉(zhuǎn)機械中預測轉(zhuǎn)子系統(tǒng)的固有頻率從而避開轉(zhuǎn)子系統(tǒng)固有頻率具有重要作用，轉(zhuǎn)子系統(tǒng)的固有頻率主要受軸承剛度影響，軸承剛度主要受球與內(nèi)、外圈的接觸狀態(tài)影響。分析滾動軸承支承的轉(zhuǎn)子系統(tǒng)時，將軸承假設為鉸鏈支承或固定支承。這樣的假定在軸承轉(zhuǎn)速較低時是允許的，隨著轉(zhuǎn)子系統(tǒng)轉(zhuǎn)速的提高或轉(zhuǎn)子柔性化，轉(zhuǎn)子系統(tǒng)往往需要跨越一階或二階臨界轉(zhuǎn)速，這種考慮是不完善的或?qū)⒁鸷艽笳`差，只有考慮軸承的動力學特性才能得到更加精確的分析結(jié)果。

本文考慮軸承轉(zhuǎn)速、軸承載荷和潤滑狀態(tài)等建立了轉(zhuǎn)子的5×5自由度振動模型，利用軸承計算軟件Cobra、三維建模軟件Solidworks和轉(zhuǎn)子動力學軟件Rotor分析軸承轉(zhuǎn)速和軸承載荷對轉(zhuǎn)子系統(tǒng)固有頻率的影響，通過試驗驗證分析的正確性。

二、建立有限元模型

對于實際問題，零件結(jié)構(gòu)往往比較復雜，通過理論分析得到轉(zhuǎn)子系統(tǒng)的剛度和軸承剛度比較困難，本文通過轉(zhuǎn)子動力學軟件Samcef Rotor和軸承計算軟件Cobra軟件求出軸承-轉(zhuǎn)子系統(tǒng)的臨界轉(zhuǎn)速分析并經(jīng)試驗驗證。

1.基本假設

在建立有限元模型前，先對轉(zhuǎn)子系統(tǒng)進行以下基本假設：（1）轉(zhuǎn)子的振動位移量較小，軸承處的變形很小即軸承的徑向游隙等基本參數(shù)改變很小，可以忽略不計；（2）各零件之間的連接為剛性連接，如端蓋與軸承端面之間的連接；（3）系統(tǒng)的阻尼可忽略不計；（4）軸承滾動體和保持架對轉(zhuǎn)子系統(tǒng)的固有頻率影響很小，可忽略不計。

2.實體模型建立

在轉(zhuǎn)子動力學Samcef Rotor分析類型“Rotor Dynamics”（ 轉(zhuǎn) 子 動 力 學） 和“Critical Speed &Stability”（臨界轉(zhuǎn)速和穩(wěn)定性分析）模塊下建立轉(zhuǎn)子系統(tǒng)實體模型。Cobra高級軸承計算軟件坐標系不能調(diào)整，為了保證Samcef Rotor和Cobra坐標系的統(tǒng)一，在Samcef Rotor建立系統(tǒng)實體模型時，應將轉(zhuǎn)子的軸線與Z軸重合。另外，實體模型中滾動軸承僅建立內(nèi)圈和外圈實體模型，不考慮滾動體和保持架部分。

建立Samcef Rotor實體模型后（圖1），在“Material”模塊中賦予轉(zhuǎn)子系統(tǒng)的材料特性。

圖1 滾動軸承-轉(zhuǎn)子實體模型

3.滾動軸承剛度矩陣計算

Cobra是高級高速滾動軸承和滾子軸承系統(tǒng)優(yōu)化分析軟件，是由美國航空航天局資助的軸承系統(tǒng)專用分析軟件，在NASA多個軸承系統(tǒng)中都得到了很好的應用，該軟件已廣泛應用在國外各大軸承公司的軸承系統(tǒng)設計分析方面。

在Cobra軟件中，輸入軸承載荷（主要包括軸向載荷和徑向載荷）、軸承結(jié)構(gòu)參數(shù)等，具體如表1所示，通過計算得到軸承的剛度矩陣，具體如表2所示。

表1 軸承結(jié)構(gòu)參數(shù)及軸承工況

表2 Cobra分析得B7000軸承5×5剛度矩陣

4.施加邊界條件及劃分網(wǎng)格

滾動體和保持架質(zhì)量很小，其影響系統(tǒng)固有頻率主要通過軸承剛度來影響，因此忽略滾動體和保持架對轉(zhuǎn)子系統(tǒng)固有頻率的影響，但軸承剛度不能忽略將軸承剛度導入到Samcef Rotor軟件中，即通過Samcef Rotor中“Bearing”約束來代替軸承的幾何模型。

將軸承內(nèi)、外圈約束類型選擇為“Bearing”，將滾動軸承剛度導入”Bearing”中“Stiffness Matrix”，利用Samcef Rotor對模型劃分網(wǎng)格，有限元模型如圖2所示。

圖2 某滾動軸承—轉(zhuǎn)子系統(tǒng)有限元模型

三、計算結(jié)果

根據(jù)模態(tài)理論，多自由系統(tǒng)自由振動時，其固有頻率個數(shù)等于自由度數(shù)。某滾動軸承-轉(zhuǎn)子系統(tǒng)被離散為數(shù)十萬個單元，由于高階模態(tài)對振動系統(tǒng)影響不大。一般情況下，分析滾動軸承轉(zhuǎn)子系統(tǒng)0～2000Hz的固有頻率可滿足工程應用的精度要求，計算結(jié)果如表3所示，各階模態(tài)圖如圖3所示。

表3 某滾動軸承轉(zhuǎn)子系統(tǒng)0～2000Hz固有頻率

四、試驗結(jié)果及對比分析

為了驗證分析結(jié)果的正確性，對轉(zhuǎn)子系統(tǒng)進行隨機振動，振動量級為0.25g/Hz，因為輪緣外徑最大，對振動比較敏感在輪緣部分和振動臺基座上各放置一個傳感器，振動現(xiàn)場圖如圖3所示，系統(tǒng)振動結(jié)果如圖4所示。通過分析可得系統(tǒng)固有頻率值主要有148Hz、413 Hz 、876Hz、974Hz、1298Hz和1650Hz。當考慮滾動軸承時，通過Samcef Rotor分析得到系統(tǒng)固有頻率，各階固有頻率均與實測值比較接近，最大誤差僅8.1%，具體如表4所示。

圖3 某滾動軸承轉(zhuǎn)子系統(tǒng)0～2000Hz固有頻率模態(tài)圖

圖4 某滾動軸承轉(zhuǎn)子系統(tǒng)振動現(xiàn)場圖

表4 考慮滾動軸承，利用Samcef Rotor分析系統(tǒng)固有頻率與實測值對比結(jié)果

根據(jù)多自由度機械振動理論可知，當系統(tǒng)的兩固有頻率比較接近時，系統(tǒng)實際振動僅顯示出一個疊加峰值，例如Samcef Rotor分析909Hz與913Hz存在該種情況。因此，分析時將909Hz和913Hz同時與實測876Hz進行對比。

另外，利用傳統(tǒng)方法，將滾動軸承設置為剛性，分析得到系統(tǒng)固有頻率。與考慮滾動軸承方法相比，各階固有頻率均與實測值相差較大，最大誤差僅75.6%，具體結(jié)果如表5所示。

圖5 隨機振動結(jié)果

表5 將滾動軸承設置為剛體時，利用Samcef Rotor分析系統(tǒng)固有頻率與實測值對比結(jié)果

五、結(jié)語

經(jīng)過分析，傳統(tǒng)滾動軸承支承轉(zhuǎn)子系統(tǒng)將將滾動軸承簡化為剛性時，分析得得系統(tǒng)固有頻率誤差較大?？紤]滾動軸承柔性時，利用軸承分析軟件Cobra分析與轉(zhuǎn)子動力學軟件Samcef Rotor耦合分析，能精確得到軸承-轉(zhuǎn)子系統(tǒng)固有頻率。通過試驗對比，驗證了本文分析方法的有效性和正確性。