ANSYS高级接触问题.pptx
《ANSYS高级接触问题.pptx》由会员分享,可在线阅读,更多相关《ANSYS高级接触问题.pptx(190页珍藏版)》请在咨信网上搜索。
1、分析中常常需要确定两个或多个相互接触物体的位移、接触区域的大小和接触面上的应力分布。接触分析存在两大难点:在求解之前,你不知道接触区域的范围;表面之在求解之前,你不知道接触区域的范围;表面之间是接触还是分开是未知的;表面之间突然接触间是接触还是分开是未知的;表面之间突然接触或突然不接触会导致系统刚度的突然变化。或突然不接触会导致系统刚度的突然变化。大多数接触问题需要计算摩擦。摩擦是与路径有关的现象,摩擦响应还可能是杂乱的,使问题求解难以收敛。ANSYS 高级接触问题 第1页/共190页1 接触分类刚柔一个表面是完全刚性的除刚体运动外无应变、应力和变形,另一表面为软材料构成是可变形的。只在一个表
2、面特别刚硬并且不关心刚硬物体的应力时有效。柔柔两个接触体都可以变形。第2页/共190页2 接触单元ANSYS采用接触单元来模拟接触问题:跟踪接触位置;保证接触协调性(防止接触表面相互穿透);在接触表面之间传递接触应力(正压力和摩擦)。接触单元就是覆盖在分析模型接触面上的一层单元。在 ANSYS 中可以采用三种不同的单元来模拟接触:面一面接触单元;点一面接触单元;点一点接触单元。第3页/共190页不同的单元类型具有完全不同的单元特性和分析过程。1.面一面接触单元用于任意形状的两个表面接触 不必事先知道接触的准确位置;两个面可以具有不同的网格;支持大的相对滑动;支持大应变和大转动。例如:面一面接触
3、可以模拟金属成型,如轧制过程。2 接触单元第4页/共190页2 接触单元第5页/共190页2.点一面接触单元用于某一点和任意形状的面的接触 可使用多个点面接触单元模拟棱边和面的接触;不必事先知道接触的准确位置;两个面可以具有不同的网格;支持大的相对滑动;支持大应变和大转动。例:点面接触可以模拟棱边和面之间的接触2 接触单元第6页/共190页2 接触单元第7页/共190页3.点点接触单元用于模拟单点和另一个确定点之间的接触。建立模型时必须事先知道确切的接触位置;多个点点接触单元可以模拟两个具有多个单元表面间的接触;每个表面的网格必须是相同的;相对滑动必须很小;只对小的转动响应有效。例如:点一点接
4、触可以模拟一些面的接触。如地基和土壤的接触2 接触单元第8页/共190页2 接触单元第9页/共190页3 关于耦合和约束方程的应用如果接触模型没有摩擦,接触区域始终粘在一起,并且分析是小挠度、小转动问题,那么可以用耦合或约束方程代替接触。使用耦合或约束方程的优点是分析还是线性的 第10页/共190页接触问题的一般特性1 接触刚度1、所有的 ANSYS 接触单元都采用罚刚度(接触刚度)来保证接触界面的协调性第11页/共190页在数学上为保持平衡,需要有穿透值然而,物理接触实体是没有穿透的 分析者将面对困难的选择:小的穿透计算精度高,因此接触刚度应该大;然而,太大的接触刚度会产生收敛困难:模型可能
5、会振荡,接触表面互相跳开。接触刚度是同时影响计算精度和收敛的最重要的参数。你必须选定一个合适的接触刚度。除了在表面间传递法向压力外,接触单元还传递切向运动(摩擦)。采用切向罚刚度保证切向的协调性。(图12)作为初值,可采用:Ktangent=0.01 Knormal切向罚刚度与法向罚刚度以同样的方式对收敛性和计算精度产生影响。第12页/共190页2、接触刚度的选取 选定一个合适的接触刚度值需要一些经验。对于面一面接触单元,接触刚度通常指定为基体单元刚度的一个比例因子。开始估计时,选用 FKN=1.0 大面积实体接触 FKN=0.01-0.1 较柔软(弯曲占主导的部分)另外,也可以指定一个绝对刚
6、度值,单位:(力/长度)/面积。点一点(除 CONTA178)和点面接触单元需要为罚刚度 KN 输入绝对值:初始估计时:对于大变形:0.1*E KN 1.0*E 对于弯曲:0.01*E KN 0.1*E E 为弹性模量第13页/共190页3、选取接触刚度的指导:Step 1.开始采用较小的刚度值Step 2.对前几个子步进行计算Step 3.检查穿透量和每一个子步中的平衡迭代次数在粗略的检查中,如以实际比例显示整个模型时就能观察到穿透,则穿透可能太大了,需要提高刚度重新分析。如果收敛的迭代次数过多(或未收敛),降低刚度重新分析。注意:罚刚度可以在载荷步间改变,并且可以在重启动中调整。牢记:接触
7、刚度是同时影响计算精度和收敛性的最重要的参数。如果收敛有问题,减小刚度值,重新分析在敏感的分析中,还应该改变罚刚度来验证计算结果的有效性。在分析中减小刚度范围,直到结果(接触压力、最大SEQV 等)不再明显改变。第14页/共190页2 摩擦1、两个接触体的剪切或滑动行为可以是无摩擦的或有摩擦的 无摩擦时允许物体没有阻力地相互滑动;有摩擦时,物体之间会产生剪切力。2、摩擦消耗能量,并且是路径相关行为。为获得较高的精度,时间步长必须小(图2-1)第15页/共190页式中:摩擦系数 一旦所受剪力超过 FT,两物体将发生相对滑动。4、弹性库仑摩擦模型:允许粘着和滑动。图213、ANSYS 中,摩擦采用
8、库仑模型,并有附加选项可处理复杂的粘着和剪切行为。库仑法则是宏观模型,表述物体间的等效剪力 FT 不能超过正压力 FN 的一部分:FT Preprocessor Modeling Create Elements Surf/Contact Surf to Surf(ESURF)对于直接生成刚性目标面,在建立目标单元之前需要要指定附加的单元属性 TSHAP第24页/共190页第25页/共190页刚性目标面的自动划分不需要 TSHAP。ANSYS 能根据实体模型确定合适的目标单元形状。划分线(LMESH)2D 刚性目标面 划分面(AMESH)3D 刚性目标面 创建关键点(KMESH)控制节点(Pil
9、ot)刚性目标面能与控制点联系起来。Pilot 实际上是只有一个节点的单元,通过这个节点的运动可以控制整个目标面的运动。ANSYS 只在 Pilot 节点上检查边界条件而忽略其它节点的约束。对可变形体目标面建立目标单元的步骤是:1.先选择可变形体表面上的节点 2.然后在可变形体上建立目标单元 Main Menu Preprocessor Modeling Create Elements Surf/Contact Surf to Surf第26页/共190页ANSYS 将根据基体的网格确定目标单元形状和外法线方向。检查外法线方向(这在自动划分刚性目标面时非常重要)图 3-3 打开单元坐标系标志并
10、重绘单元 /PSYMS,ESYS,1 目标单元外法线方向应该指向接触面。如果单元法向不指向接触面,用命令使之反转:ESURF,REVE第27页/共190页例:Seal.dat(图3-3)Step 5.建立接触面单元 设置接触单元属性、选择可变形体表面节点,并在可变形体上建立接触单元(过程与在可变形体上建立目标单元相同)Main Menu Preprocessor Modeling Create Elements Surf/Contact Surf to Surf(ESURF)这些接触单元与基体有同样的阶数(低阶或高阶)。注意,在壳或梁单元上建立目标单元或接触单元时,可以选择要在梁或壳单元的顶层
11、还是底层建立单元。图图3-33-3第28页/共190页 在选择柔体表面上的节点时,如果你确定某一部分节点永远不会接触到目标面,可以忽略它,以减少计算时间。接触面的外法向应指向目标面。如果发现外法线方向不正确,用下列命令修改之 ESURF,REVEStep 6.在有限元模型上施加边界条件 如果目标面是刚性面,目标面将会自动固定。定义了 Pilot 点 ANSYS 只检查该点的边界条件,忽略目标面上其它节点约束。控制点能控制目标面的运动。对 Seal.dat 施加的边界条见图 33。第29页/共190页Step 7.定义求解选项和载荷步,以下是默认设置 推荐使用N.L求解自动控制 使用不带自适应下
12、降的 full Newton-Raphson 法求解 时间步必须足够小。使用自动时间步。子步数的最大值(NSBMX)应较大,最小值(NSBMIN)应较小Step 8.求解Step 9.后处理 结果包括位移、应力、应变和接触等信息。接触压力、摩擦应力、总应力、接触侵入、接触间隙距离、滑动距离和接触状态都可以从/POST1 或/POST26 中得到。第30页/共190页面一面接触分析实例(建立接触对不通过接触向导)实例1:弹性环装配第31页/共190页轴对称施加位移载荷:在 L45 线上施加 0.4的 Y 向位移打开几何 N.L 开关(NLGEOM,on)设 Time=0.4 并为自动时间步给出子
13、步数(20,500,10)给出输出控制(要求输出每一子步结果)求解并查看输出和监视文件重启动分析亦将 Y 向偏移量增加到 0.55使咬接装配的第 2 个齿咬合。问题以弯曲为主,设 FKN=0.1第32页/共190页GUI 方式:Step 1.恢复数据库文件 Snap.db(Snap.db 数据库包含此有限元模型的完整几何模型、材料、边界条件。但未定义目标单元与接触单元)。Utility Menu File Resume from 选 Snap.db 【OK】基体单元:Plane42,1 轴对称材料:MAT,1 EX=0.175e6;NUXY=0.35;MU=0.0 MAT,2 EX=90000
14、;NUXY=0.35 约束节点 163 UY=0 182 UY=0第33页/共190页Step 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 Preprocessor Element Typer Add/Edit/Delete【Add】contact 2nd Surf 171 Element t
15、yper reference number=3【OK】【close】或命令:ET,3,conta171第34页/共190页Step 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 tolerance FTOLN=0.1(不用此,不收敛)或命令:R,1,0.1,0.1第35页/共190页确定罚刚度 FKN
16、 值通常在 0.01-10 之间,对于体积变形问题用 1.0(默认),对弯曲问题用 0.1确定侵入容差:侵入容差(FTOLN)是与接触单元下面的实体单元深度(h)相乘的比例因子。若此值太小会引起收敛困难,绝对不要用太小的容差!增大罚刚度(FKN)将减少侵入。将 FKN 增大 100 倍会相应地减少侵入,但是接触压力只改变 5。如不收敛可调整 FKN 或 FTOLN 值重新运行。检查侵入和每个子步的平衡迭代数,如果收敛受侵入容差的驱使可能是FKN 值估计不足或 FTOLN 值太小。如果需要多次迭代才能使残值收敛而不是侵入。FKN 值可能估计得太高。第36页/共190页Step 5.创建目标单元(
17、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第37页/共190页(3).设置单元属性 Main Menu Preprocessor Create Element Elem Attributes或命令:Type,2 Mat,1
18、 Real,1第38页/共190页(4).创建目标单元 Main Menu Preprocessor Modeling Create Elements Surf/Contact Surf to Surf注意:如果基体单元是 2D 平面或 3D 实体单元“Tlab”无效,如果基体单元是 2D 单元,“Shape”无效。选【pick all】或命令:ESURF(图3-7)图图3-73-7第39页/共190页Step 6.创建接触单元(1)为接触面选线 Utility Menu Select Entities Lines 图3-8 By Num/pick From Full 【OK】选线(图3-8)
19、【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第40页/共190页(3)设置单元属性 Main Menu Preprocessor Create Element Elem Attributes Element type number=3 ConTA 171 Material number=1 Real const
20、ant 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)或命令:ESURF第41页/共190页Step 7.选择所有选项并画单元 Utility Menu Select Everything Utility Menu Plot Elements或命令:Allsel,all Eplot(图3-11)Step 8.求解 /solu an
21、type,static 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.后处理第42页/共190页第43页/共190页4 应用接触向导创建接触对所有面面接触的单元选项和参数都可以通过接触向导来控制。使用接触向导进行接触分析方便快捷:自动定义单元类型和实常数 快速得到接触选项和参数 快速显示和反转接触法向使用接触向导必须首先对基体进行网格划分,否则不能激活接触向导。下面仍以弹性环装配接触问题为例。采用接触向导完成接触对的创建
22、来说明接触向导的使用方法。Step 1.恢复数据文件 Utility Menu File Resume from 选 Snap.db 【OK】(图4-1)此数据库 Snap.db 包括此有限元模型的几何、基体单元、分网、材料和边界条件。但不包括接触单元与目标单元及接触有关的参数。第44页/共190页Step 2.启动接触向导Main Menu Preprocessor Modeling Create Contact Pair第45页/共190页Step 3.创建接触对 图4-3 图4-4【Pick Target】(图43),拾取线(见图35)【OK】,再定义接触面 返回图44【Next】(图4
23、5)第46页/共190页图45【opening setting】图46第47页/共190页图46第48页/共190页Step 4.设置接触参数【OK】返回图45【create】【Finish】返回图42第49页/共190页退出Step 5.求解求解过程同前图42第50页/共190页5 实例 2 平面密封圈计算 验证为接触刚度估计一个合适的初始值的重要性 验证接触分析中摩擦的应用 弹性库仑摩擦 验证基于接触单元临界状态变化的时间步预测如何会对计算效率有害模型描述:2D 超弹平面密封圈 对称,密封圈压缩模拟第51页/共190页分析此模型采用两种方法:1.不使用接触向导,建立命流文件 Seal.in
24、p2.使用接触向导创建接触对,并用 GUI 方式创建 Pilot 节点。本例具有:几何 N.L(大应变与大变形)材料 N.L(超弹)接触 N.L 数据文件:Seal.inp第52页/共190页Step 1.Step 1.建模:单元类型、实常数、建模:单元类型、实常数、材料特性、基体分网材料特性、基体分网/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
25、,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$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第53页/共190页!*基体分网*l
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- ANSYS 高级 接触 问题
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【快乐****生活】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【快乐****生活】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。