CN106934165B - 一种交叉圆锥滚子轴承静刚度的快速提取方法 - Google Patents

一种交叉圆锥滚子轴承静刚度的快速提取方法 Download PDF

Info

Publication number
CN106934165B
CN106934165B CN201710166615.5A CN201710166615A CN106934165B CN 106934165 B CN106934165 B CN 106934165B CN 201710166615 A CN201710166615 A CN 201710166615A CN 106934165 B CN106934165 B CN 106934165B
Authority
CN
China
Prior art keywords
row
tapered
rollers
tapered rollers
tapered roller
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.)
Active
Application number
CN201710166615.5A
Other languages
English (en)
Other versions
CN106934165A (zh
Inventor
牛荣军
徐金超
倪艳光
邓四二
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Henan University of Science and Technology
Original Assignee
Henan University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Henan University of Science and Technology filed Critical Henan University of Science and Technology
Priority to CN201710166615.5A priority Critical patent/CN106934165B/zh
Publication of CN106934165A publication Critical patent/CN106934165A/zh
Application granted granted Critical
Publication of CN106934165B publication Critical patent/CN106934165B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Rolling Contact Bearings (AREA)

Abstract

一种交叉圆锥滚子轴承静刚度的快速提取方法，本发明首先基于Hertz接触理论和滚动轴承设计方法建立交叉圆锥滚子轴承的精确数学模型，分出第一列圆锥滚子和第二列圆锥滚子进行单独计算，采用切片单独计算载荷计算法向接触载荷的方法，并采用数值方法对力学平衡方程进行求解，将求解值代入刚度矩阵得出交叉圆锥滚子轴承静刚度值。本发明给出了一种联合外载工况下快速精确计算交叉圆锥滚子轴承刚度的方法，以辅助工程应用，指导交叉圆锥滚子轴承的设计选型。

Description

一种交叉圆锥滚子轴承静刚度的快速提取方法
技术领域
本发明涉及交叉圆锥滚子轴承领域，具体说的是一种交叉圆锥滚子轴承静刚度的快速提取方法。
背景技术
轴承的刚度是在负荷方向上轴承内、外圈产生单位的相对弹性位移量所需加的外载荷，轴承的刚度反映了其在承受负载时抵抗弹性变形的能力，是轴承结构对弹性变形难易程度的表征。轴承刚度能判断出交叉圆锥滚子轴承的结构和选型是否合理，对于轴承的设计选型具有重要的指导意义。交叉滚子轴承因其高刚性、高旋转精度、高转速等特点而广泛应用在精密旋转工作台、立式车床、卧式车床、立式磨床等场合。交叉圆锥滚子轴承受力情况复杂，快速准确获取其刚度特性十分困难。目前，国内外学者对交叉圆锥滚子轴承研究较少，对交叉圆锥滚子轴承刚度的计算也没有给出具体的方法。
发明内容
为解决上述技术问题，本发明提供一种交叉圆锥滚子轴承静刚度的快速提取方法，能够得到辅助工程应用、指导交叉圆锥滚子轴承设计的交叉圆锥滚子轴承静钢度。
为实现上述技术目的，所采用的技术方案是：一种交叉圆锥滚子轴承静刚度的快速提取方法，该方法包括以下步骤：
步骤一、建立交叉圆锥滚子轴承的力学模型，设置位置角
Figure BDA0001250003980000011
为零度且与其方向相同的滚子为第一列圆锥滚子，与第一列圆锥滚子交叉垂直的滚子为第二列圆锥滚子；
步骤二、将第一列圆锥滚子和第二列圆锥滚子沿其中心线方向切分成n片，得到第一列圆锥滚子每个切片上的法向变形量δn1(j)与内圈相对外圈在轴向力Fa作用下产生轴向位移δa、在径向力Fr作用下产生径向位移δr、在倾覆力矩M作用下产生倾角位移θ的关系式，以及得到第二列圆锥滚子每个切片上的法向变形量δn2(j)与内圈相对外圈在轴向力Fa作用下产生轴向位移δa、在径向力Fr作用下产生径向位移δr、在倾覆力矩M作用下产生倾角位移θ的关系式；
步骤三、根据Hertz线接触理论，由圆锥滚子每个切片上的法向变形量δn(j)和圆锥滚子每个切片上法向接触载荷qi(j)关系式，得到圆锥滚子的法向接触载荷为所有切片法向接触载荷qi(j)的和，根据第一列圆锥滚子的法向接触载荷得到在第一列圆锥滚子位置角
Figure BDA0001250003980000022
处的轴向分力
Figure BDA0001250003980000023
径向分力
Figure BDA0001250003980000024
轴向分力产生的转矩
Figure BDA0001250003980000025
径向分力产生的转矩
Figure BDA0001250003980000026
同理，得到第二列圆锥滚子的法向接触载荷
Figure BDA0001250003980000027
在第二列圆锥滚子位置角
Figure BDA0001250003980000028
处的轴向分力
Figure BDA0001250003980000029
径向分力
Figure BDA00012500039800000210
轴向分力产生的转矩
Figure BDA00012500039800000211
径向分力产生的转矩
Figure BDA00012500039800000212
步骤四、轴承内圈在外部载荷和滚子反作用力下处于平衡状态，根据第一列圆锥滚子位置角
Figure BDA00012500039800000213
处的轴向分力
Figure BDA00012500039800000214
径向分力
Figure BDA00012500039800000215
轴向分力产生的转矩径向分力产生的转矩
Figure BDA00012500039800000217
以及第二列圆锥滚子的法向接触载荷在第二列圆锥滚子位置角
Figure BDA00012500039800000219
处的轴向分力
Figure BDA00012500039800000220
径向分力
Figure BDA00012500039800000221
轴向分力产生的转矩
Figure BDA00012500039800000222
径向分力产生的转矩
Figure BDA00012500039800000223
建立以内圈相对外圈在轴向力Fa作用下产生轴向位移δa、在径向力Fr作用下产生径向位移δr、在倾覆力矩M作用下产生倾角位移θ为未知数的内圈的三元力学平衡方程，运用Newton-Raphson迭代法求解力学平衡方程得到δa、δr、θ的值；
步骤五、根据轴承内、外圈在负荷方向上产生单位的相对弹性位移量所需的外加载荷为滚动轴轴承的钢度，建立第一列圆锥滚子的刚度矩阵和第二列圆锥滚子的刚度矩阵，将δa、δr、θ的值代入第一列圆锥滚子的刚度矩阵和第二列圆锥滚子的刚度矩阵，刚度矩阵为第一列圆锥滚子的刚度矩阵和第二列圆锥滚子的刚度矩阵之和为交叉圆锥滚子轴承静刚度。
本发明所述的每个切片上的法向变形量δn1(j)与内圈相对外圈在轴向力Fa作用下产生轴向位移δa、在径向力Fr作用下产生径向位移δr、在倾覆力矩M作用下产生倾角位移θ的关系式为
Figure BDA00012500039800000224
其中，dm1(j)为第一列圆锥滚子每个切片的中心节圆直径，
Figure BDA00012500039800000225
为滚子母线与其中心线的夹角，h为每片圆锥滚子的厚度，leff1(j)为第一列滚子每个切片的中心位置；
leff1(j)＝0.5h(1+2(j-1))。
本发明所述的第二列圆锥滚子每个切片上的法向变形量δn2(j)与内圈相对外圈在轴向力Fa作用下产生轴向位移δa、在径向力Fr作用下产生径向位移δr、在倾覆力矩M作用下产生倾角位移θ的关系式为
Figure BDA0001250003980000031
其中，dm2(j)为第二列圆锥滚子每个切片的中心节圆直径，为滚子母线与其中心线的夹角，h为每片圆锥滚子的厚度，leff2(j)为第二列滚子每个切片的中心位置；
leff2(j)＝0.5h(1+2(j-1))。
本发明所述的圆锥滚子每个切片上的法向变形量δn(j)和圆锥滚子每个切片上法向接触载荷qi(j)关系式为
Figure BDA0001250003980000033
Kne(j)＝6.24×104lwe 0.82Dw(j)0.11/(kcount*h)0.18[1+ci 0.9cos(αe-αi)]-1.11
式中，
Figure BDA0001250003980000034
kcount表示接触的圆锥滚子的片数，δn(j)＞0时，表示该片圆锥滚子接触，lwe表示滚子有效接触长度，Dw(j)表示每片圆锥滚子的直径，αi为圆锥滚子与内圈的接触角，αf为圆锥滚子端部与外圈的接触角，αe为圆锥滚子母线与外圈的接触角。
本发明所述的第一列圆锥滚子的法向接触载荷
Figure BDA0001250003980000035
的表达式为
Figure BDA0001250003980000036
其中，
Figure BDA0001250003980000037
为第一列圆锥滚子每个切处所受法向接触载荷。
本发明所述的第一列圆锥滚子的法向接触载荷
Figure BDA0001250003980000038
得到在第一列圆锥滚子位置角处的轴向分力
Figure BDA00012500039800000310
径向分力
Figure BDA00012500039800000311
轴向分力产生的转矩
Figure BDA00012500039800000312
径向分力产生的转矩
Figure BDA00012500039800000313
的表达式为
Figure BDA0001250003980000041
Figure BDA0001250003980000042
其中，lwe表示滚子有效接触长度，Dw(j)表示每片圆锥滚子的直径，leff1(j)为第一列滚子每个切片的中心位置；
leff1(j)＝0.5h(1+2(j-1))。
本发明所述的第二列圆锥滚子的法向接触载荷
Figure BDA0001250003980000045
的表达式为
Figure BDA0001250003980000046
其中，
Figure BDA0001250003980000047
为第二列圆锥滚子每个切片所受法向接触载荷。
本发明所述的第二列圆锥滚子的法向接触载荷
Figure BDA0001250003980000048
在第二列圆锥滚子位置角
Figure BDA0001250003980000049
处的轴向分力径向分力
Figure BDA00012500039800000411
轴向分力产生的转矩
Figure BDA00012500039800000412
径向分力产生的转矩
Figure BDA00012500039800000413
的表达式为
Figure BDA00012500039800000414
Figure BDA00012500039800000416
Figure BDA00012500039800000417
其中，lwe表示滚子有效接触长度，Dw(j)表示每片圆锥滚子的直径，leff2(j)为第二列滚子每个切片的中心位置；
leff2(j)＝0.5h(1+2(j-1))。
本发明所述的三元力学平衡方程为
Figure BDA0001250003980000051
Figure BDA0001250003980000052
本发明所述的第一列圆锥滚子的刚度矩阵为
本发明所述的第二列圆锥滚子的刚度矩阵为
Figure BDA0001250003980000055
本发明所述的位置角i1＝1,2,3,…,z1,z1为第一列圆锥滚子的个数，z1＝z/2，z为圆锥滚子轴承所有滚子的个数，
Figure BDA0001250003980000061
i2＝1,2,3,…,z2，,z2为第二列圆锥滚子的个数，z2＝z/2。
本发明的有益效果是：本发明首先基于Hertz接触理论和滚动轴承设计方法建立交叉圆锥滚子轴承的精确数学模型，分出第一列圆锥滚子和第二列圆锥滚子进行单独计算，采用切片单独计算载荷计算法向接触载荷的方法，并采用数值方法对力学平衡方程进行求解，将求解值代入刚度矩阵得出交叉圆锥滚子轴承静刚度值。本发明给出的获取交叉圆锥滚子轴承静刚度的快速提取方法，给出了一种联合外载(轴向载荷、径向载荷和倾覆力矩)工况下快速精确计算交叉圆锥滚子轴承刚度的方法，以辅助工程应用，指导交叉圆锥滚子轴承的设计选型。
附图说明
图1为本发明交叉圆锥滚子轴承受力与内圈位移图；
图2为本发明交叉圆锥滚子轴承滚子位置角示意图；
图3为本发明交叉圆锥滚子轴承接触角示意图；
图4为本发明基于切片法的圆锥滚子变形和受载关系示意图；
图5为本发明轴向刚度Kaa随轴向力Fa的变化图；
图6为径向刚度Krr随径向力Fr的变化图；
图7为角刚度Kθθ随倾覆力矩M的变化图。
具体实施方式
本发明首先基于Hertz接触理论和滚动轴承设计方法建立交叉圆锥滚子轴承的精确数学模型，交叉圆锥滚子轴承内圈在外力(轴向载荷、径向载荷和倾覆力矩)和滚子反作用力下处于平衡状态，根据这种平衡状态列出以轴向位移δa、径向位移δr、角位移θ为未知量的平衡方程组，运用Newton-Raphson迭代法进行运算求解，得出δa、δr、θ的值，将所计算出的δa、δr、θ和轴承参数带入轴承刚度的计算模型中，即可得到轴承的刚度。
一种交叉圆锥滚子轴承静刚度的快速提取方法，包含以下步骤：
(1)建立交叉圆锥滚子轴承的力学模型：
假设轴承外圈固定，按照静力学的方法来建立交叉圆锥滚子轴承的力学模型，模型如图1，Fa、Fr、M分别表示交叉圆锥滚子轴承承受的轴向力、径向力和倾覆力矩，δa为轴向位移，δr为径向位移，θ为倾角位移。交叉圆锥滚子轴承有两列相互垂直的圆锥滚子及内外圈组成，记位置角为零度且与其方向相同的滚子为第一列圆锥滚子，为第一列圆锥滚子的位置角，
Figure BDA0001250003980000072
z1为第一列圆锥滚子的个数，z1＝z/2，z为圆锥滚子轴承所有滚子的个数；记与第一列圆锥滚子交叉垂直的滚子为第2列圆锥滚子，
Figure BDA0001250003980000073
为第二列圆锥滚子的位置角，
Figure BDA0001250003980000074
Figure BDA0001250003980000075
z2为第二列圆锥滚子的个数，z2＝z/2。位置角关系定义如图2所示。
(2)计算交叉圆锥滚子轴承的变形量：
内圈相对外圈在轴向力Fa作用下产生轴向位移δa，在径向力Fr作用下产生径向位移δr，在倾覆力矩M作用下产生倾角位移θ。记轴向位移引起的轴向变形为δa，径向位移引起的径向变形分量为
Figure BDA0001250003980000076
倾角位移引起的轴向变形分量为
Figure BDA0001250003980000077
(dm(j)代表每个切片的中心节圆直径)，倾角位移引起的径向变形分量为
Figure BDA0001250003980000078
(lwe表示滚子有效接触长度，
Figure BDA0001250003980000079
为滚子母线与其中心线的夹角，leff(j)为滚子每个切片的中心位置，β为内圈滚道母线与其中心线之间的夹角)。内圈相对外圈位移减小的力的方向为正，位移增加的力的方向为负。接触角定义如图3所示，αi为圆锥滚子与内圈的接触角，αf为圆锥滚子端部与外圈的接触角，αe为圆锥滚子母线与外圈的接触角。由于圆锥滚子沿母线接触位置处的直径不相同，因此产生的法向变形量也不同。为了更精确的分析圆锥滚子的载荷分布，用切片法对圆锥滚子接触变形关系进行研究分析。切片法示意图如图4，把圆锥滚子沿其中心线方向切分成n片，qi(j)代表每个切片上的法向载荷，dm(j)表示圆锥滚子每个切片的中心节圆直径，Dw(j)表示每片圆锥滚子的直径。
第一列圆锥滚子轴承变形量：
每个切片上的轴向变形量为
每个切片上的径向变形量为
Figure BDA0001250003980000081
每个切片上的法向变形量为
δn1(j)＝δa1sinαe+δr1cosαe (3)
式中dm1(j)为第一列圆锥滚子每个切片的中心节圆直径，h表示每片圆锥滚子的厚度，leff1(j)为第一列滚子每个切片的中心位置。
leff1(j)＝0.5h(1+2(j-1)) (4)
第二列圆锥滚子轴承变形量：
每个切片上的轴向变形量为
Figure BDA0001250003980000082
每个切片上的径向变形量为
Figure BDA0001250003980000083
每个切片上的法向变形量为
δn2(j)＝-δa2sinαe+δr2cosαe (7)
dm2(j)为第二列圆锥滚子每个切片的中心节圆直径，leff2(j)为第二列滚子每个切片的中心位置。
leff2(j)＝0.5h(1+2(j-1)) (8)
(3)计算交叉圆锥滚子轴承的力和力矩分量：
根据Hertz线接触理论，由圆锥滚子轴承法向变形量和法向接触载荷的关系，计算每片圆锥滚子在位置角处的法向接触载荷，每片圆锥滚子法向接触载荷为：
Figure BDA0001250003980000084
式中Kne(j)表示滚动体与内外圈之间总的负荷-变形常数：
Kne(j)＝6.24×104lwe 0.82Dw(j)0.11/(kcount*h)0.18[1+ci 0.9cos(αe-αi)]-1.11 (10)
式中
Figure BDA0001250003980000085
kcount表示接触的圆锥滚子的片数(δn(j)＞0时，表示该片圆锥滚子接触)。
第一列圆锥滚子轴承的力和力矩：
法向接触载荷为
Figure BDA0001250003980000091
法向接触载荷的轴向分力为
法向接触载荷的径向分力为
Figure BDA0001250003980000093
法向接触载荷轴向分力产生的转矩为
Figure BDA0001250003980000094
法向接触载荷径向分力产生的转矩为
Figure BDA0001250003980000095
Figure BDA0001250003980000096
为第一列圆锥滚子每片所受载荷。
第二列圆锥滚子轴承的力和力矩：
法向接触载荷为
Figure BDA0001250003980000097
法向接触载荷的轴向分力为
Figure BDA0001250003980000098
法向接触载荷的径向分力为
Figure BDA0001250003980000099
法向接触载荷轴向分力产生的转矩为
法向接触载荷径向分力产生的转矩为
Figure BDA0001250003980000102
Figure BDA0001250003980000103
为第二列圆锥滚子每片所受载荷。
(4)建立交叉圆锥滚子轴承力学平衡方程组：
交叉圆锥滚子轴承内圈在外力和滚子反作用力下处于平衡状态，受力简图如图5，规定内圈受到向下的径向力为正，向左的轴向力为正，逆时针的力矩为正。内圈的力学平衡方程为：
Figure BDA0001250003980000104
Figure BDA0001250003980000105
(5)获取交叉圆锥滚子轴承的载荷分布：
式(21)-(23)构成了以δa、δr、θ为未知量的三元非线性方程组，可运用Newton-Raphson迭代法进行数值运算求解，得出δa、δr、θ的值。然后，分别带入到式(1)-(3)和(5)-(7)计算两列交叉圆锥滚子轴承的法向变形量，再将法向变形量带入到式(11)和(16)，即可获取交叉圆锥滚子轴承两列的载荷分布。
(6)计算交叉圆锥滚子轴承的刚度矩阵：
滚动轴承的刚度是轴承内、外圈在负荷方向上产生单位的相对弹性位移量所需的外加载荷，刚度矩阵可用式(20)表达。
Figure BDA0001250003980000111
式中为主刚度，
Figure BDA0001250003980000113
为辅刚度，
Figure BDA0001250003980000114
为轴向刚度，用Kaa表示，
Figure BDA0001250003980000115
为径向刚度，用Krr表示，为角刚度，用Kθθ表示。
第一列圆锥滚子的刚度矩阵：
Figure BDA0001250003980000117
第二列圆锥滚子的刚度矩阵：
Figure BDA0001250003980000118
交叉圆锥滚子轴承总的刚度矩阵为：
K＝K1+K2 (23)
把求解的δa、δr、θ带入到刚度矩阵(21)-(23)中，即可得到交叉圆锥滚子轴承的刚度。
下面以某一具体的交叉圆锥滚子轴承为例进行说明：
某型号交叉圆锥滚子轴承，结构参数为：轴承的节圆直径dm＝536.4mm，滚动体大段直径DW＝26.99mm，滚动体小段直径DW1＝25.533mm，滚动体个数z＝52，两列滚动体个数z1＝z2＝z/2＝26，滚子有效长度lwe＝22，圆锥滚子与内圈的接触角αi＝41.2°，圆锥滚子端部与外圈的接触角αf＝46.9°，圆锥滚子母线与外圈的接触角αe＝45°。运行时受到轴向力141kN、径向力178kN、倾覆力矩170kN·m。
第一步：建立交叉圆锥滚子轴承的力学模型：
交叉圆锥滚子轴承有两列相互垂直的圆锥滚子及内外圈组成，Fa、Fr、M分别表示交叉圆锥滚子轴承承受的轴向力、径向力和倾覆力矩，δa为轴向位移，δr为径向位移，θ为倾角位移，模型如图1。交叉圆锥滚子轴承有两列相互垂直的圆锥滚子及内外圈组成，记位置角为零度且与其方向相同的滚子为第一列圆锥滚子，
Figure BDA0001250003980000121
为第一列圆锥滚子在位置角
Figure BDA0001250003980000122
处的法向载荷，
Figure BDA0001250003980000123
为第一列圆锥滚子的位置角，
Figure BDA0001250003980000124
z1为第一列圆锥滚子的个数，z为圆锥滚子轴承所有滚子的个数，z1＝z/2；
Figure BDA0001250003980000125
为第二列圆锥滚子在位置角
Figure BDA0001250003980000126
处的法向载荷，
Figure BDA0001250003980000127
为第二列圆锥滚子的位置角，z2为第二列圆锥滚子的个数，z2＝z/2，位置角关系定义如图2所示。
由于圆锥滚子在各个位置处的直径不相同，各个位置处的法向变形量、法向载荷、转矩也不同，为了更精确的分析研究圆锥滚子，用切片法对圆锥滚子进行研究分析。把圆锥滚子切成n片，qi(j)表示每个切片上的法向载荷。dm(j)为圆锥滚子每个切片的中心节圆直径，Dw(j)表示每片圆锥滚子的直径。切片法示意图如图3，β为内圈滚道母线与其中心线之间的夹角，为滚子母线与其中心线的夹角。图4为交叉圆锥滚子轴承接触角示意图，αi为圆锥滚子与内圈的接触角，αf为圆锥滚子端部与外圈的接触角，αe为圆锥滚子母线与外圈的接触角。
第二步：建立交叉圆锥滚子轴承静力学平衡方程组：
第一列圆锥滚子轴承法向变形量为：
Figure BDA0001250003980000131
式中dm1(j)为第一列圆锥滚子每个切片的中心节圆直径，lwe表示有效接触长度
h表示每片圆锥滚子的厚度，leff1(j)为第一列滚子每个切片的中心位置。
leff1(j)＝0.5h(1+2(j-1)) (2)
第二列圆锥滚子轴承法向变形量为：
式中dm2(j)为第二列圆锥滚子每个切片的中心节圆直径，leff2(j)为第二列滚子每个切片的中心位置。
leff2(j)＝0.5h(1+2(j-1)) (4)
根据Hertz线接触理论，由圆锥滚子轴承法向变形量和法向接触载荷的关系，计算两列圆锥滚子轴承在位置角处的法向接触载荷,圆锥滚子轴承法向接触载荷为：
Figure BDA0001250003980000133
式中Kne(j)表示滚动体与内外圈之间总的负荷-变形常数
Kne(j)＝6.24×104lwe 0.82Dw(j)0.11/(kcount*h)0.18[1+ci 0.9cos(αe-αi)]-1.11 (6)
式中
Figure BDA0001250003980000134
kcount表示接触的圆锥滚子的片数(δn(j)＞0时，该片圆锥滚子接触)。
第一列圆锥滚子轴承的力和力矩：
法向接触载荷为
Figure BDA0001250003980000135
法向接触载荷的轴向分力为
Figure BDA0001250003980000141
法向接触载荷的径向分力为
Figure BDA0001250003980000142
法向接触载荷轴向分力产生的转矩为
Figure BDA0001250003980000143
法向接触载荷径向分力产生的转矩为
Figure BDA0001250003980000144
Figure BDA0001250003980000145
为第一列圆锥滚子每片所受载荷。
第二列圆锥滚子轴承的力和力矩：
法向接触载荷为
Figure BDA0001250003980000146
法向接触载荷的轴向分力为
Figure BDA0001250003980000147
法向接触载荷的径向分力为
Figure BDA0001250003980000148
法向接触载荷轴向分力产生的转矩为
Figure BDA0001250003980000149
法向接触载荷径向分力产生的转矩为
Figure BDA00012500039800001410
为第二列圆锥滚子每片所受载荷。
交叉圆锥滚子轴承内圈在外力和滚子反作用力下处于平衡状态，力学平衡方程为：
Figure BDA0001250003980000151
Figure BDA0001250003980000152
Figure BDA0001250003980000153
把轴向力Fa＝141kN、Fr＝径向力178kN、M＝倾覆力矩170kN·m带入方程组式(17)、(18)、(19)里面，用Newton-Raphson迭代法进行运算求解，得出δa、δr、θ的值。
第三步：计算交叉圆锥滚子轴承的刚度矩阵：
第一列圆锥滚子的刚度矩阵：
Figure BDA0001250003980000154
第二列圆锥滚子的刚度矩阵：
Figure BDA0001250003980000161
交叉圆锥滚子轴承总的刚度矩阵为：
K＝K1+K2 (22)
把求解的δa、δr、θ带入到刚度矩阵公式(20)～(22)中，即可得到交叉圆锥滚子轴承的刚度。求得的刚度矩阵为：
式中：轴向刚度、径向刚度的单位为N/mm，角刚度的单位为N·mm/rad。
图5是当Fr＝0kN、M＝0kN·m时轴向主刚度Kaa与轴向力Fa的变化关系图；
图6是当Fa＝0kN、M＝0kN·m时径向主刚度Krr与径向力Fr的变化关系图；
图7是当Fa＝0kN、Fr＝0kN时角主刚度Kθθ与倾覆力矩M的变化关系图。

Claims (2)

1.一种交叉圆锥滚子轴承静刚度的快速提取方法，其特征在于：该方法包括以下步骤：
步骤一、建立交叉圆锥滚子轴承的力学模型，设置位置角
Figure FDA0002232975080000016
为零度且与其方向相同的滚子为第一列圆锥滚子，与第一列圆锥滚子交叉垂直的滚子为第二列圆锥滚子；
步骤二、将第一列圆锥滚子和第二列圆锥滚子沿其中心线方向切分成n片，得到第一列圆锥滚子每个切片上的法向变形量δn1(j)与内圈相对外圈在轴向力Fa作用下产生轴向位移δa、在径向力Fr作用下产生径向位移δr、在倾覆力矩M作用下产生倾角位移θ的关系式：
Figure FDA0002232975080000011
其中，dm1(j)为第一列圆锥滚子每个切片的中心节圆直径，
Figure FDA0002232975080000012
为滚子母线与其中心线的夹角，h为每片圆锥滚子的厚度，leff1(j)为第一列滚子每个切片的中心位置；
leff1(j)＝0.5h(1+2(j-1))
以及得到第二列圆锥滚子每个切片上的法向变形量δn2(j)与内圈相对外圈在轴向力Fa作用下产生轴向位移δa、在径向力Fr作用下产生径向位移δr、在倾覆力矩M作用下产生倾角位移θ的关系式：
Figure FDA0002232975080000013
其中，dm2(j)为第二列圆锥滚子每个切片的中心节圆直径，
Figure FDA0002232975080000014
为滚子母线与其中心线的夹角，h为每片圆锥滚子的厚度，leff2(j)为第二列滚子每个切片的中心位置；
leff2(j)＝0.5h(1+2(j-1))；
步骤三、根据Hertz线接触理论，由圆锥滚子每个切片上的法向变形量δn(j)和圆锥滚子每个切片上法向接触载荷qi(j)关系式，
Figure FDA0002232975080000015
Kne(j)＝6.24×104lwe 0.82Dw(j)0.11/(kcount*h)0.18[1+ci 0.9cos(αe-αi)]-1.11
式中，kcount表示接触的圆锥滚子的片数，δn(j)＞0时，表示该片圆锥滚子接触，lwe表示滚子有效接触长度，Dw(j)表示每片圆锥滚子的直径，αi为圆锥滚子与内圈的接触角，αf为圆锥滚子端部与外圈的接触角，αe为圆锥滚子母线与外圈的接触角；得到圆锥滚子的法向接触载荷为所有切片法向接触载荷qi(j)的和；
根据第一列圆锥滚子的法向接触载荷
Figure FDA0002232975080000022
得到在第一列圆锥滚子位置角处的轴向分力
Figure FDA0002232975080000024
径向分力
Figure FDA0002232975080000025
轴向分力产生的转矩径向分力产生的转矩
Figure FDA0002232975080000027
Figure FDA0002232975080000028
Figure FDA0002232975080000029
Figure FDA00022329750800000210
Figure FDA00022329750800000211
其中，lwe表示滚子有效接触长度，Dw(j)表示每片圆锥滚子的直径，leff1(j)为第一列滚子每个切片的中心位置；
同理，得到第二列圆锥滚子的法向接触载荷
Figure FDA00022329750800000212
在第二列圆锥滚子位置角
Figure FDA00022329750800000213
处的轴向分力径向分力
Figure FDA00022329750800000215
轴向分力产生的转矩
Figure FDA00022329750800000216
径向分力产生的转矩
Figure FDA00022329750800000217
Figure FDA00022329750800000218
Figure FDA00022329750800000219
Figure FDA00022329750800000220
Figure FDA00022329750800000221
其中，lwe表示滚子有效接触长度，Dw(j)表示每片圆锥滚子的直径，leff2(j)为第二列滚子每个切片的中心位置；
步骤四、轴承内圈在外部载荷和滚子反作用力下处于平衡状态，根据第一列圆锥滚子位置角
Figure FDA0002232975080000031
处的轴向分力
Figure FDA0002232975080000032
径向分力
Figure FDA0002232975080000033
轴向分力产生的转矩
Figure FDA0002232975080000034
径向分力产生的转矩
Figure FDA0002232975080000035
以及第二列圆锥滚子的法向接触载荷
Figure FDA0002232975080000036
在第二列圆锥滚子位置角
Figure FDA0002232975080000037
处的轴向分力
Figure FDA0002232975080000038
径向分力
Figure FDA0002232975080000039
轴向分力产生的转矩
Figure FDA00022329750800000310
径向分力产生的转矩
Figure FDA00022329750800000311
建立以内圈相对外圈在轴向力Fa作用下产生轴向位移δa、在径向力Fr作用下产生径向位移δr、在倾覆力矩M作用下产生倾角位移θ为未知数的内圈的三元力学平衡方程，所述的三元力学平衡方程为：
Figure FDA00022329750800000312
Figure FDA00022329750800000313
Figure FDA00022329750800000314
运用Newton-Raphson迭代法求解力学平衡方程得到δa、δr、θ的值；
步骤五、建立第一列圆锥滚子的刚度矩阵和第二列圆锥滚子的刚度矩阵；
所述的第一列圆锥滚子的刚度矩阵为
Figure FDA00022329750800000315
所述的第二列圆锥滚子的刚度矩阵为
Figure FDA0002232975080000041
将内圈相对外圈在轴向力Fa作用下产生轴向位移δa、在径向力Fr作用下产生径向位移δr、在倾覆力矩M作用下产生倾角位移θ的值代入第一列圆锥滚子的刚度矩阵和第二列圆锥滚子的刚度矩阵，第一列圆锥滚子的刚度矩阵和第二列圆锥滚子的刚度矩阵之和为交叉圆锥滚子轴承静刚度。
2.如权利要求1中所述的一种交叉圆锥滚子轴承静刚度的快速提取方法，其特征在于：所述的位置角
Figure FDA0002232975080000042
为第一列圆锥滚子的个数，z1＝z/2，z为圆锥滚子轴承所有滚子的个数，为第二列圆锥滚子的个数，z2＝z/2。
CN201710166615.5A 2017-03-20 2017-03-20 一种交叉圆锥滚子轴承静刚度的快速提取方法 Active CN106934165B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710166615.5A CN106934165B (zh) 2017-03-20 2017-03-20 一种交叉圆锥滚子轴承静刚度的快速提取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710166615.5A CN106934165B (zh) 2017-03-20 2017-03-20 一种交叉圆锥滚子轴承静刚度的快速提取方法

Publications (2)

Publication Number Publication Date
CN106934165A CN106934165A (zh) 2017-07-07
CN106934165B true CN106934165B (zh) 2020-02-14

Family

ID=59433403

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710166615.5A Active CN106934165B (zh) 2017-03-20 2017-03-20 一种交叉圆锥滚子轴承静刚度的快速提取方法

Country Status (1)

Country Link
CN (1) CN106934165B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107704652B (zh) * 2017-08-29 2021-05-04 明阳智慧能源集团股份公司 一种风力发电机组轴承刚度计算工具
CN112131769B (zh) * 2020-09-14 2023-12-12 东南大学 交叉滚柱直线导轨五自由度静刚度计算方法
CN113392544B (zh) * 2021-05-28 2022-08-26 东北林业大学 一种基于变形协调理论的行星螺纹滚柱轴承接触载荷计算方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104732068A (zh) * 2015-02-27 2015-06-24 河南科技大学 通过力学建模设置交叉圆柱滚子转盘轴承刚度矩阵的方法
CN106021790A (zh) * 2016-06-01 2016-10-12 河南科技大学 球-柱联合转台轴承动、静承载能力计算及曲线绘制方法
CN106066915A (zh) * 2016-06-01 2016-11-02 河南科技大学 一种获取球‑柱联合转盘轴承刚度矩阵的方法及装置
CN106089993A (zh) * 2016-08-27 2016-11-09 河南科技大学 一种获取球锥混合双列轮毂轴承载荷分布的方法及装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104732068A (zh) * 2015-02-27 2015-06-24 河南科技大学 通过力学建模设置交叉圆柱滚子转盘轴承刚度矩阵的方法
CN106021790A (zh) * 2016-06-01 2016-10-12 河南科技大学 球-柱联合转台轴承动、静承载能力计算及曲线绘制方法
CN106066915A (zh) * 2016-06-01 2016-11-02 河南科技大学 一种获取球‑柱联合转盘轴承刚度矩阵的方法及装置
CN106089993A (zh) * 2016-08-27 2016-11-09 河南科技大学 一种获取球锥混合双列轮毂轴承载荷分布的方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
一种交叉圆柱滚子转盘轴承的滚子偏载分析及修形;李云峰等;《机械工程学报》;20150131;第161-166页 *
精密交叉圆柱滚子转盘轴承刚性分析;程亚兵;《中国优秀硕士论文全文数据库》;2015021;第12-17页，第37页 *

Also Published As

Publication number Publication date
CN106934165A (zh) 2017-07-07

Similar Documents

Publication Publication Date Title
CN107066699B (zh) 一种获取交叉圆锥滚子轴承载荷分布的方法
CN106934165B (zh) 一种交叉圆锥滚子轴承静刚度的快速提取方法
CN107704652B (zh) 一种风力发电机组轴承刚度计算工具
CN108266459B (zh) 基于轴承磨配间隙的机床主轴径向跳动计算方法
JP6115123B2 (ja) アンギュラ玉軸受の組立装置、及びアンギュラ玉軸受の組立方法
CN102819635B (zh) 特大型双排四点接触球轴承静力学模型建立方法
CN107153734B (zh) 一种定位预紧下角接触球轴承非线性刚度的计算系统及方法
CN104008274A (zh) 一种预测高速主轴系统切削稳定性区域的方法
JP2014046389A (ja) 旋削方法および旋削装置
CN107563081B (zh) 高速空心圆柱滚子轴承的参数优化方法
CN104220767A (zh) 止推凸缘及止推板组合件
CN116822283A (zh) 考虑推进器激励的推进轴系轴承动态载荷计算方法
Cui et al. Effect of roller dynamic unbalance on cage stress of high-speed cylindrical roller bearing
CN108984936A (zh) 高速双联滚动轴承电主轴转子系统动态设计方法
CN112131769B (zh) 交叉滚柱直线导轨五自由度静刚度计算方法
CN116579104A (zh) 内外圈同时转动的中介轴承内部动态载荷获取方法及系统
CN116011258A (zh) 一种考虑套圈变形的高速球轴承动态特性分析方法
CN112231854B (zh) 双v密珠型直线导轨五自由度静刚度计算方法
CN109580223B (zh) 一种考虑芯轴变形影响的组合式轴连轴承载荷分析方法
Wang et al. Effect of inner ring centrifugal displacement on the dynamic characteristics of high-speed angular contact ball bearing
CN113946919A (zh) 一种含组合角不对中的深沟球轴承准静态解析模型分析方法
CN108372264A (zh) 一种有利于超大环辗轧稳定成形的四个抱辊布局的方法
Li et al. Study on nonlinear stiffness of rolling ball bearing under varied operating conditions
JP2010054256A (ja) 転がり軸受ユニット用荷重測定装置
CN105004293A (zh) 一种轮毂轴承套圈零件沟道轴向综合位置变差测量方法

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
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20170707

Assignee: LUOYANG HONGYUAN BEARING TECHNOLOGY Co.,Ltd.

Assignor: HENAN University OF SCIENCE AND TECHNOLOGY

Contract record no.: X2022980003705

Denomination of invention: A fast extraction method of static stiffness of cross tapered roller bearing

Granted publication date: 20200214

License type: Common License

Record date: 20220402

EE01 Entry into force of recordation of patent licensing contract