裂缝计算没有显示计算完成是为什么?

浏览:1044 回答:1

549799863016195029.png

如图所示的结构进行有限元模拟计算,按照下面的命令流进行计算,结果点击计算后最后未出现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



邀请回答 我来回答

全部回答

(1)
默认 最新
技术邻问答
@君莫
2017年10月23日
评论 点赞

没解决?试试专家一对一服务

换一批