譚秀峰，張國偉，謝里陽，何雪浤

(東北大學(xué)機械工程與自動化學(xué)院，遼寧沈陽110004)

0 引言

齒輪傳動系統(tǒng)動力學(xué)是研究齒輪傳動系統(tǒng)在動態(tài)激勵作用下動力學(xué)行為的一門工程科學(xué)［1］。齒輪傳動系統(tǒng)在車輛、電力系統(tǒng)、礦山設(shè)備等行業(yè)的重要設(shè)備中發(fā)揮著關(guān)鍵的作用。這些設(shè)備的安全、穩(wěn)定、可靠運行關(guān)系到經(jīng)濟效益及社會生活的各個方面。大型化工企業(yè)如果停產(chǎn)一天，其損失可達上百萬元［2］。我國水泥行業(yè)中應(yīng)用較廣泛的水泥磨齒輪箱，其故障可造成水泥每年至少減產(chǎn)200萬噸［3］。

齒輪傳動系統(tǒng)有3個明顯的特點:

(1)系統(tǒng)的轉(zhuǎn)速高。有的轉(zhuǎn)速甚至達到了數(shù)10萬轉(zhuǎn)/分，這會產(chǎn)生嚴(yán)重的振動，因此，振動問題一直是系統(tǒng)動力學(xué)的研究重點;

(2)建模困難。所涉及到的機械零部件有多級齒輪副、轉(zhuǎn)子和軸承等，從傳動結(jié)構(gòu)上來分有變速箱、原動機和負載等;

(3)求解困難。齒輪的間隙、時變剛度和齒面摩擦等，軸承的間隙和不同結(jié)構(gòu)形式帶來的不同剛度形式等，轉(zhuǎn)子(軸)的不同支撐形式和間隙帶來的松動等，這些非線性因素使得系統(tǒng)的求解變得非常困難，如果再考慮到齒輪傳動系統(tǒng)中可能存在的各種故障，則模型更為復(fù)雜［4］。

這些問題給齒輪傳動系統(tǒng)動力學(xué)分析、設(shè)計帶來了巨大的困難。

本研究主要根據(jù)國內(nèi)外的研究成果，闡述齒輪傳動系統(tǒng)非線性動力學(xué)理論的基本框架，對齒輪傳動系統(tǒng)的非線性動力學(xué)研究方法進行綜述;對含多級齒輪傳動的系統(tǒng)整體動力學(xué)研究進行分析總結(jié)，特別對典型的車用變速器動力學(xué)的問題進行分析;最后，針對變速器動力學(xué)，提出其在動態(tài)特性、振動和噪聲、故障診斷以及可靠性方面可能的研究熱點和方向。

1 齒輪傳動系統(tǒng)的非線性動力學(xué)研究

在齒輪傳動系統(tǒng)動力學(xué)的研究中，主要包括的3方面內(nèi)容［5］，齒輪傳動系統(tǒng)理論體系如圖1所示。

圖1 齒輪傳動系統(tǒng)理論體系

近30年來，非線性動力學(xué)的發(fā)展使人們不斷把新觀點、新方法引入到齒輪傳動系統(tǒng)的動力學(xué)分析中。許多線性方法無法解釋的分岔、跳躍、極限環(huán)、混沌等現(xiàn)象終于被揭示出來［6］。

非線性問題包括兩種模型:

(1)非線性時不變模型。齒輪傳動過程中，由于制造、安裝以及潤滑的需要，間隙是必然存在的，間隙可以是齒側(cè)間隙或滾動軸承間隙［7-10］，這種非線性模型沒有考慮嚙合剛度的時變性［11-13］;

(2)非線性時變模型。這種模型在考慮間隙非線性的同時，也考慮了嚙合剛度的時變性，從而把齒輪系統(tǒng)作為非線性的參數(shù)振動問題進行研究［14］。

國內(nèi)外許多學(xué)者已針對齒輪-轉(zhuǎn)子系統(tǒng)的動力學(xué)特性進行了大量的研究，研究結(jié)果表明:影響齒輪-轉(zhuǎn)子-滑動軸承系統(tǒng)振動特性的非線性因素主要包括齒輪嚙合力和滑動軸承的油膜力等［15］。另外，轉(zhuǎn)子和聯(lián)軸器等系統(tǒng)的材料及幾何非線性、滾動軸承的間隙、內(nèi)阻尼、內(nèi)摩擦力，密閉環(huán)中的流體作用力等，也是非線性的重要來源。

目前，關(guān)于齒輪傳動系統(tǒng)的非線性動力學(xué)研究主要有3種方法:即近似解析法、數(shù)值法和實驗法，現(xiàn)將10年來各學(xué)者的研究成果綜述如下。

1.1 解析法

解析法難以研究初值對系統(tǒng)非線性行為的影響，這是近似解析法的共同缺點，一般僅能用來考慮齒輪非線性嚙合力的單自由度和兩個自由度的系統(tǒng)［16］。

1.1.1 分段技術(shù)

分段技術(shù)是將齒輪非線性系統(tǒng)按照區(qū)間進行分塊，在每一個區(qū)間形成一個時變線性系統(tǒng)再進行求解。

柴山等［17］發(fā)現(xiàn)只有分段段數(shù)足夠多時，用分段法求得的結(jié)果才與積分法達到一致的結(jié)果。王雷［18］建立了一般的綜合考慮時變剛度、慣量的非圓齒輪扭轉(zhuǎn)振動動力學(xué)微分方程組，提出了逐段線性近似的求解方法，并采用方程的解析解來代替數(shù)值解法。劉曉寧［19］分別建立了具有分段線性特征的單自由度非線性動力學(xué)模型、綜合考慮齒側(cè)間隙和時變嚙合剛度的三自由度非線性動力學(xué)模型。

1.1.2 能量法

振動量的大小直接取決于能量輸入的強度，用振動功率流對結(jié)構(gòu)振動特性進行評價，可較傳統(tǒng)的動力響應(yīng)分析更好地揭示出振動系統(tǒng)的綜合動態(tài)特性，因而受到越來越多的重視。

邵毅敏［20］提出了能量保持因子概念，用來描述其能量通過多界面的損耗關(guān)系，并通過實驗獲得了多界面的能量關(guān)系與能量保持因子的系數(shù)表。馮婧［21］在動力學(xué)建模的基礎(chǔ)上探討了應(yīng)用能量分析研究系統(tǒng)非線性振動特性的方法，隨后推導(dǎo)了系統(tǒng)運動方程和振動功率流函數(shù)的時域仿真算法。陳思雨［22］提出了用改進的能量法求解其周期響應(yīng)的方法，詳細推導(dǎo)了該法的求解過程，得到了系統(tǒng)振動幅頻特性曲線。

1.1.3 增量諧波平衡法(IHB)

增量諧波平衡法是在諧波平衡法的基礎(chǔ)上發(fā)展起來的。增量諧波平衡法可以彌補諧波平衡法解的精度取決于諧波項個數(shù)的缺點，對于一般的非線性系統(tǒng)可以方便地求取任意階近似解。

楊紹普等［23-24］利用IHB，對直齒輪副的非線性動力學(xué)進行了研究，考慮了時變嚙合剛度和間隙，得到了這類模型統(tǒng)一形式的解，同時研究了外激勵幅值、阻尼比以及系統(tǒng)的分岔特性對系統(tǒng)幅頻曲線的影響。劉振皓［25］運用IHB對包含時變嚙合剛度、齒側(cè)間隙與綜合嚙合誤差的Ravigneaux式復(fù)合行星齒輪傳動系統(tǒng)純扭轉(zhuǎn)動力學(xué)微分方程組進行求解，得到了系統(tǒng)的基頻穩(wěn)態(tài)響應(yīng)。目前還出現(xiàn)了一些改進的增量諧波平衡法［26-27］。

1.2 數(shù)值法

數(shù)值方法是通過數(shù)值法求解非線性微分方程，得到非線性系統(tǒng)在特定初始條件和參數(shù)下的運動規(guī)律。但是該方法不能給出解的表達式，故無法對系統(tǒng)的整體作定性分析。但數(shù)值方法可作為一種驗證非線性振動問題近似解析解或者發(fā)現(xiàn)新現(xiàn)象的主要方法。

尚志勇［28］采用Newmark方法求解了含有松動與碰摩耦合故障的轉(zhuǎn)子-軸承系統(tǒng)系統(tǒng)的響應(yīng)，用Poincare映射、軸心軌跡和頻譜圖分析各個轉(zhuǎn)子-軸承系統(tǒng)在特定參數(shù)下的運動特征。陳會濤［29-30］等基于齒輪動力學(xué)理論和Lagrange方程，通過數(shù)值仿真得到了載荷和參數(shù)隨機變異時系統(tǒng)各響應(yīng)量和齒輪副間動態(tài)嚙合力的統(tǒng)計特征。陳小安［31］提出了一種基于有限單元法的多間隙耦合齒輪傳動系統(tǒng)的非線性動態(tài)特性分析方法。

1.3 實驗法

實驗方法作為一種重要的輔助方法，不僅可用來驗證理論研究的正確性，也可用來修正理論分析模型。但因為受試驗中各種誤差的影響，通過試驗來精確研究齒輪系統(tǒng)的非線性動力學(xué)行為也具有一定的困難。

馬銳等［32］通過試驗提取含有裂紋故障齒輪的振動特征，驗證了理論分析的結(jié)果。周長江［33］從理論和實驗2個方面對復(fù)雜潤滑狀態(tài)下齒面摩擦因數(shù)的計算方法作了深入的系統(tǒng)研究。韓志武等［34］以工程仿生學(xué)和有限元理論為基礎(chǔ)，并通過激光雕刻技術(shù)將仿生表面形態(tài)加工在齒面上，進行實際的臺架試驗。

2 含多級齒輪傳動的系統(tǒng)整體動力學(xué)研究

對于多級齒輪傳動系統(tǒng)，除了動力學(xué)模型的自由度數(shù)量增加以外，還存在整個齒輪系統(tǒng)包括各級齒輪軸、各級齒輪、軸承、箱體及原動機等多種振動耦合的綜合作用以及各種內(nèi)部激勵的非線性特性，這些都給系統(tǒng)的動態(tài)分析帶來了困難。

林騰蛟等［35］建立了錐-平行軸-行星多級齒輪傳動系統(tǒng)包含時變嚙合剛度、嚙合阻尼等因素的18自由度彎-扭-軸耦合非線性動力學(xué)模型，采用4～5階變步長Runge-Kutta法對動力學(xué)微分方程進行了求解。崔亞輝等［36］提出一種新的求解方法，可用該方法求出多級齒輪系統(tǒng)的動態(tài)嚙合力，采用數(shù)值仿真方法求解了系統(tǒng)的動態(tài)頻率響應(yīng)。

在現(xiàn)階段，多級齒輪傳動系統(tǒng)整體模型的降階處理，系統(tǒng)中各連接部位、邊界處、軸承處的剛度和阻尼以及動載荷處理，尋找出其動力學(xué)計算的一般方法，具有重要的理論意義和現(xiàn)實意義。

特別需要指出的是，汽車變速器作為一種典型的多級齒輪傳動系統(tǒng)，學(xué)者們已經(jīng)對其動力性能進行了一些研究［37-39］。除了需要考慮上述問題外，汽車變速器還需要考慮換擋過程中的沖擊力對變速器動力性能的影響［40］。換擋過程還涉及到操縱機構(gòu)、同步器和離合器的配合，研究它們自身參數(shù)對換擋力與換擋時間的影響同樣增加了問題的復(fù)雜性［41］。隨著變速器檔位數(shù)的增多，如何降低噪聲、減少震動和加強換擋連續(xù)性以提高駕駛員換擋舒適感是動力學(xué)亟需解決的問題。

3 結(jié)束語

本研究根據(jù)對國內(nèi)外齒輪傳動系統(tǒng)動力學(xué)研究成果的分析總結(jié)，闡述了齒輪傳動系統(tǒng)非線性動力學(xué)理論的基本框架;綜述了10年來齒輪傳動系統(tǒng)的非線性動力學(xué)的研究方法，包括解析法、數(shù)值法和實驗法;分析總結(jié)了含多級齒輪傳動的整體系統(tǒng)的動力學(xué)研究現(xiàn)狀，其重點問題在于模型的降階、內(nèi)部激勵和邊界條件的處理;特別對典型的車用變速器動力學(xué)問題進行了分析，如何降低噪聲、減少震動和加強換擋連續(xù)性以提高駕駛員換擋舒適感是其動力學(xué)亟需解決的問題。

隨著齒輪傳動系統(tǒng)非線性動力學(xué)研究的深入，以下幾個方面可能成為變速器動力學(xué)研究的熱點問題:

(1)變速器動態(tài)性能的全面分析。由于整個變速器的振動是軸、齒輪、軸承、箱體及原動機等多種振動耦合的綜合作用，研究者可利用有限元、實驗?zāi)B(tài)分析技術(shù)等方法進行整個變速器系統(tǒng)的建模及結(jié)構(gòu)設(shè)計，對系統(tǒng)動力性能進行優(yōu)化。

(2)變速器系統(tǒng)的振動和噪聲控制。研究者通過對系統(tǒng)輸出響應(yīng)的研究，了解振動機理，從而設(shè)計出平穩(wěn)性好和噪聲低的變速器系統(tǒng)，提高系統(tǒng)的有效壽命。

(3)變速器系統(tǒng)非線性動力學(xué)中的故障診斷問題。研究者通過研究變速器帶故障的齒輪系統(tǒng)，以及齒輪嚙合剛度、傳遞誤差和嚙合沖擊等與輪齒破壞的關(guān)系，對其頻譜機理進行分析，可對故障診斷進行指導(dǎo)，在理論和工程上均具有重要的價值。

(4)變速器系統(tǒng)非線性動力學(xué)中的可靠性問題。復(fù)雜零件受力的概率分布函數(shù)的確定、對系統(tǒng)考慮載荷多次作用、強度退化的可靠性模型建立和對具有失效相關(guān)性模式的系統(tǒng)可靠度的計算等問題都是進行可靠性分析亟待解決的問題。研究者可在分別建立軸承系統(tǒng)、齒輪傳動系統(tǒng)、離合器或聯(lián)軸器以及齒輪箱體的可靠性計算模型的基礎(chǔ)上，開展變速器整體的可靠性研究。

［1］李潤方，王建軍.齒輪系統(tǒng)動力學(xué)—振動、沖擊、噪聲［M］.北京:科學(xué)出版社，1997.

［2］李明，孫濤，胡海巖.齒輪傳動轉(zhuǎn)子—軸承系統(tǒng)動力學(xué)的研究進展［J］.振動工程學(xué)報，2002，15(3):249-256.

［3］陳克興，李川奇.設(shè)備狀態(tài)監(jiān)測與故障診斷技術(shù)［M］.北京:科學(xué)技術(shù)文獻出版社，1991.

［4］申永軍，楊紹普，李偉.齒輪系統(tǒng)非線性動力學(xué)研究進展及展望［J］.石家莊鐵道學(xué)院學(xué)報，2005，18(4):5-10.

［5］王建軍，李其漢，李潤方.齒輪系統(tǒng)非線性振動研究進展［J］.力學(xué)進展，2005，35(1):37-51.

［6］陳予恕.非線性振動［M］.北京:高等教育出版社，2002.

［7］KAHRAMAN A，SINGH R.Interactions between time-varying mesh stiffness and clearance non-linearities in a geared system［J］.Journal of Sound and Vibration，1991，146(1):135-156.

［8］CAI Y.Simulation on the rotaional 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.

［9］KAHRAMAN A，BLANKENSHIP G W.Interactions between commensurate parametric and forcing excitations in a system with clearance［J］.Journal of Sound and Vibration，1996，194(3):317-336.

［10］THEODOSSIADES S，NATSIAVAS S.Non-linear dynamics of gear-pair systems with periodic stiffness and backlash［J］.Journal of Sound and Vibration，2000，229(2):287-310.

［11］KAHRAMAN A，SINGH R.Nonlinear dynamics of a spur gear pair［J］.Journal of Sound and Vibration，1990，142(1):49-75.

［12］KAHRAMAN A，SINGH R.Nonlinear dynamics of a geared rotor-bearing system with multiple clearances［J］.Journal of Sound and Vibration，1991，144(1):469-506.

［13］唐增寶，朱毅芳.齒輪傳動的振動分析與動態(tài)優(yōu)化設(shè)計［M］.武漢:華中理工大學(xué)出版社，1994.

［14］王立華，李潤方，林騰蛟，等.齒輪系統(tǒng)時變剛度和間隙非線性振動特性研［J］.中國機械工程，2003，14(13):1143-1146.

［15］OZGUVEN H N.A non-linear mathematical model for dynamic analysis of spur gears including shaft and bearing dynamics［J］.Journal of Sound and Vibration，1991，145(2):239-260.

［16］崔亞輝.齒輪—轉(zhuǎn)子—滑動軸承系統(tǒng)非線性動力學(xué)特性的理論和試驗研究［D］.哈爾濱:哈爾濱工業(yè)大學(xué)能源科學(xué)與工程學(xué)院.2009.

［17］柴山，高連勇.用積分法計算直齒輪輪齒的彈性變形［J］.機械設(shè)計，2005，22(2):48-50.

［18］王雷.非圓齒輪動力學(xué)特性研究［D］.鄭州:鄭州大學(xué)機械工程學(xué)院，2008.

［19］劉曉寧.齒輪系統(tǒng)的混沌控制及仿真［D］.西安:西北工業(yè)大學(xué)機電學(xué)院，2007.

［20］邵毅敏，陳再剛，周曉君，等.沖擊振動能量通過“齒輪-軸-軸承-軸承座”多界面?zhèn)鬟f損耗研究［J］.振動與沖擊，2009，28(6):2-9.

［21］馮婧.齒輪非線性系統(tǒng)的數(shù)值仿真與功率流分析［D］.濟南:山東大學(xué)機械工程學(xué)院，2009.

［22］陳思雨.一種新的齒輪非線性振動數(shù)學(xué)模型建模與分析求解研究［D］.長沙:中南大學(xué)機電工程學(xué)院，2007.

［23］SHEN Yong-jun，YANG Shao-pu，PAN Cun-zhi，et al.Nonlinear dynamics of a spur gear pair with time-varying stiffness and backlash［J］.Journal of Low Frequency Noise，Vibration and Active Control，2004，23(3):178-187.

［24］申永軍，楊紹普，潘存治，等.參外聯(lián)合激勵下直齒輪副的非線性動力學(xué)［J］.北京交通大學(xué)學(xué)報，2005，29(1):69-73.

［25］劉振皓，巫世晶，王曉筍，等.基于增量諧波平衡法的復(fù)合行星齒輪傳動系統(tǒng)非線性動力學(xué)［J］.振動與沖擊，2012，31(3):117-122.

［26］高陽，王三民，劉曉寧.一種改進的增量諧波平衡法及其在非線性振動中的應(yīng)用［J］.機械科學(xué)與技術(shù)，2005，24(6):663-665.

［27］嚴(yán)剛峰，黃顯核.基于遺傳算法的改進諧波平衡算法［J］.電子測量與儀器學(xué)報，2009，23(10):96-100.

［28］尚志勇.轉(zhuǎn)子—軸承系統(tǒng)的非線性動力學(xué)分析［D］.哈爾濱:哈爾濱工業(yè)大學(xué)能源科學(xué)工程學(xué)院，2006.

［29］陳會濤.風(fēng)力發(fā)電機齒輪傳動系統(tǒng)隨機振動分析及動力可靠性概率優(yōu)化設(shè)計［D］.重慶:重慶大學(xué)機械傳動國家重點實驗室，2012.

［30］周志剛.隨機風(fēng)作用下風(fēng)力發(fā)電機齒輪傳動系統(tǒng)動力學(xué)及動態(tài)可靠性研究［D］.重慶:重慶大學(xué)機械傳動國家重點實驗室，2012.

［31］陳小安.基于有限單元法的多間隙耦合齒輪傳動系統(tǒng)非線性動態(tài)特性分析［J］.振動與沖擊，2010，29(2):46-49.

［32］馬銳，陳予恕.含裂紋故障齒輪系統(tǒng)的非線性動力學(xué)研究［J］.機械工程學(xué)報，2011，47(21):84-90.

［33］周長江，唐進元，鐘志華，等.齒輪傳動齒面摩擦因數(shù)計算方法的研究［J］.潤滑與密封，2006(10):185-191.

［34］韓志武，呂尤，馬榮峰，等.仿生表面形態(tài)對齒輪動力學(xué)性能的影響［J］.北京工業(yè)大學(xué)學(xué)報，2011，37(6):806-810.

［35］林騰蛟.多級齒輪傳動系統(tǒng)耦合非線性振動特性分析［J］.振動與沖擊，2013，32(17):1-7

［36］崔亞輝，劉占生，葉建槐，等.復(fù)雜多級齒輪-轉(zhuǎn)子-軸承系統(tǒng)的動力學(xué)建模和數(shù)值仿真［J］.機械傳動，2009，33(6):44-48.

［37］閔海濤，高娟，馬天飛.汽車變速器箱體結(jié)構(gòu)強度分析與優(yōu)化設(shè)計［J］.中國機械工程，2012，23(20):2514-2519.

［38］李強軍.重型汽車變速器直齒行星副變速器動力學(xué)研究［D］.重慶:重慶大學(xué)機械傳動國家重點實驗室，2012.

［39］寇海江.汽車變速器齒輪傳動系統(tǒng)非線性動力學(xué)特性研究［D］.沈陽:東北大學(xué)理學(xué)院，2011.

［40］陳玉祥，臧孟炎，陳勇，等.基于虛擬樣機技術(shù)的手動變速器換擋力分析［J］.中國機械工程，2012，23(8):996-1000.

［41］崔傳寶.機械式變速器換擋性能評價方法的初步研究［D］.長春:吉林大學(xué)汽車工程學(xué)院，2011.