4. 单元类型与网格策略

在 LS-DYNA 显式动力学分析中，单元类型与网格策略直接决定计算规模、稳定时间步长、接触搜索效率以及最终结果的物理可信度。显式求解的核心约束是条件稳定：时间步长必须小于临界值，而临界值又由最小单元尺寸、材料波速和积分方式共同决定。因此，单元选择不是“越精细越好”，而是在几何特征、变形模式、计算成本和精度之间做工程权衡。

一、四类单元的适用边界

实体单元（*SECTION_SOLID）适用于三维连续体，如碰撞吸能盒、橡胶衬套、混凝土试件、锻造坯料。显式分析中常用 ELFORM=1（常应力单点积分六面体）和 ELFORM=2（单点积分四面体）。单点积分计算快，但会产生沙漏模式，必须配合 *HOURGLASS 控制；ELFORM=-1 或 -2 为全积分，无沙漏但更易发生剪切自锁，且临界步长更小。工程上优先用六面体映射网格，四面体仅用于过渡区或复杂几何填充。

壳单元（*SECTION_SHELL）用于薄壁结构，如车身覆盖件、管件、电池包上盖。ELFORM=2 是 Belytschko-Tsay 单点积分壳，效率最高，是碰撞分析默认选择；ELFORM=16 为全积分壳，适合需要精确弯曲应力的局部区域。壳单元必须定义厚度，并注意沿厚度方向的积分点数目（NIP），通常 3～5 个积分点足以捕捉弯曲。

梁单元（*SECTION_BEAM）用于细长构件，如防滚架、螺栓、加强筋。ELFORM=1 为 Hughes-Liu 梁，ELFORM=2 为 Belytschko-Schwer 梁，后者效率更高。梁单元需要定义截面形状和方向节点，否则会因扭转方向不确定导致结果异常。

离散单元（*SECTION_DISCRETE）用于弹簧、阻尼器、安全带滑环、点焊连接等零维或一维连接。它不参与网格划分，而是通过节点对定义力—位移或力—速度关系，常用于简化局部连接，避免为螺栓或焊点建立精细实体网格。

二、积分算法与沙漏控制

显式求解普遍采用单点积分以降低计算量，但单点积分单元只有一个积分点，无法感知某些变形模式，从而产生零能模式，即沙漏。沙漏表现为网格出现锯齿状、交替变形的非物理振荡，能量曲线中沙漏能异常增长。判断标准是：沙漏能应小于总内能的 5%，超过 10% 则结果不可信。

控制沙漏的常用手段包括：在 *HOURGLASS 中设置 IHQ=4（Flanagan-Belytschko 刚度型），并配合 QM=0.05～0.1 的系数；对关键区域改用全积分单元；细化网格使沙漏模式难以形成。需要注意的是，沙漏控制系数过大会人为增加刚度，导致响应偏硬，因此应在沙漏能与物理真实性之间折中。

三、网格密度与质量的实际影响

网格密度首先影响时间步长。以钢为例，E=210 GPa，ρ=7850 kg/m³，波速 c≈5170 m/s。若最小单元边长为 2 mm，则 Δt≈0.39 μs；若局部细化到 0.5 mm，Δt 降至约 0.097 μs，计算步数增加约 4 倍。因此，局部细化必须评估其对整体计算成本的影响，必要时用质量缩放将小单元的时间步长提升到可接受范围，但需监控质量增加比例，通常控制在 1% 以内。

网格质量同样关键。翘曲壳单元会高估弯曲刚度，长宽比过大的单元会降低波传播精度，负雅可比单元直接导致计算终止。工程实践中，壳单元长宽比宜小于 5，翘曲角宜小于 15°，实体单元雅可比宜大于 0.6。对于接触区域，主从面网格密度应匹配，避免从节点穿透主面。

四、选择依据与工程决策

实际建模时，可遵循以下判断顺序：先看几何特征——一个方向尺寸远小于另外两个方向时用壳，两个方向都小时用梁，三个方向相近时用实体；再看变形模式——以弯曲为主用壳，以压缩/剪切为主用实体；最后看连接与简化需求——螺栓、焊点、弹簧可用离散单元替代。对于混合结构，实体与壳之间可用 *CONSTRAINED_SHELL_TO_SOLID 或共节点方式连接，但需注意壳的转动自由度与实体平动自由度之间的不匹配。

总之，单元类型与网格策略是显式动力学建模中“牵一发而动全身”的环节。合理的做法是：在概念设计阶段用较粗网格快速验证趋势，在详细设计阶段对关键区域局部细化，并始终监控沙漏能、质量增加比例和能量平衡，确保计算结果既高效又可信。