CN106547980B - 一种滚动轴承套圈固有频率计算方法 - Google Patents
一种滚动轴承套圈固有频率计算方法 Download PDFInfo
- Publication number
- CN106547980B CN106547980B CN201610969618.8A CN201610969618A CN106547980B CN 106547980 B CN106547980 B CN 106547980B CN 201610969618 A CN201610969618 A CN 201610969618A CN 106547980 B CN106547980 B CN 106547980B
- Authority
- CN
- China
- Prior art keywords
- rolling bearing
- thin slice
- value
- lasso
- calculation method
- 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
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
Landscapes
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Support Of The Bearing (AREA)
- Rolling Contact Bearings (AREA)
Abstract
本发明公开了一种滚动轴承套圈固有频率计算方法,包括基本参数确定及频率精度设置、薄片圆环单元划分、基础数据计算、薄片圆环单元固有频率计算、薄片圆环单元刚度计算、轴承套圈总刚度计算、计算轴承套圈的第n阶固有频率、计算结果判断和计算结果确认的步骤。其中,薄片圆环单元划分数量为100~10000;轴承套圈总刚度为各个薄片圆环单元刚度的累加求和,所得频率为套圈的第n阶径向弯曲频率。本发明的有益效果是,采用切片和积分的方法,对轴承套圈进行离散计算后。再进行积分处理,不仅简化了滚动轴承套圈固有频率的计算方法,而且还提高了计算效率和精度,有效确保了滚动轴承套圈固有频率计算的准确性、高效性和高精度等要求。
Description
技术领域
本发明涉及一种固有频率计算方法,特别是涉及一种滚动轴承套圈固有频率计算方法。
背景技术
滚动轴承作为一种重要的关键基础部件,其动态性能直接影响整机的工作性能。随着机械装备向低振动、低噪声、高精度、高速、重载等方向发展,对圆柱滚子轴承的动态性能、振动噪声水平和可靠性的要求也日益提高。而滚动轴承套圈的固有频率与轴承的振动特性直接相关。因此,在设计阶段对滚动轴承套圈的固有频率进行准确计算对轴承振动特性的预测具有重要作用。
虽然目前针对滚动轴承套圈的固有频率的计算,出现了经验公式法、解析法和有限元计算方法。但这些计算方法都存在各种不足和问题,经验公式法具有较大误差,解析法过程十分复杂,而有限元方法严重依赖网格精度且计算效率低。为此,需要开发一种计算效率和计算精度均高的滚动轴承套圈固有频率计算方法。
发明内容
本发明的目的就是针对滚动轴承套圈固有频率传统计算方法的不足,提出一种滚动轴承套圈固有频率计算新方法,其采用切片和积分的方法,对轴承套圈进行离散计算后,再进行积分处理,不仅简化了滚动轴承套圈固有频率的计算方法,而且还提高了计算效率和精度,有效确保了滚动轴承套圈固有频率计算的准确性、高效性和高精度等要求。
为实现上述目的,本发明采用如下技术方案。
一种滚动轴承套圈固有频率计算方法,包括以下步骤:
第一步,基本参数确定及频率精度设置:根据待计算的滚动轴承套圈的规格和材质的物理特性,获取滚动轴承套圈的几何尺寸参数和物理性能参数,并设置固有频率精度;其中,几何尺寸参数包括滚动轴承套圈的外径D、内径D1、套圈宽度B、沟道形状和位置,如沟道曲率半径r、沟底直径d1等;物理性能参数包括滚动轴承套圈包括弹性模量E和密度ρ;
第二步,薄片圆环单元划分:将滚动轴承套圈按宽度B方向划分成N个等厚度的薄片圆环单元,其中,N为大于或等于100的自然数;
第三步,基础数据计算:将单个薄片圆环单元的纵向截面视为矩形,按矩形截面的计算方法,计算每个薄片圆环单元的截面惯性矩Ii、质量Mi和横截面面积Ai,并计算中线半径Ri,该中线半径Ri等于薄片圆环单元内外半径之和的一半;
第四步,薄片圆环单元固有频率计算:根据物理性能参数和基础数据计算结果,按圆环的固有频率计算方法,计算第i个薄片圆环单元的n阶径向弯曲频率ωni;
第五步,薄片圆环单元刚度计算:根据第n阶径向弯曲频率ωn与刚度Kn的关系计算第i个薄片圆环单元的刚度Kni;
第六步,轴承套圈总刚度计算:对N个薄片圆环单元的刚度累加求和计算轴承套圈总刚度Knz;
第七步,计算轴承套圈的第n阶固有频率:根据第n阶径向弯曲频率ωn与刚度Kn的关系,以轴承套圈的总刚度Knz计算轴承套圈的第n阶固有频率ωn:
第八步,计算结果判断:将本次与紧前一次的ωn值计算结果比较,当二者差值的绝对值大于固有频率精度的设定值时,则返回第二步,改变所述N的取值,重复进行第三步~第八步;当二者差值的绝对值小于或等于频率精度的设定值时,则进入第九步;
第九步,计算结果确认:结束计算,并将本次计算结果的ωn值确定为该滚动轴承套圈的固有频率。
采用前述技术方案的本发明,采用切片和积分的方法,对轴承套圈进行离散计算后,再进行积分处理。相对于传统计算方法,不仅简化了滚动轴承套圈固有频率的计算方法,而且还提高了计算效率和精度,有效确保了滚动轴承套圈固有频率计算的准确性、高效性和高精度等要求。
在本方案中的中线半径Ri计算时,薄片圆环单元内外半径的值由薄片圆环单元所处的具体位置确定;在轴承外圈频率计算中,外圈外径为D不变,内径为D1或d1综合r;在轴承内圈频率计算中,内圈内径为D1不变,外径为D或d1综合r。
优选的,所述第i个薄片圆环单元的n阶径向弯曲频率ωni按下式计算:
以此获得每个薄片圆环单元的固有频率值。
优选的,所述第i个薄片圆环单元的刚度Kni按下式计算:
以此获得每个薄片圆环单元的固有频率值。
优选的,所述滚动轴承套圈总刚度Knz按下式计算:
以此获得滚动轴承套圈总刚度值。
优选的,所述轴承套圈的第n阶固有频率按下式计算:
其中,
以此获得滚动轴承套圈的径向弯曲频率。
进一步优选的,在改变所述N的取值时,新一次取值按初始值的约2~4倍递增或按初始值的1/4~1/2倍递减;或者按紧前一次值的约2~5倍或按紧前一次值的约1/5~1/2取值,且N的最大取值不大于10000。以通过较少的计算次数获得准确的计算结果。
本发明的有益效果是,通过采用切片和积分的方法,将滚动轴承套圈划分为有限数目的薄片圆环单元,每片圆环单元的切面视为矩形截面,采用矩形截面圆环的固有频率计算方法,获得单片圆环单元的频率、刚度和质量,将各个圆环单元的刚度考虑为并联方式,求解整个滚动轴承套圈的刚度和质量,再利用圆环结构固有频率计算方法获取套圈的固有频率。该计算方法克服了传统计算方法计算复杂截面惯性矩的难题,简化了滚动轴承固有频率的计算方法,且保证了计算结果的准确性。具有简单、快速、精度高等优点。
附图说明
图1是本发明计算方法的计算流程图。
图2是应用本发明计算方法的计算对象滚动轴承套圈中外圈的典型结构示意图。
图3是本发明计算对象滚动轴承套圈划分为N个薄片圆环单元的结构示意图。
图4是本发明中计算对象为滚动轴承外圈的薄片圆环单元的结构示意图。
图5是应用本发明计算方法的计算对象滚动轴承套圈中内圈的典型结构示意图。
图6是本发明中计算对象为滚动轴承内圈的薄片圆环单元的结构示意图。
具体实施方式
下面结合附图对本发明作进一步的说明,但并不因此将本发明限制在所述的实施例范围之中。
参见图1、图2、图3,一种滚动轴承套圈固有频率计算方法,包括以下步骤:
第一步,基本参数确定及频率精度设置:根据待计算的滚动轴承套圈的规格和材质的物理特性,获取滚动轴承套圈的几何尺寸参数和物理性能参数,并设置固有频率精度;其中,几何尺寸参数包括滚动轴承套圈的外径D、内径D1、套圈宽度B、沟道形状和位置,如沟道曲率半径r、沟底直径d1等;物理性能参数包括滚动轴承套圈包括弹性模量E和密度ρ;
第二步,薄片圆环单元划分:将滚动轴承套圈按宽度B方向划分成N个等厚度的薄片圆环单元,其中,N为大于或等于100的自然数;
第三步,基础数据计算:将单个薄片圆环单元的纵向截面视为矩形,按矩形截面的计算方法,计算每个薄片圆环单元的截面惯性矩Ii、质量Mi和横截面面积Ai,并计算中线半径Ri,该中线半径Ri等于薄片圆环单元内外半径之和的一半;
其中,薄片圆环单元内外半径的值由薄片圆环单元所处的具体位置确定;在轴承外圈频率计算中,参见图2并结合图4的深沟球轴承外圈外径为D不变,内径为D1或d1综合r;在轴承内圈频率计算中,参见图5并结合图6的深沟球轴承内圈内径为D1不变,外径为D或d1综合r;
第四步,薄片圆环单元固有频率计算:根据物理性能参数和基础数据计算结果,按圆环的固有频率计算方法,通过下式计算第i个薄片圆环单元的n阶径向弯曲频率ωni:
第五步,薄片圆环单元刚度计算:根据第n阶径向弯曲频率ωn与刚度Kn的关系,计算第i个薄片圆环单元的刚度Kni:
第六步,轴承套圈总刚度计算:按下式对N个薄片圆环单元的刚度累加求和计算轴承套圈总刚度Knz:
第七步,计算轴承套圈的第n阶固有频率:根据第n阶径向弯曲频率ωn与刚度Kn的关系,以轴承套圈的总刚度Knz,计算轴承套圈的第n阶固有频率ωn:
其中,
第八步,计算结果判断及确认:将本次与紧前一次的ωn值计算结果比较,当二者差值的绝对值大于固有频率精度的设定值时,则返回第二步,改变所述N的取值,重复进行第三步~第八步;当二者差值的绝对值小于或等于频率精度的设定值时,则进入第九步;
第九步,计算结果确认:结束计算,并将本次计算结果的ωn值确定为该轴承套圈的固有频率。
其中,在改变N的取值时,新一次取值按初始值的约2~4倍递增或按初始值的1/4~1/2倍递减;或者按紧前一次值的约2~5倍或按紧前一次值的约1/5~1/2倍取值,且N的最大取值不大于10000。以通过较少的计算次数获得准确的计算结果。在通常情况下,第一次取值接近下限值100,此时,第二次及之后的取值采用2~4倍的第一次取值递增,或按紧前一次值的约2~5倍取值;当第一次取值较高如1000以上时,第二次及之后的取值采用1/4~1/2倍的第一次取值递增,或按紧前一次值的约1/5~1/2倍取值。
为了进一步说明本方法的效果,现以圆环构件和实际深沟球轴承6306和6310外圈为分析对象,进行计算结果对比分析。
三者的结构参数如表1所示,深沟轴承外圈的结构示意图参见图2、图3和图4。对三者按本发明提供的计算方法与传统计算方法的计算结果进行对比,计算结果如表2所示。其中,f(2)、f(3)、f(4)分别为二阶频率、三阶频率和四阶频率,在计算径向弯曲频率时,频率代号f与ω代表相同的含义;分别计算分片数N的三次取值分别为102、502和1002的计算结果。表中原方法是指刊载于1998年第5期的“轴承”期刊第22-23页,由刘春浩、宋春磊、周晓辉等所发表的“球轴承径向弯曲振动的研究”中的所提供的传统计算方法。
表1三零件外形尺寸表(单位:mm)
表2三零件两种计算方法对比
表2的结果显示,对于圆环构件,两种方法的计算结果完全一致;对于深沟球轴承6306和6310的外圈,两种方法的计算差值均小于0.5%。上述结果证明了本发明提出方法的计算精度,且本发明提出方法在计算截面惯性矩和截面积上的效率更具有优势。
以上详细描述了本发明的较佳具体实施例。应当理解,本领域的普通技术人员无需创造性劳动就可以根据本发明的构思作出诸多修改和变化。因此,凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案,皆应在由权利要求书所确定的保护范围内。
Claims (6)
1.一种滚动轴承套圈固有频率计算方法,其特征在于,包括以下步骤:
第一步,基本参数确定及频率精度设置:根据待计算的滚动轴承套圈的规格和材质的物理特性,获取滚动轴承套圈的几何尺寸参数和物理性能参数,并设置固有频率精度;其中,几何尺寸参数包括滚动轴承套圈的外径D、内径D1、套圈宽度B、沟道形状和位置;物理性能参数包括滚动轴承套圈包括弹性模量E和密度ρ;
第二步,薄片圆环单元划分:将滚动轴承套圈按宽度B方向划分成N个等厚度的薄片圆环单元,其中,N为大于或等于100的自然数;
第三步,基础数据计算:将单个薄片圆环单元的纵向截面视为矩形,按矩形截面的计算方法,计算每个薄片圆环单元的截面惯性矩Ii、质量Mi和横截面面积Ai,并计算中线半径Ri,该中线半径Ri等于薄片圆环单元内外半径之和的一半;
第四步,薄片圆环单元固有频率计算:根据物理性能参数和基础数据计算结果,按圆环的固有频率计算方法,计算第i个薄片圆环单元的n阶径向弯曲频率ωni;
第五步,薄片圆环单元刚度计算:根据第n阶径向弯曲频率ωn与刚度Kn的关系计算第i个薄片圆环单元的刚度Kni;
第六步,轴承套圈总刚度计算:对N个薄片圆环单元的刚度累加求和计算轴承套圈总刚度Knz;
第七步,计算轴承套圈的第n阶固有频率:根据第n阶径向弯曲频率ωn与刚度Kn的关系,以轴承套圈的总刚度Knz计算轴承套圈的第n阶固有频率ωn:
第八步,计算结果判断及确认:将本次与紧前一次的ωn值计算结果比较,当二者差值的绝对值大于固有频率精度的设定值时,则返回第二步,改变所述N的取值,重复进行第三步~第八步;当二者差值的绝对值小于或等于频率精度的设定值时,则进入第九步;
第九步,计算结果确认:结束计算,并将本次计算结果的ωn值确定为该滚动轴承套圈的固有频率。
2.根据权利要求1所述的滚动轴承套圈固有频率计算方法,其特征在于,所述第i个薄片圆环单元的n阶径向弯曲频率ωni按下式计算:
3.根据权利要求1所述的滚动轴承套圈固有频率计算方法,其特征在于,所述第i个薄片圆环单元的刚度Kni按下式计算:
4.根据权利要求1所述的滚动轴承套圈固有频率计算方法,其特征在于,所述滚动轴承套圈总刚度Knz按下式计算:
5.根据权利要求1所述的滚动轴承套圈固有频率计算方法,其特征在于,所述滚动轴承套圈的第n阶固有频率按下式计算:
其中,
6.根据权利要求1~5中任意一项所述的滚动轴承套圈固有频率计算方法,其特征在于,在改变所述N的取值时,新一次取值按初始值的2~4倍递增或按初始值的1/4~1/2倍递减;或者按紧前一次值的2~5倍或按紧前一次值的1/5~1/2取值,且N的最大取值不大于10000。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610969618.8A CN106547980B (zh) | 2016-11-01 | 2016-11-01 | 一种滚动轴承套圈固有频率计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610969618.8A CN106547980B (zh) | 2016-11-01 | 2016-11-01 | 一种滚动轴承套圈固有频率计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106547980A CN106547980A (zh) | 2017-03-29 |
CN106547980B true CN106547980B (zh) | 2019-07-16 |
Family
ID=58394066
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610969618.8A Active CN106547980B (zh) | 2016-11-01 | 2016-11-01 | 一种滚动轴承套圈固有频率计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106547980B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113378319B (zh) * | 2021-06-24 | 2022-04-12 | 宁波蓝海量子精工轴承制造有限公司 | 一种双端面不对称轴承套圈端面磨削余量的设计方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105243215A (zh) * | 2015-10-12 | 2016-01-13 | 西安交通大学 | 基于参数化有限元模型的电主轴结构优化方法 |
CN105653890A (zh) * | 2016-04-07 | 2016-06-08 | 东北大学 | 一种基于轴向载荷的提升机轴承疲劳寿命模型 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160153426A1 (en) * | 2013-07-04 | 2016-06-02 | Jens GROENAGER | Bearing and gear unit for wind turbines |
-
2016
- 2016-11-01 CN CN201610969618.8A patent/CN106547980B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105243215A (zh) * | 2015-10-12 | 2016-01-13 | 西安交通大学 | 基于参数化有限元模型的电主轴结构优化方法 |
CN105653890A (zh) * | 2016-04-07 | 2016-06-08 | 东北大学 | 一种基于轴向载荷的提升机轴承疲劳寿命模型 |
Non-Patent Citations (2)
Title |
---|
Multi-DOF rotor model based measurement of stiffness and damping for active magnetic bearing using multi-frequency excitation;Jiang K;《Mechanical Systems & Signal Processing》;20151231;第358-374页 |
滚动轴承缺陷非线性激励机理与建模研究;刘静;《中国博士学位论文全文数据库 工程科技Ⅱ辑》;20141115;第C029-10页 |
Also Published As
Publication number | Publication date |
---|---|
CN106547980A (zh) | 2017-03-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111797567A (zh) | 一种基于深度学习网络的轴承故障分类方法及系统 | |
CN104484489B (zh) | 一种点蚀损伤圆柱壳的四边形有限元网格自动生成方法 | |
CN103047959B (zh) | 一种面向精密装配的基于熵理论的平面形状误差检测方法 | |
CN102346745A (zh) | 单词的用户行为数的预测方法和装置 | |
CN104199942B (zh) | 一种Hadoop平台时序数据增量计算方法及系统 | |
Cao et al. | An improvement of time-reassigned synchrosqueezing transform algorithm and its application in mechanical fault diagnosis | |
Li et al. | Automated and adaptive ridge extraction for rotating machinery fault detection | |
CN104504120A (zh) | 一种相似图片检测方法及装置 | |
WO2019200624A1 (zh) | 一种泵流动诱导振动性能综合评价方法 | |
CN106814029A (zh) | 一种多带积分盘式摩擦副摩擦性能预测方法 | |
CN106547980B (zh) | 一种滚动轴承套圈固有频率计算方法 | |
CN105004274A (zh) | 一种基于三维视觉的圆柱面半径测量方法 | |
CN105893327A (zh) | 基于fft的深沟球轴承和角接触球轴承弹性变形快速计算方法 | |
CN102269803B (zh) | 基于时间延迟的离散频谱低频成分的校正方法 | |
CN109409026B (zh) | 电机效率优化方法、装置及电机设计参数确定方法 | |
Wang et al. | Novel method for evaluating surface roughness by grey dynamic filtering | |
CN104134013A (zh) | 一种风力机叶片模态分析方法 | |
JP6090000B2 (ja) | 周波数解析装置 | |
KR20150142241A (ko) | 구조물의 피로수명을 고려한 위상최적설계방법 및 그 방법을 수행하기 위한 프로그램이 기록된 기록매체 | |
DE102021200210A1 (de) | Frühanzeichenerfassungsvorrichtung und Frühanzeichenerfassungsverfahren | |
CN105740225A (zh) | 一种融合句子局部上下文与文档领域信息的词义消歧方法 | |
Chen et al. | An improved local mean decomposition method and its application for fault diagnosis of reciprocating compressor | |
CN104135510A (zh) | 基于模式匹配的分布式计算环境性能预测方法及系统 | |
CN101988964A (zh) | 高效二次雷达目标周期相关处理工作方法 | |
CN105243279B (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 |