The novel coupling recrystallization model is proposed in this study. First, the deformation microstructure was simulated by the finite element method based on the strain gradient crystal plasticity theory. The calculated dislocation density and crystal orientation were transferred to the recrystallization phase-field simulation. The initial subgrain structures used in phase-field simulation were determined by a relationship between dislocation density and subgrain size with the dislocation density distribution calculated by crystal plasticity simulation. The so-called KWC phase-field model, which can introduce both subgrain rotation and grain boundary migration, was employed, and spontaneous nucleation and grain growth were simulated simultaneously.