一种球轴承接触刚度计算方法技术

本发明专利技术涉及一种球轴承接触刚度计算方法，首先根据球轴承的内部尺寸参数和球轴承载荷，解析确定球轴承内部接触区域参数；然后求解轴承内部载荷分布，确定内滚道和外滚道的接触角、轴承整体轴向、径向位置偏移、以及内、外滚道的接触变形；最后根据轴承的外部载荷和轴承整体的轴向、径向位置偏移，最终确定球轴承接触刚度。能够考虑球轴承不同的转速、载荷工作条件，精确确定球轴承接触刚度。

【技术实现步骤摘要】
一种球轴承接触刚度计算方法
本专利技术球轴承优化设计
，具体涉及球轴承接触刚度计算方法。
技术介绍
轴承是数控装备进给轴和工业机器人关节的主要支撑部件，在外部载荷作用下，轴承内部界面的接触刚度是影响数控装备和工业机器动力学特性的主要因素之一。在理论上关于轴承分析时，一般是假定轴承内外圈是不发生弹性变形的刚体，弹性变形只发生在滚动体上，滚动体与内外圈的接触只会产生局部的接触变形而不会改变内外圈的整体形状和尺寸。这种处理方法虽然为内外圈的位移和轴承弹性变形分析带来了很大的方便，但是它没有反映出轴承内外圈的变形因素，从而无法准确预测轴承的接触刚度。现在也有采用有限元离线分析的方法，通过确定球轴承的刚度矩阵以及接触角度来确定接触刚度，但此种方法难以实现在线监测球轴承在不同工作条件下的接触刚度实时变化情况，无法实际工业应用。因而，传统的球轴承接触刚度计算方法更多的是具有理想模型与经验意义，其精度上不够准确、实时应用效率上也不够高效。
技术实现思路
有鉴于此，本专利技术拟公开一种球轴承接触刚度计算方法，能够本文档来自技高网...

【技术保护点】
1.一种球轴承接触刚度计算方法，其特征在于包括如下步骤：/n1)首先根据球轴承的内部尺寸参数和球轴承载荷，解析确定球轴承内部接触区域参数；/n2)根据球轴承的内部接触区域参数，求解轴承内部载荷分布，确定内滚道和外滚道的接触角、轴承整体轴向、径向位置偏移、以及内、外滚道的接触变形；/n3)根据轴承的外部载荷和轴承整体的轴向、径向位置偏移，最终确定球轴承接触刚度。/n

【技术特征摘要】
1.一种球轴承接触刚度计算方法，其特征在于包括如下步骤：
1)首先根据球轴承的内部尺寸参数和球轴承载荷，解析确定球轴承内部接触区域参数；
2)根据球轴承的内部接触区域参数，求解轴承内部载荷分布，确定内滚道和外滚道的接触角、轴承整体轴向、径向位置偏移、以及内、外滚道的接触变形；
3)根据轴承的外部载荷和轴承整体的轴向、径向位置偏移，最终确定球轴承接触刚度。


2.根据权利要求1所述的球轴承接触刚度计算方法，其特征在于步骤1)中所述球轴承的内部尺寸参数包括：内滚道半径ri，外滚道半径ro、内滚道的曲率半径ri’、外滚道的曲率半径ro’、内滚道的曲率中心、外滚道的曲率中心、轴承的节径rm、滚动体的半径rb；内滚道曲率中心和外滚道曲率中心之间的距离BD＝r′o+r′i-2rb＝(fi+fo-1)2rb、球轴承的径向间隙Pd＝2ro-2ri-4rb、球轴承的初始接触角和各滚动体角度位置ψj＝2π(j-1)/Z；其中，fi＝ri’/2rb,fo＝ro’/2rb，j表示各滚动体的编号，j＝1,2,…,Z，Z表示该轴承的滚动体个数。


3.根据权利要求1所述的球轴承接触刚度计算方法，其特征在于步骤1)中首先建立无外部载荷情况下接触区域参数计算模型，设定滚动体和内、外滚道之间的弹性接触是两个具有不同曲率的椭球体的接触，各椭球体各自经过接触区域的两个主平面各自确定一个点，两个椭球体在所述的点上实现单点接触，也即滚动体和内滚道之间、以及滚动体和外滚道之间均为单点接触；由滚动体和内滚道接触的接触曲率半径关系、滚动体和外滚道接触的接触曲率半径关系联合确定两个椭球体接触时的接触曲率半径。


4.根据权利要求3所述的球轴承接触刚度计算方法，其特征在于步骤1)在无外部载荷情况下确定两个椭球体接触时的接触曲率半径后，建立外部载荷情况下的接触区域参数计算模型：在外部载荷的作用下，滚动体和内、外滚道之间分别形成椭圆接触区域，椭圆的短轴位于滚动体的运动方向；滚动体与内、外滚道椭圆接触区域的长、短轴和内外滚道接触变形的关系如下组合式：












组合式中，rI1,rI2表示其中一个椭球体在经过接触区域的两个主平面上的接触曲率半径；rII1和rII12表示另一个椭球体在经过接触区域的两个主平面上的接触曲率半径；α，b分别表示椭圆接触区域椭圆的长、短轴，δ表示内、外滚道的法向接触变形，Q表示滚动体所受的法向载荷，E1，E2分别表示两个接触椭球体的弹性模量，v1，v2分别表示两个接触椭球体的泊松比，E1＝E2，v1＝v2；根据滚动体和内、外滚道的尺寸参数和上述组合式可以确定A、B的数值；
结合接触区域参数参量α*，b*，δ*与接触椭球体的曲率半径关系设定方程组如下：



上述方程组中，e表示接触区域椭圆的偏心率，即e2＝1-(b/a)2,F(e,π/2)和E(e,π/2)分别表示第一、二类完全椭圆积分：



通过求解上述方程组前两式确定接触区域椭圆的偏心率e；并最终得到α*，b*和α，b的取值。


5.根据权利要求4所述的球轴承接触刚度计算方法，其特征在于步骤1)中结合滚动体和内外滚道关系确定A>B，根据A/B的不同取值，确定e，α*，b*，δ*的变化趋势。


6.根据权利要求1所述的球轴承接触刚度计算方法，其特征在于步骤2)中根据步骤1)确定的内部接触区域参数，由内、外滚道变形...

【专利技术属性】
技术研发人员：赵青，金超，
申请(专利权)人：熵智科技深圳有限公司，
类型：发明
国别省市：广东;44