用瑞利_李兹法求解瞬时非稳态滑动轴承油膜力的新算法.pdfVIP

第22卷第6期 　 计 算 力 学 学 报 　 V o l. 2 2 , N o . 6 2005 年 12 月 Ch in ese Journa l of Com puta t iona l M echan ics D ecem b er 2005 文章编号:(2005) 用瑞利李兹法求解瞬时非稳态滑动轴承 油膜力的新算法 肖忠会, 　沈光琰, 　郑铁生 , 　张　文 (复旦大学 力学与工程科学系, 上海 200433) 摘　要: 提出了一个利用瑞利李兹方法求解 R eyno ld s 边界条件下非稳态滑动轴承油膜力的近似算法, 充分利用 油膜力分布的特性, 用双曲余弦函数来表示油膜力的轴向压力分布, 而用多项式函数插值法来求解油膜力的周向 压力分布, 并同时计算出油膜力的破裂边界。算例表明本算法达到了相当高的精度, 可用于转子系统的非线性数 值分析, 能大大降低数值求解瞬态油膜力的计算时间。 关键词: 有限轴承; 油膜力; R eyno ld s 边界条件; 瑞利- 李兹法 中图分类号: TH 133. 3　　　文献标识码: A 法, 则系统动力相应的计算量几乎全部由油膜力的 1　引 言 计算所占据, 即使计算机的运行速度已经极大地提 刚性支承的转子, 由于其支承部分不能向系统 高, 这一计算量也是难以容忍的。因此 目前多数非 提供足够的阻尼, 致使转子无法平稳地越过临界转 [4 ] 线性分析模型多采用短轴承理论 或长轴承理 速而只能限制在低速范围内工作, 所以对于一大类 [ 5 ] 论 来简化油膜力的计算。这些理论都是基于一些 要求在超临界转速以上运行的高速转子, 或重载转 假设, 略去了计算油膜力的雷诺方程中的某些项, 子, 其支承大多数都必须实现流体润滑。滑动轴承 从而可以简化雷诺方程, 求出油膜力的解析解。根 中的油膜力在旋转机械中主要起着提供支承力, 抑 据文献[ 6 ] 只有当轴承的长径比小于 0. 2 时才能采 制振动的作用。油膜力的计算是通过求解一定边界 用短轴承假设, 当长径比大于 2. 5 时才能用长轴承 条件下的雷诺方程得到的。传统的轴承转子系统 假设, 在实际情况下, 很大一部分轴承并不满足这 的动力学分析都是借助于线性化的描述, 基于小扰 些假设, 此时如果还是运用短轴承或长轴承理论下 动的假设将非线性系统简化为一个用八个刚度, 阻 的油膜力, 就会给系统的分析带来很大误差, 甚至 [ 1 ] 导致分析的彻底失败。而且实际的轴承都是多瓦有 尼系数表征的线性模型 , 这时只需要计算轴心在 静平衡位置的油膜力和八个动特性参数, 因此可以 限长的, 想要获得油膜力的具有足够精度的解析解 [2 ] 是非常困难的。本文考虑的边界条件是R eyno ld s 用有限元法或有限差分法配合充零算法 或互补 [3 ] 边界条件, R eyn