基于完全圣维南方程的跨河桥梁设计研究.docxVIP

基于完全圣维南方程的跨河桥梁设计研究 在确定桥梁与桥下的设计流量和桥下的设计速度时，可以发现三种情况。首先，如果支流中有设计流量，则河流没有溢流。其次，如果支流中有正常流动，河流中会发生洪水，并将流入支流。第三，河流回流后，开始撤退，支流发生设计洪水。首先，确定设计桥梁的设计速度和桥下的设计速度与一般河流相同。桥的高度取决于河流回流形成的最大洪峰值，这比普通桥更大。在大洪水的情况下，平台的磨损不同于不倒水。在其次，由于水流回流，河底没有沉积物和沙量运动，桥下无法承受下游沉积物的补充。对于桥下的设计水流，应采用土壤许可和不清洁平均水流。在第三种情况下，由于滨水区内的速度逐渐减少，桥前河底的泥沙运动相对减少或完全停止。桥梁的设计速度应采用无漏水自然速度和土壤允许速度之间的速度。 可见,跨越支流的桥梁,在遇到大河倒灌的情况时,水力计算的任务是确定支流的壅水曲线和大河水位上涨及下落过程中的最大绝对值流速. 1 数学模型 1.1 力要素的等分差差方程 支流受大河倒灌时的水流运动属于明渠缓变不稳定流,其运动规律可用一维圣维南方程组来描述: 式中,Q为流量;Z为水位;B为水面宽;q为单宽旁侧入流量;A为过水断面面积;R为水力半径;n为糙率;α为动量校正系数,一般情况下取为1;g为重力加速度;x为沿河长的距离变量;t为时间变量. 用N+1个断面可将计算河段长L划分成N个子河段,每个子河段的长度为Δxi(i=1, 2, …,N),取时间步长为Δt,则应用Preissmann四点隐式差分格式,并令时段初的水力要素为时段平均水力要素,就可将式(1)和(2)离散成差分方程.对第i个子河段,该差分方程可写为 {?Qj+1i+Qj+1i+1+CiZj+1i+CiZj+1i+1=DiEiQj+1i+GiQj+1i+1?FiZj+1i+1+FiZj+1i+1=φi(3){-Qij+1+Qi+1j+1+CiΖij+1+CiΖi+1j+1=DiEiQij+1+GiQi+1j+1-FiΖi+1j+1+FiΖi+1j+1=φi(3) 式中, 其中,各变量的下标表示断面编号,上标表示时段数;θ为差分格式权重.由式(4)～(9)可知,式(3)中的系数Ci,Di,Ei,Fi,Gi,ue001φi仅与河槽几何参数、糙率和初始条件有关,故它是一个线性代数方程组.N个河段可列出包含2(N+1)个未知变量的2N个方程式,如再加上河道两端的边界条件,就可构成一个闭合的代数方程组.该闭合代数方程组的系数矩阵是一个四对角稀疏矩阵,因此,可用追赶法进行求解,从而求得各断面的水位和流量过程. 式(3)和式(4)～(9)即构成了大河倒灌情况下,确定支流壅水水面曲线和流速的数学模型.该模型包括的参数为河道横断面几何参数、糙率和差分格式权重.河道横断面几何参数一般由地形测量或断面测量成果求得.差分格式权重一般由理论分析确定,已经证明,当θ≥0.6时,差分格式是稳定的.糙率一般根据实测水文资料予以率定,当缺乏实测水文资料时,也可采用近似方法来确定. 对于为数众多的中小河流而言,一般皆缺乏必要的水文和地形资料,使用该法常常遇到困难.根据水力学原理,只需较少资料,就可以近似计算出河道壅水水面曲线. 1.2 流流场物理模型 物理模型设置在河海大学工业冷却水实验室内.该物理模型的支流总长90 m,其中试验段长70 m,矩形断面,宽0.4 m;干流长16 m,宽2.0 m;支流与干流正交,汊口处采用圆弧过渡连接,物理模型的入流有干流和支流两部分;支流采用90°三角堰控制施放不同流量,干流采用微机自动控制尾门调节,以实现干流洪水水位涨落变化的要求.在该物理模型中布置了4个水位测点,分别测量干流及支流离汊口18 m,40 m和65 m处的水位变化.流速测点设置在支流离汊口处1.5 m处,用于测量汊口附近的流速变化.图1是该物理模型实测和数学模型计算的水位过程比较图,从中可看出,两者较为吻合. 1.3 建立模型及结果 尽管可利用上述的水力学方法来推求支流任意位置的流速,但为了便于实际应用以及由于资料限制等原因,这里只采用简单的经验公式方法.因任意位置处的最大流速与空间位置有关,难以用解析式表达,故只能得到干支流交汇处的流速经验公式,对于其他位置处的设计流速,该公式具有重要的参考作用. 为确定大河倒灌情况下,大河水位涨落时支流在交汇处的最大流速,应对影响最大流速的因子进行分析.根据水力学原理,影响最大流速的因子除断面形状外,可归纳为大河水位涨落率、水位涨落历时、支流单宽流量、支流河底比降和支流糙率等.在此基础上,即可建立它们之间的经验关系式.如果断面形状已定,此种经验关系一般可具体表达为 式(10)和(11)中,(Vm)涨,(Vm)落分别为大河水位涨、落时最大流速,大河水位上涨时,(Vm)涨指向支流上游,大河水位