扭切式鼻坎挑射水流作用机理分析.docxVIP

? ? 扭切式鼻坎挑射水流作用机理分析 ? ? 吴锦钢，邱 勇，周鑫宇，葛亚飞，黄梓涵 (云南农业大学水利学院，云南 昆明 650201) 0 引 言 泄水建筑物末端挑坎体型的选择是水利水电工程中消能防冲设计中的关键问题[1]。随着水利工程开发条件日趋复杂，传统的连续式鼻坎难以满足工程布置要求，需要在此基础上对体型作出相应的修改或者优化，得到能够改变出坎水舌空间形态或水股落点位置的异型鼻坎[2]。为此，许多学者在鼻坎体型选择以及水力特性研究方面做了大量的工作。练继建等[3]以提高消能率和减轻雾化现象为目的，探究了水舌在平面上对称横向扩散的舌型挑坎与连续坎水舌的液滴粒径与下游降雨强度分布规律；孙颖等[4]为防止下泄水流冲刷河岸，采用双曲型挑坎，使得出坎水舌在对冲后仍沿对称轴线方向扩散；Deng等[5]在针对大流量、窄河谷需要将水舌纵向拉伸情况下，优化出一种新型窄缝挑坎，结合数值模拟和模型试验，研究了水舌形态的形成机理；刘昉等[6]采用正交试验，对研究了扭曲式鼻坎中的参数变化对出坎水舌和下游河床的敏感性分析。 针对狭窄弯曲型河谷布置泄水建筑物，文献[7- 8]依托水工模型试验成果，采用扭切式挑流鼻坎，通过边墙的偏折，迫使水舌在平面上转向，给出了出坎水舌形态及水股落点，解决了水流平顺归河问题。本文采用连续性方程、动量方程、弯道超高经验公式结合冲击波理论，对扭切式鼻坎挑射水流作用机理进行分析，研究鼻坎内轴线立面挑射水流的水深和流速沿程变化规律、边墙偏转侧平面偏折水流水面超高经验公式以及底部斜射水流在出坎处的出射偏折角度。 1 反弧鼻坎水流结构分区 根据水工模型试验[7- 8]和数值模拟成果[9]，扭切式鼻坎水流出坎后能够在空间上形成大面积的拉开，同时促使水股平面偏折，运动较为复杂。本文结合反弧段挑坎水流运动状态(见图1)的不同，将其概化为立面挑射水流、平面偏折水流和底部斜射水流[9]，选取微元体进行受力分析。 图1 反弧段出坎水流平面形态 2 立面挑射水流 轴线附近直线前行的水流进入挑坎段后，水面沿反弧底板不断爬升、水深增加、流速降低。 2.1 连续方程 假定水流在轴线横向上呈均匀分布[10]，流速分布沿水深满足幂指数分布规律[11]，则水流的连续方程为 (1) 式中，Q为流量；ve为沿微元体法向距离鼻坎底部高度z的流速，m/s；v0为水流表面的速度，m/s；H为鼻坎内水流深度，m；m为幂指数的次数。 忽略指数m的沿程分布变化，连续方程可化简为 (2) 2.2 动量方程 选取反弧段鼻坎水流作为控制体，如图2所示，通过动量方程推求沿轴线方向的反弧段水流沿程水深和流速。 图2 反弧段动量方程控制体 将动量变化量和外力作用项分别代入等式两边，得到 (3) (4) 等式右侧，由于压强P沿程变化，其表达式仍需推导，同时假定水流流线沿水深分布均匀且不考虑水流在行进过程中边壁对水体产生的摩阻力[14]。 2.2.1 反弧段动水压强沿程分布 在上述反弧段上任取一段水深为H、长度为ds、高度为dh的单宽微元体，微元体轴线距离反弧段曲面底部高度为h，如图3所示。 图3 动水压强微元体示意 根据牛顿第二定律，微元体沿轴线的运动方程，可表示为 (5) 式中，ρ为水的密度，kg/m3；g为重力加速度，m/s2；R为反弧段结构半径，m；v为水流速度，m/s。 对于常数C，由边界条件(h=R-H)时，P=0，解出C=-ρv2lnH-pg(R-H)cosα，将C代入式(5)，整理后得到反弧段动水压强沿程分布表达式，即 (6) 2.2.2 反弧段沿程水深和流速 基于上述分析，将式(4)、式(6)带入式(3)中，整理化简后得 (7) (8) 式(2)和式(8)构成反弧段封闭方程组，即 在此基础上，利用龙格-库塔法[15]，可以求得反弧段的沿程水深H及流速变化v0。 基于该封闭方程组，结合挑流消能公式[16]，可进一步计算出轴线附近的立面挑射水流挑距。 3 平面偏折水流 扭切式挑流鼻坎的一侧边墙以弧形向轴线方向偏转，上游端与边墙相切，下游端与反弧段末端轴线相接，其特殊体型使得水流结构发生变化，出坎水舌在平面上发生横向偏转[7]，从而影响水舌落点的分布。 图4 弯道环流条形水体 (9) 式中，v为水流速度，m/s；r为左侧圆弧边墙的半径，m；B为鼻坎宽度，m。 弯道超高计算公式(9)在绝大多数情况下，其底坡为平坡或坡度很小，而扭切式鼻坎的底板曲率较大，水流爬升效应明显[18]，对于超高的影响不能忽略。 考虑底板曲率的影响，计算该计算点(见图5)与前一点的高差Δz，Δz=R(cosα1-cosα2)。 图5 反弧段等分 由于底板曲率的存在，使得爬升的水流超高值远低于底板为平坡时的数值[19]。若想要计算某一点的超高值，可将弯道超高理论计算值ΔH乘以这一点与等分后前一点