圆柱形墩局部冲刷的三维数值模拟研究.pdfVIP

第十四届全国结构风工程学术会议论文集 圆柱形墩局部冲刷的三维数值模拟 祝志文刘震卿陈政清 (湖南大学风工程研究中心长沙410082) 摘要桥墩局部冲刷由于现场观测的困难性和试验研究的复杂性，难以获得其冲刷坑深度的合理估计．本文 以床面瞬时切应力作为泥沙起动及输运的水动力学机制建立了一个新的冲刷数值模型，在通用计算软件Fluent中使 用标准盘-￡湍流模型对圆柱形墩周围复杂流场做了数值模拟，并采用本文提出的数学模型计算河床底面坐标的实时变 化，使用C++实时修改网格文件，计算得到了圆柱墩局部冲刷坑的演化过程．通过与相关实验结果的比较，证明了 本文方法的可行性． 关键词桥墩，局部冲刷、数值模拟、最大冲刷深度． 1引言 桥墩局部冲刷是河道水流因桥墩对其强烈的扰动作用，所 形成的强烈漩涡淘刷桥墩周围的泥沙，形成的基础削弱现象。 局部从刷发生在水面以下，深度不大时，可通过涉水检查，但 当水深较大时，需潜水探测，因而难以直接观察。特别是洪水 季节，上游来流速度和水量大幅增加，桥墩局部冲刷坑发展迅 速，如设计时对最大冲刷坑深度估计，将导致桥墩偏斜或移位， 从而导致桥梁倒塌，交通中断，国内外发生过大量因局部冲刷 而导致桥梁破坏的例子【¨。因此，准确研究局部冲刷显得尤为 重要。 虽然近年来国内外一些学者提出了不少数值计算模型，但 是依然有一些不足之处需要改进：1)由于桥墩周围流场具有很 图i床面桥墩附近的流动情况 强的三维特性(见图1)，流场中有复杂的漩涡系统，采用二 维模型难以对其正确模拟；2)随着冲刷坑的发展，河床底面倾角逐渐增加，倾角对临界切应力有着 显著的影响，因此在计算过程中需要考虑河床倾角对临界切应力的影响；3)采用一般动网格技术会 造成网格质量不高，并很难保证在每个计算时间步下的网格满足使用湍流模型时的壁面函数的要求； 4)基于非通用计算程序对局部冲刷坑的数值模拟往往限制了其计算方法的推广应用。 2数值计算模型 2．1控制方程 连续性方程： 挚：o (1) dX, N．S方程： 导+专c眦吩，=专p考]-鼍+s (2)其中： =@l她殉)沏)，在本文中勋为垂直于床底方向的坐标；卜一时间O)；P+一包括静压在内的总压印却笛 动力粘度(Pa·s)。 917 第十四届全国结构风工程学术会议论文集 湍流动能k方程： c3， 掣+ib(pku,)=未[p+箦)芸]+G+G一班一肠 湍流耗散g方程： ck丢c…妊№乇妒譬(4) 挈+掣=亳陋E¨)Ea“j 。 2 ￡ 矿l a聩Jl J 的脉动值∞居)；K广可压湍流中脉动扩张的贡献，在本文中将液态水视为不可压流体，故％=0； G^『-平均速度梯度引起的萨生项，Gk：肛f竺+竺I罢；瓯一由于浮力引起的七的产生项，本文中 xt ＼Oxj j axl 取G=0；标准缸￡湍流模型中的cl。、Q。、G。、0以及与湍动能乏和湍动耗散率对应的Pmdtl数—吼， ％为经验常数，取值如下：Cl。=1．44、c2。=1．92、c3。=O、巳=O．09、砚=1．0，％=1．3． 3．河床变形模型 3．2推移质输沙