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