调压室水力计算.docx

调压室的水力计算1. 调压室断面计算当上游死水位，下游为最低水位，最小水位=188.9m，三台机满发，引水道糙率取最小值，压力管道糙率取最大值，通过水轮机的流量为57，则此时的引水隧洞水头损失的计算如表格1，压力钢管水头损失的计算如表格2。引水道应选可能的最小糙率0.012，压力管道应选择可能的最大糙率0.013。表格1引水隧洞水头损失表流量（m3/s）沿程水头损失(m)局部水头损失(m)平均糙率(n=0.014)最大糙率(n=0.016)最小糙率(n=0.012)5724.23131.64917.8020.2963810.76914.0667.9110.132192.6923.5161.9780.033表格2压力钢管水头损失表流量（m3/s）沿程水头损失(m)局部水头损失(m)平均糙率(n=0.012)最大糙率(n=0.013)最小糙率(n=0.011)572.6503.1102.2262.805381.4971.7571.2571.962190.7290.8560.6120.824托马断面面积：其中——最小水头损失，；——引水隧洞损失，=17.802+0.296=18.098；——压力管道水头损失，=3.110+2.805=5.915m；L——引水隧洞长度，12662m；g——重力加速度，g=9.81m/——引水隧洞面积，。——引水道阻力系数为了保证大波动的稳定，一般要求调压室断面大于托马斯断面，初步分析时可取（1.01.1），作为调压室的设计断面。这里选取D=7.8m，则系数k为：47.784=1.052. 最高涌波水位计算按正常蓄水位时共用同一调压室的三台机组全部满载运行瞬时丢弃全部负荷（即流量由减至流量）作为设计工况。引水隧洞的糙率取尽可能的最小值（能耗少，涌波高）。n=0.012引水道损失由表格1和表格2得：为时段开始时管中流速；f为引水隧洞断面面积。F为调压井断面面积，145.267；引水隧洞长L=12662m，g=9.81得引水道—调压室系统的特性系数。令，，则要求最高涌波水位，只需要求出即可。的符号在静水位以上为负，以下为正。运用牛顿切线法求解方程的根令牛顿迭代公式为：取迭代初值，计算结果见下表表格3迭代计算结果k1234-0.43085-0.41928-0.41899-0.41899由表格3可以看出精确到0.001，m根据《水电站调压室设计规范》调压室最高涌波水位以上的安全超高不宜小于1m。所以调压室的顶高程：Z=1279+61.128+1.5=1341.63m3. 最低涌波水位计算3.1 增加负荷时的最低涌波水位上游水位取死水位，下游取最低水位。引水隧洞糙率取最大值，取n=0.016。水轮发电机组由2台增至3台满负荷运行。由表1和表2引水隧洞水头损失为：式中，，解得：所以当增加负荷时水位降到静水位下23.74m3.2 丢弃全负荷时产生的第二波动振幅解之47.70m所以当丢弃全负荷时，水位降到静水位下47.70m综上所述调压室的最低涌波水位应为静水位下50.44m根据《水电站调压室设计规范》上游调压室最低涌波水位与调压室处压力引水道顶部之间的安全高度应不小于2～3m。所以压力引水道顶部高程为：Z=1274-50.44-2.5=1221.06m4. 四阶龙格库塔法4.1调压室水位的基本微分方程已知压力管道的流量变化规律，则调压室流量等于隧洞流量减压力管道流量，即：。上述方程两个未知量，Q和Z，采用四阶龙格库塔法进行逐步求解。4.2四阶龙格库塔法计算公式如果已知t时刻的，值，采用以下公式求解t+t时刻的+t，+t。4.3程序框图4.4 VBA程序代码定义函数：Public Function f1(t As Double, z As Double, q As Double)Dim a As Doublea = 5 a为导叶直线关闭时间If t a Thenf1 = (q - (57 - 57 * t / a)) / 47.784Elsef1 = q / 47.784End IfEnd Function定义函数：Public Function f2(t As Double, z As Double, q As Double)Dim r As Doubler = 0.0055 r为隧洞的沿程损失和局部损失系数f2 = 9.81 * 16.61 * (1279 - z - r * Abs(q) * q) / 12662End Function调压室甩负荷时水位波动程序Sub tyssw()Dim t As DoubleDim z As DoubleDim q As Double 将t,z,q定义为实数Dim k1, k2, k3, k4 As DoubleDim l1, l2, l3, l4 As Double 将k