考慮內(nèi)部激勵隨機性的兩級分流式人字齒輪傳動動力學(xué)特性

考慮內(nèi)部激勵隨機性的兩級分流式人字齒輪傳動動力學(xué)特性

廖映華1,2，秦大同1，劉長釗1

(1. 重慶大學(xué)機械傳動國家重點實驗室,重慶400044； 2. 四川理工學(xué)院機械工程學(xué)院,四川自貢643000)

摘要：根據(jù)某特種裝備用兩級分流式人字齒輪傳動系統(tǒng)的構(gòu)型特點，考慮軸承變形、嚙合剛度、嚙合誤差和齒側(cè)間隙等因素，建立了兩級分流式人字齒輪傳動系統(tǒng)的軸承-齒輪耦合非線性動力學(xué)模型。采用Runge-Kutta逐步積分法求解系統(tǒng)的非線性動力學(xué)微分方程，從而獲得隨機嚙合剛度和嚙合誤差激勵作用下兩級分流式人字齒輪傳動系統(tǒng)的動態(tài)嚙合力和動態(tài)支承力及其頻譜，采用Monte Carlo仿真獲得動態(tài)嚙合力和動態(tài)支承力的統(tǒng)計特征，研究了嚙合剛度和嚙合誤差隨機性對動態(tài)嚙合力和動態(tài)支承力的影響規(guī)律，為兩級分流式人字齒輪傳動系統(tǒng)動力學(xué)優(yōu)化以及動態(tài)可靠性優(yōu)化奠定基礎(chǔ)。

關(guān)鍵詞：隨機內(nèi)部激勵；人字齒輪傳動；動力學(xué)特性；統(tǒng)計特征

中圖分類號:TH132.4文獻標志碼：A

Dynamic characteristics of two-stage split herringbone gear trains considering the randomness of internal excitations

LIAOYing-hua1,2,QINDa-tong1,LIUChang-zhao1(1. State key laboratory of mechanical transmission, Chongqing University, Chongqing 400044, China；2. School of Mechanical Engineering, Sichuan University of Science & Engineering, Zigong 643000, China)

Abstract:According to the configurational characteristics, a bearing-gear coupled non-linear dynamics model of two-stage split herringbone gear trains used in special equipment was presented. In the model, bearing stiffness, mesh stiffness，mesh errors and gear backlash were all included. Then the Runge-Kutta step-by-step integration method was used to solve the non-linear differential equations of motion, the time history curves and frequency spectrums of the dynamic meshing force and dynamic bearing force were obtained considering the random mesh stiffness and mesh error excitations, and further, their statistical characteristics were determined by Monte Carlo simulation to analyze the influences of the randomness of mesh stiffness and mesh errors on the dynamic mesh force and dynamic bearing force. The research results lay a foundation for the dynamic optimization and dynamic reliability optimization of two-stage split herringbone gear trains.

Key words:random internal excitation; herringbone gear trains; dynamic characteristic; statistical characteristic

人字齒輪由左右完全對稱的斜齒輪構(gòu)成，人字齒輪傳動無軸向力，承載能力高，廣泛用于起重機械、艦船、礦山機械和工程機械等重載設(shè)備。由于制造安裝誤差、時變嚙合剛度和齒側(cè)間隙等因素的影響，人字齒輪傳動系統(tǒng)工作時各傳動件承受著復(fù)雜的動態(tài)載荷，對系統(tǒng)振動噪聲、疲勞壽命和可靠性有重要影響，因此人字齒輪傳動系統(tǒng)動力學(xué)特性分析是低噪聲和高可靠性設(shè)計的基礎(chǔ)。

已有大量文獻對人字齒輪傳動系統(tǒng)的動力學(xué)特性進行了研究。吳新躍等[1]在考慮齒輪軸的軸向變形，忽略齒輪摩擦的基礎(chǔ)上，建立了單級人字齒輪副嚙合耦合振動模型。王成等[2-3]建立了單級人字齒輪傳動的彎-扭-軸耦合動力學(xué)模型，討論了系統(tǒng)的內(nèi)部激勵和輪齒修形對人字齒輪傳動動態(tài)特性的影響，并通過輪齒接觸分析(TCA)和承載接觸分析(LTCA)計算內(nèi)部激勵。Sondkar[4]建立了單級人字齒輪行星齒輪動力學(xué)分析模型，并分析了系統(tǒng)的固有特性及動態(tài)響應(yīng)。郭家舜等[5]將齒廓修形作為一種時變齒側(cè)間隙，建立了單級人字齒輪傳動的彎-扭-軸耦合動力學(xué)模型，分析了修形參數(shù)對動態(tài)系數(shù)幅值的影響。卜忠紅等[6]建立了人字行星齒輪傳動系統(tǒng)的扭轉(zhuǎn)-橫向振動模型，得到了系統(tǒng)的五種固有振動模式，分析了科氏慣性力、離心慣性力以及轉(zhuǎn)速對系統(tǒng)固有頻率的影響。朱增寶等[7-8]建立了考慮支撐彈性變形、時變嚙合剛度、誤差激勵等因素的差動人字齒輪傳動系統(tǒng)動力學(xué)模型，分析了安裝誤差、偏心和齒輪誤差對系統(tǒng)均載特性的影響。上述研究將各種激勵視為確定性激勵，但實際上由于制造、安裝誤差存在的隨機性導(dǎo)致齒輪系統(tǒng)內(nèi)部激勵的不確定性，其對齒輪系統(tǒng)動力學(xué)特性的影響不應(yīng)忽視[9-13]。因此，為了研究制造誤差隨機性對兩級分流式人字齒輪傳動系統(tǒng)動力學(xué)特性的影響規(guī)律，本文以某特種裝備用兩級分流式人字齒輪傳動系統(tǒng)為對象，考慮軸承變形、嚙合剛度、嚙合誤差和齒側(cè)間隙等因素，建立了人字齒輪傳動系統(tǒng)的齒輪-軸承耦合非線性動力學(xué)模型，分析了由隨機制造誤差引起的嚙合剛度和嚙合誤差等內(nèi)部激勵的隨機性對人字齒輪傳動系統(tǒng)動力學(xué)特性的影響規(guī)律。

圖1為某種特種裝備用兩級分流式人字齒輪傳動系統(tǒng)的運動簡圖，該傳動系統(tǒng)的中間軸軸向固定，高速軸和低速軸軸向自由。將齒輪的輪體視為剛體，輪齒、軸和軸承視為線彈性體， 考慮傳動軸扭轉(zhuǎn)變形和軸向

變形，忽略粘性阻尼、重力和輪齒間的摩擦力，考慮人字齒輪制造誤差引起的嚙合誤差、軸向位移激勵和齒側(cè)間隙，則可建立圖2所示的兩級人字齒輪傳動系統(tǒng)的動力學(xué)模型。圖2中，在齒輪i(i=1,…,6)的理論中心Oi建立固定坐標系OiYiZi，坐標軸Yi平行于系統(tǒng)靜態(tài)齒輪嚙合線，坐標軸Zi位于軸承孔的理想中心線上；齒輪i的質(zhì)量為mi，轉(zhuǎn)動慣量為Ji；Td為驅(qū)動轉(zhuǎn)矩，TL為負載轉(zhuǎn)矩。

圖1　某特種裝備用兩級分流式 人字齒輪傳動系統(tǒng)的運動簡圖 Fig.1 Kinematic diagram of two-stage split herringbone gear trains used on special equipment

圖2　兩級分流式人字齒輪傳動系統(tǒng)的動力學(xué)模型 Fig.2 Dynamic model of two-stage split herringbone gear trains

1兩級分流式人字齒輪傳動系統(tǒng)的動力學(xué)模型

根據(jù)人字齒輪副的嚙合狀態(tài)以及位置關(guān)系，人字齒輪嚙合副k(k=I, Ⅱ, Ⅲ)的左(L)、右(R)側(cè)輪齒的嚙合變形δjk(j=L、R)為

(11)

式中:yi為齒輪i的質(zhì)心在Yi方向上偏離其理論位置的線位移；θi為齒輪i的角位移；βbk為構(gòu)成齒輪副k的齒輪的基圓柱螺旋角；ejk為齒輪副k的左、右側(cè)的嚙合誤差；rbi為齒輪i的基圓半徑。

考慮齒側(cè)間隙時，人字齒輪副k的左、右側(cè)的嚙合力Fjk為

(2)

式中:bk為齒輪副k的法向齒側(cè)間隙。kjk為齒輪副k的左、右側(cè)的綜合嚙合剛度。

人字齒輪副k在切向和軸向上的總嚙合力Fvyk、Fvzk分別為

(3)

取系統(tǒng)廣義坐標列陣為

q=

[y1,θ1,y2,θ2,y3,z3,θ3,y4,z4,θ4,y5,z5,θ5,y6,θ6](4)

根據(jù)牛頓第二定律，圖2所示的兩級分流式人字齒輪傳動系統(tǒng)的運動微分方程為

(5)

2人字齒輪傳動系統(tǒng)內(nèi)部激勵分析

2.1隨機嚙合剛度引起的內(nèi)部激勵

由于齒輪制造安裝誤差隨機性的影響，嚙合剛度既存在周期性分量，又存在隨機性分量，嚙合剛度激勵可表示為

(6)

式中:km為平均嚙合剛度；kai、φki分別為嚙合剛度周期性分量的幅值和相位；N為諧波次數(shù)；fm為嚙合頻率；kr(t)為嚙合剛度隨機性分量。

嚙合剛度隨機性分量kr(t)一般某一范圍內(nèi)變化，工程上通常假設(shè)其服從截尾正態(tài)分布，其概率密度函數(shù)為

(7)

2.2隨機嚙合誤差引起的內(nèi)部激勵

如圖3所示，影響齒輪系統(tǒng)傳動平穩(wěn)性的誤差主要有齒輪的偏心誤差er、齒距偏差fpt、漸開線齒廓偏差ffα。圖2中，O為齒輪的理論中心，O′為齒輪的實際中心，坐標軸y與嚙合線平行。pt為端面齒距，δ為偏心誤差的位置角，fc為軸頻。由于制造、安裝過程中的不確定性，這些誤差均包含隨機成分，則嚙合誤差激勵表示為

e(t)=ersin(2πfct+δ)+

(8)

式中，eai、φei分別為齒距和齒廓偏差引起的周期性分量的幅值和初始相位角；er(t)為齒距和齒廓偏差引起的隨機分量，它服從截尾正態(tài)分布；fm為嚙頻，且fm=fcz，z為齒輪的齒數(shù)。

圖3　影響齒輪系統(tǒng)傳動平穩(wěn)性的誤差 Fig.3 Errors affecting the drive stationarity of a gear system

3仿真與分析

某特種裝備用兩級分流式人字齒輪傳動系統(tǒng)的基本參數(shù)見表1，系統(tǒng)的輸入額定轉(zhuǎn)速nd=1 750 r/min，額定轉(zhuǎn)矩Td=2 225 N·m。以其為例分析嚙合剛度和嚙合誤差隨機性對兩級分流式人字齒輪傳動系統(tǒng)的動力學(xué)特性的影響。

3.1內(nèi)部激勵隨機性對系統(tǒng)動力學(xué)特性的影響

考慮側(cè)隙非線性和隨機內(nèi)部激勵時，兩級分流式人字齒輪傳動系統(tǒng)為多自由度非線性系統(tǒng)。本文采用Runge-Kutta數(shù)值逐步積分求解系統(tǒng)響應(yīng)、動態(tài)嚙合和動態(tài)支承力，并通過傅里葉變換得到它們的頻譜。

表1　某特種裝備用兩級分流式人字齒輪傳動系統(tǒng)的基本參數(shù)

忽略齒側(cè)間隙和嚙合誤差的情況下，嚙合剛度隨機性對齒輪副Ⅰ、Ⅲ的動態(tài)嚙合力的影響如圖4所示，嚙合剛度隨機性對齒輪1的動態(tài)支承力(Fbyk=kykyk)的影響如圖5所示。由圖4可知，考慮嚙合剛度隨機性的動態(tài)嚙合力(圖4(b)、(d))的幅值明顯比不考慮嚙合剛度隨機性的動態(tài)嚙合力(圖4 (a)、(c))大；由于齒輪副Ⅰ的嚙頻高、傳遞的載荷小，嚙合剛度隨機性引起的動態(tài)嚙合力隨機波動較大，因此齒輪副I的動態(tài)嚙合力增幅比齒輪副Ⅲ的增幅更大。由圖4(c)、(d)和圖5(a)、(b)可知，由于軸承-齒輪的耦合作用，動態(tài)支承力的幅值隨動態(tài)嚙合力的增加而增加。由圖4(e)~(h)可知，考慮嚙合剛度隨機性時，動態(tài)嚙合力的各諧波成分的幅值明顯增加，且在圖4(f)所示的齒輪副Ⅰ的動態(tài)嚙合力頻譜中出現(xiàn)了2 668 Hz、2 943 Hz的新諧波成分，這是由于嚙合剛度隨機性引起了動態(tài)嚙合力隨機波動。

圖4　嚙合剛度隨機性對齒輪副Ⅰ、Ⅲ的動態(tài)嚙合力的影響(忽略嚙合誤差和齒側(cè)間隙) Fig.4 Influences of random mesh stiffness on dynamic meshing force of gear pair Ⅰ and Ⅲ (Regardless of the meshing error and backlash)

圖5　嚙合剛度隨機性對齒輪1的動態(tài)支承力的影響(忽略嚙合誤差和齒側(cè)間隙) Fig.5 Influences of random mesh stiffness on dynamic bearing force of gear 1 (Regardless of the meshing error and backlash)

圖6　嚙合剛度隨機性對齒輪副Ⅰ的動態(tài)嚙合力的影響 (齒側(cè)間隙b 1=b 2=15 μm，b 3=25 μm，不考慮嚙合誤差) Fig.6 Influences of random mesh stiffness on dynamic mesh force of gear pair Ⅰ (the backlashes b 1=b 2=15 μm，b 3=25 μm, Regardless of the meshing error)

圖7　嚙合誤差隨機性對齒輪副Ⅰ的動態(tài)嚙合力的影響(忽略嚙合剛度變化和齒側(cè)間隙) Fig.7 Influences of random mesh errors on dynamic mesh force of gear pair Ⅰ (Regardless of the variation of meshing error and backlashes)

圖8　嚙合誤差隨機性對齒輪副Ⅰ的動態(tài)嚙合力的影響 (齒側(cè)間隙b 1=b 2=15 μm，b 3=25 μm，不考慮嚙合剛度變化) Fig.8 Influences of random mesh errors on dynamic mesh force of gear pair Ⅰ (the backlashes b 1=b 2=15 μm，b 3=25 μm，Regardless of the variation of meshing error)

在考慮齒側(cè)間隙，忽略嚙合誤差的情況下，嚙合剛度隨機性對齒輪副Ⅰ的動態(tài)嚙合力的影響如圖6所示。對比圖6(a)、(b)所示的時域曲線可知，嚙合剛度隨機性使齒輪副Ⅰ的動態(tài)嚙合力的隨機波動更加明顯、變化幅值明顯增大，齒輪副不僅發(fā)生脫嚙，導(dǎo)致嚙合沖擊，還發(fā)生了非工作齒面接觸。對比圖6(c)、(d)所示的頻譜可知，考慮嚙合剛度隨機性時，動態(tài)嚙合力和動態(tài)支承力的頻率成分明顯增多，各頻率成分對應(yīng)的幅值也發(fā)生了明顯變化，且在低頻段出現(xiàn)了連續(xù)的頻率成分，這是由于嚙合沖擊和系統(tǒng)內(nèi)部激勵隨機性使動態(tài)嚙合力中隨機成分增加所致。

忽略嚙合剛度變化和齒側(cè)間隙的情況下，嚙合誤差隨機性對齒輪副Ⅰ的動態(tài)嚙合力的影響如圖7所示。由圖7(a)、(b)可知，考慮嚙合誤差隨機性的動態(tài)嚙合力的幅值明顯比不考慮嚙合剛度隨機性的動態(tài)嚙合力大。由圖7(c)、(d)可知，考慮嚙合誤差隨機性時，齒輪Ⅰ的動態(tài)嚙合力出現(xiàn)了2 958 Hz和4 394 Hz的新頻率成分。對比圖7(a)、(b)和圖4(a)、(b)可知，嚙合誤差隨機性對動態(tài)嚙合力幅值的影響比嚙合剛度隨機性更顯著。

在考慮齒側(cè)間隙，忽略嚙合剛度變化的情況下，嚙合誤差隨機性對齒輪副Ⅰ的動態(tài)嚙合力的影響如圖8所示。由圖8(a)、(b)可知，嚙合剛度隨機性使動態(tài)嚙合力的幅值顯著增大，且出現(xiàn)了明顯的脫嚙和反向嚙合。對比圖8(c)、(d)可知，低頻部分的頻率基本不變，高頻部分頻率成分更加豐富。對比圖8(c)、(d)和圖6(c)、(d)可知，嚙合誤差隨機性對動態(tài)嚙合力高頻部分影響較大，而嚙合剛度隨機性對動態(tài)嚙合力的低頻部分和高頻部分均有顯著的影響。

3.2內(nèi)部激勵隨機性對系統(tǒng)動態(tài)力的統(tǒng)計特征的影響

多自由度非線性系統(tǒng)的隨機響應(yīng)的統(tǒng)計特征常采用隨機等效線性化[14]與MCS[15-16](Monte Carlo仿真)獲得，本文采用MCS獲得隨機嚙合剛度和隨機嚙合誤差作用下系統(tǒng)動態(tài)嚙合力和動態(tài)支承力的統(tǒng)計特征。

圖9　系統(tǒng)動態(tài)力的統(tǒng)計特征與齒輪副Ⅲ的隨機嚙合 誤差的標準差的關(guān)系 (齒側(cè)間隙b 1=b 2=15 μm，b 3=25 μm) Fig.9 Relations between the statistical characteristic of dynamic force and the standard deviation of random mesh error on gear pair Ⅲ (The gear backlashes b 1=b 2=15 μm，b 3=25 μm)

圖9為系統(tǒng)的動態(tài)力的統(tǒng)計特征與齒輪副Ⅲ的隨機嚙合誤差的標準差的關(guān)系。由圖9可知，隨機嚙合誤差的標準差較小時，隨機嚙合誤差的標準差對動態(tài)嚙合力和動態(tài)支承力的均值和變異系數(shù)(=標準差/均值)的影響較小；當(dāng)隨機嚙合誤差的標準差較大時，系統(tǒng)振動增大，齒輪副頻繁出現(xiàn)脫嚙和非工作齒面接觸，使得動態(tài)支承力的均值、動態(tài)嚙合力和動態(tài)支承力的變異系數(shù)隨隨機嚙合誤差的標準差的增大而顯著增大。

隨機嚙合剛度的標準差對系統(tǒng)動態(tài)力的統(tǒng)計特征的影響規(guī)律與隨機嚙合誤差的標準差相似。在不考慮外部載荷變化的情況下，增大隨機嚙合剛度和隨機嚙合誤差的標準差，系統(tǒng)的動態(tài)嚙合力的均值基本不變。

4結(jié)論

(1)建立了考慮軸承變形，嚙合剛度、嚙合誤差和齒側(cè)間隙的兩級分流式人字齒輪傳動系統(tǒng)的軸承-齒輪耦合動力學(xué)模型。

(2)采用Runge-Kutta數(shù)值逐步積分法求得隨機嚙合剛度或嚙合誤差激勵作用下系統(tǒng)的動態(tài)嚙合力和動態(tài)支承力，結(jié)果表明：嚙合剛度或嚙合誤差的隨機性使得系統(tǒng)的動態(tài)嚙合力和動態(tài)支承力發(fā)生隨機波動，諧波成分增多，甚至出現(xiàn)連續(xù)頻率成分；嚙合剛度或嚙合誤差的隨機性使動態(tài)嚙合力和動態(tài)支承力明顯增大，當(dāng)轉(zhuǎn)速高、載荷輕時，嚙合剛度或嚙合誤差的隨機性對動態(tài)嚙合力和動態(tài)支承力影響更明顯；考慮嚙合剛度或嚙合誤差的隨機性時，齒側(cè)間隙對系統(tǒng)動力學(xué)特性的影響更加顯著。

(3)采用Monte Carlo仿真獲得隨機嚙合剛度或嚙合誤差激勵作用下系統(tǒng)的動態(tài)嚙合力和動態(tài)支承力的統(tǒng)計特征，并分析了其對動態(tài)嚙合力和動態(tài)支承力的統(tǒng)計特征的影響。結(jié)果表明：隨機嚙合剛度或嚙合誤差的標準差增大，動態(tài)嚙合力和動態(tài)支承力的標準差以及動態(tài)支承力的均值也增大，但動態(tài)嚙合力的均值基本不變；當(dāng)隨機嚙合剛度或嚙合誤差的標準差增大至某一值時，嚙合副將發(fā)生脫嚙和非工作齒面接觸，引起嚙合沖擊，使動態(tài)嚙合力和動態(tài)支承力的波動幅值急劇增大。

參考文獻

[1]吳新躍，朱石堅. 人字齒輪傳動的振動理論分析模型[J]. 海軍工程大學(xué)學(xué)報, 2001, 13(05): 13-19.

WU Xin-yue, ZHU Shi-jian. Theoretical vibration pattern of herringbone gearing[J]. Journal of Naval University of Engineering, 2005, 13 (5):13-19.

[2]王成，方宗德，張墨林，等. 人字齒輪傳動的動態(tài)特性分析[J]. 哈爾濱工業(yè)大學(xué)學(xué)報, 2011, 43(07): 122-126.

WANG Cheng, FANG Zong-de, ZHANG Mo-lin, et al. Analysis of dynamics behavior for double helical gears transmissions[J]. Journal of Harbin Institute of Technology, 2011, 43 (7):122-126.

[3]王成，高常青，崔煥勇. 基于嚙合特性的人字齒輪動力學(xué)建模與分析[J]. 中南大學(xué)學(xué)報(自然科學(xué)版), 2012, 20(8): 3019-3024.

WANG Cheng, GAO Chang-qing, CUI Huan-yong.Dynamic modeling and analysis of double helical gears based on meshing characters[J]. Journal of Central South University (Science and Technology), 2012, 43 (8):3019-3024.

[4]Prashant S. Dynamic modeling of double-helical planetary gear sets[D]. Columbus: Ohio State University, 2012.

[5]郭家舜，王三民，王穎. 基于解析解的人字齒輪齒廓修形動態(tài)分析[J]. 航空動力學(xué)報, 2013, 28(3): 613-620.

GUO Jia-shun, WANG San-min, WANG Ying. Dynamic analysis of tooth profile modification to herringbone gears based on analytical solutions[J]. Journal of Aerospace Power, 2013, 28 (3): 613-620.

[6]卜忠紅，劉更，吳立言，等. 封閉式人字齒輪行星傳動的自由振動特性分析[J]. 中國機械工程, 2009, 20(20): 2420-2424.

BU Zhong-hong, LIU Geng, WU Li-yan, et al. Free vibration characterization analysis of encased differential herringbone gear train[J]. China Mechanical Engineering, 2009, 20 (20):2420-2424.

[7]朱增寶，朱如鵬，鮑和云，等. 偏心與齒頻誤差對封閉差動人字齒輪傳動系統(tǒng)動態(tài)均載特性的影響分析[J]. 航空動力學(xué)報,2011, 26(11): 2601-2609.

ZHU Zeng-bao, ZHU Ru-peng, Bao He-yun, et al. Impact of run-out and meshing-frequency errors on dynamic load sharing for encased differential herringbone train[J].Journal of Aerospace Power, 2011, 26 (11): 2601-2609.

[8]朱增寶，朱如鵬，李應(yīng)生，等. 安裝誤差對封閉差動人字齒輪傳動系統(tǒng)動態(tài)均載特性的影響[J]. 機械工程學(xué)報, 2012, 48(3): 16-24.

ZHU Zeng-bao,ZHU Ru-peng, LI Ying-sheng, et al. Impact of installation error on dynamics load sharing characteristic for encased differential herringbone train[J]. Journal of Mechanical Engineering, 2012，48 (3):16-23.

[9]陳會濤，秦大同，吳曉鈴，等. 考慮隨機制造誤差的風(fēng)力機行星齒輪系統(tǒng)動力學(xué)特性[J]. 機械工程學(xué)報, 2012, 48(21): 77-83.

CHEN Hui-tao, QIN Da-tong, WU Xiao-ling, et al. Dynamic characteristics of planetary gear transmission system of wind turbine with random manufacturing error[J]. Journalof Mechanical Engineering, 2012, 48 (11):77-83.

[10]魏永祥，陳建軍，馬洪波. 隨機參數(shù)齒輪系統(tǒng)的非線性動力響應(yīng)分析[J]. 工程力學(xué), 2012,29(11): 319-324.

WEI Yong-xiang, CHEN Jian-jun, MA Hong-bo. Analysis of nonlinear dynamic response of gear-rotor with random parameters [J]. Engineering Mechanics, 2012, 29(11): 319-324.

[11]盧劍偉，曾凡靈，楊漢生，等. 隨機裝配側(cè)隙對齒輪系統(tǒng)動力學(xué)特性的影響分析[J]. 機械工程學(xué)報, 2010, 46(21): 82-86.

LU Jian-wei, ZENG Fan-ling, YANG Han-sheng, et al.Influence of stochastic assembling backlash on nonlinear dynamic behavior of transmission gear pair[J].Journal of Mechanical Engineering, 2010, 46(21): 82-86.

[12]盧劍偉，劉夢軍，陳磊，等. 隨機參數(shù)下齒輪非線性動力學(xué)行為[J]. 中國機械工程, 2009, 10(3): 330-333.

LU Jian-wei, LIU Meng-jun, CHEN Lei, et al.Nonl inear dynamics behavior of gear system with stochastic parameters[J].China Mechanical Engineering, 2009, 10(3): 330-333.

[13]魏永祥，陳建軍，王敏娟. 隨機參數(shù)齒輪-轉(zhuǎn)子系統(tǒng)的扭轉(zhuǎn)振動分析[J]. 航空動力學(xué)報, 2010, 25(11): 2637-2642.

WEI Yong-xiang, CHEN Jian-jun, WANG Min-juan. Dynamic response of torsional vibration of gear-rotor with random parameters[J]. Journal of Aerospace Power,2010, 25(11): 2637-2642.

[14]何勇,金偉良. 多自由度柔性結(jié)構(gòu)非線性隨機振動響應(yīng)分析方法[J]. 計算力學(xué)學(xué)報, 2009, 26(6): 778-783.

HE Yong, JIN Wei-liang.Nonlinear random vibration of multiple-degree-of-freedom flexible structures [J]. Chinese Journal of Computational Mechanics, 2009, 26(6): 778-783.

[15]Schu?ller G I. Computational stochastic mechanics-recent advances[J]. Computers & Structures, 2001, 79(22-25): 2225-2234.

[16]Pradlwarter H J, Schu?ller G I. On advanced monte carlo simulation procedures in stochastic structural dynamics[J]. International Journal of Non-Linear Mechanics, 1997, 32(4): 735-744.