Mooney-Rivlin超弹性模型
Mooney-Rivlin超弹性的理论基础
超弹性简介
老师，「超弹性」是橡胶类材料的模型吗？
是的。超弹性（hyperelasticity）是大变形后能恢复的材料。橡胶、弹性体、生体组织。应力由应变能函数 $W$ 的导数决定。
Mooney-Rivlin模型
Mooney-Rivlin是应用最广泛的超弹性模型。应变能为：
$I_1, I_2$ 是变形张量的不变量。$C_{10}, C_{01}$ 是材料常数。
只用2个常数就能表示橡胶的大变形吗？
应变100%以内精度较好。200%以上时用Ogden模型或Arruda-Boyce模型更精确。Mooney-Rivlin是「最简单、应用最广泛」的模型。
Abaqus
```
*HYPERELASTIC, MOONEY-RIVLIN
C10, C01, D1
```
$D1$ 是体积模量（非压缩性情况下 $D1 \to 0$）。
Nastran
```
MATHE, 1, MOONEY
, C10, C01, , , , , D1
```
材料试验
Mooney-Rivlin常数由以下试验确定：
- 单轴拉伸试验 — 最基本
- 双轴拉伸试验 — 更精确
- 平面拉伸（纯剪切） — 补充
Abaqus中通过*HYPERELASTIC, TEST DATA输入试验数据，自动拟合。
总结
要点：
- $W = C_{10}(I_1-3) + C_{01}(I_2-3)$ — 2常数的应变能
- 应变100%以内精度较好 — 更大变形用Ogden等
- 注意非压缩性 — 必须用混合元素（C3D8RH等）
- 由试验数据拟合 — *HYPERELASTIC, TEST DATA
Mooney的1940年论文和Rivlin的1948年扩展
Melvin Mooney在1940年发表论文，首次用2参数（C₁、C₂）表述橡胶等向超弹性时，有限变形的数学理论体系还在建设阶段。1948年Ronald Rivlin给出了用不变量I₁、I₂、I₃的通用应变能函数展开理论，证明了Mooney函数是其最低阶近似。这就是「Mooney-Rivlin」双人名模型的由来。
Mooney-Rivlin超弹性的数值计算方法
超弹性的FEM实现
超弹性材料FEM的注意事项：
1. NLGEOM=YES必须 — 橡胶是大变形
2. 混合元素 — 橡胶近似非压缩（$\nu \approx 0.5$）。体积锁定对策
3. $C_{10}, C_{01}$ 的拟合 — 由试验数据
用哪种元素？
总结
C₁、C₂拟合的实际应用
Mooney-Rivlin参数C₁和C₂通常通过单轴拉伸、平面应变、等双轴这3种试验数据的同时最小二乘拟合来确定。天然橡胶（NR）的典型值为C₁≈0.16 MPa、C₂≈0.04 MPa，对伸长比λ=3左右能很好拟合。Abaqus的"Evaluate"功能可一屏显示各试验模式的预测曲线，还能自动检查参数是否稳定（Drucker稳定性）。
Mooney-Rivlin超弹性的实务应用
超弹性的实务
O形圈、轮胎、橡胶衬套、防振支座、医疗器械等。
材料常数的确定
稳定性检查
Mooney-Rivlin在所有变形模式（拉伸、压缩、剪切）下必须具有正刚度。拟合结果可能出现不稳定（负刚度）。用Abaqus的STABILITY CHECK检查。
实务检查清单
轮胎侧壁分析的主力模型
汽车轮胎侧壁（天然橡胶基复合材料）的大变形分析中，Mooney-Rivlin 2参数模型至今仍是业界标准的首选。普利司通、米其林的FEM分析部门从1990年代开始采用Abaqus+Mooney-Rivlin组合分析接地压分布和变形形状，已建立能将误差控制在5%以内的拟合流程。对于λ≧4的大伸长范围，推荐切换到Ogden模型。
Mooney-Rivlin超弹性的软件比较
超弹性的工具
选择指南
Marc和MSC Nastran的超弹性实现
MSC Marc与Abaqus并列是橡胶分析的老字号，Mooney-Rivlin从Marc第1代（1971年发布）就已实现。Marc 2020及以后版本新增了"Hyperelastic Fitting Tool"，可从CSV读入试验数据自动拟合。MSC Nastran传统上被认为不适合橡胶分析，但2016年起强化了非线性求解器SOL 400的超弹性材料支持，现在可在Marc混合环境中使用。
Mooney-Rivlin超弹性的先端研究
超弹性的先端研究
非压缩约束与罚函数法
橡胶几乎不产生体积变化（泊松比≈0.499）是非压缩性材料，用位移法FEM会产生体积锁定。对策是用「近似非压缩（Nearly Incompressible）」表述，把体积模量K作为罚项加入，通常设K/μ≈1,000~10,000。Abaqus的混合元素（C3D8H等）采用独立自由度方式处理静水压，这样就避免了这个问题，因此橡胶分析强烈推荐使用"H"后缀元素。
Mooney-Rivlin超弹性故障排查
超弹性的故障
不稳定参数的陷阱
Mooney-Rivlin模型中可能出现 C₂<0 的拟合结果。当C₂为较大负值时，高伸长域的应力-伸长曲线无法单调增加（Drucker不稳定），分析会出现反转、发散危险。运行Abaqus的"Material Evaluation"的Stability check功能会显示不稳定应变范围。实务中的对策是只在不含不稳定域的伸长范围内使用，或改用Ogden、Yeoh模型。
相关话题
更详细
错误