陳銳博, 張建杰, 周建星, 孫文磊

(新疆大學(xué) 機械工程學(xué)院，烏魯木齊 830047)

含軸向偏載的行星齒輪傳動系統(tǒng)動態(tài)特性研究

陳銳博, 張建杰, 周建星, 孫文磊

(新疆大學(xué) 機械工程學(xué)院，烏魯木齊 830047)

軸向偏載是影響行星輪傳動系統(tǒng)工作狀態(tài)的重要因素，以行星輪傳動系統(tǒng)為研究對象，運用牛頓動力學(xué)理論建立了行星傳動平移-扭轉(zhuǎn)動力學(xué)模型，討論了軸向偏載對嚙合剛度的影響。在考慮輪齒齒面接觸特性基礎(chǔ)上，揭示了軸向偏載對行星傳動系統(tǒng)動態(tài)特性的影響。研究表明：太陽輪發(fā)生軸向偏載后，齒面接觸區(qū)域產(chǎn)生偏移，其嚙合剛度也隨之發(fā)生改變。隨著偏載的增大，嚙合剛度逐漸減小，對行星傳動系統(tǒng)的高階固有頻率影響較為明顯。太陽輪與行星輪的外嚙合動載荷將出現(xiàn)轉(zhuǎn)頻成分，并且當(dāng)偏載增大時，其轉(zhuǎn)頻成分將明顯增大，但其他頻率成分基本保持不變。當(dāng)存在軸向偏載時，浮動軌跡由圓形分布變?yōu)闄E圓形分布，浮動軌跡中的變形量增加，從而系統(tǒng)振動更加明顯，嚙合沖擊更激烈，影響系統(tǒng)的使用壽命。

行星輪傳動；嚙合剛度；軸向偏載；接觸區(qū)域；動態(tài)特性

行星齒輪傳動具有體積小、重量輕、效率高、工作平穩(wěn)、噪聲小等優(yōu)點，被廣泛應(yīng)用于航空、船艦、風(fēng)力發(fā)電、汽車等許多工業(yè)領(lǐng)域。為提高行星傳動均載性能，太陽輪采用浮動細(xì)桿相連接，相對于其他構(gòu)件的而言，其彎曲剛度較小，則易發(fā)生彎曲情況，并產(chǎn)生較大的振動，從而縮短行星系統(tǒng)的使用壽命。

國內(nèi)外許多學(xué)者針對行星齒輪傳動系統(tǒng)動態(tài)特性問題進行了大量研究，內(nèi)容普遍涉及固有頻率、動態(tài)響應(yīng)，以及浮動軌跡等多個方面。在行星齒輪傳動系統(tǒng)的研究中需要考慮的因素較多，主要有：時變嚙合剛度、傳遞誤差、齒側(cè)間隙等。早期在分析系統(tǒng)動態(tài)特性的研究中，其剛度大多使用平均嚙合剛度，但轉(zhuǎn)動過程為單齒與多齒嚙合的交替形式，存在時變嚙合剛度，因此Lin等[1]研究了行星齒輪由嚙合剛度變化引起系統(tǒng)參數(shù)的不穩(wěn)定性。王春光[2]探討了時變嚙合剛度對系統(tǒng)不穩(wěn)定的影響，并得出調(diào)節(jié)重合度和嚙合相位可以減小特定的不穩(wěn)定區(qū)域。孫濤等[3]考慮時變嚙合剛度、嚙合誤差以及尺側(cè)間隙等因素的基礎(chǔ)上，采用解析諧波平衡法求得系統(tǒng)的非線性頻響特征。齒廓的磨損與加工誤差會使輪副產(chǎn)生齒側(cè)間隙以及齒廓周期誤差，并在齒輪間會產(chǎn)生接觸、脫離、再接觸的重復(fù)沖擊現(xiàn)象，張鎖懷等[4-5]考慮該情況研究了單級齒輪系統(tǒng)的拍擊振動現(xiàn)象。隨后劉國華等[6]學(xué)者對齒輪進行修形處理，有效的改善了齒輪機構(gòu)的動態(tài)特性。為有效引入外部激勵的影響，秦大同等[7]通過改變轉(zhuǎn)矩的方法，進行了風(fēng)力發(fā)電機行星齒輪傳動系統(tǒng)動態(tài)嚙合力的研究。楊軍[8]通過振動學(xué)方程計算了風(fēng)力發(fā)電機行星傳動系統(tǒng)的固有頻率，同時討論了剛度對系統(tǒng)各階固有頻率的影響。其中Li等[9]通過討論齒廓修形、齒面負(fù)載變化等情況下對于齒面接觸區(qū)域分布情況的影響。Lin等[10]研究了當(dāng)行星輪不均勻分布時，行星齒輪傳動系統(tǒng)的固有頻率以及振動模態(tài)的變化情況。Parker等[11]采用有限元/接觸力學(xué)模型分析了直升機行星齒輪系統(tǒng)在一定運行速度和轉(zhuǎn)矩范圍內(nèi)的動態(tài)響應(yīng)，結(jié)合了傳統(tǒng)有限元分析和Bousinesq解法的優(yōu)勢?，F(xiàn)如今國內(nèi)外學(xué)者針對行星系統(tǒng)中太陽輪發(fā)生軸向偏載情況下的研究相對較少。

本文計入了太陽輪軸向偏載的影響，建立了行星齒輪傳動系統(tǒng)動力學(xué)模型，從齒面接觸特性的角度，系統(tǒng)分析軸向偏載對嚙合剛度以及行星傳動系統(tǒng)動態(tài)特性的影響。

1 模型構(gòu)建

本文所分析的行星齒輪傳動系統(tǒng)如圖1(a)所示，其中標(biāo)識R代表齒圈，S代表太陽輪，P代表行星輪，C為行星架。行星齒輪傳動系統(tǒng)由3個行星輪、太陽輪、行星架和內(nèi)齒圈所構(gòu)成。太陽輪為輸入端，并且采用較細(xì)長的彈性細(xì)軸連接，起彈性浮動的作用。齒圈固定，行星架為輸出端[12]。

圖1 行星傳動系統(tǒng)結(jié)構(gòu)簡圖與三維模型Fig.1 The planetary gear transmission system sketch and three-dimensional model

本文分析模型基本參數(shù)如表1所示，其中轉(zhuǎn)動慣量以及質(zhì)量等均是通過UG建模得出，傳動系統(tǒng)三維裝配模型如圖1(b)所示。

表1 行星齒輪傳動系統(tǒng)參數(shù)Tab. 1 The planetary gear transmission system parameters

由于行星齒輪系統(tǒng)的復(fù)雜性，對其進行動力學(xué)研究時必須考慮零件或運動副的彈性。根據(jù)動力學(xué)建立原理的不同，可以將行星齒輪傳動的動力學(xué)模型分為兩大類，即集中參數(shù)模型和有限元模型[13]。

本文選用集中參數(shù)模型，將行星傳動的各構(gòu)件簡化成集中質(zhì)量，各構(gòu)件間以彈簧作為連接形式。將其運動看成剛體運動以及彈性變形的疊加，本次行星傳動系統(tǒng)為16自由度的平移-扭轉(zhuǎn)動力學(xué)模型。模型中各個構(gòu)件的支撐剛度，內(nèi)嚙合、外嚙合剛度以及扭轉(zhuǎn)剛度均用彈簧代替，現(xiàn)取行星架中心位置為系統(tǒng)坐標(biāo)原點，選用水平方向為X方向，豎直方向為Y方向來構(gòu)建模型。行星輪采用沿行星架均勻分布方式，并對其進行約束。

所施加載荷足夠大時，太陽輪與行星輪齒輪齒廓始終保持嚙合關(guān)系，并且兩齒面所產(chǎn)生的嚙合力與齒面彈性變形成線性關(guān)系，則可以把齒輪副剛度用線性彈簧來模擬，太陽輪與行星輪動力學(xué)模型如圖2所示，太陽輪與行星輪的扭轉(zhuǎn)和橫向微位移均會對嚙合作用產(chǎn)生影響，將其微位移轉(zhuǎn)化至嚙合線上。定義當(dāng)嚙合彈簧為壓縮狀態(tài)時取正值，反則拉伸為負(fù)值，則有：

δspi=xScosφspi+ySsinφspi+uS-xpisinα-
ypicosα-upi-espi(t)

(1)

式中：φspi=φc-α+φpi+π/2，

uS=θSrS,upi=θpirpi,φc=wct

其中,α為嚙合角，φpi是相位角 (φi=2π(i-1)/3)。rS,rpi是太陽輪與行星輪基圓半徑,espi(t)是嚙合誤差。

此時，齒輪嚙合力為

(2)

圖2 太陽輪與行星輪動力學(xué)模型Fig.2 The dynamic model of sun gear and planet gear

構(gòu)建行星傳動平移-扭轉(zhuǎn)動力學(xué)模型如圖3所示，模型中Ks，Ksθ分別為太陽輪支撐剛度與扭轉(zhuǎn)剛度；kpi為行星輪支撐剛度，其中i=1,2,3；kspi與krpi分別為太陽輪與行星輪外嚙合剛度及行星輪與齒圈內(nèi)嚙合剛度；內(nèi)齒圈在其4個對稱位置實施約束，取Kr代表內(nèi)齒圈支撐剛度。

圖3 行星傳動平移-扭轉(zhuǎn)動力學(xué)模型Fig.3 The planetary gear transmission system dynamic model

根據(jù)各個零件受力關(guān)系構(gòu)建系統(tǒng)振動微分方程：

將式(3)寫成矩陣形式得

(4)

式中：m為質(zhì)量矩陣；ka(t)為支撐剛度矩陣；km(t)為時變嚙合剛度矩陣；kn(t)為扭轉(zhuǎn)剛度矩陣；T(t)為外部激勵；系統(tǒng)各個構(gòu)件的位移向量q表示為

q=(uc,xs,ys,us,xr,yr,ur,ξpi,ηpi,upi)T

(5)

系統(tǒng)中各支撐剛度與扭轉(zhuǎn)剛度分別為：太陽輪支撐剛度Ks=5×106N/m；太陽輪扭轉(zhuǎn)剛度Ksθ=1.5×105N·m/rad；行星輪支撐剛度Kpi=2×107N/m；行星架扭轉(zhuǎn)剛度Kpθ=3.3×108N·m/rad；內(nèi)齒圈支撐剛度Kr=3×108N/m；

2 參數(shù)計算

2.1 嚙合剛度計算

本文通過有限元法，進行系統(tǒng)內(nèi)外嚙合剛度的求解，同時考慮太陽輪軸向偏載的情況下，太陽輪與行星輪之間外嚙合剛度的變化情況。在求解外嚙合剛度時，采用太陽輪施加轉(zhuǎn)矩，行星輪全約束的形式進行求解，如圖4(a)所示。在求解內(nèi)嚙合時，采用行星輪施加純扭矩，內(nèi)齒圈施加全約束的形式進行求解，如圖4(b)所示。

圖4 齒輪副分析模型Fig.4 Gear pair analysis model

在施加邊界條件時，對從動輪內(nèi)圈施加全約束，并且約束主動輪內(nèi)圈的橫向自由度，施加轉(zhuǎn)矩在主動輪上[14]。在施加轉(zhuǎn)矩時，需要將坐標(biāo)系轉(zhuǎn)換成圓柱坐標(biāo)系，以節(jié)點的形式將轉(zhuǎn)矩施加在主動輪內(nèi)圈處，計算如式(6)所示

(6)

式中：Ta為單個節(jié)點施加轉(zhuǎn)矩；n為節(jié)點數(shù)目；rn為主動輪內(nèi)圈半徑；rb為主動輪基圓半徑。

通過變形云圖計算出主動輪內(nèi)圈變形量μ1，利用式(7)求得主動輪變形角度δ1，最終求得嚙合線位置處的等效變形μ：

(7)

式中：δ1為主動輪變形角度；rn為主動輪內(nèi)圈半徑；μ1為主動輪內(nèi)圈變形；μ為嚙合線處等效變形。

求解時，選用罰函數(shù)法進行求解，該方法允許接觸目標(biāo)存在一定量的滲透，并可通過改變接觸剛度值來控制表面滲透的程度[15]。最終使用剛度計算式(8)進行剛度求解

(8)

通過上述計算方法，我們依次算出一個嚙合周期內(nèi)，10種嚙合狀態(tài)下的齒輪嚙合剛度。并通過數(shù)值計算的插值法，擬合成一個周期內(nèi)的曲線，進行求解內(nèi)外嚙合時剛度的變化情況，并繪制曲線圖，如圖5(a)所示。

圖5 嚙合剛度變化圖以及應(yīng)力分布圖Fig.5 Meshing stiffness variation and stress distribution diagram

從圖5中可以清楚的看出，嚙合剛度為時變嚙合剛度，為拱形分布。雙齒嚙合剛度大于單齒嚙合剛度，雙齒嚙合狀態(tài)下，輪齒嚙合變形減小，從而剛度增大。齒輪的嚙合周期為雙齒-單齒-雙齒嚙合，最后退出嚙合的過程，為時變嚙合剛度，因此存在剛度的突變，容易產(chǎn)生較大剛度激勵。

現(xiàn)將接觸位置的輪齒輪廓曲線長度延展后的長度定義為輪廓曲線長度，為圖5(b)中的縱坐標(biāo)。從圖5(b)中可以清楚的看出接觸區(qū)域應(yīng)力成帶狀分布，并且越接近理論接觸位置等效應(yīng)力值越大。由于齒面摩擦的原因，接觸線沿齒高方向的兩側(cè)并未呈現(xiàn)對稱性的情況，接觸線上方向著齒頂?shù)姆较虻刃?yīng)力變化較為劇烈，反之，接觸線下方向著齒根的方向等效應(yīng)力變化較為平緩。接觸線區(qū)域為主要的承載區(qū)域，該區(qū)域載荷分布較大，并且分布較為均勻，但由于存在邊緣效應(yīng)的原因，齒輪接觸區(qū)域的兩端存在較大的應(yīng)力集中，接觸范圍變化最明顯。本次所繪的嚙合位置應(yīng)力分布圖與文獻9中基本相同，但本文所繪圖在齒端兩側(cè)出現(xiàn)了明顯邊緣效應(yīng)區(qū)域。

2.2 系統(tǒng)的固有特性分析

當(dāng)系統(tǒng)外部激勵T(t)=0時，式(4)可得自由振動方程為

(9)

將系統(tǒng)固有特性變化為求解特征值問題，即：

(10)

(11)

將計算的剛度結(jié)果帶入行星傳動平移-扭轉(zhuǎn)動力學(xué)模型中，通過求解即可得到系統(tǒng)固有頻率，如表2所示，可以看出該行星傳動系統(tǒng)的固有頻率范圍420～16 703 Hz，共16階。由于行星傳動動力學(xué)模型的對稱性，會引起系統(tǒng)固有頻率產(chǎn)生重頻現(xiàn)象。本文分析對象中心構(gòu)件平移振動模式出現(xiàn)6對重頻(分別為ω1=ω2、ω3=ω4、ω6=ω7、ω12=ω13、ω14=ω15和ω8=ω9)與系統(tǒng)6構(gòu)件橫向自由度數(shù)相等；行星輪振動模式共有3階固有頻率等于系統(tǒng)3個行星輪扭轉(zhuǎn)自由度數(shù)[16]。所得結(jié)論與文獻[16]基本相同。

表2 固有頻率Tab.2 The natural frequency

行星傳動系統(tǒng)振型分為3種，分別為中心構(gòu)件平移振動模式、中心構(gòu)建扭轉(zhuǎn)振動模式、以及行星輪振動模式，其中中心構(gòu)件包括行星架，太陽輪，內(nèi)齒圈。如圖6所示。

圖6 行星傳動系統(tǒng)振型Fig.6 Planetary transmission system vibration mode

3 偏載對行星系統(tǒng)的影響

3.1 軸向偏載對于剛度的影響

對于主要承受彎曲的梁和軸，撓度和轉(zhuǎn)角過大會影響構(gòu)件或零件的正常工作。齒輪軸的撓度過大會影響齒輪的嚙合，或增加齒輪的磨損并產(chǎn)生噪聲，同時當(dāng)軸承支撐的軸在支撐處的轉(zhuǎn)角如果過大時，會增加軸承的磨損。

連接太陽輪的軸存在較大彈性浮動，在工作狀態(tài)下，易發(fā)生彎曲，從而受載不均勻，產(chǎn)生彎曲，形成軸向偏載的情況，最終導(dǎo)致齒面載荷分布不均勻。現(xiàn)定義γ為軸偏載位移量，如圖7所示本文討論γ為正常情況下，1μ，3μ的情況下，齒輪剛度的變化曲線如圖8所示。

圖7 軸向偏載示意圖Fig.7 Diagrammatic sketch when the axis bending

圖8 軸向偏載時剛度的變化曲線Fig.8 Stiffness variation curve when the axis bending

當(dāng)軸發(fā)生軸向偏載時，對齒輪剛度的分布狀況是有影響的，隨著偏載的增大，齒輪嚙合剛度逐漸降低。從圖8中可以清楚的看出雙齒嚙合區(qū)域降低幅度遠(yuǎn)大于單齒嚙合區(qū)域。剛度隨著軸向偏載的變化原因是由于齒輪在嚙入嚙出的過程中發(fā)生了齒頂?shù)母缮媲闆r，造成齒頂處等效應(yīng)力較為集中，從而產(chǎn)生相應(yīng)的影響。

從圖9中可以清楚的看出接觸區(qū)域的改變，結(jié)合剛度隨軸向偏載變化曲線如圖8所示。如圖9，可以看出隨著軸向偏載的增大，最大接觸區(qū)域逐漸向軸向偏載位置移動，局部位置呈現(xiàn)較大的應(yīng)力分布情況，并呈現(xiàn)不規(guī)則分布。如圖9(b)，接觸線兩側(cè)的應(yīng)力分布范圍同樣存在較大改變，接觸線上方沿著齒頂方向接觸應(yīng)力變化更加劇烈，同時，接觸線下方沿著齒根方向接觸應(yīng)力變化也較為劇烈。綜上所述，當(dāng)軸向偏載發(fā)生時，會導(dǎo)致接觸區(qū)域的應(yīng)力分布變化，齒輪產(chǎn)生變形，進而影響齒輪剛度的變化。

圖9 嚙合區(qū)域應(yīng)力分布圖Fig.9 Stress distribution of meshing area

3.2 軸向偏載對固有頻率的影響

現(xiàn)探討軸向偏載對行星齒輪傳動系統(tǒng)的固有頻率的影響，當(dāng)連接太陽輪的彈性細(xì)桿發(fā)生彎曲時，行星系統(tǒng)的結(jié)構(gòu)將發(fā)生改變，轉(zhuǎn)動軌跡近似于橢圓形的方式轉(zhuǎn)動，故而探討偏載狀態(tài)下，固有頻率的變化。通過改變偏載量分別為1μ，3μ時，觀察固有頻率的變化規(guī)律。

如圖10所示，通過觀察可以看出，在1～9階時，其固有頻率基本保持不變，由此可得結(jié)論：軸向偏載對行星傳動系統(tǒng)的低階固有頻率影響較小。但10階、12階、13階、14階、15階、16階情況下，其行星系統(tǒng)的固有頻率隨著偏載的增大，在逐漸減小，因此軸向偏載對行星系統(tǒng)的高階固有頻率的影響較明顯。隨軸向偏載的增大，系統(tǒng)低階固有頻率基本不變，高階固有頻率減小，但各階振型基本不變。由于本文只是改變了其連接太陽輪的細(xì)桿的軸向偏載，故其主要影響高階扭轉(zhuǎn)振形的固有頻率，以及與太陽輪嚙合的行星輪平移振動模式的固有頻率。

圖10 偏載對系統(tǒng)固有頻率變化圖Fig.10 System natural frequency change when the axis bending

3.3 軸向偏載對動態(tài)特性的變化

當(dāng)太陽輪發(fā)生軸向偏載時，行星系統(tǒng)的動態(tài)特性將發(fā)生改變，由于偏載量較小，動態(tài)特性的變化情況并不明顯，為了能充分研究偏載對系統(tǒng)動態(tài)特性的影響，將偏載量進行增大。偏載量為1μ，3μ，5μ，10μ，15μ，20μ，研究動態(tài)特性的變化。

在一定負(fù)載作用下，保證齒輪副始終處于接觸狀態(tài)，通過式(12)計算嚙合頻率，其轉(zhuǎn)速為700 r/min，轉(zhuǎn)矩為500 N·m時計算行星傳動系統(tǒng)的動載荷，如圖11所示。

(12)

式中：nc為行星架轉(zhuǎn)速；ns為太陽輪轉(zhuǎn)速；np為行星架轉(zhuǎn)速；zs為太陽輪齒數(shù)；zp為行星輪齒數(shù)；fp為嚙合頻率。

圖11 正常情況下太陽輪與行星輪動載荷 Fig.11 The dynamic force of sun gear and planet gear in the normal condition

太陽輪與行星輪外嚙合時域歷程如圖11(a)所示，可以看出，外嚙合的動載荷成周期變化，周期為T，其動載荷均值為778 N，波動幅值為194 N；太陽輪與行星輪的外嚙合動載荷的主要頻率成分為嚙合頻率(292.6 Hz)，2倍頻，3倍頻及4倍頻，其中嚙合頻率成分能量最大。

當(dāng)太陽輪軸向偏載為20μ時，通過計算行星輪繞太陽輪公轉(zhuǎn)一周的各個位置處的時變嚙合剛度，引入至系統(tǒng)動力學(xué)模型，采用時域積分方法求解模型，即可得到系統(tǒng)動載荷，如圖12(a)所示，由于轉(zhuǎn)速保持700 r/min不變，則系統(tǒng)的嚙合頻率保持一致為292.6 Hz，并且嚙合頻率的成分能量最大。由于軸向偏載的緣故，外嚙合的動載荷呈周期變化，其曲線圖為波動情況，近似于正弦曲線圖，其動載荷均值為783 N，波動幅值高達(dá)230 N。此時各頻率成分的幅值代表該時域狀態(tài)下所具有的能量。由于太陽輪發(fā)生偏載后，太陽輪類似于凸輪的形式進行轉(zhuǎn)動，其輪齒嚙合區(qū)域發(fā)生了明顯變化，類似于偏心輪接觸的方式，則其動載荷也變成正弦的波動形式。行星系統(tǒng)正常情況下太陽輪無轉(zhuǎn)頻成分，但偏載為20μ時，太陽輪出現(xiàn)轉(zhuǎn)頻成分，其幅值為21 N，即該時域狀態(tài)下轉(zhuǎn)頻所具有的能量。由此可見軸向偏載的增大，改變了太陽輪的頻率成分，出現(xiàn)了轉(zhuǎn)頻成分。由于軸向偏載情況的存在，太陽輪發(fā)生不平衡或傾斜的現(xiàn)象，當(dāng)太陽輪存在軸向偏載時，太陽輪與其相嚙合的3個行星輪的接觸區(qū)域均發(fā)生了改變，并且接觸區(qū)域均向不同方向發(fā)生了偏移，由于接觸區(qū)域的變化，則太陽輪將存在轉(zhuǎn)頻成分。

通過研究偏載狀態(tài)對嚙合頻率幅值以及轉(zhuǎn)頻幅值影響如圖13所示，由于轉(zhuǎn)速固定為700 r/min的情況下，觀察各種偏載情況下的嚙合頻率幅值并未存在明顯的變化，同時，其各倍頻也無明顯變化，其幅值基本穩(wěn)定，所以偏載的存在對其嚙合頻率的影響并不顯著。通過觀察太陽輪的轉(zhuǎn)頻曲線，正常狀態(tài)下的動載荷頻域圖，太陽輪無轉(zhuǎn)頻成分，但隨著偏載的增大，其轉(zhuǎn)頻能量明顯增大，當(dāng)偏載量增大到20μ的時候，其轉(zhuǎn)頻的幅值為21 N，太陽輪轉(zhuǎn)頻幅值隨偏載增大而逐漸增大，呈線性遞增的方式，轉(zhuǎn)頻成分的能量明顯增大。由于嚙合齒輪在轉(zhuǎn)動過程中，齒頂位置存在間隙，當(dāng)發(fā)生軸向偏載時，對齒輪嚙合頻率等相關(guān)成分影響較小。

圖12 偏載情況下太陽輪與行星輪動載荷 Fig.12 The dynamic force of sun gear and planet gear when the axis bending

圖13偏載對于嚙合頻率幅值以及轉(zhuǎn)頻幅值變化圖 Fig.13 The change of meshing frequency amplitude and rotating frequency amplitude

3.4 軸向偏載對中心輪浮動軌跡影響

通過研究偏載狀態(tài)下，太陽輪的浮動軌跡變化如圖14所示，正常情況下太陽輪浮動軌跡，當(dāng)轉(zhuǎn)速為700 r/min時，太陽輪浮動軌跡近似為圓形，浮動范圍為-10.9 μm～10.9 μm，浮動軌跡的內(nèi)圈為太陽輪所運動的軌跡，正常情況下，載荷施加均勻并無誤差存在時，太陽輪的運動軌跡為圓形。由于太陽輪的扭轉(zhuǎn)固有頻率作用，太陽輪產(chǎn)生了振動。在不同偏載情況下，其太陽輪的浮動軌跡也發(fā)生了相應(yīng)變化。觀察下圖14(a)中太陽輪浮動軌跡可得，隨著偏載的增加，其X方向以及Y方向的位移范圍也在相應(yīng)增大，當(dāng)偏載達(dá)到5μ時，其X方向上的位移量為-12 μm～12 μm，Y方向上的位移浮動范圍為-11.2 μm～11.2 μm，近似為橢圓形分布。通過計算X方向上的變形量為6.3 μm，Y方向上的變形量為4.6 μm。當(dāng)偏載達(dá)到20μ時，浮動軌跡呈橢圓形分布。其X方向上的變形量為21 μm，Y方向上的變形量為19.5 μm。軸向偏載越大，其浮動軌跡圖呈橢圓趨勢越明顯，其偏載方向X的浮動位移量相對較大，并且由于X方向為主要傾斜方向，所以X方向的變形比較大。Y方向浮動位移量相對較小，且Y方向上的變形量略小于X方向，其主要變形在X方向上。通過研究表明，當(dāng)軸向偏載存在時，其太陽輪的浮動軌跡發(fā)生了相應(yīng)變化，隨著軸向偏載的增大，太陽輪的浮動量越大，太陽輪的變形也在增大，其運動軌跡也逐漸趨近于橢圓形。由于偏載的情況，偏載越大，行星系統(tǒng)的載荷分布越不均勻，需要越大的浮動量來平衡載荷。因此偏載方向處的浮動位移量較大，也同樣會產(chǎn)生較大的沖擊力，進而增大該方向的變形量。

圖14 各偏載下中心輪與內(nèi)齒圈的浮動軌跡Fig.14 Floating track of the center wheel and the ring gear

圖14(b)所示，當(dāng)內(nèi)齒圈處于正常狀態(tài)下時，其內(nèi)齒圈浮動軌跡近似于圓形分布，X、Y方向上的變形量均為0.36 μm，偏載越大，內(nèi)齒圈的浮動軌跡同樣越趨近于橢圓形分布，偏載為20μ時，其浮動軌跡成發(fā)散狀。通過研究表明，偏載的增大，內(nèi)齒圈的變形量也在增大，但增加幅度較小，則偏載對內(nèi)齒圈的變形量變化影響并不顯著。主要影響其浮動軌跡的形狀，由于內(nèi)齒圈的浮動能夠完成構(gòu)件誤差對內(nèi)嚙合輪齒的法向位移要求。內(nèi)齒圈固定在支架上，其支撐剛度較大，則內(nèi)齒圈的浮動軌跡較不明顯，內(nèi)齒圈的變形量也相對較小。

4 結(jié) 論

(1)軸向偏載會導(dǎo)致齒輪在嚙合過程中嚙合剛度的降低，單齒嚙合區(qū)域中的嚙合剛度降低幅度比雙齒嚙合區(qū)域中嚙合剛度降低幅度小，接觸位置也發(fā)生偏移。

(2)軸向偏載變化時，對系統(tǒng)高階固有頻率影響較為顯著，低階固有頻率基本穩(wěn)定。隨著軸向偏載的增大，其高階固有頻率降低。

(3)對系統(tǒng)動態(tài)特性而言，太陽輪與行星輪的外嚙合動載荷將出現(xiàn)轉(zhuǎn)頻成分，并且當(dāng)偏載增大時，其轉(zhuǎn)頻成分將明顯增大，其他頻率成分基本保持不變。

(4)通過研究發(fā)現(xiàn)，當(dāng)隨著軸向偏載的改變，其太陽輪的浮動軌跡也發(fā)生了相應(yīng)變化，由于軸向偏載的增大，太陽輪浮動軌跡變大，其變形量變大。偏載的改變，對內(nèi)齒圈的變形量影響較小，只影響其浮動的軌跡。

[1] LIN J，PARKER R G. Planetary gear parametric instability caused by mesh stiffness variation[J]. Journal of Sound and Vibration, 2002,249(1): 129-145。

[2] 王春光. 行星齒輪傳動動態(tài)特性的研究[D].北京：中國艦船研究院,2005.

[3] 孫濤,沈允文,孫智民,等. 行星齒輪傳動非線性動力學(xué)方程求解與動態(tài)特性分析[J]. 機械工程學(xué)報,2002,38(3):11-15.

SUN Tao，SHEN Yunwen，SUN Zhimin, et al. Nonlinear dynamic equation and dynamic characteristics analysis[J]．Journal of Mechanical Engineering，2002，38(3) : 11-15．

[4] 張鎖懷,沈允文,董海軍,等. 單級齒輪系統(tǒng)拍擊特性研究[J]. 機械工程學(xué)報,2003(3): 28-31.

ZHANG Suohuai, SHEN Yunwen, DONG Haijun, et al. Research on rattling dynamics of one-stage gear system[J]. Journal of Mechanical Engineering, 2003(3): 28-31.

[5] 張鎖懷,沈允文,董海軍,等. 齒輪拍擊系統(tǒng)的動力響應(yīng)[J]. 振動工程學(xué)報,2003(1): 66-70.

ZHANG Suohuai, SHEN Yunwen, DONG Haijun, et al. Dynamic response of a gear rattling system[J]. Journal of Vibration Engineering, 2003(1): 66-70.

[6] 劉國華,李亮玉. 含間隙和時變嚙合剛度的漸開線齒輪副齒廓修形研究[J]. 中國機械工程,2007(12): 1409-1413.

LIU Guohua, LI Liangyu. Research on conjugating tooth profile of involute gear transmission with backlash and time—varying stiffness[J]. China Mechanical Engineering, 2007(12): 1409-1413.

[7] 秦大同,田苗苗,楊軍. 變風(fēng)載下風(fēng)力發(fā)電機齒輪傳動系統(tǒng)動力學(xué)特性研究[J]. 太陽能學(xué)報,2012(2):190-196.

QIN Datong, TIAN Miaomiao, YANG Jun. Study on dynamic characteristics of gear transmission system of wind generator under varying wind load [J]. Acta Energiae Solaris Sinica,2012(2): 190-196.

[8] 楊軍. 風(fēng)力發(fā)電機行星齒輪傳動系統(tǒng)變載荷激勵動力學(xué)特性研究[D].重慶：重慶大學(xué),2012.

[9] LI Shuting. Effects of misalignment error,tooth modifications and transmitted torque on tooth engagements of a pair of spur gears[J].Mechanism and Mechine Theory, 2015, 83: 125-136.

[10] LIN J, PARK R G. Structured vibration characteristics of planetary gears with unequally spaced planets[J]. Journal of Sound and Vibration,2000,233(5): 921-928.

[11] PARKER R G, AGASHE V, VIJAYAKAR S M. Dynamic response of a planetary gear system using a finite element/contact mechanics mode[J]. Journal of Mechanical Design, 2000, 122: 304-310.

[12] ZHOU Jianxing, SUN Wenlei, YUAN Liang. Nonlinear vibroimpact characteristics of a planetary gear transmission system[J]. Shock and Vibration, 2016(4): 1-11.

[13] 李潤方，王建軍. 齒輪系統(tǒng)動力學(xué)——振動、沖擊、噪聲[M].北京: 科學(xué)出版社, 1997: 380-405.

[14] 周建星,孫文磊,曹莉,等. 行星齒輪傳動系統(tǒng)碰撞振動特性研究[J]. 西安交通大學(xué)學(xué)報,2016(3): 16-21.

ZHOU Jianxing，SUN Wenlei，CAO Li，et al. Vibro-impact characteristics of planetary gear transmission[J]. Journal of Xi’an Jiaotong University,2016(3): 16-21.

[15] 卜忠紅,劉更,吳立言. 行星齒輪傳動動力學(xué)研究進展[J]. 振動與沖擊,2010,29(9): 161-166.

BU Zhonghong, LIU Geng, WU Liyan. Research advances in planetary gear trains dynamics [J]. Journal of Vibration and Shock,2010,29(9): 161-166.

[16] 許華超,孫文磊,周建星,等. 風(fēng)電機組行星傳動系統(tǒng)固有特性靈敏度研究[J]. 太陽能學(xué)報,2016(1): 201-207.

XU Huachao，SUN Wenlei，ZHOU Jianxing，et al. Sensitivity of planetary natural characteristics used in wind turbine [J]. Acta Energiae Solaris Sinica，2016(1): 201-207.

[17] 邱星輝,韓勤鍇,褚福磊. 風(fēng)力機行星齒輪傳動系統(tǒng)動力學(xué)研究綜述[J]. 機械工程學(xué)報,2014(11): 23-36.

QIU Xinghui，HAN Qinkai，CHU Fulei. Review on dynamic analysis of wind turbine geared transmission systems [J]. Journal of Mechanical Engineering, 2014(11): 23-36.

[18] WU S, ZUO M J, PAREY A. Simulation of spur gear dynamics and estimation of fault growth[J]. Journal of Sound and Vibration 2008,317 (3/4/5): 608-624.

[19] PEETERS J L M, VANDEPITTE D, SAS P. Analysis of internal drive train dynamics in a wind turbine[J]. Wind Energy, 2006(9): 141-164.

[20] AMBARISHA V K, PARKER R G. Nonlinear dynamics of planetary gears using analytical and finite element models[J]. Journal of Sound and Vibration, 2007, 302: 577-595.

[21] 石萬凱,劉敬,龔建春. 同軸對轉(zhuǎn)行星齒輪傳動系統(tǒng)動態(tài)特性分析[J]. 振動與沖擊,2014,33(7): 99-104.

SHI Wankai, LIU Jing, GONG Jianchun. Dynamic characteristics analysis for a co-axial counter-rotating planetary gear transmission system[J].Journal of Vibration and Shock, 2014,33(7): 99-104.

[22] 劉輝,蔡仲昌,項昌樂. 兩級行星齒輪傳動非線性嚙合力頻率耦合與動態(tài)特性研究[J]. 振動與沖擊,2015,34(19): 13-23.

LIU Hui, CAI Zhongchang, XIANG Changle. Frequency coupling and dynamic characteristics pf nonlinear meshing force for two-stage planetary gears[J]. Journal of Vibration and Shock, 2015,34(19): 13-23.

Astudyonthedynamiccharacteristicsofaplanetarygearsystemwithconsideringcontacttoothsurface

CHENRuibo,ZHANGJianjie,ZHOUJianxing,SUNWenlei

(SchoolofMechanicalEngineering,XinjiangUniversity,Urumqi830047,China)

Axial partial load is an important factor influencing the stability of a planet wheel drive system. Based on the planet wheel drive system and the Newtonian dynamics theory, a planetary transmission dynamics model was established. The effect of axial load for meshing stiffness was studied. At the same time, considering the influence of the tooth surface contact characteristics of mesh stiffness, the axial partial load impact on the natural frequency of the planetary system, the influence of the axial partial load on the system dynamic characteristics was studied. The research shows that: the sun wheel connects thin rod axial after partial load, the meshing stiffness decreases with the increase of partial load. Axial deviation of the existence of a planet wheel has obviously influence on high order natural frequency. When the partial load increases, the high order natural frequency decreases. With the increase of axial partial load round the sun and the planets round of external gear dynamic load changes into a sine curve, frequency domain figure appeared switch frequency. And partial load increases, the frequency amplitude increases obviously. Meshing frequency amplitude is basically stable, does not affect the frequency doubling amplitude.

planetary gear transmission; meshing stiffness; axial eccentric load; contact area; dynamic characteristics

國家自然基金(51665054；51565055)

2016-10-10 修改稿收到日期： 2016-11-17

陳銳博 男，碩士生， 1993年生

張建杰 男，副教授，碩士生導(dǎo)師，1973年生

TH113

A

10.13465/j.cnki.jvs.2017.20.028