考慮波紋度和不平衡力的薄壁軸承-轉(zhuǎn)子系統(tǒng)非線性分析

康鋒，張耀強，楊茹萍，牛青波

(1. 河南科技大學(xué) 土木工程學(xué)院，河南 洛陽 471003；2. 洛陽軸研科技股份有限公司，河南 洛陽 471039)

薄壁軸承因輕量化等特點,在機器人手臂關(guān)節(jié)、航空發(fā)動機等結(jié)構(gòu)中廣泛應(yīng)用。在軸承加工制造過程中，由于加工工藝等原因，滾道表面不可避免會產(chǎn)生波紋度，在長期磨損后也會出現(xiàn)質(zhì)量偏心，從而產(chǎn)生不平衡力，引起強迫振動。在軸承-轉(zhuǎn)子系統(tǒng)中，由于這些非線性激勵源的存在，會引起諸多用線性理論難以解釋的動力學(xué)現(xiàn)象和一些不可預(yù)測的破壞后果[1]。近年來，眾多學(xué)者開始用非線性動力學(xué)理論來分析這些現(xiàn)象。文獻[2]考慮了波紋度等非線性因素對滾子軸承-轉(zhuǎn)子系統(tǒng)的影響，其研究表明波紋度引起的振動頻率會使系統(tǒng)產(chǎn)生強烈的振動。文獻[3]研究了帶有波紋度和不平衡轉(zhuǎn)子的高速球軸承的非線性動力學(xué)特性，由于不平衡轉(zhuǎn)子和滾動體波紋度的作用，振動頻譜圖的最大振幅出現(xiàn)在波通過頻率、旋轉(zhuǎn)頻率及組合頻率處。文獻[4]提出了一種在考慮滾動零件波紋度、離心力和陀螺力矩的影響下計算球軸承特性的方法，分析表明由于波紋度的影響，軸承振動頻率成分發(fā)生了變化。文獻[5]在考慮非線性因素和不平衡力的基礎(chǔ)上建立動力學(xué)方程，對滾動軸承-轉(zhuǎn)子系統(tǒng)進行了分析，結(jié)果表明增大不平衡力，強迫振動會增強，并誘發(fā)系統(tǒng)產(chǎn)生混沌振動。文獻[6-7]所做的研究僅考慮了強迫振動，沒有涉及軸承內(nèi)外圈波紋度等相關(guān)參數(shù)。以上研究均未充分考慮內(nèi)外圈波紋度和強迫振動等因素，且關(guān)于薄壁軸承的研究較少。

現(xiàn)以薄壁軸承-轉(zhuǎn)子系統(tǒng)為研究對象，在綜合考慮內(nèi)外圈波紋度、彈性恢復(fù)力、變?nèi)岫鹊确蔷€性因素和不平衡力的基礎(chǔ)上，建立軸承-轉(zhuǎn)子系統(tǒng)非線性動力學(xué)方程，采用RK4算法求解，結(jié)合分叉圖、Poincaré映射圖和頻譜圖對薄壁軸承-轉(zhuǎn)子系統(tǒng)的混沌、分岔及振動響應(yīng)頻率等非線性動力學(xué)特性進行分析。

1 薄壁軸承-轉(zhuǎn)子系統(tǒng)的動力學(xué)分析模型

薄壁軸承-轉(zhuǎn)子系統(tǒng)可以簡化為一個剛性平衡轉(zhuǎn)子兩端由2個相同的薄壁軸承的對稱支承。薄壁軸承-轉(zhuǎn)子系統(tǒng)動力學(xué)模型如圖1所示，假設(shè)只承受徑向載荷[8]，其內(nèi)外圈分別與轉(zhuǎn)軸、機架剛性接觸，滾動體均勻分布在保持架當中做純滾動[9]，內(nèi)外滾道表面具有正弦波的波紋度[10]，設(shè)表面波紋度“凸出”為正，“凹進”為負[11]，不計潤滑作用，O,Oe分別為轉(zhuǎn)子的形心和質(zhì)心，e為轉(zhuǎn)子質(zhì)心的偏移量。薄壁軸承-轉(zhuǎn)子系統(tǒng)內(nèi)部關(guān)系示意圖如圖2所示，A為波紋度最大幅值，p為波紋度幅值，F(xiàn)e為質(zhì)量偏心引起的不平衡力，ω為轉(zhuǎn)速，φj為第j個滾動體在t時刻的角位置。

圖1 薄壁軸承-轉(zhuǎn)子系統(tǒng)動力學(xué)模型

圖2 薄壁軸承-轉(zhuǎn)子系統(tǒng)內(nèi)部關(guān)系示意圖

1.1 薄壁軸承內(nèi)外圈波紋度

在t時刻，與第j個滾動體接觸的內(nèi)外滾道表面波紋度幅值pij和pej為

j=1，2，3…Z，

(1)

式中：Nw為波紋度波數(shù)；β為波紋度初始相位角，(°)；Z為滾動體個數(shù)；ω0,ω1分別為保持架、內(nèi)圈的角速度，r/min；D為滾道直徑，mm；下標i，e分別表示軸承內(nèi)、外圈。

1.2 彈性恢復(fù)力

根據(jù)Hertz接觸理論可知，第j個滾動體與內(nèi)外滾道接觸的局部彈性恢復(fù)力Fj與彈性接觸變形量uj的關(guān)系為

(2)

uj=[(xcosφj+ysinφj)cosα-Gr+pij-pej]+，

式中：Kb為系統(tǒng)總接觸剛度，N/m1.5；τ=3/2；x,y分別為內(nèi)圈中心在豎直和水平方向上的位移，mm；α為接觸角；Gr為徑向游隙，μm；下標“+”表示如果括號內(nèi)的值為負或0，取uj=0。

在薄壁軸承-轉(zhuǎn)子系統(tǒng)中，彈性恢復(fù)力F為所有局部彈性恢復(fù)力Fj之和，即

(3)

1.3 薄壁軸承-轉(zhuǎn)子系統(tǒng)的非線性振動微分方程

根據(jù)Lagrange方程，薄壁軸承-轉(zhuǎn)子系統(tǒng)的非線性振動微分方程為

(4)

式中：M為轉(zhuǎn)子系統(tǒng)內(nèi)圈、轉(zhuǎn)子和滾動體的總質(zhì)量，kg；C為軸承運轉(zhuǎn)時的等效阻尼系數(shù)，N·s/m；Fr為施加在轉(zhuǎn)子上的恒定垂直力，N。

軸承在運轉(zhuǎn)過程中，其轉(zhuǎn)子的振動為變?nèi)岫日駝樱晦D(zhuǎn)子質(zhì)量偏心也會引起系統(tǒng)發(fā)生強迫振動；滾動體每次通過內(nèi)圈波紋度波峰時，系統(tǒng)會產(chǎn)生波通過振動WPV(Wave Passage Vibration)[12]。這些振動頻率對系統(tǒng)的非線性動力學(xué)特性都會有影響，但影響程度會隨著工況條件的改變而改變。變?nèi)岫日駝宇l率fvc、強迫振動頻率ffv、波通過振動頻率fwpv分別為

(5)

變?nèi)岫日駝宇l率fvc與強迫振動頻率ffv的關(guān)系為

fvc=ffvBN，

(6)

式中：BN為與薄壁軸承相關(guān)的系數(shù)，其值取決于軸承尺寸。

2 系統(tǒng)非線性動力學(xué)特性分析

以某機器人薄壁軸承單元為例，利用MATLAB對薄壁軸承-轉(zhuǎn)子系統(tǒng)進行動態(tài)數(shù)值仿真分析。薄壁軸承-轉(zhuǎn)子系統(tǒng)的主要參數(shù)：M=6.0 kg，F(xiàn)r=10 N，C=220 N·s/m[11]，Z=42，Di=213.808 mm，De=236.284 mm，α=0°，Gr=10 μm。

波紋度和不平衡力是影響系統(tǒng)非線性特性的重要因素。由于滾動軸承-轉(zhuǎn)子系統(tǒng)振動微分方程的強非線性特性，采用RK4算法進行求解，積分步長取一個激勵周期的1/300，計算結(jié)果有選擇性的采用分叉圖、Poincaré映射圖或頻譜圖分析系統(tǒng)的非線性動力學(xué)特性。

2.1 僅考慮外圈波紋度和不平衡力的非線性分析

為便于分析，參考文獻[13-14]，取Ai=0 μm，Ae=1 μm，不平衡力Fe=0.05Fr，x方向位移隨轉(zhuǎn)速變化的分叉圖和不同轉(zhuǎn)速下x方向的Poincaré映射圖分別如圖3和圖4所示。

圖3 x方向位移隨轉(zhuǎn)速變化的分叉圖(Fe=0.05Fr)

由圖可以看出：在轉(zhuǎn)速nrotor=820～960 r/min時，系統(tǒng)呈4周期振動;在nrotor=970～1 420 r/min時，系統(tǒng)呈2周期振動;在nrotor=630～810 r/min和1 430～2 010 r/min時，系統(tǒng)處于混沌狀態(tài);在nrotor=3 070～4 450 r/min時，系統(tǒng)呈擬周期振動。

其他參數(shù)不變，取不平衡力Fe=0.25Fr。x方向位移隨轉(zhuǎn)速變化的分叉圖如圖5所示，由圖可知，系統(tǒng)響應(yīng)混沌區(qū)由2個增加到3個，分別為轉(zhuǎn)速nrotor=520～860 r/min，nrotor=1 430～ 2 520 r/min和nrotor=3 150～3 990 r/min。由此可見，增大不平衡力會使系統(tǒng)混沌區(qū)增大。

圖5 x方向位移隨轉(zhuǎn)速變化的分叉圖(Fe=0.25Fr)

2.2 僅考慮內(nèi)圈波紋度和不平衡力的非線性分析

取內(nèi)圈波紋度最大幅值A(chǔ)i=1 μm，Ae=0 μm，不平衡力Fe=0.05Fr。nrotor=740,2 050,3 680 r/min時的振動響應(yīng)的振幅頻譜圖分別如圖6所示。

圖6 振幅頻譜圖(Ai=1 μm,Ae=0 μm)

由圖可知,系統(tǒng)在nrotor=740 r/min時，豎直x方向主要是波通過振動頻率fwpv分量，水平y(tǒng)方向主要是強迫振動頻率ffv倍頻分量，y方向比x方向振動復(fù)雜。當nrotor=2 050 r/min時，x方向fwpv分量減少，ffv出現(xiàn)，y方向僅有1倍、2倍ffv分量。在nrotor= 3 680 r/min時，x方向ffv分量明顯增大，fwpv減小，y方向只有1倍ffv。由此可見，隨轉(zhuǎn)速增大，強迫振動頻率在系統(tǒng)振動響應(yīng)中逐漸增強。

其他參數(shù)不變，取不平衡力Fe=0.25Fr。以系統(tǒng)非線性動力學(xué)特性最顯著的轉(zhuǎn)速為例分析，nrotor=3 680 r/min時的振動響應(yīng)的振幅頻譜圖如圖7所示。

圖7 轉(zhuǎn)速nrotor=3 680 r/min的振幅頻譜圖

由圖6和圖7可知，增大不平衡力，x方向振動頻譜圖變?yōu)橐贿B續(xù)的寬頻譜，頻率成分復(fù)雜，y方向出現(xiàn)新的峰值頻率，即波通過頻率分頻分量。同時x,y方向最大振幅也增大。x,y方向最大振幅隨不平衡力的變化趨勢如圖8所示，由圖可知，系統(tǒng)y方向最大振幅大于x方向，并隨不平衡力的增大呈波浪式遞增。由此可見，增大不平衡力，系統(tǒng)的振動響應(yīng)頻率變得復(fù)雜，最大振幅呈波浪式遞增。

圖8 x,y方向最大振幅隨不平衡力的變化

2.3 考慮內(nèi)外圈波紋度和不平衡力的非線性分析

取軸承內(nèi)外圈波紋度最大幅值A(chǔ)i=Ae=1 μm，不平衡力為Fe=0.05Fr和Fe=0.25Fr時，轉(zhuǎn)速nrotor=3 680 r/min時的振動響應(yīng)的振幅頻譜圖如圖9所示。

圖9 振幅頻譜圖(Ai=Ae=1 μm)

由圖可知，x方向的振動響應(yīng)比y方向復(fù)雜，其響應(yīng)有強迫振動頻率、波通過振動頻率和變?nèi)岫日駝宇l率。不平衡力較大時，系統(tǒng)振動會出現(xiàn)變?nèi)岫日駝宇l率與強迫振動頻率的組合頻率及波通過振動頻率與強迫振動頻率的組合頻率。不平衡力增大會使系統(tǒng)的振動頻率組成更加復(fù)雜，同時振幅也會增大。由此可見，不平衡力增大會使系統(tǒng)的振動響應(yīng)頻率變得復(fù)雜，其對水平y(tǒng)方向的振動影響要大于豎直x方向。

3 結(jié)束語

通過文中建立的模型來分析不平衡力和轉(zhuǎn)速對薄壁軸承-轉(zhuǎn)子系統(tǒng)非線性動力學(xué)特性的影響。結(jié)果表明：增大不平衡力會使系統(tǒng)的混沌振動范圍增大；隨轉(zhuǎn)速增加，強迫振動頻率在系統(tǒng)振動響應(yīng)中逐漸占據(jù)主要地位；不平衡力對系統(tǒng)水平方向振動的影響要遠大于豎直方向。