弹簧单元与连接器
理论与物理
什么是弹簧单元
老师，FEM中的“弹簧单元”在什么时候使用呢？
弹簧单元是在两点之间进行弹性连接的最简单单元。不仅用于物理弹簧（如螺旋弹簧等）的建模，还广泛用作连接部位或支撑条件的简化模型。
具体有哪些场景呢？
弹簧单元的刚度矩阵
线性弹簧的刚度矩阵是2×2：
其中 $k$ 是弹簧常数。这是FEM所有单元中最简单的刚度矩阵。
和梁单元的桁架部分（$EA/L$）形式一样呢。
桁架单元等价于 $k = EA/L$ 的弹簧单元。弹簧单元是其一般化形式，可以设定任意的 $k$。
弹簧单元的种类
| 种类 | 自由度 | 用途 |
|---|---|---|
| 平移弹簧（SPRING1/2） | 平移1方向 | 地基弹簧、轴向连接 |
| 旋转弹簧 | 旋转1方向 | 半刚性连接 |
| 6自由度弹簧（BUSHING） | 全部6自由度 | 衬套、弹性支座 |
| 接地弹簧（grounded spring） | 仅1个节点 | 地基支撑、弹性支撑 |
接地弹簧是一端固定在“地面”上的弹簧吗？
是的。也称为单节点弹簧。典型应用是用“地基反力系数 × 面积”的弹簧来表示地基上的基础。在Nastran中是CELAS1/CELAS2，在Abaqus中是SPRING1。
非线性弹簧
也有非线性的弹簧吗？
有的。力-位移关系可以用表格（折线）定义。
- 双线性弹簧 — 具有屈服点的弹塑性弹簧
- 间隙弹簧 — 仅在超过一定间隙时才产生力
- 非线性弹性 — 任意的F-δ曲线
Abaqus的*CONNECTOR ELEMENT是可以自由定义非线性力-位移、力矩-旋转、摩擦、阻尼的通用连接器。
各求解器对应的单元名称
| 种类 | Nastran | Abaqus | Ansys |
|---|---|---|---|
| 标量弹簧 | CELAS1/2 | SPRING1/2 | COMBIN14 |
| 衬套 | CBUSH | *CONNECTOR | COMBIN40 |
| 非线性弹簧 | CBUSH1D(NL) | *CONNECTOR(NL) | COMBIN39 |
Abaqus的CONNECTOR ELEMENT是最通用的吗？
是的。Abaqus的*CONNECTOR 可以在一个单元中定义“弹簧”、“阻尼器”、“摩擦”、“锁定”、“限位器”。Nastran的CBUSH也是多自由度的，但在非线性方面Abaqus更灵活。
总结
我来整理一下弹簧单元的理论。
要点：
- 在两点之间进行弹性连接的最简单单元 — 刚度矩阵为2×2
- 广泛用于连接部位・支撑条件的简化 — 螺栓、地基弹簧、半刚性连接
- 非线性弹簧 — 力-位移表格、间隙、摩擦
- CONNECTOR（Abaqus）是最通用的 — 整合了多自由度、非线性、摩擦
- 弹簧常数的合理性支配着结果 — 必须设定有物理依据的 $k$
最后一点很重要呢。如果随意设定弹簧常数，结果也会变得随意。
弹簧单元设置简单，但弹簧常数的物理依据就是一切。地基反力系数、螺栓的轴向刚度、连接部位的旋转刚度…能否正确计算这些，是工程师能力的体现。
弹簧单元刚度矩阵的推导
一维弹簧单元的刚度矩阵是k[1,-1;-1,1]（k：弹簧常数），是最简单的有限元。这个2×2矩阵直接从胡克定律F=kδ推导而来，出现在1956年Turner・Clough・Martin・Topp发表的有限元法原型论文《复杂结构的刚度和挠度分析》中。成为了所有结构FEM单元教科书的起点。
各项的物理意义
- 惯性项（质量项）：$\rho \ddot{u}$，即“质量×加速度”。您有过急刹车时身体被向前甩出去的经历吗？那种“被带走的感觉”正是惯性力。物体越重越难启动，一旦启动也越难停止。地震时建筑物摇晃，也是因为地面突然移动，而建筑物的质量“被落下”。静力分析中此项设为零，这是“因为缓慢施加力所以加速度可以忽略”的假设。在冲击载荷或振动问题中绝对不能省略。
- 刚度项（弹性恢复力）：$Ku$ 或 $\nabla \cdot \sigma$。拉弹簧时会感觉到“想要恢复的力”吧？那就是胡克定律 $F=kx$，也是刚度项的本质。那么提问——铁棒和橡皮筋，用相同的力拉，哪个会伸长？当然是橡皮筋。这种“不易伸长性”就是杨氏模量 $E$，它决定了刚度。常见的误解：“刚度高＝强度高”是不对的。刚度是“不易变形”，强度是“不易破坏”，是不同的概念。
- 外力项（载荷项）：体积力 $f_b$（重力等）和表面力 $f_s$（压力、接触力等）。可以这样想——桥上卡车的重量是“作用在整个内部上的力”（体积力），轮胎压路面的力是“只作用在表面上的力”（表面力）。风压、水压、螺栓的紧固力…全都是外力。这里容易犯的错误：弄错载荷的方向。本想“拉伸”却变成了“压缩”——听起来像笑话，但在3D空间中坐标系发生旋转时确实会发生。
- 阻尼项：瑞利阻尼 $C\dot{u} = (\alpha M + \beta K)\dot{u}$。试着弹一下吉他的弦。声音会一直持续吗？不，会逐渐变小。因为振动能量通过空气阻力或弦的内部摩擦变成了热。汽车的减震器也是同样的原理——特意吸收振动能量来改善乘坐舒适性。如果阻尼为零会怎样？建筑物在地震后会一直摇晃下去。实际上不会这样，所以设定适当的阻尼很重要。
假设条件与适用范围
量纲分析与单位制
| 变量 | SI单位 | 注意事项・换算备忘 |
|---|---|---|
| 位移 $u$ | m（米） | 输入mm时，载荷・弹性模量也要统一为MPa/N系 |
| 应力 $\sigma$ | Pa（帕斯卡）= N/m² | MPa = 10⁶ Pa。与屈服应力比较时注意单位制不一致 |
| 应变 $\varepsilon$ | 无量纲（m/m） | 注意工程应变与对数应变的区别（大变形时） |
| 弹性模量 $E$ | Pa | 钢: 约210 GPa，铝: 约70 GPa。注意温度依赖性 |
| 密度 $\rho$ | kg/m³ | mm系中是tonne/mm³（钢为 = 10⁻⁹ tonne/mm³） |
| 力 $F$ | N（牛顿） | mm系用N，m系也用N统一 |
数值解法与实现
弹簧单元的实现细节
实现弹簧单元时有什么需要注意的吗？
看似简单，却有意想不到的陷阱。
坐标系的问题
弹簧的方向如何定义？
弹簧单元的刚度作用于全局坐标系的特定方向。如果想在倾斜方向设置弹簧，需要定义局部坐标系或指定方向向量。
注意事项：
- Nastran的CELAS1以网格为单位指定方向（自由度编号）
- Abaqus的SPRING单元默认使用连接节点间的方向，但也可指定任意方向
- Ansys的COMBIN14默认是全局轴方向。可通过KEYOPT更改为局部方向
如果方向弄错，弹簧就会在意想不到的方向上起作用呢。
这是最常见的错误。添加了弹簧单元但结果几乎没变时，多半是弹簧方向错了。
接地弹簧的建模
地基弹簧（接地弹簧）如何建模？
以桩基础的地基弹簧为例。用弹簧表示深度方向的地基反力：
其中 $k_s$ 是地基反力系数（kN/m³），$D$ 是桩径，$\Delta z$ 是单元长度。
地基反力系数如何确定？
根据地基勘察（钻孔、标准贯入试验等）决定。
| 地基 | $k_s$ 的参考值 (kN/m³) |
|---|---|
| 软弱粘土 | 2,000 〜 5,000 |
| 中等粘土 | 10,000 〜 30,000 |
| 硬粘土 | 30,000 〜 100,000 |
| 砂土（松散） | 5,000 〜 15,000 |
| 砂土（密实） | 30,000 〜 100,000 |
弹簧常数的数量级会差2位数！ 地基的准确评估很重要呢。
地基弹簧常数的不确定性直接关系到结构的响应。必须进行灵敏度分析（上下改变 $k_s$ 观察响应的变化）。
CONNECTOR单元（Abaqus）
请告诉我Abaqus的CONNECTOR单元的设置方法。
この記事の評価ご回答ありがとうございます！参考に
なったもっと
詳しく誤りを
報告
请告诉我Abaqus的CONNECTOR单元的设置方法。
なった
詳しく
報告