Simulation for THz coherent undulator radiation from combination of velocity bunchings