2. 橡胶材料试验数据获取与处理
超弹性本构模型的精度，几乎完全取决于输入它的试验数据。Abaqus 本身不会“猜”材料行为，它只是把你在材料评估器中拟合出的应变能函数参数，代入到有限元方程里。因此，本章的核心不是“怎么点按钮”，而是“拿到什么样的数据、怎么修正、怎么判断它能不能用”。
一、三类基础试验的标准做法
对于各向同性超弹性材料，Abaqus 推荐至少提供单轴拉伸（Uniaxial）、等双轴拉伸（Biaxial）和平面拉伸（Planar）三种变形模式的数据。原因在于：单轴拉伸主要约束 I₁ 方向的行为，等双轴拉伸对 I₂ 更敏感，而平面拉伸（纯剪切）能同时提供两者的信息。只做单轴拉伸就拟合，往往在压缩或双轴工况下严重失真。
试验标准方面，橡胶类材料常参考 ISO 37、ASTM D412（单轴）、ASTM D575（压缩）、ISO 1827（剪切）等。试样通常为哑铃形或薄片，标距段长度 L₀ 一般取 20～50 mm，厚度 1～2 mm。加载速率建议控制在准静态范围（如 5～50 mm/min），因为橡胶有明显的黏弹性，速率过快会把黏弹性效应混入超弹性拟合。
关键记录量是：名义应力 = 力 / 初始截面积，名义应变 = 位移 / 初始标距。Abaqus 的超弹性模型默认使用名义应力—名义应变（工程应力—工程应变）数据，这一点必须与试验输出一致，否则拟合会整体偏移。
二、数据修正：从原始曲线到可用数据
原始试验曲线通常包含三类需要处理的问题：
- 初始刚度段（Toe region）：夹持松弛、试样对中不良会造成低应变段斜率异常。一般建议剔除
ε < 0.02的数据点，或按线性外推修正。 - 零应力偏移：力传感器零点漂移会让曲线在应变零点处不为零，需做基线平移。
- 循环软化（Mullins 效应）：首次加载曲线高于后续循环。若只做超弹性分析，应取首次加载曲线；若考虑 Mullins 效应，则需在 Abaqus 中启用相应模型。
修正后的数据应整理成两列：第一列名义应变，第二列名义应力，单位统一为 MPa。等双轴和平面拉伸的应力换算需注意：等双轴试验中，若用压力 p 加载，名义应力为 p·A/A₀；平面拉伸中，若试样宽度远大于厚度，可近似认为宽度方向应变近似为零。
三、体积压缩试验与不可压缩假设
橡胶通常被视为近似不可压缩，泊松比 ν ≈ 0.4995～0.4999。但在 Abaqus 中，若使用杂交单元（Hybrid element），可以显式输入体积模量 K；若使用普通单元，则需通过 D 参数控制可压缩性。
体积压缩试验（如 ASTM D575）用于测定体积模量。典型做法是将圆柱试样置于刚性约束中，施加静水压力，记录压力—体积应变曲线。体积模量 K = Δp / (ΔV/V₀)。对于大多数橡胶，K 在 1000～3000 MPa 量级。若试验数据缺失，可按 K = 2G(1+ν)/(3(1-2ν)) 估算，其中 G 为剪切模量。
四、数据导入 Abaqus 的完整路径
推荐流程是：Property 模块 → Material → Evaluate → Hyperelastic。在材料评估器中依次选择本构模型（如 Neo-Hookean、Mooney-Rivlin、Yeoh、Ogden），导入三类试验数据，观察拟合曲线与试验点的重合度，并查看稳定性检查（Drucker stability）是否通过。
导入时注意：
- 数据文件为纯文本，两列，制表符或空格分隔，不要带表头。
- 等双轴与平面拉伸数据需分别指定对应的变形模式，不能混用。
- 若某类试验缺失，可在评估器中勾选“允许外推”，但需谨慎，外推区域误差可能很大。
- 拟合完成后，将参数写入 Material 的 Hyperelastic 卡片，并检查单位制一致性。
五、拟合结果怎么判断
在材料评估器中，重点看三件事：
- 曲线重合度：拟合曲线应穿过三类试验点，尤其在
ε = 0.1～1.0的工作区间内不能有明显偏离。 - 稳定性检查：Abaqus 会给出 Drucker 稳定性提示。若出现不稳定，通常说明参数组合导致刚度矩阵非正定，需换模型或增加约束。
- 参数合理性：例如 Neo-Hookean 的
C₁₀应约为G/2，若拟合出负值或量级异常，说明数据或模型选择有问题。
最后，把评估器中确认的参数填入 Material 的 Hyperelastic 卡片，并在 Job 提交前用一个小变形算例做一次“单单元验证”：施加已知应变，检查反力是否与试验曲线一致。这一步能提前暴露单位、数据方向或模型选择错误，比整模型跑失败后再排查要高效得多。