3、计算机在管网水力计算中的应用-20.ppt

计算机在管网水力计算中的应用 有向线性图的概念 水力计算基本方程组的矩阵表达式 水力计算节点线性逼近解法 管网平差计算框图 管网优化设计模型及框图 （一）有向线性图的概念 在管网平差时，常用一个由点、线和回路构成的计算草图来表达所计算的网络。图中每根线有两个端点，线上标有箭头代表方向．并且假定线的正方向是由始端指向末端。这样的一个草图，图论中叫做有向线性图。 （一）有向线性图的概念 图论中还说明，一个有向线性图的结构可以用二个基本矩阵来表达，一个叫连接矩阵，另一个叫回路矩阵。如图所示为一个简单有向线性图。 1、连接矩阵 一个有向线性图的全连接矩阵就是以节点编号i代表行，枝的编号j代表列，用数值为+1、-1和0的元素所排列起来的一个n×b阶的全矩阵。它的各个元素aij是按照点和枝间的下列关系来确定的： 1、连接矩阵A 2、回路矩阵 B 由图论可知： 这二个矩阵之间还存在着一个十分重要的关系，就是A的行向量和B的行向量的内积总是等于零，这在数学上叫做A阵与B阵之间的正交性。用矩阵代数式来表达，则为： 证明BAT=0 (二)水力计算基本方程组的矩阵表达式 1、满足节点各流量代数和等于零 平差计算的二个基本方程组的矩阵代数式 AQ=q B△P=0 式中： A——A矩阵； B——回路矩阵； Q——管段流量列向量（未知量）； q——节点流量列向量（已知量）； △P——管段压力降列向量。 （三）水力计算节点线性逼近解法 应用线性逼近法解水力计算非线性方程组可以有二个途径 从求解回路流量着手，称回路法： 从求解节点压差（低压管网为节点与基准点间的压力差，中压和高压管网为节点与基准点间的压力平方差）着手，称节点法。 节点法不必提供管段初始流量就能运算，因此用节点法比较方便，现将水力计算的节点法简介如下： 节点线性逼近解法 节点线性逼近解法 （四）管网平差计算框图 （五）管网优化设计模型及框图 所谓管网系统的优化设计，就是对已知管网的网络结构，要求选用合适的管径，在满足节点流量和压力的前提下，使管网的投资和运行费用最小。 为把优化设计问题抽象为一数字规划问题，可把节点流量方程、回路压差方程和节点压力要求作为约束条件，而把管网的投资费用最小作为目标函数，归纳为一个非线性规划问题，即： 管网优化设计框图 * * b=(n-1)+c 式中：b——枝（即管段）的根数； n——点（即节点）的个数； c——独立回路（环）的个数。 Ⅰ 4 3 6 1 2 Ⅱ 5 n=6 c=2 b=6-1+2=7 并假定所示的管段的方向为正。 回路方向顺时针为正。 Ⅰ 4 3 6 1 2 Ⅱ 5 Ⅰ 4 3 6 1 2 Ⅱ 5 5 4 1 2 7 3 节点 管段 可以看出，全连接矩阵的每一列中都含有1和-1二个非零元素，这二个非零元素就是各枝的起终点。有了这个矩阵，就能画出相应的有向线性图。 如果从上面的的全矩阵中任意划去一行，根据每个枝必有两个端点的论点，仍然能够画出同样的图。这说明有一行是多余的。 就是n个节点方程中，其中一个可由其余的n-1个方程相加得到，所以只需要n-1个方程。 对应于被划去的一行所代表的节点在图论中称为“基准点”，为了适应程序上编制的方便，常将“基准点”的编号列在最后（即编号最大）。 从全连接矩阵划去基准点一行后所得的矩阵称它为连接矩阵，简称A矩阵 6 Ⅰ 4 3 6 1 2 Ⅱ 5 节点 管段 它是一个n0×b阶的矩阵，其中n0=n-1。A矩阵除了表达点和枝的连接关系，还能将枝上的量（管段流量）转化为点上的量（节点流量）。这一作用将在下面证明。 5 4 1 2 7 3 6 Ⅰ 4 3 6 1 2 Ⅱ 5 回路矩阵是以行代表有向图中的独立回路ci，列代表枝线ej，是一个c×b 阶的矩阵。它的元素bij是根据回路和枝的关系按照下列定义来确定的： 环路 枝 这个矩阵具有与A矩阵相似的作用，既是反映了枝和回路的关系，又能将枝上的量（管段压力降）转化为回路上的量（回路闭合差）。因此，A矩阵与B矩阵是同等重要的二个基本矩阵。 5 4 1 2 7 3 6 Ⅰ 4 3 6 1 2 Ⅱ 5 5 4 1 2 7 3 6 Ⅰ 4 3 6 1 2 Ⅱ 5 为了使该有向线性图能表达一般管网平差所用的计算草图，选用气源节点作为压力基准点（以下简称基准点），图中的节点6，并在其他节点上引出指向外方的箭头表示节点流量qi，各管段上注有相应的阻耗值Sj（为综合流量公式中各有关参数的常数）及管段流量Qj，这样便得到如图所示的网络计算草图。 5 4 1 2 7 3 6 q5 q4 q3 q1 q2 Ⅰ 4 3