ABAQUS 粘弹性边界及地震荷载施加的简单实现
粘弹性边界因为能够考虑地基辐射阻尼而使得结构抗震的计算结果更趋于合理，所以在需要考虑结构地基
相互作用的结构抗震计算时，是较为常用的地基边界处理和地震荷载施加方法。而 ABAQUS 软件是经常用
来进行结构响应分析的有限元软件。下面介绍一种在 ABAQUS 中实现粘弹性边界及地震荷载施加的方法。
粘弹性边界是通过在有限元模型的地基边界节点上施加弹簧阻尼器实现的，在 ABAQUS 中的实现有以下几
种方法：第一种，通过 ABAQUS 自有的弹簧单元 spring 单元和阻尼单元 dashpot 实现，具体的单元参数可
以参考文献[1]，这种较为精确；第二种是通过 ABAQUS 的 UEL 子程序实现，可以看下文献[2]；还有一种是
等效单元替代的方法，就是在地基周围加一层单元，然后设置近似的材料参数，参考文献[3]，这一种精度
较差，但实现起来较为简单。我采用的是第一种方法，但操作起来较为繁琐，具体程序及过程后面介绍。
采用粘弹性边界，其配套的地震荷载输入方法就是在已知输入地震位移和速度的情况下，计算各个时刻地
基边界各个结点上应当施加的集中力荷载，然后施加荷载，一步一步的进行计算。地震荷载的施加在 ABAQUS
中也有两种不同的思路，文献[2]中的方法是通过 ABAQUS 的 DLOAD 和 UTRACLOAD 两个子程序实现。DLOAD
子程序用于施加边界面的法向荷载，UTRACLOAD 用于施加边界面的切向荷载。而文献[1]中则是将边界上每
一个节点每个时刻的力都计算出来，然后导入到 ABAQUS 中作为幅值数据，对每个对应节点施加。
我最初的想法是两篇文章的思路各取一半，用文献[1]的方法实现粘弹性边界，用文献[2]的方法施加地震荷
载。然而尝试了很久，发现这样做的效果并不是太好，可能我编的程序哪儿还是有问题吧。最后放弃了，
统一采用文献[1]的方法实现，具体实现采用 MATLAB 语言生成 ABAQUS 的 input 文件，然后将生成的 input
文件在模型文件的指定位置插入，用 ABAQUS 运行即可。
MATLAB 程序与运行前的准备工作
首先需要准备一些必要的数据文件（上图中红色框内的文件），其余黄色框内为模型文件，蓝色框内的文件
为程序运行后的生成文件，
boundary1~5.rpt 是从 ABAQUS 反力文件中提取的反力文件，其值代表某一节点的控制面积，可以通过在地
基边界（5 个面）施加值为 1 的压力荷载，即可提取得到这些反力文件；
coord_point.rpt 为 5 个边界面上节点的坐标文件，提取方法可以在很容易的百度到；
DIS.txt 是三个方向地震波的位移文件；
VEL.txt 是三个方向地震波的速度文件；
job-996.inp 是模型的 inp 文件；
Amplitude.inp 里面是计算过程中边界节点上随时间要施加的所有集中力荷载，文件较大；
load.inp 是将 Amplitude.inp 里面的幅值施加到对应节点的荷载命令；
springs&dashpot.inp 是模型地基边界施加的弹簧阻尼器文件；
三个 input 文件在模型 inp 文件中的插入位置：
记住 springs&dashpot.inp 在 Assembly 部分，所以搜索到关键字*End Assembly，把 springs&dashpot.inp 放在
*End Assembly 之前，Amplitude.inp 放在*End Assembly 之后，load.inp 放在 load 部分即可，如下所示
·································
·································
剩余10页未读，继续阅读
评论0