The crystallization, corresponding to the fcc structure (with packing density p ≈ 0.74), of smooth equal hard spheres under batch-wised feeding and three-dimensional interval vibration is numerically obtained by using the discrete element method. The numerical experiment shows that the ordered packing can be realized by proper control of the dynamic parameters such as batch of each feeding § and vibration amplitude A. The radial distribution function and force network are used to characterize the ordered structure. The defect formed during vibrated packing is characterized as well The results in our work fill the gap of getting packing density between random close packing and fcc packing in phase diagram which provides an effective way of theoretically investigating the complex process and mechanism of hard sphere crystallization and its dynamics.