We introduce an extension of the density functional perturbation theory (DFPT) that allows self-consistent linear-response calculations from a DFT + U ground state. Using this scheme, the full phonon dispersion of strongly correlated materials, whose ground state can be captured with Hubbard-corrected functionals, can be accessed with unprecedented accuracy and numerical efficiency. The tool is applied to the study of MnO and NiO in their antiferromagnetic (AFII) ground state. Our results confirm the highly noncubic behavior of these systems and show a strong interplay between features of the phonon spectrum and the occupation of specific d states, suggesting the possibility to investigate the electronic structure of these materials through the analysis of their phonon spectrum.