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

第30卷第12期 长 江科学 院 院 报 V01．30No．12 201 3年12月 Journalof RiverScientificResearchInstitute Yangtze Dec．20 1 3 DOI：10．3969／j．issn。1001-5485．2013．12．009 8’26 杨云浩1，陈鸿杰2虬拍，王伟2 (1．中国水电顾问集团成都勘测设计研究院，成都610072； 2．河海大学a．岩土力学与堤坝工程教育部重点实验室；b．岩土工程研究所，南京210098) 摘要i探讨适用于硬岩的各向异性损伤模型在FLAC3D中的实现方法。首先，基于不可逆热力学框架下的弹塑性损伤 本构基本理论，结合弹塑性力学塑性势与非关联塑性流动法则，导出了弹塑性各向异性损伤本构模型的增量形式及 损伤演化速率的计算公式，在此基础上，采用复合M-C准则并考虑黏聚力随损伤的逐步弱化，详细给出了该损伤模型 在FLAC”中的实现流程并探讨了相关细节处理方法，编写了相应的DLL文件以供调用。进而，通过圆柱试样的压缩 数值试验验证其正确性与合理性，结果表明该损伤模型能够较好地描述岩石材料的非线性力学行为；以一个概化的 地下洞室分层开挖实例初步检验了模型的应用效果，计算所得洞周损伤较大部位与实际工程中直接或间接观测到的 围岩损伤情况较为一致，由此表明以拉应变度量的损伤及其演化方程能够合理地反映硬岩的损伤本质。 关键词：不可逆热力学；弹塑性；各向异性损伤；FLAC”；本构模型 中图分类号：TU452 文献标志码：A 文章编号：1001—5485(2013)12—0048—06 胀这一典型特征，但该模型只能考虑各向同性损伤。 文献[9]通过大理岩三轴卸荷试验发现试样主破裂面 1 研究背景 周边常有较多的近垂直于卸荷方向的微张裂隙，由此 在岩体地下工程领域，施工期围岩的稳定性长 表明试样轴向与径向的损伤必然是不同的。因此，对 期以来都是采用常规弹塑性本构模型进行分析。对 于硬岩，有必要考虑损伤的各向异性。文献[10—11] 于硬岩地下洞室，开挖卸荷引起的原岩应力重分布 分别针对花岗岩和混凝土所建立的损伤模型考虑各 导致围岩内微裂纹大量萌生，微裂纹的逐步扩展使 向异性损伤，将二阶损伤张量的损伤主值表达为拉应 得围岩的刚度和强度属性发生劣化¨≈J，从而在宏 变的幂函数，但该模型仅考虑弹性损伤且其损伤演化 观的应力应变曲线上表现出非线性特征，而传统的 方程在理论上不够严谨。文献[12]所提出的弹塑性 弹塑性模型通常不能反映这些特征，相应的稳定性 各向异性损伤模型将二阶损伤张量与应变张量的正 分析成果的可靠性也就大打折扣。 锥直接关联并假定二者共轴。该模型建立在不可逆 近年来，不少学者引入损伤的概念来描述岩石受 热力学基础上，理论依据严谨，所采用的损伤准则及 载过程中变形和强度属性的劣化并建立了相应的损 损伤演化方程也能很好地反映硬岩损伤的物理机制， 伤本构模型。文献[3]研究了弹脆塑性损伤本构模 模型参数少且物理意义较明确，有良好的应用前景。 型，但未明确损伤的启动与演化准则。文献[4—6]所 借助大型商业软件提供的二次开发接口，编制 提出的弹塑性损伤本构模型将损伤变量取为标量，以 可供其调用的用户自定义本构模块，是将最新理论 应变能或等效塑性应变作为损伤演化参量。然而，文 研究成果快速应用于工程实践的有效途径。目前，