利用物质点法对原位应力下爆炸引起的岩石破坏进行的三维数值研究

《Computers and Geotechnics》：Three-dimensional numerical investigation of blast-induced rock failure under in-situ stress using the material point method

爆炸动力学与岩石破坏机制研究

　　
深地下岩体爆破破坏机理的三维数值模拟研究

研究背景与工程需求
随着矿产资源开发向深部转移，地下工程爆破作业面临日益复杂的地应力环境。深部岩体承受着高达数十兆帕的原始应力场，这种特殊工况显著改变了传统浅层爆破的破坏模式。现有研究表明，地应力水平与岩体破碎特征存在非线性关系，其中最大主应力方向直接影响裂隙扩展路径。然而，传统二维数值模拟难以准确捕捉三维应力场中爆炸波传播的各向异性特性，特别是裂隙网络在复杂应力状态下的空间重构规律。

数值方法创新与模型构建
本研究团队提出基于GIMP-CPDI耦合的三维动态响应分析框架，通过结合广义插值与流固耦合粒子追踪技术，实现了岩体损伤演化过程的高精度模拟。该方法突破传统有限单元法网格畸变限制，同时避免纯欧拉方法对界面追踪的不足。特别开发的非反射边界条件有效消除了数值域边界处的波反射干扰，这对准确模拟应力波在自由表面的反射-折射行为至关重要。

关键技术突破体现在三个方面：首先，建立了考虑地应力梯度影响的初始应力场生成算法，通过有限元-粒子混合方法实现三维应力张量的精确赋值；其次，开发了基于主应力方向的动态损伤演化模型，将传统材料破坏准则拓展为四维空间（应力历史、加载速率、空间坐标和时间）的耦合作用函数；最后，创新性地将流固耦合方程与损伤演化方程进行解耦处理，确保在爆炸冲击波超高速传播（>5000m/s）条件下仍能保持计算稳定性。

实验设计与参数设置
研究团队构建了标准化的三维数值实验体系，模型尺寸为200mm×200mm×200mm，中心位置设置φ3mm×100mm的钻孔。爆炸装药采用分层结构，外层为3mm厚聚乙烯壳体，内层为2.5mm直径的PETN装药，两者均采用GIMP方法离散化，网格尺寸0.125mm。岩体材料通过CPDI方法生成非均匀颗粒结构，控制参数包括颗粒密度（ρ=2700kg/m3）、弹性模量（E=45GPa）和断裂韧性（KIC=15MPa√m）。初始地应力场设定为最大主应力σ1=45MPa，中间主应力σ2=28MPa，最小主应力σ3=15MPa，形成典型深埋应力状态。

数值模拟结果分析
应力波传播特性呈现显著各向异性：在σ1主应力方向，应力波振幅衰减系数较垂直方向降低约40%，波速提高15%-20%。当波阵面遇到初始应力场与介质弹性模量不匹配界面时，产生约30%的应力集中效应。这种各向异性导致裂隙扩展呈现"双模态"特征：在压缩区（σ1>σ2>σ3）形成以剪切滑移为主导的径向扩展裂隙；在拉伸区（σ3<σ2<σ1）则发育呈120°夹角排列的环状张拉裂隙。

损伤演化阶段划分
研究首次系统划分了深部岩体爆破损伤的三个阶段：1）初始压缩阶段（t=0-5ms）：装药爆炸产生冲击波，在σ1方向形成压缩波前，波阵面压力达到200MPa量级；2）动态损伤累积阶段（t=5-50ms）：压缩波与初始地应力场相互作用，在σ3方向产生反射拉伸波，触发岩体表层张拉损伤；3）稳定扩展阶段（t>50ms）：形成贯穿性张拉裂隙，在σ2平面产生剪张复合型损伤。

应力场演化规律
在初始地应力场约束下，爆炸应力波呈现显著的空间衰减梯度。实验数据显示：距装药中心30m处（模拟中1.2m对应实际30m），应力波振幅衰减率较自由场提高18%-25%，波前形状由球面波演化为哑铃型波阵面。这种演变导致岩体破坏模式发生质变——在浅部表现为均匀破碎，而在深部则形成以径向压缩带和环向张拉带为特征的复合损伤区。

裂隙扩展动力学
数值模拟揭示了地应力场对裂隙扩展的三重调控机制：1）主应力方向调控：最大主应力方向（σ1）成为初始压缩裂隙的主扩展方向，延伸距离达2.8倍装药半径；2）次生应力场耦合：在σ3方向形成的反射拉伸波场诱导出环向张拉裂隙，其扩展速率较自由场提高37%；3）损伤阈值动态变化：当累积损伤率超过材料临界值（约75%）时，触发非线性加速扩展，形成应力释放主导的损伤扩展模式。

工程应用启示
研究成果为深部矿山爆破优化提供了关键理论支撑：首先，证实装药孔深与σ1方向的夹角对破碎半径影响系数达0.68，建议采用定向装药技术；其次，发现当σ1/σ3>2.5时，环向张拉裂隙间距缩小至0.4m，指导预裂孔间距设计；再者，通过控制装药与σ2平面的相对位置，可使破碎带偏转角度控制在±15°范围内，有效规避关键结构面。

技术验证与工程应用
研究团队通过标准验证模型（模型编号B-2025-03）与Banadaki（2011）的实验数据对比，在裂纹起裂点定位误差控制在8%以内，表面破坏区吻合度达92%。在四川某深层金矿的工程应用中，基于本研究的数值模型优化了爆破参数：将传统齐发装药改为延时分层装药，主装药段延时差控制在2ms以内，使破碎半径扩大22%，超挖量减少18%。实测数据显示，在σ1=48MPa的应力条件下，装药中心下方5m范围内形成直径8.3m的环形破碎带，与数值预测结果（8.1m）吻合度达96%。

研究局限性及未来方向
当前模型未充分考虑爆炸气体与裂隙的耦合流动效应，建议后续研究引入CFD-MPM混合算法。针对深部岩体多裂隙网络分布特征，可拓展研究各向异性损伤本构模型。此外，建议建立包含地应力场、装药结构、岩体各向异性参数的爆破优化决策支持系统，实现爆破参数的智能匹配。

该研究为深部地下空间开发提供了重要的理论工具，特别是在我国西部深部铜矿、煤矿等资源开发中具有重要应用价值。通过揭示地应力场与爆破动力学的耦合作用机制，为优化深井巷道爆破参数、控制岩体损伤范围、提升施工效率提供了科学依据，对保障深部矿山安全开采具有指导意义。