Main points
☑ 常见超弹性本构模型详解
<
☑ 模型对比与选型决策
<
☑ 参数标定的技术要点
<
☑ 常见误区与注意事项
<
导语
做橡胶有限元分析，本构模型选错了，仿真结果可能差出30%甚至更多。但很多工程师在选模型时，大都凭经验拍脑袋，很少系统梳理过各模型的适用边界。
Mooney-Rivlin、Ogden、Yeoh、Neo-Hookean、Arruda-Boyce……这些模型背后的数学形式是什么？各有什么适用范围？选错了会出什么问题？
这篇文章从应变能密度函数出发，介绍常见几种超弹性本构模型的理论基础、数学表达式和选用逻辑。
01
应变能密度函数：
大变形本构模型的数学基础
橡胶材料的宏观力学行为表现为显著的非线性和大变形，其基础力学本构关系并非直接通过应力-应变函数来定义，而是由一个更为根本的热力学势函数——应变能密度函数W所决定。该函数代表了单位参考构型体积内因变形而储存的可恢复弹性能量。
对于各向同性的超弹性材料，根据物质客观性原理（即材料响应不应依赖于观察者的刚体运动），应变能密度函数 W 必须是客观应变度量的标量函数。通常选择右柯西-格林变形张量 C=FTF或其等价形式左柯西-格林张量 B=FFT作为自变量，其中 F为变形梯度张量。
进一步地，由于橡胶材料在未变形状态下无特定方向偏好，即具有各向同性，W 可以简化为 C（或 B）的三个主不变量I₁ 、I₂和 I₃的函数，即W=W(I₁ ,I₂,I₃)。其中，
第一不变量：
反应总拉伸程度。
第二不变量：
第三不变量：
则直接关联体积变化率J=det（F）
其中λ₁、λ₂、λ₃是三个主方向的伸长比，与应变ɛ的关系为：λ=ɛ+1
鉴于橡胶在变形过程中体积变化极小（通常小于1%），工程分析中普遍采用不可压缩假设，即J=1，从而I₃=1。这一约束将应变能函数简化为仅依赖于前两个不变量的形式：W=W(I₁ ,I₂)。
不同超弹性本构模型的区别，本质上是应变能密度函数W(I₁, I₂)的数学形式不同。
02
常见本构模型详解
01
Neo-Hookean 模型
最简单的超弹性模型，应变能密度函数仅依赖于应变第一不变量，具体数学形式为：
其中C₁₀是材料常数，与剪切模量G相关：G = 2C₁₀。
特点：
只有一个参数，标定简单，计算稳定
适用于小到中等应变（<30%）
无法描述填充型橡胶的应力软化和各向异性
对单轴拉伸拟合较好，双轴拉伸偏差较大
适用场景：
初步设计分析、变形较小的密封件和减震件、概念验证阶段。
02
Mooney-Rivlin 模型
经典双参数模型，应变能密度函数同时包含I₁和I₂项：
C₁₀和C₀₁是两个独立的材料常数。初始剪切模量G = 2(C₁₀ + C₀₁)。
特点：
两个参数，拟合能力优于Neo-Hookean
中等应变（30%-100%）区间表现稳定
多数有限元软件默认支持，工程应用成熟
大应变（>100%）时预测精度下降
适用场景：
轮胎、橡胶衬套、悬置、一般减震橡胶，应变范围30%-100%，且用于拟合的试验数据必须包含等双轴拉伸应力应变数据，否则容易出现拟合外推失真。
03
Yeoh 模型
多项式模型的简化形式，也称缩减多项式，仅用第一不变量I1的高阶项描述大应变行为：
通常取三项（三参数Yeoh模型），C₁₀、C₂₀、C₃₀为材料常数。其中C₂₀通常为负值，体现中等应变下的软化过程；C₃₀控制大应变下的曲线上扬，捕捉橡胶的刚化效应。
特点：
大应变拟合能力强，可描述S型应力-应变曲线（刚化效应）
计算效率高，收敛性好，适合大变形问题
中等应变区间可能不如Mooney-Rivlin稳定
适用场景：
悬置、大变形密封件、橡胶衬套大载荷工况，应变超过100%时优先考虑，对等双轴拉伸数据的依赖程度比Mooney-Rivlin模型小。
04
Ogden 模型
直接用主伸长比表示的应变能密度函数：
μᵢ和αᵢ是材料常数，N为项数（通常取N=3或N=6）。
特点：
数学形式最灵活，拟合能力最强
可处理极大应变（>300%）和复杂加载路径
参数多（6-12个），标定难度大
计算成本高，可能出现数值不稳定
适用场景：
需要高精度拟合的工程分析，如动力总成悬置、轮胎大变形分析。当其他模型拟合精度不足时，Ogden N=3或4是可靠的备选方案，但如果用于拟合模型参数的试验数据中缺少等双轴拉伸试验数据或者等双轴拉伸数据的应变范围不足以覆盖有限元仿真的完整应变范围，可能会出现拟合失真的情况。
05
Arruda-Boyce 模型
基于分子网络统计力学推导的物理模型，考虑分子链有限伸展性：
μ为初始剪切模量，λₘ为网络伸展极限，Cᵢ为级数系数（固定值，非拟合参数）。
特点：
参数有明确物理意义（μ为剪切模量，λₘ为分子链伸展极限）
大应变预测能力强，外推性比经验模型好
参数少（2个），标定相对简单
适用场景：
大变形工况推荐模型，轮胎、衬套、密封件均适用。，但如果用于拟合模型参数的试验数据中缺少等双轴拉伸试验数据，或者等双轴拉伸数据的应变范围不足以覆盖有限元仿真的完整应变范围，可能会出现拟合失真的情况。
03
模型对比与选型决策
各模型的核心差异在于应变能密度函数的构造方式，直接决定了适用范围。
模型参数与适用范围对比：
选型决策流程
按变形程度划分推荐优先级：
小应变（<30%）：
Neo-Hookean → Mooney-Rivlin中等应变（30%-100%）：
Mooney-Rivlin 为主大变形（>100%）：
Ogden、Arruda-Boyce 或 Yeoh极大应变或复杂加载：
Ogden、Arruda-Boyce（材料本构关系测试的应变范围必须覆盖产品最大应变）
实际工程中，建议同时用多个模型拟合，对比刚度曲线与实测曲线的吻合程度，选误差最小的模型，经过对比验证的本构模型可用于其他同规格产品的仿真。
04
参数标定的技术要点
选对模型只是第一步，参数标定精度决定仿真结果的可信度。
实验数据覆盖目标应变范围
标定应变范围必须覆盖或超过产品实际工况的最大应变。用25%应变的数据去标定大应变模型，预测结果不可信。下图所示L1应变等级下的拉伸曲线与L4的拉伸曲线存在显著差异。
下图所示在材料试验数据范围之外，部分模型拟合结果发散不可用，假设使用了发散的本构模型进行了仿真计算，仿真误差可达50%以上
多种加载模式联合标定
单轴拉伸+等双轴拉伸+平面拉伸数据联合标定，比单一模式可靠。如下图所示，仅使用单轴拉伸数据无法准确表征应变能密度函数的空间分布
Mullins效应的处理
对于填充型橡胶材料，Mullins效应表现为：在经历首次加载后，后续加载曲线低于初始加载曲线，呈现宏观上的软化现象，且材料的软化程度取决于其此前经历过的最大变形量（即“峰值拉伸历史”）。通常采用循环加载消除Mullins效应后，取稳定加载或卸载段应力应变数据标定超弹本构模型参数。
但如果要在仿真中考虑材料的Mullins效应（例如在疲劳仿真中，需要准确计算不同位置的橡胶单元应力应变）还需额外拟合Mullins效应模型参数（Ogden-Roxburgh模型），此种情况下，通常先利用材料的主拉伸试验曲线数据拟合超弹性本构模型参数，再利用不同应变条件下的循环拉伸曲线拟合材料的Mullins参数。
05
常见误区与注意事项
复杂模型参数多，若实验数据不足或质量不高，反而容易过拟合，预测能力下降，尤其是某些超弹性本构模型，在缺少多轴变形试验数据，或者试验数据应变范围不足的情况下，可能会出现拟合结果在大应变下发散的情况，导致仿真不收敛或者仿真结果严重失真。选择模型的原则是：在满足精度要求的前提下，优先选择参数少的模型。
同一材料在不同应变范围、不同加载模式下，最优本构模型可能不同。建议根据具体工况选择模型。
橡胶测试离散性大，建议材料试验重复3-5次，取均值并剔除异常点。数据质量差是标定失败的主要原因之一。
06
写在最后
本构模型选择没有标准答案，只有合适与否。
理解应变能密度函数的物理意义，明确产品工况需求，用高质量实验数据标定，才能做出合理判断。
如果你正在做橡胶仿真，或者在本构模型选择、参数标定上遇到具体问题，欢迎点击文末“阅读原文”联系我们。我们长期为轮胎、衬套、减震件等橡胶制品企业提供材料测试和仿真支持。
END
