【讨论】请问如何改进？

wangliang

这是个隧道开挖的例子，在别人做的基础上改的，其中有许多不足之处
请各位朋友多指点！
/com  
fini  
/cle  
*set,x1,-20  
*set,y1,-30  
*set,w,40
*set,h,50
*set,th,0.5  
  
/prep7  
  
k,,0,0 $k,,0,3.85  
k,,0.88,5.5 $k,,2.45,6.15  
k,,4.02,5.5 $k,,4.9,3.85 $k,,4.9,0  
  
larc,1,2,6,8.13 $larc,2,3,6,3.21  
larc,3,4,6,2.22 $larc,4,5,2,2.22  
larc,5,6,2,3.21 $larc,6,7,2,8.13 $larc,7,1,4,6  
line,2,6
  
a,1,2,6,7  
a,2,3,4,5,6
blc4,x1,y1,w,h  
aovl,3,2,1  
  
lsel,s,line,,1,6,
lcomb,all,, 　　　　　　　　　　
lsel,all,  
  
et,1,beam3  
et,2,plane42  
keyopt,2,3,2  
r,1,th,th*th*th/12,th, ,  
  
mp,ex,1,2.6e7  
mp,prxy,1,0.2  
mp,dens,1,25  
  
mp,ex,2,1.3e6  
mp,prxy,2,0.32  
tb,dp,2  
tbdata,1,20,30,  
mp,dens,2,22  
  
mp,ex,3,1.3e6  
mp,prxy,3,0.32  
tb,dp,3  
tbdata,1,20,30,  
mp,dens,3,22  
  
mp,ex,4,1.3e6  
mp,prxy,4,0.32  
tb,dp,4  
tbdata,1,20,30,  
mp,dens,4,22
  
lsel,s,line,,1
lesize,all,,,20,  
latt,1,1,1  
lmesh,all  
alls  
  
lsel,s,line,,6,8
lesize,all,,,20,  
latt,1,1,1  
lmesh,all  
alls  
  
lsel,s,line,,2,5
lesize,all,,,10,  
latt,1,1,1  
lmesh,all  
alls  
  
lsel,s,line,,9,12  
lesize,all,,,40
alls  
  
asel,s,area,,4  
aatt,4,,2  
amesh,all  
alls  
  
asel,s,area,,2  
aatt,3,,2  
amesh,all  
alls  
  
asel,s,area,,1  
aatt,2,,2  
amesh,all  
alls  
  
nsel,s,loc,x,x1  
nsel,a,loc,x,x1+w  
d,all,ux,0  
alls                        !select everthing
  
nsel,s,loc,y,y1  
d,all,uy,0  
alls  
  
acel,,9.8
  
fini
                               
  
/solu  
antype,static  
  
autots,on !使用自动时间步  
pred,on !打开时间步长预测器  
lnsrch,on !打开线性搜索  
nlgeom,on !打开大位移效果  
nropt,full !设定牛顿-拉普森选项  
cnvtol,f,,0.05,0  
solve
  
esel,s,mat,,2
ekill,all  
esel,all
time,3
nsubst,5
solve
  
time,100
solve
  
esel,s,mat,,3
ekill,all  
time,105
nsubst,5
solve
  
第一，把time t 中的时间t当成是真实的时间以及中间的时间间隔处理不知正确否？
第二，是否就必须要考虑流变，输入粘塑性指数，洞顶的位移才随时间变化？
第三，杀死的单元是否需要继续施加约束？