›› 2011, Vol. 32 ›› Issue (7): 2211-2218.

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

Hoek-Brown准则的主应力回映算法及其二次开发

陈培帅1,陈卫忠1,贾善坡1, 2,杨建平1   

  1. 1. 中国科学院武汉岩土力学研究所 岩土力学与工程国家重点实验室,武汉 430071;2. 长江大学 城市建设学院,湖北 荆州 434023
  • 收稿日期:2010-02-01 出版日期:2011-07-10 发布日期:2011-06-30
  • 作者简介:陈培帅,男,1985年生,硕士研究生,主要从事隧道及地下工程施工及围岩稳定方面的工作
  • 基金资助:

    国家高技术研究发展专项经费(863计划)(No. 2007AA11Z108);中国科学院“西部行动计划高新技术”项目(No. KGCX2-YW-506)

Stress return mapping algorithm of Hoek-Brown criterion in principal stress space and its redevelopment

CHEN Pei-shuai1,CHEN Wei-zhong1,JIA Shan-po1, 2,YANG Jian-ping1   

  1. 1. State Key Laboratory of Geomechanics and Geotechnical Engineering, Institute of Rock and Soil Mechanics, Chinese Academy of Sciences, Wuhan 430071, China; 2. School of Urban Construction, Yangtze University, Jingzhou, Hubei 434023, China
  • Received:2010-02-01 Online:2011-07-10 Published:2011-06-30

摘要: 在弹塑性模型数值实现的过程中,需要进行应力更新的回映算法。针对三维应力空间回映算法在奇异点收敛性方面的不足,提出主应力空间的回映算法,讨论了算法实现过程的应力空间转化问题,分析了应力更新过程中确定回映区域的方法,建立了相应的一致性刚度矩阵。基于大型有限元软件ABAQUS提供的用户材料子程序UMAT接口,利用 Fortran 编程语言,实现了Hoek-Brown准则主应力空间的隐式积分算法,利用开发的模型,进行了岩石常规三轴压缩试验的数值模拟,通过与ABAQUS内部的Mohr-Coulomb准则计算结果的对比,验证了模型和程序的可行性和准确性

关键词: 霍克-布朗本构模型, ABAQUS, 用户材料子程序, 主应力空间回映算法

Abstract: In the course of implementing elastoplastic model in the FEM, it is necessary to carry out the process of stress return mapping algorithm. It may be not converging at singular points in three-dimensional space. The transforming method in principal stress space is put forward to overcoming this shortcoming. Also, the process of determining the region of return mapping is analyzed. At last, the consistent stiffness matrix is deduced. Based on the implicit algorithm in principal stress space, the UMAT subroutine of the Hoek-Brown model is developed in ABAQUS with the Fortran programming language. Through the example of triaxial compression test, it is easy to verify the feasibility and accuracy of this method by comparing the result with Mohr-Coulomb criterion in ABAQUS.

Key words: Hoek-Brown criterion, ABAQUS, UMAT, return mapping algorithm in principal stress space

中图分类号: 

  • TU 45
[1] 刘庆彬, 潘懋, 刘洁, 郭艳军, 张小双, 姚健鹏, 李芳玉, . 基于ParaView的Abaqus有限元输出结果的 可视化与虚拟现实[J]. 岩土力学, 2019, 40(12): 4916-4924.
[2] 周佳锦,龚晓南,严天龙,张日红, . 软土地区填砂竹节桩抗压承载性能研究[J]. , 2018, 39(9): 3425-3432.
[3] 崔 旋,董威信,周汉民,孙书伟,. 基于广义塑性理论的尾砂本构模型在ABAQUS中的二次开发[J]. , 2018, 39(2): 745-752.
[4] 郭 洋,李 清,徐文龙,钱 路,田 策. 条形药包爆破预制贯通裂纹动态断裂过程研究[J]. , 2018, 39(10): 3882-3890.
[5] 王仁超,曹婷婷,刘严如. 亚塑性模型不同积分算法的数值实现[J]. , 2017, 38(5): 1510-1516.
[6] 龚迪光,曲占庆,李建雄,曲冠政,曹彦超,郭天魁. 基于ABAQUS平台的水力裂缝扩展有限元模拟研究[J]. , 2016, 37(5): 1512-1520.
[7] 李 明 ,郭培军 ,梁 力 ,李 鑫,. 含有硬包裹体分布的非均质岩石水力压裂特性研究[J]. , 2016, 37(11): 3130-3136.
[8] 刘 润 ,刘文彬 ,洪兆徽 ,王 乐,. 海底管线整体屈曲过程中土体水平向阻力模型研究[J]. , 2015, 36(9): 2433-2441.
[9] 魏丽敏,冯胜洋,何 群,杨 奇. 改进Koppejan蠕变模型及其应用[J]. , 2014, 35(6): 1762-1767.
[10] 周佳锦 ,王奎华 ,龚晓南 ,张日红 ,严天龙 ,许远荣,. 静钻根植竹节桩承载力及荷载传递机制研究[J]. , 2014, 35(5): 1367-1376.
[11] 戚玉亮 ,大塚久哲,. ABAQUS动力无限元人工边界研究[J]. , 2014, 35(10): 3007-3012.
[12] 付强,丁选明 ,刘汉龙 ,孔纲强,. 列车激振荷载下PCC桩复合地基动力分析[J]. , 2013, 34(S2): 413-420.
[13] 杨 淼 ,张忠苗 ,刘念武 ,施茂飞 ,张日红 . 新型螺旋成孔根植注浆竹节管桩抗压性状数值模拟研究[J]. , 2013, 34(7): 2119-2126.
[14] 陈培帅 ,陈卫忠 ,庄 严 . 基于断裂力学的岩爆破坏形迹两级预测方法研究[J]. , 2013, 34(2): 575-584.
[15] 刘春原,李光宏,李 兵. 预应力管桩振动台试验的数值分析[J]. , 2012, 33(S1): 265-269.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 吴昌瑜,张 伟,李思慎,朱国胜. 减压井机械淤堵机制与防治方法试验研究[J]. , 2009, 30(10): 3181 -3187 .
[2] 陈红江,李夕兵,刘爱华. 矿井突水水源判别的多组逐步Bayes判别方法研究[J]. , 2009, 30(12): 3655 -3659 .
[3] 和法国,谌文武,韩文峰,张景科. 高分子材料SH固沙性能与微结构相关性研究[J]. , 2009, 30(12): 3803 -3807 .
[4] 雷永生. 西安地铁二号线下穿城墙及钟楼保护措施研究[J]. , 2010, 31(1): 223 -228 .
[5] 肖 忠,王元战,及春宁,黄泰坤,单 旭. 波浪作用下加固软基上大圆筒结构稳定性分析[J]. , 2010, 31(8): 2648 -2654 .
[6] 柴 波,殷坤龙,陈丽霞,李远耀. 岩体结构控制下的斜坡变形特征[J]. , 2009, 30(2): 521 -525 .
[7] 赵洪波,茹忠亮,张士科. SVM在地下工程可靠性分析中的应用[J]. , 2009, 30(2): 526 -530 .
[8] 徐 扬,高 谦,李 欣,李俊华,贾云喜. 土石混合体渗透性现场试坑试验研究[J]. , 2009, 30(3): 855 -858 .
[9] 章定文,刘松玉,顾沉颖. 各向异性初始应力状态下圆柱孔扩张理论弹塑性分析[J]. , 2009, 30(6): 1631 -1634 .
[10] 邓华锋,张国栋,王乐华,邓成进,郭 靖,鲁 涛. 导流隧洞开挖施工的爆破振动监测与分析[J]. , 2011, 32(3): 855 -860 .