二00四年月
1 渗流根理
11 水工渗流危害渗流分析计算务
流体孔介质中运动称渗流水常见流体水利水电工程中广泛建 造堤坝围堰水闸等挡水建筑物形成水头差建筑物基通常透水孔介质水工渗流现象格外普遍
水工渗流造成方面危害渗流造成水库渠道水量损失渗流堤坝围堰土体饱降低坝体效容重抗剪强度导致坝坡失稳建筑物基渗流建筑物底部产 生扬压力利建筑物稳定渗流流速时造成坝体建筑物基土体 颗粒流失发生渗透变形堤坝倒塌建筑物滑移倾覆水库渗流引起游水位升导致农田冷浸渍害盐碱化作物减产拦污坝渗流造成水环境污染
水工渗流分析计算务争水渗流区域渗流流速流量水头分布浸润线等实行合理渗流掌握措施供避开减缓渗流危害
12 达西定律
19 世纪 50 年月法国工程师亨利·达西〔HDarcy〕通装圆筒中均质砂土进
行渗透试验觉察通两渗流断面间均渗流流速正两断面间水头差△h反 渗径长度L土粒构造流体性质关著名达西定律公式表达：
Dh dh
v k
k kJ (121)
L ds
h z +
p +a v 2
式中 h—测压水头
g 2g
z 位置高度p 压强γ水容重
渗流流速般流速水头a
v2 视h z + p
2g g
k—反映土粒构造流体性质系数渗透系数某具体流体〔方水〕言k 值仅土粒构造关
J—渗透坡降 J dh
ds
式中负号表示水总流水头减方
应留意达西定律中流速全断面均流速v土体孔隙中流速v 两种流速存关系：
v nv (122)
式中 n 体积孔隙率见达西流速土体孔隙中流速
应留意达西定律适层流状态渗流运动水利工程中堆石坝堆石排水体等孔隙介质中渗流 流外绝数渗流属层流达西定律适非层流渗流流淌规律公式形式表达
1
v kJ m
〔123〕
式中m1 时层流渗流m2 时完全紊流渗流1＜m＜2 时层流紊流渡区
式〔121〕等号两边xyz 轴投影便空间直角坐标系中达西公式：
v k
x
¶h
x ¶x
¶h
y ¶y
¶h
k J
x x
v k
y
k J
y y
〔124〕
v k
z z
¶z kz J z
123 渗流运动连续性方程
z
dz
x
dy
dx
y
图 121 微元体进出流量示意
渗流场中取某微分单元体〔图121〕体积V dxdydz 流入左侧面水体
质量速率r v 进出流量差
dydz 流出右侧面水体质量速率(r v
x z
¶
+ ¶z vz
dz)dxdy 左右面
¶ r v
¶x x
dxdydz
样前面面作流进流出流入量计算终累加净流入量微元体总流入量
÷
æ ¶ r ¶ r ¶ r ö
ç ¶x
vx + ¶y
vy + ¶z vz V
è ø
式开放
æ ¶v ¶v ¶v ö æ ¶r ¶r ¶r ö
÷
rç ¶ x + y + z ÷V ç v
+ v + v V
x
y z
è x ¶y ¶z ø
è ¶x
¶y ¶z ø
式中项前项相视改写成
æ ¶v
r ç ¶ x
+ ¶vy
+ ¶vz
ö
÷V 〔125〕
è x ¶y ¶z ø
式水体质量微元体积存速率质量守恒原理应等微元体水体
质量 M 时间变化速率：
¶M
¶(nrV ) nr ¶V
+ rV
¶n + nV
¶r
〔126〕
¶t ¶t ¶t ¶t ¶t
式中n 土体孔隙率ρ 水密度V 微元体体积式〔126〕右边三项分代表土体骨架孔隙体积流体密度转变速率
引入弹性压缩理导出〔126〕式中等价表达：
¶M r 2g(a + nb )V
¶h
〔127〕
¶t ¶t
式〔127〕中 α ——土颗粒骨架压缩性〔压缩模量〕
β ——水压缩性
h——渗流测压水头〔m〕 h z + p
r g
质量守恒原理式〔125〕式〔127〕相等
æ ¶v
¶v ¶v ö
y
( )¶h
ç ¶ x + + z ÷ rg a + nb
〔128〕
è x ¶y ¶z ø ¶t
假定水土行压缩时式变：
x
¶v + ¶vy
¶
+ ¶vz
0 〔129〕
x ¶y ¶z
式〔129〕行压缩流体刚体介质中流淌连续性方程
124 渗流微分方程
达西公式代入渗流连续性方程〔128〕：
¶ æ ¶h ö ¶ æ ¶h ö ¶ æ ¶h ö
k
r (a
b )¶h ¶h
¶x ç
x ¶x ÷ +
¶y ç k y
¶y ÷ +
ç k ÷ g + n
¶z z ¶z
S
¶t s ¶t
〔1210〕
è ø è ø è ø
式〔1210〕非稳定渗流微分方程般形式适合承压含水层适合压渗流式中 S rg(a + nb )称单位储存量〔尺度 1L〕含义：单位体积饱土体
s
水头降 1m 时土体压缩〔 rg a 〕水体膨胀〔 rgn b 〕释放出储存水量
介质均质性〔k k k
x y z
s
¶2h ¶2h ¶2h S ¶h
k 〕时式〔1210〕变
¶x2
+ ¶y2
+ ¶z2
k ¶t
〔1211〕
假定水土行压缩时式〔1210〕〔1211〕分变成
¶ æ ¶h ö ¶ æ ¶h ö ¶ æ ¶h ö
ç
÷
ç k ÷ + k
¶ x y
+ ç k ÷ 0 〔1212〕
z
x è ¶x ø ¶y è ¶y ø ¶z è ¶z ø
¶2h +
¶x2
¶ h
¶y2
+ ¶2h 0
¶z2
〔1213〕
式〔1212〕〔1213〕分异性性时三维稳定渗流微分方程
堤坝渗流建筑物基渗流状况简化垂直剖面二维渗流问题式〔1212〕〔1213〕变成
¶ æ ¶h ö ¶ æ ¶h ö
¶
ç k ÷ + ç k ÷ 0 〔1214〕
x z
x è ¶x ø ¶z è ¶z ø
¶2h +
¶x2
¶2h 0
¶z2
(1215)
125 面变动渗流微分方程
土坝〔堤防〕水库〔河道〕水位降时非稳定渗流状况面降 引起土体压缩弹性释放水量面降时排出水量相令式〔1210〕
中 S 0 面三维非稳定渗流微分方程
s
¶ æ ¶h ö ¶ æ ¶h ö ¶ æ ¶h ö
ç k ÷ +
¶ x
ç k ÷ + ç k
y
z
÷ 0 〔1216〕
x è ¶x ø ¶y è
¶y ø
¶z è
¶z ø
¶ æ
¶
x
ç
è
¶h ö
k
x
¶
x
÷
ø
+
¶ æ
¶
z
ç
è
¶h ö
k
z
¶
z
÷
ø
0
〔1217〕
面二式然形式稳定渗流方程式〔1212〕〔1214〕完全样结合
面变动边界条件水头分布空间座标时间函数象稳定渗流方程式解答空间座标函数
126 水工渗流问题定解条件
发生限空间流场渗流运动仅受渗流微分方程支配受流场边界条件初始渗流状态条件支配面土坝非稳定渗流状态介绍水工渗流定解条件种类型
二维非稳定渗流微分方程
图 122 土坝非稳定渗流定解条件
初始条件： h边界条件：
t 0
h (x zt) 〔仅非稳定渗流〕 〔1218〕
0
⑴水头边界 h f (x zt) 〔1219〕
G1 1
图 122 中水头边界包括：
游水头边界： h H (t)
ABC 1
游水头边界： h H
EFG 2
游坡渗出段水头边界： h z
DE
⑵流量边界： q k
n
f
¶h
¶n
G 2 2
(x zt) (1220)
图 122 中面AD 降时面流进单宽流量：
q m ¶h cosq 〔A〕
¶t
式中 θ —面法线铅直线夹角
μ —水度单位体积土体饱含水状况面降 1m 排出空隙水量称土体排水空隙率效空隙率
q k ¶h
n ¶n
k
n
¶h cosq 〔B〕
¶z
式〔A〕〔B〕面流量边界条件改写成：
m ¶h ¶h 〔C〕
k ¶t ¶z
n
坝体土异性时 k 取 k 值
n z
渗流面气相通压力水头零面应满足：
h z *
2 SEEP V30 功特色
（1） 适种工况〔水位固定游水位快降〕简单边界〔包括种防渗排水设施 〕简单土层条件：堤防土石坝渗流计算类透水材料坝〔混凝土坝浆砌石坝〕 水闸泵房基渗流计算
（2） 承受限元法计算渗流量浸润线等水头线等计算结果
（3） 承受视化界面建立原始数文件记事中建立原始数文件
（4） 计算成果印成表格动绘制出浸润线等水头线图
3 SEEP V30 解题点
31 坐标系选取计算范围确实定
图 311 示水坐标x 取右正竖直坐标高程体系全
游入渗段长度L1 游渗出段长度L2般取05~1 倍水头〔取实测值〕
H1
L1
O
L2
H2
Y
图 311
32 单元剖分
a) 单元动剖分分组分块结点编号
图 321
b) 四条边剖分点布置
软件进展渗流限元计算根微单元三角形单元需述四边形块单 元剖分三角形单元图 322 示渗流场中取出意块图中 n1n2n3n4 分表示相应边剖分点数〔剖分点数算三角形单元数目〕时针方编号n1 n3n2 n4 分两组边两边剖分点数相等相等避开出错边点数相等相差剖分点疏密应水力条件确定 部位剖分单元应加密：水力坡降部位坝轴线四周截水槽游铺盖进水处 游出口排水沟四周外剖分点数拟定必需满足三点：
① 三角形单元总数 < 500
② n n
3 1
< n n
2 4
n n
4 2
< n n
1 3
X
渗流场外形土层分布等整渗流场划分四边形子域
图 321 示子域称块结点优化满足面调整迭代需划分块分成假设干组组假设干块组成具体规定：组块首尾相接贯穿整渗流场：组块逐块编号左右组组进展块结点标号意 〔建议左右挨次〕图 321 中结点 37111215 连成折线事先假定浸润线
H1 3〔2〕〔2〕
〔3〕
7
11
4
〔3〕
［2
］ 8［5］12
〔3〕〔3〕
〔3〕
［3］
〔3〕〔3〕 3
［6］〔〕
1
〔4〕
〔3〕 ［1］
5
〔3〕
6
〔4〕
［4］
〔5〕
9 〔3〕 13
〔3〕［7］〔3〕
1〔0 2〕14
〔5〕
［8］
〔5〕
［9］
〔5〕
15
〔3〕
〔3〕16［10］〔4〕
18
〔3〕19 〔4〕
H2
21
2
〔4〕
〔3〕
17
［11］
〔4〕
〔3〕 ［12］ 〔3〕
20
〔4〕
22
③浸润线变动范围n n
1 3
4
〔5〕
n1
7
n2
〔3〕
〔3〕n4
5
n3
(6)
8
图 322
c) 剖分挨次
剖分挨次次：划分块单元分组编组号编块号核块单元结点数编结点号块边布置剖分点数
块渗透系数必需样渗透系数土层划分假设干块
块必需四边形块渗流区域边界时行避开消灭三角形块时该块边界边取点作结点样三角形块成蜕化四边形块
4 SEEP V30 操作步骤文件容介绍
41 DQB 文件产生
具体操作步骤：
第步：运行CCEDEXE 文件弹出工程视图相关参数界面(图 411)：
图 411 工程视图参数界面
工程相关参数：
(1) 填入分块组数土层种类数
(2) 等势线间距 5 10两数供选择分表示渗流场分布图中相邻两条等势线间水头势相差 5 10
(3) 渗流状况：假设计算非稳定渗流需填写水位降落受时间初始水位降水位假设计算稳定渗流需选择否计算渗流量
视图相关参数：
(1) 坐标坐标：您绘图程中输入坐标值便利绘图建议断面XY 坐标值分加 30断面XY 坐标值分减 30 输入
(2) 1m( )P表示绘图时米长度显示象素数假设图形中较短线条
填入较值方 10假设图形中线条较长建议承受默认值4
(3) 点选取半径 默认 3 素设置更捕获点会影响线条选取建议承受默认值字体般默认 10
图 412 图形绘制界面
次步：单击工程视图参数界面中确实定钮进入图形绘制界面(图 412)：
默认状态 进展图形绘制：
(1) 图形绘制：右角Line From面直接输入节点坐标回车键
(2) 假设节点坐标错误点击
进展修改
(3) 点击切换状态
左键选中线条单击右键删线条进展单元剖分：
(1) 节点编号：编号挨次直接点击该节点假设选中该节点点击 加
点选取半径
(2) 剖分掌握：先选中块四条边〔假设选中该边点击
减点选取
半径者图形放〕然点击
块剖分数填写文件保存：
进展剖分填写终块编号挨次进展
(1) 填写土层渗透系数：点击 输入层土水垂直渗透系数 K K 〔cms〕
x y
留意：输入层土渗透系数点击确定钮进展数保存
(2) 剖分块信息修改：假设剖分信息错误点击编号单击修改钮〕
进展修改〔先选中Reg 某块
(3) 设置相关参数：点击 翻开相关参数设置话框检查分组数等势线百分
数间距渗流状况否正确留意：计算工况稳定渗流计算渗流量时需填写渗流断面穿块号
(4) 保存dqbdat 文件：点击弹出 dqb 文件保存话框选择保存路径点击保存(S)
钮留意：计算时文件名必需dqbdat
(5) 生成断面轮廓图：假设需断面轮廓图直接图中点击鼠标左键选中全部保存
轮廓线点击 弹出保存话框填入文件名保存
42 运行计算
点击 钮弹出计算操作界面(图 42)：
图 42 计算操作界面
(1) 点击查 DQB 钮检查产生 DQB 文件否正确(参考说明书中 43 DQB
DAT 文件容格式)
(2) 点击剖分单元钮生成渗流计算需原始数文件G·DAT 剖分图NET0BLN
(3) 点击显示剖分图钮进展单元剖分时生成网格图NETBLN
(4) 点击编辑数钮翻开计算原始数文件G·DAT留意：该文件末尾起行输入 0 保存
(5) 点击计算钮 生成结果文件SEEP·OUT 渗流场图形文件GS·BLN
(6) 点击等势线图启动 SURFER绘制渗流场分布图：
a单击MAP 菜单中LOAD BASE MAP 消灭 IMPORT FILE 视窗双击BOUND·BLN 单击 OK消灭土坝横断面图
b单击MAP 菜单中LOAD BASE MAP 消灭 IMPORT FILE 视窗双击GS·BLN单击 OK消灭渗流场图形
c土坝横断面图渗流场图形进展叠加：鼠标渗流场图形拖堤防横断面 图鼠标时框选堤防横断面图单渗流场图形单击 MAP 消灭MAP 菜单 单击 MAP 菜单中OVERLAY MAPS 完成图形叠加
d进展文字编辑坐标轴设定等终完善渗流场分布图
(7) 假设想查计算结果SEEP·OUT 点击计算报告钮
43 DQBDAT 文件容格式
产生数文件DQBDAT 容格式
第行：SEEP1 SEEP2 SEEP3
次行：ZUKUAIDIANBIANTUCENGSLDMSLKSDBFS
第三行：KXX〔1：KUAI1：9〕 第四行：ZB〔1：DIAN1：2〕
第五行：SLCS〔1：TUCENG1：3〕 第六行：BJXX〔1：BIAN1：3〕 第七行：SLKH〔1：SLK〕
第八行：SLDMXX〔1：SLDM1：3〕 第九行：SJSW〔1：21：2〕
数说明：
（1） SEEP 1—判否保存渗流绘图数Y’确定N’否认 SEEP 2—判否动剖分数Y’确定N’否认 SEEP 3—设定N’
（2） ZU—分块组数〔计算区域剖分块次左右分假设干组组间应关系第 I 块边〔n3〕第 I+1 块边〔n1〕应边组挨次编号〕
KUAI—四边形块数目
DIAN—块角点数 BIAN—含类边界边数 TUCENG—渗透系数土层种类数
SLDM—需计算渗流量渗流断面数假设零表示计算渗流量 SLK—需计算渗流量断面穿总计块数假设零表示计算渗流量 SD—计算非稳定渗流时水库水位降落时段数假设零表示计算稳定渗流 BFS—计算等势线百分数间距值默认10
（3） KXX(1ME19)块信息数组
块 9 重量块号挨次填写第1—4 重量表示块左角起时针方四角点点号第 5—8 重量次表示块边起时针方4 条边剖分点数第9 重量块渗透系数区号
（4） ZB〔1：DIAN1：2〕块结点坐标数组
块结点序号挨次表示结点XY 坐标
（5） SLCS〔1：TUCENG1：3〕
次层土水渗透系数Kx〔单位：cms〕垂直渗透系数 Ky〔单位：cms〕水度m m 01177 k 〔该式中渗透系数k 单位：mday〕
（6） BJXX(1BN13) 边界信息数组
条边界三重量第重量边界边块块号次重量该边界边块中边序号第三重量边边界类型号规定：
1——渗流面边界〔斜心墙坝仅部面〕
3——半压流边界
5——斜心墙坝游坝壳面边界
6——均质土坝游坝坡排水体渗出段边界
7——游水头边界
8——游水头边界
9——斜心墙土坝斜墙心墙坡渗出段边界
（7） SLKH〔1：SLK〕
需计算渗流量流断面块号
（8） SLDMXX〔1：SLDM1：3〕渗流断面块单元信息
第重量某流断面块数次重量该流断面块边界相交信息 水工渗流计算般 1第三重量该流断面块中位置节点链序号般1 2
（9） SJSW〔1：21：2〕起降时间〔天〕水位降时间〔天〕水位
44 GDAT 数容
软件四边形块单元动剖分三角形单元 时 DQBDAT 动产生
GDATGDAT 软件进展渗流计算原始数文件数容格式
第行：SEEP1SEEP2SEEP3
次行：NENNNNUNNDNNM1NNM2NNFTUCENGAAAA1AA2
BB1BBBB2NNPNQSDNWNGNNM3LULDMA
第三行：H H1H2H3EPSEPS1DTY X
max
AA AA
第十四行：SW〔1：SD〕 第十五行：SJJG〔1：SD〕
第十六行：SLCS〔1：TUCENG1：3〕 第十七行：QTE(1：NQ)
第十八行：QEN(1：Q )
E
解释：
（1） SEEP1SEEP2SEEP3—DQBDAT
（2） NE—三角形单元总数
〔3〕H
—游高水位
max
H1—游水位 H2—游水位
H3—井水位 EPS—面迭代允许误差(01~0001) EPS1—井四周迭代允许误差(05~0005) DT—水位变化总历时〔单位：天〕
（14） SW〔1：SD〕
时段末游水位
（15） SJJG〔1：SD〕 时段时间间隔天数
（16） SLCS〔1：TUCENG1：3〕
次填写层土水渗透系数K 〔单位：mday〕垂直渗透系数K 〔单位：mday〕
x y
水度m m 01177 k 〔式中k 渗透系数单位：mday〕
（17） QTE(1：NQ)
计算渗流量断面单元总数
（18） QEN(1：Q )
E
计算渗流量断面全部单元号
5 算例
例 均质土坝稳定渗流计算
某均质土坝横断面尺寸游水位渗透系数图51 示求计算稳定渗流场单位坝长渗漏流量
图 51 均质坝横断面
解题步骤：
(1) 运行CCEDEXE 文件弹出工程视图相关参数界面分块组数中填入5土层种类数 中填入 2选择稳定渗流选择计算渗流量
(2) 单击视图相关参数坐标中 X 填写 170Y 填写 100坐标中 X 填写30Y 填写 10默认值
(3) 单击确定钮进入图形绘制界面参考41 DQB 文件产生中具体操作步骤建立均质坝横断面剖分图(图 52)
图 52 均质坝横断面剖分图
(4) 保存dqbdat 文件容：
Y Y N
5 8 16 7 2 1 2 0 10
1 4 5 2 3 2 3 2 1
3 6 7 4 2 3 4 4 2
4 7 8 5 4 2 4 2 1
6 9 10 7 4 4 4 3 2
7 10 11 8 4 2 4 2 1
9 12 13 10 8 4 8 4 2
10 13
14 11 8 2 8 2 1
13 15
16 14 4 2 4 2 1
0000
45000
0000
40000
40000 55000
10000 45000
10000
40000
42000
55000
42000
45000
42000
40000
67000
65000
67000
45000
67000
40000
121000
48000
118000
45000
118000
40000
140000
45000
140000
40000
0000002023 0000002023 0047157335
0000010000 0000010000 0059347456
1 1 7
2 4 7
2 1 7
4 1 7
6 1 1
6 2 8
8 1 8
6 7
2 1 2
(5) 点击图形绘制界面中 钮弹出计算操作界面：
a 点击查 DQB 钮检查产生DQB 文件否正确(参考说明书中 43 DQBDAT 文件容格式)
b 点击剖分单元钮生成渗流计算需原始数文件G·DAT 剖分图NET0BLN
c 点击显示剖分图钮进展剖分时生成网格图NETBLN
d 点击编辑数钮翻开计算原始数文件G·DAT留意：该文件末尾起行输入0 保存
e 点击计算钮 生成结果文件SEEP·OUT 渗流场图形文件GS·BLN
f 点击等势线图启动 SURFER观察渗流场分布图：
选择 MAP 菜单中LOAD BASE MAP 消灭 IMPORT FILE 视窗双击BOUND·BLN单击OK 消灭土坝横断面图
选择 MAP 菜单中LOAD BASE MAP 消灭 IMPORT FILE 视窗双击GS·BLN单击 OK 消灭渗流场图形
土坝横断面图渗流场图形进展叠加：鼠标渗流场图形拖堤防横断面图时框选堤防横断面图单渗流场图形选择MAP 菜单中OVERLAY MAPS 完成图形叠加
650
9
0
8
0
7
0
6
0
650
进展文字编辑坐标轴设定等终完善渗流场分布图〔图53〕
图 53 均质坝渗流场分布图
翻开计算报告〔SEEP·OUT〕单位坝长日渗漏流量 QQ100313〔m3daym〕
例二 堤防非稳定渗流计算
堤防堤基横断面尺寸渗透系数图54 示游水位 310m河水位 380m
45 天降 360m 35 天降 340m求计算河水位降落程中堤防渗流场
图 54 堤防横断面
解题步骤参考例建立堤防横断面剖分图〔图55〕
图 55 堤防横断面剖分图
产生问题DQB·DAT 文件容：
Y Y N
3 5 12 7 2 0 0 2 10
1 5 6 2 4 5 4 5 2
3 7 8 4 13 3 13 3 1
4 8 9 5 13 3 15 3 1
5 9 10 6 15 5 15 5 2
9 11 12 10 4 5 4 5 2
0000 30000
0000 20000
30000
38000
20000
34000
10000
30000
10000
20000
53000
33500
58000
31000
60000
30000
60000
20000
70000
30000
70000 20000
0000010000 0000010000 0059347456
0000002023 0000002023 0047157335
1 1 7
3 4 7
2 4 7
2 1 1
2 2 6
3 2 8
5 1 8
0 38000 8000 34000
留意：GSBLN 中包含初始水位 380降瞬时水位 360降瞬时水位 340 3 水位时渗流场绘图数应GSBLN 文件拆分成 3 绘图数文件 GS1BLNGS2BLNGS3BLN分 GS1BLNGS2BLNGS3BLN 堤防横断面图 LUNKE2BLN
叠加然进展文字编辑坐标轴设定等终完善渗流场分布图〔图56〕
40
380
35
310
30
0
01
0 0
6 5
0
4
0
3
2
40
45
50
55
60
65
70
25
0 0 0
9 8 7
20
0 5 10 15 20 25 30 35
i 初始水位 380 时渗流场分布图
3 8 0
90
80
3 1 0
7
0
40
35
30
0
6
0
5
0
4
0
3
0
2
01
40
45
50
55
60
65 70
25
20
0 5 10 15 20 25 30 35
40
380
35
ii
瞬时水位 360 时渗流场分布图
80
30
0 0
25
70
0
310
0
0
0 01
6 6 5 4 3 2
20
0 5 10 15 20 25 30 35 40
45 50
55 60 65 70
iii 瞬时水位 340 时渗流场分布图
图 56 水位降落程中堤防渗流场分布图
例三 心墙坝稳定渗流计算
某心墙土坝横断面尺寸渗透系数游水位图57 示求计算稳定渗流状况土坝渗流场
图 57 心墙坝横断面
解题步骤参考例建立心墙坝横断面剖分图〔图58〕：
图 58 心墙坝横断面剖分图
产生问题DQB·DAT 文件容：
Y Y N
6 12 22 10 3 1 2 0 10
1 5 6 2 4 3 4 3 3
3 7 8 4 2 3 3 3 1
4 8 9 5 3 3 4 3 1
5 9 10 6 4 3 5 3 3
7 11 12 8 2 3 3 3 2
8 12 13 9 3 3 3 3 2
9 13 14 10 3 3 2 3 2
12
15
16
13
5
3
5
3
1
13
16
17
14
5
3
5
3
3
15
18
19
16
3
3
4
3
1
16
19
20
17
4
3
4
3
3
19
21
22
20
4
3
4
3
3
0000 45000
0000 40000
69000 63000
42000 54000
15000
45000
15000
40000
75000
63000
72750
54000
70500
45000
74500
40000
76500
63000
79000
54000
81000
45000
77000
40000
118500
51000
118500
45000
118500
40000
126750
48000
135000
45000
135000
40000
150000
45000
150000
40000
0001000000
0001000000
0114581999
0000020230
0000020230
0065524904
0000500000
0000500000
0103779628
1
1
7
3
4
7
2
4
7
2
1
1
5
1
1
5
2
9
8
1
5
10
1
5
10
2
8
12
1
8
8 9
2 1 2
心墙坝渗流场分布图图 59：
630
0
480
2
0
1
图 59 心墙坝渗流场分布图
文档香网(httpswwwxiangdangnet)户传
《香当网》用户分享的内容，不代表《香当网》观点或立场，请自行判断内容的真实性和可靠性！
该内容是文档的文本内容，更好的格式请下载文档
二00四年月
1 渗流根理
11 水工渗流危害渗流分析计算务
流体孔介质中运动称渗流水常见流体水利水电工程中广泛建 造堤坝围堰水闸等挡水建筑物形成水头差建筑物基通常透水孔介质水工渗流现象格外普遍
水工渗流造成方面危害渗流造成水库渠道水量损失渗流堤坝围堰土体饱降低坝体效容重抗剪强度导致坝坡失稳建筑物基渗流建筑物底部产 生扬压力利建筑物稳定渗流流速时造成坝体建筑物基土体 颗粒流失发生渗透变形堤坝倒塌建筑物滑移倾覆水库渗流引起游水位升导致农田冷浸渍害盐碱化作物减产拦污坝渗流造成水环境污染
水工渗流分析计算务争水渗流区域渗流流速流量水头分布浸润线等实行合理渗流掌握措施供避开减缓渗流危害
12 达西定律
19 世纪 50 年月法国工程师亨利·达西〔HDarcy〕通装圆筒中均质砂土进
行渗透试验觉察通两渗流断面间均渗流流速正两断面间水头差△h反 渗径长度L土粒构造流体性质关著名达西定律公式表达：
Dh dh
v k
k kJ (121)
L ds
h z +
p +a v 2
式中 h—测压水头
g 2g
z 位置高度p 压强γ水容重
渗流流速般流速水头a
v2 视h z + p
2g g
k—反映土粒构造流体性质系数渗透系数某具体流体〔方水〕言k 值仅土粒构造关
J—渗透坡降 J dh
ds
式中负号表示水总流水头减方
应留意达西定律中流速全断面均流速v土体孔隙中流速v 两种流速存关系：
v nv (122)
式中 n 体积孔隙率见达西流速土体孔隙中流速
应留意达西定律适层流状态渗流运动水利工程中堆石坝堆石排水体等孔隙介质中渗流 流外绝数渗流属层流达西定律适非层流渗流流淌规律公式形式表达
1
v kJ m
〔123〕
式中m1 时层流渗流m2 时完全紊流渗流1＜m＜2 时层流紊流渡区
式〔121〕等号两边xyz 轴投影便空间直角坐标系中达西公式：
v k
x
¶h
x ¶x
¶h
y ¶y
¶h
k J
x x
v k
y
k J
y y
〔124〕
v k
z z
¶z kz J z
123 渗流运动连续性方程
z
dz
x
dy
dx
y
图 121 微元体进出流量示意
渗流场中取某微分单元体〔图121〕体积V dxdydz 流入左侧面水体
质量速率r v 进出流量差
dydz 流出右侧面水体质量速率(r v
x z
¶
+ ¶z vz
dz)dxdy 左右面
¶ r v
¶x x
dxdydz
样前面面作流进流出流入量计算终累加净流入量微元体总流入量
÷
æ ¶ r ¶ r ¶ r ö
ç ¶x
vx + ¶y
vy + ¶z vz V
è ø
式开放
æ ¶v ¶v ¶v ö æ ¶r ¶r ¶r ö
÷
rç ¶ x + y + z ÷V ç v
+ v + v V
x
y z
è x ¶y ¶z ø
è ¶x
¶y ¶z ø
式中项前项相视改写成
æ ¶v
r ç ¶ x
+ ¶vy
+ ¶vz
ö
÷V 〔125〕
è x ¶y ¶z ø
式水体质量微元体积存速率质量守恒原理应等微元体水体
质量 M 时间变化速率：
¶M
¶(nrV ) nr ¶V
+ rV
¶n + nV
¶r
〔126〕
¶t ¶t ¶t ¶t ¶t
式中n 土体孔隙率ρ 水密度V 微元体体积式〔126〕右边三项分代表土体骨架孔隙体积流体密度转变速率
引入弹性压缩理导出〔126〕式中等价表达：
¶M r 2g(a + nb )V
¶h
〔127〕
¶t ¶t
式〔127〕中 α ——土颗粒骨架压缩性〔压缩模量〕
β ——水压缩性
h——渗流测压水头〔m〕 h z + p
r g
质量守恒原理式〔125〕式〔127〕相等
æ ¶v
¶v ¶v ö
y
( )¶h
ç ¶ x + + z ÷ rg a + nb
〔128〕
è x ¶y ¶z ø ¶t
假定水土行压缩时式变：
x
¶v + ¶vy
¶
+ ¶vz
0 〔129〕
x ¶y ¶z
式〔129〕行压缩流体刚体介质中流淌连续性方程
124 渗流微分方程
达西公式代入渗流连续性方程〔128〕：
¶ æ ¶h ö ¶ æ ¶h ö ¶ æ ¶h ö
k
r (a
b )¶h ¶h
¶x ç
x ¶x ÷ +
¶y ç k y
¶y ÷ +
ç k ÷ g + n
¶z z ¶z
S
¶t s ¶t
〔1210〕
è ø è ø è ø
式〔1210〕非稳定渗流微分方程般形式适合承压含水层适合压渗流式中 S rg(a + nb )称单位储存量〔尺度 1L〕含义：单位体积饱土体
s
水头降 1m 时土体压缩〔 rg a 〕水体膨胀〔 rgn b 〕释放出储存水量
介质均质性〔k k k
x y z
s
¶2h ¶2h ¶2h S ¶h
k 〕时式〔1210〕变
¶x2
+ ¶y2
+ ¶z2
k ¶t
〔1211〕
假定水土行压缩时式〔1210〕〔1211〕分变成
¶ æ ¶h ö ¶ æ ¶h ö ¶ æ ¶h ö
ç
÷
ç k ÷ + k
¶ x y
+ ç k ÷ 0 〔1212〕
z
x è ¶x ø ¶y è ¶y ø ¶z è ¶z ø
¶2h +
¶x2
¶ h
¶y2
+ ¶2h 0
¶z2
〔1213〕
式〔1212〕〔1213〕分异性性时三维稳定渗流微分方程
堤坝渗流建筑物基渗流状况简化垂直剖面二维渗流问题式〔1212〕〔1213〕变成
¶ æ ¶h ö ¶ æ ¶h ö
¶
ç k ÷ + ç k ÷ 0 〔1214〕
x z
x è ¶x ø ¶z è ¶z ø
¶2h +
¶x2
¶2h 0
¶z2
(1215)
125 面变动渗流微分方程
土坝〔堤防〕水库〔河道〕水位降时非稳定渗流状况面降 引起土体压缩弹性释放水量面降时排出水量相令式〔1210〕
中 S 0 面三维非稳定渗流微分方程
s
¶ æ ¶h ö ¶ æ ¶h ö ¶ æ ¶h ö
ç k ÷ +
¶ x
ç k ÷ + ç k
y
z
÷ 0 〔1216〕
x è ¶x ø ¶y è
¶y ø
¶z è
¶z ø
¶ æ
¶
x
ç
è
¶h ö
k
x
¶
x
÷
ø
+
¶ æ
¶
z
ç
è
¶h ö
k
z
¶
z
÷
ø
0
〔1217〕
面二式然形式稳定渗流方程式〔1212〕〔1214〕完全样结合
面变动边界条件水头分布空间座标时间函数象稳定渗流方程式解答空间座标函数
126 水工渗流问题定解条件
发生限空间流场渗流运动仅受渗流微分方程支配受流场边界条件初始渗流状态条件支配面土坝非稳定渗流状态介绍水工渗流定解条件种类型
二维非稳定渗流微分方程
图 122 土坝非稳定渗流定解条件
初始条件： h边界条件：
t 0
h (x zt) 〔仅非稳定渗流〕 〔1218〕
0
⑴水头边界 h f (x zt) 〔1219〕
G1 1
图 122 中水头边界包括：
游水头边界： h H (t)
ABC 1
游水头边界： h H
EFG 2
游坡渗出段水头边界： h z
DE
⑵流量边界： q k
n
f
¶h
¶n
G 2 2
(x zt) (1220)
图 122 中面AD 降时面流进单宽流量：
q m ¶h cosq 〔A〕
¶t
式中 θ —面法线铅直线夹角
μ —水度单位体积土体饱含水状况面降 1m 排出空隙水量称土体排水空隙率效空隙率
q k ¶h
n ¶n
k
n
¶h cosq 〔B〕
¶z
式〔A〕〔B〕面流量边界条件改写成：
m ¶h ¶h 〔C〕
k ¶t ¶z
n
坝体土异性时 k 取 k 值
n z
渗流面气相通压力水头零面应满足：
h z *
2 SEEP V30 功特色
（1） 适种工况〔水位固定游水位快降〕简单边界〔包括种防渗排水设施 〕简单土层条件：堤防土石坝渗流计算类透水材料坝〔混凝土坝浆砌石坝〕 水闸泵房基渗流计算
（2） 承受限元法计算渗流量浸润线等水头线等计算结果
（3） 承受视化界面建立原始数文件记事中建立原始数文件
（4） 计算成果印成表格动绘制出浸润线等水头线图
3 SEEP V30 解题点
31 坐标系选取计算范围确实定
图 311 示水坐标x 取右正竖直坐标高程体系全
游入渗段长度L1 游渗出段长度L2般取05~1 倍水头〔取实测值〕
H1
L1
O
L2
H2
Y
图 311
32 单元剖分
a) 单元动剖分分组分块结点编号
图 321
b) 四条边剖分点布置
软件进展渗流限元计算根微单元三角形单元需述四边形块单 元剖分三角形单元图 322 示渗流场中取出意块图中 n1n2n3n4 分表示相应边剖分点数〔剖分点数算三角形单元数目〕时针方编号n1 n3n2 n4 分两组边两边剖分点数相等相等避开出错边点数相等相差剖分点疏密应水力条件确定 部位剖分单元应加密：水力坡降部位坝轴线四周截水槽游铺盖进水处 游出口排水沟四周外剖分点数拟定必需满足三点：
① 三角形单元总数 < 500
② n n
3 1
< n n
2 4
n n
4 2
< n n
1 3
X
渗流场外形土层分布等整渗流场划分四边形子域
图 321 示子域称块结点优化满足面调整迭代需划分块分成假设干组组假设干块组成具体规定：组块首尾相接贯穿整渗流场：组块逐块编号左右组组进展块结点标号意 〔建议左右挨次〕图 321 中结点 37111215 连成折线事先假定浸润线
H1 3〔2〕〔2〕
〔3〕
7
11
4
〔3〕
［2
］ 8［5］12
〔3〕〔3〕
〔3〕
［3］
〔3〕〔3〕 3
［6］〔〕
1
〔4〕
〔3〕 ［1］
5
〔3〕
6
〔4〕
［4］
〔5〕
9 〔3〕 13
〔3〕［7］〔3〕
1〔0 2〕14
〔5〕
［8］
〔5〕
［9］
〔5〕
15
〔3〕
〔3〕16［10］〔4〕
18
〔3〕19 〔4〕
H2
21
2
〔4〕
〔3〕
17
［11］
〔4〕
〔3〕 ［12］ 〔3〕
20
〔4〕
22
③浸润线变动范围n n
1 3
4
〔5〕
n1
7
n2
〔3〕
〔3〕n4
5
n3
(6)
8
图 322
c) 剖分挨次
剖分挨次次：划分块单元分组编组号编块号核块单元结点数编结点号块边布置剖分点数
块渗透系数必需样渗透系数土层划分假设干块
块必需四边形块渗流区域边界时行避开消灭三角形块时该块边界边取点作结点样三角形块成蜕化四边形块
4 SEEP V30 操作步骤文件容介绍
41 DQB 文件产生
具体操作步骤：
第步：运行CCEDEXE 文件弹出工程视图相关参数界面(图 411)：
图 411 工程视图参数界面
工程相关参数：
(1) 填入分块组数土层种类数
(2) 等势线间距 5 10两数供选择分表示渗流场分布图中相邻两条等势线间水头势相差 5 10
(3) 渗流状况：假设计算非稳定渗流需填写水位降落受时间初始水位降水位假设计算稳定渗流需选择否计算渗流量
视图相关参数：
(1) 坐标坐标：您绘图程中输入坐标值便利绘图建议断面XY 坐标值分加 30断面XY 坐标值分减 30 输入
(2) 1m( )P表示绘图时米长度显示象素数假设图形中较短线条
填入较值方 10假设图形中线条较长建议承受默认值4
(3) 点选取半径 默认 3 素设置更捕获点会影响线条选取建议承受默认值字体般默认 10
图 412 图形绘制界面
次步：单击工程视图参数界面中确实定钮进入图形绘制界面(图 412)：
默认状态 进展图形绘制：
(1) 图形绘制：右角Line From面直接输入节点坐标回车键
(2) 假设节点坐标错误点击
进展修改
(3) 点击切换状态
左键选中线条单击右键删线条进展单元剖分：
(1) 节点编号：编号挨次直接点击该节点假设选中该节点点击 加
点选取半径
(2) 剖分掌握：先选中块四条边〔假设选中该边点击
减点选取
半径者图形放〕然点击
块剖分数填写文件保存：
进展剖分填写终块编号挨次进展
(1) 填写土层渗透系数：点击 输入层土水垂直渗透系数 K K 〔cms〕
x y
留意：输入层土渗透系数点击确定钮进展数保存
(2) 剖分块信息修改：假设剖分信息错误点击编号单击修改钮〕
进展修改〔先选中Reg 某块
(3) 设置相关参数：点击 翻开相关参数设置话框检查分组数等势线百分
数间距渗流状况否正确留意：计算工况稳定渗流计算渗流量时需填写渗流断面穿块号
(4) 保存dqbdat 文件：点击弹出 dqb 文件保存话框选择保存路径点击保存(S)
钮留意：计算时文件名必需dqbdat
(5) 生成断面轮廓图：假设需断面轮廓图直接图中点击鼠标左键选中全部保存
轮廓线点击 弹出保存话框填入文件名保存
42 运行计算
点击 钮弹出计算操作界面(图 42)：
图 42 计算操作界面
(1) 点击查 DQB 钮检查产生 DQB 文件否正确(参考说明书中 43 DQB
DAT 文件容格式)
(2) 点击剖分单元钮生成渗流计算需原始数文件G·DAT 剖分图NET0BLN
(3) 点击显示剖分图钮进展单元剖分时生成网格图NETBLN
(4) 点击编辑数钮翻开计算原始数文件G·DAT留意：该文件末尾起行输入 0 保存
(5) 点击计算钮 生成结果文件SEEP·OUT 渗流场图形文件GS·BLN
(6) 点击等势线图启动 SURFER绘制渗流场分布图：
a单击MAP 菜单中LOAD BASE MAP 消灭 IMPORT FILE 视窗双击BOUND·BLN 单击 OK消灭土坝横断面图
b单击MAP 菜单中LOAD BASE MAP 消灭 IMPORT FILE 视窗双击GS·BLN单击 OK消灭渗流场图形
c土坝横断面图渗流场图形进展叠加：鼠标渗流场图形拖堤防横断面 图鼠标时框选堤防横断面图单渗流场图形单击 MAP 消灭MAP 菜单 单击 MAP 菜单中OVERLAY MAPS 完成图形叠加
d进展文字编辑坐标轴设定等终完善渗流场分布图
(7) 假设想查计算结果SEEP·OUT 点击计算报告钮
43 DQBDAT 文件容格式
产生数文件DQBDAT 容格式
第行：SEEP1 SEEP2 SEEP3
次行：ZUKUAIDIANBIANTUCENGSLDMSLKSDBFS
第三行：KXX〔1：KUAI1：9〕 第四行：ZB〔1：DIAN1：2〕
第五行：SLCS〔1：TUCENG1：3〕 第六行：BJXX〔1：BIAN1：3〕 第七行：SLKH〔1：SLK〕
第八行：SLDMXX〔1：SLDM1：3〕 第九行：SJSW〔1：21：2〕
数说明：
（1） SEEP 1—判否保存渗流绘图数Y’确定N’否认 SEEP 2—判否动剖分数Y’确定N’否认 SEEP 3—设定N’
（2） ZU—分块组数〔计算区域剖分块次左右分假设干组组间应关系第 I 块边〔n3〕第 I+1 块边〔n1〕应边组挨次编号〕
KUAI—四边形块数目
DIAN—块角点数 BIAN—含类边界边数 TUCENG—渗透系数土层种类数
SLDM—需计算渗流量渗流断面数假设零表示计算渗流量 SLK—需计算渗流量断面穿总计块数假设零表示计算渗流量 SD—计算非稳定渗流时水库水位降落时段数假设零表示计算稳定渗流 BFS—计算等势线百分数间距值默认10
（3） KXX(1ME19)块信息数组
块 9 重量块号挨次填写第1—4 重量表示块左角起时针方四角点点号第 5—8 重量次表示块边起时针方4 条边剖分点数第9 重量块渗透系数区号
（4） ZB〔1：DIAN1：2〕块结点坐标数组
块结点序号挨次表示结点XY 坐标
（5） SLCS〔1：TUCENG1：3〕
次层土水渗透系数Kx〔单位：cms〕垂直渗透系数 Ky〔单位：cms〕水度m m 01177 k 〔该式中渗透系数k 单位：mday〕
（6） BJXX(1BN13) 边界信息数组
条边界三重量第重量边界边块块号次重量该边界边块中边序号第三重量边边界类型号规定：
1——渗流面边界〔斜心墙坝仅部面〕
3——半压流边界
5——斜心墙坝游坝壳面边界
6——均质土坝游坝坡排水体渗出段边界
7——游水头边界
8——游水头边界
9——斜心墙土坝斜墙心墙坡渗出段边界
（7） SLKH〔1：SLK〕
需计算渗流量流断面块号
（8） SLDMXX〔1：SLDM1：3〕渗流断面块单元信息
第重量某流断面块数次重量该流断面块边界相交信息 水工渗流计算般 1第三重量该流断面块中位置节点链序号般1 2
（9） SJSW〔1：21：2〕起降时间〔天〕水位降时间〔天〕水位
44 GDAT 数容
软件四边形块单元动剖分三角形单元 时 DQBDAT 动产生
GDATGDAT 软件进展渗流计算原始数文件数容格式
第行：SEEP1SEEP2SEEP3
次行：NENNNNUNNDNNM1NNM2NNFTUCENGAAAA1AA2
BB1BBBB2NNPNQSDNWNGNNM3LULDMA
第三行：H H1H2H3EPSEPS1DTY X
max
AA AA
第十四行：SW〔1：SD〕 第十五行：SJJG〔1：SD〕
第十六行：SLCS〔1：TUCENG1：3〕 第十七行：QTE(1：NQ)
第十八行：QEN(1：Q )
E
解释：
（1） SEEP1SEEP2SEEP3—DQBDAT
（2） NE—三角形单元总数
〔3〕H
—游高水位
max
H1—游水位 H2—游水位
H3—井水位 EPS—面迭代允许误差(01~0001) EPS1—井四周迭代允许误差(05~0005) DT—水位变化总历时〔单位：天〕
（14） SW〔1：SD〕
时段末游水位
（15） SJJG〔1：SD〕 时段时间间隔天数
（16） SLCS〔1：TUCENG1：3〕
次填写层土水渗透系数K 〔单位：mday〕垂直渗透系数K 〔单位：mday〕
x y
水度m m 01177 k 〔式中k 渗透系数单位：mday〕
（17） QTE(1：NQ)
计算渗流量断面单元总数
（18） QEN(1：Q )
E
计算渗流量断面全部单元号
5 算例
例 均质土坝稳定渗流计算
某均质土坝横断面尺寸游水位渗透系数图51 示求计算稳定渗流场单位坝长渗漏流量
图 51 均质坝横断面
解题步骤：
(1) 运行CCEDEXE 文件弹出工程视图相关参数界面分块组数中填入5土层种类数 中填入 2选择稳定渗流选择计算渗流量
(2) 单击视图相关参数坐标中 X 填写 170Y 填写 100坐标中 X 填写30Y 填写 10默认值
(3) 单击确定钮进入图形绘制界面参考41 DQB 文件产生中具体操作步骤建立均质坝横断面剖分图(图 52)
图 52 均质坝横断面剖分图
(4) 保存dqbdat 文件容：
Y Y N
5 8 16 7 2 1 2 0 10
1 4 5 2 3 2 3 2 1
3 6 7 4 2 3 4 4 2
4 7 8 5 4 2 4 2 1
6 9 10 7 4 4 4 3 2
7 10 11 8 4 2 4 2 1
9 12 13 10 8 4 8 4 2
10 13
14 11 8 2 8 2 1
13 15
16 14 4 2 4 2 1
0000
45000
0000
40000
40000 55000
10000 45000
10000
40000
42000
55000
42000
45000
42000
40000
67000
65000
67000
45000
67000
40000
121000
48000
118000
45000
118000
40000
140000
45000
140000
40000
0000002023 0000002023 0047157335
0000010000 0000010000 0059347456
1 1 7
2 4 7
2 1 7
4 1 7
6 1 1
6 2 8
8 1 8
6 7
2 1 2
(5) 点击图形绘制界面中 钮弹出计算操作界面：
a 点击查 DQB 钮检查产生DQB 文件否正确(参考说明书中 43 DQBDAT 文件容格式)
b 点击剖分单元钮生成渗流计算需原始数文件G·DAT 剖分图NET0BLN
c 点击显示剖分图钮进展剖分时生成网格图NETBLN
d 点击编辑数钮翻开计算原始数文件G·DAT留意：该文件末尾起行输入0 保存
e 点击计算钮 生成结果文件SEEP·OUT 渗流场图形文件GS·BLN
f 点击等势线图启动 SURFER观察渗流场分布图：
选择 MAP 菜单中LOAD BASE MAP 消灭 IMPORT FILE 视窗双击BOUND·BLN单击OK 消灭土坝横断面图
选择 MAP 菜单中LOAD BASE MAP 消灭 IMPORT FILE 视窗双击GS·BLN单击 OK 消灭渗流场图形
土坝横断面图渗流场图形进展叠加：鼠标渗流场图形拖堤防横断面图时框选堤防横断面图单渗流场图形选择MAP 菜单中OVERLAY MAPS 完成图形叠加
650
9
0
8
0
7
0
6
0
650
进展文字编辑坐标轴设定等终完善渗流场分布图〔图53〕
图 53 均质坝渗流场分布图
翻开计算报告〔SEEP·OUT〕单位坝长日渗漏流量 QQ100313〔m3daym〕
例二 堤防非稳定渗流计算
堤防堤基横断面尺寸渗透系数图54 示游水位 310m河水位 380m
45 天降 360m 35 天降 340m求计算河水位降落程中堤防渗流场
图 54 堤防横断面
解题步骤参考例建立堤防横断面剖分图〔图55〕
图 55 堤防横断面剖分图
产生问题DQB·DAT 文件容：
Y Y N
3 5 12 7 2 0 0 2 10
1 5 6 2 4 5 4 5 2
3 7 8 4 13 3 13 3 1
4 8 9 5 13 3 15 3 1
5 9 10 6 15 5 15 5 2
9 11 12 10 4 5 4 5 2
0000 30000
0000 20000
30000
38000
20000
34000
10000
30000
10000
20000
53000
33500
58000
31000
60000
30000
60000
20000
70000
30000
70000 20000
0000010000 0000010000 0059347456
0000002023 0000002023 0047157335
1 1 7
3 4 7
2 4 7
2 1 1
2 2 6
3 2 8
5 1 8
0 38000 8000 34000
留意：GSBLN 中包含初始水位 380降瞬时水位 360降瞬时水位 340 3 水位时渗流场绘图数应GSBLN 文件拆分成 3 绘图数文件 GS1BLNGS2BLNGS3BLN分 GS1BLNGS2BLNGS3BLN 堤防横断面图 LUNKE2BLN
叠加然进展文字编辑坐标轴设定等终完善渗流场分布图〔图56〕
40
380
35
310
30
0
01
0 0
6 5
0
4
0
3
2
40
45
50
55
60
65
70
25
0 0 0
9 8 7
20
0 5 10 15 20 25 30 35
i 初始水位 380 时渗流场分布图
3 8 0
90
80
3 1 0
7
0
40
35
30
0
6
0
5
0
4
0
3
0
2
01
40
45
50
55
60
65 70
25
20
0 5 10 15 20 25 30 35
40
380
35
ii
瞬时水位 360 时渗流场分布图
80
30
0 0
25
70
0
310
0
0
0 01
6 6 5 4 3 2
20
0 5 10 15 20 25 30 35 40
45 50
55 60 65 70
iii 瞬时水位 340 时渗流场分布图
图 56 水位降落程中堤防渗流场分布图
例三 心墙坝稳定渗流计算
某心墙土坝横断面尺寸渗透系数游水位图57 示求计算稳定渗流状况土坝渗流场
图 57 心墙坝横断面
解题步骤参考例建立心墙坝横断面剖分图〔图58〕：
图 58 心墙坝横断面剖分图
产生问题DQB·DAT 文件容：
Y Y N
6 12 22 10 3 1 2 0 10
1 5 6 2 4 3 4 3 3
3 7 8 4 2 3 3 3 1
4 8 9 5 3 3 4 3 1
5 9 10 6 4 3 5 3 3
7 11 12 8 2 3 3 3 2
8 12 13 9 3 3 3 3 2
9 13 14 10 3 3 2 3 2
12
15
16
13
5
3
5
3
1
13
16
17
14
5
3
5
3
3
15
18
19
16
3
3
4
3
1
16
19
20
17
4
3
4
3
3
19
21
22
20
4
3
4
3
3
0000 45000
0000 40000
69000 63000
42000 54000
15000
45000
15000
40000
75000
63000
72750
54000
70500
45000
74500
40000
76500
63000
79000
54000
81000
45000
77000
40000
118500
51000
118500
45000
118500
40000
126750
48000
135000
45000
135000
40000
150000
45000
150000
40000
0001000000
0001000000
0114581999
0000020230
0000020230
0065524904
0000500000
0000500000
0103779628
1
1
7
3
4
7
2
4
7
2
1
1
5
1
1
5
2
9
8
1
5
10
1
5
10
2
8
12
1
8
8 9
2 1 2
心墙坝渗流场分布图图 59：
630
0
480
2
0
1
图 59 心墙坝渗流场分布图
文档香网(httpswwwxiangdangnet)户传
《香当网》用户分享的内容，不代表《香当网》观点或立场，请自行判断内容的真实性和可靠性！
该内容是文档的文本内容，更好的格式请下载文档