地下水揚水量 井底が不透水層に達しない場合 井底だけから流入する場合CALCULATION OF WATER VOLUME FLOWING IN WELL. CASE:UNCONFINED GROUNDWATER,H>0M,FROM BOTTOM.

Equation
Q=4kr(H-h)

000178

Scilab

000178sc01000178sc02

 -->clear;funcprot(0);
-->deff('[Q]=f1(H,h)','Q=4*k/100*r*(H-h)');
-->//地下水揚水量:Q(m^3/sec)
-->//patternf1
-->k=10^-1;//透水係数(cm/sec)
-->r=1;//(m)
-->H=linspace(5,20);//(m)
-->h=linspace(0.5,3);//(m)
-->xset("window",1);clf();xgrid();
-->condition='k='+string(k)+',r='+string(r);
-->graphtitle='地下水揚水量Q(m^3/sec).' + condition;
-->xtitle(graphtitle,'H(m)','h(m)','Q(m^3/s)',boxed=1);
-->color(240,255,240);
-->idcolor=color(240,255,240);
-->fplot3d(H,h,f1,flag=[idcolor , 2 , 4]);
-->currentaxes=get("current_axes");
-->currentaxes.font_size=6;
-->currentaxes.title.font_size=6;
-->currentaxes.x_label.font_size=6;
-->currentaxes.y_label.font_size=6;
-->currentaxes.z_label.font_size=6;
-->currentaxes.tight_limits="on";
-->currentaxes.cube_scaling="on";
-->currentaxes.auto_clear="on";
-->currentaxes.auto_scale="on";
-->xset("window",2);clf();xgrid();
-->xtitle(graphtitle,'H(m)','h(m)','Q(m^3/s)',boxed=1);
-->contour(H,h,f1,30,flag=[2,2,4]);
-->currentaxes=get("current_axes");
-->currentaxes.font_size=6;
-->currentaxes.title.font_size=6;
-->currentaxes.x_label.font_size=6;
-->currentaxes.y_label.font_size=6;
-->currentaxes.z_label.font_size=6;
-->currentaxes.tight_limits="on";
-->currentaxes.cube_scaling="on";
-->currentaxes.auto_clear="on";
-->currentaxes.auto_scale="on";
-->f1(5,3)
 ans  =
 
    0.008  
-->f1(20,0.5)
 ans  =
 
    0.078  

SciNotes

clear;funcprot(0);

deff('[Q]=f1(H,h)','Q=4*k/100*r*(H-h)');
//地下水揚水量:Q(m^3/sec)

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

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

color(240,255,240);
idcolor=color(240,255,240);
fplot3d(H,h,f1,flag=[idcolor , 2 , 4]);

currentaxes=get("current_axes");
currentaxes.font_size=6;
currentaxes.title.font_size=6;
currentaxes.x_label.font_size=6;
currentaxes.y_label.font_size=6;
currentaxes.z_label.font_size=6;
currentaxes.tight_limits="on";
currentaxes.cube_scaling="on";
currentaxes.auto_clear="on";
currentaxes.auto_scale="on";
  
xset("window",2);clf();xgrid();

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

currentaxes=get("current_axes");
currentaxes.font_size=6;
currentaxes.title.font_size=6;
currentaxes.x_label.font_size=6;
currentaxes.y_label.font_size=6;
currentaxes.z_label.font_size=6;
currentaxes.tight_limits="on";
currentaxes.cube_scaling="on";
currentaxes.auto_clear="on";
currentaxes.auto_scale="on";
  
f1(5,3)
  
f1(20,0.5)

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

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