本文用松弛法直接迭代求解弹性流体动力润滑(EHL)椭圆接触问题的一组非线性微分方程,获得了具有典型EHL特征的完全数值解。为提高收敛精度,采用了变松弛因子,并且自行推导了较高阶的有限差分公式。在计算弹性变形时,对二维的压力分布作双二次插值逼近,并通过极坐标变换克服了积分的奇异性。 In this paper, a set of nonlinear differential equations for elastohydrodynamic lubrication (EHL) elliptic contact problems are directly iteratively solved by relaxation method, and a complete numerical solution with typical EHL features is obtained. In order to improve the accuracy of convergence, a relaxation factor was adopted and a higher order finite difference formula was derived. In the calculation of elastic deformation, the two-dimensional pressure distribution for bi-quadratic interpolation approximation, and overcome the integral singularity by polar transformation.