水击基本方程的改善.pdfVIP

水 利 学 报 年 月 第 卷 第 期 文章编号 水击基本方程的改进 杨玲霞 李树慧 侯咏梅 范如琴 郑州大学环水学院 河南郑州 黄河勘测规划设计有限公司 河南郑州 摘要 本文首先指出了当前用于水击计算的数学模型中的连续性方程式在恒定流时不能成立的问题 并通过进一 步的理论推导 得到了新的水击波速计算公式和新的连续性方程 新的水击波速计算公式表明水击波速并非是一 个常量 它的大小随流体动能沿程变化率与压强势能沿程变化率之比的变化而变化 同时也随管道断面相对变化 量与压强相对变化量之比的变化而变化 新的连续性方程改善了老方程存在的不合理现象 本文还通过对水击 现象的计算比较 进一步说明了老方程的确存在不合理现象 而新方程的计算结果是合理的 最后 分析指出了当 前用于水击计算的连续性方程式在恒定流时不能成立的原因所在 关键词 水击压强 水击方程 水击理论 中图分类号 文献标识码 问题的提出 目前几乎在所有的资料和教科书中所介绍的用于水击计算的数学模型都是由意大利学者阿列维 于 年首次以数学方法建立的如下形式的运动方程和连续性方程 运动方程 连续性方程 式中 为管道某断面平均流速 为管道相应断面的测压管水头 为管道轴线与水平面的夹角 为 水击波速 为重力加速度 为管道内径 为沿程阻力系数 水击现象是管道水流从一种恒定状态过渡到另一种恒定状态的非恒定流 描述它的数学模型应该 既能体现水击现象的非恒定特性 也要满足初始和终了时刻的恒定流条件 分析式 和式 可知 运 动方程式 在恒定流时变为 显然 对上式积分后即可得到恒定流的能量方程 因此 运动方程式 在恒定流时自动成立 也即 运动方程式 能够满足恒定流条件 而连续性方程式 在恒定流时则变为如下形式 显然上式不能成立 且当前用于水击计算的基本方程组中的连续性方程不满足恒定流条件 关于 这一点 文献 都提到 但没有解决或没有彻底解决此问题 文献 给出的用于计算均质固液两 相流水击现象的连续性方程在恒定流时也不能成立 收稿日期 作者简介 杨玲霞 女 河南沁阳人 教授 主要从事有压非恒定流理论及计算研究 因此 从理论的正确性和严谨性方面考虑 可以认为当前的水击连续性方程存在着缺陷 这就有必 要进一步分析研究 重新推导出更加完善的用于水击计算的连续性方程 新的水击连续性方程 当前用于水击计算的连续性方程是在一维非恒定流连续性方程的基础上 设法消去变量 管道横 断面积 和 液体密度 引入变量 压强 而得到的 在推导过程中要用到水击压强的计算公式和水 击波速的计算公式 经过分析 认为目前的推导过程不甚严谨 下面进行新的推导 水击压强的计算公式 取图 所示的流段应用动量 定理 经过时段 水击波由下游 断面传到上游 断面 为水击波速 流段内的流速由 增至 压强由 增至 密度由 增至 管道横 断面积由 增至 则动量变化为 作用在 流段上流动方向的外力有 两端断面压力差 图 水击压强计算 重力 摩擦阻力 其中 为管轴线与水平线的夹角 为管壁摩擦切应力 为管道断面湿周 在 时段内外力的冲量为 由动量定理得 取 其中 为管道断面水力半径 为水力坡降 代入 上式得 当前广泛应用的直接水击压强的计算公式为 这相当于在式 中忽略了重力项 阻 力项 以及 和 的影响所得 水击波速