岩体爆破工程中单次爆破的振动与损伤已经很难在数值模型里准确复现当同一区域在短周期内经历多次爆破时每次爆炸产生的应力波都会在上一次损伤基础上继续传播振动时程和损伤场会变得高度非线性。要实现“多次爆破岩体全时域振动损伤累积响应”分析Hypermesh 新界面加 LS-DYNA 显式求解是一套很常用的工作流。Hypermesh 负责几何清理、网格划分、Part 和材料卡片组织LS-DYNA 负责高应变率下应力波传播、材料损伤演化、爆炸载荷与整个时域的时间积分。这篇文章围绕一个可落地的多次爆破算例从材料准备、网格质量、载荷时序、求解控制和后处理验证逐步展开并给出高频问题与排查路径。1. 先理解“多次爆破全时域振动损伤累积”模拟要完成什么1.1 这类模拟在工程中的实际场景隧道掘进、露天台阶爆破、基坑开挖、矿山多次分段爆破都会遇到同一个问题第一次爆破后岩体已经出现微裂纹和损伤第二次爆破的应力波传播路径、衰减规律、振动强度都和新鲜岩体不同。如果只做单次爆破分析即使单次结果和实测 PPV 吻合也无法回答“连续爆破后岩体是否损伤加重、振动是否放大”这类工程问题。多次爆破全时域振动损伤累积响应就是把多个爆源、多次起爆时序放在同一个求解过程中让岩体材料从初始完整状态开始经历第一次爆破、损伤演化、振动衰减再经历第二次爆破、损伤继续累积直到全时域结束。这样得到的结果既包含每一段爆破的振动速度时程也包含最终岩体的累积损伤分布还能看到中间卸载和重新加载过程中材料的力学状态变化。1.2 全时域累积和单次爆破的本质差别单次爆破模型通常只在 t0 附近加一个爆炸载荷计算几毫秒到几百毫秒就结束。多次爆破全时域模型至少包含两个差别第一载荷是时间序列。第二次爆破可能发生在 50 ms、200 ms 甚至更晚前一次爆破产生的余振尚未完全消失第二次爆炸应力波就叠加进来。这个叠加过程只能用全时域显式积分完成不能把两次爆破拆成两个独立文件简单相加。第二岩体力学状态不能重置。LS-DYNA 中如果材料模型是带损伤变量的模型比如 RHT、JHC 等那么同一套单元在第一段计算结束后已经保存了损伤、塑性应变和残余应力。全时域求解时第二段爆破会在这个状态下继续计算这才叫“累积”。用表格可以更直观地看两者的差异对比项单次爆破模型多次爆破全时域模型爆源数量通常一个多个可不同位置、不同起爆时间载荷时序只有一次爆炸压力多个压力脉冲按时间排队材料状态从初始状态开始保留上一段爆破的损伤和塑性振动评价关心单次 PPV 和主频关心各时段 PPV、持续时间、叠加效应最终输出单次振动和损伤全时域振动曲线和累积损伤场计算成本相对低时间步小、时长长成本成倍上升1.3 Hypermesh 新界面与 LS-DYNA 在本流程中的分工Hypermesh 新界面基于任务型的 Ribbon 布局相比传统菜单几何和网格入口更集中。它在本流程中负责切换 LS-DYNA User Profile建立符合 LS-DYNA 关键字的模型树创建材料卡片、Section、Part清理 CAD 几何划分六面体或高质量四面体网格检查网格质量修复负雅可比单元施加边界条件、约束、阻尼、监测点导出 K 文件并在导出前检查关键字卡片。LS-DYNA 则负责真正的求解显式中心差分时间积分爆炸压力在单元面上转化为节点力应力波在岩体网格中传播材料损伤累积和单元失效输出全场 d3plot 和节点历史时程文件。整条流程的技术主线可以概括为在 Hypermesh 中把模型准备好在 LS-DYNA 中跑完整时间域在后处理中同时读振动和损伤。任何一步如果只做“看起来差不多”后续排错都会非常痛苦。2. 模型准备阶段单位制、材料模型和网格质量必须一次性过关2.1 单位制materials 里没有“单位”按钮全靠数值换算在 Hypermesh 的 Materials 设置里很多人第一反应是找一个单位下拉框。实际没有。LS-DYNA 本身不关心单位它只按你填入的数值计算。所谓单位制是你自己约定长度、质量、时间三个基本量纲后把密度、弹性模量、波速、压力全部换算为同一组单位。常用单位制包括基本量纲常见组合 1常见组合 2常见组合 3长度mmmcm质量kgkgg时间smsus力NkNdyne 或 N?应力/压力PaMPa10^5 Pa 量级需统一密度kg/m3kg/mm3g/cm3速度m/smm/mscm/us这里最容易错的是密度和弹性模量之间的量级关系。以组合 1m-kg-s-Pa为例岩体密度可以填 2600单位是 kg/m3弹性模量可以填 50e9单位是 Pa如果误把密度按 g/cm3 填成 2.6整个模型的应力波速度会差出几十倍结果几乎必错。推荐做法在一个电子表格里先写好换算关系所有材料参数都过一遍再填进 Hypermesh。尤其是使用 *MAT_RHT、*MAT_JHC 这类高参数模型时单位换算错一处计算发散或振动频谱完全失真都正常。2.2 在 Hypermesh 新界面中创建岩体材料卡片Hypermesh 新界面下先确认 User Profile 是 LS-DYNA。如果还在 General 或其他求解器模板下很多 LS-DYNA 关键字卡片不会出现。创建材料的标准路径通常是启动 Hypermesh选择 LS-DYNA 用户配置在模型树中打开 Materials右键选择 Create输入材料名选择卡片模板。比如 *MAT_RHT 或 *MAT_PLASTIC_KINEMATIC进入 Card Edit按单位制换算后的数值填入密度、弹性模量、剪切模量、强度等字段保存材料。给材料起名时不要用中文避免 K 文件导出出现编码问题。新界面里模型树、材料卡片、Component 和 Load Collector 的层级关系比传统界面更直观。材料卡片不是存在网格上而是存在*MAT卡片中Component 通过 Section 引用材料Part 再引用 Component 和 Section。2.3 以 MAT_RHT 为案例拆解岩体材料参数RHT 模型适合岩石和混凝土类材料可以描述拉伸损伤、压缩损伤、应变率强化和压力相关强度适合爆破这类高应变率问题。它的完整参数很多这里给出一组用于案例分析的基础参数实际项目必须通过岩石室内试验或参考已发表文献校核。参数含义示例值m-kg-s-Pa 单位制RO密度2600.0G剪切模量20.0e9FC单轴抗压强度100.0e6FT单轴抗拉强度8.0e6D1压缩损伤参数0.04D2压缩损伤参数1.0EPF失效塑性应变0.01 左右需做敏感性测试K 文件中的材料卡片示意如下*MAT_RHT $ MID RO G ... 1 2600.0 20.0e9 $ ... 100.0e6 8.0e6这里必须强调RHT 模型不是所有字段都必须手填。Hypermesh 材料卡片会列出很多参数有些有默认值有些必须根据试验设定。最容易犯的错是只填密度和弹性模量其他都用默认值然后计算不收敛。建议在正式计算前先用单单元或单炮孔小模型验证材料在典型爆炸载荷下的破坏和损伤规律。对于强调脆性拉裂破坏的岩体还可以配合 *MAT_ADD_EROSION 设置拉伸失效、主应变失效等准则控制单元删除。但不要一开始就加很激进的侵蚀否则可能出现“岩体像炸药一样飞散”的伪结果。2.4 3D 网格质量检查方法与标准Hypermesh 中搜索“如何检查 3D 网格质量”时通常指的是 Quality 面板或 Check Elements 面板。新界面下可以从 Mesh 任务里找到 Check Elements或通过质量指标面板检查Jacobian雅可比越低表示单元畸变越严重。爆破显式分析建议六面体单元 Jacobian 不小于 0.7。Warpage翘曲角建议小于 10 度热点区更严。Aspect Ratio长宽比尽量小于 5。Min Length最小单元边长直接决定 LS-DYNA 最小时间步。单元最大最小角六面体尽量接近 90 度避免尖角单元。操作时可以先使用 Quality 面板设定阈值再逐个查看失效单元。质量差的单元可以用网格重划分或手动移动节点修复。需要注意在 Hypermesh 中质量图表不断变化是正常现象关键是失效单元不要集中在炮孔附近和监测路径上。从爆破应力波传播角度看炮孔附近网格尺寸直接决定压力曲线能否被分辨。一般要求最小单元尺寸与爆炸压力上升时间对应的波长匹配。如果压力上升时间是 0.5 ms岩体波速约 4000 m/s特征波长为 2 m那么热点区网格至少应小于波长四分之一即 0.5 m。更精细时可以采用 0.1 m 到 0.2 m。太粗会滤波太细则时间步小计算量爆炸性增长。2.5 节点不显示小球的显示问题排查很多人会碰到一个很让人困惑的现象模型看起来只有单元没有节点小球。这不是网格丢失绝大多数是显示设置问题。在 Hypermesh 新界面中检查顺序如下确认当前 Component 没有隐藏。Component 隐藏后节点和单元都会一起消失。查看显示控制面板中的 Node 开关。新界面的 View 或 Display 设置里节点实体可以被单独关闭。如果节点显示为很小的点而不是小球可以在显示设置里调整 Node Size。检查集盒和视图范围。有时模型离原点太远显示裁剪后节点不可见。如果只是某个集合的节点不显示确认是否开了按集合过滤显示。不要一上来就重建模型。先确认节点开关再看 Component 显示状态最后才考虑是否误删节点。3. 在 Hypermesh 新界面中搭建多次爆破模型并导出 K 文件3.1 几何清理与边界条件设置岩体几何可以简化成规则长方体。以一个 10 m x 10 m x 10 m 岩块为例把两个炮孔放在靠近中心区域。先清理 CAD 导入时的重复面、小圆角和小缝隙。边界条件方面爆破问题最怕人工边界反射。可以把外侧边界设为无反射边界使用 *BOUNDARY_NON_REFLECTING 模拟无限域如果模型尺寸足够大也可以把边界放到分析关注区很远的地方但这样网格数量太大。生产模型中更推荐在外部边界加无反射边界。底部或周围也可以施加位移约束但要明白约束位置离爆源过近会改变应力波传播导致 PPV 偏高。工程中一般只对模型最外侧约束并在后处理中验证监测点波形是否明显受反射影响。3.2 创建 Part、Section 和单元算法岩体网格在 Hypermesh 中建议先创建六面体网格。如果几何复杂可以先切分几何再用 Map 或 Solid 网格功能生成六面体。炮孔附近的网格需要局部加密。在 LS-DYNA 中岩体单元使用 *SECTION_SOLID算法一般用单点积分即 ELFORM1。它计算快但容易出现沙漏所以需要配合沙漏控制。也可以选择 ELFORM2 全积分计算稳定但成本高。对于爆破算例如果用 RHT 材料并关注损伤单点积分加沙漏控制是常见折中。Section、Part 的对应关系如下对象示例名称作用Componentrock存放岩体网格Sectionsec_rock指定单元算法和积分规则Materialmat_rht存放岩体材料模型Partpart_rock把 Component、Section、Material 关联起来在 Hypermesh 中创建 Part 后最好进入 Card Edit 检查 Part ID 是否唯一。K 文件里 *PART 的 PID 如果重复LS-DYNA 会直接报错。3.3 实现多次爆破载荷多个爆炸事件的时间序列加载多次爆破的载荷实现是“全时域累积”最关键的一步。常见有两种方案。方案一是 LS-DYNA 的 CONWEP 类爆炸载荷即 *LOAD_BLAST_ENHANCED。它可以定义多个爆炸事件每个事件独立设置炸药质量、爆源坐标和起爆时间 TBO。多个事件在同一个求解中依次起爆结构受载面持续收到压力作用岩体材料损伤自然累积。K 文件示意如下具体字段以当前 LS-DYNA 版本手册为准*LOAD_BLAST_ENHANCED $ blastid mass x0 y0 z0 tbo unit 1 30.0 2.0 2.0 2.0 0.0 2 2 30.0 6.0 6.0 2.0 0.05 2这里第一个爆炸在 0 ms 起爆第二个在 50 ms 起爆单位制通过 unit 字段说明。两个不同坐标的爆源可以让应力波在岩体内部发生叠加也能更真实地模拟前后分段爆破。方案二是直接施加压力时程曲线。在 Hypermesh 中创建 *DEFINE_CURVE用时间-压力数据逼近爆炸压力曲线再通过 *LOAD_SEGMENT 或 *LOAD_SURFACE 施加到炮孔壁单元。多个爆破可以通过多条曲线或多个 Segment Set 实现。这个方案更可控但需要自己准备压力曲线工作量大。无论哪种方案都要在 K 文件中确认爆炸载荷的作用面是与岩体单元共用的 Segment Set。作用面丢失或作用面法向方向反了爆炸载荷会完全失效表现为“岩体纹丝不动”。3.4 设置振动监测点、阻尼和求解时间振动监测点通常选在距炮孔不同距离的节点。Hypermesh 中可以在模型上记录节点 ID后续在 K 文件里用 *DATABASE_HISTORY_NODE 输出这些节点的速度、加速度和位移。一个典型的监测点配置如下*DATABASE_HISTORY_NODE $ id1 id2 id3 1001 1002 1003 *DATABASE_NODOUT $ DT 1.0e-4*DATABASE_HISTORY_NODE负责把节点加入历史输出*DATABASE_NODOUT控制输出时间间隔。这里 DT 取 0.1 ms对于爆破振动几十到几百 Hz 的频率足够。阻尼方面LS-DYNA 中可以用 *DAMPING_GLOBAL 施加质量阻尼但要谨慎。阻尼过大会把高频振动压平PPV 明显降低阻尼过小则模型长时间余振浪费计算时间。可以先用无阻尼模型做短时分析观察波形衰减情况再决定阻尼参数。求解时间用 *CONTROL_TERMINATION 设置。多次爆破模型如果二次起爆在 50 ms至少要计算到 150 ms 才能看到完整波动衰减。时间步长通过 *CONTROL_TIMESTEP 控制必要时使用质量缩放但质量缩放比例不能太大。3.5 导出 K 文件后要重点检查的卡片段在 Hypermesh 中导出 K 文件后不要直接提交 LS-DYNA。先用文本编辑器打开 K 文件按顺序检查几个关键卡片段。需要检查的内容包括*PART 是否包含全部岩体和受载面*MAT_RHT 参数的单位是否统一*SECTION_SOLID 的 ELFORM 是否匹配*LOAD_BLAST_ENHANCED 的爆源坐标、时间和炸量是否合理*DATABASE_HISTORY_NODE 中的节点 ID 是否真实存在*BOUNDARY_NON_REFLECTING 是否施加在外表面*CONTROL_TERMINATION 结束时间是否覆盖所有爆破事件。检查一个容易忽略的点如果把无限域边界和爆炸载荷同时施加在同一组外表面单元上可能导致爆炸压力被边界吸收计算出来的振动偏小。生产模型中要注意受载面和吸收边界不要重叠。4. LS-DYNA 求解与后处理如何验证全时域振动和损伤累积4.1 命令行提交求解与日志关注点K 文件准备好后用命令行提交求解ls-dyna iblasting_multi.k ncpu8 memory200m如果使用的是 Ls-Dyna 的 mpp 版本命令格式可能不同。提交后观察 message 或 d3hsp 文件重点看材料 ID 是否重复最小时间步是否太小是否存在负体积警告沙漏能是否过大是否出现单元侵蚀侵蚀比例是否合理。如果模型几分钟内内存暴涨可能是 K 文件里出现了未建模的接触或过多的输出频率。输出频率过密会让 d3plot 文件很大但不一定会导致计算错误。4.2 从二进制结果中读取振动速度时程计算完成后用 LS-PrePost 打开 d3plot再读 nodout 文件得到监测点速度时程。具体做法LS-PrePost 中打开 d3plot在 History 或 NODOUT 中加载 nodout选择监测节点 ID提取 X、Y、Z 方向速度绘制速度-时间曲线找到每个爆破事件后的 PPV把多次爆破中的每个峰值按时间排列形成整个时域的振动包络。读取时要注意速度方向要参考实际工程。评价爆破振动时通常用三个方向的合成速度不能只取一个方向的最大值。4.3 读取损伤变量并绘制时程曲线RHT 模型在计算中会更新损伤变量。LS-PrePost 中可以查看单元历史变量具体编号取决于 LS-DYNA 版本对 RHT 模型的输出定义。比较稳妥的流程是在 d3plot 后处理中查看 Element History选择岩体 Part逐个查看历史变量找到随爆破明显增大的变量如果历史变量编号不明确可以在 *DATABASE_EXTENT_BINARY 中打开额外历史变量输出根据主应力、塑性应变、压力变化综合判断损伤分布。这里要避免的误区是不要只盯一个变量。RHT 模型的损伤可能体现在多个历史变量中最终判断要结合岩体表面上是否出现拉伸破坏、监测点振动波形是否出现突然衰减、单元是否删除等证据。4.4 结果合理性判断单次爆破和累积规律一个合理的多次爆破全时域结果通常会出现以下特征第一次爆破后靠近爆源的单元开始出现损伤振动波形在一定时间内衰减第二次爆破起爆时如果时间间隔足够短监测点的振动不是从零开始而是叠加在第一次爆破的余振上PPV 不一定单调增大因为应力波叠加有相位问题但远离爆源的损伤会逐步向更远区域扩展单元删除或侵蚀后应力波传播路径会改变后续爆破的振动场和新鲜岩体不同。如果计算结果显示每次爆破都完全独立第二次起爆时波形从零开始说明载荷时序或模型状态重置出了问题。如果第二次爆破后整个岩体瞬间炸飞通常是材料参数强度过低、侵蚀阈值过大或载荷单位错误。5. 高频问题与排查路径从 Hypermesh 到 LS-DYNA 的常见坑5.1 材料量级错误导致的飞散或发散现象计算一开始岩体节点速度远超正常波速单元迅速飞散或者出现负体积。原因通常是单位制换算错误。密度、弹性模量、应力强度三个量中只要有一个量级不对材料波速就会异常。例如把 50e9 Pa 填成 50e6岩体刚度和波速会低很多。排查顺序检查密度和弹性模量是否处于同一单位制下计算岩体纵波波速对比理论值检查压力载荷的单位是 Pa、MPa 还是 GPa用小模型单爆破测试先不分析损伤只看波形是否合理。5.2 爆炸载荷没按设定时间触发现象模型前一段没有任何振动某个时刻突然出现压力或者第二次爆破根本没发生。原因*LOAD_BLAST_ENHANCED 的 TBO 设置错误或者多个爆炸卡片没有同时被求解器识别作用面丢失。排查路径在 K 文件中确认每个 blast 的 TBO 和坐标确认作用 Segment Set 与岩体单元面重合在 message 中搜索 blast 相关输出看看有几个爆炸事件被读取用单爆源去掉时间间隔先验证压力能否正确施加。5.3 负体积、网格畸变和单元侵蚀现象LS-DYNA 在某个时间步报出负体积位置通常在炮孔附近。原因炮孔附近单元畸变材料失效参数太激进或者网格尺寸与爆炸压力不匹配。处理建议细化炮孔附近网格但不要只缩尺寸还要重新正交化调整 *MAT_ADD_EROSION 中的失效准则避免单元过早删除减少时间步长取消过大的质量缩放检查沙漏能若沙漏能占比超过内能 5% 到 10%需要加强沙漏控制。5.4 节点不显示、K 文件导入异常等操作问题在 Hypermesh 新界面中节点不显示往往是显示开关问题。先检查 View 或 Display 菜单里的节点开关如果 K 文件导入 Hypermesh 后很多节点消失则可能是模型的容差设置过小导入时节点被合并或剔除。K 文件导入异常的另一个常见原因是 LS-DYNA 关键字版本与当前 Hypermesh 模板不匹配。建议用文本编辑器打开 K 文件找到报错的关键字在 Hypermesh 中确认是否支持该版本。5.5 多次爆破应该用单次多事件还是重启动如果两次爆破间隔很长比如第一次爆完 1 秒后才第二次起爆单次求解需要算 1 秒空窗期成本很高。这时可以考虑重启动或者采用分段求解再映射应力。但重启动的缺点是操作复杂而且材料损伤、历史变量在重启动文件中是否正确传递需要验证。对大多数研究性算例只要间隔在几十到几百毫秒量级用单次多事件的 *LOAD_BLAST_ENHANCED 更稳定能自然实现全时域累积。如果间隔真的很长则采用以下思路第一段算到振动基本结束后输出重启动文件第二段模型保留第一段的应力、损伤、单元状态第二段只从第二次起爆时间开始计算。这个方案要额外验证重启动后的材料状态没有丢失否则累积就不成立。6. 生产级模型落地建议与扩展方向6.1 提交求解前检查清单整理一份可以反复用的检查清单避免每次都在相同位置踩坑单位制是否统一密度、波速、压力是否按同一组量纲换算岩体材料参数是否来自试验或可靠文献是否做过材料敏感性测试网格质量是否通过检查炮孔附近最小单元尺寸是否满足压力曲线分辨率Part、Section、Material 是否正确关联有无重复 ID爆炸载荷是否包含全部爆源坐标、炸量、起爆时间是否正确作用 Segment Set 是否与岩体单元面一致法向是否指向岩体内部监测点节点 ID 是否存在是否能在 nodout 中输出无反射边界是否与受载面重叠求解结束时间是否覆盖全部爆破事件输出间隔是否合理沙漏能、质量缩放、负体积需要在计算刚开始时就盯住。6.2 参数校核、网格敏感性和阻尼设置生产级分析至少要做三件事。第一材料参数校核。用单单元或简单梁、柱试件验证材料模型在单轴拉伸、单轴压缩、围压下的响应是否符合试验。不能直接拿别人的参数套进自己的单位制。第二网格敏感性分析。粗网格和细网格各算一次对比监测点 PPV 和损伤范围。如果网格尺寸变化后结果差异超过 10%网格还需要继续加密或改变单元算法。第三阻尼和边界敏感性。无反射边界和阻尼对振动时程影响很大。建议分别计算无阻尼、低阻尼、高阻尼三组观察 PPV 和波形包络的变化选择与现场实测接近的一组。6.3 从均匀岩体到节理岩体的扩展思路本文案例把岩体处理为连续介质适用于完整或弱节理岩体。实际工程中岩体通常包含节理、层理和裂隙。要更真实地模拟多次爆破后节理张开、滑移和累积损伤可以在 Hypermesh 中建立离散节理面或使用材料模型中的初始损伤场。扩展方向包括在岩体单元中引入 Weibull 分布让材料强度具有随机性用 *MAT_JHC 或 *MAT_HJC 分析高围压冲击问题把现场实测振动波形反演成炮孔压力曲线替换 CONWEP 爆炸载荷联合模型降阶方法降低多次长间隔爆破的计算成本在后处理中增加 PPV 衰减曲线与萨道夫斯基公式对比标定数值模型。多次爆破全时域振动损伤累积并不是一个高不可攀的分析但它很依赖基本功。单位制、网格质量、材料参数、载荷时序、损伤输出每一环都会决定最终结果是否可信。对新手来说建议先从两次爆破、两个监测点的小模型开始把“第二段爆破是否在上一段损伤基础上继续”的物理过程看明白再扩展到更复杂的地质条件和爆破方案。