单源环型管道水力计算的递推算法.pdfVIP

第 20卷 第 l期 2007年 3月 苏州科技学院学报(工程技术版) J．of University of Science and Technology of Suzhou (Engineering and Technology) Vo1．20 No．1 Mar． 2o07 单源环型管道水力计算的递推算法 周 虹， 刘 敏， 高洁羽 (苏州科技学院 计算中心，江苏 苏州 215011) 摘 要：针对单源环型管道这种特定的网络结构情况，提出了一种水力计算的递推算法，可直接计算出环型管道的 流量分布。文章阐述了递推算法的数学依据 ，证明了递推算法的可行性，并给出了算法的源程序代码。 关键词：水力计算；递推算法；单源环型管道 中图法分类号 ：TP 301 文献标识码 ：A 文章编号 ：1672—0679(2007)01-0068—04 1 问题的提出 管道化运输是传输流体介质的最经济方便的方法。利用管道边界的流量信息和水力学模型，通过管道 的水力计算，可以描述出管道内部各个管段的流量、压力等信息，从而为管道的系统设计及运行管理提供依 据。目前，最常用的水力学模型是采用公式 h．= ·g 表示管段的水头与流量的关系川，辅以管道中各结点的流 量平衡关系，形成管道水力计算的数学模型。由于这是以管段水头、流量为变量的一组非线性方程和线性方 程混合的方程组，没有普遍适用的直接求解方法。通过线性化和迭代的方法得到近似解，是国际上通用的做 法【2。l，这种方法通用性强、模型概括性好，也是目前各种水力计算软件采用的基本算法。但由于迭代算法需 要指定一组较为合适的初始值，且由于迭代收敛速度的不确定，造成计算过程耗时的不确定性，所以并不适 合应用于实时监控 、仿真等实时性要求较高的场合，在这些场合，通常根据具体情况使用一些近似算法，如 文献【4】、【5】中的方法。 在实际应用情况下，管道网络的拓扑结构是多样化的，当管道网络是枝叉型结构(即内部无环)时，结点 个数P与管段个数 之间的关系为 P=L+I，由于存在 尸-1个线性无关的结点流量平衡方程，通过直接求解 这 个线性方程组，即可获得每个管段中的流量，进而可依次计算出每个管段的水头，并不需要迭代算法Il1。 当管道网络是环型管道时(如图 1所示)，结点个数P与管段个数 之间的关系为L=P，用环平衡方程与尸-1 个线性无关的结点流量平衡方程联立，就形成了环型管道水力计算的数学模型。环平衡方程是基于环型管 道中两个方向的管段总水头相同的事实，若管段的水头与流量采用公式 】-： ·g ，则环平衡方程的基本形式 为∑【k·q ]=0，其中方扩号中的符号视累加方向与流量 g 的方向是否相同而定。 在该模型中，流量的变化会引起环平衡方程中累加集合的变化，造成平衡方程是一个分段二次方程，所 以模型的解也应是一个分段函数，很难有一个解的通用解析表达式。笔 者针对这种模型，给出了一个解的递推算法，使用该算法可直接计算出 单源环型管道中的流量分布。 2 单源环型管道的计算方法 设源点为尸0，环中其它结点顺时针排列分别为 尸l，?，尸n，相应结点 的流出量分别为Q 一，Q ，每个管段的阻力系数分别为 ，? 。设管 段内流量变量分别为g。，g 一，g ，沿顺时针方向流动用正值表示，逆时 针方向流动用负值表示。则该环的水力平衡方程组共有 n+1个方程，其 中有 n个流量平衡方程(线性方程)和 1个能量平衡方程(2次方程)。 [收稿日期】2006—10—27 【基金项目】苏州科技学院基金项 目 【作者简介】周 虹(1958-)，女，云南海通人，副教授。 @ 图 1 单源环形供水管道 @ 维普资讯 ：单源环型管道水力{ qi=q~-l-Q。i=1，2，?，n ={Vi，qi10l (1) (2) 2．1 递推公式 注意到 Q。的非负性(单源)和q 的递减性，所以存在整数 k使得q~O(i-O，?，k)和 qi≤0( =JI}+l，?，n)， 即n={0，l，?，k}，将公式(1)变形为公式(3)和(4) g。 ∑ =l，2，?，七一l (3) qi=qk+l- ∑QJ =k+2，k+3，?，n j=k+2 将公式(3)和(4)分别代人公式(2)的第 1项和第2项，则有公式(5)和(6) ∑ ．g ∑ ．g ：∑ ．( ∑QJ) =∑ · 2 2 q iEn i=0 i-9 J +1 ‘-0 QJ) ：． ∑ ．g =∑ · +。 itn 1 +1 2‘q +l ， =2‘ ∑ ∑ 』 +2 ’ ~f／+2 ．(∑ +∑ i=k+2 J 1 QJ．∑ i-O ·(q 一 )2 ∑ ．，．2 q + J= 。‘ 一‘I= +l ∑ ．(∑ i=O J=l 一 1 )