王學(xué)智，于天彪，袁娜，王宛山

(1.東北大學(xué)機械工程與自動化學(xué)院，遼寧沈陽 110819;2.沈陽工學(xué)院，遼寧撫順 113122)

超高速磨削用液體動靜壓軸承承載特性分析

王學(xué)智1，于天彪1，袁娜2，王宛山1

(1.東北大學(xué)機械工程與自動化學(xué)院，遼寧沈陽 110819;2.沈陽工學(xué)院，遼寧撫順 113122)

針對超高速磨削用液體動靜壓軸承支撐形式復(fù)雜的問題，應(yīng)用FLUENT軟件對液體動靜壓軸承進行了CFD分析，獲得了不同輸入狀態(tài)下的油膜壓力場分布，并在此基礎(chǔ)上對承載特性進行分析，得出了軸承承載力隨輸入?yún)?shù)變化的規(guī)律，運用正交試驗分析了各輸入?yún)?shù)的顯著性。結(jié)果表明:承載力與偏心率和供油壓力呈近似線性關(guān)系;承載力與轉(zhuǎn)速呈非線性關(guān)系，低速時隨轉(zhuǎn)速的提高，承載力基本保持不變，中速時隨轉(zhuǎn)速的提高承載力增大，高速時隨轉(zhuǎn)速的進一步提高，承載力有下降的趨勢;供油壓力對承載力的影響比轉(zhuǎn)速更為顯著。

CFD;壓力場;承載特性;正交試驗

超高速加工是德國 Carl.J.Salomon博士提出的，他預(yù)言了超越Talor切削方程式的非切削工作區(qū)域的存在，當(dāng)大幅度提高切削速度進行加工時，可以越過切削過程產(chǎn)生的高溫死谷而使刀具在超高速區(qū)進行高速切削，這時的切削力變小，從而使切削溫度大幅度降低，超高速磨削技術(shù)就是在此基礎(chǔ)上發(fā)展起來的。超高速磨削是一種高效、高精的現(xiàn)代加工技術(shù)，一般以磨削速度的大小來劃分，當(dāng)磨削速度vS≥150 m/s時，認(rèn)為達到了超高速磨削狀態(tài)。速度是實現(xiàn)和影響超高速磨削的直接因素，速度的快慢直接影響著加工的效率和精度，對于超高速磨削機床而言，其性能的優(yōu)劣很大程度上取決于的主軸系統(tǒng)的性能，而主軸系統(tǒng)的性能很大程度上由軸承性能決定。超高速磨床主軸系統(tǒng)多采用陶瓷軸承、磁浮軸承、空氣靜壓軸承和液體動靜壓軸承進行支撐。在高速重載機床多采用液體動靜壓軸承，它綜合了靜壓軸承和動壓軸承的特點，具有磨損小、承載能力大、使用壽命長、速度范圍寬、動態(tài)特性好和剛度高等突出優(yōu)點。本文作者以東北大學(xué)超高速磨削實驗臺裝配的液體動靜壓軸承為研究對象，分析動靜壓軸承的承載特點及承載特性，更好的掌握其承載規(guī)律，為進一步提高動靜壓承載能力，提高磨削試驗臺性能打下良好基礎(chǔ)［1－2］。

應(yīng)用FLUENT軟件，在考慮了潤滑油溫黏效應(yīng)的基礎(chǔ)上，采用湍流模型對軸承進行了仿真分析，獲得了不同輸入狀態(tài)下的油膜的壓力場分布，得出了主軸轉(zhuǎn)速、軸承偏心率、供油壓力和供油溫度對油膜承載的影響規(guī)律，并結(jié)合正交試驗，對輸入?yún)?shù)的顯著性進行了進一步分析。

1 理論分析

1.1 連續(xù)性方程

所有流動過程都必須滿足連續(xù)性方程。它與能量方程、動量方程或運動方程相結(jié)合，可求解各種流動問題。連續(xù)方程也就是質(zhì)量守恒方程。任何流動的問題都必須滿足質(zhì)量守恒定律。該定律可表述為:單位時間內(nèi)流體微元中質(zhì)量的增加等同于同一時間間隔內(nèi)流入該微元體的凈質(zhì)量。按照這一定律，可以得出質(zhì)量守恒方程:

式中:ρ為流體的密度;u、v、w是速度矢量在x、

y、z方向上的分量。

1.2 N-S方程

動量守恒定律也是任何流動系統(tǒng)都必須滿足的定律。該定律可表述為:微元體中流體的動量對時間的變化率等于外界作用在該微元上的各種力之和。該定律實際是是牛頓第二定律。按照這一定律，可導(dǎo)出x，y和z3個方向的動量守恒方程。

該方程首先是由Navier在1827年提出的，只考慮了不可壓縮流體的流動。Poisson在1831年提出可壓縮流體的運動方程。Saint-Venant在1843年，Stokes在1845年獨立地提出黏性系數(shù)為一常數(shù)的形式，現(xiàn)在都成為Navier-Stokes方程，簡稱N-S方程［5］。

可壓黏性流體的N-S方程:

式中:μ表示剪切變形率;ρfi代表微元體上的體力，i可取x、y、z時，分別表示在微元體x、y、z方向上所受體力。

1.3 能量方程

能量守恒定律是包含有熱交換的流動系統(tǒng)必須滿足的基本定律。該定律可表述為:微元體中能量的增加率等于進入微元體的凈熱流量加上體力與面力對微元體所做的功。該定律實際上是熱力學(xué)第一定律。

恒方程，不是很好用，所以一般扣除勢能的變化，得到包含內(nèi)能和動能的能量守恒方程。

式中:Cp表示比熱容，T表示溫度，K表示流體的傳熱系數(shù)，ST表示流體的內(nèi)熱源及由于黏性作用流體動能轉(zhuǎn)換為熱能的部分，簡稱為黏性耗散項。

1.4 溫黏方程

軸承在運行的過程中，潤滑油膜受到剪切作用和壓縮作用而產(chǎn)生熱量，由于這些熱量不能以對流和導(dǎo)熱的方式完全排走，而造成油膜溫度的升高，并產(chǎn)生徑向、軸向和油膜厚度方向的三維不均勻分布。油膜溫度升高會直接影響潤滑油的密度和黏度，從而使軸承的支撐特性發(fā)生改變。下面就從潤滑油的密度和黏度兩個方面來分析油膜的溫升效應(yīng)［6］。

Rost溫黏關(guān)系式:

式中:μx=0.18。

Rost溫黏關(guān)系式是在Vogel關(guān)系式的基礎(chǔ)上推導(dǎo)出來的，它分別以50℃和40℃的動力黏度為基礎(chǔ)，計算t溫度下的動力黏度。

2 軸承熱特性的CFD分析

2.1 建立油膜模型和劃分網(wǎng)格

以某大學(xué)超高速磨削試驗臺主軸裝配的五腔液體動靜壓混合動力軸承為研究對象，此軸承采用五腔供油、小孔節(jié)流、小腔承載。其結(jié)構(gòu)圖如圖1所示。

圖1 動靜壓軸承結(jié)構(gòu)示意圖

該動靜壓軸承結(jié)構(gòu)參數(shù)為:軸承直徑D=80 mm，軸承寬B=80 mm，平均油膜厚度h0=0.03 mm;小孔直徑dc=1.0 mm，高度為hc=3 mm;油腔深度為c=0.25 mm，軸向?qū)挾葹閎=36 mm，周向圓心角為δ=52°;節(jié)流油腔直徑為dr=2.5 mm，高度為hr=3 mm。

首先采用Gambit軟件建立油膜的有限元模型，然后對模型進行網(wǎng)格劃分及邊界類型的設(shè)定，從而建立油膜的網(wǎng)格模型。油膜三維模型和網(wǎng)格模型如圖2所示［3］。

圖2 三維模型和網(wǎng)格模型

2.2 選擇計算模型和設(shè)置參數(shù)

采用隱式定常模型進行計算，物理模型設(shè)置為紊流。流動介質(zhì)為2號主軸油，該主軸油的參數(shù)性能指標(biāo)分別為:潤滑油的密度ρ20=810 kg/m3;潤滑油的比熱容Cp=2 000 J/(kg·K);潤滑油的導(dǎo)熱系數(shù)λ= 0.37 W/(m·K);潤滑油的動力黏度μ的設(shè)置考慮了溫度的影響，應(yīng)用式 (4)進行計算求得不同的溫度下主軸油的動力黏度值，如表1所示［4］。

表1 潤滑油的動力黏度

2.3 分析轉(zhuǎn)速對軸承承載的影響

對于動靜壓軸承而言，轉(zhuǎn)速是其產(chǎn)生動壓承載的原動力，工作過程中主軸的高速旋轉(zhuǎn)，將具有一定黏度的潤滑油帶入楔形軸瓦中，使軸承產(chǎn)生了動壓效應(yīng)，從而提高了動靜壓軸承的支撐性能。在偏心率ε=0.5，供油溫度和供油壓力保持不變的情況下，改變軸頸的旋轉(zhuǎn)速度獲得軸承在不同轉(zhuǎn)速下油膜壓力場分布規(guī)律。其不同轉(zhuǎn)速下的壓力場分布如圖3所示，承載力隨轉(zhuǎn)速的變化曲線如圖4所示。

圖3 不同轉(zhuǎn)速下的油膜壓力場分布

圖4 軸承承載力隨主軸轉(zhuǎn)速變化的曲線

由圖3可以看出高壓區(qū)主要集中在偏心位置(最小油膜厚度)，并且隨著轉(zhuǎn)速的提高，高壓區(qū)向主軸轉(zhuǎn)動方向移動，這也反映出隨著主軸轉(zhuǎn)速的提高，軸承的動壓效應(yīng)逐漸增強;由圖4可知，轉(zhuǎn)速對軸承承載力的影響較為明顯，低速時隨轉(zhuǎn)速的提高，承載力趨于恒定，中速時，承載力與轉(zhuǎn)速近似呈線性關(guān)系，當(dāng)轉(zhuǎn)速進入高速區(qū)，隨著轉(zhuǎn)速的進一步提高，承載力有下降的趨勢。

2.4 分析偏心率對軸承承載的影響

動靜壓軸承在工作過程中，由于主軸系統(tǒng)受重力作用，使得軸承產(chǎn)生偏心，偏心的出現(xiàn)使軸承的周向?qū)ΨQ結(jié)構(gòu)發(fā)生改變，這時軸承內(nèi)部的壓力和溫度分布開始不均勻，從而產(chǎn)生了新的動態(tài)特性。在相同的條件下，改變軸承的偏心率，對模型分別進行求解計算，分別獲得軸承在偏心率ε為0.1～0.5時的，壓力場分析結(jié)果。軸承壓力場分布如圖5所示，承載力變化規(guī)律如圖6所示。

圖5 不同偏心率下的油膜壓力場分布

圖6 軸承承載力隨偏心率變化曲線

由圖5可知偏心的出現(xiàn)使得壓力分布出現(xiàn)了不均勻，最小油膜厚度處的壓力要高于其他區(qū)域，最大油膜厚度處的壓力要低于其他區(qū)域，并且隨著偏心率的逐漸增大油膜壓力分布的不平衡性逐漸增大;由圖6可知，偏心率對軸承的承載的作用明顯，在偏心較小時，軸承的承載與軸承偏心率近似呈線性關(guān)系，軸承的承載力隨偏心率的增大而增大。

2.5 分析供油壓力對軸承承載的影響

在其他因素不變的情況下，供油壓力的改變會導(dǎo)致軸承流量的變化，隨著壓力的提高，軸承流量增大，這樣由潤滑油所帶走的熱量增大，從而使軸承溫度減小。通過仿真分析獲得不同供油壓力下的軸承壓力場分布如圖7所示，其承載力變化規(guī)律如圖8所示。

圖7 不同供油壓力下的油膜壓力場分布

圖8 軸承承載力隨供油壓力變化的曲線

隨著供油壓力的增大，油膜壓力場普遍增大，壓力最高點仍集中在偏心位置;隨著供油壓力的增大，軸承承載力逐漸增大，與供油壓力呈近似線性關(guān)系。

3 基于正交試驗的軸承承載特性分析

3.1 正交試驗

軸承工作的輸入?yún)?shù)，對軸承的承載特性有很大影響。主要影響因子包括:主軸轉(zhuǎn)速、軸承的偏心率和供油壓力。下面將從這幾個主要影響因子著手，對軸承的承載特性進行正交實驗分析。

表2 因素水平表

表3 實驗結(jié)果

3.2 方差分析

為了了解各因素對實驗結(jié)果影響的顯著程度，對正交試驗結(jié)果進行方差分析。

總偏差平方和與總自由度ST為:

各列偏差平方和與自由度Sj為:

誤差偏差平方和與自由度Se為:

3.3 顯著性檢驗

應(yīng)用P值檢驗來確定哪個因子對結(jié)果影響顯著。設(shè)原假設(shè)為H0，T是檢驗統(tǒng)計量，其觀測值為t，H0的拒絕域為W，則稱如下定義的概率P為原假設(shè)的檢驗P值。

在統(tǒng)計實踐中，計算出P值后，將給定的顯著性水平α與P值比較，就可得出檢驗結(jié)論:

如果α＞P值，則在顯著性水平α下拒絕原假設(shè)。

如果α≤P值，則在顯著性水平α下接受原假設(shè)。

P值越小表明該影響因素越顯著。

應(yīng)用Matlab軟件對正交實驗結(jié)果進行方差分析，獲得計算結(jié)果如表4所示。

表4 承載力方差分析表

若取顯著性水平α=0.05，對于油膜承載力而言，其顯著性影響因子為偏心與供油壓力，供油壓力的作用效果比轉(zhuǎn)速的作用效果顯著。

4 結(jié)論

(1)通過對軸承輸入?yún)?shù)的仿真分析可知:承載力與偏心率和供油壓力呈近似線性關(guān)系，隨著偏心率和供油壓力的增大，承載力逐漸增大;承載力與轉(zhuǎn)速呈非線性關(guān)系，低速時隨轉(zhuǎn)速的提高，承載力基本保持不變，中速時隨轉(zhuǎn)速的提高承載力增大，高速時隨轉(zhuǎn)速的進一步提高，承載力有下降的趨勢。

(2)通過正交試驗的方差分析可知:對承載力的影響，各因素安從大到小的順序為:偏心率＞供油壓力＞主軸轉(zhuǎn)速。

(3)對于此動靜壓主軸而言，其靜壓效果要高于動壓效果，要進一步提高軸承的支撐效果需從提高供油壓力著手。

［1］沈桂蘭.超高速磨削加工技術(shù)發(fā)展及研究［J］.價值工程，2012(23):14－15.

［2］侯亞麗，李長河，盧秉恒.超高速磨削相關(guān)技術(shù)與工業(yè)應(yīng)用［J］.軸承，2009(3):51－56.

［3］于天彪，王學(xué)智，關(guān)鵬，等.超高速磨削機床主軸系統(tǒng)模態(tài)分析［J］.機械工程學(xué)報，2012，48(17):183－188.

［4］于天彪，王學(xué)智，關(guān)鵬，等.基于FLUENT的液體動靜壓軸承的動態(tài)特性分析［J］.潤滑與密封，2012，37(6): 1－5.

［5］XU Shu.Analytical Investigation on Static Hydrostatic Bearing by Fluent［J］.Advances in Electrical Engineering＆Electrical Machines，2011，134:497－504.

［6］姜繼海，張冬泉.參考溫度t=40℃時的粘溫方程［J］.潤滑與密封，1997(5):47.

Bearing Characteristics Analysis of Liquid Hybrid Bearing on Ultra-high Speed Grinding

WANG Xuezhi1，YU Tianbiao1，YUAN Na2，WANG Wanshan1
(1.School of Mechanical Engineering＆Automation，Northeastern University，Shenyang Liaoning 110819，China; 2.Shenyang Institute of Technology，F(xiàn)ushun Liaoning 113122，China)

Aimed at complex problem of the carrying of liquid hybrid bearing for ultra-high speed grinding，the oil film pressure field distributions in different input state were obtained through CFD analysis on the bearing by FLUENT software.On this basis，the bearing characteristics were analyzed，laws of bearing capacity changing along with the input parameters were obtained，and the significance of each input parameters was analyzed through orthogonal experiment.The results show that the relationships between bearing capacity and eccentricity as well as oil supply pressure are approximately linear，and the relationship between bearing capacity and speed is non-linear.In lower speed，along with speed increasing，bearing capacity remains unchanged basically，in moderate speed，along with speed increasing，bearing capacity is increased.In high speed，along with speed further improving，there is a downward trend in bearing capacity.The effect of oil supply pressure on the bearing capacity is more significant than that of speed.

CFD;Pressure field;Bearing characteristics;Orthogonal test

TH133.3

A

1001－3881(2014)9－061－5

10.3969/j.issn.1001－3881.2014.09.017

2013－04－17

國家自然科學(xué)基金資助項目 (51275084);遼寧省自然科學(xué)基金資助項目 (201102073)

王學(xué)智 (1985—)，男，博士研究生，研究方向為數(shù)字化設(shè)計與制造。E－mail:wangxuezhineu@126.com。