• CN:11-2187/TH
  • ISSN:0577-6686

机械工程学报 ›› 2018, Vol. 54 ›› Issue (23): 56-62.doi: 10.3901/JME.2018.23.056

• 机构学及机器人 • 上一篇    下一篇

基于解析有限元的齿根裂纹时变啮合刚度计算方法

吴家腾, 杨宇, 程军圣   

  1. 湖南大学汽车车身先进设计制造国家重点实验室 长沙 410082
  • 收稿日期:2017-10-27 修回日期:2018-06-01 出版日期:2018-12-05 发布日期:2018-12-05
  • 通讯作者: 杨宇(通信作者),女,1971年生,博士,教授。主要研究方向为机械设备状态监测及故障诊断。E-mail:yangyu@hnu.edu.cn
  • 作者简介:吴家腾,男,1990年生,博士研究生。主要研究方向为设备状态监测及故障诊断。E-mail:Jiat_wu@126.com;程军圣,男,1968年生,博士,教授。主要研究方向为机械设备状态监测及故障诊断。E-mail:signalp@tom.com
  • 基金资助:
    国家重点研发计划(2016YFF0203400)、国家自然科学基金(51575168,51375152)、智能型新能源汽车国家2011协同创新中心和湖南省绿色汽车2011协同创新中心资助项目。

The Time-varying Mesh Stiffness Calculation for Gear Tooth Crack based on Analytical-finite Element Method

WU Jiateng, YANG Yu, CHENG Junsheng   

  1. State key Laboratory of Advanced Design and Manufacture for Vehicle Body, Hunan University, Changsha 410082
  • Received:2017-10-27 Revised:2018-06-01 Online:2018-12-05 Published:2018-12-05

摘要: 当齿轮发生故障时,时变啮合刚度的变化能够反映齿轮故障特征大小。因此,时变啮合刚度在齿轮传动过程中是一个重要的动力学参数。提出一种新的齿根裂纹啮合刚度计算方法,即解析有限元法(Analytical-finite element method,A-FM)。考虑到齿轮发生故障时,啮合刚度解析模型计算精度较低,将应力强度因子引入裂纹齿轮的啮合刚度计算过程。首先定义应力强度因子与啮合刚度之间的关系,通过建立齿轮接触模型计算裂纹尖端附近的应力强度因子,然后将计算结果替代解析模型中故障刚度部分。由于应力强度因子能够敏感地识别齿根裂纹的局部微小变化,故该方法相比于解析法具有更高的计算精度,相比于有限元法具备更快的计算效率。同时,建立6自由度动力学模型,通过对其振动响应进行分析,仿真结果验证了所提方法的可行性。

关键词: 动力学建模, 解析有限元法, 时变啮合刚度, 应力强度因子

Abstract: The time-varying meshing stiffness variation can reflect the size of gear fault feature when the gear fault occurs. Therefore, the time-varying meshing stiffness is an important dynamic parameter in the process of gear transmission. A novel tooth crack meshing stiffness calculation method namely analytical-finite element method (A-FM) is presented. Considering the calculation accuracy of the analytical model of the meshing stiffness is low when the gear fault occurs, the stress intensity factor is introduced into the calculation of the meshing stiffness of the cracked gear in this paper. First of all, define the relationship between the stress intensity factors and the meshing stiffness, calculate the stress intensity factor near the crack tip by means of establish the gear contact model, then replace the part of fault stiffness in the analytical model with calculation result. This method is more accurate compare with analytical method and more effective compare with finite element method due to the stress intensity factor is sensitive to small local changes in the tooth root crack. Simultaneously, the six-degree of freedom dynamic model is established. Through the analysis of vibration response, the simulation results verify the feasibility of the proposed method.

Key words: analytical-finite element method, dynamical modeling, stress intensity factor, time-varying mesh stiffness

中图分类号: