行星輪系漸開(kāi)線變位齒輪時(shí)變嚙合剛度數(shù)值計(jì)算

岳喜錚，丁問(wèn)司，丁　康， 曾智杰

(華南理工大學(xué) 機(jī)械與汽車工程學(xué)院, 廣州　510641)

?

行星輪系漸開(kāi)線變位齒輪時(shí)變嚙合剛度數(shù)值計(jì)算

岳喜錚，丁問(wèn)司，丁康， 曾智杰

(華南理工大學(xué) 機(jī)械與汽車工程學(xué)院, 廣州510641)

摘要：齒輪嚙合剛度的時(shí)變性是齒輪產(chǎn)生振動(dòng)的主要原因，對(duì)齒輪時(shí)變剛度進(jìn)行求解計(jì)算具有重要的意義?；谡兾恍行禽喯低恺X輪齒形幾何參數(shù)的計(jì)算方法，根據(jù)內(nèi)外齒輪齒形的不同，推導(dǎo)出內(nèi)齒輪齒形幾何參數(shù)計(jì)算過(guò)程。利用石川公式中輪齒變形計(jì)算公式和剛度求解公式，通過(guò)Matlab編程求得正變位行星輪系內(nèi)外嚙合時(shí)變嚙合剛度曲線，并與Weber能量法計(jì)算結(jié)果進(jìn)行對(duì)比，結(jié)果表明：靜態(tài)剛度最大誤差為3.7%，不考慮幅值為很小值時(shí)的動(dòng)態(tài)剛度，前4階諧波幅值最大誤差約為9.1%，精度滿足要求，可以為行星輪系運(yùn)動(dòng)微分方程提供較為準(zhǔn)確的嚙合剛度。

關(guān)鍵詞：變位齒輪；石川公式；外嚙合剛度；內(nèi)嚙合剛度

齒輪傳動(dòng)因結(jié)構(gòu)緊湊、效率高、能傳遞大的作用力等優(yōu)點(diǎn)而成為機(jī)械傳動(dòng)中應(yīng)用最廣泛的一種傳動(dòng)形式。齒輪傳動(dòng)系統(tǒng)是一種彈性的機(jī)械系統(tǒng)，在動(dòng)態(tài)激勵(lì)作用下會(huì)產(chǎn)生動(dòng)態(tài)響應(yīng)。動(dòng)態(tài)激勵(lì)分為外部激勵(lì)和內(nèi)部激勵(lì)兩類，即使外部激勵(lì)為0，齒輪系統(tǒng)也會(huì)因同時(shí)嚙合齒對(duì)數(shù)和嚙合接觸點(diǎn)位置的改變，以及輪齒的誤差等引起內(nèi)部激勵(lì)而產(chǎn)生振動(dòng)[1]。通常嚙合齒輪對(duì)的重合度不是整數(shù)，嚙合過(guò)程中單雙齒嚙合交替出現(xiàn)，隨時(shí)間做周期性變化，輪齒的綜合嚙合剛度也隨之變化。嚙合剛度隨時(shí)間變化產(chǎn)生剛度激勵(lì)，是產(chǎn)生齒輪嚙合內(nèi)部激勵(lì)的主要原因。因此，為了研究齒輪內(nèi)部激勵(lì)對(duì)齒輪系統(tǒng)動(dòng)態(tài)響應(yīng)的影響，對(duì)整個(gè)齒輪嚙合過(guò)程進(jìn)行分析，計(jì)算得到嚙合過(guò)程中嚙合剛度的變化規(guī)律是很有必要的，這是齒輪故障診斷、系統(tǒng)動(dòng)態(tài)特性分析以及壽命預(yù)測(cè)等研究的基礎(chǔ)。

隨著計(jì)算機(jī)的快速發(fā)展，采用有限元方法計(jì)算齒輪時(shí)變嚙合剛度比較普遍。通過(guò)建立齒輪對(duì)三維接觸有限元分析模型可得到嚙合輪齒間的時(shí)變剛度曲線[2-3]。但建立有限元模型需要考慮的因素較多，模型比較復(fù)雜，而且計(jì)算不同齒輪的嚙合剛度需要建立不同的有限元模型，效率比較低。而用數(shù)值計(jì)算方法求輪齒的嚙合剛度則省去了建模的過(guò)程，在已有的數(shù)學(xué)模型的基礎(chǔ)上，通過(guò)改變參數(shù)數(shù)值即可求解，計(jì)算速度快，效率高，精度也滿足要求。數(shù)值求解方法主要有日本的石川公式法[4]和Weber能量法。Weber能量法在國(guó)外應(yīng)用較多[5-6]，國(guó)內(nèi)也有學(xué)者應(yīng)用Weber能量法對(duì)齒輪變形量進(jìn)行推導(dǎo)[7-8]，但石川公式因固定系數(shù)較少、計(jì)算過(guò)程簡(jiǎn)單、運(yùn)算量小、適用范圍更廣而更多地被應(yīng)用[9]。

目前，隨著汽車工業(yè)快速發(fā)展，尤其是新能源汽車的發(fā)展，行星輪系因體積小、質(zhì)量輕、傳動(dòng)比大等優(yōu)點(diǎn)，越來(lái)越多地被應(yīng)用在汽車自動(dòng)變速器中。行星輪系在AT變速器中可與離合器、制動(dòng)器搭配作為主變速構(gòu)件，在CVT自動(dòng)變速器中可作為副變速構(gòu)件增大傳動(dòng)比范圍，在新能源汽車上可使用ECVT自動(dòng)變速器與電動(dòng)機(jī)共同實(shí)現(xiàn)變速。行星輪系齒輪嚙合與定軸輪系齒輪嚙合不同，行星輪系中行星架帶動(dòng)行星輪公轉(zhuǎn)使得行星輪與太陽(yáng)輪、行星輪與內(nèi)齒圈嚙合線不斷變化，外嚙合與內(nèi)嚙合同時(shí)存在。在現(xiàn)有的文獻(xiàn)中，利用有限元方法求解行星輪系時(shí)變嚙合剛度的研究較多，而在利用數(shù)值計(jì)算方法求解行星輪系時(shí)變嚙合剛度的文獻(xiàn)中以標(biāo)準(zhǔn)行星輪與太陽(yáng)輪外嚙合剛度求法為主，基本上沒(méi)考慮變位齒輪外嚙合的情況，更沒(méi)有考慮變位行星輪與內(nèi)齒圈內(nèi)嚙合的嚙合剛度。鑒于上述情況，本文對(duì)正變位行星輪與太陽(yáng)輪外嚙合進(jìn)行分析，并在此基礎(chǔ)上推導(dǎo)了正變位行星輪與內(nèi)齒圈的內(nèi)嚙合剛度。

1石川公式

在齒輪嚙合傳動(dòng)過(guò)程中，嚙合齒對(duì)數(shù)與嚙合接觸點(diǎn)位置的變化會(huì)導(dǎo)致齒輪綜合嚙合剛度隨時(shí)間周期性變化。輪齒間嚙合作用可以簡(jiǎn)化為嚙合線方向上的時(shí)變彈簧，某一接觸點(diǎn)的嚙合剛度即為此時(shí)的彈簧剛度。GB/T 3480—1997中提出了齒輪嚙合剛度的計(jì)算方法，將齒輪嚙合剛度定義為使一對(duì)或幾對(duì)同時(shí)嚙合輪齒在1 mm齒寬上產(chǎn)生1 μm變形量所需要的嚙合線上的載荷力。齒輪嚙合剛度的數(shù)學(xué)表達(dá)式為[11]

(1)

式中：Fn為嚙合線法向載荷；b為輪齒的齒寬；δ為輪齒在嚙合線方向上單位齒寬綜合彈性變形量。

嚙合剛度為輪齒的固有特性，與所受力的大小無(wú)關(guān)，通過(guò)齒輪嚙合每一時(shí)刻單位力作用下輪齒的綜合變形量即可求得此時(shí)刻的齒輪嚙合剛度。本文采用石川公式法計(jì)算齒輪嚙合線方向上輪齒嚙合的綜合彈性變形量，并在此基礎(chǔ)上求輪齒嚙合剛度。

1.1標(biāo)準(zhǔn)齒輪石川公式

石川公式作為一種常見(jiàn)的齒輪輪齒綜合彈性變形的計(jì)算方法，將復(fù)雜的齒輪漸開(kāi)線齒廓簡(jiǎn)化為由矩形與梯形組合而成的當(dāng)量齒形結(jié)構(gòu)，矩形底部所在的圓為齒輪的齒根圓，矩形與漸開(kāi)線齒廓的交點(diǎn)為嚙合齒輪對(duì)中主動(dòng)輪嚙合開(kāi)始點(diǎn)、被動(dòng)輪嚙合結(jié)束點(diǎn)，該點(diǎn)在漸開(kāi)線齒廓上所對(duì)應(yīng)的圓稱為有效齒根圓，簡(jiǎn)化后的外齒輪齒廓石川公式當(dāng)量齒形如圖1所示。

圖1　外齒輪齒廓石川公式當(dāng)量齒形

由圖1中幾何關(guān)系可得外齒輪齒廓石川公式當(dāng)量齒形幾何參數(shù)[12]。

1) 齒頂圓齒厚

(2)

式中：da為齒頂圓直徑；r為分度圓半徑；s為分度圓齒厚；α為漸開(kāi)線齒廓分度圓壓力角；αa為漸開(kāi)線齒廓齒頂圓壓力角。

2) 有效齒根圓齒厚

(3)

式中：dF為有效齒根圓直徑；αF為漸開(kāi)線齒廓有效齒根圓壓力角。

3) 齒根圓齒厚

(4)

式中：df為齒根圓直徑；αf為漸開(kāi)線齒廓齒根圓壓力角。

4) 矩形高

(5)

5) 齒形高

(6)

6) 嚙合點(diǎn)所在齒輪圓半個(gè)齒厚角

(7)

式中αx為嚙合點(diǎn)所在齒輪圓壓力角。

7) 載荷作用角

(8)

8) 嚙合點(diǎn)高

(9)

式中dx為嚙合點(diǎn)所在齒輪圓直徑。

9) 輔助尺寸

(10)

將輪齒簡(jiǎn)化之后，在求解齒輪嚙合變形量之前進(jìn)行一定的假設(shè)：

1) 在外嚙合直齒圓柱齒輪傳動(dòng)中將齒輪輪體視為剛體，輪體部分變形為0；

2) 在內(nèi)嚙合直齒圓柱齒輪傳動(dòng)中內(nèi)齒圈輪體為很薄的一部分，可近似忽略輪體變形的影響；

3) 忽略齒輪系統(tǒng)傳動(dòng)過(guò)程中的熱變形。

此時(shí)，輪齒嚙合綜合彈性變形依賴于3部分的變形量：① 輪齒漸開(kāi)線齒廓的彎曲變形δB(矩形部分的彎曲變形量δBr+梯形部分的彎曲變形量δBt)、嚙合剪切變形量δs；② 輪齒齒根與輪體相連過(guò)渡曲線圓角部分的變形量δt；③ 由齒輪對(duì)嚙合接觸應(yīng)力引起的局部變形量δc。

嚙合齒輪對(duì)中單個(gè)輪齒變形量用數(shù)學(xué)公式表示為：

(11)

嚙合齒輪對(duì)中1對(duì)輪齒的綜合彈性變形量用數(shù)學(xué)公式表示為

(12)

式中：δp為主動(dòng)輪輪齒嚙合變形量；δg為從動(dòng)輪輪齒嚙合變形量；δc為輪齒接觸應(yīng)力引起的局部變形量。

根據(jù)簡(jiǎn)化圖形的幾何參數(shù)求解各個(gè)變形量的具體公式見(jiàn)文獻(xiàn)[4]。

1.2變位齒輪石川公式

在齒輪傳動(dòng)系統(tǒng)中，為了改善標(biāo)準(zhǔn)齒輪的不足，突破標(biāo)準(zhǔn)齒輪的限制，通常對(duì)齒輪進(jìn)行修正?，F(xiàn)在最廣泛應(yīng)用的是通過(guò)改變刀具與齒輪輪坯的相對(duì)位置來(lái)切制齒輪，此時(shí)刀具的分度線與輪坯的分度圓不再相切，切制成的齒輪分度圓上齒厚s與齒槽寬e不再相等，這已不是標(biāo)準(zhǔn)齒輪，稱其為變位齒輪。變位齒輪與標(biāo)準(zhǔn)齒輪相比，模數(shù)、分度圓壓力角、分度圓直徑、齒距并沒(méi)有發(fā)生變化，發(fā)生變化的為分度圓齒厚、齒頂圓直徑、齒根圓直徑。

正變位外齒輪分度圓齒厚為

(13)

正變位內(nèi)齒輪分度圓齒厚為

(14)

式中：m為模數(shù)；x為徑向變位系數(shù)。

2行星輪系嚙合剛度

2.1行星輪與太陽(yáng)輪外嚙合剛度

在角度變位齒輪傳動(dòng)中，正傳動(dòng)可以減小齒輪機(jī)構(gòu)的尺寸，較大地提高齒輪機(jī)構(gòu)的承載能力，在體積相對(duì)較小和承載能力高的行星輪系中應(yīng)用較多。當(dāng)行星輪與太陽(yáng)輪均為正變位齒輪時(shí)，外嚙合過(guò)程中分度圓與節(jié)圓不再重合，此時(shí)要保證齒輪傳動(dòng)時(shí)無(wú)側(cè)隙正常嚙合，其中一齒輪在節(jié)圓上的齒槽寬要等于另一齒輪在節(jié)圓上的齒厚，由此可得：

(15)

式中：α′為變位后的漸開(kāi)線齒廓節(jié)圓壓力角；xs為太陽(yáng)輪變位系數(shù)；xp為行星輪變位系數(shù)。

正變位行星輪與正變位太陽(yáng)輪外嚙合齒輪傳動(dòng)的齒輪幾何尺寸計(jì)算公式如表1所示。根據(jù)表1齒輪幾何尺寸計(jì)算公式，可以得到按石川公式簡(jiǎn)化后的當(dāng)量齒形參數(shù)中除嚙合點(diǎn)高h(yuǎn)x之外的其他參數(shù)。在外嚙合過(guò)程中，齒輪副嚙合點(diǎn)不斷發(fā)生變化。當(dāng)嚙合點(diǎn)高h(yuǎn)x隨之發(fā)生變化，并且嚙合點(diǎn)確定后可得到嚙合點(diǎn)高h(yuǎn)x，即可利用石川公式的當(dāng)量齒形各部分變形量計(jì)算公式求解出嚙合過(guò)程中綜合彈性變形量，然后通過(guò)嚙合剛度表達(dá)式得到在一個(gè)嚙合周期內(nèi)齒輪副嚙合剛度隨時(shí)間變化的曲線。嚙合周期為

(16)

式中：nh為行星架轉(zhuǎn)速；np為行星輪轉(zhuǎn)速；zp為行星輪齒數(shù)。

表1　行星輪與太陽(yáng)輪幾何尺寸計(jì)算公式

2.2行星輪與內(nèi)齒輪內(nèi)嚙合剛度

內(nèi)嚙合是內(nèi)齒輪與外齒輪進(jìn)行嚙合的傳動(dòng)，外齒輪是圓盤(pán)形的，外圈有齒,而內(nèi)齒輪則是環(huán)形的，外圈光滑，輪齒分布在輪圈的內(nèi)圓上，內(nèi)齒輪齒廓如圖2所示。

由圖2內(nèi)齒輪齒廓可知：內(nèi)齒輪齒形與外齒輪齒形不同，內(nèi)齒輪的輪齒與外齒輪齒槽部分相當(dāng)，內(nèi)齒輪的齒頂圓和齒根圓相當(dāng)于外齒輪的齒根圓與齒頂圓。內(nèi)齒輪輪齒簡(jiǎn)化為梯形與矩形組合成的石川公式當(dāng)量齒形如圖3所示。

圖2　內(nèi)齒輪齒廓

圖3　內(nèi)齒輪輪齒石川公式的當(dāng)量齒形

此時(shí)，內(nèi)齒輪輪齒石川公式當(dāng)量齒形中幾何參數(shù)為:

齒頂圓齒厚

(17)

齒根圓齒厚

(18)

有效齒根圓齒厚

(19)

嚙合點(diǎn)所在齒輪圓半個(gè)齒厚角

(20)

嚙合點(diǎn)處法向載荷與輪齒中心線垂線的夾角

(21)

由于內(nèi)齒輪的齒頂圓和齒根圓的位置剛好與外齒輪相反，因此簡(jiǎn)化后的內(nèi)齒輪齒廓當(dāng)量齒形中矩形高h(yuǎn)r、齒形高h(yuǎn)與嚙合點(diǎn)高h(yuǎn)x的計(jì)算公式與式(5)、(6)、(9)相差一個(gè)負(fù)號(hào)，而輔助尺寸hi的計(jì)算公式與式(10)相同。

對(duì)于正變位的行星輪與內(nèi)齒圈內(nèi)嚙合，無(wú)側(cè)隙正常嚙合傳動(dòng)時(shí)嚙合角計(jì)算公式為

(22)

實(shí)際中心距a′及中心距變動(dòng)系數(shù)y計(jì)算公式與外嚙合時(shí)相同，而齒頂高降低系數(shù)不同，其計(jì)算公式為

(23)

由于正變位行星輪與內(nèi)齒圈內(nèi)嚙合傳動(dòng)時(shí)行星輪幾何尺寸計(jì)算方法與外嚙合時(shí)相同，因此僅列出內(nèi)齒輪的幾何尺寸計(jì)算公式如表2所示。

表2　內(nèi)齒圈的幾何尺寸計(jì)算公式

注：xr為內(nèi)齒圈變位系數(shù)。

對(duì)行星輪與內(nèi)齒圈內(nèi)嚙合過(guò)程進(jìn)行分析，畫(huà)出齒輪嚙合過(guò)程簡(jiǎn)圖，如圖4所示。圖4中：齒輪1為行星輪；齒輪2為內(nèi)齒輪；N1N2為嚙合過(guò)程的理論嚙合線，在內(nèi)嚙合時(shí)其為兩嚙合齒輪基圓的外公切線。與外嚙合時(shí)理論嚙合線為兩齒輪的內(nèi)公切線所不同，在內(nèi)嚙合過(guò)程中，行星輪為主動(dòng)輪，以行星輪有效齒根圓與內(nèi)齒圈齒頂圓接觸嚙合點(diǎn)B2開(kāi)始嚙合，行星輪的凸齒面與內(nèi)齒圈的凹齒面逐漸接觸，以行星輪齒頂圓與內(nèi)齒圈有效齒根圓接觸嚙合點(diǎn)B1結(jié)束嚙合，B1B2為實(shí)際嚙合線。與齒輪外嚙合時(shí)實(shí)際嚙合線在理論嚙合線內(nèi)不同，內(nèi)嚙合的實(shí)際嚙合線在理論嚙合線的延長(zhǎng)線上。在實(shí)際嚙合線B1B2中，B2C、B1D為雙齒嚙合區(qū)，CD為單齒嚙合區(qū)。在單雙齒嚙合區(qū)內(nèi)利用石川公式的當(dāng)量齒形各部分變形量計(jì)算公式求得綜合變形量，從而得到行星輪與內(nèi)齒輪內(nèi)嚙合時(shí)一個(gè)嚙合周期內(nèi)的嚙合剛度曲線，其嚙合周期計(jì)算公式與式(16)相同。

圖4　行星輪與內(nèi)齒輪內(nèi)嚙合過(guò)程簡(jiǎn)圖

3算例與結(jié)果分析

3.1行星輪系嚙合剛度

一個(gè)試制的等比例縮小風(fēng)電行星齒輪箱正變位行星輪系各齒輪參數(shù)如表3所示。

表3　行星輪系齒輪參數(shù)

用Matlab編程計(jì)算，設(shè)內(nèi)齒圈固定，當(dāng)輸入軸行星架轉(zhuǎn)速為600r/min時(shí)，計(jì)算得到行星輪相對(duì)行星架的轉(zhuǎn)速為1 577.78r/min。行星輪與太陽(yáng)輪外嚙合剛度時(shí)頻圖如圖5所示。外嚙合的重合度為1.187 8，雙齒嚙合區(qū)占嚙合周期比例為0.187 8。行星輪與內(nèi)齒圈內(nèi)嚙合剛度時(shí)頻圖如圖6所示。內(nèi)嚙合的重合度為1.527，雙齒嚙合區(qū)占嚙合周期比例為0.527。

圖5　外嚙合剛度時(shí)頻

圖6　內(nèi)嚙合剛度時(shí)頻

3.2與Weber能量法計(jì)算結(jié)果對(duì)比

根據(jù)Weber能量法，通過(guò)Matlab編程得到齒輪嚙合過(guò)程中時(shí)變嚙合剛度。因Weber能量法未對(duì)輪齒進(jìn)行簡(jiǎn)化，計(jì)算結(jié)果精度相對(duì)較高，ISO齒輪剛度計(jì)算標(biāo)準(zhǔn)也是以Weber能量法公式為基礎(chǔ)的，故以Weber能量法計(jì)算結(jié)果為標(biāo)準(zhǔn)，對(duì)石川公式計(jì)算得到的靜剛度和前4階動(dòng)剛度幅值進(jìn)行比較，比較數(shù)據(jù)如表4和表5所示。

表4　行星輪與太陽(yáng)輪外嚙合剛度

注：數(shù)據(jù)四舍五入，取小數(shù)點(diǎn)后4位

表5　行星輪與內(nèi)齒圈內(nèi)嚙合剛度

注：數(shù)據(jù)四舍五入，取小數(shù)點(diǎn)后4位

數(shù)據(jù)分析結(jié)果表明：內(nèi)嚙合剛度相對(duì)誤差整體要比外嚙合要小，這是由于內(nèi)嚙合時(shí)內(nèi)齒輪和外齒輪的當(dāng)量齒形簡(jiǎn)化形成的誤差相互補(bǔ)償與抵消，石川公式計(jì)算結(jié)果更接近實(shí)際值。而無(wú)論內(nèi)嚙合與外嚙合，其2、4階嚙合剛度幅值均比1、3階嚙合剛度幅值小一個(gè)數(shù)量級(jí)，相對(duì)誤差要大些，但總體誤差均在可以接受的范圍內(nèi)，可以驗(yàn)證上述分析中利用石川公式計(jì)算變位行星輪系嚙合剛度的正確性與可行性。

4結(jié)束語(yǔ)

基于常見(jiàn)的外嚙合輪齒石川公式當(dāng)量齒形，詳細(xì)推導(dǎo)了行星齒輪嚙合的正變位外齒輪齒形幾何參數(shù)的計(jì)算方法。根據(jù)內(nèi)齒輪齒形與外齒輪齒形不同，做出內(nèi)齒輪石川公式當(dāng)量齒形，推導(dǎo)出正變位內(nèi)齒輪齒形幾何參數(shù)計(jì)算公式。利用石川公式計(jì)算方法，通過(guò)Matlab編程求解行星輪系內(nèi)嚙合與外嚙合綜合彈性變形量，得到行星輪系時(shí)變嚙合剛度曲線，并將計(jì)算結(jié)果與Weber能量法計(jì)算結(jié)果進(jìn)行對(duì)比分析。結(jié)果表明：靜態(tài)剛度最大誤差為3.7%，不考慮小幅值動(dòng)態(tài)剛度時(shí)前4階諧波幅值最大誤差約為9.1%，能滿足精度要求，驗(yàn)證了本文所提計(jì)算方法的正確性，可為行星輪系動(dòng)力學(xué)特性分析提供較準(zhǔn)確的嚙合剛度，較有限元法更為有效和方便，計(jì)算效率更高。

參考文獻(xiàn)：

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

[2]LIR,YANGC,LINT,etal.Finiteelementsimulationofthedynamicalbehaviorofaspeed-increasegearbox[J].Journalofmaterialsprocessingtechnology,2004,150(1):170-174.

[3]DELRINCONAF,VIADEROF,IGLESIASM,etal.Amodelforthestudyofmeshingstiffnessinspurgeartransmissions[J].MechanismandMachineTheory,2013,61:30-58.

[4]日本機(jī)械學(xué)會(huì).齒輪強(qiáng)度設(shè)計(jì)資料[M].李茹貞,趙清慧,譯.北京：機(jī)械工業(yè)出版社,1984:28-35.

[5]TAVAKOLIMS,HOUSERDR.Optimumprofilemodificationsfortheminimizationofstatictransmissionerrorsofspurgears[J].JournalofMechanicalDesign,1986,108(1):86-94.

[6]YANGDCH,LinJY.Hertziandamping,toothfrictionandbendingelasticityingearimpactdynamics[J].JournalofMechanicalDesign,1987,109(2):189-196.

[7]楊長(zhǎng)輝,徐濤金,許洪斌,等.基于Weber能量法的直齒輪時(shí)變嚙合剛度數(shù)值計(jì)算[J].機(jī)械傳動(dòng),2015,39(2):59-62.

[8]顏海燕,唐進(jìn)元,宋紅光.直齒輪輪齒變形計(jì)算的數(shù)值積分法[J].機(jī)械傳動(dòng),2005,29(2):7-9.

[9]李亞鵬,孫偉,魏靜,等.齒輪時(shí)變嚙合剛度改進(jìn)計(jì)算方法[J].機(jī)械傳動(dòng),2010,34(5):22-26.

[10]陳再剛,邵毅敏.行星輪系齒輪嚙合非線性激勵(lì)建模和振動(dòng)特征研究[J].機(jī)械工程學(xué)報(bào),2015(7):23.

[11]GB/T3480—1997.漸開(kāi)線圓柱齒輪承載能力計(jì)算方法[S].中國(guó):國(guó)家技術(shù)監(jiān)督局,1997.

[12]黃祝慶.基于多體動(dòng)力學(xué)的齒輪副建模與實(shí)驗(yàn)研究[D].南京：南京航空航天大學(xué),2013.

(責(zé)任編輯楊文青)

Numerical Calculation of Planetary Gear Train Involute ModifiedGear Time-Varying Meshing Stiffness

YUE Xi-zheng, DING Wen-si, DING Kang, ZENG Zhi-jie

(School of Mechanical and Automotive Engineering,South China University of Technology, Guangzhou 510641, China)

Abstract:The time-varying gear meshing stiffness is the main cause of gear vibration. To caculate the time-varying meshing stiffness has vital significance. Based on geometric parameters calculation method of external gear tooth profile in positive modified planetary gear train, according to the differences between external and internal gear tooth profiles, the geometric parameters calculation method of internal gear tooth profile was deduced in detail. Using Ishikawa tooth deformation calculation formula and stiffness calculation formula, positive modified planetary gear train internal and external time-varying meshing stiffness curves can be got by Matlab programming. Compared with the calculation result of Weber energy method, maximum static stiffness error is 3.7% and the first four order harmonic amplitude maximum error is 9.1% without considering the small amplitude value of the dynamic stiffness. Calculation accuracy meets the requirements. It can provide relatively accurate meshing stiffness for motion differential equation of the planetary gear train.

Key words:modified gear; Ishikawa formula; external meshing stiffness; internal meshing stiffness

收稿日期：2016-03-09

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

作者簡(jiǎn)介：岳喜錚(1989—)，男，碩士研究生，主要從事行星輪系動(dòng)力學(xué)特性分析研究； 通訊作者 丁問(wèn)司(1968—)，男，教授，主要從事機(jī)械系統(tǒng)動(dòng)力學(xué)研究。

doi:10.3969/j.issn.1674-8425(z).2016.06.007

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

文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：1674-8425(2016)06-0038-07

引用格式：岳喜錚，丁問(wèn)司，丁康， 等.行星輪系漸開(kāi)線變位齒輪時(shí)變嚙合剛度數(shù)值計(jì)算[J].重慶理工大學(xué)學(xué)報(bào)(自然科學(xué))，2016(6):38-44.

Citation format：YUE Xi-zheng, DING Wen-si, DING Kang, et al.Numerical Calculation of Planetary Gear Train Involute Modified Gear Time-Varying Meshing Stiffness[J].Journal of Chongqing University of Technology(Natural Science)，2016(6):38-44.