郭向東,王燕霜
(1.洛陽軸研科技股份有限公司 軌道交通軸承事業部,河南 洛陽 471039;2.河南科技大學機電工程學院,河南 洛陽 471003)
現代航空發動機主軸軸系多采用雙轉子結構,即將低壓轉子和高壓轉子形成機械耦合,而連接低壓轉子和高壓轉子的軸承為中介軸承[1-2]。當今航空發動機追求大推重比,從而要求主軸軸承在更高的dm·N值下工作[3]。隨著dm·N值的提高,將帶來中介軸承保持架和滾子的打滑、接觸表面損傷和保持架不穩定等一系列問題[4-5],造成軸承壽命的降低。可見,軸承動態特性直接決定了中介軸承的壽命。所以,對中介軸承進行動態特性分析具有重要意義。
滾動軸承的動態特性研究經歷了不斷發展的過程。Jones[6]首先提出了滾道控制理論,建立了滾動軸承擬靜力學分析模型,但該模型沒有考慮彈流潤滑的作用和保持架對軸承力矩的影響。Gupta求出了軸承運動和力的瞬態解,并且編制了大型應用程序ADORE[7]。Svenska Kullager-Fabriken公司開發的大型計算機分析軟件“SHARBERTH”,采用擬動力學對滾子和保持架受力問題進行了分析和討論。但該軟件分析時考慮的因素比較多,程序過于龐大,求解時間較長,在工程應用中有一定的局限性。Aramaki[8]開發了滾動軸承動力學分析軟件BRAIN,是目前能分析軸承類型最多的軟件,但對潤滑和摩擦阻力的研究仍有待完善。Stacke等[9]開發了BEAST三維分析軟件,但對于部分膜潤滑和擠壓阻尼的計算還有待完善。
雖然動力學模型分析方法考慮的因素比較全面,但由于高速下滾動軸承各元件間動態特性非常復雜,考慮的因素較多,必然帶來數學上求解的困難,因此滾動軸承動力學理論迄今為止仍不完善[10]。鑒于此,采用擬靜力學分析的方法,對中介軸承動態特性進行研究。
因中介軸承工作時內、外圈同時同向旋轉,其受力和運動情況比外圈固定或內圈固定的軸承相對復雜,故求解中介軸承載荷分布時假設:
(1)滾子與內、外滾道之間均形成彈流潤滑油膜,不存在乏油問題。
(2)滾子與滾道的接觸中,只有滾子發生變形,忽略內、外滾道的變形。
(3)滾子與內、外圈的接觸載荷對稱于徑向載荷作用線均勻分布。
建立中介軸承在同向旋轉時的擬靜力學模型,分析軸承內部載荷的分布特征和運動特性。
簡化的中介軸承滾子運動情況如圖1所示。圖中,Di為內滾道直徑;De為外滾道直徑;Dr為滾子直徑;ωi為內圈角速度;ωe為外圈角速度;ωc為保持架角速度;ωrj為第j個滾子的角速度;ωruj為非承載區第j個滾子的角速度(未示出)。

圖1 中介軸承運動分析
滾子與內、外滾道之間的平均速度為
(1)
(2)

非承載區滾子與外圈滾道接觸表面平均速度為
(3)
滾子與內、外滾道之間的相對滑動速度為
(4)
(5)
圖2給出了滾子在承載區和非承載區的受力情況。承載區滾子受力如圖2a所示,力和力矩的平衡方程組為

圖2 滾子受力分析
(6)
式中:Qij,Qej分別為各個滾子與內、外圈滾道之間的接觸載荷;Pij,Pej,Peu為作用于滾子的切向動壓力;Tij,Tej,Teu為作用于滾子上的切向摩擦力;Fm為單個滾子的離心力;Foa,Fou為滾子受油氣混合物的繞流阻力;Fcj,Fcu為承載區和非承載區滾子與保持架之間的作用力;f為摩擦因數。承載區內當保持架推動滾子運動時取正號,反之取負號。
將(6)簡化合并得
(f-1)Tej+f(Pej-Pij)+(f+1)Tij-fFoa=0,j=0,1,2,…,N,
(7)
軸承正常工作時,保持架以勻角速度ωc平穩運轉,因此保持架處于平衡狀態, 其力矩平衡方程為
(8)
式中:Fc0為第0滾子與保持之間的作用力,z為總滾子數,2N+1為承載區內滾子數。
非承載區滾子受力分布如圖2b所示,其力平衡方程為
Fou-Fcu-Peu-Teu=0。
(9)
(6)~(9)式中,各參數意義及求解方程如下。
2.2.1 滾子接觸載荷Qij和Qej
中介軸承的動態載荷分布如圖3所示,根據平衡條件可得

圖3 中介軸承動態載荷分布
(10)
式中:Qi0為滾子與內圈之間的接觸載荷。
內圈相對于外圈的位移量為
(11)
式中:Δ0為徑向載荷作用下軸承內、外圈的相對位移;δj為第j個滾子與內、外套圈之間的接觸變形;ur為軸承的徑向游隙。
其中,第j個滾子總彈性變形量為
δj=δij+δej。
(12)
由文獻[11]可得
(13)
其中,Kn為載荷變形系數,可參考文獻[12]計算。
由(11)~(13)式可得承載區的變形協調方程
1.13(hij+hej)]=0,j=1,2,…,N。
(14)
聯立(10)和(14)式,可求解N+1個未知數Qi0,Qi1,…,QiN。承載區滾子數為M=2N+1,考慮對稱分布,為簡化計算,故只求解N+1個未知數。
2.2.2 切向動壓力Pij,Pej和Peu
軸承工作時,內、外圈之間的潤滑油受到擠壓對滾子產生切向壓力[11],承載區切向動壓力為
(15)
(16)
非承載區切向動壓力為
(17)
其中,
(18)
(19)
(20)
2.2.3 切向摩擦力Tij,Tej和Teu
由文獻[11]可知,內、外圈滾道通過油膜作用于滾子的切向摩擦力為
(21)
(22)
非承載區作用于滾子的切向摩擦力
(23)
其中,
(24)
(25)
(26)
無量綱積分
(27)
(28)
其中,
(29)
(30)
G=αE0;
(31)

2.2.4 單個滾子離心力Fm
由文獻[13]可知,單個滾子離心力計算式為
(32)
式中:mb為滾子質量。
2.2.5 油氣混合物對滾子的繞流阻力Foa和Fou
由文獻[13]可知,滾子受油氣混合物繞流阻力的計算式為
(33)
式中:Cd為繞流阻力系數;ρ為油氣混合物的密度。
2.2.6 滾子與保持架之間的作用力Fcj和Fcu
承載區和非承載區滾子與保持架之間的作用力,由(7)式和(8)式組成方程組,采用Newton-Raphson法進行求解,具體過程略。
以某型號中介軸承為算例,其主要參數見表1。

表1 中介軸承主要參數
圖4為當內圈轉速為5 000 r/min,外圈轉速為10 000 r/min時,內圈與滾子之間的接觸載荷隨徑向力變化圖。可以看出隨著徑向力的增加,內圈與滾子之間的接觸載荷隨之增加。并且在1 000 N增至2 000 N時,承載滾子數增加。

圖4 內圈與滾子之間的接觸載荷分布
圖5所示為低壓轉子轉速為5 000 r/min,高壓轉子轉速為10 000 r/min時,中介軸承滾子自轉角速度與徑向力的變化情況。其中,0號滾子為受力最大的滾子,1,2,3號滾子位于0號滾子一側,與0號滾子的距離等量遞增。從圖中可以看出,隨著徑向力的增加,承載區內各滾子的自轉角速度隨之減小,且滾子越偏離徑向力作用線,自轉角速度越大。 隨著徑向力的增加,受載區滾子與內圈間的法向接觸載荷隨之增大。滾子與外圈間的法向接觸載荷也隨之增大。但是由于離心力的作用,滾子與內圈接觸載荷增加量小于滾子與外圈之間接觸載荷的增加量,這使得內圈對滾子的牽引力增加量小于外圈對滾子的牽引力的增加量,從而使內外圈滾道作用在滾子上,使滾子自轉的摩擦力矩減小,造成滾子自轉角速度隨之減小。在同一徑向力的作用下,越接近徑向力作用線的滾子,內、外圈滾道作用在其上的摩擦力矩越小,滾子自轉角速度越小。在1 700~2 000 N時,因承載滾子數增加,內圈對滾子接觸力和牽引力的減小量小于外圈對滾子接觸力和牽引力的減小量,使得內、外圈作用在滾子上的摩擦力矩增加,致使滾子自轉角速度隨之增加。

圖5 滾子自轉角速度隨徑向力的變化情況
圖6和圖7所示為徑向載荷為2 000 N、內圈轉速為5 000 r/min、外圈轉速為13 000 r/min時,中介軸承滾子自轉角速度分別隨內、外圈轉速變化情況。從圖中可以看出:當中介軸承外圈轉速增加時,各滾子自轉角速度隨之增加,而當內圈轉速增加時,各滾子的自轉角速度隨之減小。這是因為當內、外圈轉速增加時,其作用于滾子的牽引力減小,當外圈和滾子間牽引力減小時,內、外圈作用在滾子上的摩擦力矩增加,使得滾子自轉速度增加;內圈和滾子間牽引力減小時,內、外圈作用在滾子上的摩擦力矩減小,致使滾子自轉速度減小。

圖6 滾子自轉角速度隨外圈轉速的變化情況

圖7 滾子自轉角速度隨內圈轉速的變化情況
表2為內圈轉速為5 000 r/min,外圈轉速為10 000 r/min,套圈與滾子間油膜厚度與載荷的關系。從表2中可以看出,在承載區內滾子不變時,隨著徑向力的增加,油膜厚度減小,并且內圈滾道與滾子之間的油膜厚度小于外圈滾道與滾子之間的油膜厚度。

表2 軸承油膜厚度與載荷的關系
通過雙轉子試驗機可以測得不同徑向載荷下中介軸承保持架的轉速。算例中的計算結果與試驗結果比較,最大誤差為4.28%,可見文中建立的中介軸承動態特性分析模型是正確的。
(1)建立的中介軸承動態特性分析模型誤差較小,可為中介軸承的設計提供理論指導。
(2)對于內、外圈同向旋轉的中介軸承,隨著徑向力的增加,軸承承載區內各滾子的自轉角速度隨之減小。中介軸承承載區內各滾子越偏離徑向力作用線,滾子自轉角速度越大。
(3)對于內、外圈同向旋轉的中介軸承,外圈轉速增加時,各滾子自轉角速度隨之增加,而當內圈轉速增加時,各滾子的自轉角速度隨之減小。
(4)徑向載荷的大小對彈流潤滑油膜的形成影響較大,適當地減小徑向載荷有利于形成完全彈流潤滑油膜,保證高速滾子軸承工作時能夠達到良好的潤滑狀態。