CN103810354B - 圆柱滚子轴承对数修形曲线的优化设计方法 - Google Patents

圆柱滚子轴承对数修形曲线的优化设计方法 Download PDF

Info

Publication number
CN103810354B
CN103810354B CN201410088720.8A CN201410088720A CN103810354B CN 103810354 B CN103810354 B CN 103810354B CN 201410088720 A CN201410088720 A CN 201410088720A CN 103810354 B CN103810354 B CN 103810354B
Authority
CN
China
Prior art keywords
logarithm
modification curve
contact
curve
equation
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.)
Expired - Fee Related
Application number
CN201410088720.8A
Other languages
English (en)
Other versions
CN103810354A (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.)
Dalian Jiaotong University
Original Assignee
Dalian Jiaotong 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 Dalian Jiaotong University filed Critical Dalian Jiaotong University
Priority to CN201410088720.8A priority Critical patent/CN103810354B/zh
Publication of CN103810354A publication Critical patent/CN103810354A/zh
Application granted granted Critical
Publication of CN103810354B publication Critical patent/CN103810354B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Rolling Contact Bearings (AREA)

Abstract

本发明公开了一种圆柱滚子轴承对数修形曲线的优化设计方法,首先将对数修形曲线方程简化变形,针对简化变形的对数修形曲线方程进行处理,使优化设计转化为对简化变形的对数修形曲线方程式中一个参数的优化问题,然后,应用有限元物理仿真方法获得的接触应力的分布规律曲线图,通过对多组对数曲线方程模型的计算分析和接触应力分布规律曲线图对比,优化对数修形曲线方程公式中一个参数,也就是对对数曲线的斜率和修形的凸度量进行优化设计,从而设计出与所受载荷和被修形元件几何形状相对应的最佳修形曲线。

Description

圆柱滚子轴承对数修形曲线的优化设计方法
技术领域
本发明属于轴承技术领域,尤其涉及圆柱滚子轴承设计方法。
背景技术
修形滚子轴承正在逐渐在很多重要领域取代传统的直母线滚子轴承。传统的直母线滚子轴承的滚动体与滚道间的早期接触疲劳点蚀常常发生在滚子或滚道上靠近滚子端部的区域,这是因为直母线滚子轴承在受载后滚动体两端存在边界应力集中﹐即“边缘效应”。“边缘效应”的产生使轴承的疲劳寿命大大降低,因为研究表明轴承的寿命与应力的7次方成反比。为了克服这种“边缘效应”,人们进行了大量的理论分析和实验研究。早在十九世纪30年代末Lundberg就提出了母线修形的基本理论,直至二十世纪60年代SKF轴承公司进一步发展了滚子轴承的修形技术。通过使用特殊的滚子外廓曲面可以避免或降低滚动体和内外圈接触引起的边界应力集中,目前,工程中采用的修形曲线主要有:圆弧曲线;直线两端加圆弧——滚子母线中部为直线,两端是修缘圆弧;对数曲线等。对数曲线修形被公认为最佳的修形曲线。现有对数曲线修形方法的不足是:同样采用对数曲线修形,对数曲线方程的具体参数不同,修形效果不同,影响了轴承的承载能力和使用寿命;同时,对数修形曲线存在优化设计问题;对数曲线方程具体参数的确定和修形效果的验证难以解决。
发明内容
本发明的目的在于提出一种圆柱滚子轴承对数修形曲线的优化设计方法,克服现有技术的不足。
本发明的技术方案是:圆柱滚子轴承对数修形曲线的优化设计方法,包括根据接触力学原理建立圆柱滚子轴承对数修形曲线方程和对该圆柱滚子轴承对数修形曲线方程的处理,根据接触力学原理建立圆柱滚子轴承对数修形曲线方程为:
y = k 0 P 2 l l n { 1 - ( x / l ) 2 }
式中
k 0 = 1 - γ 1 2 πE 1 + 1 - γ 2 2 πE 2
上述式中,x是以滚子接触线的中点为原点沿接触线长度方向的坐标值,y是坐标值x所对应的径向坐标值,l为接触半长,P为接触全长2l上的载荷,γ1和γ2分别为接触元件1和接触元件2的泊松比,E1和E2分别为接触元件1和接触元件2的弹性模量,其特征在于,所述对圆柱滚子轴承对数修形曲线方程的处理包括以下步骤:
1)将方程
y = k 0 P 2 l l n { 1 - ( x / l ) 2 }
简写成
y = k 2 l n ( 1 1 - k 1 x 2 )
式中k1和k2是修形参数,其中k2由轴承的材料、所受载荷和几何尺寸确定,
2)根据几何尺寸关系初步确定k1
取不同的k1值,分别命名为k1n,n为正整数,以方程式
y = k 2 l n ( 1 1 - k 1 n x 2 )
建立有限元模型,用通用商业化的有限元分析软件对修形元件进行物理仿真计算,根据有限元分析结果优选一个最佳修形参数k1值;
3)该最佳修形参数k1值对应的对数修形曲线方程式
y = k 2 l n ( 1 1 - k 1 x 2 )
即是最佳的对数修形曲线方程。
本发明的方法是将对数修形曲线的优化设计问题转化为对数修形曲线方程式中一个参数的优化,然后,应用有限元物理仿真方法获得的接触应力的分布规律曲线图,通过对多组对数曲线方程模型的计算分析和接触应力分布规律曲线图对比,优化对数修形曲线方程式中一个参数,也就是对对数曲线的斜率和修形的凸度量进行优化设计,从而设计出与所受载荷和被修形元件几何形状相对应的最佳修形曲线。
附图说明
图1为某轴承滚子的对数修形曲线族,即,相同的k2,不同的k1时的对数曲线族。
图2为未修形和4个修形轴承在载荷为20%Cr时,滚子与外圈接触线上最大接触应力沿轴向分布曲线比较图。
图3为未修形和4个修形轴承在载荷为20%Cr时,滚子与内圈接触线上最大接触应力沿轴向分布曲线比较图。
具体实施方式
以下结合附图和具体实施例对本发明作进一步说明。
根据接触力学原理建立圆柱滚子轴承对数修形曲线方程为:
y = k 0 P 2 l l n { 1 - ( x / l ) 2 }
式中
k 0 = 1 - γ 1 2 πE 1 + 1 - γ 2 2 πE 2
上述式中,x是以滚子接触线的中点为原点沿接触线长度方向的坐标值,y是坐标值x所对应的径向坐标值,l为接触半长,P为接触全长2l上的载荷,γ1和γ2分别为接触元件1和接触元件2的泊松比,E1和E2分别为接触元件1和接触元件2的弹性模量,对圆柱滚子轴承对数修形曲线方程作以下处理:
将方程
y = k 0 P 2 l l n { 1 - ( x / l ) 2 }
简写成
y = k 2 l n ( 1 1 - k 1 x 2 )
式中k1和k2是修形参数,其中k2由轴承的材料、所受载荷和几何尺寸确定。
取不同的k1值,分别命名为k1n,n为正整数,以方程式
y = k 2 l n ( 1 1 - k 1 n x 2 )
建立有限元模型,用通用商业化的有限元分析软件对修形元件进行物理仿真计算,根据有限元分析结果优选最佳修形参数k1值;
3)该最佳修形参数k1值对应的对数修形曲线方程式
y = k 2 l n ( 1 1 - k 1 x 2 )
即是最佳的对数修形曲线方程。
具体实施例:
某圆柱滚子轴承的修形设计。轴承的基本参数为:额定动载荷Cr=860(kN),额定静载荷C0r=1350(kN);轴承内圈内径130,内圈外径150,单列内圈宽62;轴承外圈内径202,外圈外径220单列外圈宽60;单列滚子数18,滚子直径26,滚子长度40,滚子两端有0.6±0.2的圆角,圆角与滚子表面的对数曲线圆滑过渡。
图2和图3是用有限元分析软件对该轴承滚子几个不同k1值的对数修形曲线处理的修形效果图,图2是未修形和修形(4个)轴承在载荷为20%Cr(Cr额定动负荷)时,滚子与外圈接触线上最大接触应力沿轴向分布曲线比较图;图3是未修形和修形(4个)轴承在载荷为20%Cr时,滚子与内圈接触线上最大接触应力沿轴向分布曲线比较图。根据以上图确定最佳的k1值,也就是确定最佳的对数修形曲线,最佳对数曲线方程为:k1=0.002657,k2=0.001585,与其对应的滚子与外圈滚道接触应力沿滚子轴向的分布规律曲线为图2中的CP-lg1-e;滚子与内圈滚道接触应力沿滚子轴向的分布规律曲线为图3中的CP-lg1-i。图2和图3中CP-line-e和CP-line-i为未修形的直母线滚子与内、外圈道之间接触应力沿轴向的分布规律曲线,其他三条曲线是不同对数曲线修形的滚子与内、外圈道之间接触应力沿轴向的分布规律曲线。由图可见,最佳对数曲线修形的滚子与内、外圈道之间接触应力沿轴向的分布规律曲线最均匀,最大接触应力最小。

Claims (1)

1.圆柱滚子轴承对数修形曲线的优化设计方法,包括根据接触力学原理建立圆柱滚子轴承对数修形曲线方程和对该圆柱滚子轴承对数修形曲线方程的处理,根据接触力学原理建立圆柱滚子轴承对数修形曲线方程为:
y = k 0 P 2 l l n { 1 - ( x / l ) 2 }
式中
k 0 = 1 - γ 1 2 πE 1 + 1 - γ 2 2 πE 2
上述式中,x是以滚子接触线的中点为原点沿接触线长度方向的坐标值,y是坐标值x所对应的径向坐标值,l为接触半长,P为接触全长2l上的载荷,γ1和γ2分别为接触元件1和接触元件2的泊松比,E1和E2分别为接触元件1和接触元件2的弹性模量,其特征在于,所述对圆柱滚子轴承对数修形曲线方程的处理包括以下步骤:
2)将方程
y = k 0 P 2 l l n { 1 - ( x / l ) 2 }
简写成
y = k 2 l n ( 1 1 - k 1 x 2 )
式中k1和k2是修形参数,其中k2由轴承的材料、所受载荷和几何尺寸确定,
2)根据几何尺寸关系初步确定k1
根据取不同的k1值,分别命名为k1n,n为正整数,以方程式
y = k 2 l n ( 1 1 - k 1 n x 2 )
建立有限元模型,用通用商业化的有限元分析软件对修形元件进行物理仿真计算,根据有限元分析结果优选出一个最佳修形参数k1值;
3)该最佳修形参数k1值对应的对数修形曲线方程式
y = k 2 l n ( 1 1 - k 1 x 2 )
即是最佳的对数修形曲线方程。
CN201410088720.8A 2014-03-11 2014-03-11 圆柱滚子轴承对数修形曲线的优化设计方法 Expired - Fee Related CN103810354B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410088720.8A CN103810354B (zh) 2014-03-11 2014-03-11 圆柱滚子轴承对数修形曲线的优化设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410088720.8A CN103810354B (zh) 2014-03-11 2014-03-11 圆柱滚子轴承对数修形曲线的优化设计方法

Publications (2)

Publication Number Publication Date
CN103810354A CN103810354A (zh) 2014-05-21
CN103810354B true CN103810354B (zh) 2016-08-17

Family

ID=50707118

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410088720.8A Expired - Fee Related CN103810354B (zh) 2014-03-11 2014-03-11 圆柱滚子轴承对数修形曲线的优化设计方法

Country Status (1)

Country Link
CN (1) CN103810354B (zh)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104636596B (zh) * 2014-12-26 2017-12-05 中国北方车辆研究所 一种特定载荷下圆柱滚子轴承非对称修形方法
CN106991202A (zh) * 2016-01-21 2017-07-28 宝山钢铁股份有限公司 一种热轧矫直机支承辊组的辊面曲线优化方法
CN105893327B (zh) * 2016-03-31 2018-06-05 重庆大学 基于fft的深沟球轴承和角接触球轴承弹性变形快速计算方法
CN106438707A (zh) * 2016-11-01 2017-02-22 重庆大学 一种深穴蜂窝小孔系空心圆柱滚子轴承
CN107358011B (zh) * 2017-09-12 2020-12-01 大连交通大学 基于圆柱滚子轴承载荷分布计算的轴承参数优化方法
CN107688716B (zh) * 2017-09-12 2021-02-26 大连交通大学 基于载荷分布和疲劳寿命的空心圆柱滚子轴承参数优化方法
CN110502765B (zh) * 2018-05-17 2023-04-07 河南科技大学 一种圆锥滚子轴承的修形方法及滚子轴承
CN108846249A (zh) * 2018-08-29 2018-11-20 中国航发哈尔滨轴承有限公司 一种面向高速重载工况下的圆柱滚子的设计方法
CN109871662B (zh) * 2019-03-27 2023-06-09 大连交通大学 高分子弹性元件缓冲器用曲面隔片内凹曲面的设计方法
CN111475895B (zh) * 2020-04-10 2023-03-24 洛阳Lyc轴承有限公司 一种球面滚子的端部圆弧修形方法
CN112580218A (zh) * 2020-12-25 2021-03-30 中国航发哈尔滨轴承有限公司 一种圆柱滚子轴承母线凸度优化设计方法
CN112727929A (zh) * 2020-12-31 2021-04-30 甘肃海林中科科技股份有限公司 一种特定工况下圆锥滚子轴承滚子表面凸度值的设计方法
CN113468695B (zh) * 2021-07-20 2023-08-25 人本股份有限公司 滚子的凸度优化设计方法
CN115139158B (zh) * 2022-06-22 2023-10-10 洛阳理工学院 一种用于双列调心球面滚子轴承的滚子修型方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6327552B2 (en) * 1999-12-28 2001-12-04 Intel Corporation Method and system for determining optimal delay allocation to datapath blocks based on area-delay and power-delay curves
CN201696494U (zh) * 2010-04-16 2011-01-05 浙江大学 一种对数修形孔的滑动轴承
CN201934514U (zh) * 2010-11-09 2011-08-17 大连交通大学 空心滚子轴承

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6327552B2 (en) * 1999-12-28 2001-12-04 Intel Corporation Method and system for determining optimal delay allocation to datapath blocks based on area-delay and power-delay curves
CN201696494U (zh) * 2010-04-16 2011-01-05 浙江大学 一种对数修形孔的滑动轴承
CN201934514U (zh) * 2010-11-09 2011-08-17 大连交通大学 空心滚子轴承

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
圆柱滚子轴承滚动体新的修形曲线;魏延刚等;《机械设计与制造》;20030430;第81-83页 *
高速客车轴箱轴承的偏载分析及轴承滚子非对称修形;魏延刚等;《润滑与密封》;20020228(第2期);第9-11页 *

Also Published As

Publication number Publication date
CN103810354A (zh) 2014-05-21

Similar Documents

Publication Publication Date Title
CN103810354B (zh) 圆柱滚子轴承对数修形曲线的优化设计方法
Kania Modelling of rollers in calculation of slewing bearing with the use of finite elements
CN107545107B (zh) 一种用于梯形截面轴承环轧制的环坯及环坯设计方法
Qian et al. Finite element modeling of power spinning of thin-walled shell with hoop inner rib
CN102644663A (zh) 圆柱滚子多圆弧变曲率凸型工程模拟方法
CN107688716A (zh) 基于载荷分布和疲劳寿命的空心圆柱滚子轴承参数优化方法
CN107563082B (zh) 基于圆柱滚子轴承接触变形量和载荷分布的轴承参数优化方法
CN107563081A (zh) 高速空心圆柱滚子轴承的参数优化方法
CN110502765B (zh) 一种圆锥滚子轴承的修形方法及滚子轴承
Qiuying et al. Spatiotemporal differentiation and the factors influencing eco-efficiency in China
QIAN et al. FE analysis for radial spread behavior in three-roll cross rolling with small-hole and deep-groove ring
CN106402162A (zh) 非圆滚道滚动轴承
CN107391893B (zh) 基于空心圆柱滚子接触变形量和载荷分布的轴承参数优化方法
CN204025307U (zh) 一种小排量内燃机的轻量化曲轴
CN205745093U (zh) 双频阻尼器
CN104268311B (zh) 汽车座椅导轨直线滚子结合部间隙评价方法
Zhao et al. Influence analysis of key design parameters on the roller safety factor of three-row and four-column roller slewing bearings based on analytical method and FEM
CN206344505U (zh) 一种绿色造纸技术用压花胶辊
CN105160179B (zh) 高速轨道车辆二系横向悬置最优阻尼比的解析计算方法
CN110826262A (zh) 一种高强双角钢组合截面构件承载力计算方法
CN201405231Y (zh) 一种对数凸度圆柱滚子精研油石
Han et al. Analysis of Residual Stress and Deformation of Rolling Strengthen Crankshaft Fillet
CN203335663U (zh) 一种用于十字轴万向节的圆弧凸度滚针
CN113468695A (zh) 滚子的凸度优化设计方法
CN107358011B (zh) 基于圆柱滚子轴承载荷分布计算的轴承参数优化方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160817

Termination date: 20210311

CF01 Termination of patent right due to non-payment of annual fee