A mesostructure-based damage model for thermal cracking analysis and application in granite at elevated temperatures

Thermal stress within rock subjected to thermal load is induced due to the different expansion rates of mineral grains, resulting in the initiation of new intergranular cracking and failure at elevated temperatures. The heterogeneity resulting from each constituent of rock should be taken into account in the study of rock thermal cracking, which may aid the better understanding of the thermal cracking mechanisms in rock. In this paper, a mesostructure-based numerical model for the analysis of rock thermal cracking is proposed on the basis of elastic damage mechanics and thermal–elastic theory. In the proposed model, digital image processing (DIP) techniques are employed to characterize the morphology of the minerals in the actual rock structure to build a numerical specimen for the rock. In addition, the damage accumulation induced by thermal (T) and mechanical (M) loads is considered to modify the elastic modulus, strength and thermal properties of individual elements with the intensity of damage. The proposed model is implemented in the well-established rock failure process analysis (RFPA) code, and a DIP based RFPA for the analysis of thermally induced stress
and cracking of rock (abbreviated as RFPA-DTM) is developed. The model is then validated by comparing the simulated results with the well-known analytical solutions. Finally, taking an image from a granite specimen as an example, the proposed model is used to study the thermal cracking process of the granite at elevated temperatures and the effects of temperature on the physical–mechanical behaviors of the granite are discussed. It is found that thermal cracks mostly initiate at the location of mineral grain boundaries and propagate along them to form locally closed polygons at the elevated temperatures. Moreover, the effects of temperature on the uniaxial compressive strength and elastic modulus of the granite are quite different. The uniaxial compressive strength decreases consistently with increasing temperature, but there exists a threshold temperature for elastic modulus which starts to decrease as the temperature increases after it exceeds the threshold.