資源描述:
《裂隙巖體滲流場的無單元法模擬及滲流特性研究.pdf》由會員上傳分享,免費在線閱讀,更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。
1、JournalofEngineeringGeology工程地質(zhì)學(xué)報1004-9665/2012/20(4)一0570—06裂隙巖體滲流場的無單元法模擬及滲流特性研究劉昌軍丁留謙寧保輝張順福(中國水利水電科學(xué)研究院北京100038)摘要裂隙巖體因含有發(fā)育程度不同的裂隙、節(jié)理和斷層等不連續(xù)面,致其滲透性具有各向異性、不連續(xù)性等特點,因此傳統(tǒng)的有限元法對分布密集的裂隙巖體滲流場求解有一定的難度。本文提出了采用無單元Glaerkin法求解有自由面裂隙滲流問題,并推導(dǎo)了無單元法求解滲流場的基本方程和積分格式,給
2、出了應(yīng)用罰函數(shù)法處理滲流邊界條件和自由面處理方法。采用IDL語言編制了二維無單元法計算軟件LIDAREFM。文中以北京懷柔橋梓鎮(zhèn)某裂隙巖體邊坡滲流場計算為例,研究了復(fù)雜裂隙共同作用下滲流場特性和自由面分布,討論了不同開度、不同連通程度的裂隙對滲流場的影響。研究結(jié)果表明:無單元法可以較好地解決有密集裂隙的巖體滲流場的求解問題,實現(xiàn)了裂隙處結(jié)點任意加密以及積分網(wǎng)格的獨立布置,避免了對有自由面和裂隙穿越的子域的重新處理,簡化了滲流問題的求解過程。關(guān)鍵詞裂隙巖體無單元法各向異性罰函數(shù)法自由面中圖分類號: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,F(xiàn)reesurface收稿13期:2