长江科学院院报 ›› 2013, Vol. 30 ›› Issue (12): 48-53.DOI: 10.3969/j.issn.1001-5485.2013.12.009

• 岩土工程 • 上一篇    下一篇

弹塑性各向异性损伤模型的FLAC3D开发与数值验证

杨云浩1,陈鸿杰2a,2b,王伟2a,2b   

  1. 1.中国水电顾问集团成都勘测设计研究院,成都610072;
    2.河海大学a.岩土力学与堤坝工程教育部重点实验室;b.岩土工程研究所,南京210098
  • 收稿日期:2012-12-26 修回日期:2013-12-10 出版日期:2013-12-10 发布日期:2013-12-10
  • 作者简介:杨云浩(1977-),男,山西黎城人,工程师,博士,从事岩石力学与地下工程方面的研究,(电话)15928690461(电子信箱)haoyunyang2001@163.com。

Implementation and Numerical Verification ofElastoplastic-Anisotropic Damage Constitutive Model in FLAC3D

YANG Yun-hao1,CHEN Hong-jie2,3,WANG Wei2,3   

  1. 1.Hydro China Chengdu Engineering Corporation,Chengdu 610072,China;
    2.Key Laboratory of Ministry of Education for Geomechanics and Embankment Engineering,Hobai University,Nanjing 210098,China;
    3.Research Institute of Geotechnical Engineering,Hohai University ,Nanjing 210098,China
  • Received:2012-12-26 Revised:2013-12-10 Online:2013-12-10 Published:2013-12-10

摘要: 探讨适用于硬岩的各向异性损伤模型在FLAC3D中的实现方法。首先,基于不可逆热力学框架下的弹塑性损伤本构基本理论,结合弹塑性力学塑性势与非关联塑性流动法则,导出了弹塑性各向异性损伤本构模型的增量形式及损伤演化速率的计算公式,在此基础上,采用复合M-C准则并考虑黏聚力随损伤的逐步弱化,详细给出了该损伤模型在FLAC3D中的实现流程并探讨了相关细节处理方法,编写了相应的DLL文件以供调用。进而,通过圆柱试样的压缩数值试验验证其正确性与合理性,结果表明该损伤模型能够较好地描述岩石材料的非线性力学行为;以一个概化的地下洞室分层开挖实例初步检验了模型的应用效果,计算所得洞周损伤较大部位与实际工程中直接或间接观测到的围岩损伤情况较为一致,由此表明以拉应变度量的损伤及其演化方程能够合理地反映硬岩的损伤本质。

关键词: 不可逆热力学, 弹塑性, 各向异性损伤, FLAC3D, 本构模型

Abstract: The implementation of an anisotropic damage constitutive model in FLAC3D for hard rock is presented in this paper. An incremental form of elastoplastic-anisotropic damage constitutive equation is first deduced using the plastic mechanics theory and the elastoplastic-damage constitutive theory within the framework of thermodynamics of irreversible process. Meanwhile the concrete form of evolution rate of damage tensor is determined by the consistence condition of damage criterion. On the basis of these formulas,the detailed flowchart of numerical calculation by using the elastoplastic-anisotropic damage model and M-C yield criterion is established and some details the numerical calculation are also discussed. A DLL(Dynamic Link Library) file corresponding to the elastoplastic-anisotropic damage model is written to be called by the FLAC3D software. The uniaxial and triaxial compression numerical tests of cylinder sample are then conducted to test the validity of this user-defined model. The results show that this model is capable of describing the nonlinear behavior of hard rock. The validity of elastoplastic-anisotropic damage model in practice is further tested by applying it to a numerical simulation of staged excavation of a fictitious underground opening. Result show that the highly damaged zone around the opening determined by the simulation is in agreement with the actual observation in project area. The result also show that using the tensile strain tensor to measure damage and establish damage criterion grasps the nature of hard rock damage.

Key words: thermodynamics of irreversible process, elastoplastic, anisotropic damage, FLAC3D, constitutive modelabout

中图分类号: