完整版本流体力学雷诺方程的推导.docVIP

1-1 为等距网格。 主要参数 R= 20mm, L=40 mm, n=1000 rpm, =0.3, c=2 mm. 各种流体润滑问题都涉及在狭小间隙中的流体粘性流动， 描写这种物理现象的基本方程为雷诺方程，他的普遍形式是 ( h2 p ) y ( h2 p ） 6(U h V h 2 h ) x x y x y t 这个椭圆形的偏微分方程仅仅对于特殊的间隙形状才可能求得解析解，而对于复杂的几何形状或者工况条件下的问题，无法用解析方法求得精确解。随着迅速发展的点算技术，数值算法成为求解润滑问题的有效途径。 数值法师讲偏微分方程转化为代数方程组的变换方法。它的一般原则是：首先将求解域划分成有限个数的单元，并使每一个单元充分的微小。以至于可以认为在各单元内的未知量 （本人毕业设计中设油膜压力为 P）相等或者依照线性变化，而不会造成很大的误差。然后，通过物理分析或数学变换方法，将求解的偏微分方程写成离散形式，即使将它转化成一组线 性代数方程。该代数方程组表示了各个单元的待求未知量于周围各单元未知量的关系。最后根据消去法或者迭代法求解代数方程组，从而求得整个求解域上的未知量。 用来求解雷诺方程的数值方法很多，最常用的是有限元差分方法、有限元法和边界元法，这些方法都是将求解域划分成许多个单元，但是处理方法各不相同。在有限差分法和有限元法中，代替基本方程的函数在求解域内是近似的，但完全满足边界条件。而边界元法所用的函数在求解域内完全满足基本方程，但是在边界上则近似的满足边界条件。 一、雷诺方程的数值解法 根据边界条件求解雷诺方程，这在数学上称为边值问题。 首先将所求解的偏微分方程无量纲化。这样做的目的是减少自变量和因变量的数目，同 时用无量纲参数表示的解具有通用性。 然后，将求解域划分成等距的或者不等距的网格，如图 图 1-1 沿轴向将 Y 划分为 8 个等距区间，沿周向从 0到 2 划分为 12 个等距区间。这样在 Y 方向有 13 个节点， 方向有 9 个节点，总计 13 9 117 个节点。则 1 ， Y 1 。 6 8 有限差分法 如果用 P 代表所求的未知量例如油膜压力， 则变量 P 在整个域中的分布可以用各节点的 P 值来表示。 根据差分原理， 任意节点 O(i, j) 的一阶和二阶偏导数都可以由其周围的节点变量值来表示。 如图 1-2 所示，如果采用中差分公式，则变量 P 在 O(i, j) 点的偏导数为 图 .1-2 p pi 1, j pi 1, j （ ）i,j 2 （1-1） p pi, j 1 pi,j 1 （ ） 2 y y i,j 2 p pi 1, j pi 1, j 2pi, j （ 2 )i, j ( )2 （1-2） 2 pi, j 1 pi, j 1 2pi, j （ p2)i, j ( y) 2 y 以 P 为润滑膜压力，雷诺方程的二维二阶偏微分方程的标准形式为： A 2P B 2 P C P D P E （1-3） 2 Y2 Y 其中 A,B,C,D 和 E 都为已知量。然后将上述方程应用到各个节点， 根据中差分公式 （1-1） 和（ 1-2）用差商代替偏导数，即可求得各个节点的变量 pi . j 于相邻各个节点变量的关系。这 种关系可以写成： pi , jC N pi , j 1 C S pi , j 1 CE pi 1 , j CW pi 1, j G （1-4） 其中 C N B ( y C S B ( y C E A ( C W A ( G E K  2 2 2 2  D ) / K y D y ) / K C K （1-5） ) / 2 C K ) / 2 K 2( A B ） 2 y 2 式（ 1-4）中各系数值随节点位置而改变。 方程（ 1-4）是有限差分法的计算方程，对于每个节点都可以写出一个方程，而在边界上的节点变量应满足边界条件，它们的数值是已知量。这样，就可以求得一组线性代数方程。方程与未知量数目相一致，所以可以求解。采用消去法或者迭代法求解代数方程组，并使计算结果满足一定的收敛精度，最终求得整个求解域上各节点的变量值。 求解代数方程使用迭代法求解。 1、雷诺方程的无量纲化 定常雷诺方程 ( h 3 p ) ( h 3 p ) 6 u h (2-1) x x y y x 将轴承表面沿平面展开，如图 1-1 所示，并代入 x R ,dx Rd . 得 ( h 3 p ) 2 ph 3 6 u Rd h R Rd y 2 等式两边同时乘以 R2 则雷诺方程变为 ( h 3 p ) 2 p2 h 3 6 u R dh (2-2) y d 若令 y YL/ 2, (2R/ L