李云，刘霁

(湖南城市学院 土木工程学院，湖南 益阳，413000)

在岩土工程爆破时，由于成本和环境、设备等因素的限制，不可能对爆破施工中的各个参数全部监测，从而限制了对爆破机理的深入研究。目前，人们对岩石中爆破力学特性进行了大量研究，如浦锡锋等[1]运用 BCM模型对岩石中爆炸应力波传播、岩石破碎过程进行了模拟；龚敏等[2]运用FRLAC/ FLAC3D程序动力分析模块研究爆破作用下岩体累积损伤效应及其稳定性；王礼立[3]利用ADINA计算程序对岩体内爆炸应力波能进行计算分析；宇文惠鑫[4]使用ABAQUS/Explicit软件对岩石爆破的动力学行为进行了三维弹塑性碎裂有限元数值模拟；夏祥[5]用离散元软件UDEC计算得到爆破荷载作用下的岩体振动特征和衰减规律等。这些研究大部分集中于对爆破后的效果进行分析，简化了装药结构，数值模拟结果没有得到现场验证，影响了计算结果的可靠性。目前，国内使用较多的爆破模拟程序是LS-DYNA程序，它是1个显示非线性动力分析通用的有限元程序，可以求解各种二维、三维爆炸动力响应问题，其优势在于可以从时间和空间2方面描述爆破力学的特征[6-9]。在此，本文作者利用LS-DYNA程序，按实际工程的柱状、多孔装药模拟岩石中爆破力学的特性分析爆破的本质。

1 数值计算模型

1.1 岩体力学模型

根据试验现场情况，将岩体本构模型简化成连续介质模型。岩石的应力-应变关系呈现复杂的非线性特征，在达到屈服极限之前将介质近似地视为……