法兰盖螺栓衬垫

1、ANASYS有限元分析题目：螺栓预紧法兰盘衬垫有限元分析姓名：老师：学号：学院：机电学院专业：机械工程螺栓预紧法兰盘衬垫有限元分析西南石油大学 ，成都 610500摘要：法兰是使管子与管子或阀门相互连接的零件，连接于管端。法兰上有孔眼，螺栓使两法兰紧连，法兰间用衬垫密封。法兰联接的螺栓预紧力是保证法兰联接点不发生泄漏的重要影响因素。给定预紧力，通过anasys 分析得出螺栓最适用于法兰螺栓连接系统的预紧力。关键词：法兰;螺栓预紧;衬垫;1. 引言1.1. 法兰连接概念法连接就是把两个管道、管件或器材，先各自固定在一个法兰盘上，两个法兰盘之间，加上法兰垫，用螺栓紧固在一起，完成了连接。有的管件和

2、器材已经自带法兰盘，也是属于法兰连接【1】。法兰盖，就是盖的边缘有法兰缘，可以用螺栓固定的法兰连接是管道施工的重要连接方式。 1.2. 螺栓预紧法兰盖螺栓法兰联接是压力容器，石油化工设备及管道应用极为广泛的一种可拆式静密封联接结构。法兰螺栓连接系统主要失效形式是泄露，而螺栓预紧和加装衬垫是保证连接面不发生泄露的重要环节之一。在螺栓连接中，螺栓在安装的时候必须拧紧，即在连接承受工作载荷之前，预受到力的作用，这个预加的力称为预紧力，预紧的目的在于增强连接的可靠性和紧密性，以防止受载后被连接件出现缝隙或发生相对滑移。螺栓预紧力不能过大或者过小。预紧力过小，达不到垫片压紧并实现初始密封条件；预紧力过大

3、，使垫片压坏或挤出。所以知道预紧力的大小值就变得尤其重要【2】。本文通过给定不同预紧力，通过anasys有限元分析，得出螺栓，衬垫，法兰的应力图。2. 理论分析及其有限元分析2.1. 有限元分析的目的和意义法兰为圆形盘状，设计螺栓组布置时要满足：在布置螺栓时要求螺栓的列应有合适的间距，边距；螺栓的布置应使各螺栓的受力合理；分布在同一圆周上螺栓数目，应取4,6,8等偶数；避免螺栓承受附加的弯曲载荷。螺栓工作时承受工作载荷和预紧力的联合作用，这种螺栓连接在紧螺栓连接中比较常见，因此也是最重要的一种。衬垫是介于法兰盖与法兰盖之间的填封介质。衬垫的材料以及预紧力不同，会对法兰盖的工作性能产生极大的影响

4、。螺栓预紧法兰盖衬垫分析，选4个螺栓对称布置，1个衬垫，2两法兰盖组成。螺栓和法兰盖由186非线性单元模拟，衬垫用194非线性单元模拟。约束边界为约束法兰内圈的径向和周向位移。与施加在4个螺栓上的力都为1000N。在anasys进行分析室，在法兰内径的径向施加约束。2.2. 模型理论分析有限单元法的基本思想是将连续的求解区域离散为一组有限个、且按一定方式相互联结在一起的单元的组合体。由于单元能按不同的联结方式进行组合，且单元本身又可以有不同形状，因此可以对复杂的模型进行求解。有限单元法作为数值分析方法的另一个重要特点是利用在每一个单元内假设的近似函数来分片地表示全求解域上待求的未知场函数。单元

5、内的近似函数通常由未知场函数或其导数在单元的各个节点的数值和其插值函数来表达。这样一来一个问题的有限元分析中，未知场函数或及其导数在各个节点上的数值就成为新的未知量，从而使一个连续的无限自由度间题变成离散的有限自由度问题。一经求解出这些未知量，就可以通过插值函数计算出各个单元内场函数的近似值，从而得到整个求解域上的近似解。显然随着单元数目的增加，也即单元尺寸的缩小，或者随着单元自由度的增加及插值函数精度的提高，解的近似程度将不断改进。如果单元是满足收敛要求的，近似解最后将收敛于精确解。目前，有限单元法的应用已由弹性力学平面间题扩展到空间问题、板壳问题，由静力平衡问题扩展到稳定间题、动力问题和波

6、动问题。分析的对象从弹性材料扩展到塑性、粘弹性、粘塑性和复合材料等，从固体力学扩展到流体力学、传热学等连续介质力学领域.在工程分析中的作用已从分析和校核扩展到优化设计并和计算机辅助设计技术相结合。2.3. Anasys分析步骤Anasys分析主要步骤：定义参数；选择单元；定义材料属性；定义材料压缩曲线；定义衬垫的线性卸载曲线；创建几何体；划分网格；定义预拉伸单元；定义接触单元节点组；定义目标单元节点组；通过接触向导建立接触对；约束法兰盘内径的径向位移；约束法兰盘内径的周向位移；施加预拉伸力；求解；后期处理。3. ANSYS计算的具体流程3.1. 分析模型在开始有限元计算前，需要对计算的工程问题

7、，进行认真的分析。这其中包括：模型的简化，能否忽略几何不规则性，能否把三维问题简化为平面问题；能否把一些载荷看作是集中载荷；能否把某些支撑看做固定的；判断模型材料的应力-应变关系，即采用线性材料分析还是非线性材料分析等。3.2. 选择单元根据第一步分析的结果选择满足条件的单元，ANSYS中有多种单元，点、面、实体等均有，而且每种单元又有较多类型，因此在具体选择时应结合使用者对软件的熟悉程度结合所分析研究模型进行选择。在单元选择时，可结合具体选项适当讲解各种典型单元的用途和常用的一些单元【4】。螺栓法兰盖单元设置步骤:Preprocessor ElementTypeAdd/Edit/Delete

8、Addsolid20node 186OK如图3-1。设置衬垫单元设置步骤：Preprocessor ElementTypeAdd/Edit/DeleteAddGasket 16node 194OKClose。两种单元设置完以后如图32所示。31设置186单元体32设置完两种单元体3.3. 设置材料属性材料属性需要结合所分析研究对象的特性来设置，主要包括泊松比、杨氏模量，属性材料包括屈服强度和切线模量（双线型材料），疲劳分析则有应力-工作次数关系曲线等，因此对构件的材料掌握非常关键。本出采用非线性线型材料来研究材料的特性。1）螺栓和法兰盖材料设置步骤：PreprocessorMaterialMo

9、delsStructural ElasticIsotropic OK,然在在弹出的对话框里输入弹性模量1.05E11，泊松比0.3。如图3-3所示。2）衬垫材料的设置步骤：PreprocessorMaterialModelsStructuralSpecialized MaterialsGasketCmpresionAdd point,如图3-4所示。把数值输入表中得到衬垫材料店压缩曲线如图3-5，再点Linaer unloain得到卸载曲线，击如图3-6,7所示。33螺栓法兰盖材料设置34衬垫材料设置35衬垫压缩曲线 36衬垫材料卸载数据点 37衬垫材料卸载特性曲线3.4. 建立模型ANSYS

10、中有四种建立模型方法，即实体建模，有限元建模，从其他CAD软件中导入和参数化建模。前三种建模方法适合除优化设计，可靠性分析的各种分析。建模过程：首先设置参数变量如图3-8所示。H1=0.01 !法兰盘1的厚度。H2=0.002 !衬垫厚度。H3=0.01 !法兰盘2的厚度。H4=0.008 !螺帽的厚度。R1=0.05 !法兰盘内半径。R2=0.08 !法兰盘外半径。R3=(R2-R1)/2+R1 !螺栓中心位置。R4=0.5*(R2-R1)/2!螺栓半径。R5=0.7*(R2-R1)/2!螺帽半径。38引入建模参数1）建立法兰盖和衬垫模型，通过圆柱命令建立不同高度的圆柱， 3-9所示图，输入

11、数字，建立1衬垫，2个法兰的模型。操作步骤：PreprocessorModelingCreateVolumesCylinderBy DimensionsOk 39法兰盘与衬垫三维模型2）通过布尔减操作，建立法兰盘和衬垫上的四个圆孔，如图3-10所示。操作步骤：PreprocessorModelingOperateBooleansSubtract选择建立的衬垫与法兰盖apply选择4个小的圆柱OK.3）通过创建圆柱操作，对建立的3个圆柱体进行布尔相加，得到螺栓，如图3-11所示。操作步骤：PreprocessorModelingCreateVolumesCylinderBy Dimension

12、4）把的螺栓通过平移到法兰盖的孔处，如图3-12。再把螺栓通过复制粘贴到其他三个法兰孔的位置，如图3-13。操作步骤：WorkplaneChange active CS toGloable Cylidric PreprocessorModelingCopyVolumes再在DY栏输入粘贴角度90,180,270。 310创建孔 311创建螺栓 312螺栓偏移 313四个螺栓创建完成3.5. 网格划分ANSYS中主要有两种网格划分方法，即自由网格和映射网格。自由网格划分的成功率高，但在动力学计算中精度稍差，用户根据需要进行选择使用。1）划分一个法兰盖网格。法兰盖网格采用映射网格。操作步骤:Pre

13、processorMeshingMeshtoolSet186单元Set网格尺寸0.005HexSweep选择一个法兰盖，如图3-14，3-15。 314法兰盖网格参数设置 2)划分衬垫网格。衬垫单元在两个体的中间，所以先划分的一个法兰，再衬垫，再划分另外一个法兰盖。操作步骤: PreprocessorMeshingMeshtoolSet186单元 Meshing Mesh Interface3D Interface选择41,42号面再在如图力输入值，完成网格划分，如图3-16,3-17所示。4）划分另一个法兰盖 另一个法兰盖的划分方法与第一法兰盖的划分网格方法一样，网格划分完如图3-18所示。

14、5)划分螺栓网格。螺栓网格采用自由网格划分。操作步骤:PreprocessorMeshingMeshtoolSet186单元TetFreeMesh选择4个螺栓Ok，如图3-19所示。 315 划分法兰盖 316衬垫网格参数设置317 衬垫单元网格划分完成法 318兰螺栓划分完成 319划分螺栓3.6. 确定分析类型这一步用户可以选择需要的分析类型：静力学分析，模态分析，谐响应分析，瞬态动力学分析，谱分析，特征值屈曲分析，子结构分析等，3.7. 施加边界条件这一步用户根据模型的实际工况定义边界条件，如对称边界，完全约束等。1）定义预拉伸单元操作步骤：:PreprocessorMoleingEle

15、mentspretensionPretensn MeshWith optionDivide at voluElements in volue。如图3-20 320定义预拉伸单元2）定义接触单元。螺栓底面与法兰盖表面由接触。创法兰表面的建节点组操作步骤：拾取法兰表面节点SelectEntiliesAreasBy PickFrom AllSele AllApplyOk,如图3-21所以。再SelectEntiliesNodesAttached toAreas AllFrom AllPlotOk,如图3-22所示。创建法兰表面节点组CN1，如图322所示。用同样创法兰给另一面的节点组CN2。 321

16、创建节点组322创建节点组创建螺栓底面节点组操作步骤：拾取法兰表面节点SelectEntiliesAreasBy PickFrom AllSele AllApply选择四个螺栓底面Ok如图3-23,再 Select Entilies Nodes Attached toAreas AllFrom AllPlotOk。 323拾取螺栓节点组 324螺栓和法兰节点组 创建接触步骤如下：Contact ManagerContace WizrdNodal ComponentTN1Nodal NexeCompentCN1Option SettingConcat OperitiesCreatFinsh.以同

17、样的方法创建第二接触。建立过程如图3-25,3-26,3-27，3-28所示。325创建接触326选择接触目标 327设置不对称 328 接触设置完成3.8. 求解用户在完成以上操作后，就进入求解阶段，ANSYS有直接求解，载荷步求解和自适应求解方法，根据需要适时选择。通过约束法兰内径的轴向和轴向位移。所以约束R1一周的节点。选取法兰内径节点：WokplaneChange WorkplaneChange active CS toGloable Cylidric SelectEntiliesAreasBy locationR1OK.过程设置如图。329选取法兰内环节点约束法兰内径节点：Prepr

18、ocessorSolutionApplyStrucuralOn NodesPick all Ux,Uy。如图3-30,3-31所示。 330设置步 331施加约束. 332约束后的显示 定义预拉伸载荷SolutionDefine LoadapplyStructuralPretnsn secth. 333设置载荷 334载荷设置过程求解SolveCurrent Ls3.9. 后处理在该步用户可以使用两种后处理方法，观察结果，校核计算结果是否符合工程结构设计要求，如果满足要求，则保存结果，如果不满足要求，则需要修改模型重新计算直到符合要求为止。 观察最大应力点：General PostprocPo

19、lt ResulContour PoltElemnet Solution。如图3-35所示。335整体变形图观察螺栓应力步骤：SelectEntiliesVolumBy pick螺栓OkSelectEntiliesElementAttacth toVolumOk,如图所示3-36所示。用同样方法观察法兰，衬垫。如图3-37,3-38所示。336螺栓应力图337法兰应力图338衬垫应力图4. 结论通过对螺栓衬垫法兰的的力学分学模型分析，再借助anasys的有限元计算，得出：衬垫在孔的位置处，应力最大。螺栓，在有形状急剧变的话的位置处，应力集中，应力最大。本模型还可以施加不同的预紧力，去观察螺栓，

20、衬垫，法兰的变形与应力状况。因此可以对施加的预紧力进行控制，得到理想的预紧力。参考文献1 王平.衬垫的压应力分析.石油矿厂机械，20012年2月，20(1):5051 2 2张宏英，李东阳.螺栓预紧力分析， 机械装备，2011，40(3)：2023.3 肖文生，ANASYS13.0基础与实例教程.2011年11月.4 秦宇.ANSYS11.0基础与实例教程，化学工业出版社，2009，21（2）：48.命令流通过对求解的模型导出命令流。通过修改命令流里面的参数，可以计算同模型下，不同问题的求解。命令流如下：/BATCH ! /COM,ANSYS RELEASE 14.5 UP20120918 0

21、9:30:28 01/03/2014 /input,start145,ans,D:ANSYS Incv145ANSYSapdl !* /FILNAME,xieshuai2,0*SET,H1,0.01 ! *SET,H2,0.002 ! *SET,H3,0.01 ! *SET,H4,0.008 ! *SET,R1,0.05 ! *SET,R2,0.08 ! *SET,R3,(R2-R1)/2+R1 ! *SET,R4,0.5*(R2-R1)/2! *SET,R5,0.7*(R2-R1)/2! /PREP7 !* ET,1,SOLID186 !* ET,2,INTER194 !* !* MPTEM

22、P, MPTEMP,1,0 MPDATA,EX,1,1.05E11MPDATA,PRXY,1,0.3 TB,GASKET,2,13,COMP ! TBPT,0.508000E-04, 0.161226E+07TBPT,0.101600E-03, 0.520884E+07TBPT,0.152400E-03, 0.113134E+08TBPT,0.203200E-03, 0.200499E+08TBPT,0.254000E-03, 0.259960E+08TBPT,0.304800E-03, 0.290345E+08TBPT,0.355600E-03, 0.357453E+08TBPT,0.406

23、400E-03, 0.440064E+08TBPT,0.457200E-03, 0.563189E+08TBPT,0.508000E-03, 0.748254E+08TBPT,0.558800E-03, 0.967287E+08TBPT,0.609600E-03, 0.129001E+09TBPT,0.683260E-03, 0.157147E+09TB,GASKET,2,5,LUNL ! TBPT,0.152400E-03, 2.430000E+11TBPT,0.304800E-03, 3.565000E+11TBPT,0.406400E-03, 5.923000E+11TBPT,0.558

24、800E-03, 1.088000E+12TBPT,0.683260E-03, 1.490000E+12!* CYLIND,R1,R2,0,H1,0,360,CYLIND,R1,R2,H1,H1+H2,0,360,CYLIND,R1,R2,H1+H2,H1+H2+H3,0,360, NUMMRG,ALL CYLIND,R4,0,0,H1+H2+H3,0,360, ! VGEN,4,R3,1 ! CSYS,1 ! ! VGEN,2,4,90,0 VGEN,2,4,180,0 VGEN,2,4,270,0 ! FLST,2,3,6,ORDE,2 FITEM,2,1 FITEM,2,-3 VSBV,

25、P51X,4 FLST,2,3,6,ORDE,2 FITEM,2,8 FITEM,2,-10 VSBV,P51X,5 FLST,2,3,6,ORDE,2 FITEM,2,1 FITEM,2,-3 VSBV,P51X,6 FLST,2,3,6,ORDE,3 FITEM,2,4 FITEM,2,-5 FITEM,2,8 VSBV,P51X,7 CYLIND,R4*0.95,0,0,H1+H2+H3,0,360, CYLIND,R5,0,H1+H2+H3,H1+H2+H3+H4,0,360, CYLIND,R5,0,-H4,0,0,360,! /VIEW,1,1,1,1 ! /ANG,1 ! /RE

26、P,FAST ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /AN

27、G,1,30,XS,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /ANG,1,30,XS,

28、1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST FLST,2,3,6,ORDE,2 FITEM,2,4 FITEM,2,-6 VADD,P51X CSYS,0 ! VGEN,7,R3,1 CSYS,1 VGEN,2,7,90,0 VGEN,2,7,180,0 VGEN,2,7,270,0 ! SAVE, xieshuai2,db, ! SAVE, xieshuai2,db, ! SAVE, xieshuai2,db, ! SAVE, 36,db,C:UsersAdministrator ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAS

29、T ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST TYPE, 1 MAT, 1REAL, ESYS, 0 SECNUM, !* ESIZE,0.005,0, /UI,MESH,OFFESIZE,0.005 TYPE,1 MAT,1 VSWEEP,2TYPE,2 MAT,2 IMESH,AREA,19,25, ,0,0,H2,0.0005TYPE,1 MAT,1 VSWEEP,3MSHKEY,0MSHAPE,1,3D VMESH,4 VMESH,5 VMESH,6 VMESH,7 PSMESH,

30、 ,P1, ,VOLU,4, 0,Z,(H1+H2+H3)/2, PSMESH, ,P2, ,VOLU,5, 0,Z,(H1+H2+H3)/2, , , , , PSMESH, ,P3, ,VOLU,6, 0,Z,(H1+H2+H3)/2, , , , , PSMESH, ,P4, ,VOLU,7, 0,Z,(H1+H2+H3)/2, , , , , ! /DIST,1,0.924021086472,1! /REP,FAST ASEL,ALLASEL,S, , , 48 NSEL,ALLNSLA,S,1! NPLOT NSLA,S,1CM,CN1,NODE ! /DIST,1,0.924021

31、086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08

32、222638492,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST ! vlist, all ! elist,all,0,0 ! /DIST,1,0.92402108

33、6472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.0822

34、2638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.

35、924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST ASEL,S, , ,28

36、 NSLA,S,1CM,CN2,NODE ALLSEL,ALL ! /MREP,EPLOT ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ASEL,ALLFLST,5,4,5,ORDE,4 FITEM,5,50 FITEM,5,54 FITEM,5,64 FITEM,5,74 ASEL,S, , ,P51X NSEL,ALLNSLA,S,1! NPLOT ! NPLOT NSLA,S,1! /MREP,EPLOT

37、CM,TN1,NODE FLST,5,4,5,ORDE,4 FITEM,5,70 FITEM,5,80 FITEM,5,53 FITEM,5,60 ASEL,S, , ,P51X NSEL,ALLNSLA,S,1CM,TN2,NODE ALLSEL,ALL ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,1.0822263

38、8492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /MREP,EPLOT ! /COM, CONTACT PAIR CREATION - START C

39、M,_NODECM,NODE CM,_ELEMCM,ELEM CM,_KPCM,KP CM,_LINECM,LINE CM,_AREACM,AREA CM,_VOLUCM,VOLU ! /GSAV,cwz,gsav,tempMP,MU,1,0.25MAT,1 MP,EMIS,1,7.88860905221e-031R,3 REAL,3 ET,4,170ET,5,174R,3,1.0,0.1,0,RMORE,1.0E20,0.0,1.0, RMORE,0.0,0,1.0,1.0,0.5RMORE,0,1.0,1.0,0.0,1.0KEYOPT,5,4,0KEYOPT,5,5,0NROPT,UNS

40、YM KEYOPT,5,7,0KEYOPT,5,8,0KEYOPT,5,9,0KEYOPT,5,10,2 KEYOPT,5,11,0 KEYOPT,5,12,0 KEYOPT,5,2,0KEYOPT,4,5,0! Generate the target surface NSEL,S,TN1CM,_TARGET,NODE TYPE,4 ESLN,S,0ESURF CMSEL,S,_ELEMCM ! Generate the contact surface NSEL,S,CN1CM,_CONTACT,NODETYPE,5 ESLN,S,0ESURF ALLSEL ESEL,ALLESEL,S,TY

41、PE,4 ESEL,A,TYPE,5 ESEL,R,REAL,3 ! /PSYMB,ESYS,1 ! /PNUM,TYPE,1! /NUM,1 ! EPLOT ESEL,ALLESEL,S,TYPE,4 ESEL,A,TYPE,5 ESEL,R,REAL,3 CMSEL,A,_NODECM CMDEL,_NODECM CMSEL,A,_ELEMCM CMDEL,_ELEMCM CMSEL,S,_KPCM CMDEL,_KPCM CMSEL,S,_LINECM CMDEL,_LINECM CMSEL,S,_AREACM CMDEL,_AREACM CMSEL,S,_VOLUCM CMDEL,_V

42、OLUCM ! /GRES,cwz,gsav CMDEL,_TARGET CMDEL,_CONTACT ! /COM, CONTACT PAIR CREATION - END ! /COM, CONTACT PAIR CREATION - START CM,_NODECM,NODE CM,_ELEMCM,ELEM CM,_KPCM,KP CM,_LINECM,LINE CM,_AREACM,AREA CM,_VOLUCM,VOLU ! /GSAV,cwz,gsav,tempMP,MU,1,0.25MAT,1 MP,EMIS,1,7.88860905221e-031R,4 REAL,4 ET,6

43、,170ET,7,174R,4,1.0,0.1,0,RMORE,1.0E20,0.0,1.0, RMORE,0.0,0,1.0,1.0,0.5RMORE,0,1.0,1.0,0.0,1.0KEYOPT,7,4,0KEYOPT,7,5,0NROPT,UNSYM KEYOPT,7,7,0KEYOPT,7,8,0KEYOPT,7,9,0KEYOPT,7,10,2 KEYOPT,7,11,0 KEYOPT,7,12,0 KEYOPT,7,2,0KEYOPT,6,5,0! Generate the target surface NSEL,S,TN2CM,_TARGET,NODE TYPE,6 ESLN,

44、S,0ESURF CMSEL,S,_ELEMCM ! Generate the contact surface NSEL,S,CN2CM,_CONTACT,NODETYPE,7 ESLN,S,0ESURF ALLSEL ESEL,ALLESEL,S,TYPE,6 ESEL,A,TYPE,7 ESEL,R,REAL,4 ! /PSYMB,ESYS,1 ! /PNUM,TYPE,1! /NUM,1 ! EPLOT ESEL,ALLESEL,S,TYPE,6 ESEL,A,TYPE,7 ESEL,R,REAL,4 CMSEL,A,_NODECM CMDEL,_NODECM CMSEL,A,_ELEM

45、CM CMDEL,_ELEMCM CMSEL,S,_KPCM CMDEL,_KPCM CMSEL,S,_LINECM CMDEL,_LINECM CMSEL,S,_AREACM CMDEL,_AREACM CMSEL,S,_VOLUCM CMDEL,_VOLUCM ! /GRES,cwz,gsav CMDEL,_TARGET CMDEL,_CONTACT ! /COM, CONTACT PAIR CREATION - END ! /MREP,EPLOT ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST CSYS,1 NSEL,ALLNSEL,S,LOC,X,R1 ! N

46、PLOT NSEL,S,LOC,X,R1 ! /REPLOT,RESIZE FLST,2,1024,1,ORDE,12 FITEM,2,2473FITEM,2,-2536 FITEM,2,2809FITEM,2,-3124 FITEM,2,3753FITEM,2,-3884 FITEM,2,6297FITEM,2,-6424 FITEM,2,9017FITEM,2,-9332 FITEM,2,9913FITEM,2,-9980 NROTAT,P51X ! /REPLOT,RESIZE FLST,2,1024,1,ORDE,12 FITEM,2,2473FITEM,2,-2536 FITEM,2

47、,2809FITEM,2,-3124 FITEM,2,3753FITEM,2,-3884 FITEM,2,6297FITEM,2,-6424 FITEM,2,9017FITEM,2,-9332 FITEM,2,9913FITEM,2,-9980 NROTAT,P51X ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST FINISH /SOLFLST,2,1024,1,ORDE,12 FITEM,2,2473FITEM,2,-2536 FITEM,2,2809FITEM,2,-3124 FITEM,2,3753FITEM,2,-3884 FITEM,2,6297FITEM

48、,2,-6424 FITEM,2,9017FITEM,2,-9332 FITEM,2,9913FITEM,2,-9980 !* /GO D,P51X, , , , , ,UX,UY, , , , ! /VIEW,1,1,1,1 ! /ANG,1 ! /REP,FAST ALLSEL,ALL CSYS,0 ! vlist, all ! elist,all,0,0 ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /D

49、IST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST

50、! vlist, all ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,0.924021086472,1! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /ANG,1,30,XS,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,1 ! /REP,FAST ! /DIST,1,1.08222638492,