超弹性 — CAE术语解释
超弹性（Hyperelasticity）
老师，我要做橡胶衬套的分析，被告知要选择"超弹性"材料模型。这和普通弹性有什么区别吗？
超弹性的理论基础
超弹性的定义和特征
"超弹性"和普通橡胶弹性有什么区别？看起来是一样的啊。
根本的区别在于变形的机制和变形的大小。金属等线性弹性体是通过晶格结构的歪曲产生变形，通常应变小于1%。而橡胶等超弹性体是通过聚合物链的绞缠解开、链本身伸缩来产生变形，应变可以超过100%，这是常见的。应力-应变关系呈非线性，而且有一个强约束条件"体积几乎不变（不可压缩性）"。
不可压缩性是说体积应变始终为0？怎样用公式表示呢？
正是这样。变形前后的体积比
应变能密度函数具体是什么样子？最著名的是什么？
实际应用中最常用的是"Mooney-Rivlin模型"和"Ogden模型"。Mooney-Rivlin的2项模型形式为
超弹性的数值计算手法
不可压缩性的处理和要素技术
有限元法中怎样处理不可压缩的约束？直接这样好像解不了吧？
眼光敏锐。严格的不可压缩约束会导致体积应变为0的条件成为压力p的约束方程，这在标准位移法中会引起"锁定"现象，数值不稳定。解决这个问题的主要方法有两种："混合定式化"和"减积分"。混合定式化中，位移和压力作为独立自由度处理，例如Abaqus的"C3D8H"单元（H代表Hybrid）增加了压力自由度，用弱形式满足条件。
"减积分"具体是做什么的？
计算体积应变（=位移微分的迹）时使用的积分点数少于计算变形时的积分点数。例如，Abaqus的"C3D8RH"单元在计算变形时用8个点（低减积分），计算体积应变时只用1个点（中心）。这样可以缓解过度的体积约束，防止锁定。但这样会引发另一种不稳定模式"沙漏模式"，需要稳定化技术。
超弹性分析中求解器设置需要特别改变吗？
必须的。大变形和强材料非线性需要迭代求解法（Newton-Raphson法）。但初始刚度不恰当会导致发散。因此"步长控制"很重要。Ansys有"Stabilization"功能，Abaqus有"Automatic stabilization"，在不稳定的接触和大变形初期阶段加入虚拟阻尼力来帮助收敛。通常要把荷载或位移分成小增量，每个增量的最大迭代次数设为15~20次。
超弹性的实务应用
材料试验数据的参数同定
材料常数（C10, C01等）怎样确定？制造商手册里有吗？
很少有。通常需要从实验数据进行"参数同定（拟合）"。最少需要"单轴拉伸试验"的数据。要追求更高精度，可以加入"单轴压缩""平面拉伸（双轴拉伸）""剪切"的数据。实验数据覆盖的变形模式越多，所选材料模型的预测精度保证的范围就越广。比如，硅胶O形圈的分析，用单轴拉伸和平面拉伸数据，用Mooney-Rivlin 2项模型拟合，是一种定式做法。
如果只有单轴拉伸的实验数据，该用什么模型？
这种情况下，Neo-Hookean模型（Mooney-Rivlin的特殊形式，C01=0）或Yeoh模型较为现实。Yeoh模型
拟合时应变范围该取多大？
至少要覆盖分析中预计的最大应变，并留出余量。例如，密封部件的名义最大应变预计为50%，实验数据至少应取到应变100%。因为材料模型只在拟合数据的范围内具有高可靠性，范围外是外推，有很大风险偏离预测。另外，为了准确捕捉应变0%附近的初始刚性，建议在应变10%以下的数据点密集采集。
超弹性的软件比较
主要软件的超弹性对应
Ansys、Abaqus、COMSOL在超弹性分析功能上有大差别吗？
基本的材料模型（Mooney-Rivlin, Ogden, Yeoh等）都具备，但实现细节和易用性有差异。Abaqus在轮胎分析等橡胶制品工业应用中很强，超弹性单元（C3D8H, C3D10H等）和稳定化算法都非常成熟。Ansys Mechanical也有同样的功能，但材料数据输入界面更直观，可以直接输入实验数据，用多个模型同时进行拟合比较，这一点工具优势明显。
COMSOL是什么位置？
COMSOL Multiphysics的最大优势不是"超弹性"单独的物理，而是它易于与其他物理场（热、电、流体）耦合的特点。例如，发热伴随的橡胶滞后现象（Mullins效应）或电场引起的介电弹性体大变形等问题容易处理。材料模型丰富，可以直接输入用户自定义的应变能函数，灵活性很高。但在纯粹大规模结构接触问题的求解效率上，比Abaqus略逊一筹。
免费或低价CAE软件（CalculiX, Code_Aster等）能做超弹性分析吗？
可以。例如，支持Abaqus兼容输入格式的CalculiX支持超弹性材料模型（HYPERELASTIC）和混合单元（C3D8H等）。Code_Aster也实现了Mooney-Rivlin和Ogden模型。但商用软件那样的高级自动收敛控制和图形化材料拟合工具基本没有。从实验数据求材料常数通常需要用Python脚本或专门工具（如MCalibration）预先处理。学习成本高，但对深入理解原理很有帮助。
超弹性的故障对应
收敛不良和对策
超弹性分析中，开始几步就突然"收敛失败"报错。为什么？
最多见的原因有两个："初始条件不适当"和"接触不稳定"。超弹性体初始状态几乎没有刚性。部件间已经有初始接触穿透、或有大的初始荷载，求解器找不到合适的位移增量就会发散。首先要精确调整初始接触状态，第一步采用很小的位移（比如0.1mm）或荷载，让部件稳定接触。
出现了"负体积错误"。这是什么原因？
单元严重扭曲，体积在数学上变成负值（单元翻转）的状态。在大变形分析中低阶单元特别容易发生。直接原因有：①网格过粗无法追随变形，②接触条件太急陡单元被压扁，③材料模型不稳定（例如Ogden模型参数α为负，某应变处应力下降），等等。对策是细化网格、提高单元阶数（切换到2阶单元）、恰当设置接触刚性、验证材料模型稳定性。
分析完成了，但反力和接触压力结果在振荡。这也是网格问题吗？
是的。这是"沙漏模式"的典型症状，特别是用低减积分单元时发生。单元内部变形模式有零能量（不产生应变能的变形），在数值中被放大了。对策是：①切换到完全积分单元（例如从Abaqus的C3D8H换到C3D8），②启用"沙漏控制"（Ansys或Abaqus选项），③细化网格（现象缓解但难以根除），等等。完全积分单元计算成本增加，但结果稳定性大幅提升。
把材料模型改成Ogden（N=3）后，超过某个应变应力急剧下降，这是物理不合理的结果。是bug吗？
不是bug，是"材料不稳定性"。Ogden模型的参数
相关主题
更多细节
错误