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

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

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

水利工程论文-黄河下游平面二维水沙运动模拟的有限元方法.doc

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

水利工程论文-黄河下游平面二维水沙运动模拟的有限元方法.doc

水利工程论文黄河下游平面二维水沙运动模拟的有限元方法摘要本文以水流和泥沙运动规律的研究成果为基础,建立了平面二维水沙模型。用有限元方法导出了本模型的离散方程式,用质量集中和预估校正法处理、迭代求解方程组。以黄河济南河段1976年汛期的洪水演进和河床演变为例从水位过程、流量过程、断面平均流速和最大流速、全域流速场和河床断面冲淤形态等方面,对模型进行了验证模拟计算。数值模拟计算结果与实测或物理模型试验结果符合较好。从而证明了本模型可靠性。关键词黄河下游有限元方法验证1前言对于河道上修建桥涵等跨河工程、引水工程等建筑物后,人们所关心的河势变化、流速分布、河床局部冲淤形态和壅水等问题,一维模型是无能为力的,但可以用二维模型解决。在进行黄河水沙运动的模拟过程中,一些模型往往需要对方程组中的许多参数进行经验处理,而且在处理复杂的河道边界以及方程的离散和求解等方面还有许多问题值得研究。本模型以黄河水利科研究院在泥沙运动规律方面的研究成果为基础,建立模型框架。并选用适用于浑水的群体沉速公式计算泥沙沉速,引入适应于从清水到高含沙的水流挟沙能力计算公式和动床阻力计算公式,克服数学模型参数过多而不能通用的缺陷。对模型研究区域的离散和方程的求解是数值模拟的关键。在离散方法方面,各种离散方法都有其优缺点,而且对每种方法又可分为多种形式,根据各种离散方法的特点和河道形态,为使对区域的离散既能很好地拟合长宽比很大,而且弯曲复杂的河道边界,又能根据河势、主流和水深的变化对局部区域加密细划,本模型选用了有限元法。由于离散后形成的方程组的计算量巨大,用一般计算方法对线性方程组求解不能满足要求,本模型应用质量集中的方法划系数矩阵为三对角矩阵,并用预估校正法处理、迭代求解方程组。大大减少了计算量。2基本方程和定解条件2.1基本方程水流运动的基本方程为12式中Ui为垂线平均流速H为水深Z为水位C为谢才系数。C1/nR1/6水力半径R=H,g,ρ,ε分别为重力加速度,水密度和粘滞系数。f为科氏力系数f2ωsinφ,ω为地球自转角速度,φ为地理纬度。为系数矩阵泥沙运动方程为134其中U为流速ωi为泥沙沉速S和S分别为水流含沙量和挟沙力,f1为泥沙非饱和系数K1为考虑紊流脉动在水平方向产生的扩散作用及泥沙存在产生的附加影响而引入的修正系数,简称为附加系数α为平衡含沙量分布系数,详见1。水流挟沙力计算,采用2中的公式计算水流挟沙力。河床糙率计算,应用3中的糙率计算公式,可以描述水力泥沙因子的变化对摩阻特性的影响。2.2定解条件边界条件。对入流边界,给出水流流速或单宽流量和含沙量过程对出口边界,给出水位过程线或流速单宽流量过程线对固壁边界,法向流速为0,水流沿切线方向流速非0。初始条件。给出在计算的初始时刻地形、流速、水位和含沙量等物理量的初始值。3有限元离散模式的建立和方程求解3.1有限元离散模式的建立有限元网格的选择。根据需要,本次选用三角形常应变单元类型离散研究区域。设整个区域共划分为NE个单元,I个节点,单元节点总体编号为i,i1,2,I。对水沙方程1~4,用Galerkin加权余量法逼近对研究区域进行剖分对单元节点和整体节点分别编号,建立局部节点编号系统和整体节点编号系统,并确定两个编号系统的关系在离散区域的基础上,求出未知量在每个单元上的形函数把形函数代入Galerkin积分表达式进行单元分析,建立局部有限元方程式对所有局部有限元求和、总体合成建立总体有限元方程式,加上本质边界条件,即可得到本问题的有限元方程式。对问题进行有限元分析式δH,δU,δV,δS表示求解变量的变分,Ω表示计算域。对以上二阶导数项利用Green公式分部积分公式,并设任一单元第i节点的平均流速,水位和含沙量及对应的加权函数的形函数值分别为Ui,Vi,Hi,Si和Ui,Vi,Hi,Si。则有将上述表达式代入方程可得有限元方程式MijdZsi/dtPijHUjPzijHVjMijdUj/dtNijUjgP1ijZsjMij/ρεQijRijUjMijdVj/dtNijVjgP2ijZsjMijτyi/ρεQijRijVjMijdSj/dtNijSjεsQijRijSjK1αωMij\f1SS/Hji,j1,2,3,I其中i,j,k1,2,I其中对以上各式,当在同一项中含有两个相同的下标时,就意味着该项表示在单元体内全体同节点编号项的叠加根据以上积分式先对各个单元分析,再叠加全域各单元方程式并使之满足边界条件即得到整体有限元方程式。对上述有限元方程在时间上对变量Z、U、V和S用向前差分格式离散可得到一个II常微分方程组。用常规方法对II阶方程组求解,计算量大,本文用质量集中方法4化方程组的系数矩阵成对角矩阵,方程组就可以容易解出。3.2方程的求解用于预估校正法迭代求解方程组第一时段的计算采用欧拉格式,第二时段采用预估校正格式,即对函数值fn进行预估、校正。用|fnfn|≤ε进行判断,若上式成立,则fnfn,否则令fnfnfn1/2,再由上式求fn进行继续迭代,直到满足精度要求。当求出Ui,Vi,Hi和Si后代入有关方程,即求出冲淤变形的河床高程,完成方程的求解。4模型的验证河段基本概况黄河下游北展滞洪区以南的北店子至后张庄河段,长约30km,是受工程控制的弯曲性河道,河道纵比降约为1,其中北店子至泺口铁桥河段长20km,堤距一般0.7km至1.5km。泺口铁桥至后张庄长约10km,堤距1.5km~2.5km,汛期平均悬移质中值粒径d500.02mm~0.027mm,河床质中值粒径D500.07mm~0.11mm,河床糙率n0.013~0.016。模型计算区域的选择和网格划分。选萨口断面为距泺口水文站上游约20km的北店子处,进口断面距其下游的曹家圈断面、郑家店断面分别约3km和9km出口断面选择在距泺口断面下游约10km的后张庄断面附近。计算区域内共有曹家圈、郑家店、洛口和后张庄四个大断面,两岸边界选择在险工、护滩控导工程和大堤等较稳固的工程所连接的边线上。整个计算区域划分成128个小断面,3728个三角网格单元体,共有2000个节点,三角形沿河宽方向最小长度为30m。一般沿河宽方向网格边长为主槽40m~50m左右,滩地100m左右。网格图如图1所示。验证时段及水沙条件选用1976年汛期8月2日至10月15日及对应的水沙过程,共75天。初始条件假设初始地形、流场、水位和含沙量为一定的常数。边界条件给出入流断面各结点的流速单宽流量和含沙量过程和出口图1计算河段网格图gridsofcalculatedriverreach断断面的水位过程线对两岸边界,根据需要分别设定滑动和不滑动边界条件。验证结果及其分析。图2给出了郑家店断面和泺口断面的水位随时间的变化过程。从图中可以看出,数模计算结果与实测值符合较好。图3给出了泺口断面流量过程与实测流量过程的比较图,可以看出,计算值与实测值基本一致。图2水位比较图Waterstagescomparisons图3流量比较图Comparisonofdischarge表1计算流速与实测流速比较表Compasionofcalculatedandfieldvelocities日期8.118.239.49.69.149.2810.2流量m3/s1550351066107410524038502120平均流速实测2.062.622.642.832.662.412.08m/s计算1.862.502.502.692.532.332.00最大流实测3.043.613.843.863.653.262.86

注意事项

本文(水利工程论文-黄河下游平面二维水沙运动模拟的有限元方法.doc)为本站会员(奋斗不息)主动上传,人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知人人文库网([email protected]),我们立即给予删除!

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

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