A coupled hydro-thermo-mechanical model based on TLF-SPH for simulating crack propagation in fractured rock mass

In this paper, a hydro-thermo-mechanical coupling model based on the smoothed particle hydrodynamics with total Lagrangian formula (HTM-TLF-SPH) is proposed to simulate the crack propagation and instability process of fractured rock mass. TLF-SPH uses the Lagrangian kernel approximation, that is, the kernel function and its gradient need only be calculated once in the initial configuration, which is much more efficient than the smoothed particle hydrodynamics (SPH) based on the Euler kernel approximation. In TLF-SPH, particles interact with each other through virtual link, and the crack propagation path of rock mass is tracked dynamically by capturing the fracture of virtual link. Firstly, the accuracy and robustness of the HTM-TLF-SPH coupling model are verified by a reference example of drilling cold shock, and the simulation results agree well with the analytical solutions. Then, the crack propagation law of surrounding rock and the evolution characteristics of physical fields (displacement, seepage and temperature fields) after excavation and unloading of deep roadway under the coupling condition of hydro-thermo-mechanical are investigated. In addition, the seepage and heat transfer processes of the surrounding rock of Daqiang coal mine under different coupling conditions are successfully simulated. Meanwhile, the effect of the boundary water pressure difference on the temperature and seepage fields under the hydro-thermal coupling condition is quantitatively analyzed. A HTM-TLF-SPH model is proposed to simulate crack propagation of fractured rock massThe HTM-TLF-SPH model can effectively simulate fissure seepage and thermal diffusionTemperature diffusion has little effect on seepage in low permeable rock massThe effect of pore water on heat transfer rate mainly depends on the seepage directionThe pore water pressure at 2 meters in front of the roadway is most affected by temperature