第5章-给水管网水力分析和计算.ppt

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,第,5,章给水管网水力分析,5.1,给水管网水力特性分析,5.2,树状管网水力分析,5.3,管网环方程组水力分析和计算,5.4,管网节点方程组水力分析和计算,节点流量方程,（根据质量守恒定律）,含义：流入某一节点的流量等于流出该节点的流量。,注意：,1,管段流量方向（指向节点为负，离开为正）,2,节点流量方向（流入为负，流出为正）,回顾,（,1,）,（,2,）,（,3,）,（,4,）,（,5,）,（,6,）,（,7,）,8,）,1,9,8,7,6,5,4,3,2,Q,7,Q,3,Q,2,Q,1,Q,4,Q,5,Q,6,Q,8,q,1,h,1,q,6,h,6,q,5,h,5,q,2,h,2,q,3,h,3,q,7,h,7,q,8,h,8,q,9,h,9,q,4,h,4,环状网节点流量方程组,回顾,线性变换,求出包括节点（,5,）和（,6,）的大节点连续性方程。,（,1,）,（,2,）,（,3,）,（,4,）,（,5,）,（,6,）,（,7,）,（,8,）,1,9,8,7,6,5,4,3,2,Q,7,Q,3,Q,2,Q,1,Q,4,Q,5,Q,6,Q,8,q,1,h,1,q,6,h,6,q,5,h,5,q,2h,2,q,3,h,3,q,7,h,7,q,8,h,8,q,9,h,9,q,4,h,4,所有节点流量方程相加：,回顾,树状网节点流量方程组,（,1,）,（,2,）,（,3,）,（,4,）,（,5,）,（,6,）,（,7,）,（,8,）,1,7,6,5,4,3,2,Q,7,Q,3,Q,2,Q,1,Q,4,Q,5,Q,6,Q,8,q,1,h,1,q,6,h,6,q,5,h,5,q,2,h,2,q,3,h,3,q,7,h,7,q,4,h,4,回顾,节点流量方程组经过线性变化，可得到：,可以看出：树状网中，各管段流量,q,i,可以用节点流量,Q,j,表示出来。,回顾,管段压降方程,（根据能量守恒定律管段两端节点水头之差等于该管段的压降：,H,Fi,H,Ti,=h,i,i-1,2,M,H,Fi,管段,i,的上端点水头；,H,Ti,管段,i,的下端点水头；,h,i,管段,i,的压降；,M,管段模型中的管段总数。,注意：,判断上下端点时按管段设定的方向，而非实际流向。,M,个管段，可以列出,M,个方程。,回顾,（,1,）,（,2,）,（,3,）,（,4,）,（,5,）,（,6,）,（,7,）,（,8,）,1,9,8,7,6,5,4,3,2,Q,7,Q,3,Q,2,Q,1,Q,4,Q,5,Q,6,Q,8,q,1,h,1,q,6,h,6,q,5,h,5,q,2,h,2,q,3,h,3,q,7,h,7,q,8,h,8,q,9,h,9,q,4,h,4,环状网管段压降方程组,回顾,h,i,可以通过管段的水力特性表示,h,i,=s,i,q,i,n,h,i,=s,i,q,i,|q,i,|,n-1,h,i,=s,i,q,i,|q,i,|,n-1,-h,ei,i=1,2,M,h,i,管段压降，,m,；,q,i,管段流量，,m/s;,s,i,管段阻力系数，应为管段上管道、管件、阀门、泵站所以设施阻力之和；,h,ei,管段扬程，如管段上未设泵站，则,h,ei,=0,；,n,管段阻力指数。,管段流向和设定方向一致，为正，即,s,i,q,i,|q,i,|,n-1,=s,i,q,i,n,回顾,管段压降方程组的变换,如果一些管段首尾相连，形成一条路径，将这些管段的能量守恒方程相加或相减，得到路径能量方程。,例：将管段,1,、,2,、,3,的能量方程相加，再减去管段,4,的能量方程，可导出从节点（,7,）到节点（,8,）之间一条路径的能量方程，即：,H,7,-H,8,=h,1,+h,2,+h,3,-h,4,（,1,）,（,2,）,（,3,）,（,4,）,（,5,）,（,6,）,（,7,）,（,8,）,1,9,8,7,6,5,4,3,2,Q,7,Q,3,Q,2,Q,1,Q,4,Q,5,Q,6,Q,8,q,1,h,1,q,6,h,6,q,5,h,5,q,2,h,2,q,3,h,3,q,7,h,7,q,8,h,8,q,9,h,9,q,4,h,4,回顾,思考题：由节点能量方程导出：,1,）由节点（,1,）、（,2,）、（,3,）、（,6,）、（,5,）、（,4,）组成的回路能量方程。,2,）由节点（,1,）、（,2,）、（,5,）、（,4,）组成的环能量方程。,可以证明，对于任意环状管网，环能量方程的一般形式为：,回顾,恒定流基本方程组,水力分析的数学含义就是解恒定流方程组。,水力分析的工程意义就是已知给水管网部分水力学参数，求其余水力参数。,回顾,5.1,给水管网水力特性分析,（,1,）给水管网水力分析的数学含义：,求解管网恒定流方程组。,（,2,）管网水力分析命题：,在满足供水需求（用水量分布、供水压力和水质）条件下，确定给水管网的科学设计方案（管网布置、管径计算、造价经济、运行安全）。,（,3,）给水管网的水力特性参数：,1,）节点：节点流量、节点压力、节点标高、自由水头；,2,）管段：,管段流量、管径、长度、摩阻系数、管段压降；,3,）环：管网供水保证率、安全可靠性。,管网水力分析条件和目的,:,1,、已知条件：,（,1,）管网布置：枝状管网、环状管网；,（,2,）节点：节点流量、地面标高、服务压力；,（,3,）管段：长度、管径、经济流速、摩阻系数；,2,、管网水力分析求解内容：,（,1,）计算管段流量；,（,2,）计算节点压力；,（,3,）确定水泵流量、扬程；,3,、管网水力分析目的,满足安全供水目标：,（,1,）设计方案水力状态流量、压力分布和变化；,（,2,）管网事故、消防、转输流量工况校核。,5.1.1,管段水力特性,1,、节点方程：求解节点压力,H,i,。方程数,N-1,。,2,、环方程：方程数环数。,3,、管段方程：方程数管段数。,5.1.2,管网恒定流方程组求解条件,（,1,）节点流量与节点水头必须有一个已知,已知节点水头而未知节点流量的节点称为,定压节点,。,已知节点流量而未知节点水头的节点成为,定流节点,。,（,2,）管网中中至少有一个定压节点,5.1.3,管网恒定流方程组求解方法,（,1,）树状管网水力计算,各管段的流量是唯一确定的，与管段流量对应的管段水头损失、管段流速、节点压力可以一次计算完成,（,2,）环状管网水力计算,将节点流量方程组和环能方程组转换成节点压力方程组或环校正流量方程组，通过求解方程组得到环状管网的水力参数。,解环方程组,解节点方程组,解管段方程组,解环方程水力分析方法,解环方程的基本思想,：先进行管段流量初分配，使节点流量连续性条件得到满足，然后，在保持节点流量连续性不被破坏的前提下，通过施加,环校正流量,，设法使各环的能量方程得到满足。,在管网水力计算时，根据求解的未知数是管段流量还是节点水压，可以分为,解环方程、解节点方程和解管段方程三类。,5.2,树状管网水力分析,特点：,（,1,）不存在环方程；,（2,）管段流量,q,i,不变化，管段水头损失,h,i,不,变化，节点方程组系数矩阵元素值为常数，未知节点压力存在直接解。即直接求解线性化节点压力方程组。,枝状管网直接算法,1,、管段流量：采用逆推法。,从树枝末端节点流量开始，用节点流量连续性方程，向前逐一累加，每一管段下游所有节点流量的和即为该管段的管段流量；,2,、节点压力（水头）：采用顺推法。,从已知压力节点出发，用管段能量方程求节点水头，可立即解出。,例,5.1,某城市树状给水管网系统如图所示，节点（,1,）处为水厂清水池，向整个管网供水，管段,1,上设有泵站，其水力特性为：,s,p1,=311.1(,流量单位,m,3,/s,，水头单位,m),h,e1,=42.6m,，,n=1.852,。根据清水池高程设计，节点（,1,）水头为,H,1,=7.8m,，各节点流量、各管段长度与直径如图所示，各节点地面标高见表,5.1,，试进行水力分析，计算各管段流量与流速、各节点水头与自由水头。,节点编号,(1),(2),(3),(4),(5,）,（,6,）,（,7,）,（,8,）,（,9,）,（,10,）,地面标高,(m),9.8,11.5,11.8,15.2,17.4,13.3,12.8,13.7,12.5,15.00,单定压节点树状管网水力分析,【,例,5.1】,计算结果,泵站扬程按水力特性公式计算：,0.89,解环方程水力分析方法,解环方程的基本思想,：先进行管段流量初分配，使节点流量连续性条件得到满足，然后，在保持节点流量连续性不被破坏的前提下，通过施加,环校正流量,，设法使各环的能量方程得到满足。,由于初分流量时是严格按照节点流量平衡来进行的，所以连,续性方程能够满足，但是能量方程就有可能不满足，即环内,正反两个方向的水头损失不相等。环内正反两个方向的水头,损失之差称作,闭合差,。调整管段流量，减少闭合差到一定精,度范围的过程就叫,管网平差,。,5.3,管网环方程水力分析和计算,管网自然环：单一闭合回路。,自然环数：,L,M,N,1,。（,M,L,N,1,）,管网水力环方程：,（,1,）管段方程线性化：,对于管网中管段，给定初始工况点，对式（,5.1,）微分得该点的切线方程：,管段,i,的阻尼系数,（,2,）环方程线性化,环方程转换：未知管段流量 未知环校正流量。,式中，,q,i,(0),管段初始分配流量；,q,k,环校正流量。,环校正流量方程：,环方程线性化泰勒展开式,在环校正流量初值点处 ，,将环水头函数用泰勒公式展开，忽略高次项，取线性项，得线性方程组：,在初步分配流量的基础上，逐步调整管段流量以满足能量方程。,L,个非线形的能量方程：,环方程线性化泰勒展开式,初步分配的流量一般不满足能量方程：,初步分配流量与实际流量的的差额为,q,，实际流量应满足能量方程：,将函数在分配流量上展开，并忽略高阶微量方程组的第一部分称为闭合差：,将闭合差项移到方程组的右边，得到关于流量误差（校正流量）的线性方程组：,环方程矩阵形式,将线性方程组式（,5.19,）表示成矩阵形式：,系数矩阵：,（对角线元素）,（,1,）牛顿,-,拉夫森算法：,1,）拟定管段流量初值，给定闭合差最大允许值，,手工计算时一般取,=0.10.5m,，,计算机计算时一般取,=0.010.1m,；,2,）由式（,5.20,）计算各环水头闭合差；,3,）闭合差均小于最大允许闭合差，则解环方程组结束，,转,7,）进行后续计算；否则继续下步；,4,）计算系数矩阵，式（,5.26,）；,5,）解线性方程组式（,5.21,），得环校正流量；6,）将环校正流量施加到环内所有管段，得到新的管段流量，作为新的初值（迭代值），转第,2,）步重新计算，管段流量迭代计算公式为：,(5.29),7,）计算管段压降、流速，用顺推法求各节点水头，最后计算节点自由水压，计算结束。,（,2,）哈代,-,克罗斯算法 水头平差法,哈代,-,克罗斯（,Hardy Crose,）,1936,年提出，适合于手工计算。,系数矩阵为对称正定、主对角优势稀疏矩阵，只保留主对角元素，忽略非对角元素，直接迭代求解：,哈代,-,克罗斯平差公式,哈代,-,克罗斯法,Q,1,Q,2,Q,3,Q,4,Q,5,Q,6,4,5,6,3,2,1,Q,1,Q,2,Q,3,Q,4,Q,5,Q,6,4,5,6,3,2,1,忽略相邻环校正流量和二阶微量的影响：,校正流量的符号与水头损失闭合差的符号相反,环方程组解法（管网平差）哈代,-,克罗斯法解环方程组步骤：,绘制管网平差计算图，标出各计算管段的长度和各节点的地面标高。,计算比流量、管段流量和节点总流量。,根据城镇供水情况，拟定环状网各管段的水流方向，按每一节点满足,Q,i,+q,ij,=0,的条件，并考虑供水可靠性要求分配流量，得出分配的管段流量,q,ij,(0),。,根据经济流速或查界限流量表选用各管段的管径。,计算各管段水头损失,h,ij,。,假定各环内水流顺时针方向管段中的水头损失为正，逆时针方向管段水头损失为负，计算该环内各管段的水头损失代数和,h,ij,，如,h,ij,0,，其差值为第一次闭合差,h,(0),。,【,如果,h,(0),0,，说明顺时针方向各管段中初步分配的流量多了些，反之，如,h,(0),0,，说明逆时针方向管段中的流量多些。,】,计算各环内各管段的,S,ij,q,ij,其总和,S,ij,q,ij,，按下式求出校正流量。,如闭合差为正，校正流量为负，反之则校正流量为正。,设图上的校正流量符号以,顺时针为正，逆时针方向为负,，凡流向和校正流量方向相同的管段，加上校正流量，否则减去校正流量。据此得第一次校正的管段流量,式中，,q,s,为本环的校正流量；,q,n,为临环的校正流量。,按此流量再行计算，如闭合差尚未达到允许的精度，在从第二步起按每次调整后的流量反复计算，直到每环的闭合差达到要求为止。手工计算，每环闭合差要求小于,0.5m,，大环闭合差小于,1.0m,。,单水源管网平差,泵站,水塔,0,Q,t,Q,p,多水源管网平差,1,2,3,4,8,12,9,11,10,6,5,7,Q,-H,t,-H,p,对置水塔系统,最高用水时,供水分界线,泵站,水塔,0,最大传输时,1,2,3,4,8,12,9,11,10,6,5,7,Q,Q,p,Q,t,-H,p,H,t,对置水塔系统,虚环,P95,应用虚环的观念，将多水源管网转化成为单水源管网。,虚环,将各水源节点与虚节点用虚管段连接成的环虚管段,水源节点与虚节点相连接的管段。虚管段中的流量等于水源,(,定压,),节点的供水量，管段流量方向是从虚节点流向水源节点。虚管段的水头损失等于各水源节点水压，方向是水源节点指向虚节点。虚管段无阻力，虚拟设一泵站，也无阻力。,虚节点水头,H,0,=0,定压节点流量,=0,。,例,5.3 p100,多定压节点管网,5.3.3,哈代,-,克罗斯平差算法的改进,1,）改各环同时平差为每次只平差一个环,每次只平差一个环，平差后立即更新管段流量，后面的平差用新的值计算闭合差，则平差的收敛速度可以得到提高。,2,）优先平差闭合差较大的环,每次只平差一个环时，可以选择闭合差绝对值最大或较大的环进行平差因为这样传递出去的闭合差最大，计算效率最高。,3,）改自然环平差为回路（大环）平差,由多个相邻环合并可以组成回路（大环），同样可以列出回路平差公式，计算效率得以提高。选择合并环的原则，一是闭合差方向相同，二是闭合差绝对值比较大。,5.4,解节点方程水力分析方法,未知量：节点水头（定压节点除外）；,方程：非线性节点流量连续性方程组；,由,得,节点方程：可以列出以,H,i,为未知数的独立方程。,管段流量系数,节点方程组举例：,节点方程通式：,节点方程求解方法：,已知条件：,节点流量、管段直径、摩阻系数、定压点压力；,计算程序：,1,）设定未知压力节点（定流点）初始压力,H,i,(0),；,2,）计算系数矩阵,A,（,0,）,；,3,）解线性化节点方程组，得新的节点压力,H,i,(1),；,4,）检验,|H,i,(1),H,i,(0),|,？,若否，计算系数矩阵,A,（,1,）,，返回,3,；,若是，计算完成，得,H,i,；,5,）计算管段流量、水头损失等。,节点方程的牛顿,-,拉夫森解法,由管段,i,水头损失,微分,变换,管段流量系数,节点方程的牛顿,-,拉夫森解法（续,1,）,节点流量函数：,式中，,H,j,为节点的水头增量；,G,j,节点流量函数，各定流节点水头增量的非线性函数，只与本节点及相邻环节点水头增量有关。,可以写出：,式中，,Gj(0,0,0),为给定节点水头初值下的节点流量闭合差：,节点方程的牛顿,-,拉夫森解法（续,2,）,由上二式，可得,近似于,G,(0),为一系数矩阵，,节点方程的牛顿,-,拉夫森解法（续,3,）,定流节点,流量方程组：,牛顿,-,拉夫森算法计算步骤：,1,）拟定定流节点水头初值，给定闭合差的最大允许值，,手工计算时一般取,=0.1L/s,，,计算机计算时一般取,=0.010.1L/s,；,2,）计算各定流节点流量闭合差,Q,；,3,）判断,Q,是否均小于最大允许闭合差，,如满足，计算结束，转（,7,）进行后续计算，,否则继续下步；,4,）计算系数矩阵算；,5,）解线性方程组式，得定流节点水头增量；,6,）将定流节点水头增量施加到相应节点上，得到新的节点水头，作为新的初值（迭代值），转第,2,）步重新计算，节点水头迭代计算公式为：,7,）计算管段流速、节点自由水压，计算结束。,节点方程的牛顿,-,拉夫森解法（续,4,）（,节点流量手工平差算法）,节点压力平差算法,。,将系数矩阵的全部非主对角元素忽略后，可导出下式：,计算步骤同上述牛顿,-,拉夫森解法步骤。,详见,P108,解节点方程,在假定每一节点水压的条件下，应用连续性方程以及管段压降方程，通过计算求出每一节点的水压。节点水压已知后，即可以从任一管段两端节点的水压差得出该管段的水头损失，进一步从流量和水头损失之间的关系算出管段流量。,原理：在初步拟订压力的基础上，逐步调整节点水压以满足连续性方程。,节点方程组具体步骤：,根据已知控制点的水压标高，假定各节点的初始水压。,由,h,ij,=H,i,-H,j,=sq,n,的关系式求出管段流量。,验证每一节点的管段流量是否满足连续性方程，即进入该节点的流量代数和是否等于零，如不等于零，则按下式求出校正水压值。,除了水压已定的节点外，按校正每一节点的水压，根据新的水压，重复上列步骤计算，直到所有节点的进出流量代数和达到预定的精确度为止。,