裂隙岩体渗流场的无单元法模拟及渗流特性研究.pdf

《裂隙岩体渗流场的无单元法模拟及渗流特性研究.pdf》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、JournalofEngineeringGeology工程地质学报1004-9665／2012／20(4)一0570—06裂隙岩体渗流场的无单元法模拟及渗流特性研究刘昌军丁留谦宁保辉张顺福(中国水利水电科学研究院北京100038)摘要裂隙岩体因含有发育程度不同的裂隙、节理和断层等不连续面，致其渗透性具有各向异性、不连续性等特点，因此传统的有限元法对分布密集的裂隙岩体渗流场求解有一定的难度。本文提出了采用无单元Glaerkin法求解有自由面裂隙渗流问题，并推导了无单元法求解渗流场的基本方程和积分格式，给

2、出了应用罚函数法处理渗流边界条件和自由面处理方法。采用IDL语言编制了二维无单元法计算软件LIDAREFM。文中以北京怀柔桥梓镇某裂隙岩体边坡渗流场计算为例，研究了复杂裂隙共同作用下渗流场特性和自由面分布，讨论了不同开度、不同连通程度的裂隙对渗流场的影响。研究结果表明：无单元法可以较好地解决有密集裂隙的岩体渗流场的求解问题，实现了裂隙处结点任意加密以及积分网格的独立布置，避免了对有自由面和裂隙穿越的子域的重新处理，简化了渗流问题的求解过程。关键词裂隙岩体无单元法各向异性罚函数法自由面中图分类号：TU4

3、43文献标识码：AMESHLESSSIMULATIoNoFSEEPAGEFIELDANDSEEPAGECHARAC．TERISTICSFoRFRACTUREDRoCKLIUChangjunDINGLiuqianNINGBaohuiZHANGShunfu(ChinaInstituteofWaterResourcesandHydropowerResearch，Beijing100038)AbstractThediscontinuitiessuchascracks，jointsandfaultscanhav

4、ethevaryingdegreesofdevelopments．Thepermeabilityoffracturedrockcanshowthecharacteristicsofanisotropyanddiscontinuity．Thetraditionalfiniteele—mentmethodencountersdifficultyinthesolutionoftheseepagefieldforintensivefracturedrocks．Thispaperpro—．posesthesol

5、utionofseepageproblemswithfreesurface．ItisobtainedthroughmeshlessGalerkinmethod．thecorre—spondingbasicequationsandintegrationschemesarederived．Applicationsofpenaltyfunctionmethodtothepro2cessingprocedureofseepageboundaryconditionandfreesurfaceboundaryar

6、ealsoputforward．Two—dimensionalmeshlessmethodsoftwareLIDAREFMisdevelopedwithIDLlanguage．Researchonseepagefieldcharacteristicsanddistributionoffreesurfaceunderdifferentdiscontinuitycombinationsarecarriedout，takingthefracturedrockslopeatQiaozitown，Huairou

7、CountyinBeijingasanexample．Theseepagefieldfeaturesunderdifferentdegreesofopeningandconnectivityarediscussed．Resultsshowthatthemeshlessmethodhascertainadvantagesinsolvingseepagefieldofintensivefractureproblem．Thelayoutindependenceofgriddingnodesonintegra

8、tionpointsisa—chieved．There—processingoffreesurfaceandsub-domainwithcrackscrossingisavoided．Theprocessofsolvingtheseepageproblemissimplified．KeywordsFracturedrock，Meshlessmethod，Anisotropy，Penaltyfunctionmethod，Freesurface收稿13期：2