辽河砂堤的渗流特性模拟分析.docxVIP

辽河砂堤的渗流特性模拟分析 摘要：渗流是堤防安全的主要问题之一，通过研究和分析渗流问题，可以有效地减小堤防被渗流破坏的可能性，也可以根据理论依据为实际堤防渗流问题提出解决措施。本文主要研究辽河干流典型河段粉细砂筑堤在汛期高水位运行时，堤防渗透变形的基本规律。以达西定律为基础，采用有限体积法进行建模、网格划分、计算，对不同边坡比率的渗流规律进行分析研究。 渗流问题是破坏水利工程安全的最大原因之一，破坏的后果对人类影响极大，故而引起建造者的广泛关注，开始研究和分析渗流问题，以期可减小渗流破坏的可能性。Fluent软件目前是CFD软件中较为领先的一款计算流体动力学软件，Fluent支持多种网格，用户使用起来更加简便。本文以达西定律为基础，通过Fluent软件对1∶2和1∶3比率的边坡进行渗流规律分析。1 研究方法1.1 渗流运动基本方程在实际工程中渗流问题是非常复杂的，如渗流区域有规则和不规则之分，当渗流区域不规则时，不能确认它是一元流或渐变流(1）连续方程。根据质量守恒原理，水流流动时，渗流场中的水在某个单元体积内进出流量速率之差等于这个单元体内水的增减速率。水体质量m随时间变化率如式（1）所示：式中：n为土体孔隙率；ρ为水的密度，g/cm(2）运动方程。根据流体力学中不可压缩黏性液体的运动微分方程可以模拟出地下水运动方程。该方程可以模拟黏性不可压缩流体动量守恒的运动过程。流体使用假定方程是连续的，并且涉及到的渗流场必须是可微的。该方程是由表示剪切应力和黏度引起的原始应力的流量梯度导出的。纳维—斯托克斯方程的矢量形式如式（2）所示：式中：V为速度矢量，m/s;g为重力加速度，m/s上述公式的物理意义是表示压力与黏性力、加速度力、质量力之间的关系。当考虑黏性时，它是实际黏性流体的微分形式的动量方程。当它不考虑黏性时，是理想流体运动的欧拉方程。(3）边界条件。渗流的边界条件有三类。给定水头的边界：当在任意一处的边界（设为W1）上的每个点的每时每刻水头都是已知的；给定流量的边界：当知道任意一处边界（设为W2）的单位面积上流出或流入的流量时（流入时用正值，流出时用负值）；混合边界条件：指流入或流出的流量与水层边界的内外水头差之间有一定的线性关系1.2 有限体积法及Fluent(1）有限体积法。一种广泛应用于CFD领域的离散化方法有限体积法也称为控制体积法，其基本思想(2)Fluent软件。Fluent软件是世界上目前最流行的计算流体动力学软件，是流体流动和热交换领域中用于模拟和分析复杂几何区域中的专用软件，它可以解决涉及流体、传热和化学反应的任何技术问题。Fluent支持多种网格，之所以选择使用非结构化的网格来划分杂乱的不规则区域，是因为非结构化网格尺寸更容易处理杂乱的形状进而计算区域。用户还可以简便、灵活的划分网格，这要得益于Fluent软件为客户提供了满足特定要求的各种人性化的功能。2 辽河堤防基于Fluent的渗流分析2.1 工程概况辽河是我国七大江河之一，它起源于河北省七老图山脉的光头山，流经河北省、内蒙古自治区、吉林省、辽宁省四省，至盘山入渤海。根据《辽河流域防洪规划》和《辽河流域综合规划（2012—2030年）》防洪工程总体布局，辽河的防洪体系已基本建成，防洪骨干以水库、拦河分洪枢纽与下游河道堤防相结合。辽河干流福德店以下河段，以石佛寺水库及支流水库构建骨干防洪工程，同时规划加固河道堤防、整治河道等，完善铁岭、盘山市城市防洪工程体系。现状干流两岸堤防共计646.286 km，其中左岸堤防长度328.553km，右岸堤防长度317.733km，除去石佛寺水库副坝，总长度608.863km。2.2 水头分析两种坡比的水头分析，坡比分别为1∶2和1∶3。首先，在Workbench上创建Fluent，进行建模分析，绘制砂堤横断面，砂堤顶宽为5m，底宽为40m，高为5m，断面如图1。网格划分，先给各边界条件命名，在进行网格划分，节点数为71 672；网格数为66 429，并进行物理性质、边界条件等设置。进行计算，输出结果如图2。将坡比改为1∶3创建模型，按照同样的步骤得到结果如图3。根据两种不同的坡比模拟出的水头分布可以看出当水头越大时，水由高势能处渗向低势能处，所以渗透的范围越大。水头差越高，渗流速度越快。3 结论本文运用有限体积法对砂堤进行了准确的数值特性模拟分析，为辽河砂堤治理提供了基本理论依据，并得出以下结论：（1）对砂堤的不同坡比的水头进行非稳定渗流视为稳定渗流来研究，在水头差越低的地方，渗流速率越小，水头差越高的位置，渗流速率越大，如果渗流物两侧没有水头差时，将不发生渗流。（2）对于砂堤不同配合比的砂料条件经行渗流场模拟分析，不