大模数磨损齿轮齿条静力学分析

代绍雄，贡宗星，陈保家，方俊豪

大模数磨损齿轮齿条静力学分析

代绍雄，贡宗星，陈保家*，方俊豪

（三峡大学 机械与动力学院，湖北 宜昌 443002）

针对磨损故障形式下的齿轮齿条起升机构，通过分析计算相同载荷下1/12磨损、1/6磨损、1/4磨损下的齿轮齿条基本参数变化情况，利用将齿轮齿条磨损问题量化为ISO相关基本参数改变的方法，对正常和不同磨损故障形势下的弯曲强度进行理论计算，得到其弯曲应力理论计算值。再通过有限元方法对正常和1/12磨损、1/6磨损、1/4磨损下的大模数齿轮齿条起升机构进行静力学分析，研究不同磨损状况下齿轮齿条的弯曲应力和接触应力应变分布规律，并将静力学分析得出的弯曲应力分析结果与利用ISO公式计算的弯曲强度计算结果进行对比分析，验证了将齿轮齿条磨损问题转化为相关基本参数改变问题的可行性，并通过对大模数故障齿轮齿条起升机构的静力学分析，得到其故障形式下的静力学特性，为大模数齿轮齿条起升机构的故障诊断和运行健康评估提供必要的理论基础。

齿轮齿条；建模；静力学分析；磨损齿轮

随着制造业的不断发展，大模数齿轮齿条起升机构由于其起升重量大、提升高度大、体积结构大、起升平稳性好、制动平稳、操作简单的特点，广泛应用于自升式海工升降平台、升船机驱动系统等大型机械设备中[1]。大模数齿轮齿条式起升机构作为核心驱动装置，制造工艺复杂、产品更换困难，在其不断啮合的同时，齿轮会出现不同程度的磨损，长期反复作用下可能会导致其传动失效。所以对磨损齿轮齿条起升机构进行理论强度计算校验和静力学分析研究就显得尤为重要。

分析大模数磨损齿轮齿条静力学特性是常见的研究方法。邵博睿等[2]研究改变压力角和齿条宽度对齿轮齿条机构接触力的影响。张俊等[3]利用有限元软件对无磨损、轻微磨损、严重磨损三种状况下的齿轮进行固有特性分析，磨损使得齿轮固有频率增大，对振型影响不大。张越[4]以标准渐开线直齿圆柱齿轮为对象，研究考虑齿面闪温的粘着磨损数值计算方法，分析齿形参数和工作参数对齿面磨损量的影响，结果表明，增大模数、压力角、传动比、齿宽或者减小输入扭矩均可使得齿面磨损量减小。冯松等[5]计算分析存在均匀磨损、微点蚀和宏观点蚀等齿面磨损时的齿轮啮合刚度。李学志等[6]针对斜齿轮啮合更替过程受到载荷突变影响，引起齿面疲劳点蚀和冲击的问题，对斜齿轮静力学接触和修形进行分析研究。邱林等[7]对渐开线弧齿圆柱齿轮的振动特性进行研究，从而分析获得该齿轮在低阶模态下的主振型和固有频率。王兆山等[8]进行了齿轮齿条磨损量数据检测方案的设计以及应用多项式拟合磨损量数据分布曲线，进而探讨磨损规律。

本文将磨损齿轮的问题转换成齿轮基本参数变化，计算与IOS公式相关变量的值，利用ISO公式[9]计算磨损齿轮弯曲应力的理论值，与静力学分析结果对比，研究其静力学特性，验证方法的可行性。

1 故障齿轮齿条起升机构装配

为了研究正常齿轮齿条起升机构和故障齿轮齿条的静力学特性，利用SolidWorks软件对齿轮齿条进行三维建模和装配，齿轮齿条的材料及设计参数如表1所示。

表1 齿轮齿条的材料及设计参数

对正常齿轮齿条进行三维建模后对齿轮的一个齿形进行均匀磨损处理，如图1所示。

对齿轮进行1/12、1/6、1/4均匀磨损处理，然后与齿条进行装配，得到装配体如图2所示。

图1 齿轮均匀磨损模型

图2 磨损齿轮齿条装配体简化模型

2 磨损齿轮齿条弯曲应力计算

对齿轮齿条弯曲强度的计算主要运用ISO公式，对于磨损齿轮来说，一些参数的变化会导致齿轮齿条弯曲应力的变化。本节根据磨损齿轮基本参数的变化，应用ISO公式对磨损齿轮弯曲应力进行计算。齿轮弯曲应力计算的相关变量符号说明如表2所示。

表2 齿根应力计算相关变量符号说明

引入[10]，在额定外载荷F＝6.664E4 N作用下齿轮弯曲应力计算公式为：

式中：当齿轮未发生磨损时，引入Y，查图[10]

可取1.53；Y取值0.789；Y与尺制、变位系数和齿数有关，查图可取2.96。

当齿轮发生均匀磨损时，把齿轮当做变位齿轮来计算齿轮的弯曲应力，、都会发生变化，从而导致、ε、Y、Y发生改变，求出磨损齿轮的相关变量的值，从而求出磨损齿轮的弯曲应力。相关参数的计算如下[7]：

相关参数计算结果如表3所示。

表3 齿根应力计算相关参数计算值

由表3可得1/12、1/6、1/4磨损齿轮的在相同载荷的作用下的弯曲应力相比较于相同载荷下的正常齿轮的弯曲应力192.083 MPa分别增加了3.8%、11.05%、16.5%。说明齿轮磨损会导致齿轮弯曲应力的增大。

3 磨损齿轮齿条静力学分析

3.1 建立分析模型，施加边界条件及载荷

将齿轮齿条转装配体文件导入软件中对其进行网格划分，图3为部分齿轮齿条模型的网格划分图，采用六面体的网格单元，对其接触部位的网格进行局部细化，得到的总划分单元为299688个，总划分节点为185969个。如图4所示，对齿条底部施加固定约束，对齿轮轴孔的自由度进行约束，为了便于添加载荷，轴孔添加圆柱型约束，对齿轮轴径施加扭矩1。

3.2 计算结果及分析

对模型进行边界条件和载荷的设置以后，对其进行静力学求解和后处理，得到正常齿轮齿条和1/12、1/6、1/4磨损齿轮弯曲应力及接触应力静力学分析结果如图5、图6所示。

由图5可知正常齿轮和不同磨损程度的齿轮通过静力学分析得到的齿轮弯曲应力，且最大弯曲应力的位置不断向齿端靠近。与理论计算值的对比如表4所示。

由表4可知，正常齿轮和1/12、1/4、1/6磨损齿轮理论计算的弯曲应力与静力学分析得出的值相差都在10%左右，说明根据磨损齿轮基本参数的变化利用ISO公式计算磨损齿轮的弯曲应力的方法具有一定的可行性。

表4 齿根弯曲应力计算与静力学分析对比表

图3 模型的网格划分

图4 齿轮齿条有限元分析模型

图5 弯曲应力静力学分析结果

图6 接触应力静力学分析结果

由图6可知，正常齿轮和不同磨损程度的齿轮通过静力学分析得到的齿轮接触应力为248.18 MPa、255.6 MPa、260.77 MPa、267.56 MPa，相比于正常齿轮，1/12、1/4、1/6磨损齿轮的接触应力分别增加了2.9%、5.07%、7.8%。

从图5、图6可以看出，齿轮轮齿的接触应力一般大于弯曲应力，弯曲应力一般分布在齿轮齿根处，接触应力一般分布在轮齿的表面，可以通过对齿面进行强化处理和对齿轮轮齿圆弧化处理来增加齿轮的接触强度和弯曲强度，为大模数齿轮齿条起升机构的故障诊断和运行健康评估提供必要的理论基础。

4 结论

（1）结合ISO公式，把磨损齿轮的磨损量转换成基本参数的变化，推导出ISO公式中的变量值的变化，从而计算出磨损齿轮的齿根弯曲应力。

（2）利用静力学分析软件对正常齿轮齿条和不同程度的磨损齿轮齿条进行静力学分析，得出其静力学特性，与理论计算值进行对比分析，验证了将齿轮齿条磨损问题量化为相关基本参数改变的方法的正确性，为以后的磨损齿轮弯曲应力计算和静力学分析提供理论基础。

（3）通过对正常和磨损齿轮齿条的理论弯曲应力计算和静力学分析，得出了磨损齿轮的静力学特性，为大模数齿轮齿条起升机构的故障诊断和运行健康评估提供必要的理论基础。

[1]向正新. 大模数齿轮齿条啮合仿真研究[D]. 荆州：长江大学，2018：1-16.

[2]邵博睿，包振明，谭祚炜，等. 自升式平台齿轮齿条磨损的有限元分析[J]. 机械工程与自动化，2019（5）：93-94，96.

[3]张俊. 直齿圆柱齿轮磨损建模及其动态响应分析[D]. 赣州：江西理工大学，2019.

[4]张越. 考虑齿面温升的渐开线直齿轮粘着磨损计算方法[D]. 长沙：湖南大学，2018.

[5]冯松，毛军红，谢友柏. 齿面磨损对齿轮啮合刚度影响的计算与分析[J]. 机械工程学报，2015，51（15）：27-32.

[6]李学志，陈蒙，陈绍康，等. 变位斜齿轮静力学分析和修形研究[J]. 机械，2019，46（1）：12-15，20.

[7]邱林，宋爱平，彭云，等. 渐开线弧齿圆柱齿轮有限元模态分析[J]. 机械，2016，43（10）：6-10.

[8]王兆山. 在役自升式平台起升齿轮齿条系统弯曲强度及疲劳损伤度研究[D]. 天津：天津大学，2018.

[9]王建敏. 大模数渐开线直齿圆柱齿轮弯曲强度研究[D]. 北京：机械科学研究总院，2006：9-17.

[10]GB/T 348－1997，渐开线圆柱齿轮承载能力计算方法[S]. 北京：中国标准出版社，1997.

Static Analysis of Gear Rack with Large Modulus Wear

DAI Shaoxiong，GONG Zongxing，CHEN Baojia，FANG Junhao

( College of Mechanical and Power Engineering, China Three Gorges University, Yichang 443002, China )

In this paper, focusing on the rack and pinion lifting mechanism under the wear failure mode, through the analysis and calculation of the basic parameters of the gear and rack under the same load of 1/12 wear, 1/6 wear, 1/4 wear, by quantifying the gear rack wear problem as the change of ISO related basic parameters, the bending strength under normal and different wear failure situations is calculated theoretically, and the theoretical calculation value of bending stress is obtained. Then, the static analysis of the large module rack and pinion lifting mechanism under normal and 1/12 wear, 1/6 wear and 1/4 wear is carried out by finite element method, and the distribution of bending stress and contact stress of the gear rack under different wear conditions is studied. The bending stress analysis results from static analysis are compared with the bending strength calculation results calculated by ISO formula, which verifies the feasibility of transforming the gear rack wear problem into the related basic parameter change problem. Through the static analysis of the large module fault rack and pinion lifting mechanism, the static characteristics under the fault mode are obtained, which provides the necessary theoretical basis for the fault diagnosis and operation health assessment of large module gear and rack lifting mechanism.

rack and pinion；modeling；static analysis；worn gear

TH132.412

A

10.3969/j.issn.1006-0316.2021.03.004

1006－0316 (2021) 03－0023－05

2020－09－01

国家自然基金面上项目：多系统耦合作用下的大模数齿轮齿条式起升机构故障预示及健康评估研究（51975324）

代绍雄（1996－），男，湖北洪湖人，硕士研究生，主要研究方向为大模数齿轮齿条动力学分析。*通讯作者：陈保家（1977－），男，湖北孝感人，博士，教授、博士生导师，主要研究方向为机械装备状态监测、故障诊断及可靠性评估与寿命预测、大模数齿轮齿条动力学分析，E-mail：1550603416@qq.com。