CN113792396B - 一种基于直接刚度法的行星螺纹滚柱轴承载荷分布规律计算方法 - Google Patents
一种基于直接刚度法的行星螺纹滚柱轴承载荷分布规律计算方法 Download PDFInfo
- Publication number
- CN113792396B CN113792396B CN202110982809.9A CN202110982809A CN113792396B CN 113792396 B CN113792396 B CN 113792396B CN 202110982809 A CN202110982809 A CN 202110982809A CN 113792396 B CN113792396 B CN 113792396B
- Authority
- CN
- China
- Prior art keywords
- thread
- deformation
- contact
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 46
- 238000004364 calculation method Methods 0.000 claims abstract description 25
- 238000006073 displacement reaction Methods 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 10
- 239000011159 matrix material Substances 0.000 claims description 9
- 230000005489 elastic deformation Effects 0.000 claims description 4
- 238000005452 bending Methods 0.000 claims description 3
- 230000008602 contraction Effects 0.000 claims description 3
- 239000000463 material Substances 0.000 claims description 3
- 238000006467 substitution reaction Methods 0.000 claims 1
- 230000009471 action Effects 0.000 abstract description 3
- 230000007246 mechanism Effects 0.000 description 3
- 238000005096 rolling process Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
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)
- Support Of The Bearing (AREA)
- Rolling Contact Bearings (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计算公式如下:
式中,N为单位宽度上螺纹接触点在轴截面上受到的法向力,β为螺纹牙的牙侧角,εa为弯矩引起的轴向变形,εb为剪切力引起的轴向变形,εc为牙根倾斜引起的轴向变形,εd为牙根剪切变形产生的轴向变形,εe为轴承内圈径向收缩产生的轴向变形,εf为滚柱径向收缩产生的轴向变形,εg为轴承外圈径向扩张产生的轴向变形;
轴承外圈1、轴承内圈3与滚柱2三者间的螺纹牙啮合的接触区域的接触变形为ρ为接触点处的主曲率,Q为接触区域的法向载荷,eε可依据ρ查表获得;
轴承外圈1和轴承内圈3与滚柱2接触侧线性化后的轴向赫兹接触刚度Rs和Rn计算公式如下:
式中,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如下:
进一步地,在步骤4中,直接刚度法计算过程如下:
依据已有等效弹簧网络模型,假设节点位移的正方向,以及等效弹簧变形量的正负判定依据,获得等效弹簧的变形矩阵与其节点位移矩阵/>之间的关系,即/>
其中为邻接矩阵,由节点位移与等效弹簧变形之间的关系式中的常系数组成;
将等效弹簧的变形量与其刚度相乘,得到等效弹簧内力:
为等效弹簧的刚度组成的对角矩阵;/>为未知的等效弹簧内力,对于两相互啮合螺纹牙之间的等效弹簧,其值代表了螺纹牙上的轴向载荷值;
对于原本的线性等效弹簧,对于已经线性化的原本非线性的等效弹簧,即赫兹接触刚度,/>是关于QavenR/2的列向量;将式(8)代入式(9)可得
可以将等效弹簧内力/>转化为已知的施加于等效弹簧节点上的外力/>即轴承内圈所承受的总轴向载荷,故
节点的位移向量可依据式(8)计算,因此行星螺纹滚柱轴承线性化等效弹簧网络中任意等效弹簧的内力及其等效变形:
其中,中与等效弹簧SRTn、NRTn相对应的内力及其变形即为螺纹牙的轴向接触载荷与等效变形。
进一步地,在步骤5中,螺纹牙的赫兹接触变形可由计算,赫兹接触应力可由如下公式计算:
式中,π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计算公式如下:
式中,N为单位宽度上螺纹接触点在轴截面上受到的法向力,β为螺纹牙的牙侧角,εa为弯矩引起的轴向变形,εb为剪切力引起的轴向变形,εc为牙根倾斜引起的轴向变形,εd为牙根剪切变形产生的轴向变形,εe为轴承内圈径向收缩产生的轴向变形,εf为滚柱径向收缩产生的轴向变形,εg为轴承外圈径向扩张产生的轴向变形;
轴承外圈1、轴承内圈3与滚柱2三者间的螺纹牙啮合的接触区域的接触变形为ρ为接触点处的主曲率,Q为接触区域的法向载荷,eε可依据ρ查表获得;
轴承外圈1和轴承内圈3与滚柱2接触侧线性化后的轴向赫兹接触刚度Rs和Rn计算公式如下:
式中,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如下:
在步骤4中,直接刚度法计算过程如下:
首先依据已有等效弹簧网络模型,假设节点位移的正方向,以及等效弹簧变形量的正负判定依据,获得等效弹簧的变形矩阵与其节点位移矩阵/>之间的关系,即/>
其中为邻接矩阵,由节点位移与等效弹簧变形之间的关系式中的常系数组成;
将等效弹簧的变形量与其刚度相乘,得到等效弹簧内力:
为等效弹簧的刚度组成的对角矩阵;/>为未知的等效弹簧内力,对于两相互啮合螺纹牙之间的等效弹簧,其值代表了螺纹牙上的轴向载荷值;
对于原本的线性等效弹簧,对于已经线性化的原本非线性的等效弹簧,即赫兹接触刚度,/>是关于QavenR/2的列向量;将式(8)代入式(9)可得
可以将等效弹簧内力/>转化为已知的施加于等效弹簧节点上的外力/>即轴承内圈所承受的总轴向载荷,故
节点的位移向量可依据式(8)计算,因此行星螺纹滚柱轴承线性化等效弹簧网络中任意等效弹簧的内力及其等效变形:
其中,中与等效弹簧SRTn、NRTn相对应的内力及其变形即为螺纹牙的轴向接触载荷与等效变形。
在步骤5中,螺纹牙的赫兹接触变形可由计算,赫兹接触应力可由如下公式计算:
式中,πeaeb为依据赫兹接触理论得到的赫兹接触参数,Pmax为最大赫兹接触应力,Q的值为Fs/cosβ。
通过上述计算,在已知行星螺纹滚柱轴承的结构尺寸参数、安装方式以及所受外部轴向载荷的条件下,即可获得行星螺纹滚柱轴承内部螺纹接触区域的轴向接触载荷、法向接触变形以及接触应力。
行星螺纹滚柱轴承尺寸参数如表1所示,选取所加载的外部轴向载荷为50kN,则行星螺纹滚柱轴承内部不同接触区域螺纹牙的接触载荷分布如图5所示;行星螺纹滚柱轴承内部不同接触区域螺纹牙的法向赫兹接触变形分布如图6所示;行星螺纹滚柱轴承内部不同接触区域螺纹牙的赫兹接触应力分布如图7所示。
表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计算公式如下:
式中,N为单位宽度上螺纹接触点在轴截面上受到的法向力,β为螺纹牙的牙侧角,εa为弯矩引起的轴向变形,εb为剪切力引起的轴向变形,εc为牙根倾斜引起的轴向变形,εd为牙根剪切变形产生的轴向变形,εe为轴承内圈径向收缩产生的轴向变形,εf为滚柱径向收缩产生的轴向变形,εg为轴承外圈径向扩张产生的轴向变形;
轴承外圈(1)、轴承内圈(3)与滚柱(2)三者间的螺纹牙啮合的接触区域的接触变形为ρ为接触点处的主曲率,Q为接触区域的法向载荷,eε可依据ρ查表获得;
轴承外圈(1)和轴承内圈(3)与滚柱(2)接触侧线性化后的轴向赫兹接触刚度Rs和Rn计算公式如下:
式中,eε_s、eε_n分别为轴承内圈(3)和轴承外圈(1)的赫兹接触参数;Fave为单个螺纹牙理论平均轴向载荷值;ρs与ρn分别为滚柱与轴承内圈(3)、轴承外圈(1)接触区域的主曲率。
5.根据权利要求4所述方法,其特征在于:在步骤3中,
所述行星螺纹滚柱轴承一维等效弹簧网络模型只考虑行星螺纹滚柱轴承的弹性变形;εSn和εNn分别为行星螺纹滚柱轴承的轴承内圈(3)与轴承外圈(1)螺纹牙轴段变形,εRn为滚柱螺纹牙轴段变形,SεTn和NεTn分别为轴承内圈(3)、轴承外圈(1)与滚柱(2)螺纹牙的自身变形及接触区域接触变形之和;上述变形与弹簧刚度一一对应,所对应的等效弹簧的等效刚度R如下:
6.根据权利要求5所述方法,其特征在于:
在步骤4中,直接刚度法计算过程如下:
依据已有等效弹簧网络模型,假设节点位移的正方向,以及等效弹簧变形量的正负判定依据,获得等效弹簧的变形矩阵与其节点位移矩阵/>之间的关系,即/>
其中为邻接矩阵,由节点位移与等效弹簧变形之间的关系式中的常系数组成;
将等效弹簧的变形量与其刚度相乘,得到等效弹簧内力:
为等效弹簧的等效刚度组成的对角矩阵;/>为未知的等效弹簧内力,对于两相互啮合螺纹牙之间的等效弹簧,其值代表了螺纹牙上的轴向载荷值;
对于原本的线性等效弹簧,
对于已经线性化的原本非线性的等效弹簧,即赫兹接触刚度;将式代入式(6)可得:
可以将等效弹簧内力/>转化为已知的施加于等效弹簧节点上的外力/>即轴承内圈所承受的总轴向载荷,故
节点的位移向量可依据式(8)计算,因此行星螺纹滚柱轴承线性化等效弹簧网络中任意等效弹簧的内力及其等效变形:
其中,中与等效弹簧SRTn、NRTn相对应的内力及其变形即为螺纹牙的轴向接触载荷与等效变形。
7.根据权利要求6所述方法,其特征在于:
在步骤5中,螺纹牙的赫兹接触变形可由计算,赫兹接触应力可由如下公式计算:
式中,πeaeb为依据赫兹接触理论得到的赫兹接触参数,Pmax为最大赫兹接触应力,Q的值为Fs/cosβ。
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 CN113792396A (zh) | 2021-12-14 |
CN113792396B true 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 (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 | 华北水利水电大学 | 机器人用摆线针轮减速机的动态传动精度数值计算方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB201520361D0 (en) * | 2015-11-18 | 2015-12-30 | Romax Technology Ltd | Driveline modeller |
-
2021
- 2021-08-25 CN CN202110982809.9A patent/CN113792396B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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)
Title |
---|
行星滚柱丝杠副弹塑性接触分析;姚琴;吴立言;刘更;马尚君;佟瑞庭;;机械传动(第11期);第106-116页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113792396A (zh) | 2021-12-14 |
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 | |
Zhen et al. | Analysis of stress and fatigue life of ball screw with considering the dimension errors of balls | |
CN106649982B (zh) | 风力发电机大锥角圆锥主轴承摩擦力矩计算方法 | |
Hwang et al. | Development of automatic variable preload device for spindle bearing by using centrifugal force | |
CN107729597A (zh) | 一种主轴轴承滚道校核工具 | |
CN109753723B (zh) | 一种向心滚动轴承疲劳寿命计算方法 | |
Keller et al. | Gearbox reliability collaborative investigation of gearbox motion and high-speed-shaft loads | |
Mitrović et al. | Effects of operation temperature on thermal expansion and main parameters of radial ball bearings | |
Śpiewak | Methodology for calculating the complete static carrying capacity of twin slewing bearing | |
CN113792396B (zh) | 一种基于直接刚度法的行星螺纹滚柱轴承载荷分布规律计算方法 | |
CN107563082B (zh) | 基于圆柱滚子轴承接触变形量和载荷分布的轴承参数优化方法 | |
Liu et al. | Load distribution calculation of a four-point-contact slewing bearing and its experimental verification | |
Park et al. | Characteristic analysis of wind turbine gearbox considering non-torque loading | |
Song et al. | Computational investigation of off-sized bearing rollers on dynamics for hypoid gear-shaft-bearing coupled system | |
CN104675947B (zh) | 一种行星滚柱丝杠及其载荷均布方法 | |
Zheng et al. | A novel planetary thread roller bearing: Design and analysis of load characteristic | |
Gao et al. | Optimization and experimental research on a new-type short cylindrical cup-shaped harmonic reducer | |
CN107577881B (zh) | 基于当量弹性模量和载荷分布的圆柱滚子轴承参数优化方法 | |
Yang et al. | Dynamics numerical simulation of planetary gear system for wind turbine gearbox | |
Blinov et al. | On unevenness of load distribution between mating thread ridges of parts of roller screw gears | |
Yang et al. | The Dynamic analysis of Planetary roller screw in Tree-climbing robot | |
Bourdon et al. | Static rolling bearing models in a CAD environment for the study of complex mechanisms: Part II—complete assembly model | |
Yu et al. | Load analysis and deformation research of the flexible bearing based on a three-force ring superposition method | |
Qin et al. | Design and structural responses of a 38-meter sectional wind turbine blade under extreme static loads | |
Zhang et al. | Load distribution and axial static rigidity for planetary roller screw mechanism |
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 |