求隧道断面流量计算fish

yxf4171

如题：
       前辈们谁做过隧道断面流量计算方面的东西？论坛上说有fish，我找遍了都没找到，谁有给发一个吧！

yxf4171

独自等待中……

yxf4171

关于gp_flow()计算得到的流量不太明白，我的理解是节点每秒不平衡流量，也就是节点每秒的流量，不知道这样理解对不对？
请路过的高手指点指点……

xiaoshi80

title
   Calculation for Water Gush of Tunnel
def setup
h=166.0            ;thickness of water bearing layer
h1=101.02          ;height from cicle center to impermeable layer
r0=3.52            ;cicle radius
l=500.0            ;length of tunnel
c_perm=0.109e-13    ermeability (m/d)/(Pa/m)
pp0=h*1e4          ;pore pressure
r00=-1*r0
r1=r0-0.001
r2=r0+0.001
end
setup

config fluid
gen zone radcyl p0 0 0 h1 p1 add 150 0 0 p2 add 0 1 0 p3 0 0 h &
                dim r0 r0 r0 r0 size 3 1 6 30 rat 1 1 1 1.1
gen zone radcyl p0 0 1 h1 p1 add 150 0 0 p2 add 0 -1 0 p3 0 1 0 &
                dim r0 r0 r0 r0 size 3 1 6 30 rat 1 1 1 1.1
                ;p4 add 150 -1 0 p5 0 0 0 p6 150 1 0 p7 150 0 0 &
                ;p8 add r0 0 0 p9 add 0 0 r00 p10 add r0 -1 0 p11 add 0 -1 r00 &
                 
gen zone reflect dd 270 dip 90 ori 0 0 0


model fl_iso
prop perm c_perm
set fl biot off
ini fmod 2e9
ini fdens 1000
ini ftens 0
ini sat 1

ini pp pp0 grad 0 0 -1e4   
fix pp 0 range z 165.99 166.01
fix pp 0 range cyl end1 0 0 h1 end2 0 1 h1 rad r1 &
           not cyl end1 0 0 h1 end2 0 1 h1 rad r2              
set grav 0,0,-10

set mech off
set fl on

def qflac
    qval = 0.0
    pnt = gp_head
    loop while pnt # null
         fval=sqrt((gp_xpos(pnt))^2+(gp_zpos(pnt)-h1)^2)-0.001  ;fval=sqrt(gp_xpos(pnt)^2+(gp_ypos(pnt)-h1)^2)-0.001
           if fval < r0 then
              qval =qval + gp_flow(pnt)
           end_if
           pnt = gp_next(pnt)
    end_loop
    qflac=qval*l
    qsol=c_perm*(h-h1+r0)*(0.676-0.06*c_perm*1e4)*1e4*l   ;empirical calculation
end

;qflac
;print qflac qsol

plot cont pp outline on

plot add flow

solve ratio 1e-2
qflac
print qflac qsol
save ysq.sav

xiaoshi80

参考一下这个吧

yxf4171

谢谢了 …… 很有帮助
我感觉这个FISH算出来的是每秒每延米的流量吧
我想算一天每延米的流量，直接乘以3600x24不对吧
是不是要用个循环，把每秒的流量叠加起来，呵呵这个还没实现

T-MACLeo

这个FISH算出来的是每秒每延米的流量，是算一天每延米的流量，可直接乘以3600x24，因为你计算的时候是达到稳定渗流的一个状态，没必要用循环把每秒的流量叠加起来你也算不过来！！！

T-MACLeo

应该不是每延米的，而是隧道整个长度的，因为计算的时候长度方向是300，理论值是乘了l的，而在数值解的记录里面是只管了竖向平面即x和z而没有管y方向的，说明这个方向也是全部记录了的，所以渗流量是整个隧道长度每秒的流量

T-MACLeo

不过我还是看不懂这个是什么：qsol=c_perm*(h-h1+r0)*(0.676-0.06*c_perm*1e4)*1e4*l，这个理论计算是怎么来的呢？求大虾指教，初始应该是根据达西定律Q=kiA来的，里面的c_perm*1e4应该是渗透系数k并且换算成了m/s为单位的,而其它那些部分中(h-h1+r0)是隧道顶点距含水层底的距离，它和i还有A有什么关系呢？而(0.676-0.06*c_perm*1e4)又是什么呢？求指教啊。。。谢过了先

xu1793543284

xiaoshi80 发表于 2009-4-27 09:41
title
   Calculation for Water Gush of Tunnel
def setup

这个su隧道是一个圆形我想问一下，如果是五心圆该如何解决？半径不是唯一的，fish语言该怎么编写