隧道开挖，isfile命令导入初始应力文件，后开挖岩土体，不收敛

tjwxx

各位高手，我在做隧道开挖时，分为两步：
（1）先加边界条件和重力，计算重力作用下的应力和位移，此时iswrite,on，将重力产生的初始应力保存为ist文件。
（2）退出ansys，重新启动，加相同的边界条件和重力，读入初始应力文件，同时杀死开挖土体单元，约束住非活动单元的节点，再进行求解。
第二步总是不能收敛，请教各位是什么原因呢？附上命令流：
/COM, Structural                             
/TITLE, Tunnel Construct Modeling Analysis      
/FILNAM,support,1                          
!（2）建立模型
！进入前处理器
/PREP7
*afun,deg
!定义单元类型         
ET,2,PLANE42                                
KEYOPT,2,3,2           
!定义材料属性
!围岩材料
MP,EX,2,1E9  
MP,PRXY,2,0.45   
MP,DENS,2,1700
TB,DP,2
TBDATA,1,0.05E6,20
!挖去土体材料
MP,EX,3,1E9  
MP,PRXY,3,0.45   
MP,DENS,3,1700
TB,DP,3
TBDATA,1,0.05E6,20
k,1,-5.7,0,0
k,2,0,5.7,0
k,3,5.7,0,0
larc,1,3,2,5.7
k,4,5.6582,-0.8269,0
k,5,5.5332,-1.6454,0
larc,3,5,4
k,6,5.3048,-2.1047,0
k,7,4.8788,-2.3904,0
larc,5,7,6
k,8,0,-3.2060,0
k,9,-4.8788,-2.3904,0
larc,7,9,8
k,10,-5.3048,-2.1047,0
k,11,-5.5332,-1.6454,0
larc,9,11,10
k,12,-5.6582,-0.8269,0
larc,11,1,12


k,60,-100,-100
k,61,-10,-100
k,62,10,-100
k,63,100,-100
k,64,100,-10
k,65,10,-10
k,66,-10,-10
k,67,-100,-10
k,68,-100,10
k,69,-10,10
k,70,10,10
k,71,100,10
k,72,100,100
k,73,10,100
k,74,-10,100
k,75,-100,100
k,76,-100,0,0
k,77,100,0,0
k,78,0,100,0
k,79,0,-100,0
l,60,61                                          
l,61,62
l,62,63
l,63,64
l,64,65
l,65,66
l,66,67
l,67,68
l,68,69
l,69,70
l,70,71
l,71,72
l,72,73
l,73,74
l,74,75
l,75,68
l,74,69
l,73,70
l,69,66
l,70,65
l,71,64
l,67,60
l,66,61
l,65,62
l,76,77
l,78,79

a,60,63,72,75                                 
asbl,all,all  
adele,1                                    
NUMCMP,AREA  
save,Tunel-geom.db

lcomb,5,6
lcomb,5,64
lcomb,2,3
lcomb,2,62
lsel,s,,,37,41,1
lsel,a,,,2,5,3
lsel,a,,,63
lesize,all,,,16
mat,3
type,2
amap,1,14,13,1
amap,2,14,13,3
amap,15,14,3,19
amap,16,14,19,1

lsel,inve
lesize,all,,,8,,,,,1
mat,2
type,2
amap,17,13,1,17,16
amap,18,13,3,17,15
amap,19,3,19,15,18
amap,20,1,19,16,18

mshkey,0
mshape,0
asel,s,,,3,14,1
amesh,all
save,Tunel-grid,db

csys,0
nsel,s,loc,x,-100                              
nsel,a,loc,x,100                                 
d,all,ux                                       
allsel
nsel,s,loc,y,-100                              
d,all,uy  
allsel,all,all
acel,,9.8
                              
/solu
antype,static
solcontrol,on
autots,on
cnvtol,f,0.05,,,
cnvtol,u,0.05,,,                                 
nsubst,200,                              
pred,on                                          
nropt,full                                       
nlgeom,on                                    
lnsrch,on                                      
outres,all,all                                                                     
time,1  
iswrite,on                                                            
allsel   
solve

退出后重启，
前面命令不变，从time，1后命令为
isfile,read,support,ist,,1
allsel                                   
esel,s,mat,,3                                      
ekill,all                                                                             
esel,all
esel,s,live                                         
nsle,s
nsel,inve
d,all,all,0                                         
allsel
solve

tjwxx

如果重启后只是导入初应力文件不进行开挖，此时可以看到应力为重力作用下的应力，位移基本为0，说明isfile命令可以清除重力产生的初始位移。
1# tjwxx

tjwxx

另外如果将iswrite 和isfile命令改为新的inistate命令，发现其不能去除重力产生的初始位移。
命令如下：建模命令和之前相同
/solu
antype,static                                    
nsubst,100                                 
pred,on                                          
nropt,full                                       
nlgeom,on                                    
lnsrch,on                                      
outres,all,all                                                                     
time,1  
inistate,write,1,                                                                    
allsel   
solve
finish

/solu
inistate,read,support,ist
allsel                                   
esel,s,mat,,3                                      
ekill,all                                                                             
esel,all
esel,s,live                                         
nsle,s
nsel,inve
d,all,all,0                                         
allsel
solve
以上命令在ansys13.0下运行，拱顶和拱底均下沉，拱顶下沉14cm左右，感觉没有去除重力产生的初始位移。

tjwxx

有人遇到过这种情况吗？

c7225237

一个人也没回复吗？

cumt410

单独查看轴向的位移，最近我也在做这个