鱼骨坝工程水流结构和水毁机理研究.pdfVIP

鱼骨坝工程水流结构与水毁机理研究 刘怀汉①潘美元②付中敏⑦曹民雄③ 摘要鱼骨坝工程是鱼嘴工程中较为典型的布置形式。本文通过结合水槽概化试验和二维水流数值计算成果， 研究了鱼骨坝修建于平整滩面或心滩上、来流与脊坝夹角为0。与15。、淹没与非淹没等工况下鱼骨坝周围的水位、流速、 流场和动水压强等水力特性，并对其水毁原因进行了分析。研究成果对鱼骨坝工程设计具有一定的理论价值。 关键词鱼骨坝工程；水力特性；水毁；水槽试验；数值计算 1 引言 鱼骨坝是一种布置形式较为复杂的鱼嘴建筑物形式，一般由顺水流方向的脊坝和垂直于脊坝轴线 的多条刺坝组成。脊坝主要用于分流、分沙和归顺水流方向，刺坝调节环流的运动，并增强坝体的稳定。 目前，鱼骨坝工程已在国内湘江下摄司滩航道整治中取得了较为成功的整治效果。在长江航道整治中， 鱼骨坝工程也有了初步的尝试，如在建的东流水道鱼骨坝工程；此外，长江陆溪口新洲整治工程中也对 鱼骨坝方案进行了讨论。然而，由于鱼骨坝周围水流条件复杂，对其周围水流特性和水毁的理论和试验 研究较少。为此，本文结合长江东流水道整治工程，利用数学模型，并根据水槽试验数据进行验证，研究 鱼骨坝的水位、流速、流场和动水压强等水力特性，并对水毁原因进行分析。 2数学模型的建立和验证 2．1平面二维水流基本方程 采用平面二维浅水方程作为水流计算的控制方程，主要方程如下： 连续方程： 照+地丝+地盟：0 Ot Ox Oy x方向运动方程： Ou批讹Ou孤¨Ou砂+g(警+鲁)一fv一詈’≯02u+cg：u^∥万=o y方向运动方程： Ov秕+u缸Ov+移ayOv+g(丛Oy+鲁)+丘一詈。等+cg：v危“雨=o 其中：x、l，是水平坐标轴，H、秽为x、y轴向流速，t是时间变量，g为重力加速度，h为水深，eto为河床泥 面高程，P为水流密度／为科氏力参数(f=2wsin6，∞为地球旋转角速度，咖为纬度)，因计算区域较小， 可不考虑科氏力的影响，计算时利用曼宁公式进行转换。 ②长江航道规划设计研究院。 ③南京水利科学研究院。 178 2．2计算方法和求解过程 利用流体力学的有限单元法米计算控制方程的数值解，其基本求解过程为1”：(1)控制方程分别在 时间和空间上进行离散。其中时间离散采用差分法，空间离散采用有限单元法；(2)通过伽辽金加权余 晕法把拧制方程从偏微分方程转变成代数方程组；(3)数值求解代数方程组。 2．3数模的建立与边界条件 x4m×0 工程方案采用概化后的鱼骨坝，放置在35m 4m(长x宽×高)的矩形水槽中进行试验研 究，槽底坡度为05％，模型}{算段长度为23．02m。根据地形特征对计算区域进行了网格剖分(图1)， 计算网格由三角形六结点和四边形八结点相结合的等参变网格单元构成，共有网格结点27063个，单元 lm，鱼骨坝坝体附近区域采用网 9048个。单元网格尺寸沿水流方向一般为04m．垂直水流方向为0 03m×0 格加密技术，网格尺寸一般为0 03m(长x宽)。 糙率按照不同的区域与材料共分四个区：水槽边壁区、坝体、水泥床底与工程区域，糙率分别取 00135、00120、00115和0010。 模型的计算边界：上游为流量边界、下游为水位边界。为了研究鱼骨坝的水流特性，共进行了淹没 和非淹没情况下6组定床水槽试验进行计算(表1)。 模型计算边界条件 表 非淹＆情况 淹没情况 I 计算工况 I 2 3 ．