华南理工大学学报(自然科学版) ›› 2016, Vol. 44 ›› Issue (10): 70-80.doi: 10.3969/j.issn.1000-565X.2016.10.011

• 力学 • 上一篇    下一篇

横观各向同性岩土材料应变局部化现象的有限元模拟

常江芳1,2 徐远杰1† 楚锡华1   

  1. 1. 武汉大学 工程力学系,湖北 武汉 430072; 2. 石家庄铁道大学 工程力学系,河北 石家庄 050043
  • 收稿日期:2015-01-19 修回日期:2016-01-14 出版日期:2016-10-25 发布日期:2016-09-01
  • 通信作者: 徐远杰(1956-),男,博士,教授,主要从事工程结构破坏数值仿真研究. E-mail:xyj9781@163.com
  • 作者简介:常江芳(1988-),女,博士,讲师,主要从事岩土材料强度与破坏的数值研究. E-mail:cjf881024@163. com
  • 基金资助:
    国家自然科学基金资助项目(11172216);国家重点基础研究发展计划项目(2010CB731502);湖北省自然科学基金资助项目(2013CFB287);河北省杰出青年科学基金资助项目(E2015210040)

Finite Element Simulation of Strain Localization in Transversely Isotropic Geomaterials

CHANG Jiang-fang1,2 XU Yuan-jie1 CHU Xi-hua1   

  1. 1.Engineering Mechanics Department,Wuhan University,Wuhan 430072,Hubei,China; 2.Mechanics Engineering Department,Shijiazhuang Tiedao University,Shijiazhuang 050043,Hebei,China
  • Received:2015-01-19 Revised:2016-01-14 Online:2016-10-25 Published:2016-09-01
  • Contact: 徐远杰(1956-),男,博士,教授,主要从事工程结构破坏数值仿真研究. E-mail:xyj9781@163.com
  • About author:常江芳(1988-),女,博士,讲师,主要从事岩土材料强度与破坏的数值研究. E-mail:cjf881024@163. com
  • Supported by:
    Supported by the National Natural Science Foundation of China(11172216),the National Program on Key Basic Research Project of China(2010CB731502)and the Natural Science Foundation of Hubei Province(2013CFB287)

摘要: 岩土材料的各向异性性质与其微结构紧密相关,文中考虑到内摩擦角是应力状态、组构张量及材料主方向的函数,发展了适用于横观各向同性岩土材料的修正的 Druck-
er-Prager 屈服准则. 基于 Cosserat 连续体理论,推导了该准则的一致性映射返回算法,形成了一致性切线模量矩阵,并利用有限元软件 Abaqus 的用户单元子程序(UEL)进行了数值实现. 通过将积分点材料强度随材料主方向及各向异性程度的变化关系与理论结果进行比较,验证了程序开发的正确性. 数值算例重点分析了材料主方向和各向异性参数对结构极限承载力及破坏模式的影响,且与经典连续体的结果进行了比较. 结果表明,文中方法能够较好地模拟具有横观各向同性的岩土材料的应变局部化现象.

关键词: 岩土材料, 岩土力学, 横观各向同性, Cosserat 连续体, 屈服准则, 应变局部化

Abstract:

The anisotropic properties of geomaterials are significantly related to their inherent microstructures.In this paper,a modified Drucker-Prager yield criterion for transversely isotropic materials is developed by evaluating the internal friction angle with the stress state,the microstructure tensor and the material principal direction.Then,based on the Cosserat continuum theory,a consistent return mapping algorithm for the modified criterion is formulated,and a consistent tangent modulus matrix is achieved.Moreover,the codes are implemented through the user defined element subroutine (UEL) in the finite element software Abaqus,and the correctness of the pro- gram is verified by comparing the theoretical results with the relationships of the material strength to the principal direction and anisotropic degree of the material in the integration points.Finally,the influences of the principal direction and anisotropic degree of the material on the bearing capacity and the failure mode of the structure are em- phatically analyzed by numerical examples,which are then compared with the results based on the classical contin- uum theory.It is found that the above-mentioned method is effective in simulating the strain localization of trans- versely isotropic geomaterials.

Key words: geomaterials, geomechanics, transverse isotropy, Cosserat continuum, yield criterion, strain locali- zation