LS-DYNA壳单元冲击仿真中的5个关键设置：以方管碰撞为例（含沙漏控制参数）

LS-DYNA壳单元冲击仿真中的5个关键设置：以方管碰撞为例

最近在帮一个做汽车防撞梁设计的团队排查仿真问题，他们反复计算一个方管的侧向冲击，结果总是不稳定，要么能量不守恒，要么变形模式看起来“怪怪的”。折腾了几轮之后，我们发现，问题往往不是出在模型有多复杂，而是几个看似基础的参数设置上。LS-DYNA这类显式动力学软件，就像一个对初始条件极其敏感的物理实验，你给的“实验装置”参数稍有偏差，出来的“现象”就可能天差地别。尤其是对于方管、薄壁梁这类大量使用壳单元的冲击仿真，计算精度和效率的平衡点，就藏在几个关键设置里。今天，我们就抛开泛泛而谈，深入聊聊在方管冲击仿真中，直接影响结果可信度的五个核心设置点。这些内容适合已经上手LS-DYNA，但在处理具体工程问题时，希望知其然更知其所以然的中级用户。

1. 沙漏控制：为何IHQ=4成为方管冲击的“默认推荐”？

沙漏现象是显式积分中一个老生常谈但又无法回避的问题。简单来说，它是因为单点积分（为了计算效率）导致的一种零能变形模式，单元像沙漏一样扭曲，却不产生任何应变能。在方管受到冲击时，大变形区域（如撞击点、褶皱产生处）的单元极易激发这种非物理的变形，严重时会导致结果完全失真。

很多教程会告诉你，在*HOURGLASS卡片里把IHQ设为4。但为什么是4？它背后的逻辑是什么？

IHQ这个参数定义了沙漏控制的算法类型。LS-DYNA提供了多种选项，从1到8各有侧重。IHQ=4对应的是Flanagan-Belytschko粘性沙漏控制。它的核心思想是引入与沙漏变形速率成正比的粘性阻尼力来抑制这种模式。对于像方管冲击这类涉及金属大塑性变形、材料有显著应变率效应的仿真，IHQ=4表现出几个独特的优势：

• 与材料模型的兼容性好：它能够较好地配合弹塑性材料模型（如MAT_24），在塑性流动阶段依然提供稳定的控制力。
• 对变形模式的适应性：在方管压溃过程中，既有弯曲褶皱，也有局部拉伸，IHQ=4算法对此类复杂变形场的适应性较强。
• 计算开销相对合理：相比一些更复杂的算法（如IHQ=5或IHQ=6的刚度型控制），它的计算成本增加不多，在精度和效率间取得了不错的平衡。

但这绝不意味着可以无脑设置IHQ=4。另一个关键参数是沙漏系数QH，它控制着阻尼力的大小。默认值（对于实体单元是0.1，对于壳单元通常是0.03）是一个安全的起点，但并非最优。

注意：过大的QH值会过度抑制沙漏能，但同时也可能“僵化”结构的真实变形，特别是影响褶皱的形态；过小的QH则可能无法有效控制沙漏。一个实用的技巧是，在完成初步计算后，查看GLSTAT文件中的沙漏能（Hourglass Energy）与内能（Internal Energy）的比值。通常要求这个比值小于5%，对于精度要求高的分析，最好能控制在1%-2%以内。如果比值过高，可以尝试微调QH，比如从0.03调整到0.05或0.01，观察变形模式和能量曲线是否改善。

下表对比了在方管冲击仿真中几种常见IHQ设置的潜在影响：