›› 2014, Vol. 35 ›› Issue (S2): 662-666.

• 数值分析 • 上一篇    下一篇

六面体覆盖的高阶数值流形方法的探讨

刘登学,张友良,谭 飞,张礼仁   

  1. 中国科学院武汉岩土力学研究所 岩土力学与工程国家重点实验室,武汉 430071
  • 收稿日期:2014-03-27 出版日期:2014-10-31 发布日期:2014-11-12
  • 作者简介:刘登学,男,1991年生,硕士研究生,主要从事数值计算方面的研究。
  • 基金资助:

    国家自然科学基金面上项目(No. 11272330); 国家重点基础研究计划(973计划)项目(No. 2014CB047100)。

Discussion on high-order numerical manifold method based on hexahedron element

LIU Deng-xue, ZHANG You-liang, TAN Fei, ZHANG Li-ren   

  1. State Key Laboratory of Geomechanics and Geotechnical Engineering, Institute of Rock and Soil Mechanics, Chinese Academy of Sciences, Wuhan 430071, China
  • Received:2014-03-27 Online:2014-10-31 Published:2014-11-12

摘要: 采用高阶近似位移覆盖函数,基于六面体数学覆盖网格建立了三维数值流形方法分析格式,给出了相应的子矩阵。利用MATLAB编制了与之对应的计算程序,对简单的地下洞室模型进行了计算,并将计算结果与其他数值分析方法结果进行了比较,证明了分析格式及相应程序的正确性和有效性。结果表明, 当数值流形方法的覆盖函数推广到高阶情况时,其求解精度会有相应的提高。最后,对该方法在隧道及地下工程的应用前景作了展望。

关键词: 数值流形方法, 高阶, 隧道及地下工程

Abstract: By introducing two covers, namely, the mathematical cover and physical cover, the numerical manifold method can solve continuous and discontinuous problems in geotechnical engineering in unified way. At present, the application of NMM is mainly limited to two-dimensional analysis which is not realistic for many underground engineering. In this paper, three-dimensional numerical manifold method based on hexahedron element is proposed; and two kinds of high-order displacement functions are used as the cover function. The stiffness matrix, loading matrix, displacement resistance matrix are given according to the cover function; and the methods of handling rank deficiency induced by high degree polynomials as local approximation is discussed. Corresponding program is made by Matlab. One simple example is calculated by the program and the numerical results, being compared with the result of other methods, proves its correctness and effectiveness. It shows that the solution accuracy of numerical manifold method will be improved, when high-order displacement function is used. An example of underground chamber is also presented. Finally, the prospects of the method in tunnel and underground engineering are given.

Key words: numerical manifold method, high-order, tunnel and underground engineering

中图分类号: 

  • O 302
[1] 孙锐, 杨峰, 阳军生, 赵乙丁, 郑响凑, 罗静静, 姚捷, . 基于二阶锥规划与高阶单元的 自适应上限有限元研究[J]. 岩土力学, 2020, 41(2): 687-694.
[2] 陈栋, 王恩元, 李楠, . 不同煤岩介质模型的波场特征研究[J]. 岩土力学, 2019, 40(S1): 449-458.
[3] 刘登学, 张友良, 丁秀丽, 黄书岭, 裴启涛, . 数值流形法中基于适合分析T样条的 局部网格加密算法[J]. 岩土力学, 2019, 40(4): 1584-1595.
[4] 魏晓明, 郭利杰, 李长洪, 张立新, 罗文冲, 刘 仁, . 高阶段胶结充填体强度空间变化规律研究[J]. 岩土力学, 2018, 39(S2): 45-52.
[5] 韩智铭,乔春生,朱 举. 含两组交叉贯通节理岩体的强度及破坏特征分析[J]. , 2018, 39(7): 2451-2460.
[6] 邢浩洁,李鸿晶,杨笑梅,. 基于切比雪夫谱元模型的成层场地地震反应分析[J]. , 2017, 38(2): 593-600.
[7] 杨石扣,任旭华,张继勋. 基于布尔交运算的三维流形单元生成研究[J]. , 2016, 37(9): 2706-2711.
[8] 陈灯红 ,杜成斌,. 结构-地基动力相互作用的时域模型[J]. , 2014, 35(4): 1164-1172.
[9] 刘泉声,刘学伟. 多场耦合作用下岩体裂隙扩展演化关键问题研究[J]. , 2014, 299(2): 305-321.
[10] 刘 建,陈 佺. 岩体黏弹性蠕变计算的高阶数值流形法研究[J]. , 2012, 33(7): 2174-2180.
[11] 刘 建,陈 佺. 一种模拟岩石蠕变的数值流形方法[J]. , 2012, 33(4): 1203-1209.
[12] 姜清辉,邓书申,周创兵. 有自由面渗流分析的三维数值流形方法[J]. , 2011, 32(3): 879-884.
[13] 焦 健,乔春生,徐干成. 开挖模拟在数值流形方法中的实现[J]. , 2010, 31(9): 2951-2957.
[14] 周小义,邓安福. 六面体有限覆盖的三维数值流形方法的非线性分析[J]. , 2010, 31(7): 2276-2282.
[15] 刘红岩,王贵和. 节理岩体冲击破坏的数值流形方法模拟[J]. , 2009, 30(11): 3523-3527.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 张文杰,陈云敏. 垃圾填埋场抽水试验及降水方案设计[J]. , 2010, 31(1): 211 -215 .
[2] 宫伟力,安里千,赵海燕,毛灵涛. 基于图像描述的煤岩裂隙CT图像多尺度特征[J]. , 2010, 31(2): 371 -376 .
[3] 万 智,董 辉,刘宝琛. 基于正交设计下SVM滑坡变形时序回归预测的超参数选择[J]. , 2010, 31(2): 503 -508 .
[4] 孙曦源,栾茂田,唐小微. 饱和软黏土地基中桶形基础水平承载力研究[J]. , 2010, 31(2): 667 -672 .
[5] 王明年,郭 军,罗禄森,喻 渝,杨建民,谭忠盛. 高速铁路大断面黄土隧道深浅埋分界深度研究[J]. , 2010, 31(4): 1157 -1162 .
[6] 胡勇刚,罗 强,张 良,黄 晶,陈亚美. 基于离心模型试验的水泥土搅拌法加固斜坡软弱土地基变形特性分析[J]. , 2010, 31(7): 2207 -2213 .
[7] 谭峰屹,姜志全,李仲秋,颜惠和. 附加质量法在昆明新机场填料压实密度检测中的应用研究[J]. , 2010, 31(7): 2214 -2218 .
[8] 柴 波,殷坤龙,肖拥军. 巴东新城区库岸斜坡软弱带特征[J]. , 2010, 31(8): 2501 -2506 .
[9] 杨召亮,孙冠华,郑 宏. 基于潘氏极大值原理的边坡稳定性的整体分析法[J]. , 2011, 32(2): 559 -563 .
[10] 王光进,杨春和,张 超,马洪岭,孔祥云,侯克鹏. 超高排土场的粒径分级及其边坡稳定性分析研究[J]. , 2011, 32(3): 905 -913 .