三维河道水力分析计算系统在桥梁优化设计的应用研究

HEC-RAS三维河道水力分析计算系统在桥梁优化设计的应用研究摘要：基于计算流体动力学的方法，运用三维河道水力学模型HEC-RAS（RiverAnalyseSystem）分析阿尔及利亚贝贾亚高速公路在建项目2号互通桥梁在百年一遇的最大洪峰流量下，不同桥梁孔跨方案在相同地质条件下对桥梁最大冲刷深度的影响，对不同桥梁孔跨方案下三维河道水力学模型计算结果进行分析研究，确定经济合理的桥梁设计方案，为今后国内外同类桥梁设计的方案研究提供了新的研究方法。关键词：HEC-RAS；冲刷计算；桥梁孔跨方案；研究HEC-RASthree-dimensionalriverhydrauliccalculationandanalysissystemintheresearchandapplicationofbridgedesignYangLeilei(ChinaRailwayFirstSurveyandDesignInstituteGroupLtd.,Xi’an710043,China)Abstract:Methodsbasedoncomputationalfluiddynamics,using3DhydraulicsmodelHEC-RAS(RiverAnalyseSystem)analysisofAlandLeighBbeJayainHighwayProjectNo.2interchangebridgeinthemaximumpeakdischargeofahundredyears,differentbridgespanofinfluenceonbridgemaximumscourdepthinthesamegeologicalconditions,thedifferentbridgespan3Dhydraulicsmodelforcalculationresultswereanalyzed,determinethebridgeeconomicrationaldesignscheme,providesnewmethodsfortheresearchofdomesticandforeignsimilarbridgedesigninthefuture.Keywords:HEC-RAS,Calculationofscour,Bridgespanschemes,Research桥梁孔跨布置是影响墩台冲刷主要因素，墩台冲刷一直以来被认为是导致桥梁破坏的重要原因之一，合理的桥梁孔跨布置是桥梁安全运营的重要保障。桥梁一般冲刷坑形态及墩台局部冲刷坑深度的正确预测可以为桥梁的设计、监控及养护提供重要的依据。虽然从上世纪50年代开始，开展了众多针对桥梁墩台局部冲刷的研究，并取得了较大进展，但是由于桥梁冲刷问题的极其复杂性，因冲刷导致的桥梁垮塌事件依然经常发生，给社会和财产带来了巨大的损失。1桥梁冲刷危害及其影响因素桥梁一般冲刷坑形态及墩台最大局部冲刷坑深度的正确预测对桥梁的设计及维护有重要的影响。如果最大局部冲刷坑被过高估计，可能造成桥梁建设成本的大幅提高；而如果局部冲刷坑被过小估计，墩台冲刷将导致基础掏空，桥梁下部结构的承载力将会大幅下降，可能造成桥梁倾斜甚至桥梁垮塌。影响桥梁一般冲刷坑形态及墩台最大局部冲刷坑深度的主要因素为河道断面宽度，桥梁孔数跨径[1,2]，桥梁墩台与水流速度的夹角，桥位处河流的流量及流速，河床地质成份的组成粒径等。河道断面宽度及河流流量、流速为自然条件下所形成，桥梁孔数、跨径及桥梁墩台与水流速度的夹角为人为所修建，因此合理的设计桥梁长度及墩、台与水流夹角成为影响一般冲刷坑形态及墩台最大局部冲刷坑深度的最主要因素。2国内桥梁冲刷现状2.1一般冲刷我国计算桥梁一般冲刷的公式，目前普遍采用1964年我国“桥渡冲刷计算学术会议”上推荐的64-1公式[3]。，64-1公式是利用我国各地桥梁实测资料并根据冲止流速的概念建立起来的。桥下一般冲刷停止时的垂线平均流速称为冲止流速[3]。桥下断面内任意垂线在一般冲刷过程中垂线平均流速降低到该垂线的冲止流速流速时，冲刷即停止，此时达到最大一般冲刷的垂线水深，河槽一般冲刷的基本形式为：536135dEhhLQAhmPpm式中，hpm为桥下断面一般冲刷后的最大水深；Qp为设计流量；L为桥孔净长；hm为设计断面上桥孔部分的最大垂线水深；h为设计断面上桥孔部分的平均水深；E为汛期含沙量有关的参数；d为河床土壤平均粒径；A为单宽流量集中系数。2.2局部冲刷我国计算桥梁局部冲刷的公式，目前普遍采用1964年我国“桥渡冲刷计算学术会议”上推荐的65-1公式[3]，65-1公式是以我国自己进行的大量室内试验资料和汛期实桥观测资料为依据建立的。当0vv时，'06.011vvBKKhB；当0vv时，'00'0'6.011vvvvvvBKKhB式中，hB为桥下局部冲刷坑深度；K为墩型系数；1B为桥墩计算宽度；v为一般冲刷后的垂线平均流速；0v为河床泥沙起动流速；0'v为墩前始冲流速；1nk为河床粒径的影响系数；E为汛期含沙量有关的参数；1n为指数。3国外桥梁冲刷计算模型HEC-RAS（RiverAnalyseSystem）是由美国陆军工程兵团水文工程中心开发的水面线计算软件包，适用于河道稳定和非稳定流三维水力计算，其功能强大，可进行各种涉水建筑物(如桥梁、涵洞、防洪堤、堰、水库、块状阻水建筑物等)的水面线分析计算，同时可生成横断面形态图、流量及水位过程曲线、复式河道三维断面图等各种分析图表，使用起来十分方便简捷，目前已成为使用最广泛的三维水力计算软件之一。3.1三维河道水力学计算简介河道水力分析模型——HEC-RAS是一个针对三维恒定/非恒定流的水力模型[2]，主要用于明渠河道流动分析和洪泛平原区域的确定。模型所得结果可以用于洪水区域管理以及洪水安全研究分析，用以评价洪水淹没区域的范围及危害程度。如在进行河道整治以及新建桥梁等工程的时候，就要分析考虑河道壅水高度、流速变化、桥涵冲刷[4,5,8]等这些因素对河流输水，城市防洪的影响。HEC-RAS主要有以下4部分组成：恒定流水面线计算、非恒定流模拟、可运动边界泥沙输移计算、水质分析[7,12]。本文主要利用HEC-RAS对非恒定流模拟、河道冲刷计算功能进行影响桥梁墩、台冲刷深度因素分析，经分析比较找出桥梁孔跨布置的最优方案。HEC-RAS的非恒定流模拟式基于连续方程和动量方程，其中连续方程[2]为：0iixut式中，为流体密度；u为流速；下标遵守爱因斯坦求和约定。动量方程[4]为：ijijijifixxuvxpfxuutu21式中，f为质量力；p为压力；v为压力流体运动粘滞系数。HEC-RAS的一般冲刷计算是基于“清水劳尔森”方程[4]进行计算的，其方程为：0732232502225.140ywDQysG式中，Gys为一般冲刷深度；2Q为河流流量；50D为河床底部沉淀物的平均直径；2W河床底部的宽度；0y为对应流量2Q时的河床平均吃水深度。HEC-RAS中桥梁墩局部冲刷计算是基于“科罗拉多州立大学”研究方程[3]进行计算的，其方程为：43.0135.0143210.2FrAYKKKKaysL式中，Lys为桥墩的局部冲刷深度；1y为桥墩上游的吃水深度；a为桥墩宽度；1k为桥墩前端形状参数，作为循环障碍时为100，作为线性障碍时为0.75，作为矩形障碍时为1.03；2k为和桥墩中心线相比与水侵蚀角度有关的参数；3k为一个可以根据河床或可见沙丘来增加冲刷深度的系数；3k为河床铺设可能的系数；1Fr为桥墩上游的弗劳德数量。HEC-RAS中桥梁台局部冲刷计算：第一种情况（L/y125）是基于“USAGE”方程[3]进行计算的，其方程为：33.0121155.04FrKKyys第二种情况（L/y125）是基于“弗勒利希”方程[3]进行计算的，其方程为：161.057.0143.02127.2yFryLKKys式中，sy为桥台的局部冲刷深度；1y为桥台基础上游的吃水深度；1k为桥台形状的校对因素；2k为水流侵蚀角度的校对因素，13.00290k；Fr为弗劳德数量；L侵蚀长度：桥台阻碍了水流的流动长度。3.2三维河道水力学模型实例分析3.2.1项目背景阿尔及利亚贝贾亚港口至东西高速公路连接线北起港口城市贝贾亚，南与东西高速公路相连接，道路全长100公里，途经5座城市，设计为双向六车道，共分3个标段，S1标段长度为21km，S2标段长度为26km，S3标段长度为53km。S1标段穿越贝贾亚省的贝贾亚、TalaHamza、吉尔河、Amizour主要城镇，S1标段位于苏马姆河的入海口处，多次穿越苏马姆河河流，并通过立交桥与国道RN09和RN75连接。苏马姆河谷位于阿尔及利亚的东北中央，卡比利亚，阿尔及尔和君士坦丁的中部。其东北—西南方向狭长，介于祖赫祖拉山西部和其延伸到阿格巴尔—古拉亚支脉东北，比班南部和延伸到巴博尔东部中间。自布塞勒哈姆河和萨赫勒河的汇流处的上游部分是苏马姆的开端——位于阿克布西南2公里处，苏马姆口是下游的终点——位于贝贾亚的东郊。3.2.2基本情况阿尔及利亚贝贾亚高速公路2号互通位于苏马姆河弯曲河道上，其桥梁设计采用百年一遇洪水频率[6,10]，我国桥梁冲刷计算基于60年代的试验公式，该试验公式仅适用于直线河道上水力计算，无法解决弯曲河道上水流方向多变性的需要。HEC-RAS三维河道水力分析计算系统精确的适用于各种空间变化河流段，本文通过对拟建桥梁河道上，在相同地质条件下，不同孔跨布置对桥位处一般冲刷、局部冲刷的影响分析，从而找出经济合理的桥梁孔跨布置方案[9,11]。2号互通位于S1标段PK6+900处，其桥梁段横跨苏马姆河主河槽区，根据收集到苏马姆河谷水文站对河流的统计资料可知，苏马姆河百年一遇最大洪水流量为2532m3/s，河床经常年冲刷与沉积，表层为含沙类粉质粘土，表层以下为沙砾石，河床沉积物的直径取值范围为0.1~100mm，人工对桥墩换填材料防护粒径为300mm。3.2.3模拟分析利用河道水力分析模型（HEC-RAS）建立水力模型，其中应注意按河道实际情况对有桥状态无效水流区域的模拟。第一种情况，苏马姆河无桥自然状态：图1河道平面图图2河道横断面图图3河道空间立体图表1无桥状态计算结果第二种情况，苏马姆河布设3-36m桥状态：图4桥位平面布置图图5桥梁墩、台冲刷线图图6桥梁纵断面图图7桥梁空间布置图表2布设3-36m桥状态计算结果(1)表2布设3-36m桥状态计算结果(2)第三种情况，苏马姆河布设4-36m桥状态：图8桥位平面布置图图9桥梁墩、台冲刷线图图10桥梁纵断面图图11桥梁空间布置图表3布设4-36m桥状态计算结果(1)表3布设4-36m桥状态计算结果(2)对比分析以上三种情况计算结果，第一种情况为河道无桥自然状态，河床底部冲刷粒径取值为10mm，百年一遇洪峰流量对应最大水位高程为10.23m，水流流速为2.97m/s，河床底部无一般冲刷现象；第二种情况为河道布设3-36m桥状态，河床底部冲刷粒径取值为10mm，百年一遇洪峰流量对应最大水位高程为10.90m，水流流速为4.61m/s，河床底部一般冲刷深度为0.81m，桥墩最大局部冲刷深度为3.98m，桥台最大局部冲刷深度为3.73m，桥墩对应总冲刷深度为4.79m，桥台对应总冲刷深度为4.54m；第三种情况为河道布设4-36m桥状态，河床底部冲刷粒径取值为10mm，百年一遇洪峰流量对应最大水位高程为10.49m，水流流速为3.80m/s，河床底部无明显一般冲刷现象，桥墩最大局部冲刷深度为1.7m，桥台最大局部冲刷深度为3.32m。运用河道水力分析模型（HEC-RAS）得出河流水力特征各项数据后，采用国内桥梁冲刷64-1、65-1公式对第三种情况河道布设4-36m桥状态进行计算