三维非稳定饱和非饱和渗流有限元法及其改进.pdfVIP

三维非稳定饱和－非饱和渗流有限元法及其改进 江沆，沈振中，邱乾勇 河海大学水利水电工程学院，江苏南京（210098 ） E-mail ：j 。 摘 要：阐述了非稳定饱和－非饱和渗流的基本方程、定解条件和有限元迭代计算过程。针 对有限元计算过程中经常出现的数值弥散现象，研究了其出现的规律和处理方法，提出了容 水度的改进算法，并编制了有限元程序。最后详细模拟复合土工膜缺陷渗漏量室内物理模型 试验，给出了试验过程中饱和－非饱和渗流场的变化过程。计算表明改进后的计算程序能有 效地消除数值弥散现象，计算成果与试验结果基本一致。 关键词：非稳定饱和－非饱和渗流；容水度；有限元法 饱和－非饱和渗流在计算方法上把饱和－非饱和区作为一个整个流场研究，不需要进行 自由面的调整，自由面仅是负压区正压区的分界面，不再是像饱和渗流场计算时流场的一个 边界，避免了饱和渗流计算中试找自由面和选取给水度补给自由面边界的麻烦。但是，通常 ∂θ k −h 和θ−h 的函数关系是由离散的点线性插值得来，θ−h 的线性插值易造成 的非连 ∂h 续性，在计算中产生振荡，造成非饱和区孔隙负压计算结果分布规律较乱、与实际物理现象 不符的现象，即所谓的数值弥散现象[ 3] [ 4 ] [ 3] 。为消除这一现象，黄康乐 采用变坐标的特征 有限元法，朱学愚等[ 4 ]采用在Galerkin有限元法的权函数中加上一个扰动的SUPG有限元法， 这些方法较好地改善了数值弥散现象，但是均较复杂。吴梦喜[ 5 ]提出了类似式（10）函数 来计算容水度C ，能较好的解决这一问题，但表达式不严密，没有分析表达式中分母等于 零时的情况。这里在文献[ 5 ]中对计算容水度C 的处理方法加以改进。 1. 非稳定饱和－非饱和土渗流定解问题 对于服从达西定律的非均质各向异性饱和－非饱和土，在忽略水的压缩性和孔隙气体对 水流运动的影响情况下渗流基本微分方程为[ 1] ： ∂ ⎡ ∂h ⎤ ∂h s c + s − ⎡ =+ ⎤ c （1） k k h k k h Q C h βS ⎢ ( ) ( )⎥ ( ) ij r c i3 r c ⎣ c s ⎦ ∂xi ⎢⎣ ∂x j ⎥⎦ ∂t 式中：h 为总水头；h 为压力水头；k s 为饱和渗透系数张量；k 为相对透水率，是非 c ij r 饱和土的渗透系数与同一种土饱和时的渗透系数的比值，在非饱和区0 k 1 ，在饱和区 r ∂θ C β k = 1 ；C 为容水度，θ为介质含水率，在正压力区 = 0