会员注册 | 登录 | 微信快捷登录 QQ登录 微博登录 | 帮助中心 人人文库renrendoc.com美如初恋!
站内搜索 百度文库

热门搜索: 直缝焊接机 矿井提升机 循环球式转向器图纸 机器人手爪发展史 管道机器人dwg 动平衡试验台设计

   首页 人人文库网 > 资源分类 > DOC文档下载

水利工程论文-图论法用于供水管网水力计算的研究.doc

  • 资源星级:
  • 资源大小:51.32KB   全文页数:5页
  • 资源格式: DOC        下载权限:注册会员
您还没有登陆,请先登录。登陆后即可下载此文档。
  合作网站登录: 微信快捷登录 支付宝快捷登录   QQ登录   微博登录
友情提示
2:本站资源不支持迅雷下载,请使用浏览器直接下载(不支持QQ浏览器)
3:本站资源下载后的文档和图纸-无水印,预览文档经过压缩,下载后原文更清晰   

水利工程论文-图论法用于供水管网水力计算的研究.doc

水利工程论文图论法用于供水管网水力计算的研究摘要图论理论是网络分析的主要工具,现用于管网的水力平衡计算,既充分发挥了图论理论的优势,使计算变得简便、迅捷,又可将管网附件加入计算,使结果更准确、更符合实际。文中采用峰阵输入管网结构,使输入数据的工作量大大减少,易于编制程序,计算大型的复杂管网。关键词供水管网水力计算图论法前言供水管网的水力平衡计算是供水系统规划设计、经济评价和运行管理的基础。水力平衡计算的目的就是在确定管径的情况下求出满足连续方程和能量方程的各节点压力水头和各管段流量。目前常用的水力平衡计算方法有哈代克罗斯法HardyCross,牛顿莱福逊法NewtonRaphson,线性理论法LinearTheory,有限元法FiniteElement等等。所有这些方法各有所长,适用范围各不相同,有的还需人工假设管段流量,使输入数据工作量增大,且未考虑管网附件的影响。本文介绍的图论法将复杂的管网处理为相应的网络图,并建立相应的数学模型,用峰阵输入原始数据来描述管网结构,输入的数据量最少,不易出错,易于计算大型的复杂管网。其计算过程可同时考虑管网附件,如控制阀、加压泵、逆止阀、减压阀等,使计算结果更符合实际。1图论原理将供水管网中的管段概化成一条线段即图中的边,将有附件的管段看成图中的特殊管段,边与边由节点相连。这样,一个供水系统的管网图就转化为图论中的网络图。而且管道中的水流是有方向的,所以管网图是有向图。根据以上所述原则,可将图1所示管网系统,转化为图2所示的网络图。图1图2图1中有一水库A,三个给水点B、C、D,Q1表示水库节点供水量,Q2\,Q3\,Q4分别表示B、C、D节点的用水量。管段视为网络图中的对应边,管段的直径、管长、管道流量、摩损系数等作为管段对应边的权。至此,与管网同构的网络图生成了。图中箭头表示各条边的方向,即管段中水流方向。网络图中节点与边的关联函数可以用完全关联矩阵I45表示如式1所示。顶点边的编号1式中Iij{1,表示j管段与i节点相连,且管内水流流离该节点0,表示此管段不与该节点关联1,表示j管段与i节点相连,且管内水流流入该节点。完全关联矩阵与管段流量列向量q以及节点流量列向量Q可组成管网节点方程即连续方程IijqQ0,qq1,q2,q3,q4,q5T,QQ1,Q2,Q3,Q4T。网络图的生成树全涉及树可以有很多种,在计算时可以任选一种。在本例中,选1、2、4这3条边为图的生成树,则补树余树的各边弦为3、5.各弦将与枝构成基本回路,一个基本回路中有且仅有1条弦。用基本回路矩阵Bf表示则如式2所示。枝124弦35Bf11010211101式中每一行表示一个基本回路环。环的方向以该环对应弦的方向为准。1表示管段中的流向与环中弦的方向相反,1表示相同,0表示该管段不在此环内。Bf可用矩阵B和单位阵U表示为式3。Bf[B|U],其中B1103111环阵与管段摩损列向量hf构成环方程如式4所示。摩损向量的元素顺序与Bf中每行元素所对应的管段顺序相同。Bfhf0。其中hfh1,h2,h4,h3,h5T4图论理论中,连续方程用割方程代替。每个割方程只含一根枝,并和相关的弦构成割集,将图2分割成互不连通的脱离体。这样,图中就有3个割集。割集和割集阵Af如式5所示割集K割阵Af枝124弦355K1e1,e3,e510011K2e2,e3,e501011K3e4,e500101割阵Af中,每一行表示一个割集。图中有3根枝,所以就有3个割集。割阵中,1表示该管段在此割集内,且管段流向与此割集内的枝中的流向相同,1表示流向相反,0表示该管段不在此割集内。式5的割阵Af和割集K一一对应。割阵Af可用一个矩阵A和一个单位阵U表示为Af[U|A],其中A[11]1101割阵与流量列向量可构成割方程。根据图论理论,割阵的行向量与环阵的行向量正交,这种关系可用式6表示。B|UU|AT0或者U|AB|UT06所以有BAT或者ABT。这样,环阵可以由割阵求出,反之亦然。关联矩阵通过选主元初等行变换即可得到割阵先选关联阵第一行中一非零枝元素为主元,并使其为1,消去其它各行中此主元再选第二行、第三行、的主元,最后即得割阵Af。因此,可以由关联矩阵导出割阵和环阵。2图论法模型任何管道的水力计算都可以用管段流量q\,水头损失h\,管径D\,管长L和管壁条件C等5个因素来描述。一般D、L和C为已知条件,只有q和h未知。因此,求解一个管网的水力平衡问题,可从两方面考虑一是利用q和h的关系,消去h,以q为未知量计算,求出q后,反求h二是首先消去q,以h为未知量计算解出h之后,再反求。图论法也可从这两方面入手,即求弦流量式和求枝摩损式。前者只适用于环状网,而后者则适用于所有类型的管网,所以本文着重介绍后者。设一管网有J个节点,P条管段,L个环,则三者满足LPJ1的关系。管网的每一管段都有q和h两个未知量,因而未知量的个数为2P。但管网环方程有L个,线性无关的连续性方程有J1个,总数为LJ1P个,不能求解2P个未知量[1]。因此,必须借助P个管段摩损方程式。管段摩损方程式线性化后的通式如7和8所示。系数R称为阻尼系数,Y称为传导系数。R和Y的具体形式与所选用的摩损公式有关,是D、C、L的函数。摩损公式线性化后,R还是q的函数,Y还是h的函数。不过,在求解过程中,总是把R和Y当作已知量来对待。阻尼式hRq7传导式qYh8式中R和Y是阻尼系数和传导系数矩阵。如果摩损公式采用HazenWilliam公式,则有hRq10.68q1.852L/C1.852D4.8710.68L|q|0.852/C1.852D4.87q9R10.68L|q|0.852/C1.852D4.8710Y1/RC1.852D4.87/10.68L|q|0.852C1.852D4.87/10.68L|q|0.85211用h向量表示管段摩损h表示枝摩损,h′表示弦摩损用q向量表示管段流量q枝管段流量,q′表示弦管段流量。割方程的右端项Q为脱离体所含节点流量之和。方环程Bfh0,即BU[h]012h′割方程AfqQ,即UA][q]Q13q′传导式[q][Y0][h]14q′0Y′h′求枝摩损式以管段摩损为未知量首先将传导式14代入割方程13得[UA][Y0][h]Q150Y′h′由环方程12可得Bhh′0,即h′Bh,代入式15得[UA][Y0][h]Q160Y′Bh即h[YAY′B]Q17根据正交定理得h[YAY′AT]Q18这就是图论法的求枝摩损式计算公式。h即为枝管段的摩损向量。解得枝摩损值h后,其余变量可由相应的公式求出。由环方程可得h′Bh,即可求出弦摩损向量h′,q、q′向量可以由式14求得。式11中C1.852D4.87/10.68L对某一管段来说是个常数,可用W表示。则传导系数Y可以表示为YW|q|0.85219在迭代计算时,第一次可以直接用W代替Y进行计算,求出h\,q后计算Y,再求新的q值,如此反复计算,直至前后两次的q值符合给定的误差标准为止。为了避免可能出现的数值摆动现象,在第三次迭代时,用前两次迭代结果的流量平均值作为初始流量值[2],即qq(1)q2220求得q3,,这样收敛速度加快。3管网附件实际管网中,有许多控制、安全、量测设施,如加压泵、控制阀、逆止阀、减压阀等附件,对管网运行产生重要影响。传统计算方法都未涉及到管网附件问题,不仅使计算准确性受损,而且其计算程序无法用于日常管理工作。图论法处理管网附件时,将附件所在管段视为特殊管段,这些管段的摩损式要根据其附件的水力学特征计算摩损值,再加入到管网中进行水力平衡计算。本文给出几种较常见管网附件的处理方法。对于其它附件,具体问题具体处理,在此就不一一详述了。3.1普通阀门闸板式阀门是用得最多的一种阀门,在一般的水力计算过程中,闸板式阀门的水头损失计算一般引用公式hfξv2/2g,ξ值见文献3。其中,a表示管段中过水断面的高度,d表示管段直径,a/d表示阀门开关。当开度为0时,阀门完全关闭,没有流量通过当开度为1时,阀门完全打开,对水流不产生影响。将阀门水头损失公式用流量表示为hfξv2/2gξ2q2/π2gD2则阻尼系数R为R2ξq/π2gD2传导系数为Yπ2gD2/2ξq1计算时只需将闸板式阀门的R或Y值加入,即可计算。蝶阀的计算方法与闸板式阀门类似3.2逆止阀逆止阀是管网中最常见的设备之一,是水流方向控制设备,只允许水流单向通过。

注意事项

本文(水利工程论文-图论法用于供水管网水力计算的研究.doc)为本站会员(奋斗不息)主动上传,人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知人人文库网([email protected]),我们立即给予删除!

温馨提示:如果因为网速或其他原因下载失败请重新下载,重复下载不扣分。

copyright@ 2015-2017 人人文库网网站版权所有
苏ICP备12009002号-5