力学

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

展开
  • 1. 武汉大学 工程力学系,湖北 武汉 430072; 2. 石家庄铁道大学 工程力学系,河北 石家庄 050043
常江芳(1988-),女,博士,讲师,主要从事岩土材料强度与破坏的数值研究. E-mail:cjf881024@163. com

收稿日期: 2015-01-19

  修回日期: 2016-01-14

  网络出版日期: 2016-09-01

基金资助

国家自然科学基金资助项目(11172216);国家重点基础研究发展计划项目(2010CB731502);湖北省自然科学基金资助项目(2013CFB287);河北省杰出青年科学基金资助项目(E2015210040)

Finite Element Simulation of Strain Localization in Transversely Isotropic Geomaterials

Expand
  • 1.Engineering Mechanics Department,Wuhan University,Wuhan 430072,Hubei,China; 2.Mechanics Engineering Department,Shijiazhuang Tiedao University,Shijiazhuang 050043,Hebei,China
常江芳(1988-),女,博士,讲师,主要从事岩土材料强度与破坏的数值研究. E-mail:cjf881024@163. com

Received date: 2015-01-19

  Revised date: 2016-01-14

  Online published: 2016-09-01

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)进行了数值实现. 通过将积分点材料强度随材料主方向及各向异性程度的变化关系与理论结果进行比较,验证了程序开发的正确性. 数值算例重点分析了材料主方向和各向异性参数对结构极限承载力及破坏模式的影响,且与经典连续体的结果进行了比较. 结果表明,文中方法能够较好地模拟具有横观各向同性的岩土材料的应变局部化现象.

本文引用格式

常江芳 徐远杰 楚锡华 . 横观各向同性岩土材料应变局部化现象的有限元模拟[J]. 华南理工大学学报(自然科学版), 2016 , 44(10) : 70 -80 . DOI: 10.3969/j.issn.1000-565X.2016.10.011

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.

文章导航

/