本发明公开了一种齿轮接触疲劳可靠性分析方法。该方法首先建立了齿轮的接触应力模型，并基于疲劳S‑N曲线获取了一定周期后的疲劳剩余强度，进而构建了齿轮接触疲劳失效的功能函数。然后，根据随机参数的分布特性，运用蒙特卡洛方法产生点集，并从中随机选出N个点作为实验设计点构建了初始Kriging模型，再通过主动学习策略，获得了最接近极限状态面或误差最大的样本点，并将其作为最佳样本点增加到原有实验设计点中，更新了Kriging模型，直至收敛。根据收敛的Kriging模型，获取了所有样本点集对应的功能函数近似值，并采用蒙特卡洛方法计算了其失效概率。最后，通过一具体算例验证了该方法的可行性。该发明的齿轮接触疲劳可靠性分析方法更具高效性和实用性。