From the theory of damage mechanics, based on the Hoek-Brown strength criterion and Weibull distribution law of rock microelement strength, a statistical constitutive model of rock thermal damage is established by using equivalent strain hypothesis, and the theoretical model is modified by considering the compression coefficient and residual strength correction coefficient. The rationality of the modified model is verified by experimental data. The results show that the stress-strain curves of rock can be divided into four stages: initial compaction, stable damage propagation, damage strengthening expansion, and damage failure according to the characteristics of rock damage evolution. The peak stress of rock increases exponentially with the increase of confining pressure, and the maximum damage evolution rate decreases exponentially with the increase of confining pressure, which indicates that confining pressure delays the development of cumulative damage. The peak stress and maximum damage evolution rate of rock decrease exponentially with the increase of temperature, which accelerates the damage of rock. The initial damage of rock is thermal damage caused by temperature, and the damage value increases with the increase of temperature. The revised theoretical curve reflects the characteristics of rock compaction stage and residual strength and improves the coincidence with the experimental curve. The research results provide a reference for the establishment of thermal damage constitutive model of rock in deep engineering.