程序使用说明


建立程序输入文件与DOS程序放在同一个目录下
文件名必须为:
in.dat
数据的意义参看输出文件说明
运行程序会生成输出文件result.dat，如上

这两个程序下载地址：
http://water.changlai.cn/土堤渗流逸出高度计算程序.zip
或者从附件下载：

我是刚刚接触堤防渗流计算，很多地方还不是很明白
把计算经验共享，仅供参考
老资格的高工们看来这只是很简单的东西
有什么问题回帖告诉我
双层堤基的计算程序很快就编好了
以后再上传来共享

今天新算的：
附录E.5 背水坡渗流出口比降计算

项目 计算式 数值 单位 备注

1.计算参数
下游坡度m2 = 5 = 5 / 1:m2
堤底高程▽底 = 17 = 17 m
下游水位▽2 = 18 = 18 m
逸出高度h0 = 2.54 = 2.54 m
下游水深H2 = 18-17 = 1 m

2.不透水地基逸出坡降计算
下游无水时
渗出点A坡降Jo = 1/SQRT(1+5^2) = 0.196 /
堤坡与不透水面交点B坡降Jo = 1/5 = 0.2 / 1/m2
下游有水时
渗出段AB内的点
计算点高程 = 19 = 19 m AB段(18~19.54)内
y = 19-17 = 2 / AB段(2.54~1)内
渗出点A坡降Jo = 1/SQRT(1+5^2) = 0.196 /
n = 0.25×1/2.54 = 0.098 /
计算点坡降J = 0.196×((2.54-1)/(2-1))^0.098 = 0.204 / 在下游坡面渗出
浸没段BC内的点 适用于y/H2<0.95
计算点高程 = 17.7 = 17.7 m BC段(17~19.413)内
y = 17.7-17 = 0.7 / BC段(0~0.95)内
背坡与水平面夹角απ = atan(1/5) = 0.197 rad 11.287度
α = 0.197/π = 0.063 /
a0 = 1/(2×(0.063×(5+0.5)×SQRT(1+5^2))) = 0.283 /
b0 = 5/(2×(5+0.5)^2) = 0.083 /
1/2α-1 1/(2×0.063)-1 = 6.937 /
计算点坡降J = 0.283×(0.7/1)^6.937/(1+0.083×1/(2.54-1)) = 0.023 /

3.透水地基逸出坡降计算
下游无水时
渗出段AB内的点
计算点高程 = 19 = 19 m AB段(17~19.54)内
y = 19-17 = 2 / AB段(0~2.54)内
计算点坡降J = 1/SQRT(1+5^2)×(2.54/2)^0.25 = 0.208 /
地基段BC内的点
计算点与堤脚距离x = 3 = 3 m
计算点坡降J = 1/(2×SQRT(5))×SQRT(2.54/3) = 0.206 /
下游有水时
渗出段AB内的点
计算点高程 = 19 = 19 m AB段(18~19.54)内
y = 19-17 = 2 / AB段(1~2.54)内
计算点坡降J = 1/SQRT(1+5^2)×(2.54/2)^0.25 = 0.208 /
浸没坡段BC内的点
计算点高程 = 17.7 = 17.7 m BC段(0.283~2.823)内
y = 17.7-17 = 0.7 / BC段(0~2.54)内
背坡与水平面夹角απ = atan(1/5) = 0.197 rad 11.287度
α = 0.197/π = 0.063 /
α1 = 1/(1+0.063) = 0.941 / 1/(1+α)
l1 = 2.54×5 = 12.7 / h0×m2
l2 = 1×5 = 5 / H2×m2
r = 0.7×5 = 3.5 / y×m2
计算点坡降J = 0.234 = 0.234 / 在下游坡面渗出
浸没地基面段CD内的点
计算点距堤脚距离x = 2 = 2 m
计算点坡降J = 0.108 = 0.108 / 在下游坡面渗出

好家伙，把你的设计计算书都传上来了，不过这样也好，可以多得积分呢。

今天新算的：
附录E.6 水位降落时均质土堤的浸润线

项目 计算式 数值 单位 备注

1.计算参数
上游坡度m1 = 3 = 3 / ( 1:m1 )
堤底高程▽底 = 17 = 17 m
堤身土料渗透系数k = 1.00E-05 = 1.00E-05 m/s ＝ 0.001cm/s ＝ 0.864m/d
土体的孔隙率n = 39.70% = 39.70% /
降前水位▽1 = 27 = 27 m ( 0.0003125m/s )
降后水位▽0 = 20 = 20 m ( 0.00023148m/s )
水位降落所需时间T = 22 = 22 h ( = 79200 秒 = 0.9167 天 )
计算上游浸润线时间t = 15 = 15 h ( = 54000 秒 = 0.625 天 )

2.h0(t)计算
降前水深H1 = 27-17 = 10 m ( ▽1-▽底 )
降后水深H0 = 20-17 = 3 m ( ▽0-▽底 )
降距H = 10-3 = 7 m ( H1-H0 )
水位降落的速度V = 7/(22/24) = 7.636 m/d ( = 0.00009m/s= 0.00212m/h )
百分数α = 113.7×(0.0001175)^(0.607^(6+log(0.001)))/100 = 15.03% / ( 公式一 )
给水度μ = 0.1503×0.397 = 0.0597 /
k/(μV) = 0.864/(0.0597×7.636) = 1.895 / ( 1/10<k/(μV)<60,需要计算浸润线 )
0.31H(1/T)(k/(μV))1/4 = 0.31×7×1/79200×(1.895)^(1/4) = 0.0000321 /
h0(t) = 7-0.0000321t
h0(t) = 7-0.0000321×54000 = 5.2666 m ( t = 15h = 54000s )

3.试算法计算t时刻渗流量q(t)和上游坡出渗点高度he(t)
he(t)1 = 3+0.01 = 3.01 m ( H0+Step )
he(t)2 = 3+5.2666 = 8.267 m ( H0+h0(t) )
步长Step = 0.01 = 0.01 m
精度 = 0.001 = 0.001 m
L = 10×3 = 30 m ( H1×m1 )
试算得到he(t) = 5.34 = 5.34 m ( 程序计算 )
q/k1 = ((3+5.2666)^2-5.34^2)/(2×(30-3×5.34)) = 1.424 m
q/k2 = (5.34-3)/3×(1+ln(5.34/(5.34-3))) = 1.424 m
平均q/k = (1.424+1.424)/2 = 1.424 m
平均渗流量q = 1.424×0.00001 = 1.42E-05 m3/s/m

4.计算 t = 15h 时刻上游段浸润线
[H0+h0(t)]2 = (3+5.2666)^2 = 68.3367 m2
2q/k = 2×1.424 = 2.848 m
浸润线方程 y=SQRT(68.3367-2.848x)

今天新算的：

附录E.7 双层堤基渗流计算和覆盖的计算
附录E.7.1 背水侧无限长双层堤基渗流计算和覆盖的计算

项目 计算式 数值 单位 备注

1.计算参数
上游坡度m1 = 3 = 3 / ( 1:m1 )
下游坡度m2 = 5 = 5 / ( 1:m2 )
堤顶宽度B = 6 = 6 m
堤顶高程▽顶 = 27 = 27 m
堤底高程▽底 = 17 = 17 m
弱透水地基底高程▽弱 = 7 = 10 m
强透水地基底高程▽强 = 7 = 5 m
上游水位▽1 = 24.8 = 24.8 m
下游水位▽2 = 18 = 18 m
堤底高程▽底 = 17 = 17 m
迎水侧有限长度L = 100 = 100 m
堤身土料渗透系数k = 1.00E-05 = 1.00E-05 m/s ( ＝ 0.001cm/s ＝ 0.864m/d )
弱透水堤基渗透系数k1 = 5.00E-05 = 5.00E-05 m/s ( ＝ 0.005cm/s ＝ 4.32m/d
强透水堤基渗透系数k0 = 7.00E-03 = 7.00E-03 m/s ( ＝ 0.7cm/s ＝ 604.8m/d

2.越流系数A计算
k0/k1 = 0.007/0.00005 = 140 / ( 属于双层堤基 )
堤身高度Hd = 27-17 = 10 m ( ▽顶-▽底 )
上游水深H = 24.8-17 = 7.8 m ( ▽1-▽底 )
弱透水层厚度T1 = 17-10 = 7 m ( ▽底-▽弱 )
强透水层厚度T0 = 10-5 = 5 m ( ▽弱-▽强 )
堤底宽度b = 6+10×(3+5) = 86 m ( B+Hd×(m1+m2) )
越流系数A = SQRT(0.00005/(0.007×7×5)) = 0.014 m-1

3.弱透水层CD段承压水头计算(x以下游堤脚为原点，向下游为正)
计算点坐标x = 10 = 10 m ( x>0 )
h = 7.8×e^(-0.014×10)/(1+0.014×86-th(0.014×100)) = 5.142 m

4.弱透水层CD段承压水头计算(x’以下游堤脚为原点，向上游为正)
计算点坐标x’ = 10 = 10 m ( x’>0 )
h = 7.8×(1+0.014×10)/(1+0.014×86-th(0.014×100)) = 6.743 m

今天下午刚刚算的：

附录E.7 双层堤基渗流计算和覆盖的计算
附录E.7.2 透水地基上弱透水层等厚有限长，强透水层无限长双层堤基渗流计算

项目 计算式 数值 单位 备注

1.计算参数
上游坡度m1 = 3 = 3 / ( 1:m1 )
下游坡度m2 = 3 = 3 / ( 1:m2 )
堤顶宽度B = 6 = 6 m
堤顶高程▽顶 = 27 = 27 m
堤底高程▽底 = 17 = 17 m
弱透水地基底高程▽弱 = 7 = 10 m
强透水地基底高程▽强 = 7 = 5 m
上游水位▽1 = 24.8 = 24.8 m
下游水位▽2 = 18 = 18 m
迎水侧有限长度L1 = 30 = 30 m ( 弱透水地基 )
背水侧有限长度L2 = 100 = 100 m ( 弱透水地基 )
堤身土料渗透系数k = 1.00E-05 = 1.00E-05 m/s ( ＝ 0.001cm/s ＝ 0.864m/d )
弱透水堤基渗透系数k1 = 5.00E-05 = 5.00E-05 m/s ( ＝ 0.005cm/s ＝ 4.32m/d )
强透水堤基渗透系数k0 = 7.00E-03 = 7.00E-03 m/s ( ＝ 0.7cm/s ＝ 604.8m/d )

2.越流系数A计算
k0/k1 = 0.007/0.00005 = 140 / ( >100,属于双层堤基 )
弱透水层厚度T1 = 17-10 = 7 m ( ▽底-▽弱 )
强透水层厚度T0 = 10-5 = 5 m ( ▽弱-▽强 )
堤身高度Hd = 27-17 = 10 m ( ▽顶-▽底 )
上游水面距不透水地基H = 24.8-5 = 19.8 m ( ▽1-▽强 )
堤底距不透水地基H1 = 17-5 = 12 m ( ＝ ▽底-▽强 ＝ T1+T0 )
弱透水层底距不透水地基H0 = 10-5 = 5 m ( ▽弱-▽强 )
堤底宽度b = 6+10×(3+3) = 66 m ( B+Hd×(m1+m2) )
越流系数A = SQRT(0.00005/(0.007×7×5)) = 0.014 m-1
A(0.441T0) = 0.014×0.441×5 = 0.031 / ( <1,公式适用 )
d’ 1/0.014×arth(0.031) = 2.215 m ( = 1/Aarth(A(0.441T0)) )

3.试算法计算ζ以确定出逸段与非出逸段的分界点,假设ζ的试算范围ζ1～ζ2，计算的步长以及精度
ζ1 = 0.0001 = 0 m ( = precision )
ζ2 = 100 = 100 m ( = L2 )
步长STEP = 0.01 = 0.01 m
精度precision = 0.0001 = 0.0001 m
试算得到ζ = 64.46 = 64.46 m ( 手动输入 )
方程左边 = 0.105 = 0.105 m ( 手动输入 )
方程右边 = 0.105 = 0.105 m ( 手动输入 )

4.出逸段AB透水层水位计算(x以下游堤脚为原点，向下游为正)
计算点坐标x = 10 = 10 m ( 0 < x < 35.54 = L2-ζ )
h = 1.482 = 1.482 m

5.非出逸段BC透水层水位计算(x’以弱透水层最下游点为原点，向上游为正)
计算点坐标x’ = 10 = 10 m ( 0 < x’ < 64.46 = ζ )
x’/T0 = 10/5 = 2 m
Δ0 = (12-5)×0.441×5/(64.46+0.441×5) = 0.232 m
Δx’/Δ0 = 0 = 0 m
Δx’ = 0.232×0 = 0 m
h = 5+(12-5)×(10+0.441×5)/(64.46+0.441×5)-0 = 6.282 m

5.非出逸段BC透水层水位计算(x’以弱透水层最下游点为原点，向上游为正)
计算点坐标x’ = 2.1 = 2.1 m ( 0 < x’ < 64.46 = ζ )
x’/T0 = 2.1/5 = 0.42 m
Δ0 = (12-5)×0.441×5/(64.46+0.441×5) = 0.232 m
Δx’/Δ0 = 0.246 = 0.246 m ( 查表E.7.2 )
Δx’ = 0.232×0.246 = 0.057 m
h = 5+(12-5)×(2.1+0.441×5)/(64.46+0.441×5)-0.057 = 5.395 m

等我把这个发完，我会分享经验的

附录E.7 双层堤基渗流计算和覆盖的计算
附录E.7.2 透水地基上弱透水层不等厚或不均质，强透水层无限长双层堤基渗流计算

项目 计算式 数值 单位 备注

1.计算参数
上游坡度m1 = 3 = 3 / ( 1:m1 )
下游坡度m2 = 5 = 5 / ( 1:m2 )
堤顶宽度B = 6 = 6 m
堤顶高程▽顶 = 27 = 27 m
堤底高程▽底 = 17 = 17 m
弱透水地基底高程▽弱 = 10 = 10 m
强透水地基底高程▽强 = 7 = 7 m
上游水位▽1 = 24.8 = 24.8 m
下游水位▽2 = 18 = 18 m
弱透水层 无限长
强透水堤基渗透系数k0 = 7.00E-03 = 7.00E-03 m/s ( ＝ 0.7cm/s ＝ 604.8m/d )

2.递推计算下游等效长度S下
下游弱透水层分段数N = 3 = 3 段 ( 以不等厚或不均质分 )

2.1递推计算
段次 = 0 = 0
强透水层厚度T0 = 10-7 = 3 m ( ▽弱-▽强 )
等效长度S0 = 0 = 0 m ( 取0 )
D0 = (1/0.022+0)/(1/0.022-0) = 1 /
段次 = 1
第1段渗透系数k1 = 5.00E-05 = 5.00E-05 m/s ( ＝ 0.005cm/s ＝ 4.32m/d )
第1段厚度t1 = 5 = 5 m
第1段长度L1 = 6 = 6 m
越流系数A1 = SQRT(0.00005/(0.007×3×5)) = 0.022 m-1
β1 = 2×0.022×6 = 0.264 / ( 2AiLi )
等效长度S1 = 1/0.022×(1×e^0.264-1)/(1×e^0.264+1) = 5.965 m
D1 = (1/0.028+5.965)/(1/0.028-5.965) = 1.401 m
段次 = 2
第2段渗透系数k2 = 5.00E-05 = 5.00E-05 m/s ( ＝ 0.005cm/s ＝ 4.32m/d )
第2段厚度t2 = 3 = 3 m
第2段长度L2 = 10 = 10 m
越流系数A2 = SQRT(0.00005/(0.007×3×3)) = 0.028 m-1
β2 = 2×0.028×10 = 0.56 / ( 2AiLi )
等效长度S2 = 1/0.028×(1.401×e^0.56-1)/(1.401×e^0.56+1) = 15.026 m
D2 = (1/0.035+15.026)/(1/0.035-15.026) = 3.219 m ( ▽底-▽弱 )
段次 = 3
第3段渗透系数k3 = 5.00E-05 = 5.00E-05 m/s ( ＝ 0.005cm/s ＝ 4.32m/d )
第3段厚度t3 = 2 = 2 m
第3段长度L3 = 30 = 30 m
越流系数A3 = SQRT(0.00005/(0.007×3×2)) = 0.035 m-1
β3 = 2×0.035×30 = 2.1 / ( 2AiLi )
等效长度S3 = 1/0.035×(3.219×e^2.1-1)/(3.219×e^2.1+1) = 26.477 m

2.2下游等效长度S下 = 26.477 = 26.477 m ( 即等效长度S3 )

3.递推计算下游等效长度S上
上游弱透水层分段数N = 3 = 3 段 ( 以不等厚或不均质分 )

3.1 递推计算
段次 = 0 = 0
强透水层厚度T0 = 15.026-3.219 = 11.807 m ( ▽弱-▽强 )
等效长度S0 = 0 = 0 m ( 取0 )
D0 = (1/0.018+0)/(1/0.018-0) = 1 /
段次 = 1
第1段渗透系数k1 = 5.00E-05 = 5.00E-05 m/s ( ＝ 0.005cm/s ＝ 4.32m/d )
第1段厚度t1 = 7 = 7 m
第1段长度L1 = 10 = 10 m
越流系数A1 = SQRT(0.00005/(0.007×3×7)) = 0.018 m-1
β1 = 2×0.018×10 = 0.36 / ( 2AiLi )
等效长度S1 = 1/0.018×(1×e^0.36-1)/(1×e^0.36+1) = 9.893 m
D1 = (1/0.022+9.893)/(1/0.022-9.893) = 1.556 m
段次 = 2
第2段渗透系数k2 = 5.00E-05 = 5.00E-05 m/s ( ＝ 0.005cm/s ＝ 4.32m/d )
第2段厚度t2 = 5 = 5 m
第2段长度L2 = 6 = 6 m
越流系数A2 = SQRT(0.00005/(0.007×3×5)) = 0.022 m-1
β2 = 2×0.022×6 = 0.264 / ( 2AiLi )
等效长度S2 = 1/0.022×(1.556×e^0.264-1)/(1.556×e^0.264+1) = 15.413 m
D2 = (1/0.028+15.413)/(1/0.028-15.413) = 2.518 m ( ▽底-▽弱 )
段次 = 3
第3段渗透系数k3 = 5.00E-05 = 5.00E-05 m/s ( ＝ 0.005cm/s ＝ 4.32m/d )
第3段厚度t3 = 3 = 3 m
第3段长度L3 = 50 = 50 m
越流系数A3 = SQRT(0.00005/(0.007×3×3)) = 0.028 m-1
β3 = 2×0.028×50 = 2.8 / ( 2AiLi )
等效长度S3 = 1/0.028×(2.518×e^2.8-1)/(2.518×e^2.8+1) = 34.03 m

3.2 上游等效长度S上 = 34.03 = 34.03 m ( 即等效长度S3 )

4. 背水侧弱透水层下各点的承压水头计算(x以下游堤脚原点，向下游为正)
计算点坐标x = 10 = 10 m ( x > 0 )
上下游水头差H = 24.8-18 = 6.8 m
堤身高度Hd = 27-17 = 10 m ( ▽顶-▽底 )
堤底宽度b = 6+10×(3+5) = 86 m ( B+Hd×(m1+m2) )
计算点的承压水头h = (26.477-10)×6.8/(34.03+86+26.477) = 0.765 m ( x = 10m )