王 星, 方宗德, 牟彥銘， 杜進(jìn)輔， 崔艷梅

(1. 西安科技大學(xué) 機(jī)械工程學(xué)院， 西安 710054； 2. 西北工業(yè)大學(xué) 機(jī)電學(xué)院， 西安 710072；3. 西安理工大學(xué) 機(jī)械與精密儀器工程學(xué)院， 西安 710048； 4. 鄭州航空工業(yè)管理學(xué)院 機(jī)電學(xué)院， 鄭州 450015)

HGT準(zhǔn)雙曲面齒輪承載傳動(dòng)誤差的優(yōu)化設(shè)計(jì)

王 星1, 方宗德2, 牟彥銘2， 杜進(jìn)輔3， 崔艷梅4

(1. 西安科技大學(xué) 機(jī)械工程學(xué)院， 西安 710054； 2. 西北工業(yè)大學(xué) 機(jī)電學(xué)院， 西安 710072；3. 西安理工大學(xué) 機(jī)械與精密儀器工程學(xué)院， 西安 710048； 4. 鄭州航空工業(yè)管理學(xué)院 機(jī)電學(xué)院， 鄭州 450015)

HGT加工方法； 準(zhǔn)雙曲面齒輪； 加工參數(shù)設(shè)計(jì)； 承載傳動(dòng)誤差； 優(yōu)化設(shè)計(jì)

準(zhǔn)雙曲面齒輪被廣泛用作汽車主減速器齒輪，是齒輪傳動(dòng)中最為復(fù)雜的一種錐齒輪，用來傳遞兩相錯(cuò)軸之間的運(yùn)動(dòng)和動(dòng)力。具有重合度大、傳動(dòng)效率高、傳動(dòng)平穩(wěn)、噪聲小以及傳動(dòng)比大(10∶1，60∶1 以至于100∶1)等優(yōu)點(diǎn)，且沿齒長和齒高方向都有相對(duì)滑動(dòng)，易于跑合，熱處理后便于研磨，有利于改善齒面接觸區(qū)、提高齒面光潔度和降低噪聲等一系列優(yōu)點(diǎn)。另外，偏置距為零的弧齒錐齒輪只是準(zhǔn)雙曲面齒輪的一個(gè)特例，因此，對(duì)準(zhǔn)雙曲面齒輪進(jìn)行研究具有一般性的意義。

楊宏斌等[1]研究了基于局部綜合法的HFT方法加工高齒準(zhǔn)雙曲面齒輪的加工參數(shù)設(shè)計(jì)。方宗德等[2]通過局部綜合法對(duì)準(zhǔn)雙曲面齒輪的優(yōu)化切齒設(shè)計(jì)進(jìn)行了前期理論探索，并沒有通過算例給以定量計(jì)算和驗(yàn)證。王星等[3]曾對(duì)HGT準(zhǔn)雙曲面齒輪進(jìn)行了工作齒面的理論推導(dǎo)，并進(jìn)行了TCA計(jì)算。SIMON[4-5]研究了如何在數(shù)控機(jī)床上加工準(zhǔn)雙曲面齒輪。NISHINO[6]對(duì)面滾式準(zhǔn)雙曲面齒輪的加載接觸性能通過數(shù)值方法進(jìn)行了模擬。SIMON[7]針對(duì)面滾式準(zhǔn)雙曲面齒輪提出了一種優(yōu)化方法，用來系統(tǒng)地定義刀盤參數(shù)和機(jī)床加工參數(shù)。

隨著工業(yè)技術(shù)的發(fā)展以及對(duì)齒輪轉(zhuǎn)速和承載能力要求的提高，齒輪傳動(dòng)過程中的振動(dòng)噪聲已日益成為影響產(chǎn)品質(zhì)量的瓶頸問題，齒輪傳動(dòng)的動(dòng)態(tài)設(shè)計(jì)與動(dòng)態(tài)質(zhì)量控制已是現(xiàn)代工業(yè)必須解決的問題。承載傳動(dòng)誤差對(duì)準(zhǔn)雙曲面齒輪的動(dòng)態(tài)特性有著直接的影響，承載傳動(dòng)誤差曲線的波動(dòng)程度可反映出齒輪副的動(dòng)態(tài)性能，波動(dòng)幅值越大，噪音越大，反之噪音越小，傳動(dòng)越平穩(wěn)，而重合度則決定了工作載荷下的傳動(dòng)誤差，當(dāng)重合度較大時(shí)，有利于改善齒輪傳動(dòng)的平穩(wěn)性和齒間載荷分配，提高齒輪動(dòng)態(tài)和強(qiáng)度性能[8-9]。WANG等[10]以TCA和LTCA為工具，以承載傳動(dòng)誤差曲線幅值等最小為優(yōu)化目標(biāo)，來確定人字齒輪最佳修行齒面，達(dá)到了減振的目的。王峰等[11]對(duì)齒輪傳動(dòng)系統(tǒng)(齒輪和箱體)進(jìn)行了減振設(shè)計(jì)。LIN等[12-15]分析了齒面修形對(duì)齒輪承載傳動(dòng)誤差、齒面沖擊的影響，并結(jié)合動(dòng)力特性分析了減振效果。LIU等[16]通過建立考慮動(dòng)載荷齒間分配、時(shí)變剛度和接觸損失的非線性分析模型，研究了齒廓修形對(duì)降低多級(jí)齒輪傳動(dòng)系統(tǒng)振動(dòng)的作用。

論文首先基于局部綜合法對(duì)HGT準(zhǔn)雙曲面齒輪進(jìn)行加工參數(shù)設(shè)計(jì)，然后對(duì)齒輪副進(jìn)行承載傳動(dòng)誤差[17]幅值最小的加工參數(shù)優(yōu)化設(shè)計(jì)，即優(yōu)化承載傳動(dòng)誤差曲線幅值最小，以提高齒輪副的動(dòng)態(tài)特性。

1 加工參數(shù)設(shè)計(jì)

近年來，隨著汽車技術(shù)的不斷發(fā)展，對(duì)準(zhǔn)雙曲面齒輪嚙合質(zhì)量的要求越來越高，而這需要通過正確的切齒參數(shù)的調(diào)整來實(shí)現(xiàn)。本文以局部綜合法為基礎(chǔ)對(duì)HGT準(zhǔn)雙曲面齒輪進(jìn)行加工參數(shù)設(shè)計(jì)，以達(dá)到預(yù)控齒面嚙合性能的目的。這里僅對(duì)加工參數(shù)設(shè)計(jì)過程做簡(jiǎn)要敘述。

(1)

圖1 加工參數(shù)設(shè)計(jì)流程圖

(a)

(b)

大輪加工參數(shù)計(jì)算相對(duì)簡(jiǎn)單，主要包括:刀盤齒形角、刀盤半徑、刀頂距、垂直刀位、水平刀位、垂直輪位、軸向輪位、床位、輪坯安裝角、滾比和刀尖圓角半徑。

小輪加工參數(shù)包括：刀盤齒形角、刀尖半徑、刀傾角、刀轉(zhuǎn)角、垂直刀位、水平刀位、垂直輪位、軸向輪位、床位、輪坯安裝角、滾比和刀尖圓角半徑。

2 優(yōu)化方法

準(zhǔn)雙曲面齒輪承載傳動(dòng)誤差的優(yōu)化過程是通過不斷改變小輪加工參數(shù)，來改變齒面接觸狀況的迭代過程，是一個(gè)非線性接觸問題，因此在實(shí)際中很難建立優(yōu)化變量和目標(biāo)函數(shù)的解析表達(dá)式，而且在優(yōu)化過程中存在多個(gè)局部最優(yōu)解，傳統(tǒng)算法不適用于該問題的求解，因此本章采用遺傳算法(Genetic Algorithm)對(duì)承載傳動(dòng)誤差進(jìn)行優(yōu)化設(shè)計(jì)。

遺傳算法的主要特點(diǎn)是：① 廣泛的適用性，可直接對(duì)結(jié)果對(duì)象進(jìn)行操作；② 群體搜索，易于并行化處理；③ 求解時(shí)使用特定問題的信息極少，容易形成通用算法程序；④ 遺傳算法屬于自適應(yīng)概率搜索技術(shù)。選擇、交叉和變異都是隨機(jī)操作，而不是確定的精確規(guī)則，選擇體現(xiàn)了向最優(yōu)解迫近，交叉體現(xiàn)了最優(yōu)解的產(chǎn)生，變異體現(xiàn)了全局最優(yōu)解的覆蓋。

遺傳算法優(yōu)化過程如圖3所示。

圖3 遺傳算法優(yōu)化流程圖

3 優(yōu)化設(shè)計(jì)模型

3.1 設(shè)計(jì)變量

設(shè)計(jì)變量可表示為

(2)

3.2 目標(biāo)函數(shù)

準(zhǔn)雙曲面齒輪動(dòng)態(tài)性能與承載傳動(dòng)誤差有著直接的關(guān)系，因此準(zhǔn)雙曲面齒輪的振動(dòng)可以用承載傳動(dòng)誤差Tl來描述，在齒輪傳動(dòng)中，當(dāng)Tl恒定時(shí)，齒輪就不會(huì)產(chǎn)生振動(dòng)，因此優(yōu)化的目標(biāo)就是設(shè)法將齒輪傳動(dòng)中Tl的變化量(即承載傳動(dòng)誤差曲線的幅值)控制在最小的范圍內(nèi)，即：

(3)

式中，Tli為第i個(gè)載荷對(duì)應(yīng)的承載傳動(dòng)誤差。

3.3 約束條件

(4)

(2)由于HGT準(zhǔn)雙曲面齒輪大輪和小輪都采用展成法加工，齒面曲率特性好，接觸跡線為直線，所以為實(shí)現(xiàn)齒輪副的高重合度設(shè)計(jì)提供了保障，因此，為了保證齒輪副傳動(dòng)的連續(xù)性與平穩(wěn)性以及防止發(fā)生邊緣接觸時(shí)齒輪副運(yùn)動(dòng)的強(qiáng)振動(dòng)性, 約束齒輪副的設(shè)計(jì)重合度大于等于2.0，η2的取值范圍設(shè)定為

(5)

綜上所述，承載傳動(dòng)誤差優(yōu)化設(shè)計(jì)流程圖如圖4所示。

圖4 承載傳動(dòng)誤差優(yōu)化過程

4 計(jì)算實(shí)例與分析

由圖4可以看出，優(yōu)化過程較為復(fù)雜，為了縮短程序運(yùn)行時(shí)間，采用小種群遺傳，種群大小取10條染色體(傳統(tǒng)遺傳算法要求種群規(guī)模為50～100條染色體)。另外為了減小不必要的重復(fù)計(jì)算，當(dāng)連續(xù)迭代5次求得的適應(yīng)度值不發(fā)生改變時(shí)，即可認(rèn)為達(dá)到最優(yōu)結(jié)果，算法終止；否則，當(dāng)遺傳算法連續(xù)迭代20次時(shí)，迭代終止。停滯代數(shù)為5，交叉概率為0.7，變異概率為0.01，適應(yīng)度函數(shù)的終值公差為1E-100。

4.1 承載傳動(dòng)誤差幅值最小的優(yōu)化設(shè)計(jì)

由于汽車用準(zhǔn)雙曲面齒輪多處于正轉(zhuǎn)狀態(tài)，因此只對(duì)齒輪副工作面(大輪凸面和小輪凹面)進(jìn)行承載傳動(dòng)誤差幅值最小的優(yōu)化設(shè)計(jì)。

表1 優(yōu)化前加工參數(shù)

由圖5可以看出，接觸跡線接近直線，且優(yōu)化后接觸跡線長于優(yōu)化前，優(yōu)化前設(shè)計(jì)重合度[18]為2.15，優(yōu)化后設(shè)計(jì)重合度為2.25。

表2 優(yōu)化后加工參數(shù)

分別在大輪常用加載扭矩800 N·m和1 500 N·m下，比較優(yōu)化前后的承載傳動(dòng)誤差曲線，承載傳動(dòng)誤差曲線如圖6(a)和圖6(b)所示。

由圖6(a)可以看出，大輪加載扭矩為800 N·m時(shí)，優(yōu)化前承載傳動(dòng)誤差幅值為3.560 Sec，優(yōu)化后承載傳動(dòng)誤差幅值為2.210 Sec，因此優(yōu)化后幅值降低了37.92%。

由圖6(b)可以看出，大輪加載扭矩為1 500 N·m時(shí)，優(yōu)化前承載傳動(dòng)誤差幅值為4.950 Sec，優(yōu)化后承載傳動(dòng)誤差幅值為4.130 Sec，因此優(yōu)化后幅值降低了16.57%，這一點(diǎn)說明文中雖然是在800 N·m載荷下對(duì)齒面進(jìn)行優(yōu)化設(shè)計(jì)，但是在一定的載荷范圍內(nèi)，齒輪副的承載傳動(dòng)誤差曲線幅值可能也會(huì)相應(yīng)減小。

因此，在大輪加載扭矩分別為800 N·m和1 500 N·m時(shí)，相同載荷下承載傳動(dòng)誤差曲線優(yōu)化后較優(yōu)化前幅值減小，波動(dòng)趨于平緩，證明了優(yōu)化方法的正確性和可行性。

4.2 多載荷下的承載傳動(dòng)誤差曲線

針對(duì)汽車工況多變的特點(diǎn)，對(duì)多載荷下的承載傳動(dòng)誤差進(jìn)行計(jì)算，優(yōu)化前后(大輪加載扭矩800 N·m時(shí)優(yōu)化)的計(jì)算結(jié)果如圖7所示。

由圖7可知，優(yōu)化后承載傳動(dòng)誤差較優(yōu)化前波動(dòng)減小，且在大輪加載扭矩為5 N·m時(shí)，承載傳動(dòng)誤差曲線與幾何接觸傳動(dòng)誤差曲線基本重合，證明當(dāng)載荷很小時(shí)，齒輪副的實(shí)際重合度為1，此時(shí)齒面接觸跡線上只有中間一段參與嚙合。另外，當(dāng)載荷小于600 N·m時(shí)，在一個(gè)嚙合周期中齒輪副單齒和雙齒交替嚙合；當(dāng)載荷大于600 N·m時(shí)，在一個(gè)嚙合周期中齒輪副雙齒和三齒交替嚙合，說明隨著載荷的增大，相應(yīng)的變形補(bǔ)償了第二對(duì)齒間的間隙，齒輪副的實(shí)際重合度增大，當(dāng)載荷增大到小輪或大輪的齒頂達(dá)到接觸時(shí)(即承載傳動(dòng)誤差曲線位于幾何接觸傳動(dòng)誤差曲線的末端)，實(shí)際重合度達(dá)到設(shè)計(jì)最大重合度，隨載荷繼續(xù)增大，則出現(xiàn)邊緣接觸(即承載傳動(dòng)誤差曲線低于幾何接觸傳動(dòng)誤差曲線的末端)，此時(shí)重合度略有增加，但嚴(yán)重的邊緣接觸易引起強(qiáng)度破壞和振動(dòng)噪音，應(yīng)予以避免。

(a) 優(yōu)化前

(b) 優(yōu)化后

承載傳動(dòng)誤差幅值對(duì)比圖如圖8所示。

(a) 800 N·m

(b) 1 500 N·m

由圖8可以看出，優(yōu)化后承載傳動(dòng)誤差幅值和優(yōu)化前相比，在載荷為100 N·m時(shí)降低了50.62%。且優(yōu)化前后曲線各有兩個(gè)局部最優(yōu)點(diǎn)，優(yōu)化前局部最優(yōu)點(diǎn)分別出現(xiàn)在300 N·m和1 000 N·m載荷下，優(yōu)化后局部最優(yōu)點(diǎn)分別出現(xiàn)在100 N·m和800 N·m載荷下。

4.3 承載傳動(dòng)誤差幅值隨載荷的變化

5 結(jié) 論

(1) 論文以HGT準(zhǔn)雙曲面齒輪幾何接觸傳動(dòng)誤差曲線轉(zhuǎn)換點(diǎn)的幅值和接觸跡線與根錐的夾角為優(yōu)化變量，以承載傳動(dòng)誤差波動(dòng)幅值最小為目標(biāo)函數(shù)，通過遺傳算法對(duì)小輪加工參數(shù)進(jìn)行了優(yōu)化設(shè)計(jì)。在大輪加載扭矩分別為800 N·m和1 500 N·m時(shí)，相同載荷下承載傳動(dòng)誤差曲線優(yōu)化后較優(yōu)化前波動(dòng)趨于平緩，優(yōu)化后承載傳動(dòng)誤差幅值分別降低了37.92%和16.57%，證明了優(yōu)化方法的正確性和可行性。

(a) 優(yōu)化前

(b) 優(yōu)化后

圖8 承載傳動(dòng)誤差幅值對(duì)比圖

圖9 設(shè)計(jì)重合度小于等于2時(shí)承載傳動(dòng)誤差幅值隨載荷的變化

Fig.9 Variation with load of loaded transmission error amplitude when the design coincidence degree is less than or equal to 2

圖10 設(shè)計(jì)重合度大于2小于3時(shí)承載傳動(dòng)誤差幅值隨載荷的變化

Fig.10 Variation with load of loaded transmission error amplitude when the design coincidence degree is greater than 2 less than 3

(2) 針對(duì)汽車工況多變的特點(diǎn)，對(duì)多載荷下的承載傳動(dòng)誤差進(jìn)行計(jì)算，研究發(fā)現(xiàn)當(dāng)設(shè)計(jì)重合度小于等于2時(shí)，承載傳動(dòng)誤差幅值曲線有1個(gè)局部最優(yōu)點(diǎn)，當(dāng)設(shè)計(jì)重合度大于2小于3時(shí)，承載傳動(dòng)誤差幅值曲線有2個(gè)局部最優(yōu)點(diǎn)，為了保持齒輪副較好的振動(dòng)特性，應(yīng)使其盡量在局部最小幅值對(duì)應(yīng)的載荷附近工作。

[1] 楊宏斌, 范明, 周彥偉，等. 高齒準(zhǔn)雙曲面齒輪的研究[J]. 中國機(jī)械工程,2000,11(8):897-900.

YANG Hongbin, FAN Ming, ZHOU Yanwei. Research on high tooth hypoid gear[J]. China Mechanical Engineering, 2000,11(8):897-900.

[2] 方宗德,楊宏斌. 準(zhǔn)雙曲面齒輪的優(yōu)化切齒設(shè)計(jì)[J].汽車工程,1998,20(5):302-307.

FANG Zongde, YANG Hongbin. Optimal cutting design for hypoid gears[J]. Automotive Engineering,1998,20(5):302-307.

[3] 王星, 方宗德, 李聲晉. HGT準(zhǔn)雙曲面齒輪傳動(dòng)的輪齒接觸分析[J]. 西北工業(yè)大學(xué)學(xué)報(bào), 2014, 32(3): 475-480.

WANG Xing, FANG Zongde, LI Shengjin. Tooth contact analysis of HGT hypoid gear drives[J]. Journal of Northwestern Polytechnical Uniersity, 2014, 32(3): 475-480.

[4] SIMON V V. Advanced manufacture of spiral bevel gears on CNC hypoid generating machine[J]. Journal of Mechanical Design, 2010, 132(3): 031001.

[5] SIMON V V. Generation of hypoid gears on CNC hypoid generator[J]. Journal of Mechanical Design, 2011, 133(12): 121003.

[6] NISHINO T. Generation and curvature analysis of face-hobbed hypoid gears[C]//Proceedings of JSME International Conference on Motion and Power Transmissions. Sendai,2009: 64-69.

[7] SIMON V V. Optimization of face-hobbed hypoid gears[J]. Mechanism and Machine Theory, 2014, 77: 164-181.

[8] GAO Jianping, FANG Zongde. Non-linear dynamic analysis of hypoid gear drive ang spiral bevel gears[C]∥The Fourth World Gear and Power Transmission Conference. France, 1999.

[9] 曹雪梅, 周彥偉, 鄧效忠,等. 弧齒錐齒輪的傳動(dòng)誤差、重合度和噪聲的關(guān)系[J]. 機(jī)械傳動(dòng), 2003, 27(2):43-45.

CAO Xuemei, ZHOU Yanwei, DENG Xiaozhong. The relationship between transmission error coincidence degree and noise for the spiral bevel gear[J]. Journal of Mechanical, 2003, 27(2):43-45.

[10] WANG C, FANG Z D, JIA H T. Investigation of a design modification for double helical gears reducing vibration and noise[J]. Journal of Marine Science and Application, 2010, 9(1): 81-86.

[11] 王峰, 方宗德, 李聲晉，等. 多工況人字齒輪傳動(dòng)系統(tǒng)結(jié)構(gòu)減振分析及優(yōu)化[J]. 振動(dòng)與沖擊, 2014, 33(22):127-130.

WANG Feng, FANG Zongde, LI Shengjin， et al. Structural vibration analysis and optimization of herringbone gear transmission system under multiple loads[J]. Journal of Vibration and Shock, 2014, 33(22):127-130.

[12] LIN H H, OSWALD F B, TOWNSEND D P. Dynamic loading of spur gears with linear or parabolic tooth profile modifications[R]. National Aeronautics and Space Administration Cleveland Ohlewis Research Center, 1989.

[13] WAGAJ P, KAHRAMAN A. Impact of tooth profile modifications on the transmission error excitation of helical gear pairs[C]//Proceedings of ESDA2002: 6th Biennial Conference on Engineering Systems Design and Analysis. Istanbul, 2002.

[14] UMEZAWA K. Low vibration design on helical gear pair[M]. AGMA, 1998.

[15] UMEYAMA M, KATO M, INOUE K. Effects of gear dimensions and tooth surface modifications on the Loaded transmission error of a helical gear pair[J]. Journal of Mechanical Design, 1998, 120(1): 119-125.

[16] LIU G, PARKER R G. Dynamic modeling and analysis of tooth profile modification for multimesh gear vibration[J]. Journal of Mechanical Design, 2008, 130: 121402.

[17] 方宗德. 齒輪輪齒承載接觸分析(LTCA)的模型和方法[J]. 機(jī)械傳動(dòng)，1998, 13(9):1-3.

FANG Zongde. Model and method of loaded tooth contact analysis (LTCA)[J]. Journal of Mechanical, 1998, 13(9):1-3.

[18] 鄧效忠，方宗德，張金良,等. 弧齒錐齒輪的高重合度設(shè)計(jì)[J]. 中國機(jī)械工程，2002, 22(2):791-795.

DENG Xiaozhong, FANG Zongde, ZHANG Jinliang, et al. High degree design for Spiral Bevel Gears[J]. China Mechanical Engineering, 2002, 13(9):791-795.

Optimization design of loaded transmission error for HGT hypoid gear drives

WANG Xing1, FANG Zongde2, MU Yanming2, DU Jinfu3, CUI Yanmei4

(1. School of Mechanical Engineering, Xi’an University of technology, Xi’an 710054, China;2. School of Mechanical Engineering, Northwestern Polytechnical University, Xi’an 710072, China;3. School of Mechanical and Precision Instrument Engineering, Xi’an University of Technology, Xi’an 710048, China;4. School of Mechatronics Engineering, Zhengzhou Institute of Aeronautical Industry Management, Zhengzhou 450015, China)

The fluctuation degree of the loaded transmission error curve can reflect the dynamic performance of a gear pair. The greater the fluctuation amplitude, the bigger the noise; the smaller the fluctuation amplitude, the smaller the noise and the transmission is more stable. Firstly, based on the local synthesis method (Synthesis Local), and according to the processing principle of the Gleason hypoid gear, the processing parameters of the HGT hypoid gear were designed. On this basis, the first order derivative of transmission ratio function and the angle between the contact trace and the root cone were treated the optimization variables, and the amplitude of the loaded transmission was treated as the objective function, then the optimization design for the processing parameters was done through the genetic algorithm, for improving the dynamic characteristics of the gear pair. It is found that: ① when the gear load torque is 800 N·m and 1 500 N·m, the amplitude of the loaded transmission error is reduced by 37.92% and 16.57% respectively; ② in order to keep the desirable vibration characteristics of gear pair, it should be to work near the load of local minimum amplitude; and with the increase of, the local minimum amplitude is moved to the direction of large load; it illustrates that, in order to make the gear pair with a smaller load bearing capacity, the amplitude of the loaded transmission error is larger.

HGT processing method; hypoid gear; design of machining parameters; loaded transmission error; optimization design

國家自然科學(xué)基金資助項(xiàng)目(51375384; 51175423); 航空科學(xué)基金(2015ZB55002)

2015-12-07 修改稿收到日期：2016-03-21

王星 女，博士生，1982年生

方宗德 男，教授，博士生導(dǎo)師，1948年生 E-mail: fauto@nwpu.edu.cn

TH132.41

A

10.13465/j.cnki.jvs.2017.08.006