1. 橡胶超弹性仿真工程导论

橡胶密封圈、减震支座、轮胎胎面、医用硅胶导管——这些零件有一个共同点：受力后会发生大幅可恢复变形，卸载后几乎回到原形。金属结构仿真中常用的线弹性本构在这里完全失效，因为橡胶的应力-应变关系既非线性、又几乎不可压缩，还强烈依赖加载历史。Abaqus 之所以成为橡胶件工程仿真的主流工具，核心就在于它内置了完整的超弹性本构族、混合单元技术和稳定的接触求解能力。本章不急于建模，而是先把"橡胶为什么特殊、超弹性理论在算什么、Abaqus 在其中承担什么角色"这三件事讲清楚，并给出整门课的学习地图。

1.1 橡胶材料的三重特殊性

理解橡胶仿真，先要理解它与金属在力学行为上的根本差异。第一是大变形与几何非线性：橡胶件工作应变常在 100%~500%，例如一个压缩型减震块在额定载荷下高度可能压缩 20%~30%，此时应变与位移不再近似线性，必须使用有限变形理论，Abaqus 中体现为打开几何非线性（NLGEOM）。第二是近乎不可压缩：橡胶体积模量远大于剪切模量，泊松比接近 0.5（典型 0.499~0.4999），体积在变形中近似守恒。这带来一个数值难题——纯位移单元会因体积自锁而刚度过高，因此 Abaqus 推荐使用杂交单元（如 C3D8H、C3D10H）或缩减积分单元配合沙漏控制。第三是本构非线性与历史依赖：单轴拉伸曲线呈明显的 S 形，小应变时较软、中等应变时变硬、大应变时再次变硬；同时存在 Mullins 效应（首次加载比后续加载更硬）和黏弹性松弛。工程上通常先用超弹性描述平衡态弹性，再按需叠加黏弹性或 Mullins 修正。

这三重特殊性直接决定了仿真流程的每一个选择：单元类型、材料参数标定方式、加载步长、收敛判据，都与金属件仿真不同。

1.2 超弹性理论在算什么

超弹性（Hyperelasticity）的核心假设是：材料存在一个应变能密度函数 W，应力可由 W 对变形的偏导直接得到，因此加载路径无关、无能量耗散。工程中最常用的是以不变量表示的模型：

• Neo-Hookean：只有一个剪切模量参数，形式最简单，仅适合小到中等应变（约 100% 以内）的粗略估算。
• Mooney-Rivlin：两参数（C10、C01），能较好拟合中等应变，是早期工程常用模型，但大应变段偏差明显。
• Yeoh：三参数（C10、C20、C30），只依赖第一不变量，对单轴拉伸数据拟合能力强，且能捕捉大应变硬化，是橡胶件仿真的主力模型之一。
• Ogden：以主伸长率表示，阶数可调，适合多轴加载和复杂工况，但参数标定需要更丰富的基础试验数据。
• Arruda-Boyce：基于分子链统计，物理意义清晰，适合需要外推到不同加载模式的场景。

需要强调的是：超弹性模型本身不描述永久变形、松弛和 Mullins 效应。如果零件长期受压后出现残余变形，或需要评估密封件的长期密封力衰减，就必须在超弹性基础上叠加黏弹性（如 Prony 级数）或 Mullins 模型。这是初学者最容易忽略的一点。

1.3 Abaqus 在橡胶仿真中的工程价值

Abaqus 对橡胶仿真的支撑体现在四个层面。其一是本构库完整：上述模型全部内置，且支持从试验数据直接拟合参数，用户只需输入单轴、双轴、平面剪切试验曲线，Abaqus 即可用最小二乘拟合出材料常数并给出稳定性检查。其二是单元技术成熟：杂交单元、杂交缩减积分单元、以及针对大变形优化的 C3D10H 等，能有效缓解体积自锁。其三是接触与自接触能力强：橡胶密封件往往涉及大变形下的自接触、密封面接触压力分布，Abaqus 的通用接触和面-面接触配合自适应网格可稳定处理。其四是求解器鲁棒：隐式 Standard 适合静态密封分析，显式 Explicit 适合冲击、跌落和复杂接触，两者可联合使用。

1.4 课程学习地图与实战目标

本课程围绕"能独立完成一个橡胶件仿真并判断结果是否可信"这一目标展开，整体分为四个阶段：材料参数获取与标定、几何与网格策略、接触与加载设置、结果判读与工程决策。下图给出课程主线与各阶段的关键交付物，帮助你在学习过程中随时定位自己处在哪一环。