刘胜超，王东峰，喻炜

(1.洛阳轴研科技股份有限公司 精密轴承制造部，河南 洛阳 471039；2.北京市三一重机有限公司，北京 102206)

符号说明

a—— 接触椭圆的长半轴

b—— 接触椭圆的短半轴

Dpw—— 球组节圆直径，mm

Dw—— 球径，mm

fi—— 内沟曲率系数

fe—— 外沟曲率系数

F(ρ)——主曲率差函数

Fa—— 轴向载荷，N

k—— 接触椭圆的椭圆率

K(e)——第一类椭圆积分

L(e)——第二类椭圆积分

Z—— 球数

ρ1Ⅰ——球在轴向平面的主曲率

ρ1Ⅱ——球在径向平面的主曲率

ρ2Ⅰ——套圈在轴向平面的主曲率

ρ2Ⅱ——套圈在径向平面的主曲率

α——初始接触角，(°)

φ——椭圆积分变量

γ——无量纲参数，γ=(Dw/Dpw)cosα

Hertz接触参数求解过程中，K(e)，L(e)和k是由循环迭代的方法求得的，迭代次数非常多，每次迭代都需要求解K(e)和L(e)，计算量大。这里利用Matlab曲线拟合工具箱，对K(e)，L(e)，k分别与F(ρ)进行拟合，得到它们和主曲率差函数F(ρ)的表达式，用拟合的结果求解给定型号轴承的其他Hertz接触参数，并与Hertz接触理论值以及最小二乘法线性回归法(回归法)得到的简化式[1]求得的结果进行对比。

1 Hertz接触参数

下面主要介绍与椭圆积分以及椭圆率计算相关的一些Hertz参数，其他参数不再赘述，可见文献[2-3]。

1.1 接触点的主曲率

在接触点，球与沟道表面的主曲率描述了球和沟道表面在接触点的几何特征，其影响到接触应力与变形及轴承动态性能[4]。球和沟道表面的主曲率计算方法见表1。

表1 球轴承接触点主曲率计算公式

主曲率差及主曲率和分别为

(1)

∑ρ=ρ1Ⅰ+ρ1Ⅱ+ρ2Ⅰ+ρ2Ⅱ。

(2)

1.2 椭圆积分

第一类完全椭圆积分为

(3)

第二类完全椭圆积分为

(4)

椭圆率为

k=b/a。

(5)

F(ρ)可以由K(e)，L(e)和k表示为[5]

(6)

2 拟合法及回归法

Matlab中提供的曲线拟合工具箱(Curve Fit Toolbox)可以对工作空间中的2组参数的数据点进行拟合，根据不同的拟合规则得到2个参数的表达式。……