解节点方程管网平差程序的开发与应用.doc

1解节点方程管网平差程序的开发与应用摘要：给水管网力计算是以解管段方程、解环方程和解节点方程为基础，对连续性方程、能量方程和压降方程应用近似优化处理方法和数值计算方法进行计算，旨在求解管段流量或节点水压，为管网设计，改扩建及运行管理提供依据 关键词：节点方程 管网平差 开发与应用 1 引言 给水管网力计算是以解管段方程、解环方程和解节点方程为基础，对连续性方程、能量方程和压降方程应用近似优化处理方法和数值计算方法进行计算，旨在求解管段流量或节点水压，为管网设计，改扩建及运行管理提供依据随着供水事业的发展，给水管网的规模不断增大，管段数和环数不断增多众所周知，传统的解环方程法是在手算基础上发展而成的，计算前需要初分管段流量对于大型复杂管网，初分流量相当繁琐，人工工作量较大，且初分值不合理会导致迭代算法不收敛为此，本文基于解节点方程的算法原理及管网数据结构的特征，研究了正定稀疏矩隈的变带宽紧缩贮存技术，运用 FORTRAN 语言编制了程序，并结合实例进行了应用和验算 2 解节点方程的有理与方法 2.1 节点方程2根据管段压降方程，Hi-Hj=Sijq2ij,将管段流量用水压表示，qij=sign(Hi-Hj)(│Hi-Hj│÷Sij）1/2，代入连续性方程，即得出节点方程，Q+Σsign(Hi-Hj)(│Hi-Hj│÷Sij）1/2=0式中 Qi—i 节点的耗水量或水源供水量（即节点流量）；HiHj——i,j 节点的水压；Sij——i,j 管段的摩阻。

若管网节点数为 M，则独立的节点方程数为 M-12.2 节点方程的线性化节点方程是以节点压力未知量的非线性方程组，令 Cij=1/(Sij│qij│),qij 的初值可用程序中所示的经验公式确定，则节点方程可化为，Q+ΣCij(Hi-Hj)=0，这是一个线性方程组，可用迭代法或牛顿法求解，程序中采用的迭代法2.3 线性方程系数矩阵的存贮根据管网图形拓扑结构可知，以上线性方程的系数矩阵为对称正定稀疏矩阵，矩阵元素中大部分为 0，节点数越多，稀疏性越明显对于 M 个节点的管网，矩阵元素共（M-1）2 个，按一般矩阵存贮需要（M-1）2 个存贮单元对称矩阵只需要存贮一半元素（上三角或下三角矩阵）即可对于稀疏矩阵，依照一定次序用一维数组紧缩存贮每行的第一个非零元素到对角线上的元素，再用指标数组存放各对角线元素在一维数组中的位置序号，这种变带宽紧缩方式可以进一步有效地节省存贮单元2.4 节点方程的计算步骤⑴读取数据，按照经验公式计算初分流量，初定管径，计算摩阻；⑵计算初始系数矩阵参数；⑶解线性方程组，求节点点压，利用压差计算管段流量，高速管径及摩阻返回；⑶重新生成系数矩阵；⑷迭代至前后两次管段流量之差在允许精度范3围内；⑸进一步计算节点自由水压，管段流速，水头损失等；⑹输出计算结果。

3 解节点方程程序的应用 解节点方程的 FORTRAN 源程序及说明从略应用程序前，需绘制计算简图，按要求将节点、管段编号，将基础数据输入文件 input.dat 中，结果文件 output.dat 中节点编号原则：已知压力节点编号；未知坟力节点编号尽可能与相邻节点编号差值小，以利于紧缩存贮某城区给水管网最大用量 822L/s，由两个泵站和水塔联合供水城区地形平坦，地面标高均按 0 米计节点要求的最小服务水头为 24m如图 1其它参数见文件 input.dat 输入文件 input.dat 的格式为： 15 21 3 0.01 24.0036.20 0.0036.80 0.0082.50 0.00436.40 0.0048.70 0.0081.50 0.00198.70 0.0066.10 0.0050.60 0.0043.20 0.005105.80 0.0035.50 0.0034.65 0.0030.36 0.0027.40 0.001270.00 0.50 0.013 2 11350.00 0.30 0.013 3 2650.0060.50 0.013 4 3620.00 0.60 0.013 5 11150.00 0.40 0.013 6 21390.00 0.40 0.013 7 31670.00 0.40 0.013 8 4760.00 0.50 0.013 6 51130.00 0.30 0.013 7 61040.00 0.30 0.01378 71730.00 0.40 0.013 9 55480.00 0.30 0.013 10 61140.00 0.20 0.013 11 71510.00 0.20 0.013 12 81500.00 0.30 0.013 10 91020.00 0.30 0.013 11 10760.00 0.20 0.013 12 118150.00 0.40 0.013 15 12225.00 0.50 0.013 13 1225.00 0.50 0.013 13 1240.00 0.50 0.013 14 4结果文件 output.dat 的输出格式为：DATA OF NODES F（1）= 33.726 Z(1 )=33.726F（2）= 31.772 Z(2 )=31.772F（3）= 28.276 Z(3 )=28.276F（4）= 29.1629Z(4 )=29.162F（5）= 32.246 Z(5 )=32.246F（6）= 30.749 Z(6 )=30.749F（7）= 24.658 Z(7 )=24.658F（8）= 25.976 Z(8 )=25.976F（9）= 29.455 Z(9)=29.455F（10）= 27.695 Z(10)=27.695F（11）= 22.779 Z(11)=22.779F（12）= 27.231 Z(12)=27.231F（13）= 34.651 Z(13)=34.651F（14）= 30.361 Z(14)=30.361F（15）= 27.40010Z(15)=27.400DATA OF PIPES Q（1）=148.292 V(1)=.756 h(1)=1.955 i(1)=1.539Q（2）=49.308 V(2)=.698 h(2)=3.496 i(2)=2.590Q（3）=-139.608 V(3)=.712 h(3)=-.887 i(3)=1.364Q（4）=300.235 V(4)=1.062 h(4)=1.480 i(4)=2.387Q（5）=62.184 V(5)=.495 h(5)=1.023 i(5)=.889Q（6）=106.415 V(6)=.847 h(6)=3.618 i(6)=2.603Q（7）=91.113 V(7)=.726 h(7)=3.187 i(7)=1.909Q（8）=167.753 V(8)=.855 h(8)=1.49711i(8)=1.970Q（9）=71.145 V(9)=1.007 h(9)=6.092 i(9)=5.391Q（10）=-34.494 V(10)=.488 h(10)=-1.318 i(10)=1.268Q（11）=83.784 V(11)=.667 h(11)=2.792 i(11)=1.614Q（12）=77.292 V(12)1.094 h(12)=3.055 i(12)=6.363Q（13）=13.353 V(13)=.426 h(13)=1.880 i(13)=1.649Q（14）=-9.482 V(14)=.302 h(14)=-1.256 i(14)=.832Q（15）=33.184 V(15)=.470 h(15)=1.760 i(15)=1.173Q（16）=67.276 V(16)=.952 h(16)=4.917 i(16)=4.821Q（17）=-25.172 V(17)=.80212h(17)=-4.453 i(17)=5.859Q（18）=-70.153 V(18)=.559 h(18)=-.170 i(18)=1.132Q（19）=-242.364 V(19)=1.235 h(19)=-.925 i(19)=4.111Q（20）=-242.364 V(20)=1.235 h(20)=-.925 i(20)=4.111Q（21）=-267.119 V(21)=1.361 h(21)=-1.199 i(21)=4.993 4 结论 本文编制的管网平差程序，运行速度快，处理信息量大，不用人工初分流量，输入简单，容错性好，这些已经多个算例得以验证。

但同时必须指出，该程序仍需要在大型复杂管网水力计算的应用中不断改进和提高 参考文献1 严煦世，范瑾初. 给水工程（第三版）. 北京：中国建筑工业出版社，1995.2 彭永臻，崔福义 . 给水排水工程计算机程序设计. 北京：中国建筑工业出版社，19943 谭浩强，田淑清 . FORTRAN77 结构化语言. 北京：清华大学出版社，19904 严煦世，赵洪宾. 给水管网理论和计算. 北京：中国建筑工业出版社，19865 吴学伟 给水管网水力计算程序的编制和应用.管道技术与设备，2000.4:31-34 。