LS-DYNA壳单元冲击仿真5个关键设置：方管碰撞沙漏与积分点全解析

跑过方管碰撞仿真的人大概率都见过这么一幕模型提交到LS-DYNA 里跑了大半天后处理一翻沙漏能占比飙到 12% 以上变形模式从规则的折叠压溃变成了歪歪扭扭的锯齿状峰值力曲线毛刺一大堆。这时候有人直接去调沙漏控制参数有人回头换单元公式但真正的问题往往在一开始建模型时的几个关键设置上。这篇文章要聊的就是用 LS-DYNA 做壳单元冲击仿真时必须要盯住的 5 个关键设置全部以方管碰撞为例展开把沙漏控制参数单独拎出来讲透顺便把“壳单元厚度与积分点个数到底什么关系”这个问题一起解决掉。方管碰撞是个特别适合当教材的案例。方管是工程里最基础的薄壁吸能结构碰撞过程中弯曲、压溃、折叠、甚至局部撕裂全都会出现几乎覆盖了壳单元冲击仿真能遇到的所有典型问题。把这 5 个设置搞明白再去做 B 柱、防撞梁、电池包护板这类真实结构思路完全可以迁移。全文会给参数推荐取值和背后的逻辑也会说明哪些情况下默认值就够用哪些情况必须动手干预。1. 为什么拿方管碰撞当案例一个设置不当就“翻车”的典型场景1.1 方管碰撞里有哪些“壳单元杀手”先说结论方管碰撞看着简单实际上比很多车身钣金件更挑设置。原因有三点。第一方管压溃是弯曲主导的问题。碰撞时壁板发生大幅弯曲塑性铰沿着管壁依次形成折叠模式是不是合理直接取决于壳单元能不能准确传递弯矩。一旦单元公式或者积分点个数不到位弯矩-曲率关系就会失真压溃力峰值和折叠波长全都会跑偏。第二方管存在复杂的自接触。折叠的管壁会层层叠在一起不同壁面之间发生摩擦和挤压这意味着接触算法必须稳定否则节点穿透、初始穿透这两个大坑会接踵而来。节点穿透会导致局部应力集中初始穿透又会让力曲线上来就是一根尖峰怎么看怎么假。第三方管压溃后单元扭曲严重。壳单元在经历大转动和大变形之后原来的平面可能被扭曲成马鞍形状如果翘曲刚度控制不好网格一旦扭曲就会“锁死”或者“软化”仿真结果自然不可信。1.2 这个案例的迁移价值方管碰撞模型不算大几千到几万个单元就能跑出像样的结果没有整车模型那么长的建模周期很适合用来摸清各种参数对结果的影响。我平时做防撞梁、副车架这类结构时遇到拿不准的壳单元设置都是先拿一个方管或者帽型梁模型做验证趋势对了再上真实结构这个方法强烈推荐。如果你发现自己的方管碰撞仿真出现了下面任何一种现象这篇博文里的 5 个设置大概率能帮你定位问题压溃力峰值异常偏高曲线毛刺多变形模式是歪的该对称的不对称沙漏能占比一直压不下来超过内能的 5%计算中途报错或者结果对网格尺寸极其敏感。2. 关键设置一壳单元公式怎么选BT 还是全积分2.1 ELFORM2 和 ELFORM16 的本质区别LS-DYNA 的壳单元公式很多日常冲击仿真里真正需要关注的其实就是两类默认的 Belytschko-TsayELFORM2和全积分壳ELFORM16。它们之间的区别用大白话说就是“计算一个单元截面内力时采样几个点”。Belytschko-Tsay 是单点积分也就是在单元中心取一个积分点。它的核心优势是快而且是又快又稳尤其适合显式动力学里动辄百万单元的碰撞模型。但单点积分有个天然缺陷单元会产生一种“变形了却不产生应变能”的变形模式也就是沙漏。这个后文会重点展开。全积分壳在每个单元里布置多个积分点理论上不会出现沙漏计算精度看似更高但代价也摆在那里计算量成倍增加而且在弯曲主导的问题中容易出现剪切锁死就是单元过于刚硬导致结构整体刚度偏大。锁死会让压溃力虚高需要加额外的假设应变方法来修正。2.2 方管碰撞选型逻辑方管碰撞这类问题我的推荐是默认用 ELFORM2同时把沙漏控制做好。理由很直接——方管碰撞模型少则几千单元多则几十万单元显式分析每增加一点点单元开销都会反映到实际计算时长上。BT 单元的计算效率高配合合理的沙漏控制参数精度完全够用。只有在下面两种情况下才建议考虑局部换成 ELFORM16管壁扭曲特别严重沙漏控制参数已经调到合理范围沙漏能仍然压不住你关心的是某一小区域的精细应力分布比如焊点附近的疲劳热点分析需要更准确的局部应力场。实际工程中我会做的是整个模型保持 ELFORM2把可能发生严重折叠的局部区域单独设置成 ELFORM16同时注意两种单元之间的连接和过渡。这个做法能兼顾效率与精度但也提醒一句——局部切换后一定要检查过渡区的应力是否连续。2.3 壳单元公式选择的参考对照单元公式积分方式沙漏风险弯曲精度计算成本适用场景ELFORM2Belytschko-Tsay单点积分较高需沙漏控制一般依赖网格质量低大模型、冲击、碰撞ELFORM16全积分多点积分基本无较高高局部精细区、压力集中区ELFORM8/20多点积分低高更高特殊工况使用较少需要特别提醒的是ELFORM16 不是万能的。有的工程师一看到沙漏就把全部单元改成全积分结果计算时间翻倍压溃力反而偏高最后发现是剪切锁死。全积分壳需要配合 LS-DYNA 的假设应变选项来缓解锁死不能无脑使用。3. 关键设置二厚度方向积分点个数2 个还是 5 个3.1 壳单元厚度方向积分点到底在算什么这是很多新手最容易忽略的一个参数。壳单元虽然在几何上是二维的但它要代表一个带厚度的三维板壳。厚度方向上的应力分布就是通过积分点来还原的。默认情况下LS-DYNA 的 BT 壳单元厚度方向只取 2 个积分点位置大约在厚度上下表面附近这只能还原一种线性的应力分布。当结构以纯弯曲为主时线性应力分布勉强够用。但方管碰撞不是纯弯曲——压溃过程中壁板一侧受拉一侧受压同时还有膜力叠加。中性层位置会随着塑性流动发生移动应力分布会变成明显的非线性。这时候 2 个积分点就像拿一把两齿的叉子去叉牛排能叉住但不稳结果要么是弯矩给的偏保守要么是塑性变形推进速度不对。3.2 压溃过程对应力分布的要求方管折叠时每一道塑性铰都是“先外侧屈服、再向内侧扩展”的过程。外侧纤维应变很大内侧纤维应变小甚至为负中性层一直在移动。厚度方向上至少要 3 到 5 个积分点才能捕捉到这个渐变过程。我做过一个对比测试同一个方管模型分别设置 NIP2、3、5其余条件完全一致。结果是 2 个积分点的模型最先出现压溃形态偏移折叠间距明显偏大峰值力比 5 积分点模型高了大约 8%5 积分点模型的变形模式对称性更好力-位移曲线也更平滑。这个差距在做整车碰撞安全件分析时是不可接受的。3.3 积分点个数与厚度的关系这里正好回应一下热搜里那个问题“ansys 壳单元厚度与积分点个数关系”。严格来说积分点个数与物理厚度本身不是绑定关系——厚度多大都可以只取 2 个积分点。但有效应力和应变的精度与积分点个数相关薄壳在弯曲时上下表面距离近应力梯度不算极端厚壳或者多层复合材料壳厚度方向应力变化剧烈积分点不足就会明显丢失精度。工程上推荐这样做普通金属薄壁件碰撞仿真NIP3 起步方管这类以弯曲压溃为主的结构NIP5费用增加不多涉及复合材料层合板每一层至少 2-3 个积分点。NIP 从 2 提高到 5计算量会增加 20%-30% 左右相对接触和总体的计算开销来说完全可以接受。为了省那 30% 的单元计算时间去冒精度风险非常不划算。4. 关键设置三网格尺寸与翘曲刚度控制4.1 网格密度如何改变折叠模式网格尺寸对碰撞仿真的影响远远超出很多人的直觉。方管碰撞中折叠波长的形成和网格密度强相关。网格太粗的时候模型没有足够的自由度去表达塑性铰的连续弯曲变形会被迫集中在个别单元上压溃力虚高折叠模式会变成“多边形的棱”而不是光滑的弧线。对于方管这类结构我建议沿管壁宽度方向至少划分 6 到 10 个单元单元长宽比尽量控制在 3:1 以内最好接近正方形。比如一根 50mm×50mm 的方管用 5mm 左右的单元尺寸就差不多如果重点关注压溃模式的细节进一步加密到 3mm 也有必要。网格大小同时决定了显式分析的最小时间步长这个后面会再提到。做个简单测试就能理解网格敏感性的“可怕”把网格从 10mm 加密到 5mm峰值力可能下降 10%-15%变形模式也可能完全不同。如果两个网格密度下结果差异很大先别急着怀疑参数大概率是网格还没收敛。碰撞仿真至少要做一个“网格加密一倍结果基本不变”的验证。4.2 翘曲刚度 warping stiffness 的含义与调整壳单元在发生大变形后原本的平面会变成空间曲面这时单元不能再用简单的平面应力公式来描述了。翘曲刚度warping stiffness就是描述单元抵抗这种“离开平面”变形能力的一个量。LS-DYNA 中通过 *CONTROL_SHELL 里的 WARPP 参数控制默认值是 1.0 左右。这个参数一般在网格质量好的时候几乎感知不到它的“存在感”恰恰在网格质量差的时候刷满。网格一旦被压溃挤歪翘曲刚度会被过度激活导致整个单元的刚度被人为抬高结果就是压溃力虚高。反过来如果网格一直保持规则的四边形翘曲刚度对结果的影响就很小。所以处理翘曲刚度的正确方式是先改善网格形态而不是先调 WARPP。尝试把 WARPP 在 0.5 到 1.5 之间做敏感性测试如果结果对 WARPP 异常敏感说明网格质量不行加密网格才是正解。4.3 壳单元厚度建模的常见错误壳单元厚度与积分点个数这个热搜词背后隐藏着一个更基础的问题——很多人建壳模型的时候厚度定义就是随手一填。壳单元的厚度通过 *SECTION_SHELL 的 T 字段定义如果厚度方向和实际结构不一致或者把同一个厚度分配给整个零件局部刚度就错了。方管建模里常见的错误有三个把厚度填成 0导致膜刚度和弯曲刚度全错忽略了方管弯角处的壁厚变化真实管材在弯角处厚度往往略小惯性释放时用集中质量代替真实质量分布导致冲击响应偏差。这些错误不会报错但是结果就是“看起来差不多细节全不对”。建议每次建模都单独检查一遍零件厚度而不是依赖默认值。5. 关键设置四沙漏控制参数方管碰撞最容易踩的坑5.1 沙漏的本质单点积分单元的零能模式沙漏这个词听起来玄乎原理其实不难理解。单点积分单元只在中心取一个积分点当单元发生某种“交替变形”时——比如顶角往相反方向翘——所有积分点处的应变都恰好为零于是单元变形了却不产生任何应变能。这种没有刚度抵抗的变形模式就是零能模式因为它的变形云图看起来像一个沙漏所以叫沙漏。方管碰撞最容易激发沙漏因为压溃过程中单元要经历剧烈的弯曲和扭曲非常容易触发这种“不花力气”的变形模式。沙漏一旦蔓延开来单元彼此之间就像没有连接一样结果就是结构看起来变形了但实际上没有承受力力的峰值会严重失真。判断沙漏有三个方法任何一个都可以用看 GLSTAT 文件里的沙漏能与内能比值控制在 5% 以内算合格严格一些可以按 1%-2% 控制看后处理中 *DATABASE_BINARY_D3PLOT 的沙漏能云图如果沙漏能集中的区域和主要塑性铰重合说明该区域单元需要处理看力-位移曲线如果出现异常锯齿或者力突然下降也要回头检查沙漏。5.2 沙漏控制卡片参数详解LS-DYNA 的沙漏控制参数集中在 *HOURGLASS 和 *CONTROL_HOURGLASS 里。这里先把核心参数拆开讲。IHQ 是沙漏控制类型对壳单元来说常用的有几种IHQ类型特点1标准粘性沙漏控制默认值稳定性好但刚度型问题误差可能偏大2Flanagan-Belytschko 粘性适合高速冲击耗散沙漏能3Flanagan-Belytschko 刚度适合低频弯曲但过大容易刚化4刚度粘性组合兼顾两种机制通用性好6Belytschko-Bindeman 假设应变对壳单元冲击仿真最推荐的组合沙漏控制精度高QM 是沙漏控制系数取值范围通常在 0.03 到 0.15 之间默认 0.05 或 0.1。QM 越大抵抗沙漏的刚度/粘性越强但过大也会让结构整体变刚压溃力虚高。QB 是壳单元弯曲刚度放大系数典型值为 1.5它主要用于提高弯曲主导问题中壳单元的沙漏抵抗能力。QH 是粘性系数常用于高速冲击取 0.1 左右。如果要在 *HOURGLASS 卡片里写全参数可以这样参考某种流行的配置*HOURGLASS $ HGIHQ IHQ QM QB QH QI QL QM 1 6 0.05 1.5 0.10 2.0 2.0 1.5这里 HGIHQ 表示沙漏控制的组编号后续通过 *SECTION_SHELL 的 HGID 字段把特定的沙漏控制类型应用到对应的壳单元上。比如部分单元用 IH6部分单元用 IH这样可以分区控制。5.3 方管碰撞实战调参流程我对方管碰撞的沙漏控制套路是这样的先跑一个什么都不调的基线模型看两个数据沙漏能占比和变形模式。如果沙漏能占比在 5% 以内、压溃模式正常那就不用动沙漏参数默认 IHQ1 就够了。如果沙漏能超标第一步不是加系数而是检查网格。网格太粗、长宽比太大是沙漏激发的首要原因。把关键区域网格加密后沙漏能往往能直接降下来一大截。如果网格没问题再动手换控制类型。优先试 IHQ6配合 QM0.05、QB1.5。原因很简单IHQ6 属于假设应变法它通过应变场投影来消除零能模式在抑制沙漏的同时对结构刚度的附加影响最小非常适合壳单元冲击仿真。对比 IHQ1IHQ6 在弯曲主导问题中能有更好的表现。如果 IHQ6 之后沙漏能仍然偏高再适度提高 QM 到 0.08 或者 0.1同时观察压溃力变化。如果峰值力明显抬高说明控制过强了需要回退。现场经验是沙漏能偏高的模型十有七八是网格质量的问题只有两三成需要动沙漏参数。一句话总结——先用网格治病再用参数调整最后才能用全积分救场。6. 关键设置五接触定义与时间步长缩放6.1 自接触与摩擦参数方管压溃时管壁会一层层叠在一起自己和自己发生接触这属于自接触问题。用 *CONTACT_AUTOMATIC_SINGLE_SURFACE 是标准选择它能自动识别所有外表面处理壳单元的自接触时比较省心。接触参数里最容易影响结果的三个值是静摩擦系数 FS、动摩擦系数 FD 和接触粘性阻尼 VDC。金属方管碰撞时推荐 FS 取 0.15 到 0.2FD 取 0.1 左右VDC 取 10 到 20。摩擦系数影响接触面上的能量耗散FS 太大压溃力会偏高FS 太小管壁之间会过于顺滑折叠模式看起来不自然。接触算法里还有一个 SOFT 参数值得重视。SOFT0 是默认的罚函数法SOFT1 或 SOFT2 用于处理初始穿透和接触稳定性问题。方管压溃时单元会发生大翻转网格节点可能瞬时穿透到对面壁内开启 SOFT2segment-based contact能明显提升接触稳定性。实际做下来SOFT2 配合单面接触可以有效减少尖峰力和初始穿透引起的伪振荡。6.2 时间步长缩放与质量增加LS-DYNA 是显式积分求解器时间步长受最小单元尺寸控制。理论上要满足 Courant 稳定性条件时间步长不能超过声波穿过最小单元所需的时间。实际计算中我们很少手动去算声速直接看 *CONTROL_TIMESTEP 里的 DT2MS 和 TSSFAC。TSSFAC 是时间步缩放系数默认 0.9。冲击仿真里我通常建议压到 0.6 到 0.7。原因很简单壳单元在碰撞过程中会扭曲扭曲后的实际临界时间步要小于初始值默认系数偏高时容易在变形后期出现数值不稳定。如果时间步太小导致计算太慢第一时间想的不是降低 TSSFAC而是去查最小单元在哪、为什么那么小。把局部细小网格通过网格划分策略优化掉往往能把整体时间步立刻提上去。另一个必须盯住的指标是质量增加。LS-DYNA 在时间步不足时会自动在节点上增加质量来维持稳定体现在 GLSTAT 文件里的增加质量占比。业界通行规则是控制在 5% 以内超过这个值说明模型里存在过小单元或者不合理的时间步缩放需要返工。质量增加超过 10% 的结果基本不可信。6.3 冲击问题的时间稳定性一点类比做电路仿真的同事有时候会跟我聊 LTspice 里的冲击电流仿真——步长设大了电流峰值会被平滑掉步长设小了运算时间又激增。LS-DYNA 冲击仿真面对的问题本质上是一样的时间步长是精度和效率的平衡点。结构冲击过程里的峰值力对应着应力波传播的瞬间步长过大就会把这个瞬间抹平。虽然两个软件解决的物理问题完全不同但这套“数值步长”的逻辑是共通的。所以当你发现力-位移曲线明显失真时除了看接触和沙漏也要回头确认一下时间步长设置。方管这类小模型时间步长即使压紧一点也花不了多少时间但到了整车级大模型这个取舍就得精打细算了。回到开头那个问题如果你的方管碰撞仿真出现了沙漏能超标、压溃模式不对、尖峰力乱跳不要急着把责任推给 LS-DYNA更不要盲目堆沙漏刚度。按这个顺序排查网格质量是否足够、单元公式是否匹配、厚度方向积分点是否太少、翘曲刚度是否有异常、沙漏控制参数是否合理、接触定义是否稳定、时间步长是否太大。大多数情况下问题出在最容易被忽略的那两三个设置上。我做这些方管碰撞算例最大的体会是显式动力学仿真的参数不是越多越好而是每个参数都得知道它在控制什么。壳单元公式、积分点个数、网格尺寸、沙漏控制、接触与时间步长这 5 个设置像一根链条任何一环出了问题最终都会在力-位移曲线和变形模式上暴露出来。建议你拿手头任何一个方管碰撞模型把网格加密一档、把 NIP 改成 5、把沙漏控制切到 IHQ6 各跑一遍对比一下峰值力和压溃模式的变化直观感受会比任何理论都深刻。