波動負載對齒輪系統(tǒng)振動特性及邊頻調(diào)制影響研究

袁運博　劉震　何濤　 陳亞輝　郭宜斌　王東華

摘要： 齒輪傳動系統(tǒng)在工作時常承受復(fù)雜多變的外部波動載荷，導(dǎo)致輪齒嚙合特性和系統(tǒng)振動響應(yīng)頻率特征復(fù)雜。基于時變嚙合剛度的能量法合成模型，建立考慮系統(tǒng)扭振和橫振響應(yīng)影響的時變嚙合剛度動態(tài)修正模型。建立單級齒輪傳動系統(tǒng)的彎扭耦合模型，用Newmark法求解系統(tǒng)的振動響應(yīng)。利用嚙合剛度動態(tài)修正模型和齒輪系統(tǒng)彎扭耦合模型，通過數(shù)值算例分析波動負載對嚙合剛度和系統(tǒng)振動響應(yīng)的影響。結(jié)果表明，在波動負載作用下，嚙合剛度和系統(tǒng)振動響應(yīng)均存在明顯的以波動負載頻率為調(diào)制頻率的邊頻調(diào)制現(xiàn)象，且被調(diào)制的中心諧波頻率越高，調(diào)制現(xiàn)象越明顯;外部波動負載的幅值越大，嚙合剛度和系統(tǒng)振動響應(yīng)的調(diào)制現(xiàn)象越明顯，且當(dāng)波動負載幅值較小時，表現(xiàn)為窄帶調(diào)頻和調(diào)幅的疊加，嚙合頻率兩側(cè)僅各有一條明顯的邊頻譜線。

關(guān)鍵詞： 齒輪傳動系統(tǒng); 外部波動負載; 邊頻調(diào)制; 嚙合剛度; 頻率特征

中圖分類號：TH132.4;TH113.1? 文獻標(biāo)志碼： A? 文章編號： 1004-4523（2019）03-0526-08

DOI：10.16385/j.cnki.issn.1004-4523.2019.03.018

引 言

齒輪傳動系統(tǒng)因其傳動平穩(wěn)、工作可靠、傳遞效率高等特點，被廣泛應(yīng)用于交通運輸、能源動力、石油化工和機械制造等領(lǐng)域。齒輪傳動系統(tǒng)的運行工況及其所承受的外部載荷復(fù)雜多變，導(dǎo)致齒輪傳動系統(tǒng)的動態(tài)響應(yīng)復(fù)雜多變，難以準(zhǔn)確預(yù)測。另外，由復(fù)雜多變的外部載荷引起的振動響應(yīng)與由齒輪內(nèi)部故障引起的振動響應(yīng)彼此相互耦合，導(dǎo)致診斷齒輪內(nèi)部故障的困難程度增大[1]。

不同內(nèi)部故障對嚙合剛度變化和系統(tǒng)響應(yīng)等的影響機理研究對齒輪傳動系統(tǒng)的故障監(jiān)測與振動具有重要意義。Chaari等[2]建立了輪齒脫落和輪齒局部斷裂兩種典型故障對時變嚙合剛度的影響模型，研究表明這兩種故障會造成時變嚙合剛度的衰減，從而造成嚙合頻率及其倍頻處系統(tǒng)響應(yīng)幅值的降低，并且會產(chǎn)生邊頻現(xiàn)象。Saxena等[3]在計入摩擦力影響的條件下，研究了輪齒脫落的形狀、大小以及位置三個因素對時變嚙合剛度和系統(tǒng)動態(tài)響應(yīng)的影響，發(fā)現(xiàn)輪齒脫落的三個因素對嚙合剛度和系統(tǒng)響應(yīng)具有顯著影響。Jia等[4]分別建立了包含齒輪局部脫落和裂紋的動力學(xué)模型，對比分析了兩種故障下系統(tǒng)振動響應(yīng)信號調(diào)制邊頻帶的特征，發(fā)現(xiàn)相對于齒根裂紋故障，局部脫落引起的信號調(diào)制邊頻帶幅值變化可以忽略。

目前很多學(xué)者集中研究了各種單一齒輪內(nèi)部故障下嚙合剛度和系統(tǒng)振動響應(yīng)的變化規(guī)律[2-7]，李永焯等[8]建立了平穩(wěn)型故障下齒輪系統(tǒng)振動響應(yīng)的通用求解模型，直觀地解釋了系統(tǒng)響應(yīng)的調(diào)制邊頻的特征規(guī)律和誘發(fā)原因。

雖然目前有關(guān)各種齒輪內(nèi)部故障下動力學(xué)特性的研究有著較為堅實的積累，但是針對外部波動負載作用下齒輪傳動系統(tǒng)動力學(xué)特性的研究卻相對匱乏。Sika和Velex[9]以及Qiu等[10]分別研究了在外部轉(zhuǎn)速波動作用下，由于扭振而引起的嚙合剛度調(diào)頻作用對系統(tǒng)參數(shù)穩(wěn)定性的影響。Kim等[11]提出了一種考慮橫振響應(yīng)影響的齒輪動力學(xué)模型，研究了由于橫振響應(yīng)而引起的時變壓力角和時變重合度對系統(tǒng)振動響應(yīng)的影響。但是目前尚未有文獻研究波動負載作用下，同時計入橫振和扭振響應(yīng)對齒輪嚙合剛度和系統(tǒng)振動響應(yīng)的影響。然而研究振動響應(yīng)與載荷波動之間的潛在關(guān)系，建立對載荷波動不敏感的齒輪箱故障診斷方法，是解決齒輪箱故障診斷精度低的一條有效途徑[1]。

本文針對上述問題，提出了一種同時計入系統(tǒng)橫振和扭振響應(yīng)影響的時變嚙合剛度修正模型，以研究在外部波動負載作用下，齒輪嚙合剛度和系統(tǒng)振動響應(yīng)的時域變化特性及其邊頻調(diào)制特征規(guī)律。

3.1 波動負載作用下嚙合剛度變化規(guī)律

不考慮波動負載作用和考慮波動負載作用下的齒輪時變嚙合剛度時域?qū)Ρ葓D如圖4所示，頻域?qū)Ρ葓D如圖5所示。時域圖4（a）顯示考慮波動負載時嚙合剛度伴有明顯的調(diào)幅作用，時域圖4（b）顯示考慮波動負載時嚙合剛度伴有一定的調(diào)頻作用。

在波動負載作用下，齒輪傳動系統(tǒng)由于齒輪副的嚙合作用，使得系統(tǒng)在產(chǎn)生扭振響應(yīng)的同時，也會產(chǎn)生橫振響應(yīng)。由1.3節(jié)的理論模型可知，由于齒輪副的橫振響應(yīng)，導(dǎo)致主從動輪間的中心距發(fā)生變化，進而導(dǎo)致齒輪副的嚙合角和重合度發(fā)生變化，從而引起1.1節(jié)中嚙合剛度求解公式（1）-（3）中積分下限βΙ發(fā)生變化，導(dǎo)致嚙合剛度產(chǎn)生調(diào)幅現(xiàn)象;另外，橫振引起的重合度變化會影響嚙合剛度單雙齒嚙合區(qū)占比，從而引起調(diào)頻和調(diào)幅綜合調(diào)制現(xiàn)象。由1.3節(jié)的理論模型也可知，由于齒輪副的扭振響應(yīng)導(dǎo)致嚙合周期發(fā)生變化，使嚙合剛度產(chǎn)生調(diào)頻現(xiàn)象。

頻域圖5（b）顯示考慮波動負載時嚙合剛度出現(xiàn)了以波動負載頻率13 Hz為調(diào)制頻率的邊頻帶，由其中心頻率（一倍嚙合頻率676.67 Hz）幅值大幅降低可知，在波動負載作用下，嚙合剛度不僅產(chǎn)生了調(diào)幅作用，而且伴有明顯的調(diào)頻作用。頻域圖5（a）顯示，中心頻率越高，由波動負載引起的調(diào)制作用越明顯，這主要是因為中心載波頻率越高，等效扭振效應(yīng)越大，則齒輪的嚙合周期變化就越明顯，進而導(dǎo)致嚙合剛度的調(diào)頻作用越明顯。

3.2 波動負載作用下系統(tǒng)振動響應(yīng)變化規(guī)律

不考慮波動負載作用和考慮波動負載作用下的系統(tǒng)的振動響應(yīng)時域?qū)Ρ葓D如圖6所示，頻域?qū)Ρ葓D如圖7和8所示。其中，圖7為橫振響應(yīng)對比，圖8為扭振響應(yīng)對比。本部分所給出的振動響應(yīng)均為主動輪的振動響應(yīng)，從動輪的振動響應(yīng)具有相似規(guī)律，本文不再贅述。

時域圖6（a）顯示在波動負載作用下主動輪的橫振（圖3中y向）響應(yīng)存在明顯的調(diào)幅現(xiàn)象，時域圖6（b）顯示主動輪的扭振響應(yīng)也存在明顯的調(diào)幅現(xiàn)象。另外，雖然主動輪的橫振和扭振響應(yīng)均伴有調(diào)頻現(xiàn)象，但在時域圖中卻不易觀察。對比圖6（a）和6（b）可知，主動輪的扭振響應(yīng)明顯大于其橫振響應(yīng)，這是由于加載的波動負載為扭轉(zhuǎn)方向的激勵。

圖7和8分別為主動輪的橫振和扭振響應(yīng)的頻譜圖。由圖7和8可知，相對于不考慮波動負載對嚙合剛度影響的情況，考慮波動負載對嚙合剛度影響時，主動輪的橫振和扭振響應(yīng)均出現(xiàn)了明顯的邊頻調(diào)制現(xiàn)象。對比圖5可以發(fā)現(xiàn)，主動輪上的橫振和扭振響應(yīng)的邊頻調(diào)制規(guī)律與嚙合剛度的相似。但同時也應(yīng)注意到，在齒輪嚙合頻率高次諧波附近，響應(yīng)的各個邊頻成分的相對幅值與嚙合剛度的差別較大，這主要是因為齒輪的振動響應(yīng)不僅僅受嚙合剛度的影響，同時還受到傳遞誤差的影響。

3.3 波動負載幅值影響分析

齒輪箱在實際運行過程中可能承受不同幅值的波動負載，因此有必要分析在不同幅值波動負載作用下，齒輪嚙合剛度和系統(tǒng)振動響應(yīng)的影響規(guī)律特征，明晰齒輪傳動系統(tǒng)振動響應(yīng)頻譜特征，為不同載荷條件下齒輪箱的故障診斷和在線監(jiān)測提供理論支撐。

本節(jié)假定小波動負載、中波動負載和大波動負載所對應(yīng)的波動扭矩幅值分別為30，100和200 N·m。圖9為嚙合剛度一倍嚙合頻區(qū)頻譜對比圖，圖10和11分別為主動輪的橫振和扭振響應(yīng)頻譜對比圖。由圖9-11可知，波動負載的幅值越大，嚙合剛度和系統(tǒng)響應(yīng)的調(diào)制越明顯，邊頻越豐富。而且當(dāng)外部波動負載較小時，調(diào)制現(xiàn)象表現(xiàn)為窄帶調(diào)頻和調(diào)幅特性，如圖9-11中小波動負載條件下，嚙合剛度以及系統(tǒng)響應(yīng)的嚙合頻率兩側(cè)僅各有一根明顯的邊頻譜線。小波動負載作用時，系統(tǒng)的振動響應(yīng)較小，因此由橫振響應(yīng)引起的重合度的改變量以及由扭振響應(yīng)引起的嚙合周期的改變量較小，則由橫振和扭振響應(yīng)引起的調(diào)頻作用較小，故呈現(xiàn)出窄帶調(diào)頻現(xiàn)象。

4 結(jié) 論

本文基于齒輪時變嚙合剛度的能量法合成模型，建立了同時考慮系統(tǒng)橫振和扭振響應(yīng)影響的時變嚙合剛度動態(tài)修正模型，同時計入原動機和負載作用，建立了單級齒輪傳動系統(tǒng)的彎扭耦合模型，并通過Newmark法求解了系統(tǒng)的振動響應(yīng)。利用嚙合剛度動態(tài)修正模型和齒輪系統(tǒng)彎扭耦合模型，通過數(shù)值算例分析了外部波動負載對齒輪嚙合剛度和系統(tǒng)振動響應(yīng)的影響，并得出如下結(jié)論。

（1）與不考慮外部波動負載所引起的橫振和扭振響應(yīng)對嚙合剛度的影響相比，考慮外部波動負載時，齒輪嚙合剛度存在明顯的以波動負載頻率為調(diào)制頻率的邊頻調(diào)制現(xiàn)象，在嚙合頻率及其諧波倍頻附近均出現(xiàn)了明顯的調(diào)制邊頻帶，且中心諧波頻率越高，調(diào)制現(xiàn)象越明顯。

（2）考慮外部波動負載所引起的橫振和扭振響應(yīng)對嚙合剛度的影響時，系統(tǒng)的振動響應(yīng)存在與齒輪嚙合剛度相似的邊頻調(diào)制現(xiàn)象，但是在齒輪嚙合頻率高次諧波附近，響應(yīng)的各個邊頻成分的相對幅值與嚙合剛度的差別較大。

（3）外部波動負載的幅值越大，嚙合剛度和系統(tǒng)振動響應(yīng)的調(diào)制現(xiàn)象越明顯，并且當(dāng)波動負載幅值小于一定值時，表現(xiàn)為窄帶調(diào)頻和調(diào)幅的疊加，嚙合頻率兩側(cè)僅各有一根明顯的邊頻譜線。

參考文獻：

[1] 雷亞國， 何正嘉， 林 京，等. 行星齒輪箱故障診斷技術(shù)的研究進展[J]. 機械工程學(xué)報， 2011， 47（19）： 59-67.

Lei Yaguo， He Zhengjia， Lin Jing， et al. Research advances of fault diagnosis technique for planetary gearboxes[J]. Journal of Mechanical Engineering， 2011， 47（19）： 59-67.

[2] Chaari F， Baccar W， Abbes M S， et al. Effect of spalling or tooth breakage on gearmesh stiffness and dynamic response of a one-stage spur gear transmission[J]. European Journal of Mechanics-A/Solids， 2008， 27（4）： 691-705.

[3] Saxena A， Parey A， Chouksey M. Time varying mesh stiffness calculation of spur gear pair considering sliding friction and spalling defects[J]. Engineering Failure Analysis， 2016， 70： 200-211.

[4] Jia S， Howard I. Comparison of localised spalling and crack damage from dynamic modelling of spur gear vibrations[J]. Mechanical Systems and Signal Processing， 2006， 20（2）： 332-349.

[5] Ma H， Song R， Pang X， et al. Time-varying mesh stiffness calculation of cracked spur gears[J]. Engineering Failure Analysis， 2014， 44（6）： 179-194.

[6] Ma H， Li Z， Feng M， et al. Time-varying mesh stiffness calculation of spur gears with spalling defect[J]. Engineering Failure Analysis， 2016， 66： 166-176.

[7] Ma H， Zeng J， Feng R， et al. Review on dynamics of cracked gear systems[J]. Engineering Failure Analysis， 2015， 55： 224-245.

[8] 李永焯， 丁 康， 何國林，等. 齒輪系統(tǒng)振動響應(yīng)信號調(diào)制邊頻帶產(chǎn)生機理[J]. 機械工程學(xué)報， 2018， 54（5）： 105-112.

Li Yongzhuo， Ding Kang， He Guolin， et al. Modulation sidebands of the vibration signal of gearbox[J]. Journal of Mechanical Engineering， 2018， 54（5）： 105-112.

[9] Sika G， Velex P. Instability analysis in oscillators with velocity-modulated time-varying stiffness—Applications to gears submitted to engine speed fluctuations[J]. Journal of Sound and Vibration， 2008， 318（1）： 166-175.

[10]Qiu X， Han Q， Chu F. Analytical investigation on unstable vibrations of single-mesh gear systems with velocity-modulated time-varying stiffness[J]. Journal of Sound and Vibration， 2014， 333（20）： 5130-5140.

[11]Kim W， Hong H Y， Chung J. Dynamic analysis for a pair of spur gears with translational motion due to bearing deformation[J]. Journal of Sound and Vibration， 2010， 329（21）： 4409-4421.

[12]Liang X， Zuo M J， Patel T H. Evaluating the time-varying mesh stiffness of a planetary gear set using the potential energy method[J]. Proceedings of the Institution of Mechanical Engineers， Part C： Journal of Mechanical Engineering Science， 2014， 228（3）： 535-547.

[13]Bartelmus W. Mathematical modelling and computer simulations as an aid to gearbox diagnostics[J]. Mechanical Systems and Signal Processing， 2001， 15（5）： 855-871.

[14]Cai Y. Simulation on the rotational vibration of helical gears in consideration of the tooth separation phenomenon （a new stiffness function of helical involute tooth pair）[J]. Journal of Mechanical Design， 1995， 117（3）： 460-469.

Abstract： Gear transmission is often subjected to complex fluctuating external torques， which usually results in complex characteristics of gear meshing and complicated frequency feature of dynamic response. Based on the potential energy method model of time-varying mesh stiffness， a modified model of time-varying mesh stiffness is proposed， which takes the influence of the torsional and translational vibration into account. The translational-torsional coupling model of a one-stage gear set is developed， and the dynamic vibration response is computed by the Newmark method. Numerical simulations are implemented to investigate the influence of the fluctuating external torque on the dynamic response and modulation sidebands. The results indicate that distinct sidebands appear in both the spectra of the mesh stiffness and the dynamic response when the effect of the fluctuating external torque is considered， and the modulation frequency of those sidebands is no other than the frequency of the fluctuating external torque. What's more， the higher the central frequency is， the more obvious the modulation phenomenon is. The results also show that more obvious and richer sidebands appear with a larger fluctuating external torque. Meanwhile， a superposition of narrow-band frequency modulation and amplitude modulation is discovered， resulting in only one significant sideband on either side of the mesh frequency.

Key words： gear transmission; fluctuated external torque; modulation sidebands; mesh stiffness; frequency feature

作者簡介：袁運博（1993-），男，博士研究生。電話：13206553155;E-mail：yuanyunbo163@163.com

通訊作者：王東華（1981-），男，副教授。電話：13633606817;E-mail：hitwdh@163.com