工程尺度岩体破裂全过程模拟的PD-FEM数值计算方法及系统

摘要

本发明提供了一种工程尺度岩体破裂全过程模拟的PD‑FEM数值计算方法及系统，根据是否发生损伤，将岩体计算模型划分为近场动力学区、有限元区和耦合区；在耦合区中进行力与位移的传递，将上覆岩层重力及构造应力转化为有限元区域力边界条件，约束转化为位移边界条件，计算出有限元区的结点位移，将此作为近场动力学区的边界条件；在上述边界条件下，迭代计算物质点位移以及损伤，求解岩体的变形特征，实现岩体破裂的工程尺度模拟。本发明能够有效减少工作量，保证计算精度的同时大幅度提高计算效率，对于解决岩体破裂工程尺度问题具有重要意义。