In crack growth life, uncertainties are caused by variance of geometry, applied loads and material properties. Therefore, the reliability estimation for these uncertainties is required to keep the robustness of calculated life. The stress intensity factors are the most important variable in crack growth life calculation, but its equation is hard to know for complex geometry, therefore they are processed by the finite element analysis which takes long time. In this paper, the response surface is considered to increase efficiency of the reliability analysis for crack growth life of a turbine wheel. The approximation model of the stress intensity factors is obtained by the regression analysis for FEA data and the response surface of crack growth life is generated for selected factors. The reliability analysis is operated by the Monte Carlo Simulation for the response surface. The results indicate that the response surface could reduce computations that need for reliability analysis for the turbine wheel, which is hard to derive stress intensity factor equation, successfully.