A coupled thermo-mechanical model based on the combined finite-discrete element method for simulating thermal cracking of rock

Abstract Based on a combined finite-discrete element method (FDEM), this study builds a coupled thermo-mechanical model (termed FDEM-TM) to simulate thermal cracking of rock. The coupled thermo-mechanical model consists of two major parts. In the first part the temperature distribution of the system is analyzed based on the heat conduction equation. In the second part the thermal stress caused by temperature change is added to perform mechanical fracture calculation. Based on these two parts, we can model rock fracture driven by thermo-mechanical coupling. Three examples with analytic solutions are used to verify the correctness of the model in dealing with the problems of steady-state heat conduction, unsteady-state conduction, and thermal-mechanic coupling, respectively. In addition, an example of thermal cracking is also given and compared with the experimental results. The simulation results are in excellent agreement with the analytical solutions or experimental results, verifying the correctness of the coupled thermo-mechanical model to simulate thermal cracking. The proposed method provides a new tool for thermal-mechanical coupling problems in geothermal exploitation.

[1]  O. K. Mahabadi,et al.  Three-dimensional FDEM numerical simulation of failure processes observed in Opalinus Clay laboratory samples , 2014 .

[2]  Hong Zheng,et al.  FDEM-flow3D: A 3D hydro-mechanical coupled model considering the pore seepage of rock matrix for simulating three-dimensional hydraulic fracturing , 2017 .

[3]  Hany F. Abdalla,et al.  Concrete cover requirements for FRP reinforced members in hot climates , 2006 .

[4]  D. Owen,et al.  A combined finite‐discrete element method in transient dynamics of fracturing solids , 1995 .

[5]  Hong Zheng,et al.  Combined Finite-Discrete Element Method for Simulation of Hydraulic Fracturing , 2016, Rock Mechanics and Rock Engineering.

[6]  Y. T. Feng,et al.  Discrete thermal element modelling of heat conduction in particle systems: Basic formulations , 2008, J. Comput. Phys..

[7]  David R. Owen,et al.  Discrete thermal element modelling of heat conduction in particle systems: Pipe-network model and transient analysis , 2009 .

[8]  A. Munjiza The Combined Finite-Discrete Element Method , 2004 .

[9]  R. P. Young,et al.  Bonded-particle modeling of thermally fractured granite , 2008 .

[10]  H. Zheng,et al.  Three-Dimensional Hydromechanical Model of Hydraulic Fracturing with Arbitrarily Discrete Fracture Networks using Finite-Discrete Element Method , 2017 .

[11]  Jianchun Li,et al.  A coupled thermo-mechanical discontinuum model for simulating rock cracking induced by temperature stresses , 2015 .

[12]  Watson L. Vargas-Escobar Discrete Modeling of Heat Conduction in Granular Media , 2002 .

[13]  S. Bradshaw,et al.  Bonded-particle modelling of microwave-induced damage in ore particles , 2010 .

[14]  Tang Chun Investigation on coupled thermo-hydro-mechanical-damage(THMD) model for rock failure process , 2008 .

[15]  B. Hobbs,et al.  Particle simulation of thermally-induced rock damage with consideration of temperature-dependent elastic modulus and strength , 2014 .

[16]  S. Bradshaw,et al.  Confined particle bed breakage of microwave treated and untreated ores , 2011 .

[17]  Tang Chun NUMERICAL INVESTIGATION ON ROCK FAILURE PROCESS INDUCED BY THERMAL STRESS , 2006 .

[18]  Giovanni Grasselli,et al.  Continuum–discontinuum analysis of failure mechanisms around unsupported circular excavations in anisotropic clay shales , 2014 .

[19]  Y. Mahmutoğlu Mechanical Behaviour of Cyclically Heated Fine Grained Rock , 1998 .

[20]  Hong Zheng,et al.  A two-dimensional coupled hydro-mechanical finite-discrete model considering porous media flow for simulating hydraulic fracturing , 2016 .

[21]  M. Xia,et al.  Thermo-mechanical coupled particle model for rock , 2015 .

[22]  A. Ibrahimbegovic,et al.  Softening behavior of quasi-brittle material under full thermo-mechanical coupling condition: Theoretical formulation and finite element implementation , 2014 .

[23]  Jiansheng Xiang,et al.  Implementation of an Empirical Joint Constitutive Model into Finite-Discrete Element Analysis of the Geomechanical Behaviour of Fractured Rocks , 2016, Rock Mechanics and Rock Engineering.

[24]  Daniel Moos,et al.  Numerical modeling of injection, stress and permeability enhancement during shear stimulation at the Desert Peak Enhanced Geothermal System , 2015 .

[25]  B. S. A. Tatone,et al.  A calibration procedure for two-dimensional laboratory-scale hybrid finite–discrete element simulations , 2015 .

[26]  O. K. Mahabadi,et al.  Numerical simulation of acoustic emission in brittle rocks by two-dimensional finite-discrete element analysis , 2013 .

[27]  W. Cho,et al.  The influence of an excavation damaged zone on the thermal-mechanical and hydro-mechanical behaviors of an underground excavation , 2008 .

[28]  A. Munjiza,et al.  Large Strain Finite Element Method: A Practical Course , 2015 .

[29]  Antonio Munjiza,et al.  A generalized anisotropic deformation formulation for geomaterials , 2016, 1805.06024.

[30]  John-Paul Latham,et al.  Challenges of a coupled combined finite-discrete element approach to explosive induced rock fragmentation , 1999 .

[31]  Antonio Munjiza,et al.  Validation of a three-dimensional Finite-Discrete Element Method using experimental results of the Split Hopkinson Pressure Bar test , 2014 .

[32]  Antonio Munjiza,et al.  Combined single and smeared crack model in combined finite-discrete element analysis , 1999 .