王歡歡

(中車青島四方機(jī)車車輛股份有限公司，山東 青島 266000)

0 引言

行星齒輪在傳動(dòng)過程中，由于承受較大的載荷，輪齒溫度會(huì)明顯升高。如果通過實(shí)驗(yàn)方法對(duì)齒輪溫度場(chǎng)進(jìn)行測(cè)定，耗時(shí)長(zhǎng)且設(shè)備昂貴，不適宜廣泛應(yīng)用。因此，建立一種較為精確的行星齒輪溫度場(chǎng)分析模型，對(duì)其溫度場(chǎng)的主要影響因素進(jìn)行分析，將有助于提升行星齒輪傳動(dòng)系統(tǒng)、冷卻散熱系統(tǒng)的優(yōu)化設(shè)計(jì)，具有十分重要的工程應(yīng)用價(jià)值。

嚙合齒輪的溫度場(chǎng)主要由穩(wěn)態(tài)溫度場(chǎng)及瞬態(tài)溫度場(chǎng)兩部分組成。對(duì)于齒輪穩(wěn)態(tài)溫度場(chǎng)，國(guó)內(nèi)外學(xué)者已進(jìn)行了大量研究。Thyla等建立了齒輪的有限元模型，計(jì)算了齒面熱流量分布，分析了輪齒的穩(wěn)態(tài)溫度場(chǎng)。肖望強(qiáng)等探討了非對(duì)稱齒輪的穩(wěn)態(tài)溫度場(chǎng)，研究表明該類型齒輪穩(wěn)態(tài)溫度低于對(duì)稱齒輪。羅彪等建立了齒輪單齒的熱力耦合有限元分析模型，通過引入修正系數(shù)，分析了輪齒的穩(wěn)態(tài)溫度場(chǎng)及熱變形狀況。方特等以彈流潤(rùn)滑狀況下的高速直齒輪為研究對(duì)象，利用熱網(wǎng)絡(luò)法計(jì)算了輪齒的穩(wěn)態(tài)溫度場(chǎng)，并通過齒輪試驗(yàn)臺(tái)進(jìn)行了驗(yàn)證。樊智敏等利用ANSYS有限元軟件，分析了雙漸開線齒輪的穩(wěn)態(tài)溫度場(chǎng)，以及齒輪階梯參數(shù)對(duì)溫度場(chǎng)的影響。另外，姚陽(yáng)迪等學(xué)者也對(duì)齒輪的穩(wěn)態(tài)溫度場(chǎng)進(jìn)行了研究。以往文獻(xiàn)采用不同方法對(duì)行星齒輪穩(wěn)態(tài)溫度場(chǎng)進(jìn)行了研究，但關(guān)于轉(zhuǎn)矩、轉(zhuǎn)速等因素對(duì)溫度場(chǎng)的作用及影響研究較少。

因此本文以2K-H(NGW)型三級(jí)行星齒輪減速器為研究對(duì)象，首先分析輪齒間的相對(duì)滑動(dòng)速度、齒面應(yīng)力、摩擦系數(shù)，得到輪齒表面的摩擦熱流量分布；然后基于有限單元法得到行星齒輪的溫度場(chǎng)分布，分析轉(zhuǎn)矩等因素對(duì)行星齒輪穩(wěn)態(tài)溫度場(chǎng)的影響。

1 行星齒輪基本參數(shù)

本文所研究的2K-H(NGW)型三級(jí)行星齒輪減速器傳動(dòng)簡(jiǎn)圖如圖1所示。

輪齒的溫度在傳動(dòng)過程中是變化的，但閃溫(輪齒接觸時(shí)的瞬時(shí)溫度)出現(xiàn)時(shí)間極短、影響范圍很小，故一般認(rèn)為穩(wěn)定傳動(dòng)時(shí)輪齒溫度為定值。不同輪齒在傳動(dòng)過程中所獲得熱流密度是相同的，所以研究齒輪溫度場(chǎng)時(shí)可以取單個(gè)輪齒來進(jìn)行研究。表1為第一級(jí)太陽(yáng)輪和行星輪基本參數(shù)。

圖1 2K-H型三級(jí)行星減速器機(jī)構(gòu)傳動(dòng)簡(jiǎn)圖

輸入扭矩T(N.m)輸入轉(zhuǎn)速ns(rpm)齒數(shù)Z1Z2模數(shù)m(mm)壓力角α(°)14891145.625 20420變位系數(shù)齒高系數(shù)ha/hf齒寬b(mm)密度ρ1 ρ2(kg/m3)泊松比μ1 μ20 0.261/1.46078500.3彈性模量E1E2(Gpa)中心距a(mm)熱導(dǎo)率λ1λ2W/(m.K)比熱容C1C2J/(kg.K)粗糙度Ra(um)20691504650.8

2 摩擦熱流量計(jì)算

在傳動(dòng)過程中，齒輪溫度場(chǎng)熱源主要來自于嚙合過程中的摩擦熱流。該摩擦熱流量由輪齒間的相對(duì)滑動(dòng)速度、齒面接觸應(yīng)力及摩擦系數(shù)等共同決定。

齒輪嚙合時(shí)的瞬時(shí)熱流密度q1C可以分別表示為：

q1C=βγμCpnVt

(1)

式中γ——熱能轉(zhuǎn)換系數(shù)；

μC——齒面摩擦系數(shù)；

β—— 熱分配系數(shù)；

pn—— 齒面接觸應(yīng)力；

Vt—— 相對(duì)滑動(dòng)速度。

根據(jù)上式，可求得齒面在傳動(dòng)過程中的平均熱流密度為：

(2)

由式(2)推導(dǎo)出的齒輪齒面摩擦熱流量分布如圖2所示。

圖2 齒面摩擦熱流量分布

3 齒輪溫度場(chǎng)影響因素分析

3.1 齒輪穩(wěn)態(tài)溫度場(chǎng)

在穩(wěn)定工況下，潤(rùn)滑油油溫約為70℃，在分析時(shí)將此溫度設(shè)置為環(huán)境溫度。齒面、齒根及端面等部位的對(duì)流換熱系數(shù)，根據(jù)之前學(xué)者的研究進(jìn)行確定。通過建立輪齒有限元三維模型，在對(duì)應(yīng)邊界施加摩擦熱流密度及對(duì)流換熱系數(shù)，利用ANSYS熱分析模塊進(jìn)行了齒輪穩(wěn)態(tài)溫度場(chǎng)的求解。輪齒穩(wěn)態(tài)溫度場(chǎng)如圖3所示。

圖3 輪齒穩(wěn)態(tài)溫度場(chǎng)

3.2 輸入轉(zhuǎn)矩影響分析

當(dāng)輸入轉(zhuǎn)矩為原轉(zhuǎn)矩50%時(shí)，齒輪本體溫度場(chǎng)和未修形齒輪本體溫度場(chǎng)如圖4所示。齒輪本體溫度場(chǎng)最高溫度為82.8℃，和原工況條件下的本體溫度場(chǎng)將比較，最高溫度下降15.4℃。不同輸入轉(zhuǎn)矩下，輪齒最高溫度如圖5所示。當(dāng)轉(zhuǎn)矩為原轉(zhuǎn)矩的2倍時(shí)，輪齒最高溫度為131.9℃，相比原工況升高33.7℃。

圖5 不同轉(zhuǎn)矩下的輪齒最高溫度

3.3 輸入轉(zhuǎn)速影響分析

當(dāng)轉(zhuǎn)速為573rpm，即穩(wěn)定工況下轉(zhuǎn)速的50%時(shí)，齒輪穩(wěn)態(tài)溫度場(chǎng)如圖6所示。最高溫度為89.6℃，和原工況條件下的最高溫度比較下降8.6℃。不同輸入轉(zhuǎn)速下，輪齒最高溫度如圖7所示。當(dāng)轉(zhuǎn)速為2291rpm時(shí)，輪齒最高溫度為112.2℃，和原工況相比升高14.0℃。

圖6 輪齒穩(wěn)態(tài)溫度場(chǎng) (573rpm)

圖7 不同轉(zhuǎn)速下的輪齒最高溫度

3.4 環(huán)境油溫影響分析

當(dāng)穩(wěn)定狀態(tài)下潤(rùn)滑油溫度為60℃時(shí)，齒輪本體溫度場(chǎng)如圖8所示，最高溫度為88.2℃，和初始工況相比下降10.0℃。不同環(huán)境油溫下，輪齒最高溫度如圖9所示，可知潤(rùn)滑油溫度變化與輪齒溫度變化呈線性關(guān)系。

圖8 輪齒穩(wěn)態(tài)溫度場(chǎng) (60℃)

圖9 不同油溫下的輪齒最高溫度

4 結(jié)論

(1)對(duì)于穩(wěn)態(tài)溫度場(chǎng)，轉(zhuǎn)矩的影響明顯要高于轉(zhuǎn)速的影響。因?yàn)楫?dāng)轉(zhuǎn)速增加時(shí)，明顯改善了對(duì)流換熱狀況。

(2)環(huán)境油溫對(duì)穩(wěn)態(tài)溫度場(chǎng)的影響呈線性關(guān)系，當(dāng)潤(rùn)滑油溫度升高或降低一定溫度時(shí)，齒輪的溫度變化與該值大致相同。

(3)在滿足工作條件的前提下，通過降低轉(zhuǎn)速和轉(zhuǎn)矩，有助于降低齒輪的穩(wěn)態(tài)溫度場(chǎng)。

[1] Thyla P R,Rudramoorthy R.Investigations On Gear Tooth Surface And Bulk Temperatures Using ANSYS[J].Psg College of Technology,2004.

[2] Fernandes C M C G,Rocha D M P,Martins R C,et al.Finite element method model to predict bulk and flash temperatures on polymer gears[J]. Tribology International, 2018, 120.

[3] 肖望強(qiáng),李威,韓建友.非對(duì)稱齒廓漸開線齒輪傳動(dòng)的熱分析[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2006,37 (12):164-167.

[4] 羅彪,葉江,李威,李必文.齒輪穩(wěn)態(tài)溫度場(chǎng)及熱變形研究[J].南華大學(xué)學(xué)報(bào)(自然科學(xué)版),2016,30(04):38-43.

[5] 方特,劉少軍.彈流潤(rùn)滑下高速直齒輪穩(wěn)態(tài)溫度場(chǎng)研究[J].機(jī)械傳動(dòng),2017,41(05):132-137.

[6] 樊智敏,王瑞雪,周萬峰.雙漸開線齒輪傳動(dòng)穩(wěn)態(tài)溫度場(chǎng)模擬分析[J].機(jī)械強(qiáng)度,2017,39(01):198-203.

[7] 姚陽(yáng)迪,林騰蛟,何澤銀.高速斜齒輪傳動(dòng)穩(wěn)態(tài)溫度場(chǎng)仿真分析[J].機(jī)械研究與應(yīng)用,2009,22(06):9-12+16.

[8] 沈允文,王彤,王三民,張永紅,劉志全.弧齒錐傳動(dòng)的穩(wěn)態(tài)本體溫度場(chǎng)分析[J].機(jī)械傳動(dòng),2001，(03):1-4+0.

[9] 徐家強(qiáng),楊磊,蔡光達(dá).少齒數(shù)齒輪溫度場(chǎng)研究[J].機(jī)械研究與應(yīng)用,2015,28(04):51-53+56.

[10]郁晗.齒輪接觸應(yīng)力和溫度場(chǎng)分析[J].機(jī)械工程與自動(dòng)化,2013(05):50-52.

[11]龍慧. 高速齒輪傳動(dòng)輪齒的溫度模擬及過程參數(shù)的敏感性分析[D].重慶：重慶大學(xué),2001.