ANSYS高级接触分析

ANSYS高级接触问题,接触问题概述在工程中会遇到大量的接触问题，如齿轮的啮合、法兰联接、机电轴承接触、卡头与卡座、密封、板成形、冲击等等。接触是典型的状态非线性问题，它是一种高度非线性行为。接触例子如下图：,分析中常常需要确定两个或多个相互接触物体的位移、接触区域的大小和接触面上的应力分布。接触分析存在两大难点：在求解之前，你不知道接触区域的范围；表面之间是接触还是分开是未知的；表面之间突然接触或突然不接触会导致系统刚度的突然变化。大多数接触问题需要计算摩擦。摩擦是与路径有关的现象，摩擦响应还可能是杂乱的，使问题求解难以收敛。,ANSYS高级接触问题,1接触分类,刚柔一个表面是完全刚性的除刚体运动外无应变、应力和变形，另一表面为软材料构成是可变形的。只在一个表面特别刚硬并且不关心刚硬物体的应力时有效。柔柔两个接触体都可以变形。,2接触单元,ANSYS采用接触单元来模拟接触问题：跟踪接触位置；保证接触协调性（防止接触表面相互穿透）；在接触表面之间传递接触应力（正压力和摩擦）。接触单元就是覆盖在分析模型接触面上的一层单元。在ANSYS中可以采用三种不同的单元来模拟接触：面一面接触单元；点一面接触单元；点一点接触单元。,不同的单元类型具有完全不同的单元特性和分析过程。1.面一面接触单元用于任意形状的两个表面接触不必事先知道接触的准确位置；两个面可以具有不同的网格；支持大的相对滑动；支持大应变和大转动。例如:面一面接触可以模拟金属成型，如轧制过程。,2接触单元,2接触单元,2.点一面接触单元用于某一点和任意形状的面的接触可使用多个点面接触单元模拟棱边和面的接触；不必事先知道接触的准确位置；两个面可以具有不同的网格；支持大的相对滑动；支持大应变和大转动。例：点面接触可以模拟棱边和面之间的接触,2接触单元,2接触单元,3.点点接触单元用于模拟单点和另一个确定点之间的接触。建立模型时必须事先知道确切的接触位置；多个点点接触单元可以模拟两个具有多个单元表面间的接触；每个表面的网格必须是相同的；相对滑动必须很小；只对小的转动响应有效。例如:点一点接触可以模拟一些面的接触。如地基和土壤的接触,2接触单元,2接触单元,3关于耦合和约束方程的应用,如果接触模型没有摩擦，接触区域始终粘在一起，并且分析是小挠度、小转动问题，那么可以用耦合或约束方程代替接触。使用耦合或约束方程的优点是分析还是线性的,接触问题的一般特性,1接触刚度1、所有的ANSYS接触单元都采用罚刚度（接触刚度）来保证接触界面的协调性,在数学上为保持平衡，需要有穿透值然而，物理接触实体是没有穿透的分析者将面对困难的选择：小的穿透计算精度高，因此接触刚度应该大；然而，太大的接触刚度会产生收敛困难：模型可能会振荡，接触表面互相跳开。接触刚度是同时影响计算精度和收敛的最重要的参数。你必须选定一个合适的接触刚度。除了在表面间传递法向压力外，接触单元还传递切向运动（摩擦）。采用切向罚刚度保证切向的协调性。（图12）作为初值，可采用：Ktangent=0.01Knormal切向罚刚度与法向罚刚度以同样的方式对收敛性和计算精度产生影响。,2、接触刚度的选取选定一个合适的接触刚度值需要一些经验。对于面一面接触单元，接触刚度通常指定为基体单元刚度的一个比例因子。开始估计时，选用FKN=1.0大面积实体接触FKN=0.01-0.1较柔软（弯曲占主导的部分）另外，也可以指定一个绝对刚度值，单位：（力/长度）/面积。对于点一点（除CONTA178）和点面接触单元需要为罚刚度KN输入绝对值：初始估计时：对于大变形:0.1*EKN1.0*E对于弯曲:0.01*EKNModelingCreateElementsSurf/ContactSurftoSurf(ESURF)对于直接生成刚性目标面，在建立目标单元之前需要要指定附加的单元属性TSHAP,刚性目标面的自动划分不需要TSHAP。ANSYS能根据实体模型确定合适的目标单元形状。划分线(LMESH)2D刚性目标面划分面(AMESH)3D刚性目标面创建关键点(KMESH)控制节点（Pilot）刚性目标面能与控制点联系起来。Pilot实际上是只有一个节点的单元，通过这个节点的运动可以控制整个目标面的运动。ANSYS只在Pilot节点上检查边界条件而忽略其它节点的约束。对可变形体目标面建立目标单元的步骤是：1.先选择可变形体表面上的节点2.然后在可变形体上建立目标单元MainMenuPreprocessorModelingCreateElementsSurf/ContactSurftoSurf,ANSYS将根据基体的网格确定目标单元形状和外法线方向。检查外法线方向（这在自动划分刚性目标面时非常重要）图3-3打开单元坐标系标志并重绘单元/PSYMS,ESYS,1目标单元外法线方向应该指向接触面。如果单元法向不指向接触面，用命令使之反转：ESURF,REVE,例：Seal.dat(图3-3)Step5.建立接触面单元设置接触单元属性、选择可变形体表面节点，并在可变形体上建立接触单元（过程与在可变形体上建立目标单元相同）MainMenuPreprocessorModelingCreateElementsSurf/ContactSurftoSurf(ESURF)这些接触单元与基体有同样的阶数（低阶或高阶）。注意，在壳或梁单元上建立目标单元或接触单元时，可以选择要在梁或壳单元的顶层还是底层建立单元。,图3-3,在选择柔体表面上的节点时，如果你确定某一部分节点永远不会接触到目标面，可以忽略它，以减少计算时间。接触面的外法向应指向目标面。如果发现外法线方向不正确，用下列命令修改之ESURF,REVEStep6.在有限元模型上施加边界条件如果目标面是刚性面，目标面将会自动固定。定义了Pilot点ANSYS只检查该点的边界条件，忽略目标面上其它节点约束。控制点能控制目标面的运动。对Seal.dat施加的边界条见图33。,Step7.定义求解选项和载荷步，以下是默认设置推荐使用N.L求解自动控制使用不带自适应下降的fullNewton-Raphson法求解时间步必须足够小。使用自动时间步。子步数的最大值（NSBMX）应较大，最小值（NSBMIN）应较小Step8.求解Step9.后处理结果包括位移、应力、应变和接触等信息。接触压力、摩擦应力、总应力、接触侵入、接触间隙距离、滑动距离和接触状态都可以从/POST1或/POST26中得到。,面一面接触分析实例（建立接触对不通过接触向导）实例1：弹性环装配,轴对称施加位移载荷：在L45线上施加0.4的Y向位移打开几何N.L开关(NLGEOM,on)设Time=0.4并为自动时间步给出子步数(20,500,10)给出输出控制（要求输出每一子步结果）求解并查看输出和监视文件重启动分析亦将Y向偏移量增加到0.55使咬接装配的第2个齿咬合。问题以弯曲为主，设FKN=0.1,GUI方式：Step1.恢复数据库文件Snap.db（Snap.db数据库包含此有限元模型的完整几何模型、材料、边界条件。但未定义目标单元与接触单元）。UtilityMenuFileResumefrom选Snap.db【OK】基体单元:Plane42,1轴对称材料:MAT,1EX=0.175e6;NUXY=0.35;MU=0.0MAT,2EX=90000;NUXY=0.35约束节点163UY=0182UY=0,Step2.创建目标面单元类型MainMenuPreprocessorElementTyperAdd/Edit/Delete【Add】contact2Dtarge169Elementtyperreferencenumber=2【OK】或命令:/PREP7ET,2,Targe169Step3.建立接触面单元类型MainMenuPreprocessorElementTyperAdd/Edit/Delete【Add】contact2ndSurf171Elementtyperreferencenumber=3【OK】【close】或命令:ET,3,conta171,Step4.指定接触法向刚度MainMenuPreprocessorRealConstantsAdd/Edit/Delete【Add】Type3CONTA171【OK】RealConstantSetNo.=1NormalPenaltyStiffnessFKN=0.1(对弯曲问题采用初值0.1)PenetrationtoleranceFTOLN=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值太小。如果需要多次迭代才能使残值收敛而不是侵入。FKN值可能估计得太高。,Step5.创建目标单元(1)为目标面选择线UtilityMenuSelectEntities拾取线(图3-5)【OK】或命令：LSEL,S,2,4LSEL,A,15,18图3-6LSEL,A,63(2).选择附于线上的全部节点(图3-6)UtilityMenuSelectEntitiesNodesAttachedtolines,all【OK】或命令：NSLL,S,1NPlot,图3-5,图3-6,(3).设置单元属性MainMenuPreprocessorCreateElementElemAttributes或命令：Type,2Mat,1Real,1,(4).创建目标单元MainMenuPreprocessorModelingCreateElementsSurf/ContactSurftoSurf注意：如果基体单元是2D平面或3D实体单元“Tlab”无效，如果基体单元是2D单元，“Shape”无效。选【pickall】或命令：ESURF(图3-7),图3-7,Step6.创建接触单元(1)为接触面选线UtilityMenuSelectEntitiesLines图3-8ByNum/pickFromFull【OK】选线(图3-8)【OK】或命令:LSEL,s,33,34LSEL,a,43,44(2)选择附于选定线上的全部节点图3-9UtilityMenuSelectEntitiesNodesAttachedtoLines,allFromFull【OK】或命令:NSLL,S,1Nplot(图3-9),图3-8,图3-9,(3)设置单元属性MainMenuPreprocessorCreateElementElemAttributesElementtypenumber=3ConTA171Materialnumber=1Realconstantsetnumber=1【OK】或命令：Type,3Mat,1Real,1(4)创建接触面单元图3-11MainMenuPreprocessorModelingCreateElementsSurf/ContactSurftoSurf图3-10【OK】【PICKALL】(图3-10)或命令:ESURF,Step7.选择所有选项并画单元UtilityMenuSelectEverythingUtilityMenuPlotElements或命令:Allsel,allEplot(图3-11)Step8.求解/soluantype,staticnlgeom,onsolc,ontime,0.4nsubst,20,500,10outres,all,alldl,45,uy,-0.4solvedl,45,uy,-0.55solvefiniStep9.后处理,4应用接触向导创建接触对,所有面面接触的单元选项和参数都可以通过接触向导来控制。使用接触向导进行接触分析方便快捷：自动定义单元类型和实常数快速得到接触选项和参数快速显示和反转接触法向使用接触向导必须首先对基体进行网格划分，否则不能激活接触向导。下面仍以弹性环装配接触问题为例。采用接触向导完成接触对的创建来说明接触向导的使用方法。Step1.恢复数据文件UtilityMenuFileResumefrom选Snap.db【OK】(图4-1)此数据库Snap.db包括此有限元模型的几何、基体单元、分网、材料和边界条件。但不包括接触单元与目标单元及接触有关的参数。,Step2.启动接触向导MainMenuPreprocessorModelingCreateContactPair,Step3.创建接触对,图4-3图4-4【PickTarget】(图43),拾取线(见图35)【OK】，再定义接触面返回图44【Next】(图45),图45,【openingsetting】图46,图46,Step4.设置接触参数【OK】返回图45【create】,【Finish】返回图42,退出Step5.求解求解过程同前,图42,5实例2平面密封圈计算,验证为接触刚度估计一个合适的初始值的重要性验证接触分析中摩擦的应用弹性库仑摩擦验证基于接触单元临界状态变化的时间步预测如何会对计算效率有害模型描述：2D超弹平面密封圈对称，密封圈压缩模拟,分析此模型采用两种方法：1.不使用接触向导，建立命流文件Seal.inp2.使用接触向导创建接触对，并用GUI方式创建Pilot节点。本例具有：几何N.L（大应变与大变形）材料N.L（超弹）接触N.L数据文件:Seal.inp,Step1.建模：单元类型、实常数、材料特性、基体分网,/prep7et,1,56!HYPER562D4nodeU_PHyperelasticSolidmp,nuxy,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命令重复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.0k,28,1.1,1.0lstr,26,25!L9lstr,27,28!L10,图5-2,!*基体分网*lesize,8,2$lesize,13,4$lesize,14,4esize,0.035type,1mat,1amesh,allsave,seal,dbStep2.创建接触对，先定义单元和接触参数et,2,169!Target1692DTargetSegmentet,3,171!Conta1712DSurfacetoSurfaceContactmp,mu,1,0!用材料特性定义摩擦，本例先无摩擦计算然后再有摩擦计算，比较结果r,1,0.1,0.1!用实常数定义接触高级选项对于弯曲为主的问题，采用接触刚度FKN=0.1作为初始估计值,FTOLN=0.1lagrange穿透比例系数（缺省值）r,2,0.1,0.1!如果FKN=1.0则不收敛,!*创建接触对1（定义实常数和分网）*,type,2!目标面1Target169real,1!Target169的实常数mat,1!接触面1定义材料与超弹单元同lesize,15,1!目标面1（L15）分网尺寸(图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!选择和所选线相连系的节点type,3!Conta171这里未发布Real,1，实常数未重新定义就是前面的Real,1，MAT也是前面的MAT,1esurf!生成接触单元lsel,all!*检查外法线方向【OK】*,图53,!*创建接触对2*type,2!目标面2Target169real,2!MAT未重新定义，就是前面的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】*,!*创建接触对2*,!刚性面1被约束!刚性面2随Pilot点移动!L4对称约束dl,4,ux,0n_load=node(kx(27),ky(27),0)!为定义刚性面2的位移做准备finish,Step3.施加边界条件,/solunlgeom,onsolc,on!N.L求解自动控制打开（缺省）time,0.85!载荷步、结束时间d,n_load,uy,-0.85nsubst,25,500,10outres,all,allmonitor,var3,n_load,fy,Step4.定义求解选项和载荷步,Step5.求解SolveStep6.查看结果/post1pldsp,2!变形图plnsol,s,eqv,0,1!VonMises应力云图plnsol,cont,pres,0,!接触压力图plnsol,epto,eqv!绘等效总应变图54save,seal,db,定义反力变量、绘载荷变形图,UtilityMenuPlotElementsMainMenuTimeHistPostproDefineVariables【Add】Reactionforce【OK】拾取节点263(控制点)【OK】User-SpecifiedLoad=FORCE:FY【OK】【close】或命令：/POST26RFORCE,2,263,F,Y,FORCEMainMenuTimeHistPostproGraphVariables1stVariabletogragh=2【OK】或命令：PLVAR,2(图54),图54,2、用接触向导创建接触对，用GUI方式创建Pilot节点（刚性面控制点），然后求解。Step1.恢复数据库文件Seal.db（包括基体的几何、单元、分网；没有选接触单元与目标单元，未定义接触对）UtilityMenuFileResumefrom选Seal.db【OK】Step2.启动接触向导MainMenuPreprocessorModelingCreateContactPair图42【contactWizard】,刚性目标面-使用Pilot节点,创建接触对1,图55,选线(1)(接触对1的目标面)【OK】返回图55【next】,图56选线(2)(接触对1的接触面)【OK】；返回图56【next】,Step3.设置接触参数,图57,图58图59,定义摩擦(图59)【OK】返回图57Coefficientfriction0.2【Create】,511图512接触对1(图511),图510图,Step5.创建目标单元控制点,图514,(1).设置单元属性用命令流：Type,4MAT,1REAL,4TSHAP,PIL10GUI:MainMenuPreprocessorModelingCreateElementElemAttributes,(2).绘关键点UtilityMenuPlotKeypointsKeypoints(3).建立目标单元控制点MainMenuPreprocessorMeshingMeshKeypoints选K27【OK】或命令流：KMESH,27Step6.施加位移约束控制点27上施加UY=-0.85MainMenuSolutionDefineLoadsApplyStructuralDisplacementOnKeypoints选K27【OK】Lab2DOFstobeconstrained:UYValue:-0.5【OK】,Step5.创建目标单元控制点,施加对称边界条件：MainMenuSolutionDefineLoadsApplyStructuralDisplacementSymmetryB.C.OnLines选线L4【OK】Step7.求解控制、求解MainMenuSolutionAnalysisTypeSolnControls,Step5.创建目标单元控制点,MainMenuSolutionSolveLS_CurrentStep8.后处理同1无摩擦(MU=0)VonMises=145.096有摩擦(MU=0.2)VonMises=142.038接触刚度取1.0不收敛0.1收敛/POST26MainMenuTimeHistPostpro选控制点27，绘图(见图5-4)。,Step5.创建目标单元控制点,面一面接触具有20个可用实常数，2个材料属性和30个可用单元选项设置。能够模拟特殊的效果和处理困难的收敛情况。通常分析开始先用缺省值，只指定罚刚度和穿透容差及子步数。只有在缺省设置遇到困难时才用高级选项。所有的单元选项和参数都可以通过接触向导来控制；也可以通过实常数和单元选项来指定。Conta178接触单元、实常数选项见图6-1及6-2。,图6-1,6面一面接触单元高级接触选项,6面一面接触单元高级接触选项,图6-2,进入接触选项菜单用下列过程：,MainMenuPreprocessorModelingCreateContactPair选【property】1、第1个重要选项：如何保证接触协调性：使用修正的拉朗日法（缺省）纯罚函数方法在接触向导中BasicContactalgorithm:AugmentedLagrangeMethodPenaltyMethod对大多数模型，修正的拉格朗日方法能很好地工作。罚函数法推荐用于具有变形很大的单元，很大的摩擦系数和/或用修正的拉格朗日方法时收敛性很差的情况。,1、最重要的选项是法向罚刚度或接触刚度,对于大面积接触采用起始值（因子）1.0，对于柔性接触采用0.1。大值对应较高精度；小值对应较好的收敛性。有时最好先以较小的接触刚度进行分析，然后在一系列载荷步中逐渐增大刚度“渐变”的接触刚度，提高收敛性图6-2在最后的载荷步逐渐提高到一个刚硬的值将提高计算精度。可通过向导的基本表或单元特性菜单设置以允许程序更新接触刚度接触向导BasicNormalPenaltyStiffness1.0factorconstant向导Basiccontactstiffnessupdate:Eachloadstep(PAIRIDbased)noneEachloadstep(允许用户指定刚度变化)Eachsubstep(允许自动和用户指定变化),2、第2个重要选项穿透容差。,缺省情况下，穿透容差是一个因子乘以基体单元厚度。对于变化很大的网格密度，采用因子会在接触表面的某些部分产生太小的容差，这时采用绝对值可能更好。不要使用太小的容差，因为它总是对收敛性有害。,3、对于临界接触状态变化的自动时间步控制,或在单元选项中控制：Contacttime/loadpredictionK7:Noprediction,1不控制：不影响自动时间步，对静力问题自动时间步打开时此选项一般是足够的。2自动二分：如果接触状态变化明显，时间步长将二分。对于动力问题自动二分通常是足够的。3合理值：比自动细分更耗时的算法。4最小值：此选项为下一步预测最小时间增量（很耗计算时间，不推荐）,4、Pinball区域影响接触状态的确定和其它许多接触特性,Pinball区域是环绕接触单元的园（2D）或球（3D），描述接触单元周围“远”和“近”区域的边界(图6-3)。在缺省情况下，Pinball区域半径是4基体单元厚度（刚柔）或2基体单元厚度（柔一柔）可以为Pinball半径指定一个不同的值。,图6-3,也可用实常数PINB调整球形区（对于初始值侵入大的问题是必要的）。,5、几种不同的接触模式这些选项使你能够模拟特殊的物理现象。,或单元选项:BehaviorcontactsurfaceK12:standard这些选项包括：标准：正常的接触闭合和打开行为，具有正常的粘着/滑动摩擦行为。粗糙：正常接触闭合和打开行为，但不发生滑动（类似于具有无限摩擦系数）不分离（滑动）：目标面和接触面一旦建立接触就不再分离（允许滑动）绑定：目标面和接触面一旦接触就粘在一起不分离（永远）：初始位于Pinball区域内或已经接触的接触检查点在法向不分离（允许滑动）绑定接触（永远）：初始位于Pinball区域内或已经接触的接触检查点在剩余的分析过程中绑定在一起（DesignSpace缺省值）绑定接触（初始接触）：只在初始接触的地方采用绑定，初始分开的地方保持分开。计算实例：悬臂梁端部旋转（不分离行为）,6、影响某些表面行为的选项,ContactOpeningstiffness（分开时的间隙刚度）保证不分离和绑定行为，它通过使用当存在间隙也具有非零刚度的弹簧来连接表面。缺省情况下，此弹簧刚度等于法向罚刚度，其效果类似于法向罚刚度刚度太小精度低；刚度太大会引起收敛问题。,摩擦系数影响基本摩擦行为：,Contactcohesion表示当没有法向压力时开始滑动的摩擦应力值。摩擦导致非对称刚度阵。因为非对称矩阵很难计算（因此导致求解变慢），程序自动控制执行对称求解，利用此算法可以解决多数含摩擦接触问题。有时，采用非对称矩阵能获用更好的收敛性。如果遇到收敛缓慢问题可以用不对称求解选项。记住：这种情况必须使用稀疏或波前求解器。对于每个支持非对称矩阵的单元，此选项也可以由下列菜单激活：MainMenuSolutionUnabridgedMenuAnalysisOptions设置Newton-Raphson选项为FullN_Runsymm,7、初始穿透,有几种技术可以模拟初始穿透接触问题（如过盈装配）。可以使用初始几何穿透，或指定偏移量，或二者皆有。(图6-5)指定偏移量（CNOF）,或在实常数中指定偏移量（CNOF）contactsurfaceoffsetCNOF:0.025,正的CNOF加大初始穿透负的CNOF减小初始穿透或导致间隙CNOF可与几何穿透组合自动CNOF调整允许ANSYS基于初始穿透自动给定CNOF值。导致“刚好接触”配置ICONT缺省为0,或单元选项:AutoCNOF/ICONTadjustmentK5:No.Auto.Adjust,7、初始穿透,初始穿透选项包括：Includeeverything：包括由几何模型和指定偏移量（如果有的话）引起的初始穿透Excludeeverything：忽略所有初始穿透效应。Includewithrampedeffects：渐变初始穿透，以提高收敛性。Includeoffsetonly：只包括由偏移量指定的基本初始穿透。Includeoffsetonlyw/ramp：只包括由偏移量指定的基本初始穿透，且渐变初始穿透以提高收敛性,7、初始穿透,如果模型包含初始几何穿透，接触力将立即“阶跃”到一个大值。载荷突变经常导致收敛困难，期望有一种机制能够将初始穿透效应渐变到零。Includewithrampedeffects和Includeoffsetonlyw/ramp选项通过在第一载荷步，将初始穿透渐变为零克服收敛困难。为求得好的结果，在第一载荷步不应施加其它载荷(图6-6)。计算实例：初始穿透,图6-6,7、初始穿透,初始不接触的两个（或多个）物体的静力分析中，在接触建立前可能产生刚体运动(图6-7)。此例中圆柱体没有施加位移约束，面由力控制。圆柱体的约束由圆柱体和平板之间的接触建立。求解过程中两个物体分离，刚度矩阵奇异。ANSYS将产生一个负主元警告。有几个选项可以解决由于初始不相连物体引起的刚体模式：,图6-7,8、刚体模式,三个高级接触特性允许调整初始接触条件以防止刚体模式：(1).自动CONF调整程序计算CNOF以清除间隙。(2).初始接触环（ICONT）将调整带内接触表面上的节点移到目标面上(3).初始允许穿透范围(PMINMat,1TYPE,2MAT,2REAL,2LSEL,LOC,X,1.9NSLL,1!1.9线上的节点ESURF!创建目标单元Target169!*TYPE,3!Mat,2;Real,2LSEL,LOC,X,2NSLL,S,1ESURF!创建接触面conta172ALLSEL,ALL,/SOLUTIOND,NODE(1.5,0,0),UY,0D,NODE(2.4,0,0),UY,0NLGEOM,ONNSUBST,10,50,5FINISH/PBC,U,1/NUMBER,1/PNUM,MAT,1EPLOTSAVE,interference,db!包含有限元模型、接触对、材料、边界条件，求解控制,例3.紧配合环（初始穿透）,分析步骤：Step1.恢复数据库文件UtilityMenuFileResumefrom选interference.db【OK】Step2.重定位输出文件UtilityMenuFileSwitchOutputtoFileInterference.output【OK】或命令：/Output,interference,OutputStep3.求解接触分析/SolusolveStep4.重定位输出到输出窗口UtilityMenuFileSwitchOutputtoOutputWindow或命令：/OUTPUT,TERMStep5.绘制径向应力MainMenuGeneralPostprocPlotResultsContourPlotNodalSoluStressX-directionSX【OK】,例3.紧配合环（初始穿透）,或命令:/POST1PLNSOL,S,X!结果：径向应力0Step6.列表显示接触结果MainMenuGeneralPostprocListResultsNodalSolution选NonlinearitemscontactCONT【OK】或命令：PRNSOL,CONTSTAT=0“张开”几乎不接触STAT=1“张开”但几乎接触STAT=2“闭合”并滑动STAT=3“闭合”并粘着,例3.紧配合环（初始穿透）,Step7.检查输出文件内容（interference.output）*NOTE*Nocontactwasdetectedforacontactpairspecifiedbyrealconstantset2大初始穿透问题可能导致无法检测到接触或检测到接触面求解存在收敛困难。在此情况下推荐在第一个载荷步中包含渐进化效应。Step8.为初始穿透选项指定渐进化效应可通过conta172单元选项设置可通过接触向导中InitialAdjustment的设置(1)MainMenuPreprocessorElementTypeAdd/Edit/Delete选Type3conta172【Option】initialpenetration/gapK9=include_ramp【OK】【close】或命令:/prep7KEYOPT,3,9,2,例3.紧配合环（初始穿透）,(2)或向导InitialAdjustmentInitialpenetration:Includeeverythingwithrampedeffect【OK】,例3.紧配合环（初始穿透）,Step9.重定位输出到文件UtilityMenuFileSwitchOutputtoFileRamped.output【OK】Step10.分析/SoluSolveStep11.重定位输出到输出窗口UtilityMenuFileSwitchOutputtoOutputWindow命令：/output,TERMStep12.绘径向应力MainMenuGeneralPostprocPlotResultsContourPlotNodalSoluStressX-directionSX【OK】或命令:/POST1PLNSOL,S,X,例3.紧配合环（初始穿透）,Step13.将轴对称模型扩展为整个模型显示UtilityMenuPlotCtrlsStyleSymmetryExpansion2DAxi-SymmetricSelectexpansionamount=Fullexpansion【OK】或命令：/EXPAND,36,Axis,10/RepxzlotISO显示/View,1,1,1/ReplotStep14.检查输出文件ramped.outputInitialpenetrationwillberampedduringthefirstloadstep在第一载荷步，初始穿透渐进化。*NOTE*Max.Initialpenetration0.1wasdetectedbetweencontactelement425andtargetelement405specifiedbyreadconstantset2检测到大的初始穿透,例3.紧配合环（初始穿透）,刚体模式应用初始接触环ICONT）目标：验证使用初始接触环调整带来克服接触分析中的刚体模式模型描述：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平面应力、平面应变和轴对称单元、大变形、大应变,图7-6,图7-7,例4.两刚性面间夹圆柱体分析,ET,2,169!Target169ET,3,171!Conta171MP,EX,1,1000MP,EX,2,30E6MP,MU,3,0ESIZE,0.1AMESH,1!type,1;mat,1MAT,2AMESH,2!*创建目标目单元*TYPE,2MAT,3REAL,2LSEL,LOC,Y,1.003NSLL,1ESURF!*创建接触单元*TYPE,3CSYS,1LSEL,LOC,X,1NSLL,1ESURFCSYS,0,LLSEL,ALL/SOLUTIONLSEL,LOC,X,0DL,ALL,UX,0LSEL,LOC,Y,0DL,ALL,UY,0LSEL,LOC,Y,1.5SFL,ALL,PRES,200ALLSEL,ALL!*求解控制*NLGEOM,ONNSUBST,10,100,5SBCTRAN/PBC,U,1/PSF,PRES,NORM,2/NUMBER,1/PNUM,MAT,1EPLOTSAVE,ICONT,db,其中ICONT.db数据文件包括此模型的有限元模型、接触对、边界条件、载荷、材料定义等。,例4.两刚性面间夹圆柱体分析,例4.两刚性面间夹圆柱体分析,分析步骤：,Step1.恢复数据库UtilityMenuFileResumefrom选ICONT.db【OK】Step2.重定位输出到文件UtilityMenuFileSwitchOutputtoFileICONT.output【OK】Step3.求解接触分析/Solusolve产生错误，不能求解UtilityMenuListFilesErrorFileError:两部分之间的小间隙产生刚体模式，导致“自由度超限”错误。DOF(e.gDisplacement)limitexceededattime3.5e-02检查输出文件ICONT.outputANSYS提供接触对信息用于诊断。ICONT.OUTPUT*NOTE*MinInitialgap4.022432169e-3wasdetectedbetweencontactelement264andtargetelement224specifiedbyrealconstantset2.此例中采用最小初始间隙估计初始接触环的值，间隙大小为0.004in，因此采用绝对值ICONT=-0.005in,Step4.指定初始接触环绝对值-0.005in,在单元实常数中赋MainMenuPreprocessorRealConstantAdd/Edit/Delete,或用向导赋值,或命令：/Prep7RMODIF,2,5,-0.005!ICONT=-0.005注意:推荐只用小值,ICONT值太大将导致模型不连续,MainMenuPreprocessorModelingCreateContactPair选【property】,Step5.求解接触分析/SolusolveStep6.重定位输出到输出窗口UtilityMenuFileSwitchOutputtoOutputWindow或命令：/OUTPUT,TERMStep7.后处理：绘位移总矢量MainMenuGeneralPostprocPlotResultsContourPlotNodalSoluDOFSolutionUSUM【OK】或命令:/POST1PLNSOL,U,SUMStep8.检查输出文件ICONT.output*NOTE*AbsoluteinitialclosureICONT0.5000e-02指定的初始接触环调整带成功的防止了刚体运动。,例4.两刚性面间夹圆柱体分析,例4.两刚性面间夹圆柱体分析,本例也可以用向导AutomaticContactadjustment:closegap或单元设置:AutomaticCNOFadjustmentKSclosegap防止刚体运动而不输入ICONT=-0.005或命令：/prep7KEYOPT,3,5,1,例4.两刚性面间夹圆柱体分析,（刚体模式-应用初始穿透范围PMIN-PMAX）,目标：验证采用初始穿透范围来克服接触分析中的刚体模式采用合适的初始穿透范围绝对值（-PMINMat,1MAT,2AMESH,2!*创建目标单元*TYPE,2REAL,2LSEL,LOC,Y,1!选L4NSLL,1ESURF!*创建接触面单元*TYPE,3!Mat,2;Real,2CSYS,1LSEL,LOC,X,1,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!*求解控制*NLGEOM,ONNSUBST,10,100,5SBCTRAN/PBC,U,1/PSF,PRES,NORM,2/NUMBER,1/PNUM,MAT,1EPLOTSAVE,PMIN-PMAX.db,Step1.恢复数据库PMIN-PMAX.dbUtilityMenuFileResumefrom选PMIN-PMAX.db.db【OK】Step2.重定位输出到文件UtilityMenuFileSwitchOutputtoFilePMIN-PMAX.output【OK】或命令：/output,PMIN-PMAX,outputStep3.求解接触分析/Solusolve出错ErrorDofe.gDisplacementlimitexceededettime0.1【OK】,分析步骤：,Step4.检查输出文件PMIN-PMAX.output*NOTE*Nocontactwasdetectedforacontactpairspecifiedbyrealconstantset2具有“大”的初始间隙的载荷控制问题可能引起检测不到接触和或自由度超限错误(图7-9)。包含小的初始穿透范围将目标面具有初始接触状态，以阻止刚体运动。Step5.指定初始穿透范围绝对值在实常数中指定在接触向导中指定MainMenuPreprocessorRealConstantAdd/Edit/Delete选Set2【Edit】大间隙，产生刚体模式，引起自由度超限错误选Type,3conta171【OK】Ini.allow.penetrationPMIN=-0.001【OK】【close】或命令：/prep7RMODIF,2,7,-0.002,-0.001,或向导MainMenuPreprocessorModelingCreateContactPairproperty,例5.两刚性卡间夹一圆球分析,Step6.求解接触分析/SolusolveStep7.重定位输出到输出窗口UtilityMenuFileSwitchOutputtoOutputWindow或命令：/OUTPUT,TERM,例5.两刚性卡间夹一圆球分析,Step8.绘轴线位移（Y向）MainMenuGener