CN114692323B - 一种无保持架球轴承滚动体打滑的数值分析方法 - Google Patents

Abstract

本发明公开了一种无保持架球轴承滚动体打滑的数值分析方法，以无保持架球轴承为对象，通过建立无保持架球轴承滚动体的动力学模型并求解，在不同轴承载荷、不同轴承转速，对无保持架球轴承滚动体出现打滑进行分析，进而判定滚动体的打滑现象。解决了目前尚无法通过试验手段获取准确的滚动体滑移速度，进而无法判定滚动体的打滑现象的问题。

Description

一种无保持架球轴承滚动体打滑的数值分析方法

技术领域

本发明涉及无保持架球轴承技术领域，尤其是一种无保持架球轴承滚动体打滑的数值分析方法。

背景技术

无保持架球轴承是磁悬浮轴承系统中重要的零件之一，是磁悬浮轴承中的辅助轴承，起着保护主轴的作用；无保持架球轴承内圈与主轴连接，无保持架球轴承外圈与轴承座连接，内外圈之间有滚动体，滚动体在无保持架球轴承不工作时主要堆积在轴承的下半圈，滚动体之间相互接触；在紧急情况下(磁悬浮轴承断电、故障)，主轴会下落到无保持架球轴承的内圈并带动其旋转，无保持架球轴承内外圈之间的滚动体会平衡主轴下落的冲击力和主轴旋转产生的力矩，不会使主轴直接掉落在轴承座上，对主轴起到保护作用；无保持架球轴承在工作时滚动体受力时刻发生变化，非常容易出现打滑现象；其中无保持架球轴承出现打滑是指滚动体和轴承的内圈或外圈在圆周方向上发生相对位移；由于不存在保持架分离滚动体，无保持架球轴承的滚动体在打滑时不仅会导致滚动体与内外圈的润滑油膜失效，使滚动体与内外圈直接接触摩擦，也会导致相邻滚动体发生非连续碰撞与高速摩擦，这都会加剧轴承发热和滚动体磨损进而影响轴承的使用寿命；因此，对无保持架球轴承滚动体的打滑分析有着重要意义。目前尚无法通过试验手段获取准确的滚动体滑移速度，进而无法判定滚动体的打滑现象。

基于上述问题，提出一种无保持架球轴承滚动体打滑的数值分析方法。

发明内容

本发明要解决的技术问题是：为了解决目前尚无法通过试验手段获取准确的滚动体滑移速度，进而无法判定滚动体的打滑现象的问题，本发明提供了一种无保持架球轴承滚动体打滑的数值分析方法，以无保持架球轴承为对象，通过建立无保持架球轴承滚动体的动力学模型并求解，在不同轴承载荷、不同轴承转速，对无保持架球轴承滚动体出现打滑进行分析，进而判定滚动体的打滑现象。

本发明解决其技术问题所采用的技术方案是：

一种无保持架球轴承滚动体打滑的数值分析方法，包括如下步骤：

1)建立滚动体的公转角加速度和自转角加速度的动力学模型，具体为：

式中，ψ_j为滚动体绕轴承回转中心角度，

为公转角加速度，由滚动体绕轴承回转中心角度ψ_j的二阶导数表示；θ_j为滚动体绕质心旋转角度，

为自转角加速度，由滚动体绕质心旋转角度θ_j的二阶导数表示；

F_ij为滚动体与内圈的摩擦力；F_oj为滚动体与外圈的摩擦力；F_bf(j，j+1)为滚动体与后一个滚动体的摩擦力；F_bf(j-1，j)为滚动体与前一个滚动体的摩擦力；

为滚动体半径；J_r为滚动体惯性矩；

为轴承内圈沟底半径；

为外圈沟底半径；G_r为滚动体重力；

为轴承节圆半径；F_{impact(j，j+1)}为滚动体与后一个滚动体的碰撞力；F_{impact(j-1，j)}为滚动体与前一个滚动体的碰撞力；

为

为

2)求解滚动体的公转角加速度和自转角加速度的动力学模型，求得公转角速度和自转角速度；

3)构建滚动体与内圈的的滑移速度公式，即：

式中：V_(ij)rel为滚动体相对轴承内圈的滑移速度；ω_i为轴承内圈转速；ω_cj为滚动体公转角速度；ω_bj为滚动体自转角速度；

4)构建滚动体与外圈的滑移速度公式，即：

式中，V_(oj)rel为滚动体相对于轴承外圈的滑移速度；

5)将公转角速度和自转角速度带入公式(3)和(4)中，得到滚动体与轴承内、外圈的滑移速度；

6)根据滚动体与轴承内、外圈的滑移速度随滚动体位置角的变化曲线评价分析滚动体的打滑情况。

具体的，滚动体与内圈的摩擦力F_ij具体为：

式中，

N_i为滚动体在和内、外圈接触点的法线方向上受到的内圈的压力。

具体的，所述滚动体在和内、外圈接触点的法线方向上受到的内圈的压力N_i具体为：

式中，

F_r为轴承所受径向载荷，Z为滚动体的个数，α为轴承接触角，ε为载荷分布系数，ψ为滚动体位置角，第j个滚动体的位置角为ψ_j。

具体的，所述滚动体与外圈的摩擦力F_oj具体为：

式中，

N_o为外圈对滚动体的反作用力。

具体的，所述外圈对滚动体的反作用力N_o具体为：

式中，F_cj为滚动体离心力，G_r为滚动体的自身重力，ψ_j为为第j个滚动体位置角，F_bf(j,j+1)＝-μF_{impact(j,j+1)}，F_bf(j-1,j)＝-μF_{impact(j-1,j)}，

为

为

所述滚动体离心力F_cj为：

式中ω_cj为滚动体公转角速度，m为滚动体的质量。

具体的，所述滚动体与后一个滚动体的碰撞力F_{impact(j,j+1)}为：

式中，δ为滚动体碰撞后的变形量，

K为刚度系数，

c为阻尼(为刚度的1％)；

式中，η为两滚动体综合弹性常数，

E₁，E₂为滚动体材料的弹性模量(在这里滚动体材料一样则E₁＝E₂)，

为与接触点主曲率相关的系数(经验值为1)，ρ为接触点主曲率(在这里为定值4/4.7625)。

本发明的有益效果是：本发明提供的一种无保持架球轴承滚动体打滑的数值分析方法，可以得到不同轴承转速、载荷下滚动体与轴承内外圈的滑移速度，如图5、图6所示，当滑移速度不等于0时即表示滚动体出现打滑现象，滚动体处于(2π/3,3π/2)的角位置时，滚动体相对于轴承内、外圈的滑移速度为0，滚动体处于纯滚动状态即没有发生打滑现象；当滚动体处于(0,2π/3)(3π/2,0)的角位置时，滚动体相对于轴承内、外圈的滑移速度不恒为0，此时滚动体产生打滑现象；滑移速度小于0时表示滚动体相对于内外圈，往运动的相反方向滑移，滑移速度大于0时表示滚动体相对于轴承内、外圈，往运动的相同方向滑移(由滚动体重力或滚动体碰撞作用产生)。

附图说明

图1：滚动体的公转角加速度和自转角加速度的动力学模型求解分析图；

图2：滚动体受力与运动示意图；

图3：滚动体公转角速度；

图4：滚动体自转角速度；

图5：极坐标系下不同转速下滚动体与轴承内外圈的滑移速度；

图6：极坐标系下不同载荷下滚动体与轴承内外圈的滑移速度。

具体实施方式

现在对本发明作进一步详细的说明。

首先，分析滚动体的受力和运动，如图2所示，首先分析单个滚动体(第j个球)的受力，在工作中轴承内圈会受到来自主轴所施加的载荷，受载的内圈会向下对滚动体施加力，因此滚动体在和内、外圈接触点的法线方向上会受到内圈的压力(N_i)和外圈对滚动体的反作用力(N_o)。由于滚动体在轴承中相对于内圈和外圈是运动的，并且滚动体与内外圈是相互接触的，因此滚动体在与内外圈接触点的切线方向会受到来自于内、外圈的摩擦力(F_o和F_i)。单个滚动体有可能会受到来自相邻滚动体的碰撞力(F_impact)，当相邻滚动体发生碰撞时，由于滚动体是运动的且在碰撞点的运动速度方向相反，会由此产生滚动体间的摩擦力(F_bf)。由于滚动体绕质心公转，这会产生离心力(F_c)。滚动体的自身重力为(G_r)。

再分析单个滚动体的运动，轴承外圈固定不动，内圈逆时针旋转驱动所有滚动体运动，滚动体的运动包括：滚动体绕轴承回转中心的公转(ω_c)和绕滚动体质心的自转(ω_b)。其公转角加速度可由滚动体绕轴承回转中心角度(ψ_j)的二阶导数

表示，其自转角加速度可由滚动体绕质心旋转角度(θ_j)的二阶导数

表示。

一种无保持架球轴承滚动体打滑的数值分析方法，如图1所示，包括如下步骤：

1)建立滚动体的公转角加速度和自转角加速度的动力学模型，具体为：

式中，ψ_j为滚动体绕轴承回转中心角度，

为公转角加速度，由滚动体绕轴承回转中心角度ψ_j的二阶导数表示；θ_j为滚动体绕质心旋转角度，

为自转角加速度，由滚动体绕质心旋转角度θ_j的二阶导数表示；

F_ij为滚动体与内圈的摩擦力；F_oj为滚动体与外圈的摩擦力；F_bf(j-1，j)为滚动体与后一个滚动体的摩擦力；F_bf(j-1，j)为滚动体与前一个滚动体的摩擦力；

为滚动体半径；J_r为滚动体惯性矩；

为轴承内圈沟底半径；

为外圈沟底半径；G_r为滚动体重力；

为轴承节圆半径；F_{impact(j，j+1)}为滚动体与后一个滚动体的碰撞力；F_{impact(j-1，j)}为滚动体与前一个滚动体的碰撞力；

为

为

2)求解滚动体的公转角加速度和自转角加速度的动力学模型，求得公转角速度和自转角速度；

3)构建滚动体与内圈的的滑移速度公式，即：

式中：V_(ij)rel为滚动体相对轴承内圈的滑移速度；ω_i为轴承内圈转速；ω_cj为滚动体公转角速度；ω_bj为滚动体自转角速度；

4)构建滚动体与外圈的滑移速度公式，即：

式中，V_(oj)rel为滚动体相对于轴承外圈的滑移速度；

5)将公转角速度和自转角速度带入公式(3)和(4)中，得到滚动体与轴承内、外圈的滑移速度；

6)根据滚动体与轴承内、外圈的滑移速度随滚动体位置角的变化曲线评价分析滚动体的打滑情况。

具体的，滚动体与内圈的摩擦力F_ij具体为：

式中，

N_i为滚动体在和内、外圈接触点的法线方向上受到的内圈的压力。

具体的，所述滚动体在和内、外圈接触点的法线方向上受到的内圈的压力N_i具体为：

式中，

F_r为轴承所受径向载荷，在本具体实施例中取值为1000N，2000N，3000N，4000N，Z为滚动体的个数，α为轴承接触角，在本具体实施例中取值为0，ε为载荷分布系数，在本具体实施例中取值为1，ψ为滚动体位置角。

具体的，所述滚动体与外圈的摩擦力F_oj具体为：

式中，

N_o为外圈对滚动体的反作用力。

具体的，所述外圈对滚动体的反作用力N_o具体为：

式中，F_cj为滚动体离心力，G_r为滚动体的自身重力，ψ_j为第j个滚动体位置角，F_bf(j,j+1)＝-μF_{impact(j,j+1)}，F_bf(j-1,j)＝-μF_{impact(j-1,j)}，

为

为

滚动体离心力F_cj为：

式中ω_cj为滚动体公转角速度，m为滚动体的质量。

具体的，滚动体与后一个滚动体的碰撞力F_{impact(j,j+1)}为：

式中，δ为滚动体碰撞后的变形量，

K为刚度系数，

c为阻尼(为刚度的1％)；

式中，η为两滚动体综合弹性常数，

E₁，E₂为滚动体材料的弹性模量(在这里滚动体材料一样则E₁＝E₂)，

为与接触点主曲率相关的系数(经验值为1)，ρ为接触点主曲率(在这里为定值4/4.7625)。

在本具体实施例中，无保持架球轴承模型拟选取的参数值如表1-1所示：

表1-1无保持架球轴承模型参数

滚动体的公转角加速度和自转角加速度方程程序代码为：

function[dy]＝Test804(t,y)

m＝0.0035519；％滚动体质量

G＝0.0348；％滚动体重量

Dw＝0.009525；％滚动体直径(m)

R＝Dw/2；滚动体半径(m)

Jr＝m*Dw^2/10；％滚动体转动惯量

dm＝0.04605；％节圆半径(m)

di＝0.03648；％内滚道沟底半径(m)

do＝0.05553；％外滚道沟底半径(m)

Rm＝dm/2；％节圆半径(m)

Ri＝di/2；％内滚道沟底半径(m)

Ro＝do/2；％外滚道沟底半径(m)

dy＝zeros(60,1)；

dy(1)＝y(2)；

dy(2)＝(mui(y(2),y(4))*Ni(y(1))*Ri*DeltaVi(y(2),y(4))-muo(y(2),y(4))*DeltaVo(y(2),y(4))*(Ni(y(1))+m*Rm*y(2)*y(2)-G*cos(y(1))-0.002*Contactball(y(1),y(5))*cos((y(1)-y(5))/2)+0.002*Contactball(2*pi+y(57),y(1))*cos((2*pi+y(57)-y(1))/2)+Contactball(y(1),y(5))*sin((y(1)-y(5))/2)+Contactball(2*pi+y(57),y(1))*sin((2*pi+y(57)-y(1))/2))*Ro-G*sin(y(1)+pi)*Rm+Contactball(y(1),y(5))*cos((y(1)-y(5))/2)*(Rm-R*sin((y(1)-y(5))/2))-Contactball(2*pi+y(57),y(1))*cos((2*pi+y(57)-y(1))/2)*(Rm-R*sin((2*pi+y(57)-y(1))/2)))/m/Rm/Rm；％mui代表滚动体与内圈的摩擦系数；muo代表滚动体与外圈的摩擦系数；DeltaVo是判断滚动体相对于内外圈滑移的方向，和运动方向同向为正，反向为负；Contactball表示滚动体间的碰撞力；pi代表π。

dy(3)＝y(4)；

dy(4)＝(mui(y(2),y(4))*Ni(y(1))*DeltaVi(y(2),y(4))-Contactball(y(1),y(5))*0.002-Contactball(2*pi+y(53),y(1))*0.002)*Dw/2/Jr+(muo(y(2),y(4))*DeltaVo(y(2),y(4))*(Ni(y(1))+m*Rm*y(2)*y(2)-G*cos(y(1))+0.002*Contactball(y(1),y(5))*cos((y(1)-y(5))/2)-0.002*Contactball(2*pi+y(57),y(1))*cos((2*pi+y(57)-y(1))/2)+Contactball(y(1),y(5))*sin((y(1)-y(5))/2)+Contactball(2*pi+y(57),y(1))*sin((2*pi+y(57)-y(1))/2)))*R/Jr；

end

程序2

微分方程求解程序

w_i＝100*pi；％轴承内圈转速

Dw＝0.009525；％滚动体直径

R＝Dw/2；滚动体半径(m)

dm＝0.04605；％节圆半径(m)

di＝0.03648；％内滚道沟底半径(m)

do＝0.05553；％外滚道沟底半径(m)

Rm＝dm/2；％节圆半径(m)

Ri＝di/2；％内滚道沟底半径(m)

Ro＝do/2；％外滚道沟底半径(m)

w_r＝R_i*R_o*w_i/(2*R*Rm)；％滚动体自转角速度

w_m＝R_i*w_i/(2*Rm)；％滚动体公转角速度

xspan＝[0,1]；％设置求解时间

s＝2*asin(R/Rm)；％相邻滚动体接触时的质心到回转中心的夹角

y0＝[14*s,w_m,0,w_r,13*s,w_m,0,w_r,12*s,w_m,0,w_r,11*s,w_m,0,w_r,10*s,w_m,0,w_r,9*s,w_m,0,w_r,8*s,w_m,0,w_r,7*s,w_m,0,w_r,6*s,w_m,0,w_r,5*s,w_m,0,w_r,4*s,w_m,0,w_r,3*s,w_m,0,w_r,2*s,w_m,0,w_r,s,w_m,0,w_r,0,w_m,0,w_r]；

％微分方程初始值的确定(数字代表滚动体编号从0开始)

h＝1e-6；％设置积分步长(每间隔10^-6s计算一次)

[t,y]＝Adams4(@Test804,xspan,y0,h)；％Adams算法；Test804表示在此处调用程序1。

获得滚动体公转角速度如图2所示，滚动体自转角速度如图3所示。

程序3

不同轴承内圈转速下(1500r/min、2000r/min、3000r/min、4500r/min)，极坐标系下，滚动体滑移速度曲线随滚动体位置角变化的绘制程序

w_i1＝50*pi；

w_i2＝200/3*pi；

w_i3＝100*pi；

w_i4＝150*pi；

％代表不同转速1500r/min、2000r/min、3000r/min、4500r/min

figure('color',[1 1 1])

g_1＝y1500(:,17)；

x_1＝find(g_1>14.5*2*pi&g_1<15.5*2*pi)；％滚动体转动到第14.5-15.5圈的数据。

g_1(x_1)；

t＝mod(g_1(x_1),2*pi)；％将滚动体时间坐标系下的角度变化转换成以2π为周期的极坐标下的角度变化。

a＝x_1(1:1)；

b＝x_1(end:end)；

I_1＝(R_m-R)*(w_i1-y1500(a:b,18))-R*(y1500(a:b,20))；％转速1500r/min下滚动体的滑移速度

polarplot(t,I_1,'b','linewidth',1.2)

hold on

g_2＝y2000(:,17)/(2*pi)；

x_2＝find(g_2>14.5&g_2<15.5)；

g_2(x_2)；

c＝x_2(1:1)；

d＝x_2(end:end)；

I_2＝(R_m-R)*(w_i2-y2000(c:d,18))-R*(y2000(c:d,20))；％转速2000r/min下滚动体的滑移速度

polarplot(g_2(x_2)*2*pi,I_2,'r','linewidth',1.2)

hold on

g_3＝y3000(:,17)/(2*pi)；

x_3＝find(g_3>14.5&g_3<15.5)；

g_3(x_3)；

e＝x_3(1:1)；

f＝x_3(end:end)；

I_3＝(R_m-R)*(w_i3-y3000(e:f,18))-R*(y3000(e:f,20))；％转速3000r/min下滚动体的滑移速度

polarplot(g_3(x_3)*2*pi,I_3,'g','linewidth',1.2)

hold on

g_6＝y4500(:,17)/(2*pi)；

x_6＝find(g_6>14.5&g_6<15.5)；

g_6(x_6)；

u＝x_6(1:1)；

w＝x_6(end:end)；

I_6＝(R_m-R)*(w_i4-y4500(u:w,18))-R*(y4500(u:w,20))；％转速4500r/min下滚动体的滑移速度

polarplot(g_6(x_6)*2*pi,I_6,'k','linewidth',1.2)

pax＝gca；

pax.ThetaAxisUnits＝'radians'；

pax.ThetaZeroLocation＝'top'；

pax.FontName＝'Times New Roman'；

pax.FontSize＝12；

F＝legend('1500r/min','2000r/min','3000r/min','4500r/min')；

set(F,'FontName','Times New Roman','FontSize',12)

legend('boxoff')

rlim([-0.1 0.3])

程序4

％不同径向载荷下(1000N，2000N，3000N，4000N)，极坐标系下，滚动体滑移速度曲线随滚动体位置角变化的绘制程序：

figure('color',[1 1 1])

g_4＝y1000(:,17)/(2*pi)；

x_4＝find(g_4>14&g_4<15)；

g_4(x_4)；

m＝x_4(1:1)；

n＝x_4(end:end)；

I_4＝(R_m-R)*(w_i-y1000(m:n,18))-R*(y1000(m:n,20))；％载荷1000N下滚动体的滑移速度

polarplot(g_4(x_4)*2*pi,I_4,'b','linewidth',1)

hold on

g_5＝y2000(:,17)/(2*pi)；

x_5＝find(g_5>14&g_5<15)；

g_5(x_5)；

o＝x_5(1:1)；

p＝x_5(end:end)；

I_5＝(R_m-R)*(w_i-y2000(o:p,18))-R*(y2000(o:p,20))；％载荷2000N下滚动体的滑移速度

polarplot(g_5(x_5)*2*pi,I_5,'r','linewidth',1)

hold on

g_6＝y3000(:,17)/(2*pi)；

x_6＝find(g_6>14&g_6<15)；

g_6(x_6)；

u＝x_6(1:1)；

w＝x_6(end:end)；

I_6＝(R_m-R)*(w_i-y3000(u:w,18))-R*(y3000(u:w,20))；％载荷3000N下滚动体的滑移速度

polarplot(g_6(x_6)*2*pi,I_6,'g','linewidth',1)

hold on

g_7＝y4000(:,17)/(2*pi)；

x_7＝find(g_7>14&g_7<15)；

g_7(x_7)；

u＝x_7(1:1)；

w＝x_7(end:end)；

I_7＝(R_m-R)*(w_i-y4000(u:w,18))-R*(y4000(u:w,20))；％载荷4000N下滚动体的滑移速度

polarplot(g_7(x_7)*2*pi,I_7,'k','linewidth',1)

pax＝gca；

pax.ThetaAxisUnits＝'radians'；

％pax.ThetaDir＝'clockwise'；

pax.ThetaZeroLocation＝'top'；

pax.FontName＝'Times New Roman'；

pax.FontSize＝12；

F＝legend('1000N','2000N','3000N','4000N')；

set(F,'FontName','Times New Roman','FontSize',12)

legend('boxoff')

rlim([-0.02 0.1])

得到不同轴承转速、载荷下滚动体与轴承内外圈的滑移速度，如图5和图6所示，当滑移速度不等于0时即表示滚动体出现打滑现象，滚动体处于(2π/3,3π/2)的角位置时，滚动体相对于轴承内、外圈的滑移速度为0，滚动体处于纯滚动状态即没有发生打滑现象；当滚动体处于(0,2π/3)(3π/2,0)的角位置时，滚动体相对于轴承内、外圈的滑移速度不恒为0，此时滚动体产生打滑现象。滑移速度小于0时表示滚动体相对于内外圈，往运动的相反方向滑移，滑移速度大于0时表示滚动体相对于轴承内、外圈，往运动的相同方向滑移(由滚动体重力或滚动体碰撞作用产生)。

以上述依据本发明的理想实施例为启示，通过上述的说明内容，相关工作人员完全可以在不偏离本项发明技术思想的范围内，进行多样的变更以及修改。本项发明的技术性范围并不局限于说明书上的内容，必须要根据权利要求范围来确定其技术性范围。

Claims (7)

1.一种无保持架球轴承滚动体打滑的数值分析方法，包括如下步骤：

1)建立滚动体的公转角加速度和自转角加速度的动力学模型，具体为：

式中，ψ_j为第j个滚动体位置角，

为公转角加速度，由第j个滚动体位置角ψ_j的二阶导数表示；θ_j为滚动体绕质心旋转角度，

为自转角加速度，由滚动体绕质心旋转角度θ_j的二阶导数表示；

F_ij为滚动体与内圈的摩擦力；F_oj为滚动体与外圈的摩擦力；F_bf(j，j+1)为滚动体与后一个滚动体的摩擦力；F_bf(j-1，j)为滚动体与前一个滚动体的摩擦力；

为滚动体半径；J_r为滚动体惯性矩；

为轴承内圈沟底半径；

为外圈沟底半径；G_r为滚动体重力；

为轴承节圆半径；F_{impact(j，j+1)}为滚动体与后一个滚动体的碰撞力；F_{impact(j-1，j)}为滚动体与前一个滚动体的碰撞力；

为

为

m为滚动体的质量；

2)求解滚动体的公转角加速度和自转角加速度的动力学模型，求得公转角速度和自转角速度；

3)构建滚动体与内圈的滑移速度公式，即：

式中：V_(ij)rel为滚动体相对轴承内圈的滑移速度；ω_i为轴承内圈转速；ω_cj为滚动体公转角速度；ω_bj为滚动体自转角速度；

4)构建滚动体与外圈的滑移速度公式，即：

式中，V_(oj)rel为滚动体相对于轴承外圈的滑移速度；

5)将公转角速度和自转角速度带入公式(3)和(4)中，得到滚动体与轴承内、外圈的滑移速度；

6)根据滚动体与轴承内、外圈的滑移速度随滚动体位置角的变化曲线评价分析滚动体的打滑情况。

2.如权利要求1所述的一种无保持架球轴承滚动体打滑的数值分析方法，其特征在于：滚动体与内圈的摩擦力F_ij具体为：

式中，

N_i为滚动体在和内、外圈接触点的法线方向上受到的内圈的压力。

3.如权利要求2所述的一种无保持架球轴承滚动体打滑的数值分析方法，其特征在于：所述滚动体在和内、外圈接触点的法线方向上受到的内圈的压力N_i具体为：

式中，

F_r为轴承所受径向载荷，Z为滚动体的个数，α为轴承接触角，ε为载荷分布系数，ψ_j为滚动体位置角。

4.如权利要求1所述的一种无保持架球轴承滚动体打滑的数值分析方法，其特征在于：所述滚动体与外圈的摩擦力F_oj具体为：

式中，

N_o为外圈对滚动体的反作用力。

5.如权利要求4所述的一种无保持架球轴承滚动体打滑的数值分析方法，其特征在于：所述外圈对滚动体的反作用力N_o具体为：

式中，F_cj为滚动体离心力，G_r为滚动体的自身重力，ψ_j为第j个滚动体位置角，F_bf(j,j+1)＝-μF_{impact(j,j+1)}，F_bf(j-1,j)＝-μF_{impact(j-1,j)}，

为

为

N_i为滚动体在和内、外圈接触点的法线方向上受到的内圈的压力。

6.如权利要求5所述的一种无保持架球轴承滚动体打滑的数值分析方法，其特征在于：所述滚动体离心力F_cj为：

式中ω_cj为滚动体公转角速度，m为滚动体的质量。

7.如权利要求1所述的一种无保持架球轴承滚动体打滑的数值分析方法，其特征在于：所述滚动体与后一个滚动体的碰撞力F_{impact(j,j+1)}为：

式中，δ为滚动体碰撞后的变形量，

K为刚度系数，

c为阻尼；

式中，η为两滚动体综合弹性常数，

E₁，E₂为滚动体材料的弹性模量，

为与接触点主曲率相关的系数，ρ为接触点主曲率，v₁、v₂均为泊松比。

Images