弹流润滑数值计算方法.docx

这是流场的求解的方法,求解点接触弹流润滑的迭代子程序ITER中采用Jacobi迭代方法与线接触差别较大,具体过程如下: (1)迭代求解是沿着X轴(流速)方向进行的,当对该行的节点迭代完毕后,在Y轴方向前进一列; (2)由于Jacobi迭代时前后两个节点上分别加、减相同的压力增量,见式(4-6),因此,同一行的压力修正是相关联的。因此,每一行的压力修正不能单独完成,而是要把全行的压力增量作为相关联的未知量进行求解。在子程序ITER中,从标号20下两行到标号50的语句是用来建立这一联立方程的。由于该联立方程是三对角的,为了节约空间,方程的系数和右端等均放在缩减的矩阵A中。 (3)然后通过调用