高级接触问题.pptx

第一章第一章 接触问题概述接触问题概述v在工程中常会遇到大量的接触问题，如齿轮的啮合、法兰联接、机电轴承接触、卡头与卡座、密封、板成形、冲击等等。接触是典型的状态非线性问题，它是一种高度非线性行为。接触例子见图1.1。图1.1 轮齿接触分析中常常需要确定两个或多个相互接触物体的位移、接触区域的大小和接触面上的应力分布。接触分析存在两大难点：1.在求解之前，你不知道接触区域、表面之间是接触或分开是未知的，表面之间突然接触或突然不接触会导致系统刚度的突然变化。2.大多数接触问题需要计算摩擦。摩擦是与路径有关的现象，摩擦响应还可能是杂乱的，使问题求解难以收敛。第一节 接触分类刚柔一个表面是完全刚性的除刚体运动外无应变、应力和变形，另一表面为软材料构成是可变形的。只在一个表面特别刚硬并且不关心刚硬物体的应力时有效。柔柔两个接触体都可以变形。第二节 接触单元ANSYS采用接触单元来模拟接触问题：跟踪接触位置保证接触协调性（防止接触表面相互穿透）在接触表面之间传递接触应力（正压力和摩擦）接触单元就是覆盖在分析模型接触面上的一层单元。在ANSYS中可以采用三种不同的单元来模拟接触：一面一面接触单元一点一面接触单元一点一点接触单元不同的单元类型具有完全不同的单元特性和分析过程。1.面面接触单元用于任意形状的两个表面接触不必事先知道接触的准确位置；两个面可以具有不同的网格；支持大的相对滑动；支持大应变和大转动。例如:面一面接触可以模拟金属成型，如轧制过程，见图1.2。图1.2 面一面接触2.点面接触单元用于某一点和任意形状的面的接触可使用多个点面接触单元模拟棱边和面的接触；不必事先知道接触的准确位置；两个面可以具有不同的网格；支持大的相对滑动；支持大应变和大转动。例：点面接触可以模拟棱边和面之间的接触，见图1.3。图 1.3 点面接触3.点点接触单元用于模拟单点和另一个确定点之间的接触。点点接触应用较少,我们不做介绍.第三节 关于耦合和约束方程的应用如果接触模型没有摩擦，接触区域始终粘在一起，并且分析是小挠度、小转动问题，那么可以用耦合或约束方程代替接触。使用耦合或约束方程的优点是分析还是线性的示例1：利用耦合模拟接触问题，见图。数据库文件：cpnorm.db1条件：接触模型无摩擦，接触区域始终粘在一起，小挠度，小转动问题，可以用耦合代替接触（耦合分析是线性的）。图 模拟接触第二章第二章 接触问题的一般特性接触问题的一般特性在进入接触分析的实施和实例讨论之前，我们先对一般条件下的通用内容进行讨论，以便能正确理解和使用相关选项。第一节 接触刚度1.所有的ANSYS接触单元都采用罚刚度（接触刚度）来保证接触界面的协调性，见图2.1。图2.1 切向罚刚度在数学上为保持平衡，需要有穿透值然而，物理接触实体是没有穿透的分析者将面对困难的选择：小的穿透计算精度高，因此接触刚度应该大；然而，太大的接触刚度会产生收敛困难：模型可能会振荡，接触表面互相跳开。接触刚度是同时影响计算精度和收敛的最重要的参数。你必须选定一个合适的接触刚度。除了在表面间传递法向压力外，接触单元还传递切向运动（摩擦）。采用切向罚刚度保证切向的协调性。见图2.2作为初值，可采用：切向罚刚度以同样的方式对收敛性和计算精度产生影响。图2.2 切向罚刚度2.接触刚度的选取选定一个合适的接触刚度值需要一些经验。对于面一面接触单元，接触刚度通常指定为基体单元刚度的一个比例因子。开始估计时，选用FKN=1.0 大面积实体接触FKN=0.010.1 较柔软（弯曲占主导的）部分另外，也可以指定一个绝对刚度值，单位：（力/长度）/面积。点面接触单元需要为罚刚度KN输入绝对值：初始估计时：对于大变形:0.1*E KN 1.0*E对于弯曲:0.01*E KN Preprocessor Modeling Create Elements Surf/Contact Surf to Surf(ESURF)对于直接生成刚性目标面，在建立目标单元之前需要要指定附加的单元属性TSHAP。TSHAP定义目标面单元形状刚性目标面刚性目标面的自动划分不需要TSHAP。ANSYS能根据实体模型确定合适的目标单元形状。划分线(LMESH)2D刚性目标面划分面(AMESH)3D刚性目标面创建关键点(KMESH)控制节点（Pilot）刚性目标面能与控制点联系起来，Pilot实际上是只有一个节点的单元，通过这个节点的运动可以控制整个目标面的运动。ANSYS只在Pilot节点上检查边界条件而忽略其它节点的约束。对可变形体目标面对可变形体目标面建立目标单元的步骤是：1.先选择可变形体表面上的节点 2.然后在可变形体上建立目标单元GUI:Main Menu Preprocessor Modeling Create Elements Surf/Contact Surf to Surf(ESURF)ANSYS将根据基体的网格确定目标单元形状和外法线方向。检查外法线方向（这在自动划分刚性目标面时非常重要）打开单元坐标系标志并重绘单元/PSYMS,ESYS,1目标单元外法线方向应该指向接触面。如果单元法向不指向接触面，用命令使之反转：ESURF,REVE例：Seal.dat Seal.dat文件中的有限元模型图3-7Step 5 建立接触面单元建立接触面单元设置接触单元属性、选择可变形体表面节点，并在可变形体上建立接触单元（过程与在可变形体上建立目标单元相同）。GUI:Main Menu Preprocessor Modeling Create Elements Surf/Contact Surf to Surf这些接触单元与基体有同样的阶数（低阶或高阶）。注意：在壳或梁单元上建立目标单元或接触单元时，可以选择要在梁或壳单元的顶层还是底层建立单元。在选择柔体表面上的节点时，如果你确定某一部分节点永远不会接触到目标面时，可以忽略它，以减少计算时间。接触面的外法向应指向目标面。如果发现外法线方向不正确，用下列命令修改之ESURF,REVEStep 6 在有限元模型上施加边界条件如果目标面是刚性面，目标面将会自动固定。定义了Pilot点ANSYS只检查该点的边界条件，忽略目标面上其它节点约束。控制点能控制目标面的运动。对Seal.dat施加的边界条件，见图3.7。Step 7 定义求解选项和载荷步，以下是默认设置：推荐使用非线性求解自动控制使用不带自适应下降的full Newton-Raphson法求解时间步必须足够小。使用自动时间步。子步数的最大值（NSBMX）应较大，最小值（NSBMIN）应较小Step 8 求解Step 9 后处理结果包括位移、应力、应变和接触等信息；接触压力、摩擦应力、总应力、接触侵入、接触间隙距离、滑动距离和接触状态都可以从/POST1或/POST26中得到。面一面接触分析实例（建立接触对不通过接触向导）实例1：弹性环装配，见图3.8。图3.8 弹性环装配轴对称施加位移载荷：在L45线上施加0.4的-Y向位移打开几何非线性开关(NLGEOM,on)设Time=0.4并为自动时间步给出子步数(20,500,10)给出输出控制（要求输出每一子步结果）求解并查看输出和监视文件重启动分析并将-Y向偏移量增加到0.55使咬接装配的第2个齿咬合。问题以弯曲为主，设FKN=0.1L45GUI方式：Step 1 恢复数据库文件 Snap.db（Snap.db数据库包含此有限元模型的完整几何模型、材料、边界条件。但未定义目标单元与接触单元）。GUI：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;NUXY=0.35约束节点 163 UY=0 182 UY=0Step 2 创建目标面单元类型GUI：Main Menu Preprocessor Element Typer Add/Edit/Delete【Add】contact 2D target 169Element typer reference number=2【OK】或命令流:/PREP 7ET,2,Target 169Step 3 建立接触面单元类型GUI:Main Menu Preprocessor Element Typer Add/Edit/Delete【Add】contact 2nd Surf 171Element typer reference number=3【OK】【close】或命令:ET,3,conta 171Step 4 指定接触法向刚度GUI:Main Menu Preprocessor Real Constants 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图3.14 建立接触面单元类型确定罚刚度 FKN值通常在0.01-10之间，对于体积变形问题用1.0（默认），对弯曲问题用0.1确定侵入容差：侵入容差(FTOLN)是与接触单元下面的实体单元深度（h）相乘的比例因子。若此值太小会引起收敛困难，绝对不要用太小的容差!增大罚刚度（FKN）将减少侵入。将FKN增大100倍会相应地减少侵入，但是接触压力只改变5。如不收敛可调整FKN或FTOLN值重新运行。检查侵入和每个子步的平衡迭代数，如果收敛受入侵容差的驱使可能是FKN值估计不足或FTOLN值太小。如果需要多次迭代才能使残值收敛而不是侵入。FKN值可能估计得太高。或命令：R,1,0.1,0.1(3)设置单元属性GUI:Main Menu Preprocessor Create Element Elem Attributes或命令流：Type,2Mat,1Real,1Step 5 创建目标单元(1)为目标面选择线GUI:Utility Menu Select Entities拾取线，见图3.9【OK】图3.9 拾取线或命令流文件：LSEL,S,2,4LSEL,A,15,18LSEL,A,63(2)选择附于线上的全部节点，见图3.10。GUI:Utility Menu Select EntitiesNodesAttached to lines,all图3.10 选择附于线上的全部节点【OK】或命令流：NSLL,S,1NPlot图3.9 拾取线图3.10 选择附于线上的全部节点创建目标单元GUI:Main Menu Preprocessor Modeling Create Elements Surf/Contact Surf to Surf注意：如果基体单元是2D平面或3D实体单元“Tlab”无效，如果基体单元是2D单元，“Shape”无效。选【Pick All】或命令：ESURF，见图3.11。图3.11 ESURFStep 6 创建接触单元(1)为接触面选线GUI:Utility Menu Select Entities LinesBy Num/pickFrom Full【OK】选线(见图3.12)【OK】图3.12 为接触面选线或命令流:LSEL,s,33,34LSEL,a,43,44(2)选择附于选定线上的全部节点GUI:Utility Menu Select Entities NodesAttached to Lines,allFrom Full【OK】或命令流:图3.13 选择线上的全部节点NSLL,S,1NPLOT，见图3.13图3.12 为接触面选线图3.13 选择线上的全部节点(3)设置单元属性GUI: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.14 选择所有单元GUI:Main Menu Preprocessor Modeling Create Elements Surf/Contact Surf to Surf【OK】【Pick All】见图3.14或命令:ESURF图3.14 创建接触单元图3.15 选择所有单元Step 7 选择所有选项并画单元GUI:Utility Menu Select EverythingUtility Menu Plot Elements或命令流:Allsel,allEplot 见图3.15Step 8 求解/soluantype,staticnlgeom,onsolc,ontime,0.4nsubst,20,500,10outres,all,alldl,45,uy,-0.4solvedl,45,uy,-0.55solvefiniStep 9 后处理第四节第四节 应用接触向导创建接触对应用接触向导创建接触对所有面面接触的单元选项和参数都可以通过接触向导来控制。使用接触向导进行接触分析方便快捷：自动定义单元类型和实常数快速得到接触选项和参数快速显示和反转接触法向使用接触向导必须首先对基体进行网格划分，否则不能激活接触向导。下面仍以弹性环装配接触问题为例。采用接触向导完成接触对的创建来说明接触向导的使用方法。Step 1 恢复数据库文件GUI:Utility Menu File Resume from选 Snap.db恢复单元【OK】见图3.16此数据库Snap.db包括此有限元模型的几何、基体单元、分网、材料和边界条件。但不包括接触单元与目标单元及接触有关的参数。图图3.16 恢复单元恢复单元Step 2 启动接触向导GUI:Main Menu Preprocessor Modeling Create Contact Pair图3.17 启动接触向导按Contact WizardStep 3 创建接触对图3.18 创建目标面【Pick Target】拾取线(见图3.19)【OK】，再定义接触面返回到图【Next】见图3.20。图3.19/SOLsolveDL,45,-UY,-0.4$第五节 平面密封圈计算实例验证为接触刚度估计一个合适的初始值的重要性验证接触分析中摩擦的应用弹性库仑摩擦验证基于接触单元临界状态变化的时间步预测如何会对计算效率有害模型描述：2D超弹平面密封圈1/2对称，密封圈压缩模拟本例具有：几何非线性（大应变与大变形）材料非线性（超弹）接触非线性数据文件:Seal.inpStep 1 建模：单元类型、实常数、材料特性、基体分网/prep7et,1,56 !HYPER56 2D 4node U_P Hyperelastic Solidmp,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.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$k,26,1.1,0$k,27,-0.6,1.0$k,28,1.1,1.0lstr,26,25 !L9lstr,27,28 !L10!*基体分网*lesize,8,2$lesize,13,4$lesize,14,4esize,0.035type,1mat,1amesh,allsave,seal,db用接触向导创建接触对用接触向导创建接触对，用GUI方式创建Pilot节点（刚性面控制点），然后求解。Step 1 恢复数据库文件Seal.db（包括基体的几何、单元、分网；没有选接接触单元与目标单元，未定义接触对）GUI：Utility Menu File Resume from选 Seal.db【OK】Step 2 启动接触向导GUI：Main Menu Preprocessor Modeling Create Contact Pair【Contact Wizard】创建接触对1。Step 5 创建目标单元控制点(1)设置单元属性 GUI:Main Menu Preprocessor Modeling Create Element Elem Attributes(2)绘关键点GUI:Utility Menu Plot Keypoints Keypoints(3)建立目标单元控制点GUI:Main Menu Preprocessor Meshing Mesh Keypoints选K27 【OK】或命令流：KMESH,27Step 6 施加位移约束控制点27上施加 UY=-0.85GUI:Main Menu Solution Define Loads Apply Structural Displacement On Keypoints选K27 【OK】Lab2 DOFs to be constrained:UYValue:-0.5【OK】施加对称边界条件：GUI:Main Menu Solution Define Loads Apply Structural Displacement Symmetry B.C.On Lines 选线L4 【OK】Step 7 求解控制、求解GUI:Main Menu Solution Analysis Type Soln Controls GUI:Main Menu Solution Solve LS CurrentStep 8 后处理/POST26GUI:Main Menu TimeHist Postpro选控制点269第五节第五节 面面接触单元高级接触选项面面接触单元高级接触选项面一面接触具有20个可用实常数，2个材料属性和30个可用单元选项设置。能够模拟特殊的效果和处理困难的收敛情况。通常分析开始先用缺省值，只指定罚刚度和穿透容差及子步数。只能在缺省设置遇到困难时才用高级选项。所有的单元选项和参数都可以通过接触向导来控制；也可以通过实常数和单元选项来指定。Conta171接触单元和Targe169目标单元的单元选项、实常数见图20。进入接触选项菜单用下列过程：GUI：Main Menu Preprocessor Modeling Create Contact Pair选【Property】我们将按照建立模型时的逻辑顺序讨论这些选项（在接触向导中并无必要按此顺序）。图3.20 单元选项、实常数1.第1个重要选项：如何保证接触协调性。使用修正的拉朗日法（缺省）纯罚函数方法在接触向导中BasicContact algorithm:Augmented Lagrange Method Penalty Method对大多数模型，修正的拉格朗日方法能很好地工作。罚函数法推荐用于具有变形很大的单元，很大的摩擦系数和/或用修正的拉格朗日方法时收敛性很差的情况。最重要的选项是法向罚刚度或接触刚度。对于大面积接触采用起始值（因子）1.0，对于柔性接触采用0.1。大值对应较高精度；小值对应较好的收敛性。有时最好先以较小的接触刚度进行分析，然后在一系列载荷步中逐渐增大刚度“渐变”的接触刚度，提高收敛性在最后的载荷步逐渐提高到一个刚硬的值将提高计算精度。可通过向导的基本表或单元特性菜单设置以允许程序更新接触刚度。向导Basiccontact stiffness update:Each load step(PAIR ID based)none Each load step(允许用户指定刚度变化)Each substep(允许自动和用户指定变化)2.第2个重要选项穿透容差穿透容差缺省情况下，穿透容差是一个因子乘以基体单元厚度。对于变化很大的网格密度，采用因子会在接触表面的某些部分产生太小的容差，这时采用绝对值可能更好。不要使用太小的容差，因为它总是对收敛性有害。3.对于临界接触状态变化的自动时间步控制不控制：不影响自动时间步，对静力问题自动时间步打开时此选项一般是足够的。自动二分：如果接触状态变化明显，时间步长将二分。对于动力问题自动二分通常是足够的。合理值：比自动细分更耗时的算法。最小值：此选项为下一步预测最小时间增量（很耗计算时间，不推荐）或在单元选项中控制：Contact time/load prediction K7:No prediction 对于临界接触状态变化的自动时间步控制4.Pinball区域影响接触状态的确定和其它许多接触特性Pinball区域是环绕接触单元的园（2D）或球（3D），描述接触单元周围“远”和“近”区域的边界，见图3.21在缺省情况下，Pinball区域半径是4基体单元厚度（刚柔）或2基体单元厚度（柔-柔）。描述接触单元周围“远”和“近”区域的边界可以为Pinball半径指定一个不同的值。Pinball赋值也可用实常数PINB调整球形区（对于初始值侵入大的问题是必要的）。5.几种不同的接触模式这些选项使你能够模拟特殊的物理现象。模拟特殊的物理现象或单元选项:Behavior contact surface K12:standard 这些选项包括：标准：正常的接触闭合和打开行为，具有正常的粘着/滑动摩擦行为。粗糙：正常接触闭合和打开行为，但不发生滑动（类似于具有无限摩擦系数）不分离（滑动）：目标面和接触面一旦建立接触就不再分离（允许滑动）绑定：目标面和接触面一旦接触就粘在一起不分离（永远）：初始位于Pinball区域内或已经接触的接触检查点在法向不分离（允许滑动）绑定接触（永远）：初始位于Pinball区域内或已经接触的接触检查点在剩余的分析过程中绑定在一起（Design Space缺省值）绑定接触（初始接触）：只在初始接触的地方采用绑定，初始分开的地方保持分开。计算实例：悬臂梁端部旋转（不分离行为）6.影响某些表面行为的选项Contact Opening stiffness（分开时的间隙刚度）保证不分离和绑定行为，它通过使用当存在间隙也具有非零刚度的弹簧来连接表面。缺省情况下，此弹簧刚度等于法向罚刚度，其效果类似于法向罚刚度刚度太小精度低；刚度太大会引起收敛问题。影响某些表面行为的选项摩擦系数影响基本摩擦行为：MU=0指定无摩擦行为MU 0确定滑动摩擦力最大摩擦应力表示可以在摩擦接触面上传递的最大应力（基体材料的剪切屈服应力）摩擦系数影响基本摩擦行为图3.22 摩擦系数影响基本摩擦行为Contact cohesion表示当没有法向压力时开始滑动的摩擦应力值。摩擦导致非对称刚度阵。因为非对称矩阵很难计算（因此导致求解变慢），程序自动控制执行对称求解，利用此算法可以解决多数含摩擦接触问题。有时，采用非对称矩阵能获用更好的收敛性。如果遇到收敛缓慢问题可以用不对称求解选项。记住：这种情况必须使用稀疏或波前求解器。对于每个支持非对称矩阵的单元，此选项也可以由下列菜单激活：GUI：Main Menu Solution Unabridged Menu Analysis Options设置Newton-Raphson选项为 full N_R unsymm计算实例：不离接触如果使用梁或壳单元建模，为了考虑梁和壳的厚度可以移动接触表面。对于Shell181单元，由于大应变变形引起的厚度变化进行也在考虑之中，见图3.23。图3.20 Shell181单元7.初始穿透有几种技术可以模拟初始穿透接触问题（如过盈装置匹配）。可以使用初始几何穿透，或指定偏移量，或二者皆有。指定偏移量（CNOF）指定偏移量（CNOF）或在实常数中指定偏移量（CNOF）contact surface offset CNOF:0.025 正的CNOF，加大初始穿透 负的CNOF，减小初始穿透或导致间隙CNOF可与几何穿透组合，见图 自动CNOF调整允许ANSYS基于初始穿透自动给定CNOF值。导致“刚好接触”配置ICONT缺省为0 CONF与几何穿透组合 自动CNOF调整或单元选项:Auto CNOF/ICONT adjustment K5:No.Auto.Adjust初始穿透选项包括：初始穿透选项包括：Include everything：包括由几何模型和指定偏移量（如果有的话）引起的初始穿透Exclude everything：忽略所有初始穿透效应。Include with ramped effects：渐变初始穿透，以提高收敛性。Include offset only：只包括由偏移量指定的基本初始穿透。Include offset only w/ramp：只包括由偏移量指定的基本初始穿透，且渐变初始穿透以提高收敛性 初始穿透选项如果模型包含初始几何穿透，接触力将立即“阶跃”到一个大值。载荷突变经常导致收敛困难，期望有一种机制能够将初始穿透效应渐变到零。Include with ramped effects和Include off set only w/ramp选项通过在第一载荷步，将初始穿透渐变为零克服收敛困难。为求得好的结果，在第一载荷步不应施加其它载荷，见图 图 第一载荷步计算实例：初始穿透图 初始不接触两个物体图图3.633.638.刚体模式刚体模式初始不接触的两个（或多个）物体的静力分析中，在接触建立前可能产生刚体运动，见图。此例中圆柱体没有施加位移约束，面由力控制。圆柱体的约束由圆柱体和平板之间的接触建立。求解过程中两个物体分离，刚度矩阵奇异。(1)刚好接触，见图 图 刚好接触这需要了解何处是刚好接触位置。如果表面弯曲或不规则这将很困难。由于有限元网格划分时的数值误差，在物体之间可能会存在小间隙或小穿透。这可能引起不收敛或接触物体弹开。三个高级接触特性允许调整初始接触条件以防止刚体模式。三个高级接触特性允许调整初始接触条件以防止刚体模式。(1)自动CONF调整程序计算CNOF以清除间隙。(2)初始接触环（ICONT）将调整带内接触表面上的节点移到目标面上(3)初始允许穿透范围(PMIN&PMAX）：将刚体表面移动到接触面上调整初始接触条件（ICONT），见下图。实常数ICONT可用于指定目标面上的“调整环”。位于调整环内的任何接触点都要移动到目标面上。推荐只进行很小的修正，如果ICONT值太大会产生不连续。如果未指定常数ICONT，ANSYS根据模型尺寸为ICONT提供一个小的缺省值。关闭ICONT，必须将其设置成非常小的值(如1e-20)0值代表非0的缺省值计算实例：刚体ICONT调整初始接触条件（PMIN和PMAX），见图3.68。实常数PMIN和PMAX指定初始穿透范围。ANSYS把整个目标面（连同变形体）移到由PMIN和PMAX指定的穿透范围内。调整初始接触条件（ICONT）如果目标面无约束，采用PMIN和PMAX的初始调整将不被执行。初始调整是一个迭代过程，ANSYS最多使用20个迭代步把目标面调整到PMIN和PMAX范围内。当目标面被移动，接触体之间不再存在间隙，而成为闭合接触的初始状态。调整初始接触条件（PMIN和PMAX）计算实例：刚体PMIN&PMAX第六节第六节 几种典型接触问题分析实例几种典型接触问题分析实例例2.悬臂梁端部旋转（“不分离”行为）目标：采用“不分离”行为选项施加大转动模拟特定物理效果。不分离：目标面和接触面一旦建立接触就不再分离（允许滑动）。模型描述：悬臂梁，端部施加旋转，见图文件：no_separation.inp采用Solid185单元：3D 8 node structural Solid单元，具有大变形能力。图 悬臂梁端部旋转/view,1,1,1/prep7k,1,10,0,0BLOCK,0,10,-0.5,0.5,-0.5,0.5RECTNG,9.5,10,-1,1AGEN,7,-0.5,1ARSYM,Z,7et,1,185mp,ex,1,10000esize,0.5vmesh,1ESIZE,2 !为目标面分网用/solutionnsel,loc,x,0d,all,all,0/pbc,u,1finisheplotSAVE,no-separation,db以下用接触向导建立接触对:例2.紧配合环（初始穿透）目标：验证接触初始穿透选项。采用缺省值和渐进化穿透选项求解一个大的初始穿透接触分析，比较结果。模型模述：紧配合环、轴对称，见图 图 紧配合环（初始穿透）文件：Interference.inp/PREP7ET,1,183,1 !轴对称，大应变ET,2,169 !Target 169 2D目标单元ET,3,172 !Conta 172 2D 3节点高阶抛物线接触单元MP,EX,1,10E3MP,NUXY,1,0MP,MU,2,0 !无摩擦RECTNG,1.5,2,0,1RECTNG,1.9,2.4,0,1ESIZE,0.05AMESH,ALL !Type,1;Mat,1TYPE,2MAT,2REAL,2LSEL,LOC,X,1.9NSLL,1 !1.9线上的节点ESURF !创建目标单元Target169!*TYPE,3 !Mat,2;Real,2LSEL,LOC,X,2NSLL,S,1 ESURF !创建接触面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,db1 !包含有限元模型、接触对、材料、边界条件，求解控制例例3.两刚性面间夹圆柱体分析两刚性面间夹圆柱体分析（刚体模式应用初始接触环ICONT）目标：验证使用初始接触环调整带来克服接触分析中的刚体模式模型描述：2D平面应力、1/4对称、施加压力，见图 文件名:ICONT.inp 有限元分析模型两刚性面间夹圆柱体分析/prep7CYL4,0,0,0,0,1,90RECTNG,0,2.5,1.003,1.5ET,1,182!Plane182 2D Structural Solid 4节点UX,UY平面应力、平面应变和轴对称单元、大变形、大应变MP,EX,1,1000MP,EX,2,30E6mp,nuxy,2,0.3ESIZE,0.1AMESH,1MAT,2AMESH,2SAVE,ICONT,db以下用接触向导创建接触对第四章第四章 点面接触点面接触点面接触单元是90年代普遍使用的接触单元。由于点面接触单元理论上的限制，使它们被更新更好的面一面接触单元取代。点一面接触单元可以用来模拟一个表面和一个节点的接触；也可以把表面指定为一组节点，用点面接触单元来模拟面一面的接触。面一面接触单元处理角点接触有困难，因为它们采用高斯点作为接触检查点，在角点处会呈现过渡穿透。在此情况下，可以混合使用面一面接触单元和点面接触单元，见图 面一面接触单元和点面接触单元混合点面接触单元不必知道接触面的位置。允许大变形，大的相对滑动，库仑摩擦滑动；接触面间可用不同的网格划分。点面接触是通过跟踪一个表面（接触面）上的点相对于另一表面（目标面）上的线或面的位置来表示的，程序使用接触单元来跟踪两个面的相对位置。接触单元形状为三角形、四面体或椎体，其底面由目标面上的节点组成，而顶点为接触面上的节点，见图 点面接触单元在节点传递力（面面接触单元在高斯点传递力）此特性使其只能用于低阶单元（角节点）这是由于中间节点的单元节点上的反力不均匀，见图 点面接触单元第一节 接触刚度点面接触单元（conta48、49）要求给出罚刚度。可以通过实验来确定一个合适的接触刚度，使求解收敛而且侵入量可以接受。选择接触刚度：对于块状实体，通常赫芝接触刚度适用于罚刚度，可以这样来估算：K=fE式中：f=0.110系数E=较软的接触体材料的弹性模量设f=1通常是一个较好的起始值。对于柔性体（梁和壳模型），系统的刚度可以比赫芝接触刚度低很多。此时可以将单位载荷施加到要接触的面上，先运行一个静态分析来确定模型的局部刚度，接触刚度可以这样来估算：上式适用于柔体接触，f=1100系数，设f=1是一个比较好的起始值。第二节 点面接触分析步骤建模与分网识别接触对生成接触单元设置单元关键字（Keyopt）和实常数给定边界条件定义求解选项求解查看结果NSEL,S,NODE,!在接触面上选择一组节点CM,CONTACT,NODE !将所有节点定义成组元“CONTACT”NSEL,S,NODE,!在目标单元上选择一组节点CM,TARGET,NODE,!将所选节点定义成组元“TARGET”NSEL,ALL !选中所有节点E,!设置单元类型R,!选择实常数!*生成接触单元*可采用接触向导施加载荷、设定边界条件建模时使接触体处于恰好的接触位置使用给定的位移将它移到某个位置接触分析中加载、设定边界条件方法与步骤和其它非线性分析相同。定义求解选项点面接触分析中常用求解设置及注意事项：时间步长必须足够小，如果时间步长太大，接触力的光滑传递将被破坏。为确保结果的准确性，可以打开自动步长（Autots,on）设置一个合适的平衡迭代次数：NEQIT,2575许多接触分析不收敛是因为设置的接触刚度太大（实常数KN取值太大）造成的，这时需要减小接触刚度重新进行分析。后处理后处理接触分析的结果主要包括位移、应力、应变和接触信息。接触信息包括：接触压力、单元的现在和过去状态：分开（没有接触）接触粘合状态接触滑动状态粘合=1；滑动=2或-2；分开=3或4两个表面间的距离，如果是正值，两表面是分开的（STAT=3或4），如果是负值代表穿透量（STAT=1或2）通过动画显示接触结果随时间的变化规律是接触分析有效的、常用的处理方法。第三节 点面接触问题分析实例例：梁端部接触目标：验证采用点-面接触单元模拟梁端部接触。建立2D点-面接触单元，求解大变形接触分析并进行后处理。模型描述：悬臂梁施加端部位移(见图)。文件:node_to_surface.inpPREP7ET,1,BEAM3B=0.5H=0.5R,1,B*H,B*(H*3)/12,HMP,EX,1,30e6K,1,0,0,0$K,2,100,0,0K,3,0,1,0$K,4,75,-10$K,5,175,-10L,1,2$L,4,5LESIZE,ALL,20LATT,1,1,1,3LMESH,ALLFINISH/SOLUTIONDK,1,ALL,0$DK,5,ALL,0$DK,2,UY,-30/PBC,U,1/PBC,ROT,1FINISHEPLOTSAVE,node-to-surface,dbnode-to-surface.db 包括有限元模型几何，材料，边界条件、加载、分网刚度估计按式K=f(bend)E/10F(bend)=1则K=1*30e6/10=30e5 选择接触点GUI：Utility Menu Select EntitiesNodesBy Num/Pick选择“From Full”【OK】拾取节点2【OK】或命令：Nsel,2建立接触节点组元GUI：Utility Menu Select Comp/Assembly Create ComponentComponent name=CONTACTComponent is made of=nodes【OK】命令：CM,CONTACT,node选择目标面节点GUI：Utility Menu Select EntitiesNodesBy Num/Pick 选择“From Full”【OK】选下面梁上节点（共21个）【OK】或命令：Nsel,22,42建立目标节点组元GUI：Utility Menu Select