流体力学雷诺方程的推导

主要参数R= 20mm, L=40 mm, n=1000 rpm,

=0.3, c=2 mm.

各种流体润滑问题都涉及在狭小间隙中的流体粘性流动，描写这种物理现象的基本方程为雷诺方程，他的普遍形式是

h2ph2phhh()(）6(UV2) xxyyxyt这个椭圆形的偏微分方程仅仅对于特殊的间隙形状才可能求得解析解，而对于复杂的几何形状或者工况条件下的问题，无法用解析方法求得精确解。随着迅速发展的点算技术，数值算法成为求解润滑问题的有效途径。

数值法师讲偏微分方程转化为代数方程组的变换方法。它的一般原则是：首先将求解域划分成有限个数的单元，并使每一个单元充分的微小。以至于可以认为在各单元内的未知量（本人毕业设计中设油膜压力为P）相等或者依照线性变化，而不会造成很大的误差。然后，通过物理分析或数学变换方法，将求解的偏微分方程写成离散形式，即使将它转化成一组线性代数方程。该代数方程组表示了各个单元的待求未知量于周围各单元未知量的关系。最后根据消去法或者迭代法求解代数方程组，从而求得整个求解域上的未知量。

用来求解雷诺方程的数值方法很多，最常用的是有限元差分方法、有限元法和边界元法，这些方法都是将求解域划分成许多个单元，但是处理方法各不相同。在有限差分法和有限元法中，代替基本方程的函数在求解域内是近似的，但完全满足边界条件。而边界元法所用的函数在求解域内完全满足基本方程，但是在边界上则近似的满足边界条件。

一、雷诺方程的数值解法

根据边界条件求解雷诺方程，这在数学上称为边值问题。

首先将所求解的偏微分方程无量纲化。这样做的目的是减少自变量和因变量的数目，同时用无量纲参数表示的解具有通用性。

然后，将求解域划分成等距的或者不等距的网格，如图1-1为等距网格。

图1-1

沿轴向将Y划分为8个等距区间，沿周向从0到2划分为12个等距区间。这样在Y方向有13个节点，方向有9个节点，总计139117个节点。则1，Y1。

68精品资料 欢迎下载

有限差分法

如果用P代表所求的未知量例如油膜压力，则变量P在整个域中的分布可以用各节点的P值来表示。根据差分原理，任意节点O(i, j)的一阶和二阶偏导数都可以由其周围的节点变量值来表示。

如图1-2所示，如果采用中差分公式，则变量P在O(i, j)点的偏导数为

图.1-2

（pi1,jpi1,jp） i,j2（1-1）

pi,j1pi,j1p（） i,jy2ypi1,jpi1,j2pi,j2p（2)i,j 2() （1-2）

pi,j1pi,j12pi,j2p（2)i,j y(y)2 以P为润滑膜压力，雷诺方程的二维二阶偏微分方程的标准形式为：

2P2PPPA2B2CDE （1-3） YY其中A,B,C,D和E都为已知量。然后将上述方程应用到各个节点，根据中差分公式（1-1）和（1-2）用差商代替偏导数，即可求得各个节点的变量pi.j于相邻各个节点变量的关系。这种关系可以写成：

精品资料 欢迎下载

pi,jCNpi,j1CSpi,j1CEpi1,jCWpi1,jG （1-4）

其中

CNCSCECWGKBy2B(y2A(2A(2EKA2(2(D2y)/K)/K)/K （1-5） )/KD2yC2C2By2）式（1-4）中各系数值随节点位置而改变。

方程（1-4）是有限差分法的计算方程，对于每个节点都可以写出一个方程，而在边界上的节点变量应满足边界条件，它们的数值是已知量。这样，就可以求得一组线性代数方程。方程与未知量数目相一致，所以可以求解。采用消去法或者迭代法求解代数方程组，并使计算结果满足一定的收敛精度，最终求得整个求解域上各节点的变量值。

求解代数方程使用迭代法求解。

1、雷诺方程的无量纲化 定常雷诺方程

h3ph3ph()()6u (2-1) xxyyx将轴承表面沿平面展开，如图1-1所示，并代入xR,dxRd. 得

h3p2ph3h()6u

RRdy2Rd等式两边同时乘以R2 则雷诺方程变为

2p3dh3p(h)h6uR (2-2)

y2d若令

6uRyYL/2,(2R/L),hc(1cos)Hc,pP2

c2精品资料 欢迎下载

代入后得

222336uRP3326uRP(Hc)HcR() 222ccYL6uR化简得

dHd

2R22PdH3P3(H)H（) LY2d（2R/L)2代入得 将 2dH3P3P(H)H （2-3） Y2d由

hc(1cos)Hc

得

H1cos

代入（2-3）式，得

2P2P2P3-3(sin)HHH

2Y22再次化简得无量纲雷诺方程

d(1cos)d

-3(sin)P2P2Psin-3 （2-4） 1cos2Y2（1cos）R为轴承半径，L为轴承长度，为偏心e/c率，e为偏心距，c为半径间隙，采用有限元差分法进行迭代计算。

式（1-4）为标准形式，参考标准式（1-3）可求得标准式中A,B,C,D,E的值。

A1，B,C将以上各值代入式（1-5）求得

3sinsin,D0,E 31cos(1cos)精品资料 欢迎下载

2CN（22Y2)2CS（22Y2) 2(1cos)3sinCE22(1cos)2(1cos)3sinCW22(1cos)3sin2Y2G(1cos)32(2Y2)22K2(Y)2Y2将已知值代入式（1-4）得

P2i（22Y2)P2,ji,j1（22Y2)Pi,j12(1cos)3sin22(1cos)Pi1.j 2(1cos)3sin22(1cos)Pi1,j 3sin2Y2(1cos)32(2Y2) 将

（2R/L)2(220/40)21,0.3

代入式（2-5）得迭代方程:

P2（22Y2)P2i,ji,j1（22Y2)Pi,j1 2(10.3cos)0.9sin22(10.3cos)Pi1.j

2(10.3cos)0.9sin22(10.3cos)Pi1,j 0.9sin2Y2(10.3cos)32(2Y2)

将116，Y8代入上式中，得

2-5） （

精品资料 欢迎下载

Pi,j0.9Pi,j10.9Pi,j1

2(10.3cos)0.47sinPi1.j

0.(10.3cos)2(10.3cos)0.47sinPi1.j

0.(10.3cos)0.012sin(10.3cos)3 （2-6）

上式为最终迭代方程。 边界问题：

将轴承表面沿平面展开，如图2-1

图.2-1

对于径向轴承，方程（2-4）中两个自变量的变化范围是：在轴承中间断面上Y=0：在边缘上 Y=1。而在0到2之间变化，这一问题的边界条件为：

（1）轴向方向

在边缘Y=1处，P=0；在中间断面Y=0上，（2）周向方向

P0. YP0按雷诺边界条件：油膜起点在0处，取P=0；油膜终点在发散区间内符合P=0及的地方。