急流弯道双曲底板的体形设计.doc

水 利 学 报 2003 年08 月 SHUILI XUEBAO 第8 期 文章编号：0559-9350（2003）08-0001-05 急流弯道双曲底板的体形设计 黄春花1，孙颖2，陈肇和2，冯新权 3 (1.中国水利水电科学研究院 水力学研究所，北京 100038；2.北京工业大学分部，北京 100044； 山东省水利勘测设计研究院，山东 济南 250013) 摘要：本文提出采用双曲底板的方法消除急流渠道弯段的冲击波。在柱面座标系中建立不可压缩粘性流体的运动方 程，导出无量纲形式的自由水面方程并建立了相应连续性方程。通过迭代法求解自由水面方程并计算出自由水面高 程，由连续性方程求出水流深度。从水面高程中减去水深得出底板高程。文中给出了按此方法设计的马河水库溢洪 道弯道段水流计算结果。结果证明渠道内水面平稳，没有发生冲击波。底板设计由编程计算完成。 关键词：双曲底板；弯道；急流；冲击波 中图分类号：TV133 文献标识码：A 在水利工程建设中，常因地形地质条件所限，其泄水建筑物布置成弯道形式 [1]。当急流进入弯 道段时发生冲击波。采用双曲底板可以削减急流弯道内的冲击波，有效调整弯道内的水流状态。 1 基本原理 双曲底板是指不仅具有纵向曲率，而且还具有横向曲率的双向弯曲的渠底。在弯道段采用双曲底板， 可以削减冲击波，以使自由水面趋于平缓。在这类结构中，水流呈三维 流动，水流中的压强并不服从静水压强分布规律。因此，对具有双向弯 曲底板结构的水力计算，既不能用一维水力学的方法，也不能用波高较 小的波浪迭加法[2]。计算这类结构的基本原理是从流体力学的基本理论 出发，建立三维弯道水流的数学模型，根据来流的水力参数、流态要求， 推求能实现对水流控制的底板形式。这类问题的实质是求解流体力学的 逆命题，即由给定水力要素反求固体边界。 基本方程 自由水面方程 在柱面坐标系(图1所示)里，建立不可压缩粘性 图1 坐标系s、y、z 液体的运动方程，并取其坐标轴分别为流线在水平面上的投影(曲线s)、 水平轴y和铅垂轴z，可以得到无量纲形式的自由水面方程式如下[3]： z = ? ? z C n ? zk ∫ + ? A? B ? ? ? ? + ? ? 1 s 2 z / s2 D z / s 1 k k k + ? F ( z / s ) + ? ? ∏ ? k ? ? ? ? ? ? dn (1) 收稿日期：2003-05-16 作者简介：黄春花(1976-)，女，安徽怀宁人，硕士研究生，主要研究方向为水工水力学与冷却塔。 1  水 利 学 报 2003 年08 月 SHUILI XUEBAO 第8 期 式中：zk ，z (s k,n)分别为边缘流线和任意位置上的自由水面方程；n为自由表面流线分布参数，对于给 定的流线，n为常数；sk为在平面图上量得的边缘流线长度与L0之比；A、B、C、D、F为变系数，它们的形 式取决于参数n的引入方式详见下文；Π为无量纲流速水头， h ?2(z ? z + h + P ? P ) ∏ = ∏ ， φ H? f ? H? ? v 2 0 ? , z? 和 其中下标φ为分层数，ΠHφ为某一层起点处的无量纲流速水头， H = ∏ 2gL 0 z 分别表示任一层 H? 水面上的某高点与L P 和 0之比和该层起点处的水面高程与L0之比， ? P 分别表示任一层水面上的某点压强 H? 与ρgL0之比和该层起点处的压强与ρgL0之比，hf 为摩擦损失与L 0之比；L0为特征长度，一般可取边缘流 线水平投影的长度作为特征长度。 连续性方程 在弯道的起始断面处，沿水平方向将水深划分为若干流层，沿铅垂方向将水流划分 为若干流片，并且假定流片在流动过程中，侧面都保持柱面形状，这样沿每一流片的每一流层在流动中都 满足水流连续方程。于是，可推导出无量纲弯曲水流的连续方程为： F ( z / s )2 ∏ + ? ? δ = δ (2) h h H k 0 ( / ) 2 ∏ + ? ? F F z s H k H 式中：δh 为流层的铅垂厚度与L0之比，δh0 为起点断面处流层的铅垂厚度与L 0之比；H为起点断面处的值。 压强分布 沿水深的无量纲压强分布方程为： P = P h ? 0 ? 1 ? ? 2 ? ? ? ?∏ ? ? ? ? ? ? z z F + ? ? + ? ? ? h s s s ? ? ? ? ? ? ? ∫ ? ? k ? k k ? F( z s s 2 / D z / ) 2 0 ∏ ? ? + ? ? ? k k + ? (F ( z s G(F + ? ? ? / )2 ) 2 k T z / s ? ? k + ? ? ( z / s k ) )