技术“慧” | Abaqus 对金属橡胶阻尼器本构的研究

本文将金属橡胶试件单轴压缩试验应力应变曲线，通过塑性数据转换后简化为四折线本构模型，并基于Abaqus 提供的用户自定义子程序接口，利用FORTRAN 语言编制了金属橡胶阻尼器本构子程序，将其嵌入到 Abaqus 材料库中。通过在单元模型和地震荷载作用下的 12层剪力墙结构中调用子程序，分别验证了其正确性。

扫码可下载白皮书

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

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

金属橡胶试件单轴压缩性能试验及拟合曲线

金属橡胶试件单轴压缩试验σ-s曲线

本文选取文献（毛晨曦，2010）中金属橡胶试件（成型密度为 0.27、应变幅值为 20%）在单轴压缩荷载作用下的试验o-s曲线，如图1所示：

图 1 a金属橡胶试件单轴压缩试验σ-s曲线

金属橡胶试件单轴压缩试验σ-s曲线的拟合

Abaqus中塑性数据的转换

在 Abaqus 中定义塑性数据时，必须使用真实应力和真实应变（输出的应力、应变也为真实应力和真实应变），而试验数据是用名义应力和名义应变给出的，这就必须将塑性材料的数据从名义应力和名义应变转化为真实应力和真实应变，具体转换过程如下：

1、用试件瞬时长度l表示应变，则每一瞬时应变增量为：

金属橡胶试件单轴压缩真实σ-s曲线的拟合

首先利用公式(2)和公式(5)将金属橡胶试件单轴压缩名义σ-s曲线(即试验σ-s曲线)转化为真实σ-s曲线，如图2所示；为了便于对金属橡胶阻尼器本构的开发，本文通过最小二乘法和滞回耗能相等两个准则将真实σ-s曲线拟合为四折线简化本构模型（即拟合σ-s曲线如图3 所示）。

图 2 转换前后的σ-ε曲线 图3 四折线简化本结构模型

金属橡胶阻尼器的真实σ-ε 曲线

在实际设计阻尼器时，为了使其具有较好的耗能能力，因此设计的金属橡胶阻尼器能够在同一轴线上的两个方向同时具有压缩的力学性能。根据金属橡胶阻尼器以上设计原理，由对称性可得金属橡胶阻尼器拉伸段(即反向压缩段)的真实o-s曲线（如图4所示）。

图 4.金属橡胶阻尼器的真实σ-s曲线及拟合σ-s曲线,

金属橡胶阻尼器 VUMAT 子程序开发

VUMAT简介

VUMAT 是适用于 Abaqus/Explicit 显示求解器的用户材料子程序。用户可使用 FORTRAN语言编写 VUMAT，定义所需要的本构关系，然后通过接口实现与 Abaqus 的数据交互，实现新材料的添加与使用。

金属橡胶阻尼器本构子程序(以下简称“MR-VUMAT”)中材料常数

如图 4(b)所示的金属橡胶阻尼器的拟合σ-s曲线，本文在编写和调用 MR-VUMAT 子程序时用到的 20 个材料参数如表1所示。

注:PROPS(1)~PROPS(10)的意义及数值与 PROPS(11)~PROPS(20)相对应(拉压对称)

MR-VUMAT 在单元中的验证

单元模型（如图5所示）及模型基本参数

图 5 单元模型

模型基本参数：

长度为 25mm；截面积为 625mm；边界条件：左端施加固定边界条件，右端施加沿x方向正方向位移边界条件，其幅值曲线如图6所示；分析步骤：动态显示分析步；单元个数：1个；单元类型：桁架单元T3D2(A2-node linear 3-D truss)

图6.加载的位移幅值曲线.

结果分析

由于状态变量 SDV9 与 SDV30 分别记录子程序中应变与应力的更新过程，而变量 LE11 与S11又分别为单元中积分点沿1方向的应变和应力(即为真实应变与应力)，两种方法得到的σ-s曲线具有很好的一致性(如图7所示)，验证了子程序的正确性。

图 7.两种方法得到的σ-s曲线对比.

MR-VUMAT 在 ABAOUS 中的应用

选取文献(董金芝，2010)中的剪力墙模型，将结构仅受拉压作用的部位放置的SMA棒替换成金属橡胶阻尼器，长度为 40mm；在基底的X方向施加 30s，峰值为 510gal的EL-Centro地震波(NSImperial Valley，1940年5月)，如图8所示。子程序中的材料常数的设置如表1所示。通过提取调用MR-VUMAT形状记忆合金阻尼器的真实σ-s曲线(即S11-LE11，如图9所示)，进一步验证了子程序的正确性。

图 8.加载的 EL-Centro 地震波

图 9. 1-12 层形状记忆合金阻尼器的真实σ-s曲线