针对角接触球轴承拟静力学模型非线性方程组数值求解过程中的不收敛问题,在现有求解方法的基础之上,提出了一种改进的数值求解法。该方法基于迭代变量之间的几何约束关系,根据几何关系的物理含义,由程序自动计算给出合理的迭代初值,避免了人为给定初值的盲目性,加速了拟静力学方程组的收敛性,减少了编程工作量。与现有方法对比计算,验证了改进算法的有效性。