基于abaqus的金属橡胶阻尼器本构研究

《基于abaqus的金属橡胶阻尼器本构研究》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、基于Abaqus的金属橡胶阻尼器本构研究1212王振营，毛晨曦，赵亚哥白，王大磊（1.东北林业大学土木工程学院，哈尔滨150040；2.中国地震局工程力学研究所，哈尔滨150080）摘要:本文将金属橡胶试件单轴压缩试验应力-应变曲线，通过塑性数据转换后简化为四折线本构模型，并基于Abaqus提供的用户自定义子程序接口，利用FORTRAN语言编制了金属橡胶阻尼器本构子程序，将其嵌入到Abaqus材料库中。通过在单元模型和地震荷载作用下的12层剪力墙结构中调用子程序，分别验证了其正确性。关键词：金属橡胶阻尼器；单轴压缩；四折线简化本构模型；Abaqus；塑性数据

2、转换；VUMAT1引言金属橡胶是将金属丝卷成螺旋形，经过编织、加压成型，并经热处理而成的金属材料，不仅具有较大的自回复弹性变形，而且具有阻尼大、重量轻、柔韧性好、吸收冲击能、不惧高低温作[2]用、不易老化等特点。目前由金属橡胶制成的隔振器广泛应用于军事、航空航天、工程机械等领域，但金属橡胶材料在土木工程结构领域的应用尚属空白。目前广泛应用的有限元软件Abaqus具有很强的线性和非线性求解能力，不仅为用户提供了大量的单元库和求解模型，同时提供了大量的用户自定义子程序（usersubroutine）,使用户可以通过FORTRAN程序接口来定义Abaqus材[1]

3、[3]料库中没有的材料模型。本文将金属橡胶试件单轴压缩试验应力-应变曲线（以下简称为试验“σε−曲线”），通过塑性数据转换后简化为四折线本构模型，并基于Abaqus提供的用户自定义子程序接口，利用FORTRAN语言编制了金属橡胶阻尼器本构子程序（以下简称为“MR-VUMAT”），将其嵌入到Abaqus材料库中。通过在单元模型和地震荷载作用下的12层剪力墙结构中调用子程序，分别验证了其正确性，为金属橡胶阻尼器在土木工程中的应用提供设计和分析参考。2金属橡胶试件单轴压缩性能试验及拟合曲线2.1金属橡胶试件单轴压缩试验σε−曲线[2]本文选取文献(毛晨曦，2010

4、)中金属橡胶试件(成型密度为0.27、应变幅值为20%)在单轴压缩荷载作用下的试验σε−曲线，如图1所示：图1.金属橡胶试件单轴压缩试验σε−曲线.2.2金属橡胶试件单轴压缩试验σε−曲线的拟合[3]2.2.1Abaqus中塑性数据的转换在Abaqus中定义塑性数据时，必须使用真实应力和真实应变（输出的应力、应变也为真实应力和真实应变），而试验数据是用名义应力和名义应变给出的，这就必须将塑性材料的数据从名义应力和名义应变转化为真实应力和真实应变，具体转换过程如下：①用试件瞬时长度l表示应变，则每一瞬时应变增量为：dldε=（1）l从l开始变形所积累的应变（即

5、真实应变）：0lldllll+∆0εε=d==ll==l(1+ε)（2）∫∫nnnnomll00lll00式中:l为原始长度；l为当前长度；∆l为试件变形量；0∆lε为名义应变，即ε=nomnoml0②通过考虑塑性变形的不可压缩性(即体积不变原理)，真实应力和名义应力的关系为：lA=lA（3）00式中：A为原始面积；A为当前面积0由公式（3）可得当前面积和原始面积的关系表达式：l0AA=（4）0l将当前面积A带入真实应力的定义式，得到：FFllll+∆0σ===σσ()()(==σε1)+（5）nomnomnomnomAAlll0000F式中：σ为名义应力，

6、即σ=nomnomA02.2.2金属橡胶试件单轴压缩真实σε−曲线的拟合首先利用公式（2）和公式（5）将金属橡胶试件单轴压缩名义σε−曲线（即试验σε−曲线）转化为真实σε−曲线，如图2所示；为了便于对金属橡胶阻尼器本构的开发，本文通过最小二乘法和滞回耗能相等两个准则将真实σε−曲线拟合为四折线简化本构模型（即拟合σε−曲线，如图3所示）。图2.转换前后的σε−曲线.图3.四折线简化本构模型.2.2.3金属橡胶阻尼器的真实σε−曲线在实际设计阻尼器时，为了使其具有较好的耗能能力，因此设计的金属橡胶阻尼器能够在同一轴线上的两个方向同时具有压缩的力学性能。根据金

7、属橡胶阻尼器以上设计原理，由对称性可得金属橡胶阻尼器拉伸段（即反向压缩段）的真实σε−曲线（如图4所示）。（a）(b)图4.金属橡胶阻尼器的真实σε−曲线及拟合σε−曲线.3金属橡胶阻尼器VUMAT子程序开发[1]3.1VUMAT简介VUMAT是适用于ABAQUS/Explicit显示求解器的用户材料子程序。用户可使用FORTRAN语言编写VUMAT，定义所需要的本构关系，然后通过接口实现与ABAQUS的数据交互，实现新材料的添加与使用。3.2金属橡胶阻尼器本构子程序（以下简称“MR-VUMAT”）中材料常数如图4（b）所示的金属橡胶阻尼器的拟合σε−曲线，

8、本文在编写和调用MR-VUMAT子程序时用到的20个