CN109550979A - 一种高速主轴单元轴承套圈轴向热位移确定方法 - Google Patents
一种高速主轴单元轴承套圈轴向热位移确定方法 Download PDFInfo
- Publication number
- CN109550979A CN109550979A CN201811284426.9A CN201811284426A CN109550979A CN 109550979 A CN109550979 A CN 109550979A CN 201811284426 A CN201811284426 A CN 201811284426A CN 109550979 A CN109550979 A CN 109550979A
- Authority
- CN
- China
- Prior art keywords
- bearing
- sphere
- thermal
- load
- thermal walking
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
- B23Q—DETAILS, COMPONENTS, OR ACCESSORIES FOR MACHINE TOOLS, e.g. ARRANGEMENTS FOR COPYING OR CONTROLLING; MACHINE TOOLS IN GENERAL CHARACTERISED BY THE CONSTRUCTION OF PARTICULAR DETAILS OR COMPONENTS; COMBINATIONS OR ASSOCIATIONS OF METAL-WORKING MACHINES, NOT DIRECTED TO A PARTICULAR RESULT
- B23Q1/00—Members which are comprised in the general build-up of a form of machine, particularly relatively large fixed members
- B23Q1/70—Stationary or movable members for carrying working-spindles for attachment of tools or work
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
- B23Q—DETAILS, COMPONENTS, OR ACCESSORIES FOR MACHINE TOOLS, e.g. ARRANGEMENTS FOR COPYING OR CONTROLLING; MACHINE TOOLS IN GENERAL CHARACTERISED BY THE CONSTRUCTION OF PARTICULAR DETAILS OR COMPONENTS; COMBINATIONS OR ASSOCIATIONS OF METAL-WORKING MACHINES, NOT DIRECTED TO A PARTICULAR RESULT
- B23Q17/00—Arrangements for observing, indicating or measuring on machine tools
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Rolling Contact Bearings (AREA)
Abstract
本发明公开了一种高速主轴单元轴承套圈轴向热位移确定方法,包括步骤一、建立高速下的角接触球轴承内部几何相容关系,得到位置角ψj处的轴承内外圈滚道最终的轴向和径向相对距离;步骤二、得到角接触球轴承沟道接触的变形几何相容方程;步骤三、求出稳定状态时轴承的各向位移值,建立轴承所受外载荷和位移之间的关系,轴承受到联合载荷的作用,即受轴向载荷Fa、径向载荷Fr和力矩载荷M的综合作用;得到对应的轴承内外圈的轴向相对热位移。本发明提出了一种应用于高速主轴单元轴承套圈轴向热位移的确定方法,考虑了高速下含接触角动态变化的轴承球体热位移对轴承内外圈热位移的影响,使获取的高速主轴单元支承刚度更加准确合理,并且便于设计更为合理的轴承。
Description
技术领域
本发明属于机械领域,尤其涉及一种高速主轴单元轴承套圈轴向热位移确定方法。
背景技术
角接触球轴承具有摩擦力矩小、极限转速高、结构简单、噪音小等优点,普遍应用于机床高速主轴单元。而高速下,受轴承摩擦发热的影响,轴承套圈将产生一定的热位移,它直接影响着轴承刚度的大小,轴承刚度则直接影响着机床主轴的临界转速、动刚度、固有频率和振型等动力学特性以及机床的加工精度。因此,准确的计算轴承套圈热位移对预测高速主轴单元的动力学特性和机床加工精度具有重要意义。
目前,高速主轴单元中角接触球轴承内外套圈相对热位移的计算方法主要是通过建立考虑热效应的球轴承拟静力学模型和主轴单元热传递模型,使用考虑轴承座和转轴热位移影响下的轴承套圈热位移计算公式进行计算。
因此在现有技术中,计算高速角接触球轴承套圈热位移时,通常忽略了滚动体的热膨胀对内外套圈热位移的影响,或者在计算过程中,将轴承内外圈接触角当成常量来计算滚动体产生的热位移。事实上,由于轴承的摩擦发热,主轴单元各个零部件都会产生热位移,并且轴承滚动体与内外圈的接触角时刻在发生变化。因此,要想准确求得轴承套圈热位移大小,必须要综合考虑主轴单元各个零部件热位移和轴承内外圈接触角动态变化等对轴承套圈热位移的影响,从而获取更精确合理的主轴单元支承刚度,以便更准确的预测高速主轴单元的动力学特性和机床的加工精度。高速主轴一般指转速超过10000转/min。
发明内容
为解决上述问题,本发明提供了一种高速主轴单元轴承套圈轴向热位移确定方法。本发明提出了一种应用于高速主轴单元轴承套圈热位移的确定方法,考虑了高速下含接触角动态变化的轴承球体热位移对轴承内外圈热位移的影响,使获取的高速主轴单元支承刚度更加准确合理,并且便于设计更为合理的轴承。
为达到上述技术效果,本发明的技术方案是:
一种高速主轴单元轴承套圈轴向热位移确定方法,包括如下步骤:
步骤一、建立高速下的角接触球轴承内部几何相容关系,得到位置角ψj处的轴承内外圈滚道最终的轴向和径向相对距离;
步骤二、得到角接触球轴承沟道接触的变形几何相容方程;
步骤三、求出稳定状态时轴承的各向位移值,建立轴承所受外载荷和位移之间的关系,轴承受到联合载荷的作用,即受轴向载荷Fa、径向载荷Fr和力矩载荷M的综合作用;根据预紧及配置方式得到对应的轴承内外圈的轴向相对热位移。
进一步的改进,所述步骤一包括如下步骤:
在主轴单元工作过程中,受热变形影响,轴承外圈与轴承座之间的配合关系由间隙配合转变为过盈配合,此时假设轴承外圈固定不动,即轴承外滚道曲率中心不动,综合考虑主轴单元热变形和轴承内圈离心变形,建立高速下的角接触球轴承内部几何相容关系,得到位置角ψj处的轴承内外圈滚道最终的轴向和径向相对距离为:
A1j=BDbsinα0+δa+Riθcosψj+ua
A2j=BDb cosα0+δr cosψj+ur+uc
式中,A1j为轴承内外圈轴向相对距离,A2j为轴承内外圈径向相对距离,B=fi+fo-1,fi、fo为内、外滚道曲率半径系数,Db为球体直径,α0为轴承初始接触角,Ri为内圈曲率中心圆半径,δa、δr、θ分别为轴承轴向位移、径向位移和角位移,ua为轴承内外圈沟道轴向相对热位移,ur为轴承内外圈沟道径向相对热位移,uc为轴承内圈离心位移;
在任意角位置ψj处有:
式中,X1j为球体中心与外沟道曲率中心的轴向距离,X2j为球体中心与外沟道曲率中心的径向距离,αij、αoj分别为轴承内、外圈与球体接触角,δij、δoj分别为内、外圈与球体接触的弹性变形趋近量。
进一步的改进,所述步骤二中,
角接触球轴承沟道接触的变形几何相容方程为:
(A1j-X1j)2+(A2j-X2j)2-[(fi-0.5)Db+δij]2=0
X1j 2+X2j 2-[(fo-0.5)Db+δoj]2=0
在通过球体中心和轴承轴线的平面上,角位置ψj处的球体所受载荷进行分析,根据Hertz点接触变形关系及球体受力平衡关系可得:
其中,Kij、Koj分别为内、外圈滚道与球体接触的载荷变形系数,Mgj、Fcj分别为第j个球体受到的陀螺力矩和离心力;λij和λoj分别表示内、外圈沟道控制系数,对于高速情况,设球体所受陀螺力矩完全被球体与外圈沟道的接触摩擦力抵消,这时,取λij=0,λoj=2。
进一步的改进,所述步骤三包括如下步骤:
为求出稳定状态时轴承的各向位移值,建立轴承所受外载荷和位移之间的关系,轴承受到联合载荷的作用,即受轴向载荷Fa、径向载荷Fr和力矩载荷M的综合作用,则轴承载荷与位移的关系可用下式表达;
对于定压预紧:
对于定位预紧:
式中,Z为轴承球体个数;
高速运转下,主轴单元温度升高而产生热位移,轴承在轴承座、转轴以及自身的热影响下,产生热位移,设δi、δo分别为轴承内、外圈轴向热位移,ui、uo分别为轴承内、外圈径向热位移,ub表示球体热位移;
当轴承单联O型配置时,轴承内外圈的轴向热位移计算方法如下:
式中,ΔLs=λsΔTsLs,表示转轴轴向热伸长量,ΔLh=λhΔThLh,表示轴承座轴向热伸长量,其中,λs、λh分别表示转轴和轴承座热导率,ΔTs、ΔTh分别表示转轴和轴承座平均温度变化量,Ls、Lh分别表示转轴和轴承座原始长度;
不考虑球体热膨胀影响时,单联O型配置下的内外圈沟道处的轴向相对热位移u'a为:
事实上,轴承高速运转过程中,球体自身也将发生热膨胀,这样,修正后的计算公式如下:
式中,ub=λbΔTbDb,表示球体热膨胀量,αi、αo分别表示球体与内、外圈的动态接触角,其中,λb表示球体热导率,ΔTb表示球体中心温度变化量;
轴承单联X型配置,此时轴承内外圈的轴向相对热位移计算方法如下:
本发明至少具有以下优点:
(1)提出了更合理的不同安装方式下的轴承套圈热位移计算方法。
(2)根据需要可以设计出更为合理的轴承。
附图说明
图1为本发明中建立的角接触球轴承内部几何相容关系;
图2为本发明中建立的角接触球轴承滚动体转动时的受力图;
图3为本发明中建立的角接触球轴承热膨胀下的套圈位移示意图;
图4(a)为本发明中建立的角接触球轴承O型配置(背对背安装)单联配置方式示意图;
图4(b)为为本发明中建立的角接触球轴承X型配置(面对面安装)单联配置方式示意图;
图5为脂润滑高速主轴单元热传递模型示意图;
图6为轴承运转性能分析计算流程;
图7为轴承热位移图;
图8为轴承各向刚度图。
具体实施方式
以下通过具体实施方式并且结合附图对本发明的技术方案作具体说明。
实施例1
本发明包括以下步骤:
步骤一:建立角接触球轴承拟静力学模型
在主轴单元工作过程中,受热变形影响,轴承外圈与轴承座之间的配合关系通常由间隙配合转变为过盈配合,此时可假设轴承外圈固定不动,即轴承外滚道曲率中心不动,综合考虑主轴单元热变形和轴承内圈离心变形,建立高速下的角接触球轴承内部几何相容关系(见图1),得到位置角ψj处的轴承内外圈滚道最终的轴向和径向相对距离为:
A1j=BDbsinα0+δa+Riθcosψj+ua
A2j=BDbcosα0+δrcosψj+ur+uc
式中,A1j为轴承内外圈轴向相对距离,A2j为轴承内外圈径向相对距离,B=fi+fo-1,fi、fo为内、外滚道曲率半径系数,Db为球体直径,α0为轴承初始接触角,Ri为内圈曲率中心圆半径,δa、δr、θ分别为轴承轴向位移、径向位移和角位移,ua为轴承内外圈沟道轴向相对热位移,ur为轴承内外圈沟道径向相对热位移,uc为轴承内圈离心位移。
在任意角位置ψj处有:
式中,X1j为球体中心与外沟道曲率中心的轴向距离,X2j为球体中心与外沟道曲率中心的径向距离,αij、αoj分别为轴承内、外圈与球体接触角,δij、δoj分别为内、外圈与球体接触的弹性变形趋近量。
故得到角接触球轴承沟道接触的变形几何相容方程为:
(A1j-X1j)2+(A2j-X2j)2-[(fi-0.5)Db+δij]2=0
X1j 2+X2j 2-[(fo-0.5)Db+δoj]2=0
在通过球体中心和轴承轴线的平面上,角位置ψj处的球体所受载荷如图2所示,其中λij和λoj分别表示内、外圈沟道控制系数,对于高速情况,可假设球体所受陀螺力矩完全被球体与外圈沟道的接触摩擦力抵消,这时,取λij=0,λoj=2。根据Hertz点接触变形关系及球体受力平衡关系可得:
式中,Kij、Koj分别为内、外圈滚道与球体接触的载荷变形系数,Mgj、Fcj分别为第j个球体受到的陀螺力矩和离心力。
为了求出稳定状态时轴承的各向位移值,还需要建立轴承所受外载荷和位移之间的关系,假设轴承受到联合载荷的作用,即受轴向载荷Fa、径向载荷Fr和力矩载荷M的综合作用,则轴承载荷与位移的关系可用下式表达。
对于定压预紧:
对于定位预紧:
式中,Z为轴承球体个数。
步骤二:确定不同安装方式下的轴承套圈热位移和离心位移计算公式
高速运转下,主轴单元温度升高而产生热位移,轴承在轴承座、转轴以及自身的热影响下,产生一定的热位移如图3所示,其中δi、δo分别为轴承内、外圈轴向热位移,ui、uo分别为轴承内、外圈径向热位移,ub表示球体热位移。
轴承在不同配置方式下的轴向热伸长量计算方法不同,这里针对两种单联配置方式下的轴承热伸长量计算方法进行说明,双联或多联配置方式下的计算方法依此类推。如图4所示,图(a)表示轴承单联O型配置,此时轴承内外圈的轴向热位移计算方法如下:
式中,ΔLs=λsΔTsLs,表示转轴轴向热伸长量,ΔLh=λhΔThLh,表示轴承座轴向热伸长量,其中,λs、λh分别表示转轴和轴承座热导率,ΔTs、ΔTh分别表示转轴和轴承座平均温度变化量,Ls、Lh分别表示转轴和轴承座原始长度。
不考虑球体热膨胀影响时,单联O型配置下的内外圈沟道处的轴向相对热位移u'a为:
事实上,轴承高速运转过程中,球体自身也将发生热膨胀,这样,修正后的计算公式如下:
式中,ub=λbΔTbDb,表示球体热膨胀量,αi、αo分别表示球体与内、外圈的动态接触角,其中,λb表示球体热导率,ΔTb表示球体中心温度变化量。
图(b)表示轴承单联X型配置,此时轴承内外圈的轴向相对热位移计算方法如下:
对于O型或X型的轴承配置方式,考虑转轴热膨胀影响的轴承内圈沟道径向热位移计算公式相同,如下所示:
式中,λi表示轴承内圈热导率,ΔTi表示轴承内圈温度变化量,di表示轴承内圈沟道直径,μs表示转轴泊松比,d表示轴承内圈直径。
对于O型或X型的轴承配置方式,考虑轴承座热膨胀影响的轴承外圈沟道径向热位移计算公式相同,如下所示:
uo=λhΔTh(1+μh)Do
式中,λh表示轴承座热导率,μh表示轴承座泊松比,Do表示外圈沟道直径。
故修正后的轴承内外圈沟道径向相对热位移的计算公式如下:
ur=ui-uo-ubcosαi-ubcosαo
对于O型或X型的轴承配置方式,高速下,轴承内圈产生的离心位移计算公式如下:
式中,ρi表示轴承内圈密度,ωi表示轴承内圈转速,Ei表示轴承内圈弹性模量,μi表示轴承内圈泊松比。
步骤三:选取算例,建立主轴单元热传递方程
这里采用一般性的脂润滑高速主轴单元热传递模型对本发明进行说明(见图5),依据热传递模型可得主轴单元热传递方程如下:
式中,Hi、Ho分别为轴承内、外圈滚道接触区的摩擦热,TLi、TLo分别为轴承内、外圈滚道接触区温度,T∞为环境温度,Tb为球体中心温度,Ri、Ro分别为轴承内、外圈的热传导热阻,Rs、Rh分别为转轴和轴承座的复合换热热阻,即同时包含了零件自身热传导热阻和零件表面与周围环境的对流换热热阻,R1、R2为并联热阻且 Rb为球体热传导热阻,RLi、RLo分别为内外圈滚道润滑脂热传导热阻。
对于转轴和轴承座的复合换热热阻,计算方法如下:
Rs=Rrads+Raxs
式中,Rrads为转轴径向热传导热阻,Raxs为转轴轴向热传导热阻与轴向热对流换热热阻之和,Rradh为轴承座径向热传导热阻与径向热对流换热热阻之和,Raxh为轴承座轴向热传导热阻与轴向热对流换热热阻之和。
下面给出本算例中各个热阻的计算公式。
内圈径向热传导热阻Ri:
式中,Bi为轴承内圈宽度。
外圈径向热传导热阻Ro:
式中,D为轴承外径,λo为轴承外圈热导率,Bo为轴承外圈宽度。
转轴径向热传导热阻Rrads:
转轴轴向复合换热热阻Raxs:
式中,hs为转轴与环境对流换热系数,且hs=23(Ts-T∞)0.25。
轴承座径向复合换热热阻Rradh:
式中,Dh为轴承座外径,hh为轴承座与环境对流换热系数,且hh=23(Th-T∞)0.25。
轴承座轴向复合换热热阻Raxh:
内滚道的润滑脂热传导阻抗RLi:
式中,λL为润滑脂热导率。
外滚道的润滑脂热传导阻抗RLo:
球体热传导热阻为:
下面给出本算例中轴承产热量的计算方法。
轴承高速运转过程中,可将轴承总摩擦力矩等额分成内外圈沟道分量,对于位置角ψj处的球体有:
式中,Mij、Moj分别为第j个球体的内、外圈沟道摩擦力矩分量,f0、f1是与轴承类型和润滑、负载有关的系数,η0为润滑油的运动粘度,ωcj为第j个球体的公转角速度,Qij、Qoj分别为第j个球体与内、外滚道的接触应力,Qimax、Qomax分别为球体与内、外滚道接触的最大应力。
对于外滚道控制,仅需考虑球体与内滚道的自旋摩擦力矩,即:
式中,Msij为球体与内滚道的自旋摩擦力矩,μsi为球体与内滚道的摩擦系数,aij为第j个球体与内滚道的赫兹接触椭圆长半轴,Σij为第j个球体与内滚道赫兹接触椭圆的第二类完全积分。
故内、外滚道接触区的摩擦发热分别为:
Hij=ωcj.Mij+ωbj.Msij
Hoj=ωcj.Moj
式中,Hij、Hoj分别为第j个球体与内、外滚道接触区的摩擦热,ωbj为球体的自旋角速度。
在本算例中,有如下关系:
步骤四:拟定轴承运转性能的计算流程
为准确求得高速主轴单元轴承热位移的大小并了解轴承套圈热位移对轴承运转性能的影响,设计了一个轴承运转性能分析计算的流程,如图6所示。针对轴承模型的强非线性,应用Newton-Raphson迭代法对其进行求解。轴承运转性能计算过程如下:输入已知参数:材料、几何参数、预紧方式、载荷和转速后,设定轴承轴向位移δa、径向位移δr和角位移θ的初始值,然后再设定中间变量X1j、X2j、δij、δoj的初始值。通过牛顿迭代法(Newton-Raphson)对建好的轴承拟静力学模型进行求解,直到中间变量X1j、X2j、δij、δoj满足设定的精度才进行下一步,否则更改中间变量值,反复迭代,中间变量达到收敛精度后,根据主轴单元热传递模型求出主轴单元温度分布,进而求得轴承内外圈沟道的热位移和离心位移,然后对轴承拟静力学模型进行修正并迭代求解,直到轴承各向位移值达到设定的收敛精度,否则返回第一步,改变轴承初始位移值δa、δr和θ,再次重复整个计算流程,当整个计算过程达到设定的收敛精度后,输出计算结果。
步骤五:根据拟定的计算流程,求解得到高速主轴单元的轴承运转性能。
本发明以一种脂润滑主轴单元为算例,根据拟定的计算流程,求解得到主轴单元轴承运转性能参数,包括轴承热位移和轴承支承刚度等随转速的变化规律,将利用现有技术分析得到的轴承性能参数与依据本发明得到的轴承性能参数作对比,得到了一些新的结论。本算例中的轴承选用角接触混合陶瓷球轴承,具体参数示于表1,轴承座材料为HT300,转轴材料为40Cr,预紧方式为定压预紧,O型配置,预紧载荷为50N,环境温度为25℃。
表1轴承参数
图7示出了本算例中定压预紧下的轴承热位移变化曲线,图中不考虑球体热位移曲线是依据现有技术得到的,考虑球体热位移是依据本发明方法得到的。从图7中可以看到,不管是否考虑球体热位移,轴承套圈轴向相对热位移均随转速上升而不断增大,且两种情况下的轴承轴向相对热位移相差不大,而考虑球体热位移与否对轴承径向相对热位移的影响比较明显,不考虑球体热位移时,轴承径向相对热位移随转速上升而增加缓慢,考虑球体热位移时,轴承径向相对热位移随转速上升而增加较快,这是由于考虑球体热位移时,球体受热膨胀,使得轴承内圈沟道热位移减小,外圈沟道热位移增大,所以轴承内外圈沟道径向相对热位移变大。
图8示出了本算例中定压预紧下的轴承各向刚度变化曲线,从图8中可以看到,考虑球体热位移与否对轴承各向刚度的影响较小,轴向刚度和角刚度均随转速上升而迅速下降,这是因为转速上升,轴承轴向相对热位移和径向相对热位移均增大,导致轴承内圈接触角减小,又由于径向相对热位移明显小于轴向相对热位移,所以外圈接触角变化不大,导致轴承轴向接触刚度下降,角刚度受轴承接触刚度的影响也随之下降。转速上升,轴承径向刚度增加,这是因为内圈接触角减小,球体与内外圈的径向接触热变形增大,导致球体与轴承内外圈的径向接触载荷增加,从而使得轴承径向刚度增加。
上述仅为本发明的一个具体导向实施方式,但本发明的设计构思并不局限于此,凡利用此构思对本发明进行非实质性的改动,均应属于侵犯本发明的保护范围的行为。
Claims (4)
1.一种高速主轴单元轴承套圈轴向热位移确定方法,其特征在于,包括如下步骤:
步骤一、建立高速下的角接触球轴承内部几何相容关系,得到位置角ψj处的轴承内外圈滚道最终的轴向和径向相对距离;
步骤二、得到角接触球轴承沟道接触的变形几何相容方程;
步骤三、求出稳定状态时轴承的各向位移值,建立轴承所受外载荷和位移之间的关系,轴承受到联合载荷的作用,即受轴向载荷Fa、径向载荷Fr和力矩载荷M的综合作用;根据预紧及配置方式得到对应的轴承内外圈的轴向相对热位移。
2.如权利要求1所述的高速主轴单元轴承套圈轴向热位移确定方法,其特征在于,所述步骤一包括如下步骤:
在主轴单元工作过程中,受热变形影响,轴承外圈与轴承座之间的配合关系由间隙配合转变为过盈配合,此时假设轴承外圈固定不动,即轴承外滚道曲率中心不动,综合考虑主轴单元热变形和轴承内圈离心变形,建立高速下的角接触球轴承内部几何相容关系,得到位置角ψj处的轴承内外圈滚道最终的轴向和径向相对距离为:
A1j=BDbsinα0+δa+Riθcosψj+ua
A2j=BDb cosα0+δr cosψj+ur+uc
式中,A1j为轴承内外圈轴向相对距离,A2j为轴承内外圈径向相对距离,B=fi+fo-1,fi、fo为内、外滚道曲率半径系数,Db为球体直径,α0为轴承初始接触角,Ri为内圈曲率中心圆半径,δa、δr、θ分别为轴承轴向位移、径向位移和角位移,ua为轴承内外圈沟道轴向相对热位移,ur为轴承内外圈沟道径向相对热位移,uc为轴承内圈离心位移;
在任意角位置ψj处有:
式中,X1j为球体中心与外沟道曲率中心的轴向距离,X2j为球体中心与外沟道曲率中心的径向距离,αij、αoj分别为轴承内、外圈与球体接触角,δij、δoj分别为内、外圈与球体接触的弹性变形趋近量。
3.如权利要求2所述的高速主轴单元轴承套圈轴向热位移确定方法,其特征在于,所述步骤二中,
角接触球轴承沟道接触的变形几何相容方程为:
(A1j-X1j)2+(A2j-X2j)2-[(fi-0.5)Db+δij]2=0
X1j 2+X2j 2-[(fo-0.5)Db+δoj]2=0
在通过球体中心和轴承轴线的平面上,角位置ψj处的球体所受载荷进行分析,根据Hertz点接触变形关系及球体受力平衡关系可得:
其中,Kij、Koj分别为内、外圈滚道与球体接触的载荷变形系数,Mgj、Fcj分别为第j个球体受到的陀螺力矩和离心力;λij和λoj分别表示内、外圈沟道控制系数,对于高速情况,设球体所受陀螺力矩完全被球体与外圈沟道的接触摩擦力抵消,这时,取λij=0,λoj=2。
4.如权利要求3所述的高速主轴单元轴承套圈轴向热位移确定方法,其特征在于,所述步骤三包括如下步骤:
为求出稳定状态时轴承的各向位移值,建立轴承所受外载荷和位移之间的关系,轴承受到联合载荷的作用,即受轴向载荷Fa、径向载荷Fr和力矩载荷M的综合作用,则轴承载荷与位移的关系可用下式表达;
对于定压预紧:
对于定位预紧:
式中,Z为轴承球体个数;
高速运转下,主轴单元温度升高而产生热位移,轴承在轴承座、转轴以及自身的热影响下,产生热位移,设δi、δo分别为轴承内、外圈轴向热位移,ui、uo分别为轴承内、外圈径向热位移,ub表示球体热位移;
当轴承单联O型配置时,轴承内外圈的轴向热位移计算方法如下:
式中,ΔLs=λsΔTsLs,表示转轴轴向热伸长量,ΔLh=λhΔThLh,表示轴承座轴向热伸长量,其中,λs、λh分别表示转轴和轴承座热导率,ΔTs、ΔTh分别表示转轴和轴承座平均温度变化量,Ls、Lh分别表示转轴和轴承座原始长度;
不考虑球体热膨胀影响时,单联O型配置下的内外圈沟道处的轴向相对热位移u'a为:
事实上,轴承高速运转过程中,球体自身也将发生热膨胀,这样,修正后的计算公式如下:
式中,ub=λbΔTbDb,表示球体热膨胀量,αi、αo分别表示球体与内、外圈的动态接触角,其中,λb表示球体热导率,ΔTb表示球体中心温度变化量;
轴承单联X型配置,此时轴承内外圈的轴向相对热位移计算方法如下:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811284426.9A CN109550979B (zh) | 2018-10-31 | 2018-10-31 | 一种高速主轴单元轴承套圈轴向热位移确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811284426.9A CN109550979B (zh) | 2018-10-31 | 2018-10-31 | 一种高速主轴单元轴承套圈轴向热位移确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109550979A true CN109550979A (zh) | 2019-04-02 |
CN109550979B CN109550979B (zh) | 2020-08-04 |
Family
ID=65865655
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811284426.9A Active CN109550979B (zh) | 2018-10-31 | 2018-10-31 | 一种高速主轴单元轴承套圈轴向热位移确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109550979B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111159880A (zh) * | 2019-12-26 | 2020-05-15 | 熵智科技(深圳)有限公司 | 一种球轴承接触刚度计算方法 |
CN111209686A (zh) * | 2020-01-16 | 2020-05-29 | 重庆大学 | 基于复合形法的滚动轴承多体润滑性能求解方法 |
CN111310373A (zh) * | 2020-02-11 | 2020-06-19 | 重庆大学 | 基于解析法的轴系统热特性分析方法、热误差建模方法和热误差补偿系统 |
CN113704981A (zh) * | 2021-08-11 | 2021-11-26 | 武汉理工大学 | 温升过程中高速轴承时变动力学行为分析方法 |
CN116357684A (zh) * | 2023-03-23 | 2023-06-30 | 山东智研连合轴承有限公司 | 一种配对角接触球轴承的组配方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015009351A (ja) * | 2013-07-02 | 2015-01-19 | 株式会社ジェイテクト | 主軸装置 |
CN107756128A (zh) * | 2017-08-30 | 2018-03-06 | 华中科技大学 | 热变形补偿方法 |
-
2018
- 2018-10-31 CN CN201811284426.9A patent/CN109550979B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015009351A (ja) * | 2013-07-02 | 2015-01-19 | 株式会社ジェイテクト | 主軸装置 |
CN107756128A (zh) * | 2017-08-30 | 2018-03-06 | 华中科技大学 | 热变形补偿方法 |
Non-Patent Citations (2)
Title |
---|
李芳芳: "基于轴承热位移的电主轴系统动态特性分析方法研究", 《中国优秀硕士学位论文全文数据库 工程科技I辑》 * |
陈维福: "机床进给系统轴承热特性分析及试验研究", 《中国优秀硕士学位论文全文数据库 工程科技I辑》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111159880A (zh) * | 2019-12-26 | 2020-05-15 | 熵智科技(深圳)有限公司 | 一种球轴承接触刚度计算方法 |
CN111209686A (zh) * | 2020-01-16 | 2020-05-29 | 重庆大学 | 基于复合形法的滚动轴承多体润滑性能求解方法 |
CN111209686B (zh) * | 2020-01-16 | 2024-04-19 | 重庆大学 | 基于复合形法的滚动轴承多体润滑性能求解方法 |
CN111310373A (zh) * | 2020-02-11 | 2020-06-19 | 重庆大学 | 基于解析法的轴系统热特性分析方法、热误差建模方法和热误差补偿系统 |
CN111310373B (zh) * | 2020-02-11 | 2022-02-18 | 重庆大学 | 基于解析法的轴系统热特性分析方法、热误差建模方法和热误差补偿系统 |
CN113704981A (zh) * | 2021-08-11 | 2021-11-26 | 武汉理工大学 | 温升过程中高速轴承时变动力学行为分析方法 |
CN113704981B (zh) * | 2021-08-11 | 2024-04-09 | 武汉理工大学 | 温升过程中高速轴承时变动力学行为分析方法 |
CN116357684A (zh) * | 2023-03-23 | 2023-06-30 | 山东智研连合轴承有限公司 | 一种配对角接触球轴承的组配方法 |
CN116357684B (zh) * | 2023-03-23 | 2023-11-03 | 山东智研连合轴承有限公司 | 一种配对角接触球轴承的组配方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109550979B (zh) | 2020-08-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109550979A (zh) | 一种高速主轴单元轴承套圈轴向热位移确定方法 | |
Liu et al. | A combined acoustic and dynamic model of a defective ball bearing | |
CN109489949B (zh) | 一种应用于高速主轴单元热动态性能的分析方法 | |
CN109027017B (zh) | 一种空间滚动轴承磨损状态评估方法 | |
CN110617953B (zh) | 新能源汽车高速电机轴承-转子系统动态特性分析方法 | |
Liu et al. | Dynamic stability analysis of cages in high-speed oil-lubricated angular contact ball bearings | |
Wang et al. | A comprehensive study on dynamic performance of ball bearing considering bearing deformations and ball-inner raceway separation | |
Zhou et al. | Thermal network model and experimental validation for a motorized spindle including thermal–mechanical coupling effect | |
Lei et al. | An integrated model to characterize comprehensive stiffness of angular contact ball bearings | |
Cheng et al. | Research on mechanical characteristics of fault-free bearings based on centrifugal force and gyroscopic moment | |
Wang et al. | Thermal-mechanical fully coupled analysis of high-speed angular contact ball bearings | |
CN111177898B (zh) | 基于bp神经网络的滚动轴承-转子系统耦合性能求解方法 | |
He et al. | Study on thermo-mechanical coupling characteristics of angle contact ball bearing with fix-position preload | |
CN110162874A (zh) | 配置振荡热管电主轴系统及其热-结构耦合特性建模方法 | |
Rabréau et al. | Influence of bearing kinematics hypotheses on ball bearing heat generation | |
Li et al. | Thermal Performances Prediction Analysis of High Speed Feed Shaft Bearings Under Actual Working Condition | |
Muszyński et al. | Modeling of heat distribution at the spindle bearing test rig | |
CN113775645A (zh) | 一种考虑油膜热效应的高速球轴承动态性能分析方法 | |
CN114169157A (zh) | 一种考虑界面摩擦的角接触球轴承动态特性计算方法 | |
Wang et al. | Prediction of the friction torque in grease lubricated angular contact ball bearings using grey system theory | |
Reddy | Thermal stress analysis of a ball bearing by finite element method | |
Lei et al. | A new method for computing contact angle of high speed ball bearing | |
BRECHER et al. | NEW BEARING STEEL FOR HIGH-SPEED APPLICATIONS. | |
Su et al. | Effect of temperature rise on a tandem of ball-bearings with an O-configuration in motorized spindle under oil-air lubrication | |
Liu | Kinematic investigation and fatigue life analysis of angular contact ball bearing in wide speed range based on raceway friction |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |