结构动力弹塑性与倒塌分析(Ⅰ)——滞回曲线改进、ABAQUS子程序开发与验证.pdf

第37卷 第 1期 地 震 研 究 V0L37．No．1 2014年 1月 JOURNAL OFSEISMOLOGICAL RESEARCH Jan．，2014 结构动力弹塑性与倒塌分析 (I)水 滞回曲线改进、ABAQUS子程序开发与验证 柳国环 ，练继建 ，国 巍 (1．天津大学 建筑工程学院，天津300072；2．天津大学 水利工程仿真与安全国家重点实验室，天津300072； 3．中南大学 土木工程学院，湖南 长沙410075) 摘要 ：详述并明确了钢材静态损伤滞回模型中3种类型加卸载路径 ，重点改进了混凝土滞回规则中再加载路径拐 点前的二次曲线与拐点后的直线规则 ，定义了再加载至骨架曲线的刚度退化指数型表达式。进而，开发了直接能 够与ABAQUS主程序保持数据相互调用的UMAT隐式和VUMAT显式算法子程序，并给出子程序与ABAQUS主程 序链接的框架示意图。最后，通过5种具有试验结果的构件对所开发程序加以检验。结果表明：(1)改进的混 凝土再加载滞 回规则可较好描述混凝土加载曲线 ，克服了二次型规则的再加载曲线与卸载曲线在拐点前相交的可 能性，体现了拐点后再加载至骨架曲线的刚度退化；(2)开发实现的隐式算法UMAT与显式算法双精度VUMAT 子程序不仅可行，而且具有足够的计算精度。 关键词：滞回模型；弹塑性；ABAQUS；UMAT；VUMAT 中图分类号：TV318 文献标识码：A 文章编号：1000—0666(2014)01—0123—09 满足相应的容差，而显式算法无需迭代，但一般 0 前言 需要将步长设置为比较小 (与结构体系中单元尺 寸有 关) 以得 到可 靠可信 的计算 结果。在 当今，随着经济的快速发展和社会的不断进 ABAQUS中，可以编制子程序并与主程序对接来实 步，大跨和超高层等超限结构 日益增多。根据当 现结构的钢筋混凝土动力弹塑性分析，相应隐式 前的规范要求，超限结构的地震反应分析需要作 和显式算法的接 口子程序头文件分别为 UMAT和 为一项专门的工作来进行，如结构在大震作用下 VUMAT，采用 UMAT和 VUMAT中的已有变量参 的动力弹塑性表现行为 (叶献国等，2003；宣钢 数直接与程序中需要额外定义变量共同实现相应 等，2003；杜修力等，2007)，然后对计算结果是 本构子程序的代码。为此，本文首先根据已有资 否符合相关规范和 (或)规程的相关具体要求进 料详细回顾了钢材和混凝土材料的单轴本构模型， 一 步审查。一般情况下，对于这种特殊结构的计 并改进了混凝土卸载后再加载曲线的规则，旨在 算通常需要用到非线性计算能力相对较高级的大 使滞回曲线避免理论上存在加载曲线与卸载曲线 型分析软件，例如：ABAQUS、MSC．MARC和 Ls． 在拐点前相交的可能性，并使得再加载曲线形状 Dyna。柳国环等 (2014a)针对 ABAQUS大型程 更接近于试验曲线形状 (先凹后凸的曲线状)。然 序，做了滞回规则改进、相关程序开发以及试验 后，给出了ABAQUS中UMAT与VUMAT之间的区 结果对比，并通过带有大跨一输 电实际工程计算 别与联系，并按照UMAT与VUMAT的规则开发了 进一步检验了相关程序的可行性、计算精度与现 直接能够与 ABAQUS主程序保持数据相互调用的 实可靠性 。 显式和隐式算法子程序AB