The aim of this paper is to present an integrated procedure for the optimization of dimensions of the ITI dental implant with dynamic chewing loads. The procedure is composed of uniform design of experiments, explicit dynamics finite element analysis, Kriging interpolation, and Nelder-Mead simplex local search optimization method. Firstly, uniform design method is employed to create a set of experiments. Then, explicit dynamics finite element analysis software ANSYS/LS-DYNA is used to analyze the micromotion of cortical and cancellous bones while the dynamic chewing loads acts on the implant. Next, Kriging interpolation is applied to construct the surrogate model of micromotion based on the input and output data of experiments. Finally, Nelder-Mead simplex local search method is applied to find the optimal solution of dimensions of dental implant under the goal of minimizing the micromotion. After performing the optimization procedure presented in this paper, the micromotion of the ITI dental implant system model can successfully be reduced by a rate of 29.7 %.