!半径为2m、球壳厚度为1.0e-4m的半球形膜，考虑内压为300pa结构的应力分布。该结构应该是有理论解为3mpa。如何分析设置才能使计算收敛，请高手指教，下面是我的命令流
fini
/cle
big_r =2
shell_t =0.1e-3
elastic =1.0e9 ！ETFE膜材
posion =0.4
prestr =3.0e8 !这里是我定义的预应力，不知道是否合适，可以修改
delt_t =-prestr/elastic
/prep7
k,1, 0,0,0
k,2, 0, big_r, 0
k,3, big_r, 0, 0
larc,2,3,1,big_r,
arotat,1, , , , , , 1, 2,90, ,
et,1,shell41
keyopt,1,1,2
mp,ex, 1,elastic
mp,prxy,1,posion
mp,alpx,1,1
tref,0
r,1,shell_t
lesize,all,,,20
mshape,1,2d
mshkey,1
amesh,1
arsym,x,all, , , ,0,0
arsym,z,all, , , ,0,0
alls
nummrg,all
numcmp,all
lsel,s,loc,y,0
nsll,s,1
d,all,all
alls
bfa,all,temp,delt_t
alls
eplot
/solu
antype,static
nlgeom,on
sstif,on
nsubst,10
cnvtol,f,,0.1
lnsrch,on
outres,all,all
solve
fini
/eof
!!!加载求解
/solu
alls
lsel,s,loc,y,0
nsll,u,1
ddele,all,all
alls
sf,all,pres,300
antype,static
nlgeom,on
sstif,on
nsubst,20
cnvtol,f,,0.01
lnsrch,on
outres,all,all
alls
solve
/post1
plnsol,s,eqv
prrsol
|
