The spindle is the main component in machine tools. To develop high speed machine tools, a lot of studies have been carried out for high speed spindle. Bearing is very important part in spindle. The bearing clearance is influenced by shrink fit and thermal expansion during operation. The designer must take into account the reduction of shrink fits. The aims of this study are to grasp the shrink fits and behavior of a bearing which is a deeply connected with fatigue life of bearing and performance of spindle through FEM(Finite Element Method). This paper proposed optimum value of shrink fit considering deformation of spindle and stress of fitting area using design of experiments. Thus, the proposed formula can be used to obtain bearing internal clearance.