ansys接触非线性.ppt
ANSYS workbench 接触简介
Contact Status Open p
Contact Status Open p
Penetration
Gap
Penetration
Gap
Closed 法向拉格朗日法
3-8
Closed 罚函数法
Penetration
Workbench Mechanical - Introduction to Contact
积分点探测
节点探测
3-9
Workbench Mechanical - Introduction to Contact
... 接触公式
Training Manual
• 对于特定的“绑定”和“不分离”两个面间的接触类型,可用多点约束 (MPC) 算法.
– MPC 内部添加约束方程来“联结”接触面间的位移 – 这种方法不基于罚函数法或Lagrange乘子法 .它是直接有效的关联绑定接触面 方式。 – MPC算法基础的绑定接触也支持大变形效应
– “N “Normal lL Lagrange” ” 法是因为拉格朗日乘子公式用于法线方向而罚函 数 方法用于切线方向而得名的.
3-12
Workbench Mechanical - Introduction to Contact
… 公式对比
Training Manual
• 下表列出不同接触算法的比较，用“+”表优势，“-”表劣势:
F
当接触协调性不被强制时会发生渗透.
Contact
Target
F
3-4
Workbench Mechanical - Introduction to Contact
B. 接触公式
• 对非线性接触实体表面, 可使用罚函数或增强拉格朗日公式:
ansys几何非线性+塑性+接触+蠕变
除去蠕变，这个模型的结果可靠性是不错的。
作了一系列接触问题，通过试验验证符合的很好。
模型解释：（1）一个弹性结构受压（接触）变形，到发生塑性变形。
（2）拿开压缩板，结构回弹，但不会回到原始位置。
（3）这时计算蠕变，释放掉应力。
（4）再压弹性结构到开始压缩位置。
比较这四步的接触力。
结果：第二，三步当然没有接触力，（若没有应力释放，第一、第四步接触力应一样，）有了应力释放，第四步接触力比第一步减小。
这个模型中的蠕变没用太好。
用的是隐式6号蠕变方程，蠕变是时间和应力的函数，参数是乱定的（应力释放太快）。
想请教有关蠕变方面的资料，尤其是材料蠕变方程选用及参数方面的资料。
/prep7!------------CuSn8----------ET,1,182,,,3mp,ex,1,115e9mp,prxy,1,0.3r,1,0.3TB,BKIN,1TBDA TA,1,470E6,0tm=100*SET,C1,1.5625E-14 !ASSIGN VALUE*SET,C2,1.5 !ASSIGN V ALUE*SET,C3, !ASSIGN V ALUE*SET,C4,0 !ASSIGN V ALUETB,CREEP,1,,,6 !ACTIV ATE DA TA TABLETBDA TA,1,C1,C2,C3,C4 !DEFINE DATA FOR TABLE!-----------contact-----------------ET,9,169ET,10,171R,9,,,0.1,0.1,,!RMORE,,,1.0E20,0.0,1.0,!RMORE,0.0,0,1.0,0,0,0.5!RMORE,,,1.0,0.0MP,MU,9,0.0!----------------geometryk,,2k,,2,0.2k,,,0.2k,,-0.2k,,-0.2,1.2k,,,1k,,2,1.2k,,1,1k,,1.25,1k,,2,1L,8,9,k,,1.5,1.2k,,1.75,1.45L, 1, 2 L, 1, 4 L, 4, 5 L, 5, 11larc,7,12,11,0.25 larc,11,12,7,0.25L, 7, 10 L, 10, 9 L, 8, 6 L, 6, 3 L, 3, 2LFILLT,11,10,0.3, ,!*LFILLT,4,5,0.5, ,!*LFILLT,11,12,0.3, ,!*LFILLT,4,3,0.5, ,FLST,2,16,4 FITEM,2,12 FITEM,2,15 FITEM,2,11 FITEM,2,13FITEM,2,10FITEM,2,1FITEM,2,9FITEM,2,8FITEM,2,7FITEM,2,6FITEM,2,5FITEM,2,14FITEM,2,4FITEM,2,16FITEM,2,3FITEM,2,2AL,P51Xrect,1,3,1.45+0.001,1.5type,1mat,1esize,0.05amesh,all!---------contact------------ allstype,10mat,9real,9lsel,s,,,6,7nsll,s,1esln,s,0esurf,alltype,9mat,9real,9lsel,s,,,17nsll,s,1esln,s,0esurf,all!------boundarylsel,s,,,3nsll,,1d,all,uxd,all,uylsel,s,,,19nsll,,1cp,11,uy,allcplgen,11,ux*get,nmin,node,,num,min d,nmin,uxksel,s,,,10nslk*get,ndis,node,,num,minfini/soluantype,staticnlgeom,onautots,onallssaverate,offtime,1e-8d,nmin,uy,-0.3nsub,20outres,all,allsolve*get,rf1,node,nmin,rf,fy *get,dis1,node,ndis,u,ytime,2e-8d,nmin,uy,0.0nsub,20outres,all,allsolve*get,rf2,node,nmin,rf,fy*get,dis2,node,ndis,u,y!BFUNIF,TEMP,90rate,onTIME,tm!NSUBST,10OUTPR,BASIC,10 ! PRINT BASIC SOLUTION FOR EVERY 10TH SUBSTEP OUTRES,ESOL,1 ! STORE ELEMENT SOLUTION FOR EVERY SUBSTEP SOLVE*get,rf3,node,nmin,rf,fy*get,dis3,node,ndis,u,yrate,offtime,tm+1e-8d,nmin,uy,-0.3nsub,20outres,all,allsolve*get,rf4,node,nmin,rf,fy*get,dis4,node,ndis,u,y/EOFtime,11d,nmin,uy,-0.0nsub,20outres,all,allsolve*get,rf11,node,nmin,rf,fy*get,dis11,node,ndis,u,y/eoffini/post1*get,rf2,node,nmin,rf,fy fini/eof。
ANSYS非线性
如果一个凸面与一个平面或凹面进入接触, 平面和凹面应该是目标面.
如果一个面比另一个面更硬, 较硬的面应该是目标面.
如果一个面是高阶, 另一个面是低阶, 低阶面应该是目标面.
如果一个面比另外一个面更大, 较大的面应该是目标面.
当指定柔性体-柔性体接触的接触面和目标面时, 如果一个面网格粗,
刚性表面
变形体
柔性体 -柔性体接触
花键轴过盈配合, 两个零件 都是柔体.
接触协调
实际接触体相互不穿透. 因此, 程序必须在这两个面间建立一种关
系, 防止它们在有限元分析中相互穿过. • 当程序防止相互穿透时, 称之为强制接触协调.
F
当没有强制接触协调时, 发生穿透.
接触面
目标面
F
接触协调 – 罚函数法
图示收敛信息
在图形窗口显示的 是图形化的收敛历 史。显示了时间、
迭代步数与不平衡
量的信息。在求解 过程中这一显示不 断更新。
非线性求解过程
下面列出了完成非线性分析所需的典型步骤：
1. 指定分析类型
2. 指定几何非线性打开或关闭 3. 为载荷步指定“ 时间 ”
4. 设定子步数
5. 施加载荷与边界条件 6. 指定输出控制与监视值 7. 保存数据库 8. 求解载荷步
结果。
输出文件的信息
在非线性求解过程中，输出窗口显示许多关于收敛的信息。输出 窗口包括： • 力/力矩不平衡量 {R}
FORCE CONVERGENCE VALUE
• 最大的自由度增量 {u}
MAX DOF INC
• 力收敛判据
CRITERION
• 载荷步与子步数
LOAD STEP 1 SUBSTEP 14
ANSYS讲义非线性分析
• 力/力矩不平衡量 {R}
FORCE CONVERGENCE VALUE
• 最大的自由度增量 {u}
MAX DOF INC
• 力收敛判据
CRITERION
• 载荷步与子步数
LOAD STEP
1 SUBSTEP 14
第35页/共97页
输出文件的信息(续)
输出窗口包括（续） ：
• 当前子步的迭代步数
EQUIL ITER 1 COMPLETED. NEW TRIANG MATRIX. MAX DOF INC= -0.1645E-01
FORCE CONVERGENCE VALUE = 10.35
CRITERION= 2.095
DISP CONVERGENCE VALUE = 0.2409E-01 CRITERION= 0.9406
第11页/共97页
概述
ANSYS 最常用的非线性功能
• 几何非线性
大应变，大位移与大转动
-- 结构稳定性 （前屈曲分析与后屈曲分析）
• 塑性
• 超弹性
• 接触非线性
第12页/共97页
非线性分析
第13页/共97页
主要内容
应理解非线性分析中所用到的基本术语：
1. Newton-Raphson法 2. 收敛 3. 载荷步，子步和平衡迭代 4. 自动时间步 5. 输出文件信息 6. 非线性求解过程 7. 高级求解控制 8. 重启动分析
1. 几何非线性
大应变，大位移，大旋转
2. 材料非线性
塑性，超弹性，粘弹性，蠕变
3. 状态改变非线性
接触，单元死活
第3页/共97页
几何非线性
如果一个结构承受大的变形，它改变的几何构形可导致非线性行 为。大位移、大应变和大旋转是几何非线性的例子。
ansys非线性接触分析中的接触行为
ansys非线性接触分析中接触行为接触是状态改变非线性，经典ANSYS版本中共提供了7种接触行为，每一种都有其特点及相应的应用范围，在选用的时候应该谨慎。
（1）标准接触行为（standard）该接触行为包括了法向接触闭合和分开行为，在该接触模式中既考虑粘着摩擦同时也考虑了滑动摩擦。
如图上，AB与BC本来是分开的，中间通过B点连接，当在A点施加力F，AB慢慢贴近BC，最终靠在一起。
但F撤销后，AB在恢复力的作用下慢慢回复到初始分开状态。
标准接触行为包括了分开状态→闭合状态→分开状态。
当AB与BC靠在一起时，既存在正压力，同时还有沿BC圆弧切线方向的摩擦力。
（2）粗糙接触行为（rough）该接触行为包括了法向接触闭合和分开行为，但滑动行为在此是不会发生的。
原因是所有参与接触的表面都被假定为非常粗糙，以致于可以认为摩擦力无穷大而不能够产生相对滑动。
在这种接触行为中，接触的两个物体或部件之间，除了存在正压力外，还有切向摩擦力，但是接触部分之间不可以产生相对滑动。
（3）绑定接触行为（bonded）是指一旦接触关系建立，那么目标面及接触面就被假定为粘结在一起（不可以分开）。
（4）绑定接触行为（始终）（bonded（always））任何初始时在许可接触容差范围内探测到的接触点或者是那些即将进入接触的点在后续的分析中将被绑定在一起。
这种接触行为的典型应用，如在组装分析中将两种不同网络的组件“加”在一起。
线性静态分析也可以用该种接触行为来解决，虽然由于有接触单元的存在，分析中将会提示为非线性分析，但往往只要一步迭代就完成了。
（5）绑定接触行为（初始接触）（bonded（initial））绑定仅发生在初始状态下就接触的面上，初始状态下没有接触的部分将继续保持分开。
典型的例子是通过焊接连接在一起的两个物体，焊接部分始终保持连接，没有焊接的部分保持分离状态。
（6）不分开型（no separation）一旦接触关系建立，目标面及接触面便被约束在一起了，但还是允许接触面之间有滑动。
非线性_接触分析ANSYS
接触分析接触问题是一种高度非线性行为，需要较大的计算资源，为了进行实为有效的计算，理解问题的特性和成立合理的模型是很重要的。
接触问题存在两个较大的难点：其一，在你求解问题之前，你不知道接触区域，表面之间是接触或分开是未知的，突然转变的，这随载荷、材料、边界条件和其它因素而定；其二，大多的接触问题需要计算摩擦，有几种摩擦和模型供你挑选，它们都是非线性的，摩擦使问题的收敛性变得困难。
一般的接触分类接触问题分为两种大体类型：刚体－柔体的接触，和柔体－柔体的接触，在刚体－柔体的接触问题中，接触面的一个或多个被看成刚体，（与它接触的变形体相较，有大得多的刚度），一般情况下，一种软材料和一种硬材料接触时，问题可以被假定为刚体－柔体的接触，许多金属成形问题归为此类接触，另一类，柔体－柔体的接触，是一种更普遍的类型，在这种情况下，两个接触体都是变形体（有近似的刚度）。
ANSYS接触能力ANSYS支持三种接触方式：点－点，点－面，和面－面，每种接触方式利用的接触单元适用于某类问题。
为了给接触问题建模，首先必需熟悉到模型中的哪些部份可能会彼此接触，若是彼此作用的其中之一是一点，模型的对应组元是一个结点。
若是彼此作用的其中之一是一个面，模型的对应组元是单元，例如梁单元，壳单元或实体单元，有限元模型通过指定的接触单元来识别可能的接触匹对，接触单元是覆盖在分析模型接触面之上的一层单元，至于ANSYS 利用的接触单元和利用它们的进程，下面分类详述。
点－点接触单元点－点接触单元主要用于模拟点－点的接触行为，为了利用点－点的接触单元，你需要预先知道接触位置，这种接触问题只能适用于接触面之间有较小相对滑动的情况（即便在几何非线性情况下）若是两个面上的结点一一对应，相对滑动又以忽略不计，两个面挠度（转动）维持小量，那么可以用点－点的接触单元来求解面－面的接触问题，过盈装配问题是一个用点－点的接触单元来模拟面－面的接触问题的典型例子。
ANSYS非线性接触问题分析汇总
ANSYS 分线性接触问题分析汇总接触非线性是一门复杂的学科，ANSYS 关于计算非线性接触的设置选项多只又多，很多人摸不到头脑，本文就基于ANSYS 模拟过的几个接触实例，研究了相关设置选项对接触结果的影响。
实例1：橡胶密封圈配合接触研究—非线性求解设置对结果的影响密封圈配合模型简图见图1，左右两端为刚体，中间圆部分为橡胶密封圈，将刚体2沿刚体1方面移动，从而实现橡胶圈密封作用，采用plane182单元，设置轴对称行为，建立橡胶密封圈与刚体接触模型，见图2。
图1 密封圈配合模型简图 图2 密封圈配合有限元模型图接触对采用默认设置，摩擦系数取0.10，研究非线性求解器设置对收敛方面的影响，大变形静态（Large Displacement Static ）效应打开，自动时间步长（Automatic time stepping ）打开，子步数（Number of substeps ）设置为50，线性搜索（Line search ）打开。
1 收敛准则对结果的影响此实例收敛准则默认采用力收敛结合力矩收敛准则（基于L2范数），收敛容差（Tolerance ）默认为0.001，工程上认为0.05的收敛容差足够满足要求。
表 1 收敛容差对计算结果的影响收敛容差 最大应力/ MPa报错与否？ 0.001 4.12364报错 0.05 4.12785 报错 0.14.12996报错查看报错信息，见图3，表示单元过于扭曲，建议提高子步数或降低时间步长，需要提高网格质量，也要考虑材料属性，接触对及约束方程的合理性，若在第一步迭代就如此，需要预先执行单元形状检查。
图3 报错信息刚体1刚体2密封圈橡胶密封圈配合Von Mises应力云图见图4。
图4 橡胶密封圈配合Von Mises应力2 子步数对结果的影响此实例子步数设置为50、100、200、500，收敛容差（Tolerance）默认为0.001，研究子步数对收敛的影响。
ANSYS非线性接触问题分析汇总
ANSYS 分线性接触问题分析汇总接触非线性是一门复杂的学科，ANSYS 关于计算非线性接触的设置选项多只又多，很多人摸不到头脑，本文就基于ANSYS 模拟过的几个接触实例，研究了相关设置选项对接触结果的影响。
实例1：橡胶密封圈配合接触研究—非线性求解设置对结果的影响密封圈配合模型简图见图1，左右两端为刚体，中间圆部分为橡胶密封圈，将刚体2沿刚体1方面移动，从而实现橡胶圈密封作用，采用plane182单元，设置轴对称行为，建立橡胶密封圈与刚体接触模型，见图2。
图1 密封圈配合模型简图 图2 密封圈配合有限元模型图接触对采用默认设置，摩擦系数取0.10，研究非线性求解器设置对收敛方面的影响，大变形静态（Large Displacement Static ）效应打开，自动时间步长（Automatic time stepping ）打开，子步数（Number of substeps ）设置为50，线性搜索（Line search ）打开。
1 收敛准则对结果的影响此实例收敛准则默认采用力收敛结合力矩收敛准则（基于L2范数），收敛容差（Tolerance ）默认为0.001，工程上认为0.05的收敛容差足够满足要求。
表 1 收敛容差对计算结果的影响收敛容差 最大应力/ MPa报错与否？ 0.001 4.12364报错 0.05 4.12785 报错 0.14.12996报错查看报错信息，见图3，表示单元过于扭曲，建议提高子步数或降低时间步长，需要提高网格质量，也要考虑材料属性，接触对及约束方程的合理性，若在第一步迭代就如此，需要预先执行单元形状检查。
图3 报错信息刚体1刚体2密封圈橡胶密封圈配合Von Mises应力云图见图4。
图4 橡胶密封圈配合Von Mises应力2 子步数对结果的影响此实例子步数设置为50、100、200、500，收敛容差（Tolerance）默认为0.001，研究子步数对收敛的影响。
Ansys非线性接触分析和设置
Ansys非线性接触分析和设置设置实常数和单元关键选项程序使用20个实常数和数个单元关键选项，来控制面─面接触单元的接触。
参见《ANSYS Elements Reference》中对接触单元的描述。
实常数在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 静、动摩擦衰减系数。
命令：RGUI：main menu> preprocessor>real constant对实常数FKN, FTOLN, ICONT, PINB, PMAX, PMIN, FKOP 和FKT，用户既可以定义一个正值，也可以定义一个负值。
程序将正值作为比例因子，将负值作为绝对值。
程序将下伏单元的厚度作为ICON，FTOLN，PINB，PMAX 和PMIN 的参考值。
例如ICON = 表明初始闭合因子是“*下层单元的厚度”。
然而，ICON = 则表示真实调整带是单位。
如果下伏单元是超单元，则将接触单元的最小长度作为厚度。
参见图5-8。
图5-8 下层单元的厚度在模型中，如果单元尺寸变化很大，而且在实常数如ICONT, FTOLN, PINB, PMAX, PMIN 中应用比例系数，则可能会出现问题。
ANSYS高级接触问题处理
• －开始估计时，选用 • FKN = 1.0 大面积实体接触 • FKN = 0.01-0.1 较柔软（弯曲占主导的）部分 • －另外，也可以指定一个绝对刚度值，单位：（力/长度）
/ 面积。 • ·点一点（除CONTA178）和点－面接触单元需要为罚刚
度KN输入绝对值： • －初始估计时： • 对于大变形: 0.1*E < KN < 1.0*E • 对于弯曲: 0.01*E < KN < 0.1*E • E 为弹性模量
ANSYS 高级接触问题
• 接触问题概述 • 在工程中会遇到大量的接触问题，如齿轮的啮合、法兰
联接、机电轴承接触、卡头与卡座、密封、板成形、冲 击等等。接触是典型的状态非线性问题，它是一种高度 非线性行为。接触例子如图1：
• 分析中常常需要确定两个或多个相互接触 物体的位移、接触区域的大小和接触面上 的应力分布。
触
• 3. 点－点接触单元用于模拟单点和另一个确定点 之间的接触。
• －建立模型时必须事先知道确切的接触位置；
• －多个点－点接触单元可以模拟两个具有多个单 元表面间的接触；
• ·每个表面的网格必须是相同的； • ·相对滑动必须很小； • ·只对小的转动响应有效。
• 例如: 点一点接触可以模拟一些面的接触。如地 基和土壤的接触
• 3、选取接触刚度的指导：
• Step 1.开始采用较小的刚度值 • Step 2.对前几个子步进行计算 • Step 3.检查穿透量和每一个子步中的平衡迭代次数
• ·在粗略的检查中，如以实际比例显示整个模型时就能观察到穿透， 则穿透可能太大了，需要提高刚度重新分析。
• ·如果收敛的迭代次数过多（或未收敛），降低刚度重新分析。 • 注意：罚刚度可以在载荷步间改变，并且可以在重启动中调整。 • 牢记：接触刚度是同时影响计算精度和收敛性的最重要的参数。如
ANSYS接触非线性应用技巧
Workbench 默认算法，计算速度快
Augmented Lagrange
ANSYS 经典界面默认算法，计算精度较高，收敛性好
Pure Lagrange
点对点的接触（较少使用）
Pure Lagrange on Normal; Penalty on Tangent
常用于模拟接触面之间零穿透并且存在摩擦的情况，例 如 模拟橡胶－金属之间的接触
Single contact pair
Contact pair-1
2.接触刚度不合适
缺省的接触刚度不再适用
ANSYS求解接触问题时，接触刚度对收敛性影响显著。大多数情况下缺省设置 有效，但不适于弯曲为主的接触问题。 不好收敛的信号是，收敛曲线平行于收敛准则，如下图。选择小的罚刚度值（ FKN），可以使接触容易收敛。
Error in element formulation
案例：非线性培训教材例子
项目背景
– 使用简化的2D 平面应力模型同样收敛失败
案例：非线性培训教材例子
解决方案
– 回到 DM 界面,增加倒圆角，右下角可看到小半径
案例：非线性培训教材例子
结果后处理.
案例：某机车结构件
工程背景
– 计算橡胶材料在 挤压作用下变形 和应力应变
Multipoint Constraint (MPC)
经常和Bonded联合使用，用于大模型的装配定义，可大大 降低计算时间，还或用于不同类型单元之间的连接
ANSYS接触不收敛原因及对策
遗漏接触对 接触刚度大 载荷步少 奇异 网格大变形 刚体位移 振荡 其他...
1.遗漏接触对
• 如果有些区域接触了（如自接触），但是没有定义接触对，单元会扭曲。
ansys教程完整PPT教学课件
jobname.log
文本
结果文件
jobname.rxx
二进制
图形文件
jobname.grph 二进制
ANSYS的数据库，是指在前处理、求解及后处理过程中，ANSYS保存在内存中的数据。数据库既存储 输入的数据，也存储结果数据:
输入数据 - 必须输入的信息 (模型尺寸、材料属性、载荷等).
结果数据 - ANSYS计算的数值 (位移、应力、应变、温度等).
OOPs!
Lines
Keypoints
第20页/共78页
2.布尔操作
1. ..... 2. ..... 3. .....
Procedure
要使用布尔操作: Main Menu: Preprocessor > -Modeling- Operate >
选择一种布尔操作 (例如： Add)
选择图形类型. 将弹出 选取菜 单 (见下页) 提示选择图形进行 布尔操作.
+ 加载的操作更加容易 ,尤其是在图形中直接拾取时.
第24页/共78页
加载 (续)
无论采取何种加载方式,ANSYS求解前都将载荷转化到有限元模型.因此, 加载到实体的载荷将自动转化到 其所属的节点或单元上。
沿线均布的压力
实体模型
加载到实 体的载荷 自动转化 到其所属 的节点或 单元上
均布压力转化到以线为边界的 各单元上
第21页/共78页
四、加载、求解
Objective
4-1. 列表和分类载荷
ANSYS中的载荷可分为:
• 自由度DOF - 定义节点的自由度（ DOF ） 值 (结构分析_位移、热 分析_ 温度、电磁分析_磁势等)
• 集中载荷 - 点载荷 (结构分析_力、热分析_ 热导率、电磁分析_
ANSYS接触问题
接触问题（参考ANSYS的中文帮助文件）当两个分离的表面互相碰触并共切时，就称它们牌接触状态。
在一般的物理意义中，牌接触状态的表面有下列特点：1、不互相渗透；2、能够互相传递法向压力和切向摩擦力；3、通常不传递法向拉力。
接触分类：刚性体－柔性体、柔性体－柔性体实际接触体相互不穿透，因此，程序必须在这两个面间建立一种关系，防止它们在有限元分析中相互穿过。
――罚函数法。
接触刚度――lagrange乘子法，增加一个附加自由度（接触压力），来满足不穿透条件――将罚函数法和lagrange乘子法结合起来，称之为增广lagrange法。
三种接触单元：节点对节点、节点对面、面对面。
接触单元的实常数和单元选项设置：FKN：法向接触刚度。
这个值应该足够大，使接触穿透量小；同时也应该足够小，使问题没有病态矩阵。
FKN值通常在0.1~10之间，对于体积变形问题，用值1.0（默认），对弯曲问题，用值0.1。
FTOLN：最大穿透容差。
穿透超过此值将尝试新的迭代。
这是一个与接触单元下面的实体单元深度（h）相乘的比例系数，缺省为0.1。
此值太小，会引起收敛困难。
ICONT：初始接触调整带。
它能用于围绕目标面给出一个“调整带”，调整带内任何接触点都被移到目标面上；如果不给出ICONT值，ANSYS根据模型的大小提供一个较小的默认值（＜0.03＝PINB：指定近区域接触范围（球形区）。
当目标单元进入pinball区时，认为它处于近区域接触，pinball区是围绕接触单元接触检测点的圆（二维）或球（三维）。
可以用实常数PINB调整球形区（此方法用于初始穿透大的问题是必要的）PMIN和PMAX：初始容许穿透容差。
这两个参数指定初始穿透范围，ANSYS 把整个目标面（连同变形体）移到到由PMIN和PMAX指定的穿透范围内，而使其成为闭合接触的初始状态。
初始调整是一个迭代过程，ANSYS最多使用20个迭代步把目标面调整到PMIN和PMAX范围内，如果无法完成，给出警告，可能需要修改几何模型。