裂缝计算没有显示计算完成是为什么?
如图所示的结构进行有限元模拟计算,按照下面的命令流进行计算,结果点击计算后最后未出现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