角接触球轴承

1、角接触球轴承打滑行为的非线性动态模型QinkaiHan,FuleiChu.TheStateKeyLaboratoryofTribology,TsinghuaUniversity,Beijing100084,China.摘要：用一个三维非线性动态模型来预测复合载荷组合条件下角接触球轴承的打滑行为。该模型考虑了钢球的自转和公转引起的离心力和陀螺效应、钢球与内外圈之间的赫兹接触变形、钢球与保持架之间的非连续接触以及弹流动体润滑。通过对试验结果的比较，验证了该动态模型正确性。在此基础上，讨论了在复合载荷作用下，轴承钢球滑动速度随时间和位置的变化规律。该模型表明，径向载荷的变化将使钢球在内外圈之间的的滑

2、动速度产生波动，对低负载区域的钢球影响更大。增加径向负荷将大幅增加滑移速度的幅度和范围，使打滑更加严重。当钢球在低载区时，大的滑动速度会使轴承和润滑油的温度升高，加剧轴承磨损，缩短轴承的使用寿命。因此，在旋转工件的设计和检测中应考虑径向载荷。1.导论：角接触球轴承是许多旋转机械的核心支撑部件，其动态特性对整个设备的使用性能、运行可靠性和使用寿命起着决定性的作用。轴承在运行过程中，滚道应为钢球提供足够大的摩擦力和摩擦力矩，以确保钢球处于纯滚动状态。否则，滚动体和内、外滚道之间可能会出相对滑移。随着现代旋转机械的高速化、重载化，轴承的滑动将使轴承和润滑油的温度升高，从而加速轴承磨损。如果轴承早期就

3、开始打滑，它可能会导致轴承寿命减少，甚至更严重的事故。因此，当前准确预测滚动轴承的打滑行为并提出防滑设计准则是很重要的问题。哈里斯1,2已经在这方面做了开创性的工作。基于沟道控制理论和准静态学，哈里斯1,2建立了用于高速角接触球轴承的滑行预测模型。该模型考虑了滚动体的各种受力情况（包括：接触力，摩擦力，流体力和离心力等），还考虑了轴向载荷、旋转速度、滚动体的数量对打滑的影响。LiaoandLin3在几何约束条件和受力平衡中考虑了每一个滚动体受到的接触力和每一个滚动体的接触角。希拉诺4打滑的评判标准中分析了在轴向和径向负荷下钢球和滚道的打滑。此外，他们还研究了热效应引起的钢球滑动5。基于准动态分

4、析，Jiangetal.6,Cuietal.7andYuanetal.8等提出了估计防止轴承打滑的最小轴向载荷的经验公式。最近，Chenetal.9,10提出推力球轴承在固体润滑条件下的准静态模型。该理论构想了一种用于准确界定钢球与内外圈之间相对运动的滑动比和旋滚比。它表明，钢球与滚道的密切接触引起的滑差使滑动和接触力分布不对称。Xuetal.11建议用预载分析法作为球轴承的滑动准则。实验结果表明采用最佳预紧力的轴承具有良好的温度特性。Chenetal12等人对高速旋转的轴承进行准静态打滑分析，发现钢球自转挤压油膜对打滑和轴承疲劳寿命有不利的影响。上述的大多数分析是在稳态条件下采用静态/动态模

5、型来研究打滑行为和防滑判据。实际上，滚动轴承上经常被施加动载荷，滚动体和保持架之间的接触和碰撞是不可避免的。这些因素往往造成滑移随时间和空间的变化。显然，基于静态模型的稳态分析很难准确地描述和预测滚动体打滑的行为。因此，发展动态分析方法在当前得到了广泛关注。采用动态法，古普塔13首先建立微分方程来模拟一个推力球轴承在弹流润滑条件下钢球的瞬时运动。奥斯滕森等人14建立了轴承运动的动力学模型，并对轴承运动进行了仿真模拟，并与测量结果进行比较。他们发现，卸载区的滑动主要受局部润滑剂分布状态和滚子在轴承中的位置的影响。Imado15提出了一个霍尔元件检测法来检测钢球在轴承中的运动。最近，米尼安拉尼阿多

6、哈科梅等人16利用有限元法建立了滚动轴承的动态接触模型。动态接触模型应用在钢球与滚道打滑的研究中。Jain和亨特17忽略离心力和钢球接触角引起的陀螺力的影响，提出了一种应用在风力机传动系统高速输出轴上的轴承的动力学模型。塞尔瓦拉和marappan18研发了一种圆柱滚子轴承试验台，用于测量各种工况下运行的轴承元件的速度，在实验结果的基础上讨论了工作参数对圆柱滚子轴承的影响。Tuetal.19以滚动体与滚道、保持架之间的接触力和摩擦力，以及滚动体的离心力为基础建立了研究滚动轴承加速时打滑的分析模型。对轴承打滑的动态问题的一些研究工作已在目前的文献中提出，但研究对象仅限于推力轴承（只适用于轴向载荷）

7、和深沟轴承（仅适用于径向载荷）。角接触轴承在轴向和径向载荷作用下（复合载荷条件下）的研究较少受到关注。高速度和复合载荷工况会使角接触轴承的滚动元件具有三维运动，而推力轴承和深沟轴承的运动可以简化为二维运动。因此，本论文提出分析角接触球轴承打滑行为的动力学模型。首先，考虑了内圈的五个自由度，对轴承的载荷分布和每个滚动体的接触力和内/外接触角进行了分析。根据赫兹接触理论和弹流动力润滑，确定了钢球与内/外滚道之间的摩擦力和摩擦力矩随时间的变化规律。考虑到钢球和保持架之间的不连续接触，利用欧拉方程建立了角接触滚动轴承的三维非线性动力学模型。本论文给出了模型的求解方法和过程并通过与实验结果的比较对其进行

8、了验证，在此基础上，讨论了在复合载荷作用下，轴承的滑动速度随时间和空间的变化规律。最后，得出了一些结论。2、角接触球轴承载荷分布分析通过对角接触球轴承载荷分布的分析，获得了三种类型的参数，包括接触力，赫兹接触面积的尺寸和钢球和内外圈之间的接触角。钢球摩擦力矩的大小和方向可以通过这些参数分析计算得到。因此，载荷分布分析在打滑分析之前进行。外圈固定，内圈以恒定的角速度i旋转。两个坐标，其中一个（X-Y-Z）固定在内圈的旋转中心线上，另一个（x-y-z）随Z轴旋转，如图1所示。x-y-z的中心与球的质心重合。忽略了内圈的旋转自由度，重点考虑余下的五个自由度（包括：三平移和两旋转）。自由向量用品=际,

9、看"皿脩唯亍表示。作用于内圈的外力向量用Fin=(FXin;FYin;FZin;MXin;MYin）T表示。第i个球相对于旋转坐标x-y-z有三个旋转角速度（xj;yj;zj）和一个基于固定轴Z的滚道旋转速度Wcj。在时间t,第j个球的位置角用下公式表示：r7图1。轴承坐标示意图。?cj=?0+cj+2冗/Nb(j-1),其中？0为初始位置角，cj为滚道角速度，Nb为钢球的数量。一般情况下？0=0,外力会使内圈移动。因此，第j个球内圈沟曲率中心的位移，包括沿y,z轴的移动和沿x轴的转动，用uj=(yj;zj;9xj)T表示。在小挠度的假设下，uj可以通过坐标变换得到uj=TjUin(

10、1)其中Tj是变换矩阵，可表示为0-此cos丸-COS巾tj8s祖。sin小叮0000其中Ri=0.5de+(rgi-0.5db)cosB0表示由内到槽曲率中心的内轴旋转中心的距离，de是轴承的节圆直径，rgi是沟曲率半径，db是钢球的直径，B0是公称接触角。没有负载时，钢球与内外滚道之间没有接触变形。由于偏转是零，所以球的中心和内外滚道的沟曲率中心位于一条直线，如图2虚线。钢球的内外接触角是相同的，等于公称接触角B0。施加负载后，内圈滚道将首先移动，然后钢球也因接触变形而移动。平衡状态，如图2所示的实心线，可以表示钢球和内外滚道之间的状态。因为外圈是固定的，所以它的沟曲率中心保持不变。根据公

11、式(1)内圈沟曲率中心的位移可以用(yj;zj)表示。钢球中心的位移是未知参数可用(vjy;vjz)表示。因为负载的影响，钢球与内外沟道的接触角变的不相等。外接触角减小为Boj,而内接触角增大为Bij。在偏转前，内圈沟曲率中心到钢球质心的距离为Lij。偏转后，距离变得lij,如图2所示。同样，外圈沟曲率中心到钢球质心的距离从Loj(偏转前）变为loj（偏转后），对于给定的内外沟曲率半径Lij=rgi-0.5db和Loj=rgo-0.5db。根据图2的几何关系,（rgi;rgo）,有可以获得以下公(4)_Lsing+v%一仙九荷载分布分析应满足上述的四个几何方程。图3给出了第j个钢球的受力情况。

12、在图中,Fcj和Mgj表示由于球的自转和公转引起的离心力和陀螺力。钢球质心钢球质心（偏转后）内圈沟曲率中心（偏转前）内圈沟曲率中心*（偏转后）f角iI凡外圈沟曲率中心C固定）轴承前后挠度的几何关系小微等I图3,第j钢球受显然有。=。.5口3和M田=m5gsin%时质量和惯性。s和c为自转角速度和沟道角速度,轴之间的夹角。对于纯滚动轴承，有力图。20mb和Ib为钢球的瞬aj为自转轴和Z亚=蜜（】一（*）"）助，砥=°5（1一%）叫%=由一】（丽黑海）钢球与内、外滚道之间的接触力用Qij和Qoj表示，根据赫兹接触理论，有Qij=xijKi6ij3/2和Qoj=xojKo6oj3

13、/2。Ki和Ko为钢球与内外滚道的接触刚度系数，6ij和6oj为钢球与内外滚道的变形量。从上面的几何分析中可以得到6ij=lij-Lij和6Oj=loj-Loj。当6ij,6oj>0；xij=1,xoj=1o6ij,6oj三0时；Xij=0,xoj=0o接触刚度系数通过接触区材料性能和几何尺寸的确定，即k=K=（兀在13G，3£况/号，其中E'=E/（1-V2）表示的有效弹性模量，E为材料的弹性模量，V为泊松比，R为等效曲率半径，k为椭圆率，己和？为椭圆第一类和第二类积分。这些参数的表达式从哈里斯和kotzalas20的论文中可以得到,接触刚度系数Ki和Ko可以通过以上

14、公式求解得到。在负荷分布分析，忽略了x水平的摩擦力，陀螺力矩由y-z平面的先对钢球进行这样的受力平衡分析Aij和入oj表示钢球与内外沟道之间的瞬时摩擦力产生，如图3所示。哈里斯和kotzalas20在他们的论文中首摩擦系数。在这里对入ij和入oj都设置初始值为1。从图3可以得到力的平衡方程如下：Qjjcos%-Qq85凡j+sin%一'2却sin/?可+%=。(7)Q斤sin*厂%亢n#町-甘%cos%+缗%ras/?0)=0(8)口0ufj第j个钢球与内沟道之间的接触力和摩擦力的方向与图3所示的方向相反。把这些力等效到内沟曲率中心可以得到.Q州一Q”cos%_sin%4=Q"

15、;允n#jj+峻。5/斤(9)所JQyJl旺订0=0J=1Qxj得到钢球与内滚道的力平衡方程(10)上述研究表明在未知参数负荷分布的位移分析中可以用该钢球与内滚道的模型(vjy,vjz和Xin,Yin,Zin,0Xin,OYin)代替进行计算。(公式(7)-(10)是一组非线性代数方程，可以用Newton-Raphson方法迭代求解。该非线性方程的维数是2Nb+5.得到这些未知参数后，每一个钢球的接触力(Qij,Qoj)和内外接触角(Bij;Boj)可以通过公式(1)和(3)-(6)计算得到。3.施加在球上的摩擦力矩和力矩运行中滚动轴承，内外沟道对钢球产生摩擦力和摩擦力矩。弹流润滑条件下的钢球

16、和滚道之间的接触与图4所示。图中的x',v'轴分别表示接触椭圆的主要和次要轴，Z'表示垂直于接触椭圆的轴。根据赫兹接触理论，接触应力在接触区为椭圆分布，并可以表示如下：其中，PHj为最大接触应力，a和b为接触椭圆的长半轴和短半轴。有Q=(印了.、K-FJ.b=a/k。在求出接触力，等效曲率半径和椭圆率后，可以求出接触区的面积(ain,bin和aout,bout)。一般情况下，油膜厚度hc远小于接触区尺寸,接触区域任一点的剪切应力r(x',y')可以用下公式表小：式M,V)=q(Pj,T)*(12)其中“(pj,T)是润滑剂的粘度，是接触压力Pj和油的温度

17、T的函数。相对之间的接触界面的滑动速度用表Au示，Au=u1-u2。图4.对于钢球和滚道之间的接触示意图。根据结果Hamrock和道森21,油膜厚度是由下面的经验公式计算：1=2.69七67Gli53w-Dgq(13)IJ_n&nl其中一E七-表示无量纲速度，刀0是20c大气压力下润滑油的粘度，&=d»2"土/48S川o)表示沿滚动方向的等效半径，外圈用“+”，内圈用“-"。.初一仁%是等效转速，dri,dro表示内圈和外圈沟道的直径，G=E'cqp是无量纲的弹性模量，其中cqp为粘度压力系数。无量纲载荷由计算得到0二藁。润滑油的粘度“（p

18、j,T）由经典的巴勒斯方程17来确定。然而，目前的研究表明在接触压力大的情况下使用该公式将会产生较大的误差。滚动轴承的接触压力通常大于1GPa,属于高接触压力。因此，在这种情况下，巴勒斯方程不能用。贝尔和科特基22的实验研究了润滑剂的接触压力和油温在高接触压力（高达2GPa）下的变化情况。以下的经验公式是通过数据拟合得到的。流岛-息）（14）其中B,R0,R分别润滑油的杜利特-泰特参数，r为油量随温度的线性变化，即r=1+?1（T-T0），其中？1为20c室温下体积的膨胀量。V/V0为体积随压力和温度的变化率，公式如下：(15)中a为热膨胀系数，K0为体积模数，K'0为常数用%=仁+&

19、amp;”,计算,这里的温度是Kelvin（绝对温度）。接触界面之间的相对滑动速度Au有钢球的自转和公转速度以及内外圈的旋转速度确定。图5给出了钢球相对于x,y,z轴的三个旋转速度（xj；yj;zj）和一个相对于轴心的轨道速度（cj)。I(18对=%cos/30j+a)zjsin夕oj沿短轴（xI轴）和长轴（yl轴）的相对滑动速度可以用以下公式表示:图5.钢球与内外滚道之间的运动关系图对于钢球与yo轴）的滑移速度可(16)117)当一口口ut<y<口Wt3nd-bom<X<bout时,soj为钢球的相对于zo轴的旋转速度，有以下公式内圈转速为I,钢球相对于内圈滚道的速度

20、为cj0同样的，钢球"双=牛（"s'n/一3不COS角）一T叫一伙Jfy-y8$用）一日F保电tj十04炉1fJHlMTtI-1*rIbcrMricrrvn外滚道的接触，短轴（xo,钢球的滚动方向）和长轴（以用下边的公式计算酗i*!r|国IiEIbCmiIjl-Iiwlb?iHarrr*c¥当一口inWyXQndnd_瓦门MXMbjn时，sij为钢球相对于zi轴的旋转速度有以下公式：21)o)sij-%cos时+&)ZjsinfJ®BV)的剪知道油膜厚度，润滑油的粘度和相对滑动速度后，接触区任意点切应力P(x',y')可以

21、按下公式求解得到。1/二广心道dydx(22)ykfewt/口皿/&A»a(23)(24)加=一匚二：)嘤切也摩擦力矩的推导如下：moj=-03Ky-Au叩大+外呵(/+必0八"-totfJ-dour£(25)(26)广广用pT)口4一311队十外孤炉+/)dydx(27)4,三维非线性动力学模型在旋转坐标(x-y-z)中，第j个钢球的运动包括三个旋转运动xj,yj,zj和一个沿滚道绕轴心的运动cj。钢球的旋转矢量定义如下：其中J4+T必q/=mj(28),钢球的惯性矩为Ixj=Iyj=Izj=Ib,Tvj为牵连角速度矩阵(29)Z为旋转坐标系（x-y-z

22、）相对于固定坐标系（X-Y-Z）的旋转速度。可以得到xX=wY=0和Z=wcjoThem尸啊卜M"为施加在滚动钢球上的总外力矩。这些力矩由钢球与内外沟道之间的摩擦力和摩擦力矩计算得到。根据前面的分析和图的几何关系，可以得出如下的公式：(30)(31)(32)削-f叫sin乩j+f闻sin%）+风3%+叫皿岛）9mMmSin乩j+m.sinftj对于钢球的在沟道上的运动，钢球表面的摩擦力矩与内外沟道的摩擦力有关，在钢球的轨道运动方程中也应该考虑润滑剂的阻力和钢球与保持架之间的不连续接触。哈里斯和kotzalas20给了一个计算球轴承润滑油阻力的经验公式fv=32g(33)其中，Pv为润

23、滑油的密度，cv为阻力系数，g为重力加速度钢球和保持架之间的相对位移6+,6-如图6所示。其中t为时间，？cage为保持架的角位移，cr为钢球与保持架之间的间隙。根据图6的2-2cos舒_3）下丁面）+0几何关系，可以得到计算6+,6-的公式：（均因此，保持架对第j钢球的接触力为-fcagej=-X+)(36)kcage为钢球与保持架之间的接触刚度。当6+,6-<0时，/-=1,y(/cagej+fv)X+=1；相反，X-=0,X+=0。根据动量定理，可以将钢球的轨道运动方程写成(37)其中，Icj为钢球相对于轴线(Z轴)的转动惯量，有1十。同(38)样的保持架的的运动方程也可以利用动量

24、定理求出Zj=1其中，Icage为保持架相对于轴线的转动惯量cage为保持架的旋转角速度。从公式(34)-(35)可以看出，要想求出钢球与保持架之间的相对位移就必须先求出钢球与保持架任意时间的角位移。图6.钢球与保持架之间的相互作用示意图由角速度与角位移之间的微分关系可以得到%工%/（39）小s卷=毕（以）每一个钢球的一系列方程包括（28）,（37）,（39）在内，都能够得到，还包括两个形容保持架运动的动态方程（38）,（40）。通过计算每个钢球与滚道的运动方程，可以得到角接触球轴承的三维非线性动力学模型。（41）苴/、中,=diag（1/仃11）是（5Nb+2）X（5Nb+2）的尺寸系数矩阵

25、，X=【口：叫小。;叫要心.F是（5Nb+2）X1的三维位移矢量，为（5Nb+2）X1的外力向量。由于钢球自转和公转产生的陀螺效应和钢球与保持架之间的间断接触，公式（41）为非线性耦合一阶微分方程。外部负载的大小，钢球自转和公转的初始值，保持架和内圈的旋转速度都对方程（41）的收敛过程有显著影响。显然，外部负载的不足会导致系统响应花费更多的时间收敛，不合理的旋转和轨道速度的初始值可能会导致系统响应的振荡。本研究中，用滚道控制理论20所获得的结果作为初始值。增加内圈速度意味着陀螺耦合效应变得显着，钢球与沟道之间的不连续接触可能会更频繁。在该计算中应减少时长以保证收敛，公式（41）可用MATLAB

26、ode15s程序求解。该解算程序是专用刚性微分方程，在每一个时间段，解算程序可以估计该计算的局部误差。这种误差必须小于或等于可接受的误差，这是一个指定的相对误差（reltol）与绝对容差（abstol）之间的函数。在本研究中RelTol和AbsTol分别设置为1e-3和1e-6,当总时间超过设定时间时，计算停止，详细的计算过程如图7所示。5.模型的验证PasdariandGentle23测试了角接触球轴承的保持架在轴向载荷下的旋转速度，35毫米口径的角接触球轴承在最高转速为3000转的变速直流电机的驱动下旋转，试验轴承的参数在表一中列出。轴承被安装在驱动轴上，驱动轴由由2个深槽球轴承支撑，并通

27、过皮带轮与电机相连。轴承垂直安装，可以用重力作为轴向载荷的一部分，使其更容易作用在轴承上。本文的动态模型可以预测保持架在各种轴向载荷下的作用，其速度随时间变化的规律。将计算数值与实验结果进行比较，以验证动态模型的正确性。本试验中未给出润滑剂的类型23。本动态模型使用了一种专用于滚动轴承的典型润滑油(MIL-L-23699)。拜尔和科特基22提供了MIL-L-23699润滑剂的参数，见表2,该润滑剂的温度T为30C。7,非线性动力学模型的求解过程。试验中，内圈转速2000rpm(i=209:4395rad/s),在纯滚动条件下=8=82。0994rad/s。在模拟中，初始角速度和位移被设置为叫2

28、00ng蜘=-800rad八通=-100闯/$价.0.帆，产2工岫0,L2随and例礴=0,麻也”H轴向载荷隆为40,80,120.,300N,基于非线性动力学模型和图8中的结果，可以得到保持架速度随时间的变化规律。结果表明在<120N时，cage的影响小与3t。由于轴向载荷较低，摩擦力和摩擦力矩也会较低，不能提供足够的阻力来保证钢球一直处于纯滚动状态。表1,滚动轴承的物理参数DescriptionValuedb(m)即g)加(m)dm(m)4(m)de(m)两EN/m2)iP(kg/mjMkcge(N/m)/cage(kgm2)0,01350.00700.00700.07010.043

29、10.056625°Z075ell03780011le84.096e-5le-5。(m)表2.润滑油MIL-L-23699参数DescriptionValueno3s)a(l/Oe-4)o(l/Cle-3)%p/kg/iTp)BKo(GPa)Kq(GPaK)0.046677.42-1.280.66418603.38210741-E0149570.80.05图8.随着cage在各种轴向负荷随时间的变化。因此，宏观滑移出现在钢球与内外沟道的接触区。加大以增大对钢球的摩擦力，可以使cage接近纯滚动速度。看越大，可以cage接近而：的时间就越短，如图8中的放大图所示。把在不同轴向载荷下的保

30、持架的速度cage,与PasdariandGentle23实验结果进行比较，如图9所示。在较低的轴向载荷作用下，cage出现明显的波动，这是因为严重的打滑导致测试结果不准确。随着轴向载荷的增加，cage为恒定值，即历"曲=0.3920。图9所示的数值结果与实验结果很吻合，表明预测滑移行为的动态模型是可靠的。止匕外，从图9可以看出在Fin>120N时，cage/i接近纯滚动比0.3920,表明保持架的打滑变弱。但是，无论球是否在纯滚动状态，应通过观察钢球的自转速度（xj,yj,zj）来判断打滑的程度。以j=1的钢球为例，三个旋转速度随时间和轴向载荷的变化如图10所示。它表明，随着

31、时间的增加，在>500N后，三个转动分量变为恒定值。该图表明，当沅=520N时，x1«0rad/s,y1M87rad/s,wz1«-376rad/s。公称接触角B0=25o,钢球的旋转速度可以通过公式（18）和（20）计算得到该值近似为零。因此，当>500N时可以认为钢球是在纯滚动状态。4.轴向和径向载荷下轴承的打滑在前一节中，仅在角接触球轴承的内圈上施加轴向载荷，在实际应用中，轴向和径向载荷（联合载荷）将同时作用于内圈。与仅仅轴向加载系统相比，联合负荷对打滑行为的影响有明显的不同。在这一节中，考虑了2种类型的载荷，即：恒定的径向载荷（时不变载荷）和正弦径向载荷

32、（时变载荷）。6.1. 恒定的径向负荷轴向负荷总是作用在内圈，设置=6o=600N,径向载荷为M为100,200,300,400,500N,第j=1个钢球运行速度随时间的变化如11所示。径向载荷的添加将引起钢球运行速度c1的小幅波动。径向载荷的增加仅仅减少CDC1的数值，一般情况下，径向负载对轨道运动的影响并不显著。通过公式（16）-（17）和（19）-（20）,可以得到Auox,Aoy,和Auix,Aiy,他们为钢球与内外沟道之间椭圆接触区短长轴的相对滑动速度。为了准确地计算滑动速度，滑动速度的绝对值定义如下：(42)(43)0.3g0.36米。米Pasdari&GfintlH&am

33、p;results。PresentresulEfi0.30.280.25L0100200300400Axiall(>3dFj(N)500600图9.在F】n不断增加的情况下对cage的计算与实验结果进行比较(a)充520¥OM770(b)图10.转动分量(CDXl,coyl,CDZl)随时间和轴向载荷的变化在不同的径向载荷作用下，Aui,Auo随j=1的钢球轨道运动的变化分别如图12和图13所示。从图12的曲线中可以明显的看出Aui的波动,在轨道周期Top=(2九/c1)的某些区域,Au的值明显大于0,表明严重打滑发生在钢球与内圈接触的区域。一个轨道周期Top可以被划分为四个区

34、域：区域A(0-0.14)Top和(0.88-1)Top),B区域(0.14-0.42)Top,C区域(0.42-0.68)Top,D区域(0.68-0.88)Top,如图12所示，很明显，在B和D区域，即使加更高的径向载荷，滑动速度也很小。这表明，在这2个区域钢球与内滚道之间几乎不存在打滑。而在A和C区域，随着钢球的旋转，滑动速度先增加后减小，A区域滑动速度的最大值大于C区域的最大值。芭Pol)Lump柒ds=s7.6与11188881.40、FT=200-Fz'=fit)Ci>nin+Ft=30D-Fz=60Dm''nF*=400-=600nmF"5

35、0t>弁=60。>nin一1),Q门Oft=i>nin0.511,5Numberofballorbitalperiod图11.钢球的运行速度C1在径向载荷F：下随时间的变化NumberofballorbitalperiodAui图12.在组合荷载作用下，钢球相对于内圈的绝对速度随着径向载荷的增加，A和C区域的滑动速度也相应增大，表明在这两个区域打滑变得更严重。从图13中可以发现外圈的滑动也有此现象，不同的是，即使加更高的径向载荷D区域的滑动速度也不太低，然而，与A区域和C区域的滑动速度相比，它仍然是较低的值。在较高的径向载荷下，钢球与内外圈之间的滑动速度出现波动。只要对内圈施

36、加径向载荷其接触载荷分布就会发生显著变化，如图14所示。径向载荷的引入会导致B,D区域的接触力变大，远大于只受到纯轴向载荷时。径向载荷的增加使接触力明显上升，B和D区域被称为载荷增加区域，两载荷区域是由轴向和径向负载综合作用产生的。从图15可以看出径向载荷使B区域的接触力增大，使D区域的接触力减小。在轴向载荷作用下，内圈相对于X轴旋转，如图15所示。然后，钢球和沟道在D区域再接触时将会产生更大的接触力。不过，在相同的载荷条件下，D区域接触力的峰值低于B区域的峰值。在联合载荷条件下，存在两个接触力增加区（B和D）,同时也存在两个接触力减小区（A和C）。当施加径向载荷时，这2个区域的接触力减少，加

37、大径向载荷接触力降低的更明显。在联合载荷作用下，该系统滑动速度的波动，可以用接触力增加的区域来解释。较高的接触力，可提供足够的阻力，以保证钢球处于纯滚动状态，此时滑移速度很小，几乎不存在打滑。钢球从接触力增加的区域移动到接触力减小的区域时，接触力降低，从而没有足够的摩擦力和摩擦力矩来带动钢球运动。芭E与可合口虫学Bu-dd一«anosqvOFIOQ-FaffOOmm诺fFoD-F7=6M1n十Fr3OD-F2=8Mm'mPY=400-1=6MF：=500-哈8tMNumberofballorbitalperiod图13.在联合荷载作用下，钢球相对于外圈的绝对速度A0.511.

38、&Numberofballorbilalperiod图14.第j=1个钢球与内圈之间的接触力图15.轴承在联合荷载作用下接触力增加和降低区域示意图宏观的滑移发生在接触区，严重打滑也发生在该接触区。当钢球从接触力减小区运动到接触力增加区时，因接触力的增加其滑动速度将减小，打滑减轻，滑动速度保持在较低值。从上面的分析中可以得到，当轴向载荷in=600N时，轴承几乎不会出现打滑，而径向载荷的加入会使滑动速度产生波动，在接触力减少区更明显，当径向载荷增大时，可以发现更为明显的波动。在这些区域，钢球仍然承受载荷。较大的滑动速度会增加轴承和润滑油的温度，加剧轴承的磨损，缩短轴承的使用寿命。6.2.

39、 正弦曲线的径向载荷由于旋转机械的旋转，在轴承上的径向载荷可能随时变化。这里，认为随时间变化的径向载荷是正弦变化的，其频率等于内圈沟道的旋转速度。它可以表示为F1=Fr"+Nin(QX),其中Fr为恒定的径向负荷，r为正弦载荷的相对振幅。轴向载荷保持恒定值即1i=600N,Fr=400N,r=0.0,0.2,0.4时；滑动速度随钢球沟道速度的变化如图16所示。r=0.0表示径向载荷为常数，只要径向载荷为时变载荷(Nr=0.2,0.4),滑动速度随着钢球轨道运动的波动就是非周期性的，这是因为内圈旋转速度和钢球沿轨道速度之比不是整数。在载荷减小的区域，时变载荷使滑动速度增大了，然而，由于

40、非周期性波动，当钢球再次进入相同的载荷降低区域，时变载荷会使滑动速度减小。因此，很难说，时变载荷是否能提高或减轻钢球的打滑。但可以肯定的是，时变载荷可以增大最大滑移速度，如图16所示。此外，其对研究滑移速度随时间稳态变化中的频率分量有很大意义。利用了快速傅里叶变换，得到的滑移速度Aui的频谱如图17所示。其中恒定径向载荷(Nr=0.0),钢球的轨道速度为频率整倍数即icocl(i=1,2,3.)。当径向负荷为正弦时间载荷(r=0.2,0.4)时，icoc1的振幅会相应地增加。除了钢球的轨道速度和它的整数倍外，内圈转速和钢球轨道速度的合速度，即(i±ic)(i=0,1,2，)。本段考虑

41、了恒定的正弦径向负荷，分析其对钢球滑动行为的影响。虽然该动态模型已通过与Ref23保持架速度的实验结果的比较得到验证，但本节中研究得到的结果在之前未被发现。几个动态的研究已经进行了，重点是分析深沟球轴承的钢球与保持架之间的作用力24,干摩擦的摩擦力矩25,26,润滑剂的阻力27。他们没有研究钢球在联合载荷作用下其在三维空间中的运动情况，下边我们将对钢球在不同载荷条件下的三维运动进行实验研究。NumberoFballorbitalperiodO.8了.6,54,3,2JO0.0.ao.0,.0,0.(图E)Yi<占POI&ABuadlF55w£<16 .时变的径向负

42、荷对滑动17 .在时变径向载荷和恒定的轴向载荷作用下滑移速度图Aui的频谱。5.结论该三维非线性动态模型用来预测角接触球轴承的打滑行为，具综合考虑了钢球旋转引起的离心力和陀螺效应，钢球与内外圈的赫兹接触，钢球和保持架之间的不连续接触，弹流动力润滑等。通过与试验结果的比较，验证了动态模型可靠性。在此基础上，得出了在联合载荷作用下，钢球滑动速度随时间和空间的变化规律。结果表明径向载荷将使钢球与内外圈之间的滑移速度发生波动，在载荷降低区域具波动更明显。径向载荷的增加将显著增加滑动速度的幅度和范围，使打滑更加严重。钢球在接触载荷降低区承受载荷会使滑动速度变大，较大的滑动速度会大大增加轴承和润滑油的温度

43、，加剧轴承磨损，减少轴承的使用寿命。因此，在旋转机械的设计与检测中应认真考虑径向载荷的影响。本研究中的弹流润滑模型用于重载接触。当轴承在低的外载荷下运行时，钢球与内外沟道之间存在着严重的打滑。在这种情况下，该模型得到的油膜厚度和润滑油的动态粘度可能是不准确的。未来的研究应致力于通过实验来改善严重打滑的条件下的该模型。致谢论文中所进行的研究工作由国家摩擦学重点实验室资助。SKLT2013D02和北京自然科学基金资助项目no.3131002。参考文献1 T.A.Harris,Ananalyticalmodeltopredictskiddinginhighspeedrollerbearings,Tr

44、ibologyTransactions9(1966)229-241.2 T.A.Harris,Ananalyticalmodeltopredictskiddinginthrust-loadedangular-contactballbearings,JournalofLubricationTechnology(1971)17-24.3 N.T.Liao,J.F.Lin,Ballbearingskiddingunderradialandaxialloads,MechanismandMachineTheory37(1)(2002)91-113.4 F.Hirano,Motionofaballinan

45、gular-contactballbearing,ASLETransactions8(1965)425-434.5 N.T.Liao,J.F.Lin,Rolling-slidinganalysisinballbearingconsideringthermaleffect,TribologyTransactions49(1)(2006)1-16.6 X.Jiang,J.Ma,G.Fan,Minimumaxialpreloadforhigh-speedprecisionrollingbearings,Bearing6(2001)1-4.(inChinese).7 L.Cui,L.Wang,D.Zh

46、eng,Studyoncriticalloadofskiddinginhighspeedballbearing,JournalofAerospacePower22(11)(2007)1971-1976.(inChinese).8 S.Yuan,K.Guo,J.Hu,Investigationontheminimumaxialforcetoavoidslidingforhigh-speedballbearings,TransactionsofBeijingInstituteofTechnology31(9)(2011)1027-1031.(inChinese).9 W.Chen,Z.Ma,L.G

47、ao,X.Li,J.Pan,Quasi-staticanalysisofthrust-loadedangularcontactballbearingspartI:theoreticalformulation,ChineseJournalofMechanicalEngineering25(1)(2012)71-80.10 W.Chen,Z.Ma,L.Gao,X.Li,J.Pan,Quasi-staticanalysisofthrust-loadedangularcontactballbearingspartII:resultsanddiscussion,ChineseJournalofMecha

48、nicalEngineering25(1)(2012)81-87.11 T.Xu,G.Xu,Q.Zhang,C.Hua,H.Tan,S.Zhang,A.Luo,Apreloadanalyticalmethodforballbearingsutilisingbearingskiddingcriterion,TribologyInternational67(2013)4450.12 W.Chen,J.Li,L.Zhang,P.Wang,Skiddinganalysisofhighspeedrollingbearingconsideringwhirlingofbearing,JournalofMec

49、hanicalEngineering49(6)(2013)38-43.(inChinese).13 P.K.Gupta,AdvancedDynamicsofRollingElements,Springer-Verlag,NewYork,1984.14 J.O.Ostensen,H.Astrom,EHoglund,Analysisofagrease-lubricatedrollerbearingunderarcticconditions,ProceedingsoftheInstitutionofMechanicalEngineers,PartJ:JournalofEngineeringTribology209(1995)213-220.15 K.Imado,FundamentalstudyofHal