1. 橡胶材料与超弹性仿真入门

橡胶密封圈、轮胎胎面、减震支座、医用硅胶管，这些零件有一个共同点：受力后能发生很大的可恢复变形。金属拉伸0.2%就进入屈服，而天然橡胶可以轻松拉伸到原长的5～7倍，卸载后仍能基本恢复原状。这种"大变形+可恢复"的组合，正是超弹性（hyperelasticity）要描述的对象。本章先建立对橡胶材料本质的认识，再说明Abaqus中一条完整的超弹性仿真链路长什么样。

一、从分子结构理解橡胶为什么"软而韧"

橡胶由长链高分子构成，链与链之间通过硫化形成的化学交联点连接成三维网络。每条分子链在自然状态下呈卷曲、纠缠的构象，熵值很高。当外力拉伸时，卷曲链被拉直，构象数减少、熵降低；卸载后分子热运动驱动链重新回到高熵的卷曲状态，宏观上表现为回弹。因此橡胶的弹性主要不是"内能弹性"（像金属靠原子键拉伸储能），而是"熵弹性"。

这个机理带来两个工程后果。第一，橡胶的剪切模量很低（典型值0.3～1.5 MPa量级），体积模量却接近水（约1000～2000 MPa），即几乎不可压缩，泊松比常在0.499～0.4999之间。第二，应力-应变曲线强烈非线性：小应变时较软，中等应变时因链段逐渐拉直而变硬（应变硬化），大应变时因接近极限伸长率而急剧变硬。此外还有Mullins效应（首次加载软化）、黏弹性（率相关）和温度敏感性。

二、超弹性本构的核心思想

超弹性模型不直接给出应力-应变关系，而是先定义应变能密度函数 W，再由 W 对应变求导得到应力。以常用的不可压缩形式为例，柯西应力可写为

σ = -p I + 2 (∂W/∂I₁) B + 2 (∂W/∂I₂) B⁻¹

其中：I₁、I₂ 为左 Cauchy-Green 变形张量的不变量
      B 为左 Cauchy-Green 变形张量
      p 为静水压力（由不可压缩约束确定）

不同模型只是 W 的展开形式不同。Neo-Hookean 只保留 I₁ 的一次项，参数少、小应变准，但大应变偏差明显；Mooney-Rivlin 增加 I₂ 项，能更好拟合中等应变；Yeoh 用 (I₁-3) 的三次多项式，对单轴拉伸数据拟合能力强，且能捕捉应变硬化；Ogden 以主伸长率 λ 为变量，适合大应变范围。选择模型时，关键不是"哪个更高级"，而是"你的零件实际处于什么应变区间、有哪些试验数据"。

三、Abaqus超弹性仿真流程概览

一条完整的超弹性分析链路可以概括为：试验数据 → 材料参数拟合 → 网格与单元选择 → 接触与边界 → 求解 → 结果判读。下图给出这条链路的信号与数据流向，以及每一步的关键输入和判断点。

四、一个可复现的入门示例

假设要分析一个截面直径3.5 mm、内径20 mm的丁腈橡胶O型圈，安装在矩形沟槽中，压缩率设定为20%。已知供应商提供了单轴拉伸数据（应变0～200%）和平面拉伸数据（应变0～100%）。

第一步，在Abaqus/CAE的Property模块中新建超弹性材料，选择Yeoh模型，把两组名义应力-应变数据分别以Uniaxial Test Data和Planar Test Data形式输入。点击Evaluate，检查在0～200%应变范围内模型是否满足Drucker稳定性，并观察拟合曲线与试验点的偏差。若平面拉伸段偏差超过10%，可改用Ogden（N=2或3）重新评估。

第二步，建立二维轴对称模型（O型圈和沟槽均可按轴对称处理），橡胶用CAX4H单元，沟槽用刚体解析面。设置摩擦系数0.1～0.2，采用罚函数接触。分析步打开NLGEOM，用位移控制把沟槽上盖下压到设计压缩量。

第三步，提交计算后重点看三件事：一是收敛历史，若出现大量Cutback，通常是网格畸变或接触过约束；二是支反力-位移曲线，它应与试验压缩曲线趋势一致；三是接触面上的法向接触压力分布，密封判据要求最大接触压力大于被密封介质压力，且接触宽度足够。

五、入门阶段最容易踩的三个坑

• 把泊松比设成0.5。 这会导致体积模量无穷大、求解直接失败。应设为0.499或更高，并配合混合单元。
• 只用单轴数据拟合。 单轴拉伸无法约束I₂项，用于压缩或剪切工况时误差可能很大。至少补充平面拉伸或双轴数据。
• 忽略Mullins效应。 若零件经历多次加卸载，首次循环的刚度明显高于后续循环。工程上可用Mullins Effect子选项或直接取稳定循环的试验数据拟合。