嚙合錯位對高速齒輪傳動振動響應的影響

李良祥， 魏靜， 張愛強， 倪德， 鐘團結(jié)

(1.重慶大學 機械傳動國家重點實驗室， 重慶 400044； 2.中國航發(fā)湖南動力機械研究所，湖南 株洲 412000)

航空齒輪等高速齒輪傳動系統(tǒng)中，由于軸、軸承、殼體和齒輪等受載變形以及零部件制造和裝配誤差，易造成齒輪嚙合錯位即齒輪不對中，根據(jù)國內(nèi)外相關文獻統(tǒng)計，由于齒輪錯位即不對中所引起或與之相關的故障占所有旋轉(zhuǎn)機械故障的60%以上[1]，微小的錯位可能造成系統(tǒng)轉(zhuǎn)速失穩(wěn)、振動加劇、齒輪偏載和疲勞破壞[2-5]等嚴重后果。關于齒輪嚙合錯位，國內(nèi)外學者進行了較多的研究，就研究方法而言，主要分為解析法[6-12]和仿真[13-16]。Shao等[6]研究了斜齒輪嚙合錯位對齒輪嚙合力的影響。馬輝等[7]研究了修形對嚙合錯位的影響，結(jié)果表明，鼓向修形可改善錯位對齒輪接觸應力的影響。Lin等[8]、Wang等[9]和Saxena 等[10]建立了斜齒輪副嚙合剛度計算模型，說明了嚙合錯位對斜齒輪嚙合剛度的影響，但他們的模型均為切齒嚙合模型，模型較為簡單，考慮因素較少。Guan等[11]和Wei等[12]以齒式聯(lián)軸器為對象，考慮錯位建立齒式聯(lián)軸器數(shù)學模型，研究了錯位對輪齒接觸特性的影響，研究結(jié)果表明，在錯位情況下輪齒表面由原來的線接觸變?yōu)辄c接觸。Shi等[13]建立了準雙曲面齒輪轉(zhuǎn)子系統(tǒng)數(shù)學模型，研究軸不對中對準雙曲面齒輪動態(tài)響應的影響，結(jié)果表明，軸系不對中對速度、轉(zhuǎn)矩波動等響應產(chǎn)生了較大影響。Zhan[14]建立了嚙合剛度的有限元模型，分析了齒輪錯位對嚙合剛度的影響。Jones等[15]和Li[16]運用有限元方法，建立齒輪模型研究了嚙合錯位對齒輪嚙合性能的影響，錯位改變了嚙合力、傳動誤差和輪齒接觸模式，輪齒修形可減小錯位的影響，但他們的建模過程中均未考慮系統(tǒng)各齒輪軸受載產(chǎn)生的彎扭耦合作用對嚙合錯位的影響。

國內(nèi)外學者在研究齒輪嚙合錯位時，所建立的模型多為切齒錯位模型，模型中未考慮齒輪軸在運轉(zhuǎn)過程中受載產(chǎn)生的彎扭耦合作用(各軸系柔性)對嚙合錯位的影響；研究多局限于錯位對齒輪副嚙合剛度、接觸狀態(tài)等齒輪自身嚙合性能影響，而錯位對系統(tǒng)振動特性，特別是錯位對應力造成的影響研究較少。針對以上研究不足，本文運用超單元縮聚方法對系統(tǒng)各軸系和齒輪輪齒進行縮聚創(chuàng)建縮聚主節(jié)點，實現(xiàn)軸系及齒輪柔性化，結(jié)合多體動力建模理論，在各主節(jié)點施加約束和載荷，建立了考慮錯位的全柔體單級航空高速直齒輪傳動系統(tǒng)動力學模型，研究了齒輪嚙合錯位對系統(tǒng)影響。

1 嚙合錯位齒輪傳動系統(tǒng)動力學模型

嚙合錯位主要分為3種：角度錯位、軸向錯位和徑向錯位，如圖1所示，錯位導致齒輪非正常嚙合，進而影響系統(tǒng)動態(tài)特性，文獻[12]的錯位試驗研究結(jié)果表明，錯位增大了振動位移頻譜中各主頻處幅值，其中角度錯位的增幅更大，因此本文主要研究角度錯位。

圖1 嚙合錯位類型Fig.1 Types of meshing misalignment

耦合動力學模型基于Guyan矩陣縮聚理論完成，通過子結(jié)構(gòu)分析法將結(jié)構(gòu)單元凝聚為一個超單元，大大縮減模型自由度。選取某一部件其動力平衡方程為：

(1)

式中：M、K、C分別為質(zhì)量、剛度和阻尼矩陣；F為載荷矩陣；u為位移矩陣。

u分為主自由度um和從自由度us，兩者關系為：

us=Tum

(2)

式中：T代表um和us相互關系，um和us分別用nm和ns階向量表示，T為nm×ns階矩陣，式(2)表示為：

(3)

靜態(tài)時：

Ku=F

(4)

假設部件有m個主自由度，且有s個從自由度，則式(4)可以表示為：

(5)

式中：Kmm、Kms、Ksm、Kss依次為m×m、m×s、s×m、s×s階剛度矩陣；主自由度上的載荷向量為Fmm。

由式(5)可得：

uss=-K-1Ksmumm

(6)

[Kmm-KmsKss-1Ksm]umm=Fmm

(7)

由式(2)、(6)可得：

(8)

聯(lián)立式(6)、(7)、(8)可得

(9)

通過有限元縮聚完成各子結(jié)構(gòu)分析，同時將傳遞矩陣T寫入式(1)的阻尼矩陣及質(zhì)量矩陣中：

Mm=(T*)TM(T*)，Cm=(T*)TC(T*)

(10)

主自由度代替了原本結(jié)構(gòu)的自由度，其縮減的自由度動力學方程為：

(11)

根據(jù)上述建模理論，建立傳動系統(tǒng)動力學模型，流程如圖2所示，其中輸入、輸出齒輪參數(shù)如表1所示，齒輪材料為9310合金鋼。

圖2 節(jié)點有限元法動力學建模流程Fig.2 Dynamics modeling process by node finite element method

表1 齒輪參數(shù)Table 1 Gear parameters

采用十節(jié)點二階四面體單元(Solid187單元)對兩齒輪進行網(wǎng)格劃分，劃分過程中對幾何進行簡化(忽略倒角等微小幾何特征)，保證模型精度兼顧計算效率。高速輸入軸單元數(shù)量80 673，節(jié)點數(shù)量134 110，輸出軸單元數(shù)量109 170，節(jié)點數(shù)量181 379。

在輸入、輸出軸上設置MPC(柔性多點約束)主節(jié)點用以創(chuàng)建多體耦合連接關系。其中輸入軸設置3個主節(jié)點：1個功率輸入節(jié)點與輸入扭矩建立耦合關系，2個軸承節(jié)點與箱體建立耦合關系；輸出軸設置3個主節(jié)點：1個功率輸出節(jié)點與輸出轉(zhuǎn)速建立耦合關系，2個軸承節(jié)點與箱體建立耦合關系。

運用軸承節(jié)點建立齒輪軸系與箱體之間的耦合關系，軸承滾動體與內(nèi)外圈之間的耦合關系為：

(12)

式中：

Fbody,1=[FX1FY1FZ1MX1MY1MZ1]T

(13)

Fbody,2=[FX2FY2FZ2MX2MY2MZ2]T

(14)

qbody,1=[x1y1z1θ1ρy1ρz1]T

(15)

qbody,2=[x2y2z2θ2ρy2ρz2]T

(16)

(17)

運用嚙合節(jié)點建立主動和從動齒輪嚙合關系，如圖3所示，在每個輪齒節(jié)圓位置沿齒寬方向創(chuàng)建均勻分布的縮聚主節(jié)點，實現(xiàn)輪齒柔性化建模。

圖3 考慮嚙合錯位的齒輪嚙合單元Fig.3 Gear meshing unit considering misalignment

將主節(jié)點載荷分布到齒面節(jié)點上，可計及輪齒彎曲、剪切等變形影響。另外，該方法可考慮軸系、輪齒瞬時變形、偏斜等導致的嚙合錯位狀態(tài)，進而影響系統(tǒng)動態(tài)響應。

2 嚙合錯位對系統(tǒng)動態(tài)特性影響

2.1 齒輪嚙合特性

2.1.1 齒輪副嚙合力

多體動力學模型中求解齒輪副所受圓周力和徑向力有效值隨錯位角度變化趨勢如圖4、5所示，由圖4、5可知，總體上圓周力和徑向力均隨錯位角度增大呈增大趨勢，0.05°范圍內(nèi)兩者變化趨勢較為緩和，0.1°時圓周力和徑向力大小分別為7 852 N和3 418 N相比于正常情況下的增幅分別為2.9%和21.8%，錯位角超過0.1°后圓周力和徑向力增幅顯著。

圖4 齒輪副圓周力Fig.4 Circumferential force of gear pair

圖5 齒輪副徑向力Fig.5 Radial force of gear pair

2.1.2 動態(tài)接觸應力

考慮齒輪的接觸強度。得到齒輪動態(tài)接觸應力可以分析其受力狀態(tài)、嚙合情況和沖擊程度等，對研究齒輪嚙合特性至關重要。

嚙合錯位角θ從0°變化至0.2°時，齒輪副動態(tài)接觸應力隨θ變化的時頻域響應如圖6所示，齒輪材料為9310鋼，屈服強度940 MPa，正常情況下接觸應力均值是595 MPa，滿足接觸強度要求。隨著θ增大接觸應力明顯增大，這是因為錯位后輪齒接觸面偏向一側(cè)，造成偏載。相同時間內(nèi)應力幅交變頻率更快，加速了材料的疲勞破壞；接觸應力頻譜中主要包含嚙頻及其2倍頻、3倍頻，嚙頻處幅值隨著θ增大明顯增加。

圖6 接觸應力隨錯位角變化時頻域響應Fig.6 Time frequency response of contact stress varying with misalignment angle

接觸應力均值和幅值隨錯位角θ變化如圖7所示，隨θ增大接觸應力均值和幅值呈增大趨勢，θ從0.05°至0.2°接觸應力均值分別為632、762、983、1 247 MPa，相比于正常情況下分別增加6.22%、28.1%、65.2%、109.6%；θ從0.05°至0.2°接觸應力幅值分別為620、827、1 235、1 446 MPa，相比于正常情況下分別增加3.85%、38.5%、106%、142%；總的來看，角度錯位對接觸應力影響較大，θ在0.05°內(nèi)接觸應力增幅較小，超過0.1°后急劇增加，錯位0.15°時接觸應力已超過材料的屈服強度。

圖7 接觸應力隨錯位角變化Fig.7 Variation of contact stress with misalignment angle

2.1.3 傳動誤差

傳動誤差峰峰值隨錯位角θ變化如圖8所示，傳動誤差峰峰值隨θ增大呈明顯增大趨勢，θ從0.05°至0.2°峰峰值分別為4.1、5.4、5.8、7.0 μm，相比于正常情況下分別增加5.1%、38.5%、48.7%、79.5%，錯位在0.05°內(nèi)，傳動誤差峰峰值變化較小，超過0.05°后增幅加劇。

圖8 傳動誤差峰峰值隨錯位角變化Fig.8 The peak-peak value of transmission error varies with misalignment angle

2.2 系統(tǒng)動態(tài)特性

2.2.1 振動加速度

錯位角θ從0°增至0.2°，拾振點振動加速度隨θ變化的時頻域響應如圖9所示，振動加速度隨θ增大急劇增大；加速度頻譜圖主要包含嚙頻及其2倍頻、3倍頻，嚙頻及其3倍頻處幅值隨θ增大明顯增大，隨著錯位程度增大最大幅值對應頻率向高頻轉(zhuǎn)變，θ為0.2°時最大幅值頻率由原先嚙頻的2倍頻轉(zhuǎn)變?yōu)閲ьl的3倍頻，與此同時由頻域圖可以明顯看出，隨著θ增大，各主頻處出現(xiàn)明顯的邊頻帶。

圖9 振動加速隨錯位角變化時頻域響應Fig.9 Time frequency response of vibration acceleration varies with misalignment angle

振動加速度RMS隨θ變化如圖10所示，總體來看，加速度隨θ增大呈增大趨勢，錯位小于0.05°時振動加速度變化較緩，超過0.05°時振動加劇，θ從0.05°至0.2°振動加速度分別為1 371.6、2 192.5、5 160.3、8 996.1 m·s-2，相比于正常情況下分別增加33.9%、114.1%、403.9%、778.5%，由此可見，振動加速度對嚙合錯位較為敏感，微小的錯位即可造成系統(tǒng)劇烈振動。

圖10 振動加速度隨錯位角變化Fig.10 Variation of vibration acceleration with misalignment angle

2.2.2 動應力

多柔體動力學模型結(jié)合有限元方法獲取系統(tǒng)動應力(模態(tài)應力恢復法)，該方法簡單高效適用于絕大多數(shù)的機械構(gòu)件。不同錯位角下，高速輸入齒輪軸動應力云圖如圖11所示，未錯位情況下，最大動應力出現(xiàn)在齒輪和軸連接處(應力集中)，且軸系動應力隨著錯位角增大明顯增大。

圖11 不同錯位角動應力云圖Fig.11 Dynamic stress nephogram of different misalignment angles

測點動應力時頻域響應如圖12所示，動應力隨著錯位角θ增大而增大，θ為0°和0.2°時應力峰值分別為65 MPa和165 MPa，增大近1.5倍；頻譜圖中，軸頻處幅值最大，嚙頻處幅值較小但隨θ增大其幅值逐漸增大，同時嚙頻處出現(xiàn)大量邊頻帶，θ為0.2°嚙頻處邊頻帶如圖13所示。

圖12 動應力隨錯位角變化時頻域響應Fig.12 Time frequency response of dynamic stress varying with misalignment angle

圖13 θ=0.2°時動應力邊頻帶Fig.13 θ=0.2° side band diagram of dynamic stress

由圖13可知，邊頻帶主要成分為嚙頻和軸頻及其倍頻，該特性可以為系統(tǒng)故障甄別提供思路[17-18]。

動應力均值和幅值隨錯位角θ變化如圖14所示，兩者均隨θ增大呈增大趨勢，錯位為0.2°時動應力均值和幅值較正常情況下分別增大109.7%和441.3%，動應力增大加速零件疲勞，減少零件壽命；同時由圖14可以發(fā)現(xiàn)，θ在0.05°內(nèi)動應力變化較小，θ大于0.05°后增幅較快。

圖14 動應力隨錯位角變化Fig.14 Variation of dynamic stress with misalignment angle

3 結(jié)論

1) 嚙合力、動態(tài)接觸應力和傳動誤差均隨錯位角度增大呈增大趨勢，錯位角度超過0.05°后嚙合力、動態(tài)接觸應力和傳動誤差峰峰值增幅明顯變大，嚙合性能惡化，其中通常錯位角的范圍 控制在2′～3′之內(nèi)以及可通過提高軸承的安裝精度來提高轉(zhuǎn)配的精度。

2) 振動加速度和軸系動應力隨錯位角增大呈增大趨勢，其中振動加速度對嚙合錯位較為敏感，較小的錯位即可造成系統(tǒng)劇烈的振動；由頻域響應可知，隨著錯位角增大，各主頻附近出現(xiàn)由軸頻及其諧波組成的豐富邊頻帶，該特性可為系統(tǒng)故障診斷提供思路。

3) 高速齒輪傳動系統(tǒng)對齒輪嚙合精度要求極高，較小錯位角度即可導致系統(tǒng)動態(tài)特性急劇惡化，對于本系統(tǒng)而言，嚙合錯位角度應控制在極小的錯位角度范圍以內(nèi)。