2. 橡胶材料本构理论精讲
在 Abaqus 中做橡胶结构仿真，第一步不是画网格，而是决定用什么本构模型描述材料。橡胶的力学行为与金属截然不同：它可以承受几百个百分点的可恢复大变形，刚度随变形增大而上升，卸载后几乎回到原形，体积几乎不可压缩。这些特征决定了它不能用弹性模量加泊松常数简单描述，而必须借助超弹性本构模型。本节把这条理论主线讲清楚：应变能密度函数是什么、各向同性和不可压缩性如何约束模型、常见模型如何分类，以及工程上如何选型。
一、应变能密度函数：超弹性的核心
超弹性材料假设存在一个标量函数 W，称为应变能密度函数（单位参考体积的储能）。应力不是独立变量，而是 W 对变形的导数。用右 Cauchy-Green 张量 C 表示时，第二 Piola-Kirchhoff 应力为：
S = 2 * ∂W/∂C这个关系意味着：只要写出 W 的形式，材料的全部应力-应变响应就确定了。因此“选本构”本质上是“选 W 的函数形式”。W 通常写成应变不变量的函数。三个基本不变量为：
- I₁ = tr(C) = λ₁² + λ₂² + λ₃²，反映整体拉伸程度；
- I₂ = ½[(tr C)² − tr(C²)]，反映面积变化相关的变形；
- I₃ = det(C) = (λ₁λ₂λ₃)²，直接关联体积变化。
其中 λ₁、λ₂、λ₃ 是三个主伸长比。橡胶的 W 一般只依赖 I₁ 和 I₂，而 I₃ 通过不可压缩约束被“锁住”。
二、各向同性假设与不可压缩性约束
各向同性假设指材料无方向性，未变形状态下任意方向力学响应相同。对未填充或填充均匀的橡胶，这一假设成立，于是 W 只能通过不变量 I₁、I₂、I₃ 依赖 C，而不能直接依赖某个方向。若橡胶中掺入大量定向纤维或存在明显取向，则需改用各向异性超弹性模型，这已超出本节范围。
不可压缩性约束指变形过程中体积不变，即 J = det(F) = λ₁λ₂λ₃ = 1，等价于 I₃ = 1。橡胶的体模量比剪切模量高约三个数量级，近似不可压缩是合理简化。在 Abaqus 中处理方式有两种：
- 完全不可压缩：用杂交单元（如 C3D8H、C3D10H）引入静水压力作为独立自由度，避免体积自锁；
- 近似不可压缩：给一个很大的但有限的体积模量 D（对应 K），用普通单元即可，但 D 过大会导致收敛困难。
三、模型分类与选择依据
按 W 的多项式展开阶数，常见模型可分为：
- Neo-Hookean：W = C₁₀(I₁ − 3)，只有一个参数，适合小到中等应变、单轴拉伸的粗略估算；
- Mooney-Rivlin：W = C₁₀(I₁ − 3) + C₀₁(I₂ − 3)，两参数，能较好拟合中等应变；
- Yeoh：三阶 I₁ 多项式，不含 I₂，对单轴、等双轴、平面拉伸的综合拟合能力强，是工程首选之一；
- Ogden：以主伸长比幂函数表达，参数多、拟合精度高，适合大应变和复杂载荷；
- Arruda-Boyce：基于分子链统计，物理意义清晰，适合大变形和循环载荷；
- Marlow：不显式给 W，直接从试验数据构造，适合只有单一种类试验数据的情况。
选择依据可归纳为三条：应变范围、可用试验数据、载荷类型。小应变选 Neo-Hookean 或 Mooney-Rivlin；中等应变且只有单轴数据选 Yeoh；大应变、多轴复杂载荷选 Ogden 或 Arruda-Boyce；若只有单轴拉伸数据，Marlow 可避免外推失真。
四、参数拟合与稳定性检查
选定模型后，需要用试验数据拟合参数。Abaqus 的 Material Evaluation 支持输入单轴、等双轴、平面拉伸、体积压缩等数据，自动完成最小二乘拟合。以 Yeoh 模型为例，其形式为：
W = C10*(I1-3) + C20*(I1-3)^2 + C30*(I1-3)^3拟合时要注意：不要只用一种试验数据外推到其他载荷类型。例如仅用单轴拉伸拟合的 Mooney-Rivlin，在等双轴压缩下可能给出错误甚至负刚度。判断拟合质量，除了看应力-应变曲线重合度，还要检查 Abaqus 给出的稳定性警告：若 W 的 Hessian 矩阵非正定，模型在部分变形路径上会失稳，求解器可能报错或给出非物理结果。
五、工程选型速查
| 场景 | 推荐模型 | 关键理由 |
|---|---|---|
| 密封圈小变形刚度估算 | Neo-Hookean | 参数少，快速评估 |
| 中等应变减振垫 | Mooney-Rivlin / Yeoh | 兼顾精度与数据需求 |
| 大变形多轴载荷 | Ogden / Arruda-Boyce | 拟合能力强，物理意义清晰 |
| 仅有单轴拉伸数据 | Marlow | 避免 I₂ 外推失真 |
| 循环载荷、Mullins 效应 | Yeoh + 损伤/粘弹性 | 可叠加永久软化模型 |
总结本节逻辑：W 决定应力，各向同性把 W 限制为不变量的函数，不可压缩性把 I₃ 锁为 1，模型分类本质是 W 的展开形式不同，选型则取决于应变范围、数据覆盖度和载荷类型。把这四点串起来，后续在 Abaqus 中定义材料、拟合参数、选择单元时就有了明确依据。