第三章：材料模型选择——弹性、弹塑性、超弹性与各向异性
做插拔力仿真这么多年，我踩过最大的坑，就是材料模型选错了。你想想看，明明网格画得漂亮，接触设置也对了，结果算出来的力曲线跟实测差了十万八千里。十有八九，问题出在材料上。
这一章，咱们就把材料模型这事儿彻底聊透。我会把四种最常见的材料模型——弹性、弹塑性、超弹性、各向异性——掰开揉碎了讲。每个模型适合什么场景、怎么定义参数、有哪些坑，我都会结合自己的项目经验来说。
核心观点：材料模型的选择，决定了仿真结果的“可信度”。选对了，事半功倍；选错了，一切白费。
3.1 弹性模型——最基础，也最容易用错
弹性模型，说白了就是胡克定律的有限元版本。应力跟应变成正比，卸载后完全恢复。这个模型只有两个参数：弹性模量 E 和泊松比 ν。
什么时候用？
- 材料变形很小，远没到屈服点
- 你只关心插拔力的初始阶段
- 做快速预分析，先看看趋势
我的经验：很多新手一上来就用弹性模型算插拔力，结果算出来的力比实测大两三倍。为什么？因为实际接触过程中，微凸体早就发生塑性变形了。我个人习惯是：只有做概念设计阶段的快速对比时，才用纯弹性模型。
参数怎么给？
*MAT_ELASTIC
$ MID E NU RO
1 2.1e5 0.30 7.85e-9注意单位！我见过最离谱的错误，是把 E 的单位搞混了。钢的 E 是 210 GPa，也就是 2.1e5 MPa。你要是写成 2.1e5 Pa，那算出来的变形量会大得离谱。
3.2 弹塑性模型——插拔力仿真的主力军
说实话，我做过的连接器插拔力仿真，90% 以上用的都是弹塑性模型。因为端子接触过程中，局部应力很容易超过屈服极限，产生永久变形。
弹塑性模型的核心要素：
- 屈服准则：最常用的是 von Mises 准则。什么时候屈服？等效应力达到屈服应力 σ_y 的时候。
- 硬化法则：各向同性硬化还是随动硬化？对于连接器这种单向加载的情况，各向同性硬化就够了。
- 应力-应变曲线：这是最关键的输入。别用理想弹塑性，那太粗糙了。
避坑指南：我曾经在一个项目中，直接用材料手册上的屈服强度来定义弹塑性模型。结果仿真出来的插拔力偏小很多。后来才发现，手册上的屈服强度是 0.2% 残余应变对应的值，而实际接触区域的应变可能已经到 5% 甚至 10% 了。正确的做法是：用完整的真实应力-应变曲线。
真实应力-应变曲线怎么转换？
工程应力 σ_eng 和工程应变 ε_eng 需要转换成真实应力 σ_true 和真实应变 ε_true：
σ_true = σ_eng * (1 + ε_eng)
ε_true = ln(1 + ε_eng)这个转换，我建议你写个小脚本自动算。手动算的话，十个数据点以内还行，多了容易出错。
弹塑性参数定义示例：
*MAT_PLASTIC_KINEMATIC
$ MID RO E PR SIGY ETAN BETA
2 7.85e-9 2.1e5 0.30 350 1500 0.0
$ SRC SRP FS
0.0 0.0 0.5这里 SIGY 是屈服强度，ETAN 是切线模量。BETA=0 表示各向同性硬化。
3.3 超弹性模型——橡胶和密封件的专属
连接器里经常有密封圈、防水垫圈这类橡胶件。它们的变形很大，而且是非线性的。这时候弹性模型和弹塑性模型都不好使了，得用超弹性模型。
超弹性模型的特点：
- 基于应变能密度函数来描述材料行为
- 可以承受大变形（应变超过 100%）
- 卸载后完全恢复，没有永久变形
常用的超弹性模型：
| 模型名称 | 适用场景 | 参数数量 |
|---|---|---|
| Mooney-Rivlin | 中等变形（应变 < 100%） | 2 个（C10, C01） |
| Ogden | 大变形（应变 > 100%） | N 阶（通常 N=1~3） |
| Neo-Hookean | 小到中等变形 | 1 个（C10） |
| Yeoh | 填充橡胶、炭黑橡胶 | 3 个（C10, C20, C30） |
注意：超弹性模型的参数不能随便从文献里抄。不同批次的橡胶，硬度 Shore A 一样，但应力-应变曲线可能差很多。我建议你每次做项目前，都让供应商提供该批次橡胶的单轴拉伸测试数据，然后自己拟合参数。
参数拟合示例：
! 从单轴拉伸数据拟合 Mooney-Rivlin 参数
! 测试数据：应变 0.1, 0.2, 0.3, 0.5, 0.8, 1.0
! 对应应力：0.5, 1.2, 2.0, 3.8, 6.5, 9.0 MPa
*MAT_HYPERELASTIC_RUBBER
$ MID RO
3 1.2e-9
$ C10 C01 C20 C30 ...
0.35 0.12 0.0 0.0拟合的时候，我习惯用最小二乘法。很多有限元软件自带拟合工具，比如 Ansys 的 Curve Fitting 模块，或者 Abaqus 的 Material Evaluation。别手算，那太折磨人了。
3.4 各向异性材料定义——方向对了，结果才对
连接器里有些材料不是各向同性的。比如磷青铜的轧制板材，沿着轧制方向和垂直于轧制方向的力学性能就不一样。再比如一些复合材料端子，纤维方向和非纤维方向的强度差很多。
什么时候需要考虑各向异性？
- 材料有明显的轧制纹理
- 复合材料（如碳纤维增强塑料）
- 经过冷加工硬化的材料
- 单晶材料（如某些 MEMS 连接器）
各向异性材料的定义方式：
最常用的是正交各向异性模型。需要定义三个方向的弹性模量 E1、E2、E3，三个泊松比 ν12、ν23、ν31，以及三个剪切模量 G12、G23、G31。一共 9 个独立参数。
*MAT_ORTHOTROPIC_ELASTIC
$ MID RO E1 E2 E3
4 7.85e-9 2.1e5 1.8e5 1.8e5
$ PR12 PR23 PR31 G12 G23 G31
0.30 0.35 0.35 7.5e4 6.5e4 6.5e4我的经验：定义各向异性材料时，最头疼的是获取这些参数。供应商通常只给一个方向的 E 和 ν。我的做法是：自己拿材料去做拉伸测试，至少测三个方向（0°、45°、90°）。如果条件不允许，那就保守一点，用各向同性模型，但把安全系数放大到 1.5 倍。
材料坐标系怎么对齐？
这个细节很多人会忽略。你定义了各向异性参数，但材料坐标系跟全局坐标系没对齐，那算出来就是错的。在大多数软件里，你需要通过定义材料方向（Material Orientation）来指定局部坐标系。
! 在 LS-DYNA 中定义材料方向
*DEFINE_COORDINATE_VECTOR
$ CID X1 Y1 Z1 X2 Y2 Z2
1 1.0 0.0 0.0 0.0 1.0 0.0
! 然后将该坐标系赋给单元
*SET_PART
$ PID
1
*PART
$ PID SECID MID ...
1 1 4
$ ... ... CID
1嗯，这里要注意：材料方向定义错了，仿真结果会完全偏离实际。我见过一个案例，工程师把轧制方向定义反了，结果算出来的插拔力比实测小了 40%。
3.5 材料模型选择的决策流程
说了这么多，到底怎么选？我总结了一个简单的决策流程：
- 先看变形量：应变 < 0.2%？用弹性模型。应变 > 0.2%？往下看。
- 再看材料类型：橡胶或类似橡胶？用超弹性模型。金属材料？用弹塑性模型。
- 检查各向异性：材料有没有明显的方向性？有的话，用各向异性模型。没有的话，各向同性就够了。
- 确认参数来源：参数是实测的还是估算的？实测的放心用，估算的加安全系数。
最后说一句：材料模型没有绝对的对错，只有合适不合适。同一个连接器，你用弹塑性模型算出来的插拔力是 5N，用超弹性模型算出来可能是 8N。哪个对？得看你的材料到底是什么。所以，我建议你每次做仿真前，先花点时间搞清楚材料的真实力学行为。这一步省了，后面全是白费功夫。
公众号:蓝海资料掘金营，微信deep3321