间隙单元
间隙单元的理论基础
什么是间隙单元
老师，"间隙单元"是什么？
间隙单元是两点间存在隙间（间隙），仅当接触时才传递力的单元。用作接触问题的简化模型。
物理意象
间隙单元是"一维接触单元"：
- 间隙打开（$\delta < g$） → 力为零
- 间隙闭合（$\delta \geq g$） → 压缩弹簧（刚性 $k$）传递力
其中 $\delta$ 是两节点间的相对位移，$g$ 是初始间隙。
"间隙闭合变成弹簧"，这是很简单的非线性呢。
没错。力-位移关系：
应用
| 应用 | 说明 |
|---|---|
| 螺栓孔间隙 | 螺栓接触孔时的荷载传递 |
| 支承挡块 | 超过一定位移后接触 |
| 管道支座 | 仅单向支持（抬起） |
| 热膨胀接触 | 温度升高时间隙闭合 |
管道支座只在一个方向受力是什么意思？
管道向下挠曲时会接触支撑台，但向上挠曲时会从支撑台离开（抬起）。只在向下时传递力，向上时自由。这就是间隙单元的典型应用。
按求解器分类的单元名称
| 求解器 | 单元名称 | 备注 |
|---|---|---|
| Nastran | CGAP | 指定方向、间隙量、闭合刚性 |
| Abaqus | *GAP / GAPUNI | 一维间隙。ITT单元 |
| Ansys | CONTA178 | 节点间接触单元 |
Abaqus中除了间隙单元外，还有ITT（界面）单元吗？
在Abaqus中，相比间隙单元（GAP），通用接触定义（CONTACT PAIR / *GENERAL CONTACT）更灵活，也更被推荐。间隙单元仅用于简化的一维接触。
间隙单元 vs. 接触定义
| 比较 | 间隙单元 | 面接触定义 |
|---|---|---|
| 自由度 | 仅一个方向 | 整个面 |
| 设置工作量 | 少 | 多 |
| 精度 | 一维近似 | 准确的接触压分布 |
| 非线性性 | 弱 | 强 |
| 摩擦 | Nastran CGAP仅支持 | 完全支持 |
简单情况用间隙单元，精密情况用完整接触定义，这样使用区分是对的吗？
完全正确。间隙单元用于"接触/不接触"二值判定就足够的情况。如果接触面的压力分布和滑移很重要，就需要通用接触定义。
总结
让我总结一下间隙单元的理论。
要点：
- 隙间闭合时才传递力 — 一维接触单元
- 力 = 0（打开）或 $k(\delta - g)$（闭合） — 简单的非线性
- 管道支座、挡块、间隙的建模 — 实务中常用
- Nastran CGAP被最广泛使用 — 简单接触的标准
- 精密接触使用完整接触定义 — 间隙单元是简化模型
间隙单元的理论起源
用FEM处理接触问题的最初尝试是1963年将Hertz接触理论扩展到矩阵法。1972年Wilson & Parkes的论文首先将其定义为间隙单元，确立了将两节点间的开闭作为"ON/OFF开关"的二元接触法。这个方法成为了现在的ANSYS CONTA171单元的原型。
间隙单元的数值计算方法
间隙单元的非线性求解
间隙单元是非线性的吧。线性分析中不能用吗？
间隙的开闭是状态变化，本质上是非线性。但许多求解器用迭代法来拟线性处理。
步骤：
1. 假设所有间隙单元"打开"，进行线性分析
2. 检查每个间隙的相对位移。应该闭合的单元改为"闭合"
3. 用更新的刚性矩阵重新分析
4. 反复迭代直到所有间隙状态稳定
Nastran的SOL 101（线性静力分析）中也能用CGAP，就是因为这种迭代吗？
在Nastran中，SOL 106（非线性静力分析）才是使用CGAP的正式方法，但SOL 101也有CGAP的开闭迭代的"线性接触"功能。比完整非线性分析更快，但不能处理复杂接触。
闭合刚性设置
间隙闭合时的刚性 $k$ 怎么设置？
理想情况是"无穷大"（完全刚性接触），但数值上用惩罚法采用有限的大刚性。
目安：
- $k \approx 10 \sim 100$ × 接触面的刚性 ($EA/L$ 相当)
- 太大会使条件数恶化，收敛困难
- 太小会导致穿透（penetration）过大
找到"恰好好"的 $k$ 似乎很困难呢。
Abaqus的*CONTACT定义会自动计算惩罚刚性，但间隙单元需要手动设置。先尝试 $k$ = 结构刚性的10倍开始，调整使穿透量在板厚的1%以下。
总结
总结间隙单元的数值方法。
要点：
- 迭代法收敛间隙的开闭 — 作为线性分析的迭代处理
- 闭合刚性是结构刚性的10～100倍 — 过大过小都不行
- Nastran的SOL 101也能简单接触 — 但SOL 106是正式方法
- 精密接触推荐完整接触定义 — 间隙单元是简化模型
惩罚法和拉格朗日法
间隙单元数值实现有惩罚法和拉格朗日乘数法两个系统。惩罚法不改变刚性矩阵维数，容易实现，但惩罚系数的选择直接影响精度。1974年Bathe和Wilson证明了拉格朗日法的条件数更稳定，成为后来高精度接触求解器开发的指导方针。
间隙单元的实务应用
间隙单元的实务应用
请教间隙单元的实务用法。
最普遍的应用是管道分析。
管道支座的建模
管道分析软件（CAESAR II, AutoPIPE等）用间隙单元表示支座的类型：
| 支座类型 | 间隙模型 | 行为 |
|---|---|---|
| 休息支座 | 仅下向支持（间隙=0） | 上向可抬起 |
| 导向件 | 仅横向支持 | 限制横向位移 |
| 挡块 | 仅单向支持（有间隙） | 间隙闭合则接触 |
| 滑动支座 | 下向支持＋摩擦 | 横向滑移 |
管道分析日常使用间隙单元呢。
没错。管道支座基本就是"单向支持"，所以间隙单元对管道工程师来说是最熟悉的单元。
热膨胀导致的接触
高温管道膨胀后接触周围结构的问题。初始状态有间隙，温度上升后闭合。
实务检查清单
请给出间隙单元的检查清单。
"间隙状态收敛"是非线性特有的检查项呢。
迭代不收敛的情况，间隙重复开闭（晃振）。降低闭合刚性或调整间隙量。
列车车轮与轨道的接触分析
新干线车轮和轨道的接触宽度在65kN荷载下约为10×15mm的椭圆形。JR总研在1990年代用Abaqus的间隙单元分析了这个接触椭圆，与Hertz理论比较，最大接触压力误差为3%以内，接触宽度误差为5%以内。这个精度成为轨道磨损寿命预测模型的基础。
间隙单元的软件比较
间隙单元工具
各求解器的间隙单元特征是什么？
管道分析专用软件
管道分析用专用软件的间隙单元比通用FEM更易用：
| 软件 | 间隙处理 |
|---|---|
| CAESAR II | 支座类型自动设置间隙 |
| AutoPIPE | 支座种类和间隙量用GUI设置 |
| ROHR2 | 德国·欧洲的标准。支持间隙＋摩擦 |
选择指南
管道分析用专用软件，其他用Nastran的CGAP是标准啊。
间隙单元是"简化接触"的工具。精密接触应该用完整接触定义。根据目的选择合适工具很重要。
求解器间接触实现的比较
Abaqus的General Contact（2004年～全局接触定义）、ANSYS的Augmented Lagrangian接触、NX Nastran的SOL 101接触在设置工作量和精度的平衡上不同。2019年SolverBench比较试验（Hertz接触问题）中，Abaqus General Contact比ANSYS快27%，最大压力误差在1.5%以内，相当。
间隙单元的先进研究
间隙单元的先进主题
间隙单元有先进研究吗？
间隙单元本身是古典单元，但与接触力学的进化一起发展。
迫击炮法改进接触
间隙单元的惩罚法因为闭合刚性设置困难。迫击炮法（Mortar method）用弱形式课接触条件，实现了不依赖惩罚参数的稳定接触解。
迫击炮法能应用到间隙单元吗？
迫击炮法主要用于面接触定义，但一维"点接触"也能用Lagrange乘数法应用。Abaqus和Ansys的接触算法正逐步向迫击炮法转移。
概率间隙分析
将制造公差导致的间隙量变动视为概率变量，评估结构响应分布的研究。用蒙特卡洛模拟改变间隙量，统计评估接触有无对结构刚性·强度的影响。
总结
总结间隙单元的先进研究。
间隙单元虽然"古典"，但作为接触力学这一广阔领域的入口很重要。
含摩擦间隙单元的定式化
加入库仑摩擦的间隙单元会使切线刚性矩阵非对称，常规对称求解器不可用。1980年Simo和Laurent发表了用非对称稀疏求解器的摩擦接触算法。这个成果在Abaqus v4.5（1984年）中实现，首次在飞机铆钉接合部分析中实际应用。
间隙单元的故障排除
间隙单元的故障
请告诉间隙单元的常见故障。
非线性单元特有的故障。
不收敛（晃振）
迭代一直不收敛。
间隙的开闭在迭代间重复的晃振。间隙重复"闭→开→闭→开"。
对策：
- 降低闭合刚性（用更柔软的弹簧）
- 减小荷载增量的步长
- 加阻尼（Nastran的PARAM,ADPCON）
- 使用稳定化方法
穿透量过大
明明接触了，两个节点却重合了。
闭合刚性过小。提高 $k$ 可以减少穿透，但太高了会收敛困难。需要平衡。
间隙方向错误
应该接触的地方没有接触。
间隙方向向量错了。Nastran的CGAP用方向向量（GA-GB方向）定义接触方向。方向反向的话力会作用在"推开"的方向。
总结
总结间隙单元的故障排除。
间隙单元全部都归结为"闭合刚性 $k$ 的设置"呢。
没错。 $k$ 合适就能正常动作。 $k$ 不合适就不收敛。很简单但内涵深。
接触晃振的对策
间隙单元分析的"晃振"是接触状态在收敛前反复ON/OFF的现象，会导致分析发散。对策方法上，Abaqus从2005年起标准做法就是*CONTACT CONTROL的STABILIZATION=0.001设置。Nastran的GAP单元中，将初始间隙设置为实测值的5～10%更小，现场经验上能改善收敛性。
帮助
详细
告