The invention discloses a gear contact fatigue reliability analysis method. Firstly, the contact stress model of gear is established, the residual fatigue strength after a certain cycle is obtained onthe basis of a fatigue S-N curve, and then the function of gear contact fatigue failure is constructed. Then, according to the distribution characteristics of the random parameters, the Monte Carlo method is used to generate point set, the initial Kriging model is constructed by randomly selecting N points as the experimental design points, and then the sample points which are closest to the limit state surface or have the largest error are obtained by active learning strategy, and the sample points are added to the original experimental design points as the best sample points, and the Kriging model is updated until convergence. According to the convergent Kriging model, the approximate value of the function corresponding to all the sample points is obtained, and the failure probability is calculated by Monte Carlo method. Finally, an example is given to verify the feasibility of the proposed method. The gear contact fatigue reliability analysis method of the invention is more efficient and practical.