Experimental and Numerical Simulation Investigation of the Mechanical Properties of Anhydrite Rock Under Freeze–Thaw Environment

To study the mechanical and cracking modes of anhydrite rock under the freeze–thaw weathering process, the physico-mechanical characteristics and morphology evolutions of anhydrite samples were determined by a series of laboratory tests. Then, a numerical simulation model was established through the PFC^2D program, and the types and number of cracks during the uniaxial compression conditions were analyzed. Finally, the distribution of maximum principal stress and shear stress was revealed. The results indicate that as the number of freeze–thaw cycles increases, there is a growth in the mass loss rate and macroscopic damage variables while the uniaxial compression strength and elastic modulus decrease exponentially. Under uniaxial compression stress, the proportion of tensile cracks in the anhydrite model is the highest, followed by tensile shear cracks and compressive shear cracks. As the number of freeze–thaw cycles increases, the proportion of tensile cracks increases exponentially, while the proportion of tensile shear cracks and compressive shear cracks decreases exponentially. Furtherly, it is found that the maximum principal stress and maximum shear stress extreme values decrease exponentially with the increase of freeze–thaw cycles. For example, after 120 cycles, the maximum shear stress at the peak stress point decreased by 47.3%. The research results will promote the comprehension of anhydrite rock geotechnical engineering disaster mechanisms in cold regions.