第5章-给水管网水力分析PPT课件.pptVIP

??ΔH ? ??ΔH ? ? ?G4 ??ΔH1 ΔH1 ?ΔQ1 (0)? ?G2 ? ?ΔH5 ?? ? ? ΔH2 ?ΔQ2 ? 0 ??ΔH3? = ??ΔQ3 (0)? ??ΔH ? ?ΔQ4 ? ?? 4? ?ΔQ(0)? ?? ?ΔH5? ? ?ΔH5 ? ?G5 ? ?ΔH5 ? 5.4.1 给水管网节点压力方程组 (续2) 将式（5.37）用泰勒公式展开，忽略高次项，可得节点校正压力线性方程组： (5.40) 1 1 ? ? ?? ? 0 ?G2 ?ΔH3 ?G3 ?ΔH3 0 0 ? 0 ? ? ?G4 ? ?G1 ?ΔH2 ?G2 ?ΔH2 ?G3 ?ΔH2 0 ?G5 ?ΔH2 ?G1 ?ΔH4 0 0 ?G4 ?ΔH4 ?G5 ?ΔH4 ? ?G1 ? ?G2 ? ? 0 ? ? ? ? 0 ? ? (0)? ? (0)? ? 5 ? ? ?∑( ), 系数矩阵的主对角元素; ?k∈j ns q(0) n?1 ? = ?? , 节点j和k衔接，第j行第k列元素及其对称元素; ? nsjk jk (0) n?1 ?0, 节点j和k不衔接，第j行第k列元素及其对称元素; c ?Gj ?∑c(0) = ? ?c(0) ? ? 5.4.1 给水管网节点压力方程组 (续3) 求以下一阶偏微分，可得节点校正压力方程组(5.40)的系数矩阵元素， 1 jk jk q ?Gj ?ΔHk ? 1 ? ? ? (5.44) 令， 1 (0) n?1 nsjk qjk = (0) jk 则有 jk 0 (0) ?ΔHk ?k∈j ? ? (5.45) (5.46) 5.4.1 给水管网节点压力方程组 (续4) 图5.9所示管网的节点校正压力矩阵方程为： (5.47) 从初始节点压力开始，用迭代法求解方 程组(5.47)，可以得到节点校正压力， 使节点流量闭合差收敛到趋近于0的条 件，由此可得各管段流量。 (0) (0) (0) (0) (0) (0) (0) (0) (0) (0) (0) (0) 5.4.2 节点校正压力方程组求解 （1）牛顿-拉夫森算法 牛顿-拉夫森算法直接求解线性化的方程组（5.47），应用迭代算法逐步逼近原方程 组解，其步骤如下： 1）拟定定流节点水头初值，并给定节点流量闭合差的最大允许值eq（手工计算时一 般取eq =0.1L/s，计算机计算时一般取eq =0.01~0.1L/s）； 2）由式（5.39）计算各定流节点流量闭合差； 3）判断各定流节点流量闭合差是否均小于最大允许闭合差，如果满足，则解方程 组结束，转（7），否则继续下步； 4）按式（5.46）计算系数矩阵元素； 5）解线性方程组式（5.47），得定流节点水头增量； 6）将水头增量施加到相应节点上，得到新的节点水头，作为新的初值，转第2）步 重新计算; 7）计算管段流速、节点自由水压，计算结束。 ∑c （2）节点压力平差算法 类似于哈代-克罗斯算法，可以忽略系数矩阵的全部非主对角元素，得到 一种简便的计算方法―节点压力平差法。 忽略系数矩阵的全部非主对角元素后，由式（5.47）可导出下式