摩擦学原理(第9章数值方法)探究.ppt

第九章 润滑分析常用数值方法 (common numerical methods for lubrication analysis ) 各种流体润滑问题都涉及到在狭小间隙中的流体粘性流动，描写这种物理现象的基本方程为Reynolds方程，它的普遍形式是： 这个椭圆型的偏微分方程仅仅对于特殊的间隙形状才可能求得解析解，而对于复杂的几何形状或工况条件下的润滑问题，无法用解析方法求得精确解。 数值法是将偏微分方程转化为代数方程组的变换方法。它的一般原则是：首先将求解域划分成有限个数的单元，并使每一个单元充分的微小，以至于可以认为在各单元内的未知量（例如油膜压力p）相等或者依照线性变化，而不会造成很大的误差。然后，通过物理分析或数学变换方法，将求解的偏微分方程写成离散形式，即是将它转化为一组线性代数方程。该代数方程组表示了各个单元的待求未知量与周围各单元未知量的关系。最后，根据Gauss消去法或者Gauss-Seidel迭代法求解代数方程组，从而求得整个求解域上的未知量。 用来求解Reynolds方程的数值方法很多，最常用的是有限差分法、有限元法和边界元法，这些方法都是将求解域划分成许多个单元，但是处理方法各不相同 。 在有限差分法和有限元法中，代替基本方程的函数在求解域内是近似的，但完全满足边界条件。 边界元方法所用的函数在求解域内完全满足基本方程，但是在边界上则近似地满足边界条件。 能量方程和弹性变形方程是流体润滑问题中考虑热效应和表面弹性变形时必须求解的重要方程，在本章中也将介绍它们的数值解法。 近年发展的多重网格法在润滑计算中开始得到应用，并有广泛的应用前景。本章的最后还介绍了多重网格法求解微分方程和积分方程。 9.1 Reynolds方程的数值解法numerical methods of Reynolds equation 9.1.1无量纲化 dimensionless 无量纲化也称为或归一化，是将有关方程的变量用合适的参考量(常数)进行转换，使变量的变化范围在1的数量级上。这个参考量称为相对单位。 为了判断影响润滑问题诸因素的影响大小，以便抓住主要问题和影响的主要方面，可采用归一化的处理方法，对方程组进行适当的简化。 另外，归一化还可以使所得的分析数据具有通用性和广泛性。 1. 方程的无量纲化 两个因素 微分方程中起作用的是变量的变化率，而不是其本身值的大小 结果的广泛性 1. 方程的无量纲化 固定瓦径向轴承 式中 不可压缩流体 流体动力粘度为常数 不可压缩流体 能量方程 温粘方程 2. 相似分析 1）相似条件 2）相似分析的基本步骤 算例 （2）选取无量纲单位并将方程无量纲化 （3）找出决定性和非决定性相似不变量 例如在稳定工况下流体静压润滑的油膜厚度h为常数，若不考虑相对滑动和热效应，则粘度?也是常数。这时Reynolds方程可简化为Laplace方程，即 （9.1） 将式（9.1）无量纲化，令x=XA，y=YB，A，B为几何尺寸；p=Ppr，pr为油腔压力；?=A2/B2；则无量纲化的Reynolds方程为 （9.2） 求解方程（9.2）的无量纲化边界条件为 （1）在油腔内 P=1 （2）在四周边缘上 P=0 若令 代入雷诺方程， 得无量纲化基本方程 (9.3) 这种形式的方程被称为Poisson方程。 再将P=Pc（1-Y2）代入方程（9.3），则变换成只含变量Pc和X的方程，即可求解中间断面上无量纲化压力Pc随X的变化。 即 （9.4） 广义雷诺方程方程的无量纲化 固定瓦径向轴承 式中 不可压缩流体 流体动力粘度为常数 不可压缩流体 能量方程 温粘方程 相似分析 相似条件 2）相似分析的基本步骤 算例 （2）选取无量纲单位并将方程无量纲化 （3）找出决定性和非决定性相似不变量 9.1.2 有限差分法finite difference method 在流体润滑计算中，有限差分方法的应用最为普遍。 现将有限差