3. 安全气囊建模的几何与网格策略
安全气囊仿真与整车碰撞仿真的最大差别在于：气囊是一个由织物膜构成的、内部充有高压气体的封闭腔体，其几何在展开过程中会发生剧烈的大变形。因此，气囊建模的核心矛盾不是"能不能算"，而是"用什么几何、什么单元、什么网格密度，才能既算得稳、又算得准、还算得快"。本节围绕折叠前处理、网格类型选择、壳单元与膜单元适用性、网格质量检查四个环节展开。
一、折叠前处理：从展开态到折叠态
气囊的 CAD 模型通常是展开态（flat pattern）或缝合后的三维袋体，而仿真需要的是折叠后收拢在模块内的初始构型。折叠前处理要解决三件事：几何清理、折叠建模、初始穿透消除。
几何清理阶段，需要缝合 CAD 中的自由边、删除微小圆角与工艺孔、统一法向。气囊织物厚度通常在 0.2～0.4 mm 量级，若按真实厚度建壳单元，单元长宽比会严重失衡，因此工程上普遍采用"零厚度"或"等效厚度"壳单元，厚度只用于接触和材料计算，不参与几何偏置。
折叠建模有两条主流路线。一是几何折叠法：在 CAD 或专用前处理器（如 LS-PrePost 的 Folding 模块、FoldPRO）中按实际折叠工艺（如 Z 折、卷折、星形折叠）逐步生成折叠构型，优点是几何真实、可复现工艺，缺点是耗时长、复杂折叠易失败。二是映射折叠法：先建立展开网格，再通过坐标映射把网格"压"进模块包络盒内，优点是快，缺点是初始构型与真实折叠差异大，展开过程可能失真。
折叠完成后必须做初始穿透检查。折叠层之间、气囊与模块壳体之间若存在穿透，LS-DYNA 会在计算初期产生巨大的接触力，导致气囊"炸开"或直接报负体积。常用做法是：用 *CONTACT_AUTOMATIC_SINGLE_SURFACE 配合 SOFT=1，并在关键字中设置 IGNORE=2 忽略初始穿透，同时通过 *CONTROL_TIMESTEP 的 DTINIT 控制初始时间步。
二、网格类型选择：壳、膜与实体
气囊织物的力学行为是"只能受拉、几乎不能受压、抗弯刚度极低"。这一物理特性直接决定了单元类型的选择。
壳单元（Shell）：LS-DYNA 中的 *SECTION_SHELL 支持多种壳公式。对于气囊，常用 ELFORM=2（Belytschko-Tsay）或 ELFORM=16（全积分）。壳单元能承受面内拉压和面外弯曲，若直接用于气囊，会人为引入抗弯刚度，使气囊"发硬"、展开形态偏圆。因此壳单元更适合气囊的结构件，如气体发生器壳体、模块盖板、拉带加强片。
膜单元（Membrane）：LS-DYNA 中通过 *SECTION_SHELL 设置 ELFORM=9（膜单元）或使用 *MAT_FABRIC 配合膜公式实现。膜单元没有弯曲刚度，只传递面内拉力，与织物真实力学行为高度一致，是气囊本体的首选。其代价是：膜单元对沙漏和负体积更敏感，需要配合沙漏控制（IHQ=4 或 8）和合理的网格密度。
实体单元（Solid）：仅用于气体发生器内部结构、点火具等体构件，不用于气囊织物。
一个实用的判断原则是：气囊织物用膜单元，气囊附件用壳单元，气体发生器用实体单元。若模型中同时存在膜与壳，需注意两者连接处的节点耦合方式，避免出现"膜单元被壳单元拉裂"的假象。
三、网格密度与质量检查
气囊网格尺寸没有统一标准，但有一个经验起点：折叠区域网格加密，展开区域适当放粗。折叠区曲率变化剧烈，若网格过大，折叠层之间会互相穿透；展开区若网格过密，计算量会成倍增加。工程上常见的做法是：折叠区单元边长 3～5 mm，展开区 8～12 mm，整体单元数控制在 2 万～8 万之间（以驾驶员气囊为例）。
网格质量检查应至少覆盖以下指标：
- 长宽比（Aspect Ratio）：建议小于 5。气囊膜单元在折叠后容易被压扁，长宽比过大会导致沙漏和负体积。
- 翘曲角（Warpage）：建议小于 15°。翘曲过大的四边形单元在展开时会产生虚假的褶皱。
- 雅可比（Jacobian）：建议大于 0.6。低于 0.3 的单元在 LS-DYNA 中极易出现负体积。
- 最小边长：不宜小于 1 mm，否则会显著压低显式时间步，拖慢计算。
- 沙漏能占比：计算完成后检查
glstat文件，沙漏能应小于内能的 5%。若超标，需加密网格或调整IHQ。
*CONTROL_TERMINATION 设置一个很短的终止时间（如 5 ms），快速跑一遍展开初期，观察是否有单元畸变。
四、一个可复用的检查清单
把上述内容整理成一张可执行的检查表，便于在每次建模后逐项核对：
| 检查项 | 目标值 / 操作 | 不满足时的处理 |
|---|---|---|
| 折叠层初始穿透 | 无穿透或 IGNORE=2 | 调整折叠间隙，或改用 SOFT=1 接触 |
| 单元类型 | 袋体 ELFORM=9 | 检查 *SECTION_SHELL 设置 |
| 长宽比 | < 5 | 局部重划，加密折叠区 |
| 雅可比 | > 0.6 | 删除并重划低质量单元 |
| 沙漏能 / 内能 | < 5% | 提高 IHQ，或加密网格 |
| 初始时间步 | 由最小单元控制 | 用 *CONTROL_TIMESTEP 设 DTINIT |
需要强调的是，气囊建模的几何与网格策略不是一次性的，而是与后续的充气模型、接触定义、材料参数相互耦合的。例如，若采用 *AIRBAG_WANG_NEFSKE 均匀压力法，网格密度对压力分布影响较小；但若采用 *AIRBAG_PARTICLE 或 CPM 粒子法，网格过粗会导致粒子穿透袋体。因此，网格策略应在明确充气建模方法之后再最终确定，而不是孤立地追求"网格好看"。