陳瑛琳，史文華，陳江義，秦東晨，袁峰

(鄭州大學(xué) 機械與動力工程學(xué)院，鄭州 450001)

為滿足RV減速器外形尺寸小，承受扭矩大的工作需求，其主軸承往往采用薄壁、多球結(jié)構(gòu)[1]。球數(shù)增多會導(dǎo)致軸承受力和變形計算量大幅增加，求解速度慢。如何快速計算和分析RV減速器主軸承的受力和變形是一個亟待解決的問題。

常用的軸承設(shè)計方法為數(shù)值分析法[2]，其中的擬靜力學(xué)方法適用于高速、重載場合，計算速度快，是目前應(yīng)用廣泛的方法之一[3-6]，該方法以赫茲理論[7]為基礎(chǔ)，考慮離心力、陀螺力矩作用，采用溝道控制假設(shè)，可分析球和溝道在點、線接觸下的載荷分布和接觸變形[8]。文獻[9-10]的實踐證明該方法在內(nèi)、外圈接觸角差異較大時具有良好的數(shù)值穩(wěn)定性和較高的計算精度。

擬靜力學(xué)法中求解非線性方程組多采用牛頓-拉弗森算法，求解簡單但對初值敏感，且方程組的數(shù)量隨球數(shù)增加而成倍增加，求解速度也相應(yīng)變慢。針對該問題，文獻[11]在算法中引入迭代步長松弛因子，并采用二分法對松弛因子的取值進行優(yōu)化，同時通過矩陣變換簡化了雅克比矩陣求解過程，從而縮減了迭代時間；文獻[12]對迭代公式進行改進和修正，降低了因初值偏差給收斂性帶來的影響;文獻[13]建立迭代約束條件，使迭代能在合理范圍內(nèi)進行;文獻[14]提出減少非線性方程數(shù)量和迭代變量約束的方法，提高了方程組求解的效率。然而，對于RV減速器主軸承，未知數(shù)和方程組數(shù)量較多，即使采用優(yōu)化算法，也難以快速得出結(jié)果[15]。文獻[16]提出可將滾子軸承沿受載方向視為對稱，非承載區(qū)滾子受力情況視為相同，減少了方程數(shù)量，使未知參數(shù)減少近一半，從而快速計算；但在大扭矩作用下，RV減速器軸承載荷分布復(fù)雜，在三維空間內(nèi)并不沿載荷方向?qū)ΨQ分布，該方法理論上不再適用。

鑒于上述原因，在建立RV減速器主軸承(角接觸球軸承)擬靜力學(xué)模型的基礎(chǔ)上，針對球數(shù)較多導(dǎo)致非線性方程組數(shù)量成倍增加，求解難度增大的問題，提出一種復(fù)雜工況下的簡化分析方法，并從計算速度和分析精度兩方面與原方法進行對比。

1 RV減速器主軸承擬靜力學(xué)模型

為便于表述，省略部分球，角接觸球軸承球的角位置關(guān)系如圖1所示，以軸承中心為原點，軸向中心線為z軸，建立坐標系Oxyz。球之間的夾角Δψ=2π/Z(Z為球數(shù))，設(shè)定第1個球的位置角為γ，則第j個球的位置角ψj=2π(j-1)/Z+γ。

圖1 球角位置示意圖Fig.1 Diagram of ball angle position

在力Fx，F(xiàn)y，F(xiàn)z和力矩Mx，My作用下，軸承發(fā)生非線性偏移，原有幾何接觸關(guān)系也發(fā)生改變。假定軸承外圈固定，則第j個球的中心和內(nèi)外圈溝曲率中心的位置關(guān)系如圖2所示。若內(nèi)圈溝道相對外圈溝道分別沿x，y，z軸方向產(chǎn)生偏移δx，δy，δz，繞x，y軸中心線方向產(chǎn)生角度偏轉(zhuǎn)θx，θy，則在位置角ψj處內(nèi)外圈溝曲率中心的軸向距離A1j和徑向距離A2j可表示為

圖2 載荷作用前后第j個球中心和內(nèi)外圈溝曲率中心的位置關(guān)系Fig.2 Position relationship between j-th ball center and curvature center of inner and outer groove before and after loading

(1)

Gi=Dpw/2+(fiDw-Dw/2)cosα0,

式中：A為受載前軸承內(nèi)外圈溝道中心的距離；α0為軸承初始接觸角；Gi為內(nèi)溝道轉(zhuǎn)動半徑；Dpw為球組節(jié)圓直徑；Dw為球直徑；fi為內(nèi)溝曲率半徑系數(shù)。

為簡化計算，引入變量X1j，X2j，分別表示受載后球中心到外溝曲率中心之間的軸向和徑向距離。根據(jù)勾股定理可得

(A1j-X1j)2+(A2j-X2j)2-[(fi-0.5)Dw+δij]2=0,

(2)

式中：fe為外溝曲率半徑系數(shù)；δij，δej分別為受載后內(nèi)、外圈溝道的法向變形。

用A1j，A2j，X1j，X2j分別表示受載后軸承內(nèi)、外圈的實際接觸角αij，αej,即

(3)

為減少計算量，提高迭代速度，聯(lián)立(2)，(3)式，將δij，δej，X1j，X2j用αij和αej表示，可使未知數(shù)減半，其轉(zhuǎn)換公式為

(4)

對垂直于軸承軸線并通過第j個球中心的平面進行分析，球受到內(nèi)外圈溝道接觸載荷Qij和Qej、離心力Fcj、陀螺力矩Mgj和摩擦力作用(圖3)。若非共面的摩擦力很小，可以假設(shè)陀螺力矩Mgj完全被球與外圈溝道接觸區(qū)的摩擦力所阻止，取陀螺力矩作用控制系數(shù)λij=0，λej=2，若條件不滿足，則取λij=1，λej=1。

圖3 第j個球所受載荷Fig.3 Load on j-th ball

由赫茲接觸理論可知任意位置角ψj處球的法向載荷與法向變形的關(guān)系為

(5)

式中：Qij，Qej分別為內(nèi)、外圈與球的法向接觸載荷；Kij，Kej分別為內(nèi)、外圈溝道位移-載荷系數(shù)。

單個球所受離心力Fcj為

(6)

式中：m為單個球質(zhì)量；ωmj為球公轉(zhuǎn)角速度。

陀螺力矩Mgj可表示為

Mgj=JωmjωRjsinβj，

(7)

式中：J為球轉(zhuǎn)動慣量；ωRj為球自轉(zhuǎn)角速度；βj為球姿態(tài)角。

由圖3球受載情況可知單個球在豎直方向和水平方向的受力平衡，平衡方程為

(8)

根據(jù)軸承整體所受力和力矩的平衡關(guān)系，可得軸承平衡方程組為

。(9)

2 簡化分析方法

傳統(tǒng)擬靜力學(xué)方法通過(8)，(9)式迭代求解，在此基礎(chǔ)上，針對RV主軸承的多球特性，設(shè)計了1套RV減速器主軸承擬靜力學(xué)模型簡化分析方法，計算流程(圖4)為：

圖4 RV減速器主軸承擬靜力學(xué)分析流程Fig.4 Analysis process of quasi-static model of main bearing for RV reducer

1)假定軸承處于靜止狀態(tài)，采用靜力學(xué)方法求解軸承在靜載荷下的內(nèi)外圈溝道接觸角、位移和變形等參數(shù)。以靜力學(xué)計算結(jié)果為基礎(chǔ)，求解軸承運轉(zhuǎn)時各零件的速度、陀螺力矩、離心力等參數(shù)。

2)假定軸承位移參數(shù)δx，δy，δz，θx，θy不變，且第1個球的位置角γ為0，以靜力學(xué)算法求解的位移和角度作為初值，對 (8)式進行迭代，求得軸承內(nèi)、外圈動態(tài)接觸角αij和αej。

3)將接觸角代入(9)式進行迭代，求得新的軸承位移、法向變形以及所受載荷等參數(shù)。判斷新位移參數(shù)δx，δy，δz，θx，θy是否達到平衡，若是則輸出結(jié)果，若否則更改第1個球的位置角γ，重新計算參數(shù)，不斷重復(fù)直至位移參數(shù)滿足精度要求。

該流程中迭代算法以牛頓-拉弗森算法為基礎(chǔ)，增加迭代步長松弛因子，保證每次迭代都使方程的絕對值減小，避免迭代跳出范圍，造成不收斂。同時，對角度迭代中的變量進行判斷和修正，避免出現(xiàn)死循環(huán)。

該算法中方程組數(shù)量為2Z+5，RV主軸承球數(shù)較多，則方程組求解易出現(xiàn)迭代緩慢，難以達到精度等問題，有必要對其進行簡化。因為球排列密集時，相鄰球受力相似，接觸角相近，假設(shè)相鄰k個球接觸狀態(tài)一致，在計算中將其作為完全相同的零件進行處理，進而成倍減少方程組數(shù)量。

若原球數(shù)Z可被k整除，則簡化后的球數(shù)N=Z/k，第j個球的位置角ψ′j=2π(j-1)k/Z+γ，則(9)式中x方向受力平衡方程可轉(zhuǎn)換為

(10)

若Z不能被整除，且余數(shù)為h，則球數(shù)N=[Z/k]+1，(10)式可轉(zhuǎn)換為

(11)

同理可對其他幾個力和力矩平衡方程式進行簡化，替代(9)式進行計算。

3 實例分析

以某減速器用主軸承H76/182RV角接觸球軸承為例分析，其主要結(jié)構(gòu)參數(shù)見表1。

表1 H76/182RV角接觸球軸承主要結(jié)構(gòu)參數(shù)Tab.1 Main structural parameters of H76/182RV angular contact ball bearing

軸承轉(zhuǎn)速n取40 r/min，施加徑向載荷Fy=4 000 N，軸向載荷Fz=13 720 N，力矩載荷My=2 450 N·m，按k個球進行分組，分別取k=1，2，3，4(k=1表示沒有簡化的情形)。采用上述簡化模型計算軸承第1個球與內(nèi)圈的法向載荷以及接觸角，其結(jié)果如圖5和圖6所示。由圖5可知：當k為2，3，4時，球所受內(nèi)圈的法向載荷與未簡化時的法向載荷有相同的分布趨勢；k取2時，計算結(jié)果與未簡化模型結(jié)果幾乎一致，k值增大，求解結(jié)果與未簡化模型結(jié)果差異增大，但誤差仍能滿足基本的工程需求。由圖6可知:簡化前后內(nèi)圈接觸角曲線分布幾乎相同。由此可知簡化模型是合理的。

圖5 簡化與原算法計算所得第1個球與內(nèi)圈的法向載荷Fig.5 Normal load of first ball and inner ring calculated by simplified and original method

為進一步檢驗算法是否在不同工況下都適用，根據(jù)RV減速器所能承受的最大載荷和轉(zhuǎn)速，分別取軸承轉(zhuǎn)速n為20～200 r/min，軸向載荷Fz為1 372mN，徑向載荷Fy為500mN，力矩載荷My為245mN·m(m為工況載荷系數(shù)，m=1，2，…，10)。取k為1，2，3，4，記錄不同工況下程序?qū)嶋H計算時間，結(jié)果見表2。并以內(nèi)外圈溝道在x方向的相對位移δx為參數(shù)，計算在不同工況下簡化算法與原算法的相對誤差，結(jié)果如圖7和圖8所示。

表2 不同工況下簡化算法的計算時間Tab.2 Calculation time of simplified method under different operating conditionss

圖7 轉(zhuǎn)速為40 r/min時不同載荷下內(nèi)外圈溝道相對位移δx的相對誤差Fig.7 Relative error of relative displacement δx of inner and outer raceway under different loads at rotational speed of 40 r/min

圖8 Fy=5 000 N，F(xiàn)z=13 720 N，My=2 450 N·m時不同轉(zhuǎn)速下內(nèi)外圈溝道相對位移δx的相對誤差Fig.8 Relative error of relative displacement δx of inner and outer raceway under load of Fy=5 000 N,Fz=13 720 N,My=2 450 N·m and different rotational speeds

由表2可知：采用簡化分析方法能夠有效地提升計算效率，減少計算時間。由圖7、圖8可知，不同轉(zhuǎn)速和載荷工況下簡化算法與原算法的相對誤差都在3%以內(nèi)。故可認為載荷和轉(zhuǎn)速對于簡化分析方法的求解精度影響不大，簡化算法在任意載荷和轉(zhuǎn)速工況下都具有適用性。

4 結(jié)束語

建立了RV減速器主軸承擬靜力學(xué)分析模型，并針對其多球特性導(dǎo)致的數(shù)值分析計算量大，求解緩慢的問題，提出了一種簡化分析方法，該方法假設(shè)相鄰k個球所受的接觸載荷及接觸角相同，可以簡化為1個零件進行計算。并以某減速器主軸承H76/182RV角接觸球軸承為例，對該方法進行驗證，結(jié)果表明簡化方法能有效提高計算效率。當k=2，3，4時，簡化分析方法誤差都小于3%，可以滿足工程需要。該簡化方法的計算精度受載荷和轉(zhuǎn)速影響較小，誤差在較小范圍內(nèi)波動，可以認為該簡化分析方法適用于不同的載荷和轉(zhuǎn)速工況。提出的方法不僅可以用于RV減速器主軸承的擬靜力學(xué)分析，還適用于球數(shù)較多的其他角接觸球軸承。