擺線針輪傳動線接觸彈流潤滑分析

孫章棟，朱才朝，劉懷舉，劉明勇，顧宗琳

（重慶大學機械傳動國家重點實驗室，重慶 400030）

擺線針輪傳動線接觸彈流潤滑分析

孫章棟，朱才朝，劉懷舉，劉明勇，顧宗琳

（重慶大學機械傳動國家重點實驗室，重慶 400030）

對擺線針輪線接觸彈流潤滑進行數值分析，得出嚙合過程中最小膜厚，與經驗公式對比，驗證模型的正確性。以擺線輪θ＝arccosK1處潤滑狀態作為判斷依據，研究載荷、轉速和流變指數對擺線針輪傳動潤滑狀態的影響規律。結果表明：隨著轉速提高，二次壓力峰高度降低并向入口區移動，膜厚相應增加；而隨著載荷增加，接觸區變寬，二次壓力峰增加且向出口區移動，膜厚略微減??；流變指數n增加，二次壓力峰先增加后減小，最后趨近于Hertz壓力，并向出口區移動，膜厚相應減小。討論了短幅系數k1對潤滑的影響，表明在滿足設計要求情況下，短幅系數k1減小，有利于提高潤滑性能。

擺線針輪；脂潤滑；數值解；Ostwald

擺線針輪行星傳動以其結構緊湊、傳動比大、傳動精度高和承載能力大等優點，而被廣泛應用。通常使用脂來實現潤滑。漸開線齒輪潤滑問題已進行廣泛的研究，但有關擺線針輪潤滑研究少有報道。擺線針輪潤滑不同于漸開線：首先擺線針輪傳動嚙合原理以及幾何特征有別于漸開線傳動；其次齒擺線輪與針輪同時嚙合齒數多于漸開線，同時嚙合齒數理論上可以達到一半擺線輪齒數；最后漸開線使用油進行潤滑，脂的流變特性有別于油。彈流潤滑性能的好壞對擺線針輪傳動的性能和使用壽命都有顯著影響。開展擺線針輪傳動彈流脂潤滑數值分析具有重要的意義。

國內外學者對齒輪彈流潤滑做了大量研究。Liu［2－3］對線接觸直齒輪的混合潤滑模型，討論了轉速及不同粗糙度對直齒輪潤滑的影響；Zhu等［4］討論了不同轉速工況下由混合潤滑過渡到干摩擦接觸的潤滑情況；Chu等［5］研究了接觸面的液面滑移和流體特性對潤滑的影響。Yang等［6］對斜齒輪進行彈流潤滑分析，討論了螺旋角對潤滑的影響；Liu等［7］首次分析了有限長線接觸的熱解，并與無限長熱解作對比；Zhu等［8］采用有限長線接觸模型分析斜齒輪穩態熱彈流潤滑問題，討論載荷與轉速工況變化對斜齒輪潤滑溫度及齒面剪切力的影響。劉曉玲等［9］展開對指數率非牛頓流體在線接觸條件下的彈流潤滑數值研究，應用多重網格法和多重網格積分法數值求解了Newton流體和Ree-Eyring流體線接觸等溫和熱彈流潤滑問題。于玫等［10］探討了不同工況下熱效應對脂潤滑彈流數值解的影響。李媛等［11］探討脂潤滑條件下激光加工微凹坑表面的摩擦特性，表明適當的凹坑深度可以改善表面的摩擦特性和溫度特性。白新瑞等［12］基于有限長線接觸彈流潤滑理論，求得乏油條件下圓柱滾子軸承彈流潤滑的完全數值解。楊沛然等［13］做了大量脂潤滑方面的研究工作。

Mihailidis等［14］推導出適合非牛頓流體的乏油廣義雷諾方程。由于擺線針輪傳動通常是用潤滑脂實現潤滑，論文以Ostwald模型對擺線針輪傳動線接觸彈流脂潤滑問題進行數值分析研究，研究載荷、轉速和流變指數對擺線針輪傳動潤滑狀態的影響規律，討論設計參數短幅系數k1對潤滑的影響。

1 數學模型

1.1 嚙合分析

擺線針輪傳動同時嚙合的齒數多，如圖1所示。在嚙合過程中，擺線輪曲率半徑ρ變化，針齒半徑rz為定值，從而可得嚙合過程中的當量曲率半徑為：

式中“＋”號為用于外凸齒廓嚙合；“－”號為用于內凹齒廓嚙合；ρ為有關嚙合相位角的函數。

圖1 擺線輪與針齒嚙合示意圖Fig．1 Mathematical model of the cycloid drive

由運動學分析可知，擺線輪與針齒在嚙合點處的線速度是相等的，設u1為擺線輪在嚙合點處的線速度，u2為針齒在嚙合點處的線速度，由圖1所示u1＝u2＝um＝，對潤滑油的卷吸速度：

1.2 潤滑模型

基于Ostwald模型的潤滑脂本構方程為：

式中：m為粘性函數；n為流變指數

采用廣義Reynolds方程［15］，引入等效粘度η*，對于Ostwald模型，其等效粘度為：

引入等效粘度可將潤滑脂的非牛頓流體轉換成牛頓流體，方便采用牛頓流體解法求解。利用Hertz接觸參數對潤滑模型進無量綱化，無量綱量參考量：

式中x為流動方向，z為膜厚方向，b為Hertz接觸半寬，h為油膜厚度，p為油膜壓力，PR為最大Hertz壓力，R為等效曲率半徑，η0為環境粘度，E為等效彈性模量。

無量綱Reynolds方程

式中ρ0為潤滑油環境密度。

無量綱Reynolds方程的邊界條件為：

式中H0為常數。

無量綱密度方程

無量綱載荷方程

2 線接觸彈流脂潤滑分析

將式（7）～式（9）離散化，建立擺線針輪傳動彈流脂潤滑數值仿真模型，通過數值直接迭代法對彈流潤滑中的膜厚和壓力進行數值求解。Gauss-Seidel迭代法用于求解低壓區壓力，Jacobi雙極子迭代法用于求解高壓區壓力［16］。

壓力求解無量綱計算域范圍X＝［XinXout］＝［－4，1．5］，油膜厚度無量綱計算范圍Z＝［0，1］。X方向的節點數1 025，油膜厚度Z向節點數11。數值迭代精度為：壓力相對誤差小于10－6。

選取擺線針輪行星傳動嚙合副為研究對象，其嚙合副與潤滑相關參數如表1。

表1 擺線針輪傳動副與潤滑相關參數Tab.1 The cycloid drive parameters and properties of the lubricant

圖2 單位嚙合力隨嚙合相位變化Fig．2 The variation of contact force with meshing process

圖3 等效曲率半徑和卷吸速度變化曲線Fig．3 The variation of reduced curve radius and rolling speed

擺線輪與針齒滾動線接觸嚙合，在嚙合過程中，單位嚙合力F、等效曲率半徑和卷吸速度隨著嚙合相位變化。圖2給出單位嚙合力隨嚙合相位變化曲線，圖3給出卷吸速度和等效曲率半徑隨嚙合相位變化曲線。

設定輸入轉速500 r／min，輸出力矩420 N／m。選取某一嚙合相位處進行分析，其潤滑脂膜壓力、膜厚如圖4所示。在出口區處，潤滑膜壓力出現明顯二次壓力峰，在二次壓力峰相對應處，潤滑膜開始收縮，形成出口區的頸縮現象，頸縮處的膜厚為最小潤滑膜厚度。擺線輪與針齒滾動嚙合剪切應力如圖5所示，在潤滑膜出口區的頸縮位置出現剪應力峰值，其它區域切應力趨于零值，摩擦小，進一步說明溫度對潤滑影響小。

擺線輪與針齒理論嚙合區間為0°～180°，而實際有效嚙合區間為30°～120°。計算最大壓力與最小膜厚隨著嚙合相位變化曲線如圖6。如圖7所示，與采用Dowson-Higginson潤滑油膜厚經驗公式計算結果變化趨勢相吻合，但比經驗公式偏大。由于Dowson-Higginson經驗公式用來計算潤滑油最小膜厚的，而修形齒廓最小膜厚數值解用來計算潤滑脂的，驗證數值模型仿真結果正確性。從圖2、圖3和圖7可以得出在嚙合相位角θ＝arccosK1＝42．97°處為曲線拐點，在嚙合過程中相位角θ＝arccosK1＝42．97°時油膜厚度最小，故以此相位處潤滑狀態作為整個嚙合過程中潤滑狀態的判斷依據。

圖5 擺線輪與針齒滾動嚙合切應力Fig．5 Cycloid gear and pin gear mesh rolling shear stress

圖6 中心壓力與最小膜厚隨嚙合相位變化Fig．6 The variation of center-pressure and min-film thickness

圖7 最小膜厚數值解與經驗公式比較Fig．7 Compare film thickness numerical solution with empirical formula

3 參數變化對彈流脂潤滑的影響

3.1 工況參數變化對潤滑影響

圖8～9給出了轉速變化對θ＝42．97°處潤滑性能的變化規律。隨著轉速提高，入口區壓力增加，二次壓力峰高度減小，其位置向入口區移動，膜厚相應增加，潤滑膜平行部分縮短，頸縮變得不明顯，表明擺線輪與針齒嚙合時，卷吸速度增大，有利于形成良好的彈流潤滑膜。

圖10～11給出了載荷變化對θ＝42．97°處潤滑性能的變化規律。隨著載荷增加，接觸區變寬，入口區壓力增加，二次壓力峰增加，其位置向出口區移動，膜厚略微減小，膜厚頸縮現象更加明顯。主要由于隨著載荷增加，油膜剛度變化較小，故膜厚變化不明顯。

3.2 流變指數n對潤滑影響

圖12、13是流變指數n對θ＝42．97°處潤滑性能的影響規律。當n＜1時，隨著流變指數n增加，入口區壓力降低，二次壓力峰增加，并向出口區移動，膜厚相應降低；當n＞1時，隨著流變指數n增加，入口區壓力略微降低，二次壓力峰減小并向出口區移動，最后二次壓力峰消失，接近Hertz壓力，膜厚相應減小，頸縮也越來越不明顯，最后消失，潤滑狀態變差。

圖8 轉速對油膜壓力分布的影響Fig．8 Pressure distribution for different rotational Speeds

圖9 轉速對油膜厚度的影響Fig．9 Film shape for different rotational Speeds

圖10 載荷對油膜壓力的影響Fig．10 Pressure distribution for different loads

圖11 載荷對油膜厚度的影響Fig．11 Film shape for different loads

圖12 流變指數n對油膜壓力的影響Fig．12 Pressure distribution for different n

圖13 流變指數n對膜厚的影響Fig．13 Film shape for different n

3．3 短幅系數k1對潤滑的影響

圖14～15給出了短幅系數k1對擺線輪平均半徑處潤滑性能的變化規律。隨著短幅系數k1增加，在大部分嚙合區間范圍內壓力增加，承載能力下降，膜厚減小。表明在滿足設計要求情況下，短幅系數k1減小，有利于提高潤滑性能。

圖14 中心壓力隨短幅系數k1變化Fig．14 The variation of center-pressure for different k1

圖15 最小膜厚隨著短幅系數k1變化Fig．15 The variation of min-film thickness for different k1

4 結 論

論文基于Ostwald模型對擺線針輪線接觸脂潤滑彈流問題進行研究，得出嚙合過程中最小膜厚值，并與經驗公式對比，驗證數學模型的正確性。以擺線輪θ＝arccos K1處潤滑狀態作為判斷依據，研究載荷、轉速和流變指數對擺線針輪傳動滑狀態的影響規律，表明隨著轉速提高，二次壓力峰高度降低并向入口區移動，膜厚相應增加，潤滑膜平行部分縮短；而隨著載荷增加，接觸區變寬，二次壓力峰增加，其位置向出口區移動，膜厚略微減小，膜厚頸縮現象更加明顯；隨著流變指數n變化，當n＜1時，隨著流變指數n增加，二次壓力峰增加，并向出口區移動，膜厚相應降低；當n＞1時，隨著流變指數n增加，二次壓力峰減小并向出口區移動，最后二次壓力峰消失，接近Hertz壓力，膜厚相應減小。討論了短幅系數k1對潤滑的影響，隨著短幅系數k1增加，在大部分嚙合區間范圍內壓力增加，承載能力下降，膜厚減小。在滿足設計要求情況下，短幅系數k1減小，有利于提高潤滑性能。

［1］朱才朝，陳爽，馬飛，等．輪齒修形對兆瓦級風電齒輪箱動態特性影響［J］．振動與沖擊，2013，32（7）：123－128．

ZHU Cai-chao，CHEN Shuang，MA Fei，et al．Effect of gear teeth modification on dynamic characteristics of a megawatt level wind turbine gearbox［J］．Journal of Vibration and Shock，2013，32（7）：123－128．

［2］Liu H J，Mao K，Zhu C C，et al．Mixed lubricated line contact analysis for spur gears using a deterministic model［J］．ASME Journal of Tribology，2012，134：021505－1－7．

［3］Liu H J，Mao K，Zhu C C，et al．Parametric studies of spur gear lubrication performance considering dynamic loads［J］．Proc．IMechE，PartJ：Journal of Engineering Tribology，2012，226（9）：731－737．

［4］Zhu D，Wang J X，Ren N，et al．Mixed Elastohydrodynamic lubrication in finite roller contacts involving realistic geometry and surface roughness［J］．ASME，Journal of Tribology，2012，134：011504－1－10．

［5］Chu L M，Lin J R，Li W L，et al．Model for line contact EHL problems consideration of effects of navierSlip and lubricant rheology［J］．ASME，Journal of Tribology，2012，134：031502 －1－8．

［6］Yang P，Yang P R．Analysis on the thermal elastohydrodynamic lubrication of tapered rollers in opposite orientation ［J］．Tribology International，2007，40：1627－1637．

［7］Liu X L，Yang P R．Analysis of the themal elastohydrodyanamic lubrication of a finite line contact［J］．Tribology International，2002，35：137－144．

［8］Zhu C C，Liu M Y，Liu，H J，et al．A thermal finite line contact EHL model of a helical gear pair［J］．Proc IMechE Part J：J Engineering Tribology，2013，227（4）299－309．

［9］劉曉玲，楊沛然．指數率流體熱彈流潤滑分析［J］．潤滑與密封，2007（8）：19－23．

LIU Xiao-lin，YANG Pei-ran．Analysis of power law fluid thermal elastohydrodynamic lubrication［J］．Lubrication engineering，2007（8）：19－23．

［10］于玫，黃平．線接觸脂潤滑熱彈性流體動力潤滑數值分析［J］．軸承，2011（1）：8－12

YU Mei，HUANG Ping．Numerical analysis on thermal elastohydrodynamic grease lubrication in line contacts［J］．Bearing，2011（1）：8－12．

［11］李媛，劉小君，王靜，等．脂潤滑條件下微凹坑表面的摩擦特性研究［J］．合肥工業大學學報（自然科學版），2012，35 （11）1444－1446．

LI Yuan，LIU Xiao-jun，WANG Jing，et al．Ivestigation of tribological properties of concave surfaces under grease lurication［J］．Journal of Hefei University of Technology，2012，35 （11）1444－1446．

［12］白新瑞，劉曉玲．乏油條件下圓柱滾子軸承的彈流潤滑分析［J］．潤滑與密封，2013，38（1）45－50．

BAI Xin-rui，LIU Xiao-ling．Analysis of elastohydrodynamic lubrication for cylindrical roller bearings under starved oilsupply condition［J］．Lubrication Engineering，2013，38（1）45 －50．

［13］Yang P，Wen S．A generalized reynolds equation for non-Newtonian thermal elastohydrodynamic lubrication［J］．ASME Journal of Tribology，1990，112（10）：631－636．

［14］Mihailidis A，Agouridas K，Panagiontidis K．Non-newtonian starved thermal-elastohydrodynamic lubrication of finite line contact［J］．Tribology Transactions，2013，56：88－100．

［15］Yang Peiran，Wen Shizhu．A generalized Reynolds equation based on non-newtonian flow in lubrication mechanics［J］．Acta Mechanica Sinica，1991，3（5）：283－289．

［16］Venner C H，Ten Napel W E，Bosma R．Advanced multilevel solution of the EHL line contact problem［J］．ASME Journal of Tribology，1993，112（7）：426－431．

Numerical analysis of elastohydrodynamic lubrication for cycloid drives

SUN Zhang-dong，ZHU Cai-chao，LIU Huai-ju，LIU Ming-yong，GU Zong-lin
（The State Key Laboratory of Mechanical Transmission，Chongqing University，Chongqing 400030，China）

Numerical analysis of grease elastohydrodynamic lubrication for cycloid drives was performed．The minimum film thickness of meshing process was obtained．Comparting it with that obtained from an empirical formula，the correctness of model was verified．The average radius of cycloid wheel's lubrication state was taken as the judgment basis，the effects of rotating speed，load and rheological index on grease EHL were analyzed．Results showed that with increase in rotating speed，the second pressure peak drops and moves toward the inlet，the film thickness increases；with increase in load，the contact field width increases，the second pressure peak rises and moves toward the out let，the film thickness decreases slightly；with increase in the rhoological index n，the second pressure peak firstly increases and then decreases，finally approaches Hertz pressure，and it moves toward the outlet，the film thickness decreases．The effect of the short ampliitude coefficient k1on lubrication was discussed，it was shown that under the conditions to meet design requirements；the decrease in k1can improve the lubrication performance．

Cycloid Drive；Grease Lubrication；Numerical Analysis；Ostwald

TH117

A

10．13465／j．cnki．jvs．2014．23．035

由國家自然科學基金（51175523）；教育部博士點基金（20110191110039）資助

2013－07－08 修改稿收到日期：2013－12－12

孫章棟男，博士生

朱才朝男，博士，教授，博士生導師，1967年生