在多炮爆破工程评估中单次爆破模拟已经很难满足实际需求。无论是露天矿山台阶爆破、隧道掘进连续掏槽还是城市拆除爆破中的同段多孔起爆岩体受到的从来都不是“一炮”的作用而是多次爆破载荷在时间和空间上的叠加。每一次爆破产生的振动波都会在岩体中形成新的损伤同时也会激活、扩展前一次爆破留下的微裂隙。这种振动损伤累积效应直接决定了岩体的最终破坏范围、边坡稳定性以及邻近建构筑物的安全状态。本文将基于 Hypermesh新界面完成前处理建模配合 LS-DYNA 显式动力学求解器完整梳理多次爆破岩体全时域振动损伤累积响应的分析流程。内容涵盖单位制设置、材料模型与损伤累积机制、Hypermesh 新界面网格质量检查、爆破载荷施加、K 文件导出、求解控制以及后处理提取振动速度与损伤云图最后给出常见报错排查方案和工程实践建议。如果你正在用 Hypermesh 新界面做 LS-DYNA 爆破仿真或者被“材料设置里单位怎么统一”“3D 网格质量怎么检查”“节点小球不显示”这类问题卡住这篇文章可以帮你省下不少弯路。1. 为什么要做“多次爆破”的振动损伤累积分析1.1 单次爆破模拟与多次爆破的本质区别单次爆破模拟的典型做法是建立岩体模型在炮孔附近施加爆炸载荷求解岩体在爆炸应力波作用下的破坏范围、振动速度峰值PPV和损伤分布。这种分析能够回答“这一炮会造成什么影响”但无法回答“连打三炮之后岩体是否已经进入失稳状态”。实际的爆破工程几乎都是多炮连续作业。相邻炮孔之间、前后排炮孔之间、同一区段多次起爆之间岩体的应力场和损伤场都在不断变化。第二次爆破时岩体已经不再是完整介质——它内部存在第一炮留下的损伤区、残余应变、甚至宏观裂纹。此时岩体的波阻抗、弹性模量、强度参数都和初始状态不同振动波的传播路径和衰减规律也会随之改变。因此多次爆破的全时域分析本质上是把“材料损伤的历史效应”纳入仿真过程。每一次求解都要继承上一次的应力、应变、损伤状态而不是把模型重置回初始状态重新计算。这比单次爆破模拟更接近工程真实。1.2 全时域振动响应与损伤累积的含义“全时域”在这里包含两层意思一是时间维度上的完整覆盖。从第一次爆破起爆开始到最后一次爆破的振动衰减结束整个时间历程都要纳入计算或至少纳入结果评估而不是只看每一次爆破单独的结果。二是状态变量的连续传递。损伤、塑性应变、应力状态随时间连续演化后一次爆破的初始条件正是前一次爆破的结束时刻状态。在 LS-DYNA 中实现这种累积效果常见做法有两类在同一个求解文件中按时间顺序设置多个爆破载荷曲线一次求解完成全过程使用重启动Restart技术每次求解结束后保存结果并继续下一次爆破加载。两种方式各有适用场景后文会详细展开。1.3 本文建模思路与适用范围本文的建模思路如下在 Hypermesh 新界面中完成岩体几何建模与网格划分设置统一的单位制避免参数量纲错乱选用适合岩体在爆炸高应变率下响应的材料模型并打开损伤累积相关选项在炮孔壁或自由面上施加多次爆破载荷曲线设置求解控制关键字、历史变量输出导出 K 文件在 LS-DYNA 中求解并利用 HyperView / LS-PrePost 后处理查看振动速度时程和损伤分布。适用场景包括露天台阶爆破边坡稳定性评估、隧道掘进循环爆破损伤区预测、采矿工程多孔起爆参数优化、爆破振动对邻近岩体及结构的影响分析等。2. 软件环境与版本说明2.1 软件组合与版本说明本文以 Hypermesh 新界面Ribbon 风格界面通常在 HyperMesh 2022 及之后版本中为默认界面作为前处理工具以 LS-DYNA 作为求解器。为了便于说明Hypermesh 新界面的操作路径均以 Ribbon 菜单结构描述。需要注意不同版本的 Hypermesh 菜单名称和图标位置可能存在差异。如果你的版本界面与本文描述不完全一致重点参考“功能入口的英文名称”不要卡在具体按钮位置上。2.2 单位制约定LS-DYNA 本身没有固定单位制它只是按照数值进行运算。所有输入参数的物理单位需要用户自己保证一致。单位制不一致是爆破仿真中最常见、也最难排查的错误一旦量纲出错计算结果可能相差几个数量级。爆破分析是一种典型的“大尺寸、高压力、短时程”问题。矿山或边坡岩体模型通常以米为长度单位爆炸压力以帕斯卡或兆帕为单位时间以秒或毫秒为单位。工程中常用的一组单位制如下物理量单位说明长度mm模型尺寸在毫米量级时更方便质量t吨密度按 t/mm³ 输入时间ms毫秒与显式求解时间步长匹配应力/压力N/mm²即 MPa1 N/mm² 1 MPa速度mm/ms1 mm/ms 1 m/s加速度mm/ms²或按换算关系处理在 Hypermesh 的 Materials 设置中输入密度时一定要检查当前单位制下的期望数值。例如岩体密度为 2700 kg/m³在 mm-t-ms 单位制下应输入 2.7e-9 t/mm³。这里的数量级很容易写错。更稳妥的做法是进入 Hypermesh 新界面的 Model Browser在创建材料时先明确当前模型单位制再按照对应单位制填写密度、弹性模量、强度参数。同一模型内绝对不允许混用“密度用 kg/m³、弹性模量用 Pa”这种方式。2.3 文件流转关系整个分析流程涉及三类核心文件Hypermesh 模型文件.hm保留几何、网格、载荷、材料等所有建模信息LS-DYNA 关键字文件.k由 Hypermesh 导出可直接提交求解也可用文本编辑器人工修改求解结果文件d3plot、binout、ascii 文件等供后处理软件读取。建议在 Hypermesh 建模过程中养成“每完成一个阶段就另存一个 .hm 文件”的习惯避免后期需要回退到某一版网格或载荷设置时必须重做。3. 核心原理拆解3.1 LS-DYNA 显式求解与爆破载荷施加方式LS-DYNA 采用显式中心差分法求解动力学方程。该方法的显著特点是不需要迭代求解联立方程组但要求时间步长满足稳定性条件即时间步长不能超过应力波穿过最小单元所需的时间。这在爆破模拟中意味着最小单元尺寸直接控制整体计算耗时。爆破载荷的施加方式通常有以下几种使用 *LOAD_BLAST_ENHANCED 关键字基于 CONWEP 经验公式计算空气中爆炸或表面爆炸的反射超压。这种方法适合描述自由空气中的冲击波作用在岩体内部炮孔爆破中并不完全适用使用 *LOAD_SEGMENT 或 *LOAD_NODE_POINT将炮孔壁上的压力-时间曲线映射到单元面上。这是岩体爆破前处理中最常用的方式适合模拟耦合装药或不耦合装药下的炮孔压力使用流固耦合方法将炸药和空气用 SPH 或 ALE 单元建模爆炸产物与岩体相互作用。该方法精度更高但建模复杂度和计算成本也大幅增加。对于多次爆破分析采用“炮孔壁压力曲线 时间延迟”的组合方式最为直接。每个炮孔或每组炮孔拥有独立的压力-时间曲线通过设置不同的起始时间模拟起爆延迟。3.2 岩体材料模型与损伤累积岩体在爆炸载荷下表现出高应变率敏感性、脆性破坏和拉压不对称等特征。常用的材料模型包括*MAT_PLASTIC_KINEMATIC适合弹塑性问题的快速起步支持应变率效应和失效但不包含真正的损伤累积变量*MAT_JOHNSON_HOLMQUIST_CONCRETEHJC 模型*MAT_111最初为混凝土冲击响应开发包含压力相关的强度模型、损伤累积和拉伸断裂软化被广泛用于岩体与混凝土的爆炸冲击分析*MAT_RHT在 HJC 模型基础上进一步发展包含更细致的拉伸损伤、应变率增强和剩余强度描述更适合岩体破裂过程模拟*MAT_ADD_EROSION作为辅助关键字可为任意材料模型添加失效准则。当单元达到指定等效塑性应变或最大主应力等阈值时单元被删除从而模拟宏观裂纹的生成。“损伤累积”在 LS-DYNA 中通常体现为材料模型内部的损伤变量 D 从 0 逐步增大到 1。当 D 达到失效阈值时单元失去承载能力或发生侵蚀删除。多次爆破分析中前一次求解结束时的损伤场必须传递到下一次求解的初始条件中才能形成真正的累积效应。需要注意的是不同材料模型对损伤变量的定义不同。HJC 模型中的损伤是等效塑性应变和体积应变的累积RHT 模型则分别考虑压缩损伤和拉伸损伤。选好材料之后必须理解其损伤变量对应的物理含义否则后处理看到的损伤云图可能被误读。3.3 多次爆破的实现方式这里介绍两种常用的多次爆破实现方式实际项目中可按需选择。方式一单文件多载荷时序加载在同一个 K 文件中为不同炮孔或不同起爆段位设置多条压力时间曲线通过曲线的时间起点形成起爆延迟。LS-DYNA 在一个求解进程内会连续计算所有载荷的作用损伤场自然累积。这种方式实现简单适合爆破时间间隔较短、总计算时程可控的情况。其局限性在于如果两次爆破之间的时间间隔很长例如数十毫秒甚至秒级而显式求解时间步长只有微秒级那么一个求解进程内要计算数百万步计算成本会急剧上升。方式二重启动分析LS-DYNA 支持将前一次求解结束时的应力、应变、损伤状态写入重启动文件然后通过 *RESTART 关键字在后续求解中继续加载。这样可以把“第一炮后的状态”和“第二炮的加载”分开保存和提交更加灵活也便于检查每一步的结果。实际工程中多次爆破前后通常存在较长间隔振动在间隔内已经衰减到较低水平。此时完全可以用重启动方式将前一次求解的响应截断到“振动基本衰减完成”的时刻再进入下一炮计算同时保证损伤状态连续传递。这是全时域振动损伤累积分析中兼顾精度与成本的关键做法。3.4 Hypermesh 新界面的 3D 网格质量要求爆破模拟对网格质量的要求高于一般静力分析。显式动力学求解时单元质量差会直接导致时间步长过小甚至出现负体积和计算发散。在 Hypermesh 新界面中3D 网格质量检查主要关注以下指标Jacobian雅可比反映单元畸变程度Hex 单元建议大于 0.6Tet 单元建议大于 0.5Aspect Ratio长宽比单元最长边与最短边的比值建议控制在 5:1 以内Warpage翘曲度单元面偏离平面的程度建议小于 15°Taper锥度单元形状与理想形状的偏差建议小于 0.5Skew偏斜角单元角度偏离理想角度的程度越小越好。检查步骤一般是在新界面的 3D 面板组下进入 Check Elements 或 Quality Index 面板选择 Tetra / Hex 规则设置上述阈值后执行检查。质量不合格的单元需要回到网格划分步骤重新调整尤其是炮孔附近、台阶拐角、结构面位置等应力集中区域。4. 完整实战案例Hypermesh 新界面 LS-DYNA 多次爆破岩体全时域分析下面以一个简化的岩体台阶多次爆破模型为例演示从 Hypermesh 前处理到 LS-DYNA 求解的主要步骤。模型尺寸不追求现场级精度重点展示流程。4.1 创建几何模型与网格划分在 Hypermesh 新界面中先创建或导入岩体几何模型。本文以简化台阶模型为例岩体总体尺寸为 20 m × 15 m × 10 m设置 3 个炮孔炮孔间距 2 m排距 1.5 m。在几何模型准备好之后进行网格划分炮孔附近区域单元尺寸控制在 100 mm 左右保证爆炸压力梯度能够被较好描述远离炮孔区域单元尺寸逐渐过渡到 500 mm 左右降低总单元数量单元类型优先使用六面体单元Hex复杂区域可局部使用四面体单元Tet。在 Hypermesh 中可以使用 Solid 面板下的自动网格划分功能。先划分二维面网格再通过 Solid Map 或 Drag 等方式生成三维网格。4.2 检查 3D 网格质量在 Hypermesh 新界面中执行以下操作进入 Mesh 组下的 Check Elements 面板在 3D / Tet 或 3D / Hex 规则中选择需要检查的单元类型设置各项阈值点击 Check 按钮查看不合格单元数量对不合格单元使用 Quality Index 面板进行自动优化或手动调整节点位置。批量修复时需要注意的是自动优化可能导致网格拓扑变化修复完成后应重新检查单元质量并确认模型边界没有出现穿透或缝隙。4.3 设置材料与单元属性在 Hypermesh 新界面的 Model Browser 中创建材料。这里以 HJC 模型为主材料模型并配合 *MAT_ADD_EROSION 定义失效条件。创建一个材料并分配属性时需要填写以下基础参数RO密度弹性模量或剪切模量泊松比抗压强度、抗拉强度损伤参数 D1、D2 等下面是一个示例材料参数表仅演示数值量纲实际参数必须根据岩体力学实验或文献标定参数示例值说明RO2.4e-9密度单位 t/mm³G1.1e4剪切模量MPaA0.79归一化内聚强度B1.6归一化压力硬化系数C0.007应变率系数FS0.04失效类型参数在 Hypermesh 新界面的 Materials 设置中并不存在一个专门的“单位制输入框”。单位正确与否完全取决于你对密度、模量、强度输入的数值是否按照统一单位制换算。这正是在 2.2 节强调单位制约定的原因。创建好材料后在 Model Browser 中创建 Components为每个 Component 分配单元类型*SECTION_SOLID和材料。炮孔附近、远区、边界区可以划分到不同 Component 中便于后续分别赋予材料参数。4.4 施加爆破载荷与边界条件边界条件方面岩体底面和侧面设置无反射边界*BOUNDARY_NON_REFLECTING用于模拟无限岩体边界避免应力波在人工边界处反射干扰结果。模型顶部为自由面。爆破载荷方面在每个炮孔壁面单元上施加压力-时间曲线。为模拟多次爆破为不同炮孔设置不同起始时间。例如炮孔 1起始时间 0 ms峰值压力 100 MPa炮孔 2起始时间 2 ms峰值压力 100 MPa炮孔 3起始时间 4 ms峰值压力 100 MPa。在 Hypermesh 中可以创建 *DEFINE_CURVE 曲线并将曲线施加到对应 Component 的单元面上。压力曲线的上升段通常非常短微秒量级下降段则根据装药量和炮孔参数设定。曲线数据点的数量会影响求解速度建议在保证曲线形状的前提下尽量减少数据点。施加完成后在 Model Browser 中检查每个 Load Collector 关联的 Component 是否正确避免出现压力施加到错误单元面的情况。4.5 求解控制关键字设置在 LS-DYNA 关键字控制面板中需要设置以下内容*CONTROL_TERMINATION设置总求解时间。多次爆破分析需要覆盖所有爆破载荷作用完成的时刻例如设置为 20 ms*CONTROL_TIMESTEP设置时间步缩放因子一般取 0.60.9初始时间步根据最小单元尺寸自动估算*DATABASE_BINARY_D3PLOT设置结果输出频率。爆破过程应力波传播速度快建议每 0.02 ms 输出一帧以保证后处理时能捕捉到振动波传播过程*DATABASE_HISTORY_NODE选择需要输出振动速度时程的节点。在模型表面和关注位置设置历史节点用于提取全时域振动速度*DATABASE_ASCIIOPTION根据需求输出 NODOUT、GLSTAT 等文件。在 Hypermesh 新界面的 LS-DYNA 面板组下上述关键字可以通过界面面板配置也可以在导出 K 文件后手动编辑。对于多次爆破分析建议导出 K 文件后检查一遍关键字顺序和参数因为 Hypermesh 界面配置有时无法覆盖所有高级选项。4.6 导出 K 文件并提交求解在 Hypermesh 新界面中执行 File → Export → Simulation Deck选择 LS-DYNA 模板导出 K 文件。导出前重点检查材料卡片是否完整曲线定义是否引用正确所有 Component 是否都分配了材料和属性是否有多余的未使用 Component 混入导出文件。将 K 文件提交给 LS-DYNA 求解。命令行提交方式如下ls-dyna iblasting_multi.k ncpu8 memory1000m不同版本的 LS-DYNA 求解器命令略有差异具体以你安装的 MPP 求解器版本为准。求解过程中注意查看 mesage 文件和 d3hsp 文件如果出现负体积、时间步骤降等提示需要停止求解并回到前处理中修改网格或参数。4.7 结果后处理振动速度时程与损伤云图求解完成后使用 HyperView 或 LS-PrePost 打开 d3plot 文件。查看振动速度时程在 History 面板中选择预先设置的历史节点绘制节点速度-时间曲线对比不同炮孔起爆时刻的振动速度峰值观察 PPV 随爆破次数的变化规律。查看损伤累积在 Contour 面板中选择损伤变量不同材料模型变量名不同HJC 模型中通常为损伤 D 或 fail 相关变量分别查看第 1 次爆破结束、第 2 次爆破结束和最终时刻的损伤云图对比不同时刻的损伤范围扩展情况识别损伤累积的主要方向。如果采用重启动方式则需要在每一次重启动求解完成后单独保存结果然后在后处理中合并或叠加查看。5. 常见问题与排查方法多次爆破岩体仿真过程中容易碰到的问题主要集中在网格质量、材料参数、单位制、显示设置以及求解稳定性几个方面。下面将实际项目中常见的现象、原因与解决思路整理成表。问题现象常见原因解决思路Hypermesh 新界面中所有节点都不显示显示设置被切换为仅显示单元或图形区 Graphics 选项被重置在 Display 设置中打开节点显示或按快捷键调整显示模式节点小球不显示新界面中节点显示样式为点状而非球状或被隐藏检查 Visualize 视图菜单下的节点显示样式设置材料参数设置了但 K 文件中缺少材料关键字材料没有分配给 Component或导出的 Component 没有关联属性检查 Model Browser 中每个 Component 的 Property 是否完整3D 网格质量检查失败六面体单元在复杂几何区域畸变严重拆分几何特征局部重新划分使用 Quality Index 自动优化求解开始后很快出现负体积单元过小、材料参数单位制错误或时间步长不合理先检查单位制统一性再检查最小单元尺寸和 *CONTROL_TIMESTEP 设置多炮载荷作用时间不准确曲线数据点设置错误或曲线起始时间未正确指定逐个炮孔检查 *DEFINE_CURVE 的起始时间和峰值位置振动速度峰值为 0 或异常大历史节点未正确关联或单位制换算导致速度量纲错误核对历史节点 ID、单位制换算表和后处理坐标轴标签损伤云图始终为 0材料模型未开启损伤累积或后处理中变量选择错误确认材料模型支持损伤变量并在后处理变量列表中查找 D 或 damage 类字段LS-DYNA 求解中断提示“Out of memory”模型单元过多、内存设置不足改用 MPP 并行版或在提交命令中增加 memory 参数模型中出现应力波在边界反射未设置无反射边界或无反射边界类型不匹配为模型底面和侧面设置 *BOUNDARY_NON_REFLECTING多次爆破求解时间过长两次爆破间隔过大计算总时程过长改用重启动方式分段计算合理截断振动衰减阶段在这些问题中有两个很容易被忽略第一个是单位制。很多刚接触 LS-DYNA 的工程师直接在 Materials 里填“2700”觉得这就是岩体密度结果在 mm-ms 单位制下实际被解释为 2700 t/mm³模型质量瞬间变成天文数字求解必然发散。第二个是重启动状态传递。有些人把多次爆破拆成多个独立求解文件每个文件都从初始状态开始重新计算。这样得到的不是累积损伤而是多炮各自独立的结果完全失去“累积响应”的物理意义。必须使用重启动功能或在同一求解进程中按时间序列加载多炮载荷曲线。6. 工程实践建议6.1 建模前先明确关注量爆破振动损伤分析中不同关注量对模型的要求差别很大。如果你关心的是岩体最终损伤范围和破裂形态那么网格密度、材料失效参数和损伤模型的选择优先度最高炮孔附近需要精细网格。如果你关心的是远区振动速度PPV那么模型尺寸、无反射边界、材料波阻抗匹配和阻尼设置更加重要炮孔附近的网格精度可以适当放宽。在实际项目中往往需要兼顾两种分析目标。建议分两个模型分别处理近区精细模型用于损伤演化分析远区大模型用于振动传播评估二者在边界处通过振动速度或应力波时程进行衔接。6.2 材料参数标定是精度瓶颈岩体不是均质材料其力学参数受节理裂隙、含水率、风化程度等因素影响。LS-DYNA 中任何一种高级材料模型都需要依赖准确的实验参数。直接照搬文献参数虽然方便但结果可能偏离现场实际。比较务实的做法是通过室内岩石力学实验获得密度、弹性模量、抗压强度、抗拉强度等基础参数通过 SHPB 实验或文献数据获得高应变率下的强度增强曲线通过现场爆破试验测得的振动速度衰减数据反向标定材料参数和阻尼使用 *MAT_ADD_EROSION 时失效参数应通过小规模试算和现场破坏范围对比校准。6.3 全时域分析的工程化处理严格意义上的全时域要求在几十毫秒甚至更长时间内持续计算。受限于显式求解时间步长其计算成本无法回避。这里给出几条工程化处理建议优先采用重启动分段计算每一段只计算到振动衰减到较低水平为止在两次爆破之间保留足够长的时间间隔但要合理设置阻尼避免振动长期不衰减导致计算时程过长如果只需要评估最终损伤云图和振动峰值可以适当压缩爆破间隔但要注意压力波叠加效应可能被人为改变如果爆破间期非常长且振动衰减后岩体应力场对下一炮影响不大可以只保留损伤状态传递不必全程计算振动衰减过程。这些处理方式在学术研究中可能需要严格验证但在工程评估中属于兼顾精度和效率的常用手段。6.4 安全与合规提示爆破工程仿真涉及实际工程安全评估。仿真结果并不能直接等同于工程结论必须结合现场爆破试验数据、地质勘察资料和规范要求进行综合判断。在进行边坡稳定性评价、临近建筑物安全评估等敏感用途时应由具备相应资质的专业机构完成计算和报告。使用 LS-DYNA 进行爆破模拟请确保你拥有合法的软件授权并遵循所在单位的仿真管理规范。7. 总结与下一步学习方向本文围绕 Hypermesh新界面配合 LS-DYNA 的多次爆破岩体全时域振动损伤累积响应梳理了完整的分析流程。核心要点可以归纳为多次爆破分析的关键不是把单炮结果相加而是保持损伤状态在时间上的连续传递单位制统一是 LS-DYNA 仿真中最基础也最重要的前提材料设置前先确定单位制3D 网格质量直接决定求解稳定性和结果可信度炮孔附近必须精细划分并严格检查质量多次爆破可以通过单文件多载荷时序加载实现也可以通过重启动分段实现后者更适合长时间间隔的工程场景后处理重点关注历史节点的振动速度时程和材料损伤变量分布并结合多次结果对比损伤扩展规律。如果你刚接触这个方向建议不要直接冲击大规模现场模型。先从单孔爆破模型开始跑通 Hypermesh 前处理、K 文件导出、LS-DYNA 求解、HyperView 后处理的完整链路再逐步增加炮孔数量、细化材料模型、引入重启动和多次爆破逻辑。每一步都用标准算例验证结果符合物理常识再向复杂模型推进。下一步可以继续研究的方向包括SPH 与 ALE 方法的爆炸模拟精度对比、岩体节理裂隙的离散模型建模、爆破振动信号的经验模态分解与频域特征提取、以及基于实测数据的材料参数反演与模型校准。希望这篇文章能帮你在 Hypermesh 新界面和 LS-DYNA 的爆破仿真路上少踩几个坑。如果你在配置单位制、检查网格质量或设置多次爆破载荷时遇到其他问题欢迎在评论区留言交流。