《ABAQUS施加地震力》由会员分享，可在线阅读，更多相关《ABAQUS施加地震力（3页珍藏版）》请在金锄头文库上搜索。
1、初学 ABAQUS 没多久，由于课程需要，想用实体单元建了一个五层的框架结构，要对其施加地震波。但是我只学了石老师实例的前面部分，只知道些比较基本的操作之类的，于是上网求助于论坛。: l4 E/ g 8 U% x# 6 y/ h在找完很多帖子后，建模并计算，基本算是成功的。所以首先是要衷心感谢各位在论坛上指导了我的楼主及他们的帖子们。不过感觉网上各位大侠可能都没体会到民间疾苦，只提了个大概，我们自己得总结半天。所以在此详细点写下自己查到的方法，回报下论坛上的各位，给其他一些初学的人一些帮助，也达到交流学习的目的。5 A8 N0 |9 w* F b e但同时还有许多自己不知其所以然的，想要请教各
2、位大侠（红色字体是引用别人说的，蓝色字体是我的疑问，望大家讨论或帮忙解答） 。4 z) r5 t0 # T5 a+ a# G5 W! B: # X& G5 Z在网上查了些方法：( A$ , p& e1 d/ j( amodule 选 load，在 tools-amplitude-creat 默认的 continue 在 Edit Amplitude 里面输入时间和加速度，点 OK。点 creat boundary condition，出现对话框 creat boundary condition，选择 acceleration/angular acceleration,continue-选择要施
3、加的边界-done-出现对话框edit bondary condition 对话框，在 amplitude 里选择你所定义的时间和加速度。点 ok 就完工了。这是在 CAE 里输入地震波的方式，我用的方法是直接在 inp 文件里加地震波的。4 Z( A, r( ?4 i. S; a6 k8 |% + Q! H* S. ?$ x j首先在 CAE 里建好模型，定义两个分析步。! . F+ k+ J9 N, l2 l. P1 5 m7 8 Y+ |5 |9 w; A% 2 |. P+ C8 m. . F8 . w f第一个分析步是加自重，采用线性加载的方式。. C& h- e$ N 6 H(a)
4、加载方式：ABAQUS 在施加 Gravity 时，默认为 Instantaneous（瞬时加载） ，如果把结构自重以瞬间加载方式加到结构上，相当于对结构施加了一个脉冲荷载，会引起结构在竖向的振动，在不考虑结构阻尼的情况，这种振动会一直持续下去。如果是混凝土结构，这种竖向振动也会造成混凝土受拉损伤，所以这种加载方式不太合理。(b)新建加载方式：创建一个新的 Amplitude，Type=smooth tpye，0 时刻 Am=0，然后再选择一个 0.5s1s 时刻，Am=1，在这个区间内线性插值，实现幅值从 0 到 1。这种方式加载要优于上述瞬时加载，但是在起初的 0.5s（或者 1s，即 s
5、mooth tpye 中设置的终点时间）内计算结果是不准确的，所以要把这部分的计算结果剔除，剔除方法就是，创建 2 个step，第一个 step 主要分析自重作用，待自重稳定后开始第二个 step 地震时程反应分析。) P: 5 ; S$ n以上方法也是得自论坛，我感觉说的有道理，就用了，不知道具体对不对。7 W# r 5 V I$ K# B% Q, c, F$ f第二个分析步就是加地震波。 # o9 w2 a; o1 输入地震波有两种方法：1、在如下位置加入下面加黑的字体部分。格式如下：时间,地震波, 时间,地震波,时间, 地震波,时间, 地震波 每行 8 个数据（我下到的地震波文件是不带时
6、间的，自己用 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-07 n9 c: q( : 9 |8 o 4 t0.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-05%2、在如下位置加入下面加黑的字体部分。区别就在于上面那一种是直接把
7、地震波放在 inp文件里，而且每一行的最后面是没有逗号的（我也不知道这是不是必须的，帮助文件里没找着） 。而这种方法是将地震波另外放置在一个 txt 文件里，当然，这个文件得跟 inp 文件在同一个文件夹里。4 7 m+ P( $ v* M3 v9 x5 u 9 B( A; c# s l: g+ D, c! s8 h3 Y* O; s%0 X o0 P8 5 i*End Assembly*Amplitude, name=Amp-1, definition=SMOOTH STEP0.,0.,) % N/ t, ?# z0.5,0.5,1.,$ u8 J( _2 J8 m+ f) d2 n |1.
8、! w, k0 u% v# t: U& q*! _( d0 / a% % q+ k1 *Amplitude, name=Amp-2, input=DIZHENBO-E-W.txt* C9 J6 O/ f x) 1 a$ Q4 j; R M* MATERIALS: 9 C6 X8 I# 9 F# P1 y%4 i4 o! p/ q% q) y9 m3 其中 DIZHENBO -E-W.txt 是地震波文件。3 Y. 9 x* D/ w8 n* 9 R1 q7 Z# t7 9 % v( W! c) $ t& w# k0,0.0650966,0.02,0.0485266,0.04,-0.001183
9、57,0.06,-0.094686,% b$ e7 c9 L: M0.08,-0.102971,0.1,-0.0514855,0.12,-0.0153865,0.14,0.0224879,; w c6 a+ U c7 T& O3 K. _0.16,0.0432005,0.18,0.0538527,0.2,0.0420169,0.22,-0.00532609,$ V6 p8 9 U9 b. E0.24,-0.0426087,0.26,-0.0467512,0.28,-0.0272222,0.3,-0.00769324,1 Z) _ S- I. ) Z6 P a9 U3 * e1 o6 x9 L5
10、 l, X%对于边界条件，在初始步时将柱底约束住：* F* l* C/ k0 z( o! q%1 1 x3 q! z3 ?1 L* BOUNDARY CONDITIONS0 k9 y. T) S c; F5 p! B9 h$ f0 A2 h* Y/ x) T. S% ?: * Name: BC-1 Type: Acceleration/Angular acceleration# w) v! K9 G( & 6 B Q*Boundary, type=ACCELERATION5 G) D# p q q_PickedSet1399, 1, 1_PickedSet1399, 2, 2_PickedSe
11、t1399, 3, 39 t9 H5 _& ; E4 T3 G* %但关键在于这里边界条件的类型是 Displacement 还是如上所示的 Acceleration 我不能确定（第二步中的边界得是 Acceleration，第一步好像是不可以用 Displacement 的，不知道对不对） 。 G5 p9 U& + q6 t! L4 E8 i& b在第二个分析步中，将地震波以加速度的形式加在柱底上：0 s* T; b# e3 q) s3 F) K3 O, J% q x0 h5 t l2 J4 z- K q7 K K6 E. - p) d% f0 L9 v%* BOUNDARY CONDITI
12、ONS* Name: BC-1 Type: Acceleration/Angular acceleration*Boundary, type=ACCELERATION_PickedSet1399, 2, 2_PickedSet1399, 3, 3. 0 9 l* z9 I. P*Boundary, amplitude=Amp-2, type=ACCELERATION_PickedSet1399, 1, 1, 1) W; I3 f2 A: R( c%& x K/ D/ u0 F6 : J; U& s5 J6 9 _8 n + i c) e其中最后的那个“1”参数是指将地震波数据乘以此系数得到最终
13、地震波的值。这里我有点不明白的是蓝色字体部分，这里不加这些语句可以吗？因为在初始步中已经定义过了，但我不知道这时是不是得重新把 2、3 方向的给再次定义下。 （我知道有些东西我自己改下再运算一遍比较下结果也许可以知道，但一来不太明白的地方太多，初学也不知道改的对不对，二来我觉得这样可能知道结果，却不知道原因，所以发上来希望有大侠告知）! q/ a- j2 D. W u4 D( IInp 文件添加到这里也就可以放入 Command 里进行计算了。; B4 o+ l; - p8 * ? 结果如下图所示：, U- f. z6 $ K: P- Y( D6 q变形 20 倍图 位移曲线（在 ABAQUS 中将顶点跟底点位移曲线相减得到相对位移） w# J# O, K/ 7 c+ c$ z7 T我想介绍的就是这些，希望对于需要用到加载地震波的各位能有所帮助。同时，这里面除了我希望有疑问的地方能得到大侠们解答之外，如有其他错误或不足之处，请大家指出。谢谢！O(_)O 哈哈
