鄭龍燕

（山東冶金技師學(xué)院 濟(jì)南 250109）

高速切削技術(shù)的發(fā)展，對機(jī)械加工精度提出更高的要求。機(jī)械加工中的加工精度受到“機(jī)床-夾具-刀具-工件”工藝系統(tǒng)各個環(huán)節(jié)熱變形的影響。英國伯明翰大學(xué)的 J.Peclenik[1-2]調(diào)查表明，熱變形已成為影響機(jī)床加工精度的首要因素，占機(jī)床總誤差的40%～70%。因此，控制機(jī)床加工精度的關(guān)鍵在于控制關(guān)鍵部件的熱變形。

軸承-主軸系統(tǒng)是機(jī)床的重要組件，軸承運(yùn)轉(zhuǎn)過程中的摩擦生熱會使主軸產(chǎn)生熱變形，影響機(jī)床的加工精度。本文以某型號數(shù)控車床軸承-主軸系統(tǒng)為例建立有限元模型，進(jìn)行發(fā)熱量和熱邊界條件的計算。利用有限元軟件ANSYS對其進(jìn)行溫度場仿真，理論計算了主軸系統(tǒng)的溫升和熱變形，并通過搭建軸承-主軸系統(tǒng)熱伸長試驗臺，進(jìn)行實驗驗證。

1 軸承-主軸系統(tǒng)有限元模型

依據(jù)通用有限元分析流程，在三維軟件Solidworks 中建立幾何模型，并對其進(jìn)行以下簡化處理：

(1)刪除螺紋孔、油管、注油孔等細(xì)小特征；

(2)刪除倒角和圓角；

(3)利用簡單結(jié)構(gòu)的實體模型代替復(fù)雜結(jié)構(gòu)的部件，比如；用簡單圓環(huán)來代替主軸冷卻套；

(4)在保證整體結(jié)構(gòu)不變的情況下，可修改各模型的尺寸，保證各零件裝配正確。

通過對幾何模型進(jìn)行簡化，可提高網(wǎng)格劃分的質(zhì)量，節(jié)約仿真計算的時間，將簡化后的三維模型導(dǎo)入 ANSYS，如圖1所示。對軸承-主軸系統(tǒng)各部件的材料屬性進(jìn)行設(shè)置，將主軸材料設(shè)定為38CrMoAlA，軸承材料設(shè)定為 GCr15，其余零件設(shè)定為45 號鋼。

圖1 軸承-主軸系統(tǒng)示意圖

對簡化的軸承-主軸系統(tǒng)展開網(wǎng)格劃分的過程，網(wǎng)格劃分過粗會降低計算精度，劃分過細(xì)會增加計算時間，而且不同的網(wǎng)格劃分方法會影響網(wǎng)格的優(yōu)劣，因此，需要根據(jù)主軸零件的幾何尺寸選擇合理的網(wǎng)格單元尺寸，主軸部件的尺寸在幾毫米到幾百毫米，選擇網(wǎng)格尺寸為1～10 mm，對簡化后的軸承-主軸系統(tǒng)采用自動網(wǎng)格劃分，如圖2所示。

圖2 軸承-主軸系統(tǒng)網(wǎng)格劃分

2 熱載荷及邊界條件計算

2.1 熱載荷計算

軸承-主軸系統(tǒng)高速運(yùn)轉(zhuǎn)時，其內(nèi)部軸承因接觸摩擦而產(chǎn)生大量的熱，其發(fā)熱主要與軸承的摩擦力矩有關(guān)，Palmgren[3]推導(dǎo)出發(fā)熱量和摩擦力矩的關(guān)系式：Nr=1.05×104Mn。

軸承的摩擦力矩主要由以下三部分構(gòu)成：潤滑劑粘性摩擦力矩M0、載荷引起摩擦力矩M1和軸承自旋摩擦力矩M2。

1）潤滑劑粘性摩擦力矩

2）載荷引起的摩擦力矩

M1=f1P1dm

3）自旋摩擦力矩

式中：f0是與潤滑有關(guān)的系數(shù)，油霧潤滑取1；ν為運(yùn)動溫度下的潤滑油粘度；f1為軸承類型和載荷有關(guān)的系數(shù)； μ 為球與滾道接觸區(qū)的摩擦因數(shù)；Q為球與滾道的法向接觸載荷(N)； α為球與滾道橢圓接觸長半軸 (mm)； L為第二類橢圓積分。

2.2 對流換熱系數(shù)計算

對流是指因為流體的運(yùn)動，使流體各部分間發(fā)生相對位移，冷熱流體相互摻混所產(chǎn)生的熱量傳遞過程。對流換熱系數(shù)[4]是指流體與固體表面之間的換熱能力，即物體表面與附近空氣溫差 1℃，單位時間單位面積上通過對流與附近空氣交換的熱量。

1）軸承座表面熱對流換熱系數(shù)

軸承座的熱對流參數(shù)指靜止軸承座表面與空氣的自然對流換熱系數(shù)，經(jīng)驗值取9.7 W m2?℃

2）旋轉(zhuǎn)軸表面與空氣的對流換熱系數(shù)

軸承在實際運(yùn)轉(zhuǎn)過程中，產(chǎn)生的部分熱量通過熱傳導(dǎo)的方式傳遞給主軸，主軸在旋轉(zhuǎn)過程中與周圍空氣產(chǎn)生熱量交換，這種交換方式屬于強(qiáng)迫對流換熱，可以通過下式得到：

式中：?sf是對流換熱系數(shù)(W?m-2?℃-1)； Nu為努謝爾特數(shù)；Ka為主軸的熱傳導(dǎo)系數(shù) W?m-1?℃-1； ds為主軸直徑(m)；Re為而誒雷諾數(shù)；Pr為普朗特數(shù)；ω 為主軸旋轉(zhuǎn)角速度(rad?s-1)； 空氣的運(yùn)動粘度(m-2s-1)； 相應(yīng)流體的運(yùn)動粘度(m-2s-1)； 相應(yīng)流體的密度(kg?m-3)； v 為相應(yīng)流體的等壓比熱容； K為相應(yīng)流體的熱傳導(dǎo)系數(shù)（W?m-1?℃-1）。

3）軸承組件與潤滑劑的對流換熱

軸承在使用時需要添加潤滑劑以減小摩擦，運(yùn)轉(zhuǎn)過程中，內(nèi)部組件與潤滑劑之間存在熱量傳遞現(xiàn)象，這種熱量傳遞屬于強(qiáng)迫對流換熱，Harris[5]推薦軸承組件與潤滑劑間產(chǎn)生的強(qiáng)迫對流換熱系數(shù)為：

式中：K 為潤滑劑導(dǎo)熱系數(shù)(W?m-1?℃-1)；us為冷卻表面與潤滑劑間相對速度(m?s-1)； ls為特征長度(m)；v0為潤滑油粘度（Pa?s）。

對于滾動體表面與潤滑油之間的強(qiáng)迫對流換熱，式中us取保持架的表面速度，ls為滾道節(jié)圓直徑。對于內(nèi)、外圈滾道表面與潤滑油之間的強(qiáng)迫對流換熱，則us取保持架表面速度的三分之一。

不同轉(zhuǎn)速下的有限元模型邊界條件如表1所示，以轉(zhuǎn)速10 000 r/min為例，分析該轉(zhuǎn)速下，軸承-主軸系統(tǒng)的溫度場分布情況。

表1 軸承有限元模型邊界條件計算值

3 軸承-主軸系統(tǒng)熱態(tài)特性分析

3.1 軸承-主軸系統(tǒng)溫度場分析

圖3是轉(zhuǎn)速10 000 r/min 時，該軸承-主軸系統(tǒng)連續(xù)運(yùn)轉(zhuǎn)達(dá)到熱平衡時的溫度場仿真結(jié)果，圖 4 ～圖6是各組件溫度分布?？梢钥闯觯S承座溫度最高為36.5℃，主軸最高溫度為40.05℃，軸承-主軸系統(tǒng)中最高溫度為 43.05℃，該最高溫度發(fā)生于左側(cè)軸承中內(nèi)圈滾道與滾動體接觸的部位。

圖3 軸承-主軸系統(tǒng)整體溫度場

圖4 主軸溫度場

圖5 軸承座溫度場

圖6 軸承溫度場

3.2 主軸熱變形分析

在軸承-主軸溫度場仿真模型的基礎(chǔ)上，展開熱力耦合分析，建立軸承-主軸熱伸長模型，計算主軸熱伸長。對軸承-主軸系統(tǒng)溫度仿真分析結(jié)果進(jìn)行后處理，把主軸的溫度場數(shù)據(jù)導(dǎo)入主軸的熱伸長模型中，并給主軸施加轉(zhuǎn)速邊界條件，給軸承座施加固定約束條件，求解主軸的熱伸長。

圖7是轉(zhuǎn)速10 000 r/min時主軸達(dá)到熱平衡后的主軸的熱變形云圖，由圖中可以看出穩(wěn)態(tài)時，主軸單端端面熱伸長（X向）最大為12.587μm，而Y方向和Z方向主軸的熱變形最大分別為 4.1557μm和 4.131μm。與徑向熱變形相比，主軸軸向熱伸長更為嚴(yán)重。建立計算主軸熱伸長的仿真模型時，主軸左右端面均未固定，處于自由狀態(tài)，因此，該主軸的實際熱伸長為右端面熱伸長數(shù)據(jù)與左端面熱伸長絕對值求和，可到主軸在10 000 r/min時，軸向熱伸長為ΔL ，如式：

ΔL=L++L-=12.587+9.994=22.581 um

圖7 主軸熱伸長

4 軸承-主軸系統(tǒng)熱伸長試驗驗證

搭建高精度軸承-主軸系統(tǒng)熱伸長測試實驗臺，如圖8所示。高精度激光位移傳感器用磁力表座固定于實驗臺表面上，其測試精度為0.1 μm，通過主軸端面位移來測試主軸軸向熱伸長。由于高精度激光位移傳感器固定于實驗臺表面，因此，其測出的熱伸長為主軸絕對熱伸長。

圖8 軸承-主軸系統(tǒng)熱伸長測試實驗臺

圖9為主軸轉(zhuǎn)速為10 000 r/min時，3h內(nèi)主軸熱伸長數(shù)據(jù)，圖中黑線是原始信號，紅色是趨勢項，對實驗原始信號進(jìn)行趨勢項處理，發(fā)現(xiàn)熱伸長數(shù)據(jù)先上升后緩慢下降趨于平穩(wěn)的現(xiàn)象。實驗在0時刻，高精度激光位移傳感器的數(shù)值為0.05151 mm，實驗結(jié)束時刻，高精度激光位移傳感器的示數(shù)為0.0725 mm，因此計算出實驗測得的主軸軸向熱伸長的數(shù)據(jù)為。

有限元方法計算的主軸熱伸長為ΔL ，實驗測試的主軸熱伸長為計算值與實驗值的相對誤差ε：

圖9 實驗主軸熱伸長的原始信號及趨勢項

5 結(jié)語

通過對軸承-主軸系統(tǒng)的分析計算理論值以及主軸熱伸長測試的實驗值，表明有限元模型計算的正確性。軸承-主軸系統(tǒng)的溫度場表明，主軸左側(cè)軸承的溫度較高，在運(yùn)轉(zhuǎn)的過程中產(chǎn)生較大溫升，溫度的變化會導(dǎo)致主軸發(fā)生一定的熱變形，影響高速機(jī)床工作時所達(dá)到的加工精度。

在后續(xù)的仿真計算中，可以考慮改變潤滑條件、冷卻條件，選用合適的軸承支撐系統(tǒng)、誤差補(bǔ)償?shù)确绞浇档椭鬏S溫升，控制主軸變形量，改善高速數(shù)控機(jī)床的加工精度。