桥墩局部冲刷发展过程三维动网格模拟.docVIP

桥墩局部冲刷发展过程三维动网格模拟 　　摘要：基于FLUENT软件的动网格更新技术和用户自定义函数功能实现了桥墩局部冲刷过程的三维动态模拟。以Melville经典冲刷试验为原型，建立数值模型。将河床面设置为主要的动边界，当床面结点瞬时剪应力大于临界剪应力时，结点位置下移，表现为冲刷，引入Van Rijn提出的沉积输运函数来控制河床面各结点的运动速度。数值模拟结果在流场形态，冲坑发生发展过程及冲坑形态均与试验结果较为吻合，模拟的冲坑深度略小于试验结果，误差约13%。误差产生的主要原因为基于雷诺平均N-S的湍流模型不能有效地反应钝形桥墩前端湍流脉动的影响 关键词：局部冲刷；动网格更新；数值模拟；桥墩；Melville冲刷试验 中图分类号：U442 文献标识码：A 文章编号：1672-1683（2017）02-0132-06 冲积河道在遇到阻水构筑物（桥墩、丁坝等）时，构筑物周围河床的局部冲刷对阻水构筑物的稳定有很大的影响。对于桥墩而言，水流在遇到桥墩后，由于桥墩的阻水使得过水面积减小，墩周流速增大，河床剪应力增加，墩周河床沉积物被水流搬运，墩周床面高程逐步降低，并产生冲坑，导致桥墩基础的埋深减小，进而会导致桥梁的倒塌，甚至生命和财产的损失 桥墩的局部冲刷是一个动态的发展过程，影响因素众多，空间分布具有很强的三维特性，这就使得冲刷模型试验成为以往研究冲刷问题的主要手段。但模型实验存在费用高，无法普遍应用，条件单一，存在模型尺寸效应等不确定因素，数值模拟方法的不断改进使得其作为一种研究手段越来越显示出其不可替代的作用 近年来，国内外学者针对桥墩冲刷三维性态发展开展了一系列数值模拟研究。Ehteram运用SSIIM软件对桥台的冲刷过程进行了三维模拟，得到了冲刷坑深度和形状并与试验结果进行了比较。Khosronejad对不同横截面形状的桥墩进行了三维动床模拟，采用了流固耦合曲线浸入边界的技术。Kim采用大涡模拟的方法对相邻的两个圆柱形墩的局部冲刷坑进行了模拟，得到的最大冲深位置与试验结果较为一致。韦雁机等基于Open-FOAM开源软件的动网格技术，用输沙率计算床面地形随时间的变化，构建起桩周局部冲刷的动态三维数学模型。祝志文等根据床底泥沙的单宽体积输沙率得到河床高程坐标的瞬时变化，采用边界自适应网格技术修改动边界计算域网格，得到圆柱形桥墩周围局部冲刷坑的演化过程。以上研究对局部冲刷的数值模拟起到了很好的推动作用，在实际应用中多少都存在一些不足的地方，如采用虚拟的浸入边界、地形函数等很难与实际条件一致，大涡模拟或分离涡的模拟计算消耗极大，基于单宽体积输沙率来计算河床的变形，计算过程复杂，涉及到梯度的计算，会使误差增加等 本文基于CFD计算软件FLUENT的动网格技术和用户白定义函数（UDF）功能，使用基于雷诺平均N-S模型的Realizable k-e湍流模型，将床面瞬时剪应力和临界剪应力带入Van Rijn沉积输运函数，得到床面坐标的变化，通过网格重构和弹性光顺结合的方法来不断修正变形较大的网格，实现了局部冲刷过程的动态模拟 1数值模型 1.1物理模型的选取 在已有的?_刷试验中，Melville和Raudkivi（MR）对局部冲刷坑发展的三个不同的阶段进行了相对较详细的定量描述。所以本次研究选取MR的经典冲刷试验资料建立数值模型，并进行对比分析。MR试验水槽长19 m，宽0.456 m，水深0.15m，平均来流流速为0.25 m/s。模型布置见图1。MR分别选取初始定床阶段（测定了河床面附近的流速），中间发展阶段（冲刷30 min时，冲坑深度达到0.04 m）和冲刷终止的平衡阶段进行分析，给出了详细的试验结果。试验中前30 min发展较为剧烈，而之后冲刷发展开始缓慢。30 min时的冲刷深度达到总冲刷深度的75%。数值模拟取前30 min进行研究 1.2计算域及网格划分 试验研究表明，圆柱形桥墩绕流流态以x轴基本呈对称分布，因此取其中一半作为本次模拟的计算域，以缩短计算时间。根据试验布置和模拟要求，将计算域高度设置为15 cm，宽度设置为3d即15.24cm，圆柱上游距桥墩中心为3d，出口处要求尾流充分发展，所以设置下游距桥墩中心为10d[s]即50.8 cm。取圆柱竖向为z轴方向，床面为x-y面水流方向为x轴正方向。具体见图2 本次模拟为动床模拟，桥墩周围由于冲刷作用使得局部变形较大，且变形不规则，所以选用适应变形能力较强的四面体非结构化网格。为了提高床面剪应力的模拟精度，在床面设置了0.2 cm（约0.5d50）边界层，且边界层会随着床面结点的移动而跟随移动，这样更大程度上保证了床面剪应力获取的精度。因此本次模拟在床面边界层内为三棱柱体网格，其余部分为四面体网格。由于墩周及靠