文档介绍：!ANSYS命令流学****笔记3-Solid单元的接触分析背景说明:两个半径均为100mm,的正交圆柱体发生正接触,作用在两圆柱接触体法线方向上的压力总和1000N,两圆柱体均为钢制分析两圆柱的接触情况。分析思路:对模型进行对称约束,施加载荷后要进行节点自由度耦合,使加载面的所有节点能在承载后具有相同的Y方向位移,并且在变形后仍然保持为水平平面。总结:,结构变形会产生刚性位移,造成不收敛。但是节点耦合造成了变形不符合实际情况,即忽略了施力面的泊松比造成的影响。,几何只有面,有限元模型是体,需要定义网格单元,定义拉伸的目标单元,对拉伸选项做出设置。,也是收敛的影响因素。可以把下列文字,直接粘贴到TXT文档修改,更为方便。!本次学****重点:!1、命令流的复杂建模!重点学****下建模。因为前处理尽可能用WB,所以这里也是规则模型。还是太不好用了。我选择狗带。!2、接触的定义。!Ansys依据实常数识别接触对,接触对要具有同样的实常数定义。而且命令流是根据选择的单元来定义接触单元。GUI操作,也不麻烦。!3、分析步数的设定!非线性分析时步数设置,线性搜索,自动步长,定义子步。!求解解析解finish/clearf=1000e=200r=100!命令不区分大小写,参数也不区分大小写。p=*(f*e**2*((r+r)/(r*r))**2)**(1/3)!,!正式分析命令流finish/clear/prep7!进入前处理et,1,solid185et,2,mesh200keyopt,2,1,6!查看help中关于mesh200和keyopt的解释,KEYOPT,ITYPE,KNUM,VALUE。此处意思为定义了mesh(2)为4nodes的四边形单元(1,6)。mp,ex,1,2e11mp,prxy,1,!定义et1的材料属性,Elasticmoduli和Poisson'sratios,(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