计算水力学05.pptVIP

第五章 河网水流计算 第一节 河网的基本概念 第二节 树状河网水流计算 一、树状河网计算方法 二、树状河网编程技巧 第三节 环状河网水流计算 一、外河道计算 二、内河道计算 三、节点水位方程 四、求解步骤 第四节 最优编码解法 二、编程技巧 第五节 矩阵标识法 二、编程技巧 李光炽 计算水力学 由此可见，系数矩阵A有如下关系 上式表明，系数矩阵A为主对角占优矩阵。 根据矩阵A所具有的特性，可采用迭代法求解。最常用的迭代法之一是超松弛法。 李光炽 计算水力学 则AZ=R的SOR法可写成 可见，如果能避免所有零元素运算，即可提高求解运算的效率。对非零元素随机分布的稀疏矩阵A，按行将A的非零元素依次排成一个元素序列{b}={b1,b2,…,bt}，存放于一维系数数组B［1:t］中。对应的元素列标排成序列{Dr}={Dr1,Dr2,…,Drt},存放于标识代码数组Dr［1:t］中。 李光炽 计算水力学 同时把每一行的第一个非零元素在B中的序号排成系列{Ri}={Ri1,Ri2,…,Rin}，存放于行代码指示数组Ri［1:n］中，其中t为矩阵A的非零元素的总个数，ｎ为矩阵Ａ的阶。 利用矩阵标识代码数组Dｒ，行代码指示数组Ri和系数数组Ｂ，可将求解式写成 李光炽 计算水力学 由此可见，如果对系数矩阵的非零元素进行代码标识，完全可以只对非零元素进行运算，避免所有零元素运算，从而高效率求解节点水位方程。 李光炽 计算水力学 (22) (3) (16) (17) (21) (20) (19) (18) (2) (1) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) (15) 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 图5-5 环状河网计算示意图 李光炽 计算水力学 外河道的求解是,通过边界条件确定首断面的递推关系,用追赶的方法求解各断面的水位与流量递推关系,得到汇入基本河网的末断面流量与内节点水位的关系。 （１）对于水位型边界条件有如下追赶方程 李光炽 计算水力学 （２）对于流量型边界条件有如下追赶方程 无论哪一类型边界条件的外河道，末断面的流量QL2都可表达成QL2=f(ZL2)=f(Z末) Z末为末节点水位。由河网联解求出节点水位ZL2，回代即可求得Zi、Qi。 李光炽 计算水力学 设河道的首断面号为L1,末断面号为L2,有如下差分方程： 有2(L2-L1+1)个未知量，2(L2-L1)个方程，方程的个数总比未知量个数少两个。因此，以首、末断面水位为基本未知量，可利用双追赶方程求解,具体如下： 李光炽 计算水力学 令 其中i=L2-2,L2-3,…,L1。 李光炽 计算水力学 对于i=L2-1有 李光炽 计算水力学 令： 其中i=L1+2,L2+3,…,L2。 李光炽 计算水力学 对于i=L1+1有 李光炽 计算水力学 因此，由上述递推公式可以得到  其中ZL1=Z首为首节点水位,ZL2=Z末为末节点水位,即首、末断面流量表达为首、末节点水位的线性组合。式(5—27)和(5—28)称为环状河网的河道追赶方程。与单一河道或树状河网河道的追赶方程的形式不同，每个河段具有6个需要保存的追赶系数。当首、末断面水位求得后，利用追赶方程，对同一断面上的流量有： 李光炽 计算水力学 联立求解得  求得Ｚi后代入追赶方程式即可求得Ｑi。  所以关键问题是如何求解节点水位，下面介绍 求解节点水位的方法。 李光炽 计算水力学 为了求解河网节点水位,必须建立节点水位方程。节点水位方程建立的依据是水量守衡原理,即流进某一节点的水量之和等于该节点蓄水量的变化。 李光炽 计算水力学 把流量与节点水位关系式代入得到与节点i相邻的节点水位为未知变量的线性代数方程：  其中Zｉ,ｊ为与节点ｉ相邻节点水位的集合。对河网每一个节点,都可建立上述的节点水位方程,形成以河网节点水位为基本未知变量的线性代数方程组： 其中:  李光炽 计算水力学 以图５—5节点（8）为例,同该节点有直接联系的河道有6、7、17、18四条，这四条河道在节点（8）处的入流方程为： 节点水量平衡条件为： 李光炽 计算水力学 将式代入得： 可见方程中的节点水位只包含与该节点有直接联系的节点水位。 李光炽 计算水力学 首先,确定基本河网。通常为了降低节点水位方程组的阶数,仅考虑内河道和内节点组成的基