会员注册 | 登录 | 微信快捷登录 QQ登录 微博登录 | 帮助中心 人人文库renrendoc.com美如初恋!
站内搜索 百度文库

热门搜索: 直缝焊接机 矿井提升机 循环球式转向器图纸 机器人手爪发展史 管道机器人dwg 动平衡试验台设计

水利工程论文-非恒定流泥沙数学模型原理及其应用.doc水利工程论文-非恒定流泥沙数学模型原理及其应用.doc -- 2 元

宽屏显示 收藏 分享

资源预览需要最新版本的Flash Player支持。
您尚未安装或版本过低,建议您

水利工程论文非恒定流泥沙数学模型原理及其应用摘要本文介绍了一维非恒定泥沙数学模型的原理和计算方法,并且对于泥沙数学模型中的一些关键技术等问题进行了较详细的讨论。利用已建立的一维非恒定泥沙数学模型,对黄河小北干流1981年洪水资料进行了模拟计算,计算结果表明黄河干流、渭河、北洛河计算水位、流量过程线与实测值符合良好,可用于黄河的洪水演进计算。关键词非恒定流泥沙数学模型河床变形1泥沙数学模型基本方程明渠或天然河流常被考虑作为一维流动,根据洪水波运动的圣维南方程、泥沙连续方程和泥沙扩散方程,可以简化推导出一维非恒定流泥沙数学模型的基本方程,其形式为水流连续方程1水流动量方程2泥沙连续方程3不平衡输沙方程4挟沙力方程SfQ,A,B,ω,S5其中Q为流量,A为断面面积,B为断面宽度,Z为水位,K为流量模数,S为断面平均含沙量,S为断面平均挟沙力,g为重力加速度,α为泥沙非平衡恢复饱和系数,ω为泥沙颗粒沉速,Ad为断面冲淤面积,γ′s为泥沙干容重。2计算方法简介一维非恒定泥沙数学模型的计算采用非耦合方法,首先求解水流连续方程和动量方程,然后求解水流挟沙力、泥沙不平衡输沙方程和泥沙连续方程,具体求解过程如下。2.1水流方程的求解首先利用Preissmann隐式差分格式将水流连续方程和动量方程转化为差分方程,然后对差分方程进行线性化,在线性化过程中,略去增量的乘积项,最后得到以下线性方程组A1jΔQjB1jΔZjC1jΔQj1D1jΔZj1E1j6A2jΔQjB2jΔZjC2jΔQj1D2jΔZj1E2j7其中Aij、Bij、Cij、Dij、Eij(i1,2)为第j单元河段差分方程的系数(j1,2,......,N1,其中N为断面个数)。给定边界条件ΔQ1Q1n1Qn1Q1tn1Q1n8ΔZNZNn1ZNnZNtn1ZNn9方程6、7式及边界条件,共有2N个未知数,2N个方程,可以求解。由于差分方程中的系数包含有未知数,方程求解不能直接求出未知变量,因此方程求解时必须进行迭代处理。下面给出用追赶法求解的步骤,追赶方程为ΔQjFjΔZjGj10ΔZjHjΔQj1IjΔZj1Jj11其中Hj、Ij、Jj、,Fj、Gj为追赶系数。2.2水流挟沙力、动床阻力、河宽变化、床沙级配调整采用文献2的方法2.3不平衡输沙方程求解利用迎风格式,将3式离散为差分方程,整理后得12当Q≥0时,利用上边界条件,自上而下计算各断面含沙量当Q0时,利用下边界条件由下至上计算各断面含沙量。2.4河床变形及淤积量计算由式3与式4相减可以得到河床变形方程为13将上式写成差分方程,各断面淤积面积为14第j河段的淤积量为ΔWjΔAjΔAj1Δxj/2153水流内边界的处理水流内边界是指河道的几何形状的不连续或水力特性的不连续点。例如,河流的汇合点、河流分流、局部河段内生产堤决口等等。在这些内部边界处,圣维南方程组和单一河道泥沙不平衡输沙方程等都不再适用,必须根据其水力特性作特殊处理。内边界条件通常包含两个相容条件即流量的连续条件和能量守恒条件(或动量守恒条件)。本模型主要考虑了以下几个类型的内边界处理。3.1水沙的汇入或汇出如图1所示,假设汇入或汇出点上下断面满足以下条件3.2支流从干流分流如图2所示,干流和支流上断面之间应满足连续方程和能量方程图1水沙的汇入或汇出示意图Sketchofwater/sedimentinfloworoutflow图2支流从干流分流示意图Outflowfromthemainstem图3支流汇入干流示意图Inflowintothemainstem支流从干流分流时,将干流分流断面按干支流流量比分为两部分,忽略时变项,将方程4式直接写成差分形式,求得干流和支流下游断面的含沙量为其中φ为分流系数,由实测资料确定。3.3支流汇入干流当支流汇入干流时,与支流从干流分流类似,干流和支流上断面之间也满足连续方程和能量方程将汇流断面按干支流流量比分为两部分,忽略时变项,将方程4式直接写成差分形式,求得汇流断面的含沙量为其中φ为分流系数,由实测资料确定。4泥沙数学模型的应用4.1计算区域和时段计算区域为龙门、华县、河津、状头至潼关河段,黄河干流上有渭河和汾河汇入,在渭河上有北洛河汇入。本模型同时模拟黄河干流、渭河、北洛河三条河流的洪水演进过程,各汇入点作为内边界处理,汾河仅作为已知水、沙过程线汇NFDA1控制的已知条件,下边界条件为潼关站出口控制水位。计算时段为1981年汛初第一场洪水,洪水时间为7月3号至7月14号,洪水持续时间为12天,龙门最大流量为6400m3/s、最大含沙量为298.0kg/m3华县最大流量为970m3/s、最大含沙量为117.0kg/m3。原始大断面资料采用1981年汛前实测大断面资料。4.2计算结果与实测值的比较图4为潼关站计算流量与实测值的比较、计算含沙量与实测值的比较。由图4可以看出,计算的潼关出口流量过程线与实测过程线比较符合,计算洪峰最大值和相应洪峰传播时间与实测值比较接近计算潼关出口含沙量过程线与实测值也是比较一致的。由图4分析得出,本模型能够比较好地模拟出三条河流的洪水传播过程和泥沙冲淤调整过程。图4潼关站流量、含沙量计算值与实测值的比较ComparisonofcalculatedandmeasureddischargeandconcentrationsatTongguanstation图5华阴站流量、水位计算值与实测值的比较ComparisonofcalculateddischargeandwaterlevelswithmeasuredatHuayinstation图5为华阴站计算水位与实测值的比较、计算流量与实测值的比较。由图中可以看出,计算水位、流量过程线与相应实测过程线比较符合。同时也看出,在洪水初期,华阴站出现了倒灌现象,流量出现了负值,本模型也比较好地模拟出了黄河干流倒灌渭河的现象。图6为朝邑站计算水位与实测值的比较、计算流量与实测值的比较。可以看出,计算水位、流量过程线与相应实测过程线比较一致。由于计算北洛河各断面的水位和流量,是通过一系列内边界条件计算得到的,因此也可以说明本模型所采用的内边界处理技术是比较合理的。
编号:201312122236046768    大小:43.46KB    格式:DOC    上传时间:2013-12-12
  【编辑】
2
关 键 词:
专业文献 建筑水利 精品文档 水利工程
温馨提示:
1: 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
2: 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
3.本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
5. 人人文库网仅提供交流平台,并不能对任何下载内容负责。
6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
  人人文库网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
0条评论

还可以输入200字符

暂无评论,赶快抢占沙发吧。

当前资源信息

4.0
 
(2人评价)
浏览:9次
奋斗不息上传于2013-12-12

官方联系方式

客服手机:17625900360   
2:不支持迅雷下载,请使用浏览器下载   
3:不支持QQ浏览器下载,请用其他浏览器   
4:下载后的文档和图纸-无水印   
5:文档经过压缩,下载后原文更清晰   

相关资源

相关资源

相关搜索

专业文献   建筑水利   精品文档   水利工程  
关于我们 - 网站声明 - 网站地图 - 友情链接 - 网站客服客服 - 联系我们
copyright@ 2015-2017 人人文库网网站版权所有
苏ICP备12009002号-5