肇庆市水务技术中心526040
摘要：桥梁是公路交通建设的重要组成部分。桥梁横跨河流、侵占河道行洪面积、壅高水位、改变局部流态，对沿河两岸堤防造成影响。按照河道管理范围内建设项目管理的相关规定，桥梁在建设时，需进行防洪影响评价，而对桥梁壅高水位及引起的冲淤变化计算是防洪影响评价计算的一项重要内容。为此结合肇庆市广宁县金场绥江大桥工程，利用有限元法，对桥梁建成后河流的局部水位壅高和流态变化进行分析评价。
关键词：防洪评价；平面二维计算；桥梁；数学模型。
1二维数学模型及计算内容
天然河道一般边界曲折、地形复杂，对于复杂河段的水流运动数值模拟，多采用基于曲线网格的坐标变换方法，其中正交曲线变换和一般（非正交）曲线变换方法是两种最常用的方法。本研究采用一般曲线坐标变换方法。一般曲线变换不受计算网格必须严格保证正交的限制，网格生成也较灵活。
1.1模型计算的基本原理
经一般曲线变换后的平面二维水流模型控制方程为：
求解过程中，为避免水位波动，控制体交界面上的流速采用动量插值处理；为避免计算迭代过程中出现溢出，采用了Patankar和Spalding提出的欠松弛技术，即在离散方程式中引入欠松弛因子，以改善离散方程式中系数的对角占优程度。
1.3计算边界条件
平面二维水流模型中，边界条件通常包括河道进出口边界、岸边界及动边界处理等。本模型中：
1）进口边界：根据已知进口全断面流量，给定入流单宽流量沿断面的横向分布。
2）出口边界：给定出口断面的水位。
3）岸边界：岸边界为非滑移边界，给定其流速为零。
4）动边界：本模型采用“冻结”法进行动边界处理，即根据水位结点处河底高程来判断该网格单元是否露出水面，若不露出，糙率取正常值，反之，糙率取一个接近于无穷大的正数。同时为了不影响水流控制方程的求解，在露出水面的结点处需给定一个薄水层，一般给定其厚度为0.5cm。
1.4计算过程
基于上述基本方程、求解方法和计算网格，即可进行程序的编写和调试，建立二维数学模型。在给定边界和初始条件下，进行数值模型实验，模型参数通过率定和验证后，即可进行相关计算。
1.5桥墩概化
目前对桩墩的模拟主要有两种模式：局部阻力修正法和直接模拟法。本报告采用的是局部阻力修正法，计算时将桥墩、承台及防撞柱均按实体不透水处理，这样考虑的计算结果是偏安全的。
2实例计算
2.1研究范围及网格布置
考虑到实测地形资料和本河段工程的影响范围，模型上边界取自工程上游430m处，下边界取至工程下游430m处，总模拟长度860m。
本模型采用非结构三角形网格，在研究范围内共布设了4178个三角形计算单元，网格尺寸在1~30m之间，桥墩局部区域网格加密处理（最小网格边长1m）。
2.2工程计算水文条件
金场绥江大桥位于金场水出口上游500m的位置，距离下游东乡水电站约4.69km，根据东乡水电站的出流特性和本河段的水面比降关系，模型上游进口流量直接取用东乡电站的出流量减去金场水的流量，模型下游出口水位取东乡水电站水位上延至本工程下游430m边界处的水位值。
防洪评价计算主要包括对河道水位、流速、流态等几方面，根据评价项目，选择以下水文条件进行计算：金场绥江大桥20年一遇工况
2.3计算结果分析
2.3.1壅水分析
工程后流速差等值线见图1。在二十年一遇水文条件下，受桥墩阻水影响，金场绥江大桥对上游洪水水位最大抬高值分别为0.06m，水位壅高值大于0.01m的最大影响范围为拟建大桥上游120m左右，桥址下游水位有所减小，最大减小距离桥址下游90m左右。
综上所述，工程对河道防洪（潮）水位的影响很小。
图120年一遇下工程前后水位变化图（红色：增加，绿色：减小）
2.3.1流速分析