边坡在地震荷载作用下稳定性

a4214172

最近在做边坡在地震荷载作用下的稳定性分析，现在已有各时间点所对应的的加速度值，如何用table命令添加进去？ 附上我的部分命令流，老是显示有错误，还望各位大神帮忙看看！感激不尽！
; generate free-field (6 nodes on left, 16 nodes on right)
ffield gen left y -5,0 np 6
ffield gen right y -5,10 np 16
ffield change mat=1
;
; fix bottom
ffield base xvel=0
ffield base yvel=0
;
; top left FF
hist ffyd 10 1
* top right FF
hist ffyd 10 2
*
* cycle FF (not applied to grid) to get static stresses
solve for 1e-5
;
save slopeff270.sav
;
def properties
m_dens = 0.0025
b_mod = 16667.0
sh_mod = 10000.0
end
properties
;
def convert
c_p = sqrt((b_mod + (4.0 * sh_mod / 3.0)) / m_dens)
c_s = sqrt(sh_mod / m_dens)
norm_str = -2.0 * m_dens * c_p
shear_str = 2.0*(-2.0 * m_dens * c_s)
end
convert
;
boun xfree yfree
;
table 1             &
0.02        -0.0215239        &
0.04        -0.0249263        &
0.06        -0.0166032        &
0.08        -0.00162895        &
0.1        -0.0036255        &
0.12        -0.00364837        &
..................................
4.96        -0.0980703        &
4.98        -0.082755        &
5        -0.0698302       
; bound. cond. at base of free-field
ffield base syy norm_str hist table 1
ffield base sxy shear_str hist table 1
ffield base yvisc xvisc
;
; apply dynamic input
bou stress (0,0,norm_str) hist table 1 range y -5.1 -4.9
bou stress (0,shear_str,0) hist table 1 range y -5.1 -4.9
;
bound xvisc yvisc range y 9.0 10.1
bound xvisc yvisc range y -5.1 -4.9
bound mat 1
;
; apply free-field conditions
bou ff range x -0.1 0.1
bou ff range x 21.9 22.1
;
reset time disp hist
hist table 1
hist xvel 11,10 yvel 11,10
hist xvel 8.4,5.4 yvel 8.4,5.4
hist xvel 5.3,0.4 yvel 5.3,0.4
hist xdis 11,10 ydis 11,10
hist xdis  8.4,5.4 ydis 8.4,5.4
hist xdis 5.3,0.4 ydis 5.3,0.4
;
damp 0.125 10.0
;
cy time 5

LHW080957

兄弟，能不能分享一下您的代码，我最近也在做地震方面，不知道如何在udec里输入地震波，qq邮箱：1490979892@qq.com

zx1991

请问能不能分享一下您的代码，我是一个初学者，麻烦您了！1173538743@qq.com