在ABAQUS中对框架结构施加地震波.docx

在ABAQUS 中对框架结构施加地震波（对初学者普及，同时向 大虾们求教）初学ABAQUS没多久，由于课程需要，想用实体单元建了一个五 层的框架结构，要对其施加地震波但是我只学了石老师《实例》 的前面部分，只知道些比较基本的操作之类的，于是上网求助于 论坛在找完很多帖子后，建模并计算，基本算是成功的所以首先是 要衷心感谢各位在论坛上指导了我的楼主及他们的帖子们不过感觉网上各位大侠可能都没体会到民间疾苦，只提了个大概， 我们自己得总结半天所以在此详细点写下自己查到的方法，回 报下论坛上的各位，给其他一些初学的人一些帮助，也达到交流 学习的目的但同时还有许多自己不知其所以然的，想要请教各位大侠（红色 字体是引用别人说的，蓝色字体是我的疑问，望大家讨论或帮忙 解答）在网上查了些方法：module 选 load，在 tools——amplitude——creat 默认的 continue 在 Edit Amplitude里面输入时间和加速度，点OK点creat boundary condition，出现对话框 creat boundary condition，选择 acceleration/angular acceleration,continue---选择要施加的边界 ---done----出现对话框 edit bondary condition 对话框，在 amplitude 里选择你所定义的时间和加速度。

点ok就完工了这是在CAE里输入地震波的方式，我用的方法是直接在inp文件 里加地震波的首先在CAE里建好模型，定义两个分析步第一个分析步是加自重，采用线性加载的方式a） 加载方式：ABAQUS在施加Gravity时，默认为Instantaneous （瞬时加载），如果把结构自重以瞬间加载方式加到结构上，相当于对结构施加了一个脉冲荷载，会引起结构在竖向的振动，在 不考虑结构阻尼的情况，这种振动会一直持续下去如果是混凝 土结构，这种竖向振动也会造成混凝土受拉损伤，所以这种加载 方式不太合理b） 新建加载方式：创建一个新的Amplitude，Type=smooth tpye， 0时刻Am=0，然后再选择一个0.5s~1s时亥U，Am=1，在这个区间 内线性插值，实现幅值从0到1这种方式加载要优于上述瞬时 加载，但是在起初的0.5s （或者1s，即smooth tpye中设置的终 点时间）内计算结果是不准确的，所以要把这部分的计算结果剔 除，剔除方法就是，创建2个step，第一个step主要分析自重作 用，待自重稳定后开始第二个step地震时程反应分析以上方法也是得自论坛，我感觉说的有道理，就用了，不知道具 体对不对。

1liudaomuxing助理工程师精华 积分 帖子 水位0161215第二个分析步就是加地震波输入地震波有两种方法：1、在如下位置加入下面加黑的字体部分格式如下：时间，地震 波，时间，地震波,时间，地震波，时间，地震波 每行8个数据（我下到的地震波文件是不带时间的，自己用C++处理了一下）End Assembly*Amplitude, name=Amp-10.005, -7.5e-08, 0.01, ・3.55e-07, 0.015, -7.03e-07, 0.02, -4.53e-070.025, 1.82e-06, 0.03, 7.01e-06, 0.035, 1.5e-05, 0.04, 2.49e-050.045, 3.54e・05, 0.05, 4.5e-05, 0.055, 5.2e-05, 0.06, 5.5e-052009-7-13 00:08技术分0%%%%%%%%%%%%%%%%%%%%%%2、在如下位置加入下面加黑的字体部分区别就在于上面那一 种是直接把地震波放在inp文件里，而且每一行的最后面是没有 逗号的（我也不知道这是不是必须的，帮助文件里没找着）而 这种方法是将地震波另外放置在一个txt文件里，当然，这个文 件得跟inp文件在同一个文件夹里。

End Assembly*Amplitude, name=Amp-1, definition=SMOOTH STEP0., 0., 0.5, 0.5, 1., 1.** *Amplitude, name=Amp-2, input=DIZHENBO--E-W.txt ♦*** MATERIALS%%%%%%%%%%%%%%%%%%%%%%其中DIZHENBO --E-W.txt是地震波文件0,0.0650966,0.02,0.0485266,0.04,-0.00118357,0.06,-0.094686,0.08,-0.102971,0.1,-0.0514855,0.12,-0.0153865,0.14,0.0224879,0.16,0.0432005,0.18,0.0538527,0.2,0.0420169,0.22,-0.00532609,0.24,-0.0426087,0.26,-0.0467512,0.28,-0.0272222,0.3,-0.00769324, %%%%%%%%%%%%%%%%%%%%%%2009-7-13#200:09liudaomuxing助理工程师精华0积分16帖子12水位15技术分0对于边界条件，在初始步时将柱底约束住：%%%%%%%%%%%%%%%%%%%%%%** BOUNDARY CONDITIONS**** Name: BC-1 Type: Acceleration/Angular acceleration*Boundary, type=ACCELERATION_PickedSet1399, 1, 1_PickedSet1399, 2, 2_PickedSet1399, 3, 3%%%%%%%%%%%%%%%%%%%%%%但关键在于这里边界条件的类型是Displacement还是如上所示 的Acceleration我不能确定（第二步中的边界得是Acceleration， 第一步好像是不可以用Displacement的，不知道对不对）。

在第二个分析步中，将地震波以加速度的形式加在柱底上：%%%%%%%%%%%%%%%%%%%%%%** BOUNDARY CONDITIONS**** Name: BC-1 Type: Acceleration/Angular acceleration*Boundary, type=ACCELERATION_PickedSet1399, 2, 2_PickedSet1399, 3, 3*Boundary, amplitude=Amp-2, type=ACCELERATION_PickedSet1399, 1, 1, 1%%%%%%%%%%%%%%%%%%%%%%其中最后的那个“ 1”参数是指将地震波数据乘以此系数得到最 终地震波的值这里我有点不明白的是蓝色字体部分，这里不加 这些语句可以吗？因为在初始步中已经定义过了，但我不知道这 时是不是得重新把2、3方向的给再次定义下我知道有些东 西我自己改下再运算一遍比较下结果也许可以知道，但一来不太 明白的地方太多，初学也不知道改的对不对，二来我觉得这样可 能知道结果，却不知道原因，所以发上来希望有大侠告知）Inp文件添加到这里也就可以放入Command里进行计算了。

结果如附件所示：20倍变形图位移曲线（在ABAQUS中将顶点跟底点位移曲线相减得到相对位 移）我想介绍的就是这些，希望对于需要用到加载地震波的各位能有 所帮助同时，这里面除了我希望有疑问的地方能得到大侠们解答之外， 如有其他错误或不足之处，请大家指出谢谢！ O（n_n）O哈哈~。