计算涌水量
在试图计算一个隧道的涌水量时，出现了一个问题，先看程序吧rest 耦合.sav
ini xvel 0.0 yvel 0.0 zvel 0.0
ini xd 0.0 yd 0.0 zd 0.0
ini sta 0
;注浆加固
prop bulk 3.34e8 shear 2e8 fric 30 coh 7e4 ten 50e4density 2500 range group zj
prop perm 1e-11 poros 0.01 biot_c 1 range group zj
mo null range group kw y 0 20
mo fl_null range group kw y 0 20
fix pp 0 range group kw y 0 20
step 50
sel shell id 1 range gro erch y 0.1 18
sel shell property isotropic 10e9 0.2 thickness 0.25 density 2500
setfluid on
def LiuLiang
sage=120
loop while sage < 3601 ;一个小时
p_gp=gp_head
command
solve age sage
endcommand
qval=0.0
loop while p_gp # null
ydis=gp_ypos(p_gp)
xdis=gp_xpos(p_gp)
zdis=gp_zpos(p_gp)
if ydis<20 then
if ydis>18 then
if zdis>64 then
if zdis<77 then
if xdis>71.9 then
if xdis<72.1 then
qval=qval+gp_flow(p_gp)
end_if
end_if
end_if
end_if
end_if
end_if
end_loop
n=sage/60
sage=sage+60
SUM=-60*qval+SUM
xtable(1,n)=n
ytable(1,n)=-1*qval
xtable(2,n)=n
ytable(2,n)=SUM
end_loop
end
LiuLiang
plot add table 1 both
plot add flow v plane
plot set b w
plot show
print sum
save
初始化中围岩的渗透系数为4.6e-11孔隙率为0.4，上面加固区渗透系数（程序中红色字母）应该提高才对，即采用了1e-11，这时就会出现 zero volume tet in zone 5 的错误，当加固区渗透系数改为1e-9时，就没有了错误，但是这不合理啊，加固区应该是减小了渗透系数才对，否则就失去了作用。请教高手，这个该怎么解决？ 楼主解决了 吗？可以把你的命令流发一分给我吗？谢谢1793543284@qq.com 为什么没大神回复呢:( zhuqiankun2016 发表于 2018-7-14 20:33
为什么没大神回复呢
请问您最近关于渗流量监测有什么进展？可以探讨一下吗？
有没有用flac3d 5.0 做隧道涌水量的计算的 zhaoge 发表于 2018-11-22 15:11
有没有用flac3d 5.0 做隧道涌水量的计算的
我在做，你好，你做的怎么样了:lol
你好，想请教一个问题，这if里的范围是怎么定义的？是直接用的全局坐标吗？同时，我用你的命令流格式简单的计算了一下，为啥一直好像是程序运行当中，但就是计算不出来，软件卡住了似的，这问题解决了嘛？
页:
[1]