基于LS-DYNA的单孔爆破模拟

本案例利用LS-DYNA软件进行简单的单孔爆破模型，因模型较为简单，建模、划分网格以及关键字定义均在LS-PREPOST中完成。其中炮孔为圆形，直径42mm，为保证能够观察到岩石完整的破坏过程，将岩石尺寸定义为3m×3m。计算采用流固耦合算法，空气和炸药采用ALE算法，岩石采用拉格朗日算法，通过* CONSTRAINT_LAGRANGE_IN_SOLID关键字将炸药和岩石结构进行耦合。炸药不直接建模，通过关键字进行填充，所采用的关键字为*INITIAL_VOLUME_FRACTION_GEOMETRY。为炸药提供流动空间的空气域网格尺寸定义为5mm。为了模拟无限大空气和无限大岩石，在四周通过*BOUNDARY_NON_REFLECTING关键字设置透射边界，防止应力波到达边界后发生反射。其中岩石材料使用*MAT_RHT，设置neiph为4，用于查看岩石的损伤值。

岩石的压力云图变化如下，可以看到爆轰波呈圆形往四周传播，应力波在到达岩石边界时，未发生发射现象，与理论相符合。

岩石的损伤过程如下图，岩石结构在爆轰波的冲击作用下，在炮孔附近区域发生损伤，靠近炮孔周围损伤最为严重；裂纹呈不规则形状向四周扩展。

通过*DATABASE_TRACER关键字在炮孔周围距离爆炸中心1cm处设置了固定监测点，用于监测整个计算过程中监测点处压力随时间的变化情况，该位置处的压力时程曲线见下图：