停泵水锤预测中的直接迭代格式与电算方法.pdfVIP

水 　　利 　　学 　　报 　2001 年 6 月 SHUILI 　　XUEBAO 第 6 期 文章编号 : (2001) 停泵水锤预测中的直接迭代格式与电算方法 1 1 1 2 谢水波 , 邱长军 , 陈胜兵 , 李步新 ( 1. 南华大学建工系 , 湖南衡阳 , 421001 ; 2. 珠江水利委员会设计院 , 广东广州 , 510000) 摘 　要 : 针对水锤数值计算中的特征线方程和边界条件方程 , 设计了一种简单实用的直接迭代格式 , 并编写了 计算软件 , 为停泵水锤数值计算提供了基础. 同时就泵阀计算参数进行了探讨 , 工程实践应用表明 , 其技术上 是可行的. 关键词 : 水锤计算 ; 直接迭代格式 ; 边界条件方程 中图分类号: TV136 + . 2 　　　文献标识码 : A 当水泵机组因突然失电或其它原因停车时 , 水泵及管路中水流速度发生递变而引起压力突变称为 停泵水锤. 停泵水锤危害极大 , 轻则停水 , 重则泵房被淹 , 交通中断甚至人员伤亡. 在水泵站尤其是 远距离输水泵站设计时 , 必须预估水锤的暂态过程 , 以便采取防护措施. 水锤暂态过程的基本方程为 准双曲型偏微分方程 , 结合具体管路的边界条件和初始条件 , 即使在十分简单的边界条件下 , 也不能 积分求解 , 只能求数值近似解. 国内外广泛采用特征线法 , 将运动偏微分方程和连续偏微分方程变换 为四个全微分方程 , 用差分形式来表达 , 借助计算机算出任一瞬时各断面上的水锤水压[1 、2 ] . 由于能 考虑非线性摩阻的影响 , 对复杂管路其优越性尤为突出 , 且精度较高. 本文重点探讨了电算中的迭代 格式和数值解法. 1 　特征线方程 将描述等直径管路中水锤暂态过程的基本方程进行变换 , 可得到其 4 个非线全微分方程[1 ] . 为了 通过有限差分法求其近似解 , 经进一步简化 , 得到近似积分公式 : 沿 c + : H = C - B Q ( 1) Pi p Pi 　　沿 c - : H = C + B Q (2) Pi M Pi 式中 : C 和 C 可根据计算时段开始瞬间的已知数据算出 :[2 ] p M C = H + B Q - RQ | Q |