›› 2009, Vol. 30 ›› Issue (S1): 215-219.

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

基于ABAQUS的颗粒材料亚塑性模型的数值实现及探讨

楚锡华,徐远杰,张明龙,余 村   

  1. 武汉大学 土木建筑工程学院,武汉 430072
  • 收稿日期:2009-05-27 出版日期:2009-08-10 发布日期:2011-03-16
  • 作者简介:楚锡华,男,1977年生,博士,主要从事计算固体力学与颗粒材料力学行为研究。
  • 基金资助:

    国家自然科学基金(No. 10802060)。

Numerical implementations and discussion on hypoplastic model for granular materials based on ABAQUS

CHU Xi-hua, XU Yuan-jie, ZHANG Ming-long, YU Cun   

  1. School of Civil and Archifectural Engineering, Wuhan University, Wuhan 430072, China
  • Received:2009-05-27 Online:2009-08-10 Published:2011-03-16

摘要:

颗粒材料的亚塑性模型以Jaumann应力率张量及变形率张量描述本构关系,基于Cauchy应力的Jaumann速率及变形率给出了有限元的切线刚度矩阵,由此指出在严格意义上基于ABAQUS二次开发的亚塑性模型数值实现必需借助UEL接口。为了简化程序开发,文中建议了与有限元相结合的两种近似切线模量矩阵,即基于矩阵广义逆一致性切向模量矩阵与近似切向模量矩阵,从而形成了通过UMAT接口实现亚塑性模型的数值方案。由此降低了程序开发的难度,同时也可借助ABAQUS的后处理功能,提高了工作效率。数值算例表明了所开发程序的正确性以及所建议方案的可行性。

关键词: 颗粒材料, 亚塑性模型, 切向刚度阵, 切向模量矩阵, ABAQUS

Abstract:

The constitutive relation of hypoplastic model for granular materials is represented by Jaumann stress rate and deformation rate, the tangential stiffness matrix for finite element simulation of hypoplasticity is given based on Gudehus-Bauer model with Jaumann rate of Cauchy stress and deformation rate, and it is illustrated that the numerical implementations of hypoplastic model based on ABAQUS must resort to UEL interface proper. For simplifying the of programming development, two methods for tangential modulus matrix are suggested combing with finite element method, viz. tangential modulus matrix based on general invert of matrix and approximate tangential modulus matrix, so numerical methods for implementations of hypoplastic model by UMAT interface are formed. By the approach suggested; the complexity of programming is reduced; and the post processing of ABAQUS can be utilized, so the work efficiency is improved. Numerical examples show the validity of methods suggested and fensibility of programming code.

Key words: granular materials, hypoplastic model, tangential stiffness matrix, tangential module matrix, ABAQUS

中图分类号: 

  • O347
[1] 付龙龙, 周顺华, 田志尧, 田哲侃, . 双轴压缩条件下颗粒材料中力链的演化[J]. 岩土力学, 2019, 40(6): 2427-2434.
[2] 刘庆彬, 潘懋, 刘洁, 郭艳军, 张小双, 姚健鹏, 李芳玉, . 基于ParaView的Abaqus有限元输出结果的 可视化与虚拟现实[J]. 岩土力学, 2019, 40(12): 4916-4924.
[3] 周佳锦,龚晓南,严天龙,张日红, . 软土地区填砂竹节桩抗压承载性能研究[J]. , 2018, 39(9): 3425-3432.
[4] 肖勇杰,陈福全,董译之. 基于Gudehus-Bauer模型的砂土中灌注桩护壁套管高频振动贯入速率[J]. , 2018, 39(8): 3011-3019.
[5] 崔 旋,董威信,周汉民,孙书伟,. 基于广义塑性理论的尾砂本构模型在ABAQUS中的二次开发[J]. , 2018, 39(2): 745-752.
[6] 郭 洋,李 清,徐文龙,钱 路,田 策. 条形药包爆破预制贯通裂纹动态断裂过程研究[J]. , 2018, 39(10): 3882-3890.
[7] 王仁超,曹婷婷,刘严如. 亚塑性模型不同积分算法的数值实现[J]. , 2017, 38(5): 1510-1516.
[8] 郭兴文,赵 骞,顾水涛,蔡 新, . 基于黏弹性接触的颗粒材料蠕变特性研究[J]. , 2016, 37(S2): 105-112.
[9] 傅 华,赵大海 ,韩华强,凌 华,. 不同级配粗颗粒材料动力特性试验研究[J]. , 2016, 37(8): 2279-2284.
[10] 易 颖,周 伟,马 刚,杨利福,常晓林, . 基于精确缩尺的颗粒材料流变研究[J]. , 2016, 37(6): 1799-1808.
[11] 龚迪光,曲占庆,李建雄,曲冠政,曹彦超,郭天魁. 基于ABAQUS平台的水力裂缝扩展有限元模拟研究[J]. , 2016, 37(5): 1512-1520.
[12] 李 明 ,郭培军 ,梁 力 ,李 鑫,. 含有硬包裹体分布的非均质岩石水力压裂特性研究[J]. , 2016, 37(11): 3130-3136.
[13] 刘 润 ,刘文彬 ,洪兆徽 ,王 乐,. 海底管线整体屈曲过程中土体水平向阻力模型研究[J]. , 2015, 36(9): 2433-2441.
[14] 蒋明镜 ,金树楼 ,张 宁 , . 不同胶结尺寸的粒间胶结强度统一表达式[J]. , 2015, 36(9): 2451-2457.
[15] 陈 龙,楚锡华,徐远杰. 基于CLoE亚塑性模型的密砂循环荷载下力学行为模拟[J]. , 2015, 36(6): 1598-1605.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 王 刚,李术才,王明斌. 渗透压力作用下加锚裂隙岩体围岩稳定性研究[J]. , 2009, 30(9): 2843 -2849 .
[2] 刘玉成,曹树刚,刘延保. 可描述地表沉陷动态过程的时间函数模型探讨[J]. , 2010, 31(3): 925 -931 .
[3] 刘恩龙. 岩土破损力学:结构块破损机制与二元介质模型[J]. , 2010, 31(S1): 13 -22 .
[4] 介玉新,杨光华. 基于广义位势理论的弹塑性模型的修正方法[J]. , 2010, 31(S2): 38 -42 .
[5] 杨建民,郑 刚. 基坑降水中渗流破坏归类及抗突涌验算公式评价[J]. , 2009, 30(1): 261 -264 .
[6] 周 华,王国进,傅少君,邹丽春,陈胜宏. 小湾拱坝坝基开挖卸荷松弛效应的有限元分析[J]. , 2009, 30(4): 1175 -1180 .
[7] 叶 飞,朱合华,何 川. 盾构隧道壁后注浆扩散模式及对管片的压力分析[J]. , 2009, 30(5): 1307 -1312 .
[8] 陈 林,张永兴,冉可新. 考虑剪应力作用的挡土墙主动土压力计算[J]. , 2009, 30(S2): 219 -223 .
[9] 罗 强 ,王忠涛 ,栾茂田 ,杨蕴明 ,陈培震. 非共轴本构模型在地基承载力数值计算中若干影响因素的探讨[J]. , 2011, 32(S1): 732 -0737 .
[10] 王云岗 ,章 光 ,胡 琦. 斜桩基础受力特性研究[J]. , 2011, 32(7): 2184 -2190 .