探讨隧道开挖过程模拟方法
tzz1997
tzz1997 Lv.2
2010年05月08日 22:45:27
来自于ANSYS
只看楼主

在网上看到有人编写这样一个例子隧道上下台阶开挖支护的模拟过程如下所示:!施加约束allselnsel,s,loc,x,x_zsnsel,a,loc,x,x_yxd,all,uxnsel,s,loc,y,y_yxd,all,uysave,step0,db!指定求解选项/solunlgeom,onnropt,fullpred,onlnsrch,onnsubst,10,50,10!荷载步1,自重应力场allsel

在网上看到有人编写这样一个例子隧道上下台阶开挖支护的模拟过程如下所示:
!施加约束
allsel
nsel,s,loc,x,x_zs
nsel,a,loc,x,x_yx
d,all,ux
nsel,s,loc,y,y_yx
d,all,uy

save,step0,db

!指定求解选项
/solu
nlgeom,on
nropt,full
pred,on
lnsrch,on
nsubst,10,50,10

!荷载步1,自重应力场
allsel
esel,s,type,,1,3,2
ekill,all

allsel
acel,,9.8
allsel
solve
save,step1,db
esel,s,live
eplot

!荷载步2,上台阶开挖,应力释放
!提取节点力
/post1
esel,s,mat,,3
nsle,s
nsel,r,ext
*get,n,node,,count
*dim,n_nf,array,n,3
*do,i,1,n
*get,n_nf(i,1),node,,num,min
nsel,u,,,n_nf(i,1)
*enddo
!*status,n_nf
*do,i,1,n
nsel,s,,,n_nf(i,1)
fsum
*get,n_nf(i,2),fsum,,item,fx
*get,n_nf(i,3),fsum,,item,fy
*enddo
!*status,n_nf
parsav,all,002,PARM

/solu
antype,,rest
!杀死上台阶土体单元
allsel
esel,s,mat,,3
ekill,all
!施加等效释放节点力
parres,,002,PARM
esel,s,mat,,3
nsle,s
nsel,r,ext
*do,i,1,n
f,n_nf(i,1),fx,0.65*n_nf(i,2)
f,n_nf(i,1),fy,0.65*n_nf(i,3)
*enddo

!约束死单元自由度
allsel
esel,s,live
nsle,s
nsel,inve
d,all,all

allsel
solve
save,step2,db
esel,s,live
eplot

!荷载步3
!删除节点力,激活上台阶支护和锚杆
allsel
esel,s,mat,,3
nsle,s
nsel,r,ext
fdele,all,all

esel,s,mat,,1,6,5
ealive,all

allsel
solve
save,step3,db

!荷载步4,下台阶开挖,应力释放
!提取节点力
/post1
esel,s,mat,,4
nsle,s
nsel,r,ext
*get,n,node,,count
*dim,n_nf,array,n,3
*do,i,1,n
*get,n_nf(i,1),node,,num,min
nsel,u,,,n_nf(i,1)
*enddo
!*status,n_nf
*do,i,1,n
nsel,s,,,n_nf(i,1)
fsum
*get,n_nf(i,2),fsum,,item,fx
*get,n_nf(i,3),fsum,,item,fy
*enddo
!*status,n_nf
parsav,all,002,PARM

/solu
antype,,rest
!杀死下台阶土体单元
allsel
ealive,all
esel,s,mat,,2,3,1
esel,a,mat,,4
ekill,all
!施加等效释放节点力
parres,,002,PARM
esel,s,mat,,4
nsle,s
nsel,r,ext
*do,i,1,n
f,n_nf(i,1),fx,0.65*n_nf(i,2)
f,n_nf(i,1),fy,0.65*n_nf(i,3)
*enddo

!约束死单元自由度
allsel
esel,s,live
nsle,s
nsel,inve
d,all,all

allsel
solve
save,step4,db
esel,s,live
eplot

!荷载步5
!删除节点力,激活下台阶支护
allsel
esel,s,mat,,4
nsle,s
nsel,r,ext
fdele,all,all

esel,s,mat,,2
ealive,all

allsel
solve
save,step5,db


请各位高人探讨一下,这程序才用的是多载荷步法,还是其他什么地应力自动释放方法,怎么看不明白,而且一般开挖支护是在一次求解中同步进行的,而这个程序却是分两步运行的,有什么区别吗?提取的节点力也没有采取反向啊,

[ 本帖最后由 tzz1997 于 2010-5-8 22:47 编辑 ]
免费打赏
tzz1997
2010年05月09日 16:34:21
2楼
各位大侠,能否高抬贵手,指点迷津
回复
tzz1997
2010年05月14日 10:38:34
3楼
我等到花儿也谢了,:Q :Q :Q :Q :Q :Q
回复
tzz1997
2010年05月27日 07:56:23
4楼
五十多个人看了,没有一个人说两句的,太失望了
回复
tzz1997
2010年05月27日 08:18:08
5楼
这是模型材料详图,开挖掉部分,为材料3,4,请求高人指点,将感激不尽
回复

相关推荐

APP内打开