桥梁河道三维水力与泥沙数值模型计算研究 A Study on Three Dimensional Hydraulic and Sediment Numerical Model Calculation for Bridge River.pdfVIP

HIGHWAY No．4 公路2012年4月第4期 Apr．2012 中图分类号：U442．39 文献标识码：A 文章编号：0451—0712(2012)04--0021--05 桥梁河道三维水力与泥沙数值模型计算研究 贾界峰，鲍卫刚 (中交路桥技术有限公司北京市100029) 摘要：从水动力和泥沙运动方程出发，采用雷诺时均的N—S方程，以志一e双方程模型对紊流方程进行封闭， 采用垂向分层的方法建立桥梁墩台冲刷的三维水动力模型。泥沙计算考虑悬移质和推移质泥沙输运，通过求解动网 格方程更新整个计算域的网格坐标。通过数值模拟结果与模型试验结果的比较，建立的桥梁墩台冲刷数值模型具有 较高的精度。 关键词：局部冲刷；三维；数值模型；动网格；冲刷深度 近年来，由于桥梁失稳造成的灾害越来越多，而 紊流是流体存在的一种基本形态。虽然N—S 桥墩局部冲刷是引起桥梁失稳的主要原因。因此， 方程可以描述紊流效应的一切细节，但是由于涡旋 正确预测桥墩局部冲刷深度是桥梁设计和养护的重 的复杂性和尺度量级，目前很难通过直接求解N— 要依据之一。自18世纪70年代开始，相关学者就S方程来解决紊流问题。因此引入紊流模型，对紊 开始对桥墩局部冲刷进行研究，但是由于现象复杂， 流运动采用模型化进行处理。 对其机理认识和理论研究都还有限。在科研和生产 目前常用的紊流模型主要有重整化群(RNG) 中，通常使用原型观测、模型试验和数值模拟三种方 模型、大涡模拟(LES)和k一甜模型等。重整化群模 法对桥墩局部冲刷进行研究计算。由于原型观测和 型方法认为，紊流在惯性子区域中可以用随机力作 模型试验受到费用高、耗时久、精度有限等因素的制 用下的N—S方程进行描述。使受边界条件和初始 约，以计算机应用为基础的数值模拟方法越来越引 条件约束的N—S方程经变换后与不受任何约束的 起研究者的重视。 其他Fourier变换等效。通过这种变化可以系统地 本文使用数值模拟方法建立了桥梁河道三维水 剔除高频小尺度紊流项，其他大尺度项可以通过数 力与泥沙模型，通过三维数值模型研究桥墩的局部 值方法求解[1]。大涡模拟是对小尺度涡旋运动通过 冲刷，了解冲刷坑形态及最大冲刷深度的发展规律。 滤波的方法进行处理，只研究尺寸较大涡旋的运动， 并以实际工程为依据，对水沙模型的计算结果进行 涡旋的尺度由网格尺寸进行控制，因此其精度在很 验证。 大程度上依赖于网格的尺度。k一∞模型是利用求 解k和甜的输运方程来封闭雷诺平均的N—S方程 1桥墩局部冲刷模型的建立与验证 实现对紊流的模拟。 采用的数值模型由流场模型、泥沙模型和底床 (2)紊流模型。 模型三部分组成。流场模型为时均雷诺方程，并以 由于紊流的高度复杂性和当前量测技术等方面 k—s紊流模型封闭计算流场。泥沙模型包括泥沙 的限制，无法对实际紊流的N—S方程提供适当的 的临界起动条件、下落速度和输运公式。由于冲刷 初始条件和边界条件，同时受计算水平的影响，尚无 计算过程中经常在床面出现意外的地形突起和沙 法直接求解N—S方程。Reynolds根据统计概念， 坡，使计算网格边缘产生畸变，因此在模拟计算过程 将速度、压力等运动要素表示为时均量与脉动量之 中采用底面光滑修正。