有自由面渗流分析的加密高斯点单元传导矩阵调整法

摘要

该文的加密高斯点单元传导矩阵调整法即是在已有固定网格法——单元传导矩阵调整法的基础上,对高斯点法的一种改进.首先,简单地介绍了渗流的基本原理和渗流计算有限单元方法.接着对变网格法和现有几种有影响的固定网格法,即剩余流量法、单元传导矩阵调整法和初流量法进行了分析.单元传导矩阵调整法以自由面为分界线,自由面以下的饱和区域的渗透系数采用实际数值,自由而以上的非饱和区域的渗透系数则折减一个系数.高斯点法针对有限单元法采用高斯积分公式近似计算单元传导矩阵这一特点,直接根据高斯点处水头与位置高程的关系来调整渗透系数,不用每一次都计算自由面的位置.该文通过加密复合单元的高斯积分点个数,对高斯点法做了改进,得到了加密高斯点单元传导矩阵调整法,使复合单元的单元传导矩阵计算值更趋精确.文中还对逸出面边界进行了处理,实现了逸出面边界节点的可逆转换.最后,采用C++语言编制了稳定渗流分析的有限元程序MESeep,从软件工程的角度对MESeep软件的计划和开发过程进行了介绍.运用MESeep软件对有压、无压,二维、三维,均质、非均质等情况进行了稳定渗流分析,结果表明,加密高斯点单元传导矩阵调整法和根据该法开发的稳定渗流分析有限元软件MESeep能进行各种情况的稳定渗流分析,是可行的.

目录

文摘

英文文摘

1综述

1.1渗流力学的发展

1.2渗流分析的必要性及主要任务

1.3渗流分析方法概述

1.4具有自由面的无压渗流分析方法及其发展趋势

1.5本论文的研究构想

2渗流的基本理论

2.1 概述

2.2渗流基本定律

2.3渗流连续性方程

2.4渗流基本微分方程

2.5渗流基本微分方程的定解条件

2.6渗流的能量泛函表达式

3渗流计算有限单元法

3.1概述

3.2变分原理

3.3单元剖分和插值函数

3.4等参数单元

3.5渗透特性矩阵

3.6方程组的求解

4有自由面渗流分析的现有方法及其评价

4.1确定自由面的传统方法——变网格法

4.2渗流计算固定网格法

5加密高斯点单元传导矩阵调整法

5.1基本原理

5.2自由面的迭代

5.3逸出面的处理

5.4收敛标准

5.5其它水力要素的计算

5.6解题步骤

6稳定渗流分析软件MESeep的开发

6.1 MESeep软件功能的定义

6.2方法的选择

6.3概要设计

6.4详细设计

6.5编码

6.6测式

7算例

7.1有压渗流问题

7.2均质矩形坝

7.3均质土坝

7.4非均质问题

7.5各向异性问题

7.6结果对照分析

8讨论

8.1论文的主要内容

8.2结论

参考文献

作者参与的科研工作、发表的文章及获得的奖励

声明

致谢