HEC-RAS初步教程-2012演示幻灯片.ppt

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,HEC-RAS,应用,武汉大学水利水电学院,曾玉红,1,What is HEC-RAS?,HEC-RAS,H,ydrologic,E,ngineering,C,enter,(HEC),R,iver,A,nalysis,S,ystem.(RAS),Steady flow water surface profile computation,.,水面线分析,Unsteady flow simulation.,非恒定流量水位过程,Sediment transport/Movable boundary sediment,.,泥沙输移,Water quality,水质模拟，温度，藻类，溶解氧，,BOD,Hydraulic design,桥梁、涵洞、堰、溢洪道等等,模拟河流的水流、泥沙、水质问题,可计算急流、缓流以及混合流,可考虑桥、涵、坝、闸、等构筑物,2,安装及帮助文件,Users manual,用户手册,安装、使用、界面菜单,guide to use HEC-RAS,Hydraulic reference manual,水力手册,理论，参数选择,theory and data requirements,Application guide,应用指南,实例,step by step,examples,wwwhec.usace.army.mil,3,4,5,一、数据要求,一般数据,几何数据,恒定流数据,非恒定流数据,plan data,geometric data,steady flow data,unsteady flow data,quasi-steady flow data,sediment data,water quality data,hydraulic design data,6,1.1 plan data,（一般数据）,project,*,.prj,7,在创建并保存,project,文件（*,.prj,）后，其他,data,文件均会自动以同样的名字保存，但采用不同的后缀来区分各类文件。,“,字母,+,数字,”的,后缀方案：主要的,data,命名规则如下：,One file for each Plan(.P01 to.P99),One Run file for each unsteady flow plan(.X01 to.X99),One Output file for each plan(.O01 to.O99),One file for each set of Geometry data(.G01 to.G99),One file for each set of Unsteady Flow data(.U01 to.U99),1.1 plan data,（一般数据）8,1.2 G,eometric data,Geometric Data,的获取：,在,HEC-RAS,中通过手工输入从头建立,利用,HEC-GeoRAS,等,GIS,扩展在,ArcGIS,中先建立好模型，然后导入,HEC-RAS,中,Geometric Data,的组成：,reach,junction,station,river,Cross-section,河流，河网,河段,汊点,断面,站点,9,每个河网由一个或多个河段（,reach,）相互连接而成，汇合处称为汊点（,junction,）。,Reach,以,reach name,彼此区分，每个,reach,还有一个其所在河流的,river name,。,以,river,为单位（注意不是,reach,）从下游向上游定义距离空间，每个,river,上的位置都对应于一个数值，即,station,，自上游往下游数值由大到小，标示断面所处的位置，但数值大小在具体计算中并无意义。,1.2 G,eometric data,10,11,12,横断面数据,横断面布置：尽量垂直流向。,13,cross section,每个,cross section,都有如下属性：,river name,，,river station,，,reach name,，,description,。,断面形态由一组（,y,z,）坐标给出，,y,轴的正方向为从上游向下游看由左岸到右岸。（,y,z,）坐标的顺序必须满足,y,值非减条件。左岸和右岸位置分别以,y,坐标值给出，且必须对应于某个（,y,z,）坐标对，用于标记河道左岸滩、主槽、右岸滩的分割位置。,在断面上可根据分区给定变化的糙率值,n,。,断面之间的间距是自上游向下游计量的，也就是当前断面与下游相邻断面之间的距离，最下游的断面上间距值没有意义，按零存储。,距离分为左岸距离（,LOB,），中心距离（,Channel,），右岸距离（,ROB,）。,每个断面上还可定义突缩和突扩的能量损失系数，用于计算当前断面与下一断面之间的能量损失。,14,1.3,如何对过水断面进行限定,死水区,ineffective flow area,15,1.3,如何对过水断面进行限定,堤防设置,16,17,1.3,如何对过水断面进行限定,障碍物,18,19,二、,Starting with HEC-RAS,开始一个新的,project,输入（导入）,geometric data,输入,flow data,、,boundary conditions,水力计算,后处理（结果显示与打印）,20,Projects,.Prj project,.R01 run file for steady flow plan,.X01 run file for unsteady flow plan,.O01 output file,.G01 geometry data,.F01 steady flow data,.U01 unsteady flow data,.Q01 quasi-unsteady flow data,.S01 sediment data,.SedCap01 sediment transport capacity,.W01 water quality,21,22,Geometric data,River-reach-junction-storage area-cross-section,River Station:,界定断面的位置，用数字表示，大的数字在上游，小的数字在下游。,Cross-section:River,reach,station and elevation,Hydraulic structure,23,Flow data and boundary conditions,Steady flow,Unsteady flow,24,三、恒定流水面线分析,地面,高程,水深,流速,水头,损失,基于一维能量方程，能量损失由沿程和局部两部分组成。,1,、原理,25,26,流程长,比阻,局部阻力系数,2,、能量损失,27,3,、断面过流能力计算,左岸，主槽，右岸采用不同的参数进行计算，三部分相加获得整个断面的过流能力。,计算可给出各部分流量的分配结果。,28,一般取相同水面深度下的较小流量为计算流量。,按糙率变化的节点进行划分：,边滩上按坐标节点进行划分：,29,4,、糙率系数的选择,30,当主槽边坡大于,5H:1V,，或是主槽有多个糙率的情况，按各部分湿周加权计算得出一个复合糙率：,31,6,、断面平均动能,32,7,、阻力系数,阻力坡降,收缩与扩张断面的损失,沿程：,局部：,33,8,、水面线计算过程,对于缓流，在上游断面假定一个水深；,若是急流，则在下游断面假定一个水深；,基于假定的水位，确定对应的总输水流量和速度水头；,计算损失；,求解,WS,（水面）；,计算值和估计值进行比较，满足要求即停止；,34,9,、临界水深的确定,临界水深：,确定流动是急流还是缓流，便于采用合适的计算方法，并进行流动分析。计算中一般自动计算该项。,断面的总能量：,H=WS+V,2,/(2g),临界水深的计算方法：,抛物线法 对于单个极值,割线法 对于多个极值,假定水面高程，采用上述方程,获得最小的能量值所对应的临界水深。,35,10,、动量方程的应用,能量方程适用于渐变流的情况。当水面从急流,-,缓流，以及缓流,-,急流的情况下，即在急变流的情况下，如渠道底坡的急剧变化如有水跃发生、通过桥梁、堰、河网节点等情况，能量方程不再使用，此时可采用经验公式，但很多情况下必须采用动量方程加以求解。,36,牛顿第二定律：,代入各项力，即得：,压力项,重力项,阻力项,速度变化,HEC-RAS,中所用的动量方程形式：37,11,、恒定流计算的一些假定,流动是恒定的,流动是渐变的,流动是一维的,渠道底坡比较缓，小于,1,：,10,38,四、非恒定流计算,1,、方程,依据：质量守恒，动量守恒,连续性方程,inflow=outflow,动量方程,控制体内动量变化,重力项,阻力项,39,2,、计算方法,隐式有限差分法：,隐式有限差分的一般格式：,时间离散：,空间离散：,40,3,、连续性方程的有限差分格式,根据上述分析：,41,3,、动量方程的有限差分格式,边滩：,主槽：,42,43,44,4,、边界条件,上游边界条件 流量,45,下游边界条件,水位过程线,水位随时间的变化、控制水深,流量过程线,下游有实测流量资料,水位流量过程线,水位和流量的单值函数,正常水深,给定阻力系数由曼宁公式确定,正常水深。一般按渠道底坡进 行估算,46,五、,Inline structure,堰（宽顶堰、尖顶堰以及弧形堰）,溢洪道,闸门（平板闸门，弧形闸门）,47,需要,4,个断面模拟构筑物段水流。,2,号断面反应下游的尾水水位，上游,3,号断面体现来流水头,48,W-,溢流堰的宽度,C-,流量系数,TE,BE,HE,系数,H1:Zu-Zd,A:,闸门面积,自由出流：,淹没出流：,过渡状态：,49,自由出流：,堰流：,过渡状态：,50,