基于變步長迭代算法的微軸承參量特性研究

李家賢　張海軍　劉全

摘 要： 利用變步長迭代算法在Matlab平臺上對微軸承的修正模型進行仿真模擬。通過傳統連續模型、一階滑移修正及Wu新滑移修正模型的對比發現，在最小氣膜間隙處各滑移修正模型得到的參量遠低于連續模型的結果。而且當氣膜間隙小于1 μm時，各模型之間的結果偏差隨氣膜間隙的減小而迅速放大，連續模型已不適于此條件下的特性分析，必須考慮滑移修正效應的影響，以便得到較精確的相關電子設備中的參量特性。

關鍵詞： 變步長迭代算法； Matlab； 滑移修正模型； 微軸承

中圖分類號： TN622+.1?34； TH117.2 文獻標識碼： A 文章編號： 1004?373X（2017）16?0010?03

Abstract： The variable step size iterative algorithm is used to simulate the modification model of the micro bearing in Matlab platform. It is found by comparison with the traditional continuum model， first?order slip correction model and Wu new slip model that the parameters obtained by all the slip correction models at minimum gas film clearance are far below the results got by the continuum model. When the gas film clearance is less than 1 μm， the deviation of results from each model increases rapidly with reduction of the gas film clearance. It means that the continuum model is not suitable for characteristic analysis under this condition. That′s why the influence of slip correction effects must be considered， so as to achieve accurate parametric characteristics of aerostatic bearing in the related electronic equipment.

Keywords： variable step size iterative algorithm； Matlab； slip correction model； micro bearing

0 引 言

微軸承是一種具有高回轉精度，耐低溫和高溫及輻射等優良特性的精密電子元件。目前在PCB電子元件、精密電子設備、控制系統等領域中都得到了非常成功地應用[1?6]。但是現有微氣體軸承承載的測量估算都只采用建立在無滑移修正條件上的傳統連續模型。然而當氣膜間隙尺寸達到亞微米級時，氣體分子在固壁的滑移效應變得極其重要。在此條件下傳統模型已失效，需采用經過修正的模型進行參量的分析[7?8]。

本文引入一階滑移修正和Wu新滑移修正模型，然后在Matlab仿真平臺[9]下采用變步長迭代算法[10]求解修正后的模型，并探討修正后得到的精確參量特性。

1 修正改進模型

1.1 微軸承結構

圖1所示電子設備中微軸承的轉子直徑為D，軸長為L，有兩排沿周向均勻分布的節流孔，每排節流孔數為N，距最近軸端為l，孔徑為d。轉速為n，偏心距為e，氣膜厚度為h。當e=0時，平均氣膜間隙為hm。

1.2 滑移修正模型

圖1所示的θ為周向角度，z為軸向坐標。則傳統連續、一階修正、及Wu新修正模型的無量綱化微軸承模型為：

式中：[Z=zL]為無量綱軸向坐標；[H=hhm=1+εcosθ]為無量綱氣膜厚度；[ε=ehm]為偏心率；[P=ppa]為無量綱氣壓，p為氣壓，pa為標準氣壓；[Λ=μπn（r2）（5pah2m）]為壓縮系數，μ為氣體黏度；[Kn=λh]為當地克努森數，[λ]為平均自由程；[f=min[1Kn，1]]；bs，cs為模型系數，見表1。

表1 Reynolds方程模型系數

劃分網格時略去表面曲率的影響，并把節流孔當作一個網格點處理，如圖2所示，陰影部分為離散點的控制體積單元。則求解式（1）的邊界條件為：大氣邊界，[P（θ，Z=0）=1]；周期邊界，[P（θ，Z）=P（2π+θ，Z）]；對稱邊界，[（?P?Z）（Z=0.5）=0]。

2 變步長迭代原理及變量修正計算

假設加壓氣源流經節流孔為絕熱過程，則流出節流孔的流量為：

式中：[ρs]為氣體密度；[ps]為氣源壓力；k為氣體絕熱指數；第σ個節流孔后氣壓為[p=pσ]；[pσ]需根據圖3所示的流入區域流量[Min，σ]和流出區域流量[Mout，σ]相等的流量守恒原理確定。對第σ個節流孔區域有：

[Min，σ=Mσ2+Mσ3-Mσ1-Mσ4] （3）

式中：[Min，σ=24μ?TMin，σ（p2ah3m）]為通過第σ個節流孔流入的無量綱化流量；[?]為空氣的氣體常數；[Mσ1]，[Mσ2]分別為無量綱化周向流量；[Mσ3]，[Mσ4]分別為無量綱化軸向流量，需根據滑移速度分布[7?8]修正計算為（同理可得[Mσ2]及[Mσ3]）：

在圖4所示算法流程框圖中，使用SOR迭代算法[11]得到除孔后高精度的氣壓，再用變步長迭代算法求解孔后氣壓，最終對面積求積分得到由沿連心方向的WH和沿連心垂直方向的WV組成的剛度K為：endprint

[K=ΔW2H+W2VΔeWH=-20L202πpcosθ rdθdZWV=-20L202πpsinθ rdθdZ] （6）

3 結果與討論

所研究的電子元件中的微軸承參數為轉子半徑4 mm，軸長0.8 mm，節流孔距最近軸端距離為[L4]，節流孔數N=8，孔徑為0.05 mm，氣源壓力為4×105 Pa。圖5（a）～圖5（c）分別為在平均氣膜間隙hm=0.2 μm、偏心率ε=0.90，轉速n=12 kr/min條件下沿θ=0°處的連續模型、一階滑移和Wu新滑移模型的氣壓分布。可以發現各模型的氣壓在最小氣膜間隙處遠大于其余區域。這是因為此處的氣膜間隙遠小于其他區域，動壓效應較強顯著增大了氣壓。而遠離此處氣壓主要靠氣源壓力起支承作用相對較小。通過對比發現在最小氣膜間隙處，一階滑移、Wu新滑移修正模型的氣壓明顯小于連續模型且此條件下的最高偏差率分別達27.50%，36.17%，遠超可接受的誤差范圍。在此條件下，為了得到精確的測量值，應考慮滑移修正的影響。

圖6、圖7分別為在偏心率ε=0.95、轉速n=10 kr/min條件下，隨氣膜間隙的變化連續和各滑移修正模型的承載和剛度參量分布。可以發現當平均氣膜間隙相對較大時，承載較弱、增幅較慢，承載力差異不明顯。但當平均氣膜間隙hm小于1.0 μm時，承載、剛度增幅變快且結果偏差隨氣膜間隙的減小而迅速放大。這是因為動壓效應的增強使得承載變強、增幅變快，間接影響剛度的變化，此外在滑移效應的影響下最小氣膜間隙處的氣壓顯著降低。而大部分承載、剛度靠最小氣膜間隙周圍的氣壓生成，所以偏差隨著滑移效應的增強被放大。圖8、圖9分別為在平均氣膜間隙hm=0.2 μm、偏心率ε=0.95條件下，隨轉速的變化連續與各滑移修正模型的承載和剛度參量分布。可以發現在相對較低的轉速5～10 kr/min條件下，由于動、靜壓效應的相互作用，無滑移與各滑移模型的承載和剛度增幅相對較快。但當轉速大于10 kr/min后，動壓增幅作用減弱使承載的增幅變緩，間接使得剛度增幅變緩。然而各模型之間的承載及剛度偏差在一定范圍內隨轉速的增大有顯著偏差。

4 結 語

滑移修正效應對電子設備中微軸承承載、剛度等參量特性的影響隨氣膜間隙的變化而改變顯著；在小氣膜間隙條件下，各滑移修正模型得到的參量遠低于連續模型的結果，而且參量之間的偏差隨氣膜間隙的減小而迅速放大，必須考慮滑移效應的影響；通過變步長迭代算法求解滑移修正模型可得到較精確的微軸承參量特性分布，以便為相關電子設備、控制系統的安全運行提供精確數據及為相關參量的估算提供堅實的理論基礎。

參考文獻

[1] WANG C C， YAU H T， KUO C L. Bifurcation and chaos analysis of a relative short spherical air bearing system via a novel hybrid method [C]// 2009 International Workshop on Chaos Fractals Theories and Applications. [S.l.]： IEEE， 2009： 327?330.

[2] CHEN M F， HUANG W L， CHEN Y P. Design of the aerostatic linear guide way with a passive disk?spring compensator for PCB drilling machine [J]. Tribology international， 2010， 43（1）： 395?403.

[3] DE MOURA E P， SOUTO C R， SILVA A A， et al. Evaluation of principal component analysis and neural network performance for bearing fault diagnosis from vibration signal processed by RS and DF analyses [J]. Mechanical systems and signal processing， 2011， 25（5）： 1765?1772.

[4] 喻麗華，謝慶生，李少波，等.鉆削電主軸空氣靜壓徑向軸承潤滑參數優化[J].機械設計與制造，2014（8）：124?126.

[5] 張洪飛，閆守成.實時汽車電子輔助制動控制系統的設計研究[J].現代電子技術，2016，39（1）：153?156.

[6] 巫發茂，蔣龍，王健，等.基于ANSYS Workbench某機載電子設備隨機振動響應分析[J].現代電子技術，2016，39（10）：96?99.

[7] BURGDORFER A. The influence of the molecular mean free path on the performance of hydrodynamic gas lubricated bearings [J]. Transactions on ASMR， 1959， 81： 94?100.

[8] WU L. A slip model for rarefied gas flows at arbitrary Knudsen number [J]. Applied physics letters， 2008， 93（25）： 253103.

[9] 王若楠，田野，馬玉峰，等.Matlab程序單步調試方法及樣例分析[J].現代電子技術，2015，38（15）：139?141.

[10] ZHENG S F， JIANG S Y. Improved finite different method for pressure distribution of aerostatic bearing [J]. Journal of Southeast University， 2009， 25（4）： 501?505.

[11] 劉麗華，馬昌鳳，唐嘉.求解廣義鞍點問題的一個新的類SOR算法[J].計算數學，2016（1）：83?95.endprint