袁運(yùn)博　劉震　何濤　 陳亞輝　郭宜斌　王東華

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

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

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

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

引 言

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

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

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

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

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

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

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

在波動(dòng)負(fù)載作用下，齒輪傳動(dòng)系統(tǒng)由于齒輪副的嚙合作用，使得系統(tǒng)在產(chǎn)生扭振響應(yīng)的同時(shí)，也會(huì)產(chǎn)生橫振響應(yīng)。由1.3節(jié)的理論模型可知，由于齒輪副的橫振響應(yīng)，導(dǎo)致主從動(dòng)輪間的中心距發(fā)生變化，進(jìn)而導(dǎo)致齒輪副的嚙合角和重合度發(fā)生變化，從而引起1.1節(jié)中嚙合剛度求解公式（1）-（3）中積分下限βΙ發(fā)生變化，導(dǎo)致嚙合剛度產(chǎn)生調(diào)幅現(xiàn)象;另外，橫振引起的重合度變化會(huì)影響嚙合剛度單雙齒嚙合區(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）顯示考慮波動(dòng)負(fù)載時(shí)嚙合剛度出現(xiàn)了以波動(dòng)負(fù)載頻率13 Hz為調(diào)制頻率的邊頻帶，由其中心頻率（一倍嚙合頻率676.67 Hz）幅值大幅降低可知，在波動(dòng)負(fù)載作用下，嚙合剛度不僅產(chǎn)生了調(diào)幅作用，而且伴有明顯的調(diào)頻作用。頻域圖5（a）顯示，中心頻率越高，由波動(dòng)負(fù)載引起的調(diào)制作用越明顯，這主要是因?yàn)橹行妮d波頻率越高，等效扭振效應(yīng)越大，則齒輪的嚙合周期變化就越明顯，進(jìn)而導(dǎo)致嚙合剛度的調(diào)頻作用越明顯。

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

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

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

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

3.3 波動(dòng)負(fù)載幅值影響分析

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

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

4 結(jié) 論

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

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

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

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

參考文獻(xiàn)：

[1] 雷亞國(guó)， 何正嘉， 林 京，等. 行星齒輪箱故障診斷技術(shù)的研究進(jìn)展[J]. 機(jī)械工程學(xué)報(bào)， 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] 李永焯， 丁 康， 何國(guó)林，等. 齒輪系統(tǒng)振動(dòng)響應(yīng)信號(hào)調(diào)制邊頻帶產(chǎn)生機(jī)理[J]. 機(jī)械工程學(xué)報(bào)， 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

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

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