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

第 31卷 第 3期 水 力 发 电 学 报 V01．3l No．3 2012年 6月 JOURNALOFHYDROELECTRIC ENGINEERING Jun．，2012 心墙土石坝渗流场的无单元法模拟 刘 昌军，丁留谦 ，宁保辉，张顺福 (中国水利水电科学研究院，北京 100038) 摘 要 ：在心墙土石坝渗流场分析 中，经常会遇到浸润线发生在渗透性差别较大的心墙和堆石体之间，这种 因两种介质 渗透系数差别较大而形成的渗流场采用传统的有限元方法难 以较好的解决。为此将无单元法应用于心墙土石坝工程的 渗流分析中，提出了心墙土石坝渗流场的无单元法的解决方案，并采用 IDL语言编制了可视化的二维无单元法计算软件 LIDAREFM。利用该软件计算分析 了心墙和堆石体渗透系数不 同比值 的渗流场分布 。研究结果表 明，当堆石体渗透系 数与心墙渗透系数 比值小于 10 数量级时，可 以直接进行计算 ；而超过 1O数量级 以上 ，需在粘土心墙和堆石体之间依次 设置接触面区域和过渡 区域 ，且接触面区域的渗透系数必须采用空间渗透张量 ，否则将 出现错误的结果或迭代不收 敛 问题 。 关键词 ：水工结构 ；接触面渗流 ；无单元法 ；心墙土石坝 ；LIDAREFM软件 ；渗流场 中圈分类号：Tu443 文献标识码 ：A Element-freeGalerkinmethodsimulationoftheseepagefield inembankmentdam withcorewall LIUChangjun，DINGLiuqian，NINGBaohui，ZHANGShun~ (ChinaInstituteofWaterResourcesandHydropowerResearch，Beijing 100038) Abstract：In analysisofseepagefield in embankmentdam with corewall，saturation linesoften appear betweencorewallandrock-fillasaresultofsignificantdifferenceinpermeabilityofthesetwozones，andfor thesamereasonasatisfactorysolutionofseepagedistributionishardtocalculatewithconventionalFEM ．This paperadoptsanelement-freeGalerkinmethod(EFGM)toanalyzesuchseepagefieldsandpresentsasolution scheme．Itwasimplementedthroughatwo-dimensionalvisualization software，LIDAREFM ，thatwascoded withinteractivedatalanguage(IDL)andappliedtoexaminetheeffectofdifferentpermeabilityratiosofthe twozones．Resultsshow thatiftheratioisnotabove1000，adirectcalculation can beused，and otherwise contactlayersand transverselayershave tobe placed between the twozonesin a dueorderforsmoother transition．In the lattercase，a coefficienttensormustbe used forthe pemr eability ofcontactzone， or incorrectcalculationre