2. 橡胶材料特性与仿真本构基础
橡胶成形与模压仿真能否给出可信的合模力、飞边、缩孔与残余应力，关键不在网格有多细，而在于本构模型是否真实反映了橡胶的力学行为。橡胶与金属最大的区别是：它几乎不可压缩（体积模量约为剪切模量的千倍量级），应力—应变曲线高度非线性，加载与卸载不重合，并且对温度和加载速率极其敏感。本节围绕超弹性、粘弹性、Mullins 效应、温度依赖性四条主线，说明常见本构模型的物理意义与参数获取途径。
一、超弹性：把“非线性弹性”写成应变能函数
超弹性假设存在一个应变能密度函数 W，应力由 W 对变形的偏导给出，因此加载路径不影响最终应力状态，卸载可完全回复。工程上常用不变量 I₁、I₂、I₃ 或主伸长比 λ₁、λ₂、λ₃ 表达 W。由于橡胶近似不可压缩，通常令 I₃=1（即 J=1），把 W 拆成偏量部分与体积部分：W = Wdev + Wvol，其中 Wvol 用体积模量 K 惩罚体积变化，K 一般取 1000～2000 MPa 量级，仅用于约束不可压缩性，不参与剪切刚度的标定。
常见模型的物理意义与适用区间如下：
- Mooney-Rivlin（2 参数）：W = C₁₀(I₁−3) + C₀₁(I₂−3)。C₁₀、C₀₁ 分别对应网络变形与缠结贡献，剪切模量 μ = 2(C₁₀+C₀₁)。它形式简单、参数少，适合单轴拉伸 100%～200% 以内的中等变形，但对大变形和双轴数据拟合偏差明显。
- Ogden（N 阶）：以主伸长比幂函数求和，W = Σ 2μp/αp² (λ₁αp+λ₂αp+λ₃αp−3)。αp 可非整数，能同时拟合单轴、双轴、平面拉伸数据，是模压大变形（如厚制品压缩、翻边）的首选，但 N≥3 时参数相关性强，需多工况数据约束。
- Arruda-Boyce（8 链）：基于 8 链网络统计，W 由初始剪切模量 μ 与极限伸长比 λm 控制，物理意义清晰——λm 表征分子链有限伸展导致的应变硬化。它只需单轴数据即可初步标定，适合预测大变形硬化，但对小应变段精度略逊于 Ogden。
参数获取途径：优先做单轴拉伸（哑铃试样）、平面拉伸（纯剪）与等双轴拉伸三组试验，应变率控制在准静态（如 10⁻³ s⁻¹），取加载段数据。用最小二乘拟合时，先固定 K，再拟合偏量参数；若只有单轴数据，Mooney-Rivlin 与 Arruda-Boyce 可勉强使用，Ogden 则不建议超过 2 阶。拟合后必须检查：小应变剪切模量是否与 DMA 测得的储能模量一致，大应变段是否出现非物理的应力下降。
二、粘弹性：时间与频率的“记忆”
橡胶在恒定应变下应力会缓慢衰减（应力松弛），在恒定应力下应变会持续增长（蠕变），动态加载时储能模量 G′ 与损耗模量 G″ 随频率变化。仿真中通常用 Prony 级数描述剪切松弛模量：G(t) = G∞ + Σ G_i·exp(−t/τ_i)，其中 G∞ 为平衡模量，τ_i 为松弛时间。模压过程中，粘弹性主要影响硫化后冷却阶段的回弹、脱模瞬间的尺寸回复以及厚制品内部的残余应力分布。
参数获取：做 30 min 应力松弛试验，或 DMA 频率扫描（0.1～100 Hz，温度 25～150 ℃）。将曲线按对数时间等距取 5～8 个 τ_i，用最小二乘求 G_i。注意 τ_i 应覆盖实际工艺时间尺度——模压冷却通常几十秒到几分钟，若 τ_i 全部落在 10⁻³ s 量级，则对脱模回弹几乎没有贡献。温度依赖可用 WLF 方程：log aT = −C₁(T−Tref)/(C₂+T−Tref)，参考温度 Tref 一般取玻璃化转变温度以上 50 ℃ 左右，C₁、C₂ 由多温度 DMA 曲线平移得到。
三、Mullins 效应：首循环软化不可忽略
橡胶第一次加载到较大应变后再卸载，第二次加载的应力明显低于第一次，且随循环次数增加逐渐稳定，这就是 Mullins 效应（应力软化）。其物理来源是填料网络（炭黑、白炭黑）的破坏与部分重建。模压中，胶料在充模、压缩、排气阶段已经历多次变形，若仿真忽略 Mullins 效应，会高估合模力与制品回弹。
工程上常用伪弹性模型描述：定义损伤变量 η = 1 − (1/λm)m，其中 λm 为历史最大伸长比，m 为材料参数；卸载时模量按 η 折减。参数 m、r、β 通过 3～5 次循环加载试验拟合，取稳定循环的骨架曲线作为超弹性基准，再叠加损伤演化。若只做一次加载，可近似用首循环与稳定循环的应力差估算 m。
四、温度依赖性：模压工艺的核心变量
橡胶模压温度通常在 150～170 ℃，胶料在模具内经历升温、硫化、冷却三个阶段。温度对力学行为的影响体现在两方面：一是模量随温度升高而下降，可用 WLF 或 Arrhenius 移位；二是硫化反应使交联密度增加，剪切模量随硫化度上升。仿真中若只做等温超弹性分析，会低估充模阶段的流动阻力，也会错过硫化度对最终模量的贡献。
实操建议：先用橡胶加工分析仪（RPA）测不同温度下的硫化曲线，得到 t₉₀ 与交联密度变化；再把硫化度作为状态变量，将剪切模量写成 μ(α) = μ₀ + (μ∞−μ₀)·αn，α 为硫化度。这样在模压仿真中，充模阶段用低模量、低粘度的未硫化参数，保压阶段逐步过渡到高模量硫化胶参数，合模力曲线才能与实测吻合。
五、模型选择速查
| 模型 | 核心参数 | 适用场景 | 数据需求 |
|---|---|---|---|
| Mooney-Rivlin | C₁₀, C₀₁ | 中小变形、快速估算 | 单轴拉伸 |
| Ogden (N=3) | μ_p, α_p | 大变形、多轴模压 | 单轴+平面+双轴 |
| Arruda-Boyce | μ, λ_m | 大变形硬化预测 | 单轴拉伸 |
| Prony + WLF | G_i, τ_i, C₁, C₂ | 冷却回弹、残余应力 | 应力松弛/DMA |
| 伪弹性 Mullins | m, r, β | 多次循环、排气后充模 | 循环加载 3～5 次 |
需要强调的是，本构模型不是越复杂越好。对于薄壁密封件，Mooney-Rivlin 加 Prony 已足够；对于厚壁减震块或大型模压制品，Ogden 三阶叠加 Mullins 与硫化耦合才必要。参数获取的试验成本与仿真精度之间要权衡，但无论选哪个模型，都必须用至少一组独立工况验证，否则本构参数只是“曲线拟合的装饰”。