黄河下游平面二维水沙运动模拟的有限元方法(李东风,张红武,许雨新,张俊华)

热度115票 浏览46次 【共0条评论】【我要评论 时间:1999年7月01日 14:13

黄河下游平面二维水沙运动模拟的有限元方法水利论文JO3S5YC$pp SN

李东风张红武,许雨新,张俊华

uA,S+]7HBXA/a0

(黄委会黄河水利科学研究院)

\I/ZbX7qc3}`a0

摘要:本文以水流和泥沙运动规律的研究成果为基础,建立了平面二维水沙模型。用有限元方法导出了本模型的离散方程式,用质量集中和预估校正法处理、迭代求解方程组。以黄河济南河段1976年汛期的洪水演进和河床演变为例从水位过程、流量过程、断面平均流速和最大流速、全域流速场和河床断面冲淤形态等方面,对模型进行了验证模拟计算。数值模拟计算结果水利论文SQ^*K/H Cg?$J

与实测或物理模型试验结果符合较好。从而证明了本模型可靠性。

:n+E*Hq;K0

关键词:黄河下游; 有限元方法; 验证水利论文_sv!p QQj

基金项目国家自然科学基金重大项目资助(59890200)水利论文`f:v.h1i9JEH`3l

作者简介:李东风(1964-),男,河南省商丘市人,黄河水利科学研究院高级工程师。水利论文W!l} g8f.Gy&I/G'}

1 前言

n r&jD4g6]0

  对于河道上修建桥涵等跨河工程、引水工程等建筑物后,人们所关心的河势变化、流速分布、河床局部冲淤形态和壅水等问题,一维模型是无能为力的,但可以用二维模型解决。在进行黄河水沙运动的模拟过程中,一些模型往往需要对方程组中的许多参数进行经验处理,而且在处理复杂的河道边界以及方程的离散和求解等方面还有许多问题值得研究。本模型以黄河水利科研究院在泥沙运动规律方面的研究成果为基础,建立模型框架。并选用适用于浑水的群体沉速公式计算泥沙沉速,引入适应于从清水到高含沙的水流挟沙能力计算公式和动床阻力计算公式,克服数学模型参数过多而不能通用的缺陷。水利论文}(a-J Q$`A T ~#}*r Cw

  对模型研究区域的离散和方程的求解是数值模拟的关键。在离散方法方面,各种离散方法都有其优缺点,而且对每种方法又可分为多种形式,根据各种离散方法的特点和河道形态,为使对区域的离散既能很好地拟合长宽比很大,而且弯曲复杂的河道边界,又能根据河势、主流和水深的变化对局部区域加密细划,本模型选用了有限元法。由于离散后形成的方程组的计算量巨大,用一般计算方法对线性方程组求解不能满足要求,本模型应用质量集中的方法划系数矩阵为三对角矩阵,并用预估校正法处理、迭代求解方程组。大大减少了计算量。

L x2zr5~0

2 基本方程和定解条件

tqIX0[hk0

2.1 基本方程水利论文RcIf;dSK8c

  水流运动的基本方程为水利论文)M)G _"FU"b

水利论文 q&f_F)i!Yp8F C

(1)水利论文iTdGie4M

:?;|J ? k*W'kH:v0

(2)水利论文1Z?4kp/`%m

式中Ui为垂线平均流速;H为水深;Z为水位;C为谢才系数。 C=1/nR1/6(水力半径R=H),g,ρ,ε分别为重力加速度,水密度和粘滞系数。f为科氏力系数(f=2ωsinφ,ω为地球自转角速度,φ为地理纬度)。为系数矩阵

  泥沙运动方程为[1]

aUiqN6E [C)x0

${b/h ]$QGe&ts;ay0

(3)水利论文9~cf7iq}L

;d6CJ6d'b D&i5vrC,Uo0

(4)

Zd4d#O8PF0o0

其中U为流速;ωi为泥沙沉速;S和S*分别为水流含沙量和挟沙力,f1为泥沙非饱和系数;K1为考虑紊流脉动在水平方向产生的扩散作用及泥沙存在产生的附加影响而引入的修正系数,简称为附加系数; α*为平衡含沙量分布系数,详见[1]。

8ZMIc(~0

  水流挟沙力计算,采用[2]中的公式计算水流挟沙力。水利论文-cJI\SL i

  河床糙率计算,应用[3]中的糙率计算公式,可以描述水力泥沙因子的变化对摩阻特性的影响。

As2Iv kk0

2.2 定解条件水利论文,d AV[*? Z3[7]

  边界条件。对入流边界,给出水流流速(或单宽流量)和含沙量过程;对出口边界,给出水位过程线或流速(单宽流量)过程线;对固壁边界,法向流速为0,水流沿切线方向流速非0。水利论文 f6p*E8LN

  初始条件。给出在计算的初始时刻地形、流速、水位和含沙量等物理量的初始值。

&Y-j+q oN%u-w0

3 有限元离散模式的建立和方程求解

Bb)P,C C v;o0

3.1 有限元离散模式的建立水利论文&|-aew]Y

  有限元网格的选择。根据需要,本次选用三角形常应变单元类型离散研究区域。

c(m?,YNn mp%R0

设整个区域共划分为NE个单元,I个节点,单元节点总体编号为i, i=1, 2,…I。

}n.Wf1u%gfX0

  对水沙方程(1) ~(4),用Galerkin加权余量法逼近;对研究区域进行剖分;对单元节点和整体节点分别编号,建立局部节点编号系统和整体节点编号系统,并确定两个编号系统的关系;在离散区域的基础上,求出未知量在每个单元上的形函数;把形函数代入Galerkin积分表达式进行单元分析,建立局部有限元方程式;对所有局部有限元求和、总体合成建立总体有限元方程式,加上本质边界条件,即可得到本问题的有限元方程式。

I.M }])cx4R0

  对问题进行有限元分析水利论文4PzfH1V

式δH,δU,δV,δS表示求解变量的变分,Ω表示计算域。

)D8d6b9{1}{2i&Q7{;{0

  对以上二阶导数项利用Green公式分部积分公式,并设任一单元第i节点的平均流速,水位和含沙量及对应的加权函数的形函数值分别为Ui,Vi,Hi,Si和U*i,V*i,H*i,S*i。则有

m V8T{.b$V0

水利论文![^;v/Yx5lo

e9z-@9A p$K4`eX$n Q0

f-M{p/TNs1x$V6N0

水利论文 U(YR:?p)K5a X^3YS

水利论文 G\1` BW C qY

DWb$JT|W0

水利论文om,G^Lm*l6^

E f%RI7RO:N)k0

将上述表达式代入方程可得有限元方程式

e3E/h%CO0
MijdZsi/dt=-Pij(HU)j-Pzij(HV)j
MijdUj/dt=-NijUj-gP1ijZsj-Mij/ρ-ε(Qij+Rij)Uj
MijdVj/dt=-NijVj-gP2ijZsj-Mijτyi/ρ-ε(Qij+Rij)Vj

MijdSj/dt=-NijSjs(Qij+Rij)Sj-K1α*ωMij\[(f1S-S*)/H]j i, j=1, 2, 3,…I

1p"lp Iz9d*VF0

其中水利论文pYkO6_&z ||a

^&UWz1y l br0

qo0_b~gNK i%J0

990411e7d.gif (1320 bytes)水利论文:]1ik$[i8^"|8mN

水利论文S d _:uL3RY^

水利论文U6@\h/w.@v*dA

o?5rG;q u0

g,N6@HFI"q^0

i,j,k=1,2,…I水利论文x5P z~8rT

其中

对以上各式,当在同一项中含有两个相同的下标时,就意味着该项表示在单元体内全体同节点编号项的叠加;根据以上积分式先对各个单元分析,再叠加全域各单元方程式并使之满足边界条件即得到整体有限元方程式。

R x)?"^h|_*bt0

  对上述有限元方程在时间上对变量Z、U、V和S用向前差分格式离散水利论文eFHT&e#U G-q0ft

2wyt:u*u6Yab0

水利论文 R-C6R1MOP

+bH'e5c Ly,t'Y0

水利论文6qdO5_/]5|

可得到一个I×I常微分方程组。用常规方法对I×I阶方程组求解,计算量大,本文用质量集中方法[4]化方程组的系数矩阵成对角矩阵,方程组就可以容易解出。

{#qcc8mE"n(WF2Ml0

3.2 方程的求解水利论文U9FF9Q$hK

  用于预估校正法迭代求解方程组水利论文 } z&i!}0ex&V

   第一时段的计算采用欧拉格式,第二时段采用预估校正格式,即对函数值fn进行预估、校正。用 |fn-f*n|≤ε进行判断,若上式成立,则fn=f*n,否则令fn(fn+fn-1)/2,再由上式求f*n进行继续迭代,直到满足精度要求。当求出Ui, Vi, Hi和Si后代入有关方程,即求出冲淤变形的河床高程,完成方程的求解。

4 模型的验证水利论文1m.HVLQr o'J _

  河段基本概况 黄河下游北展滞洪区以南的北店子至后张庄河段,长约30km,是受工程控制的弯曲性河道,河道纵比降约为1,其中北店子至泺口铁桥河段长20km,堤距一般0.7km至1.5km。泺口铁桥至后张庄长约10km,堤距1.5km~2.5km,汛期平均悬移质中值粒径d50=0.02mm~0.027mm,河床质中值粒径D50=0.07mm~0.11mm,河床糙率n=0.013~0.016。

/t?{xm*v6sG&D0

模型计算区域的选择和网格划分。选萨口断面为距泺口水文站上游约20km的北店子处,进口断面距其下游的曹家圈断面、郑家店断面分别约3km和9km;出口断面选择在距泺口断面下游约10km的后张庄断面附近。计算区域内共有曹家圈、郑家店、洛口和后张庄四个大断面,两岸边界选择在险工、护滩控导工程和大堤等较稳固的工程所连接的边线上。整个计算区域划分成128个小断面,3728个三角网格单元体,共有2000个节点,三角形沿河宽方向最小长度为30m。一般沿河宽方向网格边长为主槽40m~50m左右,滩地100m左右。网格图如图1所示。

D2pe$W1U;K0
    验证时段及水沙条件 选用1976年汛期8月2日至10月15日及对应的水沙过程,共75天。水利论文7^7p4t Ll:^q

    初始条件 假设初始地形、流场、水位和含沙量为一定的常数。

f;ru})z&w3?,@0

:J8g|C.@q0    边界条件 给出入流断面各结点的流速(单宽流量)和含沙量过程和出口断断面的水位过程线;对两岸边界,根据需要分别设定滑动和不滑动边界条件。水利论文&Dj1_kw

990411t1.gif (1971 bytes)

u Vl7lEV&cT g~0

图1 计算河段网格图

UW(a!g'AX0

grids of calculated river reach水利论文j5p1|%u\:A~*}C

  验证结果及其分析。水利论文M9B6K[9uU!h

  图2给出了郑家店断面和泺口断面的水位随时间的变化过程。从图中可以看出,数模计算结果与实测值符合较好。水利论文)WrD {%L-f

  图3给出了泺口断面流量过程与实测流量过程的比较图,可以看出,计算值与实测值基本一致。

D$k7gC8S)b0

990411t2.gif (7888 bytes)

q8\:I:fP)g~0

图2 水位比较图

j*UGf?0I0

Water stages comparisons水利论文+\Z(k _#};Jew9BFb

990411t3.gif (9160 bytes)水利论文^ v r9V+r0J#]%c

图3 流量比较图水利论文m sjt[ yJ

Comparison of discharge水利论文 GM/r%@%H[

表1 计算流速与实测流速比较表水利论文S7V#zI8f.u:MH6{Y i

Compasion of calculated and field velocities

L'mS^hK+A^0

日期8.118.239.49.69.149.2810.2

流量(m3/s)1550351066107410524038502120
平均流速实测2.062.622.642.832.662.412.08
(m/s)计算1.862.502.502.692.532.332.00
最大流速实测3.043.613.843.863.653.262.86
(m/s)计算2.763.403.563.723.603.192.70

990411t4.gif (3808 bytes)水利论文6[V m2N/yN0w*d

图4不同流量流场图

ku:ZQ3W/O}(|0

Flow fields with different dischsrge

4u1t6dS;]7E0

  表1给出了一些流量的泺口断面最大流速和最小流速的计算与实测值的比较,从图4给出了流量为2600m3/s和7800m3/s下的流场的计算结果; 图4给出了汛前小水、洪峰期大水全域流速分布,从这些图中可以看出,不同流量下河势、主流线的变化、水流漫滩、流速在主槽和滩地的不同分布情况。

Yb*Z)p9_[,C0

990411t5.gif (9522 bytes)水利论文1i c@a*v-h&{!d

图5 横断面冲淤变形比较图水利论文:ho~ mO f1` T3C

Cross-sectional erosion and deposition comparisons

%{v"i2e$hx0

  图5给出了郑家店、泺口重要断面的1976年汛前6月10日实测地形、汛后10月15日数学模型计算地形与实测地形及物理模型试验的对照结果。从图中可以看出,数模计算结果在其冲淤部位和冲淤面积等方面基本上与实测数据或物理模型试验数据相符合。水利论文$W J.?`i-S8oT Z

5 结语

mr kE:P_&D&F0

  1. 模型采用了黄河泥沙研究的新成果,理论基础可靠。水利论文H3XH:L2D

  2. 采用有限元法对区域进行离散,容易处理和更好地拟合不规则河道边界;网格划分灵活,可以对局部区域任意加密;有限元程序模块可移植性强,很容易进行不同网格形状单元离散模式之间转换。

vs$hH?'|[0

 3. 验证结果表明,在水位、流量、流速和河床冲淤等方面,数模计算结果与实测数据或物理模型试验结果符合较好,从而证明了模型的可靠性。

T ? G"mIX9Y0

参 考 文 献

mg_ Dy:F"l0

[1] 张红武,吕昕。 弯道水力学。 水利电力出版社,1993.水利论文[H8o\??

[2] 张红武,张清。黄河水流挟沙力计算公式。人民黄河,1992,(1).

3BXY#_fa0

[3] 赵连军,张红武。紊流阻力系数的计算。′95全国水动力学研讨会文集。海洋出版社,1995.

`7` vCr4J0Z0

[4] 李东风,赖国璋。 串列双圆柱绕流的有限元模拟。力学在实践中的应用。中国林业出版社, 1992.5.水利论文 b])p,uT)C:@ q"S mF

[5] 谢贻权,何福保。弹性和塑性力学中的有限元方法。机械工业出版社,1981.水利论文 b ?:jZ(rYw

 水利论文@o.g2nW;TV$f]

TAG: 黄河 李东风 有限元 张红武
顶:11 踩:16
【已经有90人表态】
12票
极差
8票
很差
10票
较差
12票
稍差
9票
稍好
13票
较好
16票
很好
10票
极好
下一篇:小浪底水库运用后三门峡水库运用方式中一些问题的初步研究(曲少军 钱意颖 郭全明)
上一篇:倒虹吸工程河段洪水与河床变形的数值模拟(刘培斌,孙东坡,杨明)
查看全部回复【已有0位网友发表了看法】