毕业设计（论文）-弹簧管刚度测量的有限元研究与分析.doc

大学学生毕业论文论文题目： 弹簧管刚度测量的有限元研究与分析 学 院： 年 级： 专 业： 姓 名： 学 号： 指导教师： 2007年5月28日摘要在航天、航海、航空领域常用的控制器件电液伺服阀和导向器件挠性陀螺中广泛地使用精密弹性元件作为力敏单元获得反馈信号。这些弹性元件的刚度参数对相对应的控制器件、导航器件的性能影响很大。因此，在加工过程中不仅对其尺寸精度和几何形状精度提出很高的要求，而且对其刚度特性也提出了一定的要求。由于这类弹性元件的几何形状较复杂，并且加工误差和材质变化均会影响到刚度值，因此生产中是通过测量的手段来保证刚度指标。随着制造技术水平的提高以及航空航天技术的不断发展，对精密弹性元件的加工精度和加工效率提出了更高的要求，所以研究以弹簧管为代表的该类弹性元件的刚度精密测量技术对增强国防实力有着重要意义。本文在哈工大研制的gcy-型测量仪的基础上，分析弹簧管头部与测量孔的相互作用，最后利用ansys建模重点分析二者之间的接触变形。首先，给出弹簧管刚度测量原理，然后，分析弹簧管头部的受力并将测量力所引起的接触力简化为均匀分布后建立接触模型，最后，利用ansys建模分析接触变形。关键词弹簧管； 接触变形；刚度abstractprecise elastic parts are widely used as sensitive parts to gain feed back signals in the controlling instrument such as electro-hydraulic servo valve or the navigation system such as flexibleness gyro in the field of astronavigation, navigation and aviation. their stiffness parameters directly correlated with the performances of corresponding system .therefore ,not only have a rather high standard been set for dimensional precision and geometrical precision of them ,but also certain demands are made against theirs stiffness characteristic. due to their complicated geometry, along with machining errors and changes of materials characteristic, it is impossible to gain the accurate values of the stiffness by calculation .so, measurement is the only way to assure the stiffness requirement in production. at present, as the production level is advanced with the technology of astronavigation and continuous development of aviation, higher accuracy and efficiency in measurement are demanded to meet them. therefore, studies on the stiffness measuring of spring tubes, typical plastic parts, is of great importance in enhancing strength of national defense.referring to hit measuring instruments, gcy-, analyzing the behavior between head of spring and hole, and, using ansys to analyze contact deformation between them. firstly the measurement principal of spring .secondly, the stress analysis in measurement process is discussed and the contact force caused by measurement force is simplified to be symmetrically distributed along the contact line .lastly, use ansys to analyze the contact deformation.key wordsspring tube; contact deformation; stiffness目录摘要iabstractii第一章 绪论11.1 课题的背景及意义11.2 本课题的主要研究内容1第二章 弹簧管刚度测量原理2 2.1 弹性元件刚度定义22.2 弹簧管刚度定义22.3 弹簧管测量原理3本章小结4第三章 弹簧管头部受力分析与建模53.1 弹簧管头部受力分析53.2 建立实体模型7本章小结7第四章 有限元分析弹簧管接触变形84.1 有限元概述84.1.1有限元法简介84.1.2有限元法常用术语94.2接触分析概述94.2.1接触问题分类94.2.2 ansys中关于接触问题的设置104.2.3接触分析的步骤114.2.4 ansys分析29总结36参考文献37致谢3837弹簧管刚度测量的有限元研究与分析第一章 绪论1.1 课题的背景及意义电液伺服系统以其响应速度快、输出功率大、控制精确性高等诸多优点, 被广泛的应用于航空、航天、军事、冶金、交通、工程机械等领域。电液伺服阀是电液控制系统的关键部件, 用于电液伺服系统的位置、速度、加速度和力的控制。其中, 弹簧管作为电液伺服阀中的关键部件, 其刚度特性直接影响控制系统的控制精度和灵敏度。弹簧管是一种薄壁弹性零件, 薄壁处的厚度仅有卜, 加工中的误差和材质变化等因素会直接影响其刚度特性, 无法对其刚度特性建立起精确的数学模型, 只能通过测量来保证其刚度的设计指标。因此，只能通过加工-测量-加工的方法来保证其精度。由于弹簧管的精度要求很高，而实际测量过程中存在着许多误差，如：测量原理误差，制造和装配误差，测量中各个环节的弹性变形和接触变形引起的误差等等。在大刚度值弹簧管的刚度测量中，重复测量精度较差，通过大量的实验发现弹簧管的头部与测杆孔之间的接触变形对测量结果有很大的影响。在大刚度值弹簧管的测量中因为施加的载荷较大，因此接触变形的作用就变得十分明显。接触变形属于非线性行为，根据弹簧管和测杆之间的接触特点，研究影响弹簧管和测杆之间的接触变形大小的因素，为提高测量系统的重复测量精度打下基础。1.2 本课题的主要研究内容（1）给出弹簧管刚度测量原理。（2）对弹簧管头部进行受力分析，并将有测量力引起的接触力简化为均匀线分布，建立接触分析模型。（3）利用ansys分析接触变形。第二章 弹簧管刚度测量原理 2.1 弹性元件刚度定义 弹性元件的刚度定义为在外加载荷作用下元件抵抗变形的能力。刚度值k等于外加载荷q与其作用下弹性元件的变形量d之比值，即： (2-1)对于一个简化为悬臂梁的弹性体，其外载荷分为力和力矩，而相应的端部位移有直线挠度和角位移之分，因此，其刚度也就有四种刚度形式，如图2-1所示。(a) (b)图2-1 简化悬臂梁受力,力矩示意图 对于图2-1（a）所示在简化悬臂梁的一段施加垂直于杆端初始方向的力f，图2-1（b）所示施加力矩m，由此在悬臂梁端部产生了挠度，并对应变形转角，的到了悬臂梁的四种刚度定义。 抗力位移刚度 (2-2) 抗力转角刚度 (2-3) 抗力矩位移刚度 (2-4) 抗力矩转角刚度 (2-5)2.2 弹簧管刚度定义根据力矩马达实际工作受力和变形，将弹簧管单独分析，受力如图2-3所示：由于弹簧管的水平微小移动，对伺服阀的性能基本没有影响，因此，我们只关心的是弹簧管的转角， (2-6) 图2-2弹簧管受力示意图2.3 弹簧管测量原理图2-3两点法测量原理图具体的测量方法是，将弹簧管装在卡具中卡紧，并在弹簧管头上安装测量杆，与测量反馈杆同步测出，其中。同样，根据所得到的n个测量点的力值和位移值，利用最小二乘法可以算出弹簧管的刚度值: 力矩： (2-7) 转角： (2-8) 当时，有： (2-9)用多点测采量，运用最小二乘法计算： (2-10)本章小结 本章主要是研究了弹簧管的测量原理，确定了弹簧管的刚度是什么，为什么被研究量是转角以及如何来测量这个刚度，为后文的分析打下基础。第三章 弹簧管头部受力分析与建模3.1 弹簧管头部受力分析弹簧管的刚度测量中，弹簧管头部与测杆孔之间的接触如图3-1所示。测量时，将弹簧管3的底端的法兰盘固定在夹具（图中未画出）中，将其头部装入刚度远远大于弹簧管的测杆2中，并由限位销4保证受力点的位置；用力矩扳手旋紧顶杆6，使之向左运动通过顶座7将弹簧管头部夹紧。在刚度测量前，弹簧管头部与测杆孔在顶紧力的作用下保持接触，有一定的接触变形，测量中在测量力f的作用下使二者之间产生进一步的变形。综上，导致二者之间发生接触变形的作用力包括两部分，一是测量力f，一是锁紧力。通过作用力分析，给出在f、作用下弹簧管与测杆之间的线接触力的分布如图4-2所示。由测量力f作用产生的线接触力呈线性分布，分布函数为：其中由式(3-1)给出。顶紧力引起来的线接触力为均匀分布，单位分布由式(3-2)给出。 图3-1 弹簧与测杆接触示意图 (3-1) (3-2)式中 弹簧管头部的高度； 任意测量点处的测量力，由测力传感器给出； 每次测量由拧紧力矩产生的顶紧力，由式（3-3）给出。 （3-3）式中 是顶杆螺纹中径； 是顶杆螺纹升角； 是顶杆的当量摩擦角； 是拧紧力矩。 图3-2 接触面受力分析弹簧管刚度测量中保证弹簧管的可靠夹紧是获得稳定测量的一个必要条件。理论上，设计的测杆孔的直径与被测型号弹簧管头部的直径越接近越能够保证夹紧，但是弹簧管头部在加工过程中是有误差的，为了保证同一型号的弹簧管都能够比较方便的装入测杆中，在测杆的设计中根据弹簧管头部的尺寸公差，取研磨后测杆孔径与弹簧管头部的轴径之差?r preprocessormeshmapped3 or4 sided 步骤二：识别接触对必须认识到模型在变形期间哪些地方可能发生接触，一是已经识别出潜在的接触面，该通过目标单元和接触单元来定义它们，目标和接触单元跟踪变形阶段的运动，构成一个接触对的目标单元和接触单元通过共享的实常号联系起来。接触环（区域）可以任意定义，然而为了更有效的进行计算（主要指cpu时间）可能想定义更小的局部化的接触环，但能保证它足以描述所需要的接触行为，不同的接触对必须通过不同的实常数号来定义（即使实常数号没有变化）。由于几何模型和潜在变形的多样形，有时候一个接触面的同一区域可能和多个目标面产生接触关系。在这种情况下，应该定义多个接触对（使用多组覆盖层接触单元）。每个接触对有不同的实常数号。步骤三：定义刚性目标面 刚性目标面可能是2d的或3d的。在2d情况下，刚性目标面的形状可以通过一系列直线、圆弧和抛物线来描述，所有这些都可以用tapge169来表示。另外，可以使用它们的任意组合来描述复杂的目标面。在3d情况下，目标面的形状可以通过三角面，圆柱面，圆锥面和球面来推述，所有这些都可以用tapge170来表示，对于一个复杂的，任意形状的目标面，应该使用三角面来给它建模。控制结点（pilot） 刚性目标面可能会和pilot结点联系起来，它实际上是一个只有一个结点的单元，通过这个结点的运动可以控制整个目标面的运动，因此可以把pilot结点作为刚性目标的控制器。整个目标面的受力和转动情况可以通过pilot结点表示出来，“pilot结点”可能是目标单元中的一个结点，也可能是一个任意位置的结点，只有当需要转动或力矩载荷时，“pilot结点”的位置才是重要的，如果定义了“pilot结点”ansys程序只在“pilot结点”上检查边界条件，而忽略其它结点上的任何约束。对于圆、圆柱、圆锥、和球的基本图段，ansys总是使用条一个结点作为“pilot结点”，基本原型使用基本几形状来模拟目标面，例如：“圆、圆柱、圆锥、球。直线、抛物线、弧线、和三角形不被允许、虽然不能把这些基本原型彼此合在一起，或者是把它们和其它的目标形状合在一起以便形成一个同一个实常数号的复杂目标面。但可以给每个基本原型指定它自己的实常数号。 单元类型和实常数 在生成目标单元之前，首先必须定义单元类型（targ169或targ170）。gui: main menupreprocessorelement type add/edit/delete 随后必须设置目标单元的实常数。 gui: main menupreprocessorreal constants 对targe169和targe170仅需设置实常数r1和r2，而只有在使用直接生成法建立目标单元时，才需要从为指定实常数r1、r2，另外除了直接生成法，可以使用ansys网格划分工具生成目标单元，下面解释这两种方法。用直接生成法建立刚性目标单元为了直接生成目标单元，使用下面的命令和菜单路径。 gui：main menupreprocessormodeling-createelementselem attributes 随后指定单元形状，可能的形状有： straight line (2d) parabola (2-d) clockwise arc (2-d) counter clock wise arc (2-d) circle (2-d) triangle (3-d) cylinder (3-d) cone (3-d) sphere (3-d) pilot node (2-d和3-d) 一旦指定目标单元形状，所有以后生成的单元都将保持这个形状，除非指定另外一种形状。然后就可以使用标准的ansys直接生成技术生成结点和单元。gui：main menupreprocessor modeling- create nodes 在建立单元之后，可以通过列示单元来验证单元形状 gui：utility menulistelementsnodes +attributes 可以使用标准的ansys网格划分功能让程序自动地生成目标单元。gui：main menupreprocessormeshing-meshkey points 为了生成一个2d目标单元，使用下面的命令和gui路径： ansys在每条直线上生成一条单一的线，在样条曲线上生成抛物线部分，在每条圆弧和倒角上生成圆弧部分，如果所有的圆弧形成一个封闭的圆，ansys生成一个单一的圆段。gui：main menupreprocessormeshing-meshlines 为了生成3d的目标单元，使用下面的命令或gui路径。如果实体模型的表面部分形成了一个完整的球，圆柱或圆锥，那么ansys程序自动生成一个基本的3d目标单元，因为生成较少的单元，从而使分析计算更有效率，对任意形状的表面，网格形状的质量不是重要的，而目标单元的形状是否能完成好的模拟刚性面的表面几何形状显得更重要。 gui：main menupreprocessor-meshing-mesharea ansys在所有可能的面上推荐使用三角形的映射网格划分，如果在表面的边界上没有曲率，则在网格划分时，指定那条边界分为一分，下面的命令或gui路径将尽可能的生成一个映射网格（如果不能进行映射，它将生成自由网格） gui：main menupreprocessor-meshing-mesh-ares-target surf 建模和网格划分的注意点： 一个目标面可能由两个或多个面断的区域组成，你应该尽可能地通过定义多个目标面来使接触区域局部比（每个目标面有一个不同的实常数号）刚性目标面上由的离散能足够指述出目标面的形状，过粗的网格离散可能导致收敛问题。如果刚性面有一个实的凸角，求解大的滑动问题时很难获得收敛结果，为了避免这些建模问题，在实体模型上，使用线或面的倒角来使尖角光滑比，或者在曲率突然变化的区域使用更细的网格。不能使用镜面对称技术（arsysm，lsymm）来映射圆、圆柱、圆锥或球面到对称平面的另一边，因为每个实常数的设置不能同时赋给多个基本原型段。检验目标面的接触方向目标面的结点号顺序是重要的，因为它定义了接触主向，对2d接触问题，当沿着目标线从第一个结点移向第二个结点时，变形体的接触单元必须位于目标面的右边。 对3d接触问题，目标三角形单元号应该使刚性面的外法线方向指向接触面，外法线通过右手原则来定义为了检查法线方向，显示单元坐标系 gui：utility menupilot curlssymbols 如果单元法向不指向接触面，选择单元反转表面的法向的方向。 gui：main menupreprocessorcreateelementon free surf 步骤四：定义柔性体的接触面为了定义柔性体的接触面，必须使用接触单元confa171或confa172（对2d）或conta173或conta174（对3d）来定义表面，程序通过组成变形体表面的接触单元来定义接触表面，接触单元与下面覆盖的变形体单元有同样的几何特性，接触单元与下面覆盖的变形体单元必须处于同一阶次（低阶或高阶）下面的变形体单元可能是实体单元、壳单元、梁单元或超单元，接触面可能壳或梁单元任何一边。 与目标面单元一样，你必须定义接触面的单元类型，然后选择正确的实常数号（实常数号必须与它对应目标的实常数号相同）最后生成接触单元。下面简单描述四种类型的接触单元：conta171：这是一种2d，2个结点的低附线单元，可能位于2d实体，壳或梁单元的表面 conta172：这是一个2d的，3结点的高阶抛物线形单元，可能位于有中结点的2d实体或梁单元的表面 conta173：这是一个3d的，4结点的低阶四边形单元可能位于3d实体或壳单元的表面，它可能退化成一个结点的三角形单元。 conta174：这是一个3d，8结点的高阶四边形单元，可能位于有中结点的3d实体或壳单元的表面，它可能退化成6结点的三角形单元。不能在高阶柔性体单元的表面上分成低阶接触单元，反之也不行，不能在高阶接触单元上消去中结点。gui：main menupreprocessorelement typeadd/edit/delete 实常数和材料特性 在定义了单元类型之后，需要选择正确的实常数的设置，每个接触对的接触面和目标面必须有相同的实常数号，而每个接触对必须有它自己不同的实常数号。ansys使用下面柔性体单元的材料特性来计算一个合适的接触（或罚）刚度，如果下面的单元是一个超单元，接触单元的材料的设置必须与超单元形成时的原始结构单元相同，生成接触单元。既可以通过直接生成法生成接触单元，也可以在柔性体单元的外表面上自动生成接触单元，推荐采用自动生成法，这种方法更为简单和可靠。可以通过下面三个步骤来自动生成接触单元。 1、选择结点选择已划分网格的柔性体表面，如果某一部分结点永远不会接触到目标面，可以忽略它以便减少计算时间，然而，必须保证设有漏掉可能会接触到目标面的节点。gui：main menupreprocessorcreateelementonfree surf 2、生成接触单元gui：main menupreprocessorcreateelementon free surf 如果接触单元是附在已用实体单元划分网格的面或体上，程序会自动决定接触计算所需的外法向，如果下面的单元是梁或壳单元，则必须指明哪个表面（上表面或下表面）是接触面。gui：main menupreprocessorcreateelementon free surf 使用上表面生成接触单元，则它们的外法向与梁或壳单元的法向相同，使用下表面生成接触单元，则它们的外法向与梁或壳单元的法向相反，如果下面的单元是实体单元，则top或bottom选项不起作用 3、检查接触单元外法线的方向，接触面的外法线方向是重要的，对3d单元，按结点程序号以右手定则决定单元的外法向，接面的外法向应该指向目标面，否则，程序可能会认为有面的过度渗透而很难找到初始解。gui：utility menuplot controlssymbols 当发现单元的外法线方向不正确时，必须通过倒不正确单元的结点号来改变它们。gui：main menupreprocessorcreateelements on free surf 步骤五：设置实常数和单元关键字程序使用九个实常数和好几个单元关键字来控制面面接触单元的接触行为。1、实常数： 9个实常数中，两个（r1和r2）用采定义目标面单元的几何形状，乘下的7个用来控制接触行为。 r1和r2定义目标单元几何形状 fkn 定义法向接触刚度因子 ftoln 定义最大的渗透范围 icont 定义初始靠近因子 pinb 定义“pinball区域 pmin和pmax 定义初始渗透的容许范围 taumar 指定最大的接触摩擦 gui：main menu preprocessorreal constant 对实常数fkn，ftoln，icont，pinb，pmax，和pmin，既可以定义一个正值也可以定义一个负值，将正值作为比例因子，将负值作为真实值，程序将下面覆盖原单元的厚度作为icon，ftoln，pinb，pmax和pmin的参考值，例如对icon，0.1表明初始间隙因子是0.1*下面覆盖层单元的厚度。然而，-0.1表明真实缝隙是0.1如果下面覆盖层单元是超单元，则将接触单元的最小长度作为厚度。2、单元关键字 每种接触单元都默认好几个关键字，对大多的接触问题缺省的关键字是合适的，而在某些情况下，需要改变缺省值，来控制接触行为。接触算法（罚函数+拉格郎日或罚函数）（keyopt（2） 出现超单元时的应力状态（deyopt（3） 接触方位点的位置 （keyopi（4） 刚度矩阵的选择 （keyopt（6） 时间步长控制 （keyopt（7） 初始渗透影响 （keyopt（9） 接触表面情况 （keyopt（12） gui：main menupreprocessorelement typeadd/edit/delete 3、选择接触算法： 对于面面的接触单元，程序可以使用扩增的拉格朗日算法或罚函数方法，通过使用单元关键字ketopt（2）来指定。扩张的拉格朗日算法是为了找到精确的拉格朗日乘子而对罚函数修正项进行反复迭代，与罚函数的方法相比，拉格朗日方法不易引起病态条件，对接触刚度的灵敏度较小，然而，在有些分析中，扩增的拉格朗日方法可能需要更多的迭代，特别在变形后网格变得太扭曲时。用拉格朗日算法的同时应使用实常数。ftoln为拉格朗日算法指定容许的最大渗透，如果程序发现渗透大于此值时，即使不平衡力和位移增量已经满足了收敛准则，总的求解仍被当作不收敛处理，ftlon的缺省值为0.1，如果此值太小则可能会造成太多的迭代次数或者不收敛。4、决定接触刚度 所有的接触问题都需要定义接触刚度，两个表面之间渗透量的大小取决了接触刚度，过大的接触刚度可能会引起总刚矩阵的病态，而造成收敛困难，一般来说，应该选取足够大的接触刚度以保证接触渗透小到可以接受，但同时又应该让接触刚度足够小以使不会引起总刚矩阵的病态问题而保证收敛性。程序会根据变形体单元的材料特性来估计一个缺省的接触刚度值，可以用实常数fkn来为接触刚度指定一个比例因子或指定一个真正的值，比例因子一般在0.01和10之间，当避免过多的迭代次数时，应该尽量使渗透到达极小值。为了取得一个较好的接触刚度值，可以按下面的步骤过行。1、开始时取一个较低的值，低估些值要比高估些值好，因为由一个较低的接触刚度导致的渗透问题要比过高的接触刚度导致的收敛性困难，要容易解决。2、对前几个子步进行计算 3、检查渗透量和每一子步中的平衡迭代次数，如果总体收敛困难是由过大的渗透引起的（而不是由不平衡力和位移增量引起的），那么可能低估了fkn的值，或者是将ftoln的值取得大小，如果总体的收敛困难是由于不平衡力和位移增量达到收敛值需要的迭代次数，而不是由于过大的渗透量，那么fkn的值可能被高估。4、按需要调查fkn或ftoln的值，重新分析。5、选择摩擦类型。 在基本的库仑摩擦模型中，两个接触面在开始相互滑动之前，在它们的界面上会有达到某一大小的剪应力产生，这种状态则作粘合状态（stick)库仑摩擦模型定义了一个等效剪应力。）一旦剪应力超过此值后，两个表面之间将开始相互滑动，这种状态，叫做滑动状态（sliding）粘合滑动计算决定什么时候一个点从粘合状态到滑动状态或从滑动状态变到粘合状态，摩擦系数可以是任一个非负值。程序缺省值为表面之间无摩擦，对rough或bonded接触（keyopt（2）=1（或3），程序将不管给定的mv值而认为摩擦阻力无限大。程序提供了一个不管接触压力的故而人为指定最大等效剪应力的选项，如果等效剪应力达到此值时，滑动发生。如图，4-1，为了指定接触界面上最大许可剪应力，设置常数taumax（缺省为1.0e20），这种限制剪应力的情况一般用于接触压力非常大的时候，以至于用库仑理论计算出的界面剪应力超过了材料的屈服极限。一对taumax的一个合理高估为y3（y3是材料的misses屈服应力）图4-1 摩擦模式对无摩擦rough和bonded接触，接触单元刚度矩阵是对称的，而涉及到摩擦的接触问题产生一个不对称的刚度，而在每次迭代使用不对称的求解器比对称的求解器需要更多的计算时间，因此ansys采用对称化算法。通过采用这种算法大多的摩擦接触问题能够使用对称系统的求解器来求解。如果摩擦应力在整个位移范围内有相当大的影响，并且摩擦应力的大小高度依赖于求解过程。对刚度阵的任何对称近似都可能导致收敛性的降低，在这种情况下，选择不对称求解选项（keyopt（6）=1）来改善收敛性。6、选择检查接触与否的位置 接触检查点位于接触单元的积分点上，在积分点上，接触单元不渗透进入目标面，然而，目标面能渗透进入接触面，看图4-2。 图4-2 接触检查点位于高斯积分点上ansys面接触单元使用gauss积分点作为缺省值，gauss积分点通常会比newton-cotes/robots结点积分项产生更精确的结果，newton-cotes/locator使用结点本身作为积分点，通过keyopt（4）来选择，然而，使用结点本身作为积分点仅应该用于角接触问题（看图4-3）。 图4-3 接触检查点位于高斯结点上使用结点作为接触发现点，可能会导致其它的收敛性问题，例如“滑脱”（结点滑下目标面的边界）看图4-4，对大多的点面的接触问题，推荐使用其它的点面的接触单元，例如conta26、conta48和conta49 图4-4 结点滑脱7、调整初始接触条件 在动态分析中，刚体运动一般不会引起问题，然而在静力分析中，当物体没有足够的约束时会产生刚体运动，有可能引起错误而终止计算。在仅仅通过接触的出现来约束刚体运动时，必须保证在初始几何体中，接触对是接触的，换句话说，要建立模型以便接触对是“刚好接触”的，然而这样作可能会遇到以下问题： 刚体外形常常是复杂的，很难决定第一个接触点发生在哪儿 既使实体模型是在初始接触状态，在网格划分后余于数值舍入误差；两个面的单元网格之间也可能会产生小的缝隙。接触单元的积分点和目标单元之间可能有小的缝隙。 同理，在目标面和接触面之间可能发生过大的初始渗透，在这种情况下，接触单元可能会高估接触力，导致不收敛或民接触面之间脱离开接触关系。定义初始接触也许是建立接触分析模型时最重要的方面，因此，程序提供了几种方法来调整接触对的初始接触条件。为了消除由于生成网格造成的数值舍入误差而引起的小缝隙或渗透，可以采用以下方法。1、使用实常数icont来指定一个好的初始接触环，初始接触环是指沿着目标面的“调整环”的深度，如果没有人为指定icont的值，程序会根据几何尺寸来给icont提供一个小值，同时输出一个表时什么值被指定的警告信息，对icont一个正值表示相对于下面变形体单元厚度的比例因子，一个负值表示接触环的真正值，任何落在“调整环”敬域内的接触检查点被自动移到目标面上，（看图4-5(a)）可以使用一个小的icont值否则，可能会发生大的不连续（看图4-5(b)） 图4-5 用icon进行接触面的调整(a) 调整前 (b) 调整后2、使用实常数pmin和pmax来指定初始容许的渗透范围，当指定pmax或pmin后，在开始分析时，程序会将目标面移到初始接触状态，如果初始渗透大于pmax，程序会调整目标面的减少渗透，接触状态的初始调节仅仅通过平移来实现。对给定载或给定位移的刚性目标面将会执行初始接触状态的初始调节。同样，对没有指定边界条件的目标面也可以进行初始接触的调整。当目标面上的节点，有给的零位移值时，使用pmax和pmin的初始调节将不会被执行。图46给出了一个初始接触调整迭代失败的例子。目标面的uy被约束住。因此，初始接触唯一容许的调整是在x方向，然而，在这个问题中，刚性目标面在x方向的任何运动都不会引起初始接触。 图4-6 初始调整失败3、设置keyopi（9）=1来消除初始渗透，看图4-7。 图4-7 消除初始渗透在某些情况下，例如过盈装配问题，期望有过度的渗，为了缓解收敛性困难，在第一个载荷步中设置keyoti（9）=2来使过度渗透渐进到0，看图4-8。当使用这种方法时，在第一个载步中不要给定其它任何载荷，也就是说要保证载荷是渐进的（kbc，0） 图4-8 渐进初始渗透在开始分析时，程序会给出每个目标面的初始接触状态的输出信息，（在输出窗口或输出文件中），这个信息有助于决定每个目标面的最大渗透成最小间隙。对于给定的目标面如果没有发现接触，可能是目标面离接触面太远（超出了pinball区域或者是接触/目标单元已经被杀死。8、决定接触状态和pinball区域。接触单元相对于目标面的运动和位置决定了接触单元的状态；程序检测每个接触单元并给出一种状态stat=0 未合的远区接触stat=1 未合的近区接触stat=2 滑动接触stat=3 粘合接触当目标面进入pinball区域后，接触单元就被当作未合上的近区域接触，pinball区域是以接触单元的积分点为中心的。使用实常数pinb来为pinball指定一个比例因子（正值）或其实值（负值），缺省时，程序将pinball区域定义为一个以4*变形体单元厚度为半径的圆（对2-d问题）或球（对3-d问题）。检查接触的计算时间依赖于pinball区域的大小，远区接触单元的计算是简单的且计算时间较少，近区接触计算将要接触的接触单元是较慢的,并是较复杂，当单元已经接触时，计算最为复杂。如果刚性面有好几个凸形区域，为了克服伪接触定义