在滚珠丝杠副运动过程中,滚珠与滚道之间的摩擦和磨损导致滚珠丝杠副精度衰退,降低其精度保持性.针对高速滚珠丝杠副螺旋线滚道接触区域的非对称性,采用传统的赫兹接触理论对丝杠和螺母滚道接触几何形状和关系的计算存在较大误差,提出了一种基于最小余能原理对高速滚珠丝杠副滚珠与滚道之间的法向接触应力进行分析计算的方法.首先,基于Frenet-Serret坐标转换公式建立滚珠丝杠副坐标系,对滚珠与滚道接触点处的接触情况进行几何描述,并通过2维插值算法重新进行网格划分以便于后续的接触力学分析;对滚珠与滚道之间的接触区域和控制方程进行离散化分析,提高数值计算的精确性;为了解决接触问题计算求解过程复杂、收敛速度缓慢,通过运用变分原理将接触问题转成极值问题,采用共扼梯度法进行循环迭代来实现快速收敛求解,提高接触数值解法的计算效率,求解出滚珠与滚道之间的接触应力分布;由于接触表面变形的计算相当于求影响系数矩阵与法向压应力之间的卷积,2维卷积运算和对应的2维快速傅里叶变换被用于求解接触问题的数值解.在数值求解中,主要的计算量集中在弹性变形的计算上,使用2维快速傅里叶变换(FFT)技术方便快速地计算接触区域弹性变形.将最小余能原理与赫兹接触理论分别用于求解光滑的球与平面之间的弹性接触应力分布,将两种计算方法得出的结果进行对比,计算数值重合度很好,验证了基于最小余能原理求解接触应力分布的正确性;基于最小余能原理分析螺旋升角对滚珠与滚道接触点处的接触面积和弹性变形的影响,并将该方法所得结果与赫兹接触理论计算结果进行对比分析,发现赫兹接触解与最小余能非赫兹接触精确解相比,随着滚珠丝杠副滚道螺旋升角的增大,丝杠与滚珠接触点A处的误差始终大于螺母与滚珠接触点B处的误差,丝杠滚道与滚珠接触点A处非赫兹解接触应力峰值逐渐变小,螺母滚道与滚珠接触点B处非赫兹解接触应力峰值逐渐增大,且接触点A处应力峰值始终大于接触点B处峰值,计算所得滚珠与滚道接触区域的面积和长短轴有所不同.最小余能原理非赫兹接触计算方法可以准确计算滚珠和滚道接触区域应力分布,并且可以全面和准确的计算滚道磨损带的宽度和深度.研究结果表明,采用最小余能非赫兹接触解可以有效提高滚珠丝杠副接触应力分布计算精度,保证滚珠丝杠副的精度保持性.