4. 橡胶超弹性单轴拉伸仿真

问题定位：为什么橡胶拉伸不能沿用金属弹塑性模型

橡胶单轴拉伸的典型特征是：初始刚度低、随伸长率增大而明显硬化、卸载后几乎完全恢复。用 Abaqus 模拟时，材料本构必须选择超弹性（Hyperelastic），而不是例 2 中的弹塑性双线性强化。本例采用两参数 Mooney-Rivlin 模型，其应变能密度函数为：

W = C10 * (I1 - 3) + C01 * (I2 - 3) + (1/D1) * (J - 1)^2

其中 I1、I2 为右 Cauchy-Green 变形张量的第一、第二不变量，J 为体积比。橡胶近似不可压缩，因此 D1 取极小值（或直接令 J≈1）。单轴拉伸下，名义应力与伸长比 λ 的关系可写为：

σ_nominal = 2*(λ - 1/λ^2) * (C10 + C01/λ)

该式说明：C10 主导小变形刚度，C01 主导大变形段的曲率。只要从试验曲线中取两个点即可反解出 C10 与 C01，这是本例参数标定的核心。

建模要点与操作顺序

橡胶试样按 GB/T 528 标准取哑铃形，仿真中简化为标距段长方体，利用对称性只建 1/8 模型以降低计算量。关键设置如下：

• 单元类型：C3D8H（八节点线性杂交单元）。杂交单元用于处理近似不可压缩材料的体积自锁，普通 C3D8 在 J≈1 时会产生虚假刚硬。
• 材料定义：在 Property 模块中创建 Hyperelastic → Mooney-Rivlin，输入 C10=0.28 MPa、C01=0.12 MPa、D1=0。D1=0 表示完全不可压缩，Abaqus 会自动采用混合公式处理。
• 分析步：Static, General，打开 NLGEOM=ON。初始增量步 0.01，最小 1e-6，最大 0.1。大变形下若不开 NLGEOM，刚度矩阵基于初始构形，结果严重偏硬。
• 边界条件：对称面施加法向约束；加载端施加位移 U2=标距×2.0（对应 λ=3.0），另一端固定。
• 输出请求：场输出选 S、LE、U；历史输出选 RF、U，用于绘制载荷-位移曲线。

参数标定计算示例

假设从试验曲线读取两个点：λ=1.5 时 σ=0.62 MPa；λ=3.0 时 σ=2.85 MPa。代入单轴拉伸公式：

点1: 0.62 = 2*(1.5 - 1/1.5²)*(C10 + C01/1.5)
     = 2*(1.5 - 0.4444)*(C10 + 0.6667*C01)
     = 2.1112*(C10 + 0.6667*C01)

点2: 2.85 = 2*(3.0 - 1/3.0²)*(C10 + C01/3.0)
     = 2*(3.0 - 0.1111)*(C10 + 0.3333*C01)
     = 5.7778*(C10 + 0.3333*C01)

联立解得：C10 = 0.28 MPa, C01 = 0.12 MPa

将上述参数输入 Abaqus 后，提取加载端 RF 并除以初始截面积得到名义应力，与试验曲线对比。若 λ=1.5 和 λ=3.0 两点的偏差均小于 5%，则认为参数标定有效。

结果验证清单