Design analysis and computer experiments (DACE) model is widely used to express efficiently nonlinear responses in the field of engineering design. As a DACE model, kriging model can approximately replace a simulation model that is very expensive or highly nonlinear. The kriging model is composed of the summation of a global model and a local model representing deviation from the global model. The local model is determined by correlation coefficient with the pre-sampled points, where the accuracy and robustness of the kriging model depends on the selection of proper correlation coefficients. Therefore, to achieve the robust kriging model, the range of the correlation coefficients is explored with respect to the degrees of the global model. Based on this study we propose the proper orders of the global model and range of parameters to make accurate and robust kriging model.