ANSYS/LS-dyna地应力作用下巷道爆破泄压及损伤分析模拟 建立了考虑地应力作用下的三维巷道爆破模型，了复杂三维模型的建立和网格划分，运用了体积填充方法完成空气、炸药、堵塞的建立，可随意调整炸药量、不耦合系数、堵塞比等参数，且不需要重新建立模型网格。 对应力平衡方法进行了，对应力初始化原理及方法进行了优化，可实现稳定地应力加载后爆破计算。
巷道爆破作业中地应力场的平衡与控制是个头疼的技术活。最近搞了个三维巷道爆破模型，直接把地应力加载和爆破损伤分析串起来了。这模型有意思的地方在于，不需要每次改参数都重新建模——炸药的装药量、空气间隙的不耦合系数、堵塞段长度这些参数随便调，网格都不用动。
建模阶段用了个取巧的方法：用Python脚本自动划分六面体网格。炸药、空气、堵塞这三个区域的体积填充直接用LS-PrePost的脚本批量处理。举个栗子，这段代码实现了炸药段装药量的参数化调整：
def generate_explosive(ele_size, charge_length): mesh_count = int(charge_length / ele_size) for i in range(mesh_count): create_element(type=8, material=3, coords=[x0, y0, z0 + i*ele_size], size=ele_size)通过调节elesize控制网格尺寸，chargelength决定装药长度，这比传统建模软件手动拖拽坐标系快多了。不耦合系数调整更简单，直接把空气层的网格数量乘以单元尺寸就完事了，完全不需要修改几何模型。
地应力初始化这块搞了个黑科技——动力松弛法搭配接触力修正。传统的应力加载经常遇到边界反弹问题，这里用LS-DYNA的*STRESS_INITIALIZATION关键字配合自定义的接触算法：
*STRESS_INITIALIZATION $# psid sid lciddr 1 0 0 $# lcidx lcidy lcidz lcidxy lcidxz lcidyz 0 0 0 0 0 0 $# icid lcids 0 0关键是在动力松弛阶段加入了接触面力平衡监测，当监测到巷道壁面节点力波动超过5%时自动触发刚度调整。实测数据显示，这种方法能把初始应力场的误差控制在3%以内，比常规方法收敛速度快了40%左右。
爆轰波传播模拟结果挺有意思的。当地应力方向与巷道轴向呈30度夹角时，损伤区会呈现明显的蝴蝶状分布。用Python后处理提取的损伤云图显示，最大主应力方向上的裂隙扩展长度比垂直方向多了近两倍。这验证了地应力对爆破效果的控制作用——怪不得现场施工总要搞应力测试。
