黄河下游平面二维水沙运动模拟的有限元方法

《黄河下游平面二维水沙运动模拟的有限元方法》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、黄河下游平面二维水沙运动模拟的有限元方法摘要：本文以水流和泥沙运动规律的研究成果为基础，建立了平面二维水沙模型。用有限元方法导出了本模型的离散方程式，用质量集中和预估校正法处理、迭代求解方程组。以黄河济南河段1976年汛期的洪水演进和河床演变为例从水位过程、流量过程、断面平均流速和最大流速、全域流速场和河床断面冲淤形态等方面，对模型进行了验证模拟计算。数值模拟计算结果与实测或物理模型试验结果符合较好。从而证明了本模型可靠性。关键词：黄河下游有限元方法验证1前言　　对于河道上修建桥涵等跨河工程、引水工程等建筑物后，人们所关心的河势变化、流速分布、河床局部冲淤形态和壅水等问题，一维模型是无能为

2、力的，但可以用二维模型解决。在进行黄河水沙运动的模拟过程中，一些模型往往需要对方程组中的许多参数进行经验处理，而且在处理复杂的河道边界以及方程的离散和求解等方面还有许多问题值得研究。本模型以黄河水利科研究院在泥沙运动规律方面的研究成果为基础，建立模型框架。并选用适用于浑水的群体沉速公式计算泥沙沉速，引入适应于从清水到高含沙的水流挟沙能力计算公式和动床阻力计算公式，克服数学模型参数过多而不能通用的缺陷。　　对模型研究区域的离散和方程的求解是数值模拟的关键。在离散方法方面，各种离散方法都有其优缺点，而且对每种方法又可分为多种形式，根据各种离散方法的特点和河道形态，为使对区域的离散既能很好地拟合

3、长宽比很大，而且弯曲复杂的河道边界，又能根据河势、主流和水深的变化对局部区域加密细划，本模型选用了有限元法。由于离散后形成的方程组的计算量巨大，用一般计算方法对线性方程组求解不能满足要求，本模型应用质量集中的方法划系数矩阵为三对角矩阵，并用预估校正法处理、迭代求解方程组。大大减少了计算量。2基本方程和定解条件2.1基本方程　　水流运动的基本方程为(1)(2)式中Ui为垂线平均流速；H为水深；Z为水位；C为谢才系数。C=1/nR1/6(水力半径R＝H)，g，ρ，ε分别为重力加速度，水密度和粘滞系数。f为科氏力系数(f=2ωsinφ，ω为地球自转角速度，φ为地理纬度)。为系数矩阵　　泥沙运动方

4、程为[1](3)(4)其中U为流速；ωi为泥沙沉速；S和S*分别为水流含沙量和挟沙力，f1为泥沙非饱和系数；K1为考虑紊流脉动在水平方向产生的扩散作用及泥沙存在产生的附加影响而引入的修正系数，简称为附加系数；α*为平衡含沙量分布系数，详见[1]。　　水流挟沙力计算，采用[2]中的公式计算水流挟沙力。　　河床糙率计算，应用[3]中的糙率计算公式，可以描述水力泥沙因子的变化对摩阻特性的影响。2.2定解条件　　边界条件。对入流边界，给出水流流速(或单宽流量)和含沙量过程；对出口边界，给出水位过程线或流速(单宽流量)过程线；对固壁边界，法向流速为0，水流沿切线方向流速非0。　　初始条件。给出在计

5、算的初始时刻地形、流速、水位和含沙量等物理量的初始值。3有限元离散模式的建立和方程求解3.1有限元离散模式的建立　　有限元X格的选择。根据需要，本次选用三角形常应变单元类型离散研究区域。设整个区域共划分为NE个单元，I个节点，单元节点总体编号为i,i=1,2,…I。　　对水沙方程(1)～(4)，用Galerkin加权余量法逼近；对研究区域进行剖分；对单元节点和整体节点分别编号，建立局部节点编号系统和整体节点编号系统，并确定两个编号系统的关系；在离散区域的基础上，求出未知量在每个单元上的形函数；把形函数代入Galerkin积分表达式进行单元分析，建立局部有限元方程式；对所有局部有限元求和

6、、总体合成建立总体有限元方程式，加上本质边界条件，即可得到本问题的有限元方程式。　　对问题进行有限元分析式δH，δU，δV，δS表示求解变量的变分，Ω表示计算域。　　对以上二阶导数项利用Green公式分部积分公式，并设任一单元第i节点的平均流速，水位和含沙量及对应的加权函数的形函数值分别为Ui，Vi，Hi，Si和U*i,V*i,H*i,S*i。则有将上述表达式代入方程可得有限元方程式MijdZsi/dt=-Pij(HU)j-Pzij(HV)jMijdUj/dt=-NijUj-gP1ijZsj-Mij/ρ-ε(Qij+Rij)UjMijdVj/dt=-NijVj-gP2ijZsj-Mijτy

7、i/ρ-ε(Qij+Rij)VjMijdSj/dt=-NijSj-εs(Qij+Rij)Sj-K1α*ωMij[(f1S-S*)/H]ji,j=1,2,3,…I其中i,j,k=1,2,…I其中对以上各式，当在同一项中含有两个相同的下标时，就意味着该项表示在单元体内全体同节点编号项的叠加；根据以上积分式先对各个单元分析，再叠加全域各单元方程式并使之满足边界条件即得到整体有限元方程式。　　对上述有限元方程在时间上对变量Z、