動壓潤滑—徑向滑動軸承的動力特性計算
2018年07月20日油膜的動力特性,包括油膜剛度和阻尼,反應(yīng)了當(dāng)軸頸偏離了靜平衡位置并在此位置附近作變位運(yùn)動時油膜力相應(yīng)的變化情況。
油膜不僅承受載荷,減輕摩擦,消除磨損等作用,從動力學(xué)角度來說,油膜是轉(zhuǎn)子軸承系統(tǒng)中的一個環(huán)節(jié)。油膜的特性,對于整個轉(zhuǎn)子系統(tǒng)的動力性態(tài)有很大影響。油膜會影響轉(zhuǎn)子系統(tǒng)的穩(wěn)定性、轉(zhuǎn)子系統(tǒng)的臨界轉(zhuǎn)速、轉(zhuǎn)子不平衡所引起的振幅、過臨界轉(zhuǎn)速時的共振放大倍數(shù)等等。
計算油膜動力特性必須以不定常工況的Reynolds方程作為分析、計算的基礎(chǔ)(雷諾方程的求解可參考:徑向滑動軸承靜特性計算(一))),其無量綱形式:
求解剛度、阻尼常采用兩種方法:
(1) 用差分概念計算剛度、阻尼
(2) 用偏導(dǎo)數(shù)概念求剛度、阻尼
其中,用差分概念計算剛度、阻尼:
① 在平衡位置取位移微小擾動,由幾何關(guān)系計算新位置下的偏心率和偏位角;
② 在新位置下直接求解雷諾方程,Simpson積分計算垂直、水平分力;
③ 由求出的力的變化量和位移擾動量求出剛度系數(shù);
④ 同理,在平衡位置去微小速度擾動,求解雷諾方程,積分得垂直、水平分力;
⑤ 由求出的力的變化量和速度擾動量求出阻尼系數(shù)。
以Kxx為例,在x方向給以+△x、-△x位移,由幾何關(guān)系,新位置下的偏心率、偏位角:
圖1 差分法求油膜剛度、阻尼
+△x擾動得Fx1、Fy1;-△x擾動得Fx2、Fy2;
同理,+△y擾動得Fx3、Fy3;-△y擾動得Fx4、Fy4;
+△Vx擾動得Fx5、Fy5;-△Vx擾動得Fx6、Fy6;+△Vy擾動得Fx7、Fy7;-△Vy擾動得Fx8、Fy8。
據(jù)上,用差分概念計算剛度、阻尼實質(zhì)上是求解八次雷諾方程,借助Matlab數(shù)值計算軟件進(jìn)行編程可計算出無量綱剛度系數(shù)阻尼系數(shù)。
取不同寬徑比、偏心率,得到無量綱剛度系數(shù)(相對單位為μωL/ψ^8)、無量綱阻尼系數(shù)(相對單位為μL/ψ^8)隨偏心率、寬徑比的變化曲線。
圖2 無量綱油膜剛度系數(shù)-偏心率-寬徑比變化關(guān)系
圖2、3分別為圓柱瓦的無量綱剛度、阻尼系數(shù)隨偏心率和寬徑比的變化曲線,由圖可見剛度、阻尼系數(shù)隨偏心率增大而增大,其增大速度隨軸承結(jié)構(gòu)、參數(shù)而異。軸承的剛度、阻尼系數(shù)是研究軸系振動和穩(wěn)定性的重要性能參數(shù),它不僅是計算軸承失穩(wěn)轉(zhuǎn)速的依據(jù),而且是計算軸承臨界轉(zhuǎn)速的邊界條件。