ANSYS—接触单元说明.docx

参考ANSYS勺中文帮助文件

接触问题（参考ANSYS勺中文帮助文件）

当两个分离勺表面互相碰触并共切时，就称它们牌接触状态。在一般勺物理意义中，牌接触状态勺表面有下列特点：

1、不互相渗透；

2、能够互相传递法向压力和切向摩擦力；

3、通常不传递法向拉力。

接触分类：刚性体－柔性体、柔性体－柔性体

实际接触体相互不穿透，因此，程序必须在这两个面间建立一种关系，防止它们在有限元分析中相互穿过。

――罚函数法。接触刚度

―― lagrange 乘子法，增加一个附加自由度（接触压力），来满足不穿透条件―― 将罚函数法和lagrange 乘子法结合起来，称之为增广lagrange 法。三种接触单元：节点对节点、节点对面、面对面。接触单元勺实常数和单元选项设臵：

FKN法向接触刚度。这个值应该足够大，使接触穿透量小；同时也应该足够小，使问题没有病态矩阵。FKN值通常在0.1~10之间，对于体积变形问题，用值1.0 （默认），对弯曲问题，用值0.1 。

FTOLN最大穿透容差。穿透超过此值将尝试新的迭代。这是一个与接触单元下面的实体单元深度（h）相乘的比例系数，缺省为0.1。此值太小，会引起收敛困难。

ICoNT初始接触调整带。它能用于围绕目标面给出一个“调整带”，调整带内任何接触点都被移到目标面上；如果不给出ICONT fi, ANSYSB据模型的大小提供一个较小的默认值（V 0.03 =

PINB 指定近区域接触范围（球形区）。当目标单元进入pinball 区时，认为它处于近区域接触，pinball 区是围绕接触单元接触检测点的圆（二维）或球（三维）。可以用实常数PINB调整球形区（此方法用于初始穿透大的问题是必要的）PMIN和PMAX初始容许穿透容差。这两个参数指定初始穿透范围，ANSYS把整个目标面（连同变形体）移到到由PMIN和PMAX旨定的穿透范围内，而使其成为闭合接触的初始状态。初始调整是一个迭代过程，ANSYS最多使用20个迭代步

把目标面调整到PMIN和PMAXE围内，如果无法完成，给出警告，可能需要修改几何模型。TAUMAX接触面的最大等效剪应力。给出这个参数在于，不管接触压力值多大，只要等效剪应力达到最大值TAUMAX就会发生滑动。该剪应力极限值通常用于接触压力会变得非常大的情况。

CNOF指定接触面偏移。+CNO增加过盈、-CNOF减少过盈或产生间隙、CNO能与几何穿透组合应用。

FKOP接触张开弹簧刚度。针对不分离或绑定接触模型，需要设臵实常数FKOP 该常数为张开接触提供了一个刚度值。FKOPH止接触面的分离；FKOP默认为1.0，用于建立粘结模型，用一个较小值（1e-5）去建立软弹簧模型。

FKT切向接触刚度。作为初值，可以采用—FKT= 0.01*FKN,这是大多数ANSYS 接触单元的缺省值。

COHE 粘滞力。即没有法向压力时开始滑动的摩擦应力值。

FACT，DC 定义摩擦系数变化规律

MU= MUK*（1+（FACT-I） EXP（- DC*vtfs∕deltaT ））

式中：MU U动摩擦系数（用户自己定义）

FACT= MUS/MUK用户自己定义）

MU U静摩擦系数

DC=衰减系数（用户自己定义）

Vtfs/deltaT U表面间的相对速度

注意：动摩擦系数由被指定为材料属性（MU ,由MP命令或GUI定义。缺省值：

FACT U 1，MUS U MUK U 0，DC U 0

Keyopt 的介绍,以Target170,Conta173 为例：

首先介绍170 的Keyopt

KEYoP（3）:定义接触行为

KEYoP（I）:单元阶数（是否含有中节点）

KEYoP（T1）U 0：低阶单元（不含中节点）

KEYoP（T1）U 1：高阶单元（含中节点）

KEYOP T2）:刚体目标面约束条件

U 0 时,自动约束选项,每一个载荷步的末尾,程序内部将是性面重新设臵约束。

满足以下条件,刚性面则缺省为自动约束没有明确定义边界条件；目标面与其它单元没有联系；没有定义耦合或约束方程。

U 1 时,用户定义选项。

Conta173 的Keyopt：

KEYOPT1）:自由度选项。U 0时，结构：UX UY和UZ; U 1时，结构和热；U 2时，TEMP（用于纯热接触问题）

KEYOP i2）：选择接触算法。U 0时，增广的拉格朗日法（缺省选项），推荐于一般应用，它对罚刚度不太敏感，但是也要求给出一个穿透容差。U 1时，罚函数法。它推荐应用于单元非常扭曲、大摩擦系数和用增广的拉格朗日法收敛行为不好的问题。

KEYOP T4）:选择接触检查点。U 0时，高斯点（缺省选项，推荐）；=1 时，节点；ANSY面对面单元默认用高斯积分点作为接触检查点。

KEYOPT5）:自动CNO调整。允许ANSYSS于初始状态自动给定CNOF S―― 导致“刚好接触”配臵。U 0时，不进行自动调整；U 1时，闭合间隙；U 2时, 减小穿透；U 3时，闭合间隙/减小穿透。

KEYOPT7）:时间步控制选项。（只有在Solution COntrOl中打开基于接触状态变化的时间步预测，此选项才起作用。Solution>Unabridged Menu>Load Step Opts>Solution Ctrl

U 0时，不控制，不影响自动时间步长。对静力问题自动时间步打开时此选项一般是足够的。

U 1 时，自动二分，如果接触状态变化明显，时间步长将二分，对于动力问题自动二分通常是足够的；

U 2时，合理值。比自动细分更耗时的算法；

U 3时，最小值。此选项为下一子步预测最小时间增量（很耗机时，不推荐）。

KEYoP（8）:防止伪接触选项。=0时，不防止；=1时，检测并忽略伪接触。KEYoPT9）:初始穿透间隙控制。=0时，包括几何穿透/间隙和CNOF = 1时, 忽略几何穿透/间隙和CNOF = 2时，包括几何穿透/间隙和CNOF且在第一个载荷步中渐变；=3时，忽略几何穿透/间隙，包括CNOF = 4时，忽略几何穿透/间隙，包括CNO，F 且在第一个载荷步中渐变。KEYOP T 10）:接触刚度更新控制。=0,闭合状态的接触刚度不进行任何更新； =1,每一载荷步更新闭合状态的接触刚度（FKN或FKT,由用户指定）；=2,与 =1同，此外，在每一子步，程序自动更新接触刚度（根据变形后下伏单元的刚度）。

KEYOP l 11）:壳、梁单元厚度影响。