fini
/clear
*set,RR,2
/title,cube
/config,nres,1000000
/prep7
et,1,plane182
KEYOPT,1,3,1 !轴对称问题
mp,ex,1,210000
mp,nuxy,1,0.3
/PNUM,KP,1
/PNUM,LINE,1
/PNUM,AREA,1
!建立外筒的模型
rectng,0,10,,30
rectng,9.8,39.8,-50,-0.5
wpoff,,15
cyl4,,,5
asba,1,3
lfillt,2,1,RR
al,4,10,11
asba,4,1
lsel,s,,,6,8,2
lesize,all,,,20
lsel,s,,,5,7,2
lesize,all,,,10
amesh,2
allsel
lsel,s,line,,12,14,1
lsel,a,line,,9
lesize,all,,,5
lsel,s,line,,1,3,2
lesize,all,,,10
lsel,s,line,,2
lesize,all,,,15
lsel,s,line,,4 !圆弧部分
lesize,all,,,5
allsel
mshkey,1 !定义使用映射网格划分
amap,3,1,4,3,13
save,conta,db
!接触单元设置
allsel
et,2,conta172 !接触单元,O环
et,3,targe169 !目标单元,上下法兰
mp,mu,1,0.2
r,10,,,0.1
real,10
keyopt,2,2,0
keyopt,2,9,0
keyopt,2,12,0
lsel,s,line,,2,4,2
cm,contact,line
type,2
mat,1
nsll,s,1
esurf,all
allsel
lsel,s,line,,8
cm,target,line
type,3
mat,1
nsll,s,1
esurf,all
allsel
fini
/solu
allsel
dl,5,2,uy,0
dl,5,2,ux,0
!nsel,s,loc,x,0
!d,all,ux,0
nsel,s,loc,y,30
d,all,uy,-40
antype,static
nlgeom,on
nsubst,20
outres,basic,all
allsel
solve
fini
/post1
plnsol,s,eqv |