一个用penalty method做的接触例子

mututu_gui

加载方法比较特殊一些，采用循环方式。
/prep7
et,1,185,,1
et,2,173,,1
et,3,170
keyopt,2,10,1
mp,ex,1,2e5
mp,nuxy,1,0.3
tb,biso,1
tbdata,1,200,1000
r,2,5,,0.1
block,-1,1,-1,1,0,90
esize,,3
lesize,2,,,2
lesize,9,,,90
vmesh,all
n,10000,-3,-6,0
n,10001,20,-6,0
type,3
real,2
tshap,cyli
e,10000,10001
asel,s,,,3
nsla,s,1
type,2
esurf

/solu
nlgeom,on
solcontrol,on
csys
asel,s,,,1
nsla,s,1
d,all,all,0
alls
dl,8,5,ux,15
*do,i,1,90,3
nsel,s,loc,y,1
sf,all,pres,0
nsel,r,loc,z,i-2,i
sf,all,pres,100+i
allsel
nsub,3,100,1
solve
*enddo

lixia_pla

再次向达人学习