《ANSYS高级接触问题1372397实用教案》由会员分享,可在线阅读,更多相关《ANSYS高级接触问题1372397实用教案(124页珍藏版)》请在金锄头文库上搜索。
1、1接触(jich)分类刚柔一个表面是完全刚性的除刚体运动外无应变、应力和变形,另一表面为软材料构成是可变形的。只在一个表面特别刚硬并且不关心(gunxn)刚硬物体的应力时有效。柔柔两个接触体都可以变形。第1页/共123页第一页,共124页。2接触(jich)单元ANSYS采用接触单元来模拟接触问题:跟踪(gnzng)接触位置;保证接触协调性(防止接触表面相互穿透);在接触表面之间传递接触应力(正压力和摩擦)。接触单元就是覆盖在分析模型接触面上的一层单元。在ANSYS中可以采用三种不同的单元来模拟接触:面一面接触单元;点一面接触单元;点一点接触单元。第2页/共123页第二页,共124页。不同的单
2、元类型具有完全不同的单元特性和分析过程。1. 面一面接触单元用于任意形状的两个表面接触不必事先知道接触的准确位置;两个面可以具有不同的网格;支持大的相对滑动(hudng);支持大应变和大转动。例如: 面一面接触可以模拟金属成型,如轧制过程。第3页/共123页第三页,共124页。第4页/共123页第四页,共124页。2.点一面接触单元用于某一点和任意形状的面的接触可使用多个点面接触单元模拟棱边和面的接触;不必事先(shxin)知道接触的准确位置;两个面可以具有不同的网格;支持大的相对滑动;支持大应变和大转动。例:点面接触可以模拟棱边和面之间的接触第5页/共123页第五页,共124页。第6页/共1
3、23页第六页,共124页。3.点点接触单元用于模拟单点和另一个确定点之间的接触。建立模型时必须事先(shxin)知道确切的接触位置;多个点点接触单元可以模拟两个具有多个单元表面间的接触;每个表面的网格必须是相同的;相对滑动必须很小;只对小的转动响应有效。例如:点一点接触可以模拟一些面的接触。如地基和土壤的接触第7页/共123页第七页,共124页。第8页/共123页第八页,共124页。3关于(guny)耦合和约束方程的应用如果接触模型(mxng)没有摩擦,接触区域始终粘在一起,并且分析是小挠度、小转动问题,那么可以用耦合或约束方程代替接触。使用耦合或约束方程的优点是分析还是线性的第9页/共123
4、页第九页,共124页。接触问题的一般(ybn)特性1接触刚度1、所有的ANSYS接触单元都采用(ciyng)罚刚度(接触刚度)来保证接触界面的协调性第10页/共123页第十页,共124页。在数学上为保持平衡,需要有穿透值然而,物理接触实体是没有穿透的 分析者将面对困难的选择:小的穿透计算精度高,因此接触刚度应该大;然而,太大的接触刚度会产生收敛困难:模型可能会振荡,接触表面互相跳开。接触刚度是同时影响(yngxing)计算精度和收敛的最重要的参数。你必须选定一个合适的接触刚度。除了在表面间传递法向压力外,接触单元还传递切向运动(摩擦)。采用切向罚刚度保证切向的协调性。(图12)作为初值,可采用
5、:Ktangent=0.01Knormal切向罚刚度与法向罚刚度以同样的方式对收敛性和计算精度产生影响(yngxing)。第11页/共123页第十一页,共124页。2、接触刚度的选取选定一个合适的接触刚度值需要一些经验。对于(duy)面一面接触单元,接触刚度通常指定为基体单元刚度的一个比例因子。开始估计时,选用FKN=1.0大面积实体接触FKN=0.01-0.1较柔软(弯曲占主导的)部分另外,也可以指定一个绝对刚度值,单位:(力/长度)/面积。点一点(除CONTA178)和点面接触单元需要为罚刚度KN输入绝对值:初始估计时:对于(duy)大变形:0.1*EKN1.0*E对于(duy)弯曲:0.
6、01*EKN0.1*EE为弹性模量第12页/共123页第十二页,共124页。3、选取接触刚度的指导:Step1.开始采用较小的刚度值Step2.对前几个子步进行计算Step3.检查穿透量和每一个(y)子步中的平衡迭代次数在粗略的检查中,如以实际比例显示整个模型时就能观察到穿透,则穿透可能太大了,需要提高刚度重新分析。如果收敛的迭代次数过多(或未收敛),降低刚度重新分析。注意:罚刚度可以在载荷步间改变,并且可以在重启动中调整。牢记:接触刚度是同时影响计算精度和收敛性的最重要的参数。如果收敛有问题,减小刚度值,重新分析在敏感的分析中,还应该改变罚刚度来验证计算结果的有效性。在分析中减小刚度范围,直
7、到结果(接触压力、最大SEQV等)不再明显改变。第13页/共123页第十三页,共124页。2摩擦(mc)1、两个接触体的剪切或滑动行为可以是无摩擦的或有摩擦的无摩擦时允许物体没有阻力地相互滑动;有摩擦时,物体之间会产生剪切力。2、摩擦消耗能量,并且(bngqi)是路径相关行为。为获得较高的精度,时间步长必须小(图2-1)第14页/共123页第十四页,共124页。图213、ANSYS中,摩擦采用库仑模型,并有附加选项可处理复杂(fz)的粘着和剪切行为。库仑法则是宏观模型,表述物体间的等效剪力FT不能超过正压力FN的一部分:FTPreprocessorModelingCreateElementsS
8、urf/ContactSurftoSurf(ESURF)对于直接生成刚性目标面,在建立目标单元之前需要要指定附加的单元属性TSHAP第23页/共123页第二十三页,共124页。第24页/共123页第二十四页,共124页。刚性目标面的自动划分不需要TSHAP。ANSYS能根据实体模型确定合适的目标单元形状。划分线(LMESH) 2D刚性目标面划分面(AMESH) 3D刚性目标面创建关键点(KMESH)控制节点(Pilot)刚性目标面能与控制点联系起来,Pilot实际上是只有一个(y )节点的单元,通过这个节点的运动可以控制整个目标面的运动。ANSYS只在Pilot节点上检查边界条件而忽略其它节点
9、的约束。对可变形体目标面建立目标单元的步骤是:1.先选择可变形体表面上的节点2.然后在可变形体上建立目标单元Main Menu Preprocessor Modeling Create Elements Surf/Contact Surf to Surf第25页/共123页第二十五页,共124页。ANSYS将根据基体的网格确定目标单元形状和外法线方向。检查外法线方向(这在自动划分刚性目标面时非常重要)图3-3打开单元坐标系标志(biozh)并重绘单元/PSYMS,ESYS,1目标单元外法线方向应该指向接触面。如果单元法向不指向接触面,用命令使之反转:ESURF,REVE第26页/共123页第二
10、十六页,共124页。例:Seal.dat(图3-3)Step5.建立接触面单元设置接触单元属性、选择可变形体表面节点,并在可变形体上建立接触单元(过程与在可变形体上建立目标单元相同)MainMenuPreprocessorModelingCreateElementsSurf/ContactSurftoSurf这些接触单元与基体有同样的阶数(低阶或高阶)。注意,在壳或梁单元上建立目标单元或接触单元时,可以选择要在梁或壳单元的顶层(dncn)还是底层建立单元。 图3-33-3第27页/共123页第二十七页,共124页。在选择柔体表面上的节点时,如果你确定某一部分节点永远不会接触到目标面时,可以忽略
11、它,以减少计算(j sun)时间。接触面的外法向应指向目标面。如果发现外法线方向不正确,用下列命令修改之ESURF,REVEStep 6.在有限元模型上施加边界条件如果目标面是刚性面,目标面将会自动固定。定义了Pilot点ANSYS只检查该点的边界条件,忽略目标面上其它节点约束。控制点能控制目标面的运动。对Seal.dat施加的边界条见图33。第28页/共123页第二十八页,共124页。Step7.定义求解选项和载荷步,以下是默认设置推荐使用N.L求解自动控制使用不带自适应下降的fullNewton-Raphson法求解时间步必须足够小。使用自动时间步。子步数的最大值(NSBMX)应较大,最小
12、值(NSBMIN)应较小Step8.求解Step9.后处理结果包括位移、应力(yngl)、应变和接触等信息。接触压力、摩擦应力(yngl)、总应力(yngl)、接触侵入、接触间隙距离、滑动距离和接触状态都可以从/POST1或/POST26中得到。第29页/共123页第二十九页,共124页。面一面接触分析实例(建立接触对不通过接触向导)实例1:弹性(tnxng)环装配第30页/共123页第三十页,共124页。轴对称施加位移载荷:在L45线上施加0.4的-Y向位移打开几何N.L开关(kigun)(NLGEOM,on)设Time=0.4并为自动时间步给出子步数(20,500,10)给出输出控制(要求
13、输出每一子步结果)求解并查看输出和监视文件重启动分析亦将-Y向偏移量增加到0.55使咬接装配的第2个齿咬合。问题以弯曲为主,设FKN=0.1第31页/共123页第三十一页,共124页。GUI方式:Step1.恢复数据库文件(wnjin)Snap.db(Snap.db数据库包含此有限元模型的完整几何模型、材料、边界条件。但未定义目标单元与接触单元)。UtilityMenuFileResumefrom选Snap.db【OK】基体单元:Plane42,1轴对称材料:MAT,1EX=0.175e6;UXY=0.35;MU=0.0MAT,2EX=90000;NUXY=0.35约束节点163UY=0182
14、UY=0第32页/共123页第三十二页,共124页。Step2.创建目标面单元(dnyun)类型MainMenuPreprocessorElementTyperAdd/Edit/Delete【Add】contact2Dtarget169Elementtyperreferencenumber=2【OK】或命令:/PREP7ET,2,Target169Step3.建立接触面单元(dnyun)类型MainMenuPreprocessorElementTyperAdd/Edit/Delete【Add】contact2ndSurf171Elementtyperreferencenumber=3【OK】【
15、close】或命令:ET,3,conta171第33页/共123页第三十三页,共124页。Step4.指定接触法向刚度MainMenuPreprocessorRealConstantsAdd/Edit/Delete【Add】Type3CONTA171【OK】RealConstantSetNo.=1NormalPenaltyStiffnessFKN=0.1(对弯曲问题采用初值0.1)PenetrationtoleranceFTOLN=0.1(不用(byng)此,不收敛)或命令:R,1,0.1,0.1第34页/共123页第三十四页,共124页。确定罚刚度FKN值通常在0.01-10之间,对于体积变
16、形问题用1.0(默认),对弯曲(wnq)问题用0.1确定侵入容差:侵入容差(FTOLN)是与接触单元下面的实体单元深度(h)相乘的比例因子。若此值太小会引起收敛困难,绝对不要用太小的容差!增大罚刚度(FKN)将减少侵入。将FKN增大100倍会相应地减少侵入,但是接触压力只改变5。如不收敛可调整FKN或FTOLN值重新运行。检查侵入和每个子步的平衡迭代数,如果收敛受入侵容差的驱使可能是FKN值估计不足或FTOLN值太小。如果需要多次迭代才能使残值收敛而不是侵入。FKN值可能估计得太高。第35页/共123页第三十五页,共124页。Step5.创建目标单元(1).为目标面选择(xunz)线Utili
17、tyMenuSelectEntities拾取线(图3-5)【OK】或命令:LSEL,S,2,4LSEL,A,15,18图3-6LSEL,A,63(2).选择(xunz)附于线上的全部节点(图3-6)UtilityMenuSelectEntitiesNodesAttachedtolines,all【OK】或命令:NSLL,S,1NPlot图3-53-5图3-63-6第36页/共123页第三十六页,共124页。(3).设置(shzh)单元属性MainMenuPreprocessorCreateElementElemAttributes或命令:Type,2Mat,1Real,1第37页/共123页第
18、三十七页,共124页。(4).创建目标单元MainMenuPreprocessorModelingCreateElementsSurf/ContactSurftoSurf注意(zhy):如果基体单元是2D平面或3D实体单元“Tlab”无效,如果基体单元是2D单元,“Shape”无效。选【pickall】或命令:ESURF(图3-7)图3-73-7第38页/共123页第三十八页,共124页。Step6.创建接触单元(1).为接触面选线UtilityMenuSelectEntitiesLines图3-8ByNum/pickFromFull【OK】选线(图3-8)【OK】或命令:LSEL,s,33,
19、34LSEL,a,43,44(2).选择(xunz)附于选定线上的全部节点图3-9UtilityMenuSelectEntitiesNodesAttachedtoLines,allFromFull【OK】或命令:NSLL,S,1Nplot(图3-9)图3-83-8图3-93-9第39页/共123页第三十九页,共124页。(3).设置(shzh)单元属性MainMenuPreprocessorCreateElementElemAttributesElementtypenumber=3ConTA171Materialnumber=1Realconstantsetnumber=1【OK】或命令:Ty
20、pe,3Mat,1Real,1(4).创建接触面单元图3-11MainMenuPreprocessorModelingCreateElementsSurf/ContactSurftoSurf图3-10【OK】【PICKALL】(图3-10)或命令:ESURF图3-103-10图3-113-11第40页/共123页第四十页,共124页。Step7.选择所有选项并画单元(dnyun)UtilityMenuSelectEverythingUtilityMenuPlotElements或命令:Allsel,allEplot(图3-11)Step8.求解/soluantype,staticnlgeom,
21、onsolc,ontime,0.4nsubst,20,500,10outres,all,alldl,45,uy,-0.4solvedl,45,uy,-0.55solvefiniStep9.后处理第41页/共123页第四十一页,共124页。第42页/共123页第四十二页,共124页。4应用接触(jich)向导创建接触(jich)对所有面面接触的单元选项和参数都可以通过接触向导来控制。使用接触向导进行接触分析方便快捷:自动定义单元类型和实常数快速(kuis)得到接触选项和参数快速(kuis)显示和反转接触法向使用接触向导必须首先对基体进行网格划分,否则不能激活接触向导。下面仍以弹性环装配接触问题为
22、例。采用接触向导完成接触对的创建来说明接触向导的使用方法。Step1.恢复数据文件UtilityMenuFileResumefrom选Snap.db【OK】(图4-1)此数据库Snap.db包括此有限元模型的几何、基体单元、分网、材料和边界条件。但不包括接触单元与目标单元及接触有关的参数。图4-14-1第43页/共123页第四十三页,共124页。Step2.启动(qdng)接触向导MainMenuPreprocessorModelingCreateContactPair第44页/共123页第四十四页,共124页。Step3.创建(chungjin)接触对图4-3图4-4【PickTarget】
23、(图43),拾取(shq)线(见图35)【OK】,再定义接触面返回图44【Next】(图45)第45页/共123页第四十五页,共124页。图45【opening setting】图46第46页/共123页第四十六页,共124页。图46第47页/共123页第四十七页,共124页。Step4.设置接触参数(cnsh)【OK】返回图45【create】【Finish】返回(fnhu)图42第48页/共123页第四十八页,共124页。退出(tuch)Step5.求解求解过程同前图42第49页/共123页第四十九页,共124页。5实例2平面(pngmin)密封圈计算验证为接触刚度估计一个合适的初始值的重
24、要性验证接触分析中摩擦的应用弹性(tnxng)库仑摩擦验证基于接触单元临界状态变化的时间步预测如何会对计算效率有害模型描述:2D超弹平面密封圈1/2对称,密封圈压缩模拟第50页/共123页第五十页,共124页。分析此模采用两种方法:1.不使用接触向导,建立(jinl)命流文件Seal.inp2.使用接触向导创建接触对,并用GUI方式创建Pilot节点。本例具有:几何N.L(大应变与大变形)材料N.L(超弹)接触N.L数据文件:Seal.inp第51页/共123页第五十一页,共124页。Step1.建模:单元类型(lixng)、实常数、材料特性、基体分网/prep7et,1,56!HYPER56
25、2D4nodeU_PHyperelasticSolidmp,nuxy,1,0.49tb,mooney,1tbdata,1,80!C10=80tbdata,2,20!C01=20k,1$k,2,0.333,0$k,3,0.867,0.867k,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命令(mnglng)重复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,
26、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第52页/共123页第五十二页,共124页。!*基体分网*lesize,8,2$lesize,13,4$lesize,14,4esize,0.035type,1mat,1amesh,allsave,seal,dbStep2.创建接触对,先定义单元和接触参数et,2,169!Target1692DTargetSegmentet,3,171!Conta1712DSurfacetoSurfac
27、eContactmp,mu,1,0!用材料特性定义摩擦,本例先无摩擦计算然后再有摩擦计算,比较结果r,1,0.1,0.1! 用实常数定义接触高级选项对于弯曲为主的问题,采用接触刚度FKN=0.1作为(zuwi)初始估计值,FTOLN=0.1lagrange穿透比例系数(缺省值)r,2,0.1,0.1! 如果FKN=1.0则不收敛第53页/共123页第五十三页,共124页。!*创建接触对1(定义(dngy)实常数和分网)*type,2!目标面1Target169real,1!Target169 的实常数mat,1!接触面1定义材料与超弹单元同lesize,15,1!目标面1(L15)分网尺寸(图
28、53)lmesh,15!目标面1(L15)分网,采用自动分网,此刚性面自动约束。检查外法线方向【OK】!*接触面1:实常数(与目标面1一致)、分网lsel,s,line,1,3,1!选L1,L2,L3lsel,a,line,9,10,1!再加L9,L10lsel,a,line,14!再加L14nsll,s,1!选择(xunz)和所选线相连系的节点type,3!Conta171这里未发布Real,1,实常数未重新定义就是前面的Real,1,MAT也是前面的MAT,1esurf!生成接触单元lsel,all!*检查外法线方向【OK】*图5 53 3第54页/共123页第五十四页,共124页。!*创
29、建接触对2*type,2!目标面2Target169real,2!MAT未重新(chngxn)定义,就是前面的MAT,1lesize,16,1!目标面2lmesh,16!L16分网(目标2分网)kmesh,27!指定Pilot点!*检查外法线方向【OK】*type,3!接触面2Conta171Real,2;Mat,1lsel,s,line,5,7,1!选L5,L6,L7lsel,a,line,11,13,1!再加L11,L12,L13nsll,s,1esurf!生成接触单元,其Real,2;Mat,1lsel,allnsel,all!*检查外法线方向【OK】*第55页/共123页第五十五页,共
30、124页。Step3.施加(shji)边界条件!刚性面1被约束!刚性面2随Pilot点移动!L4对称约束dl,4,ux,0n_load=node(kx(27),ky(27),0)!为定义刚性面2的位移做准备finishStep4.定义求解选项和载荷步/solunlgeom,onsolc,on!N.L求解自动控制打开(缺省)time,0.85!载荷步、结束时间d,n_load,uy,-0.85nsubst,25,500,10outres,all,allmonitor,var3,n_load,fy第56页/共123页第五十六页,共124页。Step5.求解solveStep6.查看结果/post1
31、pldsp,2!变形(binxng)图plnsol,s,eqv,0,1!VonMises 应力云图plnsol,cont,pres,0,!接触压力图plnsol,epto,eqv!绘等效总应变图54save,seal,db第57页/共123页第五十七页,共124页。第58页/共123页第五十八页,共124页。定义(dngy)反力变量,绘载荷变形图UtilityMenuPlotElementsMainMenuTimeHistPostproDefineVariables【Add】Reactionforce【OK】拾取(shq)节点263(控制点)【OK】User-SpecifiedLoad=FOR
32、CE:FY【OK】【close】或命令:/POST26RFORCE,2,263,F,Y,FORCEMainMenuTimeHistPostproGraphVariables1stVariabletogragh=2【OK】或命令:PLVAR,2( 图54)图5 54 4第59页/共123页第五十九页,共124页。2、用接触向导(xingdo)创建接触对,用GUI方式创建Pilot节点(刚性面控制点),然后求解。Step1.恢复数据库文件Seal.db(包括基体的几何、单元、分网;没有选接接触单元与目标单元,未定义接触对)UtilityMenuFileResumefrom选Seal.db【OK】S
33、tep2.启动接触向导(xingdo)MainMenuPreprocessorModelingCreateContactPair图42【contactWizard】第60页/共123页第六十页,共124页。创建(chungjin)接触对1图55第61页/共123页第六十一页,共124页。选线(1)(接触( jich)对1的目标面)【OK】返回图55【next】 图56选线(2)(接触( jich)对1的接触( jich)面)【OK】;返回图56【next】第62页/共123页第六十二页,共124页。Step3.设置(shzh)接触参数图5 57 7第63页/共123页第六十三页,共124页。图
34、58图59定义摩擦(mc)(图59)【OK】返回图57 Coefficient friction 0.2 【Create】第64页/共123页第六十四页,共124页。511图512接触( jich)对1(图511)图510图第65页/共123页第六十五页,共124页。Step5.创建目标单元控制点(1).设置(shzh)单元属性用命令流:Type,4MAT,1REAL,4TSHAP,PIL10GUI:MainMenuPreprocessorModelingCreateElementElemAttributes图514第66页/共123页第六十六页,共124页。(2).绘关键点UtilityMe
35、nuPlotKeypointsKeypoints(3).建立目标单元控制点MainMenuPreprocessorMeshingMeshKeypoints选K27【OK】或命令流:KMESH,27Step6.施加位移(wiy)约束控制点27上施加UY=-0.85MainMenuSolutionDefineLoadsApplyStructuralDisplacementOnKeypoints选K27【OK】Lab2DOFstobeconstrained:UYValue:-0.85【OK】第67页/共123页第六十七页,共124页。施加对称(duchn)边界条件:MainMenuSolutionD
36、efineLoadsApplyStructuralDisplacementSymmetryB.C.OnLines选线L4【OK】Step7.求解控制、求解MainMenuSolutionAnalysisTypeSolnControls第68页/共123页第六十八页,共124页。MainMenuSolutionSolveLS_CurrentStep8.后处理同1无摩擦(MU=0)VonMises=145.096有摩擦(MU=0.2)VonMises=142.038接触刚度(nd)取1.0不收敛0.1收敛/POST26MainMenuTimeHistPostpro选控制点27,绘图(见图5-4)。
37、第69页/共123页第六十九页,共124页。面一面接触具有( jyu)20个可用实常数,2个材料属性和30个可用单元选项设置。能够模拟特殊的效果和处理困难的收敛情况。通常分析开始先用缺省值,只指定罚刚度和穿透容差及子步数。只能在缺省设置遇到困难时才用高级选项。所有的单元选项和参数都可以通过接触向导来控制;也可以通过实常数和单元选项来指定。Conta178接触单元、实常数选项见图6-1及6-2。图6-1 图6-26 面一面接触单元(dnyun)高级接触选项第70页/共123页第七十页,共124页。进入接触选项菜单用下列(xili)过程:MainMenuPreprocessorModelingCr
38、eateContactPair选【property】1、第1个重要选项:如何保证接触协调性:使用修正的拉朗日法(缺省)纯罚函数方法在接触向导中BasicContactalgorithm:AugmentedLagrangeMethodPenaltyMethod对大多数模型,修正的拉格朗日方法能很好地工作。罚函数法推荐用于具有(jyu)变形很大的单元,很大的摩擦系数和/或用修正的拉格朗日方法时收敛性很差的情况。第71页/共123页第七十一页,共124页。最重要(zhngyo)的选项是法向罚刚度或接触刚度。对于大面积接触采用起始值(因子)1.0,对于柔性接触采用0.1。大值对应较高精度;小值对应较好
39、的收敛性。有时最好先以较小的接触刚度进行分析,然后在一系列载荷步中逐渐增大刚度“渐变”的接触刚度,提高收敛性图6-2在最后的载荷步逐渐提高到一个刚硬的值将提高计算精度。可通过向导的基本表或单元特性菜单设置(shzh)以允许程序更新接触刚度接触向导BasicNormalPenaltyStiffness1.0factorconstant向导Basiccontactstiffnessupdate:Eachloadstep(PAIRIDbased)noneEachloadstep(允许用户指定刚度变化)Eachsubstep(允许自动和用户指定变化)第72页/共123页第七十二页,共124页。2、第2
40、个重要(zhngyo)选项穿透容差。缺省情况下,穿透容差是一个因子乘以基体单元厚度。对于变化很大的网格密度,采用因子会在接触表面的某些部分产生太小的容差,这时采用绝对值可能更好。不要使用太小的容差,因为(ynwi)它总是对收敛性有害。第73页/共123页第七十三页,共124页。3、对于临界(lnji)接触状态变化的自动时间步控制不控制:不影响自动时间步,对静力问题自动时间步打开时此选项一般是足够的。自动二分:如果接触状态变化明显,时间步长将二分。对于动力问题自动二分通常是足够的。合理值:比自动细分更耗时的算法。最小值:此选项为下一步预测最小时间增量(很耗计算时间,不推荐)或在单元(dnyun)
41、选项中控制:Contacttime/loadpredictionK7:Noprediction第74页/共123页第七十四页,共124页。4、Pinball区域影响接触状态的确定和其它(qt)许多接触特性。Pinball区域是环绕接触单元的园(2D)或球(3D),描述接触单元周围“远”和“近”区域的边界(图6-3)。在缺省情况(qngkung)下,Pinball区域半径是4基体单元厚度(刚柔)或2基体单元厚度(柔一柔)可以为Pinball半径指定一个不同的值。图6-3第75页/共123页第七十五页,共124页。也可用实常数PINB调整球形区(对于初始值侵入大的问题(wnt)是必要的)。第76页
42、/共123页第七十六页,共124页。5、几种不同的接触模式这些选项使你能够模拟(mn)特殊的物理现象。或单元选项:BehaviorcontactsurfaceK12:standard这些选项包括:标准:正常的接触闭合和打开行为,具有正常的粘着/滑动摩擦行为。粗糙:正常接触闭合和打开行为,但不发生滑动(类似于具有无限摩擦系数)不分离(滑动):目标(mbio)面和接触面一旦建立接触就不再分离(允许滑动)绑定:目标(mbio)面和接触面一旦接触就粘在一起不分离(永远):初始位于Pinball区域内或已经接触的接触检查点在法向不分离(允许滑动)绑定接触(永远):初始位于Pinball区域内或已经接触的
43、接触检查点在剩余的分析过程中绑定在一起(DesignSpace缺省值)绑定接触(初始接触):只在初始接触的地方采用绑定,初始分开的地方保持分开。计算实例:悬臂梁端部旋转(不分离行为)第77页/共123页第七十七页,共124页。6、影响(yngxing)某些表面行为的选项ContactOpeningstiffness(分开时的间隙刚度)保证不分离和绑定行为(xngwi),它通过使用当存在间隙也具有非零刚度的弹簧来连接表面。缺省情况下,此弹簧刚度等于法向罚刚度,其效果类似于法向罚刚度刚度太小精度低;刚度太大会引起收敛问题。第78页/共123页第七十八页,共124页。摩擦(mc)系数影响基本摩擦(m
44、c)行为:Contactcohesion表示当没有法向压力时开始滑动的摩擦应力值。摩擦导致非对称刚度阵。因为非对称矩阵很难计算(因此导致求解变慢),程序自动控制执行对称求解,利用此算法可以解决多数含摩擦接触问题。有时,采用非对称矩阵能获用更好的收敛性。如果遇到收敛缓慢问题可以用不对称求解选项。记住:这种情况必须使用稀疏或波前求解器。对于每个支持非对称矩阵的单元,此选项也可以由下列(xili)菜单激活:MainMenuSolutionUnabridgedMenuAnalysisOptions设置Newton-Raphson选项为FullN_Runsymm第79页/共123页第七十九页,共124页
45、。7、初始(chsh)穿透有几种技术可以(ky)模拟初始穿透接触问题(如过盈装置配)。可以(ky)使用初始几何穿透,或指定偏移量,或二者皆有。图6-5指定偏移量(CNOF)或在实常数(chngsh)中指定偏移量(CNOF)contact surface offset CNOF: 0.025 第80页/共123页第八十页,共124页。正的CNOF加大初始穿透负的CNOF减小初始穿透或导致间隙CNOF可与几何穿透组合自动(zdng)CNOF调整允许ANSYS基于初始穿透自动(zdng)给定CNOF值。导致“刚好接触”配置ICONT缺省为0或单元(dnyun)选项:AutoCNOF/ICONTadj
46、ustmentK5:No.Auto.Adjust第81页/共123页第八十一页,共124页。初始穿透选项包括:Include everything:包括由几何模型和指定偏移量(如果有的话)引起的初始穿透Exclude everything:忽略所有初始穿透效应。Include with ramped effects:渐变(jinbin)初始穿透,以提高收敛性。Include offset only:只包括由偏移量指定的基本初始穿透。Include offset only w/ramp:只包括由偏移量指定的基本初始穿透,且渐变(jinbin)初始穿透以提高收敛性第82页/共123页第八十二页,共
47、124页。如果模型包含初始几何穿透,接触力将立即“阶跃”到一个大值。载荷突变经常导致收敛困难,期望有一种机制(jzh)能够将初始穿透效应渐变到零。Includewithrampedeffects和Includeoffsetonlyw/ramp选项通过在第一载荷步,将初始穿透渐变为零克服收敛困难。为求得好的结果,在第一载荷步不应施加其它载荷(图6-6)。计算实例:初始穿透 图6-6第83页/共123页第八十三页,共124页。8、刚体模式初始不接触的两个(或多个)物体的静力分析中,在接触建立前可能(knng)产生刚体运动(图6-7)。此例中圆柱体没有施加位移约束,面由力控制。圆柱体的约束由圆柱体和
48、平板之间的接触建立。求解过程中两个物体分离,刚度矩阵奇异。ANSYS将产生一个负主元警告。有几个选项可以解决由于初始不相连物体引起的刚体模式: 图6-7第84页/共123页第八十四页,共124页。三个高级接触特性允许调整初始接触条件以防止刚体模式:(1).自动CONF调整程序计算( jsun)CNOF以清除间隙。(2).初始接触环(ICONT)将调整带内接触表面上的节点移到目标面上(3).初始允许穿透范围(PMIN&PMAX):将刚体表面移动到接触面上调整初始接触条件(ICONT)(图6-11):实常数ICONT可用于指定目标面上的“调整环”。位于调整环内的任何接触点都要移动到目标面上。推荐只
49、进行很小的修正,如果ICONT值太大会产生不连续。如果未指定常数ICONT,ANSYS根据模型尺寸为ICONT提供一个小的缺省值。关闭ICONT,必须将其设置成非常小的值(如1e-20)0值代表非0的缺省值计算( jsun)实例:刚体ICONT图6-11第85页/共123页第八十五页,共124页。如果目标面有0约束,采用PMIN和PMAX的初始调整将不被执行。初始调整是一个迭代过程,ANSYS最多使用20个迭代步把目标面调整到PMIN和PMAX范围内。当目标面被移动,接触体之间不再存在间隙,而成为闭合(b h)接触的初始状态。计算实例:刚体PMIN & PMAX图6-12第86页/共123页第
50、八十六页,共124页。7几种典型(dinxng)接触问题分析实例例2.悬臂梁端部旋转目标:采用“不分离”行为选项施加大转动模拟特定物理效果。不分离:目标面和接触(jich)面一旦建立接触(jich)就不再分离(允许滑动)。模型描述:悬臂梁,端部施加旋转(图7-4)文件(wnjin):no_separation.inp图7-4第87页/共123页第八十七页,共124页。采用Solid185单元:3D8nodestructuralSolid单元,具有大变形能力。/view,1,1,1/prep7k,1,10,0,0BLOCK,0,10,-0.5,0.5,-0.5,0.5RECTNG,9.5,10,
51、-1,1AGEN,7,-0.5,1AGEN,ITIEM,NA1,NA2,NINC,DX,DY,DZ,KING,MOELEM,IMOVEIMOVE=1moveoriginalareastonewpositionARSYM,Z,7ARSYM,Ncomp,NA1复制一组面积,对称于某轴et,1,185mp,ex,1,10000esize,0.5vmesh,1ESIZE,2!为目标面分网用/solutionnsel,loc,x,0d,all,all,0/pbc,u,1finisheplotSAVE,no-separation,db1以下用接触向导(xingdo)建立接触对:Step1.恢复no_sep
52、aration.dbUtilityMenuFileResumefrom选no_separation.db【OK】第88页/共123页第八十八页,共124页。Step2.采用接触向导创建接触对MainMenuPreprocessorModelingCreateContactPair【contactwizard】指定目标面:Targetsurface:AreasTargetType:Rigid【Picktarget】选A7,A8【OK】指定接触面:【next】contactsurface:Areas【Picktarget】选A1,A2【OK】指定选项设置(shzh):【next】【Optional
53、Settings】选择表BASICcontactsurfaceBehavior=NoSeparation【OK】建立接触对:CoefficientofFriction=0【Create】【Finish】第89页/共123页第八十九页,共124页。Step3.创建目标(mbio)单元控制点设置单元属性:MainMenuPreprocessorModelingCreateElementElemAttributesElementtypenumber=2Targe170Materialnumber=1Realconstantsetnumber=3TargetElementShape=Pilotnode
54、或用命令:Type,2MAT,1REAL,3TSHAP,PILO绘关键点:Kplot选K1(节点198)第90页/共123页第九十页,共124页。建立目标单元控制点MainMenuPreprocessorMeshingMeshKeypoints拾取(shq)关键点1【OK】或命令:KMESH,1Step4.施加45转动(0.7854弧度)MainMenuPreprocessorLoadsDefineLoadsApplyStructuralDisplacementOnKeypointsROTXDisplacementValue=0.7854【OK】或命令:/SoluDK,1,ROTX,0.785
55、4Step5.指定“基本”求解选项MainMenuSolutionAnalysisTypeSolnControls选BASICAnalysisOption=LargeDisplacementStaticTimeatendofloadstep=45NumberofSubsteps=45MaxNo.ofSubsteps=180MinNo.ofSubsteps=5AllSolutionItemsWriteeverySubstep【OK】或命令:NLGOM,onTime,45NSUBST,45,180,5Outres,all,all第91页/共123页第九十一页,共124页。Step6.求解(qiji
56、)SolveStep7.读入最后子步结果,后处理/POST1set,last选实体单元进行后处理:UtilityMenuSelectEntitiesElementsByElemNameElementname=solid185Choose=FromFull【OK】或命令:ESEL,S,Ename,Solid185绘位移总矢量:MainMenuGeneralPostprocPlotResultsNodalSoluDOFSolutionUSUM【OK】或命令:PLNSOL,U,SUMUSUMDMX=0.855708SMX=0.855708第92页/共123页第九十二页,共124页。第93页/共123
57、页第九十三页,共124页。例3.紧配合环(初试(chsh)穿透)目标:验证接触初始穿透选项。采用缺省值和渐进化穿透选项求解一个大的初始穿透接触分析,比较结果。模型模述:紧配合环、轴对称(图7-5)文件:Interference.inp图7-5第94页/共123页第九十四页,共124页。/PREP7ET,1,183,1!轴对称大应变ET,2,169!Target1692D 目标单元ET,3,172!Target1722D3 节点( jidin)高阶抛物线接触单元MP,EX,1,10E3MP,NUXY,1,0MP,MU,2,0! 无摩擦RECTNG,1.5,2,0,1RECTNG,1.9,2.4,
58、0,1ESIZE,0.05AMESH,ALL!Type,1;Mat,1TYPE,2MAT,2REAL,2LSEL,LOC,X,1.9NSLL,1!1.9线上的节点( jidin)ESURF!创建目标单元Target169!*TYPE,3!Mat,2;Real,2LSEL,LOC,X,2NSLL,S,1ESURF!创建接触面conta172ALLSEL,ALL第95页/共123页第九十五页,共124页。/SOLUTIOND,NODE(1.5,0,0),UY,0D,NODE(2.4,0,0),UY,0NLGEOM,ONNSUBST,10,50,5FINISH/PBC,U,1/NUMBER,1/PN
59、UM,MAT,1EPLOTSAVE,interference,db!包含有限元模型、接触(jich)对、材料、边界条件,求解控制第96页/共123页第九十六页,共124页。分析步骤:Step1.恢复数据库文件UtilityMenuFileResumefrom选interference.db【OK】Step2.重定位输出文件UtilityMenuFileSwitchOutputtoFileInterference.output【OK】或命令:/Output,interference,OutputStep3.求解接触分析/SolusolveStep4.重定位输出到输出窗口UtilityMenuFi
60、leSwitchOutputtoOutputWindow或命令:/OUTPUT,TERMStep5.绘制(huzh)径向应力MainMenuGeneralPostprocPlotResultsContourPlotNodalSoluStressX-directionSX【OK】第97页/共123页第九十七页,共124页。或命令:/POST1PLNSOL,S,X!结果:径向应力0Step6.列表显示(xinsh)接触结果MainMenuGeneralPostprocListResultsNodalSolution选NonlinearitemscontactCONT【OK】或命令:PRNSOL,C
61、ONTSTAT=0“ 张开”几乎不接触STAT=1“ 张开”但几乎接触STAT=2“ 闭合”并滑动STAT=3“ 闭合”并粘着第98页/共123页第九十八页,共124页。Step7.检查输出文件内容(interference.output)*NOTE*Nocontactwasdetectedforacontactpairspecifiedbyrealconstantset2大初始穿透问题(wnt)可能导致无法检测到接触或检测到接触面求解存在收敛困难。在此情况下推荐在第一个载荷步中包含渐进化效应。Step8.为初始穿透选项指定渐进化效应可通过conta172单元选项设置可通过接触向导中Initi
62、alAdjustment的设置(1)MainMenuPreprocessorElementTypeAdd/Edit/Delete选Type3conta172【Option】initialpenetration/gapK9=include_ramp【OK】【close】或命令:/prep7KEYOPT,3,9,2第99页/共123页第九十九页,共124页。(2)或向导(xingdo)InitialAdjustmentInitialpenetration:Includeeverythingwithrampedeffect【OK】第100页/共123页第一百页,共124页。Step9.重定位输出到文
63、件(wnjin)UtilityMenuFileSwitchOutputtoFileRamped.output【OK】Step10.分析/SoluSolveStep11.重定位输出到输出窗口UtilityMenuFileSwitchOutputtoOutputWindow命令:/output,TERMStep12.绘径向应力MainMenuGeneralPostprocPlotResultsContourPlotNodalSoluStressX-directionSX【OK】或命令:/POST1PLNSOL,S,X第101页/共123页第一百零一页,共124页。Step13.将轴对称模型扩展为整
64、个模型显示UtilityMenuPlotCtrlsStyleSymmetryExpansion2DAxi-SymmetricSelectexpansionamount=Fullexpansion【OK】或命令:/EXPAND,36,Axis,10/REPlotISO显示/View,1,1,1/ReplotStep14.检查输出文件ramped.outputInitialpenetrationwillberampedduringthefirstloadstep在第一载荷步,初始穿透渐进化( jnhu)。*NOTE*Max.Initialpenetration0.1wasdetectedbetwe
65、encontactelement425andtargetelement405specifiedbyreadconstantset2检测到大的初始穿透第102页/共123页第一百零二页,共124页。例4.两刚性面间夹圆柱体分析(刚体模式应用初始接触环ICONT)目标:验证使用初始接触环调整带来克服(kf)接触分析中的刚体模式模型描述:2D平面应力、1/4对称、施加压力(图7-6)文件名:ICONT.inp/prep7CYL4,0,0,0,0,1,90RECTNG,0,2.5,1.003,1.5!图7-7ET,1,182!Plane1822DStructuralSolid4节点UX,UY平面应力、
66、平面应变和轴对称单元、大变形、大应变图7-6 图7-7第103页/共123页第一百零三页,共124页。ET,2,169!Target169ET,3,171!Conta171MP,EX,1,1000MP,EX,2,30E6MP,MU,3,0ESIZE,0.04AMESH,1!type,1;mat,1MAT,2AMESH,2!*创建目标(mbio)目单元*TYPE,2MAT,3REAL,2LSEL,LOC,Y,1.003NSLL,1ESURF!*创建接触单元*TYPE,3CSYS,1LSEL,LOC,X,1NSLL,1ESURFCSYS,0LSEL,ALL/SOLUTIONLSEL,LOC,X,0
67、DL,ALL,UX,0LSEL,LOC,Y,0DL,ALL,UY,0LSEL,LOC,Y,1.5SFL,ALL,PRES,20ALLSEL,ALL!*求解(qiji)控制*NLGEOM,ONNSUBST,10,100,5SBCTRAN/PBC,U,1/PSF,PRES,NORM,2/NUMBER,1/PNUM,MAT,1EPLOTSAVE,ICONT,db其中ICONT.db数据文件包括此模型的有限元模型、接触对、边界条件、载荷、材料(cilio)定义等。第104页/共123页第一百零四页,共124页。第105页/共123页第一百零五页,共124页。分析(fnx)步骤:Step1.恢复数据库U
68、tilityMenuFileResumefrom选ICONT.db 【OK】Step2.重定位输出到文件UtilityMenuFileSwitchOutputtoFileICONT.output【OK】Step3.求解接触分析/Solusolve产生错误,不能求解UtilityMenuListFilesErrorFileError:两部分之间的小间隙产生刚体模式,导致“自由度超限”错误。DOF(e.gDisplacement)limitexceededattime3.5e-02检查输出文件ICONT.outputANSYS提供接触对信息用于诊断(zhndun)。ICONT.OUTPUT*NOT
69、E*MinInitialgap4.022432169e-3wasdetectedbetweencontactelement264andtargetelement224specifiedbyrealconstantset2.此例中采用最小初始间隙估计初始接触环的值,间隙大小为0.004in,因此采用绝对值ICONT=-0.005in第106页/共123页第一百零六页,共124页。Step4.指定初始(chsh)接触环绝对值-0.005in在单元(dnyun)实常数中赋在向导中赋值Main Menu Preprocessor Real Constant Add/Edit/Delete第107页/共
70、123页第一百零七页,共124页。或用向导(xingdo)赋值MainMenuPreprocessorModelingCreateContactPair选【property】或命令:/Prep7RMODIF,2,5,-0.005!ICONT=-0.005注意(zhy):推荐只用小值,ICONT值太大将导致模型不连续第108页/共123页第一百零八页,共124页。Step5.求解接触分析/SolusolveStep6.重定位(dngwi)输出到输出窗口UtilityMenuFileSwitchOutputtoOutputWindow或命令:/OUTPUT,TERMStep7.后处理:绘位移总矢量
71、MainMenuGeneralPostprocPlotResultsContourPlotNodalSoluDOFSolutionUSUM【OK】或命令:/POST1PLNSOL,U,SUMStep8.检查输出文件ICONT.output*NOTE*AbsoluteinitialclosureICONT0.5000e-02指定的初始接触环调整带成功的防止了刚体运动。第109页/共123页第一百零九页,共124页。第110页/共123页第一百一十页,共124页。本例也可以用向导AutomaticContactadjustment:closegap或单元(dnyun)设置:AutomaticCNO
72、FadjustmentKSclosegap防止刚体运动而不输入ICONT=-0.005或命令:/prep7KEYOPT,3,5,1第111页/共123页第一百一十一页,共124页。例5.两刚性卡间夹一圆球分析(刚体模式-应用初始穿透(chuntu)范围PMIN-PMAX)目标:验证采用初始穿透范围来克服(kf)接触分析中的刚体模式采用合适的初始穿透范围绝对值(-PMIN&PMAX)求解接触分析。模型描述:轴对称、1/2对称施加压力(图7-8)文件:PMIN-PMAX.INP/prep7k,1,0,2.5$k,2,2,2.5k,3,2,0.5k,4,1.25,0.5$k,5,0,1.5l,1,2
73、$l,2,3$l,3,4$l,4,5l,5,1图7-8第112页/共123页第一百一十二页,共124页。al,1,2,3,4,5!A1k,6,0,1$k,7,1,0$k,8,0,0larc,6,7,8,1.0l,7,8$l,8,6al,6,7,8ET,1,182,1!Keyopt(3)=1AxismetricET,2,169$ET,3,171MP,EX,1,30E6$MP,EX,2,1000MP,MU,3,0ESIZE,0.1AMESH,1!type,1;Mat,1MAT,2AMESH,2!*创建(chungjin)目标单元*TYPE,2REAL,2LSEL,LOC,Y,1! 选L4NSLL,
74、1ESURF!*创建(chungjin)接触面单元*TYPE,3!Mat,2;Real,2CSYS,1LSEL,LOC,X,1第113页/共123页第一百一十三页,共124页。NSLL,1ESURFCSYS,0ALLSEL,ALL!*加约束和载荷*/SOLUTIONLSEL,LOC,X,0DL,ALL,UX,0LSEL,LOC,Y,0DL,ALL,UY,0LSEL,LOC,Y,2.5SFL,ALL,PRES,50ALLSEL,ALL!*求解(qiji)控制*NLGEOM,ONNSUBST,10,100,5SBCTRAN/PBC,U,1/PSF,PRES,NORM,2/NUMBER,1/PNUM
75、,MAT,1EPLOTSAVE,PMIN-PMAX.db第114页/共123页第一百一十四页,共124页。分析(fnx)步骤:Step1.恢复数据库PMIN-PMAX.dbUtilityMenuFileResumefrom选PMIN-PMAX.db.db【OK】Step2.重定位输出到文件UtilityMenuFileSwitchOutputtoFilePMIN-PMAX.output【OK】或命令:/output,PMIN-PMAX,outputStep3.求解接触(jich)分析/Solusolve出错ErrorDofe.gDisplacementlimitexceededettime0.
76、1【OK】第115页/共123页第一百一十五页,共124页。Step4.检查输出文件PMIN-PMAX.output*NOTE*Nocontactwasdetectedforacontactpairspecifiedbyrealconstantset2具有“大”的初始(chsh)间隙的载荷控制问题可能引起检测不到接触和或自由度超限错误(图7-9)。包含小的初始(chsh)穿透范围将目标面具有初始(chsh)接触状态,以阻止刚体运动。Step5.指定初始(chsh)穿透范围绝对值在实常数中指定在接触向导中指定MainMenuPreprocessorRealConstantAdd/Edit/Del
77、ete选Set2【Edit】大间隙,产生刚体模式,引起自由度超限错误选Type,3conta171【OK】Ini.allow.penetrationPMIN=-0.001【OK】【close】或命令:/prep7RMODIF,2,7,-0.002,-0.001第116页/共123页第一百一十六页,共124页。或向导(xingdo)MainMenuPreprocessorModelingCreateContactPairproperty第117页/共123页第一百一十七页,共124页。Step6.求解接触分析/SolusolveStep7.重定位(dngwi)输出到输出窗口UtilityMenu
78、FileSwitchOutputtoOutputWindow或命令:/OUTPUT,TERM第118页/共123页第一百一十八页,共124页。Step8.绘轴线(zhuxin)位移(Y向)MainMenuGeneralPostprocPlotResultsContourPlotNodalSoluDOFSolutionUY【OK】或命令:/POST1PLNSOL,U,Y(图7-10)图7-10注意:此结果中等值线还包括由于目标面调整产生的刚体(gngt)位移,因此接触面未接触第119页/共123页第一百一十九页,共124页。Step9.检查输出(shch)文件PMIN-PMAX.output*A
79、bsoluteMax.allowpenetrationrangePMAX0.2000E-2AbsoluteMin.allowpenetrationrangePMIN0.1000E-2*NOTE*TheMax.Initialpenetration1.002326104e-3hasbeenadjustedtomatchallowablerange(1.e-03,2.e-03),forthetargetsurfacespecifiedbyrealconstantset2.*NOTE*Forthispurposetheentirerigidtargetsurfacespecifiedbyrealcon
80、stantset2hasbeenmovedbyX=0,Y=-0.221729606,Z=0指定的初始穿透范围成功阻止了刚体运动,由初始调整引起目标面移动列于结果文件中。Forceconvergencevalue-0.8867criterion=1.471SolutionConvergedAfterequilibriumiteration2*NOTE*Thedisplacementvectorinresultsfile(.rst)hasbeenupdatedduetorigidbodymotionoftargetsurfaceindentifiedbyrealconstantset2:X=0,Y
81、=-0.221729606,Z=0.Solutionisdone!第120页/共123页第一百二十页,共124页。ANSYS在结果文件中(PMIN-PMAX.RST)自动更新位移矢量以包含目标面刚体运动。不包含刚体位移的结果等值线图中最小Y向位移-0.189845与目标面调整值Y=-0.22173组合,总Y向位移为-0.411575:(-0.189845)+(-0.22173)=-0.41149Step10.读入最后子步的结果MainMenuGeneralPostprocReadResultsLastSet命令:Set,last结果文件包含目标面刚体运动的位移,因此,为对包含刚体位移的结果进行
82、后处理,必须用SET命令从结果文件(.rst)中读入数据Step11.绘制(huzh)轴线位移图7-11MainMenuGeneralPostprocPlotResultsContourPlotNodalSoluDOFSolutionUY【OK】或命令:PLNSOL,U,Y(图7-11)图7-11第121页/共123页第一百二十一页,共124页。第122页/共123页第一百二十二页,共124页。感谢您的观赏(gunshng)!第123页/共123页第一百二十三页,共124页。内容(nirng)总结1 接触分类。只在一个表面特别刚硬并且不关心刚硬物体的应力时有效。除了在表面间传递法向压力外,接触单元还传递切向运动(摩擦)。3最小值:该选项为下一子步、预报时间增量(zn lin)的最小值(计算费用十分昂贵,建议不用)。给出输出控制(要求输出每一子步结果)。MAT未重新定义,就是前面的MAT, 1。L16分网(目标2分网)。自动二分:如果接触状态变化明显,时间步长将二分。Step 11.绘制轴线位移图7-11第一百二十四页,共124页。