4. 接触对创建流程
在 ANSYS 接触非线性分析中，接触对（Contact Pair）是承载全部接触行为的核心对象。一个接触对由两个面集合构成：接触面（Contact Surface）与目标面（Target Surface）。接触面通常由 CONTA17x 系列单元（如 CONTA174）离散，目标面由 TARGE17x 系列单元（如 TARGE170）离散。二者通过共享同一个实常数编号（Real Constant Set ID）建立配对关系，因此实常数既是几何配对标识，也是接触算法的参数载体。
本节围绕“怎么建、怎么选、怎么设”三个问题展开：先明确目标面与接触面的分配原则，再走一遍接触对向导的完整操作路径，最后落到实常数与关键选项的具体取值。
一、目标面与接触面的分配原则
接触对是“非对称”的：求解器只检查接触面节点是否穿透目标面，而不反向检查。因此分配错误会直接导致穿透、收敛困难或结果失真。工程上遵循以下判据：
- 刚度与网格：凸面、网格较粗的一侧作目标面；凹面、网格较密的一侧作接触面。若两侧网格密度接近，则把刚度较大（材料模量高）的一侧作目标面。
- 几何阶次：高阶单元面（如 186/187 的二次面）宜作接触面，低阶单元面作目标面，避免高阶面被低阶面“切穿”。
- 运动关系：大滑动问题中，通常把运动剧烈、可能滑入对方的一侧设为接触面。
- 刚柔组合：刚体—柔体接触中，刚体面必须作目标面，柔体面作接触面。
以轴孔过盈配合为例：孔内表面为凹面且网格较密，设为接触面；轴外表面为凸面，设为目标面。若反过来设置，轴表面节点会穿透孔壁，过盈量无法正确传递。
二、接触对向导的操作路径
在 ANSYS Mechanical APDL 中，接触对向导（Contact Wizard）可通过 Main Menu → Modeling → Create → Contact Pair 启动，也可用命令流实现。向导的核心逻辑是“先选面、再定算法、后设参数”。
以轴孔过盈为例，假设孔内表面已选中为组件 HOLE_IN，轴外表面为 SHAFT_OUT，实常数编号取 3，接触算法选用增广拉格朗日法（Augmented Lagrange），命令流如下：
! 1. 定义接触面单元类型（CONTA174）
ET,2,CONTA174
KEYOPT,2,2,2 ! 接触算法：2=增广拉格朗日
KEYOPT,2,12,5 ! 接触面行为：5=绑定/不分离按需
! 2. 定义目标面单元类型（TARGE170）
ET,3,TARGE170
KEYOPT,3,2,0 ! 目标面约束方式
! 3. 定义实常数（ID=3）
R,3
RMODIF,3,3,0.3 ! FKN=0.3，过盈问题降低法向刚度
RMODIF,3,4,0.05 ! FTOLN=0.05
RMODIF,3,5,0.5 ! FKT=0.5
RMODIF,3,10,1 ! ICONT=1，允许初始接触调整
! 4. 在接触面上生成 CONTA174
TYPE,2
REAL,3
ESEL,S,COMPONENT,HOLE_IN
NSLE,S
ESLN,S,1
ESURF
! 5. 在目标面上生成 TARGE170
TYPE,3
REAL,3
ESEL,S,COMPONENT,SHAFT_OUT
NSLE,S
ESLN,S,1
ESURF
ALLSEL关键点在于：接触面与目标面必须使用同一个 REAL 编号，否则求解器无法识别配对关系。生成后应检查接触单元法向是否指向目标面，可用 /PSYMB,ESYS,1 显示单元坐标系，法向箭头应指向对方。
三、实常数与关键选项设置
实常数是接触对参数的主要入口，常用项及其工程取值如下表所示（示例值针对轴孔过盈配合，属合理假设，实际需按模型试算调整）：
| 实常数项 | 含义 | 示例取值 | 调整依据 |
|---|---|---|---|
| FKN | 法向接触刚度因子 | 0.3 | 过盈量大、收敛困难时降低 |
| FTOLN | 允许穿透容差 | 0.05 | 过盈配合可适当放宽 |
| FKT | 切向（摩擦）刚度因子 | 0.5 | 摩擦主导时取 0.01~1.0 |
| ICONT | 初始接触调整 | 1 | 过盈配合用 1 或 2 |
| PINB | 球区半径 | 1.0 | 默认即可，特殊几何再调 |
| CNOF | 接触面偏移 | 0 | 需模拟过盈量时设为负值 |
除实常数外，KEYOPT 关键选项同样决定接触行为。最常用的几项：
- KEYOPT(2) 接触算法：0=罚函数法，1=罚函数+拉格朗日，2=增广拉格朗日，3=拉格朗日乘子法。过盈配合推荐 2，兼顾精度与收敛。
- KEYOPT(5) 初始接触闭合：0=不调整，1=自动闭合间隙，2=自动减少穿透，3=同时处理。过盈问题常用 1 或 3。
- KEYOPT(9) 初始穿透/间隙：0=忽略，1=按实常数 CNOF 处理，2=按几何处理。需精确控制过盈量时选 2。
- KEYOPT(12) 接触面行为：0=标准，1=不分离，2=绑定，3=不分离+绑定，5=绑定但允许滑动。过盈配合常用 0 或 1。
CNCHECK 命令检查接触对：CNCHECK,SUMMARY 列出所有接触对及实常数，CNCHECK,DETAIL 检查初始穿透与间隙。若出现大量初始穿透，优先调整 ICONT 或 CNOF，而不是盲目降低 FKN。/PSYMB,ESYS,1 目视检查法向。完成上述设置后，接触对即具备可求解条件。后续求解中若出现“接触状态频繁切换”或“穿透量超限”，应回到实常数与 KEYOPT 逐项排查，而不是直接修改载荷步。接触对的质量，本质上取决于面分配是否合理、实常数是否匹配物理场景这两件事。