4. 气囊织物材料模型与参数获取

气囊织物是典型的柔性各向异性材料，其力学行为与金属、塑料差异极大。在 LS-DYNA 中，织物通常用壳单元（*SECTION_SHELL）离散，材料模型的选择直接决定折叠展开、充气冲击和乘员接触力的可信度。本节围绕三类常用材料模型、织物非线性行为的来源，以及从试验数据到卡片参数的标定流程展开。

一、三类材料模型的定位与适用边界

MAT_FABRIC（*MAT_024 之外的 *MAT_034/*MAT_058 系列中的织物专用模型，常用 *MAT_FABRIC，即 *MAT_034 的织物版本，实际卡片号以所用版本为准）是气囊织物的首选模型。它支持正交各向异性弹性、经/纬向独立刚度、非线性应力-应变曲线、应变率效应以及失效删除。其核心特点是：织物几乎不能承压，只能承拉，因此模型中通过“拉为正、压为零刚度”的方式处理，避免出现不真实的压缩应力。

MAT_ELASTIC（*MAT_001）只提供线弹性各向同性或正交各向异性响应，没有非线性、没有失效、没有应变率。它适合做初步方案对比、织物刚度量级估算，或作为复杂模型的调试替身。若直接用它模拟气囊展开，会明显高估或低估峰值张力，因为真实织物在大应变下刚度会急剧上升。

MAT_PIECEWISE_LINEAR_PLASTICITY（*MAT_024）是金属塑性模型，带屈服面、硬化曲线和应变率。织物没有明显屈服平台，用它模拟织物属于“借用”，仅在需要快速验证接触算法、且不关心织物真实本构时临时使用。它的塑性流动假设与织物纱线滑移、涂层剥离的物理机制不符，因此不建议作为最终气囊模型。

二、织物非线性行为的三个来源

织物拉伸曲线通常呈“先软后硬”的 J 形，而不是金属的直线加屈服平台。其非线性主要来自三方面：纱线屈曲与伸直——初始阶段纱线从波浪状态被拉直，模量很低；纱线间摩擦与滑移——经纬交织点处纱线互相挤压，载荷通过摩擦传递，产生滞后；涂层与基布耦合——涂层在低应变时先承担部分载荷，高应变时基布主导。此外，织物对应变率敏感，高速充气下（应变率可达 10~100 /s）的模量明显高于准静态。

在 MAT_FABRIC 中，这些行为通过 LCID 定义的应力-应变曲线表达。曲线应使用真实应力-真实应变（或工程应力-工程应变，但需全文统一），并保证从原点出发、单调递增、无负斜率。若曲线出现下降段，LS-DYNA 会按卸载处理，导致能量异常。

三、参数标定的实际操作

标定不是“把试验曲线抄进卡片”，而是让单轴、双轴仿真同时复现试验。推荐流程如下：

1. 整理试验数据：将经向、纬向单轴拉伸的力-位移曲线换算为应力-应变。以 500 mm×50 mm 试样、厚度 0.35 mm 为例，若某点拉力 1200 N，则应力 = 1200/(50×0.35) ≈ 68.6 MPa；应变由引伸计或夹头位移扣除柔度后得到。
2. 建立单轴标定模型：取 100 mm×20 mm 壳单元条带，一端固定，另一端施加位移，用 *MAT_FABRIC 赋初值。初值可取低应变段割线模量，例如经向 200 MPa、纬向 180 MPa。
3. 提取仿真曲线并对比：输出拉伸端反力与位移，换算为应力-应变，与试验曲线叠图。若低应变段偏高，说明初始模量给大；若高应变段偏低，说明曲线后段刚度不足。
4. 修正 LCID 与 RATIO：按分段调整曲线控制点。例如在 ε=0.05 处试验应力 25 MPa、仿真 32 MPa，则将该点下调约 20%。应变率效应通过 RATIO（或 *DEFINE_CURVE 的缩放）引入，通常以准静态曲线为基准，高速曲线乘以 1.2~1.5。
5. 双轴验证：用 1:1 和 2:1 双轴试验检查经纬耦合。若双轴下经向应力明显高于试验，需检查 PRBA、PRCA 泊松比是否过大，或曲线是否只按单轴拟合。
6. 失效标定：用 FAIL 或 *MAT_ADD_EROSION 定义失效应变，通常取单轴断裂应变的 0.8~0.9 倍，避免单元过早删除。

四、参数示例与判断依据

下表给出一组用于教学演示的 MAT_FABRIC 参数示例（非实测值，仅用于说明字段含义与量级）。

标定完成的判据是：单轴经向、单轴纬向、双轴 1:1 三条仿真曲线与试验曲线的峰值误差均小于 10%，且初始段斜率偏差不超过 15%。若只有单轴数据，至少应保证经、纬向单轴复现良好，并在报告中注明双轴未验证。对于折叠气囊，还需检查展开过程中织物是否出现异常压缩应力或单元负体积，必要时降低时间步缩放系数或细化网格。

掌握上述模型差异、非线性来源和标定闭环后，后续气囊折叠、充气加载与乘员接触的仿真结果才具备物理可信度。参数获取的核心不是一次拟合，而是让模型在单轴、双轴、不同应变率下都能复现织物真实响应。