王建，羅善明，蘇德瑜

(1. 廈門理工學(xué)院 機械與汽車工程學(xué)院，福建 廈門，361024；2. 湖南科技大學(xué) 機電工程學(xué)院，湖南 湘潭，411201)

重合度是衡量齒輪副的傳動連續(xù)性、平穩(wěn)性、傳遞載荷均勻性的重要度量指標，其定義是齒輪副的一對輪齒從開始嚙合到終止嚙合過程中，某一齒輪所轉(zhuǎn)過的角度與該齒輪單個齒距圓心角的比值。它所表達的是齒輪副的平均嚙合齒對數(shù)的信息[1-2]。齒輪傳動的重合度不僅與傳動的回轉(zhuǎn)誤差有關(guān)，還對齒輪強度有著重要的影響，而且為了保證能進行連續(xù)傳動，必須滿足重合度大于1 的條件。為了提高齒輪傳動的承載能力和動態(tài)性能，目前采用重合度大于2 的高重合度齒輪傳動[3-5]。與常規(guī)的齒輪傳動相比，高重合度齒輪傳動的同時嚙合輪齒對數(shù)較多，因此承載能力大[6]。高重合度直齒圓柱齒輪的設(shè)計，一般是通過增大齒頂高系數(shù)、增大齒數(shù)、減小壓力角等方法來實現(xiàn)[7]。增大齒頂高系數(shù)可以顯著提高齒輪副的重合度，但會導(dǎo)致齒頂厚減小使齒頂變尖，容易引起靠近齒頂?shù)牟课话l(fā)生崩齒。齒輪副齒數(shù)增大重合度也會增加，但會使齒輪傳動體積增大。減小壓力角可使重合度增大，且能增大齒厚，減小輪齒變尖的可能，但是在滿足重合度要求的基礎(chǔ)上，要盡量選擇較大的壓力角，以減小彎曲和接觸應(yīng)力，提高齒輪傳動的強度。上述高重合度齒輪的設(shè)計方法都屬于被動設(shè)計方法，無法實現(xiàn)通過改變預(yù)先給定的齒輪傳動重合度主動控制齒廓形狀的目的。為了滿足航空和汽車領(lǐng)域中齒輪傳動的高平穩(wěn)性和高重合度需求，本文作者提出了一種基于壓力角的高重合度齒輪主動設(shè)計方法，即在建立的壓力角為一次函數(shù)時齒廓方程描述方法的基礎(chǔ)上，首先給定齒輪傳動的重合度，根據(jù)齒輪嚙合原理，推導(dǎo)齒廓的形狀，從而達到通過改變重合度主動控制齒廓形狀的目的。

1 基于壓力角的齒廓方程的數(shù)學(xué)描述

齒輪嚙合坐標系如圖1 所示。齒輪1 與坐標系Σ 1(O1,x1,y1)相固聯(lián)，齒輪2 與坐標系Σ2(O2,x2,y2)相固聯(lián)，通過節(jié)點P 的曲線為嚙合線。假設(shè)兩齒廓接觸點從P 移動到P1位置，則P1位置的接觸點在坐標系Σ 0(O0,x0,y0)下的坐標可表示為

其中：λ 為極距；α 為壓力角；r1為齒輪1 的節(jié)圓半徑。

若齒輪1 順時針旋轉(zhuǎn)一個角度φ，坐標系Σ1(O1,x1,y1)也隨著齒輪1 旋轉(zhuǎn)一個角度到達圖示位置，則P1點在坐標系Σ1(O1,x1,y1)下的坐標可表示為

同樣地，將接觸點在坐標系Σ1(O1,x1,y1)中的坐標變換到坐標系Σ2(O2,x2,y2)下，可得齒輪2 的齒廓方程為

其中：u 為齒輪傳動比。

當φ 取正值時，式(2)和式(3)分別代表齒輪1 齒廓齒頂部分和齒輪2 齒廓齒根部分的坐標方程，而當φ取負值時，分別代表齒輪1 齒廓齒根部分和齒輪2 齒廓齒頂部分的坐標方程，同時，極距λ 和壓力角α 都分別是轉(zhuǎn)角φ 的函數(shù)。

根據(jù)齒輪嚙合原理，齒廓上接觸點的法向量和切向量應(yīng)相互垂直，即

其中：τ11為接觸點P1在齒輪1 齒廓上的切向量；n11為接觸點P1在齒輪1 齒廓上的法向量。

根據(jù)圖1 和式(2)，式(4)可化簡整理為

要使兩齒輪正確嚙合，齒輪2 齒廓上接觸點的法向量和切向量也應(yīng)相互垂直，同樣可計算得式(5)成立。

由式(2)，(3)和(5)可知，若壓力角α 也能表示成φ的函數(shù)，則式(2)和(3)分別表示相互嚙合的兩齒輪即齒輪1 和齒輪2 的齒廓方程。

2 以壓力角描述的齒輪嚙合特性的定量表達式

2.1 曲率

齒輪1 上接觸點的曲率為[8]

同理可得，齒輪2 上接觸點的曲率為

由式(6)和(7)可知，當曲率為正值時，齒廓為凸齒廓，因此，齒廓上尖點存在的條件為曲率半徑為0，即

2.2 滑動系數(shù)

根據(jù)齒輪嚙合原理，齒輪1 和2 的滑動系數(shù)可表示為[9]

其中，

將式(10)和(11)代入式(9)可得，齒輪1 和2 的滑動系數(shù)為

2.3 重合度

齒輪的重合度可表示為[10-11]

其中：z1為齒輪1 的齒數(shù)。

假設(shè)齒輪的α 可表示為φ 的一次函數(shù)，即

其中：C0和C1分別為一次函數(shù)的系數(shù)； αa為節(jié)圓和齒頂圓之間齒廓上點的壓力角；αd為節(jié)圓和齒根圓之間齒廓上點的壓力角。

當兩齒輪的一對輪齒開始嚙合時，必為主動輪的齒根推動從動輪的齒頂。因此，開始時嚙合點是從動輪的齒頂與嚙合線的交點(如圖1(a)所示)，隨著嚙合傳動的進行，結(jié)束時的嚙合點是主動輪的齒頂與嚙合線的交點(如圖1(b)所示)，為了方便計算，以節(jié)點處為參考點。

由圖1 可知，下式成立：

其中：r1和r2分別表示齒輪1 和齒輪2 的分度圓半徑；λa表示節(jié)圓和齒頂圓之間齒廓上的點對應(yīng)的極距；λd表示節(jié)圓和齒根圓之間齒廓上的點對應(yīng)的極距。

當齒輪1 的模數(shù)m，齒數(shù)z1，傳動比u，齒頂高ha，壓力角的函數(shù)表達式已知時，根據(jù)式(14)和式(16)可求得接觸點位于齒廓頂部時對應(yīng)的φa(φa＞0)的大小，根據(jù)式(15)和式(17)可求得接觸點位于齒廓根部時對應(yīng)的φd(φd＜0)。將其代入式(13)，可計算得齒輪傳動的重合度。

3 高重合度齒廓主動設(shè)計數(shù)學(xué)模型的建立

一元二次方程(16)和(17)根的存在條件為

由圖1 可知：壓力角范圍為[0,π/2]，而當φ=0 時，α=C0±C1φ 表 示 分 度 圓 上 的 壓 力 角，因 此C0滿 足0≤C0≤π/2。

當C0=0，將C0=0 代入式(14)，(15)，(5)和(8)可知，當φ=0 時，不管C1取何值，齒廓都會存在尖點。因此，C0=0 這種情況不適合設(shè)計齒輪齒廓。

當C1在上述范圍內(nèi)，且C1＜0 時， φa和 φd的取值范圍分別為

而當C1＞0 時， φa和 φd的取值范圍分別為

將式(22)和式(23)代入式(13)，即可計算得齒輪的重合度。同理，當齒輪的模數(shù)m，齒數(shù)z1，傳動比u，齒頂高系數(shù)ha*，頂隙系數(shù)C*，系數(shù)C0已知時，根據(jù)式(22)、式(23)和式(13)，可畫出系數(shù)C1和重合度ε的關(guān)系曲線。根據(jù)該曲線關(guān)系，當重合度確定后，可選擇系數(shù)C0和C1，將C0和C1分別代入式(14)，(15)和(5)，可計算得α 和λ，將上述數(shù)值代入式(8)驗證是否出現(xiàn)尖點，若不出現(xiàn)尖點，則可將上述數(shù)值分別代入式(2)和(3)可得主動齒廓和從動齒廓的表達式；若出現(xiàn)尖點，則需重新選擇系數(shù)C0和C1。

4 設(shè)計實例

4.1 齒廓方程

為了驗證上述設(shè)計方法的正確性，以重合度為3.0的齒輪傳動為例，說明其具體設(shè)計過程。齒輪的其他參數(shù)為模數(shù)m=2.5，齒頂高系數(shù)ha*=1，頂隙系數(shù)C*=0.25，齒數(shù)z1=25，z2=39。根據(jù)式(22)、式(23)和式(13)，若選取C0為π/8 和π/16，則系數(shù)C1和重合度ε的關(guān)系曲線如圖2 所示。當齒輪的重合度為3.0 時，本文選取C0=π/8，從而可得C1=0.539 9。將C0=π/8 和C1=0.539 9 代入式(14)，(15)和(5)，可計算得α 和λ，將上述數(shù)值代入式(8)可知，齒廓不會出現(xiàn)尖點。因此，將上述數(shù)值分別代入式(2)和式(3)可得主動齒廓和從動齒廓的方程分別為

圖2 系數(shù)C1 和重合度ε的關(guān)系曲線Fig.2 Relationship between contact ratio ε and coefficient C1

圖3 所示為高重合度齒輪與漸開線齒輪齒廓對比圖，其中，虛線為漸開線齒廓，實線為高重合度齒廓。由圖3、式(24)和式(25)可知：

1) 當φ≥0 時，式(24)和式(25)的第1 個式子分別表示主動齒輪分度圓和齒頂圓之間的齒廓、從動齒廓齒根圓和分度圓之間的齒廓；而當φ≤0 時，式(24)和式(25)的第2 個式子分別表示主動齒輪分度圓和齒根圓之間的齒廓以及從動齒廓齒頂圓和分度圓之間的齒廓。

2) 與相同參數(shù)的漸開線齒廓相比，高重合度齒輪在分度圓以下部分的齒寬大于漸開線齒輪在分度圓以下部分的齒寬，而在分度圓以上部分則相反。

圖3 齒廓對比圖Fig.3 Tooth profiles of gears

4.2 曲率

根據(jù)式(6)和(7)，利用MATLAB 軟件，可得出齒廓的曲率隨轉(zhuǎn)角φ的變化曲線，如圖4 所示。其中實線表示主動齒廓的曲率，虛線表示從動齒廓的曲率。綜合曲率半徑是衡量齒輪接觸強度的重要指標，與齒廓接觸應(yīng)力和接觸強度有著密切的關(guān)系，因此有必要找出高重合度齒輪綜合曲率半徑的變化曲線。

圖4 齒廓曲率隨齒輪轉(zhuǎn)角的變化曲線Fig.4 Curvature of tooth profiles with different rotation angles

綜合曲率半徑ρ可表示為[12]：

當接觸部分曲率方向相同時取“-”號，相反則取“+”號。

圖5 所示為高重合度齒輪綜合曲率半徑隨轉(zhuǎn)角φ的變化曲線。由圖5 可知：高重合度齒輪的綜合曲率半徑很大，根據(jù)赫茲公式，兩輪齒在接觸處的最大接觸應(yīng)力與綜合曲率半徑成反比[13]，故高重合度齒輪在接觸點處的最大接觸應(yīng)力較小，對應(yīng)的接觸強度會相應(yīng)較大。

4.3 滑動系數(shù)

圖5 高重合度齒輪的綜合曲率半徑Fig.5 Composite curvature radius of proposed gear

圖6 高重合度齒輪的滑動系數(shù)Fig.6 Sliding coefficients of proposed gear drive

根據(jù)式(12)繪制高重合度齒輪的滑動系數(shù)如圖6所示。當φ＞0 時，可繪制主動齒輪分度圓與齒頂圓之間齒廓的滑動系數(shù)以及從動齒輪分度圓和齒根圓之間的滑動系數(shù)，而當φ＜0 時，可繪制主動齒輪分度圓與齒根圓之間齒廓的滑動系數(shù)以及從動齒輪分度圓和齒頂圓之間的滑動系數(shù)。

當壓力角為常數(shù)時，齒廓為漸開線齒廓[14]，若漸開線齒輪的壓力角為20°，其他參數(shù)和高重合度齒輪相同，則C0=π/9，C1=0。根據(jù)式(12) 繪制漸開線齒輪的滑動系數(shù)，如圖7 所示。

圖7 漸開線齒輪的滑動系數(shù)Fig.7 Sliding coefficients of involute gear drive

由圖6 和圖7 可知：節(jié)圓處滑動系數(shù)為0，而在遠離節(jié)圓處，滑動系數(shù)逐漸增大，表示齒面磨損量增大；高重合度齒輪主動輪的滑動系數(shù)變化區(qū)間為[-2.33, 0.61]，從動輪的滑動系數(shù)變化區(qū)間為[-1.55,0.70]，而漸開線齒輪主動輪的滑動系數(shù)變化區(qū)間為[-5.98, 0.68]，從動輪的滑動系數(shù)變化區(qū)間為[-2.12,0.86]，因此，高重合度齒輪的滑動系數(shù)遠遠小于漸開線齒輪的滑動系數(shù)，有利于改善齒輪的磨損狀況，提高齒輪的壽命。

5 動力學(xué)分析

為了分析重合度對齒輪動力學(xué)特性的影響，根據(jù)Hertz 碰撞理論，本文利用Adams 軟件，對設(shè)計實例中的高重合度齒輪進行了動力學(xué)分析。分別在兩齒輪上施加旋轉(zhuǎn)副，齒輪之間加實體-實體碰撞力，輸入軸上加轉(zhuǎn)速驅(qū)動，輸出軸上加負載轉(zhuǎn)矩，建立的完整動力學(xué)仿真模型如圖8 所示。

圖8 動力學(xué)仿真模型Fig.8 Dynamic simulation model

小齒輪施加的轉(zhuǎn)速為300 r/min，為了避免由于開始的沖擊造成的速度突變，轉(zhuǎn)速驅(qū)動以階躍函數(shù)施加，為STEP(time，0，0d，0.2，1 800d)。在被動齒輪的轉(zhuǎn)動軸上施加的負載函數(shù)為STEP(time，0，0，0.01，150 000)，大小為150 N·m 的階躍函數(shù)。接觸碰撞函數(shù)可表示為IMPACT(x，x，x1，k，e，Cmax，d)[15]，其中，x 表示用來計算碰撞函數(shù)里的距離變量；x相對于位移的碰撞速度；x1表示位移x 自由長度的正實數(shù)變量，若x＜x1則力賦一個正值，否則力為0 N，x1可以被定義為一個實數(shù)、函數(shù)或者變量；e 表示力的變形特性指數(shù)的一個正實變量；Cmax表示最大阻尼系數(shù)的一個非負實變量；d 表示在施加全阻尼時的邊界穿透量。根據(jù)設(shè)計實例中齒輪參數(shù)可計算得阻尼系數(shù)C=50 N·s/mm，穿透深度d=0.1mm，k=6.4×105N/mm，e=1.5。本文取仿真時間為0.5 s，步長為6 000。圖9所示為高重合度齒輪嚙合力仿真曲線。為了便于比較，圖10 給出了相同參數(shù)的漸開線齒輪的嚙合力仿真曲線。由圖9 和圖10 可知：

圖9 高重合度齒輪嚙合力曲線Fig.9 Meshing force of high contact ratio gear

圖10 漸開線齒輪嚙合力曲線Fig.10 Meshing force of involite gear

1) 與漸開線齒輪傳動相同，在仿真開始的瞬間，齒輪突然加速，嚙合力隨齒輪嚙合的進行發(fā)生波動，但波峰間距較大，即齒輪嚙合周期較長，隨著傳動的進行，嚙合力在一個均值附近以一定的幅值和周期上下波動，周期和幅值趨于穩(wěn)定，這是齒輪周期性嚙入嚙出沖擊的表現(xiàn)。

2) 漸開線齒輪傳動的徑向嚙合力在均值處的波動幅值為758.43 N，而高重合度齒輪傳動的徑向嚙合力在均值處的波動幅值為459.16 N，比漸開線齒輪傳動減少了39.5%。

3) 漸開線齒輪傳動的切向嚙合力在均值處的波動幅值為399.24 N，而高重合度齒輪傳動的切向嚙合力在均值處的波動幅值為351.05 N，比漸開線齒輪傳動減小了12.1%。因此，從仿真結(jié)果來看，高重合度齒輪傳動的嚙入嚙出沖擊相對較小，傳動更加平穩(wěn)。

6 結(jié)論

1) 利用基于壓力角的齒輪主動設(shè)計方法，可以在設(shè)計之前就預(yù)知齒輪傳動的重合度，還可以直接計算其曲率、滑動系數(shù)等嚙合特性，為方便快捷地設(shè)計滿足不同重合度需求的特殊齒輪奠定了基礎(chǔ)。

2) 與相同參數(shù)的漸開線齒輪相比，利用本文提出的方法設(shè)計的重合度為3 的齒輪，具有最大接觸應(yīng)力小，滑動系數(shù)小的優(yōu)點，有利于改善齒輪的磨損狀況，提高齒輪壽命。

3) 在相同條件下，高重合度齒輪傳動的徑向嚙合力在均值處的波動幅值比漸開線齒輪減少了39.5%，切向嚙合力在均值處的波動幅值比漸開線齒輪減小了12.1%，因此，其嚙入嚙出沖擊和傳動的平穩(wěn)性均優(yōu)于漸開線齒輪。

[1] Tsai M H,Tsai Y C.Design of high-contact-ratio spur gear using quadratic parametric tooth profiles[J]. Mechanism and Machine Theory,1998,33(5):551-564.

[2] TAN Weiming. Common definition for end-surface contact ratio of gears and its applications[J]. Chinese Journal of Mechanical Engineering,2004,17(4):595-597.

[3] LI Shuting. Effect of addendum on contact strength, bending strength and basic performance parameters of a pair of spur gears[J]. Mechanism and Machine Theory, 2008, 43(12):1557-1584.

[4] 尹剛. 高重合度齒輪應(yīng)力場有限元分析[J]. 重慶大學(xué)學(xué)報(自然科學(xué)版),2010,33(7):53-57.YIN Gang. Finite element analysis of high contact-ratio gear stress field[J].Journal of Chongqing University (Natural Science Edition),2010,33(7):53-57.

[5] 鄧效忠, 方宗德, 魏冰陽, 等. 高重合度弧齒錐齒輪加工參數(shù)設(shè)計與重合度測定[J]. 機械工程學(xué)報,2004,40(6):95-99.DENG Xiaozhong, FANG Zongde,WEI Bingyang,et al.Design and experiment of spiral bevel gears with high contact ratio[J].Chinese Journal of Mechanical Engineering,2004,40(6):95-99.

[6] 王三民, 紀名剛. 高速大重合度直齒圓柱齒輪的齒廓最佳修形研究[J]. 航空學(xué)報,1996,17(1):121-124.WANG Sanmin, JI Minggang. Study on the best profile modification of high-speed spur gears with high-contact-ratio[J].Acta Aeronautica et Astronautica Sinica,1996,17(1):121-124.

[7] 渠珍珍, 楊杭旭. 高重合度齒輪的實現(xiàn)方法及其參數(shù)設(shè)計[J].機電技術(shù),2012(3):103-105.QU Zhenzhen,YANG Hangxu.Realization and parameter design of high contact ratio gears[J]. Mechanical & Electrical Technology,2012(3):103-105.

[8] 林菁. 基于嚙合角函數(shù)的平面圓齒輪共軛齒廓求解及幾何特性研究[J]. 機械傳動,2004,28(3):6-8.LIN Jing. A direct-profile-design method for generating noncircular gear pairs[J]. Journal of Mechanical Transmission,2004,28(3):6-8.

[9] Yeh T, Yang D C H, Tong S H, Design of new tooth profile for high-load capacity gears[J]. Mechanism and Machine Theory,2001,36(10):1105-1120.

[10] Pedrero J I, Artés M, García-Prada J C. Determination of the addendum modification factors for gears with pre-established contact ratio[J]. Mechanism and Machine Theory, 1996, 31(7):937-945.

[11] WANG Jian, HOU Liang, LUO Shanming, et al. Active design of tooth profiles using parabolic curve as the line of action[J].Mechanism and Machine Theory,2013,61(9):47-63.

[12] 孟兆明, 許基琴, 汪惠琴. 漸開線齒輪點蝕位置的確定及其計算綜合曲率半徑的誤差修正[J]. 青島科技大學(xué)學(xué)報,2004,25(4):339-342.MENG Zhaoming, XU Jiqin,WANG Huiqin.The definition and calculation of the location of spot erosion for involute gear correction of the errors for the composite curve radius[J].Journal of Qingdao University of Science and Technology, 2004, 25(4):339-342.

[13] 劉令濤, 焦永和, 張春林, 等. 內(nèi)平動全擺線齒輪曲率半徑的分析[J]. 北京理工大學(xué)學(xué)報,2010,30(2):158-160,173.LIU Lingtao, JIAO Yonghe, ZHANG Chunlin, et al. Analysis of the curvature radius of internal parallel moving gears with whole cycloidal profiles[J]. Journal of Beijing Institute of Technology,2010,30(2):158-160,173.

[14] 林菁. 共軛曲線嚙合角函數(shù)理論及其應(yīng)用[M]. 北京: 科學(xué)出版社,2005:31-35.LIN Jing. Theory and application of working pressure angle function for conjugate curves[M]. Beijing: Science Press, 2005:31-35.

[15] 李三群, 賈長治, 武彩崗, 等. 基于虛擬樣機技術(shù)的齒輪嚙合動力學(xué)仿真研究[J]. 系統(tǒng)仿真學(xué)報,2007,19(4):901-904.LI Sanqun, JIA Changzhi, WU Caigang, et al. Dynamic simulation study of gear meshing based on virtual prototyping technology[J]. Journal of System Simulation, 2007, 19(4):901-904.