基于多轴疲劳准则的齿轮点蚀寿命预测  

Gear Pitting Life Prediction Based on Multi-Axial Fatigue Criterion

作　　者：王雄 董庆兵[1,2] 史修江 魏静 秦毅[1,2] 冯成程[1,2] 罗振涛 WANG Xiong;DONG Qingbing;SHI Xiujiang;WEI Jing;QIN Yi;FENG Chengcheng;LUO Zhentao(State Key Laboratory of Mechanical Transmissions,Chongqing University,Chongqing 400044,China;College of Mechanical and Vehicle Engineering,Chongqing University,Chongqing 400044,China;College of Power and Energy Engineering,Harbin Engineering University,Heilongjiang Harbin 150001,China)

机构地区：[1]重庆大学机械传动国家重点实验室,重庆400044 [2]重庆大学机械与运载工程学院,重庆400044 [3]哈尔滨工程大学动力与能源工程学院,黑龙江哈尔滨150001

出　　处：《摩擦学学报》2023年第1期92-103,共12页Tribology

基　　金：重点研发计划(2020YFB2010101);自然科学基金(51905051);两机重大专项(J2019-Ⅳ-0018-0086);重庆大学机械传动重点实验室项目(SKLMT-ZZKT-2021M06);中央高校基本科研业务费专项(2020CDCGJX058,2021CDJKYJH018)资助。

摘　　要：综合齿轮动力学和弹性流体动力润滑理论,建立基于临界平面法多轴疲劳寿命预测模型.首先根据齿轮啮合特性获取齿面接触的时变参数,采用平均滤波方法模拟齿面磨合后的粗糙状态,并将齿面粗糙形貌带入油膜厚度计算,并基于量纲化差分方法建立齿轮的热弹流动力润滑模型;随后,通过润滑界面压力和摩擦力的分布计算近表面应力状态,确定接触近表面任意平面的应力与应变幅;最终,采用临界平面方法计算Smith-Watson-Topper(SWT)参数和最易萌生裂纹的平面,最终确定齿轮疲劳点蚀寿命,并试验验证模型有效性.结果表明:粗糙表面造成压力、油膜厚度和温度等波动较大,最大应力集中分布在表面,疲劳点蚀的微裂纹首先在表面萌生;齿轮疲劳点蚀数值模型可有效预测不同润滑条件下的疲劳点蚀寿命.Pitting is one of the primary failure modes of gears resulting from crack initiation and propagation caused by the friction contact at the power transmission interface.The corresponding fatigue is essentially induced by multi-axial loadings,and the fatigue life is somehow related to the stress/strain states at a specific plane during the meshing process.The uniaxial fatigue methods are thus not applicable any longer to predict the lifetime.In this paper,a multi-axial fatigue life prediction model was established based on the critical plane method coupling with the gear meshing dynamics and elastohydrodynamic lubrication theory.The contact parameters at each time step in a meshing cycle,such as geometrical curvatures,sliding/rolling velocities and loadings,were determined based on the analysis of gear dynamic behaviors and then introduced into the lubrication analysis.A two-dimensional moving average filter was used to take the running-in effects on the surface roughness into account,which had proven to be capable of evaluating the morphology change during the running-in period while remaining the dominant characteristics of the initial morphology.Normally,the roughness peak height was decreased by 20%~50%by such a method,which had been previously proven to be reliable for research purposes on rough surfaces after the running-in process.The sectional morphology was introduced into the calculation of the lubrication film thickness,and a transient thermal elastohydrodynamic lubrication model in line contact was then developed by solving the governing equation system based on the non-normalized discretization method to calculate the surface pressure and tractions.The fast Fourier transform was employed to improve the calculation efficiency.The moving coordinate method based on the analysis of the motion distance in each time interval was used to determine the stress/strain at plane strain conditions in a fixed computational domain,and the near-surface elastic fields in a random plane were properly obtained based on t

关 键 词：齿轮啮合 疲劳点蚀 平均滤波法 热弹流润滑 临界平面法 

分 类 号：TH132.413[机械工程—机械制造及自动化]