水位暴风雨条件下密集排水孔三维渗流场精细模拟.docxVIP

水位暴风雨条件下密集排水孔三维渗流场精细模拟 事实上，大多数节水项目都会出现带密集排水孔和自由面的非稳定渗流分析。如大坝初期蓄水、水库应急放空、洪水期堤坝水位快速上涨和回落等。由于非稳定渗流场产生的渗透压力及岩土体工程性质的迅速变化,将给工程稳定性增加新的不安全因素。 密集排水孔(幕)是降低坝体浸润线和坝基扬压力比较有效的方法之一。近年来我国越来越多的水利工程采用密集排水孔作为主要排水措施。而在一些高寒区,为了防止坝体冰冻,采用了上游坝坡进行混凝土护坡,同时混凝土护坡预制块布置一定数量的排水孔的方案,以便库水位在应急放空时,坝体内水能够快速地排出。而排水孔孔距、孔径直接影响坝体内渗流场分布,特别是在水位骤降下,不同间排距、不同孔径排水效果明显不同。因此需要对密集排水孔进行精细模拟。对于密集排水孔作用下渗流精细模拟是工程界和学术界关心的问题,排水孔径向孔径较小,数量众多,且空间分布复杂,对每一个排水孔渗流行为进行精细准确模拟是十分困难的事情。目前模拟排水孔的方法主要有以沟代井法、解析解法和有限元结合方法、排水子结构法和改进的排水子结构法。目前国内外采用固定网格的有限元法求解有自由面的渗流场方法主要有结点虚流量法、初流量法、截止负压法和改进的截止负压法。由于不用在每次迭代中确定自由面近似位置和判别自由面单元相交的实际情形等诸多优点,改进的截止负压法是目前固定网格中求解非稳定渗流自由面问题较为有效的方法。 针对大坝上游坝坡布置密集排水孔的护坡工程,有必要采用更加精确的数值方法研究水位骤降下密集排水孔排水效果和坝体内非稳定渗流场变化过程,进而分析非稳定渗流产生的渗透力对大坝稳定性的影响。本文以下坂地沥青混凝土心墙坝为例,采用改进的截止负压法并引入子结构方法和排水孔开关器方法,对坝体上游护坡密集排水孔在水库应急放空情况下的渗流场分布进行了有限元数值模拟研究,详细分析了不同孔距和孔径的排水孔布置方案对大坝渗流场影响,并根据计算结果提出了坝上游混凝土护坡的布置排水孔的合理方案。 1 非稳定渗透理论及渗透解决方案 1.1 分控制方程 非均质各向异性多孔介质中有自由面非稳定渗流连续微分控制方程为: 式中:Q———源或汇项; 自由面边界条件: 初始条件: 式中:h(x,y,z,t)———已知水头函数; θ———自由面外法线方向与垂线的交角。 1.2 lerkin法求解 对上述数学模型在空间域上合理离散,并采用Galerkin法和等参八结点六面体单元,可得到相应的有限元求解格式: 在此基础上,对时间域进行合理离散,可得到有限元时间隐式差分迭代格式: 1.3 计算重复格式 由于自由面位置是未知的,式(6)无法一次求解,这里采用文献和介绍的增量迭代法: 可推导得如下适于计算的迭代格式: 其中: 1.4 成立后的期流场求解方法 本文采用改进的截止负压法对有自由面非稳定渗流场进行求解。文献首先提出了有自由面稳定渗流场的求解方法-截止负压法,文献在文献基础上又对截止负压法进行了改进。详细公式推导见文献。 1.5 用汉字法处理农业生产中的问题 (1)对于已知水头结点,其结点压力为已知,因而其每步迭代求得的压力增量为0,故可将渗透矩阵中该结点对应的对角元素乘以一个大数即可。 (2)对于逸出面上的结点,如果发现其结点压力大于0,可以采用处理截止负压的方法进行类似处理,即将渗透矩阵中该结点对应的对角线元素置为一个大数λ,同时将其对应的结点不平衡力赋为-λp,因此可以迫使该结点的压力等于0。 (3)自由面边界计算时,由于第k步迭代时渗压场是未知的,无法计算当前自由面边界,故用第k-1步迭代时求得的渗压场通过差值得到自由面,将此近似看作第k步计算自由面边界积分时的自由面,随着迭代的进行,这样求得的自由面将逐步逼近真实自由面。对于第1步迭代,自由面边界则为该时段初的自由面。 2 计算排水孔子结构的原则 2.1 设置开关器、定位孔,使其成为排水孔的顶排型 目前水利工程中排水孔有2类排水孔[9~11],一类是顶排型,即渗出水流沿排水孔顶口排出。另一类是底排型,即渗出水流沿孔壁向下流动,从排水孔底口流出。而对于大坝上游坡的密集排水孔在水位骤降时,排水孔由顶排型排水孔转变为底排型并最终逐渐失效。 针对大坝上游坡密集排水孔的具体渗流行为,结合改进截止负压法求解技术,对水位骤降下密集排水孔的处理引入开关器概念,即:在每一个排水孔顶口处设置一个开关器K,在排水孔计算时,先假设开关器打开,孔内边界全部用孔顶口高程作为已知水头边界(即一类边界),每一步迭代后对孔内的结点流量进行甄别。其表达式如下: 式中:kejl———孔周边单元e传导矩阵中第j行第l列元素,并且单元e中局部编号为j的孔周边结点对应结点xi; it———渗流场迭代次数。 此类排水孔计算第k步Δt时段每个排水孔端口结点流量记