基于ABAQUS弹脆塑性本构模型二次开发及应用.pdfVIP

第29卷第1期 红水河 Vd．29．No．1 River Feb．2010 2010年2月 HongShui 基于ABAQUS的弹脆塑性本构 模型二次开发及其应用 刘 明1，章青1，刘仲秋2 (1．河海大学工程力学系，江苏南京210098；2．山东农业大学水利土木工程学院，山东泰安271018) 摘要：在复杂应力状态下，大理岩、花趣岩等脆性岩石呈弹脆塑性破坏形式。利用ABAQlUS提供的用户材井子 程序IMAT的二次开发平台，编制完成了基于D—P准则的弹脆塑性模型的UMAT接口程序。论文简要介绍了 ABAQUS有限元软件及其IMAT用户材料予程序，重点对基于D—P屈服准则的弹脆塑性本构模型的应力跌落 公式进行了推导，并讨论了有关的应力更新算法。最后应用该程序对某引水隧洞工程进行了数值分析，并与理想 弹塑性情况下的计算结果进行对比研究。结果表明：与理想弹塑性模型相比，弹脆塑性模型得到的洞周围岩变形 和塑性区范围较大，与实际情况更为接近。 关键词：弹脆塑性本构模型；D—P准则；UMAT子程序；围岩稳定 中图分类号：Tv223．3文献标识码：A 文章编号：1001—408x(20lO)Ol一0047一05 包括结构分析、热传导等多方面的问题，具有强大的 1引言 非线性分析功能。ABAQI您另一个重要的特点就 很多脆性岩石的单轴压缩试验的结果显示：当 在于它良好的开放性，提供了多种二次开发的接口， 岩石的应力强度达到一定峰值后其应力值会迅速降 可以方便用户根据自己的需要扩展ABAQUs的功 MaterialMe． 低，然后在一较低应力水平上做塑性流动【lJ，基本 能。用户材料子程序(User—defined 呈脆性破坏形式，尤其在复杂应力条件下这种脆性 chanical 破坏形式更为明显【2j。其应力应变关系可用图1所 用户定义自己的材料属性的Fortran程序接口。用 示的弹脆塑性模型[3]来表示。因此对该模型的研 户可以通过该程序中定义自己的本构模型，相当于 究对实际工程具有重要的意义，而目前的很多通用 定义了自己的材料模型，因此称做用户材料子程序。 商业软件的本构模型库中并没有这种模型，致使很 UMAT子程序不同与独立的有限元计算程序，它需 多采用通用商业软件计算的实际工程问题只能使用 要同ABAQUS进行数据传递，甚至共享一些变量， 理想弹塑性模型进行模拟，大大影响了计算结果的 因此UMAT子程序有其特定的书写格式[51。 准确度。如文献[4]中对三峡船闸开挖进行的理想 弹塑性与弹脆塑性对比分析发现：弹脆塑性的塑性 区比理想弹塑性的大60％以上。本文利用 ABAQI玲提供的用户材料子程序UMAT的二次开 发平台，编制基于D—P准则的三维弹脆塑性模型 的UMAT用户材料子程序，并应用该程序对某引 应变 水隧洞围岩进行数值分析。 图1弹脆塑性模型图 2 ABAQUS软件及其用户材料子程 序IMAT 3基于D-P准则的弹脆塑性本构模型 删Us是功能强大的有限元软件，可以模拟 弹脆塑性本构模型可以看做是由线弹性增长、 收稿日期：2009一09—14；修回日期：2009一09一” 基金项目：国家自然科学基金重点项目(No