Ansys接触分析

1、- 1 -Ansys 接触分析接触问题是一种高度非线性行为， 需要较大的计算资源， 为了进行实为有效的计算， 理 解问题的特性和建立合理的模型是很重要的。接触问题存在两个较大的难点： 其一，在你求解问题之前，你不知道接触区域，表面之 间是接触或分开是未知的， 突然变化的， 这随载荷、 材料、边界条件和其它因素而定； 其二， 大多的接触问题需要计算摩擦， 有几种摩擦和模型供你挑选， 它们都是非线性的， 摩擦使问 题的收敛性变得困难。一般的接触分类 接触问题分为两种基本类型： 刚体柔体的接触， 和柔体柔体的接触， 在刚体柔体 的接触问题中，接触面的一个或多个被当作刚体， （与它接触的变形体相比，有

2、大得多的刚 度），一般情况下， 一种软材料和一种硬材料接触时， 问题可以被假定为刚体柔体的接触， 许多金属成形问题归为此类接触，另一类，柔体柔体的接触，是一种更普遍的类型， 在这 种情况下，两个接触体都是变形体（有近似的刚度） 。ANSY漲触能力ANSYS 支持三种接触方式：点点，点面，和面面，每种接触方式使用的接触单元 适用于某类问题。为了给接触问题建模， 首先必须认识到模型中的哪些部分可能会相互接触， 如果相互作 用的其中之一是一点， 模型的对应组元是一个结点。 如果相互作用的其中之一是一个面， 模 型的对应组元是单元， 例如梁单元， 壳单元或实体单元， 有限元模型通过指定的接触单元来 识

3、别可能的接触匹对，接触单元是覆盖在分析模型接触面之上的一层单元，至于ANSYS使用的接触单元和使用它们的过程，下面分类详述。点点接触单元点点接触单元主要用于模拟点点的接触行为， 为了使用点点的接触单元， 你需要 预先知道接触位置， 这类接触问题只能适用于接触面之间有较小相对滑动的情况 （即使在几 何非线性情况下）如果两个面上的结点一一对应， 相对滑动又以忽略不计， 两个面挠度 （转动）保持小量， 那么可以用点点的接触单元来求解面面的接触问题， 过盈装配问题是一个用点点的接 触单元来模拟面面的接触问题的典型例子。点面接触单元 点面接触单元主要用于给点面的接触行为建模，例如两根梁的相互接触。 如果

4、通过一组结点来定义接触面， 生成多个单元， 那么可以通过点面的接触单元来模 拟面面的接触问题， 面即可以是刚性体也可以是柔性体，这类接触问题的一个典型例子是插头到插座里。使用这类接触单元， 不需要预先知道确切的接触位置，接触面之间也不需要保持一致的网格，并且允许有大的变形和大的相对滑动。Contact48 和 Contact49 都是点面的接触单元， Contact26 用来模拟柔性点刚性面 的接触， 对有不连续的刚性面的问题， 不推荐采用 Contact26 因为可能导致接触的丢失， 在 这种情况下， Contact48 通过使用伪单元算法能提供较好的建模能力。面面的接触单元ANSYS 支持

5、刚体柔体的面面的接触单元， 刚性面被当作 “目标”面，分别用 Targe169 和 Targe170 来模拟 2D 和 3D 的“目标”面，柔性体的表面被当作“接触”面，用 Conta171,Conta172,Conta173,Conta174 来模拟。一个目标单元和一个接单元叫作一个“接 触对”程序通过一个共享的实常号来识别“接触对” ，为了建立一个“接触对”给目标单元 和接触单元指定相同的实常的号。- 2 -与点面接触单元相比，面面接触单元有好几项优点，支持低阶和高阶单元支持有大滑动和有摩擦的大变形，协调刚度阵计算，单元提供不对称刚度阵的选项。提供为工程目的采用的更好的接触结果，例如法向压

6、力和摩擦应力。 没有刚体表面形状的限制，刚体表面的光滑性不是必须的，允许有自然的或网格离 散引起的表面不连续。与点面接触单元比，需要较少的接触单元，因而造成需要较小的磁盘空间和CPU时间。允许多种建模控制，例如：绑定接触渐变初始穿透目标面自动移动到补始接触平移接触面（老虎梁和单元的厚度）支持死活单元使用这些单元，能模拟直线（面）和曲线（面） ，通常用简单的几何形状例如圆、抛物 线、球、圆锥、圆柱采模拟曲面，更复杂的刚体形状能使用特殊的前处理技巧来建模。 执行接触分析不同的接触分析类型有不同的过程，下面分别讨论 面面的接触分析我们可以使用面面的接触单元来模拟刚体柔体和柔体柔体的面面接触。对大多数

7、接触问题 ， 接 触 向 导 为 建 立 接 触 对 提 供 了 一 种 简 单 的 方 法 （ 通 过选 择菜 单路 径 PreprocessorCreatContact pairContact Wizard 进入接触向导 ） 。接触向导将会引导我 们方便地创建一个接触对。接触向导对话框的帮助系统将解释对话框及其选择的用法。如果模型还没有被划分网格，接触向导便会保持灰色（不起作用） 。如果希望建立一个 柔体的接触模型， 那么在进入接触向导之前应给可能成为接触面的所有部分划分网格。如果希望建立刚体柔体的接触模型， 那么在进入接触向导之前， 只需给可能成为柔体接触面的部 分划分网格。在涉及到两个

8、边界的接触问题中， 很自然把一个边界作为 “目标”面而把另一个作为 “接 触”面，对刚体柔体的接触， “目标”面总是刚性的， “接触”面总是柔性面，这两个面合起 来叫作“接触对”使用 Targe169 和 Conta171 或 Conta172 来定义 2-D 接触对，使用 Targe170 和 Conta173 或 Conta174 来定义 3-D 接触对，程序通过相同的实常收号来识别“接触对” 。 面面接触分析的步骤：执行一个典型的面面接触分析的基本步骤列示如下：1 建立模型，并划分网格2 识别接触对3 指定接触面和目标面4 定义刚性目标面5 定义柔性接触面6 设置单元关键字和实常数7 定

9、义控制刚性目标面的运动8 给定必须的边界条件9 定义求解选项和载荷步10求解接触问题11查看结果步骤一：建立模型，并划分网格- 3 -在这一步中， 你需要建立代表接触体几何形状的实体模型。 与其它分析过程一样， 设置 单元类型，实常的，材料特性。用恰当的单元类型给接触体划分网格。命令： AMESHVMESHGUI ： Main MenuPreprocessormeshMapped3 or4 SidedMain MenuPreprocessormeshmapped4 or 6 sided步骤二：识别接触对必须认识到，模型在变形期间哪些地方可能发生接触，一旦已经识别出潜在的接触面， 你应该通过目标

10、单元和接触单元来定义它们， 目标和接触单元跟踪变形阶段的运动， 构成一 个接触对的目标单元和接触单元通过共享的实常数号联系起来。接触环(区域)可以任意定义，然而为了更有效的进行计算(主要指CPU时间)你可能想定义更小的局部化的接触环， 但能保证它足以描述所需要的接触行为， 不同的接触对必 须通过不同的实常数号来定义(即使实常数号没有变化) 。由于几何模型和潜在变形的多样性，有时候一个接触面的同一区域可能和多个目标面 产生接触关系。在这种情况下，应该定义多个接触对(使用多组覆盖层接触单元)。每个接触对有不同的实常数号。步骤三：指定接触面和目标面 接触单元被约束住，不能穿透目标面，然而目标单元可以

11、穿透接触面，对刚体柔体的接 触，目标面的指定是显而易见的， 对柔体柔体的接触问题， 目标面接触面的不同选择可能 引起穿透量的不同，从而影响求解精度，当指定目标面和接触面时应考虑以下准则。(1)、当凸面和平面或凹面接触时，应指定平面或凹面作为目标面。(2)、如果一个面上的网格较细， 而另一个面上的网格较粗， 应指定细网格所在面作为接触 面，粗网格所在面作为目标面。(3)、如果两个面的刚度不同时，应当较硬的面为目标面，教软的面作为接触面。(4)、如果一个面上的基础单元是高阶单元， 而另一个面上的基础单元是低阶单元， 应将基 础单元为高阶单元的面作为接触单元。(5)、如果两个面的大小明显不同，应将大

12、面作为目标面。 步骤四：定义刚性目标面刚性目标面可能是 2D的或3 D的。在2D情况下，刚性目标面的形状可以通过一一 系列直线、圆弧和抛物线来描述，所有这些都可以用TAPGE169来表示。另外，可以使用它们的任意组合来描述复杂的目标面。在3D情况下，目标面的形状可以通过三角面，圆柱面，圆锥面和球面来推述，所有 这些都可以用 TAPGE17C来表示，对于一个复杂的，任意形状的目标面，应该使用三角面来 给它建模。控制结点( Pilot )刚性目标面可能会和“ pilot 结点”联系起来，它实际上是一个只有一个结点的单元， 通过这个结点的运动可以控制整个目标面的运动， 因此可以把 pilot 结点作

13、为刚性目标的控 制器。整个目标面的受力和转动情况可以通过 pilot 结点表示出来， “pilot 结点”可能是 目标单元中的一个结点，也可能是一个任意位置的结点，只有当需要转动或力矩载荷时，“pilot 结点”的位置才是重要的，如果你定义了pilot 结点” ANSYS程序只在pilot结点”上检查边界条件，而忽略其它结点上的任何约束。对于圆、圆柱、圆锥、和球的基本图段，ANSYS总是使用条一个结点作为 “pilot结点”。基本图元 能够使用基本几形状来模拟目标面， 例如：圆、圆柱、圆锥、球。直线、抛物线、 弧线、 和三角形不被允许、 虽然你不能把这些基本原型彼此合在一起， 或者是把它们和其

14、它的目标 - 4 -形状合在一起以便形成一个同一实常数号的复杂目标面。 但可以给每个基本图元指定它自己 的实常数号。单元类型和实常数在生成目标单元之前，首先必须定义单元类型(TARG169或TARG170。命令： ETGUI: main menupreprocessorElement Type Add/Edit/Delete 随后必须设置目标单元的实常数。命令 :RealGUI: main menupreprocessorreal constants对TARGE16刖TARGE17C仅需设置实常数 R1和R2,而只有在使用直接生成法建立目标 单元时，才需要从为指定实常数R1、R2,另外除了直接

15、生成法，也可以使用ANSYS网格划分工具生成目标单元，下面解释这两种方法。使用直接生成法建立刚性目标单元 为了直接生成目标单元，使用下面的命令和菜单路径。命令： TSHAPGUI: main menupreprocessormodeling-createElementsElem Attributes 随后指定单元形状，可能的形状有：straight line (2D)parabola (2-D)clockwise arc(2-D)counterclockwise arc (2-D)circle(2-D)Triangle (3-D)Cylinder (3-D)Cone (3-D)Sphere (

16、3-D)Pilot node (2-D 和 3-D)一旦指定目标单元形状， 所有以后生成的单元都将保持这个形状， 除非指定另外一种形 状。然后就可以使用标准的 ANSYS直接生成技术生成结点和单元。命令： NEGUI ：main menupreprocessor modeling- create nodes main menupreprocessor modeling- createElements在建立单元之后，你可以通过列示单元来验证单元形状命令： ELISTGUI ：utility menulistElementsNodes+Attributes使用ANSYS网格划分工具生成刚性目标单元

17、你也可以使用标准的 ANSYS网格划分功能让程序自动地生成目标单元，ANSYS程序将会以实体模型为基础生成合适的目标单元形状而忽略TSHAP命令的选项。为了生成一个“ PILOT结点”使用下面的命令或GUI路径：命令： KmeshGUI ： main menupreprocessormeshing-meshkeypoints注意：KMESH总是生成“ PILOT结点”为了生成一个 2D 目标单元，使用下面的命令和GUI 路径：ANSYS 在每条直线上生成一条单一的线，在样条曲线上生成抛物线部分，在每条圆弧和 倒角上生成圆弧部分，如果所有的圆弧形成一个封闭的圆，ANSYS生成一个单一的圆。- 5

18、 -命令： LMESHGUI ： main menupreprocessormesling-meshlines为了生成3 D的目标单元，使用下面的命令或GUI路径。如果实体模型的表面部分形成了一个完整的球，圆柱或圆锥，那么ANSY皐序自动生成一个基本的3 D目标单元，因为生成较少的单元，从而使你分析计算更有效率，对任意形 状的表面，应该使用Amesh命令来生成目标单元，在这种情况下，网格形状的质量不是重要的，而目标单元的形状是否能完成好的模拟刚性面的表面几何形状显得更重要。命令： AMESHGUI ： main menupreprocessor-meshing-meshAreaANSYS 在所

19、有可能的面上推荐使用三角形的映射网格划分，如果在表面的边界上没有曲率，则在网格划分时，指定那条边界分为一分， 下面的命令或 GUI路径将尽可能的生成一个 映射网格（如果不能进行映射，它将生成自由网格）命令： MSHKF，Y 2GUI ： main menupreprocessor-meshling-mesh-Ares-Target Surf 建模和网格划分的注意点：一个目标面可能由两个或多个间断的区域组成， 你应该尽可能地通过定义多个目标面来 使接触区域局部化 （每个目标面有一个不同的实常数号） 刚性目标面上由的离散能足够指述 出目标面的形状， 过粗的网格离散可能导致收敛问题。 如果刚性面有一

20、个尖的凸角， 求解大 的滑动问题时很难获得收敛结果， 为了避免这些建模问题， 在实体模型上， 使用线或面的倒 角来使尖角光滑化，或者在曲率突然变化的区域使用更细的网格。注意：不能使用镜面对称技术（ARSYSM LSYMM来映射圆、圆柱、圆锥或球面到对称平面的另一边，因为每个实常数的设置不能同时赋给多个图元。检验目标面的接触方向。目标面的结点号顺序是重要的，因为它定义了接触主向，对2D接触问题，当沿着目标线从第一个结点移向第二个结点时，变形体的接触单元必须位于目标面的右边。对3D接触问题，目标三角形单元号应该使刚性面的外法线方向指向接触面，外法线 通过右手原则来定义为了检查法线方向，显示单元坐标

21、系命令： /PSYMS， ESYS， 1GUI ： Utility menuplotctrlssymbols 如果单元法向不指向接触面，选择单元反转表面的法向的方向。命令： ESURF， REVEGUI ： main menupreprocossorcreateElementon free surf 步骤 4：定义柔性体的接触面为了定义柔性体的接触面，必须使用接触单元CONFA171或CONFA172（对2 D）或CONTA17或CONTA174对3 D）来定义表面程序通过组成变形体表面的接触单元来定义接触表面， 接触单元与下面覆盖的变形体单 元有同样的几何特性，接触单元与下面覆盖的变形体单元

22、必须处于同一阶次（低阶或高阶） 下面的变形体单元可能是实体单元、 壳单元、 梁单元或超单元， 接触面可能壳或梁单元任何 一边。与目标面单元一样， 你必须定义接触面的单元类型， 然后选择正确的实常数号 （实常数 号必须与它对应目标的实常数号相同）最后生成接触单元。单元类型： 下面简单描述四种类型的接触单元CONTA171 :这是一个2 D， 2个结点的低阶线单元，可能位于2 D实体，壳或梁单兀的表面CONTA172 :这是一个 2 D 的，3结点的高阶抛物线形单元，可能位于有中结点的2D实体或梁单元的表面- 6 -CONTA173:这是一个 3 D 的，4 结点的低阶四边形单元可能位于3D实体或

23、壳单兀的表面，它可能退化成一个结点的三角形单元。CONTA174 :这是一个3 D, 8结点的高阶四边形单元，可能位于有中结点的3-D实体或壳单元的表面，它可能退化成 6 结点的三角形单元。 不能在高阶柔性体单元的表面上分成低阶接触单元， 反之也不行， 不能在高阶接触单元 上消去中结点。命令: ETGUI : main menupreprocessorElement typeAdd/Edit/Delete 实常数和材料特性 在定义了单元类型之后， 需要选择正确的实常数的设置， 每个接触对的接触面和目标面 必须有相同的实常数号，而每个接触对必须有它自己不同的实常数号。ANSYS使用下面柔性体单元

24、的材料特性来计算一个合适的接触（或罚）刚度，如果下面 的单元是一个超单元。接触单元的材料的设置必须与超单元形成时的原始结构单元相同。生成接触单元 。我们既可以通过直接生成法生成接触单元， 也可以在柔性体单元的外表面上自动生成接 触单元，我们推荐采用自动生成法，这种方法更为简单和可靠。可以通过下面三个步骤来自动生成接触单元1 、选择结点 选择已划分网格的柔性体表面的结果，如果确定某一部分结点永远不会接触到目标面， 你可以忽略它以便减少计算时间，然而，你必须保证设有漏掉可能会接触到目标面的结点。命令: NSELGUI : main menupreprocessorcreateElementonfr

25、ee surf2、生成接触单元命令: ESURFGUI : main menupreprocessorcreateElementon free surf 如果接触单元是附在已用实体单元划分网格的面或体上，程序会自动决定接触计算所需 的外法向， 如果下面的单元是梁或壳单元， 则必须指明哪个表面 （ 上表面或下表面 ）是接触面 命令: ESURF， TOP OR BOTIOMGUI : main menupreprocessorcreateElementon free surf 使用上表面生成接触单元， 则它们的外法向与梁或壳单元的法向相同， 使用下表面生成 接触单元，则它们的外法向与梁或壳单元的

26、法向相反，如果下面的单元是实体单元，则TOP或BOTTO选项不起作用3 、检查接触单元外法线的方向，当程序进行是否接触的检查时，接触面的外法线方向 是重要的，对3D单元，按结点程序号以右手定则来决定单元的外法向，接面的外法向应- 7 -KEYOP（I 4）KEYOP（T 6）KEYOP（T 7）KEYOP（T 9）KEYOP（T 12）该指向目标面， 否则， 在开始分析计算时， 程序可能会认为有面的过度穿透而很难找到初始 解。在此情况下，程序一般会立即停止执行，你可以检查单元外法线方向是否正确。命令： /PSYMBGUI ： Utility menuplotctrlssymbols 当发现单元

27、的外法线方向不正确时，必须通过反转不正确单元的结点号来改变它们。 命令： ESURF， REVEGUI ： main menupreprocessorCreateElements on free surf 步骤六：设置实常数和单元关键字 程序使用九个实常数和好几个单元关键字来控制面面接触单元的接触行为。 实常数9 个实常数中，两个（R1和R2）用采定义目标面单元的几何形状，乘下的7个用来控制接触行为。R1 和R2定义目标单元几何形状FKN 定义法向接触刚度因子FTOLN 定义最大的穿透范围ICONT 定义初始靠近因子PINB 定义“ Pinball 区域PMIN和PMAX定义初始穿透的容许范围

28、TAUMAR 指定最大的接触摩擦CNOF给接触面指定一个正的或负的偏移值FKOP指定接触发生时所给的刚度因子 命令： RGUI ： main menu preprocessorreal constant对实常数 FKN FTOLN ICONT，PINB, PMAX PMIN,和FKOP既可以定义一个正值也可 以定义一个负值， 程序将正值作为比例因子， 将负值作为真实值， 程序将下面基本变形体单 元的厚度作为ICON, FTOLN PINB, PMAX和 PMIN的参考值，例如对 ICON, 0.1表明初始间 隙因子是 0.1* 下面基本变形体单元的厚度。然而， -0.1 表明真实缝隙是 0.1

29、 ，如果下面的 基本变形体单元是超单元，则将接触单元的最小长度作为厚度。单元关键字每种接触单元都有好几个关键字， 对大多的接触问题缺省的关键字是合适的， 而在某些 情况下，可能需要改变缺省值，来控制接触行为。接触算法（罚函数 +拉格郎日或罚函数）（KEYOP （T 2）出现超单元时的应力状态（ DEYOP（T 3）接触方位点的位置 刚度矩阵的选择 时间步长控制 初始穿透影响 接触表面情况 命令： KEYOPTETGUI ： main menupreprocessorElemant TypeAdd/Edit/Delete 选择接触算法：对面面的接触单元， 程序可以使用扩增的拉格朗日算法或罚函数方

30、法， 通过使用单元 关键字 KETOP（T 2）来指定。扩展的拉格朗日算法是为了找到精确的拉格朗日乘子而对罚函数修正项进行反复迭代， 与罚函数的方法相比，拉格朗日方法不易引起病态条件，对接触刚度的灵敏度较小，然而， 在有些分析中， 扩增的拉格朗日方法可能需要更多的迭代， 特别是在变形后网格变得太扭曲 时。使用拉格朗日算法的同时应使用实常数 FTOLNFTOLN 为拉格朗日算法指定容许的最大穿透，如果程序发现穿透大于此值时，即使不平 - 8 -衡力和位移增量已经满足了收敛准则，总的求解仍被当作不收敛处理，FTLON的缺省值为0.1 ,你可以改变这个值，但要注意如果此值太小可能会造成太多的迭代次数

31、或者不收敛。决定接触刚度所有的接触问题都需要定义接触刚度， 两个表面之间穿透量的大小取决了接触刚度， 过 大的接触刚度可能会引起总刚矩阵的病态， 而造成收敛困难， 一般来说， 应该选取足够大的 接触刚度以保证接触穿透小到可以接受， 但同时又应该让接触刚度足够小以使不会引起总刚 矩阵的病态问题而保证收敛性。程序会根据变形体单元的材料特性来估计一个缺省的接触刚度值，你能够用实常数 FKN来为接触刚度指定一个比例因子或指定一个真正的值，比例因子一般在 0.01 和 10之间， 当避免过多的迭代次数时，应该尽量使穿透到达极小值。为了取得一个较好的接触刚度值，又可需要一些经验，你可以按下面的步骤过行。1

32、、开始时取一个较低的值，低估些值要比高估些值好因为由一个较低的接触刚度导 致的穿透问题要比过高的接触刚度导致的收敛性困难，要容易解决。2、对前几个子步进行计算3、检查穿透量和每一子步中的平衡迭代次数，如果总体收敛困难是由过大的穿透引起的（而不是由不平衡力和位移增量引起的），那么可能低估了 FKN的值或者是将FTOLN的值取得大小，如果总体的收敛困难是由于不平衡力和位移增量达到收敛值 需要过多的迭代次数，而不是由于过大的穿透量导致的，那么FKN的值可能被高估。4、按需要调查FKN或FTOLN的值，重新分析。选择摩擦类型。在基本的库仑摩擦模型中， 两个接触面在开始相互滑动之前， 在它们的界面上会有

33、达到 某一大小的剪应力产生， 这种状态则作粘合状态 （ stick ）。库仑摩擦模型定义了一个等效剪 应力。一旦剪应力超过此值后，两个表面之间将开始相互滑动，这种状态，叫作滑动状态（ Sliding ）粘合 滑动计算决定什么时候一个点从粘合状态到滑动状态或从滑动状态变到粘 合状态，摩擦系数可以是任一非负值。程序缺省值为表面之间无摩擦，对rough 或 bonded接触（KEYOPT（2）=1 （或3），程序将不管给定的 MU直而认为摩擦阻力无限大。程序提供了一个不管接触压力的故而人为指定最大等效剪应力的选项，如果等效剪应 力达到此值时，滑动发生。见图41 ，为了指定接触界面上最大许可剪应力，设

34、置常数TAUMAX缺省为1.0E20 ）,这种限制剪应力的情况一般用于接触压力非常大的时候，以至于 用库仑理论计算出的界面剪应力超过了材料的屈服极限。一对TAUMAX勺一个合理高估为（ 是材料的 mises 屈服应力）。-9 -Gauss integration point图2 接触检查点位于高斯积分点上对无摩擦rough和bonded接触，接触单元刚度矩阵是对称的，而涉及到摩擦的接触问 题产生一个不对称的刚度，而在每次迭代使用不对称的求解器比对称的求解器需要更多的计 算时间，因此ANSYS序采用对称化算法。通过采用这种算法大多的摩擦接触问题能够使用 对称系统的求解器来求解。如果摩擦应力在整个

35、位移范围内有相当大的影响，并且摩擦应力的大小高度依赖于求解过程。对刚度阵的任何对称近似都可能导致收敛性的降低，在这种情况下，选择不对称求解选项( KEYOPT(6) =1)来改善收敛性。选择检查接触与否的位置接触检查点位于接触单元的积分点上，在积分点上，接触单元不穿透进入目标面，然 而，目标面能穿透进入接触面，看图ANSYS面接触单元使用GAUSS积分点作为缺省值，GAUSS积分点通常会比 Newto n-Cotes/robatto结点积分项产生更精确的结果，Newto n-cotes/lobatto使用结点本身作为积分点，通过 KEYOPT(4)来选择，你想使用的方法，然而，使用结点本身作为

36、积分 点仅应该用于角接触问题(看图4 3 )。-10 -图3 接触检查点位于高斯结点上然而，使用结点作为接触发现点，可能会导致其它的收敛性问题，例如“滑脱”（结点滑下目标面的边界）看图4 4，对大多的点一面的接触问题，我们推荐使用其它的点一面的接触单元，例如 CONTA26 CONTA4和口 CONTA49图4 结点滑脱调整初始接触条件在动态分析中，刚体运动一般不会引起问题，然而在静力分析中，当物体没有足够的约束时会产生刚体运动，有可能引起错误而终止计算。在仅仅通过接触的出现来约束刚体运动时，必须保证在初始几何体中， 接触对是接触的，换句话说，你要建立模型以便接触对是“刚好接触”的，然而这样作

37、可能会遇到以下问题：刚体外形常常是复杂的，很难决定第一个接触点发生在哪儿既使实体模型是在初始接触状态，在网格划分后余于数值舍入误差；两个面的单元 网格之间也可能会产生小的缝隙。接触单元的积分点和目标单元之间可能有小的缝隙。同理，在目标面和接触面之间可能发生过大的初始穿透，在这种情况下，接触单元可能会高估接触力，导致不收敛或民接触面之间脱离开接触关系。定义初始接触也许是建立接触分析模型时最重要的方面，因此，程序提供了几种方法来调整接触对的初始接触条件。注意：下面的技巧可以在开始分析时独立执行成几个联合起来执行，它们是为了消除由于生成网格造成的数值舍入误差而引起的小缝隙或穿透，而不是为了改正网格或

38、几何数据的错误。1、使用实常数ICONT来指定一个好的初始接触环，初始接触环是指沿着目标面的“调 整环”的深度，如果没有人为指定ICONT的值，程序会根据几何尺寸来给ICONT提供一个小值，同时输出一个表时什么值被指定的警告信息，对ICONT 一个正值表示相对于下面变形体单元厚度的比例因子，一个负值表示接触环的真正值，任何落在“调整环”敬域内的接触检查点被自动移到目标面上，(看图5(a)建议使用一个小的ICONT值否则，可能会发生大的不连续(看图4 5(b)Deformedintegraflori pointTarget surfa.ceR ic|id body-11 -(b)图5用ICON进

39、行接触面的调整(a)调整前(b)调整后2、使用实常数PMIN和PMAX来指定初始容许的穿透范围，当指定PMAX或 PMIN后，在 开始分析时，程序会将目标面移到初始接触状态，如果初始穿透大于PMAX程序会调整目标面的减少穿透，接触状态的初始调节仅仅通过平移来实现。对给定载荷或给定位移的刚性目标面将会执行初始接触状态的初始调节。同样，对没有指定边界条件的目标面也可以进行初始接触的调整。当目标面上的节点有给的零位移值时，使用PMAX和PMIN的初始调节将不会被执行。注意：ANSYS程序独立地处理目标面上节点的自由度，例如：如果你指定自中度UX值为“ 0”，那么，沿着X方向就没有初始调查，然而，在

40、Y和Z方向仍然会激活 PMAX和PMIN 选项。初始状态调整是一个迭代过程，程序最多进行20次迭代，如果目标面不能进入可接受的穿透范围，程序会给出一个警告信息，你可能需要调整你的初始几何模型。图46给出了一个初始接触调整迭代失败的例子。目标面的UY被约束住。因此，初始接触唯一容许的调整是在 X方向，然而，在这个问题中，刚性目标面在X方向的任何运动 都不会引起初始接触。-12 -Target sirface图6一个初始调整失败的例子3、使用实常数CNOF来指定接触面的偏移。正值表示将整个接触面向目标面移动，负值表示让整个接触面远离目标面。图7 消除初始始穿透真正的初始穿透包括两部分：?由几何形状

41、所导致的?由用户指定的接触表面偏移量(CNOF所引起的穿透KEY0PT(9)有如下几种选择?包括由几何形状和接触表面偏移量所引起的穿透，设置KEY0PT(9)= 0?忽略上述两方面所引起的穿透，设置KEY0PT(9)= 1?包括由用户指定的接触表面偏移量，但是忽略由几何形状所引起的初始穿透，设置 KEY0PT(9 )= 3在某些情况下，例如过盈装配问题，常常期望有过度的穿透，为了缓解收敛性困难，在第一个载荷步中设置 KEY0TI(9) =2来使过度穿透渐进到 0,看图4 &当使用这种方法 时，在第一个载荷步中不要给定其它任何载荷，也就是说要保证载荷是渐进的(KBC 0)。下面的KEY0PT选择

42、提供了渐变的能力。?让总的穿透量渐变为零(CN0F几何形状)，设置KEY0PT(9 )= 2?渐变用户指定的接触面过盈量而忽略由几何形状所引起的穿透，设置KEY0PT(9)=4?当使用上面的KEY0PT(9的设置时，在第一个荷载步中不要给定其它任何载荷，并且应设置KBC 0。同时应保证 Pin ball区域足够大，可以捕捉到初始过盈量。-13 -Target surfaceTarget surface图8 渐进初始穿透在开始分析时，程序会给出每个目标面的初始接触状态的输出信息，（在输出窗口或输出文件中），这个信息有助于决定每个目标面的最大穿透成最小间隙。对于给定的目标面如果没有发现接触，可能是

43、目标面离接触面太远（超出了Piaball区域或者是接触/目标单元已经被杀死。决定接触状态和 Pin ball区域。接触单元相对于目标面的运动和位置决定了接触单元的状态；程序检测每个接触单元并给出一种状态STAT=O未合的远区接触STAT=1未合的近区接触STAT=2滑动接触STAT=3粘合接触当目标面进入pin ball区域后，接触单元就被当作未合上的近区域接触，pin ball区域是以接触单元的积分点为中心的。使用实常数PINB来为pin ball指定一个比例因子（正值）或其实值（负值），缺省时，程序将pin ball区域定义为一个以 4*变形体单元厚度为半径的 圆（对2-D问题）或球（对

44、3-D问题）。检查接触的计算时间依赖于pin ball区域的大小，远区接触单元的计算是简单的且计算时间较少，近区接触计算将要接触的接触单元是较慢的，并是较复杂，当单元已经接触时，计算最为复杂。如果刚性面有好几个凸形区域，为了克服伪接触定义，设置一个合适的pin ball区域是 - 14 -程序够用KEYOP（I 11）来考虑壳（2-D和3-D）和梁（2-D）的厚度。缺省时，程序不考虑单元厚度，用中面来表示它。当设置有用的，名而对大多数问题，缺省值是合适的选择表面作用模式 通过设置 keyopt（12） 来选择下面的某种作用模式法问单边接触（KEYOPT（12）=0）.粗糙接触，用来模拟无滑动的

45、，表面相当粗糙的摩擦接触问题，这种设置对 应于摩擦系数无限大（MU ,因此用户定义的摩擦系的（MU被忽略KEYOPT12）=1）. 不分开的接触，用来模拟那种一旦接触就再不分开的问题，这种不分开是指对法向接触而言，允许有相对滑动。（ KEYOP （T 12） =2）. 绑定接触 用来模拟那种接触一旦发生表面在所省方向都被绑定的问题。 一旦接触就再也不能脱开也不允许有相对滑动（KEYOP（T 12） =3）用超单元建立接触模型面一面的接触单元能模拟刚体和另一个有运动的线弹性体的接触，而线弹性体又以体用超单元来建模， 这大大降低了进行接触迭代的自由度数， 记住任荷接触结定都必须是超单元 的主自由度

46、。既然超单元仅仅由一组保留的结点自由度组成，它没有用来定义接触的表面几何形状， 因此，必须在形成超单元之前在单元表面上形成接触单元， 来自超单元的信息包括结点连结 和组合刚度，但是没有材料特性和应力状态， （是否轴支称，平面应力或平面应变） ，一个限 制是接触单元的材料特性设置必须与形成超单元之前的原始单元的材料特性相同。使用KEY0PT（3）来提供接触分析的信息，对2-D单元（C0NTA171 CONTA17）关键字选项如下所示：. 不使用超单元（ KET0PT（ 3） =0）.轴对称（ KEY0P（T 3） =1）. 平面应变或单位厚度的平面应力（KEY0P（T 3） =2）-.需要厚度输

47、X的平面应力（KEYOPT 3） =3）,对这种情况使用实常数的 R2来 指定指定厚度对3 D单元（CONTA173 C0NTA17）关键字选项如下示:使用 H 单元（ KET0P（I 3） =0）使用超单元（ KEYOP（I 3） =1）考虑厚度影响KFTOPI（11） =1 时则考虑梁或壳的厚度，从底面 或顶面来计算接触距离， 建模时要考虑到厚度， 记住刚性目标面会向任一边移动半个梁或壳 单元的厚度，当使用壳单元 181 号时，在变形期间厚度的变化也将被考虑。使用时间步长控制 时间步长控制是一个自动时间步长特征，这个特征预测什么时间接触单元的状态，将发生变化或者需要二分当前的时间步长，使用

48、KEYOP（T 7）来选择下列四种行为之一来控制时间步长。 KEYOPT（7） =0 时不提供控制， KEYOP（T 7） =3 提供最多的控制。 KEYOP（I 7） =0，没有控制，时间步长的大小不受预测影响，当自动时间步长被激 活且允许一个很小的时间步长时，这个设置是合适的。KETOP（I 7） =1 如果一次迭代期间有太大的穿透发生或者接触状态突然变化，则进 行时间步长二分。KEYOP（I 7） =2对下一个子步预测一个合理的时间增量KETOP（I 7） =3对下一个子步，预测一个最小的时间增量- 15 -pilot 结点。pilot 结点上给定边界条件pilot 结点上定义一个质pi

49、lot ”结点，你不能另外CP来结点，那么使用死活单元选项面面的接触单元允许激活或杀死单元， 能够在分析的某一阶段中杀死这个单元而在以 后的阶段再重新激活它， 这个特征对于模拟复杂的金属成形过程是有用的、 在此过程的不同 分析阶段有多个目标需要和接触面相互作用， 回弹模拟常常需要在成形过程的后期移走刚性 工具。步骤七： 控制刚性目标的运动。 刚性目标面是按照物体的原始外形来建立的且整个面的运动是通过“pilot ”结点的给定的节点位移来定义的， （如果没有定义 “ pilot ”结点， 则通过刚性目标面上的不同结点。 ） 为了控制整个目标面的运动，在下面的任何情况下都必须使用目标面上作用着给定

50、的外力 目标面发生旋转目标面和其它单元相连（例结构质量单元） pilot 结点的厚度代表着整个刚性面的运动，你可以在 （位移、初速度 ） 集中载荷，转动等等，为了考虑刚体的质量，在 量单元。当使用 pilot 结点时，记住下面的几点局限性 每个目标面只能有一个“ Pilot 的结点.圆、圆锥、圆柱、球的第一个结点（结点工）是” 定义或改变 pilot 结点.程序忽略不是 pilot 结点的所有其它结点上的边界条件。. 只有“ pilot ”结点能与其它单元相连.当定义了“ pilot ”结点后，不能使用约束方程（ CF或结点来耦合（控制目标面的自由度，如果在刚性面上给定任意载荷或者约束，必须定

51、义“ pilot ”结点，只能在 pilot 结点上加载，如果没有使用“ pilot ” 则只能有刚体运动。在每个载荷步的开始， 程序检查每个目标面的边界条件， 如果下面的条件都满足， 程序将目标面作为固定处理：在目标面结点上没有明确定义边界条件或给定力. 目标面结点没有和其它单元相连. 没有目标面结上使用约束方程或结点来耦合 在每个载荷步的末尾，程序将会放松被内部设置的约束条件步骤八：给变形体单元加必要的边界条件现在可以按需要加上任你边界条件。加载过程与其它的分析类型相同 步骤九：定义求解和载荷步选项接触问题的收敛性随问题不同而不同， 下面列式了一些典型的在大多数面面的接触分 析中推荐使用的

52、选项时间步长必须足够小以描述适当的接触 。如果时间步长太大， 则接触力的光滑传递 会被破坏，设置精确时间步长的可信赖的方法是打开自动时间步长。命令： Autots, onGUI ： Main MenuSolution-load step opts-Time/FrequenceTime&Time step/Time& substeps. 如果在迭代期间接触状态变化， 可能发生不连续，为了避免收敛太慢， 使用修改的刚度阵，将牛顿一拉普森选项设置成 FULL命令： NROPT,FULL,OFFGUI ： Main MenuSolutionAnalysis options不要使用自下降因子，对面一面的

53、问题，自适应下降因子通常不会提供任何帮助，因 此我们建议关掉它。- 16 -设置合理的平衡迭代次数，一个合理的平衡迭代次数通常在25和50之间命令： NEQITGUI： Main MenuSolution-load step opts-NonlinearEquilibriwm iter因为大的时间增量会使迭代趋向于变得不稳定，使用线性搜索选项来使计算稳定化。命令： LNSRCHGUI： Main menusolution-load step opts-Nonlinearlins search 除非在大转动和动态分析中，打开时间步长预测器选项 命令： PREDGUI： main meansola

54、rion-load step opis-nonlinearpredictor 在接触分析中许多不收敛问题是由于使用了太大的接触刚度引起的，（实常数 FKN）检验是否使用了合适的接触刚度。步骤十：求解 现在可以对接触问题进行求解，求解过程与一般的非线问题求解过程相同步骤十一：检查结果 接触分析的结果主要包括位移、 应力、应变， 支反力， 和接触信息 （接触压力、 滑动等） 你可以在一般的后处理器（ post1 ）或时间历程后处理器（ post26 ）中查看结果。1.为了在 post1 中查看结果， 数据库文件所包含的模型必须与用于求解的模型相同。2.必须存在结果文件在 post1 中查看结果1.

55、从输出文件中查看分析是否收敛。 如果不收敛，你可能不想后处理，而更在乎为什么不收敛。 如果已经收敛，继续后处理。2.进入postl如果你的模型不在当前的数据库中，使用恢复命令（resume）来恢复它。命令： /post1GUI： main menuGeneral postproc3.读入所期望的载荷步和子步的结果，这可以通过载荷步和子步数也可以通过时间 来实现。命令： SETGUI： main menugenerad postproc4.使用下面的任何一个选项来显示结果选项：显示变形形状态命令： PLDISPGUI： main menugeneral postprocplot resnlt d

56、eformed shape选项：等值显示命令： PLNSOLPLESOLGUI： main menugeneral postprocplot resultcontour plot-noded solu或element和 solu使用这个选项来显示应力， 应变或其它项的等值图， 如果相邻的单元有不同的材料行为 （例如塑性或多弹性材料特性， 不同的材料类型， 或不同的死活属性） 则在结果显示时应避 免结点应力平均错误。也可以将定的接触信息用等值图显示出来，对2D 接触分析，模型用灰色表示，所要求显示的项将沿着接触单元存在的模型的边界以梯型面积表示出来，对3D 接触分析，模型将用灰色表示，而要求的项

57、在接触单元存在的2D 表面上等值显示。还可以等值显示单元表的数据和线性化单元数据。命令： PLETAB PLLSGUI： main menugeneral postprocElement TablePlot Element Table main menu General - 17 -PostprocPlot Results-Contour plot-line Elem Res 选项：列表显示 命令： PRNSOLPRESOLPRRSOLPRETABPRITERNSORTESORTGUI ： Main menuGeneral PostprocList ResultsNoded SolutionM

58、ain menuGeneral PostprocLost ResultsElement SolutionMain menuGeneral PostprocList ResultsReaction Solution在列表显示它们之前，可以用NSOR和ESOR俅对它们进行排序选项：动画 可以动画显示接触结果随时间的变化 命令： ANIMEGUI： Wility mennPlotctrls AnimateContours lner fime post26 中查看结果可以使用 post26 来查看一个非线性结构对加载历程的响应使用 post26 ，可以比较一个 变量陡另一个变量的变化关系， 例如， 可

59、以画出某个结点位移随给定载的曲线关系， 某个结 点的塑性应变与时间的关系，一个典型的 post26 后处理过程需要分以下几个步骤1.从输出文件中检查是否分析已经收敛2.求解已收敛，进入 post26 ，如果模型不在当前数据库中恢复它 命令： /Post26GUI:Main menuTimehist Postpro3.定义变量 命令： NSOLESOLRFORCEGUI： Main menuTime Hist PostproDefine Variable4.画曲线或列表显示 命令： PLVARPRVAREXTREMGUI： Main menuTime Hist PostprocGraph Var

60、iableMain menuTime Hist PostprocList Variarle Main menuTime Hist PostprocList Extremes点面接触分析我们能使用点面接触单元来模拟一个表面和一个结点的接触， 另外， 可以通过把表面 指定为一组结点，从而用点面接触来代表面面的接触。ANSYS 程序的点面接触单元允许下列非线性行为：有大变形的面-面接触分析接触和分开库仑摩擦滑动热传递点面的接触是一种在工程应用中普遍发生的现象， 例如：夹子、金属成形等等，工程 技术人员对由于结构之间的接触而产生的应力变形为和温度改变是感兴趣的。 使用点面的接触单元- 18 -在ANS