岩体p型自适应块体单元法研究.pdfVIP

第24 卷 第4 期 岩石力学与工程学报 Vol.24 No.4 2005 年2 月 Chinese Journal of Rock Mechanics and Engineering Feb. ，2005 岩体p 型自适应块体单元法研究 汪卫明，陈胜宏 (武汉大学 水资源与水电工程科学国家重点试验室，湖北 武汉 430072) 摘要：初步研究了岩体p 型自适应块体单元法的基本理论和数值实现方法，以解决块体单元法中计算精度的控制 问题。根据当前的块体单元法数值解，利用广义刚度矩阵的阶谱特性，估计各更高阶广义节点形函数的引入对提 高计算精度所起的作用；然后引入作用较大的高阶广义节点形函数，扩充广义刚度矩阵，并求解新的整体平衡方 程，得到更高精度的数值解。经过几次自适应升阶，即可得到满足精度要求的计算结果。通过算例分析，验证了 所提方法的准确性和高效性。 关键词：岩石力学；数值分析；块体单元法；p 型自适应；阶谱；岩体 中图分类号：TU 452 文献标识码：A 文章编号：1000 –6915(2005)04 –0559 –06 STUDY ON p-VERSION ADAPTIVE BLOCK ELEMENT METHOD FOR ROCK MASSES WANG Wei-ming ，CHEN Sheng-hong (State Key Laboratory of Water Resources and Hydropower Engineering Science ，Wuhan University，Wuhan 430072 ，China) Abstract ：A p-version adaptive block element method for rock masses is presented to control the computational accuracy of block element method. According to current numerical solution of block element method ，the contribution of higher-order shape function on each generalized node to the improvement of computational accuracy is evaluated by using the hierarchical characteristics of generalized stiffness matrix. Then ，those useful generalized nodes are selected ，and the generalized stiffness matrix is enlarged. A more accurate numerical solution can be obtained after the new equilibrium equations are solved. The adaptive processes are repeated until the results are accurate enough. The methodology proposed is proved to be valuable and efficient by a numerical example. Key words ：rock mechanics ；numerical analysis ；block element method ；p-version adaptation ；hierarchy ；rock mass