馬俊，司東宏，薛玉君，，董永祥

（1.河南科技大學(xué) 機(jī)電工程學(xué)院，河南 洛陽(yáng) 471003；2.河南省機(jī)械設(shè)計(jì)及傳動(dòng)系統(tǒng)重點(diǎn)實(shí)驗(yàn)室，河南 洛陽(yáng) 471003）

主軸系統(tǒng)的動(dòng)態(tài)特性是直接影響高速雙轉(zhuǎn)子試驗(yàn)機(jī)整體性能的主要因素。主軸軸系受力復(fù)雜，不僅承受徑向加載力，還承受轉(zhuǎn)子旋轉(zhuǎn)質(zhì)量慣性力引起的周期性振動(dòng)，而周期性振動(dòng)會(huì)使軸系的疲勞壽命降低；因此，分析試驗(yàn)機(jī)主軸系統(tǒng)的動(dòng)態(tài)特性對(duì)于提高試驗(yàn)機(jī)整體穩(wěn)定性、可靠性和延長(zhǎng)使用壽命具有重要的意義［1］。

文獻(xiàn)［2］分別分析了雙轉(zhuǎn)子軸承試驗(yàn)機(jī)左右端主軸的固有頻率，計(jì)算出單根軸的臨界轉(zhuǎn)速，由此得出單根軸的固有頻率滿足使用要求。然而，雙轉(zhuǎn)子軸承試驗(yàn)機(jī)單根軸的模態(tài)變化不能等同于雙軸耦聯(lián)時(shí)的模態(tài)變化，因此，使用有限元分析軟件ANSYSWorkbench分析雙軸耦聯(lián)時(shí)的軸系模態(tài)變化，同時(shí)考慮高轉(zhuǎn)速軸系旋轉(zhuǎn)時(shí)產(chǎn)生的預(yù)應(yīng)力對(duì)軸系模態(tài)的影響。

1 模態(tài)分析理論

模態(tài)分析是研究機(jī)械結(jié)構(gòu)動(dòng)力學(xué)特性、振動(dòng)分析和動(dòng)態(tài)優(yōu)化設(shè)計(jì)的常用方法。多自由度結(jié)構(gòu)系統(tǒng)的動(dòng)力學(xué)方程可以表示為［3-4］

式中：M，C，K分別為結(jié)構(gòu)系統(tǒng)的質(zhì)量、阻尼、剛度矩陣；F（t）為外激勵(lì)矢量；X為位移矢量。

在實(shí)際工況中，阻尼對(duì)結(jié)構(gòu)的固有頻率和振型影響不大，則無(wú)外載荷作用且忽略阻尼時(shí)，可將（1）式簡(jiǎn)化為

當(dāng)軸系旋轉(zhuǎn)時(shí)，由于離心力的作用會(huì)產(chǎn)生一定的旋轉(zhuǎn)預(yù)應(yīng)力，從而影響到軸系的模態(tài)振動(dòng)性能。因此，在考慮預(yù)應(yīng)力的情況下，（2）式方程解的形式可表示為

式中：kr為預(yù)應(yīng)力剛度矩陣；m為振幅列陣；ω為角頻率；φ為初相角。通過(guò)求解上式可以得到ω1，ω2，…，ωn等n個(gè)固有頻率，并且滿足0≤ω1＜ω2＜…＜ωn。

2 主軸系統(tǒng)有限元分析

2.1 軸系模型的建立

使用Unigraphics軟件建立試驗(yàn)機(jī)主軸系統(tǒng)結(jié)構(gòu)模型，主軸系統(tǒng)主要有內(nèi)軸、外軸、陪試軸承（角接觸球軸承、圓柱滾子軸承）、試驗(yàn)軸承（角接觸球軸承）、徑向加載軸承等組成（圖1）。試驗(yàn)軸承將內(nèi)軸和外軸耦聯(lián)在一起。陪試角接觸球軸承分別裝在內(nèi)軸左端和外軸右端，陪試圓柱滾子軸承裝在外軸左端，主要對(duì)軸系起支承作用；徑向加載軸承對(duì)軸系不起支承作用，計(jì)算時(shí)可以忽略其剛度。為了提高網(wǎng)格劃分質(zhì)量并降低有限元分析時(shí)間，建模的過(guò)程中忽略一些不影響分析結(jié)果的倒角、倒圓角、孔等特征［5-7］。

圖1 軸系模型Fig.1 Shafting model

2.2 材料特性

主軸是軸承試驗(yàn)機(jī)的關(guān)鍵部件，必須具備足夠的強(qiáng)度、剛度及耐磨性，且需要熱處理變形小，因此需選擇彈性模量、泊松比、屈服強(qiáng)度合適的材料。根據(jù)以上特點(diǎn)，主軸材料選用42CrMo，線性各向同性，詳細(xì)參數(shù)見表1。

表1 主軸材料屬性Tab.1 Material properties of spindle

2.3 有限元模型的建立

通過(guò)Unigraphics將軸系模型導(dǎo)成ANSYS Workbench分析軟件可以識(shí)別的X＿T格式。根據(jù)本次分析所用計(jì)算機(jī)的配置及對(duì)分析結(jié)果的精確度要求，選用了Sizing劃分網(wǎng)格，網(wǎng)格單元類型為CTETRA（10），單元大小為15；部件中的單元總數(shù)為34179，節(jié)點(diǎn)總數(shù)為57438，有限元模型如圖2所示。

圖2 有限元模型Fig.2 Finite element model

2.4 軸承剛度的計(jì)算

根據(jù)試驗(yàn)機(jī)主軸系統(tǒng)裝配方式可知，陪試軸承和試驗(yàn)軸承對(duì)主軸系統(tǒng)起支承作用，需要計(jì)算其支承剛度，一般將軸承簡(jiǎn)化成彈簧支承。實(shí)際工況中，軸承只有支承剛度而沒(méi)有拉剛度，因此，圓柱滾子軸承只有一個(gè)方向的支承剛度，而角接觸球軸承接觸力比較復(fù)雜，在2個(gè)互相垂直的方向各自建立一個(gè)彈簧支承，軸承支承簡(jiǎn)化效果如圖3所示。

圖3 軸系支承示意圖Fig.3 Diagram of shafting support

現(xiàn)計(jì)算出滿足設(shè)計(jì)要求的軸承剛度的最小值，支承軸承參數(shù)見表2，假設(shè)不考慮軸系上不同型號(hào)軸承的交叉剛度和阻尼的影響，只考慮軸承對(duì)軸系支承正剛度的作用［8-9］。

表2 軸承參數(shù)Tab.2 Parameters for bearing

根據(jù)經(jīng)驗(yàn)公式［10］可知，角接觸球軸承及圓柱滾子軸承的徑向剛度分別為

將表2中的參數(shù)代入（4），（5）式可得K＝176957.4369 N/mm，K*＝1120621.355 N/mm。

2.5 邊界條件

在ANSYS Workbench中，約束角接觸球軸承T2和T4的全部自由度，T1和T3的軸向移動(dòng)自由度；約束圓柱滾子軸承T6的全部自由度，如圖4所示。

圖4 邊界條件Fig.4 Boundary condition

2.6 模態(tài)分析

雙轉(zhuǎn)子軸系設(shè)計(jì)通常只做靜態(tài)下的模態(tài)分析，不考慮軸旋轉(zhuǎn)時(shí)產(chǎn)生的預(yù)應(yīng)力對(duì)軸系模態(tài)的影響，但在實(shí)際工況下，軸系旋轉(zhuǎn)時(shí)產(chǎn)生的預(yù)應(yīng)力不僅會(huì)影響軸系的剛度和強(qiáng)度，而且也影響軸承剛度；隨著軸系轉(zhuǎn)速的升高預(yù)應(yīng)力會(huì)變大，從而改變與內(nèi)外圈的接觸角，使軸承的剛度軟化，并對(duì)軸系的模態(tài)產(chǎn)生相應(yīng)影響［11］。

考慮到軸系旋轉(zhuǎn)形成的預(yù)應(yīng)力對(duì)軸承剛度的影響，分析軸系在靜態(tài)和不同轉(zhuǎn)速下的預(yù)應(yīng)力對(duì)軸系模態(tài)的影響，其趨勢(shì)如圖5所示。由圖可知：靜態(tài)下軸系模態(tài)的各階固有頻率值大于考慮預(yù)應(yīng)力的各階固有頻率值，這是由于軸系旋轉(zhuǎn)形成的預(yù)應(yīng)力弱化了軸承的支承剛度，從而降低了固有頻率；而隨著轉(zhuǎn)速的升高，軸系的固有頻率下降越來(lái)越明顯，當(dāng)轉(zhuǎn)速達(dá)到15000 r/min時(shí)，軸系1階固有頻率與靜態(tài)下的1階固有頻率相比降低了25.9%，這是因?yàn)殡S著轉(zhuǎn)速的增加，軸系產(chǎn)生的預(yù)應(yīng)力增大，使球與內(nèi)圈的接觸角增大，與外圈的接觸角減小，并改變了接觸力的大小，使軸承的剛度隨著預(yù)應(yīng)力的增加而降低，最終軸系固有頻率因軸承剛度的降低而減小。

圖5 不同轉(zhuǎn)速下軸承各階固有頻率的變化趨勢(shì)Tab.5 Variation trend of each natural frequency of bearing under different speeds

從圖中還可以看出：1階固有頻率的變化最為明顯，這說(shuō)明預(yù)應(yīng)力對(duì)1階固有頻率影響較大，為改善預(yù)應(yīng)力對(duì)主軸模態(tài)頻率的影響，軸系工作時(shí)，應(yīng)增加與其對(duì)應(yīng)的軸向力，提高軸承剛度。

3 軸系臨界轉(zhuǎn)速的計(jì)算

臨界轉(zhuǎn)速是軸系設(shè)計(jì)時(shí)工作轉(zhuǎn)速避開共振頻率點(diǎn)的主要參考依據(jù)，軸系在臨界轉(zhuǎn)速附近工作時(shí)，將產(chǎn)生劇烈的振動(dòng)，長(zhǎng)時(shí)間運(yùn)行會(huì)造成主軸部件的嚴(yán)重彎曲變形，甚至?xí)斐赊D(zhuǎn)子和軸承的破壞［12］。通過(guò)分析計(jì)算軸系的模態(tài)，求得各階固有頻率，從而確定軸系各階臨界轉(zhuǎn)速。臨界轉(zhuǎn)速與固有頻率的關(guān)系為［13］

式中：n為臨界轉(zhuǎn)速，r/min；f為固有頻率，Hz。

3.1 Campbell 圖求解

在ANSYS Workbench中，通過(guò)不同轉(zhuǎn)速下的多載荷步模態(tài)分析，可以獲得Campbell圖并能得到軸系的臨界轉(zhuǎn)速，軸系的Campbell圖如圖6所示。實(shí)際工況中軸系的工作轉(zhuǎn)速小于1階臨界轉(zhuǎn)速，因此只關(guān)注1階臨界轉(zhuǎn)速值，約為23206.2 r/min。

圖6 Campbell圖Tab.6 Campbell diagram

3.2 理論計(jì)算值

軸系由內(nèi)軸和外軸通過(guò)試驗(yàn)圓柱滾子軸承連接在一起，因軸的剛度較大，故忽略軸的彈性而將其設(shè)為質(zhì)量塊，軸系上的軸承簡(jiǎn)化為彈簧，所以整個(gè)軸系的簡(jiǎn)化模型如圖7所示。

圖7 軸系簡(jiǎn)化模型Tab.7 Simplified model for shafting

軸系簡(jiǎn)化模型中：X1，X2分別為內(nèi)、外軸的位移；M1和M2分別為內(nèi)、外軸的質(zhì)量，M1＝50.53 kg，M2＝42.72 kg；K1～K4為軸承的剛度，由（4），（5）式求得K1＝K2＝K4＝K＝176957.44 N/mm，K3＝K*＝1120621.36 N/mm。

由軸系的簡(jiǎn)化模型圖建立能量方程式，系統(tǒng)在任意時(shí)刻的動(dòng)能T和勢(shì)能V分別為

由（7）式減（8）式可得Lagrange函數(shù)為

對(duì)比理論計(jì)算與有限元分析的結(jié)果可知：理論分析值略大于有限元分析值，這是由于有限元軟件自動(dòng)考慮了主軸材料的阻尼。根據(jù)軸設(shè)計(jì)手冊(cè)可知，為了軸系的安全運(yùn)行，通常取軸系1階臨界轉(zhuǎn)速值的70%作為軸系在剛性階段的極限工作轉(zhuǎn)速。使用2種方法得到的1階臨界轉(zhuǎn)速分別為23206.2 和 24660 r/min，其 70% 分 別 為16244.34和17262 r/min，均高于主軸系統(tǒng)最高設(shè)計(jì)工作轉(zhuǎn)速15000 r/min，說(shuō)明試驗(yàn)機(jī)軸系的設(shè)計(jì)符合要求。

4 結(jié)束語(yǔ)

通過(guò)ANSYSWorkbench中Block Lanczos有限元法計(jì)算了軸系靜態(tài)下和預(yù)應(yīng)力下的2種模態(tài)頻率值，對(duì)比可知：軸系轉(zhuǎn)速在3000～15000 r/min時(shí)，預(yù)應(yīng)力下的模態(tài)值都低于靜態(tài)下的模態(tài)值，且隨著轉(zhuǎn)速的升高差值越來(lái)越大；當(dāng)轉(zhuǎn)速達(dá)15000 r/min時(shí)，1階固有頻率相差了25.9%，表明預(yù)應(yīng)力對(duì)軸系的影響明顯，因此模態(tài)分析時(shí)應(yīng)考慮預(yù)應(yīng)力的影響。

根據(jù)Campbell圖和理論分析法獲得軸系的1階臨界轉(zhuǎn)速分別為23206.2和24660 r/min，均遠(yuǎn)大于主軸的設(shè)計(jì)工作轉(zhuǎn)速，可以滿足設(shè)計(jì)要求；而有限元分析的數(shù)值小于理論值，說(shuō)明有限元分析的數(shù)值比較保守，在實(shí)際工況中，為了安全考慮采用此臨界轉(zhuǎn)速值。