俞建衛(wèi)　王禮飛　魏　巍

合肥工業(yè)大學(xué)，合肥，230009

關(guān)節(jié)軸承摩擦溫度場(chǎng)計(jì)算方法

俞建衛(wèi)王禮飛魏巍

合肥工業(yè)大學(xué)，合肥，230009

以推力型關(guān)節(jié)軸承為研究對(duì)象，分析了軸向載荷下的接觸應(yīng)力，并根據(jù)試驗(yàn)所得摩擦力矩反推出了動(dòng)態(tài)摩擦因數(shù)。結(jié)合接觸應(yīng)力及摩擦因數(shù)求出了摩擦熱流率，并將其作為第二類邊界條件進(jìn)行了軸承非定常溫度場(chǎng)有限元分析，最后利用紅外熱像儀拍攝了摩擦過程中軸承側(cè)表面溫度分布，檢驗(yàn)仿真結(jié)果的可靠性。仿真與試驗(yàn)結(jié)果表明：該計(jì)算方法解決了球面摩擦副摩擦因數(shù)的精確計(jì)算和熱流率的空間分配問題，建模時(shí)忽略球面間隙會(huì)直接影響接觸應(yīng)力分布并造成溫度場(chǎng)分布出現(xiàn)較大誤差，考慮球面間隙可以更精確地計(jì)算軸承摩擦溫度場(chǎng)。

關(guān)節(jié)軸承；球面間隙；接觸應(yīng)力；溫度場(chǎng)；紅外測(cè)量；有限元

0　引言

關(guān)節(jié)軸承又稱球面滑動(dòng)軸承，是一種摩擦副為球面的特殊滑動(dòng)軸承，常在重負(fù)荷條件下做低速擺動(dòng)或旋轉(zhuǎn)，在軍用機(jī)械、建筑及航空等眾多領(lǐng)域均有大量應(yīng)用[1-3]。潤滑型關(guān)節(jié)軸承材料多數(shù)使用GCr15鋼,采用淬火后低溫回火處理[4],當(dāng)使用溫度高于回火溫度時(shí),材料硬度急劇下降,材料磨損加劇。另外,摩擦熱效應(yīng)還會(huì)顯著降低潤滑油膜的黏度,黏度的下降又會(huì)導(dǎo)致潤滑油膜的承載能力降低,甚至導(dǎo)致潤滑失效[5]。

研究表明：摩擦副表層或亞表層材料發(fā)生變形或撕裂時(shí)產(chǎn)生的能量大部分會(huì)轉(zhuǎn)化為熱能，引起摩擦表面溫度升高并向四周傳遞[6]。軸承摩擦過程中摩擦面作為生熱面，溫度會(huì)明顯高于外表面，傳統(tǒng)的測(cè)溫方法僅在軸承表面布置熱電偶測(cè)量單點(diǎn)溫度，對(duì)摩擦面及整體溫度場(chǎng)分布缺乏研究;又由于摩擦面和軸承內(nèi)部的不可見性，其溫度值不宜直接獲得。因此有必要對(duì)軸承展開摩擦三維非定常溫度場(chǎng)的分析。

在平面滑動(dòng)摩擦副的溫度場(chǎng)研究中,假設(shè)接觸區(qū)域壓力恒定，將理論計(jì)算與實(shí)測(cè)溫度對(duì)比可知,當(dāng)考慮熱流率隨時(shí)間和旋轉(zhuǎn)半徑變化時(shí)，可以得到更符合實(shí)際的摩擦溫度場(chǎng)模型[7]。但關(guān)節(jié)軸承的摩擦副為空間球面,其接觸形式不同于赫茲接觸，需要進(jìn)行數(shù)值方法求解[8]。又由于球面間隙的存在，因此假設(shè)球面完全協(xié)調(diào)接觸也會(huì)導(dǎo)致較大的計(jì)算誤差[9]。本文先使用有限元方法分析軸承接觸應(yīng)力分布,結(jié)合實(shí)測(cè)摩擦力矩反求動(dòng)態(tài)摩擦因數(shù)，根據(jù)摩擦生熱原理計(jì)算熱流率隨時(shí)間和空間坐標(biāo)分布,并將其作為摩擦面的第二類邊界條件建立溫度場(chǎng)模型,然后使用熱像儀測(cè)試溫度進(jìn)行試驗(yàn)驗(yàn)證,最后基于上述方法分別研究建模時(shí)忽略球面間隙和考慮球面間隙對(duì)軸承溫度場(chǎng)的影響。

1　試驗(yàn)裝置與實(shí)體模型

如圖1所示，以GX25S推力型油潤滑關(guān)節(jié)軸承為試驗(yàn)試樣,基于HDM-20型號(hào)端面磨損試驗(yàn)機(jī)設(shè)計(jì)試驗(yàn)夾具。試驗(yàn)過程中,由下平臺(tái)的加載面施加向上的軸向載荷;試驗(yàn)機(jī)通過銷槽傳遞扭矩,帶動(dòng)芯軸和內(nèi)圈做順時(shí)針單向旋轉(zhuǎn)運(yùn)動(dòng);軸承外圈與下平臺(tái)采用過渡配合，裝置在下平臺(tái)的牽引螺釘連接拉力傳感器，可實(shí)時(shí)測(cè)量軸承摩擦力矩；試驗(yàn)潤滑劑為鋰基潤滑脂。使用Thermo VisionTMA40-M在線紅外熱像儀拍攝軸承表面,獲得外表面實(shí)際溫度分布。

圖1　試驗(yàn)原理圖

有限元建模時(shí)忽略芯軸的銷槽、潤滑油槽及其他細(xì)小特征的影響，則整體模型關(guān)于芯軸中心對(duì)稱,為縮短計(jì)算時(shí)間將其簡化為平面軸對(duì)稱模型。

由于試驗(yàn)工況下軸承溫升較小,故假設(shè)材料各項(xiàng)參數(shù)在溫升范圍內(nèi)未發(fā)生變化。各材料參數(shù)如表1所示。

表1　材料屬性

2　有限元模型的建立

2.1模型建立及簡化

在軸承的實(shí)際結(jié)構(gòu)中,為滿足儲(chǔ)存潤滑油及滿足轉(zhuǎn)動(dòng)靈活的需要,外圈的球面直徑會(huì)略大于內(nèi)圈的球面直徑。使用ZEISS三坐標(biāo)儀分別測(cè)量內(nèi)外球面直徑,實(shí)測(cè)值如下:外圈內(nèi)徑為68.1000 mm,內(nèi)圈外徑為67.9400 mm,兩球面的曲率半徑差值為0.008 mm。數(shù)據(jù)顯示,內(nèi)外球面在幾何上存在初始的曲率半徑差值(下文中將這個(gè)差值統(tǒng)稱為球面間隙),而并非完全相同。但由于間隙值相對(duì)整體尺寸非常小,不利于有限元的建模處理，因此有研究在分析球面接觸應(yīng)力時(shí)忽略了球面間隙的存在,假設(shè)內(nèi)外球面完全一致,但試驗(yàn)證明忽略間隙值會(huì)導(dǎo)致應(yīng)力分析結(jié)果出現(xiàn)較大誤差[9]。為研究間隙值對(duì)接觸應(yīng)力及溫度場(chǎng)的影響,分別將建模時(shí)忽略間隙(內(nèi)外球面直徑均為68.000 mm)和考慮間隙(球面直徑為實(shí)測(cè)值)作為兩種不同情況進(jìn)行有限元建模研究,則對(duì)球面接觸可作如下假設(shè)：①忽略間隙時(shí)球面完全接觸，假設(shè)內(nèi)外球面接觸位置的溫度相同;②考慮間隙時(shí),由于球面間隙極小并涂有大量潤滑脂，內(nèi)外球面存在較大的熱交換，也假設(shè)內(nèi)外圈球面接觸位置的溫度相同。

基于上述假設(shè),在溫度場(chǎng)計(jì)算中使用MPC算法將接觸時(shí)對(duì)應(yīng)球面位置的溫度自由度“系”在一起。

2.2摩擦熱量的計(jì)算

摩擦生熱量可由下式計(jì)算[10]:

q=μp vr

(1)

式中,q為單位接觸面積上產(chǎn)生的熱量,即熱流密度;μ為摩擦因數(shù);p為接觸正壓力;vr為接觸位置的相對(duì)線速度。

由式(1)可知,摩擦熱的生成與摩擦因數(shù)、接觸應(yīng)力分布及相對(duì)線速度均有直接關(guān)系。

關(guān)節(jié)軸承試驗(yàn)機(jī)上的摩擦因數(shù)計(jì)算式為[11]

(2)

式中,M為摩擦力矩;N為載荷;dk為球面直徑。

但式(1)并沒有考慮加載方式和球面間隙等參數(shù)對(duì)接觸應(yīng)力及摩擦因數(shù)的影響,引入到熱流率的計(jì)算時(shí)會(huì)導(dǎo)致誤差較大。

采用等效離散法,在結(jié)構(gòu)有限元中對(duì)模型進(jìn)行網(wǎng)格劃分,可將球面等效成一系列不同半徑的圓環(huán)面接觸。假設(shè)在某一載荷條件下共有n個(gè)圓環(huán)參與實(shí)際接觸，忽略潤滑脂的影響，則各圓環(huán)對(duì)旋轉(zhuǎn)中心的摩擦力矩之和等于實(shí)測(cè)摩擦力矩。

試驗(yàn)過程中,設(shè)各圓環(huán)摩擦因數(shù)為μi(t),則由力矩平衡可得

(3)

為簡化模型,假設(shè)各圓環(huán)摩擦因數(shù)相同,且均為μ(t),則有

(4)

式中,M(t)為實(shí)測(cè)摩擦力矩;Fi為圓環(huán)接觸總反力;ri為圓環(huán)半徑,且ri∈[17.25 mm,27.701 mm]。

根據(jù)式(1),摩擦生熱量可表示為

(5)

式中,H(r,t)為熱流率;ω(t)為旋轉(zhuǎn)角速度;Ai為圓環(huán)面積。

對(duì)任一圓環(huán)有

Fi=piAi

(6)

因此,關(guān)節(jié)軸承的球面摩擦熱流率可表示為

Hi(r,t)=μ(t)Fiω(t)ri

(7)

2.3熱量傳導(dǎo)及定解條件

基于傅里葉導(dǎo)熱定律,對(duì)各向同性的固體導(dǎo)熱有如下方程[12]：

(8)

定解條件包括初始條件和邊界條件。試驗(yàn)在室溫下進(jìn)行，初始溫度為25 ℃。

對(duì)流換熱可根據(jù)牛頓冷卻公式計(jì)算：

Φ=hA(Tw-Tf)

(9)

式中,h為表面換熱系數(shù);A為表面面積;Tw和Tf分別為表面溫度和環(huán)境溫度。

熱輻射計(jì)算為非線性,增大了收斂難度,因此將傳熱系數(shù)提高10%替代熱輻射效果[13]。

2.4有限元計(jì)算

根據(jù)式(7)所需參數(shù),使用APDL語言編寫計(jì)算控制程序，程序流程如圖2所示。

圖2　APDL程序流程圖

試驗(yàn)工況為：載荷2 kN,轉(zhuǎn)速30 r/min。按照?qǐng)D2所示流程,程序先計(jì)算出兩類模型下的接觸應(yīng)力分布(圖3),再根據(jù)實(shí)測(cè)摩擦力矩計(jì)算出實(shí)時(shí)摩擦因數(shù)(圖4),結(jié)合接觸反力、摩擦因數(shù)及轉(zhuǎn)速計(jì)算出熱流率分布(圖5),最后計(jì)算出溫度場(chǎng)(圖6)。

3　結(jié)果與分析

3.1接觸應(yīng)力分布

圖3所示分別為兩類模型的球面接觸應(yīng)力徑向坐標(biāo)的變化情況。其中無間隙模型是假設(shè)內(nèi)外球面直徑均為68.000 mm,軸承內(nèi)外圈恰好處于完全協(xié)調(diào)接觸；有間隙模型根據(jù)實(shí)測(cè)球面的曲率半徑進(jìn)行建模處理，軸承內(nèi)外圈在初始狀態(tài)下并未接觸，由于載荷的施加，初始間隙被逐漸合并，并最先在球面底部產(chǎn)生初始接觸。

(a)軸承內(nèi)外圈有間隙

(b)忽略軸承內(nèi)外圈間隙圖3　接觸應(yīng)力

由圖3可知:兩類模型均在球面最底端出現(xiàn)最大接觸應(yīng)力，且均隨徑向坐標(biāo)的增大先急劇減小，而后下降趨勢(shì)趨于平緩。無間隙模型與文獻(xiàn)[9]得到了類似的接觸應(yīng)力分布趨勢(shì)，間隙模型接觸區(qū)域的徑向坐標(biāo)范圍遠(yuǎn)遠(yuǎn)小于無間隙情況，其最大接觸應(yīng)力也遠(yuǎn)遠(yuǎn)大于無間隙狀態(tài)。這是因?yàn)殚g隙模型中球面間隙破壞了球面接觸的協(xié)調(diào)性，初始接觸時(shí)只有球面最底端的區(qū)域參與接觸，導(dǎo)致其應(yīng)力分布與無間隙模型出現(xiàn)明顯差異。此外，由式(1)可知，接觸應(yīng)力的分布會(huì)直接影響摩擦生熱量在空間上的分布;由式(4)可知,接觸應(yīng)力的分布也同時(shí)影響摩擦因數(shù)的反推。因此建模時(shí)考慮實(shí)際的幾何間隙,可以提高熱流率的計(jì)算精度。

3.2摩擦因數(shù)及熱流率

圖4中,在忽略間隙和計(jì)入間隙兩種模型中,摩擦因數(shù)與摩擦力矩的變化趨勢(shì)均保持一致,但無間隙模型反推出的摩擦因數(shù)明顯小于有間隙模型反推出的摩擦因數(shù)。由于二者采用相同計(jì)算方法由相同摩擦力矩反推得到,無間隙模型中球面接觸區(qū)域較大，軸承轉(zhuǎn)動(dòng)不靈活易造成較大摩擦力矩，因此反推出的摩擦因數(shù)較有間隙模型反推出的摩擦因數(shù)偏小。

(a)摩擦力矩

(b)摩擦因數(shù)圖4　摩擦力矩及摩擦因數(shù)

圖5中兩類模型的熱流率均隨球面徑向坐標(biāo)和時(shí)間變化。有間隙模型中接觸面上只有球面最底端附近區(qū)域有摩擦熱生成，無間隙模型的熱流率分布范圍比有間隙模型的熱流率分布范圍更廣,但其峰值也明顯偏小。根據(jù)式(1)可知，在相同的轉(zhuǎn)速下接觸應(yīng)力分布和摩擦因數(shù)直接決定摩擦熱的分布,因此兩類模型的熱流率在空間上分布與接觸應(yīng)力分布成正比，在時(shí)間上的變化與動(dòng)態(tài)摩擦因數(shù)的變化成正比。

(a)無間隙狀態(tài)下熱流率分布

(b)間隙狀態(tài)下熱流率分布圖5　熱流率分布

3.3溫度分布及驗(yàn)證

圖6為軸承外圈在1800 s時(shí)的溫度灰度圖。由圖6可知:兩類模型的總體溫升大致相同，最高溫度都出現(xiàn)在球面底部；但二者的溫度場(chǎng)分布出現(xiàn)較大差異，無間隙模型的總體溫度分布更加均勻，摩擦生熱區(qū)域與軸承表面的溫差較間隙模型更小。出現(xiàn)此類現(xiàn)象的原因主要為：①式(7)是將摩擦阻力做功轉(zhuǎn)換成熱能，因此單位時(shí)間內(nèi)輸入軸承的總熱量相同，在相同的散熱邊界條件下,兩類模型的總體溫升也大致相同；②圖6中兩類模型的熱流率都集中于球面底部，故大量熱量輸入此處形成高溫區(qū)；③無間隙模型的接觸區(qū)域大于間隙模型，有利于球面上生成的摩擦熱向低溫區(qū)擴(kuò)散并形成較均勻的溫度分布，而間隙模型的生熱區(qū)域明顯偏小，不利于溫度擴(kuò)散，因此摩擦面的最高溫度與軸承表面出現(xiàn)更大的溫差。

(a)無間隙模型

(b)有間隙模型圖6　軸承外圈溫度分布

取軸承外圈表面距上表面3 mm處點(diǎn)A(圖1所示)作為溫度驗(yàn)證點(diǎn)，實(shí)測(cè)溫度與計(jì)算溫度隨時(shí)間變化曲線如圖7所示。

在摩擦前期,理論計(jì)算值與測(cè)量值都以較快速度上升并基本保持一致；隨著時(shí)間的延長溫升速度均有所減緩，無間隙模型溫度逐漸高于間隙模型溫度和實(shí)測(cè)溫度，實(shí)測(cè)溫度低于兩類模型的計(jì)算溫度且差距逐漸增大。在第30 min計(jì)算值與實(shí)測(cè)值出現(xiàn)最大誤差，實(shí)測(cè)值與無間隙模型的誤差約為6.75℃,實(shí)測(cè)值與間隙模型的誤差約為4.53℃。

圖7　驗(yàn)證點(diǎn)溫升過程對(duì)比

分析產(chǎn)生誤差的原因?yàn)?①建模過程中將幾何模型的細(xì)小特征作了簡化,同時(shí)也忽略了潤滑油和磨屑帶走的熱量，導(dǎo)致計(jì)算溫度偏高；②有研究指出金屬摩擦生熱過程中摩擦力所做功約有95%轉(zhuǎn)化為熱量[14],而式(7)假設(shè)摩擦力做功全部轉(zhuǎn)化為熱能也會(huì)造成理論計(jì)算溫度高于實(shí)測(cè)溫度;③無間隙模型的熱量流入范圍較大，熱源距離測(cè)溫點(diǎn)的距離也較間隙模型更近,導(dǎo)致A點(diǎn)溫升高于間隙模型的溫升，同時(shí)誤差也大于間隙模型的誤差。另外，由于球面間隙較小，高溫導(dǎo)致的熱變形會(huì)改變球面的接觸面積，進(jìn)而影響熱流率在空間上的分配，當(dāng)使用溫度較高時(shí)需要進(jìn)一步考慮熱力耦合的影響。

4　結(jié)論

(1)球面間隙的存在對(duì)球面接觸應(yīng)力分布、摩擦因數(shù)及摩擦熱的生成均有重要影響，建模時(shí)忽略間隙值會(huì)產(chǎn)生較大誤差。

(2)摩擦產(chǎn)生的熱流率隨時(shí)間和空間坐標(biāo)動(dòng)態(tài)變化,其中動(dòng)態(tài)的摩擦因數(shù)決定熱流率在時(shí)間上的變化規(guī)律，接觸應(yīng)力分布決定熱流率在空間上的能量分配。

(3)根據(jù)間隙模型的分析結(jié)果，當(dāng)軸承使用溫度在50 ℃左右時(shí)，內(nèi)部高溫區(qū)比軸承表面溫度高出約5℃，因此，在試驗(yàn)中通過測(cè)量表面溫度判斷軸承實(shí)際溫升的方法存在較大誤差。

(4)依據(jù)本文建立的間隙模型,在已知軸承尺寸、材料參數(shù)及換熱條件的情況下，結(jié)合試驗(yàn)所得摩擦力矩可以反推實(shí)時(shí)摩擦因數(shù)和熱流率分布。

[1]Kim B C,Park D C,Kim H S,et al.Development of Composite Spherical Bearing[J].Composite Structures,2006,75(1):231-240.

[2]潘漢明,陳以一,趙憲忠,等.廣州新電視塔雙向鉸節(jié)點(diǎn)試驗(yàn)研究[J].工業(yè)建筑,2009(1):117-121.

Pan Hanming,Chen Yiyi,Zhao Xianzhong,et al.Experimental Study on Bidirectional Pin-hinged Joints of Guangzhou New TV Tower[J].Industrial Construction,2009(1):117-121.

[3]Germaneau A,Peyruseigt F,Mistou S,et al.3D Mechanical Analysis of Aeronautical Plain Bearings:Validation of a Finite Element Model from Measurement of Displacement Fields by Digital Volume Correlation and Optical Scanning Tomography[J].Optics and Lasers in Engineering，2010,48(6):676-683.

[4]《機(jī)械工程材料性能數(shù)據(jù)手冊(cè)》編委會(huì).機(jī)械工程材料性能數(shù)據(jù)手冊(cè)[M].北京:機(jī)械工業(yè)出版社,1995．

[5]溫詩鑄, 黃平. 摩擦學(xué)原理[M].北京：清華大學(xué)出版社有限公司, 2002.

[6]Yang Jianheng.Analysis to the Thermal Mechanism of Wear in Dry Friction Condition[J].Lubrication Engineering,2005,171(5):173-176.

[7]魏巍,俞建衛(wèi),沈持正,等.基于試驗(yàn)數(shù)據(jù)的滑動(dòng)摩擦溫度場(chǎng)仿真方法研究[J].中國機(jī)械工程,2013,24(18):2426-2430.

Wei Wei,Yu Jianwei,Shen Chizheng,et al.Research on Sliding Friction Temperature Field Simulation Method Based on Experimental Data[J].China Mechanical Engineering,2013,24(18):2426-2430.

[8]Germaneau A,Peyruseigt F,Mistou S,et al.Verification of a Spherical Plain Bearing Finite-element Model Using Scattered Light Photoelasticity Tests[J].Proceedings of the Institution of Mechanical Engineers,Part J:Journal of Engineering Tribology,2008,222(5):647-656.

[9]劉六井,高偉.大型推力關(guān)節(jié)軸承結(jié)構(gòu)有限元分析[J].船海工程,2004(2):25-27.

Liu Liujing,Gao Wei.The Structure Finite Element Analysis of Large-scale Thrust Jiont Bearing[J].Ship and Ocean Engineering,2004(2):25-27.

[10]Yang Y C,Chen W L.A Nonlinear Inverse Problem in Estimating the Heat Flux of the Disc in a Disc Brake System[J].Applied Thermal Engineering,2011,31(14):2439-2448.

[11]王國鋒. 復(fù)合擺動(dòng)式關(guān)節(jié)軸承性能試驗(yàn)機(jī)的研制及試驗(yàn)研究[D].洛陽:河南科技大學(xué), 2011.

[12]楊世銘,陶文銓.傳熱學(xué)[M].3版.北京:高等教育出版社, 1998.

[13]李奇亮. 基于界面測(cè)溫的面接觸摩擦溫度場(chǎng)研究[D].合肥:合肥工業(yè)大學(xué),2009.

[14]Majcherczak D,Dufrenoy P,Berthier Y.Tribological Thermal and Mechanical Coupling Aspects of the Dry Sliding Contact[J].Tribology International,2007,40:834-843.

(編輯王艷麗)

Calculation Method of Spherical Plain Bearing Friction Temperature Field

Yu JianweiWang LifeiWei Wei

Hefei University of Technology,Hefei,230009

Taking thrust spherical plain bearing as the research object,the contact stress under axial load was analyzed and based on the experimental friction torque the dynamic friction coefficient was solved.Then combined with contact stress and friction coefficient,heat flow rate was calculated and served as the second boundary condition for the unsteady temperature field finite element analysis of the bearing.Finally,using infrared camera to shoot the outer surface of two-dimensional temperature distribution of the bearing friction process,then reliability of simulation results was tested.Simulation and experimental results show that: the calculation method has solved the problems of accurate calculation of spherical friction coefficient and spatial distribution of the heat flow rate.Ignoring the spherical gap when modeling will affect the contact stress distribution directly and cause a larger error in temperature distribution.Considering a spherical gap a more accurate bearing friction temperature field can be obtained.

spherical plain bearing;spherical gap;contact stress;temperature field;infrared measurement;finite element

2014-06-13

國家自然科學(xué)基金資助項(xiàng)目(51075114)

TH117.1DOI：10.3969/j.issn.1004-132X.2015.10.002

俞建衛(wèi),男,1956年生。合肥工業(yè)大學(xué)摩擦學(xué)研究所及安徽省粉末冶金工程技術(shù)研究中心研究員。主要研究方向?yàn)槟Σ翆W(xué)測(cè)試技術(shù)。王禮飛,男,1991年生。合肥工業(yè)大學(xué)機(jī)械與汽車工程學(xué)院碩士研究生。魏巍，男，1986年生。合肥工業(yè)大學(xué)儀器科學(xué)與光電工程學(xué)院博士研究生。