This paper investigates the estimation of high-dimensional factor models in which factor loadings undergo an unknown number of structural changes over time. Given that a model with multiple changes in factor loadings can be observationally indistinguishable from one with constant loadings but varying factor variances, this reduces the high-dimensional structural change problem to a lower-dimensional one. Due to the presence of multiple breakpoints, the factor space may expand, potentially causing the pseudo factor covariance matrix within some regimes to be singular. We define two types of breakpoints: {\bf a singular change}, where the number of factors in the combined regime exceeds the minimum number of factors in the two separate regimes, and {\bf a rotational change}, where the number of factors in the combined regime equals that in each separate regime. Under a singular change, we derive the properties of the small eigenvalues and establish the consistency of the QML estimators. Under a rotational change, unlike in the single-breakpoint case, the pseudo factor covariance matrix within each regime can be either full rank or singular, yet the QML estimation error for the breakpoints remains stably bounded. We further propose an information criterion (IC) to estimate the number of breakpoints and show that, with probability approaching one, it accurately identifies the true number of structural changes. Monte Carlo simulations confirm strong finite-sample performance. Finally, we apply our method to the FRED-MD dataset, identifying five structural breaks in factor loadings between 1959 and 2024.