湍流工况下水导轴承承载力及其Jacobin矩阵的一维直接算法.PDFVIP

湍流工况下水导轴承承载力及其 Jacobin 矩阵 的一维直接算法 闫志勇 （复旦大学力学与工程科学系上海 200433 ） 摘要：基于变分原理，把水导轴承湍流工况下修正的 Reynolds 方程转化为一个泛函 极值问题，并提出对应的一维有限元算法。算例表明，本文算法不仅结果精度能得 到保证，而且运算时间大大减小。 关键词：水导轴承； 湍流； 变分原理 Direct solution method for load force and Jacobin matrix of water guide bearing in turbulent regime Yan Zhiyong (Department of Mechanics and Engineering Science, Fudan University, Shanghai, 200433) Abstract: Based on the variational theory, the amended Reynolds equation of water guide bearing is transformed to functional extreme value problem, and one-dimensional finite element algorithm is proposed. The numerical results show that this method can diminish the computing time remarkably with the precision guaranteed. Keywords: water guide bearing; turbulent; variational principle 1、引言 水导轴承广泛地应用于大型立式水泵机组、核主泵上。由于制造、安装误差、来流等因 素，运行时叶轮和主轴承受径向力作用，水导轴承起着稳定主轴和叶轮的作用。 对采用水这种低粘度润滑剂的轴承而言，由于粘度低，雷诺数一般较大，层流理论往往 不准确，需要采用紊流理论加以修正。在实际工程计算中，以 Peichardt 壁面定律为基础的 Ng-Pan 理论，只适用于 Couette 流为主的高速轻载轴承[1] 。对于以轴向流动为主的Poiseuille 流，Hirs 通过经验关系建立了 Hirs 理论[2] 。张运清在考虑惯性效应的情况下，提出了一种适 用于高压密封和高速重载轴承的复杂流场紊流 Reynolds 方程[3] 。Elord 等提出了可用于求解 Couette 流和 Poiseuille 流同时存在的 Elrod-Ng 理论[4] 。水导轴承紊流计算模式的求解，在实 际应用中往往采取有限差分法、有限体积法、数据库等方法，计算量巨大[3, 5] 。尤其是对于 水导轴承支撑的水轮机械系统进行非线性稳定性分析时，在求解系统的动力响应的每一步都 必须重复计算承载力，显然这一计算量是难以忍受的。 本文基于自由边值变分原理，把水导轴承湍流工况下满足的修正 Reynolds 方程等价于 一个泛函极值问题，选取权函数与线性插值函数乘积的形式作为压力函数，建立对应的一维 有限元变分方程，并采用直接法求解承载力及其对应的 Jacobin 矩阵。 2、基本原理 对于湍流工况下，文献[1]给出了如下修正的非稳态 Reynolds 方程 1  G h3 p  G h3 p  h h 2 (  )  ( z )   (1) R z z t       2   其中 G ,G 分别为无量纲湍流系数。  z 如果水导轴承内以 Poi