ANSYS高级接触分析优秀PPT.ppt
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_05.gif)
《ANSYS高级接触分析优秀PPT.ppt》由会员分享,可在线阅读,更多相关《ANSYS高级接触分析优秀PPT.ppt(191页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。
1、ANSYS 高级接触问题 接触问题概述在工程中会遇到大量的接触问题,如齿轮的啮合、法兰联接、机电轴承接触、卡头与卡座、密封、板成形、冲击等等。接触是典型的状态非线性问题,它是一种高度非线性行为。接触例子如下图:分析中常常须要确定两个或多个相互接触物体的位移、接触区域的大小和接触面上的应力分布。接触分析存在两大难点:在求解之前,你不知道接触区域的范围;表面之间是接触还是分开是未知的;表面之间突然接触或突然不接触会导致系统刚度的突然变更。大多数接触问题须要计算摩擦。摩擦是与路径有关的现象,摩擦响应还可能是杂乱的,使问题求解难以收敛。ANSYS 高级接触问题 1 接触分类刚柔一个表面是完全刚性的除刚
2、体运动外无应变、应力和变形,另一表面为软材料构成是可变形的。只在一个表面特殊刚硬并且不关切刚硬物体的应力时有效。柔柔两个接触体都可以变形。2 接触单元ANSYS 接受接触单元来模拟接触问题:跟踪接触位置;保证接触协调性(防止接触表面相互穿透);在接触表面之间传递接触应力(正压力和摩擦)。接触单元就是覆盖在分析模型接触面上的一层单元。在 ANSYS 中可以接受三种不同的单元来模拟接触:面一面接触单元;点一面接触单元;点一点接触单元。不同的单元类型具有完全不同的单元特性和分析过程。1.面一面接触单元用于随意形态的两个表面接触 不必事先知道接触的精确位置;两个面可以具有不同的网格;支持大的相对滑动;
3、支持大应变和大转动。例如:面一面接触可以模拟金属成型,如轧制过程。2 接触单元2 接触单元2.点一面接触单元用于某一点和随意形态的面的接触 可运用多个点面接触单元模拟棱边和面的接触;不必事先知道接触的精确位置;两个面可以具有不同的网格;支持大的相对滑动;支持大应变和大转动。例:点面接触可以模拟棱边和面之间的接触2 接触单元2 接触单元3.点点接触单元用于模拟单点和另一个确定点之间的接触。建立模型时必需事先知道精确的接触位置;多个点点接触单元可以模拟两个具有多个单元表面间的接触;每个表面的网格必需是相同的;相对滑动必需很小;只对小的转动响应有效。例如:点一点接触可以模拟一些面的接触。如地基和土壤
4、的接触2 接触单元2 接触单元3 关于耦合和约束方程的应用假如接触模型没有摩擦,接触区域始终粘在一起,并且分析是小挠度、小转动问题,那么可以用耦合或约束方程代替接触。运用耦合或约束方程的优点是分析还是线性的 接触问题的一般特性1 接触刚度1、全部的 ANSYS 接触单元都接受罚刚度(接触刚度)来保证接触界面的协调性在数学上为保持平衡,须要有穿透值然而,物理接触实体是没有穿透的 分析者将面对困难的选择:小的穿透计算精度高,因此接触刚度应当大;然而,太大的接触刚度会产生收敛困难:模型可能会振荡,接触表面相互跳开。接触刚度是同时影响计算精度和收敛的最重要的参数。你必需选定一个合适的接触刚度。除了在表
5、面间传递法向压力外,接触单元还传递切向运动(摩擦)。接受切向罚刚度保证切向的协调性。(图12)作为初值,可接受:Ktangent=0.01 Knormal切向罚刚度与法向罚刚度以同样的方式对收敛性和计算精度产生影响。2、接触刚度的选取 选定一个合适的接触刚度值须要一些阅历。对于面一面接触单元,接触刚度通常指定为基体单元刚度的一个比例因子。起先估计时,选用 FKN=1.0 大面积实体接触 FKN=0.01-0.1 较松软(弯曲占主导的部分)另外,也可以指定一个确定刚度值,单位:(力/长度)/面积。对于点一点(除 CONTA178)和点面接触单元须要为罚刚度 KN 输入确定值:初始估计时:对于大变
6、形:0.1*E KN 1.0*E 对于弯曲:0.01*E KN 0.1*E E 为弹性模量3、选取接触刚度的指导:Step 1.起先接受较小的刚度值Step 2.对前几个子步进行计算Step 3.检查穿透量和每一个子步中的平衡迭代次数在粗略的检查中,如以实际比例显示整个模型时就能视察到穿透,则穿透可能太大了,须要提高刚度重新分析。假如收敛的迭代次数过多(或未收敛),降低刚度重新分析。留意:罚刚度可以在载荷步间变更,并且可以在重启动中调整。牢记:接触刚度是同时影响计算精度和收敛性的最重要的参数。假如收敛有问题,减小刚度值,重新分析在敏感的分析中,还应当变更罚刚度来验证计算结果的有效性。在分析中减
7、小刚度范围,直到结果(接触压力、最大SEQV 等)不再明显变更。2 摩擦1、两个接触体的剪切或相互滑动行为可以是无摩擦的,也可以是有摩擦的 无摩擦时允许物体没有阻力地相互滑动;有摩擦时,物体之间会产生剪切力(摩擦力)。2、摩擦消耗能量,并且是路径相关行为。为获得较高的精度,时间步长必需很小(图2-1)图213、ANSYS 中,摩擦接受库仑模型,并有附加选项可处理困难的粘着和剪切行为。库仑法则是宏观模型,表述物体间的等效剪力 FT 不能超过正压力 FN 的一部分:FT Preprocessor Modeling Create Elements Surf/Contact Surf to Surf(
8、ESURF)对于干脆生成刚性目标面,在建立目标单元之前须要要指定附加的单元属性 TSHAP刚性目标面的自动划分不须要 TSHAP。ANSYS 能依据实体模型确定合适的目标单元形态。划分线(LMESH)2D 刚性目标面 划分面(AMESH)3D 刚性目标面 创建关键点(KMESH)限制节点(Pilot)刚性目标面能与限制点联系起来。Pilot 事实上是只有一个节点的单元,通过这个节点的运动可以限制整个目标面的运动。ANSYS 只在 Pilot 节点上检查边界条件而忽视其它节点的约束。对可变形体目标面建立目标单元的步骤是:1.先选择可变形体表面上的节点 2.然后在可变形体上建立目标单元 Main
9、Menu Preprocessor Modeling Create Elements Surf/Contact Surf to SurfANSYS 将依据基体的网格确定目标单元形态和外法线方向。检查外法线方向(这在自动划分刚性目标面时特殊重要)图 3-3 打开单元坐标系标记并重绘单元 /PSYMS,ESYS,1 目标单元外法线方向应当指向接触面。假如单元法向不指向接触面,用叮嘱使之反转:ESURF,REVE例:Seal.dat(图3-3)Step 5.建立接触面单元 设置接触单元属性、选择可变形体表面节点,并在可变形体上建立接触单元(过程与在可变形体上建立目标单元相同)Main Menu Pr
10、eprocessor Modeling Create Elements Surf/Contact Surf to Surf(ESURF)这些接触单元与基体有同样的阶数(低阶或高阶)。留意,在壳或梁单元上建立目标单元或接触单元时,可以选择要在梁或壳单元的顶层还是底层建立单元。图图3-33-3 在选择柔体表面上的节点时,假如你确定某一部分节点恒久不会接触到目标面,可以忽视它,以削减计算时间。接触面的外法向应指向目标面。假如发觉外法线方向不正确,用下列叮嘱修改之 ESURF,REVEStep 6.在有限元模型上施加边界条件 假如目标面是刚性面,目标面将会自动固定。定义了 Pilot 点 ANSYS
11、只检查该点的边界条件,忽视目标面上其它节点约束。限制点能限制目标面的运动。对 Seal.dat 施加的边界条见图 33。Step 7.定义求解选项和载荷步,以下是默认设置 举荐运用N.L求解自动限制 运用不带自适应下降的 full Newton-Raphson 法求解 时间步必需足够小。运用自动时间步。子步数的最大值(NSBMX)应较大,最小值(NSBMIN)应较小Step 8.求解Step 9.后处理 结果包括位移、应力、应变和接触等信息。接触压力、摩擦应力、总应力、接触侵入、接触间隙距离、滑动距离和接触状态都可以从/POST1 或/POST26 中得到。面一面接触分析实例(建立接触对不通过
12、接触向导)实例1:弹性环装配轴对称施加位移载荷:在 L45 线上施加 0.4的 Y 向位移打开几何 N.L 开关(NLGEOM,on)设 Time=0.4 并为自动时间步给出子步数(20,500,10)给出输出限制(要求输出每一子步结果)求解并查看输出和监视文件重启动分析亦将 Y 向偏移量增加到 0.55使咬接装配的第 2 个齿咬合。问题以弯曲为主,设 FKN=0.1GUI 方式:Step 1.复原数据库文件 Snap.db(Snap.db 数据库包含此有限元模型的完整几何模型、材料、边界条件。但未定义目标单元与接触单元)。Utility Menu File Resume from 选 Sna
13、p.db 【OK】基体单元:Plane42,1 轴对称材料:MAT,1 EX=0.175e6;NUXY=0.35;MU=0.0 MAT,2 EX=90000;NUXY=0.35 约束节点 163 UY=0 182 UY=0Step 2.创建目标面单元类型Main Menu Preprocessor Element Typer Add/Edit/Delete【Add】contact 2D targe169 Element typer reference number=2【OK】或叮嘱:/PREP 7 ET,2,Targe169Step 3.建立接触面单元类型Main Menu Preproces
14、sor Element Typer Add/Edit/Delete【Add】contact 2nd Surf 171 Element typer reference number=3【OK】【close】或叮嘱:ET,3,conta171Step 4.指定接触法向刚度 Main Menu Preprocessor Real Constants Add/Edit/Delete 【Add】Type 3 CONTA171 【OK】Real Constant Set No.=1 Normal Penalty Stiffness FKN=0.1(对弯曲问题接受初值0.1)Penetration tole
15、rance FTOLN=0.1(不用此,不收敛)或叮嘱:R,1,0.1,0.1确定罚刚度 FKN 值通常在 0.01-10 之间,对于体积变形问题用 1.0(默认),对弯曲问题用 0.1确定侵入容差:侵入容差(FTOLN)是与接触单元下面的实体单元深度(h)相乘的比例因子。若此值太小会引起收敛困难,确定不要用太小的容差!增大罚刚度(FKN)将削减侵入。将 FKN 增大 100 倍会相应地削减侵入,但是接触压力只变更 5。如不收敛可调整 FKN 或 FTOLN 值重新运行。检查侵入和每个子步的平衡迭代数,假如收敛受侵入容差的驱使可能是FKN 值估计不足或 FTOLN 值太小。假如须要多次迭代才能
16、使残值收敛而不是侵入。FKN 值可能估计得太高。Step 5.创建目标单元(1)为目标面选择线 Utility Menu Select Entities 拾取线(图3-5)【OK】或叮嘱:LSEL,S,2,4 LSEL,A,15,18 图3-6 LSEL,A,63(2).选择附于线上的全部节点(图3-6)Utility Menu Select Entities Nodes Attached to lines,all 【OK】或叮嘱:NSLL,S,1 NPlot图图3-53-5图图3-63-6(3).设置单元属性 Main Menu Preprocessor Create Element Ele
17、m Attributes或叮嘱:Type,2 Mat,1 Real,1(4).创建目标单元 Main Menu Preprocessor Modeling Create Elements Surf/Contact Surf to Surf留意:假如基体单元是 2D 平面或 3D 实体单元“Tlab”无效,假如基体单元是 2D 单元,“Shape”无效。选【pick all】或叮嘱:ESURF(图3-7)图图3-73-7Step 6.创建接触单元(1)为接触面选线 Utility Menu Select Entities Lines 图3-8 By Num/pick From Full 【OK】
18、选线(图3-8)【OK】或叮嘱:LSEL,s,33,34 LSEL,a,43,44 (2)选择附于选定线上的全部节点图3-9 Utility Menu Select Entities Nodes Attached to Lines,all From Full 【OK】或叮嘱:NSLL,S,1 Nplot(图3-9)图图3-83-8图图3-93-9(3)设置单元属性 Main Menu Preprocessor Create Element Elem Attributes Element type number=3 ConTA 171 Material number=1 Real constan
19、t set number=1 【OK】或叮嘱:Type,3 Mat,1 Real,1(4)创建接触面单元图3-11 Main Menu Preprocessor Modeling Create Elements Surf/Contact Surf to Surf 图3-10 【OK】【PICK ALL】(图3-10)或叮嘱:ESURFStep 7.选择全部选项并画单元 Utility Menu Select Everything Utility Menu Plot Elements或叮嘱:Allsel,all Eplot(图3-11)Step 8.求解 /solu antype,static
20、nlgeom,on solc,on time,0.4 nsubst,20,500,10 outres,all,all dl,45,uy,-0.4 solve dl,45,uy,-0.55 solve finiStep 9.后处理4 应用接触向导创建接触对全部面面接触的单元选项和参数都可以通过接触向导来限制。运用接触向导进行接触分析便利快捷:自动定义单元类型和实常数 快速得到接触选项和参数 快速显示和反转接触法向运用接触向导必需首先对基体进行网格划分,否则不能激活接触向导。下面仍以弹性环装配接触问题为例。接受接触向导完成接触对的创建来说明接触向导的运用方法。Step 1.复原数据文件 Utili
21、ty Menu File Resume from 选 Snap.db 【OK】(图4-1)此数据库 Snap.db 包括此有限元模型的几何、基体单元、分网、材料和边界条件。但不包括接触单元与目标单元及接触有关的参数。Step 2.启动接触向导Main Menu Preprocessor Modeling Create Contact PairStep 3.创建接触对 图4-3 图4-4【Pick Target】(图43),拾取线(见图35)【OK】,再定义接触面 返回图44【Next】(图45)图45【opening setting】图46图46Step 4.设置接触参数【OK】返回图45【c
22、reate】【Finish】返回图42退出Step 5.求解求解过程同前图425 实例 2 平面密封圈计算 验证为接触刚度估计一个合适的初始值的重要性 验证接触分析中摩擦的应用 弹性库仑摩擦 验证基于接触单元临界状态变更的时间步预料如何会对计算效率有害模型描述:2D 超弹平面密封圈 对称,密封圈压缩模拟分析此模型接受两种方法:1.不运用接触向导,建立命流文件 Seal.inp2.运用接触向导创建接触对,并用 GUI 方式创建 Pilot 节点。本例具有:几何 N.L(大应变与大变形)材料 N.L(超弹)接触 N.L 数据文件:Seal.inpStep 1.Step 1.建模:单元类型、实常数、
23、建模:单元类型、实常数、材料特性、基体分网材料特性、基体分网/prep7et,1,56 !HYPER56 2D 4node U_P Hyperelastic Solidmp,nuxy,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.133k,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
24、$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$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!*基体分网*lesize,8,2$lesize,13,4$lesize,14,4esize,0.035type,1mat,1amesh,allsave,seal,dbStep 2.创建接触对,先定义单元和接触参数et,2,169 !Target169 2D Target Seg
25、mentet,3,171 !Conta171 2D Surface to Surface Contactmp,mu,1,0 !用材料特性定义摩擦,本例先无摩擦计算然后再有摩擦计算,比较结果r,1,0.1,0.1 !用实常数定义接触高级选项对于弯曲为主的问题,接受接触刚度FKN=0.1作为初始估计值,FTOLN=0.1 lagrange穿透比例系数(缺省值)r,2,0.1,0.1 !假如FKN=1.0则不收敛!*!*创建接触对创建接触对1 1(定义实常数和分网)(定义实常数和分网)*type,2 !目标面 1 Target169real,1 !Target169的实常数mat,1 !接触面 1
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- ANSYS 高级 接触 分析 优秀 PPT
![提示](https://www.deliwenku.com/images/bang_tan.gif)
限制150内