土石坝静动力全面分析.ppt

渗流分析 基本变量 及单元插值 导数 位移 方程右边项 节点力 有限元方程 的物理意义 节点力平衡 应变 应力 水头 水力梯度 流速 节点流量 节点流量平衡 应力变形分析 渗流与应力变形有限元的比较 第二章 土石坝的渗流分析 2.5 稳定渗流有限元计算 流速边界条件的处理 在有限元方程中考虑了流速边界。对不透水边界（vn=0），对方程不产生影响，故在计算中不必特别给出 流速边界： 第二章 土石坝的渗流分析 2.5 稳定渗流有限元计算 水头边界条件的处理 水头边界：对如下4节点的有限元方程，设H2=H0为已知边界条件 第一个方程变为： 在实际的有限元程序中，有限元方程直接集成 第二章 土石坝的渗流分析 2.5 稳定渗流有限元计算 对如下4节点的有限元方程，设H2和H3为已知边界条件 回代到已知水头结点的方程中，可得到所需的结点流量 为了使结点具有所要求的水头值，需在结点上补充或抽出水量，在求解之前未知 求解方程组，得到其余结点的水头值 水头边界条件处理的物理意义 第二章 土石坝的渗流分析 2.5 稳定渗流有限元计算 入流边界 出流边界 i j Qi 1 2 3 4 5 6 1 2 Qi1 Qi2 j 3 4 5 6 Qj3 Qj4 Qj5 Qj6 渗流有限元方程的物理意义 第二章 土石坝的渗流分析 2.5 稳定渗流有限元计算 如何计算透过整个渗透断面的渗漏量？ 入渗节点 出流节点 第二章 土石坝的渗流分析 2.5 稳定渗流有限元计算 自由水面的处理（1） 条件： 自由水面的位置事先无法确定，须在计算中用迭代计算确定，常采用的办法有： 网格修正法 （最原始的办法） 单元传导矩阵修正法 （Bathe,1979) 剩余流量法 （Desai,1976) 初始流量法 （Gell, Wittke, 1984) 改变渗透系数： 见后介绍 第二章 土石坝的渗流分析 2.5 稳定渗流有限元计算 1）对全域进行计算 2）用H=z确定近似自由水面 3）计算近似自由水面的法向流速vn 4）在近似自由水面上叠加-vn，计算其相应节点流量， 修正方程右边项。 5）重新进行计算。并按上述的步骤进行迭代，直到近 似自由水面上近似满vn=0 自由水面处理 – 剩余流量法 h1 h2 vn 第二章 土石坝的渗流分析 2.5 稳定渗流有限元计算 取计算的hz的区域S’为考虑的域，该域不应包含在计算域之内，即不应对其它部分所形成的方程有影响。但其中一单元e的影响为： 自由水面处理 – 初始流量法 第一次迭代时hz的区域 迭代结束后的自由面 迭代过程中的自由面 S’ 影响相当于产生了结点流量 校正： 相当于施加相反的结点流量，以校正影响 第二章 土石坝的渗流分析 2.5 稳定渗流有限元计算 自由水面处理 – 初始流量法 计算高斯点上的水头Hji 求解方程，计算Hi+1 对高斯点循环 修正方程右边项： 对所有单元循环 是 否 1 2 3 4 7 高斯点1-4和7位于渗透区之外 （i=1,2,3,4,7)） 误差比较 第二章 土石坝的渗流分析 2.5 稳定渗流有限元计算 节点2为合适的渗出节点 节点1为不恰当的渗出节点， 需进行修正 1 2 可能的 渗透面 Q20 Q10 1）假定可能的渗透面上的节点全部为渗出点 ? Hi=zi 2）初步进行计算 3）对可能的渗透面上的单元，计算 其对应的节点上的流量值Qi 4）对节点流量值Qi 0的节点，取Hji+1=??Hji (?1, 如?=0.95) 计算由于水头差（ Hji+1-Hji）所产生的方程右边项的修正值 修正方程右边项 7）重新进行计算。并按上述的步骤进行迭代，直到所有的渗出点满足条件。 渗出面的处理 第二章 土石坝的渗流分析 2.5 稳定渗流有限元计算 b) 薄斜心墙坝 a) 厚心墙坝 在进行渗流稳定期坝体有限元的应力变形分析时，如何考虑水压力的作用？给出考虑水压力作用的可能方法 第二章 土石坝的渗流分析 2.5 稳定渗流有限元计算 土中水及渗流计算 非稳定渗流有限元计算 控制方程及边界和初始条件 有限元方程 时间导数项处理 边界条件处理 第二章 土石坝的渗流分析 2.6 非稳定渗流有限元计算 饱和非稳定渗流问题描述 运动方程： 边界条件： 水头边界条件 ?1上给定水头 1 2 流速(流量)边界条件 ?2上给定法向流速(量) 4 3 渗出面 ?3上h=z，vn0 自由水面* ?4上h=z，vn?0 初始条件： 第二章 土石坝的渗流分析 2.6 非稳定渗流有限元计算 在二维