目前,齿轮时变啮合刚度计算大多参考ISO实验值进行公式简化,该方法不能准确反映不同齿数齿轮啮合刚度的综合特性。针对此问题,利用APDL语言建立渐开线直齿轮的参数化有限元模型,通过接触分析得到轮齿接触位移值并计算刚度。该方法考虑了轮缘厚度对刚度的影响;通过最小二乘法得到直齿轮单齿刚度公式,并依此推导得到任意齿数啮合的时变啮合刚度公式;并与石川公式、势能法、Kuang方法对比,说明了各计算方法在计算时变刚度时的优缺点。将FEA公式与ISO计算单齿最大刚度的结果进行对比,验证了该刚度公式的可靠性。