Ansys非线性接触分析和设置.doc

1、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 指定施

2、加于接触面的正或负的偏移值。FKOP 指定在接触分开时施加的刚度系数。FKT 指定切向接触刚度。COHE 制定滑动抗力粘聚力。TCC 指定热接触传导系数。FHTG 指定摩擦耗散能量的热转换率。SBCT 指定 Stefan-Boltzman 常数。RDVF 指定辐射观察系数。FWGT 指定在接触面和目标面之间热分布的权重系数。FACT 静摩擦系数和动摩擦系数的比率。DC 静、动摩擦衰减系数。命令： R GUI：main menu preprocessorreal constant对实常数 FKN, FTOLN, ICONT, PINB, PMAX, PMIN, FKOP 和 FKT，用户既可以定

3、义一个正值，也可以定义一个负值。程序将正值作为比例因子，将负值作为绝对值。程序将下伏单元的厚度作为ICON，FTOLN，PINB，PMAX 和 PMIN 的参考值。例如 ICON = 0.1 表明初始闭合因子是“0.1*下层单元的厚度”。然而，ICON = -0.1 则表示真实调整带是 0.1 单位。如果下伏单元是超单元，则将接触单元的最小长度作为厚度。参见 图5-8 。图5-8 下层单元的厚度在模型中，如果单元尺寸变化很大，而且在实常数如 ICONT, FTOLN, PINB, PMAX, PMIN 中应用比例系数，则可能会出现问题。因为从比例系数得到的实际结果，取决于下层单元的厚度，这就可

4、能引起大、小单元之间的重大变化。如果出现这一问题，请用绝对值代替比例系数。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

5、(7)伪接触预防 KEYOPT(8)初始穿透或间隙的影响 KEYOPT(9)法向和切向接触刚度修正方法控制 KEYOPT(10)壳的厚度影响 KEYOPT(11)接触面行为(粗糙、绑定等) KEYOPT(12)命令： KEYOPT ET GUI：main menupreprocessorElemant TypeAdd/Edit/Delete5.4.9.3 选择接触算法对面面接触单元，程序可以使用增进的拉格朗日方法或罚函数方法。通过单元关键字 KEYOPT(2)来指定。增进的拉格朗日方法是为了找到精确的拉格朗日乘子（即接触力），而对罚函数进行一系列修正迭代。与罚函数的方法相比，拉格朗日方法容易得

6、到良态条件，对接触刚度的敏感性较小。然而，在有些分析中，增进的拉格朗日方法可能需要更多的迭代，特别是在变形后网格变得太扭曲时。使用拉格朗日方法的同时应使用实常数 FTOLN。FTOLN 为拉格朗日方法指定容许的最大穿透。如果程序发现穿透大于此值时，即使不平衡力和位移增量已经满足了收敛准则，总的求解仍被当作不收敛处理。FTLON 的缺省值为0.1。用户可以改变这个值，但要注意，如果此值太小，可能会造成太多的迭代次数或者不收敛。5.4.9.4 确定接触刚度所有的接触问题都需要定义接触刚度，两个表面之间穿透量的大小取决于接触刚度。过大的接触刚度可能会引起总刚矩阵的病态，从而造成收敛困难。一般来说，应

7、该选取足够大的接触刚度以保证接触穿透小到可以接受，但同时又应该让接触刚度足够小以不致引起总刚矩阵的病态而保证收敛性。ANSYS 程序根据下伏柔体单元的材料特性，来估计一个缺省的接触刚度值。用户可用实常数 FKN 来为接触刚度指定一个比例因子或指定一个绝对值。比例因子一般在0.01和10之间；对于大变形问题，选1是比较好的；而对于弯曲为主的问题，通常为0.010.1。用户应当总是检验以使穿透到达极小值，而又避免过多的迭代次数。注意 -FTOLN 和 FKN 从一个荷载步到另一个荷载步中，都可以修改。也可以在重启动中修改。这时，必须定义KEYOPT(10)=1,2。为了确定一个较好的接触刚度值，可

8、能需要一些经验。用户可以按下面的步骤来进行尝试：1、 开始时取一个较低的值。低估值要比高估值好，因为由一个较低的接触刚度导致的穿透问题，比过高的接触刚度导致的收敛性困难，要容易解决。2、 对前几个子步进行计算分析，直到最终荷载的一个比例(刚好完全建立接触)。3、 检查每一子步中的穿透量和平衡迭代次数。如果总体收敛困难是由过大的穿透引起的(而不是由不平衡力和位移增量引起的)，那么可能低估了FKN 的值，或者是将 FTOLN 的值取得大小。如果总体的收敛困难是由于不平衡力和位移增量达到收敛值时需要过多的迭代次数，而不是由于过大的穿透量引起的，那么 FKN 的值可能被高估。4、 按需要调整 FKN

9、或 FTOLN 的值，重新进行完整的分析。注意 - 如果穿透控制变成总体平衡迭代中的主因(如果为使问题收敛到穿透容差内，比收敛到不平衡力的容差内，需要更多的迭代)，用户应该增大 FTOLN 值，以允许更多的穿透，或增大 FKN。5.4.9.5 选择摩擦类型在基本的库仑摩擦模型中，两个接触面在开始相互滑动之前，在它们的界面上会有达到某一大小的剪应力产生。这种状态称为粘合状态(stick)。库仑摩擦模型定义了一个等效剪应力，在某一法向压应力p作用下剪应力达到此值时表面开始滑动 (=p+COHE，其中是摩擦系数-MU-作为材料特性定义，而 COHE 是粘聚力)。一旦剪应力超过此值后，两个表面之间将开

10、始相互滑动。这种状态，叫作滑动状态(Sliding)。粘合/滑动计算决定什么时候一个点从粘合状态到滑动状态，或从滑动状态变到粘合状态。摩擦系数可以是任一非负值，程序缺省值为表面之间无摩擦。对于粗糙或绑定接触( KEYOPT(12)=1、3、5、6)，程序将不管给定的 MU 值而认为摩擦阻力无限大。程序提供了一个人为指定最大等效剪应力的选项，不管接触压力值的大小，如果等效剪应力达到此值时，即发生滑动，见 图5-9 。为了指定接触界面上最大容许剪应力，设置常数 TAUMAX (缺省为1.0E20)。这个剪应力极限，通常用于在接触压力非常大的时候(如在某些加工过程中)的一些情况，以至于用库仑理论计算

11、出的界面剪应力超过了材料的屈服极限。TAUMAX 的一个合理上限估值为 是表面附近材料的 von Mises屈服应力)。经验数据有助于决定 TAUMAX 的值。图5-9 摩擦模式5.4.9.5.1 静、动摩擦系数摩擦系数依赖于接触面的相对滑动速度，通常静摩擦系数高于动摩擦系数。ANSYS提供了如下表示的指数衰减摩擦模型：=MU(1+(FACT-1)exp(-DCV rel )其中： 为摩擦系数。 MU动摩擦系数，用MP命令输入。 FACT是静摩擦系数与动摩擦系数之比，缺省为最小值1.0。 DC为衰减系数，缺省为0.0，单位为time/length。因此，时间在静态分析中有一些意义。 V rel

12、 是ANSYS计算的滑动速度。如果知道静、动摩擦系数和至少一个数据点（ 1 ，V rel ），则可以确定摩擦衰减系数为： 如果不指定衰减系数，且FACT大于1.0，当接触进入滑动状态时，摩擦系数会从静摩擦系数突变到动摩擦系数，这种行为类似于CONTAC46和CONTAC49单元所用的动摩擦模型，因为这会导致收敛困难，所以不建议采用。5.4.9.5.2 对称、不对称求解器对无摩擦、粗糙和绑定接触，接触单元刚度矩阵是对称的。而涉及到摩擦的接触问题产生一个不对称的刚度。在每次迭代使用不对称的求解器，比对称的求解器需要更多的计算时间。因此ANSYS程序采用对称化算法。通过采用这种算法大多数摩擦接触问题

13、，能够使用对称系统的求解器来求解。如果摩擦应力在整个位移场内有相当大的影响，并且摩擦应力的大小高度依赖于求解过程，则对刚度阵的任何对称近似都可能导致收敛性降低。在这种情况下，选择不对称求解选项( NROPT ,UNSYM)来改善收敛性。5.4.9.6 选择接触检查的位置接触检查点位于接触单元的积分点上。在积分点上，接触单元不穿透进入目标面。然而，目标面能穿透进入接触面。见 图5-10 。图5-10 接触检查点位于高斯积分点上图5-11 接触检查点位于节点上ANSYS面面接触单元使用高斯积分点作为缺省值，高斯积分点通常会比 Newton-Cotes/Lobatto 节点积分方案产生更精确的结果，

14、Newton-cotes/Lobatto 用节点本身作为积分点。通过KEYOPT(4)来选择用户想使用的方法。这一选项仅适用于低阶接触( CONTAC171 和 CONTAC173)。然而，使用节点本身作为积分点仅应该用于角接触问题(看 图5-11 )。注意，使用节点作为接触检查点可能会导致其它收敛性问题，例如“滑脱”(节点滑出目标面的边界)，见 图5-12 。对大多数的点面的接触问题，我们推荐使用其它的点面的接触单元，例如CONTAC26、CONTAC48 和 CONTAC49。见本书5.5。 图5-12 节点滑脱5.4.9.7 调整初始接触条件在动态分析中，刚体运动一般不会引起问题。然而在

15、静力分析中，当物体没有足够的约束时会产生刚体运动，有可能引起错误而终止计算。在仅仅通过接触的出现来约束刚体运动时，必须保证在初始几何体中，接触对是接触的。换句话说，用户要建立模型以便接触对是“刚好接触”的。然而这样做，可能会遇到以下问题： 刚体外形常常是复杂的，很难决定第一个接触点发生在哪儿。 既使实体模型在初始时处于接触状态，在网格划分后由于数值舍入误差，两个面的单元网格之间也可能会产生小缝隙。 接触单元的积分点和目标单元之间可能有小缝隙。同理，在目标面和接触面之间可能发生过大的初始穿透。在这种情况下，接触单元可能会高估接触力，导致不收敛或接触面之间脱开接触关系。定义初始接触也许是建立接触分

16、析模型时最重要的方面。因此，程序提供了几种方法来调整接触对的初始接触条件。注意 ：下面的技巧可以在开始分析时独立执行，或几个联合起来执行。它们是为了消除由于生成网格造成的数值舍入误差而引起的小间隙或穿透，而不是为了改正网格或几何数据的错误。1、应用实常数 CNOF 来指定一个接触面偏移。指定正的值来使整个接触面偏向目标面。 指定负的值来使接触面离开目标面。ANSYS 能够自动提供 CNOF 值到刚好闭合间隙或减少初始穿透。如下设置KEYOPT(5)：=1: 闭合间隙；=2: 减少初始穿透；=3: 闭合间隙或减少初始穿透。如果设置了 KEYOPT(5)0 ，则 ICONT 缺省值为0。2、使用实

17、常数 ICONT 来指定一个小的初始接触环，初始接触环是指沿着目标面的“调整环”的深度。如果没有人为指定 ICONT 的值，程序会根据几何尺寸来给 ICONT 提供一个小值(但有意义的值)，同时输出一个表示什么值被指定的警告信息。 ICONT 正值表示相对于下层单元厚度的比例因子；负值表示接触环的绝对值。任何落在“调整环”域内的接触检查点被自动移到目标面上，(参见 图5-13a )。建议使用一个十分小的 ICONT 值，否则可能会发生严重不连续(看 图5-13b )图5-13 用ICON进行接触面的调整。(a)调整前;(b)调整后CNOF 与 ICONT 之间的差别，是前者把整个接触面移动 C

18、NOF 的距离，而后者把所有初始分开的(刚好位于调整环 ICONT 内的)接触点向目标面移动。如果用户应用其他方法来平衡初始未约束的自由体(如 FTOLN, PINB, PMAX 和 PMIN)，基本上消除ICONT 的影响(把它关联一个很小的值，如1E-20)，是一个好办法。但是，设置 ICONT=0，并不会关闭它。而是导致 ANSYS 用缺省值代替。当与其他约束自由体的方法联用时，这反过来可能不能达到效果。3、使用实常数 PMIN 和 PMAX 来指定初始容许的穿透范围。当指定 PMAX 或 PMIN 后，在开始分析时，程序会将目标面移到初始接触状态，见 图5-14 。如果初始穿透大于 P

19、MAX，程序会调整目标面来减少穿透。接触状态的初始调节仅仅通过平移来实现。对给定载荷或给定位移的刚性目标面，将会执行初始接触状态的初始调节。对没有指定边界条件的目标面，也同样可以进行初始接触的调整。当目标面上的所有节点，有给定的零位移值时，使用 PMAX 和 PMIN 的初始调节将不会被执行。图5-14 接触面调整(PMIN，PMAX)注意 -ANSYS程序独立地处理目标面上节点的自由度。例如：如果用户指定自由度 UX 值为“0”，那么沿着X方向就没有初始调查。然而，在 Y 和 Z 方向仍然会激活 PMAX 和 PMIN 选项。初始状态调整是一个迭代过程，程序最多进行20次迭代。如果目标面不能

20、进入可接受的穿透范围(即 PMIN，PMAX 范围)，程序将在原始几何实体上操作。这时程序会给出一个警告信息，用户可能需要调整用户的初始几何模型。图5-15 给出了一个初始接触调整迭代失败的例子。目标面的 UY 被约束。因此，初始接触唯一容许的调整是在 X 方向，然而，在这个问题中，刚性目标面在 X 方向的任何运动都不会引起初始接触。对于柔体-柔体接触，这种方法不仅移动整个目标面，还同时移动与目标面相连的整个柔体。请确保没有其他接触面或目标面与柔体相连。图5-15 一个初始调整失败的例子4、设置KEYOPT(9)=1来调整初始穿透或间隙，见 图5-16 。图5-16 忽略初始穿透，KEYOPT

21、(9)=1真正的初始穿透包括两部分： 几何模型产生的穿透或间隙； 用户定义的接触面偏移(CNOF)产生的穿透或间隙；KEYOPT(9)提供下列功能： 包括由几何和接触面偏移产生的初始穿透，设置 KEYOPT(9)=0。这是缺省。 忽略上面两者引起的初始穿透，设置KEYOPT(9)=1。在KEYOPT(12)=4或5时，这一KEYOPT(9)=1，也将忽略间隙弹簧的初始力，这样，建立了一个初始的“理想的”接触面-在接触截面上没有初始力的作用。 为了包括定义的接触面偏移(CNOF)，但忽略由于几何模型引起的初始穿透，设置KEYOPT(9)=3。在KEYOPT(12)=4或5时，这一KEYOPT(9

22、)=3，也将忽略打开间隙弹簧的初始力，这样，建立了一个初始的“理想的”接触面-在接触截面上没有初始力的作用。在某些情况下，例如过盈装配问题，期望有过度的穿透。如在第一个载荷步施加阶跃初始穿透，可能造成收敛困难。为了缓解收敛性困难，在第一个载荷步中设置渐变的初始穿透来克服。见 图5-17 。 下面的KEYOPT(9)设置用来提供渐变功能： 设置KEYOPT(9)=2，来渐变施加初始穿透(CNOF+由于几何模型造成的偏移)。 设置KEYOPT(9)=4，来渐变施加接触面穿透，但忽略由于几何模型造成的穿透。对于上面两个 KEYOPT(9)设置，用户还应该设置 KBC ,0 ，并在第一个载荷步中不要给

23、定任何其它外载荷。还要确保球形（Pinball）区域足够大以捕捉到初始过盈。用户可以联合应用上面的技术。例如，用户可能希望设置十分精确的初始渗透或间隙，但有限元节点的初始坐标可能无法提供足够的精度。这时，可以： 应用 ICONT 来移动初始张开的接触点刚好碰到目标面。 应用 CNOF 来指定穿透（正值）或间隙(负值)。 应用 KEYOPT(9)=3 来在第一个子步求解初始穿透，或应用 KEYOPT(9)=4 来逐渐求解初始穿透。图5-17 渐变初始过盈在开始分析时，程序会给出每个目标面的初始接触状态的输出信息(输出窗口或输出文件中)，这个信息有助于决定每个目标面的最大穿透或最小间隙。对于给定的

24、目标面，如果没有发现接触，可能是目标面离接触面太远(超出了球形区域)，或者是接触/目标单元已经被杀死。5.4.9.8 决定接触状态和球形区域。接触单元相对于其目标面的运动和位置，决定了接触单元的状态；程序检测每个接触单元，并给出一种状态： STAT=0 未闭合的远场接触 STAT=1 未闭合的近场接触 STAT=2 滑动接触 STAT=3 粘合接触当目标面进入球形区域后，接触单元就被当作未闭合的近场接触，球形区域是以接触单元的积分点为中心的。使用实常数 PINB 来为球形区域指定一个比例因子(正值)，或其绝对值(负值)。缺省时，程序将球形区域定义为一个以“4*下层单元厚度(对于刚体-柔体接触)

25、”或“2*下层单元厚度(对于柔体-柔体接触)”为半径的圆(对2-D问题)或球(对3-D问题)。检查接触的计算时间依赖于球形区域的大小，远场接触单元计算简单，计算时间较少。近场接触计算(对于接近接触或实际接触的接触单元)较慢并且较复杂。当单元已经接触时，计算最为复杂。如果刚性面有好几个凸形区域，为了克服伪接触定义，设置一个合适的球形区域是有用的。而对大多数问题，缺省值是合适的。5.4.9.9 在自接触问题中避免伪接触在一些对称接触问题(包含自接触)中，ANSYS可能错误地假设在十分接近的几何位置上的接触面和目标面之间的接触。在自接触问题的角点会发生这一问题。它可能是由单元的初始几何位置引起的，也

26、可能在分析时通过变形而引起的。在二个面位于球形区域内，而且它们之间夹角小于90度时，会产生这一问题。在这种情况下，ANSYS程序认为发生了十分大的穿透。 图5-18 说明了这种情况。图5-18 检测伪接触用户可以通过单元 CONTA171CONTA174的KEOPT(8)=1，来防止ANSYS程序考虑这一问题。当应用这个 KEYOPT 时，ANSYS将忽略在如下情况下产生的“过度穿透”接触： 初始检测到的穿透大于接触容差(FTOLN)的20%。见 图5-18a 。 接触状态由 图5-18a 所示的远离接触突然变化到 图5-18b 所示的过度穿透。如果KEYOPT(8)激活，ANSYS各个荷步中

27、第一次检测到伪接触时会发出一个警告。如ANSYS在第1荷载步中发现这种接触，用户可看到如下的信息：“Contact element x has too much penetration related to target element y. We assume it (may be more elements) is spurious contact.”如果ANSYS检测到归类为伪接触的突变，用户看到如下信息：“Contact element x status changed abruptly with target element y. We assume it (may be more

28、elements) is spurious contact.”ANSYS在一个荷载步中仅发出一次这样的信息。在该荷载步中如果还存在其他的伪接触，ANSYS不再提醒。注意 -在应用KEYOPT(8)=1之前，请认真检查模型。ANSYS将忽略符合了过度穿透准则的真正接触。这一技术应当仅用于在改变球形区域 PINB 不能阻止伪接触的情况。5.4.9.10 修正法向和切向接触刚度在分析过程中，可以修正法向和切向接触刚度。可以自动修正(由于改变下伏单元刚度的大应变效应产生)，也可以显式地修正(由用户重新指定 FKN 或 FKT 值)。KEYOPT(10)控制法向和切向接触刚度如何修正： KEYOPT(1

29、0)=0，禁止那些已经处于“闭合”状态的单元的接触刚度修正。对于从“张开”变化到“闭合”状态的单元，将在每一个子步上修正接触刚度。 KEYOPT(10)=1，允许已处于“闭合”状态的单元的接触刚度，在荷载步之间或在重启动期间改变。对于从“张开”变化到“闭合”状态的单元，将在每一个子步上修正接触刚度。 KEYOPT(10)=2与 KEYOPT(10)=1相同，只是对所有单元(不论其状态)将在每一个子步上，由程序决定自动修正。5.4.9.11 选择表面作用模式面-面接触单元支持法向单向接触模式及其他力学表面作用模式。通过设置 KEYOPT(12) 来选择下面的某种作用模式： KEYOPT(12)=

30、0，法向单向接触，即在接触分开时，法向压力等于0。 KEYOPT(12)=1，理想粗糙接触，用来模拟无滑动的，表面完全粗糙的摩擦接触问题，这种设置对应于摩擦系数无限大，因此用户定义的摩擦系数(MU)被忽略。 KEYOPT(12)=2，“不分离”接触，接触面和目标面一旦接触，在其后的分析中就连在一起(虽然允许有相对滑动)。 KEYOPT(12)=3，绑定接触模式，目标面和接触面一旦接触，随后就在所有方向上绑定。 KEYOPT(12)=4，不分离接触，其中的接触积分点，或初始在球形区域内，或一旦接触，就总是与目标面沿接触面的法向连在一起， 但允许滑动。调整 FKOP(见下)，可用“软弹簧”把这些区

31、域联系在一起。 KEYOPT(12)=5，绑定接触模式，其中接触积分点，或者初始在球形区域内，或者一旦接触，就总是与目标面沿接触面的法向和切向绑定在一起。 KEYOPT(12)=6，绑定接触模式，其中初始接触的接触积分点保持与目标面接触，而初始处于打开状态的接触积分点，在整个分析期间保持打开状态。这个选项与在初始接触的区域应用 CEINTF 类似。对于模拟不分离或绑定接触，用户可能需要设置 FKOP 实常数。这在接触张开时，提供一个刚度系数。如果 FKOP 为正值，则真正的接触张开刚度等于 FKOP乘以接触闭合时施加的刚度。如果 FKOP 为负值，该值作为 接触张开刚度的绝对值。缺省的 FKO

32、P 值为1。不分离或绑定接触，在接触发生张开时，产生回拉力，这个力可能不足以阻止分离。为了减小分离，定义一个较大的 FKOP 值。在有些时候，希望接触面分离，但需要在接触面之间建立联系来阻止刚体运动，在这种情况下，可以指定较小的 FKOP 值，以使接触面之间保持联系(这是软弹簧效应)。5.4.9.12 用超单元建立接触模型面一面的接触单元能模拟刚体(或一个线弹性体)和另一个有小位移的线弹性体的接触。这些线弹性体可用超单元来建模，这大大降低了进行接触迭代的自由度数。记住任何接触或目标节点都必须是超单元的主自由度。由于超单元仅仅由一组保留的节点自由度组成，它没有用来定义接触面和目标面的几何形状。因

33、此，必须在形成超单元之前在原始单元表面上定义接触面和目标面。来自超单元的信息，包括节点连结和组合刚度，但是没有材料特性和应力状态(轴对称、平面应力或平面应变)。一个限制是接触单元的材料特性设置必须与形成超单元之前的原始单元的材料特性相同。使用 KEYOPT(3)来提供2D接触分析的信息。CONTA171、CONTA172 单元的选项如下： 不使用超单元(KETOPT(3)=0)。 轴对称(KEYOPT(3)=1)。 平面应变或单位厚度的平面应力(KEYOPT(3)=2)。 需要厚度输入的平面应力(KEYOPT(3)=3)。对这种情况，使用实常数R2来指定厚度。对于3D接触分析，CONTA173

34、，CONTA174单元的KEYOPT(3)选项忽略。ANSYS 将自动检查下伏单元是否为超单元。5.4.9.13 考虑厚度影响程序能够用KEYOPT(11)来考虑壳(2D和3D)、梁(2D)的厚度。对于刚体-柔体接触，ANSYS将自动移动接触面到壳/梁的底面或顶面。对于柔体-柔体接触，ANSYS将自动移动与壳/梁单元相连的接触面和目标面。缺省时，程序不考虑单元厚度，用中面来表示梁和壳，而穿透距离从中面计算。当设置 KFTOPI(11)=1 时，则考虑梁或壳的厚度。从指定的底面或顶面来计算穿透距离。注意 -仅在使用节点位于中面的壳或梁单元时（例如，KEYOPT(11)=0的SHELL91单元），

35、用KEYOPT(11)=1来考虑厚度影响。建模时如要考虑厚度，记住偏移可能来自接触面或目标面或两者。在和KEYOPT(11)=1一起指定接触偏移(CNOF)时，CNOF从壳/梁的顶面和底面计算，而不是中面。当和SHELL181一起使用时，还考虑变形过程中的厚度变化。5.4.9.14 使用时间步长控制时间步长控制是一个自动时间步长特征，这个特征预测什么时间接触单元的状态将发生变化，或者二分当前时间步。使用KEYOPT(7)来选择下列四种行为之一来控制时间步长。KEYOPT(7)=0时不提供控制(缺省)，KEYOPT(7)=3提供最多的控制： KEYOPT(7)=0，无控制。时间步大小不受预测影响

36、。当自动时间步长激活，且允许一个很小的时间步长时，这个设置对大多数情况是合适的。 KETOPT(7)=1，如果一次迭代期间产生太大的穿透，或者接触状态急剧变化，则进行时间步长二分。 KEYOPT(7)=2，对下一个子步预测一个合理的时间增量。 KETOPT(7)=3，对下一个子步，预测一个最小的时间增量。5.4.9.15 使用死活单元选项面面接触的接触单元和目标单元允许激活或杀死，而且也跟随其下伏单元的死活状态。能够在分析的某一阶段中杀死这个单元，而在以后的阶段再重新激活它。这个特征对于模拟复杂的金属成形过程是有用的，在此过程的不同分析阶段，有多个刚性目标面需要和接触面相互作用。回弹模拟常常需

37、要在成形过程的后期移走刚性工具。这一选项不能用于“不分离”或绑定接触。5.4.10 控制刚性目标面的运动(刚体-柔体接触)刚性目标面是在其原始构形上定义的。而整个面的运动，通过控制节点上的给定位移来定义(如果没有定义控制节点，则通过目标面上的不同节点来定义)。为控制整个目标面的边界条件(和运动)，在下面任何情况下，必须使用控制节点: 目标面上作用着给定的外力。 目标面发生旋转。 目标面和其它单元相连(例如结构质量单元 MASS21 等)。 目标面的运动有平衡条件调节。控制节点的自由度代表着整个刚性面的运动，包括2D中的2个平移和1个转动自由度，或3D中的3个平移和3个转动自由度。用户可以在控制

38、节点上施加边界条件(位移、初速度)、集中载荷、转动等等。为了考虑刚体的质量，可在控制节点上定义一个质量单元。当使用控制节点时，记住对目标面有如下一些限制： 每个目标面只能有一个控制节点。 ANSYS忽略除了控制外的所有节点上的边界条件。 只有控制节点能与其它单元相连。 当定义了控制节点后，不能使用约束方程( CE )或节点耦合( CP )来控制目标面的自由度。如果要在刚性面上施加任意载荷或者约束，用户必须定义控制节点，并在控制节点上加载。如果没有使用控制节点，则只能有刚体运动。注意 控制节点可以是目标单元上的一个节点，或者是任意位置的节点，但不应该是接触单元上的节点。只有在施加了转角或力矩时，

39、控制节点的位置才重要。对于每一个控制节点，ANSYS将自动定义一个内节点及一个内部约束方程。通过内部约束方程，控制节点的转动自由度与内节点的平移自由度联系了起来。缺省时，目标单元的 KEYOPT(2)=0，ANSYS 对每个目标面检查边界条件。如果下面的条件都满足，那么程序将目标面作固定处理： 在目标面节点上没有明确定义边界条件或给定力。 目标面上的节点没有和其它单元相连。 没有在目标面上的节点上使用约束方程或节点耦合。在每个载荷步的末尾，程序将会释放内部设置的约束条件。在结果文件(Jobname.RST)和数据库(Jobname.DB)中保存的约束条件可能会由于这些改变而被修改。用户应当在重启动一个分析或用交互模式重新求解之前，仔细检查当前的约束条件是否合乎要求。如果需要，用户可以通过在目标单元定义中设置KEYOPT(2)=1，来控制目标节点的约束条件。