CN113792396A - 一种基于直接刚度法的行星螺纹滚柱轴承载荷分布规律计算方法 - Google Patents

一种基于直接刚度法的行星螺纹滚柱轴承载荷分布规律计算方法 Download PDF

Info

Publication number
CN113792396A
CN113792396A CN202110982809.9A CN202110982809A CN113792396A CN 113792396 A CN113792396 A CN 113792396A CN 202110982809 A CN202110982809 A CN 202110982809A CN 113792396 A CN113792396 A CN 113792396A
Authority
CN
China
Prior art keywords
contact
thread
deformation
bearing
axial
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
Application number
CN202110982809.9A
Other languages
English (en)
Other versions
CN113792396B (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.)
Northeast Forestry University
Original Assignee
Northeast Forestry University
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 Northeast Forestry University filed Critical Northeast Forestry University
Priority to CN202110982809.9A priority Critical patent/CN113792396B/zh
Publication of CN113792396A publication Critical patent/CN113792396A/zh
Application granted granted Critical
Publication of CN113792396B publication Critical patent/CN113792396B/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/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design
    • 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
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

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)
  • Support Of The Bearing (AREA)

Abstract

本发明提出一种基于直接刚度法的行星螺纹滚柱轴承载荷分布规律计算方法，属于行星滚柱丝杠及行星螺纹滚柱轴承接触载荷计算技术领域；根据行星螺纹滚柱轴承的结构特征，对其结构进行离散化并将离散化结构简化为等效弹簧；计算螺纹牙刚度、轴段刚度及赫兹接触刚度，建立行星螺纹滚柱轴承的一维等效弹簧网络模型；利用直接刚度法，计算外部载荷作用下，弹簧网络内各个弹簧的内力与变形，得到对应的螺纹牙轴向接触载荷，从而计算接触变形与接触应力；本发明在考虑所有结构变形的基础上，通过直接刚度法精确计算了行星螺纹滚柱轴承的接触载荷分布，无需采用迭代计算，计算速度快，实用性与通用性好。

Description

一种基于直接刚度法的行星螺纹滚柱轴承载荷分布规律计算 方法
技术领域
本发明属于行星滚柱丝杠及行星螺纹滚柱轴承领域，具体地，涉及一种基于直接刚度法的行星螺纹滚柱轴承载荷分布规律计算方法。
背景技术
行星螺纹滚柱轴承是一种螺纹滚柱作为滚动体的滚动轴承，为行星滚柱丝杠的衍生机构，其主要结构包括保持架与具有环形螺纹的轴承外圈，轴承内圈和滚柱，其具有大承载和高刚度的性能优势，将其应用于电动执行机构将显著提高系统的功率体积比。由于行星螺纹滚柱轴承与行星滚柱丝杠同属于螺旋传动机构，因此其内部螺纹接触载荷分布规律是影响其力学性能指标的重要因素。行星螺纹滚柱轴承在轴向载荷作用下，其内部存在螺纹牙变形、轴段变形与接触变形，且三种变形之间存在相互作用关系。因此，为计算行星螺纹滚柱轴承各个螺纹接触点处的接触应力及应变，需要从整体的角度出发充分考虑其内部弹性变形之间的影响关系，计算其内部载荷分布。由于行星螺纹滚柱轴承中滚柱与轴承外圈、滚柱与轴承内圈的螺纹牙之间的接触，可以等效为多个角接触球轴承的组合，因此目前已有的球轴承接触载荷计算方法对于行星螺纹滚柱轴承的接触载荷计算存在一定的借鉴意义。但球轴承的接触载荷分布计算方法多数仅分析了轴向与周向上单个滚动体与轴承内圈及外圈的接触力学特性。因此，需要在获得行星螺纹滚柱轴承单个接触点接触载荷计算方法的基础上，针对行星螺纹滚柱轴承的结构和变形方式，提出其载荷分布规律计算方法，进一步分析其每个螺纹接触区域的接触力学特性，为其承载特性的研究工作提供理论支持。
发明内容
本发明提出了一种基于直接刚度法的行星螺纹滚柱轴承载荷分布规律计算方法，能够精确计算行星螺纹滚柱轴承的接触载荷分布，无需采用迭代计算，计算速度快，实用性与通用性好。
本发明是通过以下方案实现的：
一种基于直接刚度法的行星螺纹滚柱轴承载荷分布规律计算方法：
所述方法具体包括以下步骤：
步骤1：对行星螺纹滚柱轴承的结构进行离散化处理，离散为螺纹牙间的轴段、螺纹牙和螺纹牙啮合的接触区域；
步骤2：分别计算螺纹牙的赫兹接触刚度、螺纹牙自身刚度和螺纹牙的轴段刚度；
步骤3：依据行星螺纹滚柱轴承的离散化结构建立一维等效弹簧网络模型，基于各离散化结构刚度，计算各等效弹簧的等效刚度；
步骤4：利用直接刚度法，计算弹簧网络模型中各等效弹簧的内力与变形，获得对应的螺纹牙轴向接触载荷与等效变形；
步骤5：通过螺纹牙的轴向接触载荷与赫兹接触理论计算相应的赫兹接触应力与接触变形。
进一步地，在步骤1中，所述的行星螺纹滚柱轴承包括轴承外圈1、滚柱2、轴承内圈3和保持架4，所有滚柱2的两端均安装于保持架4的圆周方向的均匀分布的通孔内，滚柱2的外环槽螺纹同时与轴承内圈3的外环槽螺纹和轴承外圈1的内环槽螺纹啮合。
进一步地，所述的轴承外圈1和轴承内圈3的螺纹为标准三角形环槽螺纹，滚柱2螺纹类型的为凸弧形环槽螺纹。
进一步地，在步骤2中，
所述螺纹牙的轴段刚度为RB＝EA/p，
其中，E为材料的弹性模量，A为轴段的有效横截面积，p为螺距；
所述轴承外圈1的螺纹牙自身轴向刚度Rnt、滚柱2的螺纹牙自身轴向刚度Rrt和轴承内圈3的螺纹牙自身轴向刚度Rst计算公式如下：
Figure BDA0003229525070000021
式中，N为单位宽度上螺纹接触点在轴截面上受到的法向力，β为螺纹牙的牙侧角，εa为弯矩引起的轴向变形，εb为剪切力引起的轴向变形，εc为牙根倾斜引起的轴向变形，εd为牙根剪切变形产生的轴向变形，εe为轴承内圈径向收缩产生的轴向变形，εf为滚柱径向收缩产生的轴向变形，εg为轴承外圈径向扩张产生的轴向变形；
轴承外圈1、轴承内圈3与滚柱2三者间的螺纹牙啮合的接触区域的接触变形为
Figure BDA0003229525070000022
ρ为接触点处的主曲率，Q为接触区域的法向载荷，eε可依据ρ查表获得；
轴承外圈1和轴承内圈3与滚柱2接触侧线性化后的轴向赫兹接触刚度Rs和Rn计算公式如下：
Figure BDA0003229525070000031
Figure BDA0003229525070000032
式中，eε_s、eε_n分别为轴承内圈3和轴承外圈1的赫兹接触参数；Fave为单个螺纹牙理论平均轴向载荷值；ρs与ρn分别为滚柱与轴承内圈3、轴承外圈1接触区域的主曲率。
进一步地，在步骤3中，
所述行星螺纹滚柱轴承一维等效弹簧网络模型只考虑行星螺纹滚柱轴承的弹性变形；εSn和εNn分别为行星螺纹滚柱轴承的轴承内圈3与轴承外圈1螺纹牙轴段变形，εRn为滚柱螺纹牙轴段变形，SεTn和NεTn分别为轴承内圈3、轴承外圈1与滚柱2螺纹牙的自身变形及接触区域接触变形之和；上述变形与弹簧刚度一一对应，所对应的等效弹簧刚度R如下：
Figure BDA0003229525070000033
Figure BDA0003229525070000034
进一步地，在步骤4中，直接刚度法计算过程如下：
依据已有等效弹簧网络模型，假设节点位移的正方向，以及等效弹簧变形量的正负判定依据，获得等效弹簧的变形矩阵
Figure BDA0003229525070000035
与其节点位移矩阵
Figure BDA0003229525070000036
之间的关系，即
Figure BDA0003229525070000037
Figure BDA0003229525070000038
其中
Figure BDA0003229525070000039
为邻接矩阵，由节点位移与等效弹簧变形之间的关系式中的常系数组成；
将等效弹簧的变形量与其刚度相乘，得到等效弹簧内力：
Figure BDA00032295250700000310
Figure BDA00032295250700000311
为等效弹簧的刚度组成的对角矩阵；
Figure BDA00032295250700000312
为未知的等效弹簧内力，对于两相互啮合螺纹牙之间的等效弹簧，其值代表了螺纹牙上的轴向载荷值；
对于原本的线性等效弹簧，
Figure BDA0003229525070000041
对于已经线性化的原本非线性的等效弹簧，即赫兹接触刚度，
Figure BDA0003229525070000042
是关于QavenR/2的列向量；将式(8)代入式(9)可得
Figure BDA0003229525070000043
Figure BDA0003229525070000044
可以将等效弹簧内力
Figure BDA0003229525070000045
转化为已知的施加于等效弹簧节点上的外力
Figure BDA0003229525070000046
即轴承内圈所承受的总轴向载荷，故
Figure BDA0003229525070000047
节点的位移向量
Figure BDA0003229525070000048
可依据式(8)计算，因此行星螺纹滚柱轴承线性化等效弹簧网络中任意等效弹簧的内力及其等效变形:
Figure BDA0003229525070000049
Figure BDA00032295250700000410
其中，
Figure BDA00032295250700000411
中与等效弹簧SRTn、NRTn相对应的内力及其变形即为螺纹牙的轴向接触载荷与等效变形。
进一步地，在步骤5中，螺纹牙的赫兹接触变形可由
Figure BDA00032295250700000412
计算，赫兹接触应力可由如下公式计算：
Figure BDA00032295250700000413
式中，πeaeb为依据赫兹接触理论得到的赫兹接触参数，Pmax为最大赫兹接触应力，Q的值为Fs/cosβ。
本发明有益效果
(1)本发明根据行星螺纹滚柱轴承的结构特征，对其结构进行离散化处理；将离散化的弹性结构简化为等效弹簧，计算不同结构的轴向刚度，建立行星螺纹滚柱轴承的一维等效弹簧网络模性；依据直接刚度法，计算外部载荷作用下等效弹簧网络中各弹簧的内力与变形，得到对应的螺纹牙轴向接触载荷；依据所求得螺纹牙轴向接触载荷，利用赫兹接触理论计算相应的接触变形与接触应力；
(2)本发明可精确计算行星螺纹滚柱轴承的接触载荷分布，无需采用迭代计算，计算速度快，实用性与通用性好。
附图说明
图1为本发明中基于直接刚度法的行星螺纹滚柱轴承载荷分布规律计算方法的流程图；
图2为行星螺纹滚柱轴承的结构示意图；其中(1)为轴承外圈，(2)为滚柱，(3)为轴承内圈，(4)为保持架，(5)为弹性挡圈；
图3为螺纹牙接触界面的受力分析示意图；
图4为行星螺纹滚柱轴承的一维等效弹簧网络模型；
图5为螺纹牙接触载荷分布计算结果；其中(a)为滚柱与轴承外圈接触侧轴向接触载荷分布，(b)为滚柱与轴承内圈接触侧轴向接触载荷分布；
图6为螺纹牙赫兹接触变形分布计算结果；其中(a)为滚柱与轴承外圈接触侧的赫兹接触变形布，(b)为滚柱与轴承内圈接触侧的赫兹接触变形；
图7为螺纹牙赫兹接触应力分布计算结果；其中(a)为滚柱与轴承外圈接触侧接触应力，(b)为滚柱与轴承内圈接触侧接触应力。
具体实施方式
下面将结合本发明实施例中的附图对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例；基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
结合图1至图7，
实施例1：如图1，本实施例具体包括：
步骤1：对行星螺纹滚柱轴承的结构进行离散化处理，离散为螺纹牙间的轴段、螺纹牙和螺纹牙啮合的接触区域；
步骤2：分别计算螺纹牙的赫兹接触刚度、螺纹牙自身刚度和螺纹牙的轴段刚度；
步骤3：依据行星螺纹滚柱轴承的离散化结构建立一维等效弹簧网络模型，基于各离散化结构刚度，计算各等效弹簧的等效刚度；
步骤4：利用直接刚度法，计算弹簧网络模型中各等效弹簧的内力与变形，获得对应的螺纹牙轴向接触载荷与等效变形；
步骤5：通过螺纹牙的轴向接触载荷与赫兹接触理论计算相应的赫兹接触应力与接触变形。
其中行星螺纹滚柱轴承结构如图2所示，包括轴承外圈1、轴承内圈3、滚柱2、保持架4，多个滚柱2通过保持架4均布于轴承内圈3的周围，滚柱2的螺纹同时与轴承内圈3和轴承外圈1的螺纹啮合实现旋转支撑作用，其中轴承内圈3及轴承外圈1的螺纹是牙型角为90°的三角形螺纹，滚柱2螺纹的牙型侧面为凸弧形，轴承外圈1、滚柱2和轴承内圈3的螺纹均为导程为零的环槽螺纹；
在步骤2中，
所述螺纹牙的轴段刚度为RB＝EA/p，
其中，E为材料的弹性模量，A为轴段的有效横截面积，p为螺距；
所述轴承外圈1的螺纹牙自身轴向刚度Rnt、滚柱2的螺纹牙自身轴向刚度Rrt和轴承内圈3的螺纹牙自身轴向刚度Rst计算公式如下：
Figure BDA0003229525070000061
式中，N为单位宽度上螺纹接触点在轴截面上受到的法向力，β为螺纹牙的牙侧角，εa为弯矩引起的轴向变形，εb为剪切力引起的轴向变形，εc为牙根倾斜引起的轴向变形，εd为牙根剪切变形产生的轴向变形，εe为轴承内圈径向收缩产生的轴向变形，εf为滚柱径向收缩产生的轴向变形，εg为轴承外圈径向扩张产生的轴向变形；
轴承外圈1、轴承内圈3与滚柱2三者间的螺纹牙啮合的接触区域的接触变形为
Figure BDA0003229525070000062
ρ为接触点处的主曲率，Q为接触区域的法向载荷，eε可依据ρ查表获得；
轴承外圈1和轴承内圈3与滚柱2接触侧线性化后的轴向赫兹接触刚度Rs和Rn计算公式如下：
Figure BDA0003229525070000063
Figure BDA0003229525070000064
式中，eε_s、eε_n分别为轴承内圈3和轴承外圈1的赫兹接触参数；Fave为单个螺纹牙理论平均轴向载荷值；ρs与ρn分别为滚柱与轴承内圈3、轴承外圈1接触区域的主曲率。
在步骤3中，行星螺纹滚柱轴承一维等效弹簧网络模型如附图4所示
所述行星螺纹滚柱轴承一维等效弹簧网络模型只考虑行星螺纹滚柱轴承的弹性变形；εSn和εNn分别为行星螺纹滚柱轴承的轴承内圈3与轴承外圈1螺纹牙轴段变形，εRn为滚柱螺纹牙轴段变形，SεTn和NεTn分别为轴承内圈3、轴承外圈1与滚柱2螺纹牙的自身变形及接触区域接触变形之和；上述变形与弹簧刚度一一对应，所对应的等效弹簧刚度R如下：
Figure BDA0003229525070000071
Figure BDA0003229525070000072
在步骤4中，直接刚度法计算过程如下：
首先依据已有等效弹簧网络模型，假设节点位移的正方向，以及等效弹簧变形量的正负判定依据，获得等效弹簧的变形矩阵
Figure BDA0003229525070000073
与其节点位移矩阵
Figure BDA0003229525070000074
之间的关系，即
Figure BDA0003229525070000075
Figure BDA0003229525070000076
其中
Figure BDA0003229525070000077
为邻接矩阵，由节点位移与等效弹簧变形之间的关系式中的常系数组成；
将等效弹簧的变形量与其刚度相乘，得到等效弹簧内力：
Figure BDA0003229525070000078
Figure BDA0003229525070000079
为等效弹簧的刚度组成的对角矩阵；
Figure BDA00032295250700000710
为未知的等效弹簧内力，对于两相互啮合螺纹牙之间的等效弹簧，其值代表了螺纹牙上的轴向载荷值；
对于原本的线性等效弹簧，
Figure BDA00032295250700000711
对于已经线性化的原本非线性的等效弹簧，即赫兹接触刚度，
Figure BDA00032295250700000712
是关于QavenR/2的列向量；将式(8)代入式(9)可得
Figure BDA00032295250700000713
Figure BDA00032295250700000714
可以将等效弹簧内力
Figure BDA00032295250700000715
转化为已知的施加于等效弹簧节点上的外力
Figure BDA00032295250700000716
即轴承内圈所承受的总轴向载荷，故
Figure BDA0003229525070000081
节点的位移向量
Figure BDA0003229525070000082
可依据式(8)计算，因此行星螺纹滚柱轴承线性化等效弹簧网络中任意等效弹簧的内力及其等效变形:
Figure BDA0003229525070000083
Figure BDA0003229525070000084
其中，
Figure BDA0003229525070000085
中与等效弹簧SRTn、NRTn相对应的内力及其变形即为螺纹牙的轴向接触载荷与等效变形。
在步骤5中，螺纹牙的赫兹接触变形可由
Figure BDA0003229525070000086
计算，赫兹接触应力可由如下公式计算：
Figure BDA0003229525070000087
式中，πeaeb为依据赫兹接触理论得到的赫兹接触参数，Pmax为最大赫兹接触应力，Q的值为Fs/cosβ。
通过上述计算，在已知行星螺纹滚柱轴承的结构尺寸参数、安装方式以及所受外部轴向载荷的条件下，即可获得行星螺纹滚柱轴承内部螺纹接触区域的轴向接触载荷、法向接触变形以及接触应力。
行星螺纹滚柱轴承尺寸参数如表1所示，选取所加载的外部轴向载荷为50kN，则行星螺纹滚柱轴承内部不同接触区域螺纹牙的接触载荷分布如图5所示；行星螺纹滚柱轴承内部不同接触区域螺纹牙的法向赫兹接触变形分布如图6所示；行星螺纹滚柱轴承内部不同接触区域螺纹牙的赫兹接触应力分布如图7所示。
Figure BDA0003229525070000088
Figure BDA0003229525070000091
表1有限元仿真所用模型尺寸参数
由附图5可知，行星螺纹滚柱轴承不同接触区域的螺纹接触载荷分布不均程度不同。滚柱与轴承外圈接触的最大螺纹法向接触载荷小于滚柱与轴承内圈接触侧，且该侧的接触载荷分布总体上较为均匀，轴承内部两接触侧的螺纹接触载荷自受力端起开始下降，且滚柱与轴承内圈接触侧的接触载荷呈现几何下降的趋势，当螺纹牙序号超过30时接触载荷的变化趋于平缓。由附图6与附图7可知，行星螺纹滚柱轴承的螺纹赫兹接触变形与赫兹接触应力的分布情况与其螺纹法向接触载荷分布情况相近。
以上对本发明所提出的一种基于直接刚度法的行星螺纹滚柱轴承载荷分布规律计算方法，进行了详细介绍，对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。

Claims (7)

1.一种基于直接刚度法的行星螺纹滚柱轴承载荷分布规律计算方法，其特征在于：
所述方法具体包括以下步骤：
步骤1：对行星螺纹滚柱轴承的结构进行离散化处理，离散为螺纹牙间的轴段、螺纹牙本体和螺纹牙间啮合的接触区域；
步骤2：分别计算螺纹牙的赫兹接触刚度、螺纹牙自身刚度和螺纹牙的轴段刚度；
步骤3：依据行星螺纹滚柱轴承的离散化结构建立一维等效弹簧网络模型，基于各离散化结构刚度，计算各等效弹簧的等效刚度；
步骤4：利用直接刚度法，计算弹簧网络模型中各等效弹簧的内力与变形，获得对应的螺纹牙轴向接触载荷与等效变形；
步骤5：通过螺纹牙的轴向接触载荷与赫兹接触理论计算相应的赫兹接触应力与接触变形。
2.根据权利要求1所述方法，其特征在于：在步骤1中，所述的行星螺纹滚柱轴承包括轴承外圈(1)、滚柱(2)、轴承内圈(3)和保持架(4)，所有滚柱(2)的两端均安装于保持架(4)的圆周方向的均匀分布的通孔内，滚柱(2)的外环槽螺纹同时与轴承内圈(3)的外环槽螺纹和轴承外圈(1)的内环槽螺纹啮合。
3.根据权利要求2所述方法，其特征在于：所述的轴承外圈(1)和轴承内圈(3)的螺纹为标准三角形环槽螺纹，滚柱(2)螺纹类型的为凸弧形环槽螺纹。
4.根据权利要求2所述方法，其特征在于：在步骤2中，
所述螺纹牙的轴段刚度为RB＝EA/p，
其中，E为材料的弹性模量，A为轴段的有效横截面积，p为螺距；
所述轴承外圈(1)的螺纹牙自身轴向刚度Rnt、滚柱(2)的螺纹牙自身轴向刚度Rrt和轴承内圈(3)的螺纹牙自身轴向刚度Rst计算公式如下：
Figure FDA0003229525060000011
式中，N为单位宽度上螺纹接触点在轴截面上受到的法向力，β为螺纹牙的牙侧角，εa为弯矩引起的轴向变形，εb为剪切力引起的轴向变形，εc为牙根倾斜引起的轴向变形，εd为牙根剪切变形产生的轴向变形，εe为轴承内圈径向收缩产生的轴向变形，εf为滚柱径向收缩产生的轴向变形，εg为轴承外圈径向扩张产生的轴向变形；
轴承外圈(1)、轴承内圈(3)与滚柱(2)三者间的螺纹牙啮合的接触区域的接触变形为
Figure FDA0003229525060000021
ρ为接触点处的主曲率，Q为接触区域的法向载荷，eε可依据ρ查表获得；
轴承外圈(1)和轴承内圈(3)与滚柱(2)接触侧线性化后的轴向赫兹接触刚度Rs和Rn计算公式如下：
Figure FDA0003229525060000022
Figure FDA0003229525060000023
式中，eε_s、eε_n分别为轴承内圈(3)和轴承外圈(1)的赫兹接触参数；Fave为单个螺纹牙理论平均轴向载荷值；ρs与ρn分别为滚柱与轴承内圈(3)、轴承外圈(1)接触区域的主曲率。
5.根据权利要求4所述方法，其特征在于：在步骤3中，
所述行星螺纹滚柱轴承一维等效弹簧网络模型只考虑行星螺纹滚柱轴承的弹性变形；εSn和εNn分别为行星螺纹滚柱轴承的轴承内圈(3)与轴承外圈(1)螺纹牙轴段变形，εRn为滚柱螺纹牙轴段变形，
Figure FDA0003229525060000024
和
Figure FDA0003229525060000025
分别为轴承内圈(3)、轴承外圈(1)与滚柱(2)螺纹牙的自身变形及接触区域接触变形之和；上述变形与弹簧刚度一一对应，所对应的等效弹簧刚度R如下：
Figure FDA0003229525060000026
Figure FDA0003229525060000027
。
6.根据权利要求5所述方法，其特征在于：
在步骤4中，直接刚度法计算过程如下：
依据已有等效弹簧网络模型，假设节点位移的正方向，以及等效弹簧变形量的正负判定依据，获得等效弹簧的变形矩阵
Figure FDA0003229525060000028
与其节点位移矩阵
Figure FDA0003229525060000029
之间的关系，即
Figure FDA00032295250600000210
Figure FDA0003229525060000031
其中
Figure FDA0003229525060000032
为邻接矩阵，由节点位移与等效弹簧变形之间的关系式中的常系数组成；
将等效弹簧的变形量与其刚度相乘，得到等效弹簧内力：
Figure FDA0003229525060000033
Figure FDA0003229525060000034
为等效弹簧的刚度组成的对角矩阵；
Figure FDA0003229525060000035
为未知的等效弹簧内力，对于两相互啮合螺纹牙之间的等效弹簧，其值代表了螺纹牙上的轴向载荷值；
对于原本的线性等效弹簧，
Figure FDA0003229525060000036
对于已经线性化的原本非线性的等效弹簧，即赫兹接触刚度，
Figure FDA0003229525060000037
是关于QavenR/2的列向量；将式(8)代入式(9)可得
Figure FDA0003229525060000038
Figure FDA0003229525060000039
可以将等效弹簧内力
Figure FDA00032295250600000310
转化为已知的施加于等效弹簧节点上的外力
Figure FDA00032295250600000311
即轴承内圈所承受的总轴向载荷，故
Figure FDA00032295250600000312
节点的位移向量
Figure FDA00032295250600000313
可依据式(8)计算，因此行星螺纹滚柱轴承线性化等效弹簧网络中任意等效弹簧的内力及其等效变形:
Figure FDA00032295250600000314
Figure FDA00032295250600000315
其中，
Figure FDA00032295250600000316
中与等效弹簧SRTn、NRTn相对应的内力及其变形即为螺纹牙的轴向接触载荷与等效变形。
7.根据权利要求6所述方法，其特征在于：
在步骤5中，螺纹牙的赫兹接触变形可由
Figure FDA00032295250600000317
计算，赫兹接触应力可由如下公式计算：
Figure FDA00032295250600000318
式中，πeaeb为依据赫兹接触理论得到的赫兹接触参数，Pmax为最大赫兹接触应力，Q的值为Fs/cosβ。
CN202110982809.9A 2021-08-25 2021-08-25 一种基于直接刚度法的行星螺纹滚柱轴承载荷分布规律计算方法 Active CN113792396B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110982809.9A CN113792396B (zh) 2021-08-25 2021-08-25 一种基于直接刚度法的行星螺纹滚柱轴承载荷分布规律计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110982809.9A CN113792396B (zh) 2021-08-25 2021-08-25 一种基于直接刚度法的行星螺纹滚柱轴承载荷分布规律计算方法

Publications (2)

Publication Number Publication Date
CN113792396A true CN113792396A (zh) 2021-12-14
CN113792396B CN113792396B (zh) 2023-08-22

Family

ID=79182295

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110982809.9A Active CN113792396B (zh) 2021-08-25 2021-08-25 一种基于直接刚度法的行星螺纹滚柱轴承载荷分布规律计算方法

Country Status (1)

Country Link
CN (1) CN113792396B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20180365357A1 (en) * 2015-11-18 2018-12-20 Romax Technology Limited Driveline Modeller
CN111199013A (zh) * 2020-02-21 2020-05-26 内蒙古工业大学 一种针对行星滚柱丝杠副接触滚道的主曲率计算方法
CN111881529A (zh) * 2020-07-27 2020-11-03 大连交通大学 机器人用rv减速器主轴承受力分析计算方法
CN112464481A (zh) * 2020-12-02 2021-03-09 华北水利水电大学 机器人用摆线针轮减速机的动态传动精度数值计算方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20180365357A1 (en) * 2015-11-18 2018-12-20 Romax Technology Limited Driveline Modeller
CN111199013A (zh) * 2020-02-21 2020-05-26 内蒙古工业大学 一种针对行星滚柱丝杠副接触滚道的主曲率计算方法
CN111881529A (zh) * 2020-07-27 2020-11-03 大连交通大学 机器人用rv减速器主轴承受力分析计算方法
CN112464481A (zh) * 2020-12-02 2021-03-09 华北水利水电大学 机器人用摆线针轮减速机的动态传动精度数值计算方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
姚琴;吴立言;刘更;马尚君;佟瑞庭;: "行星滚柱丝杠副弹塑性接触分析", 机械传动, no. 11, pages 106 - 116 *

Also Published As

Publication number Publication date
CN113792396B (zh) 2023-08-22

Similar Documents

Publication Publication Date Title
Liu et al. An analytical calculation method of the load distribution and stiffness of an angular contact ball bearing
CN108984933B (zh) 弹流润滑条件下计算滚动轴承载荷和压力的边界元法
Zhen et al. Analysis of stress and fatigue life of ball screw with considering the dimension errors of balls
CN106649982B (zh) 风力发电机大锥角圆锥主轴承摩擦力矩计算方法
CN107247856B (zh) 一种单滚柱包络环面蜗杆副时变啮合刚度解析方法
CN110516345B (zh) 行星滚柱丝杠误差源耦合建模方法
Shi et al. Calculation of time-varying backlash for an involute spur gear pair
CN109753723B (zh) 一种向心滚动轴承疲劳寿命计算方法
Yongqi et al. Analysis of Stress and Strain of the Rolling Bearing by FEA method
Park et al. Load sharing and distributed on the gear flank of wind turbine planetary gearbox
Mitrović et al. Effects of operation temperature on thermal expansion and main parameters of radial ball bearings
CN107563082B (zh) 基于圆柱滚子轴承接触变形量和载荷分布的轴承参数优化方法
CN113792396A (zh) 一种基于直接刚度法的行星螺纹滚柱轴承载荷分布规律计算方法
CN113392544B (zh) 一种基于变形协调理论的行星螺纹滚柱轴承接触载荷计算方法
CN104675947B (zh) 一种行星滚柱丝杠及其载荷均布方法
Dindar et al. Experimental determination and analytical model of friction torque of a double row roller slewing bearing
Jingyang et al. The load distribution of the main shaft bearing considering combined load and misalignment in a floating direct-drive wind turbine
Göncz et al. Load capacity of a three-row roller slewing bearing raceway
Blinov et al. On unevenness of load distribution between mating thread ridges of parts of roller screw gears
Zhilnikov et al. A method of calculating the friction moment in cageless bearings
Bourdon et al. Static rolling bearing models in a CAD environment for the study of complex mechanisms: Part II—complete assembly model
CN114722660A (zh) 行星滚柱丝杠刚柔耦合动态接触特性分析方法
Shan et al. A new method for determining the preload in a wire race ball bearing
CN112949123B (zh) 一种考虑涂层影响的点接触力载与变形关系的建立方法
CN108009342B (zh) 一种w2r型轴连轴承及其载荷分布分析方法

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