第十四章 水电站的水锤8.docVIP

第十四章 水电站的水锤 第九节 水电站有压引水系统非恒定流电算法简介 ? ? 水电站有压引水系统的非恒定流计算包括水锤计算和调压室涌波计算。这两种计算各有特点而又相互联系。在负荷变化时，机组的转速变化与水锤和调压室涌波也有联系。把这三种过渡过程联系起来研究的理论虽然早已基本具备，但由于计算过于繁琐，在电子计算机应用于工程实际之前，很少有把它们联系求解的实例，一般都是用孤立的、简化的方法计算。即使对于分岔管的水锤，为了避免繁琐的计算，也往往采用很粗略的简化方法。电子计算机的应用给较精确、合理地计算上述问题开避了新途径，现简要介绍如下。? ? 一、简单管水锤计算? ? 简单管水锤计算一般不必利用电子计算机。但如欲在计算中考虑水头损失或机组特性的影响，用电子计算机能较好地处理这类问题，用一般的方法则难以解决。? ? 对于图14-12的简单管，若水锤波通过管段AP和PB的时间均为，则求解A、P、B三点压强和流量的方程为? ? A点：式中h=H/Ho=1+ζ；q=Q/Qo=v；α=/Ho，扩为对应于Qo的AP段水管的水头损失。? ? 式(a)可写成式中， 。? ? 将式(b)代人式（c），并令S=一A，得? ? 解式(d)，舍去增根，得? ? 故式(a)、式(b)可写成? ? P点：? ? 式（e）和式(f)可写成? ? 解式(g)和式(h)，得? ? B点：? ? 上式可写成? ? 根据式（14-76）、式（14-77）和式（14-78）所列的顺序，不难编出简单的程序迭代地求出A, P, B三点压强和流量的变化过程，计算可以从t=2,开始，到所要求的时刻为止。根据开度变化曲线确定。? ? 二、分岔管的水锤计算 图 14-25 分岔管示意图 ? ? 对于图14-25所示的分岔管，若P点有n个分支，水锤波通过各分支的时间均为,通过主管PB的时间为m,,参照式(14-76)，水轮机端，，……、An点的水锤压强和流量可用以下式组求出式中i=1，2,…，n,故以上式组共n个。? ? P点的压强和流量有n＋2个未知量，用下列n十2个方程求解特征方程? ? 即式中，，i=1，2,……，n,故上式有n个。? ? P-B特征方程? ? 即式中；。? ? P点连续方程? ? 以上共有n+2个方程，其中具有下标t+者为未知量，亦有n+2个。将前n+1个方程写成矩阵形式，为了避免系数矩阵主对角线的积为零，将第n+1个方程排在第二位，并以的关系代人，得? ? 式（18-80)可用高斯(Guass)消元法求解。关于高斯消元法的电算程序在一般的电算程序书中均可找到，可以借用。有的计算机中存有这种程序，只要调用即可。? ? 参照式(18-78)，B点在t+ (m+1)时刻的流量可用以下式组求出? ? 计算可从t = 2开始，到要求的时刻结束。、、……、根据、、……、An点的开度变化曲线确定。? ? 程序框图如下： ? ? 程序框图中，NH为预定的总时段，I为正在进行计算的时段数。? ? 三、具有调压室的引水系统非恒定流计算? ? 可有以下三种情况：? ? (1)对细而高的调压室（如差动调压室的升管等），可将调压室视为一个支管，把水锤．计算和调压室的涌波计算用水锤方程统一起来，调压室水面处的流速即为调压室的水位变化速度。? ? (2）若调压室比较粗短，为了避免因计算时段过小而大大增加计算工作量，可以忽略调压室和其中水体的弹性变形。在这种假定下，调压室底部的流量等于其顶部的流量，底部的压强等于调压室水柱的高度和其中水体惯性力的代数和。? ? （3)若调压室极粗短，且与压力水管的连接处断面又较大（如无连接管的简单调压室），在水锤计算过程中调压室水位变化不大，则可视调压室为水库。? ? 四、水电站的调节保证计算? ? 用计算机可以把水锤、调压室涌波和机组的转速变化三种过渡过程联系起来加以研究。其方法是在以上所述的基础上引人机组的动力方程，并根据水轮机的特性和开度变化规律，从而求出压强、水位和转速的变化过程，限于篇幅，具体方法从略。? ? 五、简单管水锤的计算程序及计算实例? ? 1.程序中符号的意义? ? 程序根据式(14-76)、式(14-77)和式(14-78)所列的顺序，用FORTRAN语言编制。程序中变量的意义如下：? ? HO—水头Ho；C-水锤波速c；D-管道直径D；N-管道糙率n；? ? QMAX-管道中最大流量；TS-阀门关闭时间Ts；? ? TMAX-水锤计算总时间；? ? VMAX管中最大流速Vmax=Qmax/ (/4)；? ? HW-管道中的水头损失；TR—水锤波通过AP和PB管段的时间；? ? P-水锤常数ρ=c