4. 单元类型选择与网格划分策略

在整车碰撞仿真中，网格是承载所有物理计算的骨架。单元类型选错，再精细的材料曲线和接触设置也无法挽救结果；网格尺寸与质量失控，则会出现负体积、沙漏能异常、时间步骤降甚至计算中断。本节围绕车身结构中三类主力单元——壳单元、实体单元、梁单元——的适用边界，以及网格尺寸与质量准则的工程判断方法展开。

一、三类单元的力学假设与适用场景

壳单元基于板壳理论，通过厚度方向的积分点描述弯曲与膜行为，适合一个方向尺寸远小于另外两个方向的结构。车身中的覆盖件、A柱/B柱加强板、门槛内板、地板、前后防撞梁薄壁结构，厚度通常在0.6～2.5 mm，而面内尺寸在几十到几百毫米量级，是壳单元的典型对象。在LS-DYNA中常用*SECTION_SHELL配合ELFORM=2（Belytschko-Tsay）或ELFORM=16（全积分）定义，前者效率高、适合大规模车身，后者在翘曲敏感区域更稳定。

实体单元适合三个方向尺寸接近、需要捕捉厚度方向应力梯度或三维应力状态的对象。车身中典型应用包括：发动机悬置支架、铸铝减震塔、副车架接头、电池包箱体加强块、螺栓连接局部、以及需要模拟压溃褶皱细节的吸能盒。常用*SECTION_SOLID配合ELFORM=1（单点积分）或ELFORM=-1（全积分），单点积分需配合沙漏控制，全积分计算成本高但不易出现沙漏。

梁单元适合一个方向尺寸远大于另外两个方向、且主要承受轴向力、剪力与弯矩的构件。车身中典型对象是车门防撞杆、座椅骨架管件、稳定杆、部分副车架管梁。LS-DYNA中常用*SECTION_BEAM配合ELFORM=1（Hughes-Liu）或ELFORM=2（Belytschko-Schwer），需定义截面形状与积分规则。梁单元计算效率极高，但无法描述截面局部变形，因此只适用于截面刚度足够、失效模式以整体弯曲或轴向压溃为主的结构。

二、网格尺寸的确定逻辑

网格尺寸不是越细越好，而是要在计算精度与时间步之间取得平衡。显式积分的时间步由最小单元尺寸控制，公式为 Δt = L_min / c，其中 c 为材料声速。以钢为例，c ≈ 5000 m/s，若最小单元尺寸为 5 mm，则 Δt ≈ 1×10⁻⁶ s；若缩小到 2 mm，Δt 降至 4×10⁻⁷ s，计算时间成倍增加。因此工程上通常采用分区策略：远离碰撞区的基础网格取 8～10 mm，A柱、B柱、门槛、前纵梁等主要吸能与变形区域取 5 mm，接触面、圆角、褶皱诱导槽等局部取 3 mm。对于薄壁吸能盒，若厚度方向需要多个实体单元，则面内尺寸可放宽到 4～6 mm，厚度方向至少 2 层。

网格过渡要平滑，相邻单元尺寸比不宜超过 2:1，避免刚度突变引起应力波反射。在LS-DYNA中可通过*CONTROL_TIMESTEP设置TSSFAC（时间步缩放因子，默认0.9）和DTMIN（最小时间步），当单元尺寸过小时触发质量缩放，但缩放比例应控制在5%以内，否则会显著改变惯性效应。

三、网格质量准则与检查方法

壳单元重点检查四项指标：翘曲角（Warpage）建议小于15°，长宽比（Aspect Ratio）小于5，雅可比（Jacobian）大于0.6，最小内角大于30°、最大内角小于150°。实体单元重点检查雅可比大于0.3、体积比（Volume Ratio）小于5、最小内角大于20°。梁单元主要检查截面定义与节点偏置是否正确，单元长度不宜小于截面高度的2倍。

在HyperMesh或ANSA中完成网格划分后，应导出LS-DYNA关键字文件，利用LS-PrePost的“Check”功能或*CONTROL_ENERGY中的沙漏能监控进行验证。若初始沙漏能超过总内能的5%，说明单点积分单元存在零能模式，需调整沙漏控制系数或改用全积分。若出现负体积，通常源于实体单元在压缩下过度畸变，应检查材料失效定义与单元删除准则。

四、单元类型与网格策略的协同

实际车身模型中，壳、实体、梁往往混合使用。例如前纵梁主体用壳单元，后端与副车架连接处用实体单元过渡，中间通过*CONSTRAINED_INTERPOLATION或*CONSTRAINED_NODAL_RIGID_BODY耦合。梁单元用于车门防撞杆时，两端需与门框壳单元建立刚性连接，避免出现“悬空”梁。混合建模的关键是保证连接区域刚度传递连续，耦合节点不应产生局部应力集中。

网格划分完成后，建议用一段10 ms的刚性墙低速冲击进行“网格体检”，观察沙漏能、接触力曲线和变形模式是否合理。若沙漏能突增或接触力振荡剧烈，优先排查网格质量与单元类型匹配问题，再调整接触参数。这一步骤能在正式碰撞工况前暴露大部分网格隐患，显著降低返工成本。