做两次初始应力场，用iswrite和isfile遇到的问题，求指点

季冰乙

这里用个简单的模型来反应问题，如图1：隧道开挖前，做了一次初始地应力，用iswrite和isfile，很好的到达目标了。想再做一次初始地应力，把隧道开挖后形成的位移场也给化成零，用同样的办法总会出问题，如图2。
求高手指点，觉得第一次应该没问题的，如果正确大家可以分享下。关于第二次 位移场请零，如果有更好的办法也请指教。还要个附带的问题，有见过说isfile的时候 ，所在的步 time 号必须是1，如果是这种做2次的，time怎么调，重启么？
坐等回复中，谢谢指导


/PREP7      
ET,1,SOLID45
et,2,mesh200
keyopt,2,1,6
mp,ex,1,7.65e6          !最外圈土
mp,prxy,1,0.35
mp,dens,1,2030
mp,ex,2,7.65e6         !开挖的
mp,prxy,2,0.35
mp,dens,2,2030.01
blc4,,,10,10,10
wpoffs,5,5
cyl4,,,,,2,360,10
vovlap,all
allsel
nummrg,all   
numcmp,all
wprota,,,90
vsbw,all
wprota,,90
vsbw,all
allsel
nummrg,all   
numcmp,all
allsel
mshape,1,3d
mshkey,,
lesize,all,0.6
vsel,s,,,1,4
mat,2
vmesh,all
vsel,s,,,5,8
mat,1
vmesh,all
save,qwe01,db
/sol
ANTYPE,0
!* the zero load step only for get initial stress file
lsclear,all
TIME,0  
AUTOTS,-1   
NSUBST,1, , ,1  
KBC,0
nropt,full
nsel,s,loc,x,0,0
nsel,a,loc,x,10,10
d,all,ux
nsel,s,loc,y,0,0
d,all,uy
nsel,s,loc,z,0,0
nsel,a,loc,z,10,10
d,all,uz
!/psf,pres,norm,2               !显示压力!(norm表示正交,2以箭头形式显示)
allsel
acel,0,10,0
allsel
iswrite,on
solve
finish
save,qwe11,db

/solu   
ANTYPE,0
!* the first load step
lsclear,all
TIME,1  
AUTOTS,-1   
NSUBST,1, , ,1  
KBC,0
nropt,full
allsel
nsel,s,loc,x,0,0
nsel,a,loc,x,10,10
d,all,ux
nsel,s,loc,y,0,0
d,all,uy
nsel,s,loc,z,0,0
nsel,a,loc,z,10,10
d,all,uz
allsel
acel,0,10,0
allsel
esel,s,live
isfile,read,qwe,ist,,1
iswrite,off
solve
save,qwe12,db


!退出重读，为了TIME 1?
/solu
lsclear,all
TIME,0
AUTOTS,-1   
NSUBST,1, , ,1  
KBC,0
nropt,full
allsel
esel,s,mat,,2
ekill,all

allsel
nsel,s,loc,x,0,0
nsel,a,loc,x,10,10
d,all,ux
nsel,s,loc,y,0,0
d,all,uy
nsel,s,loc,z,0,0
nsel,a,loc,z,10,10
d,all,uz
allsel
acel,0,10,0
allsel
iswrite,on
solve
save,qwe13,db      


/solu   
ANTYPE,0
!* the first load step
lsclear,all
TIME,1  
AUTOTS,-1   
NSUBST,1, , ,1  
KBC,0
nropt,full

allsel
nsel,s,loc,x,0,0
nsel,a,loc,x,10,10
d,all,ux
nsel,s,loc,y,0,0
d,all,uy
nsel,s,loc,z,0,0
nsel,a,loc,z,10,10
d,all,uz
allsel
esel,s,mat,,2
ekill,all
allsel
acel,0,10,0

allsel
esel,s,live
isfile,read,qwe,ist,,1
iswrite,off
solve
save,qwe14,db

季冰乙

坐等指导中

季冰乙

没人气，还是说的不清楚呀，来个人 交流下呀:)

terryy11

季冰乙 发表于 2011-10-9 20:11
没人气，还是说的不清楚呀，来个人 交流下呀

最后两个求解模块提出来编到一起，改下单元生死那，最后结果就对了，最大位移.175E-7=0.

季冰乙

terryy11 发表于 2011-10-9 21:42
最后两个求解模块提出来编到一起，改下单元生死那，最后结果就对了，最大位移.175E-7=0.

仔细研究下

cheng1002046

能介绍一下，第一次的隧道开挖模拟吗？具体是个什么流程？刚开始学ansys，谢谢！

英雄气不喘

能帮我看看出了什么问题吗？

英雄气不喘

我擦，第二次单元忘删了。。。