轴向受载高速角接触球轴承的动态特性求解方法

本资料包含pdf文件1个，下载需要1积分

Solving M ethod for Dynamic Performance of High Speed Angular ContactBal Bearings Loaded AxialyYu Xiao-kai ，Liu Kai-ge ，Zhao Chun-jiang ，Sun Bei-qi ，Xie Peng-fei(1.Luoyang Bearing Science＆Technology Co.，Ltd.，Luoyang 471039，China；2.Taiyuan University of Science and Technology，Taiyuan 030024，China)Abstract：A calculation method is introduced based O1"1 analysis of the mechanical characteristics of high speed angularcontact ball bearings only under axial load，which takes the actual contact angle as the basic variables to iterative solu-tion of dynamic performance of bearings.The method is compared with Harris method through an example，the resultsshow that the method iS effeetive and coincides with Harris result.The calculation parameters are reduced by half.re。

suhing in simple calculation and fast speed of convergence。

Key words：angular contact ball beating；axially load；contact angle；dynamic perform ance高速电动机等高速旋转机械主轴中常采用的高速角接触球轴承承受的径向力与轴向力相比很小，因此，轴承受力模型可以简化为纯轴向受载模型。高速角接触球轴承的动态特性如接触角、轴向变形量和刚度会影响高速旋转机械的性能；因收稿日期：2012-10-22；修回日期：2012-12-25此计算不同轴向力下高速角接触球轴承的实际接触角、轴向变形量和轴向刚度意义重大。其计算方法通常采用 Harris方法，但这种方法需要求解非常复杂的非线性方程，直接求解很困难，计算精度也很难达到要求。文中将实际接触角设定为迭代基本变量，分析高速角接触球轴承轴向受力工况下的动态方程组，将非线性方程组的参数由[11]徐悦，田爱梅，张振鹏，等.导弹垂直发射系统弹性多体动力学建模与仿真[J].兵工学报，2008，29(9)：1 083-1 087。

[12] Zakhariev E.Dynamics of Rigid Muhibody Systemswith Clearances in the JointsJ].Mechanics of Strnc-tures＆ Machines.1999，27：75-99。

[13] Raji NIHar，SADEGHI Farshid，RATEICK Richard G，et a1.A Numetical Model for Life Scatter in Roling E1-ement Bearings[J].Journal of Tribology，2008，130(1)：1-9。

[14]Harris T A.Rolling Bearing Analysis[M].4th ed.NewYork：John Wiley&Sons，2000。

[15]Karkainen Anti，HELFERT Marlene，AESCHLIMANNBeat.et a1.Dynamic Analysis of Rotor System with Mis·aliaigned Retainer Bearings[J].Journal of Tribology，2008，130(2)：1-10。

[16]王燕霜，邓四二，杨伯原.润滑油五参数流变模型的研究[J].摩擦学报，2007，127(5)：461-465。

(编辑：温朝杰)《轴承)2013.No.4原来的7个减少为2个，简化了计算。

1 求解方程的简化求解角接触球轴承的动态特性需要从接触变形人手，建立轴承的几何方程、变形方程、球平衡方程、套圈平衡方程，包括 6Z2(Z为球数)个方程，再建立包括离心力、陀螺力矩、接触角等在内的辅助方程组，应用 Newton-Raphson方法求解 。

由于纯轴向载荷下高速角接触球轴承只发生轴向变形，基于此建立变形方程。由于轴承承受纯轴向载荷时，各个球的受力情况相同，因此力平衡方程个数比较少，可对力平衡方程进行归-化处理。

仅轴向受载时，球心和沟道曲率中心位置变化关系如图1所示。根据承载前、后内、外沟道曲率中心距的径向分量不变，可得变形协调方程为[( -0.5) ]COS Ot。[( -0.5)D 6i]COS i(fi -1)D COS Oto， (1)式中： ， 分别为内、外沟曲率半径系数；D 为球直径，mm； ，6 为球与内、外圈间的弹性变形量，mm， ，ot 为球与内、外圈间的实际接触角，(。)； 。为原始接触角，(。)。

内圈沟道曲率中心图 1 球 中心和 沟道 曲率 中心位置因为陀螺力矩对实际内、外接触角的影响很小，所以以实际内、外接触角为迭代基本变量的求解方法中，可以忽略陀螺力矩的影响。故球所受到的载荷可以仅考虑与内、外圈间的接触载荷以及自身的离心力，如图2所示。

根据球的受力关系，可知Qi Q。 Fa， (2)式中：F 为轴承受到的轴向力，N；Z为球数；Qi ，Q 为球与内、外圈间接触载荷的轴向分量。

图 2 离心力作 用下球 的受力图球与内、外圈的接触载荷为 ， (3) 。 (4)球与内、外圈间接触载荷的径向分量为Qi Qicos i cot i， (5)Q Q。c0s cot 。。 (6)由于球径向受力平衡，则Qi F Q。 ， (7)式中：F 为球所受的离心力，N。

将(5)～(6)式代人(7)式，得力平衡方程为Fat (8)球轴承的接触载荷与接触变形间的关系为QiK1 ， (9)Q。K ， (10)式中：Ki，K 为球与内、外圈间的载荷 -变形常数，N/ram 。

将(3)～(4)式分别代入(9)～(10)式，可得6i赢 丁， ( )6 ：( )丁2。 12)( 6e 赢 )叮将(11)和(12)式代人(1)式，可得[( -。.5)。 ( Fa)丁COS Oe[ -0.5)Dw(志 )丁]c0 -1)D COS O/o。 (13)联立(8)式和(13)式即得实际接触角的求解方程组。通过求解方程组，可得高速角接触球轴承的于晓凯，等：轴向受载高速角接触球轴承的动态特性求解方法实际接触角 OLi， ，从而通过图1所示的球中心和沟道曲率中心位置关系得到轴承的轴向位移量 6为6 [( -0.5)D 。]sin Ot [( -0.5)D 6 ]sin i-(fi -1)·D sin o。 (14)2 求解方程的基本变量分析根据轴承的动态求解方程组，分析其中的基本变量 J。离心力方程为F 2.26 ×10 D。 D 3 2， (15)式中：D。 为球组节圆直径，mm；n 为球公转速度，r/min。球的公转速度 n 为(1- )COS(OL。-/3)±Oo(1 )COS(O/。-/s)” - (1 )COS(O/i-/s)(1-yi)cos( -/3) 。

(16)根据外沟道控制理论假设卢： rcta ( )， (17)COS Ote 。

Yi( )： ， (18) -- - L lpw式中：内、外圈旋转方向相反时劝-”，相同时取”； 为套圈角速度，rad/s；/3为球姿态角，(。)； 为几何参数。

接触变形常数 为K2.15×10 (∑P) 6 ， (19)式中： 为接触变形系数，可以根据主曲率差函数，(p)查表求得；∑P为主曲率和函数，mm～。

根据(15)～(19)式，依次寻找离心力、接触载荷变形常数与内、外接触角的关系，如图 3所示，图中箭头左边是箭头右边的函数。

图 3 变量传 递关 系3 求解步骤及算例纯轴向载荷下高速角接触球轴承上述方程组的求解流程如图4所示。求解过程中，采用全局收敛算法(Broyden算法)。

利用文中给出的计算方法对 7218B角接触球轴承在给定的受力和转速情况下进行编程计算 ，图4 计算流程 图得到表 1所示的计算结果。7218B角接触球轴承的几何参数为 D 22.23 mm，D 125.26 mm，fi 0.523 2，Z16，Odo40。。由表 1可知，采用文中的方法与Haris的方法得到的结果相差不大，充分说明了基于接触角迭代方法求解纯轴向受力的高速角接触球轴承动态特性的可行性 。

表 1 不同工况下的计算结果4 结束语依据纯轴向受力的高速角接触球轴承的特点，建立了轴承的变形协调方程和力平衡方程，基于实际接触角对方程组进行迭代求解，求解结果与 Haris方法的计算结果比较吻合，说明了该计算方法的有效性，同时，新计算方法中求解参数比Haris方法减少了-半，计算量较校