面齒輪輪齒剛度計(jì)算方法研究

黃鵬，李政民卿

(南京航空航天大學(xué) 江蘇省精密與微細(xì)制造技術(shù)重點(diǎn)實(shí)驗(yàn)室, 江蘇 南京 210016)

?

面齒輪輪齒剛度計(jì)算方法研究

黃鵬，李政民卿

(南京航空航天大學(xué) 江蘇省精密與微細(xì)制造技術(shù)重點(diǎn)實(shí)驗(yàn)室, 江蘇 南京 210016)

摘要：基于Buckingham的觀點(diǎn)，將面齒輪齒形看作是由沿齒長(zhǎng)方向一系列變壓力角的齒條組成，得到沿軸向和徑向都為變截面的面齒輪簡(jiǎn)化齒形，獲得了面齒輪輪齒嚙合變形的計(jì)算公式，求解出了面齒輪的輪齒剛度。與有元法的分析結(jié)果進(jìn)行比較，二者相對(duì)誤差較小，驗(yàn)證了面齒輪輪齒剛度計(jì)算方法的可行性。

關(guān)鍵詞：面齒輪；輪齒變形；輪齒剛度；計(jì)算方法

面齒輪傳動(dòng)是圓柱齒輪與面齒輪相嚙合的齒輪傳動(dòng)，具有重合度大、無(wú)需防錯(cuò)位設(shè)計(jì)、與面齒輪嚙合的直齒圓柱齒輪無(wú)軸向力作用和傳遞動(dòng)力性能優(yōu)越等優(yōu)點(diǎn)[1,3]；據(jù)波音公司報(bào)道[4]，采用面齒輪傳動(dòng)作為主減速器的阿帕奇Brock-Ⅲ型直升機(jī)已于2009年11月23日首飛成功。

在面齒輪傳動(dòng)中，由面輪齒輪齒嚙合變形引起的剛度激勵(lì)是其傳動(dòng)中的主要?jiǎng)討B(tài)激勵(lì)之一，面齒輪傳動(dòng)過(guò)程中的綜合嚙合剛度是研究動(dòng)力學(xué)的基礎(chǔ)，而確定面齒輪輪齒剛度是研究面齒輪傳動(dòng)綜合嚙合剛度的一項(xiàng)重要工作，因此有必要對(duì)面齒輪輪齒剛度的計(jì)算方法開(kāi)展研究。

圓柱齒輪和圓錐齒輪嚙合剛度的計(jì)算和測(cè)量方法從20世紀(jì)就開(kāi)始研究[5,6]，并取得了一定的研究成果。Isa Yesilyurt等人[7]對(duì)圓柱齒輪的單齒嚙合剛度的理論計(jì)算方法進(jìn)行了研究，并采用模態(tài)法來(lái)測(cè)量圓柱齒輪的單齒剛度。而目前對(duì)于面齒輪輪齒剛度多采用有限元法計(jì)算得到，其計(jì)算效率低。本文對(duì)面齒輪輪齒剛度的解析計(jì)算方法進(jìn)行研究，尋求一種高效準(zhǔn)確的面齒輪輪齒剛度計(jì)算方法。

1 面齒輪齒形分析

基于Buckingham的觀點(diǎn)[8]：面齒輪齒形可以看作是沿齒長(zhǎng)方向一系列變壓力角的齒條組成，如圖1(a)所示；其沿齒長(zhǎng)方向上某一截面的齒形分析，如圖1(b)所示。

圖1　面齒輪齒形及齒形分析圖

根據(jù)圖1(b)中的幾何關(guān)系，可求出面齒輪嚙合半徑r2處的齒厚t[8]為：

t=2r1[(T/2R1)+invφ1-invφ]

(1)

式中：r1為漸開(kāi)線直齒輪嚙合半徑；T為漸開(kāi)線直齒輪在分度圓處的齒槽寬，其值為：

T=πm/2

其中：m為模數(shù)；R1為漸開(kāi)線直齒輪分度圓半徑；φ為漸開(kāi)線直齒輪壓力角；φ1為節(jié)點(diǎn)壓力角，其值為：

φ1=arccos(Rb/r1)

其中Rb為漸開(kāi)線直齒輪基圓半徑。

齒條的齒頂高h(yuǎn)a為：

ha=r1-Rr-c

(2)

式中：Rr為漸開(kāi)線直齒輪齒根圓半徑；c為齒頂間隙。

齒條的齒根高h(yuǎn)f為：

hf=Roc-r1

(3)

式中Roc為漸開(kāi)線直齒輪齒頂圓半徑。

根據(jù)上述面齒輪齒形分析方法，以表1所列面齒輪參數(shù)在matlab中繪制無(wú)過(guò)渡圓弧的面齒輪單齒模型，如圖2所示。

表1　面齒輪傳動(dòng)的齒輪副主要幾何設(shè)計(jì)參數(shù)

圖2　面齒輪單齒模型仿真圖

2面齒輪輪齒剛度計(jì)算方法

基于面齒輪的復(fù)雜齒形，在求解面齒輪輪齒剛度時(shí)需要對(duì)面齒輪齒形進(jìn)行合理簡(jiǎn)化。根據(jù)仿真圖2，可將面齒輪齒形簡(jiǎn)化成由一系列沿齒高方向變化的梯形截面組成，如圖3所示。

圖3　面齒輪單齒模型簡(jiǎn)化示意圖

根據(jù)圖3中面齒輪齒形的簡(jiǎn)化形式，其輪齒剛度的解析法求解思想仍然基于圓柱齒輪輪齒剛度的當(dāng)量齒輪解析分析方法，即面齒輪的輪齒綜合彈性變形包括輪齒的彎曲變形、剪切變形、接觸變形以及由齒根基礎(chǔ)彈性傾斜引起的附加變形[9-11]。分析中，先通過(guò)求解輪齒各單一變形后，再將其相加，得出面齒輪輪齒在單位載荷作用下接觸點(diǎn)處的總變形量δf為：

δf=δ1+δ2+δ3+δ4

(4)

式中：δ1為彎曲變形；δ2為剪切變形；δ3為輪齒接觸變形；δ4為齒根基礎(chǔ)彈性傾斜引起的附加變形。

根據(jù)剛度和變形量-柔度的關(guān)系可知，面齒輪在接觸點(diǎn)處的輪齒剛度kf為:

kf=1/δf

(5)

2.1　面齒輪輪齒的彎曲變形和剪切變形

面齒輪輪齒的彎曲變形和剪切變形分析中，將在接觸范圍內(nèi)的輪齒分成若干小段i，并把小段i看成一個(gè)懸臂梁，如圖4(a)所示；每個(gè)小段i沿齒長(zhǎng)方向的截面可看作是隨面齒輪齒高方向變化的梯形截面，如圖4(b)所示。

圖4　面齒輪輪齒簡(jiǎn)化成的懸臂梁模型

根據(jù)截面上切應(yīng)力的分布規(guī)律作以下2個(gè)假設(shè)[12]：

1) 橫截面上各點(diǎn)切應(yīng)力的方向都平行于剪應(yīng)力的方向；

2) 切應(yīng)力沿截面寬度均勻分布。

根據(jù)上述假設(shè)，面齒輪輪齒在接觸點(diǎn)處的彎曲變形和剪切變形可通過(guò)對(duì)各小段的變形量求和獲得。

面齒輪輪齒在接觸點(diǎn)處彎曲變形包括單位長(zhǎng)度載荷作用在小段i上產(chǎn)生的彎矩變形和橫向變形。

根據(jù)圖4(a)所示，取小段i的厚度為ti，截面面積為Ai，截面模量為Ii；小段i至載荷作用點(diǎn)沿x方向的距離為h，載荷作用點(diǎn)處的半齒厚為yi，齒根處的半齒厚為ym；泊松比為v，彈性模量為E；面齒輪齒寬為b；載荷與y軸間的夾角為βj；則單位長(zhǎng)度載荷Fj作用在小段i上的彎矩變形δ11為：

(6)

由Fj對(duì)小段i產(chǎn)生的橫向變形δ12為：

(7)

因此，單位載荷作用下小段i的彎曲變形δbi為：

δbi=δ11+δ12

(8)

由Fj對(duì)小段i產(chǎn)生的剪切變形δsi為：

(9)

根據(jù)圖4(b)可知，由于各小段的截面均是隨面齒輪齒高位置變化的梯形截面，因此其截面模量Ii為：

(10)

截面面積Ai為：

Ai=b(Sn+Sw)/2

(11)

式中，沿齒高方向任意高度h處的Sn、Sw可采用圖5所示模型計(jì)算。

圖5　Sn，Sw的計(jì)算模型

根據(jù)圖5可知，利用梯形ABCD可求得:

Sn=d1+2(h-a1)tanφmin

(12)

式中：φmin為面齒輪內(nèi)半徑節(jié)線處的壓力角；d1為面齒輪內(nèi)半徑節(jié)線處的齒厚；a1為面齒輪內(nèi)半徑節(jié)線位置到齒頂?shù)木嚯x。d1和a1的計(jì)算式分別為：

d1=2r11(T/2R1-invφmin+invφ)

a1=r11-Rr-c

式中r11為面齒輪內(nèi)半徑。

利用梯形EFGH可求得:

Sw=d2+2(h-a2)tanφmax

(13)

式中：φmax為面齒輪外半徑節(jié)線處的壓力角；d2為面齒輪外半徑節(jié)線處的齒厚；a2為面齒輪外半徑節(jié)線位置到齒頂?shù)木嚯x。d2和a2的計(jì)算式分別為:

d2=2r12(T/2R1-invφmax+invφ)

a2=r12-Rr-c

式中r12為面齒輪外半徑。

利用式(8)-式(13)，可分別求出面齒輪輪齒在載荷作用點(diǎn)處沿載荷方向的彎曲變形和剪切變形，其中彎曲變形δ1為:

(14)

剪切變形δ2為:

(15)

式中n為小段i的個(gè)數(shù)。

2.2　面齒輪輪齒接觸變形和齒根基礎(chǔ)彈性傾斜引起的附加變形

面齒輪接觸變形分析是以布希涅斯克問(wèn)題的解為基礎(chǔ)，根據(jù)物體表面上有相同距離的點(diǎn)在公切面上形成類橢圓區(qū)域的結(jié)論，得出面齒輪輪齒的接觸變形δ3[13]為：

(16)

式中：J為橢圓積分系數(shù)，可通過(guò)查橢圓積分系數(shù)表獲得；ρx為接觸區(qū)域橢圓長(zhǎng)半徑，θg和θf(wàn)可由下式確定：

(17)

其中：i分別為g和f；Si為材料縱向延伸和橫向壓縮比的系數(shù)；E為材料的彈性模量。

面齒輪齒根基礎(chǔ)彈性傾斜引起的附加變形δ4[14-15]為：

(18)

式中：Lf=xj-xm-yjtanβj，Hf=2ym。

3輪齒剛度解析法的可行性驗(yàn)證

為驗(yàn)證上述面齒輪輪齒剛度計(jì)算方法的可行性，將解析法計(jì)算結(jié)果與有限元法計(jì)算結(jié)果進(jìn)行對(duì)比分析。根據(jù)表1中所列面齒輪參數(shù)，建立面齒輪有限元分析模型，如圖6所示。

圖6　面齒輪有限元分析模型

利用圖6中模型，將計(jì)算得到的變形量轉(zhuǎn)換為輪齒剛度并與上述面齒輪輪齒剛度解析計(jì)算方法進(jìn)行比較，如表2所示。

表2　解析法與有限元法的剛度計(jì)算比較表

采用有限元法和本文介紹的方法計(jì)算得到的面齒輪輪齒變形量及輪齒剛度沿齒高方向變化的曲線分別如圖7(a)和圖7(b)所示。

圖7　解析法與有限元法計(jì)算結(jié)果對(duì)比

根據(jù)圖7可知，解析法計(jì)算的面齒輪輪齒嚙合變形及輪齒剛度與有限元法分析所得結(jié)果相差較小。為進(jìn)一步分析面齒輪剛度解析法的可行性，定義解析法與有限元法的相對(duì)平均偏差η為:

(19)

式(19)中，N為接觸點(diǎn)個(gè)數(shù)。利用式(19)對(duì)表2中數(shù)據(jù)分析可知，解析法與有限元法的相對(duì)平均偏差為10.32%。因此，上述面齒輪輪齒剛度的解析計(jì)算方法是可行的。

4結(jié)論

1) 基于Buckingham的觀點(diǎn)，將面齒輪齒形看成是由一系列變壓力角齒條組成，形成了用于求解面齒輪輪齒剛度的齒形建模方法。

2) 將無(wú)過(guò)渡圓弧的面齒輪齒形看作是沿軸向和徑向均為變截面的梯形彈性懸臂梁，運(yùn)用材料力學(xué)方法，形成了面齒輪輪齒剛度的解析計(jì)算方法。

3) 分別運(yùn)用解析法和有限元法計(jì)算出了面齒輪輪齒剛度，得到了輪齒剛度沿齒高方向的變化規(guī)律，兩種計(jì)算方法得到的結(jié)果相差很小，證明了面齒輪輪齒剛度的計(jì)算方法的可行性。

參考文獻(xiàn):

[1] LITVIN F L, WANG J C, BOSSLER R B Application of face gear drives in helicopter transmissions [J]. Journal of Mechanical Design, Transactions of the ASME, 1994, 116 (3) : 672～676.

[2] HANDSCHUH R F, LEWICKI D G, BOSSLER R B. Experimental testing of prototype face gears for helicopter transmissions [J]. Journal of Aerospace Engineering, 1994, 208 (2) : 129～136.

[3] BILL R C. Summary highlights of the advanced rotorcraft transmission (ART) program [R]. American Inst of Aeronautics and Astronautics, New York, 1992.

[4] MESA Ariz. Boeing AH-64D Apache Block III structures test aircraft makes 1st flight [EB/OL]. [2009-12-4]. http://boeing.mediaroom.com/index.php?s=20295&item=977.

[5] MONSAK P, KAZEM K. Efficient evaluation of spur gear tooth mesh load using pseudo-interference stiffness estimation method [J]. Mechanism and Machine Theory. 2002, 37(8): 769~786.

[6] COMELL R W. Compliance and stress sensitivity of spur gear teeth [J]. Journal of Mechanical Design, 1981,103: 447~459.

[7] ISA Y, GU F, ANDREW D B. Gear tooth stiffness reduction measurement using modal analysis and its use in wear fault severity assessment of spur gears [J]. NDT&E International 2003,36 (5) :357~372.

[8] BUCKINGHAM E. Analytical mechanics of gears [M]. New York, Dover Publications, 1949: 302~320.

[9] BAND R V，PETERSON R E. Load and stress cycle in gear teeth [J]．Mechanical Engineering，1929，51(9)：653～662.

[10] WEBER C．The deformation of loaded gears and the effect on their load carrying capacity [J]．Sponsored Research (Germany) ，British Dept．of Scientific and Industrial Research，Report，1949 (3).

[11] ATTIA A Y． Deflection of spur gear teeth cut in thin rims [J]．ASME Eng. Ind，1964，86(4)： 333~341.

[12] 劉鴻文．材料力學(xué)[M]. 北京：高等教育出版社，2004.

[13] 李政民卿，朱如鵬．面齒輪傳動(dòng)的承載接觸分析[J]．南京航空航天大學(xué)學(xué)報(bào)，2010，42(2)：219~223.

[14] 吳俊飛，郭建章．變厚齒輪輪齒嚙合綜合剛度確定方法研究[J]．青島化工學(xué)院學(xué)報(bào)，2002，23(1)：74~76.

[15] 李潤(rùn)方，王建軍.齒輪系統(tǒng)動(dòng)力學(xué)-振動(dòng)、沖擊、噪聲[M]．北京：科學(xué)出版社，1996.

Research on Method of Calculating Stiffness of Face Gear Tooth

HUANG Peng，LI ZHENG min-qing

( Jiangsu Key Laboratory of Precision and Micro-Manufacturing Technology,

Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China)

Abstract:Based on Buckingham’s standpoint, along the length of the tooth of face gear , any one of its sections can be considered as a section of rack with varying pressure angle. According to above theory, the tooth form of the face gear is simulated successfully. The calculating method of the tooth’s stiffness of the face gear is obtained. Compared with the result of analyzing the tooth’s stiffness by FEM, both relative errors are almost same. it is proved that the calculation method for its stiffness is correct.

Keywords:face gear; gear tooth deformation; gear tooth stiffness; calculating method

中圖分類號(hào)：TH132.429

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1671-5276(2015)02-0046-04

作者簡(jiǎn)介：黃鵬(1987-)，男，湖北隨州人，碩士研究生，主要從事機(jī)械傳動(dòng)與減振降噪研究。

基金項(xiàng)目：國(guó)家自然科學(xué)基金(51105194)；國(guó)家留學(xué)基金資助(201206835011)

收稿日期：2014-11-08