4. 橡胶材料试验数据获取与处理
超弹性本构模型的精度，几乎完全取决于输入试验数据的质量与覆盖范围。Abaqus 中常用的 Neo-Hookean、Mooney-Rivlin、Yeoh、Ogden 等模型，其参数都是通过对试验应力—应变曲线拟合得到的。如果只做单轴拉伸就拟合，模型在压缩、剪切或大变形工况下往往严重失真。因此，工程上通常要求至少提供单轴拉伸（UT）、等双轴拉伸（ET）和平面拉伸（PT）三类基础试验数据，才能较完整地约束模型参数。
三类基础试验的力学状态与标准
三类试验的本质区别在于变形模式，也就是主伸长比之间的关系。设三个主方向伸长比为 λ₁、λ₂、λ₃，并约定加载方向为 1 方向：
- 单轴拉伸（UT）：λ₁ = λ，λ₂ = λ₃ = λ-1/2，一个方向受拉，另两个方向自由收缩，常用于获取拉伸段响应。
- 等双轴拉伸（ET）：λ₁ = λ₂ = λ，λ₃ = λ-2，两个方向同时等量拉伸，第三方向自由收缩，等效于压缩段响应。
- 平面拉伸（PT）：λ₁ = λ，λ₂ = 1，λ₃ = λ-1，一个方向拉伸、一个方向约束不变形，接近纯剪切状态。
试验标准方面，单轴与平面拉伸常参照 ISO 37、ASTM D412；等双轴拉伸没有统一的国际标准，多采用专用的双轴夹具或气泡膨胀法，参照 ISO 23529 的试样制备与调节要求。试样通常为哑铃形或矩形条，标距段厚度建议 1～2 mm，标距长度 20～50 mm。加载速率一般控制在 10～50 mm/min，对应名义应变率约 10⁻³～10⁻² s⁻¹，以尽量接近准静态。
需要特别强调的是预循环。橡胶首几次加载存在明显的 Mullins 效应，应力软化显著。标准做法是先做 3～5 次加卸载循环，取第 3 次或稳定后的循环作为有效数据，这样得到的曲线才具有可重复性。
数据修正：从名义量到真实量
试验机输出的是名义应力 σ_nom（力除以初始截面积）和名义应变 ε_nom（位移除以初始标距）。超弹性模型内部使用的是真实应力与真实应变，因此必须换算：
σ_true = σ_nom · (1 + ε_nom)
ε_true = ln(1 + ε_nom)以一组示例数据说明：某橡胶试样初始截面积 A₀ = 20 mm²，标距 L₀ = 25 mm。当拉力 F = 60 N、位移 ΔL = 12.5 mm 时，σ_nom = 60/20 = 3.0 MPa，ε_nom = 12.5/25 = 0.5。换算后 σ_true = 3.0 × 1.5 = 4.5 MPa，ε_true = ln(1.5) ≈ 0.405。可见大变形下两者差异明显，若直接用名义量拟合，模型会偏软。
除换算外，还需做三项修正：一是剔除预循环段，只保留稳定循环的加载支；二是去除夹持端滑移与初始松弛段，通常舍弃应变小于 0.02 的数据点；三是统一应变率，三类试验尽量在同一速率下完成，否则拟合出的参数会隐含黏弹性影响。
曲线拟合与模型校准流程
在 Abaqus/CAE 中，通过 Property 模块的 Material → Hyperelastic 进入拟合界面。以 Yeoh 模型为例，其应变能函数为 W = C₁₀(I₁−3) + C₂₀(I₁−3)² + C₃₀(I₁−3)³，只需单轴数据即可初步拟合，但强烈建议同时勾选等双轴与平面数据。操作要点如下：
- 在 Test Data 中分别录入 Uniaxial、Biaxial、Planar 三组真实应力—应变数据，注意单位统一为 MPa。
- 选择模型阶数。Mooney-Rivlin 适合应变 < 100% 的中小变形；Yeoh 适合大变形且对压缩友好；Ogden 阶数高、拟合能力强但参数多，需警惕不稳定。
- 点击 Evaluate，Abaqus 用最小二乘法求解参数，并给出各试验曲线的拟合对比图与残差。
- 检查稳定性：Drucker 稳定性要求参数满足特定不等式，若出现负刚度或曲线交叉，应降低阶数或补充数据。
常见问题与处理
实际工作中最常遇到三类问题。第一，只有单轴数据：此时可用 Yeoh 一阶或 Neo-Hookean 做粗略估算，但必须明确告知模型在压缩与剪切工况下不可信。第二，曲线出现负斜率：多因试样打滑或夹持端应力集中，应检查夹具并重新试验。第三，拟合参数为负：通常意味着数据范围过窄或噪声过大，可尝试截取有效段或改用低阶模型。
完成拟合后，建议把参数回代到一个小型验证模型中，例如单轴压缩或简单剪切工况，对比仿真与试验的力—位移曲线。只有三类基础试验数据齐全、修正规范、拟合残差可控，超弹性模型才具备工程可信度。