王 偉 王林林 吳倩文

(青島科技大學(xué)橡塑材料與工程教育部重點實驗室,山東省橡塑材料與工程重點實驗室 山東青島 266042)

橡膠油封被廣泛應(yīng)用于機械設(shè)備中，適于旋轉(zhuǎn)和往復(fù)運動密封[1]。橡膠油封在使用中涉及到固體力學(xué)、摩擦學(xué)、流體力學(xué)、邊界潤滑等多方面的理論知識，因而對橡膠油封在安裝和使用中的變形及密封界面上接觸壓力分布、界面流體壓力[2]和唇口附近溫度場的精確計算，在理論上存在較大的困難[3]。而數(shù)值模擬可對橡膠油封的上述問題進(jìn)行準(zhǔn)確分析[4-7]。

油封與高速旋轉(zhuǎn)軸之間因摩擦生熱，會導(dǎo)致唇口附近溫度很高，并會加速橡膠的老化和磨損[8]，使其服役壽命縮短[9-10]。因此，準(zhǔn)確預(yù)測油封唇口附近的溫度場具有重要實際意義。LIU等[11]研究了表面帶有不同紋理結(jié)構(gòu)的軸對油封磨耗和密封性能的影響。馬靈童和孟慶睿[12]用FLUENT建立金屬薄盤結(jié)構(gòu)摩擦副熱流固耦合模型，對其生熱過程進(jìn)行了仿真分析。張付英等[13]建立了油封唇口的混合流體潤滑模型，利用MATLAB求解得到了不同速度下唇口溫度的變化。孫煊廣等[14]用ABAQUS建立風(fēng)電高速軸制動器熱-機耦合三維有限元模型，對其摩擦生熱進(jìn)行了分析，獲得了較好的效果，可借鑒用于橡膠密封件摩擦生熱的研究。雖然許多文獻(xiàn)對油封的變形、磨耗和溫度場進(jìn)行了仿真分析和實驗研究，但采用熱力耦合方法進(jìn)行模擬并進(jìn)行實驗驗證的工作尚不多。目前相關(guān)文獻(xiàn)研究存在的問題是單獨進(jìn)行油封變形的力場分析或流固耦合溫度場的計算分析，針對干摩擦情況下油封熱力完全耦合的研究很少。本文作者僅用軸對稱模型進(jìn)行熱力耦合計算，建模方便且計算效率高。另外，考察干摩擦狀況下油封的摩擦生熱和磨損情況，對油封在環(huán)境惡劣情況下的服役壽命評估有指導(dǎo)意義。文中研究考慮了橡膠材料的非線性力學(xué)行為，油封唇口的大變形和與軸的接觸非線性邊界條件，利用力場-溫度場完全耦合的方法對油封唇口的摩擦生熱過程進(jìn)行了數(shù)值模擬，并與油封工裝試驗結(jié)果進(jìn)行了比較和驗證，二者吻合較好，為解決類似旋轉(zhuǎn)密封件的動態(tài)摩擦生熱問題提供了一種有效途徑。

1 有限元模型建立

1.1 幾何模型

根據(jù)36 mm×52 mm×12 mm耐油丁腈橡膠油封的實際斷面尺寸及相應(yīng)的裝配軸尺寸，在ABAQUS軟件中建立有限元軸對稱模型，模型共有3 630個節(jié)點，3 364個單元。軸對稱有限元模型如圖1所示。

圖1 軸對稱有限元模型Fig 1 Axisymmetric finite element model

1.2 材料模型

對于橡膠類物理非線性材料的力學(xué)行為采用Yeoh模型[15]來描述，如式(1)所示。

W=C10(I1-3)+C20(I1-3)2+C30(I1-3)3

(1)

其中Rivlin系數(shù)C10、C20、C30由橡膠的單向拉伸試驗和最小二乘法來確定。文中丁腈橡膠的參數(shù)C10、C20、C30分別為0.756、-0.206、0.073 MPa。實驗測試得到丁腈橡膠的導(dǎo)熱系數(shù)為0.35 W/(m·K)，比熱容為1 600 J/(kg·K)。金屬軸的導(dǎo)熱系數(shù)為80 W/(m· K)，比熱容為460 J/(kg·K)。

1.3 邊界條件

油封與軸的接觸采用罰函數(shù)法求解。在裝配油封時在軸表面和油封唇部涂上一層潤滑油，軸的表面粗糙度約為Ra0.4 μm，因此摩擦因數(shù)取0.18是合理的。測試中沒有油壓，也無油潤滑，油封處于干摩擦狀態(tài)。在ABAQUS中采用摩擦生熱功能，考慮油封和軸與外界的對流換熱和熱輻射邊界條件。為了便于紅外測溫準(zhǔn)確，軸端部涂抹成黑色。實驗室溫度為22 ℃，試驗初始溫度與室溫一致。試驗時在很短時間內(nèi)將軸逐漸加速到穩(wěn)定轉(zhuǎn)速1 000 r/min。

1.4 數(shù)值模擬方法

計算采用瞬態(tài)力-熱完全耦合方法，首先模擬軸安裝到油封的裝配過程，預(yù)先將金屬骨架一側(cè)固定，對軸施加一定位移到達(dá)裝配位置；然后在軸的耦合參考點(選取在旋轉(zhuǎn)軸線上)施加104.72 rad/s(相當(dāng)于1 000 r/min)轉(zhuǎn)速，根據(jù)實際情況大約10 s達(dá)到要求的轉(zhuǎn)速；最后軸在恒定轉(zhuǎn)速下連續(xù)轉(zhuǎn)動25 min，使得溫度場達(dá)到平衡，整個模擬過程完全與試驗一致。

在ABAQUS中計算油封摩擦生熱時，油封唇部與軸之間產(chǎn)生的摩擦功率作為熱傳導(dǎo)方程中的內(nèi)熱源(生熱率)，摩擦生熱功率[10]的計算則由唇口摩擦力與旋轉(zhuǎn)軸的表面線速度得到，從而實現(xiàn)力場與溫度場的耦合迭代求解，最后以溫度場達(dá)到收斂誤差0.1 ℃為準(zhǔn)。摩擦功率的50%分配給橡膠油封，余下的分配給軸作為傳熱計算的內(nèi)熱源。在傳熱過程的計算中考慮了油封、軸外表面與外界環(huán)境的對流和輻射，忽略了材料熱膨脹對力場的影響以及橡膠的黏彈性特征。

2 實驗裝置

采用油封測試工裝試驗機，通過傳動軸后側(cè)設(shè)置的扭矩傳感器可得到油封轉(zhuǎn)動的摩擦扭矩，通過功率模塊換算得到消耗的功率，通過變頻電機調(diào)節(jié)旋鈕控制軸的轉(zhuǎn)速并以數(shù)字顯示，用紅外成像儀測量軸端部邊緣的溫度。測試工裝如圖2所示，圖中顯示了轉(zhuǎn)速穩(wěn)定后的摩擦扭矩和轉(zhuǎn)速以及功率值。

圖2 油封測試工裝試驗機Fig 2 Lip seal testing device

3 結(jié)果與討論

3.1 油封主要性能參數(shù)比較及接觸壓力分布

油封裝配到軸上后模擬得到的抱軸力為22.81 N，速度達(dá)到1 000 r/min后的摩擦扭矩為452.75 N·mm。試驗測試得到的摩擦扭矩約為448 N·mm (見圖2)，利用公式T=πDμFr(式中D為軸直徑，μ為摩擦因數(shù)，F(xiàn)r為抱軸力)計算得到的摩擦扭矩T為464.35 N·mm，三者基本吻合。

在模擬油封與軸的裝配過程中，先將軸外移與油封唇口分開，然后逐漸施加位移到達(dá)裝配位置。該分析步完成后可以得到油封唇口的接觸壓力分布如圖3所示。可以看出，由于油封唇部與軸的過盈量較大導(dǎo)致尖角附件接觸壓力最高約為2.59 MPa，然后沿著接觸界面向下逐漸減小(見圖4唇口)；接觸寬度約為0.23 mm，唇口與軸接觸界面上的平均接觸壓力約為0.94 MPa，大于油封工作時的油壓0.055 MPa，說明密封可靠。因此，接觸界面上的平均壓力分布可用于油封唇部密封可靠性的評價。

圖3 油封與軸界面上的接觸壓力分布Fig 3 Distribution of contact pressure between lip seal and shaft

圖4 油封唇口的接觸壓力分布及軸裝配前后變形輪廓比較Fig 4 Distribution of contact pressure in the vicinity of lip andcomparison of deformed profile after shaft applied

圖4進(jìn)一步顯示了軸裝配前后油封的變形輪廓比較(虛線為軸裝配前的輪廓)，可以看出唇口部分變形較大，而腰部變形很小，抱軸力主要由彈簧力和腰部以上橡膠的彈性變形力組成。仿真得到的變形輪廓變化以及唇口壓力分布為合理設(shè)計和優(yōu)化油封結(jié)構(gòu)有重要指導(dǎo)意義。

3.2 摩擦生熱溫度場

油封唇部在高速旋轉(zhuǎn)軸的作用下，由于局部摩擦功較高轉(zhuǎn)化為熱量耗散，但橡膠的導(dǎo)熱系數(shù)較小，從而使得唇部溫度逐漸升高。經(jīng)過一段時間運行后，唇口生熱與外部環(huán)境的對流和輻射散熱達(dá)到平衡，整個溫度場達(dá)到穩(wěn)態(tài)，此時的溫度場分布如圖5所示。可以看到唇部溫度較高，而其他部位溫度較低。唇口最高溫度達(dá)到140.6 ℃，與紅外成像儀測試的最高溫度141 ℃接近，如圖6所示?？梢姡M得到的油封摩擦扭矩和唇口溫度與試驗結(jié)果吻合很好。

圖5 油封溫度場達(dá)到穩(wěn)態(tài)時的溫度分布Fig 5 Steady temperature field of the lip seal

圖6 紅外成像儀得到的溫度場分布Fig 6 Temperature field tested by infrared imager(a) wholetemperature distribution;(b) temperaturedistribution of local seal and shaft

4 結(jié)論

(1)通過建立橡膠油封在干摩擦狀態(tài)下力場-溫度場耦合的有限元模型，預(yù)測了油封的抱軸力和摩擦扭矩，以及穩(wěn)態(tài)下的摩擦生熱溫度場，并與相應(yīng)油封工裝測試結(jié)果和紅外成像儀得到的溫度進(jìn)行了比較，結(jié)果表明:力-熱耦合方法求解油封的摩擦生熱溫度與試驗結(jié)果吻合較好，驗證了文中提出方法的正確性；模擬得到的油封摩擦扭矩與試驗結(jié)果接近，進(jìn)一步證實所采用計算方案的合理性。

(2)文中初步建立了橡膠唇形油封的力-熱耦合有限元模型，得到了摩擦生熱的穩(wěn)態(tài)溫度場，但在計算中尚未考慮油封唇口的磨耗對接觸壓力和摩擦生熱的影響，以及溫度對橡膠老化后性能和熱膨脹的影響，進(jìn)一步的工作還需建立更加完善的模型并開展相應(yīng)的試驗來加以驗證，才能更準(zhǔn)確地預(yù)測橡膠油封的服役壽命和密封的可靠性。