本发明提供了一种粗糙表面法向接触刚度计算方法，包括：基于随机数值矩阵和粗糙表面自相关特性生成一个满足Gauss分布的粗糙表面；基于粗糙表面生成三维粗糙表面几何模型；对三维粗糙表面几何模型进行网格划分得到网格化的粗糙表面模型；获取约束边界条件和载荷边界条件并进行求解计算，得到粗糙表面上不同位移下的平均接触压力和平均变形量；基于不同位移下的平均接触压力和平均变形量，拟合生成接触面法向的平均变形量与平均接触压力的关系；基于接触面法向的平均变形量与平均接触压力的关系，得到粗糙表面的法向接触刚度。该方法能够考虑安装边接触面粗糙度的影响，定量计算出粗糙表面法向接触刚度随接触压力的变化，计算结果准确可靠。