《流体力学雷诺方程的推导汇编》由会员分享，可在线阅读，更多相关《流体力学雷诺方程的推导汇编（9页珍藏版）》请在装配图网上搜索。
1、流体力学雷诺方程的推导主要参数 R= 20mm, L=40 mm, n=1000 rpm, =0.3, c=2 mm.各种流体润滑问题都涉及在狭小间隙中的流体粘性流动，描写这种物理现象的基本方程 为雷诺方程，他的普遍形式是-p)-p)6(U V 2)x x y yx y t这个椭圆形的偏微分方程仅仅对于特殊的间隙形状才可能求得解析解，而对于复杂的几何形 状或者工况条件下的问题，无法用解析方法求得精确解。随着迅速发展的点算技术，数值算 法成为求解润滑问题的有效途径。数值法师讲偏微分方程转化为代数方程组的变换方法。它的一般原则是：首先将求解域 划分成有限个数的单元，并使每一个单元充分的微小。以至于
2、可以认为在各单元内的未知量 (本人毕业设计中设油膜压力为 P)相等或者依照线性变化，而不会造成很大的误差。然 后，通过物理分析或数学变换方法，将求解的偏微分方程写成离散形式，即使将它转化成一 组线性代数方程。该代数方程组表示了各个单元的待求未知量于周围各单元未知量的关系。 最后根据消去法或者迭代法求解代数方程组，从而求得整个求解域上的未知量。用来求解雷诺方程的数值方法很多，最常用的是有限元差分方法、有限元法和边界元 法，这些方法都是将求解域划分成许多个单元，但是处理方法各不相同。在有限差分法和有 限元法中，代替基本方程的函数在求解域内是近似的，但完全满足边界条件。而边界元法所 用的函数在求解域
3、内完全满足基本方程，但是在边界上则近似的满足边界条件。一、雷诺方程的数值解法根据边界条件求解雷诺方程，这在数学上称为边值问题。首先将所求解的偏微分方程无量纲化。这样做的目的是减少自变量和因变量的数目，同 时用无量纲参数表示的解具有通用性。然后，将求解域划分成等距的或者不等距的网格，如图1-1为等距网格。图1-1沿轴向将Y划分为8个等距区间，沿周向从0到2划分为12个等距区间。这样在Y方向有13个节点，方向有9个节点,总计13 9 117个节点。则 工丫小 有限差分法如果用P代表所求的未知量例如油膜压力，则变量 P在整个域中的分布可以用各节点的 P值来表示。根据差分原理，任意节点 O(i,j)的
4、一阶和二阶偏导数都可以由其周围的节点变量值来表示。如图1-2所示，如果采用中差分公式，则变量 P在O(i, j)点的偏导数为图.1-2（_p）Pu Pui，j 2(1-1)2(fj2()2/i,jyp,j 1 p,j 12yp 1,j p 1,j 2pFT,j(1-2)p,j1 p,j1 2p,j以P为润滑膜压力，雷诺方程的二维二阶偏微分方程的标准形式为:2PY2cD EY(1-3)其中A,B,C,D和E都为已知量。然后将上述方程应用到各个节点,根据中差分公式（1-1）和（1-2）用差商代替偏导数，即可求得各个节点的变量R.j于相邻各个节点变量的关系。这种关系可以写成:pi, j CN pi,
5、jC Spi, j 1CE pi 1 ,j CW pi 1, j G(1-4)其中CwB 不BypA2DyDy C2C2-(1-5)B_)y2式（1-4）中各系数值随节点位置而改变。方程（1-4）是有限差分法的计算方程，对于每个节点都可以写出一个方程，而在边界上 的节点变量应满足边界条件，它们的数值是已知量。这样，就可以求得一组线性代数方程。方程与未知量数目相一致，所以可以求解。采用消去法或者迭代法求解代数方程组，并使计 算结果满足一定的收敛精度，最终求得整个求解域上各节点的变量值。求解代数方程使用迭代法求解。1、雷诺方程的无量纲化定常雷诺方程将轴承表面沿平面展开，如图-(匚上) xx-(M
6、y匕) y(2-1)1-1所示，并代入R ,dxRd .等式两边同时乘以R2则雷诺方程变为3 J)hRd(h3-)2P2 yh3dh6u R d(2-2)若令y YL/2,(2R/L)2,h c(1cos) Hc, pc6u RPrc代入后得(H 3c3 6uc3R2 6u R2cdH6u Rd化简得(HdH d将（2R/L）2代入得3 P(H 3 上)dH d(2-3)h c(1cos )Hccos代入（2-3）式，得-3( sin)H2 22P22PY7再次化简得无量纲雷诺方程-3(R为轴承半径，L为轴承长度,限元差分法进行迭代计算。d(1cos ) dsin ) P1 cos为偏心2P-
7、22Psin(1 cos )3(2-4)e/c率e为偏心距，c为半径间隙,采用有式（1-4）为标准形式，参考标准式（1-3）可求得标准式中A,B,C,D,E的值。A 1, B ,C3 sinE，d QEsin(1 cos )3将以上各值代入式（1-5）求得22Y)22 (1 cos )3 sin22 (1 cos )2 (1 cos )3 sin22 (1 cos )Y 2)3 sin(1 cos 2T 2(2 Y 2 )将已知值代入式（1-4）得Pi,j2222 Y 2 ) Pi，j 122 Y 2 ) Pi，j2 (1 cos )3 sin-2 Pi 1.j2(1 cos )2 (1 co
8、s )3 sin-2 Pi 1 j22 (1 cos )(2-5)3 sin2 Y 2(1 cos2(2 Y 2 )将2_ _ 2_(2R/L)2(220 /40)21,0.3代入式(2-5)得迭彳t方程： 22Pi,j 2- Pi,j 12- Pi,j 12(2 Y 2 )2(2 Y 2 )2 (10 .3 cos )0 .9 sinZ2P i 1 . j22 (10.3 cos)j2(10.3 cos)0.9sin24 Pi 1,j2(10.3 cos )_220 .9 sin Y(10 .3 cos )32(2 Y 2 ),11, ,_将 -,Y 1代入上式中，得 68Pi,j0.9 P
9、-0.9 P-12(10 .3 cos )0 .4 7 sinZZ-ZP i 1 . j0 .5 4 (10 .3 cos )j2 (10 .3 cos )0 .4 7 sinP i 1 . j0 .5 4 (10 .3 cos )0 .0 12 sinZ3(10 .3 cos )(2-6)上式为最终迭代方程。边界问题：将轴承表面沿平面展开，如图 2-1图.2-1对于径向轴承，方程(2-4)中两个自变量的变化范围是：在轴承中间断面上丫=0:在边缘上丫=1。而 在0到2之间变化，这一问题的边界条件为：(1)轴向方向在边缘Y=1处，P=0;在中间断面 Y=0上，0.(2)周向方向按雷诺边界条件：油膜起点在0处，取P=0;油膜终点在发散区间内符合 P=0及P 0的地万。
流体力学雷诺方程的推导汇编
