计算渗流量一种简易方法.pdfVIP

计算渗流量的一种简易方法 蔡晓英1 朱军2 (1．广东省水利水电科学研究院，广州5106lo； 2．广东省水利厅，广州510150) 摘要 渗流有限元计算渗流量的常规方法是通过对法向流速在过流面上进行积分得到 通过该断面上的流量。本文提出一种不同于常规方法的渗流量计算简易方法，利用计算水头 值及渗流传导矩阵，简单、高效地求出渗流场的过流量。与常规方法相比，该方法具有简单且 计算精度高的特点。通过对两个算例的比较分析证明了此方法的正确性。 关键词渗流量，计算，简易方法 1 前言 渗流量是水工建筑物渗流分析中的一个重要指标，用它可以表示建筑物的漏水情况，为防 渗、排水设施提供可靠的设计依据。因此高效、正确地计算渗流量具有很重要的实际意义。 在渗流有限元计算中，一般先求出渗流场的水头值，再求出过流断面上的法向流速，法向 流速在过流面上的积分值即为通过该断面上的流量。 本文不按常规方法，仅利用计算出的水头值及渗流传导矩阵，便可简单、高效地求出渗流 场的过流量。 2渗流有限元及常规计算流量方法 有限元求解稳定渗流场的方程为 [K][H]=[Q] 式中，[K]为总体渗透矩阵；[Q]为等效结点流量项；[H]为待求的结点水头向量。 对于无压渗流场，则须对上式进行多次迭代，以得出真实的渗流区域，具体方法有很多，见 文献[1]、[2]。对于有压渗流场，渗流区域已知，无须迭代，直接求解方程(1)便可。算出水头 值后，利用其计算流速，由于流速在单元面上不连续，故常规方法的计算断面一般由若干单元 中截面组成(如图1)。由结点水头值求出单元中截面上各点的法向流速，再在中截面上对其 积分求单元过流量Ql，断面所穿过单元中截面过流量之和即为计算断面流量Q。 图l某土坝剖面 图2 计算断面右边的一个单元 CoS n工 r矗， ， CoSn S ， 了 d (2) Q=∑Q=∑』卜ds=∑『f[裂，并，蓦]{是， CoS扎 ： 、●，●，●。L，●●J (，＼( ， ))) 法向向量。 由(2)式可知，为求流量，必须先求出若干单元中截面上各点的法向流速，再将法向流速乘 以中断面面积。这种求解方式要求先计算法向流速，且其计算流速时要对插值函数求导一次， 这就意味插值函数的多项式将降低一次，所得出的计算结果精度也随之降低。 如果能避开法向流速的求解，而仅利用已计算出来的水头值及渗透矩阵加上较简单的运 算彳辱出流量值，将会简化计算，并得出精度较高的结果。 3 求解流量的简易方法 我们可以由若干单元面所组成的面为计算断面，该断面可以看作是由在其同边(以图1为 例，为右边或左边，不同边得出的过流量绝对值相等但正负号相反)若干单元的单元面组成，只 要求出这些单元面上的结点流量，将它们相加然后取绝对值，即为该断面的过流量。我们以图 l中1—1断面右边的若干单元面为计算对象来说明本方法。 对于式(1)中的总体渗透矩阵(真实渗流区域)，由各单元渗透矩阵叠加而成，即[K]= ∑[K]i。 对于计算断面右边的一个单元j，如图2，用下式可计算该单元对i结点(位于计算断面上) 的结点流量贡献， 4