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

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