一种考虑岩体热力耦合问题的连续-非连续数值模拟方法

本发明公开了一种考虑岩体热力耦合问题的连续‑非连续数值模拟方法，将模拟区域划分为潜在破坏区和不可破坏区，并对模拟区域进行网格划分；然后生成用于模拟计算的热计算节点系统和力计算节点系统，两者间存在映射关系；在潜在破坏区网格边界设置节理单元，通过节理单元的破坏模拟岩体破裂过程；模拟过程分解为若干计算步，在每一个计算步内热传导过程和力学计算过程是依次求解的，依据节理单元的变形和应力状态判断其是否破坏，并更新映射关系，直至最后一个计算步结束，完成岩体的热破裂过程模拟。本发明能够避免在热传导模拟过程中引入人工传热参数，有针对性地研究岩体局部区域的破裂演变规律，降低模型的计算成本，并有效地模拟岩体热力耦合问题中的动态破裂过程。