›› 2013, Vol. 34 ›› Issue (2): 526-534.

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

接触处理Munjiza方法在数值流形法中的应用

徐栋栋,孙冠华,郑 宏   

  1. 中国科学院武汉岩土力学研究所 岩土力学与工程国家重点实验室,武汉 430071
  • 收稿日期:2012-07-13 出版日期:2013-02-11 发布日期:2013-03-01
  • 作者简介:徐栋栋,男,1986年生,博士研究生,主要从事计算岩石力学方面的研究工作。
  • 基金资助:

    国家杰出青年基金资助项目(No. 50925933);国家自然科学基金-面上项目(No. 11172313)

Application of Munjiza’s algorithm to contact treatment in numerical manifold method

XU Dong-dong, SUN Guan-hua, ZHENG Hong   

  1. State Key Laboratory of Geomechanics and Geotechnical Engineering, Institute of Rock and Soil Mechanics, Chinese Academy of Sciences, Wuhan 430071, China
  • Received:2012-07-13 Online:2013-02-11 Published:2013-03-01

摘要: 数值流形法(NMM)在接触处理过程包括:接触方式判断,包括点-点和点-线接触;开-闭迭代,确定块体系统的约束状态;以及接触传递,将有效的接触对传递到下一个时步。过程稍显繁琐和耗时,且与真实的物理接触状态有异。由Munjiza所提出的NBS(no binary search)接触检测算法,将单元映射到规则格子中,以链表结构将其有效地连接在一起,只在单元所在格子以及周围格子内部进行接触判断,接触检测效率大为提高,计算量仅随单元数线性增长,内存需求也很低。在计算接触力时以所定义的势为媒介,用重叠面积来衡量接触力的大小,属于分布式接触力,更接近于实际,避免了原NMM接触处理过程的繁琐。因此,将其作为一种平行的接触处理方法引入到NMM中,并以算例验证,证实了算法的可行性。

关键词: 势接触力, 接触判断, NMM, FEM-DEM

Abstract: Contact treatment of the numerical manifold method (NMM) mainly includes contact detection, including point-point and point-line contact, open-close iteration for determining the correct constrained state, and contact transfer by which effective contact couples can be transferred to the next time step. The whole process is slightly complicated and inefficient. Besides, the contact modes adopted in the NMM are quite different from the real physical contact state. The no binary search (NBS) contact detection algorithm proposed by Munjiza maps all elements into regular cells and connects the elements by the list structure effectively. The contact detections are limited within the cell and the four neighboring cells. With low requirement of random access memory, contact detection operation is just proportional to the number of the elements, and the efficiency is accordingly increased significantly. With the defined ‘potential’, contact pressure, belonging to the distribution category, is calculated by the contact overlapped area, which is closer to the real contact state, and avoids the cumbersome process involved in the NMM. As a result, it is introduced to the NMM as an alternative contact treatment method. The effectiveness is verified by the examples.

Key words: potential contact force, contact detection, numerical manifold method (NMM), FEM-DEM

中图分类号: 

  • O 241
[1] 柯锦福, 王水林, 郑宏, 杨永涛, . 基于修正对称和反对称分解的 三维数值流形元法应用推广[J]. 岩土力学, 2020, 41(2): 695-706.
[2] 杨石扣,张继勋,任旭华,. 基于数值流形法的三维裂纹扩展研究[J]. , 2018, 39(S1): 488-494.
[3] 严成增,郑 宏,葛修润. 基于统一标定的势接触力计算[J]. , 2015, 36(1): 249-256.
[4] 刘泉声,刘学伟. 多场耦合作用下岩体裂隙扩展演化关键问题研究[J]. , 2014, 299(2): 305-321.
[5] 周先齐 ,徐卫亚 ,钮新强 ,崔玉柱,. 离散单元法研究进展及应用综述[J]. , 2007, 28(S1): 408-417.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 王川婴,胡培良,孙卫春. 基于钻孔摄像技术的岩体完整性评价方法[J]. , 2010, 31(4): 1326 -1330 .
[2] 谈云志,孔令伟,郭爱国,万 智. 压实红黏土水分传输的毛细效应与数值模拟[J]. , 2010, 31(7): 2289 -2294 .
[3] 王生新,陆勇翔,尹亚雄,郭定一. 碎石土湿陷性试验研究[J]. , 2010, 31(8): 2373 -2377 .
[4] 王云岗,熊 凯,凌道盛. 基于平动加转动运动场的边坡稳定上限分析[J]. , 2010, 31(8): 2619 -2624 .
[5] 龙 照,赵明华,张恩祥,刘峻龙. 锚杆临界锚固长度简化计算方法[J]. , 2010, 31(9): 2991 -2994 .
[6] 史旦达,周 健,贾敏才,杨永香. 考虑蠕变性状的港区软土地基参数反演和长期沉降预测[J]. , 2009, 30(3): 746 -750 .
[7] 邓宗伟,冷伍明,李志勇,岳志平. 喷混凝土边坡温度场与应力场耦合的有限元时效分析[J]. , 2009, 30(4): 1153 -1158 .
[8] 李 健 ,谭忠盛 ,喻 渝 ,倪鲁肃. 下穿高速公路浅埋大跨度黄土隧道施工措施研究[J]. , 2011, 32(9): 2803 -2809 .
[9] 应宏伟 ,郑贝贝 ,谢新宇. 狭窄基坑平动模式刚性挡墙被动土压力分析[J]. , 2011, 32(12): 3755 -3762 .
[10] 刘振平,贺怀建,白世伟,张 华,邹 丹,李火兵. 软基中500 kV输电桩板塔基现场试验与数值模拟[J]. , 2012, 33(4): 1031 -1039 .