4. 中面抽取与壳单元简化
薄壁结构是汽车、家电、航空钣金件中最常见的形态。一块2 mm厚的钢板，如果按实体单元离散，厚度方向至少要铺2～3层单元才能让弯曲刚度不失真；而按壳单元处理，厚度只是一个截面属性参数，网格只铺在板面上。两者自由度差距通常在10～50倍量级。本节要解决的就是：如何把已经建好的实体薄壁模型，可靠地转成壳单元模型，并保证转换后刚度、连接和传力路径不发生明显失真。
一、先判断“能不能抽”：厚度与特征尺度比
中面抽取的前提是结构满足薄壁假设。工程上常用的判据是：板厚 t 与面内最小特征尺寸 L 的比值 t/L 小于约 1/10，且厚度方向应力梯度可以忽略。例如某支架板厚 1.5 mm，面内最小圆角半径 8 mm，t/L≈0.19，属于偏厚区域，抽取后需要关注圆角处的刚度损失；而一块 1.2 mm 厚、面内尺寸 300 mm 的门板，t/L≈0.004，抽取后壳单元结果与实体结果通常能控制在 5% 以内。
实际操作中，不必逐块手算。可以在前处理中按厚度分组显示，把 t/L 明显偏大的区域（如厚法兰、加强凸台、实心凸缘）标记为“保留实体”，其余薄壁区域进入抽取流程。这一步决定了后续是“全壳模型”还是“壳—实体混合模型”。
二、中面抽取的三种典型方式
主流前处理工具（HyperMesh、ANSA、LS-PrePost 等）提供的中面抽取方式大致分三类：
- 偏置法（Offset）：选取实体外表面，沿法向偏置 t/2 得到中面。适合厚度均匀、表面平整的板件，速度快，但对圆角和曲面过渡区容易产生自交。
- 配对法（Midsurface / Pair）：识别实体上下两个对应面，取其中间位置生成中面。对变厚度、带拔模的钣金件更稳健，但需要面配对正确。
- 直接建模法：对已经丢失实体或本身只有曲面的模型，直接在曲面基础上补出中面并赋予厚度。适合逆向或导入模型。
抽取完成后，必须检查三件事：中面是否连续、厚度方向是否一致、自由边是否与原始实体边界吻合。常见问题是圆角被抽成尖角、翻边高度丢失、局部面片重叠。这些缺陷在壳网格划分后表现为单元畸变或初始穿透，必须在划网格前修掉。
三、厚度偏置与截面定义
壳单元本身没有厚度几何，厚度通过截面属性（*SECTION_SHELL）定义。LS-DYNA 中常用的是 ELFORM=2（Belytschko-Tsay）或 ELFORM=16（全积分），厚度方向积分点 NIP 一般取 3～5。对于需要输出厚度方向应力的场合，NIP 取 5 更合适。
偏置方向必须与中面法向一致。如果抽取时中面法向混乱，会导致壳单元正反面颠倒，接触和压力加载方向出错。处理办法是统一法向：选取一个参考面，用 Normal 工具把所有壳单元法向调整到一致，再检查自由边法向是否连续。
四、连接处理：从实体焊点到壳单元连接
实体模型中的焊接、螺栓、铆接，在壳模型中不能直接沿用。常见处理方式有：
- 共节点：两块壳在翻边处直接共享节点，适合连续焊缝。要求网格在连接线上对齐。
- 刚性连接（*CONSTRAINED_NODAL_RIGID_BODY）：用一组节点绑定到主节点，适合点焊、螺栓连接。注意不要跨越大变形区域，否则会引入虚假刚度。
- 梁单元模拟焊点：用短梁连接两层壳，可传递力和弯矩，适合需要观察焊点受力的场合。
- 接触绑定：用 *CONTACT_TIED_SHELL_EDGE_TO_SURFACE 等关键字，适合网格不匹配的装配连接。
连接处理的核心原则是：连接刚度不能高于母材，也不能低到让结构“散架”。点焊直径通常取 4～6 mm，对应刚性连接的影响半径约 3～5 mm；如果焊点间距小于 30 mm，建议改用连续绑定或梁单元，避免局部刚度过高。
五、简化后的检查与验证
壳模型建好后，不要直接投入大批量计算。先做三项检查：
- 自由模态对比：对同一结构分别做实体和壳的自由模态分析，比较前 10 阶频率。若偏差超过 10%，说明中面或连接处理有问题。
- 静刚度对比：在典型加载点施加单位力，比较位移。壳模型通常偏刚，若偏刚超过 15%，需检查圆角和连接。
- 质量核对：壳模型总质量应与实体模型一致，偏差应小于 2%。质量偏差大通常意味着厚度或密度定义错误。
下面这张图给出从实体薄壁件到壳单元模型的关键处理链，以及每一步需要关注的参数和检查点。
六、一个可复用的检查清单
把上面的流程固化成检查表，可以在每次壳简化后快速过一遍：
| 检查项 | 判断依据 | 不通过时的处理 |
|---|---|---|
| 厚度比 t/L | 薄壁区 t/L < 1/10 | 厚区保留实体，做混合模型 |
| 中面连续性 | 无自由边错位、无面片重叠 | 修补面、重新配对 |
| 法向一致性 | 所有壳单元法向同侧 | 统一法向后再划网格 |
| 截面属性 | ELFORM、NIP、厚度与设计一致 | 逐组核对 *SECTION_SHELL |
| 连接刚度 | 焊点影响半径 3～5 mm | 改用梁焊点或绑定接触 |
| 模态/刚度/质量 | 频率偏差 < 10%，位移偏差 < 15%，质量偏差 < 2% | 回修中面或连接 |
需要强调的是，壳简化不是“越简越好”。对于局部需要观察厚度方向应力、接触压溃或三维裂纹扩展的区域，应保留实体单元，形成壳—实体混合模型。混合模型的关键是壳与实体之间的过渡：常用 *CONSTRAINED_SHELL_TO_SOLID 或共用节点方式，过渡区网格密度要匹配，避免出现刚度突变。
完成上述步骤后，一个典型的薄壁支架模型可以从约 80 万实体单元降到约 3 万壳单元，自由度从约 240 万降到约 18 万，在保持前 10 阶模态频率偏差小于 8% 的前提下，显式分析时间步长也可从 0.3 μs 提升到 1.2 μs 量级。这就是中面抽取与壳单元简化在 LS-DYNA 计算加速中最直接的价值。