A phase field model is developed to simulate the grain evolution of 17-4PH steel during cyclic heat treatment (CHT). Our simulations successfully reproduce the grain morphologies of every CHT. In the process of every CHT, phase transformation recrystallization happens. The recrystallized grains appear mainly on the original grain boundaries. The average grain size of 13.2 μm obtained by 1040 ℃×1 h solution treatment for this experimental steel can be refined to 2.2 μm after five CHT's. Furthermore, the effects of phenomenological parameters in our model are discussed.