发明名称 一种基于轴向载荷的提升机轴承疲劳寿命模型
摘要 本发明公开了一种基于轴向载荷的提升机轴承疲劳寿命模型,包括如下步骤:S1、对滚动体在轴承纵向截面内倾斜偏移量进行分板;S2、分析滚动体在仅受径向载荷滚动体状态下的受力;S3、分析滚动体在同时承受径向载荷和轴向载荷下的受力;S4、建立基于滚动体的提升机轴承疲劳寿命模型,在引入修正系数η,则基于最大剪应力的轴承寿命公式为:<img file="DDA0000960095960000011.GIF" wi="329" he="98" />本发明的基于轴向载荷的提升机轴承疲劳寿命模型,在考虑轴向载荷作用的情况下,提出新的提升机轴承疲劳寿命模型,为矿井提升机轴承寿命预测分析的可靠性提供技术支持。本发明在考虑轴向载荷的情况下,对提升机轴承系统的疲劳寿命模型并对原有模型进行修正,以适用于各种复杂工况下的轴承寿命分析。
申请公布号 CN105653890A 申请公布日期 2016.06.08
申请号 CN201610221199.X 申请日期 2016.04.07
申请人 东北大学 发明人 王新刚;王宝艳;陈小辉
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 沈阳优普达知识产权代理事务所(特殊普通合伙) 21234 代理人 俞鲁江
主权项 一种基于轴向载荷的提升机轴承疲劳寿命模型,其特征在于,包括如下步骤:S1、对滚动体在轴承纵向截面内倾斜偏移量进行分板;A1、对通过轴承轴线的纵向切面上的轴承受力情况进行分析,此时发生在该平面上的滚动体的倾斜会导致滚动体在套圈滚道上的不利压力分布,并使滚动体与轴承内外圈所接触区域的承受压力增大,而保持架也会由于滚动体的倾斜发生轻微形变,使得最后滚动体位置接近于轴承径向载荷分布图中的对应位置;A2、应用Palmgren公式,引出满足导出的等式的条件,确定每个滚动体上的轴向和径向载荷的分布以及滚动体与轴承环的相对位置;A3、计算滚动体末端和套圈滚道凸缘处的接触力Q<sub>f</sub>,将接触区域分为若干个宽度微元y,在划分的宽度微元中,近似地可认为宽度微元是处于刚性状态的,在原有的变形后其不再发生位移或塑形变形;此时,在滚动体和套圈滚道的接触区域内宽度微元与套圈所呈的角度为χ角,在这部分接触区域中可以取决于滚动体在套圈滚道内的嵌入深度的横截面的形式来表达,该凸缘部分的宽度间距为:<maths num="0001"><math><![CDATA[<mrow><msub><mi>&theta;</mi><mi>f</mi></msub><mo>=</mo><mo>-</mo><mfrac><mn>2</mn><msup><mi>E</mi><mo>&prime;</mo></msup></mfrac><msubsup><mo>&Integral;</mo><mn>0</mn><msub><mi>b</mi><mi>f</mi></msub></msubsup><msub><mi>p</mi><mi>f</mi></msub><mrow><mo>(</mo><msub><mi>z</mi><mi>f</mi></msub><mo>)</mo></mrow><mi>ln</mi><mo>|</mo><mrow><msub><mi>z</mi><mi>f</mi></msub><mo>-</mo><msub><mi>s</mi><mi>f</mi></msub></mrow><mo>|</mo><msub><mi>ds</mi><mi>f</mi></msub><mo>+</mo><mi>C</mi></mrow>]]></math><img file="FDA0000960095930000011.GIF" wi="867" he="146" /></maths>式中,C为常量;A4、Q<sub>f</sub>可以通过对宽度微元y进行积分得到,同时滚动体与套圈接触产生的接触力位置由到滚动体两侧的圆心距离f表示;S2、仅受径向载荷滚动体分析;在轴承整体发生倾斜的情况,轴承在受到径向力F<sub>r</sub>,轴向力F<sub>a</sub>和力矩M<sub>y</sub>后出现了轴承体相对于静止状态下的偏转角ξ,此时同时受到轴向力和径向力作用的滚动体在轴承整体横截面部分的差角为ψ<sub>lim</sub>,而差角ψ<sub>ε</sub>所代表的是在这角度范围内的套圈滚道间距恰好等于在这其内运动的圆柱滚动体的直径;此范围内的滚动体只受到径向力F<sub>r</sub>的作用,轴向力无影响;横截面圆周部分除上述两个角度范围外,其余差角部分内的滚动体可以近似认为恰好处于和套圈滚道接触的临界点,不受任何载荷作用;该差角ψ<sub>ε</sub>可以表示为:<img file="FDA0000960095930000021.GIF" wi="830" he="97" />式中:g值为轴承的径向游隙;Δδ表示为轴承圆周在圆心处的变形偏移量;因为这个差角所对应的滚动体只是承受径向力,在径向载荷的作用下,滚动体与滚道的接触所产生的接触变形由公式可以表达为:δ=0.39(θ<sub>a</sub>+θ<sub>b</sub>)<sup>0.9</sup>P/l<sub>e</sub><sup>0.8</sup><maths num="0002"><math><![CDATA[<mrow><mo>(</mo><msub><mi>&theta;</mi><mi>a</mi></msub><mo>+</mo><msub><mi>&theta;</mi><mi>b</mi></msub><mo>)</mo><mo>=</mo><mn>4</mn><mo>&lsqb;</mo><mfrac><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><msub><mi>v</mi><mi>a</mi></msub><mn>2</mn></msup><mo>)</mo></mrow><msub><mi>E</mi><mi>a</mi></msub></mfrac><mo>+</mo><mfrac><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><msub><mi>v</mi><mi>b</mi></msub><mn>2</mn></msup><mo>)</mo></mrow><msub><mi>E</mi><mi>b</mi></msub></mfrac><mo>&rsqb;</mo></mrow>]]></math><img file="FDA0000960095930000022.GIF" wi="802" he="147" /></maths>式中:P值为径向负荷;l<sub>e</sub>为滚动体的有效长度;θ<sub>a</sub>和θ<sub>b</sub>由杨氏模量与泊松比进行表示;a滚动体的材料系数;b套圈滚道的材料系数;这里引入一个积分因子:<maths num="0003"><math><![CDATA[<mrow><mi>J</mi><mrow><mo>(</mo><mi>&epsiv;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mi>&pi;</mi></mrow></mfrac><msubsup><mrow><mo>&Integral;</mo><mo></mo></mrow><mn>0</mn><mrow><mo>+</mo><mi>&psi;&epsiv;</mi></mrow></msubsup><msup><mrow><mo>[</mo><mn>1</mn><mo>-</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mi>&epsiv;</mi></mrow></mfrac><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>cos</mi><mo>)</mo></mrow><mo>]</mo></mrow><mi>t</mi></msup><mi>cos</mi><mi>&psi;d&psi;</mi></mrow>]]></math><img file="FDA0000960095930000023.GIF" wi="956" he="147" /></maths>式中:参数t为常数,对于滚子轴承,t=1.1;J(ε)的积分值和该接触变形量δ有关,而δ值又需要由J(ε)确定,则根据差角ψ<sub>ε</sub>即可以得到对应的接触变形量,即径向载荷P为:<maths num="0004"><math><![CDATA[<mrow><mi>P</mi><mo>=</mo><mfrac><mrow><mn>07117</mn><msup><msub><mi>l</mi><mi>e</mi></msub><mrow><mn>8</mn><mo>/</mo><mn>9</mn></mrow></msup><msup><mi>&delta;</mi><mrow><mn>10</mn><mo>/</mo><mn>9</mn></mrow></msup></mrow><mrow><mo>[</mo><mfrac><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><msub><mi>v</mi><mi>a</mi></msub><mn>2</mn></msup><mo>)</mo></mrow><msub><mi>E</mi><mi>a</mi></msub></mfrac><mo>+</mo><mfrac><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><msub><mi>v</mi><mi>b</mi></msub><mn>2</mn></msup><mo>)</mo></mrow><msub><mi>E</mi><mi>b</mi></msub></mfrac><mo>]</mo></mrow></mfrac></mrow>]]></math><img file="FDA0000960095930000024.GIF" wi="496" he="187" /></maths>则最大受载的滚动体负荷为<maths num="0005"><math><![CDATA[<mrow><mi>Q</mi><mo>=</mo><mfrac><mi>P</mi><mrow><mi>Z</mi><mi>J</mi><mrow><mo>(</mo><mi>&epsiv;</mi><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000960095930000025.GIF" wi="243" he="132" /></maths>式中:Z为轴承内的滚动体数;则在该滚动体分布圆周内距离负荷作用线为ψ值的圆柱滚动体的接触载荷为:<img file="FDA0000960095930000026.GIF" wi="581" he="96" />S3、同时承受径向载荷和轴向载荷的滚动体分析;由之前所给出的滚动体和套圈凸缘部分的宽度间距为θ<sub>f</sub>和轴承圆周圆心的变形偏移量Δδ+g/2,可以确定出一个参量:<maths num="0006"><math><![CDATA[<mrow><mi>k</mi><mo>=</mo><mfrac><mrow><msub><mi>&theta;</mi><mi>f</mi></msub><msub><mi>l</mi><mi>e</mi></msub></mrow><mrow><mrow><mo>(</mo><msub><mi>&Delta;</mi><mi>&delta;</mi></msub><mo>+</mo><mfrac><mi>g</mi><mn>2</mn></mfrac><mo>)</mo></mrow><msub><mi>D</mi><mi>r</mi></msub></mrow></mfrac></mrow>]]></math><img file="FDA0000960095930000031.GIF" wi="355" he="155" /></maths>式中:l<sub>e</sub>为滚动体与套圈滚道的有效接触长度;D<sub>r</sub>为圆柱滚动体的底面直径;g值为轴承的径向间隙;根据参数k得到一个比参量m为:<maths num="0007"><math><![CDATA[<mrow><mi>m</mi><mo>=</mo><mfrac><mrow><mn>2</mn><mi>k</mi></mrow><mrow><mi>k</mi><mo>+</mo><mi>cos</mi><mi>&psi;</mi></mrow></mfrac></mrow>]]></math><img file="FDA0000960095930000032.GIF" wi="323" he="131" /></maths>式中:ψ为同时受到径向载荷和轴向载荷的滚动体在轴承受力圆周中所对应的差角;根据m值的大小范围可以确定出此时条件下一个滚动体所受径向载荷量为:<maths num="0008"><math><![CDATA[<mrow><mi>Q</mi><mrow><mo>(</mo><mi>&psi;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mn>9</mn><msup><mi>C&delta;</mi><mrow><mn>10</mn><mo>/</mo><mn>9</mn></mrow></msup></mrow><mrow><mn>19</mn><mi>m</mi></mrow></mfrac><mo>&lsqb;</mo><mn>1</mn><mo>-</mo><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>m</mi><mo>)</mo></mrow><mrow><mn>19</mn><mo>/</mo><mn>9</mn></mrow></msup><mo>&rsqb;</mo></mrow>]]></math><img file="FDA0000960095930000033.GIF" wi="730" he="131" /></maths>其中根据几何关系为:<maths num="0009"><math><![CDATA[<mrow><mrow><mo>(</mo><msub><mi>&Delta;</mi><mi>&delta;</mi></msub><mo>+</mo><mfrac><mi>g</mi><mn>2</mn></mfrac><mo>)</mo></mrow><mi>cos</mi><mi>&psi;</mi><mo>=</mo><msub><mi>&delta;</mi><mi>max</mi></msub><mo>-</mo><mn>2</mn><mi>&chi;</mi><mfrac><msub><mi>l</mi><mi>e</mi></msub><mn>2</mn></mfrac></mrow>]]></math><img file="FDA0000960095930000034.GIF" wi="675" he="126" /></maths><maths num="0010"><math><![CDATA[<mrow><mi>tan</mi><mi>&chi;</mi><mo>=</mo><mi>&chi;</mi><mo>=</mo><mfrac><msub><mi>&theta;</mi><mi>f</mi></msub><msub><mi>D</mi><mi>r</mi></msub></mfrac></mrow>]]></math><img file="FDA0000960095930000035.GIF" wi="355" he="131" /></maths>即可表示为:<maths num="0011"><math><![CDATA[<mrow><msub><mi>&delta;</mi><mi>max</mi></msub><mo>=</mo><mrow><mo>(</mo><mo>(</mo><mrow><msub><mi>&Delta;</mi><mi>&delta;</mi></msub><mo>+</mo><mfrac><mi>g</mi><mn>2</mn></mfrac></mrow><mo>)</mo><mi>cos</mi><mi>&psi;</mi><mo>+</mo><mfrac><mi>k</mi><mrow><mo>(</mo><msub><mi>&Delta;</mi><mi>&delta;</mi></msub><mo>+</mo><mfrac><mi>g</mi><mn>2</mn></mfrac><mo>)</mo></mrow></mfrac><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000960095930000036.GIF" wi="783" he="163" /></maths>此时也可根据这个m值求出的径向载荷应用点为:<maths num="0012"><math><![CDATA[<mrow><mi>e</mi><mrow><mo>(</mo><mi>&psi;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><msub><mi>l</mi><mi>e</mi></msub><mn>2</mn></mfrac><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mn>2</mn><mi>m</mi></mfrac><mo>{</mo><mfrac><mrow><mn>19</mn><mo>[</mo><mn>1</mn><mo>-</mo><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>m</mi><mo>)</mo></mrow><mfrac><mn>28</mn><mn>9</mn></mfrac></msup></mrow><mrow><mo>[</mo><mn>1</mn><mo>-</mo><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>m</mi><mo>)</mo></mrow><mfrac><mn>19</mn><mn>9</mn></mfrac></msup></mrow></mfrac><mo>-</mo><mn>1</mn><mo>}</mo><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000960095930000037.GIF" wi="931" he="211" /></maths>滚动体受到的轴向载荷为:<maths num="0013"><math><![CDATA[<mrow><msub><mi>Q</mi><mi>a</mi></msub><mo>=</mo><mfrac><mi>Z</mi><mi>&pi;</mi></mfrac><msubsup><mo>&Integral;</mo><mn>0</mn><msub><mi>&psi;</mi><mi>lim</mi></msub></msubsup><mfrac><mrow><mn>2</mn><mi>e</mi><mrow><mo>(</mo><mi>&psi;</mi><mo>)</mo></mrow><mi>Q</mi><mrow><mo>(</mo><mi>&psi;</mi><mo>)</mo></mrow></mrow><msub><mi>&theta;</mi><mi>f</mi></msub></mfrac><mi>d</mi><mi>&psi;</mi></mrow>]]></math><img file="FDA0000960095930000041.GIF" wi="630" he="155" /></maths>在上述条件下产生的滚动体相对于轴承轴线的倾斜角为η,而其上下表面相对于原滚道位置的倾斜角为η<sub>ij</sub>和η<sub>oj</sub>(以内套圈滚道和外套圈滚道进行划分),与原套圈滚道相接触的滚动体,因接触形变产生的位移标定为Δij,Δoj;S4、建立基于滚动体的提升机轴承疲劳寿命模型;Hertz接触理论表明,套圈滚道表层承受的应力分布与滚道次表层中应力分布可以利用三向应力状态进行简化分析,该分析指出次表层的最大剪应力是由于接触表面上最大赫兹接触应力引起的;根据三向应力状态的分析可知,确定了由Hertz接触应力引起的次表层最大剪应力发生在z=0.78615的深度上,从而得到最大剪应力与表面Hertz接触应力的关系为:τ<sub>max</sub>=‑0.30028σ<sub>c</sub>滚子轴承的滚动体由于是圆柱体,其与内滚道接触时相当于两个圆柱体外接触,而滚动体与外滚道接触时相当于两个圆柱体内接触,结合曲率半径,即表面Hertz接触应力为:<maths num="0014"><math><![CDATA[<mrow><msub><mi>&sigma;</mi><mi>c</mi></msub><mo>=</mo><mn>0.418</mn><msqrt><mfrac><mi>PE</mi><msub><mi>l</mi><mi>e</mi></msub></mfrac><mrow><mo>(</mo><mfrac><mn>2</mn><mi>d</mi></mfrac><mo>+</mo><mfrac><mn>1</mn><mi>r</mi></mfrac><mo>)</mo></mrow></msqrt></mrow>]]></math><img file="FDA0000960095930000042.GIF" wi="516" he="186" /></maths>式中:P为受力最大的滚动体所承受的力;L为滚动体的有效工作长度;假设由于内外圈偏斜因素下产生的应力变化量为,由于轴承内部滚动体的受力并不是都相同,一部分数量的滚动体是不受到周向应力的影响,对该部分的滚动体不加以考虑,则在该工况下的次表面剪应力应为:<maths num="0015"><math><![CDATA[<mrow><msub><msup><mi>&tau;</mi><mo>&prime;</mo></msup><mi>max</mi></msub><mo>=</mo><mo>-</mo><mn>0.30028</mn><msub><mi>&sigma;</mi><mi>c</mi></msub><mo>+</mo><mfrac><mrow><mo>(</mo><msub><mi>&sigma;</mi><mi>r</mi></msub><mo>-</mo><msub><mi>&sigma;</mi><mi>&theta;</mi></msub><mo>)</mo></mrow><mn>2</mn></mfrac></mrow>]]></math><img file="FDA0000960095930000043.GIF" wi="730" he="129" /></maths>因此,通过该工况下对最大剪应力的影响后,计算出的轴承寿命为:<maths num="0016"><math><![CDATA[<mrow><msup><mi>L</mi><mo>&prime;</mo></msup><mo>=</mo><msup><mrow><mo>(</mo><mfrac><msub><mi>&tau;</mi><mi>max</mi></msub><mrow><msub><msup><mi>&tau;</mi><mo>&prime;</mo></msup><mi>max</mi></msub></mrow></mfrac><mo>)</mo></mrow><mn>9</mn></msup><mi>L</mi></mrow>]]></math><img file="FDA0000960095930000051.GIF" wi="334" he="119" /></maths>在引入修正系数η,则基于最大剪应力的轴承寿命公式为:L″=ηL′其中修正系数η≤1。
地址 110004 辽宁省沈阳市和平区文化路3号三巷11号