接触分析实例
ABAQUS应用培训-09 接触分析及实例
15
2015/7/5
15
4
2015/7/5
4
接触 –利用接触模块定义并管理以下信息：
– 模型区域之间或模型区域同它周围模型之间的机械和热接触。
机械接触 一般和双面接触 自接触 弹性基础 激励器/传感器
其它接触 (约束) 捆绑约束 刚体 显示体 耦合 壳到实体耦合 嵌入区域 线性多点约束方程
热接触 热膜条件 对环境的辐射 热接触
1
4
5
2 7
3 10
6 11
8
9
12
3
接触操作（Interaction）
Interaction
1. 不同的接触算法：standard explicit 2. 作为一种边界条件，接触定义在不同的step中可以发生变化 3. 接触算法还分为Kinematic, Penalty，初学者不用选择 4. Surface to Surface, Surface to node 5. 主面（master, first）和从面（slace,second） 6. 接触体之间的相互作用：摩擦和热传递等设置 7. 接触面的间隙调整
k14 k24 k34 k44 k54 k64 0 0 0 0 0 0
0 0 k35 k45 k55 k65 0 0 0 0 0 0
0 0 k36 k 46 k56 k66 0 0 0 0
k116 0
0 0 0 0 0 0 k77 k87 k97 k107 0 0
0 0 0 0 0 0 k78 k88 k98 k108 0 0
有限元分析-ABAQUS的应用
接触分析
方 刚
2015年7月
2015/7/5
1
1
接触
• “为什么要定义机械接触？
Ansys workbench 螺栓接触实例操作
8例1 螺栓连接件分析如图所示为一螺栓连接的法兰连接件简图，法兰一端及内侧面固定约束。
载荷1为螺栓预应力1000N载荷2为螺栓预应力1500N载荷3为螺栓预应力2000N根据实际情况，自己设定接触类型，其中摩擦类型接触对时，摩擦系数为0.1 为方便设置，材料均取钢材，求其变形及应力。
边界条件螺栓连接件分析1 导入几何模型，进入DS模块2 材料设置选择默认的材料：Structural Steel3 设置接触螺栓与螺母的接触类型为Bonded螺栓杆与法兰的接触类型为Frictional，摩擦系数为0.1螺栓杆与垫片内壁的接触类型为Frictional，摩擦系数为0.1其余接触类型为No Separation4 网格划分5 选择分析类型·在“New Analysis”中选择结构静力学分析“Static Structural”；6 施加约束与载荷1）施加固定约束·点击“Static Structural”，在“Supports”中选择固定约束“Fixed Support”·选择法兰一端及内侧面固定约束；2）施加载荷·选择载荷1处螺栓杆表面，添加螺栓预应力“Bolt Pretension”大小为1000N ·选择载荷2处螺栓杆表面，添加螺栓预应力“Bolt Pretension”大小为1500N ·选择载荷3处螺栓杆表面，添加螺栓预应力“Bolt Pretension”大小为2000N5 设定求解类型1）求解变形·点击“solution”，点击“Deformation”选择“Total”，求解变形·点击“Stress”，选择“Equivalent (V on-Mises)”，求解等效应力6 单击“Solve”求解7 观察求解结果·点击“Total Deformation”查看变形·点击“Equivalent Stress”查看应力分布例2卡紧散热片的不锈钢扣件受力分析扣紧件是一个不锈钢的卡子，因为散热片同功率部件之间的接触力同最终的散热有很大关系，因此研究力的大小是很有意义的。
有限元实例分析
作业一：有限元分析实例实例：请对一个盘轴配合机构进行接触分析。
轴为一等直径空心轴，盘为等厚度圆盘，其结构及尺寸如图所示。
盘和轴为一种材料，材料参数为：弹性模量Ex=2.5E5，泊松比NUXY=0.3，摩擦系数MU=0.25，试采用有限元计算方法分析轴和盘在过盈配合时的应力应变分布以及将轴从盘心拔出时轴和盘的接触情况。
问题分析说明（1）本题主要分析装配过程中结构的静态响应，所以分析步选择通用静态分析步。
由于为过盈配合，属于大变形，故应考虑几何非线性的影响。
（2）模型具有轴对称性，所以可以采取轴对称模型来进行分析，先建立二维模型计算，再转换为三维模型计算，这样可以节省计算时间。
分析过程由两个载荷步组成, 第一个载荷步为过盈分析, 求解过盈安装时的情况。
第二个载荷步为将轴从盘心拔出时的接触分析, 分析在这个过程中盘心面和轴的外表面之间的接触应力。
它们都属于大变形问题, 属于非线性问题。
在分析时需要定义一些非线性选项来帮助问题的收敛。
（3）接触面之间有很大的相对滑动，所以模型要使用有限滑移。
模型建立的分析说明（1）进定义单元类型此项实例分析的问题中涉及到大变形, 故选用So li d185 单元类型来建立本实例入部件模块，的模型。
盘轴接触问题属于面面接触, 目标面和接触面都是柔性的,将使用接触单元T ARGET 170 和CO NTAT17 4来模拟接触面。
分别创建名为为part1、part2的部件。
（2）定义材料属性，在线性各向同性材料属性对话框中的EX (弹性模量) 文本框中输入 2 . 5E5,PRX Y (泊松比) 文本框中输入0 . 3,并将定义的材料属性赋予给part1和part2。
如下图所示。
（3）进入装配模块，创建两者间的装配关系。
（4）进入分析步模块定义名为step1和step2的两个分析步。
（5）进入相互作用模块，创建相互作用属性，设置摩擦系数；然后定义接触关系。
如下图所示。
（6）进入载荷模块，创建边界条件，依次定义名为BC -2（类型为：完全固定）、BC -3（类型为：位移/转角，约束U1、UR3），分析步均为Initial 。
ansysworkbench 接触实例分析
前言WokBench 是众所周知的好东西，以下是自己琢磨的一个小应用，肯定有不对的地方，欢迎指出，便于大家共同提高。
问题描述这是一个塑料小卡扣的例子，主要想使用WorkBench 了解在使用中，塑料件的变形是否足够。
模型是用ProE 制作的，为了简化，只切取了关于变形的部分，如下图：其中蓝色的部分是活动的，只有一个方向的运动，红色的部分是固定的。
大体的尺寸如下，单位是毫米：注意：在模型中，蓝色和红色部件的距离要控制好（这是由ProE 中，模型装配关系决定的），如果太近，软件将自动计算出一个接触区域，但对于这个例子，还需要手动扩大接触区域。
如果距离太远，在手动设置Pinball 类型的接触区域时，Pinball 的半径要设得很大，可能导致无法计算。
请参考上面的尺寸图纸调节两个部件之间的距离。
之后，设置接触面（2、3）：需要将两个部件在运动过程中，会接触的地方一一标出，千万不要加无用的面。
将Pinball Region 设置为Radius 方式（4），并将Radius 设置一个合适的值（5），本例设置了3 毫米（如图，会形成一个蓝色的大圆球），求解的时候软件会使用这个PinBall 自动探测接触。
还需要将接触方式设置为无摩擦的（6）。
最后将接触面计算方式设置为Adjust To Touch（7）。
也可以尝试其他的方式，不过对于这个仅研究红色部件变形的例子就无所谓了。
关于单元格WorkBench 中可以不自行划分单元格（在解算的时候，如果没有手动的设置，软件就会先自动划分），软件帮你自动产生。
如果你的其他设置正确，即便是这个自动的值也能很精确了。
添加分析这个分析用静力学就可以了（1）。
之后要设置Analysis Setting（2）。
将Nuber Of Step 设置为2（3）。
注意：1）蓝色部件在运动的过程中，先压迫红色部件，再逐渐松开，因此必须将这个过程至少分解为至少两个阶段（阶段指“Step”）。
ANSYS接触实例分析参考
ANSYS接触实例分析参考1．实例描述一个钢销插在一个钢块中的光滑销孔中。
已知钢销的半径是0.5 units, 长是2.5 units，而钢块的宽是 4 Units, 长4 Units,高为1 Units,方块中的销孔半径为0.49 units,是一个通孔。
钢块与钢销的弹性模量均为36e6，泊松比为0.3.由于钢销的直径比销孔的直径要大，所以它们之间是过盈配合。
现在要对该问题进行两个载荷步的仿真。
（1）要得到过盈配合的应力。
（2）要求当把钢销从方块中拔出时，应力，接触压力及约束力。
2．问题分析由于该问题关于两个坐标面对称，因此只需要取出四分之一进行分析即可。
进行该分析，需要两个载荷步：第一个载荷步，过盈配合。
求解没有附加位移约束的问题，钢销由于它的几何尺寸被销孔所约束，由于有过盈配合，因而产生了应力。
第二个载荷步，拔出分析。
往外拉动钢销1.7 units，对于耦合节点上使用位移条件。
打开自动时间步长以保证求解收敛。
在后处理中每10个载荷子步读一个结果。
本篇先谈第一个载荷步的计算。
下篇再谈第二个载荷步的计算。
3．读入几何体首先打开ANSYS APDL然后读入已经做好的几何体。
从【工具菜单】-->【File】-->【Read Input From】打开导入文件对话框找到ANSYS自带的文件（每个ansys都自带的）\Program Files\Ansys Inc\V145\ANSYS\data\models\block.inp【OK】后，四分之一几何模型被导入。
4．定义单元类型只定义实体单元的类型SOLID185。
至于接触单元，将在下面使用接触向导来定义。
5．定义材料属性只有线弹性材料属性：弹性模量36E6和泊松比0.36．划分网格打开MESH TOOL,先设定关键地方的网格划分份数然后在MESH TOOL中设定对两个体均进行扫略划分，在volumeSweeping中选择pick all，按下【Sweep】按钮，在主窗口中选择两个体，进行网格划分。
ansys实例命令流-简支梁接触分析命令流
/FILNAME,Simply Supported Beam,1/TITLE, Beam Analysis/PREP7 !前处理模块。
/VIEW,1,1,1,1 !调整窗口视角。
ET,1,PLANE42 !定义单元。
KEYOPT,1,3,3!定义材料属性。
MPTEMP,,,,,,,, !定义触头材料属性。
MPTEMP,1,0MPDATA,EX,1,,2.06e5MPDATA,PRXY,1,,0.35MPTEMP,,,,,,,, !定义接触单元属性(摩擦系数)。
MPTEMP,1,0MPDATA,MU,2,,0.16R,1,30, !梁板度。
! 建立模型。
BLC4, , ,960,80 !建立梁。
CYL4,40,40,20CYL4,920,40,20FLST,3,2,5,ORDE,2FITEM,3,2FITEM,3,-3ASBA, 1,P51X !建立梁完成。
CYL4,40,40,20 !建立轴。
CYL4,920,40,20 !建立轴完成。
! 网格划分。
FLST,5,20,4,ORDE,2 ! 定义网格密度。
FITEM,5,1FITEM,5,-20CM,_Y,LINELSEL, , , ,P51XCM,_Y1,LINECMSEL,,_YLESIZE,_Y1,4, , , , , , ,1 ! 定义网格密度完成。
FLST,5,3,5,ORDE,3 ! 网格划分。
FITEM,5,1FITEM,5,-2FITEM,5,4CM,_Y,AREAASEL, , , ,P51XCM,_Y1,AREACHKMSH,'AREA'CMSEL,S,_YCMDELE,_YCMDELE,_Y1CMDELE,_Y2 ! 网格划分完成。
! 建立接触。
FLST,5,8,4,ORDE,2 ! 定义接触面。
FITEM,5,5FITEM,5,-12LSEL,R, , ,P51XALLSEL,BELOW,LINELPLOTNSLL,S,1CM,contact,NODE ! 定义接触面完成。
hyperworks接触分析实例教程1
hyperworks接触分析实例教程1Pre-Processing for Pipes Impact using RADIOSS Block - RD-3520For this tutorial it is recommended that you complete the introductory tutorial, Getting Started with HyperMesh - HM-1000. Working knowledge of the creation and editing of collectors and card images are a definite pre-requisite. Familiarity with the interfaces panel, and the creation of boundary conditions are useful, although not required.ObjectiveIn this tutorial you will learn how to set up a RADIOSS input file in HyperMesh for analyzing the impact response between two pipes. The modeling steps that are covered are:Model DescriptionThe units used in this tutorial are milliseconds, millimeters and kilograms (ms, mm, kg), and the tutorial is based on RADIOSS Block 51.Pipe model ExerciseStep 1: Load the RADIOSS Block User Profile and import the modelIn the following steps, load the RADIOSS Block User Profile, then import a RADIOSS input deck containing the mesh for the two pipes.Selecting the RADIOSS User Profile sets the FE input reader to RADIOSS Block and loads theRADIOSS51 FE output template. It also loads a RADIOSS macro menu with numerous tools specific to this interface. The graphical user interface is tailored to RADIOSS users with panels such as the admas panel on the 1D page, panel names and options renamed or removed to match RADIOSS terminology.The model is loaded and is organized into two component collectors: 1 and 2 (named by HyperMesh). No materials or properties are defined at this time.Creating materials, sections, and parts for the model.?Defining the contact between the two pipes using /INTER/TYPE7.?Applying a translational initial velocity to a pipe using the /INIVEL card.?Applying local constraints to the other pipe using the /BCS card.1.Load the Radioss Block User Profile using the pull down menu Preferences > User Profile > RADIOSS(Block Format) > BLOCK90 ....2.From the files panel, select the import subpanel and import the following RADIOSS deck:/tutorials/hwsolvers/radioss/pipesd00.Notes:On import of a RADIOSS deck, any HyperMesh warning and error messages are written to a fileOn import, any RADIOSS cards not supported by HyperMesh are written to the control cardunsupp_cards . This card is accessed from the control cards panel on the BCs page and is apop-up text editor. The unsupported cards are exported with the rest of the model.Care should be taken if an unsupported card points to an entity in HyperMesh. An example of this is an unsupported material referenced by a /PART card. HyperMesh stores unsupported cards as text and does not consider pointers.On import, HyperMesh renumbers entities having the same ID as other entities. In HyperMesh, for example, all elements must have a unique ID. The message file radiossblk.msg provides a list of renumbered elements and their original and new IDs.Step 2: Understand the relationships between the /PART, /SHELL, /MAT and /PROP cards in HyperMeshA /PART shares attributes such as section properties (/PROP) and a material model (/MAT). A group of shells (/SHELL) sharing common attributes generally share a common part ID (PID).The figure below shows how these keywords are mapped to HyperMesh entities:Map to HyperMesh entities Component, property and material collectors are created and edited from the collectors panel.For the RADIOSS keyword interface, there is only one component card image and it is named Part . There are several property card images, such as P1_shell P2_truss , P14_solid . There are many material card images, such as M1_ELAST , M48_HONEYCOMB .The complete list of card images is available from the collectors panel, as you assign card images to the various types of collectors.A HyperMesh card image allows you to view the image of keywords and data lines for defined RADIOSS entities as interpreted by the loaded template. The keywords and data lines appear in the exported RADIOSS input file as you see them in the card images. Additionally, for some card images, you can define and edit various parameters and data items for the corresponding RADIOSS Block.Use the card (card editor) panel from the permanent menu to review and edit card images. Also, for many entities, their card image can be viewed and edited from the panels in which they are created.Step 3: Create a /MAT card In HyperMesh, a /MAT card is associated to a material collector. To relate it to a /PART card, the material collector needs to be assigned to a component collector.You can assign the material to the component collector as you create the component using the createsubpanel of the collectors panel. In situations where the material was not assigned to the component at the time of creation (and in this case, a dummy material is created with the same name as the componentnamed radiossblk.msg . This file is created in the folder from which HyperMesh is started. The content of the file is also displayed in a pop-up window./SHELL elem_ID part_ID> Organized into component collectors/PART part_ID prop_ID mat_ID > Component collector with a componentcard image/PROP prop_ID > Property collector with a property cardimage/MAT mat_ID>Material collector with a material cardimagecollector), update the component collector's definition by assigning the material in the update subpanel of the collectors panel.In this step, create a material collector with the M1_ELAST card image using the collectors panel. This material will be assigned to both pipes.1.Create/edit a material collector with the name elast1, and a M1_ELAST card image using the createsubpanel of the collectors panel.2.In the card previewer, click [Rho_I] to activate its field.For density, specify 7.8 E-6For Young's modulus [E], specify 208For Poisson's ratio [nu], specify 0.30Note:If you have difficulties completing any task with the creation, update or editing of collectors in this tutorial, refer to the on-line help for the collectors panel by clicking help from the permanent menu. Hint:Any collector that you mistakenly create instead of create/edit can be edited using the card image subpanel of the collectors panel.In this step, we created the material we will use for the analysis. We can now define the /PROP card that will be used to define the properties of the elements in the model.Step 4: Create a /PROP cardIn HyperMesh, the /PROP card is assigned to a property collector. To generate this card, create a property collector using either Property collector icon in the toolbar or from the pull down menu Properties > Create. The model consists of two pipes modeled with shell elements. Create a property collector witha /PROP/SHELL card that will be used for all the elements.1.Create a property collector with the name prop shell, set Type = to surface, set card image= toP1_SHELL, and Thickness = 2.5 card.2.In the card previewer for the /PROP/SHELL card, check 2.5for the shell thickness [Thick].Step 5: Assign the /PART, /MAT and /PROP cards to the elementsAssign the /PART card to the component for the coarse pipe and specify the /PROP/SHELL card ID in it.1.Load the Component Collector Update panel from either the toolbar Component collector or from thepull down menu, select Component > Create.2.Click update subpanel.3.Select the pipe 1 component.4.For card image, select PART.5.Assign the material elast1 and property prop shell.6.Click update.Repeat this procedure for pipe 2.Step 6: Create contact cardsRADIOSS contacts are created in the interfaces panel from the Analysis page or from the pull down menu, selectBCs>Create>Interfaces.A RADIOSS contact is a HyperMesh group. When you want to manipulate an /INTER card, such as delete it, renumber it, or turn it off, you need to work with HyperMesh group entities.In this step, create a contact between the two pipes using /INTER/TYPE7. The pipe with the coarser mesh (2) will be the master surface while the one with finer mesh (1) will be the slave surface. RADIOSS has multiple ways to define master andslave entity types from which to choose; in this example define the master and slave entities as components, doing this the master will be exported as /SURF/PART and the slave asa /GRNOD/PART.First create a group interface with the TYPE7card image using the create subpanel of the interfaces panel.Next, add the master and slave to the group using the add subpanel. Lastly, edit the group’s card image using the card image subpanel, and specify a friction coefficient.Automatically TYPE7 is selected for card image.In this step, we defined the contact between the two pipes as /INTER/TYPE7.Creating boundary conditions for a deck in RADIOSS Block can be efficiently carried out with the BC’s Manager available on the Utility Browser , which can be switched on the screen from the View pull down menu.Every load needs to be stored in a load collector with the corresponding card image. For example, store the velocity loads ina velocity load collector and boundary conditions in the respective collector.1.Create a group with the name contact and the TYPE7 card image using the create subpanel of the interfaces panel.2.From the BCs page, select the interfaces panel or from the pull down menu, select BCs > Create > Interfaces .3.Select the create subpanel.4.In the name= field, enter contact .5.For type=, select TYPE7.6.Optionally select a color.7.Click create .8.From the interfaces panel, select the add subpanel.9.Under master , set the entity type to comps :10.Click the yellow comps selector and select the coarser pipe component 2.11.Click update in the master: field, to the right of the yellow comps selector.12.For slave:, select comps .13.Click the yellow comps selector in the slave:line and select the finer mesh pipe component 1.14.In the slave: field, click update .15.Click review to graphically view the entities in the interface, the master entities of the interface are drawn in blue and the slave entities in red.16.Edit the definition of the group, using the card image subpanel to set the static coefficient to 0.10.17.From the interfaces panel, select the card image subpanel.18.Click edit to review the card image.19.Specify 0.10 for the static coefficient [FRIC].Step 7: Create an initial velocity (/INIVEL/TRA)In this step, we will apply a translational initial velocity to the coarse pipe using /INIVEL/TRA applied to a predefined set of nodes /GRNOD/PART.1.Click on BC’s Manager on the Utility browser.2.From the Boundary Conditions Manager, enter the name tran_vel and select the type as initialvelocity under the Create header.3.For creating a entity set of type GRNOD which is referred to in the INIVEL card, click on Parts button andselect the coarser pipe from the GUI (component ID 2). Click proceed.4.In the BC’s Manager, enter the initial velocity components as 0,0 and -30 for Vx, Vy and Vz fields.5.There is an option for creating/referring the initial velocity card to a local coordinate system. However ifnothing is specified, the global coordinate system is selected by default.6.There is an option to specify the size of the load display on the screen under Label Scale.7.Click the create button. Double check in the model browser for your reference that a load collector andan entity set is created.This completes the creation of an initial velocity for the pipe in the negative global Z direction.Step 8: Create a /BCS and constrain the finer mesh pipeIn this step, we will fully constrain the end nodes of the bottom pipe by using the Boundary Conditions Manager.1.In the BCs Manager under the Create subheading, enter the Name SPC and set Select type asBoundary Condition.2.Now specify the node set of type as GRNOD for the BCS card, switch the entity from parts to Nodesand select the end nodes of the bottom pipe, which are to be constrained.3.Under the Boundary condition components subheading (as illustrated below) activate all thetranslational and rotational check boxes. Click the create button. A load collector with a BCS card is Step 9: Create output definitions and control cardsA window pops up.created and applied the nodes as selected in the above steps. A corresponding node set is created.1.In the Utility menu of the solver browser, click the Engine button.2.Select the options, as shown below.Step 10: Export the modelThis concludes this tutorial. You may discard this HyperMesh model or save it for your own reference. In this tutorial we introduced some of the concepts that govern the HyperMesh interface to RADIOSS. We also use numerous panels that allowed us to do basic modeling in terms of RADIOSS such as defining contacts or boundary conditions.EXERCISE EXPECTED RESULTS3.Click Apply and Close to store the selected options for control and output.1.On Standard toolbar at the top of the HyperMesh window, click Export.2.For File:, click the foldericon and navigate to destination directory where you want to run.3.For Name , enter pipeimpact and click Save .4.Click the downward-pointing arrows next to Export options to expand the panel.5.Click Auto export engine file to export the engine file with the model file.6.Click Export to export both model and engine file.Final deformation and energy balance plotSee alsoSee HyperMesh Tutorials for a complete list of tutorials.。
ANSYS接触实例分析参考
ANSYS接触实例分析参考ANSYS是工程仿真领域广泛使用的一种有限元分析软件。
在实际工程中，接触问题经常出现，例如机械装配中的接触、摩擦、磨损等现象需要进行分析和优化。
本文将介绍几个ANSYS接触实例，并分析其分析方法和结果。
第一个实例是机械装配中的接触分析。
假设有一个由两个金属块组成的简单装配，要分析它们之间的接触情况。
首先需要建立两个金属块的几何模型，并进行网格划分。
然后，使用ANSYS中的接触分析模块，设置接触类型、接触参数和材料特性等。
接着，施加相应的边界条件和载荷条件，运行分析并获取接触压力和接触面积等结果。
最后，根据结果对接触情况进行评估和优化。
第二个实例是摩擦接触问题的分析。
假设有一个由摩擦带和基体组成的摩擦副，需要分析摩擦力和热量的分布。
首先需要建立摩擦带和基体的几何模型，并进行网格划分。
然后，使用ANSYS中的摩擦接触分析模块，设置摩擦带和基体的材料特性、摩擦系数和接触压力等参数。
接着，施加相应的边界条件和载荷条件，运行分析并获取摩擦力、摩擦热量和温度分布等结果。
最后，根据结果对摩擦副的性能进行评估和优化。
第三个实例是磨损接触问题的分析。
假设有一个由金属零件和砂轮组成的磨削装置，需要分析金属零件表面的磨损情况。
首先需要建立金属零件和砂轮的几何模型，并进行网格划分。
然后，使用ANSYS中的磨损接触分析模块，设置金属零件和砂轮的材料特性、初始接触压力和磨粒等参数。
接着，施加相应的边界条件和载荷条件，运行分析并获取磨损量、磨损深度和磨损形貌等结果。
最后，根据结果对磨削装置进行评估和优化。
以上三个实例只是ANSYS接触分析的一小部分应用，接触分析的对象和问题种类都非常多样。
在实际工程中，可以根据具体问题的特点选择不同的接触分析方法和技术，以获取更准确和可靠的结果。
同时，还可以通过对接触问题的分析和优化，改善产品的性能和可靠性，提高工程效率和经济效益。
总结起来，ANSYS接触实例分析主要包括机械装配中的接触分析、摩擦接触问题的分析和磨损接触问题的分析。
ANSYS高级接触分析
§2 面－面接触单元
• §1 概述
• 面－面接触单元，是模拟任意两个表面间接触的方法。表面可以具有任意形 状。是 ANSYS 中最通用的接触单元。精度高、特性丰富还可使用接触向导， 建模方便。（其它接触单元目前尚不能用向导）。
•
面－面接触单元在面的高斯点处传递压力，这种先进技术使面－面接触
单元具有很多优点：
• 库仑法则是宏观模型，表述物体间的等效剪力 FT 不能超过正压力 FN 的一部分：
FT <= μ× FN 式中： μ－ 摩擦系数 • 一旦所受剪力超过 FT，两物体将发生相对滑动。 4、弹性库仑摩擦模型：允许粘着和滑动。
§3 自动时间步、控制
接触单元的 Keyopt（7）选项控制时间步的预报。 • 0－无控制：不影响时间步尺寸。当自动时间步开关
和土壤的接触
§2 接触单元
§3 关于耦合和约束方程的应用
• 如果接触模型没有摩擦，接触区域始终粘在一起，并且分 析是小挠度、小转动问题，那么可以用耦合或约束方程代 替接触。
• 使用耦合或约束方程的优点是分析还是线性的
接触问题的一般特性
• §1 接触刚度
• 1、所有的 ANSYS 接触单元都采用罚刚度（接触刚度）来 保证接触界面的协调性
• 接触单元就是覆盖在分析模型接触面上的一层 单元。
• 在 ANSYS 中可以采用三种不同的单元来模拟 接触：
◦
面一面接触单元；
◦
点一面接触单元；
◦
点一点接触单元。
§2 接触单元
• 不同的单元类型具有完全不同的单元特性和分 析过程。
• 1. 面一面接触单元用于任意形状的两个表面接 触
• 不必事先知道接触的准确位置； • 两个面可以具有不同的网格； • 支持大的相对滑动； • 支持大应变和大转动。 • 例如: 面一面接触可以模拟金属成型，如轧制
弹塑性力学基础与有限元分析-接触分析实例
06
结论与展望
结论
1
本文通过理论分析和有限元模拟，深入研究了弹 塑性力学基础与有限元分析在接触分析中的应用。
2
研究结果表明，弹塑性力学基础与有限元分析在 接触分析中具有较高的精度和可靠性，能够有效 地模拟复杂接触问题。
3
本文所采用的有限元分析方法在处理接触问题时 具有较好的通用性和扩展性，为进一步研究复杂 接触问题提供了有力支持。
弹塑性本构模型
弹塑性本构模型的定义
弹塑性本构模型是描述弹塑性材料力学行为的数学模型，它通过应力应变关系来描述材料的弹塑性行 为。
常见的弹塑性本构模型
常见的弹塑性本构模型包括Mohr-Coulomb模型、Drucker-Prager模型、Cam-Clay模型等。这些模 型在描述材料的弹塑性行为方面各有特点，适用于不同的材料和工程问题。
接触面完全贴合，无相对运动。
滑动状态
接触面部分贴合，存在相对运动。
混合状态
接触面同时存在分离、粘结和滑动。
接触检测与跟踪
初始接触检测
确定初始状态下接触面的位置和状态。
接触状态跟踪
实时监测接触面的运动状态和相互作用。
接触面更新
根据接触状态调整接触面的几何形状和参数。
接触刚度与阻尼
1 2
接触刚度
描述接触面间的相互作用力与相对位移的关系。
求解阶段主要进行有限元 方程的求解，得到各节点 的位移和应力等结果。
ABCD
前处理阶段主要完成有限元 模型的建立和网格划分，为 求解阶段提供输入数据。
后处理阶段主要对求解结果进 行可视化、分析和评估，为工 程设计和优化提供依据。
04
接触分析原理
接触状态描述
分离状态
六起氰化物中毒实例及案例分析
六起氰化物中毒实例及案例分析Revised by Hanlin on 10 January 2021六起氰化物中毒实例及案例分析化学结构中含有氰根（CN-）的化合物均属于氰化物。
一般将其无机化合物归为氰类，有机化合物归为腈类。
氢氰酸，HCN，别名氰化氢，是一种具有苦杏仁特殊气味的无色液体。
易溶于水、酒精和乙醚。
易在空气中均匀弥散，在空气中可燃烧。
氰化氢在空气中的含量达到5.6~12.8%时，具有爆炸性。
氰化氢为气体，其水溶液称氢氰酸。
氢氰酸属于剧毒类。
其主要应用于电镀业（镀铜、镀金、镀银）、采矿业（提取金银）、船舱、仓库的烟熏灭鼠，制造各种树脂单体如丙烯酸树酯、甲基丙烯酸树酯等行业，此外也可在制备氰化物的生产过程中接触到本物质。
氰化物的毒性主要由其在体内释放的氰根而引起。
氰根离子在体内能很快与细胞色素氧化酶中的三价铁离子结合，抑制该酶活性，使组织不能利用氧。
氰化物对人体的危害分为急性中毒和慢性影响两方面。
氰化物所致的急性中毒分为轻、中、重三级。
轻度中毒表现为眼及上呼吸道刺激症状，有苦杏仁味，口唇及咽部麻木，继而可出现恶心、呕吐、震颤等；中度中毒表现为叹息样呼吸，皮肤、粘膜常呈鲜红色，其他症状加重；重度中毒表现为意识丧失，出现强直性和阵发性抽搐，直至角弓反张，血压下降，尿、便失禁，常伴发脑水肿和呼吸衰竭。
氢氰酸对人体的慢性影响表现为神经衰弱综合症，如头晕、头痛、乏力、胸部压迫感、肌肉疼痛、腹痛等，并可有眼和上呼吸道刺激症状。
皮肤长期接触后，可引起皮疹，表现为斑疹、丘疹，极痒。
在某些镀铜、镀镍工艺中可采用无氰电镀。
凡发生氰化氢的工序，应严加密闭，或放在隔离室内局部排风，室内保持负压，防止有毒气体逸出。
在进入用氢氰酸烟熏过的仓库时，必须事先通风，并戴隔离式防毒面具方可入内。
实例11983年12月22日下午7时，上海某化工厂电焊工李某（男、21岁）在该厂丙酮氰醇车间对堵塞的管道进行切割时，不慎管内余存的氢氰酸逸出，李某由此而吸入氰化氢气体，致头晕、乏力，进而呼吸困难、意识丧失，皮肤粘膜呈樱桃红色。
5.3.9 接触设置综合实例[共10页]
![5.3.9 接触设置综合实例[共10页]](https://img.taocdn.com/s3/m/cabc9f6769dc5022abea0018.png)
5.3 状态非线性分析——接触5.3.9 接触设置综合实例通过前面例子的学习，已经了解了WB中接触设置。
下面以一个2D压片弯曲挤压胶片，胶片再承受密封流体压力的例子综合描述接触分析。
本例包含刚柔接触、自接触、密封流体压力。
1．建立2D模型如图5-3-99所示，建立一个含压模板、压片、胶片的2D模型。
由于压片上端为曲线，且压片与胶片均处于相对自由状态，所以很难精确定义压模板和胶片与压片相切的位置，因此压模板距压片有微小间隙，胶片与压片呈过盈状态。
压模板在整个过程中几乎不变形，而且也不是本分析所关注的目标，所以将其定义为刚体；压片在整个过程中存在大的弯曲变形，其结果将表现为首尾相接触，将其材料定义为非线性铝合金；胶片为橡胶件，整个过程中存在大应变，且胶片内部存在自接触可能，将其本构定义为Ogden 3rd Order类型。
压模板，命名tie，刚体压片，命名Surface Body，材料本构为非线性铝合金胶片，命名rub，材料本构为Ogden 3rd Order图5-3-99 2D模型2．2D模型及材料设置调用WB默认材料库内的非线性铝合金（General Non-linear Materials→Aluminum Alloy NL），新增一个材料，命名为rub，本构选择Hyperelastic→Ogden 3rd Order，9个参数分别为：MU1=0.043438MPa，A1=1.3，MU2=8.274E−5MPa，A2=5，MU3=−0.0006895MPa，A3=−2，D1=0.029MPa^−1，D2=0MPa^−1，D3=0MPa^−1。
在Geometry→2D Behavior处定义为Plane Stress（平面应力），如图5-3-100所示。
– 435 –第5章 非线性静力学分析– 436 – 3．Virtual Topology （虚拟拓扑）设置虚拟拓扑一般用于合并几个不同平面，使其保证为一个有限元拓扑模型，除此之外，还可用于分割模型。
《接触分析多点接触》课件
齿轮的接触分析主要关注齿面接触应力和齿根弯曲应力。通过分析这些应力，可 以评估齿轮的承载能力和疲劳寿命，为齿轮的设计和优化提供依据。同时，接触 分析还可以用于研究齿轮的啮合特性，优化齿轮的几何参数和材料选择。
实例三：曲轴的接触分析
总结词
曲轴是发动机中的重要元件，其接触分析对于了解曲轴的工作状态和疲劳寿命至关重要。
识别和判断物体间的接触状态，包括接触 点的位置和方向。
研究接触点上的应力、应变和能量分布， 以及摩擦和磨损行为。
02
多点接触模型
多点接触模型的定义
定义
多点接触模型是一种描述两个或多个 物体在接触过程中相互作用的方法， 它考虑了多个接触点的影响，以更准 确地模拟物体的接触行为。
对比
与传统的单点接触模型相比，多点接 触模型能够更好地处理复杂接触情况 ，更准确地模拟物体的接触响应。
收敛判定
设定合适的收敛准则，确 保计算结果的准确性。
结果分析
后处理
对计算结果进行后处理，如绘制 应力、应变分布图等。
结果评估
根据实际需求，对计算结果进行评 估，判断其是否满足工程要求。
优化建议
基于计算结果，提出优化建议，以 改进接触性能或提高系统稳定性。
04
接触分析的实例
实例一：轴承的接触分析
适用性强
可用于各种复杂接触情况，如曲面接 触、多体接触等。
接触分析的优缺点
• 可视化效果好：能够生成高质量的图形和动画，方便 用户理解和分析接触过程。
接触分析的优缺点
计算量大
由于多点接触分析需要考虑多个 接触点之间的相互作用，计算量 较大，需要高性能计算机支持。
建模复杂
需要建立详细的几何模型和接触 约束，建模过程较为复杂。
接触分析全文 文档
参考：接触分析全文一般的接触分类 1ANSYS接触能力 1点─点接触单元 1点─面接触单元 2面─面的接触单元 2执行接触分析 3面─面的接触分析 3接触分析的步骤： 3步骤1：建立模型，并划分网格 4步骤二：识别接触对 4步骤三：定义刚性目标面 4步骤4：定义柔性体的接触面8步骤5：设置实常数和单元关键字10步骤六：21步骤7：给变形体单元加必要的边界条件22步骤8：定义求解和载步选项22第十步：检查结果23点─面接触分析25点─面接触分析的步骤26点－点的接触36接触分析实例（GUI方法）39非线性静态实例分析（命令流方式）42接触分析接触问题是一种高度非线性行为，需要较大的计算资源，为了进行实为有效的计算，理解问题的特性和建立合理的模型是很重要的。
接触问题存在两个较大的难点：其一，在你求解问题之前，你不知道接触区域，表面之间是接触或分开是未知的，突然变化的，这随载荷、材料、边界条件和其它因素而定；其二，大多的接触问题需要计算摩擦，有几种摩擦和模型供你挑选，它们都是非线性的，摩擦使问题的收敛性变得困难。
一般的接触分类接触问题分为两种基本类型：刚体─柔体的接触，半柔体─柔体的接触，在刚体─柔体的接触问题中，接触面的一个或多个被当作刚体，（与它接触的变形体相比，有大得多的刚度），一般情况下，一种软材料和一种硬材料接触时，问题可以被假定为刚体─柔体的接触，许多金属成形问题归为此类接触，另一类，柔体─柔体的接触，是一种更普遍的类型，在这种情况下，两个接触体都是变形体（有近似的刚度）。
ANSYS接触能力ANSYS支持三种接触方式：点─点，点─面，平面─面，每种接触方式使用的接触单元适用于某类问题。
为了给接触问题建模，首先必须认识到模型中的哪些部分可能会相互接触，如果相互作用的其中之一是一点，模型的对立应组元是一个结点。
如果相互作用的其中之一是一个面，模型的对应组元是单元，例如梁单元，壳单元或实体单元，有限元模型通过指定的接触单元来识别可能的接触匹对，接触单元是覆盖在分析模型接触面之上的一层单元，至于ANSTS使用的接触单元和使用它们的过程，下面分类详述。
本文介绍了ANSYS APDL命令和轴承接触分析相关实例
本文介绍了ANSYS APDL命令和轴承接触分析相关实例。
finish/clear/filename,ex5-3 !设定工作名/title,the contact analysis of Bearing !设定工作标题!************************************!进入前处理模块!************************************/prep7*set,roll_r,5.7 !定义变量参数*set,inside_r1,15*set,inside_r2,18.5*set,outside_r1,29.8*set,outside_r2,33.3*set,deep,8*set,axes_r3,15.5*set,fillet_r4,1*set,axes_deep,20*set,fillet_r5,1!定义单元类型和属性et,1,solid45 !定义单元类型solid45MP,EX,1,2.06e5 !定义弹性模量MP,PRXY,1,0.3 !定义泊松比MP,MU,1,0.3 !定义摩擦系数!建立模型sph4,,-inside_r2-roll_r,roll_r+0.9 !生成轴承滚珠cyl4,,,inside_r1,-120,inside_r2,-60,deep !生成+z向的轴承内圈cyl4,,,outside_r1,-120,outside_r2,-60,deep !生成+z向的轴承外圈cyl4,,,inside_r1,-120,inside_r2,-60,-deep !生成-z向的轴承内圈cyl4,,,outside_r1,-120,outside_r2,-60,-deep !生成-z向的轴承外圈wpoff,0,0,-deep !偏移工作平面cyl4,,,0,-120,axes_r3,-60,-axes_deep !生成装配轴lfillt,57,59,fillet_r4,, !生成倒角线al,57,62,59,54,61 !由线生成面vrotat,32,,,,,,42,37,60,1, !旋转面生成体vdele,6,,,1 !删除编号为6的体!lfillt,31,39,fillet_r5,, !生成倒角线al,34,39,53,31,37 !由线生成面vdrag,27,,,,,,35 !拖拉面生成体vdele,4,,,1 !删除编号为4的体!vgen,,7,,,,,0.27,,,1 !平移编号为7的体!wpoff,0,0,deep !偏移工作平面csys,1 !激活柱坐标系asel,s,loc,x,inside_r2 !选择x=inside_r2的面asel,a,loc,x,outside_r1 !选择x=ouside_r1的面vsba,1,all !体被面分割vdele,4,,,1 !删除编号为4的体vdele,8,,,1 !删除编号为8的体allsel,all !选择全部图元vsel,u,volu,,7 !不选编号为7的体vglue,all !粘接全部的体!以下通过一些布尔操作以方便网格划分wpoff,0,-inside_r2-roll_r,0 !偏移工作平面vsbw,1 !用工作平面分割体1wpro,,-90, !旋转工作平面vsbw,2 !用工作平面分割体2vsbw,3 !用工作平面分割体3wpro,,,-90 !旋转工作平面vsbw,1 !用工作平面分割体1vsbw,2 !用工作平面分割体2vsbw,5 !用工作平面分割体5vsbw,6 !用工作平面分割体6!voffst,2,-4 !沿面的法向平移面2生成体voffst,9,-4 !沿面的法向平移面9生成体voffst,23,-4 !沿面的法向平移面23生成体voffst,53,-4 !沿面的法向平移面53生成体!voffst,3,4 !沿面的法向平移面3生成体voffst,25,4 !沿面的法向平移面25生成体voffst,38,4 !沿面的法向平移面38生成体voffst,58,4 !沿面的法向平移面58生成体!vovlap,all !对体进行搭接操作vdele,25,,,1 !删除编号为25的体及其所属的低阶图元vdele,32,,,1 !删除编号为32的体及其所属的低阶图元vdele,33,,,1 !删除编号为33的体及其所属的低阶图元vdele,34,,,1 !删除编号为34的体及其所属的低阶图元!vdele,31,,,1 !删除编号为31的体及其所属的低阶图元vdele,35,,,1 !删除编号为35的体及其所属的低阶图元vdele,36,,,1 !删除编号为36的体及其所属的低阶图元vdele,37,,,1 !删除编号为37的体及其所属的低阶图元vglue,all !对体进行粘接操作!划分网格esize,2 !设定网格单元尺寸mshape,0,3d !设定网格形状为六面体单元mshkey,1 !设定为映射网格划分方式vsel,s,volu,,1,3,2 !选择编号为1、3 的体vsel,a,volu,,4,5 !同时选择编号为4,5的体vsel,a,volu,,9 !同时选择编号为9的体vsel,a,volu,,12,14 !同时选择编号为12、13、14的体cm,sphere,volu !生成体的组件spherevmesh,all !对体进行网格划分!esize,1 !设定网格单元尺寸!vsel,inve,volu !对当前体选择集进行反选vsel,s,volu,,6vsel,a,volu,,22,23vsel,a,volu,,26,30vsel,a,volu,,38,40vsweep,all !对体sweep网格划分esize,1.5 !设定网格单元尺寸allsel,allvsweep,8,50,49 !设定源面和目标面并进行sweep网格划分vsweep,7,32,37 !设定源面和目标面并进行sweep网格划分!!生成耦合设置cmsel,s,sphere,volu !选择名称为sphere的组件vgen,2,all,,,,,,,0 !复制该组件cmsel,s,sphere,volu !选择名称为sphere的组件vclear,all !清除该组件包含图元的网格vdele,all,,,1 !删除该组件包含的图元!csys,1 !激活柱坐标系asel,s,loc,x,inside_r2 !选择x=inside_r2的面asel,a,loc,x,outside_r1 !同时选中x=outside_r1的面asel,u,loc,y,-90 !从当前选择集中不选y=-90的面nsla,s,1 !选择面所属的节点nrotat,all !旋转节点坐标系与当前激活坐标系平齐cpintf,ux !在重合节点生成自由度ux的耦合设置cpintf,uy !在重合节点生成自由度uy的耦合设置cpintf,uz !在重合节点生成自由度uz的耦合设置!!设定接触参数/PREP7ALLSEL,ALL !选择全部图元/COM, CONTACT PAIR CREATION - START !接触对设置开始/GSA V,cwz,gsav,,temp !将当前的图形设置保存在cwz.gsav文件中!MP,MU,1,0.3 !定义摩擦系数MAT,1 !激活材料属性1R,3 !定义实常数3REAL,3 !激活实常数3ET,2,170 !定义单元类型2ET,3,174 !定义单元类型3KEYOPT,3,9,0 !设定单元类型3的关键项9KEYOPT,3,10,1 !设定单元类型3的关键向10R,3,,,0.1, !设定法向接触刚度为0.1!生成目标面ASEL,S,,,30 !选择编号为30的面ASEL,A,,,90 !同时选中编号为90的面ASEL,A,,,98 !同时选中编号为98的面ASEL,A,,,104 !同时选中编号为104的面ASEL,A,,,113 !同时选中编号为113的面ASEL,A,,,138 !同时选中编号为138的面ASEL,A,,,143 !同时选中编号为143的面CM,AREA_TARGET,AREA !生成目标面组件target TYPE,2 !激活单元类型2NSLA,S,1 !选择面所属的节点ESLN,S,0 !选择节点依附的单元ESURF !在当前选择的单元上覆盖生成单元ESEL,ALL !选择所有的单元!生成接触面ASEL,S,,,35 !选择编号为35的面ASEL,A,,,36 !同时选中编号为36的面CM,AREA_CONTACT,AREA !生成接触面组件contact TYPE,3 !激活单元类型3NSLA,S,1 !选择面所属的节点ESLN,S,0 !选择节点依附的单元ESURF !在当前选择的单元上覆盖生成单元ALLSEL !选择全部图元ESEL,ALL !选择全部单元ESEL,S,TYPE,,2 !选择单元类型为2的单元ESEL,A,TYPE,,3 !同时选中单元类型为3的单元ESEL,R,REAL,,3 !在当前选择集中选出实常数为3的单元/PSYMB,ESYS,1 !打开单元坐标系显示/PNUM,TYPE,1 !打开单元类型编号/NUM,1 !打开颜色显示EPLOT !图形显示单元ALLSEL,ALL !选择全部图元/GRES,cwz,gsav !从cwz.gsav文件中恢复图形设置/COM, CONTACT PAIR CREATION - END !接触对结束!**********************************!进入求解模块!**********************************/solu !进入求解模块csys,1 !激活柱坐标系nsel,s,loc,x,outside_r2 !选择x=outside_r2的节点d,all,all !在节点上施加全部自由度约束asel,s,loc,y,-60 !选择y=-60的面asel,a,loc,y,-120 !同时选中y=-120的面da,all,symm !施加对称边界条件!施加装配轴的移动位移da,33,uz,2*deep !在编号为33的面上施加位移约束!非线性求解设置lnsrch,on !打开线性搜索pred,on !打开预测矫正autot,on !打开自动时间步nsubst,40,100,10 !设定子步数outres,all,all !输出所有子步上的全部数据allsel,all !选择所有图元solve !开始求解!**********************************!进入后处理模块!**********************************/POST1 !进入通用后处理器PLDISP,2 !图形显示结构变形图/DSCALE,1,1.0 !设定显示比例为1.0/EXPAND,6,POLAR,FULL,0,60,0,, !将结果扩展到360度/REPLOT !重绘当前图形PLNSOL,S,Z !图形显示z方向的应力PLNSOL,S,EQV !图形显示平均等效应力PLNSOL, CONT,PRES, 0,1.0 !接触应力等值线图PLNSOL, CONT,STA T, 0,1.0 !接触状态等值线图PLDI,2,ANMODE,10,0.5, ,0 !变形前后动画效果FINISH !后处理模块结束/EXIT,ALL !退出并保存全部数据。
ANSYS齿轮接触案例讲解
齿轮的接触分析实例
c. ANSYS提示是否删除原来的线，选择【Delete】，点击 【OK】。
（18）偏移工作平面到总坐标系的原点： WorkPlane>Offset WP to>Global Origin。
（19）将工作平面与总体坐标系对齐： WorkPlane>Align WP with>Global Cartesian。
（2）在下图的列表框中选择“Solid”, “4node 182”, 单击【OK】。
齿轮的接触分析实例
(3)在下图的Element Types对话框中单击【Options】 弹出单元选项对话框，对PLANE182单元进行设 置。设置完成后点击【OK】，然后【Close】。
齿轮的接触分析实例
2.3 定义实常数
齿轮的接触分析实例
（4）偏移工作平面到给定位置。 a.从实用菜单中选择
WorkPlane>Offset WP to>Keypoints + 。 b.在ANSYS图形窗口选择110号辅 助点，点击【OK】。 （5）旋转工作平面 a.从实用菜单中选择 WorkPlane>Offset WP by Increments。 b.在“XY,YZ,ZX,ZXAngles”文本框 中输入-50,0,0,点击【OK】.
Preprocessor>Modeling>Operate>Booleans>Add>Lines。 b.分别选择齿底上的两条线，点击【OK】。 c.把齿底上的所有线加起来。 (28)把所有线粘起来。 a.从主菜单选择
Preprocessor>Modeling>Operate>Booleans>Glue>Lines。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实例1：带孔方板的接触分析实例
静力分析分为线性静力分析和非线性静力分析。
线性分析，即外载荷与系统的响应之间为线性关系。
而在真实的物理结构中，结构刚度随变形而发生变化，即所谓的“非线性分析”。
非线性问题的分类：
材料非线性（material nonlinearity）:材料的应力—应变关系为非线性，即本构关系为非线性。
几何非线性：指位移的大小对结构响应的影响，包括大位移、大转动、初始应力、几何刚性化和突然翻转等，实例：过盈装配
边界条件非线性：即边界条件在分析过程中不断的发生变化，实例：接触问题
后处理：
(1)延展平面应力单元来构造三维视图
View—odb display option—sweep&extrude
(2)显示接触压强(CPRESS)和接触状态(COPEN)
Result—field output
接触分析中的主要问题及解决收敛问题的常用方法：
1.接触面的分类：
由单元构成的柔性或刚性接触面；节点构成的接触面；解析刚体接触面
2.一个接触对最多有一个由节点构成的接触面；如果只有一个接触面，称为“自接触”
3.模拟接触问题的两种算法：
通用接触算法和接触对算法。
其中后一种算法，对接触面的类型有较多的限制，需要指定相互接触的面
4.定义接触对：
主面必须连续，且由节点组成；如果是有限滑移，主面发生接触的部位必须光滑，不能有尖角，不要让从面节点落在主面之外（尤其不能落在主面的背面）；接触面在发生接触的部位有较大尖角或凹角时，应该将其分别定义为两个面
5.有限滑移和小滑移
在有限滑移过程中，需不断的判断从面节点和主面的哪部分接触；而在小滑移中，在分析开始时就定义了接触的部位，而在分析过程中不会发生变化
小滑移也可以用于几何非线性问题，有两种算法：点对面和面对面算法
实例2：过盈装配过程模拟
边界条件给出的位移值是相对于模型初始状态的绝对位移值，而不是当前分析步中的增量值；可以控制在每个增量步中的位移均与；轴对称模型将三维模型转化为二维平面问题，只有三个自由度。
后处理：
1.延展轴对称单元来构造等效的三维视图；显示接触压强；显示压头节点的反作用力RF
接触分析中用到的输出变量：
CFS: Total force due to frictional stress (CFS n, n = 1, 2, 3).，摩擦应力的合力
CFN:接触压力的合力
CAREA:接触面积
CMN:接触压力相对于原点的合力矩
CMS:摩擦压力相对于原点的合力矩
CFT:接触压力和摩擦压力的合力
CMT:接触压力和摩擦压力相对于原点的合力矩
面节点输出变量：
STATUS:接触状态
CPRESS:接触压强
CSHEAR1：局部方向1上的摩擦剪应力
COPEN:从面节点与主面的距离
CSLIP1:在局部方向1上的相对切向滑移。