第二章：有限元基础理论

各位同学，欢迎来到第二章。说实话，这一章是整门课里最“硬”的部分。弹性力学、接触力学、非线性求解器——听起来就头大对吧？

我当年刚入行时也觉得这些理论离实际太远。直到有一次，我算出来的插拔力跟实测差了30%……嗯，从那以后我再也不敢跳过理论基础了。

2.1 弹性力学基础：别被张量吓到

弹性力学说白了就是研究“东西怎么变形、怎么受力”的学问。你想想看，连接器端子插进去，弹片被撑开，这不就是典型的弹性变形吗？

我个人习惯把弹性力学拆成三块来理解：

• 应力：单位面积上的力。说白了就是“多挤”。
• 应变：变形程度。说白了就是“多歪”。
• 本构关系：应力和应变之间的“交易规则”。

对于连接器常用的铜合金，我们一般用线弹性本构。公式很简单：σ = E·ε。E就是弹性模量，铜合金通常在110~130 GPa之间。

我在项目中遇到过一位同事，他把不锈钢的弹性模量输成了200 MPa而不是200 GPa……结果算出来的插拔力只有0.01N。这种低级错误，大家千万别犯。

2.2 接触力学基础：这才是核心

连接器插拔力仿真，核心就是接触。没有接触，就没有插拔力。

接触力学研究的是两个物体“碰在一起”时发生了什么。你想想看，端子插入时，公端和母端的接触面在干嘛？

• 法向：互相挤压，产生正压力
• 切向：相对滑动，产生摩擦力

有限元里处理接触，最常用的方法是罚函数法。它的逻辑很简单：

接触力 = 罚刚度 × 穿透量

罚刚度越大，穿透越小，但计算越容易发散。我建议新手把罚刚度设为材料弹性模量的0.1~1倍，先试算一下。

还有一个重要概念：摩擦模型。连接器里常用库仑摩擦，公式是：

摩擦力 = μ × 正压力

μ是摩擦系数。镀金端子一般在0.1~0.3之间，镀锡端子可能到0.5以上。这个值对插拔力影响很大，我建议你查一下材料手册，别瞎猜。

2.3 非线性求解器原理：为什么不是一步算完？

你可能会问：为什么不能像静力学分析那样，一步算出结果？

因为连接器插拔过程有三个非线性来源：

1. 几何非线性：弹片大变形，形状变了，刚度也变了
2. 材料非线性：可能进入塑性
3. 接触非线性：接触状态不断变化（接触/分离/滑动）

这三个“非”加在一起，求解器只能一步步来。常用的方法是牛顿-拉普森法。

具体流程是这样的：

1. 施加一小步位移（比如0.01mm）
2. 计算当前刚度矩阵
3. 求解位移增量
4. 检查内外力是否平衡
5. 如果不平衡，修正位移，重复2-4
6. 平衡了，进入下一步

我刚开始做仿真时，总想一步算完，把步长设得很大。结果求解器反复迭代就是不收敛，最后报错退出。后来我学乖了，把初始步长设成总位移的1%，稳得很。

2.4 知识体系框架

下面这张图是我自己整理的，把本章三个核心内容串起来了。你保存好，以后做仿真时对照着看。

2.5 新手常见问题

我带过不少新人，发现几个高频问题，这里统一说一下：

• 问题一：“我材料参数输对了，为什么结果不对？”——检查单位！我见过有人把mm和m混用，结果差了1000倍。
• 问题二：“接触就是不收敛怎么办？”——先减小步长，再检查接触方向有没有设反。
• 问题三：“算出来的插拔力曲线有锯齿”——正常！那是接触状态切换导致的，把输出步长加密就能平滑。

好了，这一章的内容就到这里。弹性力学、接触力学、非线性求解器——这三个东西你搞懂了，后面的仿真操作就是水到渠成的事。