简单圆柱滑动轴承油膜压力计算Matlab编程说明.pdfVIP

大卫贝克汉公 （著） 2018/8/1, NJ, SEU 典型圆柱滑动轴承油膜压力计算的 Matlab 程序编制说明 滑动轴承 图 1 圆柱型滑动轴承示意图 编程思路 步骤1：计算膜厚H 的二维分布(,) ； 步骤2：计算每一列 （等于某一角度）上的差分法计算式中的系数A, B, C, D, E, F; 步骤2：计算压力P 分布，同时导入Reynolds 边界条件，设置好迭代公式。 编程中关键问题 显然油腔是环形的，周向连续，当我们把他切开，铺平，近似为矩形，这里 我们设置方向的网格长度为1° ，网格数则为360 格，方向上分成 114 格，这 里考虑到对称性，只计算方向上的一半57 格，以节省计算量，环形网格切开铺 平后如图2. 大卫贝克汉公 （著） 2018/8/1, NJ, SEU 图 2 差分网格展开示意图 在油膜起始边切开，平铺后，则第1 列和第361 列是重合的，第58 行为 向 （轴向）上的对称线。 关键问题1：差分网格最后一行、一列的计算问题 用一个 P 矩阵存储压力初值。事实上我们只需要计算图2 中红色虚线内的 内节点，即57×359 矩阵，即，以外节点均为边界，压力P=0 ，但是我们计算一 个节点的P 值，需要用到相邻的四个节点的P 值，我们从第 （2, 2 ）节点分别向 右和向下计算是，当计算到第360 列和58 行，遇到麻烦，因为往右和往下，没 有后续节点压力值用了，所以这时编程计算到最后一列，最后一行时，需要单独 处理，破坏了编程的连续性，是程序变得稍微负责一点。这是有一些技巧可以运 用，I 方向上我们多存储一列 P 值，即将第361 列也加进去，361 列其实就是第 一列，但是由于程序是递增向前推进，计算第360 列时，不能掉头摘取第一列值 用，加入第361 列后，第360 行右边节点P 值就有了，所以保证了i 方向上程序 推进的连续性。 j 方向上的第58 行的计算却不能采用上述方法，因为第361 列时边界，P 值 在每次迭代计算中始终为0 ，而第58 行下面其实有第59 行，他的P 值等第57 行的P 值得，第57 行的P 值是随着每一次的迭代计算为变化的，所以如果在第 58 行下面添加一行，但是其P 值不好确定，所以编制程序时，不如直接先算2→57 行 ， 58 行 单 独 计 算 。 计 算 程 序 如 图 3 ： 大卫贝克汉公 （著） 2018/8/1, NJ, SEU 图 3 Matlab 部分程序 关键问题2：导入Reynolds 边界条件。 Reynolds 边界条件导入常用的有效且简便的方法是：每行均从油膜起始边向 终止边计算，当计算出某点的压力为负值时，即取为0 ，此点即为油膜自然破裂 边的近似位置，后面所有点的压力值均取为0 。编制程序时可这样处理，如图三， 在计算每一行点P 值之前，先进行判断，P 值是否≤ ，若