Analytical Solution of Ice-Rock-Model Stress Field and Stress Intensity Factors in Inhomogeneous Media

The stress distribution and fracture parameter calibration of ice-rock models are important aspects of studying rock properties at high altitudes and latitudes. However, progress in ice-rock modeling has been slow and singular, and it is limited due to the discrete nature of rocks and the applicability of fracture mechanics. In this study, a circular inhomogeneous ice-rock model is proposed for the first time, and a method is provided for calculating the stress field of the model under biaxial loading. A method for calculating the single-crack stress intensity factor of the model subjected to biaxial compressive loading is also provided. The novelty of this work is that the inhomogeneous ice-rock model is treated as a superposition of two models, namely, a circular pore plate and circular ice, according to the superposition principle. The key is that the stress field distribution law of the ice-rock model is obtained based on the basis of the displacement continuity of the ice-rock interface. The analytical and approximate solutions of the stress intensity factor of a single crack were also obtained by considering the normal phase effect of the crack surface and combining the stress distribution law of the ice-rock model. Comparison with the CAE method was made to verify the correctness of the stress field and stress intensity factor calculation methods. The evolution laws of lateral pressure coefficients, the elastic modulus ratio of ice and rock on the stress field, and the stress intensity factor were analyzed. The effects of lateral pressure coefficients, elastic modulus ratios, and crack distributions on the failure modes were investigated using the extended finite element method (XFEM). This study can provide a theoretical basis for the evaluation of mechanical properties and prediction of the failure modes of frozen rock bodies.