黃冠華， 張衛(wèi)華， 付永佩， 梁樹林， 王興宇

高速列車齒輪傳動(dòng)系統(tǒng)參數(shù)振動(dòng)穩(wěn)定性

黃冠華1， 張衛(wèi)華1， 付永佩1， 梁樹林2， 王興宇2

（1.西南交通大學(xué)牽引動(dòng)力國家重點(diǎn)實(shí)驗(yàn)室，四川成都610031；2.中國北車集團(tuán)長(zhǎng)春軌道客車股份有限公司，吉林長(zhǎng)春130024）

為了準(zhǔn)確表達(dá)參數(shù)激勵(lì)下高速列車齒輪系統(tǒng)振動(dòng)的穩(wěn)定性，利用有限元方法得到高速列車齒輪系統(tǒng)時(shí)變嚙合剛度，并用傅里葉級(jí)數(shù)展開進(jìn)行擬合.考慮齒輪嚙合誤差，建立了高速列車齒輪傳動(dòng)系統(tǒng)扭轉(zhuǎn)振動(dòng)模型.結(jié)合多尺度近似解析方法，推導(dǎo)了參激振動(dòng)下高速列車齒輪系統(tǒng)的近似解析解，得到了系統(tǒng)的穩(wěn)定性邊界曲線，并分析了影響齒輪傳動(dòng)系統(tǒng)穩(wěn)定性的相關(guān)因素.研究結(jié)果表明：齒輪系統(tǒng)的不穩(wěn)定性區(qū)域隨著列車運(yùn)行的速度降低總體呈減小趨勢(shì)，但是在發(fā)生參數(shù)共振速度處存在明顯不穩(wěn)定區(qū)域；增大阻尼有利于系統(tǒng)的穩(wěn)定性，當(dāng)阻尼系數(shù)從0.01增加到0.05時(shí)，處于穩(wěn)定區(qū)域的剛度波動(dòng)幅值從5%增加至20%；增加齒輪的重合度可以減小嚙合剛度的諧波特性，從而增強(qiáng)系統(tǒng)的穩(wěn)定性.

穩(wěn)定性；參數(shù)振動(dòng)；齒輪傳動(dòng)系統(tǒng)；多尺度法；高速列車

高速列車傳動(dòng)系統(tǒng)一般為齒輪傳動(dòng)系統(tǒng)，從動(dòng)輪直接壓裝在車軸上，主動(dòng)輪采用聯(lián)軸節(jié)與牽引電機(jī)相連，通過齒輪箱懸吊在構(gòu)架橫梁上［1］.因此，齒輪系統(tǒng)的振動(dòng)特性將直接影響著高速列車驅(qū)動(dòng)傳動(dòng)系統(tǒng)的運(yùn)行性能.在齒輪傳動(dòng)中，由于參與嚙合的輪齒對(duì)數(shù)的周期變化，使得齒輪輪齒的綜合嚙合剛度也是周期變化的，所以在動(dòng)力學(xué)模型中體現(xiàn)為周期性時(shí)變的彈性剛度.考慮這種因素后，齒輪動(dòng)力學(xué)問題在力學(xué)上是參數(shù)振動(dòng)問題，其分析的模型是參數(shù)振動(dòng)方程［2］.

齒輪系統(tǒng)是一種參數(shù)振動(dòng)系統(tǒng)，判定系統(tǒng)穩(wěn)定性以及影響穩(wěn)定性的因素是首要問題，眾多齒輪方面的學(xué)者在這方面做了大量工作.文獻(xiàn)［3-5］中利用數(shù)值方法分析了單自由度齒輪系統(tǒng)的穩(wěn)定性和穩(wěn)態(tài)響應(yīng)，模型將齒輪嚙合剛度假設(shè)為矩形波和正弦波，揭示了齒輪系統(tǒng)的諧波共振、概周期響應(yīng)和混沌響應(yīng).文獻(xiàn)［6-7］中利用直接積分算法對(duì)多自由度齒輪系統(tǒng)的穩(wěn)定性進(jìn)行了研究，并探討了通向混沌的雙周期分岔途徑.文獻(xiàn)［8-10］中運(yùn)用Floquet理論對(duì)穩(wěn)定性尤其是參數(shù)穩(wěn)定性方面也有較多的研究［8-10］.

上述研究主要集中在一般機(jī)械領(lǐng)域，針對(duì)鐵路車輛，尤其是高速列車傳動(dòng)系統(tǒng)的研究較少.

高速列車齒輪傳動(dòng)系統(tǒng)動(dòng)態(tài)響應(yīng)較為復(fù)雜，在特定頻率激勵(lì)往往出現(xiàn)超諧共振、亞諧共振等多種參數(shù)共振形式，對(duì)齒輪傳動(dòng)系統(tǒng)的服役將產(chǎn)生不利影響，嚴(yán)重的會(huì)導(dǎo)致系統(tǒng)的共振失效.為了準(zhǔn)確表達(dá)參數(shù)激勵(lì)下的高速列車齒輪系統(tǒng)振動(dòng)穩(wěn)定性，本文針對(duì)某型高速動(dòng)車組齒輪傳動(dòng)系統(tǒng)，采用多尺度解析方法對(duì)齒輪系統(tǒng)方程作近似展開，得到系統(tǒng)穩(wěn)定區(qū)的近似解析解，給出動(dòng)力穩(wěn)定性圖譜，并從系統(tǒng)穩(wěn)定性角度提出了高速列車齒輪傳動(dòng)系統(tǒng)參數(shù)選取建議.

1 齒輪嚙合動(dòng)力學(xué)模型

1.1 動(dòng)力學(xué)方程

當(dāng)忽略傳動(dòng)軸和支撐系統(tǒng)的彈性變形時(shí)，可將高速列車齒輪傳動(dòng)系統(tǒng)簡(jiǎn)化處理為齒輪副的扭轉(zhuǎn)振動(dòng)系統(tǒng)，如圖1所示，圖中：

θp、θg為主、被動(dòng)齒輪的扭轉(zhuǎn)振動(dòng)角位移；

Ip、Ig為主、被動(dòng)齒輪的轉(zhuǎn)動(dòng)慣量；

Rp、Rg為主、被動(dòng)齒輪的基圓半徑；

i為傳動(dòng)比；

e（t）為輪齒嚙合傳遞誤差；

km為嚙合綜合剛度；

cm為嚙合阻尼；

Tp、Tg為作用在主、被動(dòng)齒輪上的外載荷力矩.

動(dòng)力學(xué)方程可表示為［11］

圖1 齒輪動(dòng)力學(xué)模型Fig.1 Dynamic model of a gear pair

為了消除系統(tǒng)剛性位移，定義系統(tǒng)動(dòng)態(tài)傳遞誤差和靜態(tài)傳遞誤差的差值為

將式（1）、（2）相減，得到單自由度系統(tǒng)的動(dòng)力學(xué)方程為

式中：

me為當(dāng)量質(zhì)量，

1.2 齒輪嚙合剛度和傳遞誤差

齒輪嚙合剛度的獲取方法有很多種，通常先計(jì)算出嚙合剛度的峰值和平均值，然后按嚙合頻率將嚙合剛度簡(jiǎn)化成矩形波周期函數(shù)，略去高階項(xiàng)后再將其展開成傅里葉級(jí)數(shù)［12］，即

式中：

Ks為平均剛度；

Kj為剛度波動(dòng)幅值；

j為剛度有限諧波項(xiàng)數(shù)；

φj為相位角；

ωe為齒輪副的嚙合圓周頻率，

ωe＝2πZ1n1／60＝2πZ2n2／60，

其中，

Z1、Z2分別為主、被動(dòng)輪的齒數(shù)，

n1、n2分別為主、被動(dòng)輪的轉(zhuǎn)速.

上述方法簡(jiǎn)單易用，但是對(duì)于輪齒剛度的時(shí)變表達(dá)不夠精確.事實(shí)上，輪齒綜合嚙合剛度定義為使一對(duì)或幾對(duì)同時(shí)嚙合的輪齒在1 mm齒寬上產(chǎn)生1 μm撓度所需的載荷［13］.根據(jù)這一定義，建立高速列車輪齒三維實(shí)體接觸有限元模型，本文建立的齒輪傳動(dòng)的輪齒接觸有限元模型如圖2所示.

圖2 輪齒接觸的三維有限元模型Fig.2 3D finite element model of gear pair contact

在可能接觸區(qū)域部分進(jìn)行網(wǎng)格加密，得到的模型共有68 460個(gè)單元，86 750個(gè)節(jié)點(diǎn).將主動(dòng)輪和被動(dòng)輪的齒面定義為接觸對(duì)，在齒輪軸線上建立參考點(diǎn)，并在參考點(diǎn)和大小齒輪內(nèi)圈和剖面間建立耦合約束，將轉(zhuǎn)矩載荷、約束施加在主動(dòng)輪和被動(dòng)輪的參考點(diǎn)上.計(jì)算隨時(shí)間變化的嚙合輪齒之間彈性變形和受力，得到齒輪嚙合剛度，并采用傅里葉級(jí)數(shù)對(duì)時(shí)變剛度進(jìn)行擬合.

圖3為小齒輪在4 200 r／min轉(zhuǎn)速下，采用有限元方法和傅里葉級(jí)數(shù)擬合的齒輪時(shí)變嚙合剛度曲線.

輪齒嚙合誤差通常采用齒輪嚙合頻率的傅里葉級(jí)數(shù)表示［12］，即

式中：

e0、ej分別為齒輪誤差的常數(shù)和幅值；

θj為相位角.

1.3 運(yùn)動(dòng)方程的無量綱化

將式（4）、（5）代入式（3），令x＝bu（b為特征尺寸），對(duì)其進(jìn)行無量綱化，可得

式中：K1＝Kj／Ks；

圖3 時(shí)變嚙合剛度曲線Fig.3 Time-varying curve of mesh stiffness

2 穩(wěn)定性分析

引入小參數(shù)ε，則有：

ξ＝εμ； kj＝εK1.

式（6）的齊次形式可表示為¨u＋2εμ˙u＋

使用多尺度法［14］，討論式（7）的一次近似解.設(shè)零次近似方程的解為

式中：

T0、T1、Ti為多尺度法的時(shí)間變量；

A為待定的復(fù)函數(shù).

一次近似方程表示為

式中：

cc為前面各項(xiàng)的共軛復(fù)數(shù).

將式（10）代入到式（9）中，消除久期項(xiàng)，得

設(shè)

分離實(shí)部和虛部，得

式中：

其中，

b1、b2為常數(shù)，

λ為特征值.

特征方程為

由式（13）可知，當(dāng)λ具有正實(shí)部時(shí)，系統(tǒng)不穩(wěn)定，由此可得系統(tǒng)穩(wěn)定性邊界的臨界曲線方程為

3 實(shí)例分析

根據(jù)上述結(jié)果，對(duì)某型高速列車齒輪系統(tǒng)的穩(wěn)定性進(jìn)行分析.齒輪副的相關(guān)參數(shù)如表1所示.

表1 齒輪副參數(shù)Tab.1 Parameters of gear couples

圖4（a）為不考慮嚙合阻尼，展開項(xiàng)數(shù)j取1～6項(xiàng)時(shí)，系統(tǒng)在kj-ˉω平面上的穩(wěn)定性圖譜，V形區(qū)域內(nèi)為不穩(wěn)定區(qū)域（以下同）.

從圖4（a）可以看出，隨著項(xiàng)數(shù)j的增大，嚙合剛度的諧波特性會(huì)降低，系統(tǒng)的不穩(wěn)定性區(qū)逐漸減??；在嚙合剛度不變時(shí)，隨著參數(shù)激勵(lì)ˉω的減小，不穩(wěn)定區(qū)域也會(huì)減小，出現(xiàn)不穩(wěn)定區(qū)域的重疊.

圖4（b）為相應(yīng)的速度穩(wěn)定性圖譜.

從圖4（b）可以看出，齒輪嚙合頻率隨著列車運(yùn)行速度的降低而減小，不穩(wěn)定的區(qū)域總體呈減小趨勢(shì)，但在發(fā)生參激振動(dòng)的轉(zhuǎn)速時(shí)，不穩(wěn)定的區(qū)域明顯更大，在實(shí)際運(yùn)行中應(yīng)引起注意.

圖4 齒輪系統(tǒng)穩(wěn)定性Fig.4 Stability of gear system

圖5 為阻尼系數(shù)對(duì)穩(wěn)定性的影響.從圖5中可以看出，系統(tǒng)阻尼對(duì)穩(wěn)定性有較大的影響，阻尼可以減小系統(tǒng)的不穩(wěn)定區(qū)域，改善系統(tǒng)的動(dòng)態(tài)特性.當(dāng)阻尼系數(shù)從0.01增加到0.05時(shí)，處于穩(wěn)定區(qū)域的剛度波動(dòng)幅值從5%增加至20%.

從以上分析可以看出，可以通過以下途徑增加高速列車齒輪傳動(dòng)系統(tǒng)穩(wěn)定性：降低嚙合諧波剛度比值、合理選取系統(tǒng)的參激頻率（嚙合頻率與固有頻率的匹配）以及增大輪齒的嚙合阻尼.在實(shí)際的設(shè)計(jì)過程中，首先應(yīng)該保證列車常用的運(yùn)行速度避開固有頻率與嚙合頻率容易發(fā)生參激共振時(shí)的轉(zhuǎn)速（從圖4（b）看應(yīng)盡量避免240 km／h的常速行駛），增大嚙合阻尼主要依靠材料的選取或采用附加阻尼的方式，降低嚙合諧波剛度比值可以通過增大齒輪的嚙合重合度.

式（15）是端面重合度與嚙合剛度均值的表達(dá)式［15］，

式中：

εα為端面重合度；

c′為單對(duì)齒剛度.

圖6為端面重合度分別取1.2和1.9時(shí)，采用數(shù)值直接積分對(duì)式（6）進(jìn)行求解得到的位移隨時(shí)間變化圖，從圖中可以看出，當(dāng)端面重合度分別取1.2和1.9時(shí)，系統(tǒng)趨于不穩(wěn)定和穩(wěn)定，這也是高速列車齒輪系統(tǒng)參數(shù)設(shè)計(jì)中普遍采用高重合度的原因.

圖5 不同阻尼下齒輪系統(tǒng)參數(shù)振動(dòng)穩(wěn)定性Fig.5 Parametric vibration stability of gear system with different dampings

圖6 不同端面重合度系統(tǒng)時(shí)間歷程圖Fig.6 Dynamic response time history of gear system with different transverse contact ratios

4 結(jié) 論

本文從理論上分析了高速列車齒輪系統(tǒng)在參數(shù)時(shí)變嚙合剛度下的穩(wěn)定性問題，通過對(duì)時(shí)變嚙合剛度的傅里葉展開，運(yùn)用非線性多尺度近似解析方法得到了齒輪系統(tǒng)的穩(wěn)定性圖譜，從系統(tǒng)穩(wěn)定性的角度得到了齒輪設(shè)計(jì)時(shí)應(yīng)考慮參激頻率、剛度的諧波特性以及嚙合阻尼三方面因素，主要結(jié)論如下：

（1）齒輪系統(tǒng)的不穩(wěn)定性區(qū)域隨著列車運(yùn)行的速度降低總體呈減小趨勢(shì)，但是在參激頻率處存在明顯不穩(wěn)定區(qū)域，應(yīng)根據(jù)系統(tǒng)的固有頻率合理地制定運(yùn)營速度.

（2）系統(tǒng)的阻尼比和嚙合剛度的諧波分量對(duì)系統(tǒng)的穩(wěn)定性影響較大.增大阻尼有利于系統(tǒng)的穩(wěn)定性，通過增加齒輪嚙合的重合度可以減小嚙合剛度的諧波特性，從而減小系統(tǒng)的不穩(wěn)定區(qū)域，當(dāng)端面重合度從1.2增加到1.9，對(duì)系統(tǒng)直接數(shù)值積分也驗(yàn)證了這一結(jié)果.

（3）文中從穩(wěn)定性的角度分析了齒輪嚙合引起的參數(shù)振動(dòng)，分析模型可為高速列車驅(qū)動(dòng)傳動(dòng)系統(tǒng)動(dòng)力學(xué)等研究提供借鑒.

［1］ 張衛(wèi)華.動(dòng)車組總體與轉(zhuǎn)向架［M］.北京：中國鐵道出版社，2011：197-204.

［2］ 王建軍，洪濤，吳仁智，等.齒輪系統(tǒng)參數(shù)振動(dòng)問題研究綜述［J］.振動(dòng)與沖擊，1997，16（4）：69-73.WANG Jianjun，HONG Tao，WU Renzhi，et al.Researches on parametric vibration of gear transmission systems-review［J］.Journal of Vibration and Shock，1997，16（4）：69-73.

［3］ BENTON M，SEIREG A.Factors influencing instability and resonance in geared systems［J］.ASME Journal of Mechanical Design，1981，103：372-378.

［4］ KAHRAMAN A，SINGH R.Interactions between timevarying mesh stiffness and clearance non-linearities in a geared systems［J］.ASMEJournalofSoundand Vibration，1991，146（1）：135-156.

［5］ KAHRAMAN A.Effect of axial vibrations on the dynamics of a helical gear pair［J］.Journal of Sound and Acoustics，1993，115（1）：33-39.

［6］ RAGHOTHAMA A，NARAYANAN S.Bifurcation and chaos in geared rotor bearing system by incremental harmonic balance method［J］.Journal of Sound and Vibration，1999，226（3）：469-492.

［7］ SEYRANIAN A P，SOLEM F，PEDEREEN P.Multiparameter linear periodic systems：sensitivity analysis and applications［J］.Journal of Sound and Vibration，2000，229（1）：89-111.

［8］ 李同杰，朱如鵬，鮑和云，等.行星齒輪傳動(dòng)非線性振動(dòng)系統(tǒng)參數(shù)穩(wěn)定域的計(jì)算方法［J］.航空動(dòng)力學(xué)報(bào)，2012，27（6）：1416-1423.LI Tongjie，ZHU Rupeng，BAO Heyun，et al.Method of stability region determination for planetary gear train's parameters based on nonlinear vibration model［J］.Journal of Aerospace Power，2012，27（6）：1416-1423.

［9］ 周麗杰，王能建，張德福.艦載直升機(jī)甲板牽引系統(tǒng)運(yùn)動(dòng)穩(wěn)定性分析［J］.西南交通大學(xué)學(xué)報(bào)，2011，46（3）：409-414.ZHOU Lijie，WANG Nengjian，ZHANG Defu.Motion stability analysis of carrier helicopter traction system on deck［J］.Journal of Southwest Jiaotong University，2011，46（3）：409-414.

［10］ 王延忠，周元子，李國權(quán)，等.螺旋錐齒輪嚙合剛度及參數(shù)振動(dòng)穩(wěn)定性研究［J］.航空動(dòng)力學(xué)報(bào)，2010，25（7）：1664-1669.WANG Yanzhong，ZHOU Yuanzi，LI Guoquan，et al.Mesh stiffnessfunctionandstabilityanalysisfor parametric vibration of spiral bevel gears［J］.Journal of Aerospace Power，2010，25（7）：1664-1669.

［11］ 李潤(rùn)方，王建軍.齒輪系統(tǒng)動(dòng)力學(xué)［M］.北京：科學(xué)出版社，1997：154-162.

［12］ AMABILI M，RIVOLA A.Dynamic analysis of spur gear pairs：steady-state response and stability of the SDOF model with time-varying meshing damping［J］.Mechanical Systems and Signal Processing，1997，11（3）：375-390.

［13］ SENEY S.Torsional properties of spur gears in mesh using nonlinear finite element analysis［D］.Perth：Curtin University of Technology，1999.

［14］ 劉延柱，陳立群.非線性振動(dòng)［M］.北京：高等教育出版社，2001：83-86.

［15］ 國家技術(shù)監(jiān)督局.GB／T3480—1997漸開線圓柱齒輪承載能力計(jì)算方法［S］.北京：中國標(biāo)準(zhǔn)出版社，1997.

（中文編輯：秦 瑜 英文編輯：蘭俊思）

Stability Analysis of Parametric Vibration for Gear Transmission System in High-Speed Train

HUANG Guanhua1， ZHANG Weihua1， FU Yongpei1， LIANG Shulin2， WANG Xingyu2
（1.State Key Laboratory of Traction Power，Southwest Jiaotong University，Chengdu 610031，China；2.CNR，Changchun Railway Passenger Vehicle Company，Changchun 130024，China）

In order to study the stability of the gear transmission system in high-speed trains，a dynamic model describing the torsional vibration behaviors of the gear system was developed.In this model，the time-varying mesh stiffness of meshing teeth pairs was calculated through finite element analysis，and the mesh stiffness and transmission error were expanded using the technique of Fourier series.Based on this model，the multiple scales method was used to solve the nonlinear differential equations of gear systems，and the approximate analytical solution and transition curves that separate stable from unstable regions were obtained.In addition，the main factors that influence the stability were discussed.The results show that the unstable regions decrease with the decrease of the train＇s running speed，but an unstable region always exists at the speed where parametric resonance occurs；increasing the damping is effective to reduce the unstable regions：as the damping increases from 0.01 to 0.05，the amplitude of mesh stiffness fluctuation in stable regions increases from 5%to 20%；and，an increase in the contact ratio can help suppress the harmonic characteristics of mesh stiffness so as to improve the stability of system.

stability；parametric vibration；gear transmission system；method of multiple scales；high-speed train

U270.1

：A

0258-2724（2014）06-1010-06

10.3969／j.issn.0258-2724.2014.06.012

2013-09-18

國家自然科學(xué)基金和鐵道部高速鐵路基礎(chǔ)研究基金聯(lián)合資助項(xiàng)目（U1234208）

黃冠華（1987－），男，博士研究生，研究方向?yàn)楦咚倭熊噦鲃?dòng)系統(tǒng)動(dòng)力學(xué)，E-mail：hgh7735＠126.com

黃冠華，張衛(wèi)華，付永佩，等.高速列車齒輪傳動(dòng)系統(tǒng)參數(shù)振動(dòng)穩(wěn)定性［J］.西南交通大學(xué)學(xué)報(bào)，2014，49（6）：1010-1015.