摘要
螺栓是转盘轴承的关键紧固件。为考察紧固螺栓工作参数对转盘轴承载荷分布的影响,提出了一种针对三排滚柱式转盘轴承的有限元等效建模方法。首先,基于三排滚柱式转盘轴承的结构特点和受力情况,建立考虑紧固螺栓的转盘轴承有限元模型,并从减少计算量和保证计算效率等方面考虑,采用非线性弹簧单元模拟滚子,通过有限元仿真快速获取了该转盘轴承的载荷分布。然后,设计并开展三排滚柱式转盘轴承静态加载试验,通过有限元仿真结果和试验结果的对比验证了所构建有限元模型的准确性。最后,基于构建的有限元模型,分析了螺栓预紧力和螺栓结合面处摩擦系数对三排滚柱式转盘轴承载荷分布的影响。仿真结果表明:螺栓预紧力对三排滚柱式转盘轴承载荷分布的影响显著,而螺栓结合面处摩擦系数的影响微弱,其中上排滚子承受的最大载荷随螺栓预紧力的增大而增大,下排滚子承受的最大载荷随预紧力的增大而减小;在轴向力和倾覆力矩的联合作用下,三排滚柱式转盘轴承下排滚子承受的载荷总体上高于上排滚子,充足的螺栓预紧力可以提高转盘轴承的承载能力,而螺栓预紧力不足则会使下排滚子承受更大的载荷,导致转盘轴承提前失效。研究结果可为转盘轴承的安装和后期维护提供参考。
转盘轴承是大型旋转机械的核心部件,被广泛应用于镗床、塔式起重机、风力发电机、挖掘机和盾构机等设
为了提高转盘轴承的承载能力,确定其滚道上的载荷分布是关键。目前,国内外许多学者针对转盘轴承的载荷分布开展了相关研究。何培瑜
三排滚柱式转盘轴承的结构复杂,其滚道上的载荷分布容易受到多种因素的影响。螺栓是转盘轴承的重要紧固件,但其对转盘轴承载荷分布的影响无法直接通过理论计算得到。基于此,笔者拟采用非线性弹簧单元模拟实体滚子,通过建立三排滚柱式转盘轴承的有限元模型来分析螺栓预紧力及螺栓结合面处摩擦系数对转盘轴承载荷分布的影响,旨在为转盘轴承的实际工程应用提供参考。
以130.20.1005型三排滚柱式转盘轴承为研究对象,其结构参数如
三排滚柱式转盘轴承的结构如

图1 三排滚柱式转盘轴承结构示意
Fig. 1 Schematic diagram of three-row roller slewing bearing structure
在轴向力Fa的单独作用下,三排滚柱式转盘轴承内的载荷均匀分布,如
(1) |
(2) |
式中:Qmax 1、Qmax 2分别为上、下排滚子承受的最大载荷,N;QM max 为倾覆力矩作用下滚子承受的最大载荷,N;Qn为轴向力作用下滚子承受的载荷,N;N为滚子数量;D为滚道节圆直径,mm。

图2 轴向力作用下三排滚柱式转盘轴承的载荷分布
Fig. 2 Load distribution of three-row roller slewing bearing under the axial force

图3 轴向力和倾覆力矩联合作用下三排滚柱式转盘轴承的载荷分布
Fig. 3 Load distribution of three-row roller slewing bearing under the combined action of axial force and tilting moment
三排滚柱式转盘轴承中含有大量滚子-滚道接触对,且滚子与滚道间的接触情况复杂。此外,分析非线性接触问题时还可能出现计算不收敛或计算耗时过长等情况。因此,须对滚子-滚道接触模型进行简化,从而提高计算效率。在本文分析中,采用非线性弹簧单元来模拟滚子。通过将滚子与滚道在载荷作用下的接触变形关系赋予非线性弹簧单元,以使弹簧单元可用于表征滚子与滚道之间的接触特性。根据赫兹接触理论,滚子与滚道在载荷作用下的接触变形关系可表示
(3) |
式中:Q为载荷,N;L为滚子长度,mm;δ为滚子与滚道的总变形量,mm。
将滚子与滚道在载荷作用下的接触变形关系简化为:
(4) |
式中:k为接触刚度,N/mm。
根据

图4 三排滚柱式转盘轴承滚子的变形量—载荷曲线
Fig. 4 Curve of deformation‒load of roller of three-row roller slewing bearing

图5 三排滚柱式转盘轴承滚子-滚道接触模型
Fig. 5 Contact model of roller-raceway of three-row roller slewing bearing
考虑到三排滚柱式转盘轴承为对称结构,为了提高计算效率,只构建一半的转盘轴承有限元模型,其中1个实体滚子简化成1对非线性弹簧单元。在后续分析中,为便于观察各个滚子承受的载荷,按顺时针方向对非线性弹簧单元进行编号,如

图6 三排滚柱式转盘轴承有限元模型
Fig. 6 Finite element model of three-row roller slewing bearing
建模时考虑三排滚柱式转盘轴承所受的轴向力和倾覆力矩。为了便于施加载荷,在转盘轴承外圈的几何中心设置主节点RP1,并将RP1与外圈滚道上表面进行耦合。由于只构建了一半的三排滚柱式转盘轴承有限元模型,则施加在模型上的载荷为:轴向力Fa=59.5 kN,倾覆力矩M=4.155 kN·m。在设置边界条件时,将安装基础的下表面定义为完全约束,释放主节点RP1处倾覆力矩方向上的旋转自由度和轴向自由度。
在有限元建模时,划分网格的单元类型和疏密程度会极大地影响计算精度。参考文献[

图7 三排滚柱式转盘轴承滚道区域网格划分
Fig. 7 Mesh generation of raceway area of three-row roller slewing bearing finite element model

图8 螺栓预紧力的施加示意
Fig. 8 Schematic diagram of application of bolt preload
为了验证中间排滚子对上、下排滚子载荷分布的影响,分别建立了有限元模型1和有限元模型2。其中,模型1忽略了中间排滚子,模型2考虑了中间排滚子。基于2个有限元模型,分析得到三排滚柱式转盘轴承的载荷分布曲线,如

图9 三排滚柱式转盘轴承的载荷分布
Fig. 9 Load distribution of three-row roller slewing bearing
为验证上述有限元模型的准确性,对比由仿真分析和经验
为了进一步验证有限元模型的准确性,对三排滚柱式转盘轴承开展静态加载试验。转盘轴承静态加载试验台如

图10 转盘轴承静态加载试验台
Fig. 10 Static loading test bench for slewing bearing

图11 三排滚柱式转盘轴承内侧应变片布置
Fig. 11 Arrangement of strain gauges on the inner side of three-row roller slewing bearing
通过试验测量,获得了三排滚柱式转盘轴承上、下排滚子内圈周向上的应力分布,如

图12 三排滚柱式转盘轴承应力分布
Fig. 12 Stress distribution of three-row roller slewing bearing
为了便于分析仿真结果与试验结果的误差,采用均方根误差(root mean square error,RMSE
(5) |
式中:为应力仿真值,MPa;为应力试验值,MPa;i为测点编号;n为测点数量,本文n=18 个。
将相应的数据代入
转盘轴承通常以螺栓连接的方式固定在支承结构上,紧固螺栓的预紧力及其结合面处的摩擦系数会影响转盘轴承的整体刚度,从而影响其载荷分布。为此,基于上文构建的有限元模型,分别考察螺栓预紧力及其结合面处摩擦系数对三排滚柱式转盘轴承载荷分布的影响,旨在为转盘轴承的工程应用提供参考和借鉴。
130.20.1005型三排滚柱式转盘轴承的外圈为动圈,内圈为固定圈,内圈通过螺栓固定在安装基础上。为保证转盘轴承装配体结构的刚性和连接紧固性,须对紧固螺栓施加预紧力。但由于预紧力施加方法不同,使得实际预紧力大小与目标加载值间存在误差,从而影响转盘轴承的载荷分布。为了分析螺栓预紧力对三排滚柱式转盘轴承载荷分布的影响,在其有限元模型中对螺栓施加不同的预紧力,分析不同预紧力下该转盘轴承的载荷分布情况。
分别对三排滚柱式转盘轴承有限元模型中的螺栓施加0,100和170 kN预紧力,通过有限元分析得到该转盘轴承的载荷分布,如

图13 不同螺栓预紧力下三排滚柱式转盘轴承的载荷分布
Fig. 13 Load distribution of three-row roller slewing bearing under different bolt preloads

图14 不同螺栓预紧力下三排滚柱式转盘轴承的最大载荷
Fig. 14 Maximum load of three-row roller slewing bearing under different bolt preloads
紧固螺栓与转盘轴承结合面处的摩擦系数受制造工艺、结合面处清洁度和氧化物等因素的影响。基于上文构建的有限元模型,考察螺栓结合面处不同摩擦系数对转盘轴承载荷分布的影响。在三排滚柱式转盘轴承有限元模型中,设置不同的螺栓结合面处摩擦系数(0.05,0.10和0.15),通过有限元分析得到该转盘轴承的载荷分布,如

图15 螺栓结合面处不同摩擦系数下三排滚柱式转盘轴承的载荷分布
Fig. 15 Load distribution of three-row roller slewing bearing under different friction coefficients of bolt joint surface
从

图16 螺栓结合面处不同摩擦系数下三排滚柱式转盘轴承的最大载荷
Fig. 16 Maximum load of three-row roller slewing bearing under different friction coefficients of bolt joint surface
通过有限元方法,建立了考虑紧固螺栓的三排滚柱式转盘轴承有限元模型。建模时通过非线性弹簧单元的等效处理,减小了运算规模,使得在普通计算机上也能实现对转盘轴承载荷分布的求解。同时,通过转盘轴承的静态加载试验验证了有限元模型的准确性。
结果表明:在紧固螺栓的工作参数中,螺栓预紧力对三排滚柱式转盘轴承载荷分布的影响显著,而螺栓结合面处摩擦系数的影响微弱;当螺栓预紧力增大时,上排滚子承受的最大载荷随之增大,而下排滚子承受的最大载荷显著减小;在轴向力和倾覆力矩的联合作用下,下排滚子承受的载荷总体上高于上排滚子,则增大螺栓预紧力有利于提高三排滚柱式转盘轴承的承载能力。在转盘轴承的安装和后期维护中,应当保证充足的螺栓预紧力,避免因螺栓预紧力不足而导致轴承提前失效。
参考文献
刘荣,李玉川,刘静,等.单排四点接触球转盘轴承承载能力研究[J].轴承,2018(12):34-38. doi:10.19533/j.issn1000-3762.2018.12.009 [百度学术]
LIU Rong, LI Yu-chuan, LIU Jing, et al. Study on load capacity of single row four point contact ball slewing bearings[J]. Bearing, 2018(12): 34-38. [百度学术]
何培瑜,洪荣晶,王华.转盘轴承对数修形滚子的特性分析[J].南京工业大学学报(自然科学版),2018,40(6):70-78. doi:10.3969/j.issn.1671-7627.2018.06.011 [百度学术]
HE Pei-yu, HONG Rong-jing, WANG Hua. Movement characteristics analysis of logarithmic modification roller slewing bearings[J]. Journal of Nanjing Tech University(Natural Science Edition), 2018, 40(6): 70-78. [百度学术]
张占立,周鹏举,左朋辉,等.负游隙对四点接触转盘球轴承摩擦力矩影响研究[J].机床与液压,2018,46(1):56-58. doi:10.3969/j.issn.1001-3881.2018.01.012 [百度学术]
ZHANG Zhan-li, ZHOU Peng-ju, ZUO Peng-hui, et al. Study of negative axial clearance affected friction torque of four contact-point slewing ball bearing[J]. Machine Tool & Hydraulics, 2018, 46(1): 56-58. [百度学术]
王燕霜,袁倩倩.负游隙对特大型双排四点接触球轴承载荷分布的影响[J].机械工程学报,2012,48(21):110-115. doi:10.3901/JME.2012.21.110 [百度学术]
WANG Yan-shuang, YUAN Qian-qian. Influence of negative clearance on load distributions of large-size double row four-point-contact ball bearings[J]. Journal of Mechanical Engineering, 2012, 48(21): 110-115. [百度学术]
WANG Yan-shuang, CAO Jia-wei. Determination of the precise static load-carrying of pitch bearings based on static models considering clearance[J]. International Journal of Mechanical Sciences, 2015, 100: 209-215. doi:10.1016/j.ijmecsci.2015.06.023 [百度学术]
赵春江,熊杰,于晓凯,等.惯性载荷对高速角接触球轴承最优预紧力的影响分析[J].工程设计学报,2015,22(2):150-154. doi:10.3785/j.issn.1006-754X.2015.02.008 [百度学术]
ZHAO Chun-jiang, XIONG Jie, YU Xiao-kai, et al. Analysis on the influence of inertial load on the optimal preload of high-speed angular contact ball bearing[J]. Chinese Journal of Engineering Design, 2015, 22(2): 150-154. [百度学术]
HERAS I, AGUIRREBEITIA J, ABASOLO M, et al. Load distribution and friction torque in four-point contact slewing bearings considering manufacturing errors and ring flexibility[J]. Mechanism & Machine Theory, 2019, 137: 23-36. doi:10.1016/j.mechmachtheory.2019.03.008 [百度学术]
GAO X, HUANG X, WANG H, et al. Effect of raceway geometry parameters on the carrying capability and the service life of a four-point-contact slewing bearing[J]. Journal of Mechanical Science and Technology, 2010, 24(10): 2083-2089. doi:10.1007/s12206-010-0709-5 [百度学术]
HE P, LIU R, HONG R J, et al. Hardened raceway calculation analysis of a three-row roller slewing bearing[J]. International Journal of Mechanical Sciences, 2018, 137: 133-144. doi:10.1016/j.ijmecsci.2018.01.021 [百度学术]
PETER G, ROK P, SRECKO G. Computational model for determination of static load capacity of three-row roller slewing bearings with arbitrary clearances and predefined raceway deformations[J]. International Journal of Mechanical Sciences, 2013, 73: 82-89. doi:10.1016/j.ijmecsci.2013.04.012 [百度学术]
LUDWIK K, RAFAL P. Modification of the raceway profile of a single-row ball slewing bearing[J]. Mechanism & Machine Theory, 2018, 128: 1-15. doi:10.1016/j.mechmachtheory.2018.05.009 [百度学术]
AGUIRREBEITIA J, ABASOLO M, RAFAL A, et al. General static load-carrying capacity for the design and selection of four contact point slewing bearings: finite element calculations and theoretical model validation[J]. Finite Elements in Analysis & Design, 2012, 55: 23-30. doi:10.1016/j.finel.2012.02.002 [百度学术]
DAIDIE A, CHAIB Z, GHOSN A. 3D simplified finite elements analysis of load and contact angle in a slewing ball bearing[J]. Journal of Mechanical Design, 2008, 130(8): 082601. doi:10.1115/1.2918915 [百度学术]
王存珠,陈观慈,李肖杰,等.螺栓及其预紧力对大型转盘轴承接触载荷分布的影响[J].轴承,2016(1):24-27,33. doi:10.3969/j.issn.1000-3762.2016.01.009 [百度学术]
WANG Cun-zhu, CHEN Guan-ci, LI Xiao-jie, et al. Effects of bolts and its pre-tightening force on contact load distribution in large-scale slewing bearings[J]. Bearing, 2016(1): 24-27, 33. [百度学术]
王永全,王华,高学海,等.螺栓预紧力对转盘轴承载荷分布的影响[J].轴承,2013(12):13-16,50. doi:10.3969/j.issn.1000-3762.2013.12.004 [百度学术]
WANG Yong-quan, WANG Hua, GAO Xue-hai, et al. Influence of bolt pretension on load distribution for slewing bearings[J]. Bearing, 2013(12): 13-16, 50. [百度学术]
WANG H, HE P, PANG B, et al. A new computational model of large three-row roller slewing bearings using nonlinear springs[J]. Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, 2017, 231: 203-210. doi:10.1177/0954406217704223 [百度学术]
姬丽丽,王华,潘裕斌,等.基于非线性弹簧的三排圆柱滚子组合转盘轴承静态承载能力模型[J].轴承,2016(2):1-5. doi:10.3969/j.issn.1000-3762.2016.02.001 [百度学术]
JI Li-li, WANG Hua, PAN Yu-bin, et al. Static load capacity model of three-row roller cylindrical slewing bearing based on nonlinear spring[J]. Bearing, 2016(2): 1-5. [百度学术]
石亦平,周玉蓉.ABAQUS有限元分析实例讲解[M].北京:机械工业出版社,2006:42-51. [百度学术]
SHI Yi-ping, ZHOU Yu-rong. ABAQUS finite element analysis examples[M]. Beijing: China Machine Press, 2006: 42-51. [百度学术]
国家质量监督检验检疫总局, 中国国家标准化管理委会.六角头螺栓: GB/T5782—2000[S]. 北京:中国标准出版社,2000:1-8. [百度学术]
State General Administration of Quality Supervision, Inspection and Quarantine, Standardization Administration of China. Hexagonal head bolts: GB/T5782‒2000[S]. Beijing: Standards Press of China, 2000:1-8. [百度学术]
Association of German Engineers. Systematic calculation of highly stressed bolted joints (joints with one cylindrical bolt): VDI 2230[S]. Berlin: Beuth Press, 2003: 113-123. [百度学术]
LIU R, WANG H, PANG B. Load distribution calculation of a four-point-contact slewing bearing and its experimental verification[J]. Experimental Techniques, 2018, 42: 243-252. doi:10.1007/s40799-018-0237-2 [百度学术]