给水管网水力分析和计算.ppt

重新计算各环水头闭合差： 解线性方程组： 得环校正流量解为： ， 施加环校正流量，得到新的管段流量，重新计算有关管段数据，见表5.8。 闭合差减小，但仍不满足要求，需再次构成系数矩阵，并解线性方程组： 解此方程得： 施加该环流量，得到新的管段流量，重新计算有关管段数据，见表5.9。 重新计算各环水头闭合差： 各环水头闭合差满足要求，平差计算结束。 计算管段流速，见表5.9。 由节点（8）出发，用顺推法计算各节点水头和节点自由水压，见表5.10。 （2）哈代-克罗斯算法（1936） 系数矩阵F(0) 是一个对称正定的主对角优势稀疏矩阵，主对角元素值是较大的正值，非主对角的大多数元素为零，不为零的元素都是较小的负值，只保留主对角元素，忽略非主对角元素，则线性方程组可直接求解： (5.30) 此式称为哈代-克罗斯平差公式，哈代-克罗斯算法又称为水头平差法。 哈代-克罗斯算法水力分析的步骤与牛顿-拉夫森算法基本相同，只是计算环校正流量采用水头平差公式（5.30），代替解线性方程组。 【例5.3】多定压节点管网如图5.7所示，节点(1)为清水池，节点水头12.00m，节点(5)为水塔，节点水头为48.00m，各节点地面标高如表5.11，管段［1］上设有泵站，其水力特性如图中所示，计算各管段流量与流速、各节点水头与自由水压（水头损失采用海曾-威廉公式计算，C=110）。 【解】设虚节点（0）及从节点（0）到定压节点（1）和（5）的两条虚管段［10］和［11］，虚管段设有泵站，泵站静扬程为定压节点的节点水头，分别为12.00和48.00m，如图5.7所示。由虚管段［10］和［11］与管段［1］、［2］、［3］和［4］构成一个环，称为虚环，编码为③。 用哈代-克罗斯法进行平差计算，见表5.12。经过两次平差，各环水头闭合差均小于0.5m，最后计算管段流速和节点水头等，见表5.13，计算结果如图5.8所示。 图5.8　多定压节点管网 水力分析结果 5.4　解节点方程水力分析方法 未知量：节点水头（定压节点除外）； 方程：非线性节点流量连续性方程组； 由 得 节点方程：可以列出以 Hi 为未知数的独立方程。 管段流量系数 节点方程组举例： 节点方程通式： 节点方程求解方法： 已知条件： 节点流量、管段直径、摩阻系数、定压点压力； 计算程序： 1）设定未知压力节点（定流点）初始压力Hi(0)； 2）计算系数矩阵A（0）； 3）解线性化节点方程组，得新的节点压力Hi(1)； 4）检验| Hi(1)－ Hi(0)| ε？ 若否，计算系数矩阵A（1），返回3； 若是，计算完成，得Hi； 5）计算管段流量、水头损失等。 节点方程的牛顿-拉夫森解法 由管段i水头损失 微分 变换 管段流量系数 节点方程的牛顿-拉夫森解法（续1） 节点流量函数： 式中，△Hj――为节点的水头增量； Gj――节点流量函数，各定流节点水头增量的非线性函数，只与本节点及相邻环节点水头增量有关。 可以写出： 式中，Gj(0,0,…,0)为给定节点水头初值下的节点流量闭合差： 节点方程的牛顿-拉夫森解法（续2） 由上二式，可得 近似于 G(0)为一系数矩阵， 节点方程的牛顿-拉夫森解法（续3） 定流节点 流量方程组： 牛顿-拉夫森算法计算步骤： 1）拟定定流节点水头初值，给定闭合差的最大允许值， 手工计算时一般取=0.1L/s， 计算机计算时一般取=0.01~0.1L/s； 2）计算各定流节点流量闭合差△Q； 3）判断△Q是否均小于最大允许闭合差， 如满足，计算结束，转（7）进行后续计算， 否则继续下步； 4）计算系数矩阵算； 5）解线性方程组式，得定流节点水头增量； 6）将定流节点水头增量施加到相应节点上，得到新的节点水头，作为新的初值（迭代值），转第2）步重新计算，节点水头迭代计算公式为： 7）计算管段流速、节点自由水压，计算结束。 节点方程的牛顿-拉夫森解法（续4）（节点流量手工平差算法） 节点压力平差算法。 将系数矩阵的全部非主对角元素忽略后，可导出下式： 计算步骤同上述牛顿-拉夫森解法步骤。 详见P108 解节点方程 在假定每一节点水压的条件下，应用连续性方程以及管段压降方程，通过计算求出每一节点的水压。节点水压已知后，即可以从任一管段两端节点的水压差得出该管段的水头损失，进一步从流量和水头损失之间的关系算出管段流量。 原理：在初步拟订压力的基础上，逐步调整节点水压以满足连续性方程。 节点方程组具体步骤： ①根据已知控制点的水压标高，假定各节点的初始水压。 ②由hij=Hi-Hj=sqn的关系式求出管段流量。 ③验证每一