CN107066699A - 一种获取交叉圆锥滚子轴承载荷分布的方法 - Google Patents
一种获取交叉圆锥滚子轴承载荷分布的方法 Download PDFInfo
- Publication number
- CN107066699A CN107066699A CN201710166616.XA CN201710166616A CN107066699A CN 107066699 A CN107066699 A CN 107066699A CN 201710166616 A CN201710166616 A CN 201710166616A CN 107066699 A CN107066699 A CN 107066699A
- Authority
- CN
- China
- Prior art keywords
- row
- mrow
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000009826 distribution Methods 0.000 title claims abstract description 40
- 238000006073 displacement reaction Methods 0.000 claims description 53
- 230000014509 gene expression Effects 0.000 claims description 10
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 241000969729 Apteryx rowi Species 0.000 claims 1
- 238000005096 rolling process Methods 0.000 abstract description 8
- 238000005520 cutting process Methods 0.000 abstract description 3
- 238000013178 mathematical model Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 6
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 230000003068 static effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000004753 textile Substances 0.000 description 1
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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Rolling Contact Bearings (AREA)
Abstract
一种获取交叉圆锥滚子轴承载荷分布的方法，本发明首先基于Hertz接触理论和滚动轴承设计方法建立交叉圆锥滚子轴承的精确数学模型，分出第一列圆锥滚子和第二列圆锥滚子进行单独计算，采用切片单独计算载荷计算法向接触载荷的方法，利用法向接触载荷得到轴向分力、径向分力，轴向分力产生的转矩、径向分力产生的转矩得出三元力平衡方程，并采用数值方法对力学平衡方程进行求解，进而获取交叉圆锥滚子轴承载荷分布。本文给出的获取交叉圆锥滚子轴承载荷分布的方法，解决了目前鲜有研究的交叉圆锥滚子载荷分布的问题。为交叉圆锥滚子轴承各参数的合理设计提供依据，提高了设计的准确性与效率。
Description
技术领域
本发明涉及交叉圆锥滚子轴承领域，具体说的是一种获取交叉圆锥滚子轴承载荷分布的方法。
背景技术
交叉滚子轴承具有高旋转精度、高刚性、高转速、低转动惯量、低扭矩、低摩擦、复合承载(轴向载荷、径向载荷、倾覆力矩组合)、设计紧凑、节省空间、操作安装简化等优点，广泛应用在精密旋转工作台、立式车床、卧式车床、立式磨床、大型滚齿机、旋转转向架、加工中心转台、工业机器人、医疗设备、测量装置、IC卡制造设备和纺织机械等领域。由此可见，交叉滚子轴承具有很好的市场和前景。交叉滚子轴承分为两大类，一类是圆柱交叉滚子轴承，另一类是圆锥交叉滚子轴承。目前，国内外学者对交叉圆柱滚子轴承有所研究，但对交叉圆锥滚子轴承鲜有研究，对交叉圆锥滚子轴承的载荷分布也没有给出具体的研究方法。因此开展交叉圆锥滚子轴承载荷分布的研究具有重要的理论研究意义。
发明内容
为解决上述技术问题，本发明提供一种获取交叉圆锥滚子轴承载荷分布的方法，能够得到交叉圆锥滚子轴承载荷精确分布。
为实现上述技术目的，所采用的技术方案是：一种获取交叉圆锥滚子轴承载荷分布的方法，该方法包括以下步骤：
步骤一、建立交叉圆锥滚子轴承的力学模型，设置位置角为零度且与其方向相同的滚子为第一列圆锥滚子，与第一列圆锥滚子交叉垂直的滚子为第二列圆锥滚子；
步骤二、将第一列圆锥滚子和第二列圆锥滚子沿其中心线方向切分成n片，以第一列圆锥滚子每个切片上的轴向变型量δa1(j)和每个切片上的径向变形量δr1(j)，得到第一列圆锥滚子每个切片上的法向变形量δn1(j)与内圈相对外圈在轴向力Fa作用下产生轴向位移δa、在径向力Fr作用下产生径向位移δr、在倾覆力矩M作用下产生倾角位移θ的关系式，以第二列圆锥滚子每个切片上的轴向变型量δa2(j)和每个切片上的径向变形量δr2(j)，得到第二列圆锥滚子每个切片上的法向变形量δn2(j)与内圈相对外圈在轴向力Fa作用下产生轴向位移δa、在径向力Fr作用下产生径向位移δr、在倾覆力矩M作用下产生倾角位移θ的关系式；
步骤三、根据Hertz线接触理论，由圆锥滚子每个切片上的法向变形量δn(j)和圆锥滚子每个切片上法向接触载荷qi(j)关系式，得到圆锥滚子的法向接触载荷为所有切片法向接触载荷qi(j)的和，根据第一列圆锥滚子的法向接触载荷得到在第一列圆锥滚子位置角处的轴向分力径向分力轴向分力产生的转矩径向分力产生的转矩同理，得到第二列圆锥滚子的法向接触载荷在第二列圆锥滚子位置角处的轴向分力径向分力轴向分力产生的转矩径向分力产生的转矩
步骤四、根据轴承内圈在外部载荷和滚子反作用力下处于平衡状态，建立以内圈相对外圈在轴向力Fa作用下产生轴向位移δa、在径向力Fr作用下产生径向位移δr、在倾覆力矩M作用下产生倾角位移θ为未知数的内圈的三元力学平衡方程，运用Newton-Raphson迭代法求解力学平衡方程得到δa、δr、θ的值；
步骤五、将δa、δr、θ的值代入步骤二，得到第一列圆锥滚子每个切片上的法向变形量δn1(j)和第二列圆锥滚子每个切片上的法向变形量δn2(j)，再将第一列圆锥滚子每个切片上的法向变形量δn1(j)和第二列圆锥滚子每个切片上的法向变形量δn2(j)代入步骤三，得到交叉圆锥滚子轴承两列滚子的载荷分布。
本发明所述的每个切片上的法向变形量δn1(j)与内圈相对外圈在轴向力Fa作用下产生轴向位移δa、在径向力Fr作用下产生径向位移δr、在倾覆力矩M作用下产生倾角位移θ的关系式为
δn1(j)＝δa1sinαe+δr1cosαe
其中，dm1(j)为第一列圆锥滚子每个切片的中心节圆直径，为滚子母线与其中心线的夹角，h为每片圆锥滚子的厚度，leff1(j)为第一列滚子每个切片的中心位置；
leff1(j)＝0.5h(1+2(j-1))。
本发明所述的第二列圆锥滚子每个切片上的法向变形量δn2(j)与内圈相对外圈在轴向力Fa作用下产生轴向位移δa、在径向力Fr作用下产生径向位移δr、在倾覆力矩M作用下产生倾角位移θ的关系式为
δn2(j)＝-δa2sinαe+δr2cosαe
其中，dm2(j)为第二列圆锥滚子每个切片的中心节圆直径，为滚子母线与其中心线的夹角，h为每片圆锥滚子的厚度，leff2(j)为第二列滚子每个切片的中心位置；
leff2(j)＝0.5h(1+2(j-1))。
本发明所述的圆锥滚子每个切片上的法向变形量δn(j)和圆锥滚子每个切片上法向接触载荷qi(j)关系式为
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为圆锥滚子母线与外圈的接触角。
本发明所述的第一列圆锥滚子的法向接触载荷的表达式为
其中，为第一列圆锥滚子每个切处所受法向接触载荷。
本发明所述的第一列圆锥滚子的法向接触载荷得到在第一列圆锥滚子位置角处的轴向分力径向分力轴向分力产生的转矩径向分力产生的转矩的表达式为
其中，lwe表示滚子有效接触长度，Dw(j)表示每片圆锥滚子的直径，leff1(j)为第一列滚子每个切片的中心位置；
leff1(j)＝0.5h(1+2(j-1))。
本发明所述的第二列圆锥滚子的法向接触载荷的表达式为
其中，为第二列圆锥滚子每个切片所受法向接触载荷。
本发明所述的第二列圆锥滚子的法向接触载荷在第二列圆锥滚子位置角处的轴向分力径向分力轴向分力产生的转矩径向分力产生的转矩的表达式为
其中，lwe表示滚子有效接触长度，Dw(j)表示每片圆锥滚子的直径，leff2(j)为第二列滚子每个切片的中心位置；
leff2(j)＝0.5h(1+2(j-1))。
本发明所述的三元力学平衡方程为
本发明所述的位置角i1＝1,2,3,…,z1,z1为第一列圆锥滚子的个数，z1＝z/2，z为圆锥滚子轴承所有滚子的个数，i2＝1,2,3,…,z2，,z2为第二列圆锥滚子的个数，z2＝z/2。
本发明的有益效果是：本发明首先基于Hertz接触理论和滚动轴承设计方法建立交叉圆锥滚子轴承的精确数学模型，分出第一列圆锥滚子和第二列圆锥滚子进行单独计算，采用切片单独计算载荷计算法向接触载荷的方法，并采用数值方法对力学平衡方程进行求解，进而获取交叉圆锥滚子轴承载荷分布。本发明给出的获取交叉圆锥滚子轴承载荷分布的方法，解决了目前鲜有研究的交叉圆锥滚子载荷分布的问题。为交叉圆锥滚子轴承各参数的合理设计提供依据，提高了设计的准确性与效率。
附图说明
图1为本发明交叉圆锥滚子轴承受力与内圈位移图；
图2为本发明交叉圆锥滚子轴承滚子位置角示意图；
图3为本发明交叉圆锥滚子轴承接触角示意图；
图4为本发明基于切片法的圆锥滚子变形和受载关系示意图；
图5为本发明交叉圆锥滚子轴承内圈受力简图；
图6为不同修型交叉圆锥滚子轴承载荷分布图；
图7为不同凸型最大受载圆锥滚子每片接触半宽沿母线分布；
图8为不同凸型最大受载圆锥滚子每片接触压力沿母线分布。
具体实施方式
本发明首先基于Hertz接触理论和滚动轴承设计方法建立交叉圆锥滚子轴承的精确数学模型，并采用数值方法对力学平衡方程进行求解，进而获取交叉圆锥滚子轴承载荷分布。本发明给出的获取交叉圆锥滚子轴承载荷分布的方法，解决了目前鲜有研究的交叉圆锥滚子载荷分布的问题。
获取交叉圆锥滚子轴承载荷分布的方法，包含以下步骤：
(1)建立交叉圆锥滚子轴承的力学模型：
假设轴承外圈固定，按照静力学的方法来建立交叉圆锥滚子轴承的力学模型，模型如图1，Fa、Fr、M分别表示交叉圆锥滚子轴承承受的轴向力、径向力和倾覆力矩，δa为轴向位移，δr为径向位移，θ为倾角位移。交叉圆锥滚子轴承有两列相互垂直的圆锥滚子及内外圈组成，记位置角为零度且与其方向相同的滚子为第一列圆锥滚子，为第一列圆锥滚子的位置角，(i1＝1,2,3,…,z1),z1为第一列圆锥滚子的个数，z1＝z/2，z为圆锥滚子轴承所有滚子的个数；记与第一列圆锥滚子交叉垂直的滚子为第2列圆锥滚子，为第二列圆锥滚子的位置角，(i2＝1,2,3,…,z2),z2为第二列圆锥滚子的个数，z2＝z/2。位置角关系定义如图2所示。
(2)计算交叉圆锥滚子轴承的变形量：
内圈相对外圈在轴向力Fa作用下产生轴向位移δa，在径向力Fr作用下产生径向位移δr，在倾覆力矩M作用下产生倾角位移θ。记轴向位移引起的轴向变形为δa，径向位移引起的径向变形分量为倾角位移引起的轴向变形分量为(dm(j)代表每个切片的中心节圆直径)，倾角位移引起的径向变形分量为(lwe表示滚子有效接触长度，为滚子母线与其中心线的夹角，leff(j)为滚子每个切片的中心位置，β为内圈滚道母线与其中心线之间的夹角)。内圈相对外圈位移减小的力的方向为正，位移增加的力的方向为负。接触角定义如图3所示，αi为圆锥滚子与内圈的接触角，αf为圆锥滚子端部与外圈的接触角，αe为圆锥滚子母线与外圈的接触角。由于圆锥滚子沿母线接触位置处的直径不相同，因此产生的法向变形量也不同。为了更精确的分析圆锥滚子的载荷分布，用切片法对圆锥滚子接触变形关系进行研究分析。切片法示意图如图4，把圆锥滚子沿其中心线方向切分成n片，qi(j)代表每个切片上的法向载荷，dm(j)表示圆锥滚子每个切片的中心节圆直径，Dw(j)表示每片圆锥滚子的直径。
第一列圆锥滚子轴承变形量：
每个切片上的轴向变形量为
每个切片上的径向变形量为
每个切片上的法向变形量为
δn1(j)＝δa1sinαe+δr1cosαe (3)
式中dm1(j)为第一列圆锥滚子每个切片的中心节圆直径，h表示每片圆锥滚子的厚度，leff1(j)为第一列滚子每个切片的中心位置。
leff1(j)＝0.5h(1+2(j-1)) (4)
第二列圆锥滚子轴承变形量：
每个切片上的轴向变形量为
每个切片上的径向变形量为
每个切片上的法向变形量为
δn2(j)＝-δa2sinαe+δr2cosαe (7)
dm2(j)为第二列圆锥滚子每个切片的中心节圆直径，leff2(j)为第二列滚子每个切片的中心位置。
leff2(j)＝0.5h(1+2(j-1)) (8)
(3)计算交叉圆锥滚子轴承的力和力矩分量：
根据Hertz线接触理论，由圆锥滚子轴承法向变形量和法向接触载荷的关系，计算每片圆锥滚子在位置角处的法向接触载荷，每片圆锥滚子法向接触载荷为：
式中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)
式中kcount表示接触的圆锥滚子的片数(δn(j)＞0时，表示该片圆锥滚子接触)。
第一列圆锥滚子轴承的力和力矩：
法向接触载荷为
法向接触载荷的轴向分力为
法向接触载荷的径向分力为
法向接触载荷轴向分力产生的转矩为
法向接触载荷径向分力产生的转矩为
为第一列圆锥滚子每片所受载荷。
第二列圆锥滚子轴承的力和力矩：
法向接触载荷为
法向接触载荷的轴向分力为
法向接触载荷的径向分力为
法向接触载荷轴向分力产生的转矩为
法向接触载荷径向分力产生的转矩为
为第二列圆锥滚子每片所受载荷。
(4)建立交叉圆锥滚子轴承力学平衡方程组：
交叉圆锥滚子轴承内圈在外力和滚子反作用力下处于平衡状态，受力简图如图5，规定内圈受到向下的径向力为正，向左的轴向力为正，逆时针的力矩为正。内圈的力学平衡方程为：
(5)获取交叉圆锥滚子轴承的载荷分布：
式(21)-(23)构成了以δa、δr、θ为未知量的三元非线性方程组，可运用Newton-Raphson迭代法进行数值运算求解，得出δa、δr、θ的值。然后，分别带入到式(1)-(3)和(5)-(7)计算两列交叉圆锥滚子轴承的法向变形量，再将法向变形量带入到式(11)和(16)，即可获取交叉圆锥滚子轴承两列的载荷分布。
下面以某一具体的交叉圆锥滚子轴承为例进行说明：
某型号交叉圆锥滚子轴承，结构参数为：轴承的节圆直径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、径向力50kN、倾覆力矩110kN·m。
第一步：假设轴承外圈固定，按照静力学的方法来建立交叉圆锥滚子轴承的力学模型，模型如图1，Fa、Fr、M分别表示交叉圆锥滚子轴承承受的轴向力、径向力和倾覆力矩，δa为轴向位移，δr为径向位移，θ为倾角位移。交叉圆锥滚子轴承有两列相互垂直的圆锥滚子及内外圈组成，记位置角为零度且与其方向相同的滚子为第一列圆锥滚子，为第一列圆锥滚子的位置角，(i1＝1,2,3,…,z1),z1为第一列圆锥滚子的个数，z1＝z/2，z为圆锥滚子轴承所有滚子的个数；记与第一列圆锥滚子交叉垂直的滚子为第2列圆锥滚子，为第二列圆锥滚子的位置角，(i2＝1,2,3,…,z2),z2为第二列圆锥滚子的个数，z2＝z/2。位置角关系定义如图2所示。
第二步：内圈相对外圈在轴向力Fa作用下产生轴向位移δa，在径向力Fr作用下产生径向位移δr，在倾覆力矩M作用下产生倾角位移θ。记轴向位移引起的轴向变形为δa，径向位移引起的径向变形分量为倾角位移引起的轴向变形分量为(dm(j)代表每个切片的中心节圆直径)，倾角位移引起的径向变形分量为(lwe表示滚子有效接触长度，为滚子母线与其中心线的夹角，leff(j)为滚子每个切片的中心位置，β为内圈滚道母线与其中心线之间的夹角)。内圈相对外圈位移减小的力的方向为正，位移增加的力的方向为负。接触角定义如图3所示，αi为圆锥滚子与内圈的接触角，αf为圆锥滚子端部与外圈的接触角，αe为圆锥滚子母线与外圈的接触角。由于圆锥滚子沿母线接触位置处的直径不相同，因此产生的法向变形量也不同。为了更精确的分析圆锥滚子的载荷分布，用切片法对圆锥滚子接触变形关系进行研究分析。切片法示意图如图4，把圆锥滚子沿其母线方向切分成n片，qi(j)代表每个切片上的法向载荷。dm(j)表示圆锥滚子每个切片的中心节圆直径，Dw(j)表示每片圆锥滚子的直径。
第一列圆锥滚子轴承变形量：
每个切片上的轴向变形量为
每个切片上的径向变形量为
每个切片上的法向变形量为
δn1(j)＝δa1sinαe+δr1cosαe (3)
式中dm1(j)为第一列圆锥滚子每个切片的中心节圆直径，h表示每片圆锥滚
子的厚度，leff1(j)为第一列滚子每个切片的中心位置。
leff1(j)＝0.5h(1+2(j-1)) (4)
第二列圆锥滚子轴承变形量：
每个切片上的轴向变形量为
每个切片上的径向变形量为
每个切片上的法向变形量为
δn2(j)＝-δa2sinαe+δr2cosαe (7)
dm2(j)为第二列圆锥滚子每个切片的中心节圆直径，leff2(j)为第二列滚子每个切片的中心位置。
leff2(j)＝0.5h(1+2(j-1)) (8)
第三步：根据Hertz线接触理论，由圆锥滚子轴承法向变形量和法向接触载荷的关系，计算每片圆锥滚子在位置角处的法向接触载荷，每片圆锥滚子法向接触载荷为：
式中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)
式中kcount表示接触的圆锥滚子的片数(δn(j)＞0时，表示该片圆锥滚子接触)。
第一列圆锥滚子轴承的力和力矩：
法向接触载荷为
法向接触载荷的轴向分力为
法向接触载荷的径向分力为
法向接触载荷轴向分力产生的转矩为
法向接触载荷径向分力产生的转矩为
第二列圆锥滚子轴承的力和力矩：
法向接触载荷为
法向接触载荷的轴向分力为
法向接触载荷的径向分力为
法向接触载荷轴向分力产生的转矩为
法向接触载荷径向分力产生的转矩为
第四步：交叉圆锥滚子轴承内圈在外力和滚子反作用力下处于平衡状态。受力简图如图5，规定内圈受到向下的径向力为正，向左的轴向力为正，逆时针的力矩为正。内圈的力学平衡方程为：
其中：轴向力Fa＝141kN、径向力Fr＝50kN、倾覆力矩M＝110kN·m。
第五步：式(21)-(23)构成了以δa、δr、θ为未知量的三元方程组，运用Newton-Raphson迭代法进行运算求解，得出δa、δr、θ的值，分别带入到式(1)-(3)和(5)-(7)计算两列交叉圆锥滚子轴承的法向变形量，再将法向变形量带入到式(11)和(16)，即可获取交叉圆锥滚子轴承两列的载荷分布。
在已知外部载荷条件下，分别计算滚子母线无修型、圆弧修型、对数修形3种情况下交叉圆锥滚子的载荷分布，如图6。图7是不同凸型最大受载圆锥滚子每片接触半宽沿母线分布图，图8是不同凸型最大受载圆锥滚子每片接触应力沿母线分布图。
Claims (10)
1.一种获取交叉圆锥滚子轴承载荷分布的方法，其特征在于：该方法包括以下步骤：
步骤一、建立交叉圆锥滚子轴承的力学模型，设置位置角为零度且与其方向相同的滚子为第一列圆锥滚子，与第一列圆锥滚子交叉垂直的滚子为第二列圆锥滚子；
步骤二、将第一列圆锥滚子和第二列圆锥滚子沿其中心线方向切分成n片，以第一列圆锥滚子每个切片上的轴向变型量δa1(j)和每个切片上的径向变形量δr1(j)，得到第一列圆锥滚子每个切片上的法向变形量δn1(j)与内圈相对外圈在轴向力Fa作用下产生轴向位移δa、在径向力Fr作用下产生径向位移δr、在倾覆力矩M作用下产生倾角位移θ的关系式，以第二列圆锥滚子每个切片上的轴向变型量δa2(j)和每个切片上的径向变形量δr2(j)，得到第二列圆锥滚子每个切片上的法向变形量δn2(j)与内圈相对外圈在轴向力Fa作用下产生轴向位移δa、在径向力Fr作用下产生径向位移δr、在倾覆力矩M作用下产生倾角位移θ的关系式；
步骤三、根据Hertz线接触理论，由圆锥滚子每个切片上的法向变形量δn(j)和圆锥滚子每个切片上法向接触载荷qi(j)关系式，得到圆锥滚子的法向接触载荷为所有切片法向接触载荷qi(j)的和，根据第一列圆锥滚子的法向接触载荷得到在第一列圆锥滚子位置角处的轴向分力径向分力轴向分力产生的转矩径向分力产生的转矩同理，得到第二列圆锥滚子的法向接触载荷在第二列圆锥滚子位置角处的轴向分力径向分力轴向分力产生的转矩径向分力产生的转矩
步骤四、根据轴承内圈在外部载荷和滚子反作用力下处于平衡状态，建立以内圈相对外圈在轴向力Fa作用下产生轴向位移δa、在径向力Fr作用下产生径向位移δr、在倾覆力矩M作用下产生倾角位移θ为未知数的内圈的三元力学平衡方程，运用Newton-Raphson迭代法求解力学平衡方程得到δa、δr、θ的值；
步骤五、将δa、δr、θ的值代入步骤二，得到第一列圆锥滚子每个切片上的法向变形量δn1(j)和第二列圆锥滚子每个切片上的法向变形量δn2(j)，再将第一列圆锥滚子每个切片上的法向变形量δn1(j)和第二列圆锥滚子每个切片上的法向变形量δn2(j)代入步骤三，得到交叉圆锥滚子轴承两列滚子的载荷分布。
2.如权利要求1所述的一种获取交叉圆锥滚子轴承载荷分布的方法，其特征在于：所述的每个切片上的法向变形量δn1(j)与内圈相对外圈在轴向力Fa作用下产生轴向位移δa、在径向力Fr作用下产生径向位移δr、在倾覆力矩M作用下产生倾角位移θ的关系式为
δn1(j)＝δa1sinαe+δr1cosαe
其中，dm1(j)为第一列圆锥滚子每个切片的中心节圆直径，为滚子母线与其中心线的夹角，h为每片圆锥滚子的厚度，leff1(j)为第一列滚子每个切片的中心位置；
leff1(j)＝0.5h(1+2(j-1))。
3.如权利要求1所述的一种获取交叉圆锥滚子轴承载荷分布的方法，其特征在于：所述的第二列圆锥滚子每个切片上的法向变形量δn2(j)与内圈相对外圈在轴向力Fa作用下产生轴向位移δa、在径向力Fr作用下产生径向位移δr、在倾覆力矩M作用下产生倾角位移θ的关系式为
δn2(j)＝-δa2sinαe+δr2cosαe
其中，dm2(j)为第二列圆锥滚子每个切片的中心节圆直径，为滚子母线与其中心线的夹角，h为每片圆锥滚子的厚度，leff2(j)为第二列滚子每个切片的中心位置；
leff2(j)＝0.5h(1+2(j-1))。
4.如权利要求1所述的一种获取交叉圆锥滚子轴承载荷分布的方法，其特征在于：所述的圆锥滚子每个切片上的法向变形量δn(j)和圆锥滚子每个切片上法向接触载荷qi(j)关系式为
<mrow>
<msub>
<mi>q</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>h</mi>
<mo>*</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>n</mi>
<mi>e</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>&delta;</mi>
<mi>n</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mn>1.11</mn>
</msup>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>&ForAll;</mo>
<msub>
<mi>&delta;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>&GreaterEqual;</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>&ForAll;</mo>
<msub>
<mi>&delta;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo><</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
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为圆锥滚子母线与外圈的接触角。
5.如权利要求4所述的一种获取交叉圆锥滚子轴承载荷分布的方法，其特征在于：所述的第一列圆锥滚子的法向接触载荷的表达式为
其中，为第一列圆锥滚子每个切处所受法向接触载荷。
6.如权利要求5所述的一种获取交叉圆锥滚子轴承载荷分布的方法，其特征在于：所述的第一列圆锥滚子的法向接触载荷得到在第一列圆锥滚子位置角处的轴向分力径向分力轴向分力产生的转矩径向分力产生的转矩的表达式为
<mrow>
<msubsup>
<mi>Q</mi>
<mi>a</mi>
<mn>1</mn>
</msubsup>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msubsup>
<mi>q</mi>
<mi>i</mi>
<mn>1</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>sin&alpha;</mi>
<mi>e</mi>
</msub>
</mrow>
其中，lwe表示滚子有效接触长度，Dw(j)表示每片圆锥滚子的直径，leff1(j)为第一列滚子每个切片的中心位置；
leff1(j)＝0.5h(1+2(j-1))。
7.如权利要求4所述的一种获取交叉圆锥滚子轴承载荷分布的方法，其特征在于：所述的第二列圆锥滚子的法向接触载荷的表达式为
其中，为第二列圆锥滚子每个切片所受法向接触载荷。
8.如权利要求4所述的一种获取交叉圆锥滚子轴承载荷分布的方法，其特征在于：所述的第二列圆锥滚子的法向接触载荷在第二列圆锥滚子位置角处的轴向分力径向分力轴向分力产生的转矩径向分力产生的转矩的表达式为
<mrow>
<msubsup>
<mi>Q</mi>
<mi>a</mi>
<mn>2</mn>
</msubsup>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msubsup>
<mi>q</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>sin&alpha;</mi>
<mi>e</mi>
</msub>
</mrow>
其中，lwe表示滚子有效接触长度，Dw(j)表示每片圆锥滚子的直径，leff2(j)为第二列滚子每个切片的中心位置；
leff2(j)＝0.5h(1+2(j-1))。
9.如权利要求1所述的一种获取交叉圆锥滚子轴承载荷分布的方法，其特征在于：所述的三元力学平衡方程为
10.如权利要求1-9中任意一项所述的一种获取交叉圆锥滚子轴承载荷分布的方法，其特征在于：所述的位置角i1＝1,2,3,…,z1,z1为第一列圆锥滚子的个数，z1＝z/2，z为圆锥滚子轴承所有滚子的个数，i2＝1,2,3,…,z2，,z2为第二列圆锥滚子的个数，z2＝z/2。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710166616.XA CN107066699B (zh) | 2017-03-20 | 2017-03-20 | 一种获取交叉圆锥滚子轴承载荷分布的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710166616.XA CN107066699B (zh) | 2017-03-20 | 2017-03-20 | 一种获取交叉圆锥滚子轴承载荷分布的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107066699A true CN107066699A (zh) | 2017-08-18 |
| CN107066699B CN107066699B (zh) | 2020-02-21 |
Family
ID=59620176
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710166616.XA Active CN107066699B (zh) | 2017-03-20 | 2017-03-20 | 一种获取交叉圆锥滚子轴承载荷分布的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107066699B (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107391893A (zh) * | 2017-09-12 | 2017-11-24 | 大连交通大学 | 基于空心圆柱滚子接触变形量和载荷分布的轴承参数优化方法 |
| CN107704652A (zh) * | 2017-08-29 | 2018-02-16 | 明阳智慧能源集团股份公司 | 一种风力发电机组轴承刚度计算工具 |
| CN108050153A (zh) * | 2018-01-05 | 2018-05-18 | 成都市青白江四方非标轴承有限公司 | 一种交叉圆锥滚子轴承 |
| CN109933889A (zh) * | 2019-03-11 | 2019-06-25 | 中车大连机车研究所有限公司 | 一种铁路轴承用滚动体多曲线段修型方法 |
| CN114357621A (zh) * | 2022-01-04 | 2022-04-15 | 北京机电工程研究所 | 飞行器活动部件展开过程支反力载荷计算方法及系统 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120247315A1 (en) * | 2011-02-01 | 2012-10-04 | Corvid Technologies, Inc. | Blast Deflecting Shield for Ground Vehicles and Shielded Ground Vehicles and Methods Including Same |
| CN104732068A (zh) * | 2015-02-27 | 2015-06-24 | 河南科技大学 | 通过力学建模设置交叉圆柱滚子转盘轴承刚度矩阵的方法 |
| CN106066915A (zh) * | 2016-06-01 | 2016-11-02 | 河南科技大学 | 一种获取球‑柱联合转盘轴承刚度矩阵的方法及装置 |
| CN106089993A (zh) * | 2016-08-27 | 2016-11-09 | 河南科技大学 | 一种获取球锥混合双列轮毂轴承载荷分布的方法及装置 |
-
2017
- 2017-03-20 CN CN201710166616.XA patent/CN107066699B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120247315A1 (en) * | 2011-02-01 | 2012-10-04 | Corvid Technologies, Inc. | Blast Deflecting Shield for Ground Vehicles and Shielded Ground Vehicles and Methods Including Same |
| CN104732068A (zh) * | 2015-02-27 | 2015-06-24 | 河南科技大学 | 通过力学建模设置交叉圆柱滚子转盘轴承刚度矩阵的方法 |
| CN106066915A (zh) * | 2016-06-01 | 2016-11-02 | 河南科技大学 | 一种获取球‑柱联合转盘轴承刚度矩阵的方法及装置 |
| CN106089993A (zh) * | 2016-08-27 | 2016-11-09 | 河南科技大学 | 一种获取球锥混合双列轮毂轴承载荷分布的方法及装置 |
Non-Patent Citations (1)
| Title |
|---|
| 李云峰等: "交叉圆柱滚子转盘轴承的滚子偏载分析及修形", 《机械工程学报》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107704652A (zh) * | 2017-08-29 | 2018-02-16 | 明阳智慧能源集团股份公司 | 一种风力发电机组轴承刚度计算工具 |
| CN107704652B (zh) * | 2017-08-29 | 2021-05-04 | 明阳智慧能源集团股份公司 | 一种风力发电机组轴承刚度计算工具 |
| CN107391893A (zh) * | 2017-09-12 | 2017-11-24 | 大连交通大学 | 基于空心圆柱滚子接触变形量和载荷分布的轴承参数优化方法 |
| CN107391893B (zh) * | 2017-09-12 | 2020-09-29 | 大连交通大学 | 基于空心圆柱滚子接触变形量和载荷分布的轴承参数优化方法 |
| CN108050153A (zh) * | 2018-01-05 | 2018-05-18 | 成都市青白江四方非标轴承有限公司 | 一种交叉圆锥滚子轴承 |
| CN109933889A (zh) * | 2019-03-11 | 2019-06-25 | 中车大连机车研究所有限公司 | 一种铁路轴承用滚动体多曲线段修型方法 |
| CN109933889B (zh) * | 2019-03-11 | 2023-05-05 | 中车大连机车研究所有限公司 | 一种铁路轴承用滚动体多曲线段修型方法 |
| CN114357621A (zh) * | 2022-01-04 | 2022-04-15 | 北京机电工程研究所 | 飞行器活动部件展开过程支反力载荷计算方法及系统 |
| CN114357621B (zh) * | 2022-01-04 | 2024-06-11 | 北京机电工程研究所 | 飞行器活动部件展开过程支反力载荷计算方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107066699B (zh) | 2020-02-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107066699B (zh) | 一种获取交叉圆锥滚子轴承载荷分布的方法 | |
| CN107704652B (zh) | 一种风力发电机组轴承刚度计算工具 | |
| CN103823940B (zh) | 一种基于整体动力学模型的静压转台动态响应计算方法 | |
| CN102819635B (zh) | 特大型双排四点接触球轴承静力学模型建立方法 | |
| CN106934165B (zh) | 一种交叉圆锥滚子轴承静刚度的快速提取方法 | |
| CN109416241A (zh) | 摆动碾压装置的检查装置、检查用工具、检查方法、轴承单元的制造装置及轴承单元的制造方法 | |
| CN106089993A (zh) | 一种获取球锥混合双列轮毂轴承载荷分布的方法及装置 | |
| CN106294971A (zh) | 一种超大尺寸环件控稳轧制抱辊力设计方法 | |
| CN102890475A (zh) | 大型回转体零件表面轮廓误差测量及实时补偿方法 | |
| CN104008274A (zh) | 一种预测高速主轴系统切削稳定性区域的方法 | |
| CN103791870A (zh) | 平面二次包络环面蜗杆测量方法 | |
| CN111894973A (zh) | 一种铜保持架结构的双列圆锥滚子轴承 | |
| CN103954503B (zh) | 一种用于测量环件刚度的装置 | |
| CN106441887B (zh) | 一种转盘轴承滚动体与外圈滚道间接触载荷的测量方法 | |
| CN204286265U (zh) | 转盘轴承齿圈径向跳动测量装置 | |
| CN206709723U (zh) | 圆跳动检测机 | |
| CN202101655U (zh) | 大型球轴承测量仪 | |
| CN116579104A (zh) | 内外圈同时转动的中介轴承内部动态载荷获取方法及系统 | |
| CN103925898B (zh) | 一种直接测量面对面角接触轴承凸出量的方法 | |
| CN206399400U (zh) | 一种双半外圈交叉滚子轴承、四点接触球轴承轴向游隙检测装置 | |
| CN201792204U (zh) | 超短圆柱滚动回转工作台 | |
| CN205111328U (zh) | 一种风电轮毂类机械的高精度回转工装 | |
| CN103823247B (zh) | 一种重力加速度测量方法 | |
| 陈洋 et al. | Research on rotating motion of cylindrical roller based on double-disc straight groove grinding | |
| CN203773069U (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 |