附上命令
/prep7
et,1,solid45
mp,ex,1,1e10
mp,nuxy,1,.2
et,2,conta173
et,3,targe170
mp,mu,2,0.2
et,4,mesh200
keyopt,4,1,6
keyopt,2,9,0
!keyopt,2,2,4
keyopt,2,4,2
r,1,,,.1,0.1 !定义法向刚度因子为1
!建模
/view,,1,1,1
k,1,0,0,0
k,2,0,10,0
k,3,0,50,0
k,4,10,0
k,5,50,0,0
larc,2,4,1,10
larc,3,5,1,50
l,2,3
l,4,5
al,1,2,3,4
!用mesh200划分面
type,4
lesize,all,,,8
mshape,0,2d
mshkey,1
amesh,all
!拉伸生成体,共拉伸80
esize,,12
type,1
mat,1
vext,all,,,0,0,80
!生成小柱
*get,kmax,kp,,num,max
*get,lmax,line,,num,max
lsel,none
k,kmax+1,0,0,200
k,kmax+2,11,0,200
k,kmax+3,0,11,200
l,kmax+1,kmax+2
l,kmax+1,kmax+3
larc,kmax+2,kmax+3,kmax+1,11
asel,none
al,lmax+1,lmax+2,lmax+3
!用mesh200划分面
type,4
lesize,all,,,6
mshape,0,2d
mshkey,1
amesh,7
!拉伸生成体,共拉伸100
esize,,16
type,1
mat,1
asel,s,,,7
vext,all,,,0,0,-120
allsel,all
!生成目标面
type,3
mat,2
real,1
tshap,quad
asel,s,,,3
cm,_target,area
nsla,s,1
!esln,s,0
esurf
!生成接触面
type,2
mat,2
real,1
tshap,quad
asel,s,,,10
cm,_contact,area
nsla,s,1
!esln,s,0
esurf
allsel,all
/expand,4,polar,half,,90
/view,,1,1,1
eplot,all
fini
/solu
allsel,all
antype,static,new
outpr,all,1
outres,all,1
!NROPT,SYMm
nlgeom,on
solcontrol,on
rescontrol,define,all,-1
autots,off
nsubst,20
!施加对称边界条件
asel,s,,,4,6,2
da,all,symm
asel,s,,,9,11,2
da,all,symm
!nsel,s,loc,z,0
!d,all,all
!allsel,all
da,5,all,0
da,7,ux,0
da,7,uy,0
da,7,uz,-10
allsel,all
eplot,all
solve |