ABAQUS二次开发 三维多层土黏弹性边界+地震动输入Python脚本

       地震动输入前，需要设置黏弹性人工边界，以实现基地的弹性恢复能力。在ABAQUS中，黏弹性人工边界需要在边界节点设置并联弹簧阻尼器，手动实现工作量大，可以借助Python程序依次读取节点并设置弹簧阻尼器。该程序考虑了不同土层地震波传播速度及波在土层分界面反射、透射问题，实现三维六层图以内的三向地震动入射。模型信息如下：

1、模型尺寸

为便于计算，将模型定为长条型，模型长宽各200m，高2400m，土层分为四层，层高均为600m，网格大小40×40×40m。

2、材料参数

自下向上土层编号1、2、3、4。土层弹性模量分别为10e9, 2e9, 10e9, 2e9。泊松比均为0.25，密度2000。

3、波幅放大系数：

土层1→2：透射系数:1.38，反射系数0.62。土层2→1：透射系数:0.62，反射系数-0.38。

4、地震波位移时程曲线：

5、验证（见附件）