/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