心墙土石坝渗流场的无单元法模拟  

Element-free Galerkin method simulation of the seepage field in embankment dam with core wall

作　　者：刘昌军[1] 丁留谦[1] 宁保辉[1] 张顺福[1] 

机构地区：[1]中国水利水电科学研究院,北京100038

出　　处：《水力发电学报》2012年第3期182-187,共6页Journal of Hydroelectric Engineering

基　　金：国家十一五科技支撑课题(2008BAB42B05;2008BAB42B06);国际科技合作与交流基金项目(2010DFA74520)

摘　　要：在心墙土石坝渗流场分析中,经常会遇到浸润线发生在渗透性差别较大的心墙和堆石体之间,这种因两种介质渗透系数差别较大而形成的渗流场采用传统的有限元方法难以较好的解决。为此将无单元法应用于心墙土石坝工程的渗流分析中,提出了心墙土石坝渗流场的无单元法的解决方案,并采用IDL语言编制了可视化的二维无单元法计算软件LIDAREFM。利用该软件计算分析了心墙和堆石体渗透系数不同比值的渗流场分布。研究结果表明,当堆石体渗透系数与心墙渗透系数比值小于103数量级时,可以直接进行计算;而超过103数量级以上,需在粘土心墙和堆石体之间依次设置接触面区域和过渡区域,且接触面区域的渗透系数必须采用空间渗透张量,否则将出现错误的结果或迭代不收敛问题。In analysis of seepage field in embankment dam with core wall,saturation lines often appear between core wall and rock-fill as a result of significant difference in permeability of these two zones,and for the same reason a satisfactory solution of seepage distribution is hard to calculate with conventional FEM.This paper adopts an element-free Galerkin method（EFGM） to analyze such seepage fields and presents a solution scheme.It was implemented through a two-dimensional visualization software,LIDAREFM,that was coded with interactive data language（IDL） and applied to examine the effect of different permeability ratios of the two zones.Results show that if the ratio is not above 1000,a direct calculation can be used,and otherwise contact layers and transverse layers have to be placed between the two zones in a due order for smoother transition.In the latter case,a coefficient tensor must be used for the permeability of contact zone,or incorrect calculation results may be caused by diverged iteration.

关 键 词：水工结构 接触面渗流 无单元法 心墙土石坝 LIDAREFM软件 渗流场 

分 类 号：TU443[建筑科学—岩土工程]