Numerical Modeling of Temperature Effect on Tensile Strength of Granitic Rock

The aim of this paper is to numerically predict the temperature effect on the tensile strength of granitic rock. To this end, a numerical approach based on the embedded discontinuity finite elements is developed. The underlying thermo-mechanical problem is solved with a staggered method marching explicitly in time while using extreme mass scaling, allowed by the quasi-static nature of the slow heating of a rock sample to a uniform target temperature, to increase the critical time step. Linear triangle elements are used to implement the embedded discontinuity kinematics with two intersecting cracks in a single element. It is assumed that the quartz mineral, with its strong and anomalous temperature dependence upon approaching the alpha-beta transition at the Curie point (similar to 573 degrees C), in granitic rock is the major factor resulting in thermal cracking and the consequent degradation of tensile strength. Accordingly, only the thermal expansion coefficient of quartz depends on temperature in the present approach. Moreover, numerically, the rock is taken as isotropic except for the tensile strength, which is unique for each mineral in a rock. In the numerical simulations mimicking the experimental setup on granitic numerical rock samples consisting of quartz, feldspar and biotite minerals, the sample is first heated slowly to a target temperature below the Curie point. Then, a uniaxial tension test is numerically performed on the cooled down sample. The simulations demonstrate the validity of the proposed approach as the experimental deterioration of the tensile strength of the rock is predicted with agreeable accuracy.

Automated summary

This paper aims to predict the temperature effect on the tensile strength of granitic rock using a numerical approach. It is found that the temperature dependence of the quartz mineral near 573 degrees Celsius is the main factor leading to thermal cracking and degradation of tensile strength. Numerical simulations demonstrate that the proposed approach can accurately predict the experimental deterioration of the tensile strength of the rock.