轴承非线性油膜力的一种变分近似解ξ.pdfVIP

第 19 卷 　第 3 期 应 用 力 学 学 报 Vol. 19 　No. 3 2002 年 9 月 CHINESE JOURNAL OF APPL IED MECHANICS Sep. 2002 文章编号 (2002) 轴承非线性油膜力的一种变分近似解 陈　龙　郑铁生　张　文　马建敏 ( 复旦大学　上海　200433) 摘要 :基于自由边值理论和凸集上的变分方法 ,提出一种求解当轴颈大扰动时实际轴承瞬态油膜力 ( ) 的近似公式。公式中引入一个参数来模拟油膜破裂自由边界 雷诺边界条件 ,则凸集上的泛函极 值问题就转换为求此参数的代数极值问题。对椭圆轴承在轴颈作大范围扰动情况下的计算结果表 明 ,这一方法达到了很高的精度 ,可用于转子轴承系统的非线性动力分析 ,能大量降低数值计算瞬 态油膜力所需的计算量。 关键词 :滑动轴承 ;非定常油膜力 ;变分方法 中图分类号: TH1333 　　　文献标识码 : 　A 油膜力 ,也需要油膜力的Jacobi 矩阵 ,而由油膜力数 1 　引　言 据的插值来获得多变量偏导数构成 Jacobi 矩阵 ,其 精度的保证是一个难点。于是寻求实际轴承油膜力 转子轴承系统动力学分析需要求解轴承的油 具有良好近似性的半解析快速算法的问题又提了出 膜力。对线性系统 , 只需计算轴心在静平衡位置的 来。 油膜力和 8 个动特性系数。因此可以用有限差分或 区别于稳态线性油膜力模型 ,在进行非线性分 有限元法配合冲零算法[1 ] 或互补算法[2 ] 来求得 析时 Reynolds 方程和边界条件都是非稳态的 ,必需 Reynolds 边界条件下的油膜力 ,而后就可以计算转 同时考虑它们的影响。文献[ 5 ]提出了动载轴承非 子轴承系统的动力特性。然而对非线性系统 ,在求 稳态非线性油膜力的一般数学模型只需三个函数表 解系统动力响应的每一步都必须重复计算油膜力。 达 ,并具体给出了无限短全圆柱轴承非稳态油膜力 如果仍采用大规模的数值方法[3 ] ,则系统动力响应 的解析式。短轴承模型油膜力模型完全忽略了 的计算量几乎为油膜力的计算所占据。即使计算机 Reynolds 方程中的周向压力导数项 ,致使周向压力 的运行速度已极大的提高 ,这一计算量也是难以容 分布仅取决于 Reynolds 方程的“载荷”项 ,也就无法 忍的。因此目前多数非线性转子动力学分析仍采用 适应周向边界条件 , 更不能适应 Reynolds 边界条 短轴承模型 ,如文献[ 4 ,5 ] 。无可置疑 ,这种模型虽 件。这一情况与实际轴承的运行状况相差很大 ,必 然简化了油膜力的计算 ,但与实际轴承相去甚远 ,仅 然导致大的误差。实际轴承都是多瓦有限长的 ,要 具有定性意义。最近有文献[ 6 ]提出建立有限轴承 想获得油膜力的具有足够精度的完全解析解是非常 的非线性油膜力数据库 ,运用插值来求油膜力。但 困难的。本文运用自由边值变分原理 ,取短轴承模 非稳态多参数所需的巨大数据量使这一构思的应用 型压力轴向分布的合理成分 ,将二维变分化成一维 , 范围受到限制 ;更何况系统