›› 2014, Vol. 35 ›› Issue (7): 2064-2070.

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

基于局部单元劈裂的FEM/DEM自适应分析方法

严成增,孙冠华,郑 宏,葛修润   

  1. 中国科学院武汉岩土力学研究所 岩土力学与工程国家重点实验室,武汉 430071
  • 收稿日期:2013-06-08 出版日期:2014-07-10 发布日期:2014-07-11
  • 作者简介:严成增,男,1986年生,博士生,研究兴趣包括:离散元有限元耦合、离散元、颗粒物质物理与力学、多物理场耦合、DDA、大规模高性能并行计算以及软件工程。
  • 基金资助:

    国家自然科学青年基金项目(No. 11202223);国家重点基础研究发展计划(973)项目(No. 2011CB013505)。

Adaptive FEM/DEM analysis method based on local splitting elements

YAN Cheng-zeng,SUN Guan-hua,ZHENG Hong,GE Xiu-run   

  1. State Key Laboratory of Geomechanics and Geotechnical Engineering, Institute of Rock and Soil Mechanics, Chinese Academy of Sciences, Wuhan 430071, China
  • Received:2013-06-08 Online:2014-07-10 Published:2014-07-11

摘要: 为了模拟岩体中裂纹的萌生、扩展,Munjiza提出了有限元法/离散元法(FEM/DEM)耦合分析方法。因为裂纹是沿单元边界进行扩展的,亦即裂纹扩展具有网格依赖性,为获得较好的裂纹扩展形态,需要划分密集的初始网格。为解决上述难题,基于FEM/DEM耦合分析方法,提出了基于局部单元动态劈裂的FEM/DEM自适应分析方法,以克服裂纹扩展形态对网格的依赖性。该方法在最初建模时无需划分很密的初始网格,随着荷载的施加,对裂纹尖端附近的局部单元进行动态劈裂,为裂纹的后续扩展提供了更多可能的扩展方向,使得裂纹扩展不必沿着初始网格的单元边界扩展,即可以沿着单元内部进行扩展,裂纹扩展形态更为平滑,与实际情况更为接近。同时相对原FEM/DEM耦合分析方法一开始就划分很密的网格而言,新方法可以划分较为稀疏的初始网格,计算成本降低。最后,通过巴西劈裂算例与原FEM/DEM耦合分析方法对比,分析表明,新方法在一定程度上克服了裂纹扩展形态对初始网格的依赖性。

关键词: 有限元法/离散元法(FEM/DEM), 单元劈裂法, 自适应技术, 裂纹扩展

Abstract: In order to simulate initiation and propagation of rock crack, Munjiza proposed FEM/DEM coupling analysis method. Because the cracks extend along the cell boundary, i.e. the propagation of crack has mesh-dependency problems. For the purpose of a better crack propagation shape, dense initial mesh is needed. To solve the above problems, the FEM/DEM adaptive analysis methods of local unit dynamic splitting is proposed to overcome mesh-dependency problems of crack morphology based on FEM/DEM coupling analysis method. However, very fine elements will not be necessary in this procedure; and dynamic splitting can take place in the local elements containing crack tips with increasing loading, which provides more possible expansion direction for subsequent expansion of the crack. Cracks are not confined along edges of elements. Cracks can grow in more natural way, which can be extended inside the element and more smoother. While original FEM/DEM coupling method divides very dense mesh, the proposed method can divide initial mesh relatively sparsely; so the cost of calculation is low. Finally, the numerical example of Brazilian disc indicates that the proposed method overcomes the dependence of crack morphology on the initial grid.

Key words: finite element method/discrete element method (FEM/DEM), method of element splitting, procedure of adaptive elements, crack growth

中图分类号: 

  • O 319.56
[1] 张国凯, 李海波, 王明洋, 李晓锋, . 基于声学测试和摄像技术的单裂隙岩石 裂纹扩展特征研究[J]. 岩土力学, 2019, 40(S1): 63-72.
[2] 张传庆, 刘振江, 张春生, 周辉, 高阳, 侯靖, . 隐晶质玄武岩破裂演化及破坏特征试验研究[J]. 岩土力学, 2019, 40(7): 2487-2496.
[3] 王海军, 郁舒阳, 任然, 汤雷, 李欣昀, 贾宇, . 基于3D-ILC含内裂纹孔口脆性固体断裂特性试验[J]. 岩土力学, 2019, 40(6): 2200-2212.
[4] 杨石扣, 张继勋, 任旭华, . 基于改进数值流形法的接触裂纹问题研究[J]. 岩土力学, 2019, 40(5): 2016-2021.
[5] 张 波, 李 垚, 杨学英, 朱飘扬, 朱春帝, 刘子豪, 刘文杰, 罗志恒, . 一种用于水压致裂试验的水压供给装置研制及应用[J]. 岩土力学, 2019, 40(5): 2022-2028.
[6] 丛 怡, 丛 宇, 张黎明, 贾乐鑫, 王在泉, . 大理岩加、卸荷破坏过程的三维颗粒流模拟[J]. 岩土力学, 2019, 40(3): 1179-1186.
[7] 李铮, 郭德平, 周小平, 王允腾, . 模拟岩石中裂纹扩展连接的近场动力学方法[J]. 岩土力学, 2019, 40(12): 4711-4721.
[8] 张珂, 王海军, 任然, 汤雷, 郁舒阳, 刘鑫娜, 顾浩, . 基于3D-ILC球体45º三维双内裂纹复合断裂研究[J]. 岩土力学, 2019, 40(12): 4731-4739.
[9] 马鹏飞, 李树忱, 周慧颖, 赵世森. 岩石材料裂纹扩展的改进近场动力学方法模拟[J]. 岩土力学, 2019, 40(10): 4111-4119.
[10] 王桂林, 梁再勇, 张 亮, 孙 帆, . Z型裂隙对砂岩强度和破裂行为影响机制研究[J]. 岩土力学, 2018, 39(S2): 389-397.
[11] 伍天华, 周 喻, 王 莉, 孙金海, 赵 欢, 孙 铮, . 单轴压缩条件下岩石孔-隙相互作用机制细观研究[J]. 岩土力学, 2018, 39(S2): 463-472.
[12] 杨石扣,张继勋,任旭华,. 基于数值流形法的三维裂纹扩展研究[J]. , 2018, 39(S1): 488-494.
[13] 徐文彬,曹培旺,程世康,. 不同偏置裂纹充填体断裂特性试验[J]. , 2018, 39(5): 1643-1652.
[14] 宋义敏, 邢同振, 吕祥锋, 赵泽鑫, 邓琳琳, . 不同加载速率Ⅰ型预制裂纹花岗岩断裂特征研究[J]. 岩土力学, 2018, 39(12): 4369-4375.
[15] 郭 洋,李 清,徐文龙,钱 路,田 策. 条形药包爆破预制贯通裂纹动态断裂过程研究[J]. , 2018, 39(10): 3882-3890.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 刘小文,常立君,胡小荣. 非饱和红土基质吸力与含水率及密度关系试验研究[J]. , 2009, 30(11): 3302 -3306 .
[2] 黄建华,宋二祥. 大型锚碇基础围护工程冻结帷幕力学性态研究[J]. , 2009, 30(11): 3372 -3378 .
[3] 姚仰平,冯 兴,黄 祥,李春亮. UH模型在有限元分析中的应用[J]. , 2010, 31(1): 237 -245 .
[4] 王朝阳,许 强,倪万魁. 原状黄土CT试验中应力-应变关系的研究[J]. , 2010, 31(2): 387 -391 .
[5] 万少石,年廷凯,蒋景彩,栾茂田. 边坡稳定强度折减有限元分析中的若干问题讨论[J]. , 2010, 31(7): 2283 -2288 .
[6] 刘 嘉,王 栋. 正常固结黏土中平板锚基础的吸力和抗拉力[J]. , 2009, 30(3): 735 -740 .
[7] 徐维生,柴军瑞,陈兴周,孙旭曙. 岩体裂隙网络非线性非立方渗流研究与应用[J]. , 2009, 30(S1): 53 -57 .
[8] 赵尚毅,郑颖人,李安洪,邱文平,唐晓松,徐 俊. 多排埋入式抗滑桩在武隆县政府滑坡中的应用[J]. , 2009, 30(S1): 160 -164 .
[9] 刘振平,贺怀建,朱发华. 基于钻孔数据的三维可视化快速建模技术的研究[J]. , 2009, 30(S1): 260 -266 .
[10] 陈保国 ,孙金山 ,张 磊. 上埋式钢筋混凝土拱涵受力特性及地基处理研究[J]. , 2011, 32(5): 1500 -1506 .