基于彈性圓柱殼的柔性齒輪副建模及動力學(xué)分析

張 洋，關(guān)先磊，王青山

（中南大學(xué) 高性能復(fù)雜制造國家重點(diǎn)實(shí)驗(yàn)室，長沙 410083）

齒輪傳動系統(tǒng)廣泛應(yīng)用于汽車、飛機(jī)等各種場合，被認(rèn)為是最重要的機(jī)械部件之一。在工作過程中，傳動齒輪存在橫向與面內(nèi)的彈性振動，由此會引起結(jié)構(gòu)失效等動力學(xué)現(xiàn)象。這種結(jié)構(gòu)失效現(xiàn)象在高速輕質(zhì)重載齒輪中經(jīng)常出現(xiàn)，例如高速動車牽引齒輪傳動系統(tǒng)以及航空發(fā)動機(jī)齒輪系統(tǒng)。因此，建立考慮輪體柔性的齒輪傳動動力學(xué)模型具有十分重要的意義。

為了獲得高效齒輪，人們對齒輪系統(tǒng)的動態(tài)建模和振動特性分析進(jìn)行了大量的研究。對于齒輪的輪體結(jié)構(gòu)，通常有三種建模方法。其一是扭振模型即剛性模型，Han和Chu[1]研究了齒輪傳動轉(zhuǎn)子系統(tǒng)周期角角速度運(yùn)動下的動力學(xué)特性。Chen等[2]研究了雙螺旋齒輪的旋轉(zhuǎn)動力學(xué)，考慮到了軸承、回轉(zhuǎn)效果，其中采用鐵木辛柯梁作為軸單元，齒輪則視為剛性體。He 等[3]提出了一種新穎的剛性齒輪嚙合模型，并采用改進(jìn)的勢能法計(jì)算了齒輪偏心對時效嚙合剛度的影響。Shi等[4]完成了考慮時變參數(shù)和側(cè)隙的嚙合圓柱齒輪對的建模和分析。其二是以彈性環(huán)作為齒輪的輪體模型，Christopher和Robert[5]研究以耦合彈性環(huán)為模型，分析研究高速柔性齒輪副的振動。Chen等[6]研究了帶撓性正齒環(huán)行星齒輪組的動態(tài)仿真。Fan 等[7]根據(jù)殼理論和鐵木辛柯梁理論建立了行星齒輪內(nèi)環(huán)和轉(zhuǎn)軸模型。其三是以柔性圓盤作為齒輪輪體模型。Vinayak 和Singh[8]建立了柔性齒輪的多體動力學(xué)模型，考慮了分布齒輪嚙合剛度和齒胚模型，求解了一系列的模態(tài)問題和響應(yīng)問題。

綜合比較上述幾種模型可知，剛性模型具有模型簡單，計(jì)算量小的優(yōu)點(diǎn)，但是剛性模型沒有考慮到齒輪齒面變形，和實(shí)際齒輪貼合度存在差距，尤其是齒輪輪體厚度較薄的齒輪，剛性模型和柔性模型的差距較大。相比較于剛性模型，柔性輪體模型則更加符合工程實(shí)際，因此柔性模型的建模和分析是非常有必要的?；趶椥原h(huán)的柔性齒輪副模型只適用于齒輪的齒寬較小的情況，而以柔性圓盤為齒輪輪體的齒輪模型雖然適用范圍較廣，但是存在模型收斂性不好的缺點(diǎn)?；诂F(xiàn)有柔性齒輪副模型研究的不足，本文基于彈性圓柱殼單元，建立一種新穎的柔性齒輪輪體模型，相比于上述兩種柔性模型，基于彈性圓柱殼的齒輪副模型，具有良好的收斂性，可適用于齒輪寬度較大的情況，還可以進(jìn)行齒輪寬度方向的動力學(xué)特性研究。

綜上所述，本文首先基于微分求積有限元法和旋轉(zhuǎn)圓柱殼理論，推導(dǎo)旋轉(zhuǎn)圓柱殼的動能和勢能，完成柔性齒輪的輪體建模。兩圓柱齒輪由空間固定的具有時變嚙合剛度的彈簧耦合，表示齒輪之間的嚙合效應(yīng)。并且每個圓柱殼用彈性基礎(chǔ)支承，來表示軸以及軸承的彈性。最后分別使用拉格朗日方程和Newmark 迭代來求解柔性齒輪副的模態(tài)和動態(tài)響應(yīng)。本文以某航空航天齒輪副為例，通過收斂性分析，研究了本文模型的數(shù)值穩(wěn)定性；然后與ANSYS的有限元模型進(jìn)行模態(tài)特性的對比，驗(yàn)證本文模型的正確性。接下來，數(shù)值計(jì)算了大轉(zhuǎn)速范圍下柔性齒輪副系統(tǒng)固有頻率隨轉(zhuǎn)速的變化規(guī)律。最后通過與傳統(tǒng)的基于齒輪體剛性假設(shè)的扭振模型進(jìn)行對比，闡述齒輪柔性對系統(tǒng)動力學(xué)特性帶來的影響。

1 理論推導(dǎo)

1.1 柔性齒輪副的幾何模型

柔性齒輪副的模型如圖1所示。齒輪簡化為圓柱殼模型，齒輪嚙合副簡化為彈簧模型。分別在兩個圓柱殼的中心層上定義正交圓柱坐標(biāo)系來描述殼體的振動，坐標(biāo)軸xi，θi，zi(i=1，2)分別沿著軸向，周向以及徑向。以該坐標(biāo)系為參考，殼體的位移用ui，vi，wi(i=1，2)來表示。齒輪的轉(zhuǎn)速為Ωi，根據(jù)齒輪嚙合原理可知Ω2=-R1Ω1/R2，齒輪的基圓半徑為Ri，殼厚為hi，齒寬為li。假設(shè)齒輪的密度、泊松比、彈性模量以及剪切模量分別用ρi，μi，Ei，Gi表示。每個齒輪分別采用徑向彈簧（kri）和切向彈簧（kθi）來代表軸的彈性，嚙合彈簧的剛度為km。

圖1 柔性齒輪副模型

1.2 齒輪副系統(tǒng)的能量分析

根據(jù)1 階剪切變形理論，殼體的位移分量可由中間面的位移來表示[9]：

殼體上任意點(diǎn)的位置矢量和速度矢量分別為：

其中：i、j、k分別表示x、θ、z方向上的單位矢量。

旋轉(zhuǎn)圓柱殼的動能為：

將式（2）代入式（3）可得：

圓柱殼的彈性勢能公式為：

其中應(yīng)力應(yīng)變?yōu)椋?/p>

其中：

E為彈性模量，μ為泊松比。

兩齒輪對應(yīng)的嚙合點(diǎn)之間用彈簧連接，其彈簧的嚙合勢能為：

其中：v1和v2分別表示主、從動輪上嚙合點(diǎn)位置的切向位移，km(t)為時變嚙合彈簧的剛度。

在圓柱殼中間環(huán)內(nèi)部，分別采用徑向彈簧和切向彈簧模擬軸的彈性。支撐彈簧勢能為：

kri和kθi分別表示徑向和切向單位長度上的彈簧剛度。但是在有限元方法中，彈簧只能添加在節(jié)點(diǎn)上，因此，采用一下等效方式，得到節(jié)點(diǎn)上的等效彈簧剛度：

其中：C表示圓環(huán)周長，N表示有限元網(wǎng)格的節(jié)點(diǎn)數(shù)，keri表示施加在節(jié)點(diǎn)上的等效彈簧剛度。

因此系統(tǒng)的總應(yīng)變能為：

1.3 系統(tǒng)的離散及求解

將系統(tǒng)的的總應(yīng)變能和總動能代入拉格朗日方程可得齒輪嚙合系統(tǒng)的動力學(xué)方程為：

即可得到圓柱殼齒輪嚙合結(jié)構(gòu)的整體質(zhì)量矩陣M、剛度矩陣K、阻尼矩陣C分別為：

其中：

其中：Mk是對稱矩陣。A和B為微分權(quán)系數(shù)矩陣；C為積分權(quán)系數(shù)矩陣[10]，E為單位矩陣。

最后使用Newmark 積分迭代方法進(jìn)行動態(tài)響應(yīng)求解[11]。Newmark積分迭代方法的算法流程圖如圖2所示。

2 算例分析

本文采用實(shí)際工程中的齒輪參數(shù)作為算例對象，主從動齒輪的尺寸參數(shù)如表1所示。材料參數(shù)為：彈性模量E=2.1×1011Pa；密度ρ=7 850 kg/m3；泊松比μ=0.3。需要特別說明支撐彈簧設(shè)置和載荷設(shè)置。在研究系統(tǒng)模態(tài)時，同時考慮徑向和周向的支撐彈簧，模擬轉(zhuǎn)軸和載荷作用下的工況。而在進(jìn)行動力響應(yīng)計(jì)算時，只在主動輪上施加周向載荷（即扭矩），同時去掉周向支撐彈簧，保留徑向支撐彈簧；而從動輪上則不施加載荷，繼續(xù)保留周向和徑向支撐彈簧。

表1 齒輪副尺寸參數(shù)

根據(jù)參考文獻(xiàn)得到的齒輪嚙合剛度的解析模型[12]，可得到如圖3所示的嚙合剛度曲線，取算術(shù)平均值可得平均嚙合剛度1.350 3×109N/m，這將作為后續(xù)進(jìn)行模態(tài)特性分析的嚙合剛度。齒輪的嚙合剛度平均分配到對應(yīng)的有限元嚙合節(jié)點(diǎn)上。

圖3 時變嚙合剛度曲線

2.1 收斂性和模型驗(yàn)證

首先進(jìn)行模態(tài)收斂性分析，通過計(jì)算主從動齒輪嚙合模型的模態(tài)，來驗(yàn)證微分求積有限元法的收斂性。分別沿兩圓柱殼的軸向和周向劃分單元數(shù)為M和N，并且每個單元的微分求積節(jié)點(diǎn)數(shù)為n。由于計(jì)算量過大，先確定軸向單元數(shù)M=2；研究周向單元數(shù)和節(jié)點(diǎn)數(shù)對收斂性的影響。圖中偏差計(jì)算公式為(fi-fend)/fend%，其中fi為某單元數(shù)N和單元節(jié)點(diǎn)數(shù)n時的頻率，fend為N和n同時取最大值時得到的頻率，由圖4易得，模型的收斂性極好，根據(jù)收斂性結(jié)果，在后續(xù)的分析中選擇單個圓柱殼軸向單元數(shù)量為2，周向單元數(shù)量為2×6，單元節(jié)點(diǎn)為6進(jìn)行計(jì)算。

圖4 柔性嚙合副模型的前4階頻率收斂性

為進(jìn)一步驗(yàn)證模型結(jié)果的正確性，接下來進(jìn)行模態(tài)對比驗(yàn)證。在ANSYS中建立有限元模型，圓柱殼采用Shell單元，每個圓柱殼40×80個單元；嚙合彈簧和支撐彈簧使用COMBIN14 單元有限元模型如圖5所示。

圖5 ANSYS模型

表2展示了本模型和ANSYS 模型的前12 階自然頻率對比，可以發(fā)現(xiàn)，兩者結(jié)果吻合很好，圖6展示了齒輪嚙合模型的前4 階模態(tài)振型，其結(jié)果也與ANSYS 結(jié)果吻合。上述結(jié)論充分驗(yàn)證了模型的正確性。

圖6 柔性齒輪副前4階振型結(jié)果對比

表2 兩種模型自然頻率對比

2.2 坎貝爾圖分析

在進(jìn)行動態(tài)響應(yīng)分析之前，首先計(jì)算系統(tǒng)的自然頻率并得坎貝爾圖。在嚙合剛度為1.350 3×107N/m的情況下，計(jì)算不同轉(zhuǎn)速時系統(tǒng)的自然頻率，結(jié)果如圖7所示。在零速時，由于嚙合剛度的存在破壞了圓柱殼的軸對稱性，導(dǎo)致固有頻率不同。

圖7 柔性齒輪副坎貝爾圖

當(dāng)速度從零開始增加時，一些固有頻率增加而另一些固有頻率減少。當(dāng)自然頻率減小到0 時，出現(xiàn)臨界轉(zhuǎn)速。由于此模型的尺寸較小，導(dǎo)致在轉(zhuǎn)速增加到3 000 rad/s時，還未出現(xiàn)第一臨界轉(zhuǎn)速。

2.3 時間步長收斂性研究

從1.3的模型求解部分，可以知道進(jìn)行動態(tài)響應(yīng)計(jì)算時，需要選擇合適的時間步長，時間步長Δt的選擇會影響到動態(tài)響應(yīng)結(jié)果的收斂性。時間步長的選取依據(jù)為：所關(guān)注頻段內(nèi)的結(jié)構(gòu)主要貢獻(xiàn)的若干振型的最小周期（最大頻率）的1/10～1/20[11]。設(shè)取第20 階頻率為最大頻率，取其倒數(shù)的1/c2的值作為一種時間步。對應(yīng)的時間步長區(qū)間為：Δt=4.518×10-6～9.036×10-6s；

同時，還需考慮齒輪嚙合的頻率，不同轉(zhuǎn)速對應(yīng)的時間步長為：

需要注意的是，取兩者時間步長的較小值?，F(xiàn)在轉(zhuǎn)速為100 rad/s 的工況下，選取主動輪上的中間嚙合點(diǎn)作為研究對象，c1取200，研究c2取值對收斂性的影響。表3為c2的取值對應(yīng)的時間步長。

表3 時間步長取值

圖8是時間步長收斂性結(jié)果示意圖，其中局部放大圖的橫縱坐標(biāo)單位與整體圖的單位一致。從圖中可以發(fā)現(xiàn)，在合理區(qū)間內(nèi)，不同時間步長對模型的動態(tài)響應(yīng)曲線影響很小，所得響應(yīng)曲線幾乎完全重合。說明在時間步長足夠小的情況下，時間步長的選擇已經(jīng)不影響動態(tài)響應(yīng)分析結(jié)果，從而驗(yàn)證了本模型在計(jì)算動態(tài)響應(yīng)時時間步的收斂性很好。

圖8 動態(tài)響應(yīng)時間步長收斂性

2.4 剛?cè)崮Ｐ蛣討B(tài)特性對比

接下來，通過和傳統(tǒng)的剛性扭振模型[13]進(jìn)行動態(tài)特性結(jié)果對比，通過兩種模型的差異性，闡述輪體柔性對齒輪副系統(tǒng)的影響。兩種模型處于相同的工況之下，即載荷以及支承彈簧，時變嚙合剛度，模型尺寸都相同。圖9為不同轉(zhuǎn)速下剛?cè)崮Ｐ偷姆€(wěn)定之后的動態(tài)響應(yīng)曲線對比，其中柔性模型的響應(yīng)曲線為主動輪上中間嚙合點(diǎn)處的響應(yīng)曲線，剛性模型的響應(yīng)曲線為主動輪基圓半徑上的位移響應(yīng)曲線。

由圖可知，兩種模型的動態(tài)響應(yīng)周期完全一致，這是由于兩者的時變嚙合剛度一致。剛性模型的動態(tài)響應(yīng)曲線波動范圍明顯大于柔性模型，振動平衡位置也大于柔性模型。這是由于柔性模型中圓柱殼模型的柔性引起。在剛性模型中，僅考慮了嚙合彈簧和支撐彈簧的柔性，因此在相同載荷下，外力做功全部轉(zhuǎn)化成彈簧的勢能；而在柔性模型中，外力做功一部分轉(zhuǎn)化成為彈簧的勢能，還有一部分轉(zhuǎn)化為柔性圓柱殼的應(yīng)變能。換而言之，圓柱殼的柔性對于齒輪副系統(tǒng)有減振吸能效果，因此柔性模型中嚙合彈簧的勢能小于剛性模型中嚙合彈簧的勢能，從而表現(xiàn)為剛性結(jié)果的平衡位置和振動波動范圍都大于柔性模型。

圖10展示了不同轉(zhuǎn)速下剛?cè)崮Ｐ偷念l譜成分，圖10(a)為剛性模型，可以發(fā)現(xiàn)圖中頻率成分以一倍嚙合頻率和多倍嚙合頻率成分為主，當(dāng)一倍頻率和多倍頻率的幅值線與一條平行于轉(zhuǎn)速軸的直線相交時，頻譜圖的幅值將會達(dá)到峰值，該條直線對應(yīng)的頻率為該模型下的第1階固有頻率f1（扭轉(zhuǎn)模型的第1、2 階固有頻率分別為1 156.5 Hz 和15 368 Hz），這種現(xiàn)象叫做參數(shù)共振，這是由于嚙合頻率是該模型的激振頻率，當(dāng)激勵頻率接近固有頻率或接近固有頻率的兩倍時，就會發(fā)生這種共振。由于頻譜圖的關(guān)注范圍為0～10 000 Hz，所以只出現(xiàn)了一個共振區(qū)域。除去參數(shù)共振區(qū)域附近，一倍頻和高倍頻對應(yīng)的峰值隨轉(zhuǎn)速增加，總體上也呈現(xiàn)線性增加的趨勢，且一倍頻率的峰值明顯大于多倍頻率的峰值。

圖10(b)為柔性模型，可以發(fā)現(xiàn)圖中頻率成分仍然以一倍嚙合頻率和高倍嚙合頻率成分為主，但是峰值隨轉(zhuǎn)速的變化波動很大。和剛性模型類似，柔性模型的頻譜圖的幅值也會出現(xiàn)參數(shù)共振現(xiàn)象，但是由于柔性模型使用的圓柱殼單元，該模型的第一階固有頻率為1 230.5 Hz。除去峰值附近區(qū)域，一倍頻和高倍頻對應(yīng)的峰值隨轉(zhuǎn)速增加呈現(xiàn)波動上升的趨勢，但是在高轉(zhuǎn)速下，一倍頻率的峰值并不占有明顯優(yōu)勢。

圖10 剛?cè)崮Ｐ偷念l譜圖

從兩種模型頻譜圖的整體幅值量級上看，剛性模型的頻譜幅值大于柔性模型的頻譜幅值。圖9的動態(tài)響應(yīng)曲線也可以很好說明這個現(xiàn)象，在轉(zhuǎn)速為10 rad/s和500 rad/s的轉(zhuǎn)速下，剛性的振動范圍都遠(yuǎn)大于柔性的波動范圍，因此將時域曲線轉(zhuǎn)換為頻譜圖時，就會出現(xiàn)剛性頻譜幅值大于柔性頻譜幅值。

3 結(jié)語

(1)本文基于圓柱殼理論，考慮時變嚙合剛度以及彈性支撐的影響，建立了柔性齒輪副的等效簡化模型。該模型具有適用范圍廣以及收斂性好的優(yōu)點(diǎn)，豐富了現(xiàn)有的理論知識。

(2)基于該模型，開展了柔性齒輪副動力學(xué)特性驗(yàn)證和研究。推導(dǎo)了基于微分求積有限單元法的圓柱殼齒輪嚙合模型的線性微分方程組。在確定軸向單元數(shù)的前提下，分別研究不同周向單元數(shù)以及單元節(jié)點(diǎn)數(shù)下模型的自然頻率，討論了模型的數(shù)值收斂性。接下來，與ANSYS的有限元模型對比自然頻率，來驗(yàn)證本文模型的正確性。

(3)建立了和圓柱殼齒輪嚙合副模型等效的集中質(zhì)量剛性模型。通過兩種模型的動態(tài)響應(yīng)曲線、位移響應(yīng)頻譜圖的對比，來研究齒輪柔性的影響。

(4)在相同的工況下，柔性模型和剛性模型的動態(tài)響應(yīng)曲線周期相同，并且兩種模型在同一轉(zhuǎn)速下的響應(yīng)曲線形狀相似，但是柔性模型結(jié)果的平衡位置和振動波動范圍都小于剛性模型。這是由于柔性模型中圓柱殼模型的柔性引起。在剛性模型中，僅考慮了嚙合彈簧和支撐彈簧的柔性，因此在相同載荷下，外力做功全部轉(zhuǎn)化成彈簧的勢能；而在柔性模型中，外力做功一部分轉(zhuǎn)化成為彈簧的勢能，還有一部分轉(zhuǎn)化為柔性圓柱殼的應(yīng)變能。

(5)從位移響應(yīng)頻譜圖可以看出，兩種模型的頻率成分都以一倍嚙合頻率和高倍嚙合頻率成分為主。但是由于圓柱殼模型的柔性，導(dǎo)致了其固有頻率與剛性模型不一致，從而導(dǎo)致了參數(shù)共振區(qū)域的不一致，進(jìn)而影響了一倍嚙合頻率和高倍嚙合頻率的幅值。并且由于齒輪柔性的影響，柔性模型頻譜圖的幅值整體上小于剛性模型的幅值。