4. 网格划分规范与质量准则

整车碰撞模型的网格质量，直接决定显式积分能否稳定推进、接触能否正确建立、以及加速度与侵入量曲线是否可信。一个常见的误区是先把几何“划满”，再回头修质量；实际上，单元类型选型与尺寸策略应在划分前确定，质量指标则作为划分过程中的实时判据，而不是事后补救。

一、单元类型选型：先看受力与厚度，再看几何

壳单元适用于其中一个方向尺寸远小于另外两个方向的薄壁结构，典型对象是车身覆盖件、纵梁、门槛梁、B柱、地板、前后防撞梁等钣金件。整车碰撞中，壳单元承担了绝大部分吸能与传力路径，因此是主力单元类型。选用壳单元时需明确厚度方向积分方式：对于需要模拟弯曲与面内行为的钣金，通常采用沿厚度方向多个积分点（如 NIP=5）的壳公式，以较少的计算代价获得较合理的弯曲刚度；若厚度方向应力梯度显著或需要更精确的塑性弯曲，可考虑增加积分点或改用实体单元。

实体单元适用于三个方向尺寸相近的块体，如发动机缸体、变速器壳体、悬置支架、铸铝减震塔、电池包内部模组等。碰撞中实体单元常处于三向受压或复杂应力状态，采用单点积分实体单元时需注意沙漏控制，通常配合 *HOURGLASS 中的刚度型或粘性型控制，并检查沙漏能占总内能比例，一般建议控制在 5% 以内。若实体单元用于薄壁件，则会出现厚度方向单元数不足、弯曲刚度过高的问题，此时应改回壳单元。

梁单元适用于一个方向尺寸远大于另外两个方向、且以轴向力、剪力与弯矩为主的结构，如副车架管件、防滚架、座椅骨架管、部分加强撑杆。梁单元的关键在于截面定义与积分规则：需通过 *SECTION_BEAM 指定截面形状与尺寸，并选择与受力匹配的积分方式（如管件常用 ELFORM=2 或 ELFORM=3）。梁单元与壳、实体连接处应避免单点铰接，必要时用刚性单元或短梁过渡，防止局部应力集中导致非物理失效。

二、整车网格尺寸策略：按吸能贡献分区控制

整车碰撞模型不可能也不需要对所有区域采用同一尺寸。合理的做法是按“吸能贡献 + 接触参与度 + 变形梯度”分区。以某轿车正面100%重叠刚性壁障碰撞（示例工况，速度50 km/h）为例：前纵梁、前防撞梁、吸能盒、前舱纵梁后段属于主吸能区，单元尺寸取5～8 mm，保证压溃褶皱能由足够数量的单元捕捉；门槛梁、B柱、地板纵梁、前围板属于次吸能区，取8～12 mm；顶盖、后围、门内板非碰撞侧、内饰支架等远场区域可取15～25 mm。相邻区域之间应设置过渡带，单元尺寸过渡比控制在1.5～2以内，避免刚度突变引起局部应力波反射。

局部加密同样重要。焊点、粘胶连接、螺栓连接附近，以及圆角、翻边、接触可能发生的位置，应至少保证2～3层单元覆盖，否则接触搜索会漏检或穿透。最小单元尺寸一般不低于2～3 mm，过小单元会显著降低显式积分稳定步长，使计算时间成倍增加。

三、质量指标控制：翘曲、长宽比、雅可比与内角

壳单元质量指标中，翘曲（Warpage）描述四边形单元四个节点不共面的程度。翘曲过大时，单元面内与面外行为耦合，弯曲刚度失真。工程上目标值控制在5°以内，5°～15°为警告区，超过15°应优先调整节点位置或改用三角形单元过渡。长宽比（Aspect Ratio）是单元最长边与最短边之比，目标不超过5，警告区5～8。长宽比过大时，单元在不同方向刚度差异明显，容易在冲击下产生方向性误差。

雅可比（Jacobian）反映单元从参考坐标到实际坐标映射的畸变程度，取值0～1，越接近1越规则。壳单元目标值不低于0.7，0.5～0.7为警告区，低于0.5通常意味着单元严重扭曲，应重新划分。四边形内角建议控制在45°～135°，过尖或过钝的内角会降低积分精度。三角形单元在关键吸能区占比建议不超过5%，整体模型不超过10%，因为三角形单元刚度偏硬，且常伴随应力集中。

实体单元除上述指标外，还需关注体积比与最小高度。单点积分实体单元要检查沙漏能，若沙漏能超过总内能5%，应加密网格或调整沙漏控制系数。梁单元则需检查截面属性与单元长度匹配，避免单元长度小于截面尺寸导致局部刚度异常。

四、可执行的检查清单示例

完成划分后，建议先做一次无接触的短时（如5～10 ms）试算，观察能量曲线与沙漏能变化，确认网格没有引起初始穿透或负体积，再进入完整碰撞工况。网格质量不是一次性的验收项，而是贯穿建模、试算、正式计算与结果复核的持续控制过程。