技术领域
[0001] 本
发明属于轴承设计领域,具体涉及一种三点接触球轴承异常接触擦伤定量评估方法。
背景技术
[0002] 三点接触球轴承是一类特殊的
角接触轴承,是一种可承受以双向轴向
载荷为主、少量径 向载荷及倾覆
力矩的单列球轴承,其优点在于可承受联合载荷,高速性能好;由于其
内圈可 分所以拆装方便,并且可设计
保持架为一体式;工作过程
中轴向力变向时窜动小,所以广泛 地应用于航空
发动机、
火箭发动机、舰船柴油机以及各类高速
齿轮箱
传动轴中。所谓三点接 触轴承是指其在零游隙设计状态下轴承接触点为三点,但是一般来说其正常工作状态下,均 为两点接触。不过由于几何参数设计及工作状态变化会出现在非主承力内半圈
沟道上产生第 三点接触,该接触点为异常接触会导致滚球和内圈沟道的严重滑动,出现严重擦伤故障,使 得
钢球和沟道
精度下降,改变轴承内部的稳定运动状态,严重时会导致轴承过早损坏。所以 针对三点接触球轴承分析其工作状态下产生三点接触的条件及产生后擦伤程度评估工作对于 轴承设计、装配、运行等工作环节均具有重要的指导意义。
[0003] 为了研究三点接触球轴承的异常接触状态及擦伤定量评估方法,需要分析异常接触影响 因素并最终建立合理的轴承分析模型及定量描述擦伤磨损的评定标准。轴承分析模型方面, 徐洪江等(徐洪江,徐卫军.三点接触球轴承设计改进浅析[J].哈尔滨轴承,2007,28(3):6-7.) 根据几何关系,在设计中增加了避免多点接触的游隙验算,对于某型航空发动机轴承进行了 改进设计,提升其所涉及轴承的性能和使用寿命。李杰等(李杰,田拥胜,张华良,et al.考虑 轴向力影响的三点接触球轴承
刚度特性研究[J].推进技术,
2018.)运用商用有限元
软件建立了 三点接触球轴承模型分析在引入了异常三点接触刚度变化。Leblanc A等(Leblanc A,Nelias D, Plona D.Ball Motion and Sliding Friction in an Arched Ball Bearing[C]//Rolling Element Bearings.2007.)应用轴承滚球的运动学分析代替了
滚道控制理论建立了四点接触轴承的拟静 力学模型并对不同工况下轴承接触点内部载荷及相对滑移进行了分析。Halpin J D等(Halpin J D,Tran AN.An Analytical Model of Four-Point Contact Rolling Element Ball Bearings[J]. Journal of Tribology,2016,138(3):031404.)应用接触点内摩擦耗能最小准则计算了轴承滚球自 转的方位角替代了滚道控制理论建立了四点接触轴承拟静力学模型并推导了轴承刚度计算公 式并提出改模型可进一步改写为三点接触球轴承模型,最后分析了轴承接触情况载荷分布及 轴承寿命。Olofsson等(OLOFSSON,ANDERSSON,
Simulation of mild wear in boundary lubricated spherical roller thrust bearings[J].Wear,2000,241(2):180-185.)建立了一个迭 代磨损计算模型,利用阿查得磨损公式对接触时的磨损进行了仿真,实现了磨损后接触表面 几何形貌的估计并定性地进行了实验验证。Feng等(H Liu,C&Y Chen,X&M Gu,J&N Jiang, S&L Feng,Z.(2009).High-speed wear lifetime analysis of instrument ball bearings.Proceedings of The Institution of Mechanical Engineers Part J-journal of Engineering Tribology-PROC INST MECH ENG J-J ENG TR.223.497-510.10.1243/13506501JET496.)利用阿差得磨损公式建立了 球轴承的磨损寿命模型。
[0004] 从文献调研中可以发现,针对三点接触球轴承的建模分析中包括几何分析、有限元计算 分析以及拟静力学分析,上述模型均有对应的假设不符合轴承结构和实际工况。需要考虑轴 承每个部件的三维运动、滑动
牵引力和保持架效应,并进行动力学仿真,根据仿真结果分析 其三点接触状态并根据金属材料磨损原理引出擦伤程度评估参数具有重要工程意义及实际需 求。
发明内容
[0005] 本发明的目的在于针对上述
现有技术的不足,提供了一种三点接触球轴承异常接触擦伤 定量评估方法。
[0006] 本发明采用以下技术方案来实现的:
[0007] 一种三点接触球轴承异常接触擦伤定量评估方法,包括以下步骤:
[0008] 1)三点接触球轴承几何参数计算
[0009] 根据轴承球径和内
外圈沟道
曲率半径建立
垫片角/垫片宽度、初始接触角、
轴向游隙/径向 游隙的几何关系;
[0010] 2)三点接触球轴承异常接触状态判别
[0011] 基于步骤1)中的几何关系,由游隙和垫片角计算初始接触角,并以初始接触角小于等 于垫片角为产生异常接触判别准则获得异常接触几何边界条件,即游隙和垫片角参数化关系, 并基于
轴承外圈、双半内圈沟道曲率中心及球心几何
位置关系判断球轴承接触状态;
[0012] 3)基于数值模型的轴承运动仿真
[0013] 考虑轴承主要运动部件及
润滑剂的相互作用,建立轴承部件运动微分方程,并采用数值 积分方法进行时间积分运算,其中轴承主要运动部件包括左半内圈、右半内圈、外圈和保持 架;
[0014] 4)异常接触点内相对滑移速度场及压力场计算
[0015] 根据步骤3)计算所得轴承滚球和套圈运动,考虑滚球公转、自转、陀螺运动、自旋运 动及套圈的绕轴旋转,计算在
轴承套圈与滚球接触椭圆内相对滑动速度并计算滑移速度场, 根据赫兹接触理论计算接触椭圆内压力分布;
[0016] 5)异常接触点内PV值计算及异常接触擦伤
风险分析
[0017] 由阿查得磨损公式引出以PV值定量描述擦伤磨损,根据步骤4)中所得滚球相对轴承内 圈运动速度计算异常接触点内的相对滑动速度分布及压
应力分布,计算接触点内PV值的极 值,对轴承产生擦伤的风险进行定量评估,其中PV值为相对滑动速度和压应力的乘积。
[0018] 本发明进一步的改进在于,步骤2)中三点接触球轴承异常接触状态判别,具体方法如 下:
[0019] 根据步骤1)中分析三点接触球轴承异常接触的边界条件为垫片角αs等于初始接触角α, 其计算公式为:
[0020]
[0021] 其中s为垫片宽度αs中下标s表示垫片,即角度αs为垫片宽度量s相关角度,如图2、3 所示,Ri为内圈沟道
曲率半径,D为轴承滚球直径;
[0022]
[0023] 其中Z=Sd/D为
径向游隙系数,Sd为径向游隙,B=fi+fo-1,fi为内圈沟道曲率 半径系数,fo为外圈沟道曲率半径系数,其中下标i表示内圈,o表示外圈,Sd中下标d表示 径向方向;
[0024] 产生异常接触的条件为:
[0025] α≤αs (3)
[0026] 进一步考虑轴承受载条件对三点异常接触进一步进行判断:
[0027] 规定轴承左右半内圈固结,轴向力为正时左半内圈为主承力内圈,反之亦然,则主承力 内圈曲率中心与非主承力内圈曲率中心关系如下:
[0028]
[0029] 其中 为主承力内圈曲率中心坐标, 为非主承力内圈曲率中心坐标,FX为轴承内 圈轴向受力;
[0030] 轴承非承力内圈产生接触的条件为:
[0031]
[0032] 其中 分别为球心、主承力内半圈沟道曲率中心、非主承力内半圈沟道曲 率中心在内圈固结坐标下矢量, 为非主承力内半圈接触点在内圈固结坐标下矢量,其值为:
[0033]
[0034] 若满足上述式(5)条件则产生非承力半内圈产生三点异常接触,其接触力求解先求得接 触点趋近量并根据赫兹接触理论进行计算;
[0035] 上述式(4)至(6)中上标ir表示该向量为内圈固结
坐标系下向量,下标b表示滚球即向量 表示滚球球心在内圈固结坐标系下矢量,下标irm表示主承力内圈,下标ira表示非主承力内 圈,即 表示主承力内圈滚道曲率中心在内圈固结坐标系下矢量,同理 下标ica表示非主 承力内圈与滚球接触点,即向量 表示非主承力内圈接触点在内圈固结坐标中是啊零两位 置。
[0036] 本发明进一步的改进在于,步骤3)中基于数值模型的轴承运动仿真,如下:
[0037] 运用步骤2)中所建立的三点接触球轴承异常接触条件,考虑轴承受载、赫兹接触力、
离心力、陀螺力矩、拖动力及阻力具体工作条件下受力条件进一步分析轴承运动状态,为接 触点内相对滑移速度计算提供
基础,其中轴承受载包括轴向力、径向力和倾覆力矩。
[0038] 本发明进一步的改进在于,步骤4)中异常接触点内相对滑移速度场计算,如下:
[0039] 运用步骤3)中考虑实际工况中受力条件的所的轴承滚球和套圈的运动,考虑轴承滚球 公转、自转、自旋、陀螺运动及
轴承内圈绕轴运动,计算在轴承套圈与滚球接触椭圆内相对 滑动速度并计算滑移速度场。
[0040] 本发明进一步的改进在于,步骤5)中异常接触点内PV值计算及异常接触擦伤风险分析, 如下:
[0041] 首先由阿查得磨损公式:
[0042]
[0043] 其中Q为磨损量、K为微凸体产生磨粒概率、H为材料布氏硬度值、W为法向载荷、L 为相对滑动距离,由此引出以PV值定量描述擦伤磨损程度,根据步骤4)中所得相对滑动速 度分布及压应力分布,计算接触点内PV值的极值,对轴承产生擦伤的风险进行定量评估。
[0044] 本发明具有如下有益的技术效果:
[0045] 本发明提供的一种三点接触球轴承异常接触擦伤定量评估方法,该方法考虑轴承每个部 件的三维运动、滑动牵引力和保持架效应建立三点接触球轴承模型,对异常三点接触及擦伤 状况进行分析,获得轴承在实际运行工况下异常三点接触情况,并根据金属材料磨损理论引 入PV值定量评估轴承擦伤严重程度。本发明为三点接触球轴承不同工况下异常接触擦伤程 度评估提供现实可靠的参考准则。
附图说明
[0046] 图1是本发明中方法技术路线图。
[0047] 图2是本发明中三点接触球轴承垫片宽度及径向游隙示意图。
[0048] 图3是本发明中三点接触球轴承初始接触角及垫片角示意图,其中,图3(a)为三点接 触球轴承初始接触角示意图,图3(b)为三点接触球垫片角示意图。
[0049] 图4是本发明中案例轴承垫片角、游隙系数及初始接触角关系。
[0050] 图5是本发明中三点接触球轴承套圈与滚球接触关系。
[0051] 图6是本发明中案例轴承套圈滚球接触点内相对滑移速度分布,其中图6(a)中轴向力 900N,图6(b)中轴向力1500N。
[0052] 图7是本发明中案例轴承套圈滚球接触点内PV值分布。
具体实施方式
[0053] 以下结合附图和
实施例对本发明提供的一种三点接触球轴承异常接触擦伤定量评估方法 作进一步详细说明。整体
流程图如图1所示。
[0054] 本发明实施实例采用如表1所示参数的球轴承。
[0055] 具体实施步骤如下所示:
[0056] 表1轴承参数
[0057]
[0058] (1)轴承几何参数分析。根据轴承球径、内外圈沟道曲率半径等几何参数建立垫片角/ 垫片宽度、初始接触角、轴向游隙/径向游隙及的几何关系;
[0059] 分析三点接触球轴承异常接触的边界条件为垫片角αs等于初始接触角α,其计算公式 为:
[0060]
[0061] 其中s为垫片角宽度,Ri为内圈沟道曲率半径,D为轴承滚球直径。
[0062]
[0063] 其中Z=Sd/D为径向游隙系数,Sd为径向游隙,B=fi+fo-1,fi为内圈沟道曲率 半径系数,fo为外圈沟道曲率半径系数。
[0064] (2)三点接触球轴承异常接触状态判别。基于步骤(1)中推导关系,由游隙和垫片角 计算初始接触角,并以初始接触角小于等于垫片角为产生异常接触判别准则获得异常接触几 何边界条件,即游隙和垫片角参数化关系,并基于轴承外圈、双半内圈沟道曲率中心及球心 几何位置关系判断球轴承接触状态;
[0065] 产生异常接触的条件为
[0066] α≤αs (3)
[0067] 所得参数关系如图4所示
[0068] 规定轴承左右半内圈固结,轴向力为正时左半内圈为主承力内圈,反之亦然,则主承力 内圈曲率中心与非主承力内圈曲率中心关系如下:
[0069]
[0070] 其中 为主承力内圈曲率中心坐标, 为非主承力内圈曲率中心坐标,FX为轴承内 圈轴向受力。
[0071] 轴承非承力内圈产生接触的条件为
[0072]
[0073] 其中 分别为球心、主承力内半圈沟道曲率中心、非主承力内半圈沟道曲 率中心在内圈固结坐标下矢量, 为非主承力内半圈接触点在内圈固结坐标下矢量,其值为:
[0074]
[0075] 若满足上述式(5)条件则产生非承力半内圈产生三点异常接触,其接触力求解先求得接 触点趋近量并根据赫兹接触理论进行计算。非承力内半圈产生接触条件几何关系示意图如图 6,趋近量小于零对应于式(5)中上式,真实的接触点应该位于内圈沟道上即在图X中内半 圈沟道所在弧线实线处,对应代数关系为式(5)中下式。
[0076] 异常接触几何边界条件如图5所示,本案例选择初始接触角为20°,计算所得垫片宽度 为0.26mm,为讨论制造误差产生多点接触及擦伤风险评估,选择加工误差为0.01mm即垫片 宽度为0.27mm。
[0077] (3)基于数值模型的轴承运动仿真。运用步骤(2)中所建立的三点接触球轴承异常接 触条件,考虑轴承受载(轴向力、径向力、倾覆力矩)、赫兹接触力、离心力、陀螺力矩、拖 动力、阻力等具体工作条件下受力条件进一步分析轴承运动状态,为接触点内相对滑移速度 计算提供基础。
[0078] (4)异常接触点内相对滑移速度场计算。运用步骤(3)中考虑实际工况中受力条件的 所的轴承滚球和套圈的运动,考虑轴承滚球公转、自转、自旋、陀螺运动及轴承内圈绕轴运 动。计算在轴承套圈与滚球接触椭圆内相对滑动速度并计算滑移速度场,相对滑动速度分布 如图6所示。
[0079] (5)异常接触点内PV值计算及异常接触擦伤风险评估。由阿查得磨损公式引出以PV 值定量描述擦伤磨损程度,根据步骤(4)中所得相对滑动速度分布及压应力分布,计算接触 点内PV(相对滑动速度和压应力乘积)的极值,对轴承产生擦伤的风险进行定量评估。
[0080] 对于步骤(2)所述几何参数及装配参数轴承在轴向力900N及1300N、转速12000rpm 工况下运用步骤(3)和(4)中仿真分析方法进行计算,由图7可得异常接触时接触点内速 度场,并计算可得轴向力为100N-2000N时PV极值如线标所示,轴向力为100-400N及 1700-2000N时不产生三点接触,轴向力为900N时PV极值取得最大值,磨损程度最高,实 际工况应避免转速为12000rpm时轴向力降低至900N区间。
[0081] 表2异常接触PV值随轴向力变化
[0082]
[0083]