齒輪因其傳動(dòng)優(yōu)勢(shì)而被廣泛應(yīng)用在各類傳動(dòng)設(shè)備中，其動(dòng)力學(xué)特性直接反應(yīng)了設(shè)備運(yùn)行品質(zhì)。然而，由于齒輪內(nèi)外部激勵(lì)沖擊等因素的存在，會(huì)對(duì)傳動(dòng)系統(tǒng)的動(dòng)力學(xué)特性產(chǎn)生影響，加快齒輪磨損或其他故障的發(fā)生，引起劇烈的振動(dòng)和噪聲。因此，研究齒輪傳動(dòng)系統(tǒng)內(nèi)部沖擊對(duì)提升系統(tǒng)運(yùn)行品質(zhì)、降低振動(dòng)噪聲、延長(zhǎng)設(shè)備使用壽命具有十分重要的工程價(jià)值。

　　齒輪表面形貌直接影響齒輪傳動(dòng)動(dòng)力學(xué)特性。Majumdar 等研究表明基于分形幾何理論的 M-B 模型可以較好的模擬齒面粗糙形貌。通過(guò)大量實(shí)驗(yàn)得出分形參數(shù)與粗糙度之間的關(guān)系式為 這表明可以定量、準(zhǔn)確的實(shí)現(xiàn)分形參數(shù)與表面粗糙度之間的相互轉(zhuǎn)化。Yu 等通過(guò)分形幾何理論建立粗糙齒面模型，并修正了粗糙表面上存在微凸體時(shí)齒輪副的齒隙值。Zhao 等考慮了分形參數(shù)對(duì)齒輪接觸剛度的影響，隨著粗糙度的增大，時(shí)變嚙合剛度逐漸減小。Li 等將修正后的齒隙和嚙合剛度應(yīng)用到齒輪系統(tǒng)中，并分析了齒隙和剛度對(duì)系統(tǒng)動(dòng)力學(xué)特性的影響。然而，包括上述研究在內(nèi)，很少有基于各向異性三維粗糙齒面模型開展相關(guān)研究的。所以本研究首先基于分形幾何理論建立了三維各向異性分形粗糙齒面模型，然后對(duì)齒隙和時(shí)變嚙合剛度進(jìn)行了修正研究，同時(shí)創(chuàng)新性的引入了基于三維粗糙齒面的非均勻磨損，基于此研究了分形參數(shù)對(duì)系統(tǒng)動(dòng)力學(xué)特性的影響。本著降低振動(dòng)噪聲的原則，可以合理選擇使系統(tǒng)處于周期運(yùn)動(dòng)的分形參數(shù)。因此，從分形參數(shù)的角度實(shí)現(xiàn)了對(duì)齒輪傳動(dòng)系統(tǒng)振動(dòng)和噪聲的控制。研究思路和結(jié)果可以為齒輪加工和齒輪設(shè)計(jì)提供一定的理論指導(dǎo)。

　　一、三維分形粗糙齒面建模

　　齒輪齒面機(jī)加工后微觀凹凸不平，具有一定的表面粗糙度。同時(shí)實(shí)際加工齒面微觀形貌呈現(xiàn)出特定的方向性，因此本研究將 x 方向和 y 方向的兩個(gè)二維 M-B 分形函數(shù)疊加在一起構(gòu)建 x 、y 兩向異性的三維分形粗糙表面，建模公式為式(1)所示。

　　式中，G_x 表示 x 方向二維尺度系數(shù)，D_x 表示 x 方向二維分形維數(shù)，y 方向同理。Φ _m,n 是隨機(jī)相位，基于表面平整度和頻率分布密度的考慮，取γ=1.5。分形理論模擬的粗糙表面具有連續(xù)性、不可分性、尺度不變性和自親和性的特征，克服了表面粗糙度表征對(duì)樣本大小的依賴，使得對(duì)粗糙表面的描述更穩(wěn)定、更精確。

　　為了求解粗糙齒面對(duì)齒側(cè)間隙和嚙合剛度的影響，各向異性三維粗糙齒面建模過(guò)程如圖 1 所示。首先根據(jù)表 1 齒輪參數(shù)建立光滑齒面模型如圖 1(b)所示。本文將圖 1(a)所示的粗糙表面經(jīng)過(guò)旋轉(zhuǎn)后疊加到光滑齒面上即可得到各向異性三維粗糙齒面，如圖 1(c) 所示。其中圖 1(d)和(e)是粗糙齒面的局部細(xì)節(jié)圖。

圖 1 各向異性三維粗糙齒面建模過(guò)程

表 1 齒輪參數(shù)

　　二、齒隙的計(jì)算

　　相比于理想光滑齒面，粗糙齒面嚙合時(shí)由于微凸體的存在，齒隙值將發(fā)生改變。另外，粗糙齒面磨損率也與光滑齒面不同。因此，本研究設(shè)定初始齒隙值b₀為50 μm 。微凸體平均高度值記為bi ，粗糙齒面磨損深度值為b_w。綜合修正齒隙為b(t) = b₀ - b_i - b_w 。

　　分形參數(shù)對(duì)齒隙的影響：齒面粗糙體通過(guò)改變兩個(gè)輪齒相應(yīng)接觸區(qū)域之間的距離來(lái)改變齒隙值，如圖2所示。粗糙表面之間實(shí)際值取決于粗糙峰的平均高度。粗糙齒面導(dǎo)致的齒隙修正值應(yīng)為初始齒隙與表面微凸體高度值之差。根據(jù)分形函數(shù)(式 1)可以計(jì)算在(x₀ , y₀ ) 處單個(gè)微凸體的高度。同時(shí)沿齒廓方向?qū)⒒￠L(zhǎng)等分為n份。則第nj份區(qū)域內(nèi)微凸體的算術(shù)平均值高度為：

　　同理，可以得到主動(dòng)輪和從動(dòng)輪上所有微凸體絕對(duì)值的算術(shù)平均值bi _ p 和bi _ g 。則 微凸體導(dǎo)致的齒隙修正值為bi = bi _ p + bi _ g 。圖 3 為不同分形參數(shù)時(shí)凸體高度值對(duì)比圖。從圖中可以得出隨著分形維數(shù)增大、特征尺度系數(shù)減小，凸體平均高度逐漸減小。

圖 2 考慮三維粗糙齒面嚙合時(shí)的齒側(cè)間隙示意圖　　

圖 3 不同分形參數(shù)時(shí)齒面微凸體高度值

　　分形參數(shù)對(duì)磨損率的影響： 齒面粗糙度不同，其磨損速率必然不同。將齒面形貌特性融入到 Flodin 的齒面磨損研究中，求解不同分形參數(shù)時(shí)的齒面非均勻磨損速率。齒面磨損深度表達(dá)式為：

　　a_H 代表嚙合點(diǎn)的赫茲接觸寬度，k 是磨損系數(shù)，v是嚙合點(diǎn)處的滾動(dòng)速度，P 是嚙合點(diǎn)處的接觸載荷，具體計(jì)算方法參考文獻(xiàn)。ζ 是黏著磨損系數(shù)，ζ a_HkP 是黏著磨損量。

　　圖 4 為不同分形參數(shù)時(shí)齒面磨損深度，從圖中可以得出，隨著分形維數(shù)的減小、特征尺度系數(shù)(表面粗糙度)的增大磨損深度逐漸增大。綜上所述，考慮分形粗糙齒面的綜合齒隙如圖 5 所示。

圖 4 分形參數(shù)對(duì)齒面磨損深度的影響　　

圖 5 粗糙齒面?zhèn)鲃?dòng)時(shí)的綜合齒隙

　　三、基于三維分形齒面接觸模型的 TVMS 計(jì)算

　　時(shí)變嚙合剛度是齒輪傳動(dòng)系統(tǒng)的主要激勵(lì)因素之一，影響著系統(tǒng)動(dòng)力學(xué)特性的表現(xiàn)。勢(shì)能法是求解時(shí)變嚙合剛度使用較為廣泛的方法。它根據(jù)材料力學(xué)和彈性力學(xué)將時(shí)變嚙合剛度分解為軸向壓縮剛度k_a 、彎曲剛度k_b 、剪切剛度 k_s 、基體變形剛度k_f 以及赫茲接觸剛度kh進(jìn)行求解。在考慮三維分形齒面后，微觀上的齒面接觸實(shí)際上是微凸體之間的擠壓與滑動(dòng)，導(dǎo)致齒面實(shí)際接觸面積明顯小于理論接觸面積，所以傳統(tǒng)的赫茲接觸剛度不再適用。本文建立了三維各向異性分形齒面接觸模型，用于求解分形齒面接觸剛度。

　　在計(jì)算具有分形特征的接觸剛度時(shí)，通?？梢詫⒁粚?duì)齒輪之間的彈性接觸問(wèn)題轉(zhuǎn)化為具有等效曲率半徑和等效彈性模量的彈性圓柱與剛性平面的接觸。如圖 2(b)所示，其中 P 為剛體平面施加在橢圓微凸體上的法向接觸載荷，δ_n 為橢圓微凸體的高度，ω_n 為橢圓微凸體的實(shí)際變形量，取值在(0，δ_n ) 之間。由式(1)可得以橢圓區(qū)域的長(zhǎng)軸l _x 和短軸l _y 為基底的微凸體的輪廓曲線 z(x, y)，即:

　　單個(gè)橢圓微凸體的接觸面積a 為：

　　其中e 表示橢圓微凸體接觸的偏心量。當(dāng)微凸體變形達(dá)到臨界變形δ_c 時(shí)，臨界接觸面積 a_c 的表達(dá)式為：

　　微凸體的最大接觸面積為：

　　齒輪副分形齒面接觸面積分布函數(shù)可表示為：

　　其中ζ 表示橢圓微凸體的接觸系數(shù)，

　　當(dāng)橢圓微凸體處于完全彈性變形時(shí)，法向接觸剛度kh可表示為：

　　另外，上述公式中r_x 、r_y 、e 、δ_n 和ω_n 等可在文獻(xiàn)中查得?？倳r(shí)變嚙合剛度可表示為：

　　考慮分形齒面的 TVMS 計(jì)算結(jié)果如圖 6 所示。從圖中可以看出隨著D的增大或G 的減小，嚙合剛度逐漸增大。這是由于D的增大或G的減小會(huì)造成分形齒面接觸面積增大，總體變形量減小，最終導(dǎo)致剛度值增大。

圖 6 分形參數(shù)對(duì)嚙合剛度的影響

　　四、分形參數(shù)對(duì)系統(tǒng)特性的影響

　　齒輪傳動(dòng)系統(tǒng)動(dòng)力學(xué)模型參考文獻(xiàn)建立，不再詳細(xì)介紹。通過(guò)上述研究可以得出分形參數(shù)對(duì)齒隙和嚙合剛度值的變化具有顯著影響，這勢(shì)必會(huì)導(dǎo)致齒輪系統(tǒng)動(dòng)態(tài)特性的改變。設(shè)定模型參數(shù)值為：b₀ = 50 μm ，e₀ = 25 μm ，μ₁ = μ₂ = 0.1 ，c_px = c_gx = c_py = c_gy = 3×10³ Ns/m，k_px = k_gx = k_py = k_gy = 6×10⁸ Ns/m 。將不同分形參數(shù)時(shí)的齒隙和嚙合剛度值帶入動(dòng)力學(xué)模型中，采用 Runge‐Kutta 算法對(duì)動(dòng)力學(xué)方程求解。

　　將轉(zhuǎn)速作為變量參數(shù)，得到了不同轉(zhuǎn)速下系統(tǒng)豐富的動(dòng)力學(xué)特性，如圖 7 所示。隨著轉(zhuǎn)速的增加系統(tǒng)經(jīng)歷了 1 周期、2 周期、4 周期和混沌運(yùn)動(dòng)等多種響應(yīng)狀態(tài)，另外混沌意味著齒輪傳動(dòng)系統(tǒng)穩(wěn)定性差，不利于系統(tǒng)穩(wěn)定運(yùn)行。因此，可以通過(guò)選擇避開圖中混沌運(yùn)動(dòng)區(qū)域來(lái)增強(qiáng)系統(tǒng)穩(wěn)定性，減少振動(dòng)以達(dá)到降噪的作用。

圖 7 不同轉(zhuǎn)速時(shí)系統(tǒng)動(dòng)態(tài)行為

　　分形維數(shù)對(duì)齒輪振動(dòng)特性的影響：選取 3600r/min 和 6000r/min 兩個(gè)一般轉(zhuǎn)速工況，以分形維數(shù) D_x 和 D_y 為變量參數(shù)、特征尺度系數(shù)G_x = G_y = 5×10^-7 m 、其余參數(shù)不變時(shí)系統(tǒng)進(jìn)行動(dòng)態(tài)響應(yīng)研究，結(jié)果如圖 8 所示。

圖 8 不同分形維數(shù)( D_x 和 D_y )時(shí)系統(tǒng)動(dòng)態(tài)行為

　　從圖 8 中可以看出，隨著齒輪旋轉(zhuǎn)速度的增加(由3000r/min 提升到 6000r/min)， 齒輪傳動(dòng)系統(tǒng)動(dòng)力學(xué)響應(yīng)發(fā)生明顯變化，系統(tǒng)穩(wěn)定性顯著降低。這表明隨著齒輪系統(tǒng)旋轉(zhuǎn)速度的增加，系統(tǒng)振動(dòng)明顯加劇。仔細(xì)觀察特定轉(zhuǎn)速時(shí)的分岔圖可以發(fā)現(xiàn)，在分形維數(shù) D_x 和 D_y 均較小時(shí)，動(dòng)態(tài)傳動(dòng)誤差幅值較大，傳動(dòng)系統(tǒng)振動(dòng)劇烈。隨著 D_x 或 D_y 單因素的增加，DTE 幅值快速減小，在 Dx 和 Dy 都增大到 1.9 時(shí)，DTE 幅值最小，表明此刻系統(tǒng)傳動(dòng)最為平穩(wěn)。

　　特征尺度系數(shù)對(duì)齒輪振動(dòng)特性的影響： 以特征尺度系數(shù)Gx 和Gy 為變量參數(shù)、分形維數(shù) D_x = D_y =1.4 、其余參數(shù)不變對(duì)系統(tǒng)進(jìn)行動(dòng)態(tài)響應(yīng)研究，結(jié)果如圖 9 所示。隨著齒輪旋轉(zhuǎn)速度的增加，系統(tǒng)穩(wěn)定性明顯降低，這與 4.1 中的結(jié)論一致，均屬于工況對(duì)傳動(dòng)系統(tǒng)的固有影響，不可改變。觀察特征尺度系數(shù)G_x 和G_y 對(duì)系統(tǒng)動(dòng)力學(xué)響應(yīng)似乎影響并不明顯。仔細(xì)觀察可以得出隨著表面特征尺度系數(shù)的增大，系統(tǒng)穩(wěn)定性逐漸降低，振動(dòng)愈發(fā)劇烈。所以為了系統(tǒng)的穩(wěn)定運(yùn)行，可以在加工齒輪表面時(shí)適當(dāng)?shù)脑龃蠓中尉S數(shù)，減小特征尺度系數(shù)。而對(duì)于特殊用途(如需要振動(dòng)的場(chǎng)景)的傳動(dòng)齒輪，在齒面加工時(shí)可以適當(dāng)減小分形維數(shù)而增大特征尺度系數(shù)。

圖 9 不同特征尺度系數(shù)(Gx 和Gy )時(shí)系統(tǒng)動(dòng)態(tài)行為

　　五、結(jié)論

　　利用分形幾何理論建立了各向異性三維粗糙齒面模型，并對(duì)齒輪嚙合時(shí)的齒隙和時(shí)變剛度進(jìn)行修正。研究了齒面形貌對(duì)系統(tǒng)動(dòng)力學(xué)特性的影響。主要結(jié)論如下：

　　(1)通過(guò)改變分形參數(shù)可以得到特定的各向異性粗糙齒面模型，分形維數(shù)越小、特征尺度系數(shù)越大，表面越粗糙，粗糙凸體引起的齒隙修正值越大。

　　(2)隨著表面粗糙度增大，時(shí)變嚙合剛度值減小，齒面磨損速率增加，磨損深度在靠近節(jié)點(diǎn)的方向逐漸減小。

　　(3)在相同表面粗糙度時(shí)，隨著齒輪轉(zhuǎn)速的增大，系統(tǒng)穩(wěn)定性降低，動(dòng)態(tài)傳動(dòng)誤差幅值增大。

　　(4)隨著分形維數(shù)的增大和特征尺度系數(shù)的減小，系統(tǒng)具有更高的穩(wěn)定性，可以有效降低齒輪傳動(dòng)系統(tǒng)的振動(dòng)特性。

　　參考文獻(xiàn)略.