• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
    • / 15
    • 下载费用:30 金币  

    重庆时时彩108期开奖: 一种滚动轴承等效刚度与等效阻尼的计算方法.pdf

    关 键 词:
    一种 滚动轴承 等效 刚度 阻尼 计算方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201610831098.4

    申请日:

    2016.09.19

    公开号:

    CN106354974A

    公开日:

    2017.01.25

    当前法律状态:

    实审

    有效性:

    审中

    法律详情: 实质审查的生效IPC(主分类):G06F 17/50申请日:20160919|||公开
    IPC分类号: G06F17/50 主分类号: G06F17/50
    申请人: 湖南科技大学
    发明人: 伍济钢; 蒋勉; 石海波; 袁继广; 张双健; 王刚
    地址: 411201 湖南省湘潭市雨湖区桃源路
    优先权:
    专利代理机构: 湘潭市汇智专利事务所(普通合伙) 43108 代理人: 颜昌伟
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201610831098.4

    授权公告号:

    |||

    法律状态公告日:

    2017.03.01|||2017.01.25

    法律状态类型:

    实质审查的生效|||公开

    摘要

    本发明公开了一种滚动轴承等效刚度与等效阻尼的计算方法,包括如下步骤:1)基于弹性流体动力润滑理论,将滚动体、油膜、滚道三者之间的接触区域进行细化;2)通过压力分布双线性逼近函数以及差分法,求出细化后各区域内的接触的弹性变形刚度、油膜压力以及油膜厚度,利用线性扰动方程求解细化后各区域内的油膜动刚度与油膜阻尼;进而计算滚动体与滚道相接触的接触刚度和接触阻尼;3)在步骤2)计算得到的滚动体与滚道相接触的接触刚度和接触阻尼的基础上,计算出滚动轴承的等效刚度与等效阻尼。本发明具有计算结果准确度高的优点,可为轴承转子系统的设计提供了技术支持,降低了设计的轴承转子系统的故障率。

    权利要求书

    1.一种滚动轴承等效刚度与等效阻尼的计算方法,包括如下步骤:1)基于弹性流体动
    力润滑理论,将滚动体、油膜、滚道三者之间的接触区域进行细化;2)通过压力分布双线性
    逼近函数以及差分法,求出细化后各区域内的接触的弹性变形刚度、油膜压力以及油膜厚
    度,利用线性扰动方程求解细化后各区域内的油膜动刚度与油膜阻尼;进而计算滚动体与
    滚道相接触的接触刚度和接触阻尼;3)在步骤2)计算得到的滚动体与滚道相接触的接触刚
    度和接触阻尼的基础上,计算出滚动轴承的等效刚度与等效阻尼。
    2.根据权利要求1所述的滚动轴承等效刚度与等效阻尼的计算方法,步骤1)中滚动体、
    油膜、滚道三者之间的接触细化时,不考虑滚动体与内外滚道相对滑动,将滚动体、油膜、滚
    道三者之间的接触,根据接触弹性变形分为油膜入口区、弹性接触区和油膜出口区三个区。
    3.根据权利要求1所述的滚动轴承等效刚度与等效阻尼的计算方法,步骤2)中细化后
    各区域内的油膜动刚度与油膜阻尼求解方法为:联立弹性形变方程和Reynolds方程,并运
    用复合直接迭代法,求出油膜静态下的压力和厚度值,然后利用扰动方程建立扰动方程组
    并求解微动量和一阶微动量,求得线性近似油膜刚度值和油膜阻尼值。
    4.根据权利要求1所述的滚动轴承等效刚度与等效阻尼的计算方法,步骤2)中求解弹
    性变形时,在四节点矩形单元上使用双线性函数逼近压力分布求解弹性变形。
    5.根据权利要求3所述的滚动轴承等效刚度与等效阻尼的计算方法,在扰动方程组求
    解时,利用泰勒级数将油膜厚度与油膜压力展开,取坐标原点为静平衡位置。
    6.根据权利要求2所述的滚动轴承等效刚度与等效阻尼的计算方法,步骤2)中,忽略油
    膜出口区的油膜刚度和油膜阻尼,对油膜入口区的油膜刚度与油膜阻尼及弹性接触区滚动
    体与滚道的弹性接触刚度、油膜刚度以及油膜阻尼分别计算,然后在此基础上计算滚动体
    与滚道相接触的接触刚度和接触阻尼。
    7.根据权利要求6所述的滚动轴承等效刚度与等效阻尼的计算方法,计算油膜入口区
    油膜阻尼和油膜刚度时,假设以下条件:
    由于滚动轴承中滚动体与滚道接触面的形状是狭长的椭圆,因此滚动方向的短半轴远
    比长半轴小,所以接触面形状可以近似于同等椭圆长径与短径的矩形接触面,并忽略接触
    区域边缘上的漏油;
    滚动体与滚道之间的间隙为抛物线;
    轴承中油膜厚度大于滚动体、内外圈的表面粗糙度;
    惯性力的值小于粘性力;
    忽略重力影响;
    考虑到Reynolds方程适用性,认为粘度是恒定的;
    忽略接触过程中气穴现象产生的反压力。
    8.根据权利要求1所述的滚动轴承等效刚度与等效阻尼的计算方法,计算轴承径向等
    效刚度与等效阻尼时,考虑滚动轴承中滚动体与滚道之间存在径向游隙及在轴承径向载荷
    下轴承的内外圈套之间产生相对位移。

    说明书

    一种滚动轴承等效刚度与等效阻尼的计算方法

    技术领域

    本发明涉及一种滚动轴承的等效刚度与等效阻尼的计算方法。

    背景技术

    轴承转子系统设计需要考虑轴承的刚度和阻尼。在进行滚动轴承如临界转速等动
    力学特性分析时必须给出滚动轴承的动态特性参数:等效刚度与等效阻尼。等效刚度与等
    效阻尼的不精确直接给滚动轴承转子系统的动力学特性分析带来很大的误差,造成轴承转
    子系统设计时误差大,严重时会造成临界转速与额定转速接近,造成机毁人亡的大事故。

    现在,对滚动轴承进行动力学分析时,常将滚动轴承内部接触弹性刚度作为滚动
    轴承整体刚度而不考虑油膜的影响,同时忽略了阻尼对系统稳态下的实质影响,对于阻尼
    与刚度的获取主要依赖经验值,得到的等效刚度与等效阻尼的精确度低。方兵等对轴承特
    性进行分析时利用实验测量轴承的等效阻尼和等效刚度。何芝仙等对于轴承动力学分析时
    考虑了刚度而忽略阻尼的影响。Gupta的模型将阻尼的影响简化处理。Hagiu提出了一个动
    态理论分析模型,强调了高速滚动接触动态取决于机理Hertz接触的弹性刚度和触区入口
    处的润滑剂刚度与阻尼。刘秀海利用滚动体在流体中平移时所受的阻尼来模拟滚子所受的
    粘性阻尼,但忽略了阻尼挤压特性。Elsermans与Walford认为球轴承径向、轴向的刚度与阻
    尼实验测量结果均比预期值要大,这个结果可以通过外圈-滚动体-内圈-轴的分析来解释,
    但理论计算不够完善。Harsha与Kankar提出基于Hertz弹性变形的球轴承非线性模型,并引
    入了经验阻尼来分析其振动传递。陈斌等对油膜阻尼进行了理论计算,但缺乏有力的实验
    验证。

    发明内容

    为了解决上述技术问题,本发明提供一种滚动轴承等效刚度与等效阻尼计算结果
    精度高的滚动轴承等效刚度和等效阻尼计算方法,它为轴承转子系统的设计提供了技术支
    持,降低了设计的轴承转子系统的故障率。

    本发明解决上述技术问题的技术方案是:一种滚动轴承的等效刚度与等效阻尼的
    计算方法,包括如下步骤:1)基于弹性流体动力润滑理论,将滚动体、油膜、滚道三者之间的
    接触区域进行细化;2)通过压力分布双线性逼近函数以及差分法,求出细化后各区域内的
    接触的弹性变形刚度、油膜压力以及油膜厚度,利用线性扰动方程求解局部油膜动刚度与
    油膜阻尼;3)根据滚动轴承动态载荷分布,并结合步骤2)计算得到的弹性变形刚度、油膜刚
    度以及油膜阻尼,计算出滚动轴承的等效刚度与等效阻尼。

    上述的滚动轴承等效刚度与等效阻尼的计算方法中,步骤1)中滚动体、油膜、滚道
    三 者之间的接触细化时,不考虑滚动体与内外滚道相对滑动,将滚动体、油膜、滚道三者之
    间的接触,根据接触弹性变形分为油膜入口区、弹性接触区和油膜出口区三个区。

    上述的滚动轴承等效刚度与等效阻尼的计算方法中,步骤2)中细化后各区域内的
    油膜动刚度与油膜阻尼求解方法为:联立弹性形变方程和Reynolds方程,并运用复合直接
    迭代法,求出油膜静态下的压力和厚度值,然后利用扰动方程建立扰动方程组并求解微动
    量和一阶微动量,求得线性近似油膜刚度值和油膜阻尼值。

    上述的滚动轴承等效刚度与等效阻尼的计算方法中,步骤2)中求解弹性变形时,
    在四节点矩形单元上使用双线性函数逼近压力分布求解弹性变形。

    上述的滚动轴承等效刚度与等效阻尼的计算方法中,在扰动方程组求解时,利用
    泰勒级数将油膜厚度与油膜压力展开,取坐标原点为静平衡位置。

    上述的滚动轴承等效刚度与等效阻尼的计算方法中,步骤2)中,忽略油膜出口区
    的油膜刚度和油膜阻尼,对油膜入口区的油膜刚度与油膜阻尼及弹性接触区滚动体与滚道
    的弹性接触刚度、油膜刚度以及油膜阻尼分别计算,然后在此基础上计算滚动体与滚道相
    接触的接触刚度和接触阻尼。

    上述的滚动轴承等效刚度与等效阻尼的计算方法,计算油膜入口区油膜阻尼和油
    膜刚度时,假设以下条件:

    由于滚动轴承中滚动体与滚道接触面的形状是狭长的椭圆,因此滚动方向的短半
    轴远比长半轴小,所以接触面形状可以近似于同等椭圆长径与短径的矩形接触面,并忽略
    接触区域边缘上的漏油;

    滚动体与滚道之间的间隙为抛物线;

    轴承中油膜厚度大于滚动体、内外圈的表面粗糙度;

    惯性力的值小于粘性力;

    忽略重力影响;

    考虑到Reynolds方程适用性,认为粘度是恒定的;

    忽略接触过程中气穴现象产生的反压力。

    上述的滚动轴承等效刚度与等效阻尼的计算方法中,计算轴承径向等效刚度与等
    效阻尼时,考虑滚动轴承中滚动体与滚道之间存在径向游隙及在轴承径向载荷下轴承的内
    外圈套之间产生相对位移。

    与现有技术相比,本发明具有以下有益效果:

    (1)本发明细化了滚动体、油膜、滚道三者之间的接触,根据接触弹性变形分为三
    个区, 油膜入口区、弹性接触区、和油膜出口区,并考虑了滚动轴承中滚动体与滚道之间存
    在径向游隙,在轴承径向载荷下,轴承的内外圈套产生的相对位移,这样更加符合实际情
    况,计算出的结果更加准确;为轴承转子系统的设计提供了技术支持,降低了设计的轴承转
    子系统的故障率。

    (2)本发明在联立弹性形变方程和Reynolds方程的基础上,运用复合直接迭代法,
    求出油膜静态下的压力和厚度值,利用扰动方程求解方程组的微动量和一阶微动量,最终
    求得线性近似油膜刚度值和油膜阻尼值,具有计算速度快,计算精度高等优点。

    (3)本发明的计算不需要对滚动轴承进行各种测试,获得相关参数,而是直接可以
    应用,具有简单实用的优点。

    附图说明

    图1为本发明的简化的弹性流体润滑接触模型。

    图2为本发明的滚动轴承的油膜几何形状。

    图3为本发明的滚动轴承载荷分布与形变。

    图4为本发明的滚动轴承的单个滚动体接触的等效刚度与等效阻尼。

    图5为本发明的滚动轴承的滚动体与径向的位置夹角时外圈与滚动体之间的
    油膜压力值。

    图6为本发明的滚动轴承的滚动体与径向的位置夹角时内圈与滚动体接触的
    油膜压力值。

    图7为本发明的滚动轴承的滚动体与径向的位置夹角时外圈与滚动体接触
    的油膜压力值。

    图8为本发明的滚动轴承的滚动体与径向的位置夹角时内圈与滚动体接触
    的油膜压力值。

    图9为本发明的用于实验验证的轴-轴承-基座系统的前四阶固有频率测量结果。

    图10为本发明的用于实验验证的轴-轴承-基座系统基于Dyrobes转子-轴承简化
    分析结果。

    具体实施方式

    下面结合附图对本发明作进一步的说明。

    本发明包括如下步骤:

    第一步:建立滚动体与滚道的EHL接触模型。

    在不考虑滚动体与内外滚道相对滑动的情况下,滚动体、油膜、滚道三者之间的接
    触区域(如图1所示),根据接触弹性变形分为三个区:油膜入口区A、弹性接触区B和油膜出
    口区C。kef、cef为油膜入口区A的油膜刚度与油膜阻尼。kc、kf、cf分别为弹性接触区B滚动体
    与滚道的弹性接触刚度、油膜刚度以及油膜阻尼。由于其油膜出口区C的油膜开始拉伸其影
    响力弱化,相比整体的接触阻尼及刚度影响非常小,这里不考虑油膜出口区C的油膜刚度与
    油膜阻尼。这样主要的刚度与阻尼为kef、cef、kc、kf、cf这五个参数,为了对滚动轴承等效刚
    度与等效阻尼进行计算,首先应分析计算出这五个参数。

    (1)弹性变形接触刚度kc

    Hertz接触理论是根据完全弹性体的静态接触条件得出来的,通常被用来作为异
    向曲面接触副的弹性变形和应力场计算的依据,在滚动体与沟道接触时,其接触面的宽度
    远远小于接触表面的曲率半径,由于将接触视为点接触,则其接触面可认为是椭圆。

    通过计算得出接触椭圆的长轴a与短轴b以及接触形变δ如下:

    <mrow> <mi>a</mi> <mo>=</mo> <msup> <mi>a</mi> <mo>*</mo> </msup> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mn>3</mn> <mi>Q</mi> </mrow> <mrow> <mn>2</mn> <msup> <mi>&Sigma;&rho;E</mi> <mo>&prime;</mo> </msup> </mrow> </mfrac> <mo>)</mo> </mrow> <mrow> <mn>1</mn> <mo>/</mo> <mn>3</mn> </mrow> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>

    <mrow> <mi>b</mi> <mo>=</mo> <msup> <mi>b</mi> <mo>*</mo> </msup> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mn>3</mn> <mi>Q</mi> </mrow> <mrow> <mn>2</mn> <msup> <mi>&Sigma;&rho;E</mi> <mo>&prime;</mo> </msup> </mrow> </mfrac> <mo>)</mo> </mrow> <mrow> <mn>1</mn> <mo>/</mo> <mn>3</mn> </mrow> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>

    <mrow> <mi>&delta;</mi> <mo>=</mo> <msup> <mi>&delta;</mi> <mo>*</mo> </msup> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mn>3</mn> <mi>Q</mi> </mrow> <mrow> <mn>2</mn> <msup> <mi>&Sigma;&rho;E</mi> <mo>&prime;</mo> </msup> </mrow> </mfrac> <mo>)</mo> </mrow> <mrow> <mn>2</mn> <mo>/</mo> <mn>3</mn> </mrow> </msup> <mfrac> <mrow> <mi>&Sigma;</mi> <mi>&rho;</mi> </mrow> <mn>2</mn> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>

    <mrow> <msup> <mi>a</mi> <mo>*</mo> </msup> <mo>=</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mn>2</mn> <mi>E</mi> <mrow> <mo>(</mo> <mi>e</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mi>&pi;</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mn>2</mn> </msup> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>)</mo> </mrow> <mrow> <mn>1</mn> <mo>/</mo> <mn>3</mn> </mrow> </msup> <mo>,</mo> <msup> <mi>b</mi> <mo>*</mo> </msup> <mo>=</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mn>2</mn> <msqrt> <mrow> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mn>2</mn> </msup> </mrow> </msqrt> <mi>E</mi> <mrow> <mo>(</mo> <mi>e</mi> <mo>)</mo> </mrow> </mrow> <mi>&pi;</mi> </mfrac> <mo>)</mo> </mrow> <mrow> <mn>1</mn> <mo>/</mo> <mn>3</mn> </mrow> </msup> <mo>,</mo> <msup> <mi>&delta;</mi> <mo>*</mo> </msup> <mo>=</mo> <mfrac> <mrow> <mn>2</mn> <mi>k</mi> <mrow> <mo>(</mo> <mi>e</mi> <mo>)</mo> </mrow> </mrow> <mi>&pi;</mi> </mfrac> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mn>2</mn> </msup> <mo>)</mo> <mi>&pi;</mi> </mrow> <mrow> <mn>2</mn> <mi>E</mi> <mrow> <mo>(</mo> <mi>e</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>)</mo> </mrow> <mrow> <mn>1</mn> <mo>/</mo> <mn>3</mn> </mrow> </msup> </mrow>

    <mrow> <mfrac> <mn>1</mn> <msup> <mi>E</mi> <mo>&prime;</mo> </msup> </mfrac> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mrow> <mo>(</mo> <mfrac> <mrow> <mn>1</mn> <mo>-</mo> <msubsup> <mi>v</mi> <mn>1</mn> <mn>2</mn> </msubsup> </mrow> <msub> <mi>E</mi> <mn>1</mn> </msub> </mfrac> <mo>+</mo> <mfrac> <mrow> <mn>1</mn> <mo>-</mo> <msubsup> <mi>v</mi> <mn>2</mn> <mn>2</mn> </msubsup> </mrow> <msub> <mi>E</mi> <mn>2</mn> </msub> </mfrac> <mo>)</mo> </mrow> <mo>,</mo> <mi>&Sigma;</mi> <mi>&rho;</mi> <mo>=</mo> <mfrac> <mn>1</mn> <msub> <mi>R</mi> <mrow> <mn>1</mn> <mi>x</mi> </mrow> </msub> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <msub> <mi>R</mi> <mrow> <mn>1</mn> <mi>y</mi> </mrow> </msub> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <msub> <mi>R</mi> <mrow> <mn>2</mn> <mi>x</mi> </mrow> </msub> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <msub> <mi>R</mi> <mrow> <mn>2</mn> <mi>y</mi> </mrow> </msub> </mfrac> <mo>,</mo> </mrow>

    <mrow> <mfrac> <mn>1</mn> <msub> <mi>R</mi> <mi>x</mi> </msub> </mfrac> <mo>=</mo> <mfrac> <mn>1</mn> <msub> <mi>R</mi> <mrow> <mi>x</mi> <mn>1</mn> </mrow> </msub> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <msub> <mi>R</mi> <mrow> <mi>x</mi> <mn>2</mn> </mrow> </msub> </mfrac> <mo>,</mo> <mfrac> <mn>1</mn> <msub> <mi>R</mi> <mi>y</mi> </msub> </mfrac> <mo>=</mo> <mfrac> <mn>1</mn> <msub> <mi>R</mi> <mrow> <mi>y</mi> <mn>1</mn> </mrow> </msub> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <msub> <mi>R</mi> <mrow> <mi>y</mi> <mn>2</mn> </mrow> </msub> </mfrac> </mrow>

    <mrow> <mi>E</mi> <mrow> <mo>(</mo> <mi>e</mi> <mo>)</mo> </mrow> <mo>=</mo> <msubsup> <mo>&Integral;</mo> <mn>0</mn> <mrow> <mi>&pi;</mi> <mo>/</mo> <mn>2</mn> </mrow> </msubsup> <msqrt> <mrow> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mn>2</mn> </msup> <msup> <mi>sin</mi> <mn>2</mn> </msup> <mi>&theta;</mi> </mrow> </msqrt> <mi>d</mi> <mi>&theta;</mi> <mo>,</mo> <mi>k</mi> <mrow> <mo>(</mo> <mi>e</mi> <mo>)</mo> </mrow> <mo>=</mo> <msubsup> <mo>&Integral;</mo> <mn>0</mn> <mrow> <mi>&pi;</mi> <mo>/</mo> <mn>2</mn> </mrow> </msubsup> <mfrac> <mn>1</mn> <msqrt> <mrow> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mn>2</mn> </msup> <msup> <mi>sin</mi> <mn>2</mn> </msup> <mi>&theta;</mi> </mrow> </msqrt> </mfrac> <mi>d</mi> <mi>&theta;</mi> <mo>,</mo> </mrow>

    K(e)=1.5277+0.6023ln(Ry/Rx)、E(e)=1.0003+0.5968(Rx/Ry)

    k=1.0339(Ry/Rx)0.6360、

    曲率符号规定接触的两个面凸表面为正,凹表面为负。E(e),k(e)分别是第一类和
    第二类完全椭圆积分函数,a为椭圆接触面的长轴半径,b为椭圆接触面的短轴半径,Q为接
    触面的相互压力,δ为弹性最大变化量,∑ρ为曲率和,ν是泊松比,E为弹性模量,e为椭圆参
    数。

    利用(3)式可求出单个滚动体的Hertz理论弹性变形,将(3)式进行简化为:

    δ=GQ2/3 (4)

    其中单个滚动体与内圈或外圈的Hertz接触刚度为:

    <mrow> <msub> <mi>k</mi> <mi>c</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mi>d</mi> <mi>Q</mi> </mrow> <mrow> <mi>d</mi> <mi>&delta;</mi> </mrow> </mfrac> <mo>=</mo> <mfrac> <mn>3</mn> <mn>2</mn> </mfrac> <msup> <mi>G</mi> <mrow> <mo>-</mo> <mn>1.5</mn> </mrow> </msup> <msup> <mi>&delta;</mi> <mn>0.5</mn> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>

    从中可以看出,其刚度不是一个常量,它会随位移(或载荷)的变化而变化。

    (2)弹性接触区的油膜刚度kf和油膜阻尼cf

    如图2所示,两个弹性物体间的点接触问题可以视为具有当量主曲率半径Rz,Ry和
    当量弹性模量E'的弹性椭球体与刚性平面的接触。两个表面之间存在润滑油膜,且接触中
    心点o的实际油膜厚度为hc,在油膜压力作用下,接触面产生的弹性变形为δ(x,y),油膜厚
    度的表达式可以写成:

    <mrow> <mi>h</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>h</mi> <mi>c</mi> </msub> <mo>+</mo> <mfrac> <msup> <mi>x</mi> <mn>2</mn> </msup> <mrow> <mn>2</mn> <msub> <mi>R</mi> <mi>x</mi> </msub> </mrow> </mfrac> <mo>+</mo> <mfrac> <msup> <mi>y</mi> <mn>2</mn> </msup> <mrow> <mn>2</mn> <msub> <mi>R</mi> <mi>y</mi> </msub> </mrow> </mfrac> <mo>+</mo> <mi>&delta;</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>-</mo> <mi>&delta;</mi> <mrow> <mo>(</mo> <mn>0</mn> <mo>,</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>

    hc=h0+δ(0,0),h0为刚体中心油膜厚度。

    基于等温条件下的Reynolds方程对油膜的参数计算,等温条件下的Reynolds方程
    普遍形式(假设u2和v2不随x和y变化)如下:

    <mrow> <mtable> <mtr> <mtd> <mrow> <mfrac> <mo>&part;</mo> <mrow> <mo>&part;</mo> <mi>x</mi> </mrow> </mfrac> <mo>(</mo> <mrow> <mfrac> <mrow> <msup> <mi>&rho;h</mi> <mn>3</mn> </msup> </mrow> <mi>&eta;</mi> </mfrac> <mfrac> <mrow> <mo>&part;</mo> <mi>p</mi> </mrow> <mrow> <mo>&part;</mo> <mi>x</mi> </mrow> </mfrac> </mrow> <mo>)</mo> <mo>+</mo> <mfrac> <mo>&part;</mo> <mrow> <mo>&part;</mo> <mi>y</mi> </mrow> </mfrac> <mo>(</mo> <mrow> <mfrac> <mrow> <msup> <mi>&rho;h</mi> <mn>3</mn> </msup> </mrow> <mi>&eta;</mi> </mfrac> <mfrac> <mrow> <mo>&part;</mo> <mi>p</mi> </mrow> <mrow> <mo>&part;</mo> <mi>y</mi> </mrow> </mfrac> </mrow> <mo>)</mo> <mo>=</mo> <mn>6</mn> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>u</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>&times;</mo> <mfrac> <mrow> <mo>&part;</mo> <mrow> <mo>(</mo> <mi>&rho;</mi> <mi>h</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mo>&part;</mo> <mi>x</mi> </mrow> </mfrac> <mo>+</mo> <mn>6</mn> <mrow> <mo>(</mo> <msub> <mi>v</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>v</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mfrac> <mrow> <mo>&part;</mo> <mrow> <mo>(</mo> <mi>&rho;</mi> <mi>h</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mo>&part;</mo> <mi>y</mi> </mrow> </mfrac> <mo>+</mo> <mn>12</mn> <mfrac> <mrow> <mo>&part;</mo> <mrow> <mo>(</mo> <mi>&rho;</mi> <mi>h</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mo>&part;</mo> <mi>t</mi> </mrow> </mfrac> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>

    联立式(6)、式(7)进行求解,在油膜厚度以及油膜力的相互关系的基础上,获得接
    触区的油膜刚度和油膜阻尼。一些学者在Hamrock和Dowson推导的点接触最小油膜厚度的
    公式上直接进行求导来获取刚度值,这样存在两个问题:①油膜刚度是指微量形变下的刚
    度值,②因为在弹性形变的影响下局部油膜厚度的变化不一致,油膜刚度不能直接求导获
    取,而且油 膜阻尼值无法直接获取。因此在联立弹性形变方程和Reynolds方程的基础上,
    运用复合直接迭代法,求出油膜静态下的压力和厚度值,利用扰动方程求解方程组的微动
    量和一阶微动量,最终求得线性近似油膜刚度值和油膜阻尼值。

    在四节点矩形单元上使用双线性函数逼近压力分布求解弹性变形,在油膜压力作
    用下,两个接触表面的法向位移总和为:

    <mrow> <mi>&delta;</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>2</mn> <mrow> <msup> <mi>&pi;E</mi> <mo>&prime;</mo> </msup> </mrow> </mfrac> <mo>&Integral;</mo> <msub> <mo>&Integral;</mo> <mi>&Omega;</mi> </msub> <mfrac> <mrow> <mi>p</mi> <mrow> <mo>(</mo> <mi>&xi;</mi> <mo>,</mo> <mi>&zeta;</mi> <mo>)</mo> </mrow> <mi>d</mi> <mi>&xi;</mi> <mi>d</mi> <mi>&zeta;</mi> </mrow> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>-</mo> <mi>&xi;</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mi>y</mi> <mo>-</mo> <mi>&zeta;</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>

    Ω为求解区域,e为单元区域,pij为单元区域e的中心压力值,λ(x,y)为系数值。

    利用无量纲形式,将式(6)和式(7)化简为点接触润滑无量纲Reynolds方程:

    <mrow> <mfrac> <mo>&part;</mo> <mrow> <mo>&part;</mo> <mi>X</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <mfrac> <mrow> <mover> <mi>&rho;</mi> <mo>&OverBar;</mo> </mover> <msup> <mi>H</mi> <mn>3</mn> </msup> </mrow> <mover> <mi>&eta;</mi> <mo>&OverBar;</mo> </mover> </mfrac> <mfrac> <mrow> <mo>&part;</mo> <mi>P</mi> </mrow> <mrow> <mo>&part;</mo> <mi>X</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mn>1</mn> <msup> <mi>k</mi> <mn>2</mn> </msup> </mfrac> <mfrac> <mo>&part;</mo> <mrow> <mo>&part;</mo> <mi>Y</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <mfrac> <mrow> <mover> <mi>&rho;</mi> <mo>&OverBar;</mo> </mover> <msup> <mi>H</mi> <mn>3</mn> </msup> </mrow> <mover> <mi>&eta;</mi> <mo>&OverBar;</mo> </mover> </mfrac> <mfrac> <mrow> <mo>&part;</mo> <mi>P</mi> </mrow> <mrow> <mo>&part;</mo> <mi>Y</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>A</mi> <mn>1</mn> </msub> <mfrac> <mrow> <mo>&part;</mo> <mrow> <mo>(</mo> <mover> <mi>&rho;</mi> <mo>&OverBar;</mo> </mover> <mi>H</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mo>&part;</mo> <mi>X</mi> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow>

    油膜厚度方程:

    <mrow> <msub> <mi>H</mi> <mrow> <mi>k</mi> <mi>l</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>H</mi> <mn>0</mn> </msub> <mo>+</mo> <msub> <mi>A</mi> <mn>3</mn> </msub> <msubsup> <mi>X</mi> <mi>k</mi> <mn>2</mn> </msubsup> <mo>+</mo> <msub> <mi>A</mi> <mn>4</mn> </msub> <msubsup> <mi>Y</mi> <mi>l</mi> <mn>2</mn> </msubsup> <mo>+</mo> <mover> <mo>&Sigma;</mo> <mi>i</mi> </mover> <mover> <mo>&Sigma;</mo> <mi>j</mi> </mover> <msubsup> <mi>D</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mrow> <mi>k</mi> <mi>l</mi> </mrow> </msubsup> <msub> <mi>P</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow>

    膜厚参数载荷参数速度参数材料参数G=αE1α为粘
    压系数,为最大Hertz压力,这里w为载荷值。

    <mrow> <msub> <mi>A</mi> <mn>1</mn> </msub> <mo>=</mo> <mfrac> <mrow> <mn>12</mn> <msup> <mi>E</mi> <mn>1</mn> </msup> <msubsup> <mi>UR</mi> <mi>x</mi> <mn>3</mn> </msubsup> </mrow> <mrow> <msup> <mi>a</mi> <mn>2</mn> </msup> <msub> <mi>bp</mi> <mi>H</mi> </msub> </mrow> </mfrac> <mo>,</mo> <msub> <mi>A</mi> <mn>3</mn> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <mi>k</mi> </mrow> </mfrac> <mo>,</mo> <msub> <mi>A</mi> <mn>4</mn> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>kR</mi> <mi>x</mi> </msub> </mrow> <mrow> <mn>2</mn> <msub> <mi>R</mi> <mi>y</mi> </msub> </mrow> </mfrac> <mo>,</mo> <msub> <mi>A</mi> <mn>5</mn> </msub> <mo>=</mo> <mfrac> <mrow> <mn>12</mn> <msup> <mi>E</mi> <mn>1</mn> </msup> <msubsup> <mi>R</mi> <mi>x</mi> <mn>2</mn> </msubsup> </mrow> <mrow> <msup> <mi>a</mi> <mn>2</mn> </msup> <msub> <mi>p</mi> <mi>H</mi> </msub> </mrow> </mfrac> <mo>,</mo> <msubsup> <mi>D</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mrow> <mi>k</mi> <mi>l</mi> </mrow> </msubsup> <mo>=</mo> <mfrac> <mrow> <msubsup> <mi>C</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mrow> <mi>k</mi> <mi>l</mi> </mrow> </msubsup> <msub> <mi>R</mi> <mi>x</mi> </msub> <msub> <mi>p</mi> <mi>H</mi> </msub> </mrow> <mrow> <mi>a</mi> <mi>b</mi> </mrow> </mfrac> </mrow>

    为无量纲Roelands粘压关系式。

    为无量纲密度方程。

    将式(9)与式(10)联立,并采用三点中心差分格式来替代偏导,经整理可得如下的
    差分方程组:

    <mrow> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>E</mi> <mrow> <mi>k</mi> <mi>l</mi> </mrow> <mi>E</mi> </msubsup> <msub> <mi>P</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> <mo>,</mo> <mi>l</mi> </mrow> </msub> <mo>+</mo> <msubsup> <mi>E</mi> <mrow> <mi>k</mi> <mi>l</mi> </mrow> <mi>N</mi> </msubsup> <msub> <mi>P</mi> <mrow> <mi>k</mi> <mo>,</mo> <mi>l</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>+</mo> <msubsup> <mi>E</mi> <mrow> <mi>k</mi> <mi>l</mi> </mrow> <mi>W</mi> </msubsup> <msub> <mi>P</mi> <mrow> <mi>k</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>l</mi> </mrow> </msub> <mo>+</mo> <msubsup> <mi>E</mi> <mrow> <mi>k</mi> <mi>l</mi> </mrow> <mi>S</mi> </msubsup> <msub> <mi>P</mi> <mrow> <mi>k</mi> <mo>,</mo> <mi>l</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>+</mo> <msubsup> <mi>E</mi> <mrow> <mi>k</mi> <mi>l</mi> </mrow> <mi>o</mi> </msubsup> <msub> <mi>P</mi> <mrow> <mi>k</mi> <mi>l</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mover> <mo>&Sigma;</mo> <mi>i</mi> </mover> <mover> <mo>&Sigma;</mo> <mi>j</mi> </mover> <mrow> <mo>{</mo> <mrow> <msub> <mi>A</mi> <mn>1</mn> </msub> <mrow> <mo>&lsqb;</mo> <mrow> <msub> <mrow> <mo>(</mo> <mfrac> <mrow> <mo>&part;</mo> <mover> <mi>&rho;</mi> <mo>&OverBar;</mo> </mover> </mrow> <mrow> <mo>&part;</mo> <mi>X</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mrow> <mi>k</mi> <mi>l</mi> </mrow> </msub> <msubsup> <mi>D</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mrow> <mi>k</mi> <mi>l</mi> </mrow> </msubsup> <mo>+</mo> <mover> <msub> <mi>&rho;</mi> <mrow> <mi>k</mi> <mi>l</mi> </mrow> </msub> <mo>&OverBar;</mo> </mover> <msubsup> <mi>DX</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mrow> <mi>k</mi> <mi>l</mi> </mrow> </msubsup> </mrow> <mo>&rsqb;</mo> </mrow> </mrow> <mo>}</mo> </mrow> <msub> <mi>P</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>=</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>A</mi> <mn>1</mn> </msub> <msub> <mrow> <mo>(</mo> <mfrac> <mrow> <mo>&part;</mo> <mover> <mi>&rho;</mi> <mo>&OverBar;</mo> </mover> </mrow> <mrow> <mo>&part;</mo> <mi>X</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mrow> <mi>k</mi> <mi>l</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>H</mi> <mn>0</mn> </msub> <mo>+</mo> <msub> <mi>A</mi> <mn>3</mn> </msub> <msubsup> <mi>X</mi> <mi>k</mi> <mn>2</mn> </msubsup> <mo>+</mo> <msub> <mi>A</mi> <mn>4</mn> </msub> <msubsup> <mi>Y</mi> <mi>l</mi> <mn>2</mn> </msubsup> <mo>)</mo> </mrow> <mo>+</mo> <mn>2</mn> <msub> <mi>A</mi> <mn>1</mn> </msub> <msub> <mi>A</mi> <mn>3</mn> </msub> <msub> <mover> <mi>&rho;</mi> <mo>&OverBar;</mo> </mover> <mrow> <mi>k</mi> <mi>l</mi> </mrow> </msub> <msub> <mi>X</mi> <mi>k</mi> </msub> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow>

    求得静态油膜厚度和油膜压力值,在扰动方程基础上,用泰勒级数将油膜厚度与
    油膜压力展开,取坐标原点为静平衡位置,则油膜压力在静态值附近微动时,油膜压力和油
    膜厚度可由下列线性关系式表示:

    <mrow> <mi>p</mi> <mo>=</mo> <msub> <mi>p</mi> <mn>0</mn> </msub> <mo>+</mo> <mi>k</mi> <mi>&Delta;</mi> <mi>h</mi> <mo>+</mo> <mi>c</mi> <mi>&Delta;</mi> <mover> <mi>h</mi> <mo>&CenterDot;</mo> </mover> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>12</mn> <mo>)</mo> </mrow> </mrow>

    h=h0+Δh (13)

    式中:p0为静平衡油膜压力,h0为静平衡油膜厚度,k为油膜近似
    刚度,c为油膜近似阻尼,Δh,均为扰动参数,其量值很小。则等温条件下的Reynolds方
    程简化后的无量纲形式为:

    <mrow> <mtable> <mtr> <mtd> <mrow> <mfrac> <mo>&part;</mo> <mrow> <mo>&part;</mo> <mi>X</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <mfrac> <mrow> <mover> <mi>&rho;</mi> <mo>&OverBar;</mo> </mover> <msup> <mi>H</mi> <mn>3</mn> </msup> </mrow> <mover> <mi>&eta;</mi> <mo>&OverBar;</mo> </mover> </mfrac> <mo>)</mo> </mrow> <mfrac> <mrow> <mo>&part;</mo> <mi>P</mi> </mrow> <mrow> <mo>&part;</mo> <mi>X</mi> </mrow> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <msup> <mi>k</mi> <mn>2</mn> </msup> </mfrac> <mfrac> <mo>&part;</mo> <mrow> <mo>&part;</mo> <mi>Y</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <mfrac> <mrow> <mover> <mi>&rho;</mi> <mo>&OverBar;</mo> </mover> <msup> <mi>H</mi> <mn>3</mn> </msup> </mrow> <mover> <mi>&eta;</mi> <mo>&OverBar;</mo> </mover> </mfrac> <mo>)</mo> </mrow> <mfrac> <mrow> <mo>&part;</mo> <mi>P</mi> </mrow> <mrow> <mo>&part;</mo> <mi>Y</mi> </mrow> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <mrow> <mo>(</mo> <mfrac> <mrow> <mover> <mi>&rho;</mi> <mo>&OverBar;</mo> </mover> <msup> <mi>H</mi> <mn>3</mn> </msup> </mrow> <mover> <mi>&eta;</mi> <mo>&OverBar;</mo> </mover> </mfrac> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mrow> <mfrac> <mrow> <msup> <mo>&part;</mo> <mn>2</mn> </msup> <mi>P</mi> </mrow> <mrow> <mo>&part;</mo> <msup> <mi>X</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <msup> <mi>k</mi> <mn>2</mn> </msup> </mfrac> <mfrac> <mrow> <msup> <mo>&part;</mo> <mn>2</mn> </msup> <mi>P</mi> </mrow> <mrow> <mo>&part;</mo> <msup> <mi>Y</mi> <mn>2</mn> </msup> </mrow> </mfrac> </mrow> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>A</mi> <mn>1</mn> </msub> <mfrac> <mrow> <mo>&part;</mo> <mrow> <mo>(</mo> <mover> <mi>&rho;</mi> <mo>&OverBar;</mo> </mover> <mi>H</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mo>&part;</mo> <mi>X</mi> </mrow> </mfrac> <mo>+</mo> <msub> <mi>A</mi> <mn>5</mn> </msub> <mover> <mi>&rho;</mi> <mo>&OverBar;</mo> </mover> <mfrac> <mrow> <mo>&part;</mo> <mrow> <mo>(</mo> <mover> <mi>&rho;</mi> <mo>&OverBar;</mo> </mover> <mi>H</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mo>&part;</mo> <mi>T</mi> </mrow> </mfrac> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>14</mn> <mo>)</mo> </mrow> </mrow>

    h0与时间变量无关,将式(12)、式(13)无量纲式代入式(14)中,将Δh、
    二次方以及二次方以上的项略掉,并将同次项归一,得到如下三个等式:

    <mrow> <mfrac> <mo>&part;</mo> <mrow> <mo>&part;</mo> <mi>X</mi> </mrow> </mfrac> <mo>(</mo> <mrow> <mfrac> <mrow> <mover> <mi>&rho;</mi> <mo>&OverBar;</mo> </mover> <msup> <msub> <mi>H</mi> <mn>0</mn> </msub> <mn>3</mn> </msup> </mrow> <mover> <mi>&eta;</mi> <mo>&OverBar;</mo> </mover> </mfrac> <mfrac> <mrow> <mo>&part;</mo> <msub> <mi>P</mi> <mn>0</mn> </msub> </mrow> <mrow> <mo>&part;</mo> <mi>X</mi> </mrow> </mfrac> </mrow> <mo>)</mo> <mo>+</mo> <mfrac> <mn>1</mn> <msup> <mi>k</mi> <mn>2</mn> </msup> </mfrac> <mfrac> <mrow> <mo>&part;</mo> <mi>k</mi> </mrow> <mrow> <mo>&part;</mo> <mi>Y</mi> </mrow> </mfrac> <mo>(</mo> <mrow> <mfrac> <mrow> <mover> <mi>&rho;</mi> <mo>&OverBar;</mo> </mover> <msup> <msub> <mi>H</mi> <mn>0</mn> </msub> <mn>3</mn> </msup> </mrow> <mover> <mi>&eta;</mi> <mo>&OverBar;</mo> </mover> </mfrac> <mfrac> <mrow> <mo>&part;</mo> <msub> <mi>P</mi> <mn>0</mn> </msub> </mrow> <mrow> <mo>&part;</mo> <mi>Y</mi> </mrow> </mfrac> </mrow> <mo>)</mo> <mo>=</mo> <msub> <mi>A</mi> <mn>1</mn> </msub> <mfrac> <mrow> <mo>&part;</mo> <mrow> <mo>(</mo> <mover> <mi>&rho;</mi> <mo>&OverBar;</mo> </mover> <msub> <mi>H</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <mo>&part;</mo> <mi>X</mi> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>15</mn> <mo>)</mo> </mrow> </mrow>

    <mrow> <mtable> <mtr> <mtd> <mrow> <mfrac> <mo>&part;</mo> <mrow> <mo>&part;</mo> <mi>X</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <mrow> <mfrac> <mrow> <mover> <mi>&rho;</mi> <mo>&OverBar;</mo> </mover> <msup> <msub> <mi>H</mi> <mn>0</mn> </msub> <mn>3</mn> </msup> </mrow> <mover> <mi>&eta;</mi> <mo>&OverBar;</mo> </mover> </mfrac> <mrow> <mo>(</mo> <mrow> <mfrac> <mrow> <mo>&part;</mo> <mi>k</mi> </mrow> <mrow> <mo>&part;</mo> <mi>X</mi> </mrow> </mfrac> <mo>+</mo> <mfrac> <mn>3</mn> <msub> <mi>H</mi> <mn>0</mn> </msub> </mfrac> <mfrac> <mrow> <mo>&part;</mo> <msub> <mi>P</mi> <mn>0</mn> </msub> </mrow> <mrow> <mo>&part;</mo> <mi>X</mi> </mrow> </mfrac> </mrow> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <mfrac> <mn>1</mn> <msup> <mi>k</mi> <mn>2</mn> </msup> </mfrac> <mfrac> <mo>&part;</mo> <mrow> <mo>&part;</mo> <mi>Y</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <mrow> <mfrac> <mrow> <mover> <mi>&rho;</mi> <mo>&OverBar;</mo> </mover> <msup> <msub> <mi>H</mi> <mn>0</mn> </msub> <mn>3</mn> </msup> </mrow> <mover> <mi>&eta;</mi> <mo>&OverBar;</mo> </mover> </mfrac> <mrow> <mo>(</mo> <mrow> <mfrac> <mrow> <mo>&part;</mo> <mi>k</mi> </mrow> <mrow> <mo>&part;</mo> <mi>Y</mi> </mrow> </mfrac> <mo>+</mo> <mfrac> <mn>3</mn> <msub> <mi>H</mi> <mn>0</mn> </msub> </mfrac> <mfrac> <mrow> <mo>&part;</mo> <msub> <mi>P</mi> <mn>0</mn> </msub> </mrow> <mrow> <mo>&part;</mo> <mi>Y</mi> </mrow> </mfrac> </mrow> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>A</mi> <mn>1</mn> </msub> <mfrac> <mrow> <mo>&part;</mo> <mover> <mi>&rho;</mi> <mo>&OverBar;</mo> </mover> </mrow> <mrow> <mo>&part;</mo> <mi>X</mi> </mrow> </mfrac> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>16</mn> <mo>)</mo> </mrow> </mrow>

    <mrow> <mfrac> <mo>&part;</mo> <mrow> <mo>&part;</mo> <mi>X</mi> </mrow> </mfrac> <mo>(</mo> <mrow> <mfrac> <mrow> <mover> <mi>&rho;</mi> <mo>&OverBar;</mo> </mover> <msup> <msub> <mi>H</mi> <mn>0</mn> </msub> <mn>3</mn> </msup> </mrow> <mover> <mi>&eta;</mi> <mo>&OverBar;</mo> </mover> </mfrac> <mfrac> <mrow> <mo>&part;</mo> <mi>c</mi> </mrow> <mrow> <mo>&part;</mo> <mi>X</mi> </mrow> </mfrac> </mrow> <mo>)</mo> <mo>+</mo> <mfrac> <mn>1</mn> <msup> <mi>k</mi> <mn>2</mn> </msup> </mfrac> <mfrac> <mo>&part;</mo> <mrow> <mo>&part;</mo> <mi>Y</mi> </mrow> </mfrac> <mo>(</mo> <mrow> <mfrac> <mrow> <mover> <mi>&rho;</mi> <mo>&OverBar;</mo> </mover> <msup> <msub> <mi>H</mi> <mn>0</mn> </msub> <mn>3</mn> </msup> </mrow> <mover> <mi>&eta;</mi> <mo>&OverBar;</mo> </mover> </mfrac> <mfrac> <mrow> <mo>&part;</mo> <mi>c</mi> </mrow> <mrow> <mo>&part;</mo> <mi>Y</mi> </mrow> </mfrac> </mrow> <mo>)</mo> <mo>=</mo> <mfrac> <mrow> <mn>12</mn> <msup> <mi>E</mi> <mn>1</mn> </msup> <msubsup> <mi>R</mi> <mi>x</mi> <mn>2</mn> </msubsup> </mrow> <mrow> <msub> <mi>P</mi> <mi>H</mi> </msub> <msup> <mi>a</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mover> <mi>&rho;</mi> <mo>&OverBar;</mo> </mover> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>17</mn> <mo>)</mo> </mrow> </mrow>

    利用直接复合迭代法计算式(16)和式(17),得出油膜刚度kf和油膜阻尼cf。

    (3)油膜入口区油膜刚度kef与油膜阻尼cef

    简化分析入口区油膜阻尼和油膜刚度,假设以下条件:

    ①由于滚动轴承中滚动体与滚道接触面的形状是狭长的椭圆,因此滚动方向的短
    半轴b远比长半轴a小,所以接触面形状可以近似于同等椭圆长径与短径的矩形接触面,并
    忽略接触区域边缘上的漏油。

    ②滚动体与滚道之间的间隙为抛物线且x≥b,这里R为y方向的
    综合曲率半径。

    ③轴承中油膜厚度大于滚动体、内外圈的表面粗糙度。

    ④惯性力的值小于粘性力。

    ⑤忽略重力影响。

    ⑥考虑到Reynolds方程适用性,可认为粘度是恒定的。

    ⑦忽略接触过程中气穴现象产生的反压力。

    根据上述假设,将Reynolds方程式(7)简化为:

    <mrow> <mfrac> <mo>&part;</mo> <mrow> <mo>&part;</mo> <mi>x</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <mfrac> <mrow> <msup> <mi>&rho;h</mi> <mn>3</mn> </msup> </mrow> <mi>&eta;</mi> </mfrac> <mfrac> <mrow> <mo>&part;</mo> <mi>p</mi> </mrow> <mrow> <mo>&part;</mo> <mi>x</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mo>=</mo> <mn>6</mn> <msub> <mi>u</mi> <mi>s</mi> </msub> <mfrac> <mrow> <mo>&part;</mo> <mrow> <mo>(</mo> <mi>&rho;</mi> <mi>h</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mo>&part;</mo> <mi>x</mi> </mrow> </mfrac> <mo>+</mo> <mn>12</mn> <msub> <mi>u</mi> <mi>z</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>18</mn> <mo>)</mo> </mrow> </mrow>

    其中us=u1+u2,uz为法向挤压速度。η与ρ均为定值。对式(18)中的x积分,并考虑到
    sommerfeld条件和半sommerfeld条件可得:

    <mrow> <msub> <mi>Q</mi> <mi>e</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mn>2</mn> <mi>&eta;</mi> <mi>u</mi> <mi>R</mi> <mi>L</mi> </mrow> <msub> <mi>h</mi> <mn>0</mn> </msub> </mfrac> <mo>+</mo> <mfrac> <mn>3</mn> <msqrt> <mn>2</mn> </msqrt> </mfrac> <mfrac> <mrow> <msup> <mi>&pi;v&eta;R</mi> <mn>1.5</mn> </msup> <mi>L</mi> </mrow> <msubsup> <mi>h</mi> <mn>0</mn> <mn>1.5</mn> </msubsup> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>19</mn> <mo>)</mo> </mrow> </mrow>

    η为油膜粘度、u滚动体圆周速度、L接触椭圆区域长轴、v油膜进入速度。这样可得
    入口区的油膜阻尼:

    <mrow> <msub> <mi>c</mi> <mrow> <mi>e</mi> <mi>f</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>dQ</mi> <mi>e</mi> </msub> </mrow> <mrow> <mi>d</mi> <mi>v</mi> </mrow> </mfrac> <mo>=</mo> <mfrac> <mn>3</mn> <msqrt> <mn>2</mn> </msqrt> </mfrac> <mfrac> <mrow> <msup> <mi>&pi;&eta;R</mi> <mn>1.5</mn> </msup> <mi>L</mi> </mrow> <msubsup> <mi>h</mi> <mn>0</mn> <mn>1.5</mn> </msubsup> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>20</mn> <mo>)</mo> </mrow> </mrow>

    通过数值仿真,发现油膜入口区的形变(包括接触面的弹性变形和油膜的厚度变
    化)相比弹性接触区的形变是十分小的,因此,在载荷微变的情况下,油膜入口区的刚度kef
    可以忽略不计。将计算得到的cef、cf、kc、kf带入式(19)、式(20)最终获得滚动体与滚道相接
    触的 接触刚度k和接触阻尼c:

    <mrow> <mi>k</mi> <mo>=</mo> <mfrac> <mrow> <msub> <mi>k</mi> <mi>c</mi> </msub> <msub> <mi>k</mi> <mi>f</mi> </msub> </mrow> <mrow> <msub> <mi>K</mi> <mi>c</mi> </msub> <mo>+</mo> <msub> <mi>k</mi> <mi>f</mi> </msub> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>21</mn> <mo>)</mo> </mrow> </mrow>

    c=cef+cf (22)

    第二步:计算滚动轴承等效阻尼与等效刚度。

    通过计算得到的滚动体与滚道相接触的接触刚度k和接触阻尼c,计算轴承径向等
    效刚度krc与等效阻尼cre??悸枪龆岢兄泄龆逵牍龅乐浯嬖诰断蛴蜗?,在轴承径向载
    荷Fr下,轴承的内外圈套产生的相对位移为δr,如图3所示。

    在平衡条件下,内圈的径向载荷必须等于滚动体载荷的竖直分量之和:


    这里为单个滚动体与径向的位置夹角,
    注意积分转换后的是连续,α为修正系数。

    <mrow> <msub> <mi>Q</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> <mo>=</mo> <mfrac> <msub> <mi>F</mi> <mi>r</mi> </msub> <mrow> <msub> <mi>ZJ</mi> <mi>r</mi> </msub> <mrow> <mo>(</mo> <mi>&epsiv;</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>24</mn> <mo>)</mo> </mrow> </mrow>


    联立式(24)、式(25)式即可获得滚动体2与滚道的载荷分布,从而求出单个滚动体
    2与内圈1、外圈3的接触刚度和阻尼,如图4所示,求出单个滚动体2与内圈1、外圈3相接触的
    接触刚度与接触阻尼:

    <mrow> <msub> <mi>k</mi> <mi>e</mi> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>k</mi> <mn>1</mn> </msub> <msub> <mi>k</mi> <mn>2</mn> </msub> </mrow> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>k</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>26</mn> <mo>)</mo> </mrow> </mrow>

    <mrow> <msub> <mi>c</mi> <mi>e</mi> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>c</mi> <mn>1</mn> </msub> <msub> <mi>c</mi> <mn>2</mn> </msub> </mrow> <mrow> <msub> <mi>c</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>c</mi> <mn>2</mn> </msub> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>27</mn> <mo>)</mo> </mrow> </mrow>

    将各个滚动体2与内圈1、外圈3相接触的接触刚度与接触阻尼联立,求得滚动轴承
    的径向等效刚度krc与等效阻尼cre:



    同理,这里的为单个滚动体与径向的位置夹角。

    实验验证

    使用B&K测试系统,在轴承部件模态的基础上,利用锤击法进行测量,获取其系统
    刚体固有频率,利用公式求得轴承等效阻尼。

    考虑到轴承本身质量轻、体积小,属于部件模态测试,因此不利于传感器的布置,
    且传递函数计算困难等因素,将轴承外圈固定在轴承支座上,将一定质量的轴安装在轴承
    内圈并通过紧锁螺钉加以紧固,这样轴-轴承-基座成为系统。

    采用单点拾振法,在轴的一段布置传感器,轴与轴承座上所标记的白点为锤击点。
    将测量得到的传递函数通过模态分析获得轴承系统部件刚体固有频率。该测试方法是在轴
    承以600rpm连续工作10分钟后进行,测量部件系统所得到的固有频率如图9所示。

    从图9中我们可以得知:部件模态的前几阶固有频率值ωn以及阻尼比ζ,这里取第
    2阶模态固有频率和阻尼比(部件刚体模态)代入cb=2ζmωn中计算,计算得出的阻尼为
    483.2,这与我们理论计算的所得到的值cre=484.0255相差0.17%。

    为对刚度进行验证,通过Dyrobes转子-轴承分析软模拟出双轴承-单转轴下的一
    阶临界转速57.97Hz,如图10所示,中间加入了0.82kg的质量块

    这里对转轴进行了等效,分别获得等效质量m=1.64kg,等效轴承刚度ks=
    219450,等效阻尼ck=121.3705。

    <mrow> <mi>k</mi> <mo>=</mo> <mfrac> <mrow> <mn>2</mn> <msub> <mi>k</mi> <mi>b</mi> </msub> <msub> <mi>k</mi> <mi>s</mi> </msub> </mrow> <mrow> <mn>2</mn> <msub> <mi>k</mi> <mi>b</mi> </msub> <mo>+</mo> <msub> <mi>k</mi> <mi>s</mi> </msub> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>30</mn> <mo>)</mo> </mrow> </mrow>

    <mrow> <mi>c</mi> <mo>=</mo> <mfrac> <mrow> <msub> <mi>c</mi> <mrow> <mi>r</mi> <mi>e</mi> </mrow> </msub> <msub> <mi>c</mi> <mi>s</mi> </msub> </mrow> <mrow> <msub> <mi>c</mi> <mrow> <mi>r</mi> <mi>e</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>c</mi> <mi>s</mi> </msub> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>31</mn> <mo>)</mo> </mrow> </mrow>

    <mrow> <mn>2</mn> <mi>&pi;</mi> <mi>f</mi> <mo>=</mo> <msqrt> <mrow> <mn>1</mn> <mo>-</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mi>c</mi> <mrow> <mn>2</mn> <msqrt> <mrow> <mi>m</mi> <mi>k</mi> </mrow> </msqrt> </mrow> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mo>&CenterDot;</mo> <msqrt> <mfrac> <mi>k</mi> <mi>m</mi> </mfrac> </msqrt> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>32</mn> <mo>)</mo> </mrow> </mrow>

    联立式(30)、(31)、(32),求得轴承径向动态刚度为kb=5.4743×107,这与我们理
    论计算 值krc=5.3067×107的偏差为3.06%。

    关于本文
    本文标题:一种滚动轴承等效刚度与等效阻尼的计算方法.pdf
    链接地址://www.4mum.com.cn/p-6027307.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
  • pk10公式计划软件安卓 吉林时时快三预测 11选5模拟投注器 德国赛车人工计划 一分快三大小单双回血技巧 福少时时彩 足彩二串一最牛技巧 重庆时时五星走势图彩经网 pk10分析软件手机版 北京pk10走势规律分析 双色球大乐透软件 彩票站转让协议书 2019女篮亚洲杯赛程 千里马计划软件ios 刮刮彩 11选5三胆拖五多少钱