理正岩土使用手册-渗流分析

2.2.4 编辑原始数据
录入或选择渗流分析所需的各种原始数据，有限元法和公式法交互窗口分别如图2.2-4和2.2-5。
图2.2-4 有限元数据交互对话框
图2.2-5 公式法数据交互对话框 注意：
1. 集中的参数交互界面，即把几乎所有的参数置于一个界面上，操作简单，大大提高了人机交互的效 率，这是理正岩土系列软件的一个共性特征。
最后得：
9． 已 知 水 头 边 界 的 处 理 方 法 处理方法1为类似于强制位移的处理方法，乘大数法，此方法的优点是简单，缺点是容易造成病态矩
阵。。 处理方法2是将已知水头的节点都移到方程的右端，方程的未知变量中不出现这些节点，优点是方程不
会出现病态，未知数个数比第一种减少，缺点是处理起来比较复杂。 本系统中采用第二种方法。
5.2 一般稳定渗流计算方法
一般稳定渗流计算方法适用于以下情况： 土堤类型：
均质土堤 心墙式土堤 斜墙式土堤 地基类型： 透水地基 不透水地基 排水条件： 无排水设备或贴坡式排水 褥垫式排水 棱体排水
5.2.1 均质土堤
5.2.1.1 不透水地基 1．不透水地基均质土堤下游坡无排水设备或贴坡式排水稳定渗流计算
2.1 操作流程
第二章 快速操作指南
图2.1-1 操作流程
2.2 快速操作指南
2.2.1 选择工作路径
I:/…/06渗流分析（一二三）.htm
1/10
2010-5-26
06渗流分析（一二三）
图2.2-1 指定工作路径 注 意 ： 此处指定的工作路径是所有岩土模块的工作路径。进入某一计算模块后，还可以通过按钮【选 工程】重新指定此模块的工作路径。
图5.2.1-1 不透水堤基均质土堤的渗流计算（无排水设备） I:/工作/…/06渗流分析（五）.htm
（5.2.1.1-1）
（5.2.1.1-2） （5.2.1.1-3） （5.2.1.1-4）
1/12
2010-5-26
06渗流分析（五）
式中： q —— 单位宽度渗流量（m3/d.m）； k —— 堤身渗透系数（m/d）；
7/10
2010-5-26
06渗流分析（一二三）
3.3.3 图形查询辅助工具
1． 图 形 查 询 工 具 栏
图3.3-2 计算结果图形查询界面
2． 图 形 查 询 快 捷 菜 单 在图形结果查询窗口单击鼠标右键，弹出图形查询快捷菜单，可以方便地查看图形。
3． 【 图 形 查 询 】 菜 单
3.4 计算书的编辑修改
3.1.4 读入理正“边坡稳定”数据文件
用于读入理正边坡稳定分析软件的数据文件，与边坡稳定共享坡线、地层信息。当点此菜单命令时，将 会弹出一个对话框选择所要读入的文件名称，再点【打开】按钮则自动读入。
3.1.5 读入DXF文件自动形成坡面、节点、土层数据（复杂土层模块中）
提供了另一种交互复杂土层的方法，对于熟练掌握AutoCAD的用户，能更快捷完成复杂土层的交互。通 过【辅助功能】菜单中的“读入DXF文件自动形成坡面、节点、土层数据”把土层交互到土层参数表中。选 择命令后，弹出注意事项对话框（图3.1-1）：
公式法：
“计算简图” “计算结果图”
06渗流分析（一二三）
通过在输出选项列表中选择相应的选项来输出对应的图形结果。
图3.3-1 图形输出界面
3.3.2 通过辅助功能菜单查看图形结果
单击【辅助功能】菜单（图3.1-3）中的【查看计算图形结果】项，可查看当前例题的图形结果。
I:/…/06渗流分析（一二三）.htm
2010-5-26
06渗流分析（一二三）
第一章 功能概述
渗流分析计算软件主要分析土体中的渗流问题。适用于勘察、设计等单位进行土堤、土坝的渗流分析、 闸坝地基的渗流分析、堤防的渗流分析、基坑降水的流场分析等。并可以将流场的数据传递到稳定分析软 件，以便分析考虑流场的稳定问题。
⑴ 渗流的分析方法：公式方法和有限元方法。 ⑵ 公式方法依据《堤防工程设计规范》提供的计算公式。适用于下列情况： 一般稳定渗流计算； 双层地基稳定渗流计算； 水位上升过程中不稳定渗流计算； 水位降落过程中不稳定渗流计算。 ⑶ 有限元方法是依据非饱和土理论、根据基本的渗流理论――达西定律等，采用有限元方法分析稳定 流及非稳定流中多种边界条件、多种材料的堤坝、或土体的渗流分析。但有限元法分析渗流问题是以线性达 西定律为基础，因此不适应非线性达西定律的流场分析及不满足达西定律的流场分析。
达到该值但仍未满足计算要求时，程序也会停止计算，这时会提示计算不成功。应修改参数或边界条件后， 再重新计算。
2. 计算不收敛的可能原因 1）渗透系数的数量级相差太大－－将渗透系数较大的地层去掉，不参与分析； 2）土层区域的构造不合理－－不要形成有较小锐角的区域； 3）网格剖分单元长度不合理－－有狭窄区域时，应特别注意； 4）边界条件交互不合理－－注意边界条件是否与实际情况相吻合； 5）交互的计算迭代次数不够! 由于导致计算不成功的可能原因并不唯一，因此需要用户认真检查交互的参数，分析不成功的原因，修 改相应的参数，重新计算。
06渗流分析（一二三）
图3.1-2 读入DXF文件注意事项对话框
3.2 计算简图辅助操作菜单
在数据交互界面的左侧图形窗口单击鼠标右键，弹出图形显示快捷菜单，使用该菜单可有效的查看计算 简图，可把计算简图存为DXF格式的文件，用AUTOCAD等图形编辑器进行编辑。
3.3 快速查询图形结果
3.3.1 选择输出图形结果
2.2.2 计算项目选择
选择渗流计算所采用的方法（有限元分析法与公式法）：
图2.2-2 计算项目选择
2.2.3 增加计算项目
点击【工程操作】菜单中的【增加项目】菜单或“增”按钮来新增一个计算项目。
I:/…/06渗流分析（一二三）.htm
2/10
2010-5-26
06渗流分析（一二三） 图2.2-3 增加计算项目界面
2. 同时提供了有关参数的即时弹跳说明信息，方便用户理解参数的意义。
2.2.5 计算结果查询
I:/…/06渗流分析（一二三）.htm
3/10
2010-5-26
06渗流分析（一二三）
图2.2-6 计算结果查询窗口 计算结果查询界面分为左右两个窗口，左侧窗口用于查询图形结果，右侧窗口用于查询文字结果。
注 意 ：当采用不同的有 限 网 格 剖 分 长 度 两次计算的结果非常接近时，可认为此时的解是合理的解。
最 大 迭 代 次 数 ——指程序计算时迭代次数的上限值。 迭 代 收 敛 判 断 误 差 ——迭代计算时，两次计算结果误差的允许值。
注意： 1. 当程序计算过程中已经满足计算要求时，迭代次数即使未达到该值也会自动停止迭代；而迭代次数
10/Biblioteka Baidu0
2010-5-26
06渗流分析（四）
第四章 有限元法编制原理
4.1 编制依据
⑴《渗流计算原理及应用》顾慰慈 编著 中国建材工业出版社 2000.08； ⑵《土工原理与计算》（第二版）河海大学 钱家欢 殷宗泽 主编 中国水电出版社 1996.05； ⑶《渗流数值计算与程序应用》毛昶熙 段祥宝 李祖贻 主编 河海大学出版社 1999.01； ⑷《有限单元法原理与应用》（第二版）朱伯芳 编著 中国水利水电出版社1998.10。
图3.1-1 读入DXF文件注意事项对话框
3.1.6 生成准流网图
通过【辅助功能】菜单中的“生成准流网图”生成流网图，并通过“间隔数”的交互及“加删等势 线”和“加删流径”按钮命令编辑“等势线”和“流径”的间距。选择命令后弹出对话框（图3.1-2）
I:/…/06渗流分析（一二三）.htm
5/10
2010-5-26
4.3 非稳定流有限元分析方法
非稳定渗流问题的有限元方程是与时间有关的方程，对时间进行有限差分运算，可得到如下方程：
式中： Δt —— 时间增量； ω —— 0 ~ 1之间的一个系数； H1 —— 时间增量结束时的水头值； H0 —— 时间增量开始时的水头值； Q1 —— 时间增量结束时的流量值； Q0 —— 时间增量开始时的流量值； 其它符号的意义同前。
其中：Ω为单元内源或蒸发造成的流量m3/s，内源为正，蒸发为负。 8． 边 界 流 量 所 引 起 的 力 向 量
其中：Vn=lxVx+lyVy为单元边界的流速m/s，流入为正，流出为负。
I:/工作/…/06渗流分析（四）.htm
2/4
2010-5-26
06渗流分析（四）
而lx=c os α，ly=c os β，如上图显示。假定Vn在边界上为线性变化，在1点为q1，在2点为q2，将上述积分 运算展开：
4.2 有限元法基本计算公式
1． 渗 流 有 限 元 分 析 基 本 方 程
式中： [K] —— 透水系数矩阵； {H} —— 总水头向量； [M] —— 单元储水量矩阵； {Q} —— 流量向量； t —— 时间。
2． 透 水 系 数 矩 阵 [K]
其中：
其中：α── 最大透水系数的方向角，用于各项异性的材料中。 ３．由节点水头求单元内任意点水头
取ω=0时，上述方程变为：
根据上式，如果已知前一个时间步的水头值，就可以根据上式递。
I:/工作/…/06渗流分析（四）.htm
3/4
2010-5-26
06渗流分析（四）
I:/工作/…/06渗流分析（四）.htm
4/4
2010-5-26
06渗流分析（五）
第五章 公式法编制原理
5.1 编制依据
《堤防工程设计规范》 GB 50286-98
有限元法：
“计算简图” “压力水头各点计算值” “压力水头等值线” “压力水头彩色云图” “总水头各点计算值” “总水头等值线”
I:/…/06渗流分析（一二三）.htm
6/10
2010-5-26
“总水头彩色云图” “流速矢量图+等水头线” “流速Vx彩色云图” “流速Vy彩色云图” “流速V彩色云图” “X向水头梯度计算值” “X向水头梯度彩色云图” “Y向水头梯度计算值” “Y向水头梯度彩色云图” “给定界面数值曲线”
第三章 操作说明 3.1 关于计算例题的编辑
3.1.1 增加例题与删除当前例题
1．通过【工程操作】菜单的“增加项目”和“删除当前项目”来增加一个新的例题或删除当前的例 题。
2．“增”或“删”按钮增加一个新的例题或删除当前的例题。通过“算”计算当前选择的例题。
3.1.2 数据的读写
通过【辅助功能】菜单的“读入数据文件”可以将原来保存好的数据读进来进行计算；通过【辅助功 能】菜单的“数据存盘到文件”可以将当前例题的数据保存在磁盘上。
文字结果输出较为完整的计算书，主要包括以下信息：
“计算简图” “计算已知条件”
I:/…/06渗流分析（一二三）.htm
8/10
2010-5-26
“计算结果”
06渗流分析（一二三）
文字结果可以使用文字编辑菜单进行编辑，也可以用其它文本编辑器进行编辑。
3.5 几个参数的说明
有 限 网 格 剖 分 长 度 ——在计算时划分的单元长度，该参数影响计算的精度及速度。一般可以先 取1m～5m之间的数进行计算，查看计算结果，不满意可修改该值重新计算。
I:/…/06渗流分析（一二三）.htm
4/10
2010-5-26
06渗流分析（一二三）
3.1.3 把典型例题加入例题模板库
实际工程中会有一些具有一般代表性的典型例题，当完成该例题的数据交互后，可通过【辅助功能】菜 单中的“将此例题加入模板库”把该例题存为例题模板，从而在每次新增例题时可以重复调用该例题的数 据，在此基础上修改少量的数据进行计算。
H1 —— 上游水位（m）； H2 —— 下游水位（m）； h0 —— 下游出逸点高度（m）； m1 —— 上游坡坡率； m2 —— 下游坡坡率；
b —— 土堤顶部宽度（m）； H —— 土堤高度（m）； L —— （m）； ΔL —— （m）； L1 —— （m）。
其中：
I:/工作/…/06渗流分析（四）.htm
1/4
2010-5-26
06渗流分析（四）
其中：A为三角形面积 4． 由 节 点 水 头 求 单 元 内 任 意 点 流 速
也就是：
注 意 ： 单元内任意点的流速均相等。 5． 单 元 刚 度 矩 阵 [K]
6． 单 元 储 水 量 矩 阵 [M]
7． 内 源 或 蒸 发 所 引 起 的 力 向 量
3.6 数据和结果文件
数据和结果文件位于用户设定好的工作目录下。有限元分析法和公式法的数据文件格式分别
I:/…/06渗流分析（一二三）.htm
9/10
2010-5-26
06渗流分析（一二三）
为*.SL和*.GSL，图形文件格式*.DXF，计算书格式为*.RTF。
I:/…/06渗流分析（一二三）.htm