In the study of large-scale localized strain features, localized strain can lead to energy release and seismic effects. From the macroscopic process of localized strain, it is found that localized deformation is accompanied by structural weakening, which is due to phase changes in rocks. The phase change corresponds to the secondary phase transition process in physics. The present study illustrates this phenomenon from the perspective of secondary phase transition theory in statistical physics, combining the localized strain features on a large scale. Theoretical analysis and experimental studies were carried out using three brittle rocks, including marble, granite and red sandstone. First, two perturbation methods, i.e., Krylov–Bogoliubov method and Poincare method, are used to calculate the higher-order control equations in analytical model. The influences of control equation coefficients on the strain localization process are analyzed. Then the uniaxial compression tests are carried out on the three rocks to record the strain process. Finally, the theoretical and experimental results are compared to analyze the strain localization phenomenon. The comparison results show that both the Krylov–Bogoliubov solution and the Poincare solution in the analytical model can well describe the evolution characteristics of localized strain. It indicates that the theoretical model is valid and has high accuracy. The Poincare method better simulates the nonlinear phenomenon of strain localization. The present study provides a new theoretical method to better understand the strain localization phenomenon.