第3章材料本构模型

3. 材料本构模型：金属弹塑性模型、塑料粘弹性模型与摩擦系数

各位同学，大家好。这一章我们聊聊材料本构模型。说实话，这是插拔力仿真里最容易「翻车」的地方。我见过太多人，网格画得漂亮，接触设置也对了，结果算出来的力曲线跟实测差了十万八千里。为什么？材料参数给错了。

插拔力仿真涉及的材料，说白了就三大类：金属（端子）、塑料（外壳/绝缘体）、以及接触面的摩擦行为。咱们一个一个来拆解。

端子材料，比如磷青铜、铍铜，在插拔过程中一定会进入塑性。你想想看，插拔力仿真如果不考虑塑性，那正压力算出来会偏大，插拔力自然也不准。

我个人习惯用双线性等向强化模型。为什么？简单、稳定、参数好标定。

双线性模型核心参数：

在 Abaqus 里，定义方式是这样的：

*Material, name=Copper_Alloy
*Elastic
 115000, 0.34
*Plastic
 450, 0.0
 520, 0.02

注意看，第一行是屈服点（450MPa，塑性应变为0），第二行是某个塑性应变下的应力。双线性嘛，两点确定一条直线。

我的经验： 切线模量 E_t 一般取弹性模量的 1%~5%。我曾经遇到一个案例，某学生把 E_t 设成了 10%，结果回弹量算出来特别大，跟实测对不上。后来一查，材料实际只有 2% 的硬化率。

如果材料有更复杂的硬化行为，比如不锈钢，我会用多线性模型。说白了，就是把应力-应变曲线拆成好几段直线来拟合。

*Plastic
 300, 0.0
 380, 0.01
 450, 0.03
 500, 0.06

这里有个坑：数据点一定要单调递增。我见过有人把数据输反了，结果 Abaqus 报错「The values must be monotonically increasing」——嗯，别问我怎么知道的。

避坑指南： 我曾经用多线性模型算一个连接器，结果力曲线出现锯齿状波动。查了半天，发现是应力-应变曲线在某个点斜率突变太大。解决办法：增加数据点密度，或者改用平滑的 Ramberg-Osgood 模型。

塑料件，比如 LCP、PA9T，在插拔过程中会表现出粘弹性。什么意思？就是它的力学行为跟时间有关。你插得快和插得慢，测出来的力不一样。

我建议用线性粘弹性模型，配合 Prony 级数来定义。在 Abaqus 里，你需要定义：

举个例子：

*Viscoelastic, time=PRONY
 0.3, 0.1, 0.2
 0.2, 1.0, 0.1
 0.1, 10.0, 0.05

这里每一行代表一个 Prony 项：第一列是模量比例 g_i，第二列是松弛时间 τ_i（秒），第三列是体积模量比例 k_i。

我的经验： 塑料的粘弹性参数很难从材料手册上直接查到。我通常的做法是：做一组不同速率下的压缩/拉伸实验，然后反推 Prony 参数。别偷懒，这一步省不了。

还有一个常见问题：塑料的屈服。很多塑料在插拔过程中不会明显屈服，但会有永久变形。这时候可以用Drucker-Prager 模型或者Crushable Foam 模型。我个人更推荐前者，因为它能考虑静水压力对屈服的影响。

避坑指南： 我曾经用线弹性模型算一个 LCP 外壳，结果插拔力偏小 30%。后来发现，塑料在接触区域发生了局部屈服，导致刚度下降。改用弹塑性模型后，结果就对上了。

摩擦系数，这是插拔力仿真里最「玄学」的参数。你想想看，端子跟端子接触，表面有镀层、有润滑剂、有氧化物，摩擦系数怎么可能是一个固定值？

我建议用库仑摩擦模型，但摩擦系数要区分静摩擦和动摩擦。在 Abaqus 里：

*Surface Interaction, name=Contact_Prop
*Friction
 0.25, 0.20

第一列是静摩擦系数，第二列是动摩擦系数。注意，Abaqus 默认用指数衰减来过渡，你还可以指定衰减系数。

摩擦系数标定方法：

这里有个细节：摩擦系数跟正压力有关。你想想看，正压力大了，接触面可能发生塑性变形，摩擦系数反而会变化。所以，我建议在仿真里用压力依赖的摩擦模型。

*Friction, pressure dependency
 0.25, 0.0
 0.20, 100.0
 0.18, 200.0

每一行：第一列是摩擦系数，第二列是对应的接触压力（MPa）。这样，Abaqus 会根据实际接触压力自动插值。

我的经验： 镀金端子的摩擦系数一般在 0.15~0.25 之间，镀锡的在 0.3~0.5 之间。如果加了润滑剂，可以降到 0.1 以下。但注意，润滑剂会随着插拔次数增加而消耗，摩擦系数会逐渐增大。

为了让大家更直观地理解这一章的内容，我画了一张图：

这一章的内容，说白了就是三件事：

嗯，这里要注意：材料参数不是越复杂越好。我见过有人用 Chaboche 模型做端子，参数标定了两个月，结果精度只提高了 3%。没必要。对于插拔力仿真，双线性模型 + 合理的摩擦系数，已经能解决 90% 的问题。

最后，送大家一句话：仿真准不准，材料参数是关键。别在网格和接触设置上花太多时间，先把材料参数搞准了。