[发明专利]行星轮系非赫兹接触的接触刚度计算与动力学建模方法

说明书

本发明公开了一种行星轮系非赫兹接触的接触刚度计算与动力学建模方法。所述接触刚度计算方法的步骤如下：获取轮齿表面测试数据；采用渐开线加权补偿方法将轮齿表面测试数据与轮齿表面弧度相匹配计算接触刚度。所述动力学建模方法按照上述方法计算接触刚度；然后将轴向抗压刚度、弯曲刚度、剪切刚度、计算的接触刚度等带入行星轮系动力学模型，采用ode45计算方法对动力学方程进行求解，获得轮齿的振动数据；对轮齿的振动数据进行分析获得轮齿的受力。将轮齿表面的接触形式考虑为非赫兹接触，充分考虑轮齿表面的微观形貌，减小接触刚度的计算误差，提高动力学建模的准确性，提高接触特性分析、振动特性预估、功率损耗计算、寿命预测的准确性。

技术领域

本发明涉及仿真分析技术领域，具体涉及一种行星轮系非赫兹接触的接触刚度计算与动力学建模方法。

背景技术

行星轮系传动比大，承载能力强，传动形式多样，不同传动形式的传动比不相同，例如在相同的轮齿参数与太阳轮输入转速条件下，传动形式不同，行星架输出与齿圈输出的传动比也会不相同。行星轮系传动形式多样，故此在高速重载或者多级变速的机构中，行星轮系得到广泛的应用。对行星轮系进行动力学建模，可以对其运动特性进行分析，计算功率损耗等，可用于预测行星轮系的寿命。

目前的行星轮系的动力学建模考虑轮齿表面的微观形貌不充分，且将轮齿表面的接触形式考虑为赫兹接触，造成计算误差较大，导致进一步的接触特性分析、振动特性预估、功率损耗计算乃至寿命预测的准确性降低。

发明内容

本发明的目的在于克服现有技术中所存在的上述接触形式考虑为赫兹接触对轮齿表面的微观形貌考虑不充分的不足，提供一种行星轮系非赫兹接触的接触刚度计算与动力学建模方法，通过充分考虑轮齿表面的微观形貌，和将轮齿表面的接触形式考虑为非赫兹接触，降低了计算误差，提高后续计算仿真的准确性。

为了实现上述发明目的，本发明提供了以下技术方案：

第一方面，本发明实施例提供了一种行星轮系非赫兹接触的接触刚度计算方法，包括以下步骤：

S1，获取轮齿表面测试数据；采用渐开线加权补偿方法将轮齿表面测试数据与轮齿表面弧度相匹配，计算接触刚度。

优选地，所述步骤S1具体包括以下步骤：

S11，通过超景深显微扫描仪对轮齿表面进行显微扫描，获取轮齿表面测试数据；

S12，对轮齿表面测试数据进行补偿并计算行星齿轮齿面轮廓分形维数；

采用渐开线加权处理方法，结合轮齿表面不同接触位置的压力角和显微扫描仪摄影方向与被扫描轮齿中心之间的夹角，获得不同接触位置的补偿比例；对轮齿表面测试数据进行补偿，补偿后的轮齿表面测试数据记作第一轮齿表面测试数据；

对功率谱函数取双对数后计算最小二乘中线斜率；所述分形维数根据最小二乘中线斜率得到；

S13，轮齿表面测试数据选点，并对选点进行曲线拟合得到轮廓层数据拟合曲线：

轮齿表面测试数据的最大值对应的切面称为表面层，表面层的轮廓线记作特征啮合线；在特征啮合线上，以特征啮合线的数据峰值为中心，左右等间隔取点，进行轮齿表面测试数据选点；对选点进行曲线拟合得到轮廓层数据拟合曲线；

S14，通过轮廓层数据拟合曲线计算实际接触面积与拟合曲线接触面积；

S15，计算表面接触系数，根据变形阶段类型和分形维数的取值范围对无量纲法向载荷进行分段计算；

S16，计算非赫兹接触刚度。

优选地，所述步骤S12根据如下所示的公式，计算分形维数D：

其中，k为对功率谱函数取双对数后计算得到的最小二乘中线斜率；