多重网格法求解雷诺方程

内燃机的主轴承 、连杆轴承承受着非常高的非稳定载荷 ，即动载荷，故 内燃机 轴承是典 型动载滑 动轴承。内燃机轴承在结构设计方面有很多特殊性，随着对轴承润滑理论及设计计算法研究的深入 ，尤其是计算机技术的发展和普及 ，使 内燃机 的设计计算方法更加完 善、更 加可 靠、更 加 系统。1986年，OsborneReynolds在-定假设条件下，由流体力学原理提出动压滑动轴承的润滑基本方程，总结出流体膜建立动压力的机理，为流体力学润滑理论奠定基矗内燃机轴承油膜压力计算、轴心轨迹计算、最小油膜厚度计算都是以雷诺方程为基础的。

对于 Reynolds方程的求解，现在主要的数值计算方法包括有限差分法和有限元法 。有限差分法主要针对常用的轴承形式 ；有 限元法用在轴瓦表面结构 复杂的情况，例如开有凶、沟槽、台阶等时。对于普通常用的轴承，采用有限差分法在很短的计算时间内就能得到比较满意的结果。有限差分法使用的是单纯的超松弛迭代法求解雷诺方程，而多重网格法是加速后的超松弛迭代方法 ，它 比单纯 的超松 弛迭代法具有更短的计算时间、更高的计算精度和更快的收敛速度。本文主要研究用多重网格法求解雷诺方程。

1 多重 网格法多重网格法的实质就是求解线性方程组的-种加速收敛的手段。多重网格法 的三大基本思想是：①细网格松弛，负责消除高频振荡的误差；②粗网格校正，负责低频光滑误差；③套迭代技术，负责通过限制和延拓算子连接所有层共同求解同-个问题。多重网格法就是在不同的层上进行求解，所有层相互协调地求解同-个问题，从而使高、低频偏差分量都能很快地消除和最大限度地提高迭代效率 。由于限制和拖延的次数不同，多重网格法可分为如图 1所示的几种循环方法。

肼 - JI，-1- -JI，-2- - - 3- - Jl，-4- .--I------ -- ]-二(a)两重 V”循环 (b)多重 V”循环 (c)完全多重网格法图 1 多重网格的不同循环方法下面以两重网格V”循环为例介绍多重 网格法 的计算方法。取粗网格和细网格步长，且粗网格步长是细网格步长的 2倍。具体步骤如下：①设定初始值；②计算细网格上的亏损量 ；③从 细网格到粗 网格转移 亏损量 ；④在粗网格上精确求解修正量 ；⑤ 由粗网格到细网格转移修正量；⑥计算细网格修正后的量。多次重复步骤①～⑥ ，直至结果收敛 。

2 求解二维雷诺方程加速性验证本文通过具体实例 ，验证 多重 网格法是-种高效的加速收敛的手段。实例是基于有限差分法求解滑动轴承的雷诺方程得到油膜 压力 ，分别采用单纯的超松弛迭代法和多重网格法，最后再与使用了迭代精度要求的计算结果进行对比。

计算 参 数 如 下：径 向 滑 动 轴 承 的 内 径 D-2O mm，宽径比 B/D-1，偏心率 E-0.6，半径间隙f-0.145 5 mm，轴颈转速 n-500 r/min，润滑油在工作状态下的动力黏度 叩-0.028 Pa·s，相对间隙 -0.001 1；网格密度取 20×20；总迭代次数为 2O次。超松弛迭代法的计算过程稳定，所以第-种计算过程和多重网格法的各层计算都采用了超松弛迭代法，且取松弛收稿 El期：2012-11-13；修回日期：2012-11-20作者简介：许胜丽 (1986-)，女，河南开封人，在谈硕i：研究生，主要从事机械设计及理论的学习研究。

2013年第 2期 许胜丽，等：多重网格法求解雷诺方程 · 71 ·因子为 1.3。油膜压力的坐标原点放在轴承宽度中央。

(1)以单纯的超松弛迭代法求解雷诺方程，迭代 2O次，得到轴承宽度中央的周向油膜压力，如图 2所示。

周向角度/rad图 2 单纯超松 弛迭代 法(2)以多重网格法求解雷诺方程，得到轴承宽度中央的周 向油膜压力。计算过程是采用两重网格V”循环法，设定第-重网格迭代 88次，第二重网格迭代 4次，即在第-重细网格上迭代 8次，转移到第二重粗网格迭代 4次，再转回到第-重细网格迭代 8次，总共迭代 2O次。

第-重网格迭代 8次结果如图 3所示，第二重网格迭代 4次，迭代次数 k的值依次取 1、2、3、4，且油膜压力P的初值采用第-重网格的计算结果和误差乘以相应的完全加权限制操作数。第二重网格结点(除了网格边界点)均由第-重网格的8个点包围，故设定限制操作数由第-重的每个点的 1/8相加得到。应用限制操作数可把第-重的计算结果转移作为第二重计算的初值。

O.1O7 L. i. .7- ～ / 0 0.5 1.U 1.5 2.U 2.5 3.0 3.5周向角度/tad图 3 第-层 网格求解 结果有了第二重 网格求解结果 ，再转 回第-重网格迭代 8次，油膜压力 P的初值采取由第二重网格的计算结果乘以对应的延拓操作数得到。同样限制操作数可把第二重网格的计算结果转移作为第二次在第-重网格计算的初值。最终求解结果如图4所示。

(3)采用迭代精度要求为0.01的相对误差来限制迭代次数，得到轴承宽度中央的周向油膜压力。此次计算比前面两种计算在精度上要高，用时也长，更接近精确值，所以它可以作为前面两种计算方法相对比的标准。

增加收敛准则，用迭代精度要求来限制迭代结果，保证相对误差为 0.0l，求解结果如图 5所示。

由图2、图 4和图5可知，将迭代同样的次数分别用单纯的超松弛迭代法和用多重网格法求得的结果与使用迭代精度要求 0.01相对误差的结果比较，可知多重网格法的迭代结果更接近精确值，说明多重网格法有更高的迭代效率。在-般求解径向滑动轴承的雷诺方程时，都要求达到-定精度，采用普通的计算方法，可能要迭代的次数非常多，同时仪常耗时；而用多重网格法能在很少次数的迭代中达到很高的迭代精度，减少了计算量、节省了计算时间。

f -广 l/7/ 周向角度/rad图 4 第三层网格求解结果f i. . . .f..z.，。

/ 1- I /-/ i i /周 向角度 /rad图 5 相对误差为 0.01的迭代结果3 结束语在求解二维雷诺方程中，使用常规的单网格，如果要达到很高的精度，将需要迭代相当多的次数，整个计算过程也将耗费很多时间，但是采用多重网格法在各层网格上进行较少的迭代，便可以得到接近精确解的解，所以多重网格法求解雷诺方程将可以比较快地达到所要求的精度，是-种比较高效的计算方法。