不透水層まではあまり深くないが集水埋渠底が不透水層まで達していない場合、CALCULATION OF WATER VOLUME FLOWING IN DRAIN. CASE:RIVER-BED WATER,CIRCULAR DRAIN.

Equation
Q=\frac{kL(H^{2}-h^{2})}{R(\frac{h}{t+0.5r})^{0.5}(\frac{h}{2h-t})^{0.25}}

000185

Scilab
185sci01185sci02185sci03185sci04

 -->clear;funcprot(0);
-->deff('[Q]=f1(R,r)','Q=k/100*L*(H^2-h^2)/(R*(h/(t+0.5*r))^0.5*(h/(2*h-t))^0.25)');
-->//地下水揚水量:Q(m^3/sec)
-->//patternf1
-->k=10^-1;//透水係数(cm/sec)
-->H=5;//(m)
-->h=2.5;//(m)
-->L=20;//(m)
-->t=0.5;//(m)
-->R=linspace(5,20);//(m)
-->r=linspace(0.5,3);//(m)
-->xset("window",1);clf();xgrid();
-->condition='k='+string(k)+',H='+string(H)+',h='+string(h)+',L='+string(L)+',t='+string(t);
-->graphtitle='地下水揚水量Q(m^3/sec).' + condition;
-->xtitle(graphtitle,'R(m)','r(m)','Q(m^3/s)',boxed=1);
-->color(240,255,240);
-->idcolor=color(240,255,240);
-->fplot3d(R,r,f1,flag=[idcolor , 2 , 4]);
-->xset("window",2);clf();xgrid();
-->xtitle(graphtitle,'R(m)','r(m)','Q(m^3/s)',boxed=1);
-->contour(R,r,f1,30,flag=[2,2,4]);
-->f1(5,3)
 ans  =
 
    0.0777006  
-->f1(20,0.5)
 ans  =
 
    0.0118954  
-->clear;funcprot(0); 
-->deff('[Q]=f2(H,h)','Q=k/100*L*(H^2-h^2)/(R*(h/(t+0.5*r))^0.5*(h/(2*h-t))^0.25)');
-->//地下水揚水量:Q(m^3/sec)
-->//patternf2
-->k=10^-1;//透水係数(cm/sec)
-->R=25;//(m)
-->r=1;//(m)
-->L=20;//(m)
-->t=0.5;//(m)
-->H=linspace(5,15);//(m)
-->h=linspace(1,3);//(m)
-->xset("window",3);clf();xgrid();
-->condition='k='+string(k)+',R='+string(R)+',r='+string(r)+',L='+string(L)+',t='+string(t);
-->graphtitle='地下水揚水量Q(m^3/sec).' + condition;
-->xtitle(graphtitle,'H(m)','h(m)','Q(m^3/s)',boxed=1);
-->color(255,239,219);
-->idcolor=color(255,239,219);
-->fplot3d(H,h,f2,flag=[idcolor , 2 , 4]);
-->xset("window",4);clf();xgrid();
-->xtitle(graphtitle,'H(m)','h(m)','Q(m^3/s)',boxed=1);
-->contour(H,h,f2,30,flag=[2,2,4]);
-->f2(15,7)
 ans  =
 
    0.0627137  

SciNotes

clear;funcprot(0);

deff('[Q]=f1(R,r)','Q=k/100*L*(H^2-h^2)/(R*(h/(t+0.5*r))^0.5*(h/(2*h-t))^0.25)');
//地下水揚水量:Q(m^3/sec)

//patternf1
k=10^-1;//透水係数(cm/sec)
H=5;//(m)
h=2.5;//(m)
L=20;//(m)
t=0.5;//(m)
R=linspace(5,20);//(m)
r=linspace(0.5,3);//(m)
  
xset("window",1);clf();xgrid();

condition='k='+string(k)+',H='+string(H)+',h='+string(h)+',L='+string(L)+',t='+string(t);
graphtitle='地下水揚水量Q(m^3/sec).' + condition;
xtitle(graphtitle,'R(m)','r(m)','Q(m^3/s)',boxed=1);

color(240,255,240);
idcolor=color(240,255,240);
fplot3d(R,r,f1,flag=[idcolor , 2 , 4]);
  
xset("window",2);clf();xgrid();

xtitle(graphtitle,'R(m)','r(m)','Q(m^3/s)',boxed=1);
contour(R,r,f1,30,flag=[2,2,4]);
  
f1(5,3)
  
f1(20,0.5)
    
clear;funcprot(0); 

deff('[Q]=f2(H,h)','Q=k/100*L*(H^2-h^2)/(R*(h/(t+0.5*r))^0.5*(h/(2*h-t))^0.25)');
//地下水揚水量:Q(m^3/sec)

//patternf2
k=10^-1;//透水係数(cm/sec)
R=25;//(m)
r=1;//(m)
L=20;//(m)
t=0.5;//(m)
H=linspace(5,15);//(m)
h=linspace(1,3);//(m)
  
xset("window",3);clf();xgrid();

condition='k='+string(k)+',R='+string(R)+',r='+string(r)+',L='+string(L)+',t='+string(t);
graphtitle='地下水揚水量Q(m^3/sec).' + condition;
xtitle(graphtitle,'H(m)','h(m)','Q(m^3/s)',boxed=1);

color(255,239,219);
idcolor=color(255,239,219);
fplot3d(H,h,f2,flag=[idcolor , 2 , 4]);
  
xset("window",4);clf();xgrid();

xtitle(graphtitle,'H(m)','h(m)','Q(m^3/s)',boxed=1);
contour(H,h,f2,30,flag=[2,2,4]);
  
f2(15,7)

参考文献
川島普、篠原紀、西川泰治(1999).『新訂衛生工学』.明現社.184pp.

アプリケーション
Scilab http://www.scilab.org/