裂縫計算沒有顯示計算完成是為什么?

如圖所示的結構進行有限元模擬計算,按照下面的命令流進行計算,結果點擊計算后最后未出現solution is done。我想知道為什么?求大神們給我解答。
finish
/clear
/filname,L
/title,liefeng
/prep7
*afun,deg
*set,length,3
*set,width,3
*set,rad,0.05
*set,ahi,0.02
*set,siga,2e6
*set,ksiga,4e6
*set,p,12e6 !水壓
!材料屬性
et,1,plane182
!keyopt,1,3,2
mp,ex,1,5.2e10
mp,prxy,1,0.3
r,1,0.002
!有限元建模
wplane,,length/2,width/2,,length/2+(rad+ahi)/sqrt(2),width/2+(rad+ahi)/sqrt(2),0,length/2-(rad+ahi)/sqrt(2),width/2+(rad+ahi)/sqrt(2),0
cswpla,11,0
rectng,-.1,.1,0,0.06
rectng,-.1,.1,-0.06,0
cyl4,rad+ahi,0,0.002,180
cyl4,rad+ahi,0,0.005,180
asba,1,4,,delete,keep
asba,4,3,,delete,keep
cyl4,rad+ahi,0,0.002,-180
cyl4,rad+ahi,0,0.005,-180
asba,2,6,,delete,keep
asba,6,4,,delete,keep
ksel,s,loc,x,rad+ahi,2*(rad+ahi)
nummrg,kp
allsel,all,all
cyl4,rad,0,0.002,180
cyl4,rad,0,0.005,180
asba,5,8,,delete,keep
asba,8,6,,delete,keep
cyl4,rad,0,0.002,-180
cyl4,rad,0,0.005,-180
asba,7,10,,delete,keep
asba,10,8,,delete,keep
ksel,s,loc,x,rad,rad+ahi/2
nummrg,kp
allsel
cyl4,-(rad+ahi),0,0.002,180
cyl4,-(rad+ahi),0,0.005,180
asba,9,12,,delete,keep
asba,12,10,,delete,keep
cyl4,-(rad+ahi),0,0.002,-180
cyl4,-(rad+ahi),0,0.005,-180
asba,11,14,,delete,keep
asba,14,12,,delete,keep
ksel,s,loc,x,-(rad+ahi),-2*(rad+ahi)
nummrg,kp
allsel
cyl4,-rad,0,0.002,180
cyl4,-rad,0,0.005,180
asba,13,16,,delete,keep
asba,16,14,,delete,keep
cyl4,-rad,0,0.002,-180
cyl4,-rad,0,0.005,-180
asba,15,18,,delete,keep
asba,18,16,,delete,keep
ksel,s,loc,x,-rad,-(rad+ahi/2)
nummrg,kp
allsel
cyl4,0,0,rad,360
asba,17,18
cyl4,0,0,rad,360
asba,19,17
cyl4,0,0,rad,360
asba,5,17
cyl4,0,0,rad,360
asba,6,5
cyl4,0,0,rad,360
asba,7,5
cyl4,0,0,rad,360
asba,8,5
cyl4,0,0,rad,360
asba,13,5
cyl4,0,0,rad,360
asba,14,5
cyl4,0,0,rad,360
asba,15,5
cyl4,0,0,rad,360
asba,16,5
wplane,,0,0,,10,0,0,0,10,0
rectng,length/2-0.5,length/2+0.5,width/2-0.5,width/2+0.5
cyl4,length/2,width/2,rad,360
asba,5,16
allsel
asba,21,all,,delete,keep
rectng,0,length,0,width
cyl4,length/2,width/2,rad,360
asba,16,21
allsel
asba,22,all,,delete,keep
csys,0
!網格劃分
kscon,11,0.0005,0,8,1.0
LESIZE,1,,,2
LESIZE,7,,,2
lesize,17,,,2
kscon,24,0.0005,0,8,1.0
LESIZE,33,,,2
LESIZE,47,,,2
LESIZE,26,,,2
kscon,16,0.0005,0,8,1.0
LESIZE,15,,,2
kscon,29,0.0005,0,8,1.0
LESIZE,46,,,2
AESIZE,18,0.005
AESIZE,20,0.005
AESIZE,5,0.05
AESIZE,16,0.5
allsel
aatt,1,1,1
amesh,3
amesh,4
amap,1,9,12,13,10
amap,2,9,12,17,14
amesh,10
amesh,12
amap,9,21,25,26,22
amap,11,28,31,26,22
amesh,7
amesh,17
amap,6,7,18,47,38
amap,19,7,18,46,45
amesh,13
amesh,15
amap,8,27,32,48,33
amap,14,27,32,49,23
amesh,18
amesh,20
amesh,5
amesh,16
allsel
!加載
finish
/solu
!dk,34,all
!dk,50,all
!dk,51,all
!dk,52,all
allsel
nsel,s,loc,y,0
sf,all,pres,siga
allsel
nsel,s,loc,y,width
sf,all,pres,siga
allsel
nsel,s,loc,x,0
sf,all,pres,ksiga
allsel
nsel,s,loc,x,length
sf,all,pres,ksiga
allsel
lsel,s,,,61,62,1
lsel,a,,,64,65,1
nsll,r,1
sf,all,pres,p
allsel
nsel,s,,,14
lie1num=14
cm,lie1,node
allsel
cint,new,1
cint,type,sifs
cint,ctnc,lie1
cint,symm,off
cint,ncon,5
cint,normal,11,2
cint,list
allsel,all
nsel,s,,,277
lie2num=277
cm,lie2,node
allsel
cint,new,2
cint,type,sifs
cint,ctnc,lie2
cint,symm,off
cint,ncon,5
cint,normal,11,2
cint,list
allsel,all
antype,static
solve

技術鄰問答 


















