各位老大,你们看了,以下是我修改一新的命令流,但是无法计算成功,请示下:::
/prep7
et,1,shell93
!keyopt,1,2,1
mp,ex,1,2.55e10
mp,prxy,1,0.167
mp,dense,7850
r,1,0.04
k,1
k,2,,1000,-3
k,3,,2000,-13
k,4,,3000,-30
k,5,,4000,-53
k,6,,5000,-92
k,7,,6000,-120
bsplin,all
k,8,200,,19
k,9,400,,76
k,10,600,,171
k,11,800,,304
k,12,990,,465
bsplin,1,8,9,10,11,12
allsel,all
lplot
/pnum,line,1
/replot
adrag,1,,,,,,2
numcmp,all
nummrg,all
lsel,s,,,1,3,2
lesize,all,,,80
lsel,s,,,2,4,2
lesize,all,,,30
mshkey,1
amesh,all
fini
/solu
antype,static
nsel,s,loc,y
dsym,symm,y
nsel,s,loc,x
dsym,symm,x
nsel,s,loc,y,6000
d,all,ux,,,,,uz
allsel,all
sfa,1,1,pres,1.5e-3
nlgeom,on
outpres,,1
nsubstp,30
solve |