ANSYS高级接触问题1372397.pptx
《ANSYS高级接触问题1372397.pptx》由会员分享,可在线阅读,更多相关《ANSYS高级接触问题1372397.pptx(123页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、1 接触分类刚柔一个表面是完全刚性的除刚体运动外无应变、应力和变形,另一表面为软材料构成是可变形的。只在一个表面特别刚硬并且不关心刚硬物体的应力时有效。柔柔两个接触体都可以变形。第1页/共123页2 接触单元ANSYS采用接触单元来模拟接触问题:跟踪接触位置;保证接触协调性(防止接触表面相互穿透);在接触表面之间传递接触应力(正压力和摩擦)。接触单元就是覆盖在分析模型接触面上的一层单元。在ANSYS中可以采用三种不同的单元来模拟接触:面一面接触单元;点一面接触单元;点一点接触单元。第2页/共123页不同的单元类型具有完全不同的单元特性和分析过程。1.面一面接触单元用于任意形状的两个表面接触不必
2、事先知道接触的准确位置;两个面可以具有不同的网格;支持大的相对滑动;支持大应变和大转动。例如:面一面接触可以模拟金属成型,如轧制过程。第3页/共123页第4页/共123页2.点一面接触单元用于某一点和任意形状的面的接触可使用多个点面接触单元模拟棱边和面的接触;不必事先知道接触的准确位置;两个面可以具有不同的网格;支持大的相对滑动;支持大应变和大转动。例:点面接触可以模拟棱边和面之间的接触第5页/共123页第6页/共123页3.点点接触单元用于模拟单点和另一个确定点之间的接触。建立模型时必须事先知道确切的接触位置;多个点点接触单元可以模拟两个具有多个单元表面间的接触;每个表面的网格必须是相同的;
3、相对滑动必须很小;只对小的转动响应有效。例如:点一点接触可以模拟一些面的接触。如地基和土壤的接触第7页/共123页第8页/共123页3 关于耦合和约束方程的应用如果接触模型没有摩擦,接触区域始终粘在一起,并且分析是小挠度、小转动问题,那么可以用耦合或约束方程代替接触。使用耦合或约束方程的优点是分析还是线性的 第9页/共123页接触问题的一般特性1 接触刚度1、所有的ANSYS接触单元都采用罚刚度(接触刚度)来保证接触界面的协调性第10页/共123页在数学上为保持平衡,需要有穿透值然而,物理接触实体是没有穿透的 分析者将面对困难的选择:小的穿透计算精度高,因此接触刚度应该大;然而,太大的接触刚度
4、会产生收敛困难:模型可能会振荡,接触表面互相跳开。接触刚度是同时影响计算精度和收敛的最重要的参数。你必须选定一个合适的接触刚度。除了在表面间传递法向压力外,接触单元还传递切向运动(摩擦)。采用切向罚刚度保证切向的协调性。(图12)作为初值,可采用:Ktangent=0.01Knormal切向罚刚度与法向罚刚度以同样的方式对收敛性和计算精度产生影响。第11页/共123页2、接触刚度的选取选定一个合适的接触刚度值需要一些经验。对于面一面接触单元,接触刚度通常指定为基体单元刚度的一个比例因子。开始估计时,选用FKN=1.0 大面积实体接触FKN=0.01-0.1 较柔软(弯曲占主导的)部分另外,也可
5、以指定一个绝对刚度值,单位:(力/长度)/面积。点一点(除CONTA178)和点面接触单元需要为罚刚度KN输入绝对值:初始估计时:对于大变形:0.1*E KN 1.0*E对于弯曲:0.01*E KN 0.1*EE 为弹性模量第12页/共123页3、选取接触刚度的指导:Step 1.开始采用较小的刚度值Step 2.对前几个子步进行计算Step 3.检查穿透量和每一个子步中的平衡迭代次数在粗略的检查中,如以实际比例显示整个模型时就能观察到穿透,则穿透可能太大了,需要提高刚度重新分析。如果收敛的迭代次数过多(或未收敛),降低刚度重新分析。注意:罚刚度可以在载荷步间改变,并且可以在重启动中调整。牢记
6、:接触刚度是同时影响计算精度和收敛性的最重要的参数。如果收敛有问题,减小刚度值,重新分析在敏感的分析中,还应该改变罚刚度来验证计算结果的有效性。在分析中减小刚度范围,直到结果(接触压力、最大SEQV等)不再明显改变。第13页/共123页2 摩擦1、两个接触体的剪切或滑动行为可以是无摩擦的或有摩擦的无摩擦时允许物体没有阻力地相互滑动;有摩擦时,物体之间会产生剪切力。2、摩擦消耗能量,并且是路径相关行为。为获得较高的精度,时间步长必须小(图2-1)第14页/共123页 图213、ANSYS中,摩擦采用库仑模型,并有附加选项可处理复杂的粘着和剪切行为。库仑法则是宏观模型,表述物体间的等效剪力FT不能
7、超过正压力FN的一部分:FT Preprocessor Modeling Create Elements Surf/Contact Surf to Surf(ESURF)对于直接生成刚性目标面,在建立目标单元之前需要要指定附加的单元属性TSHAP第23页/共123页第24页/共123页刚性目标面的自动划分不需要TSHAP。ANSYS能根据实体模型确定合适的目标单元形状。划分线(LMESH)2D刚性目标面划分面(AMESH)3D刚性目标面创建关键点(KMESH)控制节点(Pilot)刚性目标面能与控制点联系起来,Pilot实际上是只有一个节点的单元,通过这个节点的运动可以控制整个目标面的运动。A
8、NSYS只在Pilot节点上检查边界条件而忽略其它节点的约束。对可变形体目标面建立目标单元的步骤是:1.先选择可变形体表面上的节点2.然后在可变形体上建立目标单元Main Menu Preprocessor Modeling Create Elements Surf/Contact Surf to Surf第25页/共123页ANSYS将根据基体的网格确定目标单元形状和外法线方向。检查外法线方向(这在自动划分刚性目标面时非常重要)图3-3打开单元坐标系标志并重绘单元/PSYMS,ESYS,1目标单元外法线方向应该指向接触面。如果单元法向不指向接触面,用命令使之反转:ESURF,REVE第26页
9、/共123页例:Seal.dat(图3-3)Step 5.建立接触面单元设置接触单元属性、选择可变形体表面节点,并在可变形体上建立接触单元(过程与在可变形体上建立目标单元相同)Main Menu Preprocessor Modeling Create Elements Surf/Contact Surf to Surf这些接触单元与基体有同样的阶数(低阶或高阶)。注意,在壳或梁单元上建立目标单元或接触单元时,可以选择要在梁或壳单元的顶层还是底层建立单元。图3-33-3第27页/共123页在选择柔体表面上的节点时,如果你确定某一部分节点永远不会接触到目标面时,可以忽略它,以减少计算时间。接触面
10、的外法向应指向目标面。如果发现外法线方向不正确,用下列命令修改之ESURF,REVEStep 6.在有限元模型上施加边界条件如果目标面是刚性面,目标面将会自动固定。定义了Pilot点ANSYS只检查该点的边界条件,忽略目标面上其它节点约束。控制点能控制目标面的运动。对Seal.dat施加的边界条见图33。第28页/共123页Step 7.定义求解选项和载荷步,以下是默认设置推荐使用N.L求解自动控制使用不带自适应下降的full Newton-Raphson法求解时间步必须足够小。使用自动时间步。子步数的最大值(NSBMX)应较大,最小值(NSBMIN)应较小Step 8.求解Step 9.后处
11、理结果包括位移、应力、应变和接触等信息。接触压力、摩擦应力、总应力、接触侵入、接触间隙距离、滑动距离和接触状态都可以从/POST1或/POST26中得到。第29页/共123页面一面接触分析实例(建立接触对不通过接触向导)实例1:弹性环装配第30页/共123页轴对称施加位移载荷:在L45线上施加0.4的-Y向位移打开几何N.L开关(NLGEOM,on)设Time=0.4并为自动时间步给出子步数(20,500,10)给出输出控制(要求输出每一子步结果)求解并查看输出和监视文件重启动分析亦将-Y向偏移量增加到0.55使咬接装配的第2个齿咬合。问题以弯曲为主,设FKN=0.1第31页/共123页GUI
12、方式:Step 1.恢复数据库文件 Snap.db(Snap.db数据库包含此有限元模型的完整几何模型、材料、边界条件。但未定义目标单元与接触单元)。Utility Menu File Resume from 选 Snap.db 【OK】基体单元:Plane42,1 轴对称材料:MAT,1 EX=0.175e6;UXY=0.35;MU=0.0 MAT,2 EX=90000;NUXY=0.35 约束节点 163 UY=0 182 UY=0第32页/共123页Step 2.创建目标面单元类型Main Menu Preprocessor Element Typer Add/Edit/Delete【A
13、dd】contact 2D target 169Element typer reference number=2【OK】或命令:/PREP 7ET,2,Target 169Step 3.建立接触面单元类型Main Menu Preprocessor Element Typer Add/Edit/Delete【Add】contact 2nd Surf 171Element typer reference number=3【OK】【close】或命令:ET,3,conta 171第33页/共123页Step 4.指定接触法向刚度Main Menu Preprocessor Real Constan
14、ts Add/Edit/Delete【Add】Type 3 CONTA 171【OK】Real Constant Set No.=1Normal Penalty Stiffness FKN=0.1(对弯曲问题采用初值0.1)Penetration tolerance FTOLN=0.1(不用此,不收敛)或命令:R,1,0.1,0.1第34页/共123页确定罚刚度 FKN值通常在0.01-10之间,对于体积变形问题用1.0(默认),对弯曲问题用0.1确定侵入容差:侵入容差(FTOLN)是与接触单元下面的实体单元深度(h)相乘的比例因子。若此值太小会引起收敛困难,绝对不要用太小的容差!增大罚刚度(
15、FKN)将减少侵入。将FKN增大100倍会相应地减少侵入,但是接触压力只改变5。如不收敛可调整FKN或FTOLN值重新运行。检查侵入和每个子步的平衡迭代数,如果收敛受入侵容差的驱使可能是FKN值估计不足或FTOLN值太小。如果需要多次迭代才能使残值收敛而不是侵入。FKN值可能估计得太高。第35页/共123页Step 5.创建目标单元(1).为目标面选择线Utility Menu Select Entities拾取线(图3-5)【OK】或命令:LSEL,S,2,4LSEL,A,15,18 图3-6LSEL,A,63(2).选择附于线上的全部节点(图3-6)Utility Menu Select
16、EntitiesNodesAttached to lines,all【OK】或命令:NSLL,S,1NPlot图3-53-5图3-63-6第36页/共123页(3).设置单元属性Main Menu Preprocessor Create Element Elem Attributes或命令:Type,2Mat,1Real,1第37页/共123页(4).创建目标单元Main Menu Preprocessor Modeling Create Elements Surf/Contact Surf to Surf注意:如果基体单元是2D平面或3D实体单元“Tlab”无效,如果基体单元是2D单元,“S
17、hape”无效。选【pick all】或命令:ESURF(图3-7)图3-73-7第38页/共123页Step 6.创建接触单元(1).为接触面选线Utility Menu Select Entities Lines图3-8By Num/pickFrom Full【OK】选线(图3-8)【OK】或命令:LSEL,s,33,34LSEL,a,43,44(2).选择附于选定线上的全部节点图3-9Utility Menu Select Entities NodesAttached to Lines,allFrom Full【OK】或命令:NSLL,S,1Nplot(图3-9)图3-83-8图3-93
18、-9第39页/共123页(3).设置单元属性Main Menu Preprocessor Create Element Elem AttributesElement type number=3 ConTA 171Material number=1Real constant set number=1【OK】或命令:Type,3Mat,1Real,1(4).创建接触面单元图3-11Main Menu Preprocessor Modeling Create Elements Surf/Contact Surf to Surf图3-10【OK】【PICK ALL】(图3-10)或命令:ESURF图3
19、-103-10图3-113-11第40页/共123页Step 7.选择所有选项并画单元Utility Menu Select EverythingUtility Menu Plot Elements或命令:Allsel,allEplot(图3-11)Step 8.求解/soluantype,staticnlgeom,onsolc,ontime,0.4nsubst,20,500,10outres,all,alldl,45,uy,-0.4solvedl,45,uy,-0.55solvefiniStep 9.后处理第41页/共123页第42页/共123页4 应用接触向导创建接触对所有面面接触的单元选
20、项和参数都可以通过接触向导来控制。使用接触向导进行接触分析方便快捷:自动定义单元类型和实常数快速得到接触选项和参数快速显示和反转接触法向使用接触向导必须首先对基体进行网格划分,否则不能激活接触向导。下面仍以弹性环装配接触问题为例。采用接触向导完成接触对的创建来说明接触向导的使用方法。Step 1.恢复数据文件Utility Menu File Resume from选 Snap.db【OK】(图4-1)此数据库Snap.db包括此有限元模型的几何、基体单元、分网、材料和边界条件。但不包括接触单元与目标单元及接触有关的参数。图4-14-1第43页/共123页Step 2.启动接触向导Main M
21、enu Preprocessor Modeling Create Contact Pair第44页/共123页Step 3.创建接触对 图4-3 图4-4【Pick Target】(图43),拾取线(见图35)【OK】,再定义接触面 返回图44【Next】(图45)第45页/共123页图45【opening setting】图46第46页/共123页图46第47页/共123页Step 4.设置接触参数【OK】返回图45【create】【Finish】返回图42第48页/共123页退出Step 5.求解求解过程同前图42第49页/共123页5 实例2平面密封圈计算验证为接触刚度估计一个合适的初始
22、值的重要性验证接触分析中摩擦的应用弹性库仑摩擦验证基于接触单元临界状态变化的时间步预测如何会对计算效率有害模型描述:2D超弹平面密封圈1/2对称,密封圈压缩模拟第50页/共123页分析此模采用两种方法:1.不使用接触向导,建立命流文件Seal.inp2.使用接触向导创建接触对,并用GUI方式创建Pilot节点。本例具有:几何N.L(大应变与大变形)材料N.L(超弹)接触N.L 数据文件:Seal.inp第51页/共123页Step 1.建模:单元类型、实常数、材料特性、基体分网/prep7et,1,56 !HYPER56 2D 4node U_P Hyperelastic Solidmp,nu
23、xy,1,0.49tb,mooney,1tbdata,1,80 !C10=80tbdata,2,20 !C01=20k,1$k,2,0.333,0$k,3,0.867,0.867 k,4,1.1,0.867k,5,1.1,1$k,6,0.8,1$k,7,0.267,0.133 k,8,0,0.133l,1,2*repeat,7,1,1 !将l,1,2命令重复7次l,8,1lfil,1,2,0.20$lfil,2,3,0.15$lfil,5,6,0.20lfil,6,7,0.15$lfil,7,8,0.05$lfil,8,1,0.05al,all !应用所有选择的线生成面k,25,-0.6,0$
24、k,26,1.1,0$k,27,-0.6,1.0$k,28,1.1,1.0lstr,26,25 !L9lstr,27,28 !L10图5-25-2第52页/共123页!*基体分网*lesize,8,2$lesize,13,4$lesize,14,4esize,0.035type,1mat,1amesh,allsave,seal,dbStep 2.创建接触对,先定义单元和接触参数et,2,169 !Target169 2D Target Segmentet,3,171 !Conta171 2D Surface to Surface Contactmp,mu,1,0 !用材料特性定义摩擦,本例先无
25、摩擦计算然后再有摩擦计算,比较结果r,1,0.1,0.1 !用实常数定义接触高级选项对于弯曲为主的问题,采用接触刚度FKN=0.1作为初始估计值,FTOLN=0.1 lagrange穿透比例系数(缺省值)r,2,0.1,0.1 !如果FKN=1.0则不收敛第53页/共123页!*创建接触对1(定义实常数和分网)*type,2 !目标面1 Target169real,1 !Target169 的实常数mat,1 !接触面1定义材料与超弹单元同lesize,15,1 !目标面1(L15)分网尺寸(图53)lmesh,15 !目标面1(L15)分网,采用自动分网,此刚性面自动约束。检查外法线方向【O
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- ANSYS 高级 接触 问题 1372397
限制150内