張西金,賈　超,候祥穎,郭　芳,翟　佳

(西北工業(yè)大學(xué) 機(jī)電學(xué)院，西安 710072)

斜齒輪具有重合度大,承載能力強(qiáng)等特點(diǎn)，在工業(yè)領(lǐng)域有廣泛的應(yīng)用.齒輪制造及安裝誤差不可避免，斜齒輪制造及安裝誤差會(huì)造成齒輪副軸線不平行，出現(xiàn)的軸線角誤差會(huì)使其傳動(dòng)性能變差，而齒輪修形可以改善其傳動(dòng)性能[1-2].因此,可以通過斜齒輪修形技術(shù)研究，以便減少軸線角誤差給其傳動(dòng)性能帶來的不利影響.

近年來，斜齒輪修形研究取得了一些進(jìn)展.文獻(xiàn)[3]提出一種新的拓?fù)湫扌锡X面，采用齒輪接觸分析(Tooth Contact Analysis,TCA)技術(shù)進(jìn)行了接觸路徑和傳動(dòng)誤差分析，并應(yīng)用有限元接觸分析方法研究了輪齒彎曲應(yīng)力和接觸應(yīng)力，但尚未研究修形齒輪的承載傳動(dòng)誤差(Loaded Transmission Errors,LTE)情況.文獻(xiàn)[2,4]推導(dǎo)了修形斜齒輪的齒面接觸與邊緣接觸的全部計(jì)算，探討了齒面修形對(duì)改善齒輪傳動(dòng)性能的作用，并在TCA基礎(chǔ)上提出了一種斜齒輪承載接觸分析(Loaded Tooth Contact Analysis,LTCA)(有限元方法和非線性數(shù)學(xué)規(guī)劃方法的LTCA)，對(duì)修形斜齒輪進(jìn)行了LTCA計(jì)算，探討了刀具修形及CNC機(jī)床修形兩種方法,對(duì)其他修形方法未做探討.文獻(xiàn)[5]研究了沿齒廓方向拋物線修形的齒面結(jié)構(gòu)，未涉及齒向修形.文獻(xiàn)[6-7]通過理論齒面疊加修形齒面，基于TCA和LTCA，建立動(dòng)力學(xué)模型，通過三維修形完成了寬斜齒輪多目標(biāo)齒面優(yōu)化，他也進(jìn)行了斜齒輪對(duì)角修形多目標(biāo)優(yōu)化設(shè)計(jì)以減小振動(dòng)與噪音，但都只進(jìn)行了標(biāo)準(zhǔn)安裝齒輪副的修形優(yōu)化.目前國內(nèi)外斜齒輪修形技術(shù)研究還不成熟，只有針對(duì)特定工況進(jìn)行拓?fù)湫扌尾拍塬@得較好的傳動(dòng)性能，因此斜齒輪修形方法還需要探索完善.

文中針對(duì)軸線角誤差的斜齒輪拓?fù)湫扌渭夹g(shù)進(jìn)行仿真研究，以減小軸線角誤差對(duì)其承載傳動(dòng)誤差(LTE)的影響，為斜齒輪修形設(shè)計(jì)提供了參考.

1　軸線角誤差的斜齒輪拓?fù)湫扌?/h2>

斜齒輪拓?fù)湫扌伟X廓修形和齒向修形，齒廓修形是通過改變齒輪刀具切削刃的方法實(shí)現(xiàn)的，齒向修形是通過數(shù)控機(jī)床的多軸聯(lián)動(dòng)實(shí)現(xiàn)的，這樣的修形方案可以提高修形效率.對(duì)于加工刀具齒條(滾刀)，可以將直線齒廓改變?yōu)閳A弧形、拋物線形等形式，由于圓弧齒廓修形研究較少，所以這里研究圓弧形齒廓修形，也就是在法平面內(nèi)刀具齒廓為圓弧形.為了減少修形工作量，采用小輪修形，大輪為標(biāo)準(zhǔn)漸開線齒輪的方案.建立如圖1所示的動(dòng)坐標(biāo)系Sc10和Sc1，它們都固定在小輪的加工刀具齒條上，Xc10Yc10為齒條齒寬中心處的法平面，采用圓弧形齒廓(如圖1所示)，其中R為圓弧半徑，R越大修形量越大；圓弧與原來的直線齒廓相切，切點(diǎn)距離節(jié)線與直線齒廓的交點(diǎn)為c，節(jié)線上下工作齒面的修形量應(yīng)該有所不同，設(shè)置參數(shù)c的目的就是為了減小靠近齒根部位工作齒面的修形量而增大靠近齒頂部位工作齒面的修形量；a表示節(jié)線上法向半齒厚.

圖1　修形刀具齒形

在圓弧形齒廓刀具展成加工齒輪的基礎(chǔ)上，齒向采用鼓形修形，中間段采用直線，兩側(cè)采用曲線，研究表明高階修形具有較好的傳動(dòng)性能，這里采用四階修形曲線.考慮軸線角誤差的影響，兩側(cè)的修形情況應(yīng)該有所不同，即采用不對(duì)稱的齒向修形方案.齒向修形如圖2所示.

圖2　齒向鼓形修形

在Sc10坐標(biāo)系中刀具齒面的四元位矢γc10表示為

(1)

式中：u1為Xc10坐標(biāo)值，v1為Zc10坐標(biāo)值.

齒向修形公式為

(2)

將γc10變換到Sc1坐標(biāo)系中需要經(jīng)過四次坐標(biāo)變換，每次的轉(zhuǎn)換矩陣可以由以下公式得出[8]：

最終的變換矩陣如下：

式中：α為壓力角;β為螺旋角.在Sc1坐標(biāo)系中齒條齒面的四元位矢和單位法矢表示為：

rc1=Dc1c10rc10

(3)

(4)

式(3)和式(4)完全表達(dá)了齒條刀具的空間曲面情況，刀具齒面在加工齒輪時(shí)按展成原理就可以加工出齒輪齒面，由于刀具齒廓已不再是直線形式，因此加工出來的齒輪齒廓也不再是標(biāo)準(zhǔn)的漸開線齒廓.

2　拓?fù)湫扌涡饼X輪齒面接觸分析

2.1　修形齒面建模

拓?fù)湫扌锡X面建立采用了通過理論齒面與法向修形曲面疊加構(gòu)造修形齒面的方法[6]，理論齒面采用刀具展成加工形成的曲面，而法向修形曲面是考慮齒向鼓形修形產(chǎn)生的修形曲面，法向修形量為式(2)所示.建立如圖3所示的展成加工坐標(biāo)系，Sa1為固定坐標(biāo)系，S1為小齒輪的動(dòng)坐標(biāo)系，Za1和Z1都是和齒輪的回轉(zhuǎn)中心重合.S1隨小齒輪一起繞Za1軸轉(zhuǎn)動(dòng)，刀具沿Yc1軸平移，刀具齒面展成過程中形成齒輪齒面，因此小齒輪齒面位矢和單位法矢可以表示為

r1=M(φ1)1c1rc1

(5)

n1=M(φ1)1c1nc1

(6)

其中M(φ1)1c1是刀具坐標(biāo)系Sc1到小輪坐標(biāo)系S1的轉(zhuǎn)換矩陣，展成過程中滿足關(guān)系式H=Rp1×φ1，Rp1為節(jié)圓半徑，H為位移，φ1為齒輪轉(zhuǎn)角.

根據(jù)微分幾何中的嚙合方程可以得出以下關(guān)系

(7)

將式(7)代入式(5)和式(6)可以確定小輪理論齒面情況，也就是可以得到刀具齒面上任何一點(diǎn)在小齒輪齒面上的對(duì)應(yīng)點(diǎn).

實(shí)際小輪齒面是由理論面與法向修形曲面疊加形成的，其位矢和法矢[6]為

r1r=δ(u1,v1)n1+r1

(8)

(9)

大輪采用標(biāo)準(zhǔn)齒條刀具展成加工，齒面建模過程在這里不再累述.

圖3　齒輪展成示意圖

2.2　齒面分析

輪齒承載接觸分析技術(shù)是在載荷作用下對(duì)齒輪的嚙合過程進(jìn)行高精度的計(jì)算機(jī)仿真，以獲得設(shè)計(jì)齒輪在擬真實(shí)工況條件下的工作性能的一個(gè)重要手段，LTCA計(jì)算得到承載傳動(dòng)誤差是分析齒輪動(dòng)態(tài)性能(振動(dòng)、噪聲)的一個(gè)參數(shù)，LTE波動(dòng)越小齒輪嚙合動(dòng)態(tài)性能一般越好.這里采用LTCA進(jìn)行斜齒輪拓?fù)湫扌畏治?，而LTCA是建立在TCA基礎(chǔ)上的，因此需要同時(shí)研究TCA技術(shù).

建立如圖4所示的TCA分析坐標(biāo)系，Sa1和Sa2是輔助坐標(biāo)系，不隨齒輪轉(zhuǎn)動(dòng)且互相平行，用于表示軸線角誤差.其中固定坐標(biāo)系Ss的Zs軸與大輪的回轉(zhuǎn)軸重合.在固定坐標(biāo)系中，大小輪配對(duì)齒面在嚙合傳動(dòng)時(shí)嚙合點(diǎn)位矢及單位法矢必定相同，由此可得到如下TCA分析方程.

M(φ1,γ1,γ2)s1φ=M(φ2)s2r2

(10)

M(φ1,γ1,γ2)s1φ=M(φ2)s2n2

(11)

式(10)～(11)中M(φ1,γ1,γ2)s1和M(φ2)s2分別表示小輪和大輪坐標(biāo)系到固定坐標(biāo)系的變換矩陣，φ1和φ2是大小輪的瞬時(shí)轉(zhuǎn)角，γ1和γ2是不同方向的軸線角誤差.由于n1r和n2為單位法矢，所以構(gòu)成了5個(gè)獨(dú)立方程，結(jié)合大小輪展成過程中的嚙合方程可以得到嚙合點(diǎn).進(jìn)一步可以獲得各接觸橢圓的長(zhǎng)軸及齒面間隙.

圖4　TCA坐標(biāo)系

式(12)中F為接觸輪齒的法向柔度矩陣；P為傳遞力矩計(jì)算得到的齒面法向力；p和d分別為齒對(duì)沿接觸橢圓長(zhǎng)軸離散點(diǎn)構(gòu)成的法向載荷向量及變形后的齒面間隙向量，w為沿接觸橢圓長(zhǎng)軸離散點(diǎn)變形前的齒面間隙向量，Z為輪齒的法向位移向量，pj和dj為沿接觸橢圓長(zhǎng)軸的離散法向載荷及變形后的齒面間隙.

在TCA基礎(chǔ)上進(jìn)行承載接觸(LTCA)分析，接觸分析方程為

(12)

在此基礎(chǔ)上可以得到承載接觸分析數(shù)學(xué)規(guī)劃模型[9]，F(xiàn)采用有限元計(jì)算獲得，w由TCA計(jì)算獲得，P由傳遞力矩計(jì)算得到，由此求解承載接觸分析數(shù)學(xué)規(guī)劃模型獲得參數(shù)p、d、Z.將Z轉(zhuǎn)化為大輪轉(zhuǎn)角，即可得到承載傳動(dòng)誤差.

3　斜齒輪齒廓修形參數(shù)

斜齒輪齒廓拓?fù)湫扌伟X廓修形和齒向修形，首先研究齒廓修形問題.齒輪的基本參數(shù)為:小輪齒數(shù)30，大輪齒數(shù)45，法面模數(shù)6，壓力角20°，螺旋角15°，傳遞力矩750 N·m-1.設(shè)計(jì)圓弧半徑R=K1×mn，mn為法向模數(shù)，為了不影響接觸強(qiáng)度，圓弧半徑不應(yīng)太小.設(shè)計(jì)c=K2×mn,其中K2=0時(shí)圓弧齒廓切點(diǎn)在齒條刀具節(jié)線位置，K2越大表示圓弧齒廓切點(diǎn)離齒條刀具節(jié)線越遠(yuǎn).

通過仿真發(fā)現(xiàn)：當(dāng)圓弧半徑R過小時(shí)會(huì)產(chǎn)生齒側(cè)邊緣接觸，而且傳動(dòng)誤差和齒面接觸應(yīng)力都比較大；當(dāng)圓弧半徑R太大時(shí)，齒廓修形量太小時(shí)，又會(huì)造成齒頂邊緣接觸.綜合考慮取R=250mn，且取c=0，使鼓形中間區(qū)的大小輪嚙合點(diǎn)都在分度圓上.

4　軸線角誤差的斜齒輪齒向修形

實(shí)際齒輪安裝時(shí)一般都會(huì)存在軸線角誤差，由于不同軸線角誤差需要的修形是不同的，所以需要研究不同軸線角誤差下的不對(duì)稱的修形問題，以便齒輪副在實(shí)際安裝工況下取得較好的傳動(dòng)性能.鑒于齒向修形參數(shù)對(duì)承載傳動(dòng)誤差影響較大，所以這里齒廓修形參數(shù)取前面的值，只研究齒向修形的修形量和修形長(zhǎng)度對(duì)傳動(dòng)誤差的影響，從而確定合理的修形參數(shù).

對(duì)于高精度齒輪ISO推薦齒向鼓形最大修形量為：10～25 μm，這里我們?nèi)∽畲笮扌瘟繛?0 μm和25 μm兩種情況進(jìn)行研究.

4.1　軸線角誤差γ1的修形

以Ya2Za2投影面內(nèi)軸線角誤差γ1取0.05′為例，b為齒寬，左右兩端最大修形長(zhǎng)度分別取0.1b、0.2b、0.3b、0.4b和0.5b，左右兩端最大修形量分別取(δLmax=0.010 mm,δRmax=0.010 mm)、(δvLmax=0.025 mm,δRmax=0.025 mm)、(δLmax=0.010 mm,δRmax=0.025 mm)和(δLmax=0.025 mm,δRmax=0.010 mm)四種情況進(jìn)行LTCA仿真，其中一種情況的仿真結(jié)果如圖5所示，四種情況的最小傳動(dòng)誤差波動(dòng)見表1.

圖5　軸線角誤差γ1為0.05′時(shí)LTE波動(dòng)量(δLmax =0.010,δRmax =0.025)

從圖5可以看出，TLE波動(dòng)量隨最大修形長(zhǎng)度及最大修形量變化關(guān)系比較復(fù)雜，但左右最大修形長(zhǎng)度都為0.5b時(shí)TLE波動(dòng)量最大，這是最不利的傳動(dòng)情況.

表1 　γ1為0.05′時(shí)LTE波動(dòng)量

從表1可以看出，四種情況所對(duì)應(yīng)的最小LTE波動(dòng)量相差很小，最小LTE波動(dòng)量對(duì)應(yīng)的最大修形長(zhǎng)度相同(左端最大修形長(zhǎng)度0.3b、右端最大修形長(zhǎng)度0.2b)，說明取最大修形長(zhǎng)度不對(duì)稱的設(shè)計(jì)可以取得最好的傳動(dòng)性能.

根據(jù)表1確定合理的修形參數(shù)為：左端最大修形長(zhǎng)度0.3b、右端最大修形長(zhǎng)度0.2b、左右兩端最大修形量均為0.025 mm，此時(shí)承載傳動(dòng)誤差最大波動(dòng)量為0.291 8，與未修形相比承載傳動(dòng)誤差最大波動(dòng)量減少了42.85%(未修形時(shí)LTE最大波動(dòng)量為0.510 6)，修形結(jié)果表明該修形方案可以提高傳動(dòng)性能.該組參數(shù)仿真得到傳動(dòng)誤差及載荷分配系數(shù)如圖6所示，從圖6中可以看出LTE波動(dòng)小，主要承載區(qū)域在齒面中間，兩側(cè)區(qū)域沒有承載，說明還有富裕的承載能力，可以傳遞更大的力矩.

4.2　軸線角誤差γ2的修形

以Xa2Za2投影面內(nèi)軸線角誤差γ2取0.05′為例，左右兩端最大修形長(zhǎng)度分別取0.1b、0.2b、0.3b、0.4b和0.5b，左右兩端最大修形量同樣采用情況進(jìn)行LTCA仿真，其中一種情況的仿真結(jié)果如圖7所示，四種情況的最小傳動(dòng)誤差波動(dòng)見表2.

圖6　軸線角誤差γ1為0.05′時(shí)齒輪嚙合性能

圖7　軸線角誤差γ2為0.05′時(shí)LTE波動(dòng)量

從圖7可以看出，TLE波動(dòng)量隨最大修形長(zhǎng)度及最大修形量變化關(guān)系比較復(fù)雜，但左右修形長(zhǎng)度都為0.5b時(shí)TLE波動(dòng)量最大，這是最不利的傳動(dòng)情況.從表2可以看出，四種情況所對(duì)應(yīng)的最小LTE波動(dòng)量各不相同，取最大修形量不對(duì)稱(左端最大修形量0.010 mm、右端最大修形量0.025 mm)的設(shè)計(jì)可以取得最好的傳動(dòng)性能.

表2　γ2為0.05′時(shí)LTE波動(dòng)量

根據(jù)表2確定合理的修形參數(shù)為：左右兩端最大修形長(zhǎng)度均為0.3b、左端最大修形量0.010 mm、右端最大修形量0.025 mm，此時(shí)承載傳動(dòng)誤差最大波動(dòng)量為0.282 0，與未修形相比承載傳動(dòng)誤差最大波動(dòng)量減少了53.47%(未修形時(shí)LTE最大波動(dòng)量為0.606 1)，修形結(jié)果表明該修形方案可以提高傳動(dòng)性能.該組參數(shù)仿真得到傳動(dòng)誤差及載荷分配系數(shù)如圖8所示.從圖8可以看出LTE波動(dòng)小，主要承載區(qū)域在齒面中間，兩側(cè)區(qū)域沒有承載，可以傳遞更大的力矩.

圖8　軸線角誤差γ2為0.05′時(shí)嚙合性能

5　結(jié) 論

1) 提出了軸線角誤差的斜齒輪拓?fù)湫扌涡路桨?，用圓弧齒廓刀具展成加工齒輪，齒向采用高階非對(duì)稱鼓形修形.

2) 當(dāng)軸線角誤差γ1為0.05′時(shí)采用最大修形長(zhǎng)度不對(duì)稱的設(shè)計(jì)承載傳動(dòng)誤差最大波動(dòng)量減少了42.85%，當(dāng)軸線角誤差γ2為0.05′時(shí)采用最大修形量不對(duì)稱的設(shè)計(jì)承載傳動(dòng)誤差最大波動(dòng)量減少了53.47%.存在軸線角誤差時(shí)采用合理的齒向不對(duì)稱修形參數(shù)能有效減少承載傳動(dòng)誤差波動(dòng)量，提高齒輪副的嚙合性能.

3) 下一步將研究修形參數(shù)的優(yōu)化設(shè)計(jì)問題，使存在軸線角誤差的斜齒輪獲得最佳的嚙合性能.

參考文獻(xiàn)：

[1]王統(tǒng),賈毅,邱良恒,等.漸開線齒輪修形方法的進(jìn)展[J].上海交通大學(xué)學(xué)報(bào),1998,32(5):133.

WANG Tong,JIA Yi,QIU Liangheng,et al.Progress Outline of Gear Tooth Profile Modification[J].Journal of Shanghai Jiaotong University,1998,32(5):133.(in Chinese)

[2]方宗德.修形斜齒輪的輪齒接觸分析[J].航空動(dòng)力學(xué)報(bào),1997,12(3):247.

FANG Zongde.Tooth Contact Analysis of Modified Helical Gears[J].Journal of Aerospace Power,1997,12(3):247.(in Chinese)

[3]LITVIN F L,GONZALEZ-PEREZ I,FUENTES A,et al.Topology of Modified Surfaces of Involute Helical Gears with Line Contact Developed for Improvement of Bearing Contact,Reduction of Transmission Errors,and Stress Analysis[J].Mathematical and Computer Modeling,2005,42(9):1063.

[4]方宗德.修形斜齒輪的承載接觸分析[J].航空動(dòng)力學(xué)報(bào),1997,12(3):251.

FANG Zongde.Loaded Tooth Contact Analysis of Modified Helical Gears[J].Journal of Aerospace Power,1997,12(3):251.(in Chinese)

[5]盛鋼,沈云波.基于拋物線修形的斜齒輪傳動(dòng)嚙合特性研究[J].西安工業(yè)大學(xué)學(xué)報(bào),2007,27(6):531.

SHENG Gang,SHEN Yunbo.The Meshing of Modified Helical Gear Drive in Profile Direction[J].Journal of Xi’an Technological University,2007,27(6):531.(in Chinese)

[6]蔣進(jìn)科,方宗德,蘇進(jìn)展.寬斜齒輪多目標(biāo)修形優(yōu)化設(shè)計(jì)[J].西安交通大學(xué)學(xué)報(bào),2014,48(8):91.

JIANG Jinke,FANG Zongde,SU Jinzhan.Multi-Objective Optinal and Modified Design for Wide Helical Gear[J].Journal of Xi’an Jiaotong University,2014,48(8):91.(in Chinese)

[7]蔣進(jìn)科,方宗德,王峰.降低斜齒輪噪聲的對(duì)角修形優(yōu)化設(shè)計(jì)[J].振動(dòng)與沖擊,2014,33(7):63.

JIANG Jinke,FANG Zongde,WANG Feng.Optimal Design with Diagonal of Modification for Reducing Helical Gear Noise[J].Journal of Vibration and Shock,2014,33(7):63.(in Chinese)

[8]FAYDOR L.LITVIN,ALFONSO FUENTES.Gear Geometry and Applied Theory[M].Cambridge:Cambridge University Press,2004.

[9]方宗德,劉更，沈允文.圓柱齒輪傳動(dòng)的三維接觸分析[J].西北工業(yè)大學(xué)學(xué)報(bào),1992,10(1):74.

FANG Zongde,LIU Geng,SHEN Yunwen.3-D Tooth Contact Analysis of Helical Gears[J].Journal of Northwestern Polytechnical University,1992,10(1):74.

(in Chinese)