论文部分内容阅读
在水利水电工程中,常常会遇到有自由面的无压渗流问题,由于该类问题的自由面边界和逸出面边界事先未知,需要在迭代的过程中逐步形成,因此属于边界非线性问题。如何准确、快捷地确定自由面和逸出面的位置历来都是无压渗流分析的一个关键,也是一个难点。目前,有自由面的渗流问题主要是通过有限单元法来求解,其中的固定网格法凭借自身的优势,在无压渗流分析中逐渐取代了变网格法的传统地位,应用日益广泛。本文的加密高斯点单元传导矩阵调整法即是在已有固定网格法——单元传导矩阵调整法的基础上,对高斯点法的一种改进。首先,简单地介绍了渗流的基本原理和渗流计算有限单元方法。接着对变网格法和现有几种有影响的固定网格法,即剩余流量法、单元传导矩阵调整法和初流量法进行了分析。单元传导矩阵调整法以自由面为分界线,自由面以下的饱和区域的渗透系数采用实际数值,自由面以上的非饱和区域的渗透系数则折减一个系数。高斯点法针对有限单元法采用高斯积分公式近似计算单元传导矩阵这一特点,直接根据高斯点处水头与位置高程的关系来调整渗透系数,不用每一次都计算自由面的位置。本文通过加密复合单元的高斯积分点个数,对高斯点法做了改进,得到了加密高斯点单元传导矩阵调整法,使复合单元的单元传导矩阵计算值更趋精确。文中还对逸出面边界进行了处理,实现了逸出面边界节点的可逆转换。最后,采用C++语言编制了稳定渗流分析的有限元程序MESeep,从软件工程的角度对MESeep软件的计划和开发过程进行了介绍。运用MESeep软件对有<WP=3>压、无压,二维、三维,均质、非均质等情况进行了稳定渗流分析,结果表明,加密高斯点单元传导矩阵调整法和根据该法开发的稳定渗流分析有限元软件MESeep能进行各种情况的稳定渗流分析,是可行的。