Ansys 非线性接触分析和设置5.4.9 设置实常数和单元关键选项程序使用 20 个实常数和数个单元关键选项，来控制面─面接触单元的接触参见《ANSYS Elements Reference》中对接触单元的描述5.4.9.1 实常数在 20 个实常数中，两个(R1 和 R2)用来定义目标面单元的几何形状剩下的用来控制接触面单元R1 和 R2 定义目标单元几何形状FKN 定义法向接触刚度因子FTOLN 是基于单元厚度的一个系数，用于计算允许的穿透ICONT 定义初始闭合因子PINB 定义“Pinball"区域PMIN 和 PMAX 定义初始穿透的容许范围TAUMAR 指定最大的接触摩擦CNOF 指定施加于接触面的正或负的偏移值FKOP 指定在接触分开时施加的刚度系数FKT 指定切向接触刚度COHE 制定滑动抗力粘聚力TCC 指定热接触传导系数FHTG 指定摩擦耗散能量的热转换率SBCT 指定 Stefan-Boltzman 常数。
RDVF 指定辐射观察系数FWGT 指定在接触面和目标面之间热分布的权重系数FACT 静摩擦系数和动摩擦系数的比率DC 静、动摩擦衰减系数命令： R GUI：main menu> preprocessor>real constant对实常数 FKN, FTOLN, ICONT, PINB, PMAX, PMIN, FKOP 和 FKT，用户既可以定义一个正值，也可以定义一个负值程序将正值作为比例因子，将负值作为绝对值程序将下伏单元的厚度作为 ICON，FTOLN，PINB，PMAX 和 PMIN 的参考值例如 ICON = 0.1 表明初始闭合因子是“0.1*下层单元的厚度”然而，ICON = -0.1 则表示真实调整带是 0.1 单位如果下伏单元是超单元，则将接触单元的最小长度作为厚度参见 图 5-8 图 5-8 下层单元的厚度在模型中，如果单元尺寸变化很大，而且在实常数如 ICONT, FTOLN, PINB, PMAX, PMIN 中应用比例系数，则可能会出现问题因为从比例系数得到的实际结果，取决于下层单元的厚度，这就可能引起大、小单元之间的重大变化。
如果出现这一问题，请用绝对值代替比例系数TCC, FHTG, SBCT, RDVF 和 FWGT 仅用于热接触分析[KEYOPT(1)=1]5.4.9.2 单元关键选项每种接触单元都包括数个关键选项对大多的接触问题，缺省的关键选项是合适的而在某些情况下，可能需要改变缺省值下面是可以控制接触行为的一些关键选项：自由度 KEYOPT(1)接触算法(罚函数+拉格朗日乘子或罚函数) KEYOPT(2)存在超单元时的应力状态(仅 2D) KEYOPT(3)接触检测点的位置(仅低阶接触单元) KEYOPT(4)CNOF 自动调整 KEYOPT(5)时间步控制 KEYOPT(7)伪接触预防 KEYOPT(8)初始穿透或间隙的影响 KEYOPT(9)法向和切向接触刚度修正方法控制 KEYOPT(10)壳的厚度影响 KEYOPT(11)接触面行为(粗糙、绑定等) KEYOPT(12)命令： KEYOPT ET GUI：main menu>preprocessor>Elemant Type>Add/Edit/Delete5.4.9.3 选择接触算法对面─面接触单元，程序可以使用增进的拉格朗日方法或罚函数方法。
通过单元关键字 KEYOPT(2)来指定增进的拉格朗日方法是为了找到精确的拉格朗日乘子（即接触力），而对罚函数进行一系列修正迭代与罚函数的方法相比，拉格朗日方法容易得到良态条件，对接触刚度的敏感性较小然而，在有些分析中，增进的拉格朗日方法可能需要更多的迭代，特别是在变形后网格变得太扭曲时使用拉格朗日方法的同时应使用实常数 FTOLNFTOLN 为拉格朗日方法指定容许的最大穿透如果程序发现穿透大于此值时，即使不平衡力和位移增量已经满足了收敛准则，总的求解仍被当作不收敛处理FTLON 的缺省值为0.1用户可以改变这个值，但要注意，如果此值太小，可能会造成太多的迭代次数或者不收敛5.4.9.4 确定接触刚度所有的接触问题都需要定义接触刚度，两个表面之间穿透量的大小取决于接触刚度过大的接触刚度可能会引起总刚矩阵的病态，从而造成收敛困难一般来说，应该选取足够大的接触刚度以保证接触穿透小到可以接受，但同时又应该让接触刚度足够小以不致引起总刚矩阵的病态而保证收敛性ANSYS 程序根据下伏柔体单元的材料特性，来估计一个缺省的接触刚度值用户可用实常数 FKN 来为接触刚度指定一个比例因子或指定一个绝对值。
比例因子一般在 0.01 和 10 之间；对于大变形问题，选 1 是比较好的；而对于弯曲为主的问题，通常为 0.01～0.1用户应当总是检验以使穿透到达极小值，而又避免过多的迭代次数注意 --FTOLN 和 FKN 从一个荷载步到另一个荷载步中，都可以修改也可以在重启动中修改这时，必须定义 KEYOPT(10)=1,2为了确定一个较好的接触刚度值，可能需要一些经验用户可以按下面的步骤来进行尝试：1、 开始时取一个较低的值低估值要比高估值好，因为由一个较低的接触刚度导致的穿透问题，比过高的接触刚度导致的收敛性困难，要容易解决2、 对前几个子步进行计算分析，直到最终荷载的一个比例(刚好完全建立接触)3、 检查每一子步中的穿透量和平衡迭代次数如果总体收敛困难是由过大的穿透引起的(而不是由不平衡力和位移增量引起的)，那么可能低估了 FKN 的值，或者是将 FTOLN 的值取得大小如果总体的收敛困难是由于不平衡力和位移增量达到收敛值时需要过多的迭代次数，而不是由于过大的穿透量引起的，那么 FKN 的值可能被高估4、 按需要调整 FKN 或 FTOLN 的值，重新进行完整的分析注意 -- 如果穿透控制变成总体平衡迭代中的主因(如果为使问题收敛到穿透容差内，比收敛到不平衡力的容差内，需要更多的迭代)，用户应该增大 FTOLN 值，以允许更多的穿透，或增大 FKN。
5.4.9.5 选择摩擦类型在基本的库仑摩擦模型中，两个接触面在开始相互滑动之前，在它们的界面上会有达到某一大小的剪应力产生这种状态称为粘合状态(stick)库仑摩擦模型定义了一个等效剪应力 τ，在某一法向压应力 p 作用下剪应力达到此值时表面开始滑动 (τ=μp+COHE，其中 μ 是摩擦系数--MU--作为材料特性定义，而 COHE 是粘聚力)一旦剪应力超过此值后，两个表面之间将开始相互滑动这种状态，叫作滑动状态(Sliding)粘合/滑动计算决定什么时候一个点从粘合状态到滑动状态，或从滑动状态变到粘合状态摩擦系数可以是任一非负值，程序缺省值为表面之间无摩擦对于粗糙或绑定接触( KEYOPT(12)=1、3、5、6)，程序将不管给定的 MU 值而认为摩擦阻力无限大程序提供了一个人为指定最大等效剪应力的选项，不管接触压力值的大小，如果等效剪应力达到此值时，即发生滑动，见 图 5-9 为了指定接触界面上最大容许剪应力，设置常数 TAUMAX (缺省为 1.0E20)这个剪应力极限，通常用于在接触压力非常大的时候(如在某些加工过程中)的一些情况，以至于用库仑理论计算出的界面剪应力超过了材料的屈服极限。
TAUMAX 的一个合理上限估值为 是表面附近材料的 von Mises 屈服应力)经验数据有助于决定 TAUMAX 的值图 5-9 摩擦模式5.4.9.5.1 静、动摩擦系数摩擦系数依赖于接触面的相对滑动速度，通常静摩擦系数高于动摩擦系数ANSYS 提供了如下表示的指数衰减摩擦模型：μ=MU×(1+(FACT-1)exp(-DC×V rel )其中：μ 为摩擦系数MU 动摩擦系数，用 MP 命令输入FACT 是静摩擦系数与动摩擦系数之比，缺省为最小值 1.0DC 为衰减系数，缺省为 0.0，单位为 time/length因此，时间在静态分析中有一些意义V rel 是 ANSYS 计算的滑动速度如果知道静、动摩擦系数和至少一个数据点（μ 1 ，V rel ），则可以确定摩擦衰减系数为： 如果不指定衰减系数，且 FACT 大于 1.0，当接触进入滑动状态时，摩擦系数会从静摩擦系数突变到动摩擦系数，这种行为类似于 CONTAC46 和CONTAC49 单元所用的动摩擦模型，因为这会导致收敛困难，所以不建议采用5.4.9.5.2 对称、不对称求解器对无摩擦、粗糙和绑定接触，接触单元刚度矩阵是对称的。
而涉及到摩擦的接触问题产生一个不对称的刚度在每次迭代使用不对称的求解器，比对称的求解器需要更多的计算时间因此 ANSYS 程序采用对称化算法通过采用这种算法大多数摩擦接触问题，能够使用对称系统的求解器来求解如果摩擦应力在整个位移场内有相当大的影响，并且摩擦应力的大小高度依赖于求解过程，则对刚度阵的任何对称近似都可能导致收敛性降低在这种情况下，选择不对称求解选项( NROPT ,UNSYM)来改善收敛性5.4.9.6 选择接触检查的位置接触检查点位于接触单元的积分点上在积分点上，接触单元不穿透进入目标面然而，目标面能穿透进入接触面见 图 5-10 图 5-10 接触检查点位于高斯积分点上图 5-11 接触检查点位于节点上ANSYS 面─面接触单元使用高斯积分点作为缺省值，高斯积分点通常会比 Newton-Cotes/Lobatto 节点积分方案产生更精确的结果，Newton-cotes/Lobatto 用节点本身作为积分点通过 KEYOPT(4)来选择用户想使用的方法这一选项仅适用于低阶接触( CONTAC171 和 CONTAC173)然而，使用节点本身作为积分点仅应该用于角接触问题(看 图 5-11 )。
注意，使用节点作为接触检查点可能会导致其它收敛性问题，例如“滑脱”(节点滑出目标面的边界)，见 图 5-12 对大多数的点─面的接触问题，我们推荐使用其它的点─面的接触单元，例如 CONTAC26、CONTAC48 和 CONTAC49见本书§5.5 图 5-12 节点滑脱5.4.9.7 调整初始接触条件在动态分析中，刚体运动一般不会引起问题然而在静力分析中，当物体没有足够的约束时会产生刚体运动，有可能引起错误而终止计算在仅仅通过接触的出现来约束刚体运动时，必须保证在初始几何体中，接触对是接触的换句话说，用户要建立模型以便接触对是“刚好接触”的然而这样做，可能会遇到以下问题：刚体外形常常是复杂的，很难决定第一个接触点发生在哪儿既使实体模型在初始时处于接触状态，在网格划分后由于数值舍入误差，两个面的单元网格之间也可能会产生小缝隙接触单元的积分点和目标单元之间可能有小缝隙同理，在目标面和接触面之间可能发生过大的初始穿透在这种情况下，接触单元可能会高估接触力，导致不收敛或接触面之间脱开接触关系定义初始接触也许是建立接触分析模型时最重要的方面因此，程序提供了几种方法来调整接触对的初始接触条件注意 ：下面的技巧可以在开始分析时独立执行，或几个联合起来执行。
它们是为了消除由于生成网格造成的数值舍入误差而引起的小间隙或穿透，而不是为了改正网格或几何数据的错误1、应用实常数 CNOF 来指定一。