地下工程渗流场有限元模拟研究-地质力学学报.PDFVIP

第 卷 第 期 地 质 力 学 学报 ’ ?@ * ’ 5? * 年 月 % ! ( 12345/6 27 8.29.:;/5:= =AB , % ! 文章编号： （ ） !#$##!# %! $%’!$’ 地下工程渗流场有限元模拟研究 ! % 李 杰 ，李 强 （ 福建省建筑科学研究院，福建 福州 ； !* )% % * 长江水利委员会汉江集团，湖北 丹江口 ++!( , ） 摘 要：在水电、石油、交通等部门的地下工程中，经常会遇到地下水的渗流问 题，如果处理不当，可能诱发应力恶化而影响结构稳定。本文以非饱和介质的渗流 模型为基础，实现了地下洞室开挖过程中的非稳定渗流场解析。 按分期施工方案，在保持渗流场变化连续性的前提下，对各施工阶段进行了准 稳定渗流分析，并在此基础上建立了地下工程非稳定渗流数值模型。其中，对于洞 室系统，基于水量平衡原理，采用了子结构模型；对于防渗帷幕，采用了隔水层模 型；对于排水井列，采用了杆单元模型。 根据上述模型，采用有限差分法作为数值模拟方法，编写了地下工程非稳定渗 流场分析的有限元程序，实现了复杂渗流区域的渗流场的动态模拟。工程实例计算 分析表明，用该模型计算的结果是合理的。 关键词：模拟；地下工程；非饱和；子结构 中图分类号：-.!% 文献标识码：/ 引言 三维非稳定渗流问题数值分析方法以往一般是以自由水面为边界，在饱和区内进行计 算。这种方法在计算中每段都要试求自由水面边界，比较麻烦，而且也无法考虑渗流场内的 初始含水量分布，因此不能全面真实地反映地下水的渗流动态。本文基于饱和0 非饱和介质 渗流分析模型，采用子结构模拟洞室系统，按照施工分期方案，对各施工阶段进行准稳定渗 流分析，较好地实现了对帷幕井列的模拟。伴随施工进度需要保持渗流场变化的连续性，这 样既能正确模拟非稳定渗流场又能反映施工方案对渗流的影响，同时，还可以大大简化非稳 定渗流的复杂计算。 ! 渗流基本理论 饱和 非饱和渗流模型考虑了地下水的非饱和流动，把饱和区与非饱和区耦合在一起进 0 收稿日期：%!$($! 作者简介：李杰 （ —），男，硕士，主要从事流固耦合数值模拟研究。 !(’) ; 地 质 力 学 学 报 ## 行整体分析，压力水头在饱和区为正值，在非饱和区为负值，零压力面就是自由水面，即饱 和区与非饱和区的分界面。 ! ! 基本方程 ［］ 假设水不可压缩，在无内源的情况下，三维饱和 ! 非饱和渗流的基本方程可表示为 ： # # # # ! ! ! ! ! ! ! （） （） （） （） （） ! ( ! # $