周鎮(zhèn)宇, 唐進(jìn)元, 董建雄

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

?

同功重比修形斜齒和直齒面齒輪性能對比研究

周鎮(zhèn)宇, 唐進(jìn)元, 董建雄

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

為深入了解同功重比修形斜齒與直齒面齒輪的性能差異，選擇更適合于高速重載工況下的面齒輪傳動.基于嚙合原理推導(dǎo)了修形斜齒與直齒面齒輪齒面方程，基于CATIA建立了修形斜齒與直齒面齒輪三維模型，采用有限元接觸分析方法，以接觸應(yīng)力、彎曲應(yīng)力和重合度為面齒輪傳動性能指標(biāo)展開研究.研究結(jié)果表明：修形斜齒面齒輪相比修形直齒面齒輪接觸應(yīng)力大幅降低，算例最大接觸應(yīng)力降低16.3%；修形斜齒面齒輪相比修形直齒面齒輪彎曲應(yīng)力大幅降低，算例最大彎曲應(yīng)力降低32.4%；修形斜齒面齒輪相比修形直齒面齒輪重合度大幅提高，算例重合度提高10.3%.所以同功重比情況下，修形斜齒面齒輪傳動性能優(yōu)于修形直齒面齒輪，前者更適合于高速重載工況下的輕量化設(shè)計(jì).

輕量化; 功重比; 接觸應(yīng)力; 彎曲應(yīng)力; 重合度; 修形

面齒輪傳動是一種相對較新的齒輪傳動方式.正交面齒輪傳動中,與直齒輪嚙合傳動的稱為直齒面齒輪,與斜齒輪嚙合傳動的稱為斜齒面齒輪.面齒輪目前主要應(yīng)用于武裝直升機(jī)主傳動中,而輕量化設(shè)計(jì)已成為飛機(jī)、汽車設(shè)計(jì)的一大趨勢[1].對直升機(jī)主減速器而言,高功率密度是其主要性能指標(biāo),因此研究面齒輪功重比(單位質(zhì)量所承載的功率)非常有必要.

斜齒面齒輪和直齒面齒輪均易出現(xiàn)邊緣接觸.為改善接觸軌跡的齒面分布,避免邊緣接觸,得到均勻的載荷分布,需對面齒輪進(jìn)行修形.Litvin等[2-8]率先提出拋物線齒廓修形方法,并在斜齒面齒輪接觸分析中提出點(diǎn)接觸形成條件和安裝誤差坐標(biāo)系,得到均勻的齒面接觸軌跡；初步給出面齒輪重合度的定義,但未給出具體計(jì)算方法.雷敦財(cái)[9]研究了有限元分析中輪轂結(jié)構(gòu)對嚙合剛度的影響.

本文基于有限元分析方法,在相同功重比工況下對修形斜齒面齒輪和修形直齒面齒輪進(jìn)行接觸應(yīng)力、彎曲應(yīng)力和重合度對比研究.

1　修形斜齒和直齒面齒輪齒面方程

齒輪嚙合傳動中主、被動齒輪的基節(jié)必須處處相等,從理論上講,精確的漸開線剛性齒輪是能夠?qū)崿F(xiàn)上述目標(biāo)的.但實(shí)際的齒輪副均為彈性體,在一定嚙合力作用下會產(chǎn)生相應(yīng)彈性變形,使處于嚙合線位置的主動輪和被動輪基節(jié)出現(xiàn)變化,嚙合過程中產(chǎn)生嚙合沖擊.為了消除輪齒嚙合沖擊,通常采用齒廓修形方法,即沿齒高方向從齒面上去除部分材料,從而改變齒廓形狀,消除嚙合沖擊,減小嚙合過程中的最大應(yīng)力,提高齒輪副的傳動平穩(wěn)性和使用壽命.而拋物線齒廓是常用的修形齒廓.

考慮到標(biāo)準(zhǔn)的斜齒面齒輪和直齒面齒輪在嚙合傳動過程中均易產(chǎn)生邊緣接觸,故對這2種面齒輪均進(jìn)行拋物線齒廓修形.修形面齒輪的齒面由修形產(chǎn)形輪齒面包絡(luò)得到[10],而修形產(chǎn)形輪齒面由修形齒條刀具包絡(luò)得到的,修形齒條齒廓直線被拋物線代替.

1.1修形斜齒面齒輪齒面方程

傳統(tǒng)加工斜齒輪的齒條刀刃∑c是直線,端面如圖1(a)所示,虛線所示∑cs和∑c1為修形拋物線齒廓,頂點(diǎn)Oa的位置由uo決定；空間齒廓如圖1(b)所示,∑i(i=cs,cl)為修形齒條齒面,li(i=cs,cl)為齒寬縱向變量參數(shù)；圖1(c)為拋物線修形齒條端面兩側(cè)齒廓,Sc(xc,yc,zc)為齒條動標(biāo)系,Sa(xa,ya,za)和Sb(xb,yb,zb)為輔助坐標(biāo)系,ui(i=s,l)為坐標(biāo)系Sa(xa,ya,za)中橫向變量參數(shù)；圖1(d)所示為修形齒條加工修形產(chǎn)形輪的坐標(biāo)系,Sm(xm,ym,zm)為靜標(biāo)系,rpi(i=s,l)為節(jié)圓半徑,Si(xi,yi,zi)為

圖1　修形斜齒輪加工坐標(biāo)系Fig.1　The processing coordinate system of the modified helical pinion

產(chǎn)形輪(i=s)和傳動輪(i=l)動標(biāo)系.

圖1(c)中修形齒條刀具左右兩側(cè)齒面∑i和∑j(i,j=cs,cl)方程為(符號?和±均按左右順序表示)

(1)

動標(biāo)系Sa到Sc的轉(zhuǎn)換矩陣[11]為Mca，得Sc中修形齒條刀具齒面∑i(i=cs,cl)方程為

rcl(ui,li)=Mcaral(ui,li)，

(2)

其齒面法向量為ncl；Sc到Ss的轉(zhuǎn)換矩陣為Msc(ψs)；得Ss中∑i(i=cs)包絡(luò)的修形產(chǎn)形輪齒面∑s1方程為

(3)

其齒面法向量為 ns2；圖2為修形斜齒面齒輪加工坐標(biāo)系；Ss到S2轉(zhuǎn)化矩陣為M2s；φs和φ2滿足以下關(guān)系,其中Z2為面齒輪齒數(shù):

(4)

圖2　修形斜齒面齒輪加工坐標(biāo)系Fig.2　The processing coordinate system of the modified face gear with helical pinion

修形產(chǎn)形輪齒面∑s1包絡(luò)的修形斜齒面齒輪齒面∑1在動標(biāo)系S2中的方程為

(5)

用同樣方法可由修形齒條右側(cè)刀刃包絡(luò)得到修形產(chǎn)形輪齒面,進(jìn)而包絡(luò)求出修形斜齒面齒輪另一齒面的方程.

1.2修形直齒面齒輪齒面方程

與修形斜齒面齒輪不同的是,修形直齒面齒輪兩側(cè)齒面是對稱的,故只需推導(dǎo)一面的方程即可.此處所述修形直齒面齒輪的修形方法跟1.1節(jié)方法一致,推導(dǎo)過程也基本相同,不贅述其過程,在圖2所示加工坐標(biāo)系中可求得動標(biāo)系S2中修形面齒輪方程為

(6)

其中rs(us,ls)為修形直齒產(chǎn)形輪方程.

2　修形斜齒和直齒面齒輪高精度幾何建模

有限元仿真對幾何模型精度要求較高.表1和表2為修形斜齒和直齒面齒輪設(shè)計(jì)參數(shù).

表1修形斜齒面齒輪設(shè)計(jì)參數(shù)

Table 1The design parameters of the modified face gear with helical pinion

設(shè)計(jì)參數(shù)傳動小輪產(chǎn)形輪修行斜齒面齒輪齒數(shù)Z2528160壓力角α/(°)2525—模數(shù)m/mm6.356.356.35旋向左旋左旋右旋螺旋角β/(°)151515外徑/mm——582內(nèi)徑/mm——497軸交角γ/(°)909090修形系數(shù)as0.0020.003—齒頂系數(shù)1.01.251.0齒底系數(shù)1.251.251.25比值so/ωo0.90.9—

表2修形直齒面齒輪設(shè)計(jì)參數(shù)

Table 2The design parameters of the modified face gear with spur pinion

設(shè)計(jì)參數(shù)參數(shù)值小齒輪齒數(shù)Z119產(chǎn)形輪齒數(shù)Zs21面齒輪齒數(shù)Z273名義壓力角α/(°)25模數(shù)m/mm4.5面齒輪內(nèi)徑/mm158面齒輪外徑/mm189拋物線位置參數(shù)Uo/mm3.035軸交角γ/(°)90產(chǎn)形輪齒頂高系數(shù)as1.25產(chǎn)形輪齒根高系數(shù)bs1小輪齒頂高系數(shù)a11小輪設(shè)計(jì)齒根高系數(shù)bm1

根據(jù)修形斜齒和直齒面齒輪方程,利用MATLAB計(jì)算功能計(jì)算得到齒面離散點(diǎn),計(jì)算流程如圖3所示.

圖3　修形斜齒面齒輪齒面離散點(diǎn)計(jì)算流程圖Fig.3　The calculation flow chart of the tooth surface discrete points of the modified face gear with helical pinion

計(jì)算得到的齒面點(diǎn)數(shù)據(jù),導(dǎo)入CATIA中進(jìn)行擬合,經(jīng)分割、陣列、去除等步驟可得到高精度修形斜齒面齒輪幾何模型裝配圖，如圖4所示[12-13].同樣方法可得修形直齒面齒輪幾何模型裝配圖，如圖5所示.

圖4　修形斜齒面齒輪Fig.4　The modified face gear with helical pinion

圖5　修形直齒面齒輪Fig.5　The modified face gear with spur pinion

3　修形斜齒和直齒面齒輪功重比對比分析

功重比為承載功率和質(zhì)量的比值,因2組齒輪采用相同材料,質(zhì)量等價于體積.本文采用有限元分析的方法,對比相同功重比的修形斜齒面齒輪和修形直齒面齒輪的承載能力.在CATIA中測得修形斜齒面齒輪大小輪體積總和為22 900 000 mm3，修形直齒面齒輪大小輪體積總和為31 320 000 mm3.

3.1修形斜齒和直齒面齒輪有限元建模

基于ABAQUS有限元軟件[14-17],對圖4和圖5的2組面齒輪進(jìn)行有限元建模,其建模過程按以下步驟進(jìn)行：1)幾何模型導(dǎo)入；2)網(wǎng)格劃分；3)裝配；4)設(shè)置材料屬性：E=2.06×105MPa,v=0.3；5)設(shè)置分析步；6)建立接觸與約束；7)定義邊界條件和載荷；8)設(shè)置網(wǎng)格類型以及創(chuàng)建并提交作業(yè).

對修形斜齒面齒輪進(jìn)行網(wǎng)格劃分時,需先在CATIA中將修形斜齒面齒輪的單齒進(jìn)行剖分,過程如圖6所示.有些情況下,若齒輪內(nèi)外端粗細(xì)不均,還需在齒寬方向進(jìn)行剖分,把分塊同時導(dǎo)入ABAQUS進(jìn)行網(wǎng)格劃分,再把劃分好網(wǎng)格的分塊進(jìn)行接合.修形斜齒面齒輪有限元網(wǎng)格如圖7所示.

圖6　CATIA中剖分單齒Fig.6　The partition of single gear in CATIA

圖7　修形斜齒面齒輪有限元網(wǎng)格模型Fig.7　The finite element mesh model of the modified face gear with helical pinion

修形直齒面齒輪有限元建模過程相比修形斜齒面齒輪更加簡單,不具體介紹.

修形斜齒面齒輪轉(zhuǎn)速為5 000 r/min，功率為524 kW，計(jì)算得大輪扭矩為1 000 N·m.在相同轉(zhuǎn)速下,根據(jù)測得的體積,可求得同功重比的修形直齒面齒輪承載功率為716.7 kW，大輪扭矩為1 368 N·m.把已建好的2組齒輪有限元模型定義載荷后即可提交計(jì)算.

3.2接觸應(yīng)力對比

2組齒輪經(jīng)有限元分析后,得到同功重比情況下修形斜齒面齒輪最大接觸應(yīng)力如圖8所示,修形直齒面齒輪最大接觸應(yīng)力如圖9所示.

圖8　修形斜齒面齒輪最大接觸應(yīng)力云圖Fig.8　The maximum contact stress cloud chart of the modified face gear with helical pinion

圖9　修形直齒面齒輪最大接觸應(yīng)力云圖Fig.9　The maximum contact stress cloud chart of the modified face gear with spur pinion

由圖8和圖9可知修形斜齒面齒輪最大接觸應(yīng)力為1 241 MPa，修形直齒面齒輪最大接觸應(yīng)力為1 482 MPa.對比分析可知：同功重比工況下，修形斜齒面齒輪相比修形直齒面齒輪最大接觸應(yīng)力降低16.3%((1482-1241)/1482×100%).

3.3彎曲應(yīng)力對比

修形斜齒面齒輪最大彎曲應(yīng)力如圖10所示,修形直齒面齒輪最大彎曲應(yīng)力如圖11所示.

圖10　修形斜齒面齒輪最大彎曲應(yīng)力云圖Fig.10　The maximum bending stress cloud chart of the modified face gear with helical pinion

圖11　修形直齒面齒輪最大彎曲應(yīng)力云圖Fig.11　The maximum bending stress cloud chart of the modified face gear with spur pinion

由圖10可知修形斜齒面齒輪最大彎曲應(yīng)力為579.9 MPa，由圖11可知修形直齒面齒輪最大彎曲應(yīng)力為857.9 MPa.對比可知：相同功重比工況下,修形斜齒面齒輪相比修形直齒面齒輪最大彎曲應(yīng)力降低32.4%((857.9-579.9)/857.9×100%).

3.4重合度對比

修形斜齒面齒輪有限元仿真結(jié)果的齒對接觸力隨時間分布如圖12所示.

圖12　修形斜齒面齒輪接觸力Fig.12　The contact force diagram of the modified face gear with helical pinion

修形直齒面齒輪齒對接觸力隨時間分布如圖13所示.

圖13　修形直齒面齒輪接觸力Fig.13　The contact force diagram of the modified face gear with spur pinion

同樣方法算得同功重比的修形直齒面齒輪傳動重合度為1.636.

對比分析可知：在上述同功重比工況下,修形斜齒面齒輪相比修形直齒面齒輪重合度提高10.3%((1.804-1.636)/1.636×100%).

4　結(jié)　論

同功重比工況下,修形斜齒面齒輪相比修形直齒面齒輪：最大接觸應(yīng)力降低16.3%,最大彎曲應(yīng)力降低32.4%,重合度提高10.3%.修形斜齒面齒輪比修形直齒面齒輪承載能力更強(qiáng),重合度更高,更加符合高速重載工況下的輕量化設(shè)計(jì)要求.

[1] 王征兵. 大型偏航變槳齒輪箱輕量化設(shè)計(jì)技術(shù)研究[D]. 鄭州：機(jī)械科學(xué)研究總院鄭州機(jī)械研究所，2012：28-30．

WANG Zheng-bing. Research on lightweight design technology for large scale yaw and pitch gearboxes [D]. Zhengzhou: China Academy of Machinery Science & Technology， Zhengzhou Research Institute of Mechanical Engineering， 2012: 28-30．

[2] LITVIN F L,ZHANG Y,WANG J C,et al.Design and geometry of face-gear drives[J].Journal of Mechanical Design,1992,114(4):642-647．

[3] LITVIN F L.Development of face-gear technology for industrial and aerospace power transmission[R].Washington D.C.:NASA/CR-2002-211320．

[4] HEATH G F,BOSSLER R B.Advanced Rotorcraft Transmission (ART) Program-Final Report[R].Washington D.C.:NASA/CR-1993-191057.

[5] LITVIN F L, FUENTES A,ZANZI C.Design,generation and stress analysis of two versions of geometry of face-gear drives[J].Mechanism and Machine Theory,2002,37(10):1179-1211．

[6] LITVIN F L,GONZALWZ-PEREZ I,FUENTES A.Design,generation and stress analysis of face-gear drive with helical pinion[J].Computer Methods in Applied Mechanics and Engineering,2005,194(36/38):3870-3902．

[7] LITVIN F L,FUENTES A,HOWKINS M.Design,generation and TCA of new type of asymmetric face-gear drive with modified geometry[J].Computer Methods in Applied Mechanics and Engineering， 2001， 190(43/44):5837-5865．

[8] LITVIN F L ,VECCHIATO D,GUROVICH E.Computerized developments in design,generation,simulation of meshing,and stress analysis of gear drives[J].Meccanica,2005,40(3):291-323．

[9] 雷敦財(cái).面齒輪時變嚙合剛度計(jì)算及動態(tài)嚙合性能研究[D].長沙：中南大學(xué)機(jī)電工程學(xué)院,2013:42-46．

LEI Dun-cai.Study of time varying mesh stiffness and dynamic performance of face-gear pair [D].Changsha: Central South University,School of Mechanical and Electrical Engineering,2013:42-46．

[10] 李特文.齒輪幾何學(xué)與應(yīng)用理論[M].國凱,葉凌云,范琳,等,譯.上海：科學(xué)技術(shù)出版社,2008:484-516．

LITVIN F L.Gear geometry and applied theory [M].Translated by GUO Kai,YE Ling-yun,FAN Lin,et al.Shanghai:Science and Technology Press,2008:484-516．

[11] 李瑰賢.空間幾何建模及工程應(yīng)用[M].北京：高等教育出版社,2007:135-137．

LI Kui-xian.Spatial geometry modeling and its application in engineering[M].Beijing:Higher Education Press,2007:135-137．

[12] 丁仁亮.CATIA V5教程[M].北京：機(jī)械工業(yè)出版社,2007:2-30．

DING Ren-liang.CATIA V5 Course[M].Beijing:China Machine Press,2007:2-30．

[13] 陳興明.含安裝誤差的修形直齒面齒輪傳動齒面接觸分析的研究[D].長沙：中南大學(xué)機(jī)電工程學(xué)院,2012:21-24．

CHEN Xing-ming.Tooth contact analysis of modified spur face gear drives with alignment errors[D].Changsha:Central South University,School of Mechanical and Electrical Sngineering, 2012:21-24．

[14] 劉艷平.直齒面齒輪加載接觸分析及彎曲應(yīng)力和接觸應(yīng)力計(jì)算方法研究[D].長沙：中南大學(xué)機(jī)電工程學(xué)院,2012:29-32．

LIU Yan-ping.The loaded tooth contact analysis and bending stress and contact stress calculating method research of face gear drive with spur involute pinion[D].Changsha:Central South University,School of Mechanical and Electrical Engineering, 2012:29-32．

[15] 石亦平,周玉蓉.ABAQUS有限元分析實(shí)例詳解[M].北京：機(jī)械工業(yè)出版社,2007:24-32．

SHI Yi-ping,ZHOU Yu-rong.Detailed answers of instances by ABAQUS finite element analysis [M].Beijing:China Machine Press,2007,24-32．

[16] 羅希年.齒輪彎曲應(yīng)力計(jì)算的優(yōu)化[D].西安：長安大學(xué)機(jī)械工程學(xué)院,2012:42-47．

LUO Xi-nian.The optimization of calculation of gear bending stress[D].Xi′an: Chang′an University,School of Mechanical Engineering,2012:42-47．

[17] 戈紅霞,呂慶軍,張志凱.關(guān)于面齒輪接觸和彎曲應(yīng)力有限元計(jì)算方法的研究[J].新技術(shù) 新工藝,2014(1):43-47．

GE Hong-xia,Lü Qing-jun,ZHANG Zhi-kai.Research on finite element method of face gear contact and bending stress calculation[J].New Technology & New Process,2014(1):43-47．

[18] 董建雄,唐進(jìn)元.基于輪齒接觸分析的修形面齒輪傳動重合度計(jì)算研究[J].機(jī)械傳動,2015,39(4):8-10．

DONG Jian-xiong,TANG Jin-yuan.Study on contact ratio calculation of modified face gears based on contact analysis[J].Mechanical Transmission,2015,39(4):8-10．

Comparative study on the performance of modified face gears withhelical pinion and spur pinion in the case of same power density

ZHOU Zhen-yu, TANG Jin-yuan, DONG Jian-xiong

(State Key Laboratory of High Performance Complex Manufacturing, Central South University, Changsha 410083, China)

In order to know the performance difference of the modified face gears with helical pinion and spur pinion under the same power density better and choose a more suitable kind of face gear drive for high speed and heavy load conditions. Based on meshing theory, the tooth surface equations of the modified face gears with helical pinion and spur pinion were derived. The three-dimensional models of the modified face gears with helical pinion and spur pinion were built by CATIA. The finite element contact analysis method was used and the contact stress, bending stress and contact ratio were considered as performance indexes of the face gears to expand the research. Researches proved that the contact stress of modified face gear with helical pinion was much lower than that of modified face gear with spur pinion. The maximum contact stress was reduced by 16.3% in the case. The bending stress of modified face gear with helical pinion was much lower than that of modified face gear with spur pinion. The maximum bending stress was reduced by 32.4% in the case. The contact ratio of modified face gear with helical pinion was much higher than that of modified face gear with spur pinion. The contact ratio was increased by 10.3% in the case. So the transmission performance of the modified face gear with helical pinion is better than that of modified face gear with spur pinion under the same power density, the former is more suitable for the lightweight design under high speed and heavy load condition.

lightweight; power density; contact stress; bending stress; contact ratio; modification

2016-01-20.

本刊網(wǎng)址·在線期刊：http://www.journals.zju.edu.cn/gcsjxb

國家自然科學(xué)基金資助項(xiàng)目(51535012，51305462，51275530).

周鎮(zhèn)宇(1991—),男,江西宜春人,碩士生,從事齒輪傳動及數(shù)字化制造研究,E-mail:zhenyuzhou@csu.edu.cn.

10.3785/j.issn. 1006-754X.2016.04.007

TH 132.41

A

1006-754X(2016)04-0338-07

http://orcid.org//0000-0001-5888-2892

通信聯(lián)系人：唐進(jìn)元(1962-),男,湖南長沙人,教授,從事齒輪傳動及數(shù)字化制造研究,E-mail:jytangcsu@163.com.http://orcid.org//0000-0001-7186-1316