裂隙岩体渗流场的无单元法模拟及渗流特性研究-工程地质学报.pdfVIP

ＪｏｕｒｎａｌｏｆＥｎｇｉｎｅｅｒｉｎｇＧｅｏｌｏｇｙ　　工程地质学报　　１００４－９６６５／２０１２／２０（４）０５７００６ 裂隙岩体渗流场的无单元法模拟及渗流特性研究 刘昌军　丁留谦　宁保辉　张顺福 （中国水利水电科学研究院　北京　１０００３８） 摘　要　裂隙岩体因含有发育程度不同的裂隙、节理和断层等不连续面，致其渗透性具有各向异性、不连续性等特点，因此传 统的有限元法对分布密集的裂隙岩体渗流场求解有一定的难度。本文提出了采用无单元Ｇｌａｅｒｋｉｎ法求解有自由面裂隙渗流 问题，并推导了无单元法求解渗流场的基本方程和积分格式，给出了应用罚函数法处理渗流边界条件和自由面处理方法。采 用ＩＤＬ语言编制了二维无单元法计算软件ＬＩＤＡＲＥＦＭ。文中以北京怀柔桥梓镇某裂隙岩体边坡渗流场计算为例，研究了复杂 裂隙共同作用下渗流场特性和自由面分布，讨论了不同开度、不同连通程度的裂隙对渗流场的影响。研究结果表明：无单元 法可以较好地解决有密集裂隙的岩体渗流场的求解问题，实现了裂隙处结点任意加密以及积分网格的独立布置，避免了对有 自由面和裂隙穿越的子域的重新处理，简化了渗流问题的求解过程。 关键词　裂隙岩体　无单元法　各向异性　罚函数法　自由面 中图分类号：ＴＵ４４３　　文献标识码：Ａ ＭＥＳＨＬＥＳＳＳＩＭＵＬＡＴＩＯＮＯＦＳＥＥＰＡＧＥＦＩＥＬＤＡＮＤＳＥＥＰＡＧＥＣＨＡＲＡＣ ＴＥＲＩＳＴＩＣＳＦＯＲＦＲＡＣＴＵＲＥＤＲＯＣＫ ＬＩＵＣｈａｎｇｊｕｎ　ＤＩＮＧＬｉｕｑｉａｎ　ＮＩＮＧＢａｏｈｕｉ　ＺＨＡＮＧＳｈｕｎｆｕ （ＣｈｉｎａＩｎｓｔｉｔｕｔｅｏｆＷａｔｅｒＲｅｓｏｕｒｃｅｓａｎｄＨｙｄｒｏｐｏｗｅｒＲｅｓｅａｒｃｈ，Ｂｅｉｊｉｎｇ　１０００３８） Ａｂｓｔｒａｃｔ　Ｔｈｅｄｉｓｃｏｎｔｉｎｕｉｔｉｅｓｓｕｃｈａｓｃｒａｃｋｓ，ｊｏｉｎｔｓａｎｄｆａｕｌｔｓｃａｎｈａｖｅｔｈｅｖａｒｙｉｎｇｄｅｇｒｅｅｓｏｆｄｅｖｅｌｏｐｍｅｎｔｓ．Ｔｈｅ ｐｅｒｍｅａｂｉｌｉｔｙｏｆｆｒａｃｔｕｒｅｄｒｏｃｋｃａｎｓｈｏｗｔｈｅｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆａｎｉｓｏｔｒｏｐｙａｎｄｄｉｓｃｏｎｔｉｎｕｉｔｙ．Ｔｈｅｔｒａｄｉｔｉｏｎａｌｆｉｎｉｔｅｅｌｅ ｍｅｎｔｍｅｔｈｏｄｅｎｃｏｕｎｔｅｒｓｄｉｆｆｉｃｕｌｔｙｉｎｔｈｅｓｏｌｕｔｉｏｎｏｆｔｈｅｓｅｅｐａｇｅｆｉｅｌｄｆｏｒｉｎｔｅｎｓｉｖｅｆｒａｃｔｕｒｅｄｒｏｃｋｓ．Ｔｈｉｓｐａｐｅｒｐｒｏ ｐｏｓｅｓｔｈｅｓｏｌｕｔｉｏｎｏｆｓｅｅｐａｇｅｐｒｏｂｌｅｍｓｗｉｔｈｆｒｅｅｓｕｒｆａｃｅ．ＩｔｉｓｏｂｔａｉｎｅｄｔｈｒｏｕｇｈｍｅｓｈｌｅｓｓＧａｌｅｒｋｉｎｍｅｔｈｏｄ．ｔｈｅｃｏｒｒｅ ｓｐｏｎｄｉｎｇｂａｓｉｃｅｑｕａｔｉｏｎｓａｎｄｉｎｔｅｇｒａｔｉｏｎｓｃｈｅｍｅｓａｒｅｄｅｒｉｖｅｄ．Ａｐｐｌｉｃａｔｉｏｎｓｏｆｐｅｎａｌｔｙｆｕｎｃｔｉｏｎｍｅｔｈｏｄｔｏｔｈｅｐｒｏ ｃｅｓｓｉｎｇｐｒｏｃｅｄｕｒｅｏｆｓｅｅｐａｇｅｂｏｕｎｄａｒｙｃｏｎｄｉｔｉｏｎａｎｄｆｒｅｅｓｕｒｆａｃｅｂｏｕｎｄａｒｙａｒｅａｌｓｏｐｕｔｆｏｒｗａｒｄ．Ｔｗｏｄｉｍｅｎｓｉｏｎａｌ ｍｅｓｈｌｅｓｓｍｅｔｈｏｄｓｏｆｔｗａｒｅＬＩＤＡＲＥＦＭｉｓｄｅｖｅｌｏｐｅｄｗｉｔｈＩＤＬｌａｎｇｕａｇｅ．Ｒｅｓｅａｒｃｈｏｎｓｅｅｐａｇｅｆｉｅｌｄｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ａｎｄｄｉｓｔｒｉｂｕｔｉｏｎｏｆｆｒｅｅｓｕｒｆａｃｅｕｎｄｅｒｄｉｆｆｅｒｅｎｔｄｉｓｃｏｎｔｉｎｕｉｔｙｃｏｍｂｉｎａｔｉｏｎｓａｒｅｃａｒｒｉｅｄｏｕｔ，ｔａｋｉｎｇｔｈｅｆｒａｃｔｕｒｅｄｒｏｃｋ ｓｌｏｐｅａｔＱｉａｏｚｉｔｏｗｎ，ＨｕａｉｒｏｕＣｏｕｎｔｙｉｎＢｅｉｊｉｎｇａｓａｎｅｘａｍｐｌｅ．Ｔｈｅｓｅｅｐａｇｅｆｉｅｌｄｆｅａｔｕｒｅｓｕｎｄｅｒｄｉｆｆｅｒｅｎｔｄｅｇｒｅｅｓ ｏｆｏｐｅｎｉｎｇａｎｄｃｏｎｎｅｃｔｉｖｉｔｙａｒｅｄｉｓｃｕｓｓｅｄ．Ｒｅｓｕｌｔｓｓｈｏｗｔｈａｔｔｈｅｍｅｓｈｌｅｓｓｍｅｔｈｏｄｈａｓｃｅｒｔａｉｎａｄｖａｎｔａｇｅｓｉｎｓｏｌｖｉｎｇ ｓｅｅｐａｇｅｆｉｅｌｄｏｆｉｎｔｅｎｓｉｖｅｆｒａｃｔｕｒｅｐｒｏｂｌｅｍ．Ｔｈｅｌａｙｏｕｔｉｎｄｅｐｅｎｄｅｎｃｅｏｆｇｒｉｄｄｉｎｇｎｏｄｅｓｏｎｉｎｔｅｇｒａｔｉｏｎｐｏｉｎｔｓｉｓａ ｃｈｉｅｖｅｄ．Ｔｈｅｒｅ－ｐｒｏｃｅｓｓｉｎｇｏｆｆｒｅｅｓｕｒｆａｃｅａｎｄｓｕｂｄｏｍａｉｎｗｉｔｈｃｒａｃｋｓｃｒｏｓｓｉｎｇｉｓａｖｏｉｄｅｄ．Ｔｈｅｐｒｏｃｅｓｓｏｆｓｏｌｖｉｎｇ ｔｈｅｓｅｅｐａｇｅｐｒｏｂｌｅｍｉｓｓｉｍｐｌｉｆｉｅｄ． Ｋｅｙ