荷作用下滑动轴承轴心轨迹计算仿真.docVIP

. . 　　 本文档下载自 HYPERLINK / 文库下载网，内容可能不完整，您可以点击以下网址继续阅读或下载：  HYPERLINK /doc/c2e87d0af78a6529647d5311.html /doc/c2e87d0af78a6529647d5311.html 荷作用下滑动轴承轴心轨迹计算仿真 很好的资料 第３０卷第１期 ２０１０年２月 振动、测试与诊断 ＪｏｕｒｎａｌｏｆＶｉｂｒａｔｉｏｎ．Ｍｅａｓｕｒｅｍｅｎｔ＆Ｄｉａｇｎｏｓｉｓ Ｖ０１．３０Ｎｏ．１ Ｆｅｂ．２０１０ 瞬变载荷作用下滑动轴承轴心轨迹计算仿真。 马金奎， 路长厚， 陈淑江 （山东大学高效洁净机械制造教育部重点实验室济南，２５００６１） 摘要为了研究动压滑动轴承在瞬变载荷作用下的润滑状况，在考虑轴颈惯性力和非线性油膜力的基础上，建立了滑动轴承轴心轨迹的运动方程，提出了采用欧拉法求解有限长滑动轴承瞬时轴心轨迹的数值方法。分别计算了在阶跃、矩形脉冲和正弦脉冲瞬变载荷作用下轴心轨迹、最大油膜压力及最小油膜厚度的变化规律。分析结果表明，在瞬变载荷作用时，轴心轨迹、最大油膜压力及最小油膜厚度都有较大的变化并呈现出一定的振荡过程。由于脉冲载荷的作用时问有限，随着其消失，轴心乃收敛于原平衡位置，而阶跃载荷则使轴心收敛于新的平衡位置。 关键词滑动轴承润滑仿真轴心轨迹瞬变载荷 中图分类号ＴＨｌ３３．３ ＴＨｌｌ７．２ 有限长滑动轴承瞬时轴心轨迹的数值方法。由初始 引言 滑动轴承的轴心轨迹直接反映轴承的工作状 条件，首先求解该瞬时轴心位置的Ｒｅｙｎｏｌｄｓ方程，根据油膜压力分布积分得到油膜合力，然后求解运动方程，得到轴心的瞬时加速度，再由欧拉方法计算出轴心的下一个位置参数，如此循环迭代，直到得到完整的轴心轨迹。该方法充分考虑了油膜力的非线性，对滑动轴承参数设计及非线性行为分析具有实用价值。 态，通过轴心轨迹可以判定轴承的稳定性［１］，了解轴承的润滑状况，验证轴承设计参数的合理性。当轴承承受静载荷稳定工作时，轴心处于静平衡位置或有小的涡动。如果轴承受到瞬态动态载荷作用，轴心轨迹运动复杂，油膜压力和最小油膜厚度发生较大变化，尤其当最小油膜厚度较小时，轴承和轴会发生摩擦或碰壁［２｛］。因此，计算瞬变动载荷作用下滑动轴承的轴心轨迹和最小油膜厚度，对滑动轴承设计及/doc/c2e87d0af78a6529647d5311.html其故障诊断有重要意义。文献［４］借助于ＣＦＤ软件ＡＤＩＮＤ８．１建立滑动轴承的液固模型，对轴施加ｚ，Ｙ方向的稳态正弦激励，研究了轴心轨迹的变化，并利用ｚ，Ｙ方向的位移响应识别了轴承的非线性特性系数。何芝仙［５３等用多体动力学软件ＡＤＡＭＳ和Ｍａｔｌａｂ联合仿真，研究了矩形冲击载荷作用下考虑轴颈倾斜时轴心轨迹的变化，但该文所选工况接近半速涡动，油膜产生自激振动，不能准确反映载荷对轴心轨迹的影响。文献［６］计算了冲击载荷对油膜压力的影响，没有计算轴心轨迹和最小油膜厚度的变化。 本文基于轴颈惯性力、非线性油膜力和瞬变载荷平衡，建立了滑动轴承轴心轨迹的运动方程。通过求解雷诺方程和运动方程，提出了采用欧拉法求解 １轴承系统模型 １．１运动方程 图ｌ为滑动轴承计算模型。Ｑ为轴承的中心，即ｘＯｂｙ坐标原点；ＯＪ为轴的中心，轴心坐标为（ｚ，ｙ）；轴以叫。旋转，Ｆ。（ｃｏｏｔ），Ｆ，（ｏＪｏｔ）为非线性油膜力在ｚ，Ｙ方向的积分分量；Ｑ为作用于轴上的瞬变载荷；Ｍ为轴的质量。轴心在坐标系ｘＯ。Ｙ中的运动方程为 鬈纂荔；嚣地 Ｍ多＝Ｆ，（∞∥）＋Ｑ， ㈩ 用ＭＣ∞２０除以方程（１）两端，则其无量纲形式为 戈＝口瓦（ｒ）＋一ｑｘ（ｒ）十ｐ 舯ｑ ２蚓ｒｌＬ引３佤属卜瀚； （２） 矿＝口Ｆ，（ｒ）＋己（ｒ） ．国家“九七三”计划资助项目（编号：２００９ＣＢ７２４４０４） 收稿日期：２００９—０８—０６；修改稿收到日期：２００９—１１—１６ 万方数据 第ｌ期马金奎，等：瞬变载荷作用下滑动轴承轴心轨迹计算仿真 图１滑动轴承计算模型 ｒ＝㈤（一ｑｚ面＝％等；ｐ＝南埘㈤＝ 等；Ｌ为轴承长度；７为润滑油黏度；ｃ为轴承半 径间隙；Ｒ为轴承半径。 Ｉ．２雷诺方程 不司压缩、层流状态。/doc/c2e87d0af78a6529647d5311.htmlＦ动载荷作用雷诺方程为】１ ａ［＾３雾］＋丕１３警］＝ ６∞。７耄＋ｌＺｒ／（ｘｃｏｓ口＋５，ｓｉｎ口） （３） 其中：０为圆周计算坐标；ｚ为轴向坐标。 油膜厚度为 ｈ（咿，ｆ）＝Ｃ＋ｘｃｏｓ０＋ｙｓｉｎＯ （４） 石＝南，（文，ｐ）＝警，则式（３）的无量纲形 令无量纲变量Ⅳ＝告，（ｘ，ｙ）＝垒半，ｚ－－一ｉＺ， 弘吐虿Ｊ 式为 ；［日３等］＋（引２曩［Ｈ３差］＝ ６等＋１２（２ｃ。ｓ０＋９ ｓｉｎ口） （５） 无量纲油膜厚度 Ｈ