2023年流体力学雷诺方程的推导.pdf

精品资料 欢迎下载 主要参数 R= 20mm, L=40 mm, n=1000 rpm, =0.3, c=2 mm. 各种流体润滑问题都涉及在狭小间隙中的流体粘性流动， 描写这种物理现象的基本方程为雷诺方程，他的普遍形式是 )2(6()(22thyhVxhUyphyxphx） 这个椭圆形的偏微分方程仅仅对于特殊的间隙形状才可能求得解析解，而对于复杂的几何形状或者工况条件下的问题，无法用解析方法求得精确解随着迅速发展的点算技术，数值算法成为求解润滑问题的有效途径 数值法师讲偏微分方程转化为代数方程组的变换方法它的一般原则是：首先将求解域划分成有限个数的单元，并使每一个单元充分的微小以至于可以认为在各单元内的未知量（本人毕业设计中设油膜压力为P）相等或者依照线性变化，而不会造成很大的误差然后，通过物理分析或数学变换方法，将求解的偏微分方程写成离散形式，即使将它转化成一组线性代数方程该代数方程组表示了各个单元的待求未知量于周围各单元未知量的关系最后根据消去法或者迭代法求解代数方程组，从而求得整个求解域上的未知量 用来求解雷诺方程的数值方法很多，最常用的是有限元差分方法、有限元法和边界元法，这些方法都是将求解域划分成许多个单元，但是处理方法各不相同。

在有限差分法和有限元法中，代替基本方程的函数在求解域内是近似的，但完全满足边界条件而边界元法所用的函数在求解域内完全满足基本方程，但是在边界上则近似的满足边界条件 一、雷诺方程的数值解法 根据边界条件求解雷诺方程，这在数学上称为边值问题 首先将所求解的偏微分方程无量纲化这样做的目的是减少自变量和因变量的数目，同时用无量纲参数表示的解具有通用性 然后，将求解域划分成等距的或者不等距的网格，如图1-1为等距网格 图 1-1 沿轴向将 Y 划分为 8 个等距区间，沿周向从20 到划分为 12 个等距区间这样在 Y方向有 13 个节点，方向有 9 个节点，总计117913个节点则8161Y， 精品资料 欢迎下载 有限差分法 如果用 P 代表所求的未知量例如油膜压力， 则变量 P 在整个域中的分布可以用各节点的 P值来表示根据差分原理，任意节点 O(i, j)的一阶和二阶偏导数都可以由其周围的节点变量值来表示 如图 1-2所示，如果采用中差分公式，则变量 P 在 O(i, j) 点的偏导数为 图.1-2 2, 1, 1,jijijippp）（ （1-1 ） yppypjijiji21,1,,）（ 2,, 1, 1,22)(2)jijijijipppp（ （1-2 ） 2,1,1,,22)(2)ypppypjijijiji（ 以 P 为润滑膜压力，雷诺方程的二维二阶偏微分方程的标准形式为： EYPDPCYPBPA2222 （1-3 ） 其中 A,B,C,D 和 E 都为已知量。

然后将上述方程应用到各个节点， 根据中差分公式 （1-1 ）和（1-2）用差商代替偏导数，即可求得各个节点的变量jip.于相邻各个节点变量的关系这种关系可以写成： 形状或者工况条件下的问题无法用解析方法求得精确解随着迅速发展的一个单元充分的微小以至于可以认为在各单元内的未知量本人毕业设计数方程该代数方程组表示了各个单元的待求未知量于周围各单元未知量精品资料 欢迎下载 GpCpCpCpCpjiWjiEjiSjiNji, 1, 11,1,, （1-4 ） 其中 ）y( 2/)2(/)2(/)2(/)2(222222BAKKEGKCACKCACKyDyBCKyDyBCWESN （1-5 ） 式（1-4 ）中各系数值随节点位置而改变 方程（1-4）是有限差分法的计算方程，对于每个节点都可以写出一个方程，而在边界上的节点变量应满足边界条件，它们的数值是已知量这样，就可以求得一组线性代数方程方程与未知量数目相一致，所以可以求解采用消去法或者迭代法求解代数方程组，并使计算结果满足一定的收敛精度，最终求得整个求解域上各节点的变量值。

求解代数方程使用迭代法求解 1、雷诺方程的无量纲化 定常雷诺方程 xhuyphyxphx6)()(33 (2-1) 将轴承表面沿平面展开，如图1-1所示，并代入.,RddxRx 得 RdhuhypRdphR6)(3223 等式两边同时乘以2R 则雷诺方程变为 ddhRuhypph6)(3223 (2-2) 若令 226,)cos1 (,)/2 (, 2/cRuPpHcchLRYLy 形状或者工况条件下的问题无法用解析方法求得精确解随着迅速发展的一个单元充分的微小以至于可以认为在各单元内的未知量本人毕业设计数方程该代数方程组表示了各个单元的待求未知量于周围各单元未知量精品资料 欢迎下载 代入后得 2222233233)2(6)6(LYPcRuRcHPcRucH ddHRu6 化简得 ddHYPLRHPH22233)2)(（ 将 2)/2LR（代入得 ddHYPHPH2233)( （2-3 ） 由 Hcch)cos1 ( 得 cos1H 代入（2-3 ）式，得 2232222)sin(3-YPHPHPH dd)cos1 (  再次化简得无量纲雷诺方程 32222cos1sin-cos1)sin( 3-）（YPPP （2-4 ） R 为轴承半径，L 为轴承长度，为偏心ce/率，e为偏心距，c为半径间隙，采用有限元差分法进行迭代计算。

式（1-4 ）为标准形式，参考标准式（1-3 ）可求得标准式中 A,B,C,D,E 的值 3)cos1 (sin, 0,cos1sin3,1EDCBA， 将以上各值代入式（1-5 ）求得 形状或者工况条件下的问题无法用解析方法求得精确解随着迅速发展的一个单元充分的微小以至于可以认为在各单元内的未知量本人毕业设计数方程该代数方程组表示了各个单元的待求未知量于周围各单元未知量精品资料 欢迎下载 22222222322222222)(2)(2)cos1(sin3)cos1(2sin3)cos1(2)cos1(2sin3)cos1(2)2)2YYKYYGCCYCYCWESN（（ 将已知值代入式（1-4 ）得 1,2221,222,)2)2jijijiPYPYP（（ jiP.12)cos1(2sin3)cos1(2 jiP,12)cos1(2sin3)cos1(2 )(2)cos1(sin322223YY （2-5 ） 将 .30, 1)40/202()/222LR（ 代入式（ 2-5 ）得迭代方程: 1,2221,222,)2)2jijijiPYPYP（（ jiP.12)cos.301(2sin.90)cos.301(2 jiP,12)cos.301(2sin.90)cos.301(2 )(2)cos.301(sin.9022223YY 将8161Y，代入上式中，得 形状或者工况条件下的问题无法用解析方法求得精确解随着迅速发展的一个单元充分的微小以至于可以认为在各单元内的未知量本人毕业设计数方程该代数方程组表示了各个单元的待求未知量于周围各单元未知量精品资料 欢迎下载 1,1,,.90.90jijijiPPP jiP.1)cos.301(4.50sin7.40)cos.301(2 jiP.1)cos.301(4.50sin7.40)cos.301(2 3)cos.301(sin12.00 （2-6 ） 上式为最终迭代方程。

边界问题： 将轴承表面沿平面展开，如图 2-1 图.2-1 对于径向轴承，方程（2-4 ）中两个自变量的变化范围是：在轴承中间断面上 Y=0：在边缘上 Y=1而在2到0之间变化，这一问题的边界条件为： （1）轴向方向 在边缘 Y=1 处，P=0；在中间断面 Y=0 上，0YP. （2）周向方向 按雷诺边界条件： 油膜起点在0处， 取 P=0； 油膜终点在发散区间内符合 P=0 及0P的地方 形状或者工况条件下的问题无法用解析方法求得精确解随着迅速发展的一个单元充分的微小以至于可以认为在各单元内的未知量本人毕业设计数方程该代数方程组表示了各个单元的待求未知量于周围各单元未知量。