提出了一种将有限元方法和接触理论相结合的内啮合齿轮副啮合刚度计算方法。该方法通过齿轮整体和局部有限元模型分离出啮合点宏观变形,利用线接触变形解析公式计算啮合点接触变形,求解非线性啮合平衡方程后得到齿轮副时变啮合刚度和载荷分布。该方法相比一般有限元法具有更高的计算效率和稳定性,同时克服了解析方法难以考虑斜齿轮和不同齿圈结构影响的缺点。最后,分析了内齿圈不同支撑数目和齿圈厚度对啮合刚度的影响。该方法可为内齿轮副强度及动态特性设计提供有效指导。