📚 课程目录第2章复杂装配体接触非线性分析

2. 复杂装配体接触非线性分析

复杂装配体与单一零件的最大区别，在于载荷要通过零件之间的接触面传递。接触是典型的边界非线性：接触状态（分离、黏着、滑移）会随载荷变化而改变，刚度矩阵也随之变化。Ansys Mechanical 中，接触设置是否合理，往往直接决定求解能否收敛、结果是否可信。本节围绕接触类型选择、接触刚度与穿透控制、摩擦模型、螺栓预紧力施加，以及收敛调试展开。

一、接触类型：绑定、不分离、无摩擦、摩擦与粗糙

Mechanical 提供五种接触类型，其法向与切向行为差异明显，选择依据是“装配面之间是否允许分离、是否允许滑移”。

Bonded（绑定）：法向不分离、切向不滑移，相当于把两个面焊在一起，是线性接触，收敛最快。适用于焊接、过盈配合的简化、螺栓与被连接件的近似连接。
No Separation（不分离）：法向不分离、切向可无摩擦滑移。适用于允许微小滑移但不允许脱开的配合面。
Frictionless（无摩擦）：法向可分离、切向摩擦系数为零。适用于润滑良好或摩擦影响可忽略的滑动面，但容易因刚体位移导致收敛困难。
Frictional（摩擦）：法向可分离、切向按库仑摩擦传递剪应力，是最贴近真实装配的模型，也是收敛最困难的类型。
Rough（粗糙）：法向可分离、切向不滑移（摩擦系数无穷大）。适用于不允许滑移的粗糙结合面。

需要强调：接触类型不是“越真实越好”，而是“与关注量匹配”。若只关心整体刚度与应力分布，绑定往往足够；若关心接触面滑移量、接触压力分布或密封性，则必须使用摩擦接触。

二、接触刚度与穿透控制

Mechanical 默认使用增广拉格朗日法（Augmented Lagrange），法向刚度由法向刚度因子 FKN 控制，默认 1.0。FKN 越大，穿透越小，但越难收敛；FKN 越小，越易收敛，但穿透可能大到影响结果。经验做法是：先取 0.1～0.5 快速获得收敛解，再逐步提高到 1.0～2.0 复核穿透量。

穿透量可通过“Contact Tool → Penetration”查看。工程上一般要求最大穿透不超过局部单元尺寸的 5%～10%，或不超过关注特征尺寸的 1%。若穿透过大，可提高 FKN、细化接触区网格，或改用 MPC 公式（适用于绑定与不分离）。

法向刚度更新方式也影响收敛：Each Iteration 每迭代更新，收敛稳健但耗时；Each Iteration, Aggressive 更激进；Never 不更新，仅适合线性接触。摩擦接触建议使用默认的每迭代更新。

三、摩擦模型

库仑摩擦的切向剪应力上限为 τ_max = μ·p + c，其中 μ 为摩擦系数，p 为接触压力，c 为黏聚力。Mechanical 中可设置摩擦系数、黏聚力，以及弹性滑移容差（Elastic Slip Tolerance）。弹性滑移容差默认 0.01，表示在达到滑移极限前允许的弹性切向变形比例；该值越小，切向刚度越大，越接近理想库仑摩擦，但收敛越难。

摩擦系数需按材料配对选取，例如钢-钢干燥面约 0.15～0.3，钢-混凝土约 0.4～0.6，润滑面约 0.05～0.1。这些是常见工程参考值，实际应以试验或手册为准。若摩擦系数过大且载荷路径依赖切向传递，可能出现“黏着-滑移”振荡，此时可适当降低弹性滑移容差或改用增广拉格朗日法。

四、螺栓预紧力施加

螺栓预紧是装配体分析的核心载荷之一。Mechanical 提供 Bolt Pretension 载荷，施加在螺栓圆柱面上，需先定义局部坐标系使轴向与螺栓轴线一致。预紧力 F_i 与拧紧扭矩 T 的近似关系为 T = K·F_i·d，其中 K 为扭矩系数（干燥钢约 0.2），d 为螺纹公称直径。例如 M16 螺栓、目标预紧力 60 kN，则 T ≈ 0.2×60000×0.016 = 192 N·m。

预紧力施加分两步：第一步在螺栓上施加预紧力并锁定（Lock），第二步施加外载荷。若直接一步加载，预紧力与外载耦合，结果失真。Mechanical 中通过“Bolt Pretension → Adjustments → Lock”实现，或在多载荷步中设置。

示例：M16 螺栓预紧
目标预紧力 F_i = 60 kN
扭矩系数 K = 0.2
拧紧扭矩 T = K·F_i·d = 0.2 × 60000 × 0.016 = 192 N·m
接触设置：螺栓头/螺母与被连接件 Bonded 或 Frictional(μ=0.15)
被连接件结合面：Frictional(μ=0.2)
法向刚度因子 FKN = 1.0，弹性滑移容差 = 0.01

五、收敛调试与工程案例

接触非线性收敛困难通常表现为：力残差震荡、接触状态反复切换、刚体位移。调试顺序建议为：先检查是否存在未约束的刚体位移；再降低 FKN 或增大弹性滑移容差；然后细化接触区网格；必要时启用“Stabilization（稳定化）”或“Contact Stabilization”阻尼；最后检查载荷步是否过大，改为多子步逐步加载。

以一个法兰-螺栓装配体为例：法兰结合面用 Frictional(μ=0.2)，螺栓用 Bolt Pretension 60 kN 并锁定，外载为内压 2 MPa。初始 FKN=1.0 时求解在第二子步发散，残差曲线在 1e-2 附近震荡。将 FKN 降至 0.3、弹性滑移容差改为 0.02 后收敛，随后用 Contact Tool 检查穿透，最大穿透 0.012 mm，小于局部单元尺寸 0.5 mm 的 5%，满足要求。最终接触压力在螺栓附近呈环状分布，法兰边缘出现分离，符合预期。

调试口诀：先约束、再降刚、后细化、加稳定、分步加载。每次只改一个参数，便于定位原因。

综上，复杂装配体接触分析的关键在于：接触类型与关注量匹配、接触刚度与穿透量平衡、摩擦参数有据可依、螺栓预紧分步施加、收敛调试有序推进。掌握这五点，才能在保证收敛的同时获得可信的接触压力、分离区与滑移量结果。