连续下降的非稳定渗流分析.docxVIP

连续下降的非稳定渗流分析 在水库主断面的运行中，由于维护等原因，应降低水库的水位，并继续降低水位（尽管下降缓慢）。上游坝板可能存在滑动风险。这是因为水库中的碱性线很高。当水库水位降低时，水库中的饱和水从上游流向上游，导致上游非稳定排水压力。为了准确计算非稳定排水流对上游坝板滑动稳定的影响，在水位下降过程中，各阶段水库的水位、润湿线的位置以及渗透水的压力的确定非常重要。在这项工作中，我们使用自由流量近似来确定每个时期的自由面和等势线，计算不同时期的渗透水压力，以便更好地分析水库上游滑动的稳定性。 1 流量分布有限单元的分析和设计 1.1 有“1+0”因素,有“n”或“n”因素,有“n”或2,3,42,3,45,45,45.2,45.2,3,45.2,45.4;5.2,4,5,45,4,5,45,4,5,45,4,45,4,5,44.5;4.2,45.2,3,45.4;5.2,3,45.2,3,45.2,3,45.2,3,45.2,45.2,4.45.2,45.2,45.2,45.2,45.2,45.4;5.4;5.4;5.4;5.4;5.2,3,45.4;5.2,3,45.4;5.2,3,45.4;5.2,3,45.4;5.2,3,45.4;5.2,3,45.4;5.2,3,45.4;5.2,3,45.2,3,45.2,3,45.2,3,45.2,3,45.2,3,45.2,3,45.2,3,45.2,3,45.2,3,45.2,3,45.2,3,45.2,3,45.2,45.2,45.2,45.2,45.2,45.2,45.2,45.2,45.2,45.4;5.2,45.4;5.2,45.4;5.2,45.4;5.2,45 单元中水头分布函数可表示为: Η=4∑i=1ΝiΗi(1)H=∑i=14NiHi(1) 其中:Νi=14(1+ξ0)(1+η0)(i=1,2,3,4)(2)Ni=14(1+ξ0)(1+η0)(i=1,2,3,4)(2) 并有:x=4∑i=1Νixi(3)x=∑i=14Nixi(3) y=4∑i=1Νiyi(4)y=∑i=14Niyi(4) 1.2 dy+2qhd5的渗流场模式 渗流场剖分成若干单元后, 非稳定渗流泛函式相应地分解为有关单元的泛函之和, 即: Ι(Η)=m∑e=1?12[kx(?Η?x)2+ky(?Η?y)2]dxdy+k∑j=1∫Γ2qΗdΓ(5)I(H)=∑e=1m?12[kx(?H?x)2+ky(?H?y)2]dxdy+∑j=1k∫Γ2qHdΓ(5) 为方便计, 以Ie表示单元e上的泛函, 即: Ιe(Ηe)=?e12[kx(?Η?x)2+ky(?Η?y)2]dxdy+∫Γ2qΗdΓ(6)Ie(He)=?e12[kx(?H?x)2+ky(?H?y)2]dxdy+∫Γ2qHdΓ(6) 对上式右边第一项记为即Ie, 对水头H求偏导并写成矩阵表达式如下: {?Ι1?Η}e=[k]e{Ηi}e(7){?I1?H}e=[k]e{Hi}e(7) 式中: kei?j=∫1-1∫1-1{Bi}Τ[J-1]Τ[Μ][J-1]{Bj}|J|dξdη(8)kei?j=∫1?1∫1?1{Bi}T[J?1]T[M][J?1]{Bj}|J|dξdη(8) 其中: {Bi}Τ=[?Ηi?ξ?Νi?y]?[J]=[?x?ξ?y?ξ?x?η?y?η]?[Μ]=[kx00ky] 对(6)式右边第二项, 表示Γ2边界的流量边界条件, 设m,p点为边界上二相邻点, 对非稳定渗流变动的渗流自由面作为流量补给边界条件时有Ιe2=∫Γ2μΗ?Η?tcosθdΓ=∫xpxmμΗ?Η?tdx, 对H求偏导并写成矩阵表达式 (式中μ为渗流自由面变动区土体给水度): {?Ι2?Η}e=[p]e{?Η?t}e(9) 式中:[p]=Μ(xp-xm)6[0000000000210012] 综合(8)式和(9)式, 对任意单元e有: {?Ι?Η}e=[k]3{h}3+[Ρ]3{?Η?t}e(10) 求泛函极值可得泛函对结点水头微分的方程组, 采用置大数法, 把已知水头的结点也加入方程中写成矩阵形式为: [Κ]{Η}+[p]{?Η?t}={f}(11) 式中: {F}是已知常数项, 对于未知水头为零, 对于已知水头则输入该水头. 对时间项取隐式有限差分则(11)式变为: [Κ]{Η}t+Δt+1Δt[Ρ]{Η}t+Δt-1Δt[Ρ]{Η}t={F}(12) 由上式可知, 已知前一时刻t的结点水头分布, 即可求出下一时刻t+Δt的水头分布. 因此只要知道初始条件下的渗流水头分布, 即可计算库水位降落等边界条件改变时的渗流场水头分布. 1.3 程序的组成 1.3.1 表面位置验算 自由面是渗流场特有的一个待定边界, 本文采用迭代法. 具体迭代步骤如下: