侯曉帥　趙則祥　于賀春　張國(guó)慶　何　鑫　姚　博

(中原工學(xué)院機(jī)電學(xué)院，河南 鄭州 451191 )

?

幾何誤差對(duì)多孔質(zhì)氣靜壓軸承性能的影響*

國(guó)家自然科學(xué)基金資助項(xiàng)目(51475485)；國(guó)家自然科學(xué)基金資助項(xiàng)目(51405523)；鄭州市科技攻關(guān)項(xiàng)目(121PPTGG363-6)

侯曉帥趙則祥于賀春張國(guó)慶何鑫姚博

(中原工學(xué)院機(jī)電學(xué)院，河南 鄭州 451191 )

針對(duì)幾何誤差對(duì)多孔質(zhì)氣靜壓軸承性能的影響，采用基于有限體積法的Fluent軟件對(duì)流場(chǎng)進(jìn)行三維建模仿真。研究了幾何誤差對(duì)多孔質(zhì)氣靜壓軸承承載力和剛度的影響，得到以下結(jié)論：幾何誤差形成的氣膜間隙變化造成的壓力分布不均，當(dāng)中凸-直線、中凸-圓弧和圓錐度誤差在3 μm以內(nèi)時(shí)，對(duì)氣靜壓軸承的性能有明顯的提高；沿徑向方向上傾斜3 μm以內(nèi)時(shí)，對(duì)氣靜壓軸承的性能影響不大；中凹和橢圓度誤差對(duì)氣靜壓軸承的性能有負(fù)作用。

氣靜壓軸承；幾何誤差；承載力；剛度

多孔質(zhì)材料作為氣靜壓軸承的節(jié)流器具有廣泛的應(yīng)用范圍，主要原因是多孔質(zhì)材料節(jié)流器供氣面積大，穩(wěn)定性能好，具有較高的承載能力和剛度等優(yōu)點(diǎn)[1]。目前，多孔質(zhì)材料有青銅、不銹鋼、陶瓷和石墨等。青銅和不銹鋼材料的塑性較大，用機(jī)械加工方法難以獲得不被堵塞的多孔質(zhì)表面。經(jīng)磨削加工的陶瓷孔隙內(nèi)極易出現(xiàn)脆裂的微粒。采用冷等靜壓的工藝方法制造出的各向同性多孔質(zhì)石墨是比較理想的材料，它不僅具備一般石墨的全部特性而且材料結(jié)構(gòu)均勻且致密、孔隙分布均勻、抗氧化能力強(qiáng)[2]。但是多孔質(zhì)石墨由于質(zhì)軟，加工過(guò)程中不利于保證幾何精度。

吳笛[3]建立了局部多孔質(zhì)氣靜壓徑向軸承的數(shù)學(xué)模型并研究了各參數(shù)對(duì)軸承靜態(tài)特性的影響；李歡歡等人[4]研究了橢圓誤差對(duì)多孔質(zhì)氣靜壓軸承特性的影響；李樹森等人[5]理論推導(dǎo)了幾何形誤差對(duì)氣靜壓圓柱軸承運(yùn)動(dòng)精度的影響關(guān)系式，提出設(shè)計(jì)時(shí)應(yīng)當(dāng)考慮幾何誤差的影響，但未給出詳細(xì)的設(shè)計(jì)指導(dǎo)；徐建民等人[6]應(yīng)用邊界元法研究了在二維流動(dòng)假設(shè)下，軸承副具有幾何形狀誤差時(shí)空氣靜壓軸承的工作特性；侯志泉[7]系統(tǒng)研究了工作參數(shù)和結(jié)構(gòu)參數(shù)以及靜壓軸承摩擦副形狀誤差對(duì)液體靜壓軸承性能的影響規(guī)律；邊新孝等人[8-9]考慮氣膜厚度誤差對(duì)小孔節(jié)流徑向軸承承載能力的影響；采用數(shù)值計(jì)算的方法，分析了小孔節(jié)流徑向軸承支點(diǎn)反力不在中點(diǎn)或在中點(diǎn)且有彎矩時(shí)偏載情況的影響；Pande S S等人[10]利用有限差分法分析了幾何形誤差和表面粗糙度對(duì)小孔節(jié)流徑向軸承靜態(tài)特性的影響；James Nevel[11]研究了定子和轉(zhuǎn)子的n棱圓誤差對(duì)氣靜壓軸承性能的影響，分析了由于氣膜的彈性作用彌補(bǔ)了n棱圓誤差的缺陷，n棱圓誤差對(duì)氣靜壓軸承性能的影響可以忽略。綜合國(guó)內(nèi)外的研究現(xiàn)狀，筆者通過(guò)分析幾何誤差形成的氣膜間隙變化造成的壓力分布情況，研究了多孔質(zhì)氣靜壓主軸的圓度、圓柱度和傾斜度等幾何誤差對(duì)氣靜壓軸承承載力和剛度的影響。

1　模型建立

1.1物理模型的建立

建立的多孔質(zhì)氣靜壓徑向軸承副模型如圖1所示。設(shè)定多孔質(zhì)徑向軸承的主要參數(shù)如表1所示。理想狀態(tài)下的具體結(jié)構(gòu)參數(shù)：直徑d=50 mm，長(zhǎng)L=80 mm，單側(cè)氣膜厚度h=10 μm。

表1主要參數(shù)表

符號(hào)說(shuō)明符號(hào)說(shuō)明Ps/MPa供氣壓力Pa/MPa環(huán)境壓力D/mm軸承外徑d/mm主軸直徑a/mm橢圓實(shí)軸長(zhǎng)b/mm橢圓虛軸長(zhǎng)L/mm軸承長(zhǎng)度Δx/mm主軸偏心量e橢圓離心率α傾斜角

根據(jù)國(guó)家標(biāo)準(zhǔn)GB/T1182-2008《產(chǎn)品幾何技術(shù)規(guī)范(GPS)幾何公差 形狀、方向、位置和跳動(dòng)公差標(biāo)注》的規(guī)定，幾何誤差包括形狀公差、方向公差、位置公差和跳動(dòng)公差。形狀公差主要包括圓度和圓柱度等[12]。在實(shí)際生產(chǎn)中，加工出來(lái)的氣靜壓徑向軸承(孔)與其相配的軸不僅存在尺寸偏差而且存在諸如圓度、圓柱度等幾何誤差，造成氣膜厚度在軸向和徑向都是不均勻的。同時(shí)，多孔質(zhì)氣靜壓軸承與主軸在裝配過(guò)程中會(huì)產(chǎn)生傾斜度誤差。

為便于理論分析，采用簡(jiǎn)化了的主軸幾何誤差模型計(jì)算多孔質(zhì)氣靜壓軸承的承載力和剛度。其中簡(jiǎn)化原則是：保證不改變幾何公差值的情況下，將多孔質(zhì)氣靜壓主軸表面不同波紋度的曲線通過(guò)最小二乘法擬合為直線或圓弧等形狀。

1.1.1圓度誤差

圓度誤差主要指橢圓度誤差包括沿實(shí)軸和虛軸兩種模型，如圖2所示。

1.1.2圓柱度誤差

主軸圓柱度誤差主要包括中凸、中凹型誤差和圓錐度誤差，其它形狀類型暫不討論。

(1)主軸中間凹、凸引起的誤差通過(guò)最小二乘法擬合為直線形和圓弧形兩種，如圖3～4所示。

(2)圓錐度誤差，造成氣膜厚度沿軸線方向的不均勻，如圖5所示。

1.1.3傾斜度誤差

裝配過(guò)程中主軸傾斜度誤差，如圖6所示。

1.2網(wǎng)格的劃分

通過(guò)分析，建立的物理模型包括多孔質(zhì)氣靜壓徑向軸承和氣膜間隙兩部分，考慮到流場(chǎng)尺寸與氣膜尺寸比例失調(diào)，計(jì)算精度及效率等問(wèn)題，采用分區(qū)劃分的方法[13]。對(duì)于氣膜間隙的厚度方向進(jìn)行加密處理，兩部分均采用結(jié)構(gòu)化網(wǎng)格。本次分析設(shè)定氣膜間隙網(wǎng)格數(shù)值為6，其余部分網(wǎng)格間距為1 mm，多孔質(zhì)石墨透氣率為18%，收斂精度為1×10-9。

1.3邊界條件的設(shè)置

計(jì)算有以下3個(gè)假設(shè)條件[2]：

(1)壁面絕對(duì)光滑，不考慮壁面粗糙度影響。

(2)不考慮轉(zhuǎn)速的影響。

(3)牛頓流體通過(guò)多孔介質(zhì)的流動(dòng)，假定為層流流動(dòng)并遵守Darcy定律，不考慮滑移邊界的影響。

根據(jù)以上假設(shè)條件，具體邊界條件設(shè)定如下：

(1)環(huán)境壓力Pa=0.1 MPa，環(huán)境溫度Ta=300 K。

(2)供氣壓力恒定Ps=0.5 MPa(絕對(duì)壓力)，供氣溫度Ts=300 K；出口壓力等于環(huán)境壓力Pa，溫度為環(huán)境溫度Ta。

1.4剛度的求解

多孔質(zhì)氣靜壓徑向軸承的承載能力隨著偏心量變化的曲線。對(duì)于靜態(tài)剛度Kw可以采用差分方法進(jìn)行計(jì)算[14]，計(jì)算公式如下：

(1)

式中：Kw為靜態(tài)剛度，N/μm；W為承載力，N；h為氣膜厚度，μm。

2　Fluent仿真結(jié)果分析

2.1幾何誤差Fluent分析

2.1.1橢圓度誤差對(duì)承載力和剛度的影響分析

不同離心率的橢圓度對(duì)承載力和剛度的影響如圖7～8所示。

2.1.2圓柱度誤差對(duì)承載力和剛度的影響分析

(1)主軸中間凹引起的誤差對(duì)承載力和剛度的影響如圖9～10所示。

在偏心量相同的情況下，隨著中凹誤差的增大承載力減?。粡姆抡娣治龅玫降臄?shù)據(jù)比較，中凹-圓弧的承載力(圖10a)下降速度比中凹-直線(圖9a)的快；

隨著中凹誤差的增大剛度減小且幅度較大。在相同中凹誤差的情況下，隨著偏心量的增大承載力增大；隨著偏心量的增大剛度小幅度波動(dòng)；當(dāng)中凹誤差大于1μm時(shí)，隨著偏心量的增大，剛度增加但幅度較小。中凹-直線(圖9)的承載力和剛度整體比中凹-圓弧(圖10)的承載力和剛度好；具有中凹-直線(圖9)和圓弧(圖10)等圓柱度誤差的氣靜壓軸承的承載力和剛度比理想情況下差。

(2)主軸中間凸引起的誤差對(duì)承載力和剛度的影響如圖11～12所示。

在偏心量相同的情況下，隨著中凸誤差的增大承載力小范圍內(nèi)增加；從仿真分析得到的數(shù)據(jù)比較，中凸-直線的承載力(圖11a)增長(zhǎng)速度比中凸-圓弧(圖12a)的快；隨著中凸誤差的增大剛度增大；當(dāng)中凸誤差大于3 μm時(shí)，剛度增長(zhǎng)幅度較小。在相同的中凸誤差情況下，隨著偏心量的增大承載力增大；隨著偏心量的增大剛度減?。划?dāng)中凸誤差大于3 μm時(shí)，剛度減小且幅度較大。中凸-直線(圖11)比中凸-圓弧(圖12)誤差整體的承載力和剛度好。具有中凸-直線(圖11)和圓弧(圖12)等圓柱度誤差的氣靜壓軸承的承載力和剛度比理想情況下的好。

(3)圓錐度對(duì)承載力和剛度的影響如圖13所示。

在圓錐度相同的情況下，隨著偏心量的增大承載力隨之增大；剛度隨之減小但幅度較小。在偏心量相同的情況下，隨著圓錐度誤差(圖13)的增大承載力和剛度均有不同程度的增加；當(dāng)圓錐度誤差大于3 μm時(shí)，剛度增長(zhǎng)發(fā)生劇變。由于氣靜壓徑向軸承副性能不受錐度(正錐或倒錐)的影響，在一定情況下可以允許有一定的圓錐度誤差，但應(yīng)盡量保持在3 μm以內(nèi)。

2.1.3傾斜度對(duì)承載力和剛度影響分析

裝配中產(chǎn)生的主軸傾斜度引起的位置誤差對(duì)承載力和剛度的影響，如圖14所示。

當(dāng)0≤α≤arctan(0.006/L)時(shí)，在傾斜度(圖14)相同的情況下，隨著偏心量的增大剛度減小但幅度較小。在偏心相同的情況下，隨著傾斜角度的變化承載力和剛度幾乎不變；沿徑向方向上傾斜在3 μm以內(nèi)時(shí)，對(duì)氣靜壓軸承的性能影響不大。

2.2原因分析

多孔質(zhì)氣靜壓主軸具有中凸圓柱度和圓錐度幾何誤差時(shí)，承載力和剛度反而增大主要是由于幾何誤差作用促使氣膜間隙沿軸向和徑向都不均勻而產(chǎn)生的壓力分布不均，如圖15所示。

圖15中，90°方向?yàn)闅饽らg隙較小一側(cè)；270°方向?yàn)闅饽らg隙較大一側(cè)；Z=0(圖15a)處為軸向中間位置；Z=-20(圖15b)和Z=20(圖15c)分別為軸向1/4和3/4處。氣體通過(guò)多孔質(zhì)氣靜壓軸承進(jìn)入孔/軸之間形成的氣膜間隙，90°方向壓力最大，270°方向壓力較小。理想模型壓力分布沿軸向在小范圍內(nèi)波動(dòng)，在軸向?qū)ΨQ位置的壓力分布相等，沿徑向形成的壓差約為供氣壓力的10%；中凸模型由于幾何誤差的作用，在軸向Z=0處時(shí)，氣膜間隙最小從而導(dǎo)致在徑向90°方向壓力最大；在軸向Z=-20和Z=20時(shí)，氣膜間隙逐漸增大，壓力分布逐漸減小，沿徑向形成的壓差約為供氣壓力的31%；圓錐模型(正錐)由于幾何誤差的作用，在軸向Z=-20處時(shí)，氣膜間隙較小，壓力分布較大；在軸向Z=20處時(shí)，氣膜間隙逐漸增大，壓力分布逐漸減小，沿徑向形成的壓差約為供氣壓力的34%；中凸和圓錐模型在軸向1/4和3/4處，沿徑向壓差比理想模型的大，所以中凸和圓錐模型承載力和剛度增加。由于圓錐模型沿徑向壓差比中凸模型的大，所以在相同圓柱度誤差值下，隨著偏心量的增大，圓錐模型剛度減小幅度比中凸模型小。最終，氣體沿著具有中凸圓柱度和圓錐度等幾何誤差的主軸線性流向兩端。

3　結(jié)語(yǔ)

本文研究了幾何誤差對(duì)多孔質(zhì)氣體靜壓軸承副的特性影響，得到如下結(jié)論：

(1)具有中凸和圓錐度誤差的氣靜壓軸承的承載力和剛度比理想情況下要好，可以允許有一定的中凸和圓錐度誤差，但應(yīng)盡量在3 μm以內(nèi)。

(2)具有橢圓度和中凹圓柱度誤差的氣靜壓軸承的承載力和剛度比理想情況下差；生產(chǎn)中應(yīng)盡量避免主軸有橢圓度和中凹誤差。

(3)沿橢圓實(shí)軸的氣靜壓軸承承載力和剛度比沿橢圓虛軸的好。

(4)沿徑向方向上傾斜3 μm以內(nèi)時(shí)，對(duì)氣靜壓軸承的性能影響不大。

[1]杜金名. 多孔質(zhì)流體靜壓軸承潤(rùn)滑技術(shù)的研究[D]. 哈爾濱：哈爾濱工業(yè)大學(xué)，2003.

[2]王云飛. 氣體潤(rùn)滑理論與氣體軸承設(shè)計(jì)[M]. 北京：機(jī)械工業(yè)出版社，1999.

[3]吳笛. 局部多孔質(zhì)氣體靜壓徑向軸承的建模與仿真[J]. 軸承，2010(10)：31-36.

[4]李歡歡，于賀春，張國(guó)慶，等. 橢圓誤差對(duì)多孔質(zhì)徑向氣體靜壓軸承特性的影響[J]. 制造技術(shù)與機(jī)床，2015(6)：84-87.

[5]李樹森，劉暾，張鵬順，等. 幾何形誤差對(duì)氣體靜壓圓柱軸承運(yùn)動(dòng)精度的影響[J]. 潤(rùn)滑與密封，2001(3)：45-48.

[6]徐建民，陳明. 形狀誤差對(duì)氣體靜壓圓柱軸承靜特性的影響[J]. 武漢化工學(xué)院學(xué)報(bào)，1998，20(1)：71-74.

[7]侯志泉. 基于CFD動(dòng)網(wǎng)格方法的液體靜壓主軸回轉(zhuǎn)精度機(jī)理研究[D]. 長(zhǎng)沙：湖南大學(xué)，2013.

[8]邊新孝，李謀渭，李威. 加工誤差對(duì)氣體靜壓徑向軸承的影響[J]. 北京科技大學(xué)學(xué)報(bào)，2005，27(3)：331-333.

[9]邊新孝，李謀渭. 靜壓空氣軸承偏載時(shí)的承載分析[J]. 航空精密制造技術(shù)，2005，41(1)：17-19.

[10]Pande S S, Somasundaram S. Effect of manufacturing errors on the performance of aerostatic journal bearings[J]. Wear, 1981，66(2)：145-156.

[11]James Nevel. The effect of stator and rotor form error on air bearing spindle accuracy[D]. US:The pennsylvania state university, 2011.

[12]毛平淮. 互換性與測(cè)量技術(shù)基礎(chǔ)[M]. 北京：機(jī)械工業(yè)出版社，2011.

[13]Liang Yingchun, Liu Jingshi. Computational fluid dynamics analysis of an aerostatic journal bearing with slot-entry restrictors[J]. Advanced Science Letters, 2011，4(8-10)：2817-2821.

[14]孫雅洲，盧澤生，饒河清. 基于FLUENT軟件的多孔質(zhì)靜壓軸承靜態(tài)特性的仿真與實(shí)驗(yàn)研究[J]. 機(jī)床與液壓，2007，35(3)：170-172.

(編輯譚弘穎)

如果您想發(fā)表對(duì)本文的看法，請(qǐng)將文章編號(hào)填入讀者意見(jiàn)調(diào)查表中的相應(yīng)位置。

Effect of geometric error on the performance of porous aerostatic bearing

HOU Xiaoshuai, ZHAO Zexiang, YU Hechun, ZHANG Guoqing, HE Xin, YAO Bo

(School of Mechanical & Electrical Engineering，Zhongyuan University of Technology, Zhengzhou 451191，CHN )

The effect of geometric error on the performance of porous aerostatic bearing is studied in this paper. The three-dimensional models of porous aerostatic bearing are built and the simulations are carried out with the Fluent software based on the finite volume method. The influence of the capacity and stiffness of porous aerostatic bearing are obtained; the distribution of uneven pressure caused by gas film clearance changes, error within 3μm of conicity and bulge in the middle can improvement the performance of porous aerostatic bearing; within 3μm of slope on the radial direction has a little influence; error of ovality and concave have a negative effect on the performance of porous aerostatic bearing.

aerostatic bearing；geometric error；loading capacity；stiffness

TH133.35

A

10.19287/j.cnki.1005-2402.2016.07.004

侯曉帥，男，1990年生，碩士研究生，主要研究方向?yàn)榫苤圃旒夹g(shù)與裝備。

2016-01-26)

160713