4. 薄壁结构有限元建模流程
薄壁吸能结构的仿真精度，往往在网格划分之前就已经被决定了。一个薄壁帽形梁、方管或双帽梁，如果几何模型带着倒角碎面、重复面、微小缝隙进入网格阶段，后续无论用多细的网格、多好的材料卡，结果都会出现刚度虚高、接触穿透或能量异常。因此，本节把建模流程拆成五个相互制约的环节：几何清理、中面抽取、网格划分策略、单元类型选择、厚度定义，并说明每个环节的输入、判断依据和常见陷阱。
一、几何清理：先让模型“干净”
薄壁件通常由钣金冲压或辊压成形，CAD 模型里常见的问题包括：圆角半径过小（R0.5 以下）、面与面之间存在微小缝隙或重叠、翻边处有非流形边、对称件只做了一半、以及大量对吸能贡献很小的安装孔和凸台。几何清理的目标不是把模型修得漂亮，而是让中面抽取和网格划分能够稳定进行。
实际操作中，建议按以下顺序处理：
- 缝合与容差检查：先以 0.01～0.05 mm 的容差缝合自由边，观察剩余自由边数量。若自由边集中在翻边末端，通常是正常的开口边；若出现在平面中部，说明存在裂缝或重复面。
- 去除微小特征：对吸能贡献可忽略的圆角（R<1 mm）、小孔（直径<3 mm）、凸台和印痕，直接删除或填充。保留这些特征会迫使网格尺寸降到 0.5 mm 以下，显著增加计算量。
- 简化翻边与搭接：点焊搭接区域可简化为共面或小台阶，避免抽取中面后出现两层中面间距过小、网格重叠。
- 对称处理：若结构几何与载荷均对称，可只保留 1/2 或 1/4 模型，在对称面上施加对称边界条件，减少网格量并提高稳定性。
二、中面抽取：从实体到壳
薄壁结构厚度方向尺寸远小于面内尺寸，用壳单元比实体单元更经济。中面抽取就是把实体壁厚的中位面提取出来，后续所有壳单元都铺在这个面上。抽取质量直接决定厚度定义是否准确。
抽取时要注意三点。第一，厚度方向必须一致：每个中面的法向应指向壳的同一侧（通常指向外侧或内侧），否则后续偏置和接触方向会混乱。第二，变厚度区域要单独处理：如加强板、拼焊板（Tailor Welded Blank）厚度从 1.2 mm 变到 1.8 mm，应在中面上划分出对应区域，分别赋厚度，而不是取平均值。第三，抽取后要检查面连续性：相邻中面之间应共享边或留有可控间隙，避免出现悬空面。
对于无法稳定抽取中面的复杂区域（如铸铝节点、厚翻边），可保留局部实体，用实体单元与壳单元过渡连接，但过渡区要避免刚度突变。
三、网格划分策略：尺寸、形状与过渡
薄壁吸能仿真的网格策略核心是“在关键区域足够细，在非关键区域足够粗”。推荐做法如下：
- 基础尺寸：以最小特征尺寸为基准。若结构最小圆角为 R3 mm，目标单元尺寸可取 1.5～2 mm；若最小圆角为 R1 mm，则取 0.5～1 mm。
- 变形集中区加密：压溃诱导槽、折角、焊点周围、翻边根部是塑性铰形成位置，单元尺寸应比基础尺寸小 30%～50%。
- 单元形状：壳单元优先用四边形，长宽比控制在 3:1 以内，翘曲角小于 15°，避免出现三角形密集区。三角形单元在弯曲中偏刚，只用于过渡。
- 厚度方向：壳单元不需要厚度方向网格，但若局部改用实体，厚度方向至少 2 层单元。
- 网格质量检查：划分后检查最小内角（>30°）、雅可比（>0.6）、重复节点和自由边，确保没有孤立单元。
四、单元类型选择：壳、实体与连接
LS-DYNA 中薄壁结构常用单元类型如下：
| 单元类型 | 适用场景 | 关键设置 |
|---|---|---|
| Belytschko-Tsay 壳（ELFORM=2） | 常规薄壁冲压件、辊压件 | 默认壳单元，效率高，需设置厚度与积分点 |
| Belytschko-Levathan 壳（ELFORM=1） | 翘曲明显、大变形区域 | 比 ELFORM=2 稍慢，但翘曲精度更好 |
| 全积分壳（ELFORM=16） | 需要避免沙漏的局部区域 | 计算量较大，慎用 |
| 实体单元（ELFORM=1 或 -1） | 厚壁节点、铸铝接头、填充泡沫 | 厚度方向至少 2 层，配合壳-实体连接 |
| 焊点（*CONSTRAINED_SPOTWELD 或梁单元） | 点焊搭接 | 焊点直径、失效准则需与试验一致 |
壳单元沿厚度方向通常取 5 个积分点（NIP=5），足以捕捉弯曲和膜力。若关注厚度方向应力分布，可增加到 7 个，但收益有限。沙漏控制建议用类型 4 或 5，系数取 0.03～0.10，并通过能量比（沙漏能/总内能）判断是否小于 5%。
五、厚度定义：从几何到属性
壳单元的厚度不是几何厚度，而是属性厚度。定义方式有两种：
- 直接赋厚度：在 *SECTION_SHELL 中设置 T1～T4 厚度值，适用于等厚度区域。
- 偏置中面：若中面抽取时已按实际中位面建立，厚度直接取板厚；若中面位置有偏差，可用偏置量修正，但偏置会影响接触和惯性，需谨慎。
对于拼焊板或变厚度区域，应建立多个 *SECTION_SHELL，分别赋给对应单元集。厚度值必须与材料卡中的密度、弹性模量一起检查，确保单位制一致。例如采用 mm-ton-s 单位制时，钢的密度为 7.85e-9 ton/mm³，弹性模量为 210000 MPa，厚度单位为 mm。
六、建模流程的整体关系
五个环节不是线性独立的：几何清理决定中面能否抽取，中面质量决定网格能否铺好，网格尺寸决定单元类型和厚度赋值的粒度，厚度定义又反过来影响接触和沙漏。下图给出这一流程的输入、判断与输出关系，便于在建模时逐项核对。
把这张关系图当作建模检查表使用：每完成一个环节，回到对应判断条件确认一次。几何清理不彻底，中面抽取就会失败；中面法向不统一，厚度偏置和接触方向就会出错；网格尺寸与最小圆角不匹配，单元质量就会拖垮显式积分步长。厚度定义完成后，建议先跑一个 5～10 ms 的短时压溃，观察沙漏能、接触能和总内能的比例，再决定是否进入完整工况。这样可以在投入大量计算资源之前，尽早发现建模层面的问题。