ANSYS中正确地模拟过盈配合

在ANSYS中正确地模拟过盈配合过盈配合在机械产品的装配中使用的相当普遍。比如轴与轴承、轴与轴瓦、汽车的制动盘等，都是通过一定的过盈量来使两个装配部件紧密地连接起来。下面讨论如何在ANSYS中正确地模拟过盈配合。过盈配合在有限元分析中是一种典型的非线性接触行为。在有限元分析中设定了接触，从本质上来讲就是对相互接触的两个部件施加了某种约束，不同的接触算法对于接触约束的处理方法有所不同。接触约束的理论算法的选择，在ANSYS中是通过设置contact单元的KEOPT(2)选项来实现的。在ANSYS中目前主要有5种接触约束算法：KEYOPT(2)=0AugmentedLagrangian-加强的拉格朗日算法，这是ANSYS的缺省选择；KEYOPT(2)=1Penaltyfunction-罚函算法；KEYOPT(2)=2Multipointconstraint(MPC)-多点约束算法；KEYOPT(2)=3Lagrangemultiplieroncontactnormalandpenaltyontangent-接触法向采用拉格朗日乘子，接触切向采用罚函数的综合算法。KEYOPT(2)=4PureLagrangemultiplieroncontactnormalandtangent-法向和切向均采用拉格朗日乘子算法。各种不同的约束算法各有其优缺点，各有各自最适用的场合，具体情况需要具体对待。大部分情况下，默认选择KEYOPT(2)=0就够用了。过盈配合所致的接触分析的难点在于如何确定初始接触状态。初始接触状态设置得不对，会导致错误的计算结果或者不准确的计算结果，下面举两个例子来说明。例1．两个圆柱体在几何上是刚好接触，划分网格后有限元模型有间隙。如图1所示。这两个圆柱体，在几何上是刚好相切的，即处于几何上刚好接触的初始状态。划分网格后，由于在圆周上用小段直线代替了弧线，两个圆柱体之间产生了一定的间隙，两个圆柱体的有限元模型的初始状态不再是接触的。此时，如果接触参数设置不当，就会因为初始约束不足，圆柱体出现刚体位移，得到错误的结果。(说明：例1本来与设置过盈量是无关的，为了说明初始接触状态的重要性顺带说说。)例2．有的人把两个接触部件的几何位置设定一定的过盈量，想用这个过盈量来模拟过盈配合，这种做法是错误的，几何上的过盈量不等于划分网格后有限元模型的实际过盈量。下面的图2中，是一个孔类零件和一个轴类零件的截面图，轴和孔在几何位置上预设了过盈量。(内圈的红色圆是孔边界，外圈的蓝色圆是轴边界，轴和孔在几何上是相互侵入的)。在几何上，图2的轴和孔有一定的过盈配合量，其大小等于两个圆的半径之差，我们的本意是想用这个几何位置上的过盈量来模拟过盈配合。不幸的是，两个部件划分网格之后，实际的过盈量应该为单元之间的距离，即图中靠得比较近的两条线段之间的距离，显然，这个距离不再等于我们预先设置的过盈量了。更何况，上面这个图还是两个部件的网格对得比较整齐的情况，如果网格对得不整齐，过盈量就和我们预设的差得更远了。对于过盈配合来讲，过盈量的数值变化对于过盈产生的应力的影响是很大的。在ANSYS中，要正确的设置过盈配合，主要分3步：第一步：设置KEYOPT(9)=4KEYOPT(9)的默认值为0，意思是既考虑两个接触部件由于初始几何位置造成的初始侵入量(或者间隙)，同时也考虑CNOF参数设置的偏移量。意即接触部件的初始接触状态是由CNOF和初始侵入量(或间隙)共同决定的。在这种情况下，两个接触部件的初始几何位置对初始接触状态是有影响的，这对于准确设置过盈量是很不利的。(前面例2已经说明了通过几何位置设置初始过盈量是不准确的)。在设置了KEOPT(9)=4之后，程序在计算初始接触状态的时候就只考虑CNOF的设置值，不考虑接触部件的几何位置造成的侵入或间隙，而且过盈量是以ramp方式施加的。(ramp施加方式即逐步施加)。第二步：通过设置Icont实常数划分网格后，通常情况下，Targetsurface和contactsurface上的单元之间会有间隙或者过盈量，如果间隙或者过盈量在Icont设定误差范围内，间隙或者过盈量会被消除掉，程序会使contactsurface和targetsurface上的单元处于刚好接触的状态。这个值的的具体设置可以参见帮助文档，本文中设置为0.2。第三步：通过设置实常数CNOF来设置过盈量在第二步中，通过Icont的设置，已经使得Contactsurface上的单元和Targetsurface上面的单元处于刚好接触的位置了，此时再设置CNOF，CNOF的值就是过盈量。(CNOF的本意并不是过盈量，只是在有了前面的设定后，它的值就是我们所要的过盈量，其具体含义请参考ANSYS的帮助文档)。下面这个例子实际上是一本ANSYS书上的一个例子，这个例子的PDF版本在网上流传甚广，但是原书上的分析结果是错误的，具体错误之处，将在后面提及。例3.一个简单的轴和带孔圆盘的过盈配合的实例。圆盘的基本尺寸为：内径Rpin=35mm(原书中此值为34mm)，外径Rpout=100mm，盘高Hp=25mm；轴的基本尺寸为：内径Rain=25mm，外径Raout=35mm，轴长La=150mm。（原书中圆盘孔内径为Rpin＝34mm，和轴在几何上形成1mm的过盈量。由于结构是完全轴对称的，故可只取四分之一模型分析之。本例分析中，取过盈量f=0.01mm，而且本例仅仅计算由于过盈配合所产生的应力。按照本例各个物理量所取的单位，最终的计算结果中，应力单位应该为MPa;完整的命令流如下：Finish/clear,start/TITLE,Contactanalysiswithinitialinterference/PREP7!带孔圆盘的基本尺寸；Rpin=35Rpout=100Hp=25!轴的基本尺寸；Rain=25Raout=35La=150!过盈量f;f=0.01!实体的单元类型为带中间节点的2阶六面体单元；ET,1,solid186MP,EX,1,2.1E5!弹性模量；MP,PRXY,1,0.3!poisson系数；!生成带孔圆盘的1/4实体模型；CYL4,0,0,Raout,0,Rpout,90,Hp!轴的1/4实体模型；CYL4,0,0,Rain,0,Raout,90,La!把轴的位置沿着轴向移动一段距离；VGEN,,2,,,,,-10,,,1!!!!************对实体划分网格****************************!*LESIZE,17,,,15,,,,,1LESIZE,19,,,15,,,,,1!*LESIZE,18,,,2,,,,,1LESIZE,20,,,2,,,,,1!*LESIZE,22,,,20,,,,,1!*LESIZE,5,,,10,,,,,1LESIZE,7,,,10,,,,,1!*LESIZE,6,,,8,,,,,1LESIZE,8,,,8,,,,,1!*LESIZE,10,,,3,,,,,1!*VSWEEP,ALL!*************网格划分完毕**********************************!!/COM,CONTACTPAIRCREATION-STARTMP,MU,1,0.2MAT,1R,3REAL,3ET,2,170ET,3,174R,3,,,0.2,0.2,0.9,0!Icont系数设置为0.2;!设置过盈量为f;RMORE,,,1.0E20,f,1.0,0RMORE,0.0,0,1.0,,1.0,0RMORE,0,1.0,1.0,0.0,,1.0KEYOPT,3,4,0KEYOPT,3,5,0NROPT,UNSYMKEYOPT,3,7,0KEYOPT,3,8,0!KEYOPT(9)=4;!不考虑初始几何位置造成的过盈或者间隙，!只考虑CNOF参数设置的值，即过盈量。KEYOPT,3,9,4KEYOPT,3,10,2KEYOPT,3,11,0KEYOPT,3,12,0KEYOPT,3,2,0KEYOPT,2,5,0!GeneratethetargetsurfaceASEL,S,,,4CM,_TARGET,AREATYPE,2NSLA,S,1ESLN,S,0ESLL,UESEL,U,ENAME,,188,189ESURFCMSEL,S,_ELEMCM!GeneratethecontactsurfaceASEL,S,,,9CM,_CONTACT,AREATYPE,3NSLA,S,1ESLN,S,0ESURFALLSELESEL,ALLESEL,S,TYPE,,2ESEL,A,TYPE,,3ESEL,R,REAL,,3/PSYMB,ESYS,1/PNUM,TYPE,1/NUM,1EPLOTESEL,ALL/COM,CONTACTPAIRCREATION-ENDFINISH/SOLFLST,2,4,5,ORDE,4FITEM,2,5FITEM,2,-6FITEM,2,11FITEM,2,-12DA,P51X,SYMM!设置对称约束；FLST,2,1,5,ORDE,1FITEM,2,3/GODA,P51X,ALL,!*ANTYPE,0NLGEOM,1!必须打开大变形效应；NSUBST,1,0,0AUTOTS,0TIME,100/STATUS,SOLUSOLVEFINISH/POST1SET,1,LAST,1,PLNSOL,s,eqv!查看单元的VonMisesstress;esel,s,ename,,174PLNSOL,CONT,PRES,0,1.0!在contact单元上查看contactpressure;SAVE计算结果：说明：修改f值即可修改过盈量。如果过盈量设置的过大，接触部件可能会产生塑性变形，如果要考虑塑性变形，则应该修改材料属性定义为塑性材料。