齒輪齒面溫度場計算

賴 鋒，林騰蛟，邢惠斌，姚陽迪，劉 文，孫建國

(1.中國船舶重工集團公司第七○三研究所，黑龍江 哈爾濱 150036;2.重慶大學，重慶 400030)

齒輪齒面溫度場計算

賴 鋒1，林騰蛟2，邢惠斌1，姚陽迪2，劉 文2，孫建國2

(1.中國船舶重工集團公司第七○三研究所，黑龍江 哈爾濱 150036;2.重慶大學，重慶 400030)

齒輪齒面溫度場的大小對齒輪傳動的性能與失效以及齒輪潤滑系統的設計等有著重要的影響。齒輪溫度的熱平衡狀態和齒輪本體溫度以及齒面瞬時溫度的變化極為復雜。本文綜合考慮齒輪轉速，潤滑油熱傳導率、密度、比熱、運動粘度以及間隙冷卻過程中的標準化總冷卻量等影響，給出了齒輪齒面溫度場計算方法。本次研究成果對理論、設計具有一定的指導意義。

齒輪;齒面溫度場計算;機械傳動

0 引言

輪齒表面溫度的大小對齒輪傳動的性能與失效以及齒輪潤滑系統的設計與操控等有著重要的影響。齒輪傳動過程中輪齒嚙合面摩擦熱流量的產生，齒輪在潤滑冷卻系統作用下的強制對流傳熱等因素的綜合影響，使得齒輪溫度的熱平衡狀態和齒輪本體溫度以及齒面瞬時溫度的變化極為復雜。因此，研究齒輪齒面的溫度場分析方法，建立精確的齒輪輪齒的溫度分析模型，開發齒輪齒面溫度場計算程序，并應用于工程實際，將有助于指導齒輪傳動裝置及潤滑冷卻系統的設計并降低設計所需的周期和費用，這對高速重載齒輪傳動設計尤其重要。本文綜合考慮齒輪轉速，潤滑油熱傳導率、密度、比熱、運動粘度以及間隙冷卻過程中的標準化總冷卻量等影響，給出了齒輪齒面溫度場計算方法。

1 漸開線齒輪參數化建模

對于漸開線圓柱齒輪，可以將嚙合運動看作是在一個平面內或平行一個平面。下面運用運動學法重點推導斜齒輪齒廓方程。

1.1 齒廓方程公式的推導

經過坐標轉換等，推導了斜齒輪齒廓方程，如下:

式中:b為齒厚;γ為滾動角;β為分度圓螺旋角;r為分度圓半徑。

1.2 漸開線齒輪建模

本文借助ANSYS的APDL語言自行編制了一套參數化建模程序。輸入齒輪的主要幾何參數便可以建立齒輪副接觸模型，并實現有限元網格自動劃分。本例建立的內嚙合齒輪副接觸有限元網格模型如圖1所示。

圖1 內嚙合斜齒輪有限元模型Fig.1 Finite element model of inner mesh inclined gear

2 齒輪運轉過程接觸分析

本例采用了赫茲接觸理論和有限元方法分析了嚙合狀態下標準齒廓及修形齒廓的漸開線齒輪齒面接觸壓力的分布和變化規律。根據漸開線圓柱齒輪的傳動特點，主、從動輪的嚙合位置一直處在1條直線上，經計算，斜齒輪嚙合時總嚙合線變化如圖2所示。

圖2 斜齒輪嚙合時總嚙合線變化Fig.2 Variational mesh-line when the inclined gear meshing

3 齒輪齒面嚙合摩擦生熱分析

的摩擦熱流量對齒輪溫度和分布規律的改變以及輪齒熱平衡產生了重要的影響。本文著重研究了摩擦系數和齒面摩擦熱的計算方法，并分析了其沿嚙合面的分布規律及相關影響因素。

圖3和圖4為未修形和修形的外嚙合斜齒輪齒面摩擦熱流量。由圖可知，外嚙合斜齒輪在節點附近熱流量較小，與主、從動輪和修形與否關系不大。在齒根與齒頂處，未修形齒輪的摩擦熱流量相對整個齒面較大而修形齒輪的摩擦熱流量較小。主、從動輪嚙合面的摩擦熱流量分配是不均的，規律基本上和直齒輪相似，但斜齒輪在齒寬方向上齒面摩擦熱流量不是均勻分布的，不像直齒輪在齒寬方向上是成對稱分布的。

齒輪參數和齒廓形狀直接決定了齒間載荷分配和接觸壓力的分布，以及接觸輪齒對嚙合齒面的相對滑動。在高速運轉條件下嚙合齒對間的相對滑動導致了顯著的摩擦熱生成，它是由齒面承受的載荷、旋轉速度的大小和齒面的摩擦系數所共同決定的。齒面上由摩擦產生的熱能一部分通過熱傳導進入齒輪本體，而另一部分則由作用于嚙合齒面的潤滑油通過強制對流傳熱而冷卻擴散。因此，齒輪嚙合齒面產生

4 齒輪嚙合的對流傳熱分析

在齒輪溫度場分析計算中必須事先確定對流傳熱系數，齒輪換熱系數主要取決于齒輪的冷卻潤滑方式及齒輪的運行條件等，沒有現有的理論可以精確計算，一般都采用經驗公式法去確定對流傳熱系數。下面給出的對流傳熱系數均是國內外文獻中提供的經驗公式。

本例在計算輪齒對流傳熱系數時，潤滑油采用L－ECD40柴油機油，40℃時其性能參數如表1所示。

表1 潤滑油性能參數Tab.1 Lubrication performance parameters

根據斜齒輪幾何參數及工況參數計算端面對流傳熱系數。端面對流傳熱系數隨轉速的變化如圖5所示。由圖可知，斜齒輪端面的對流傳熱系數隨轉速增大而增大。

圖5 斜齒輪端面對流傳熱系數隨速度變化曲線Fig.5 The curve of convection diathermanous coefficient by variational speed of the inclined gear

5 輪齒溫度場有限元分析

由于嚙合輪齒接觸滑動而產生的齒面摩擦加熱，齒輪運轉過程中輪齒的表面溫度是不斷變化的。在齒面摩擦滑動生熱和潤滑油強制對流冷卻作用的綜合影響下，齒輪輪齒本體的溫度經過一定周期的旋轉后可在固定載荷條件下達到一個熱的平衡。盡管達到熱平衡狀態后齒面的瞬時溫度(閃溫)仍然隨時間而改變，但輪齒本體的溫度在一個旋轉周期內的變化則非常小。這是因為達到平衡狀態后，輪齒在每一嚙合周期內齒面摩擦熱量的極大部分都被齒面冷卻的潤滑油帶走:另一個原因是高速運轉條件下齒面獲得的熱量由齒面向齒體傳導的時間非常短。因此，達到熱平衡狀態后輪齒的本體溫度基本上是相對于時間穩定的隨位置變化的分布，瞬時溫度則是在本體溫度基礎上隨時間和位置而改變的分布。

5.1 齒輪溫度場有限元分析模型的建立

首先建立齒輪的單齒有限元模型，在進行齒輪溫度場分析時，需要加載熱流密度及對流傳熱邊界條件。齒輪對流傳熱的加載邊界如圖6所示。

圖6 齒輪對流傳熱的加載邊界Fig.6 Convection diathermanous verge load of the gear

熱流密度的加載通常采取在齒面上直接加載平均熱流密度的方法來實現，這樣加載的熱流密度在齒面上是平均分布的，但實際上，輪齒齒面的熱流密度分布是呈一定規律變化的，為了使熱流密度加載更符合實際情況，對每一個節點加載熱流密度。本文以未修形內嚙合斜齒輪為例給出加載熱流密度后的輪齒溫度場有限元模型，并計算修形前后輪齒的溫度場。

5.2 邊界條件的處理

摩擦熱流量加載到齒面后的齒輪溫度場有限元模型如圖7所示。5.3 齒輪齒面溫度場計算

圖7 斜齒輪溫度場有限元網格及熱流量邊界Fig.7 The finite element reseau and heat flux verge of the inclined gear

未修形及修形內嚙合斜齒輪溫度場計算結果分別如圖7及圖8所示。由圖可知，斜齒輪的溫度場規律和直齒輪齒面溫度場分布規律相似，齒輪節圓附近并不是最高溫度區，齒輪的最高溫度出現在齒根附近或齒頂附近。修形后齒輪的溫度要明顯降低。與直齒輪不同的是在齒寬方向溫度分布是不對稱的，最高溫度區偏向嚙出端，上述情況分別見圖8和圖9所示。

6 結語

通過建立標準漸開線齒輪及修形齒輪實體模型和彈性接觸有限元模型并分析齒面接觸力的變化規律和相關的影響因素，給出了計算齒輪嚙合過程中主、從動齒輪齒面摩擦系數的方法，以及齒輪齒面摩擦熱流量計算方法，并分析齒輪齒面摩擦熱流量沿嚙合面的分布規律、大小以及相關的影響因素;齒輪端面的對流傳熱系數以及輪齒嚙合面的對流傳熱系數的計算方法，分析了對流傳熱系數的大小和分布及其影響因素。在此基礎上推導了輪齒本體溫度的熱平衡方程，建立了本體溫度有限元分析模型以及熱載荷和對流傳熱的邊界條件，對不同轉速和載荷工況下輪齒本體溫度進行了計算，并分析了輪齒本體溫度變化規律，得到了輪齒的溫度云圖、位移云圖以及應力云圖，建立了齒輪齒面溫度場計算方法，該計算方法對理論研究及工程應用有一定的指導意義。

［1］李潤方，龔劍霞.接觸問題數值方法及其在機械設計中的應用［M］.重慶:重慶大學出版社，1991.

LI Run-fang，GONG Jian-xia.The numerical method of the contact problem and the application of it in mechanical design［M］.Chongqing:Chongqing University Press，1991.

［2］李瑰賢.空間幾何建模及工程應用［M］.北京:高等教育出版社，2007.

LI Gui-xian.The model fouding of spacical geometry and its application in project.［M］.Beijing:Higher Education Press，2007.

［3］龔桂義，陳式椿，王永潔.漸開線圓柱齒輪強度計算與結構設計［M］.北京:機械工業出版社，1986.

GONG Gui-yi，CHEN Shi-chun，WANG Yong-jie.The calculation of intensity and the design of sructure for the involute circular gear［M］.Beijing:Mechanical Industry Press，1986.

［4］薩本佶.高速齒輪傳動設計［M］.北京:機械工業出版社，1986.

SA Ben-ji.The design of the hign-speed-transmission gear［M］.Beijing:Mechanical Industry Press，1986.

［5］李潤方.齒輪傳動的剛度分析和修形方法［M］.重慶:重慶大學出版社，1998.

LI Run-fang.Gear-driven rigidity analysis and modify method［M］.Chongqing:Chongqing University Press，1998.

［6］龍慧.高速齒輪傳動輪齒的溫度模擬及過程參數的敏感性分析［D］.重慶:重慶大學，2001.

LONG Hui.The temperature emulation of hign speed driven gear and the sensitivity analysis of parameters under processes.［D］.Chongqing:Chongqing University，2001.

［7］王統.用有限元法分析估算直齒輪的體積溫度［J］.上海交通大學學報，1981，(4):53 －71.

WANG Tong.Theanalysisofstraightgear'svolume temperature by finite element［J］.Jornal of Shanghai Jiaotong University，1981，(4):53 －71.

［8］李潤方，湯慶平.運轉過程中輪齒偶合熱彈性接觸有限元分析［J］.齒輪，1989，13(1):23 －28，57.

LI Run-fang，TANG Qing-ping.The analysis of gear tooth's coupling heat elasticity under operating by contact finite element method.［J］.Gear，1989，13(1):23 －28，57.

［9］方宗德，張永才.斜齒輪的齒廓修形［J］.汽車工程，1993，15(1):59 －64.

FANG Zong-de，ZHANG Yong-cai.The dentiform finishing of inclined gear［J］.Automobile Project，1993，15(1):59－64.

The computational research of the gear's tooth surface temperature

LAI Feng1，LIN Teng-jiao2，XING Hui-bin1，YAO Yang-di2，LIU Wen2，SUN Jian-guo2
(1.The 703 Research Institute of CSIC，Harbin 150078，China;2.Chongqing University，Chongqing 400030，China)

It is an important influence to performance and fail of gear transmission to temperature field of tooth surface.The gear lubrication system will be influenced by the temperature field.The heat balance status，tooth temperature and tooth surface instantaneous temperature is very complicated.In the paper，we consider synthetically gear rotate speed，heat exchange ratio of lubricating oil，density，specific heat，kinetic viscosity and so on;we give the computational method to the temperature field of tooth surface under the conditions.We can find that the results are significative to theory and design in the gear transmission.

gear;temperature field;computational method;mechanical drive

TH132.3

A

1672－7649(2012)03－0050－04

10.3404/j.issn.1672－7649.2012.03.010

2011－03－28;

賴鋒(1977－)，男，工程師，長期從事船舶后傳動設計工作。