基于特征线法的水锤计算程序的开发.docxVIP

基于特征线法的水锤计算程序的开发 停止泵的水锤是水锤的一种。这意味着，当水泵突然断开或其他原因发生时，由于水流突然变化，一些严重的压力，如水泵和管道系统中的站版快速增长，会导致频繁上升的水流损失。一般情况下停泵水锤最为严重, 其对泵房和管路的安全有极大的威胁, 国内有几座水泵房曾发生停泵水锤而导致泵房淹没或管路破裂的重大事故。 停泵水锤值的大小与泵房中水泵和输水管路的具体情况有关。在泵房和输水管路设计时应考虑可能发生的水锤情况, 并采取相应的防范措施避免水锤的发生, 或将水锤的影响控制在允许范围内。我院在综合国内外关于水锤的最新科研成果并结合多年工程实践的经验, 以特征线法为基础开发了水锤计算程序。这一程序可较好地模拟各种工况条件下水泵及输水管路系统的水锤状况, 为高扬程长距离输水工程提供设计依据。 1 原理和计算方法 停泵水锤的计算有多种方法:图解法、数解法和电算法。其基本原理是按照弹性水柱理论, 建立水锤过程的运动方程和连续方程, 这两个方程是双曲线族偏微分方程。 运动方程式为: ue014H/ue014x+1/g×ue014V/ue014t+V/g×ue014V/ue014x+ f/D×V2/2g=0 (1) 连续方程式为: ue014H/ue014t+Vue014H/ue014x+a2/g×ue014V/ue014x=0 (2) 式中H——管中某点的水头 V——管内流速 a——水锤波传播速度 x——管路中某点坐标 g——重力加速度 t——时间 f——管路摩阻系数 D——管径 通过简化求解得到水锤分析计算的最重要的基础方程: H-H0=F(t-x/a) +F(t+x/a) (3) V-V0=g/a×F(t-x/a) -g/a× F(t+x/a) (4) 式中F(t-x/a) ——直接波 F(t+x/a) ——反射波 在波动学中, 直接波和反射波的传播在坐标轴 (H, V) 中的表现形式为射线, 即特征线。它表示管路中某两点处在水锤过程中各自相应时刻的水头H与流速V之间的相互关系。为了方便计算机的计算, 将上述方程组变换为水头平衡方程和转速改变方程, 即成事故停泵时水泵的两个边界条件方程式: F1=PM-BQv+Hn(β2+v2) (A0+ A1x) -ΔH0v2/ (τ2) =0 (5) F2= (β2+v2) (B0+B1x) +m0-C3 (β0-β) =0 (6) 式中β——N/Nn(实际转速/额定转速) v——Q/Qn(实际流量/额定流量) 通过上述两式的联立, 采用牛顿—莱福生迭代公式, 可以解出v和β的近似数值。 将水泵的全面性能曲线改造为仅与转速和流速有关的全面性能曲线, 以便计算机在解方程时取值, 即: WH(x) =h/ (β2+v2) (7) WM(x) =m/ (β2+v2) (8) 式中h——H/Hn(实际扬程/额定扬程) m——M/Mn(实际转矩/额定转矩) 现行的水锤计算方法就是基于上述原理。 2 水泵机组基本情况 根据我国南方某城市取水泵房的水泵及输水管线的实际情况, 采用计算机程序模拟水锤情况如下 (均按最低枯水位计算) : 基本情况: 水泵机组:Qn=5 000 m3/h,Hn=55 m,Nn=741 r/min,Ns=132.4,GD2=874.7 kg·m2,Mn=932.72 kg·m, 近期单台运行, 远期两台运行。 输水管线:DN=1 400 mm,L=5 750 m, 几何扬程:35 m (近期) , 45 m (远期) 。 泵房和输水管线如图1所示: 2.1 假设采用军阀管停止泵的水锤 实际阀门关闭时具有实际流速 普通止回阀管路停泵水锤计算结果如表1所示。 表1中所列数据为假设水泵出口处的流速为零时阀门即刻关闭所产生的水锤压力值。实际工况中, 阀门关闭总要一段时间, 因此实际水锤值将与表中所列数据有出入。根据计算机模拟结果, 如果在此条件下适当增加水泵机组的转动惯量可以将水锤压力值明显降低。 单泵慢组合关闭 缓闭止回阀管路停泵水锤计算结果如表2所示。 经过计算机模拟, 当关阀时间和快慢组合与最佳模拟条件不同时, 泵前最大压力值都将有所增加。因此一个装有两阶段关闭阀门的输水系统, 其阀门的操作过程应经过计算确定, 并应在试运行中调整。此种设备定货时应向制造厂提出具体的技术要求 (快、慢关闭时间及可调性) 。 实际计算机模拟结果 此泵房出水管在穿越大堤处 (距泵出口40 m) 形成了驼峰, 经计算, 此处将发生弥合水锤。实际观测与计算机模拟的结果相近, 计算机模拟结果如表3所示。 从表3得知, 当管路中发生断流的停泵水锤 (即弥合水锤) 时, 水锤值很大, 达到几何