基于Carreau流变模型的渐开线齿轮摩擦因数研究Word格式文档下载.docx

1、齿轮传动中的齿面摩擦力将引起功率损失和表面磨损，从而降低传动效率，加剧轮齿失效，产生振动与噪声，此外还将影响齿根弯曲疲劳强度与齿面接触疲劳强度1等。因此，准确求解啮合齿面上各点的摩擦因数，进行齿面摩擦特性研究对于减少摩擦损失、增大轮齿承载能力和改善系统传动性能等具有重要的意义。目前齿面摩擦因数的计算方法主要分为两类。一类以齿面摩擦特性试验为基础2。Benedict和Kelley3采用圆盘模拟试验机测量了模拟试件的摩擦力。Rebbrchi等4设计出齿面摩擦力动测试验台，并将测试结果与相关研究结论进行了比较。周长江等2指出齿面摩擦因数动测实验的优越性，并补充了线外啮合冲击摩擦模型及其摩擦因数的计算

2、方法。另一类以弹流润滑理论为基础。润滑油的非牛顿流体模型是准确预测弹流润滑膜的油膜厚度和弹流牵曳特性的关键因素5。高速重载场合下，弹流润滑接触区的压力很高，润滑剂此时会具有极强的剪切稀化特性、黏性和弹性，用牛顿流体模型来计算弹流润滑摩擦力往往产生较大的误差。为此，人们先后提出了多种非牛顿流体模型。其中对数函数型Ree-Eyring流变模型在弹流润滑领域得到了广泛应用68。然而近些年的研究表明，相对于对数函数非牛顿流变模型，基于幂函数非牛顿流变模型的理论分析更能真实地预测润滑油膜在弹流工况下的润滑特性910。由于齿轮传动大都处于混合弹流润滑状态，Johnson等11结合粗糙表面的 Greenwo

3、od-Williamson模型，提出了载荷分担理论。Gelinck和 Schipper12将载荷分担理论用于了解决线接触弹流润滑问题。近年来，Akbarzadeh和Khonsari1314应用此理论对直齿轮进行了考虑表面粗糙度的弹流润滑分析。Han等15应用此理论研究了计及表面粗糙度的斜齿轮润滑特性。但以上研究大都基于牛顿流体模型或者Ree-Eyring流变模型。本文作者将结合载荷分担和弹流润滑理论，并基于Carreau流变模型研究润滑剂流变性以及转速、载荷、表面粗糙度、润滑剂的压黏特性等对渐开线齿轮齿面摩擦因数的影响。1 混合弹流润滑模型1.1 载荷分担理论绝大多数齿轮传动在混合弹流状态下工

4、作。Johnson等认为齿面法向载荷F d由润滑油膜和粗糙峰接触共同承担5，齿面摩擦力F f则由润滑油膜摩擦力F fh和粗糙峰接触摩擦力F fa合成，即其中:式中:1/1为润滑油膜承载比例;1/2为粗糙峰承载比例。两表面在相对滑动过程中由于润滑油膜的黏性剪切而产生滑动摩擦力，其大小为润滑膜作用在表面上的剪应力和，即:x其中剪应力和润滑剂的流变模型有关:为切应变率，等于流动速度沿流体厚度方向的变化梯度。通常认为粗糙峰峰点接触处于边界润滑状态，研究表明1617，边界润滑下的滑动摩擦因数在较小的范围内 （0.070.15）变化。若忽略载荷、卷吸速度和滑动速度的影响，在润滑剂特性和齿面粗糙度相同时，一

5、般认为边界润滑下所有的粗糙峰具有相同的摩擦因数f a，因此粗糙峰摩擦力可表示为1.2 流变模型常用的幂函数剪切稀化非牛顿流变模型为Carreau-Yasuda模型，其本构方程1819为为低压剪切黏度。黏度对压力和温度的变化均十分敏感，压力增加会使黏度急剧增大，温度增加却使黏度急剧减小。因此，在齿轮润滑研究中，必须把黏度同时作为压力和温度的函数。在弹流润滑领域被广泛使用的2个黏度模型是Baras指数黏度模型和Roelands黏度模型。而自由体积黏度模型能够精确预测润滑油的黏压关系18，20。本文作者选用Doolittle-Tait自由体积黏度模型来描述弹流工况下润滑膜的黏压关系。黏压温方程为下标

6、r代表参考状态;参数B0、R0、则对应于不同的润滑剂。润滑油体积随着压力的变化关系由Tait状态方程描述:由式（1）（3）知，要获得总摩擦因数，需要知道润滑油膜承载比例或粗糙峰承载比例。1.3 线接触弹流润滑完全光滑 （即不考虑表面粗糙度）的两滚子间中心油膜厚度可以由MOE方程21得到:上式中的量纲一参数分别定义20为如果将式 （12）中的E替换为E/1，把F d替换为F d/1，即可得到适用于等温无限长线接触混合弹流润滑即考虑表面粗糙度的中心油膜厚度公式:而式 （13）包含两个未知数1、。要想求解方程（13），还需要一个方程。1.4 粗糙线接触理论当两粗糙表面间的油膜厚度小于一定值时，表面间

7、将发生粗糙峰直接接触。此时，载荷由润滑油膜和粗糙峰共同承担，即处于混合润滑状态。Greenwood和Tripp22研究了当两粗糙表面接触时，接触中心压力即由粗糙峰承担的压力p a与润滑油的油膜厚度及接触表面的微观形貌参数之间的关系为函数F5/2通常由下式定义:n s为粗糙峰密度;s为粗糙峰半径;s为粗糙峰均方根值;d d为粗糙峰中心线距表面中心线之间的距离。Gelinck和Schipper23则研究了两粗糙表面线接触时，中心接触压力p a和最大Hertz压力p H、综合曲率半径R以及表面粗糙度微观形貌参数之间的关系，得到了很高精确度的拟合公式:式 （17）为干接触时的情况。当两粗糙表面处于混合

8、弹流润滑状态时，根据Johnson的载荷分担理论，将上式中的E/2，F d替换为F d/2，n s替换为n s2，便可得到在混合弹流润滑中由粗糙峰承担的压力为1.5 油膜承载比例求解Greenwood-Tripp接触理论与 Gelinek-Schipper接触理论是从不同的角度来计算线接触时接触区中心粗糙峰的接触压力p a的。但在某一润滑状态下，量纲一载荷参数W、综合曲率半径R、综合弹性模量E、表面粗糙度等相关参数都是定值。因此，Greenwood-Tripp接触理论与 Gelinek-Schipper接触理论两者所算得的粗糙峰压力p a值应该是相等的，即:由于MOE方程的膜厚对应于牛顿流体，

9、因此需要对牛顿流体膜厚进行修正。考虑到润滑油剪切稀化效应等非牛顿流变特性对弹流润滑油膜厚度的影响，Jang和Khonsaril24基于Carreau流变模型给出了线接触弹流中心膜厚修正公式:另外考虑到齿轮啮入端润滑油温度变化对油膜厚度的影响，根据Hsu和Lee25的研究，引入热修正系数C t，W、G分别为量纲一载荷和材料参数;L c为热载荷参数;s r为滑滚比;h ct为非牛顿流体热油膜厚度。k f为润滑剂的热传导系数;为温黏系数。因此首先假定1、2，利用式 （13）得到牛顿流体油膜厚度h Nc，考虑热效应和非牛顿效应的影响，利用式（20）、（21）得到修正后的非牛顿热油膜厚度h ct并将其代

10、入式 （19）。如果方程两边的差大于允许值，重新假定1，直到方程两边的差小于允许值，从而可以得到不同啮合位置时的1、2、h ct。1.6 温度分析在齿轮传动润滑过程中，润滑油膜因受到压缩和黏性剪切作用而发热，导致润滑油温度升高和黏度降低。当高速滚动或滚动中伴有较大的滑动时，齿轮传动润滑计算尤其是摩擦力的计算，必须考虑润滑中的热效应。发热产生的能量通过传导和对流而散失，因此，润滑油膜的温度场需要根据能量方程及界面上的边界条件来确定17。1.6.1 能量方程线接触弹流润滑的总能量方程17为f、c f、u、k f、T分别为润滑油的密度、比热容、速度、热传导系数和温度。准稳态条件下的能量能量方程可化简

11、为上式左边为热传导项，右边分别代表热耗散和对流项，如果将粗糙峰接触产生的热Q a包含进去，需在上式右侧加上一项Q a:u s为相对滑动速度，见本文2.1节。因此考虑粗糙峰接触的能量方程为将式 （5）对油膜厚度求导，可得速度方程:e=反映了剪应力沿膜厚方向的变化，其值等于p/x，由压力流引起。1.6.2 边界条件求解能量方程时，需根据齿轮的热传导性能来确定油膜与齿轮轮齿接触面上的边界条件。界面温度的确定，与齿轮的热传导性能有关。对于弹流接触问题，接触面的温度计算可以归结为具有移动热源的半无限大体的热传导问题，其基本方程26可简化为1.6.3 温度场求解能量方程 （25）和速度方程（26）相互耦合

12、，求解时采用有限差分法。首先划分网格，然后假定均匀温度场和线性速度分布，利用式 （26）得到有效黏度以及剪应变率后，采用有限差分法求解速度方程（26）得到新的速度分布，将新的速度分布代入能量方程 （25），得到新的温度场。并利用复化积分法求解界面温度方程 （27），得到两表面的温度。反复迭代即可得到最终的温度场和剪应力分布。1.7 摩擦因数求解由式 （16）可得总摩擦因数为因此根据本文1.11.5节求解得到粗糙峰承载比例1/2，并根据本文1.6节温度分析，求解得到剪应力，将其代入式 （28），即可得到总摩擦因数f。2 渐开线直齿轮啮合模型齿轮啮合是相当复杂的运动过程，考虑到接触区宽度远小于齿面

13、在接触点的曲率半径，因此对啮合齿面作适当简化，将轮齿沿啮合线上任一点的啮合情况可转换为两个当量圆柱的接触情况17。2.1 齿轮接触参数当渐开线齿轮在啮合线上任意点K啮合时，根据渐开线的性质，两个当量圆柱的中心分别在点N p和N g（见图1），则两齿轮轮齿润滑接触区内啮合点K处的瞬时曲率半径为则综合曲率半径为为啮合点K的坐标 （以啮合起始点A点为原点，x轴正方向沿啮合线方向）。两表面相对于接触点K的速度为所以两个当量圆柱的平均速度即卷吸速度为而两个当量圆柱的相对滑动速度为2.2 载荷分布基于最小弹性势能原理建立载荷分布模型，可计算得到轮齿间的载荷分布。直齿轮双齿啮合区载荷分配率在0.3340.6

14、66之间，单齿啮合时载荷分配率为1。直齿轮的载荷分配与齿廓啮合点的位置有关，但不受模数、压力角等齿轮基本参数的影响。3 分析与讨论以表1中的齿轮传动为例，对比分析研究其在表2给出的2种润滑剂下的油膜厚度、油膜承载比例及摩擦因数，并研究转速、载荷、表面粗糙度以及压黏特性等对渐开线齿轮齿面摩擦因数的影响。采用2种不同润滑剂时，油膜厚度、油膜承载比例、平均油膜温升、摩擦因数沿啮合线的变化对比分别见图25。从图中可看出，不同润滑剂时的油膜厚度、油膜承载比例不同，因此摩擦因数也不同。采用PAO润滑油的齿轮传动，总体油膜厚度厚，且油膜承载比例高，油膜温升小，因此摩擦因数小。另外啮入区油膜温升大，油膜较薄，

15、油膜承载比例较低，因此摩擦因数较大。而退出啮合段，油膜厚度增厚，油膜承载比例也增大，摩擦因数显著减小。而采用PG460润滑油的齿轮传动，总体油膜厚度相对较薄，而且一个啮合周期内油膜承载比例变化不大，整体比较低，且油膜温升较大，因此摩擦因数较大，接近0.1（粗糙峰摩擦因数）。从以上分析可看出，选用不同的润滑剂对齿轮传动的摩擦因数影响较大，因此润滑剂的选择对齿轮传动非常重要。采用PAO润滑油的齿轮传动，载荷、转速、粗糙度、压黏系数等因素对摩擦因数的影响见图69。从图中可看出:（1）载荷、转速、粗糙度和压黏系数对摩擦因数的影响都较大。（2）其中载荷对油膜厚度影响较大，因此不同载荷下的摩擦因数相差也较

16、大。随着转矩的增加，摩擦因数明显增大;而当转矩增加至一定值，随着转矩的增加啮入区的摩擦因数减小，而啮出区的摩擦因数则增大，但转矩对摩擦因数的影响已不是很大。（3）摩擦因数先是随着转速的增加而显著减小，而当转速增加至一定值时摩擦因数则又随之增大。（4）随着表面粗糙度的增大，摩擦因数随之明显增大。而随着黏压系数的增大，摩擦因数也随之明显增大。4 结论（1）不同的润滑剂的剪切稀化特性不同，油膜厚度和摩擦因数也不同，因此润滑剂的选择对于齿轮传动非常重要。（2）而同一润滑剂下，齿轮载荷、转速、表面粗糙度和润滑剂黏压系数对摩擦因数都有较大影响。（3）摩擦因数随着转矩的增大而显著增大，但超过某转矩值时，摩擦

17、因数则随着转矩的增大而缓慢变化。（4）摩擦因数随着转速的增大而减小，但超过某一转速值时，摩擦因数随着转速的增大而增大。（5）随着表面粗糙度的增大，摩擦因数随之明显增大。随着黏压系数的增大，摩擦因数也随之增大。参考文献【1】高创宽.渐开线齿轮传动的混合弹流润滑研究D.太原:太原理工大学机械工程学院，2005:1126.【2】周长江，唐进元，钟志华，等.齿轮传动齿面摩擦因数计算方法的研究J.润滑与密封，2006，31（10）:185191.Zhou Changjiang，Tang Jinyuan，Zhong Zhihua，et al.Studies on the calculation of te

18、eth surface friction coefficients in gear driveJ.Lubrication Engineering，2006，31（10）:【3】Benedict G H，Kelley B W.Instantaneous coefficients of gear tooth frictionJ.ASLE Transactions，1961，4（1）:5970.【4】Rebbrchi B，Oswald F B，Townsend D P.Measurement of gear tooth dynamic frictionC/Proceeding of 7th ASME

19、 International Power Transmission and Gears Conference.San Diego，1996:355363.【5】刘剑平.非牛顿流体弹流润滑膜失效的理论与实践研究D.北京:北京邮电大学自动化学院，2012:1124.【6】Johnson K L，Greenwood J A.Thermal Analysis of an Eyring Fluid in EHD TractionJ.Wear，1980，61:353374.【7】Conry T F，Wang S，Cusano C.A Reynolds-Eyring equation for elastoh

20、ydrodynamic lubrication in line contactsJ.Journal of Tribology，1987，109:648658.【8】Wang Youqiang，Yi Xuejuan.Non-Newtonian transient thermoelastohydrodynamic lubrication analysis of an involute spur gearJ.Lubrication Science，2010，22（10）:465478.【9】Bair S，Qureshi F.The generalized Newtonian fluid model

21、and elastohydrodynamic film thicknessJ.Journal of Tribology，2003，125:7075.【10】Bair S，Vergne P，Querry M.A unified shear-thinning treatment of both film thickness and traction in EHDJ.Tribology Letters，2005，18（2）:145152.【11】Johnson K L，Greenwood JA，Poon SY.A simple theory of asperity contact in elasto

22、hydrodynamic lubricationJ.Wear，1972，19（1）:91108.【12】Gelinck E R M，Schipper D J.Calculation of Stribeck curves for line contactsJ.Tribology International，2000，33（3）:175181.【13】Akbarzadeh S，KhonsariM M.Performance of spur gears considering surface roughness and shear thinning lubricantJ.Journal of tri

23、bology，2008，130（2）:021503.【14】Akbarzadeh S，Khonsari M M.Thermoelastohydrodynamic analysis of spur gears with consideration of surface roughnessJ.Tribology Letters，2008，32（2）:129141.【15】Han L，Zhang DW，Wang F J.Predicting film parameter and friction coefficient for helical gears considering surface ro

24、ughness and load variationJ.Tribology Transactions，2013，56（1）:4957.【16】Williams JA.Advances in themodelling of boundary lubricationJ.Tribology Series，2002，40:3748.【17】温诗铸，黄平.摩擦学原理M.北京:清华大学出版社，2002.【18】Bair SS.High pressure rheology for quantitative elastohydrodynamicsM.Amsterdam:Elsevier，2007.【19】Ya

25、suda K Y，Armstrong R C，Cohen R E.Shear flow properties of concentrated solutions of linear and star branched polystyrenesJ.Rheologica Acta，1981，20（2）:163178.【20】Doolittle A K.Studies in Newtonian flow.II.the dependence of the viscosity of liquids on free-spaceJ.Journal of Applied Physics，2004，22（12）

26、:14711475.【21】Moes H.Optimum similarity analysiswith applications to elastohydrodynamic lubricationJ.Wear，1992，159（1）:5766.【22】Greenwood J A，Tripp JH.The contact of two nominally flat rough surfacesJ.Proceedings of the institution ofmechanical engineers，1970，185（1）:625633.【23】Gelinck E R M，Schipper

27、D J.Deformation of rough line contactsJ.Journal of Tribology，1999，121（3）:449454.【24】Jang JY，Khonsaril M M.Correction factor formula to predict the central and minimum film thickness for shear-thinning fluids in EHLJ.Journal of Tribology，2008，130:024501.14.【25】Hsu C H，Lee R T.An efficient algorithm f

28、or thermal elastohydrodynamic lubrication under rolling/sliding line contactsJ.Journal of Tribology，1994，116（4）:762769.【26】Jaeger J C.Moving surfaces of heat and the temperature at sliding contactsJ.Proc R Soc N.S.W.，1942，56:203224.DOI:10.3969/j.issn.02540150.2015.03.009*基金项目:国家高技术研究发展计划 （863计划）项目（2012AA110701）.【文献来源】