本人最近在用Ansys模拟隧道开挖，现在先做一个简单的实例，二维隧道模拟。
起初我先使用施加重力加速度方法模拟自重应力，发现最后计算地面出现较大的隆起，通过查询相关帖子发现这是通有的问题，于是改用读入初始应力场的方法模拟自重应力。但由于与以前没有用过isfile命令，现在算起来不收敛，计算部分命令如下
/sol !求解
allsel
nsel,s,loc,x,-30 !施加约束
nsel,a,loc,x,30
d,all,ux,0
nsel,s,loc,y,-29.12
d,all,uy,0
acel,0,10,0 !施加重力加速度
save
nropt,full !设置求解选项
nlgeom,1
allsel
esel,s,type,,2 !杀死衬砌单元
ekill,all
allsel
iswrite,on !读出初始应力场，保存为ceshi.ist文件
solve
finish
save
/sol
antype,,new !重新求解计算
allsel
acel,0,10,0 !施加重力加速度
nsel,s,loc,x,-30 !施加约束
nsel,a,loc,x,30
d,all,ux,0
nsel,s,loc,y,-29.12
d,all,uy,0
save
nropt,full !设置求解选项
nlgeom,1
esel,s,type,,2 !杀死衬砌单元
ekill,all
allsel
esel,s,live !选择生单元
isfile,read,ceshi,ist,,1 !读入初始应力场到应力点
allsel
solve
finish
save
/sol
antype,,rest,1,1,0 !重新求解计算
allsel
esel,s,type,,2
ekill,all
asel,s,area,,18,20,2 !杀死开挖土单元
esla,r
ekill,all
allsel
esel,s,type,,2 !激活衬砌单元
ealive,all
allsel
……
希望各位高手指点，谢谢！
