ansys 导入初始应力文件 位移清零后隧道开挖,控制应力释放不收敛
采用的地层参数,在位移未清零情况下,用应力释放α进行隧道开挖不收敛,但是在位移未清零情况下,能够收敛!情况为何?是否我导入初始应力文件有什么问题?
采用地层参数:r=2100,E=6E+007,v=0.35,c=20000,φ=35
应力释放系数取α=0.2
/solu
nsubst,10 !指定荷载子步
nlgeom,on !指定大变形
nropt,full !指定newton-lapnace迭代
pred,on !启用非线性求解预测器
lnsrch,on !启用线性搜索
iswrite,on
time,1
allsel
solve
save,model4,db
finish
!导入模型model0
/solu
antype,,new
allsel
!!应力场的施加
isfile,read,DCM,ist,,2 !DCM是文件名,2代表节点应力
nsubst,10 !指定荷载子步
nlgeom,on !指定大变形
nropt,full !指定newton-lapnace迭代
pred,on !启用非线性求解预测器
lnsrch,on !启用线性搜索
time,1
allsel
solve
save,model5,db
finish
/post1
!调出节点力
/solu
antype,,rest !设定分析类型为重启动分析,重启动点默认为上次求解结束时
time,2
allsel
nlgeom,on !指定大变形
allsel
csys,1 !将当前坐标变为前面自定义的极座标
asel,s,loc,x,0,2.7
esla,s
ekill,all !杀死开挖土体单元,包括核心图、管片层、注浆层
allsel
parres,,003,PARM
csys,1
asel,s,loc,x,0,2.7
esla,s
nsle,s
nsel,r,ext
*do,i,1,n
f,n_nf(i,1),fx,0.80*n_nf(i,2)
f,n_nf(i,1),fy,0.80*n_nf(i,3)
*enddo
!cnvtol,f,,0.02,2,1
allsel !选择所有元素
solve !求解