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

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

   首页 人人文库网 > 资源分类 > DOC文档下载

水利工程论文-考虑河岸冲刷的弯曲河道水流及河床变形的数值模拟.doc

  • 资源星级:
  • 资源大小:33.85KB   全文页数:11页
  • 资源格式: DOC        下载权限:注册会员/VIP会员
您还没有登陆,请先登录。登陆后即可下载此文档。
  合作网站登录: 微信快捷登录 支付宝快捷登录   QQ登录   微博登录
友情提示
2:本站资源不支持迅雷下载,请使用浏览器直接下载(不支持QQ浏览器)
3:本站资源下载后的文档和图纸-无水印,预览文档经过压缩,下载后原文更清晰   

水利工程论文-考虑河岸冲刷的弯曲河道水流及河床变形的数值模拟.doc

水利工程论文考虑河岸冲刷的弯曲河道水流及河床变形的数值模拟摘要本文将正交曲线坐标系下的平面二维水沙数学模型和粘性河岸的冲刷模型结合,用于模拟弯道内的水流运动、悬移质泥沙的输移、河床的纵向及横向变形.用水槽试验资料验证了本文提出的水流模型,结果表明流速分布、水位等计算结果与实测值相当符合.应用建立的水沙数学模型以及河岸冲刷模型,模拟了一概化弯道在持续清水冲刷下的主流线位置、断面形态、主槽比降的变化过程,模拟结果符合弯道演变规律.关键词弯曲河道平面二维水沙数学模型河岸冲刷模型河床纵向及横向变形弯曲河道是天然河道中常见的一种河流形态,多年来人们从各个方面对弯曲河道特有的水流运动规律、河床演变规律,进行了广泛的研究1,2.如张红武2从模型试验和理论分析出发,较为深入地研究了弯道内的水面形态、环流运动以及纵向流速的沿程分布规律等.随着计算机性能的提高,以及数值计算方法的发展,很多研究者开始用数学方法模拟弯道内的水沙运动及河床变形.最早以Engelund3、Odgaard4等人提出的恒定水流模型为代表.近年来,部分研究者提出的三维数学模型,既能模拟弯道内的水沙运动,又能考虑床面的冲淤变化,较以前有了更大的进步5,6.尽管三维数学模型可以模拟弯道内较为复杂的水流运动和河床变形问题,不过在实际工程中用的最广的还是平面二维模型7,8.通常现有这些模型只考虑河床的纵向变形,而很少注意弯道内的横向变形问题,也即河岸冲刷问题.本文针对当前在弯道水流和河床变形模型的研究状况,以前人的研究结果为基础,建立了正交曲线坐标下的平面二维水沙数学模型,用于模拟弯道内的水沙运动及河床纵向变形在深入分析河岸冲刷机理的基础上,采用力学模型模拟粘性河岸冲刷、崩塌的过程.最后利用弯道模型的试验资料,验证了本模型计算的流场、岸边水位等与实测值符合较好.同时模拟了90°弯道在持续清水冲刷下的河床变形过程,对主流线、断面形态、主槽比降等计算结果随时间变化的分析表明,模拟结果符合弯道变形规律.1数学模型本文建立的模型具有如下特点1为适应不规则的岸边界,建立正交曲线坐标下的平面二维水沙模型,并用MADI方法9求解水流方程组.2改进以往对泥沙恢复饱和系数取为经验常数的做法,模型中采用了张红武10提出的不平衡输沙理论,用于研究悬移质泥沙的输移以及河床变形过程.3引入Osman11提出的河岸冲刷模型,用于模拟粘性河岸的横向冲刷过程以及在重力作用下的崩塌过程.1.1平面二维水沙数学模型1.1.1水流基本方程式中εξ、εη分别表示ξ、η方向的泥沙扩散系数S0k、Sk、Sk、ωk分别为第k粒径组泥沙的侧向输入项、分组含沙量、分组挟沙力及有效沉速α、f1、K1分别为平衡含沙量分布系数,非饱和系数以及附加系数.上述3个参数以及水流挟沙力的计算方法,详见文献10.而分组挟沙力级配以及床沙级配的计算方法,见文献12.1.1.3河床纵向变形方程由悬移质泥沙不平衡输移引起的河床纵向变形方程为式中ρ′为床沙干密度,N为非均匀泥沙的分组数.1.2几个关键问题的处理1.2.1正交曲线网格的生成设x,y为物理平面上的笛卡尔坐标系,ξ,η为相应计算平面上的变换坐标系,用以下方程可实现两个坐标系之间的变换,即求解上述方程组,即可得到物理平面x,y与计算平面ξ,η上对应点的关系.当网格正交时,β0.控制函数的具体表达式,参见文献13.1.2.2初、边界条件对初始条件,一般给出水位、流速、含沙量等变量的初值.对进口边界、给定流量、沙量过程线,以及含沙量级配.对出口边界、给定水位过程线,并认为S/ξ0.对岸边界,对水流可采用滑移或无滑移边界条件对泥沙,可取S/0为岸边界法线方向.在本文模拟的弯道中,横断面由主槽和滩地组成.假设滩地不过水,则岸边界是滩岸,同时认为滩岸可以冲刷.实际计算区域包括滩地和主槽,故采用冻结法技术14处理滩地上那些不过水的节点.一般岸坡坡顶和坡脚间的水平距离较小,可能远小于网格尺寸,在程序中不能分辨.为了能模拟河岸的冲刷过程,同时又不增加计算量,故本文在计算中用一数组记录岸边相邻两节点之间的实际地形.1.3河岸冲刷模型现有的很多泥沙数学模型,很少考虑到河岸冲刷问题.即使那些考虑了河岸冲刷的模型,往往对弯道中河岸冲刷机理、冲刷速率,做了大量的假设和简化.如Ikeda15认为弯道凹岸附近的纵向垂线平均流速大于弯道中心处的垂线平均流速时,弯道凹岸冲刷,反之则淤积.Hasegawa16认为河岸冲刷速率与近岸的剩余流速成正比,根据泥沙连续方程,得出了适用于弯曲河道通用的河岸冲刷系数.不过,这类模型中存在一个共同缺点,即很少考虑到河岸组成物质、土体特性、几何形态等对河岸冲刷的影响.众所周知,粘性河岸和非粘性河岸的冲刷机理和冲刷速率是完全不同的.为准确模拟河岸的横向冲刷过程,必须深入分析河岸冲刷机理.一般来讲,河岸冲刷后退是河岸土体和近岸水流相互作用的结果.除了岸边植被生长情况、河道内水位升降、渗流、管涌等因素影响河岸后退外,但主要是以下2种情况,是导致河岸后退的主要原因1通过水流直接横向冲刷河岸导致河岸后退.近岸水流直接作用于河岸,冲动河岸边坡上水面以下的表层土体,并被水流带走,从而导致河岸后退.2通过河岸崩塌导致河岸后退.崩塌是河岸上的一部分土体在重力作用下,沿某一滑动面发生移动的过程.一般床面发生冲刷,导致河岸高度增加,或者水流淘刷河岸坡脚,使河岸坡度变陡,都会降低河岸的稳定性,当河岸的稳定性降低到一定程度后,河岸便会发生崩塌,导致岸边界后退.图1粘性河岸冲刷的计算模式本文主要考虑粘性河岸的冲刷后退过程,采用Osman11提出的河岸冲刷模型图1.该模型首先计算河岸横向冲刷距离,然后分析河岸是否会失稳、崩塌.在Δtsec时间内,粘性河岸被水流横向冲刷后退的距离为式中γs河岸土体的容重kN/m3ΔB为Δt时间内河岸因水流横向冲刷而后退的距离mτ为作用在河岸上的水流切应力N/m2τc为河岸土体的起动切应力N/m2.Cl为横向冲刷系数,取决于河岸土体的物理化学特性.当由式7得河槽冲宽ΔB,用平面二维水沙模型算出河床冲深ΔZ后,河岸高度增加,坡度变陡,稳定性降低.根据土力学中的边坡稳定性关系,采用若干假定,可得到河岸发生初次崩塌时的临界条件.若河岸已发生初次崩塌,则假定以后的河岸崩塌方式为平行后退,即崩塌后的边坡角度恒为β,仍可用土力学的方法判断是否会发生二次崩塌.1.4计算方法和求解过程首先,采用MADI法9计算流场.这种方法采用非交错网格,将水位、流速等变量均布置在同一网格点上,对水流连续方程和动量方程均不按方向剖开,由此将基本运动方程离散而形成新的差分代数方程组,

注意事项

本文(水利工程论文-考虑河岸冲刷的弯曲河道水流及河床变形的数值模拟.doc)为本站会员(奋斗不息)主动上传,人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知人人文库网([email protected]),我们立即给予删除!

温馨提示:如果因为网速或其他原因下载失败请重新下载,重复下载不扣分。

copyright@ 2015-2017 人人文库网网站版权所有
苏ICP备12009002号-5