第四章 点－面接触单元.ppt

! *** 创建接触单元 *** 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 ! *** 求解控制、求解 *** /solu solcontrol,0 antype,static cnvtol,f,,,,-1 nlgeom,on nsubst,6 outres,,1 ! 求解结果输出设置 d,ncen,uy,-0.1 ! 对耦合节点施加载荷uy=-0.1 solve fini ! *** 后处理器 *** /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 例3 卡头与卡座的接触分析(图4-4) 目标：计算卡头插进卡座和拉出卡座所需的力。此问题属于带厚度的平面应力问题。图4-4 建模要点：卡头和卡座的底板被认为是刚性的，因此建模时不考虑。由于对称，用右半部分来建模和计算。 卡头与卡座板厚＝5 材料: E = 2.8e3, = 0.3 MU = 0.2 (摩擦系数) 图4-4 文件: IPROG.DAT Step 1.建模：定义单元模型、实常数、材料特性、分网 /prep7 et,1,plane42,,,3 ! 带厚度的平面应力 et,2,plane42,,,3 mp,ex,1,2.8e3 r,1,5 et,3,48,,,1 ! contac48 2D点－面接触单元,KEYOPT(3)=1弹性库仑摩擦 k,1,10 $k,2,20 $k,3,15,18.5 $k,4,10,20k,5,12.5,30 $k,6, 20,30 l,1,3 $l,4,3 $l,4,5 $l,5,6 $l,6,2 $l,2,1 al,all ! A1 lgen,2,1,3,1 ! L7,L8,L9 k,11,5 $k,12,5,30 l,11,7 $l,11,12 $l,12,10 ! L10, L11, L12 lsla,s ! 选L1, L2, L3, L4, L5, L6 lsel,invert al,all lsel,all recting,0,15,0,10 ! A3 asba,2,3 ! 从A2-A3→A4 recting,2,4,,,0,30,0,,,0,1 asel,s,,,4 ! 选A4 aatt,1,1,2 ! MAT=1,Real=1 ! Type=Plane42asel,all esize,4 amesh,all ! A1、A4分网 fini /solu nsel,s,loc,y,0 d,all,all nsel,s,loc,y,60 d,all,ux nsel,all fini save,IPROG,db ! *** 以下进行接触分析 *** Step 2.识别接触对 L3→L17、L9 L2→L8 ! *** 创建4个节点组元 *** ①将L9和L17线上的节点定义为组元“snapins”（卡头插入） lsel,s,,,9 ! 选L9 lsel,a,,,