HEC-RAS原理0001

1、12345678910111213141516171819202122【精品文档】HEC-RAS1序一维恒定流计算原理1 一维恒定流计算能量方程原理一维恒定流水面线可通过求解能量方程来获得，具体表达式如下:Z2丫22V222g1Vl2Yrhe2g(1)式中：Z,乙为河道底高程；Y, Y2为断面水深;V, V2为断面平均流速;2为动能修正系数；g为重力加速度;he为水头损失。两个断面间的水头损失包括沿程水头损失和局部水头损失，水头损失表达式如下:he LSf C2V222giVi22g式中：L为断面平均距离；Sf为两断面间沿程水头损失坡度；C为收缩或扩散损 失系数。断面平均距离表达式如下:L L

2、iobQebLQch LrobQrobQlob Qch Qrob式中：Llob , Lch, Lrob分别是两断面间左边滩地、主槽、右边滩地的距离；Qlob , Qch ,Qr0b分别是左边滩地、主槽、右边滩地平均流量。根据不同糙率分界点划分滩地，利用曼宁公式计算每个分区的流量，表达 式如下：【精品文档】-仅供参考学习与交流23242526272829303132333435363738394041424344454647【精品文档】Q KS1f/2(4)K1AR2/3n式中：K为流量模数；n为曼宁糙率系数；A为分区面积；R为水力半径。动能修正系数 可通过滩地和主槽流量来进行计算，表达式如下：

3、3 3 32 KlobKchKrobAt. 2. 2. 2AobAchArobK?式中：A为整个过流断面面积；Aob, Ah, Aob分别为左边滩地、主槽、右边滩地 过流面积；K为整个过流断面的流量模数；Kob, Kh, Kob分别为左边滩地、主槽、 右边滩地流量模数。沿程水头损失坡度Sf可通过下式求解:2SfQiQ2Ki K22 一维恒定流计算动量方程原理当水面线越过临界水深，能量方程已经不再适用。自然界中从缓流过渡到 急流或者从急流过渡到缓流现象十分普遍，例如在比降变化较大的渠道、桥梁 堰、河流汇合处，都可能发生这种现象。因此，需要采用动量方程对水面线进 行求解。动量方程是通过牛顿第二定律

4、推导得来的，表达式如下：Fx maP2 Pi Wx Ff Q Vx式中：R、R分别为断面1和断面2上所受的水压力；W为水重力沿X轴分量;Ff为摩擦力；Q为流量；为水的密度； Vx为流速变化沿X轴分量。【精品文档】-仅供参考学习与交流【精品文档】其中:P aYiP2 A2Y2式中： 为水的容重；A1、人分别为断面1和断面2的过水面积；Y1、Y2分别为断面1和断面2水面至形心的距离WxA1A22LS0式中：L为断面1至断面2距离沿X轴方向的分量；So为河底比降。FfA A22SfLmaQ1V1g式中：为动量修正系数。将式（10）式（14）带入式（8）可得:需 A2Y2Ai A2LS0AiA2LSf

5、Q； 1-1A1Y1gA(15)3桥梁建筑物模拟方法3.1断面布设为了准确计算水流通过桥梁建筑物的能量损失，一般在桥梁附近布设断面，断面布设示意图见图148495051525354555657585960616263646566676869【精品文档】-仅供参考学习与交流【精品文档】LeLcBA/理想水流扩散区域转变模式标准水流转变模式707172收缩区域图1桥梁附近断面布设示意图断面1。断面1应布设在桥梁下游水流不受桥梁影响的位置，扩散比73(2Le/AB CD )主要受建筑物结构、河道糙率以及河底比降的影响，表1给项目比降糙率nob/ nc=1not/ nc=2nob/ nc=4b/B=0

6、.10.000191.4-3.61.3-3.01.2-2.10.0009481.0-2.50.8-2.00.8-2.00.0018961.0-2.20.8-2.00.8-2.0b/ B=0.250.000191.6-3.01.4-2.51.2-2.00.0009481.5-2.51.3-2.01.3-2.00.0018961.5-2.01.3-2.01.3-2.0b/B=0.50.000191.4-2.61.3-1.91.2-1.474 出了不同工况的扩散比，以供参考。75 表1扩散比参照表【精品文档】-仅供参考学习与交流767778798081828384858687888990919293

7、949596979899【精品文档】0.0009481.3-2.11.2-1.61.0-1.40.0018961.3-2.01.2-1.51.0-1.4表中b/B为桥下过水宽度与河道宽度比值；S为河底比降；nob/nc为阻水建筑物断面糙率与主槽糙率比值。上表中扩散比只能作为模拟计算的初值，需要根据计算所得的弗劳德数进行修正，当河道宽度接近305m,桥下过水宽度在30.5m152.5m,流量在142nVs851nVs之间，桥下游扩散段长度可用下面的回归方程求解：Fc2八Le298 257 c20.918/ 0.00479Q（16）式中：Le为桥下游扩散段长度（ft ）； Fc2、Fc1分别为为断

8、面2处和断面1处主槽 弗劳德数；Lobs为所有阻水建筑物长度的一半（ft ）； Q为流量（cfs ）。如果河道宽度和流量都小于上述范围，Le可用下式进行计算：” 0.421 0.485 10.000018QLobsFc1(17)如果河道宽度和流量大于上述范围，Le可用下式进行计算:- 0.489 0.608 F2LobsFci（18）扩散系数范围最大不超过4:1，最小不能小于0.5:1。如果扩散系数大于3:1 , 则需要在断面1和断面2之间加一个断面，以便准确计算能量损失。断面2和断面3。断面2和断面3分别位于桥下游和桥上游距离桥梁较近的位置上，断面2应布设在水流经过桥梁后突然扩散的断面上，而

9、断面 3应布设 在水流突然收缩的断面上。断面4。断面4应布设在桥梁上游水流流线平行的断面，一般来讲，桥上游收缩段长度Lc要小于Leo桥梁上游收缩段距离Le与河底比降、糙率以及流量有关，表 2是收缩比参【精品文档】-仅供参考学习与交流100101102103104105106107108109110111112113114【精品文档】照表。当过水断面宽度和流量等级在上述范围内,Lc可用下面的回归方程求解。(19)Fc2Lc 263 38.8 *Fc10.5nob58.70.161L0bsnc式中：Qb为滩地过流流量。如果过水断面和流量等级不在上述范围内，应该使用下式来估算2Lc值。(20)- 1

10、.4 0.333 F21.86 -QobLobsFc1Q0.50.19 nc收缩比不能超过2.5:1 ,不能小于0.3:1 。3.2扩散系数和收缩系数扩散系数和收缩系数是用来描述断面形状发生变化而产生的能量损失， 般情况下，扩散系数要大于收缩系数，下表是扩散系数和收缩系数的参照表。 表3收缩和扩散系数参照表工况收缩系数扩散系数断面形状未发生变化0.00.0渐变断面0.10.3典型桥梁断面0.30.5【精品文档】-仅供参考学习与交流表2收缩比参照表比降糙率nob/ nc=1nob/ nc=2nob/ nc=40.000191.0-2.30.8-1.70.7-1.30.0009481.0-1.90

11、.8-1.50.7-1.20.0018961.0-1.90.8-1.40.7-1.2【精品文档】突然变化断面0.60.8为了准确描述能量损失，需要对两个系数进一步细化。扩散系数表达式如 下:Ce0.09 0.57 -Dob0.075 F2DcFd(21)式中：4为断面1处滩地水力半径；口为断面1处主槽水力半径。使用时可先 参照表3选择初始值，通过计算后可得到 ce的计算值，如果初始值和计算值相 差不大，则G值满足要求，如果相差过大，则重新调整该值，直至二者误差满 足要求为止。收缩系数没有自身的回归方程，但可通过建筑物的结构进行选取，下表为 推荐的收缩系数取值表。表4收缩系数推荐表压缩比推荐的收

12、缩系数0b/ B0.250.3-0.50.25 b/ B0.50.1-0.30.5b/B10.13.3死水区域划分由于引桥阻水，使得桥前桥后有部分死水区域，在计算过程中需要将其划分出去。断面3死水区域划分方法，一般是将对应桥梁过流断面两侧分别加上 断面3与桥梁之间的距离，此范围之外的都为死水区域。断面 2死水区域的划 分方法与断面3类似。死水区域划分的高程与桥梁发生堰流高程相关，对于断面2,开始时刻是不知道堰流发生的高程阈值，需要估计一个高程，这个高程一般介于桥底板高程 与桥面最低点高程之间。断面 3死水区域高程可设为桥面最低点高程。115116117118119120121122123124

13、125126127128129130131132133134135【精品文档】-仅供参考学习与交流136137138139140141142143144145146147148149150151152153【精品文档】4过桥水流一维数值模拟3 BUBD 2图2桥梁附近计算断面位置断面2与断面3之间的水流流态十分复杂，能量方程不能准确计算能量损 失，动量方程是较好的选择。动量方程的运用分三步，第一步是在断面2和断面BD之间建立动量方程，表达式如下：2 2 -AbdYbdbdQbd A2Y2 Q2 AP Ypq Ff Wxg、DgA2 BD式中：A、Abd分别为断面2和断面BD过水面积；Ap为桥梁下游侧桥墩阻水面 BD积；Ypbd为水面至桥墩阻水断面形心的距离;2、bd分别为断面2和断面BD的动量修正系数；Q、Qd为流量；第二步，在断面BD与断面BU之间建立动量平衡方程，表达式如下:AbuYbubuQ；ugABUAbdYbdQ2QDFfWxgABD最后一步，在断面BU和断面3之间建立动量方程，表达式如下:A3Y32 3Q3gA22BU Qbu1APBU Q3AbuYbuQBU%Ypbu1 Cd Ff WxgABU2gA3(25)【精品文档】-仅供