做岩土爆破动力分析很多人的困惑是从“单次爆破”转向“多次爆破”开始的。单次爆破的应力波传播、单元损伤云图、测点速度曲线按教程跑一遍基本能出来可一旦遇到隧道循环进尺、矿山多排孔爆破、硐室群反复开挖问题立刻变复杂第一爆后岩体已经损伤第二爆的载荷是在一个“伤痕累累”的介质里传播材料刚度和强度都变了如果每次计算都从零初始状态开始结果等于把两场互不相关的单次爆破拼在一起谈不上损伤累积。从更底层的角度看多次爆破模拟能不能成立取决于三件事材料模型是否具备损伤累积能力多炮孔载荷能否按时间次序生效两次爆破之间的应力、应变、损伤状态是否被完整传递。Hypermesh新界面主要承担前处理LS-DYNA负责求解两者缺一不可。这篇文章会从实际建模痛点出发讲清楚多次爆破岩体全时域振动损伤累积响应的实现路线、Hypermesh新界面下的网格与单位处理、LS-DYNA材料模型选择以及一套可以直接迁移的k文件模板。如果你正在做隧道爆破振动预测、矿山台阶爆破围岩损伤分析或者准备把“单次爆破”升级为“多次爆破”这篇文章可以帮你少走很多弯路。读完你会明白全时域不是把每一微秒都算满而是用合理的状态传递方式把多次爆破的振动波和损伤场完整连接起来。1. 多次爆破与单次爆破的仿真差异1.1 为什么不能简单延长计算时间很多新手第一次做多次爆破时会直接把两次炸药放进同一个模型然后拉长求解时间认为“算完第一次接着算第二次”就是全时域。这个想法方向没错但工程上很难落地。原因在于LS-DYNA显式积分的稳定性条件。时间步长由最小单元尺寸和材料声速决定爆炸近区的网格通常很密时间步可能小到微秒级甚至亚微秒级。如果两次爆破的间隔是分钟级或小时级按真实时间积分意味着海量的时间步单机几乎不可能完成。即便是毫秒级微差爆破全部空间单元都要在两次起爆之间继续积分计算量也会成倍放大。更关键的是爆炸载荷作用时间很短但岩体的应力重分布、损伤区扩展、振动波衰减都需要一段时间。实际工程中两次爆破之间的“静默期”对整个损伤累积结果没有直接贡献却让计算成本变得不可接受。所以多次爆破仿真的核心思路不是“把所有时间都算满”而是把每次爆破后的关键状态保存下来在下次爆破开始时恢复并继续计算。这才是“全时域”的正确含义。1.2 单次爆破与多次爆破的本质区别对比维度单次爆破多次爆破初始状态岩体应力通常按原岩应力或零应力处理第一次爆破后的应力、应变、损伤场保留材料响应材料参数基本恒定损伤区材料刚度、强度劣化波阻抗改变振动特征一个脉冲波列多个波列叠加后次爆破振动可能因损伤而放大结果输出峰值速度、单次损伤云图各测点PPV时程、累积损伤区扩展计算策略单次求解即可需要重启动、状态映射或统一时序控制单次爆破分析只看“这次爆破产生了多大振动、多少损伤”多次爆破分析还要回答另一个问题前一次爆破留下的损伤会不会让后一次爆破的振动波传得更远、损伤扩展得更快。这在工程上是极有价值的判断因为围岩稳定往往不是被某一次爆破击穿而是被多次爆破累积损伤“磨”坏的。1.3 数值实现的两个突破口从实现层面看多次爆破累积损伤需要解决两个问题。第一是材料层面的累积能力。如果材料模型本身没有损伤变量或者损伤变量在单元删除后被重置那后一次爆破就不知道前一次爆破留下了什么。常用的带损伤材料模型如LS-DYNA中的MAT_272RHT和MAT_111HJC内部都定义了损伤变量的演化与累积这是做多次爆破的重要基础。第二是求解层面的状态传递。状态传递常见做法有三种全程多炮孔统一模型、局部重启动、完全重启动。选哪种取决于两次爆破之间网格是否变化、炮孔位置是否变化、间隔时间是长还是短。这部分决定了整个项目的技术架构也是很多教程讲得最少的地方。2. 两种多次爆破技术路线2.1 路线A单模型多炮孔延时起爆如果两次爆破是同一轮开挖中的“微差爆破”炮孔位置和网格在求解过程中都不改变可以采用单模型多炮孔方案。第一次炸药先起爆计算到一定时刻第二次炸药再通过*INITIAL_DETONATION指定的起爆时间触发。这样一个求解流程就能输出完整的多波列振动过程。这种方案的优点是模型统一应力波传播、损伤演化和振动响应都是真实连续的不需要额外处理状态传递。炮孔之间的延时通常只有几十毫秒时间尺度上可以接受因此很适合研究毫秒延时对振动叠加的影响。缺点是总求解时间仍然受限。如果两次爆破间隔达到秒级单元数量又大计算成本仍然会肉眼可见地上升。而且所有炮孔位置必须在建模阶段就完全确定如果第二次爆破炮孔布置变了这个方案就用不了。2.2 路线B分步重启动工程中更常见的“多次爆破”是循环进尺爆破第一次爆破开挖一段掌子面出渣、钻孔、装药再进行第二次爆破两次间隔可能是几小时甚至几天。这种情况下炮孔位置往往每次都不一样网格甚至可能因开挖边界改变而重新构建。此时需要分步求解。第一次爆破计算完成后把应力、应变、损伤等状态写入d3dump文件第二次计算开始时使用重启动文件把状态读回来再施加新的爆破载荷。这种方案的优点是时间间隔不再受计算量限制可以精确模拟“每次爆破都是独立装药、但岩体状态是继承的”真实过程。缺点是实现难度高特别是第二次爆破需要新增炮孔或修改网格时涉及局部重启动和完全重启动的版本差异问题。2.3 路线对比与选型建议对比项路线A单模型多炮孔路线B分步重启动网格是否变化不变可不变可变化炮孔位置建模时固定每次爆破可单独设置适用时间间隔毫秒级到秒级任意间隔状态传递自动连续通过d3dump/应力映射实现难度低中高结果连续性完全连续取决于重启动方法我的建议是优先尝试路线A因为流程简单、结果可视化直观如果项目遇到“每次爆破网格重新划分”的情况再切换到路线B。先用小模型把两种路线都跑通再扩大网格规模是比较稳的工程路径。3. Hypermesh新界面下的前处理流程3.1 新界面与经典界面的差异Hypermesh新界面采用Ribbon风格与经典版的菜单栏和图标操作相比变化很大。新界面把操作按流程组织Geometry、Mesh、Setup、Analysis、Post等模块比较清晰同时依赖模型浏览器和流程模板适合新建模型时快速推进。对新用户来说最大的不习惯是“找不到旧按钮”。比如经典面板里的2D网格划分、3D网格划分在新界面中被整合到Mesh页签对应的子面板里。材料卡片编辑则在Setup或Model浏览器中完成。只要你接受了“按流程找工具”的思路上手速度其实比经典界面快。对于LS-DYNA前处理思路是一样的先在Geometry中完成几何准备再到Mesh中做2D/3D网格最后在Setup中指定材料、属性和求解控制卡片导出k文件。3.2 Hypermesh材料设置里怎么设置单位这是个高频问题也是新手最容易出错的地方。Hypermesh本身没有内置单位制你在Materials面板里填的每一个数值最终都会原样写入LS-DYNA的k文件。LS-DYNA也没有“单位自动换算”功能它只要求你保持单位一致。所以在Hypermesh的Materials设置里“设置单位”实际上是指你自己选定一套单位制然后把所有物理量换算到这套单位制下再去填表。如果模型用m-kg-s密度填2600弹性模量填3e10如果改用mm-kg-ms密度要填2.6e-9弹性模量要填30因为在这个单位制下应力单位是GPa。单位制长度质量时间密度示例应力/压力单位速度单位m-kg-smkgs2600 kg/m³Pam/smm-kg-msmmkgms2.6e-9 kg/mm³GPamm/ms即m/scm-g-μscmgμs2.6 g/cm³Mbarcm/μs建议在项目开始前把单位制写在建模说明文档里Hypermesh中每个材料卡片旁边也可以加注释避免三个月后回头看模型时忘记当初用的是哪套单位。3.3 几何导入与简化爆破模型往往由岩体、炸药、可能还有空气区组成。从CAD导入后第一步是几何清理删除重复曲面、补齐缺失面、去除对动力响应影响极小的圆角和小凸台。新界面的Geometry页签里有自动清理和手动清理工具建议优先用自动清理再检查自由边。几何清理不干净的后果是网格划分时出现畸形单元最终在LS-DYNA计算中产生负体积或沙漏。3.4 材料、属性与组件分配使用Hypermesh新界面时推荐先创建Component再创建Material和Property最后把Component关联到Property。不要把材料直接赋在单元上否则后处理和管理都会乱。在Model浏览器中创建材料后在Card Editor里按当前单位制填写参数。LS-DYNA的双线性弹塑性材料卡片、JWL状态方程卡片都在这里维护。做完这一步建议导出一次k文件用文本编辑器检查数值是否与预期一致这一步能避免很多低级错误。4. 三维网格划分与质量检查4.1 网格策略近区密、远区疏爆破模型的网格尺寸不能全场统一否则要么近区精度不够要么远区单元太多算不动。工程上通常把爆源附近5到10倍炮孔直径的区域划为加密区单元尺寸控制在应力波高频分量的波长范围以内向外逐渐过渡到稀疏网格。爆炸载荷在岩体中产生的高频成分衰减很快近区网格太粗时损伤区形状会失真振动速度峰值也会被严重低估。从网格角度看这往往比材料模型参数影响更大。4.2 Hypermesh新界面下的3D网格划分几何准备完成后进入Mesh页签的3D子面板。常用操作是先在岩体表面划分2D三角形或四边形网格再由2D网格生成3D四面体或六面体网格。对于爆破模型六面体网格在计算稳定性和单元数量上都有优势但复杂岩体几何往往难以全部六面体化。四面体网格比较灵活但注意低阶四面体单元在爆炸冲击下容易产生负体积需要严格控制网格质量。划分完成后建议先不要急着做求解先检查一遍单元质量再进入材料设置。4.3 Hypermesh如何检查3D网格质量新界面中网格质量检查主要集中在Mesh页签的Check Elements或Quality Index工具里。2D网格检查的是翘曲、长宽比、四边形角度、三角形