有限元分析及工程应用-2016第六章.ppt

《有限元分析及工程应用-2016第六章.ppt》由会员分享，可在线阅读，更多相关《有限元分析及工程应用-2016第六章.ppt（114页珍藏版）》请在三一文库上搜索。

1、有限元基本理论及应用,第六章 接触问题的有限元分析,接触问题存在两个较大的难点： 1）在求解问题之前，接触区域不确定，表面之间是否接触或分开是未知的、瞬时变化的，它由载荷、材料、边界条件和其它因素而定； 2）大多的接触问题需要计算摩擦，有多种摩擦及其模型可供挑选，但它们都是非线性的，摩擦使问题的收敛性变得困难。 接触问题基本类型：刚体柔体接触，柔体柔体接触。 在刚体柔体的接触中，接触面的一个或多个被当作刚体(与它接触的变形体相比，有大得多的刚度)，一般情况下，一种软材料和一种硬材料接触时，问题可以被假定为刚体柔体的接触，许多金属成形问题归为此类接触； 而柔体柔体的接触，是一种更普遍的类型，在这

2、种情况下，两个接触体都是变形体(有近似的刚度)。,6.1 接触边界的有限元算法,（1）直接迭代法,在用有限元位移法求解接触问题时，首先假设初始接触状态形成系统刚度矩阵，求得位移和接触力后，根据接触条件不断修改接触状态，重新形成刚度矩阵求解，反复迭代直至收敛。,每次迭代都要重新形成刚度矩阵，求解控制方程，而实际上接触问题的非线性主要反映在接触边界上，因此，通常采用静力凝聚技术，使得每次迭代只是对接触点进行， 大大提高了求解效率。,虚力法：用沿边界的虚拟等效压力来模拟接触状态，这样在每次迭代中并不重新形成刚度矩阵，所做的只是回代工作。,有限元混合法：以结点位移和接触力为未知量，并采用有限元形函数插

3、值，将接触区域的位移约束条件和接触力约束条件均反映到刚度矩阵中去，构成有限元混合法控制方程,6.1 接触边界的有限元算法,（1）直接迭代法,对弹塑性接触问题，在求解过程中接触非线性和材料非线性都需要迭代求解。 通常是利用系统刚度矩阵的变化来反映材料非线性的影响，在每次塑性修正迭代过程中都要结合对接触状态的判断进行接触迭代计算，并且，荷载增量更是受到不允许在一个增量步中出现两种非线性的限制。,6.1 接触边界的有限元算法,（2）接触约束算法,接触约束算法就是通过对接触边界约束条件的适当处理，将约束优化问题转化为无约束优化问题求解。 根据无约束优化方法的不同，可分为罚函数方法和Lagrange 乘

4、子法。,1）罚函数方法,将接触非线性问题转化为材料非线性问题。分为障碍函数法和惩罚函数法,6.1 接触边界的有限元算法,（2）接触约束算法,1）罚函数方法,障碍函数法假设接触面之间充满某虚拟物质，在未接触时其刚度趋于零，不影响物体的自由运动，在接触时其刚度变得足够大，能阻止接触物体之间的相互嵌入。 常用的间隙元等方法均属于此类，该方法处理简单，编程方便，只是在传统有限元分析中增加一种单元模式而已； 惩罚函数法对接触约束条件的处理是通过在势能泛函中增加一个惩罚势能。,惩罚因子,嵌入深度，是节点位移的函数,接触问题就等价于无约束优化问题:,6.1 接触边界的有限元算法,（2）接触约束算法,1）罚函

5、数方法,其中,由于人为假设了很大的罚因子，可能引起方程的病态。,2） Lagrange 乘子法与增广Lagrange 乘子法,约束最小化问题转化为无约束最小化问题,6.1 接触边界的有限元算法,（2）接触约束算法,2） Lagrange 乘子法与增广Lagrange 乘子法,代入,Lagrange 乘子法中接触约束条件可以精确满足,增广Lagrange 乘子法：最直接的一种方法是构造修正的势能泛函：,6.1 接触边界的有限元算法,（2）接触约束算法,2） Lagrange 乘子法与增广Lagrange 乘子法,相应的控制方程为：,考虑Lagrange 乘子的物理意义，可将其用接触对的接触应力代

6、替，通过迭代计算得到问题的正确解。 在迭代过程中，接触应力作为已知量出现，这样既吸收了罚函数方法和Lagrange 乘子法的优点，又不增加系统的求解规模，而且收敛速度也比较快。,另一种增广Lagrange 乘子法主要是为了弥补Lagrange 乘子法中控制矩阵存在零主元的弱点：,6.1 接触边界的有限元算法,（2）接触约束算法,2） Lagrange 乘子法与增广Lagrange 乘子法,解收敛于,解,系统的控制方程写为：,6.1 接触边界的有限元算法,（2）接触约束算法,3） 数学规划法,基于势能或余能原理，并利用变分不等式等现代数学方法而导出，最初该方法是针对无摩擦接触问题提出，它利用了无

7、摩擦接触问题的非穿透条件和互补条件：,经有限元离散后，无摩擦接触问题被归结为二次规划( 线性互补) 问题求解。把摩擦条件可以写成如下带导数的互补形式：,6.1 接触边界的有限元算法,（2）接触约束算法,3） 数学规划法,对摩擦条件的另一种处理方法是引进惩罚因子，然后仿照塑性力学将摩擦接触条件表示成有惩罚因子的互补形式:,有了上述摩擦接触条件的互补关系，就可以利用参变量变分原理或虚功原理建立摩擦接触问题的有限元二次规划(线性互补)模型。 对这类线性互补问题常采用Lemke 算法求解。,6.1 接触边界的有限元算法,（2）接触约束算法,3） 数学规划法,对三维摩擦接触问题，为了能利用线性互补方法求

8、解，通常以多面体棱锥近似代替Coulomb 圆锥，从而实现滑动函数的线性化，但该方法大大增加了问题的求解规模。 为了尽量减少线性化所增加的求解规模，有的学者又提出了参数二次规划迭代算法、序列线性互补方法等。 三维摩擦接触问题本质上属于非线性互补问题，由此出发可以建立非线性互补接触力法模型和非线性互补接触位移法模型。 数学规划方法在弹塑性接触问题的应用上通常用迭代法反映材料非线性特征，在每次迭代中用数学规划方法求解接触问题。 钟万勰等利用参变量变分原理将接触问题和弹塑性问题表示成具有相同形式的有限元参数二次规划问题, 很方便地实现了弹塑性接触问题的数学规划解法。,6.2 ANSYS软件中接触算法

9、,（1）显式接触边界算法,ANSYS/LS-DYNA程序中处理不同结构界面的接触碰撞和相对滑动是程序中非常重要和独特的功能，有二十多种不同的接触类型可供选择。 主要是变形体与变形体的接触、离散点与变形体的接触、变形体本身不同部分的单面接触、变形体与刚性体的接触、变形结构固连以及根据失效准则接触固连，模拟钢筋在混凝土中固连和失效滑动的一维滑动线等。,ANSYS/LS-DYNA程序处理接触碰撞面主要采用3种不同的算法，即节点约束法（仅用于固连界面）、对称罚函数法（最常用的算法）和分配参数法（仅用于滑动界面）。,不同结构可能相互接触的两个面分别称为主表面(其中的单元表面称为主片、节点称为主节点)和从

10、表面(其中的单元表面称为从片、节点称为从节点)。,6.2 ANSYS软件中接触算法,（1）显式接触边界算法,节点约束法是最早选用的接触算法。 原理：在每一时间步长修正构型之前，检查每一个没有与主表面接触的从节点，是否在当前时间步长内贯穿主表面。 如果从节点贯穿主表面，则将时间步长缩小，使那些贯穿从节点都不贯穿主表面，而其中有的刚到达主表面。 在下一个时间步长开始时，对刚到达主表面的从节点施加碰撞条件。 对所有已经与主表面接触的从节点都施加约束条件，以保持从节点与主表面接触。 此外，检查与主表面接触的从节点所属单元是否存在受拉交界面力。 如有受拉面力，则用释放条件使从节点脱离主表面。由于此算法比

11、较复杂，后来只用于固连界面。,6.2 ANSYS软件中接触算法,（1）显式接触边界算法,分配参数法在LS-DYNA程序中仅用于滑动处理 原理是：将每一个正在接触的从单元一半质量分配到被接触的主表面面积上，同时由每一个从单元的内应力确定所用在接受质量的主表面面积上的分布压力。 在完成质量和压力分配后，程序修正主表面的加速度。然后对从节点的加速度和速度施加约束，以保证从节点沿主表面运动。程序不允许从节点穿透主表面，从而避免反弹。,对称罚函数法是一种新的算法，1982年8月开始用于DYNA2D程序。 原理：每一时步先检查各从节点是否穿透主表面，没有穿透则对该从节点不做任何处理。 如果穿透，则在该从节

12、点与被穿透主表面之间引入一个较大的界面接触力，其大小与穿透深度、主片刚度成正比，称为罚函数值。,6.2 ANSYS软件中接触算法,（1）显式接触边界算法,物理意义：相当于在从节点和被穿透主表面之间放置一个弹簧，以限制从节点对主表面的穿透。 对称罚函数法是同时再对各主节点处理一遍，其算法与从节点一样。 对称罚函数法编程简单，很少激起网格沙漏效应，没有噪音，这是由于算法具有对称性，动力守恒准确，不需要碰撞和释放条件。 罚函数值大小受到稳定性限制，若计算中发生明显穿透，可以放大罚函数值或缩小时间步长来调节。 对称罚函数法：它在每一时步分别对从节点和主节点循环处理一遍，算法相同。,6.2 ANSYS软

13、件中接触算法,（1）显式接触边界算法,6.2 ANSYS软件中接触算法,（1）显式接触边界算法,6.2 ANSYS软件中接触算法,（1）显式接触边界算法,不用作任何处理,6.2 ANSYS软件中接触算法,（1）显式接触边界算法,主片的刚度因子:,6.2 ANSYS软件中接触算法,（1）显式接触边界算法,6)摩擦力的计算,界面刚度,6.2 ANSYS软件中接触算法,（1）显式接触边界算法,时间步长,衰减系数,6.2 ANSYS软件中接触算法,（1）显式接触边界算法,由库仑摩擦造成界面的剪应力，在某些情况下，超过材料承受的能力，程序采用某种限制措施，令：,7) 节点力总体坐标轴方向分量的获得。 将

14、接触力矢量和摩擦力矢量投影到总体坐标轴方向并组集到总体载荷矢量P。,对称罚函数法是将上述算法对从节点和主节点分别循环处理，如果仅对从节点循环处理，则称为“分离和摩擦滑动一次算法”。 用于金属模具压力成形问题，此时，将模具视为刚体，变形体表面定义为从表面，而刚体表面定义为主表面，它比对称算法节省机时。,6.2 ANSYS软件中接触算法,（1）显式接触边界算法,ANSYS/LS-DYNA程序对表面与表面间的接触算法又做了不少改进: 主要是为了使薄板模压成形问题计算更加精确。例如接触搜索，采用刚体近似冲压模具时有的单元长宽比很不好，搜索与从节点最靠近的主节点有时很困难，采用搜索最接近的主片位置来代替

15、搜索最接近的主节点。 由于计算接触表面位置考虑到壳单元的厚度，因而在薄板模压成形时，壳单元的厚度变化对接触表面摩擦力影响很大。 在接触计算中增加粘性接触阻尼项，以阻止薄板模压成形计算过程中垂直于接触表面的振荡。,（2）隐式接触边界算法,两个较大的难点： 1）在用户求解问题之前，用户通常不知道接触区域。 随载荷、材料、边界条件和其它因素的不同，表面之间可以接触或者分开，在很大程度上是难以预料的，并且还可能是突然变化的。,6.2 ANSYS软件中接触算法,（2）隐式接触边界算法,2）大多数的接触问题需要考虑摩擦作用，有几种摩擦定律和模型可以挑选，它们都是非线性的。 摩擦效应可能是无序的，所以摩擦使

16、问题的收敛性成为一个难点。,ANSYS隐式算法中接触问题： 刚体柔体的接触，一个或多个接触面能被当作刚体（与它接触的变形体相比，有较大的刚度）。 一般情况下，当一种软材料和一种硬材料接触时，可以假定为刚体柔体的接触（许多金属成形问题）。 柔体柔体的接触是一种更常用的类型，在这种情况下，两个接触体都是变形体(有相似的刚度)，比如螺栓法兰连接。,ANSYS接触方式：点点、点面、点点、面面、线线和线面接触，每种接触方式使用不同的接触单元集，并适用于某一特定类型的问题。,6.2 ANSYS软件中接触算法,（2）隐式接触边界算法,1）面面的接触单元,ANSYS支持刚柔体的面面接触，刚性面(“目标”面)，

17、分别用TARGE169和TARGE170单元来模拟2D和3D的“目标”面，柔性体的表面（ “接触”面），用CONTA171、CONTA172、CONTA173、CONTA174等单元来模拟。,一个目标单元和一个接触单元可称作一个“接触对”，并通过一个共享的实常数号来识别“接触对”。 为了建立一个“接触对”，必须要给目标单元和接触单元指定相同的实常数号。,面面接触非常适合于过盈装配、插入、锻造和深拉伸等问题的接触分析。,6.2 ANSYS软件中接触算法,（2）隐式接触边界算法,1）面面的接触单元,支持低阶和高阶单元。 为工程需要提供更好的接触结果，如法向压力和摩擦应力。 对目标面的形状没有限制，

18、表面可以是自然不连续或由于网格离散引起不连续。 允许流体压力渗透载荷的模拟。,面面接触单元，能模拟直线(面)和曲线(面)，通常用简单的几何形状模拟曲面，更复杂的刚体形状能使用特殊的前处理技巧来建模。,不适合于点点、点面、边缘面或3维的线线接触。仅支持一般的静态和瞬态分析、屈曲分析、谐分析、模态或随机分析、子结构分析。,6.2 ANSYS软件中接触算法,（2）隐式接触边界算法,2）点面的接触单元,点面接触单元是CONTA175，它支持大滑动、大变形以及接触面之间存在的不同网格。当单元渗透到指定目标面上的某个目标单元时接触就发生了。（模拟如扣合件的角点沿配合面滑动）。 如果接触面是由一组节点组成时

19、，也可以使用CONTA175单元来模拟面面接触，表面可也是刚体或变形体（电线插入插槽内）。 用户不需要知道接触的精确位置，接触组件之间也不需要匹配的网格。 尽管可用于模拟小滑动，但同时也允许模拟大变形和更大的相对滑动。 CONTA175单元不支持接触表面具有3维高阶单元。如果目标面出现严重不连续，这个也会出现失效。不能采用等值线来显示接触结果。,6.2 ANSYS软件中接触算法,（2）隐式接触边界算法,3）三维线线接触,线线接触单元(CONTA176)可用来模拟3维梁梁的接触，或管子在另一根管子内滑动，某些典型的应用也包括机织物和网球球拍弦。 CONTA176被附着在3维梁或管单元上，支持接触

20、表面的低阶或高阶单元，目标面可用3维线段(直线或抛物线单元)来模拟，可适用于大滑动、大位移等问题。,线面接触单元(CONTA177)可用来模拟一个3维梁或壳边缘接触实体或壳单元。 它支持接触面的低阶或高阶单元，目标面采用TRAGE170单元来模拟，这个单元也适用于大滑动、大位移等问题。,4）线面接触,6.2 ANSYS软件中接触算法,（2）隐式接触边界算法,5）点点接触单元,点点接触单元(CONTA178)主要用于模拟点点的接触行为，为了使用点点的接触单元，必须要预先知道接触位置，只能适用于接触面之间有较小相对滑动的情况(即使在几何非线性情况下)。 如果两个面上的结点一一对应，相对滑动变形又可

21、以忽略不计，两个面的挠度(转动)较小，则可用点点接触单元来求解面面的接触问题。,6.3 接触对生成与参数设置,（1）接触对的生成,接触管理器能够指定、浏览和编辑接触对，它提供了一种非常方便的方式来管理所有的接触对；接触向导能引导生成接触对，它也可从接触管理器中进入。,1）接触管理器(Contact Manager),可在开始阶段、前处理器、求解器和通用后处理器中使用。,6.3 接触对生成与参数设置,（1）接触对的生成,1）接触管理器(Contact Manager),接触向导(Contact Wizard)：允许用户采用手动方式来指定目标和接触面，可应用于2维和3维几何模型。 接触特性(Con

22、tact Properites)：对于使用的接触单元，可以通过实常数和关键字(KEYOPT)来指定接触对的性能。 删除接触对(Delete Contact Pairs)：删除接触对列表中被选的接触对。 接触对选择选项(Contact Selection Options)：指定接触单元、目标单元或两者的显示。 显示单元/结果(Plot Element/Results)：显示被选接触对的单元，如果当前处于后处理器，则可显示接触结果。 显示接触对的法向(Show Normals)：确定在显示接触对时是否显示单元的法向。,6.3 接触对生成与参数设置,（1）接触对的生成,1）接触管理器(Contact

23、 Manager),转换法向(Flip Normals)：转换被选接触对单元的法向，它仅适用于在接触对选择选项中指定的单元。 目标和接触单元转换(Switch Contact and Target)：将接触面和目标面之间相互转换，仅适用于柔性柔性的面面接触对。,列表单元/结果(List Elements/Results)：列表出被接触对的单元，如果处于通用后处理器中，则可列表出接触结果。 模型背景(Model Context)：采用半透明方式在整体模型上显示接触对，或仅显示接触对。 接触对状态(Check Contact Status)：提供被选接触对的相关信息。 接触对结果(Contact

24、Results)：显示接触结果项。,6.3 接触对生成与参数设置,（1）接触对的生成,2） 接触向导(Contact Wizard),引导用户采用手工方式创建接触对，可用于刚体柔性和柔性柔性接触，也可用于面面接触和点面接触的配置，不支持刚性目标初始化。,如果模型没有划分网格，则接触向导不可用。 如果希望创建一个刚体柔性模型，在启动向导之前，仅需要对作为柔性表面的模型划分网格； 如果想创建一个柔性柔性模型，则在启动向导之前，则必须要对作为接触表面的模型全部划分网格。,可使用线、面积、体、被选的一组节点或节点组件来生成目标和接触面，它也允许同时选择多个面来指定为目标和接触面。 在指定好目标和接触面

25、后，在创建接触对之前也可以指定接触对的属性，然后单击Create则创建了一个接触对。,6.3 接触对生成与参数设置,（2）接触分析的参数设置,1）面面接触分析接触算法的选择,面面接触单元，其接触算法有： 罚函数方法(KEYOPT(2)=1) 增进的Lagrange法(KEYOPT(2)=0) 在接触法向为Lagrange法，接触切向为罚函数法(KEYOPT(2)=3) 在接触法向和切向均为Lagrange法(KEYOPT(2)=4) 内部多点约束法(MPC)( KEYOPT(2)=2),罚函数:使用一个接触“弹簧”来建立两个接触面之间的联系，弹簧刚度就称为接触刚度 使用的实常数有：FKN和FK

26、T，如果KEYOPT(10)=2的话，则还需要有FTOLN和SLTO。,增进的Lagrange法:对罚函数进行一系列修正迭代，在平衡迭代中接触力被放大，从而使最后的渗透要比许可公差(FTOLN)要小得多。,6.3 接触对生成与参数设置,（2）接触分析的参数设置,1）面面接触分析接触算法的选择,与罚函数法相比，Lagrange法容易得到良态条件，对接触刚度的敏感性较小。 在有些分析中，增进的Lagrange法可能需要更多的迭代，特别是在变形后网格变得太扭曲时。使用Lagrange法的同时应使用实常数FTOLN（Lagrange法指定容许的最大穿透量）。,当接触关闭时，Lagrange法将施加零的

27、渗透，当出现粘性接触时，它施加一个“零的滑移”。 不需要接触刚度FKN和FKT，相反它需要颤动控制参数：FTOLN和TNOP。 将接触力作为一个附加的自由度添加到模型上，从而需要附加迭代来建立平衡条件，与增进的Lagrange法相比，它的计算时间要长一些。同时要关闭稀疏解器。,6.3 接触对生成与参数设置,（2）接触分析的参数设置,1）面面接触分析接触算法的选择,在接触法向使用Lagrange法而在其摩擦面采用罚函数，该法对于粘性接触条件将施加一个零的渗透，并允许小量的滑动，除了最大许可弹性滑动参数SLTO外，还需要颤动控制参数：FTOLN和TNOP。 内部多点约束算法与绑定接触(KEYOPT

28、(12)=5或6)和不分离接触(KEYOPT(12)=4)一起使用，可用来模拟接触装配和运动学约束等。,1）面面接触分析决定接触刚度,对于增进的Lagrange法和罚函数法需要定义法向和切向接触刚度。 在接触与目标面之间的渗透大小取决于法向刚度，而粘性接触中的滑动量的大小取决于切向刚度。,6.3 接触对生成与参数设置,（2）接触分析的参数设置,1）面面接触分析决定接触刚度,大的刚度值会减少渗透/滑动量，但会导致一个病态条件的整体刚度矩阵并引起收敛困难； 小的刚度值会导致大的渗透/滑动量并产生一个不精确的解。 理想状态是定义一个足够大的刚度以使渗透/滑动量达到可接受的小量，而一个足够小量的刚度会

29、使问题在收敛方面具有良好的状态。,ANSYS软件为接触刚度(FKN，FKT)、许可渗透(FTOLN)和许可滑移(SLTO)提供了一个缺省值，在大多数情况下，不需要定义接触刚度，同时建议使用“KEYOPT(10)=2”允许系统自动更新接触刚度。,对于某些接触问题，可选择使用实常数FKN(缺省值为1.0)来定义法向接触刚度系数，该系数介于0.011.0之间，其缺省值对于许多变形来说是合适的，如果指定了弯曲变形，则建议使用一个较小的值。,6.3 接触对生成与参数设置,（2）接触分析的参数设置,1）面面接触分析决定接触刚度,缺省的法向接触刚度会受到模型中的材料属性、单元大小和总自由度数的影响，在求解过

30、程许多因素也许会施加到实际的法向接触刚度。,实常数FTOLN(缺省值为0.1)与增进的Lagrange法一起使用，它是一个公差因子将被施加在表面的法向，它的范围要小于1.0(通常小于0.2)，取决于其下实体单元、壳单元或梁单元的厚度。该因子用来确定渗透兼容性是否满足，如果渗透是在一个许可的公差范围(FTOLN乘以其下单元的厚度)内，则接触的协调性是满足的。,如果ANSYS系统检测到任意渗透量大于这个许可范围，即使残余力和位移增量已达到收敛，则整体求解还是考虑为不收敛。通过对FTOLN指定一个负值来定义一个绝对许可渗透量。同时FTOLN也可与Lagrange法一起使用，但此时它将是一个颤动控制参

31、数。,6.3 接触对生成与参数设置,（2）接触分析的参数设置,1）面面接触分析决定接触刚度,系统会自动生成一个缺省的切向接触刚度，它与MU(摩擦系数)和法向刚度FKN成比例。 缺省的切向刚度是：FKT=1.0，正值为FKT提供一个比例因子，负值则指定绝对切向刚度值。 对于KEYOPT(10)=2或KEYOPT(2)=3的状态，系统会根据当前法向接触压力、最大许可弹性滑移(SLTO)自动更新切向接触刚度(KT=FKT*MU*PRES /SLTO)。 当FKT在每次迭代自由更新时，SLTO则控制最大滑移距离。一个大的值会提高收敛但会影响精度。,试运行的步骤：,开始时取一个较低的接触刚度值，一般来说

32、，低估其值比高估该值要好一些，因为较低的接触刚度所导致的渗透问题要比过高的接触刚度导致的收敛性困难要容易解决；,6.3 接触对生成与参数设置,（2）接触分析的参数设置,1）面面接触分析决定接触刚度,用最终载荷的一部分来运行分析(刚好使接触完全建立)。 检查每个子步中的渗透和平衡迭代的次数。 如果整体不收敛是由于渗透量过大引起，而不是由残余力和位移增量引起，这时FKN的值太大或者FTOL的值过小； 如果整体收敛需要许多平衡迭代步才能满足残余力和位移的收敛，则表明FKN和FTOL的值高估了。 根据需要调整FKN、FKT、FTOLN或SLOT的值，开始一个完整的分析过程，如果在整体平衡迭代中，用户能

33、够支配渗透量，这时可以增大FTOLN的值来获得更大的许可渗透量，或者增大FKN的值。,1）面面接触分析选择摩擦模型,粘和状态(stick）：在基本的库仑摩擦模型中，两个接触面在开始相互滑动之前，在界面上会有某一大小的剪应力产生。,6.3 接触对生成与参数设置,（2）接触分析的参数设置,1）面面接触分析选择摩擦模型,粘聚力,滑动状态(sliding): 一旦剪应力超过此值后，两个表面将开始相互滑动。 粘和/滑动计算决定什么时候一个点从粘和状态到滑动状态，或从滑动状态到粘和状态。,对于另一种系统支持的摩擦模型，可以使用USERFRIC子程序来定义用户需要的摩擦模型。,对于无摩擦、粗糙和绑定接触，接

34、触单元的刚度矩阵是对称的。,大多数摩擦接触问题采用对称求解器来求解，但如果摩擦应力对整体位移有很大的影响时会导致低的收敛率，可通过“NROPT，UNSYM”来设置非对称求解器以改善收敛精度。,6.3 接触对生成与参数设置,（2）接触分析的参数设置,1）面面接触分析选择摩擦模型,实常数TAUMAX是具有应力单位的最大接触摩擦。当摩擦应力达到这个值时，不管法向接触压力值的大小，滑动都将会发生。,另一个使用实常数是凝聚力：COHE(缺省值为0)，具有应力单位，即使具有法向压力为零，它也提供一个滑动阻力,摩擦系数依赖于接触面的相对滑动速度，通常静摩擦系数要高于动摩擦系数。ANSYS中提供了如下的指数衰

35、减摩擦模型：,6.3 接触对生成与参数设置,（2）接触分析的参数设置,1）面面接触分析选择摩擦模型,如果不指定衰减系数，且FACT1.0，当接触进入滑动摩擦状态时，摩擦系数会从静摩擦系数突变到动摩擦系数，模型会导致收敛困难。,6.3 接触对生成与参数设置,（2）接触分析的参数设置,1）面面接触分析接触检查位置的选择,接触检查点位于接触单元的积分点上。在积分点上，接触单元不能穿透进入目标面，目标面能穿透进入接触面。,ANSYS面面接触单元计算使用高斯积分点(缺省值)，比节点检测方法一般来说会提供更精确的结果，节点检测法是采用节点本身作为积分点。 接触单元CONTA175、CONTA176和CON

36、TA177总是采用节点检测法。,节点检测法需要一个光滑的接触表面(KEYOPT(4)=1)或一个光滑的目标面(KEYOPT(4) =2)，这是相当耗费时间的，仅用于处理角点、点面或边缘表面的接触问题。,6.3 接触对生成与参数设置,（2）接触分析的参数设置,1）面面接触分析接触检查位置的选择,使用节点作为接触检查点可能会导致其它收敛性问题，如“节点滑脱”(节点滑出目标面的边界)。,为了防止滑脱，当实常数TOLS的缺省值不能避免这个问题时，应设置TOLS来延长目标面，对于大多数的点面接触问题，建议使用CONTA175。,当接触节点有可能位于目标的边缘时，则TOLS是非常有用的。TOLS的单位是百

37、分数，即：1.0表示增加接触边缘长度的1.0%。 一个非常小的值也许就可以防止滑脱的发生，对于小变形问题其缺省值为10，对于大变形问题(“NLGEOM, on”)2.0。,6.3 接触对生成与参数设置,（2）接触分析的参数设置,1）面面接触分析调整初始接触条件,在动态分析中，刚体运动一般不会引起问题，然而在静力分析中，当物体没有足够的约束时会产生刚体运动。在静态分析中，出现“Zero or negative”的警告信息或不实际的、非常大的位移将表明为没有约束的运动。,刚体外形常常是复杂的，很难决定第一个接触点发生在哪儿。 既使实体模型是在初始接触状态，在网格划分后基于数值舍入误差，两个面的单元

38、网格之间也可能会产生小的缝隙。 接触单元的积分点和目标单元之间可能有小的缝隙。,接触和目标面之间会出现过多的初始渗透。在此情况下，接触单元也许高估了接触力，从而导致不收敛或偏离了接触。,目的：消除由于生成网格造成的数值舍入误差而引起的小缝隙或渗透,6.3 接触对生成与参数设置,（2）接触分析的参数设置,1）面面接触分析调整初始接触条件,实常数CNOF：指定接触表面偏移量。指定正值将使整体接触面接近目标面，负值将使接触面偏离目标面。 如果同时输入CNOF和PINB时，要确定PINB的值要大于CNOF的值，否则CNOF将会被忽略； 如果输入了CNOF，而PINB采用缺省值，则PINB将会被调节到大

39、于CONF的值。 通过设置KEYOPT(5)的值，ANSYS系统会自动调整CNOF的值。 KEYOPT(5)=1：关闭间隙。 KEYOPT(5)=2：减少初始渗透。 KEYOPT(5)=3：既关闭间隙又减少初始渗透。,6.3 接触对生成与参数设置,（2）接触分析的参数设置,1）面面接触分析调整初始接触条件,实常数ICONT：指定初始接触环。初始接触环是指沿着目标面的“调整环”深度。 输入正值表示相对于变形体单元厚度的比例因子，负值表示为接触环的真正值。 当KEYOPT(5)=0、1、2或3时，其缺省值为0；若KEYOPT(5)=4，系统将根据几何尺寸提供较小且有意义的值，并给出一个该值已指定的

40、警告信息。,任何落在“调整环”领域内的接触检查点会被自动移到目标面上。,6.3 接触对生成与参数设置,（2）接触分析的参数设置,1）面面接触分析调整初始接触条件,使用CNOF和ICONT的区别是：前者用CNOF的距离改变整个接触表面；后者移动落在ICONT调节环内所有初始开着的接触点到目标面。,实常数PMIN和PMAX：指定初始容许的渗透范围。当指定PMAX或PMIN后，在开始分析时，系统会将目标面移到一个初始接触的状态。 如果初始渗透大于PMAX，系统会调整目标面以减少渗透，如果初始渗透小于PMIN，则系统调节目标面以保证初始接触，接触状态的初始调节仅仅通过平移来实现。,对给定载荷或给定位移

41、的刚性目标面将会执行初始接触状态的初始调节。 当目标面上的节点有给定的零位移时，使用PMAX和PMIN的初始调节将不会被执行。,6.3 接触对生成与参数设置,（2）接触分析的参数设置,1）面面接触分析调整初始接触条件,ANSYS程序独立地处理目标面上节点的自由度，如：若指定自由度UX值为“0”，那沿着X方向就不会有初始调查，然而，在Y和Z方向仍然会激活PMAX和PMIN选项。,初始状态调整是一个迭代过程，程序最多进行20次迭代，如果目标面不能进入可接受的渗透范围，系统则会使用未调节的几何模型，并会给出一个警告信息，这时可能需要调整初始几何模型。,刚性目标面在X方向的任何运动都不会引起初始接触。

42、,6.3 接触对生成与参数设置,（2）接触分析的参数设置,1）面面接触分析调整初始接触条件,设置KEYOPT(9)=1来消除初始渗透,过盈装配问题，期望有过度的渗透，为了缓解收敛性困难，在第一个载荷步中设置KEYOPT(9)=2来使过度渗透渐进到0,在第一个载步中不要给定其它任何载荷，也就是说要保证载荷是渐进的(KBC，0)。 在开始分析时，程序会给出每个目标面的初始接触状态的输出信息，(在输出窗口或输出文件中)，这个信息有助于决定每个目标面的最大渗透或最小间隙。,6.3 接触对生成与参数设置,（2）接触分析的参数设置,1）面面接触分析决定接触状态和Pinball区域,接触单元相对于目标面的运

43、动和位置决定了接触单元的状态；系统检测每个接触单元并给出一种状态 STAT=0 未合的且未靠近接触 STAT=1 未合的近区接触 STAT=2 滑动接触 STAT=3 粘合接触,当它的积分点位于与相关目标面系统计算(用户指定)的距离时，接触单元就被考虑为近区域接触，这个距离就被认定为pinball区域，且为关于Gauss点的圆(二维)或球(三维)。,使用实常数PINB来为pinball指定一个比例因子(正值)或真实值(负值)，可以为PINB指定任何值。,在缺省状态，并假定大变形效应打开时，系统将pinball区域定义为一个以r为半径的圆(对2-D问题)或球(对3-D问题),6.3 接触对生成与参数设置,（2）接触分析的参数设置,1）面面接触分析决定接触状态和Pinball区域,如果大变形效应未打开，pinball区域是大变形效应的一半。如果输入了实常数CNOF，且缺省的PINB值要小于CNOF的绝对值，那么缺省的PINB值将被设置为1.1*CNOF的绝对值。,寻找接触的计算时间依赖于pinball区域的大小