ansys 接触分析.doc

《ansys 接触分析.doc》由会员分享，可在线阅读，更多相关《ansys 接触分析.doc（13页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、【精品文档】如有侵权，请联系网站删除，仅供学习与交流ansys 接触分析.精品文档.ansys 接触分析 接触问题是一种高度非线性行为，需要较大的计算资源，为了进行实为有效的计算，理解问题的特性和建立合理的模型是很重要的。 接触问题存在两个较大的难点：其一，在你求解问题之前，你不知道接触区域，表面之间是接触或分开是未知的，突然变化的，这随载荷、材料、边界条件和其它因素而定；其二，大多的接触问题需要计算摩擦，有几种摩擦和模型供你挑选，它们都是非线性的，摩擦使问题的收敛性变得困难。一般的接触分类 接触问题分为两种基本类型：刚体柔体的接触，半柔体柔体的接触，在刚体柔体的接触问题中，接触面的一个或多个

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

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

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

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

6、的表面不连续。与点面接触单元比，需要较多的接触单元，因而造成需要较小的磁盘空间和CPU时间。允许多种建模控制，例如：绑定接触渐变初始渗透目标面自动移动到补始接触平移接触面（老虎梁和单元的厚度）支持死活单元 使用这些单元，能模拟直线（面）和曲线（面），通常用简单的几何形状例如圆、抛物线、球、圆锥、圆柱采模拟曲面，更复杂的刚体形状能使用特殊的前处理技巧来建模。执行接触分析 不同的接触分析类型有不同的过程，下面分别讨论面面的接触分析 在涉及到两个边界的接触问题中，很自然把一个边界作为“目标”面而把另一个作为“接触”面，对刚体柔体的接触，“目标”面总是刚性的，“接触”面总是柔性面，这两个面合起来叫作“

7、接触对”使用Targe169和Conta171或Conta172来定义2-D接触对，使用Targe170和Conta173或Conta174来定义3-D接触对，程序通过相同的实常收号来识别“接触对”。接触分析的步骤：执行一个典型的面面接触分析的基本步骤列示如下：1 建立模型，并划分网格2 识别接触对3 定义刚性目标面4 定义柔性接触面5 设置单元关键字和实常的6 定义控制刚性目标面的运动7 给定必须的边界条件8 定义求解选项和载荷步9 求解接触问题10查看结果步骤1：建立模型，并划分网格在这一步中，你需要建立代表接触体几何形状的实体模型。与其它分析过程一样，设置单元类型，实常的，材料特性。用恰

8、当的单元类型给接触体划分网格。命令：AMESHVMESHGUI：Main MenuPreprocessormeshMapped3 or4 Sided Main MenuPneprocessormeshmapped4 or 6 sided步骤二：识别接触对你必须认识到，模型在变形期间哪些地方可能发生接触，一是你已经识别出潜在的接触面，你应该通过目标单元和接触单元来定义它们，目标和接触单元跟踪变形阶段的运动，构成一个接触对的目标单元和接触单元通过共享的实常号联系起来。接触环（区域）可以任意定义，然而为了更有效的进行计算（主要指CPU时间）你可能想定义更小的局部化的接触环，但能保证它足以描述所需要的

9、接触行为，不同的接触对必须通过不同的实常数号来定义（即使实常数号没有变化）。由于几何模型和潜在变形的多样形，有时候一个接触面的同一区域可能和多个目标面产生接触关系。在这种情况下，应该定义多个接触对（使用多组覆盖层接触单元）。每个接触对有不同的实常数号。步骤三：定义刚性目标面刚性目标面可能是2D的或3D的。在2D情况下，刚性目标面的形状可以通过一系列直线、圆弧和抛物线来描述，所有这些都可以用TAPGE169来表示。另外，可以使用它们的任意组合来描述复杂的目标面。在3D情况下，目标面的形状可以通过三角面，圆柱面，圆锥面和球面来推述，所有这些都可以用TAPGE170来表示，对于一个复杂的，任意形状的

10、目标面，应该使用三角面来给它建模。 控制结点（Pilot） 刚性目标面可能会和“pilot结点“联系起来，它实际上是一个只有一个结点的单元，通过这个结点的运动可以控制整个目标面的运动，因此可以把pilot结点作为刚性目标的控制器。整个目标面的受力和转动情况可以通过pilot结点表示出来，“pilot结点”可能是目标单元中的一个结点，也可能是一个任意位置的结点，只有当需要转动或力矩载荷时，“pilot结点”的位置才是重要的，如果你定义了“pilot结点”ANSYS程序只在“pilot结点”上检查边界条件，而忽略其它结点上的任何约束。 对于圆、圆柱、圆锥、和球的基本图段，ANSYS总是使用条一个结

11、点作为“pilot结点” 基本原型 你能够使用基本几形状来模拟目标面，例如：“圆、圆柱、圆锥、球。直线、抛物线、弧线、和三角形不被允许、虽然你不能把这些基本原型彼此合在一起，或者是把它们和其它的目标形状合在一起以便形成一个同一实常数号的复杂目标面。但你可以给每个基本原型指定它自己的实常的号。 单元类型和实常数 在生成目标单元之前，首先必须定义单元类型（TARG169或TARG170）。 命令：ET GUI:main menupreprocessorElement Type Add/Edit/Delete 随后必须设置目标单元的实常数。 命令:Real GUI:main mennpreproce

12、ssorreal constants 对TARGE169和TARGE170仅需设置实常数R1和R2，而只有在使用直接生成法建立目标单元时，才需要从为指定实常数R1、R2，另外除了直接生成法，你也可以使用ANSYS网格划分工具生成目标单元，下面解释这两种方法。 使用直接生成法建立刚性目标单元 为了直接生成目标单元，使用下面的命令和菜单路径。 命令：TSHAP GUI：main menupreprocessormodeling-createElementsElem Attributes 随后指定单元形状，可能的形状有： straight line (2D) parabola (2-D) clock

13、wise arc(2-D) counterclokwise arc (2-D) circle(2-D) Triangle (3-D) Cylinder (3-D) Cone (3-D) Sphere (3-D) Pilot node (2-D和3-D) 一旦你指定目标单元形状，所有以后生成的单元都将保持这个形状，除非你指定另外一种形状。 然后你就可以使用标准的ANSYS直接生成技术生成结点和单元。 命令：N E GUI：main menupnoprocessor modeling- create nodes main menupnoprocessor modeling- createEleme

14、nts 在建立单元之后，你可以通过列示单元来验证单元形状 命令：ELIST GUI：utility menulistElementsNodes+Attributes 使用ANSYS网格划分工具生成刚性目标单元 你也可以使用标准的ANSYS网格划分功能让程序自动地生成目标单元，ANSYS程序将会以实体模型为基础生成合适的目标单元形状而忽略TSHAP命令的选项。 为了生成一个“PILOT结点”使用下面的命令或GUI路径： 命令：Kmesh GUI：main menuproprocessormeshing-meshkeypoints 注意： KMESH总是生成“PILOT结点” 为了生成一个2D目标

15、单元，使用下面的命令和GUI路径： ANSYS在每条直线上生成一条单一的线，在样条曲线上生成抛物线部分，在每条圆弧和倒角上生成圆弧部分，如果所有的圆弧形成一个封闭的圆，ANSYS生成一个单一的圆段。 命令：LMESH GUI：main menupneprocessormesling-meshlines 为了生成3D的目标单元，使用下面的命令或GUI路径。 如果实体模型的表面部分形成了一个完整的球，圆柱或圆锥，那么ANSYS程序自动生成一个基本的3D目标单元，因为生成较少的单元，从而使你分析计算更有效率，对任意形状的表面，应该使用Amesh命令来生成目标单元，在这种情况下，网格形状的质量不是重要

16、的，而目标单元的形状是否能完成好的模拟刚性面的表面几何形状显得更重要。 命令：AMESH GUI：main menupreprocessor-meshing-meshArea ANSYS在所有可能的面上推荐使用三角形的映射网格划分，如果在表面的边界上没有曲率，则在网格划分时，指定那条边界分为一分，下面的命令或GUI路径将尽可能的生成一个映射网格（如果不能进行映射，它将生成自由网格） 命令：MSHKFY，2 GUI：main menupreprocessor-meshling-mesh-Ares-Target Surf 建模和网格划分的注意点： 一个目标面可能由两个或多个面断的区域组成，你应该尽

17、可能地通过定义多个目标面来使接触区域局部比（每个目标面有一个不同的实常数号）刚性目标面上由的离散能足够指述出目标面的形状，过粗的网格离散可能导致收敛问题。如果刚性面有一个实的凸角，求解大的滑动问题时很难获得收敛结果，为了避免这些建模问题，在实体模型上，使用线或面的倒角来使尖角光滑比，或者在曲率突然变化的区域使用更细的网格。 注意：不能使用镜面对称技术（ARSYSM，LSYMM）来映射圆、圆柱、圆锥或球面到对称平面的另一边，因为每个实常数的设置不能同时赋给多个基本原型段。 检验目标面的接触方向。 目标面的结点号顺序是重要的，因为它定义了接触主向，对2D接触问题，当沿着目标线从第一个结点移向第二个

18、结点时，变形体的接触单元必须位于目标面的右边。 对3D接触问题，目标三角形单元号应该使刚性面的外法线方向指向接触面，外法线通过右手原则来定义 为了检查法线方向，显示单元坐标系 命令：/PSYMS，ESYS，1 GUI：Utility menuplotctrlssymbols 如果单元法向不指向接触面，选择单元反转表面的法向的方向。 命令：ESURF，REVE GUI：main menupreprocossorcreateElementon free surf步骤4：定义柔性体的接触面 为了定义柔性体的接触面，必须使用接触单元CONFA171或CONFA172（对2D）或CONTA173或CON

19、TA174（对3D）来定义表面 程序通过组成变形体表面的接触单元来定义接触表面，接触单元与下面覆盖的变形体单元有同样的几何特性，接触单元与下面覆盖的变形体单元必须处于同一阶次（低阶或高阶）下面的变形体单元可能是实体单元、壳单元、梁单元或超单元，接触面可能壳或梁单元任何一边。 与目标面单元一样，你必须定义接触面的单元类型，然后选择正确的实常数号（实常数号必须与它对应目标的实常数号相同）最后生成接触单元。 单元类型： 下面简单描述四种类型的接触单元 CONTA171：这是一种2D，2个结点的低附线单元，可能位于2D实体，壳或梁单元的表面 CONTA172：这是一个2D的，3结点的高阶抛物线形单元，

20、可能位于有中结点的2D实体或梁单元的表面 CONTA173：这是一个3D的，4结点的低阶四边形单元可能位于3D实体或壳单元的表面，它可能褪化成一个结点的三角形单元。 CONTA174：这是一个3D，8结点的高阶四边形单元，可能位于有中结点的3D实体或壳单元的表面，它可能褪化成6结点的三角形单元。 不能在高阶柔性体单元的表面上分成低阶接触单元，反之也不行，不能在高阶接触单元上消去中结点。 命令：ET GUI：main menupreprocessorElement typeAdd/Edit/Delete 实常数和材料特性 在定义了单元类型之后，需要选择正确的实常数的设置，每个接触对的接触面和目标

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

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

23、法向相同，使用下表面生成接触单元，则它们的外法向与梁或壳单元的法向相反，如果下面的单元是实体单元，则TOP或BOTTOM选项不起作用 3、检查接触单元外法线的方向，当程序进行是否接触的检查时，接触面的外法线方向是重要的，对3D单元，按结点程序号以右手定则来决定单元的外法向，接面的外法向应该指向目标面，否则，在开始分析计算时，程序可能会认为有面的过度渗透而很难找到初始解。在此情况下，程序一般会立即停止执行，你可以检查单元外法线方向是否正确。 命令：/PSYMB GUI：Utility menuplotctrlssymbols 当发现单元的外法线方向不正确时，必须通过倒不正确单元的结点号来改变它们

24、。 命令：ESURF，REVE GUI：main menupreprocossorCreateElements on free surf前提：有限元模型。已识别接触面及目标面。（*可应用自由度耦合来替代接触。）定义接触单元及实常数（刚性）目标单元 TARGE169 TARGE170 ；（柔性）接触单元 CONTA171CONTA172。 *Commands* ET,K,169 !K - 指定的单元编号 ET,K+1,172实常数 一个接触对对应同一个实常数号。 TARGE单元的实常数包括：R1、R2 定义目标单元几何形状 CONTA单元的实常数包括： No. Name Description

25、1 R1 Target circle radius（刚性环半径） 2 R2 Superelement thickness（单元厚度） *3 FKN Normal penalty stiffness factor（法向接触刚度因子） *4 FTOLN Penetration tolerance factor（最大允许的穿透） *5 ICONT Initial contact closure（初始闭合因子） 6 PINB Pinball region（“Pinball”区域） *7 PMAX Upper limit of initial allowable penetration（初始穿透的最大值

26、） *8 PMIN Lower limit of initial allowable penetration（初始穿透的最小值） *9 TAUMAX Maximum friction stress（最大的接触摩擦） *10 CNOF Contact surface offset（施加于接触面的正或负的偏移值） 11 FKOP Contact opening stiffness or contact damping *12 FKT Tangent penalty stiffness factor（切向接触刚度） 13 COHE Contact cohesion（滑动抗力粘聚力） 14 TCC T

27、hermal contact conductance（热接触传导系数） 15 FHTG Frictional heating factor（摩擦耗散能量的热转换率） 16 SBCT Stefan-Boltzmann constant 17 RDVF Radiation view factor 18 FWGT Heat distribution weighing factor 19 ECC Electric contact conductance 20 FHEG Joule dissipation weight factor 21 FACT Static/dynamic ratio（静摩擦系数和

28、动摩擦系数的比率） 22 DC Exponential decay coefficient（摩擦衰减系数） 23 SLTO Allowable elastic slip 24 TNOP Maximum allowable tensile contact pressure 25 TOLS Target edge extension factor 附注： +值作为比例因子，-值作为绝对值； 带*号的实常数比较重要，关乎接触分析的收敛； 一般实常数可为缺省值。 *Commands* REAL,K !K - 指定169和172的共同的实常数编号 ET,K,169 !K - 指定的单元编号 ET,K+1

29、,172 R,K,FKN0,FTOLN0 !指定实常数16的值 RMORE,FKT0, !指定实常数712的值 * *材料参数 仅需定义摩擦系数 *Commands* MP,MU,J,0.6 !J为材料编号，切勿与此前定义的重复；0.6为接触摩擦系数。 * * 附注：当接触对数目较多，且各接触对的性质相似时，可利用*Do循环语句定义单元、实常数。 *Commands* MP,MU,J,0.6 *DO,K,8,14,2 !Do循环语句定义单元、实常数。 REAL,K ET,K,169 ET,K+1,172 !KEYOPT,K+1,5,3 KEYOPT,K+1,9,3 !KEYOPT,K+1,8,

30、0 !KEYOPT,K+1,10,1 R,K,FKN0,FTOLN0 ! ,CNOF0, RMORE,FKT0 NROPT,UNSYM *ENDDO * * 生成单元应用直接生成法（ESURF）建立单元. *Commands* REAL,K !指定实常数编号 type,K !指定TARGE单元编号 mat,J !指定材料编号 lsel,1 !选择编号为1的接触边界线 nsll,1 !选择依附于线1的节点 esurf !生成刚性目标单元 type,K+1 !指定CONTA单元编号 lsel,2 nsll,2 esurf !生成柔性接触单元 allsel 关乎收敛的若干参数（实常数及单元关键字）设置接触刚度设置 FKN FKT KEYOPT(10)初始穿透设置 CNOF KEYOPT(5) KEYOPT(9) *Commands* KEYOPT,K+1,9,3最大允许穿透 FTOLN可应用自由度耦合来替代接触 NUMMRG (命令)