1. 项目概述:台阶爆破模拟的技术价值
在矿山开采、隧道掘进等工程爆破领域,台阶爆破(也称梯段爆破)是最常用的爆破技术之一。这种爆破方式通过将岩体分成若干台阶进行逐层爆破,能有效控制爆破范围、提高破碎效果。但传统爆破方案设计主要依赖经验公式,难以精确预测爆破效果和周边岩体损伤情况。
ANSYS/LS-DYNA作为显式动力学分析领域的标杆软件,其强大的非线性计算能力和丰富的材料模型,特别适合模拟爆破这种瞬态、高应变的复杂物理过程。通过建立台阶爆破数值模型,工程师可以在施工前:
- 可视化爆破应力波传播路径
- 量化评估爆破振动对周边岩体的影响
- 优化孔网参数(孔距、排距、超深等)
- 预测爆破块度分布
提示:LS-DYNA的MAT_HIGH_EXPLOSIVE_BURN材料模型配合EOS_JWL状态方程,是目前模拟炸药爆轰过程最成熟的解决方案之一。
2. 模型搭建关键步骤
2.1 几何建模技巧
在ANSYS Workbench中创建台阶模型时,推荐使用DesignModeler的"Extrude"功能构建阶梯状几何体。实际操作中需注意:
- 基准平面创建:通过"Plane"工具建立与主平面成倾角的参考面(对应台阶坡面角)
- 布尔运算:用"Slice"工具分割岩体形成台阶,建议保留0.3-0.5m的台阶根部
- 炮孔建模:使用"Body of Revolution"旋转生成圆柱体,直径按实际装药孔径设置(常见76-102mm)
apdl复制! 典型APDL命令流片段
CYL4,0,0,0.05, , , ,0.5 ! 创建直径100mm、深0.5m的炮孔
VGEN,2,ALL, , ,2, , , ,1 ! 阵列生成多排炮孔
2.2 材料参数设置
岩体材料建议采用*MAT_JOHNSON_HOLMQUIST_CONCRETE(JH-2)模型,该模型能准确反映高应变率下的岩石力学行为。关键参数包括:
- 密度:花岗岩通常取2600-2800kg/m³
- 弹性模量:30-60GPa
- 抗压强度:100-200MPa
- 应变率系数C取0.007,P取3.0
炸药材料使用MAT_HIGH_EXPLOSIVE_BURN配合EOS_JWL:
k复制$ 乳化炸药典型参数
*MAT_HIGH_EXPLOSIVE_BURN
1 1.35 0.693 6.5e9
*EOS_JWL
1 3.712e11 3.231e9 4.15 0.95 1.0 7.0e9
2.3 网格划分策略
采用显式动力学分析特有的网格处理技术:
- 炮孔周边加密网格(5-10mm单元尺寸)
- 采用六面体主导网格(Hex Dominant Method)
- 设置沙漏控制(HOURGLASS=5)
- 添加无反射边界(*BOUNDARY_NON_REFLECTING)
经验:在台阶坡面处网格应沿45°方向偏置,避免出现极端锐角单元。
3. 求解设置与损伤评估
3.1 接触算法选择
- 炸药与岩石:*CONTACT_ERODING_SURFACE_TO_SURFACE
- 岩石破碎体:*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE
- 设置SOFT=1(软接触)防止初始穿透
3.2 求解控制参数
k复制*CONTROL_TERMINATION
0.05 ! 终止时间(秒)
*CONTROL_TIMESTEP
0.67 ! 时间步缩放因子
*DATABASE_BINARY_D3PLOT
5.0e-5 ! 结果输出间隔
3.3 损伤评估方法
- 塑性应变判据:当等效塑性应变>0.01时认为岩体破坏
- 损伤云图:通过*DATABASE_EXTENT_BINARY输出STATUS变量
- 振动监测:在关键位置设置*SET_NODE_LIST采集振动速度
k复制*DATABASE_HISTORY_NODE
1001 ! 监测节点编号
4. 常见问题解决方案
4.1 计算发散处理
- 现象:计算在初始阶段崩溃
- 解决方案:
- 检查材料单位制一致性(常用cm-g-μs制)
- 降低初始时间步缩放因子(0.3-0.5)
- 添加*CONTROL_SHELL控制翘曲
4.2 能量异常排查
- 健康能量比例:内能/总能量应在70%-90%
- 异常处理:
- 检查沙漏能占比(应<5%)
- 增加体积粘性系数(QV1=0.06)
- 调整接触刚度(SFS=0.1)
4.3 后处理技巧
- 在LS-PrePost中使用Fringe Range功能聚焦损伤区域
- 通过"History>Math"功能绘制振动速度衰减曲线
- 使用"Split Window"对比不同时刻应力云图
5. 工程应用实例
某铁矿台阶爆破参数优化案例:
- 原方案:孔距5m×排距4m,单耗0.45kg/m³
- 模拟发现:底部抵抗线不足导致根底
- 优化方案:调整为4.5m×3.8m,超深增加0.5m
- 效果:大块率降低37%,振动速度下降22%
实际建模时发现,当采用分段装药结构时,需要特别注意:
- 空气间隔段需单独建模
- 使用*INITIAL_VOLUME_FRACTION_GEOMETRY定义装药位置
- 起爆点设置采用*INITIAL_DETONATION
在最近参与的隧道掘进项目中,通过对比三种不同起爆网络(逐孔、排间、V型)的模拟结果,最终采用的V型起爆使振动峰值降低了约35%。这个过程中最关键的发现是:相邻炮孔的起爆时差应控制在3-8ms范围内才能实现最优叠加效果。
