高级接触问题.pptVIP

模型分析：本例为三维实体接触，选用HYPER58（三维超弹性实体单元）、CONTAC49（三维点一面接触单元）进行分析。由于模型对称性可以采用简化的模型进行分析：用1/4园（平面），用2D模型进行计算，选用HYPER56（二维超弹性实体单元）、CONTAC26（二维点－基础接触单元）。 目标：介绍一种二维点－基础接触单元进行点－面接触分析。用该单元建立接触点与目标面的方法不使用GCGEN命令，而采用直接法用E命令建立单元。 CONTAC26 2D Point-to-Ground Contact (图4-3) 该单元在面的法向方向具有承受压缩的能力，在切线方向具有承受剪切（库仑摩擦）的能力。有三个节点I, J, K，2个自由度UX, UY。面可以是直线或圆弧。 R,npt,STIFN, … Mp,npt,MU（Keyopt(1) = 1时，应用库仑摩擦） 图4-3 第29页，共62页，2022年，5月20日，21点49分，星期四 文件：C_N_S.dat R = 0.2 ! 半径参数变量 /prep7 et,1,HYPER56 ! 超弹性实体单元，基体分网使用 et,2,contac26 ! 二维点－基础接触单元 r,2,2000 mp,ex,1,2.82 mp,nuxy,1,0.49967 tb,mooney,1 tbdata,1,0.293 tbdata,2,0.177 csys,1 ! 设置为柱坐标系 k,1 $k,2,r,-90 $k,3,r $k,4,(0.5*r),-90 k,5,(0.6*r),-45 $k,6, (0.5*r) $k,7,r,-45 l,2,7 $l,7,3 csys,0 ! 重新设置为笛卡尔坐标系 a,2,7,5,4 $a,7,3,6,5 $a,4,5,6,1 esize,,4 amesh,all ! 基体分网 第30页，共62页，2022年，5月20日，21点49分，星期四 ! *** 创建接触单元 *** n,1001,(-2*r),-r ! 定义目标节点J n,1002,(2*r),-r ! 定义目标节点K csys,1 ! 将坐标系设置为柱坐标 type,2 ! 接触单元CONTAC26 real,2 *do,j,1,8 nodc=node(r,90*((j-1)/8-1),0) ! 生成CONTAC26的I节点证基体的圆弧共8个点 e,nodc,1001,1002 ! 生成接触单元 *enddo csys,0 nsel,s,loc,x ! *** 选x=0的节点 *** d,all,ux ! 将x=0的节点x轴方向的位移置为0 nsel,s,loc,y ! 选y=0的节点 cp,1,uy,all ! 耦合y=0的节点在 轴方向的位移 *get,ncen,node,,num,min ! 提取所选节点中最小的节点编号 nsel,all fini 第31页，共62页，2022年，5月20日，21点49分，星期四 ! *** 求解控制、求解 *** /solu solcontrol,0 antype,static cnvtol,f,,,,-1 nlgeom,on nsubst,6 outres,,1 ! 求解结果输出设置 d,ncen,uy,-0.1 ! 对耦合节点施加载荷uy=-0.1 solve fini 第32页，共62页，2022年，5月20日，21点49分，星期四 ! *** 后处理器 *** /post1 /dscale,1,1 pldisp,1 ! 绘变形图 /post26 ! 进入历程后处理器 /axlab,y,force ! 标注y轴 /axlab,x,displacement ! 标注x轴 nsol,2,ncen,u,y rforce,3,ncen,f,y prod,2,2,,,,,,,-2 prod,3,2,,,,,,,-2 xvar,2 plvar,3 ! 绘制力－位移曲线 plvar,2,3 ! 列表显示力、位移取值 fini 第33页，共62页，2022年，5月20日，21点49分，星期四 第34页，共62页，2022年，5月20日，21点49分，星期四 例3 卡头与卡座的接触分析(图4-4) 目标：计算卡头插进卡座和拉出卡座所需的力。此问题属于带厚度的平面应力问题。图4-