ANSYS命令流学习笔记3-Solid单元的接触分析.docx

!ANSYS命令流学习笔记3-Solid单元的接触分析背景说明：两个半径均为100mm，的正交圆柱体发生正接触，作用在两圆柱接触体法线方向上的压力总和1000N，两圆柱体均为钢制分析两圆柱的接触情况。分析思路：对模型进行对称约束，施加载荷后要进行节点自由度耦合，使加载面的所有节点能在承载后具有相同的Y方向位移，并且在变形后仍然保持为水平平面。总结：1. 如果不进行节点耦合，结构变形会产生刚性位移，造成不收敛。但是节点耦合造成了变形不符合实际情况，即忽略了施力面的泊松比造成的影响。2. 由面网格生成体网格，几何只有面，有限元模型是体，需要定义网格单元，定义拉伸的目标单元，对拉伸选项做出设置。3. 接触部分的细分十分重要，也是收敛的影响因素。可以把下列文字，直接粘贴到TXT文档修改，更为方便。! 本次学习重点：!1、命令流的复杂建模!重点学习下建模。因为前处理尽可能用WB，所以这里也是规则模型。还是太不好用了。我选择狗带。!2、接触的定义。!Ansys依据实常数识别接触对，接触对要具有同样的实常数定义。而且命令流是根据选择的单元来定义接触单元。GUI操作，也不麻烦。!3、分析步数的设定!非线性分析时步数设置，线性搜索，自动步长，定义子步。!求解解析解finish/clearf=1000e=200r=100 !命令不区分大小写，参数也不区分大小写。p=2.45*(f*e*2*(r+r)/(r*r)*2)*(1/3) !接触应力的解析解61.73，有限元法结果59.9MPa!正式分析命令流finish/clear/prep7 !进入前处理et,1,solid185et,2,mesh200keyopt,2,1,6 !查看help中关于mesh200和keyopt的解释，KEYOPT, ITYPE, KNUM, VALUE。此处意思为定义了mesh(2)为4 nodes的四边形单元(1,6)。mp,ex,1,2e11mp,prxy,1,0.3 !定义et 1的材料属性，Elastic moduli 和Poissons ratios，（lab，mat，co） !lab这里是系统规定的，输出结果表格时lab是可以自定义的。cyl4, , , 100, , ,90 !在原点，创建一个r100的1/4圆面cyl4, ,200,100, , ,-90 !在(0,200)，创建一个r100的1/4圆面wpoff,10 !workplane向X移动10wprota, , ,90 !workplane绕Y轴旋转90度asbw,all !用workplane切分面esize,10,0 !定义线的网格尺寸10，这里0表示均分，定义了尺寸后，此值无意义mshape,0,2D !划分为四边形的2D网格mshkey,0 !自由划分amesh,all !划分面网格ksel,s, , ,2 !选择关键点2ksel,a, , ,4 !再选关键点4，这时候选了KP2和4krefine,all, , ,3 !关键点2和4处细化，细化等级为3CSYS,5 !转化坐标系，Y轴为旋转轴的圆柱坐标系。特烦坐标系换来换去的asel,s, , ,4,6,2 !选择4、6两个面agen, ,all, , , ,90, , , ,1 ! 面2和4绕Y旋转90度。gen表示generate,CSYS0下就是平移CSYS,0 !转换坐标系为cartesian坐标系asel,s, , ,4,6,2 !选择4、6两个面arsym,z,all, , , ,0,0 !对称生成面，xy平面对称，0,0是保留原面，保留网格信息asel,s, , ,3,5,2 !选择3、5两个面arsym,x,all, , , ,0,0 !对称生成面，yz平面对称，0,0是保留原面，保留网格信息eplot/view,1,1,1,1/replotasel,s, , ,1,2asel,a, , ,4,6,2 !选择面1,2,4,6agen, ,all, , ,50, ,50, , ,1 !平移选择面，x50，y50asel,s, , ,1,2asel,a, , ,4,6,2 !选择面1,2,4,6allsel,below,areanummrg,nodenummrg,kp !合并1246面上的节点allsel,all !此处全选，为了进行下一步的选择工作asel,s, , ,3,7,2asel,a, , ,8 !选择面3,5,7,8allsel,below,areanummrg,nodenummrg,kp !合并1246面上的节点allsel,alltype,1mat,1extopt,esize,10,0 !均分10份extopt,aclear,0extopt,attr,0,0,0 !定义生成体的特性asel,s, , ,3,7,2asel,a, , ,8 !选择面3,5,7,8vext,all, , , , , 100, , , , !向Z拉伸100asel,s, , ,1,2asel,a, , ,4,6,2 !选择面1,2,4,6vext,all, , , -100, , , , , , !向-X拉伸100!建模完成，开始定义接触单元mp,mu,1, !et1的摩擦系数为0 不是两者之间的摩擦系数？r,3 !定义实常数3，是为了让170和174单元同属一个实常数，识别他们的接触对。因为ansys里同一个接触对的单元实常数相同。real,3 !选择实常数3et,3,170 !定义目标单元TARGE170et,4,174 !定义接触单元CONTA174keyopt,4,9,0 !包括由几何和接触面偏移产生的初始穿透，也是缺省值0.keyopt,4,10,0 !禁止闭合的单元进行刚度修正，也是缺省值0. 可以不用定义，此处为了方便对接触的理解。asel,s, , ,27 !ASEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP。sel的命令要熟记asel,a, , ,31,33,2 asel,a, , ,37 !选择面27,31,33,37type,3nsla,s,1 !选择面和线上的所有点，节点，关键点esln,s,0 !选择已选择节点相关的所有单元esurf !创建TARGE170asel,s, , ,10 asel,a, , ,14asel,a, , ,20 asel,a, , ,23 !选择面10,14,20,23type,4nsla,s,1 !选择面和线上的所有点，节点，关键点esln,s,0 !选择已选择节点相关的所有单元esurf !创建CONTA174allsel !方便下一步操作!定义接触完成，开始加载约束条件asel,s, , ,9,21,4 !选择面9,13,17,21da,all,uz, !约束uzasel,s, , ,12 asel,a, , ,16asel,a, , ,18 asel,a, , ,22 !选择面12,16,18,22da,all,uydl,26, ,ux !约束面uy，l26ux。DL, LINE, AREA, Lab, 比da多了一个面选择。asel,s, , ,1,2 asel,a, , ,4,6,2 !选择面1,2,4,6da,all,uxdl,47, ,uz !约束面ux，l47uz。asel,s, , ,26,38,4 !选择面26,30,34,38nsla,s,1 !选择面上所有点,为了后续耦合自由度sfa,all,1,pres,0.05 !surface loadscp,1,uy,all !耦合自由度？为啥耦合自由度？难道施加压力的方向不能作为依据？结构不对称造成模型节点异常，必须cp加点位移allsel,allfinish/solulnsrch,1 !线性搜索 Ui+1 = Ui + S*D,通过自动调整S值，提高迭代效率。nsubst,5,10,3 !number of substeps,5步，最大10，