CN106326570A - 球、锥混合双列轮毂轴承最佳预紧量的确定方法及装置 - Google Patents
球、锥混合双列轮毂轴承最佳预紧量的确定方法及装置 Download PDFInfo
- Publication number
- CN106326570A CN106326570A CN201610736007.9A CN201610736007A CN106326570A CN 106326570 A CN106326570 A CN 106326570A CN 201610736007 A CN201610736007 A CN 201610736007A CN 106326570 A CN106326570 A CN 106326570A
- Authority
- CN
- China
- Prior art keywords
- hub
- ball
- bearing
- load
- row
- 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
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
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F16—ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
- F16C—SHAFTS; FLEXIBLE SHAFTS; ELEMENTS OR CRANKSHAFT MECHANISMS; ROTARY BODIES OTHER THAN GEARING ELEMENTS; BEARINGS
- F16C43/00—Assembling bearings
- F16C43/04—Assembling rolling-contact bearings
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Mechanical Engineering (AREA)
- Rolling Contact Bearings (AREA)
Abstract
本发明涉及一种球、锥混合双列轮毂轴承最佳预紧量的确定方法及装置,属于轴承设计技术领域。本发明通过制定球、锥混合双列轮毂轴承单元载荷谱,根据球、锥混合双列轮毂轴承的力学模型计算径向位移δr、轴向位移δa和角位移θ,并以载荷谱为已知条件,以预紧量作为分析变量得到轮毂轴承单元相应的单列套圈寿命和综合寿命,建立预紧量与轮毂轴承单元寿命关系,以轮毂轴承单元单列套圈和或整个轮毂轴承单元寿命最长所对应的预紧量为最佳预紧量。本发明提供的基于载荷谱下对球、锥混合双列轮毂轴承最佳预紧量的确定方法,能够精确得出球、锥混合双列轮毂轴承单元装配中的最佳预紧量。
Description
技术领域
本发明涉及一种球、锥混合双列轮毂轴承最佳预紧量的确定方法及装置,属于轴承设计技术领域。
背景技术
轮毂轴承单元是汽车中最为关键的零部件之一,是汽车车轴处用来承重和为轮毂转动提供精确引导的零部件,一般采用双列结构。当下主流车用轮毂轴承单元分为两种:双列角接触球轮毂轴承单元和双列圆锥滚子轮毂轴承单元。在如今家用汽车制造趋向提高舒适性、多功能性并力争扩大汽车内部空间而致使家用汽车重量随之增加,从而对轮毂轴承承载能力有更高要求的背景下,现今出现了一种新型轮毂轴承结构:球、锥混合双列轮毂轴承,即一列钢球,一列圆锥滚子。由于双列轮毂轴承单元中,多数情况下只有靠近汽车转向轴侧列轴承承受主要载荷,另一列轴承受载时间较短,主要承受次要载荷,所以将球、锥混合双列轮毂轴承的圆锥滚子列安装于靠近转向轴侧,承受主要载荷。球、锥混合双列轮毂轴承兼具了双列角接触球轴承结构紧凑、摩擦力矩小、制造成本低和双列圆锥滚子轴承承载性能较强、轴承刚度大等两者共同优点,更加适合应用于质量较重、对轮毂轴承承载性能要求较高的越野汽车上,很可能会成为未来轮毂轴承的主要发展趋势。
预紧量作为轮毂轴承单元装配过程中最重要的因素之一,直接影响轮毂轴承的使用寿命和振动噪声及旋转精度等工作性能。当预紧量过大时会使滚道接触应力急剧增大而产生早期疲劳,当预紧量过小时在轴承受载后容易造成松动而形成振动,所以确定出轮毂轴承单元的最佳装配过盈量,对于轮毂轴承单元延长使用寿命和发挥最佳工作性能都具有重要意义。而针对新型球、锥混合双列轮毂轴承单元,虽然提出了其设计结构,至今还没有具体的力学建模和承载性能精确计算方法,更没有能精确确定装配过程中其最佳预紧量的方法和装置。
发明内容
本发明的目的是提供一种球、锥混合双列轮毂轴承最佳预紧量的确定方法及装置,以解决难以精确确定球、锥混合双列轮毂轴装配过程中最佳预紧量的问题。
本发明为解决上述技术问题而提供一种球、锥混合双列轮毂轴承最佳预紧量的确定方法,该方法包括以下步骤:
1)根据轮胎受力公式计算不同行驶情况下球、锥混合双列轮毂轴承单元所受的径向载荷Frn、轴向载荷Fan和倾覆力矩Mn,并根据各行驶情况所对应的时间分配率确定球、锥混合双列轮毂轴承单元相对应的载荷谱;
2)根据球、锥混合双列轮毂轴承的力学模型建立以其内圈中心径向位移δr、轴向位移δa和角位移θ三个未知量构成的内圈力学平衡方程,并以步骤1)中所确定的径向载荷Frn、轴向载荷Fan和颠覆力矩Mn作为外部负载求解该平衡方程,以得到球、锥混合双列轮毂轴承单元内圈对应的径向位移量δr、轴向位移量δa和角位移量θ;
3)基于所建立的载荷谱计算平均当量动载荷,根据轴承额定寿命计算球、锥混合双列轮毂轴承单元钢球基本额定寿命L10_1和或圆锥滚子基本额定寿命L10_2;
4)以每个轴向预紧量值作为轴向游隙,建立预紧量与轮毂轴承单元寿命关系,以轮毂轴承单元单列套圈和或整个轮毂轴承单元寿命最长所对应的预紧量为最佳预紧量。
所述步骤4)中平均当量动载荷的计算公式为:
其中Qei为第i种行驶情况下轴承的当量动载荷,ni为第i种行驶情况下轮毂轴承转速,qi为第i种行驶情况轮毂轴承单元所承受载荷的对应时间分配率,下标n为行驶情况的种类。
所述步骤4)中的综合寿命L10为:
其中e1为点接触,e2为线接触,e为e1和e2中的较小者。
步骤2)所建立的平衡方程为:
Q1i为滚道对钢球的法向载荷,α1ψi表示受载后任意位置钢球的接触角,ψi表示钢球的位置角,Z1表示钢球数,dm1为钢球节圆直径,de1为轴承内圈中心O与钢球节圆中心O1’之间距离的2倍,Q2i为滚子与外圈接触载荷,αe为圆锥滚子与外圈接触角,φi为第i个圆锥滚子的位置角,Z2表示滚子数,dm2为滚子节圆直径,de2为轴承内圈中心O与滚子节圆中心O2’之间距离的2倍。
所述滚道对钢球的法向载荷和滚子与外圈接触载荷的计算公式分别为:
Q1i=Kn1[max(δ1i,0)]3/2
Q2i=Kn2[max(δ2i,0)]10/9
δ2i=(δr-R22θ)cosφi cosαe+(δ0-0.5dm2θcosφi-δa)sinαe
其中Kn1为滚动体与内外圈之间的总负载-变形常数,α0为原始接触角,Ri为钢球内滚道沟曲率中心所在圆的半径,A为内外圈原始沟曲率中心距,R22为轴承轴承内圈中心O与钢球节圆中心O2’之间的距离,δ0为轴向游隙即轴向预紧量。
所述步骤1)中径向载荷Frn、轴向载荷Fan和倾覆力矩Mn的计算公式分别为:
Mn=Fmso+FanRs
其中fw为载荷冲击系数,W为汽车重量,g为重力加速度,ag为汽车侧向加速度率,h为汽车质心高度,l为前轴轮间距,s0为径向受力偏心距,Rs为汽车轮胎半径,tr为汽车前轴轮间距。
步骤2)中的平衡方程采用Newton-Raphson迭代法进行求解。
本发明还提供了一种球、锥混合双列轮毂轴承最佳预紧量的确定装置,该装置包括载荷确定模块、位移量确定模块、寿命计算模块和最佳预紧量确定模块,
所述载荷确定模块用于根据轮胎受力公式计算不同行驶情况下球、锥混合双列轮毂轴承单元所受的径向载荷Frn、轴向载荷Fan和倾覆力矩Mn,并根据各行驶情况所对应的时间分配率确定球、锥混合双列轮毂轴承单元相对应的载荷谱;
所述位移量确定模块用于根据球、锥混合双列轮毂轴承的力学模型建立以其内圈中心径向位移δr、轴向位移δa和角位移θ三个未知量构成的内圈力学平衡方程,并以所确定的径向载荷Frn、轴向载荷Fan和颠覆力矩Mn作为外部负载求解该平衡方程,以得到球、锥混合双列轮毂轴承单元内圈对应的径向位移量δr、轴向位移量δa和角位移量θ;
所述的寿命计算模块用于根据于所建立的载荷谱计算平均当量动载荷,根据轴承额定寿命计算球、锥混合双列轮毂轴承单元钢球基本额定寿命L10_1和或圆锥滚子基本额定寿命L10_2;
所述最佳预紧量确定模块用于以每个轴向预紧量值作为轴向游隙,建立预紧量与轮毂轴承单元寿命关系,以轮毂轴承单元单列套圈和或整个轮毂轴承单元寿命最长所对应的预紧量为最佳预紧量。
所述寿命计算模块所采用的平均当量动载荷的计算公式为:
其中Qei为第i种行驶情况下轴承的当量动载荷,ni为第i种行驶情况下轮毂轴承转速,qi为第i种行驶情况轮毂轴承单元所承受载荷的对应时间分配率,下标n为行驶情况的种类。
所述位移量确定模块中建立的平衡方程为:
Q1i为滚道对钢球的法向载荷,α1ψi表示受载后任意位置钢球的接触角,ψi表示钢球的位置角,Z1表示钢球数,dm1为钢球节圆直径,de1为轴承内圈中心O与钢球节圆中心O1’之间距离的2倍,Q2i为滚子与外圈接触载荷,αe为圆锥滚子与外圈接触角,φi为第i个圆锥滚子的位置角,Z2表示滚子数,dm2为滚子节圆直径,de2为轴承内圈中心O与滚子节圆中心O2’之间距离的2倍。
本发明的有益效果是:本发明通过制定球、锥混合双列轮毂轴承单元载荷谱,根据球、锥混合双列轮毂轴承的力学模型计算径向位移δr、轴向位移δa和角位移θ,并以载荷谱为已知条件,以预紧量作为分析变量得到轮毂轴承单元相应的单列套圈寿命和综合寿命,建立预紧量与轮毂轴承单元寿命关系,以轮毂轴承单元单列套圈和或整个轮毂轴承单元寿命最长所对应的预紧量为最佳预紧量。本发明提供的基于载荷谱下对球、锥混合双列轮毂轴承最佳预紧量的确定方法,能够精确得出球、锥混合双列轮毂轴承单元装配中的最佳预紧量,对此新型轮毂轴承单元使用寿命的延长和最佳工作性能的发挥具有重要意义,为球、锥混合双列轮毂轴承单元的装配制造和使用提供了有力的科学依据。
附图说明
图1是球、锥混合双列轮毂轴承结构及内圈受力示意图;
图2是转角θ对右列轴承产生的位移分量示意图;
图3是径向位移和轴向位移外滚道投影示意图;
图4是本发明实例中越野车行驶载荷谱示意图;
图5是本发明实例中左列钢球套圈寿命-预紧量关系图;
图6是本发明实例中右列圆锥滚子套圈寿命-预紧量关系图;
图7是本发明实例中轮毂轴承单元的综合寿命-预紧量关系图。
具体实施方式
下面结合附图对本发明的具体实施方式做进一步的说明。
本发明球、锥混合双列轮毂轴承最佳预紧量的确定方法的实施例
本发明通过制定球、锥混合双列轮毂轴承单元载荷谱,根据球、锥混合双列轮毂轴承的力学模型计算径向位移δr、轴向位移δa和角位移θ,并以载荷谱为已知条件,以预紧量作为分析变量得到轮毂轴承单元相应的单列套圈寿命和综合寿命,建立预紧量与轮毂轴承单元寿命关系,以轮毂轴承单元单列套圈和或整个轮毂轴承单元寿命最长所对应的预紧量为最佳预紧量。该方法的具体实施过程如下:
1.制定球、锥混合双列轮毂轴承单元载荷谱
载荷谱的制定便于更好地分析轮毂轴承的性能,为轮毂轴承的耐久性分析计算提供了更加科学合理的依据。由于汽车行驶环境的多样性和不确定性,对于球、锥混合双列轮毂轴承单元,本实施例中的载荷谱包含四种情况:良好路面直行、左转、颠簸路面直行、右转。
借鉴国外轮毂轴承耐久性试验标准周期,取一个实验周期为100s,基于中国路面基本情况,取其直线行驶时间占用率为94%,转弯时间占用率为6%。其中良好路面直行时间占用率和颠簸路面直行时间占用率分别为50%和44%;左转和右转时间占有率都为3%。
对汽车驾驶室侧(即汽车前轴左侧)轮毂轴承进行分析,其轮胎径向力Frn和轮胎轴向力Fan计算公式为:
其中fw为载荷冲击系数,在良好路面行驶取值为1.0、颠簸路面行驶取值为1.3;W为汽车质量;g为重力加速度,取值为9.8m/s2;ag为汽车侧向加速度率(即侧向加速度与重力加速的的比值),直行时取值为0、左转时取值为0.55、右转时取值为-0.55;h为汽车质心高度,l为前轴轮间距。
倾覆力矩计算公式Mn为:
Mn=Frnso+FanRs (3)
其中so为径向受力偏心距,Rs为汽车轮胎半径。
当已知汽车相关参数后,按照不同行驶路况下对应的载荷冲击系数fw和侧向加速度率ag,可分别计算出良好路面直行、左转、颠簸路面直行、右转等四种行驶情况下,球、锥混合双列轮毂轴承单元所受的径向载荷Frn、轴向载荷Fan和倾覆力矩Mn,(n取1、2、3、4,对应于四种工况条件)。
以时间为横坐标,求取的力和力矩值为纵坐标,取一个实验周期为100s,按照汽车四种行驶情况下所对应的时间占用率,即可获得球、锥混合双列轮毂轴承单元相对应的载荷谱。
2.建立球、锥混合双列轮毂轴承单元力学模型
假设球、锥混合双列轮毂轴承单元其外圈固定,当内圈承受径向力Fr、轴向力Fa和颠覆力矩M联合作用时,其径向位移为δr、轴向位移为δa、角位移为θ,本文定义对于每列轴承使滚动体载荷增加的位移为正位移,反之则为负位移。球、锥混合双列轮毂轴承单元结构及受力示意图如图1所示。并规定左列钢球相关表达式所用下标为1,右列圆锥滚子相关表达式所用下标为2,用i表示滚动体(左列钢球、右列圆锥滚子)。
左列钢球与滚道之间为点接触,α0为原始接触角,受载后在不同位置钢球的接触角会发生不同的变化,ψi表示钢球的位置角,用α1ψi表示受载后任意位置钢球的接触角,则:
A=(fi+fe-1)Dw1 (7)
其中,S1ψi为受负载位移之后任意位置的内外圈沟曲率中心距,A为内外圈原始沟曲率中心距,Ri为钢球内滚道沟曲率中心所在圆的半径。
任意位置钢球的变形量δ1i等于受负载之后内外圈沟曲率中心距与原始沟曲率中心距之差:
δ1i=S1ψi-A (9)
钢球为点接触,所以滚道对钢球的法向载荷与钢球的接触变形的关系为:
Q1i=Kn1[max(δ1i,0)]3/2 (10)
其中Kn1为滚动体与内外圈之间的总负载-变形常数。左列钢球数为Z1,则角接触钢球接触载荷径向分量总和为:
角接触钢球接触载荷轴向分量总和为:
由于球-锥混合轴承为双列轴承,轴承内圈中心与钢球节圆中心不同心(如图1所示两中心距离为R12),故其接触载荷径向分量也会产生抵抗力矩M1ir,所有钢球的接触载荷径向分量抵抗力矩M1ir与接触载荷轴向分量抵抗力矩M1ia相加构成左列钢球总抵抗力矩M1。
M1ir=R12Q1icosα1ψicosψi=0.5de1Q1icosα1ψicosψi (13)
M1ia=0.5dm1Q1isinα1ψicosψi (14)
其中,dm1为左列钢球节圆直径(dm1=2R11),为表述便捷引入de1(de1=2R12),R12为轴承内圈中心O与钢球节圆中心O1’之间的距离。
右列圆锥滚子与滚道之间为线接触,对于右列使滚子载荷产生位移的方向与左列有所区别。在第i个滚子处,转角θi的分量为:
θi=θcosφi (16)
其中,φi为第i个圆锥滚子的位置角。
由于力矩作用而使右列圆锥滚子产生的径向和轴向位移分量如图2所示。在第i个滚子处产生的径向位移分量δrθi与轴向位移分量δaθi为:
δrθi=R22θcosφi (17)
δaθi=0.5dm2θcosφi (18)
上式中R22为轴承内圈中心O与钢球节圆中心O2’之间的距离,dm2为圆锥滚子节圆直径,则在第i个滚子处的径向总位移δ2ri为:
δ2ri=δrcosφi-δrθi=(δr-R22θ)cosφi (19)
在第i个滚子处的轴向总位移δ2ai为:
δ2ai=δ0-δaθi-δa=δ0-0.5dm2θcosφi-δa (20)
δ0为轴向游隙(即轴向预紧量);在第i个滚子处其外滚道法向接触总位移δ2i为径向总位移δ2ri和轴向总位移δ2ai在外滚道投影之和,其中αe为圆锥滚子与外圈接触角,αi为圆锥滚子与内圈接触角,如图3所示。
δ2i=δ2ricosαe+δ2aisinαe (21)
即:
δ2i=(δr-R22θ)cosφi cosαe+(δ0-0.5dm2θcosφi-δa)sinαe (22)
则第i个滚子与外圈接触载荷为:
Q2i=Kn2[max(δ2i,0)]10/9 (23)
对于圆锥滚子,其内、外接触角之差等于锥角,若其半锥角较小,则其内、外接触角之差也随之较小,则滚子与外滚道接触处的总刚度系数Kn2可按照圆柱滚子接触刚度公式计算:
Kn2=2.89×104×l0.82×Dw2 0.11 (24)
式中:l为滚子有效接触长度,Dw2为圆锥滚子平均直径。
圆锥滚子接触载荷径向分量总和为:
圆锥滚子接触载荷轴向分量总和为:
由于球-锥混合轴承为双列轴承,轴承内圈中心与滚子节圆中心不同心(如图1所示两中心距离为R22),故其接触载荷径向分量也会产生抵抗力矩M2ir,右列滚子的接触载荷径向分量总抵抗力矩M2ir与接触载荷轴向分量总抵抗力矩M2ia相加构成圆锥滚子总抵抗力矩M2,如下式所示:
M2ir=R22Q2i cosαecosφi=0.5de2Q2i cosαecosφi (27)
M2ia=0.5dm2Q2i sinαecosφi (28)
其中dm2为右列滚子节圆直径(dm2=2R21),为表述便捷引入de2(de2=2R22),R22为轴承内圈中心O与滚子节圆中心O2’之间的距离。
轴承内圈在外部载荷和所有滚动体反作用力的合力作用下处于平衡状态,平衡方程如下:
平衡方程(30)构成了以δr、δa、θ为未知量的三个非线性方程组,可运用Newton迭代法进行运算求解。将汽车不同行驶状态下驾驶室侧轮胎所承受的径向力Frn、轴向力Fan和颠覆力矩Mn分别作为外部负载(n=1、2、3、4),即可求取球、锥混合双列轮毂轴承单元内圈对应的径向位移量δr、轴向位移量δa和角位移量θ。
3.基于载荷谱下确定球、锥混合双列轮毂轴承单元最佳预紧量。
在一个试验周期内(周期为100s)轮毂轴承单元所受载荷要随着行驶情况(良好路面直行、左转、颠簸路面直行、右转)的不同而发生变化。假设在良好路面直行、左转、颠簸路面直行、右转四种行驶情况下,每种行驶情况下轮毂轴承单元所承受载荷的对应时间分配率分别为q1、q2、q3、q4,轮毂轴承转速分别为n1、n2、n3、n4(单位r/min),此时轴承的当量动载荷分别为Qe1、Qe2、Qe3、Qe4,则基于载荷谱多工况下,轴承套圈的平均当量动载荷为:
其中:(下标n=1、2、3、4),nn为良好路面直行、左转、颠簸路面直行、右转四种行驶情况下对应的轮毂轴承主轴转速(单位r/min),un为四种行驶情况下对应的汽车行驶速度(单位km/h)。
求得基于载荷谱下轴承套圈的平均当量动载荷Qe,根据轴承额定寿命计算方法,分别计算球、锥混合双列轮毂轴承单元左列钢球基本额定寿命L10_1和右列圆锥滚子基本额定寿命L10_2(L10_1和L10_2求取方法与L10求取方法类似),然后得到整个球、锥混合双列轮毂轴承单元的综合寿命L10为:
对于点接触e1=9/8,对于线接触e2=10/9,e取两者较小值,故e=e2。
由于轮毂轴承单元均采用过盈配合进行装配,所以其轴向预紧量δ0不大于零。以0作为轴向预紧量的初始值,-0.0025作为增量,对球、锥混合双列轮毂轴承单元的轴向预紧量取一系列初值:0、-0.0025、-0.005、-0.0075、...。
分别以每个轴向预紧量值作为轴向游隙,利用上文平衡方程(式30)和综合寿命公式(式32),求取每个轴向预紧量下相对应的轮毂轴承单元的单列套圈寿命和综合寿命,以一系列轴向预紧量初值作为横坐标,求取的相对应轮毂轴承单元套圈寿命作为纵坐标,可分别绘制出:左列钢球套圈寿命-预紧量关系图、右列圆锥滚子套圈寿命-预紧量关系图、轮毂轴承单元的综合寿命-预紧量关系图。从中可以确定出轮毂轴承单元单列套圈和整个轮毂轴承单元的最佳预紧量。
下面以某一具体的球、锥混合双列轮毂轴承为例对本发明进行说明。
步骤1:制定球、锥混合双列轮毂轴承单元载荷谱
载荷谱的制定便于更好地分析轮毂轴承的性能,为轮毂轴承的耐久性分析计算提供了更加科学合理的依据。由于汽车行驶环境的多样性和不确定性,对于球、锥混合双列轮毂轴承单元,其载荷谱包含四种行驶状态:良好路面直行、左转、颠簸路面直行、右转(见表1所示,)四种状态分别对应下标n=1、2、3、4。
借鉴国外轮毂轴承耐久性试验标准周期,取一个实验周期为100s,基于中国路面基本情况,取其直线行驶时间分配率为94%,转弯时间占用率为6%。其中良好路面直行时间占用率和颠簸路面直行时间占用率分别为50%和44%;左转和右转时间占有率都为3%。对某型号越野车驾驶室侧(即汽车前轴左侧)轮毂轴承进行分析,采用公式(1)、(2)和(3)计算其轮胎径向力Frn、轮胎轴向力Fan和倾覆力矩Mn。具体参数值见表1和表2所示。
表1
表2
根据汽车的相关参数,按照不同行驶路况下对应的载荷冲击系数fw和侧向加速度ag,计算出良好路面执行、左转、颠簸路面直行、右转四种行驶情况下,球、锥混合双列轮毂轴承单元所受的径向载荷Frn、轴向载荷Fan和倾覆力矩Mn(n=1、2、3、4),具体结果见表3.
表3
以时间为横坐标,不同行驶状态载荷值为纵坐标,取一个实验周期为100s,按照表1所示越野车四种行驶情况下所对应的时间分配率,即可获得球、锥混合双列轮毂轴承单元相对应的载荷谱,如图4所示。
步骤2:建立球锥双列轮毂轴承单元静力学模型并求取其内圈中心位移量。
假设球、锥混合双列轮毂轴承单元其外圈固定,当内圈承受径向力Fr、轴向力Fa和颠覆力矩M联合作用时,其径向位移为δr、轴向位移为δa、角位移为θ,本文定义对于每列轴承使滚动体载荷增加的位移为正位移,反之则为负位移。球、锥混合双列轮毂轴承单元结构及受力示意图如图1所示。并规定左列钢球相关表达式所用下标为1,右列圆锥滚子相关表达式所用下标为2,用i表示滚动体(左列钢球、右列圆锥滚子)。
以某企业设计出的新型球、锥混合双列轮毂轴承单元为例,其轴承相关参数如下:左列钢球直径、钢球数和钢球节圆直径分别为:DW1=12.6mm、Z1=16、dm1=65mm、钢球原始接触角α0=45°,内、外沟曲率系数分别是fi=0.513、fe=0.523,轴承内圈中心O与钢球节圆中心O1’之间的距离R12=10.5mm;右列圆锥滚子平均直径、滚子数、滚子母线接触长度和滚子节圆直径分别为:DW2=9.46mm、Z2=21、L=12.3mm、dm2=63mm、圆锥滚子与外圈接触角αe=25.46°、圆锥滚子与内圈接触角αi=22.16°,轴承内圈中心O与滚子节圆中心O2’之间的距离R22=10.5mm。
以步骤1中计算出的:良好路面直行、左转、颠簸路面直行、右转等四种行驶状态下,驾驶室侧轮胎所承受的径向力Frn、轴向力Fan和颠覆力矩Mn分别作为外部负载(n=1、2、3、4),轴向预紧量δ0设置为0,运用Newton迭代法进行运算求解,即可求取球、锥混合双列轮毂轴承单元内圈对应的径向位移量δr、轴向位移量δa和角位移量θ,如表4所示。
表4
步骤3:基于载荷谱下确定球、锥混合双列轮毂轴承单元最佳预紧量。
在一个试验周期内(周期为100s)轮毂轴承单元所受载荷随着行驶情况(良好路面直行、左转、颠簸路面直行、右转)的不同而发生变化。假设在良好路面直行、左转、颠簸路面执行、右转四种行驶情况下,每种行驶情况下轮毂轴承单元所承受载荷的对应时间分配率分别为q1、q2、q3、q4;轮毂轴承转速分别为n1、n2、n3、n4(单位r/min),其中:(下标n=1、2、3、4),un为四种行驶情况下对应的汽车行驶速度(单位km/h),数值见表1所示。
此时轴承的当量动载荷分别为Qe1、Qe2、Qe3、Qe4,则基于载荷谱多工况下,由公式(31)可计算出轴承套圈的平均当量动载荷求得基于载荷谱下轴承套圈的平均当量动载荷根据轴承额定寿命计算方法,可分别计算出球、锥混合双列轮毂轴承单元左列钢球基本额定寿命L10_1、右列圆锥滚子基本额定寿命L10_2和整个球、锥混合双列轮毂轴承单元的综合寿命L10(见式32)。
由于轮毂轴承单元均采用过盈配合进行装配,所以其轴向预紧量δ0不大于零。以0作为轴向预紧量的初始值,-0.0025作为增量,对球、锥混合双列轮毂轴承单元的轴向预紧量取一系列初值:0、-0.0025、-0.005、-0.0075、...。
分别以每个轴向预紧量值作为轴向游隙,利用上文平衡方程(式30)和综合寿命公式(式32),求取每个轴向预紧量下相对应的轮毂轴承单元的单列套圈寿命和综合寿命,从中可以精确确定出轮毂轴承单元单列套圈和整个轮毂轴承单元的最佳预紧量。
以一系列轴向预紧量初值作为横坐标,求取的相对应轮毂轴承单元套圈寿命作为纵坐标,可分别绘制出:左列钢球套圈寿命-预紧量关系图(如图5所示)、右列圆锥滚子套圈寿命-预紧量关系图(如图6所示)、轮毂轴承单元的综合寿命-预紧量关系图(如图7所示)。从套圈寿命-预紧量关系图中可以更方便快捷的看出轮毂轴承单元单列套圈和整个轮毂轴承单元的最佳预紧量。从图5中可以确定左列钢球套圈最佳预紧量δ0为-0.025;从图6中可以确定右列圆锥滚子预紧量δ0在0到-0.0125之间都比较好;从图7中可以确定整个球、锥混合轮毂轴承单元的最佳预紧量δ0为-0.0175。
本发明获取球、锥混合双列轮毂轴承载荷分布的装置的实施例
本实施例中获取球、锥混合双列轮毂轴承载荷分布的装置包括载荷确定模块、位移量确定模块、寿命计算模块和最佳预紧量确定模块,载荷确定模块用于根据轮胎受力公式计算不同行驶情况下球、锥混合双列轮毂轴承单元所受的径向载荷Frn、轴向载荷Fan和倾覆力矩Mn,并根据各行驶情况所对应的时间分配率确定球、锥混合双列轮毂轴承单元相对应的载荷谱;位移量确定模块用于根据球、锥混合双列轮毂轴承的力学模型建立以其内圈中心径向位移δr、轴向位移δa和角位移θ三个未知量构成的内圈力学平衡方程,并以所确定的径向载荷Frn、轴向载荷Fan和颠覆力矩Mn作为外部负载求解该平衡方程,以得到球、锥混合双列轮毂轴承单元内圈对应的径向位移量δr、轴向位移量δa和角位移量θ;寿命计算模块用于根据于所建立的载荷谱计算平均当量动载荷,根据轴承额定寿命计算球、锥混合双列轮毂轴承单元钢球基本额定寿命L10_1和或圆锥滚子基本额定寿命L10_2;最佳预紧量确定模块用于以每个轴向预紧量值作为轴向游隙,求取每个轴向预紧量下相对应的轮毂轴承单元的单列套圈寿命和综合寿命,确定出轮毂轴承单元单列套圈和整个轮毂轴承单元的最佳预紧量。各模块的具体实现手段已在方法的实施例中进行了说明,这里不再赘述。
Claims (10)
1.一种球、锥混合双列轮毂轴承最佳预紧量的确定方法,其特征在于,该方法包括以下步骤:
1)根据轮胎受力公式计算不同行驶情况下球、锥混合双列轮毂轴承单元所受的径向载荷Frn、轴向载荷Fan和倾覆力矩Mn,并根据各行驶情况所对应的时间分配率确定球、锥混合双列轮毂轴承单元相对应的载荷谱;
2)根据球、锥混合双列轮毂轴承的力学模型建立以其内圈中心径向位移δr、轴向位移δa和角位移θ三个未知量构成的内圈力学平衡方程,并以步骤1)中所确定的径向载荷Frn、轴向载荷Fan和颠覆力矩Mn作为外部负载求解该平衡方程,以得到球、锥混合双列轮毂轴承单元内圈对应的径向位移量δr、轴向位移量δa和角位移量θ;
3)基于所建立的载荷谱计算平均当量动载荷,根据轴承额定寿命计算球、锥混合双列轮毂轴承单元钢球基本额定寿命L10_1和或圆锥滚子基本额定寿命L10_2;
4)以每个轴向预紧量值作为轴向游隙,建立预紧量与轮毂轴承单元寿命关系,以轮毂轴承单元单列套圈和或整个轮毂轴承单元寿命最长所对应的预紧量为最佳预紧量。
2.根据权利要求1所述的球、锥混合双列轮毂轴承最佳预紧量的确定方法,其特征在于,所述步骤4)中平均当量动载荷的计算公式为:
其中Qei为第i种行驶情况下轴承的当量动载荷,ni为第i种行驶情况下轮毂轴承转速,qi为第i种行驶情况轮毂轴承单元所承受载荷的对应时间分配率,下标n为行驶情况的种类。
3.根据权利要求1所述的球、锥混合双列轮毂轴承最佳预紧量的确定方法,其特征在于,所述步骤4)中的综合寿命L10为:
其中e1为点接触,e2为线接触,e为e1和e2中的较小者。
4.根据权利要求1所述的球、锥混合双列轮毂轴承最佳预紧量的确定方法,其特征在于,步骤2)所建立的平衡方程为:
Q1i为滚道对钢球的法向载荷,α1ψi表示受载后任意位置钢球的接触角,ψi表示钢球的位置角,Z1表示钢球数,dm1为钢球节圆直径,de1为轴承内圈中心O与钢球节圆中心O1’之间距离的2倍,Q2i为滚子与外圈接触载荷,αe为圆锥滚子与外圈接触角,φi为第i个圆锥滚子的位置角,Z2表示滚子数,dm2为滚子节圆直径,de2为轴承内圈中心O与滚子节圆中心O2’之间距离的2倍。
5.根据权利要求4所述的球、锥混合双列轮毂轴承最佳预紧量的确定方法,其特征在于,所述滚道对钢球的法向载荷和滚子与外圈接触载荷的计算公式分别为:
Q1i=Kn1[max(δ1i,0)]3/2
Q2i=Kn2[max(δ2i,0)]10/9
δ2i=(δr-R22θ)cosφicosαe+(δ0-0.5dm2θcosφi-δa)sinαe
其中Kn1为滚动体与内外圈之间的总负载-变形常数,α0为原始接触角,Ri为钢球内滚道沟曲率中心所在圆的半径,A为内外圈原始沟曲率中心距,R22为轴承轴承内圈中心O与钢球节圆中心O2’之间的距离,δ0为轴向游隙即轴向预紧量。
6.根据权利要求1所述的球、锥混合双列轮毂轴承最佳预紧量的确定方法,其特征在于,所述步骤1)中径向载荷Frn、轴向载荷Fan和倾覆力矩Mn的计算公式分别为:
Mn=Frnso+FanRs
其中fw为载荷冲击系数,W为汽车重量,g为重力加速度,ag为汽车侧向加速度率,h为汽车质心高度,l为前轴轮间距,s0为径向受力偏心距,Rs为汽车轮胎半径,tr为汽车前轴轮间距。
7.根据权利要求1-6中任一项所述的球、锥混合双列轮毂轴承最佳预紧量的确定方法,其特征在于,步骤2)中的平衡方程采用Newton-Raphson迭代法进行求解。
8.一种球、锥混合双列轮毂轴承最佳预紧量的确定装置,其特征在于,该装置包括载荷确定模块、位移量确定模块、寿命计算模块和最佳预紧量确定模块,
所述载荷确定模块用于根据轮胎受力公式计算不同行驶情况下球、锥混合双列轮毂轴承单元所受的径向载荷Frn、轴向载荷Fan和倾覆力矩Mn,并根据各行驶情况所对应的时间分配率确定球、锥混合双列轮毂轴承单元相对应的载荷谱;
所述位移量确定模块用于根据球、锥混合双列轮毂轴承的力学模型建立以其内圈中心径向位移δr、轴向位移δa和角位移θ三个未知量构成的内圈力学平衡方程,并以所确定的径向载荷Frn、轴向载荷Fan和颠覆力矩Mn作为外部负载求解该平衡方程,以得到球、锥混合双列轮毂轴承单元内圈对应的径向位移量δr、轴向位移量δa和角位移量θ;
所述的寿命计算模块用于根据于所建立的载荷谱计算平均当量动载荷,根据轴承额定寿命计算球、锥混合双列轮毂轴承单元钢球基本额定寿命L10_1和或圆锥滚子基本额定寿命L10_2;
所述最佳预紧量确定模块用于以每个轴向预紧量值作为轴向游隙,建立预紧量与轮毂轴承单元寿命关系,以轮毂轴承单元单列套圈和或整个轮毂轴承单元寿命最长所对应的预紧量为最佳预紧量。
9.根据权利要求8所述的球、锥混合双列轮毂轴承最佳预紧量的确定装置,其特征在于,所述寿命计算模块中采用的平均当量动载荷的计算公式为:
其中Qei为第i种行驶情况下轴承的当量动载荷,ni为第i种行驶情况下轮毂轴承转速,qi为第i种行驶情况轮毂轴承单元所承受载荷的对应时间分配率,下标n为行驶情况的种类。
10.根据权利要求8所述的球、锥混合双列轮毂轴承最佳预紧量的确定装置,其特征在于,所述位移量确定模块中建立的平衡方程为:
Q1i为滚道对钢球的法向载荷,α1ψi表示受载后任意位置钢球的接触角,ψi表示钢球的位置角,Z1表示钢球数,dm1为钢球节圆直径,de1为轴承内圈中心O与钢球节圆中心O1’之间距离的2倍,Q2i为滚子与外圈接触载荷,αe为圆锥滚子与外圈接触角,φi为第i个圆锥滚子的位置角,Z2表示滚子数,dm2为滚子节圆直径,de2为轴承内圈中心O与滚子节圆中心O2’之间距离的2倍。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610736007.9A CN106326570B (zh) | 2016-08-27 | 2016-08-27 | 球、锥混合双列轮毂轴承最佳预紧量的确定方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610736007.9A CN106326570B (zh) | 2016-08-27 | 2016-08-27 | 球、锥混合双列轮毂轴承最佳预紧量的确定方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106326570A true CN106326570A (zh) | 2017-01-11 |
CN106326570B CN106326570B (zh) | 2019-09-13 |
Family
ID=57790668
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610736007.9A Active CN106326570B (zh) | 2016-08-27 | 2016-08-27 | 球、锥混合双列轮毂轴承最佳预紧量的确定方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106326570B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107577904A (zh) * | 2017-11-01 | 2018-01-12 | 杭州电子科技大学 | 一种等寿命轮毂轴承的设计方法 |
CN108692938A (zh) * | 2017-04-06 | 2018-10-23 | 湖南南方宇航高精传动有限公司 | 一种获取滚动轴承寿命的方法 |
CN110502816A (zh) * | 2019-08-13 | 2019-11-26 | 上海应用技术大学 | 汽车轮毂轴承的载荷谱编制方法、寿命分析方法及装置 |
CN111767617A (zh) * | 2020-06-28 | 2020-10-13 | 上海羿弓精密科技有限公司 | Rv减速机主轴承预紧量选配方法及选配系统和存储介质 |
CN112033255A (zh) * | 2020-08-14 | 2020-12-04 | 常州机电职业技术学院 | 一种双列短圆柱滚子轴承预紧量的计算方法 |
CN112560168A (zh) * | 2020-11-17 | 2021-03-26 | 湖北文理学院 | 轮毂轴承单元过盈量确定方法、装置、设备及存储介质 |
CN112747884A (zh) * | 2021-01-29 | 2021-05-04 | 重庆长安汽车股份有限公司 | 一种轮毂轴承抗冲击能力评价方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3759588A (en) * | 1971-11-12 | 1973-09-18 | Nasa | High speed hybrid bearing comprising a fluid bearing & a rolling bearing connected in series |
CN103090835A (zh) * | 2013-01-31 | 2013-05-08 | 烟台光洋精密轴承有限公司 | 一种模拟式汽车轮毂轴承轴向游隙测量装置 |
CN105179471A (zh) * | 2014-05-26 | 2015-12-23 | 斯凯孚公司 | 新型圆锥滚子轴承以及应用该轴承的新型轴系统 |
CN106089993A (zh) * | 2016-08-27 | 2016-11-09 | 河南科技大学 | 一种获取球锥混合双列轮毂轴承载荷分布的方法及装置 |
-
2016
- 2016-08-27 CN CN201610736007.9A patent/CN106326570B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3759588A (en) * | 1971-11-12 | 1973-09-18 | Nasa | High speed hybrid bearing comprising a fluid bearing & a rolling bearing connected in series |
CN103090835A (zh) * | 2013-01-31 | 2013-05-08 | 烟台光洋精密轴承有限公司 | 一种模拟式汽车轮毂轴承轴向游隙测量装置 |
CN105179471A (zh) * | 2014-05-26 | 2015-12-23 | 斯凯孚公司 | 新型圆锥滚子轴承以及应用该轴承的新型轴系统 |
CN106089993A (zh) * | 2016-08-27 | 2016-11-09 | 河南科技大学 | 一种获取球锥混合双列轮毂轴承载荷分布的方法及装置 |
Non-Patent Citations (3)
Title |
---|
牛荣军 等: ""双列球型轮毂轴承单元耐久性寿命计算分析"", 《河南科技大学学报:自然科学版》 * |
董晓 等: ""轿车轮毂轴承疲劳寿命的计算与分析"", 《河南科技大学学报:自然科学版》 * |
董晓: ""汽车轮毂双列角接触球轴承动力学仿真分析"", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108692938A (zh) * | 2017-04-06 | 2018-10-23 | 湖南南方宇航高精传动有限公司 | 一种获取滚动轴承寿命的方法 |
CN107577904A (zh) * | 2017-11-01 | 2018-01-12 | 杭州电子科技大学 | 一种等寿命轮毂轴承的设计方法 |
CN107577904B (zh) * | 2017-11-01 | 2021-01-05 | 杭州电子科技大学 | 一种等寿命轮毂轴承的设计方法 |
CN110502816A (zh) * | 2019-08-13 | 2019-11-26 | 上海应用技术大学 | 汽车轮毂轴承的载荷谱编制方法、寿命分析方法及装置 |
CN111767617A (zh) * | 2020-06-28 | 2020-10-13 | 上海羿弓精密科技有限公司 | Rv减速机主轴承预紧量选配方法及选配系统和存储介质 |
CN111767617B (zh) * | 2020-06-28 | 2022-10-25 | 上海羿弓精密科技有限公司 | Rv减速机主轴承预紧量选配方法及选配系统和存储介质 |
CN112033255A (zh) * | 2020-08-14 | 2020-12-04 | 常州机电职业技术学院 | 一种双列短圆柱滚子轴承预紧量的计算方法 |
CN112560168A (zh) * | 2020-11-17 | 2021-03-26 | 湖北文理学院 | 轮毂轴承单元过盈量确定方法、装置、设备及存储介质 |
CN112560168B (zh) * | 2020-11-17 | 2022-06-24 | 湖北文理学院 | 轮毂轴承单元过盈量确定方法、装置、设备及存储介质 |
CN112747884A (zh) * | 2021-01-29 | 2021-05-04 | 重庆长安汽车股份有限公司 | 一种轮毂轴承抗冲击能力评价方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106326570B (zh) | 2019-09-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106326570A (zh) | 球、锥混合双列轮毂轴承最佳预紧量的确定方法及装置 | |
CN106089993B (zh) | 一种获取球锥混合双列轮毂轴承载荷分布的方法及装置 | |
US20220364953A1 (en) | Compilation Method for Reliability Test Load Spectrum of High-Speed Bearing of Electric Drive System | |
CN110617953B (zh) | 新能源汽车高速电机轴承-转子系统动态特性分析方法 | |
CN104354700B (zh) | 一种基于无迹卡尔曼滤波的车辆参数在线估计方法 | |
CN103703263B (zh) | 制造车轮轴承装置的方法 | |
CN102564763A (zh) | 一种机床主轴轴承动态预紧力测量方法 | |
CA2641534C (en) | Multi-row symmetrical rolling bearing | |
CN107153734A (zh) | 一种定位预紧下角接触球轴承非线性刚度的计算系统及方法 | |
CN111814259A (zh) | 一种基于轮毂电机行星齿轮减速器及设计方法 | |
CN106768543A (zh) | 一种轮胎滚动阻力的测试装置及其方法 | |
CN101233000A (zh) | 车轮轴承单元 | |
CN115406669A (zh) | 一种多轴特种车辆侧翻指标优化方法 | |
Szabó | Vehicle test based validation of a tire brush model using an optical velocity sensor | |
CN202012566U (zh) | 汽车轮毂轴承单元 | |
CN206563656U (zh) | 一种试验室用轮毂双向动载荷的加载装置 | |
Huang et al. | Research on loaded brake performance test of trucks | |
CN106066915A (zh) | 一种获取球‑柱联合转盘轴承刚度矩阵的方法及装置 | |
JP2008106904A (ja) | 車輪用軸受装置の製造方法 | |
Lu et al. | Design Optimization of an Integrated E‐Type Multilink Suspension Wheel‐Side Drive System and Improvement of Vehicle Ride Comfort | |
Joshi et al. | Analysis and topological optimization of motorcycle front wheel | |
CN207033989U (zh) | 轴承用承载拉力的平衡组合式双半变径中圈 | |
CN113946919B (zh) | 一种含组合角不对中的深沟球轴承准静态解析模型分析方法 | |
Ismailov | DEVELOPMENT OF METHODOLOGY FOR CALCULATING TRANSMISSION OF WHEEL VEHICLES WITH 4× 2 WHEEL POSITION | |
Chaudhari et al. | Design and Development of Open Differential for Transmission System of Quad Bike |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20170111 Assignee: Zhejiang Zhaofeng Mechanical and Electronic Co., Ltd. Assignor: Henan University of Science and Technology Contract record no.: X2019980000338 Denomination of invention: Method and device for determining optimal pretension quantity of sphere and cone mixed double-line hub bearing Granted publication date: 20190913 License type: Exclusive License Record date: 20191111 |
|
EE01 | Entry into force of recordation of patent licensing contract |