急求：用flac2d 地震波对边坡的影响的做法

conf dyn gw
; shaking table test for liquefaction
g 1 5
m finn
gen 0 0 0 5 50 5 50 0
fix x y j=1
fix x
set grav 10, flow=off
prop dens 2000 shear 2e8 bulk 3e8
prop fric 35 poros 0.5
water dens 1000 bulk 2e9 tens 1e10
ini pp 5e4 var 0 -5e4
ini syy -1.25e5 var 0 1.25e5
ini sxx -1e5 var 0 1e5 szz -1e5 var 0 1e5
prop ff_latency=50
;
; parameters for Martin formula
prop ff_switch = 0
prop ff_c1=0.8 ff_c2=0.79
prop ff_c3=0.45 ff_c4=0.73
;
; parameters for Byrne formula
; prop ff_switch = 1
; def _setCoeff_Byrne
; ff_c1_ = 8.7*exp(-1.25*ln(n1_60_))
; ff_c2_ = 0.4/ff_c1_
; ff_c3_ = 0.0000
; end
; set n1_60_ = 7
; _setCoeff_Byrne
; prop ff_c1=ff_c1_ ff_c2=ff_c2_
; prop ff_c3=ff_c3_
;
set ncwrite=50
def sine_wave
while_stepping
vv = ampl * sin(2.0 * pi * freq * dytime)
loop j (1,jzones)
vvv = vv * float(jgp - j) / float(jzones)
loop i (1,igp)
xvel(i,j) = vvv
end_loop
end_loop
end
def eff_stress
eff_stress = (sxx(1,2)+syy(1,2)+szz(1,2))/3.0 + pp(1,2)
settlement = (ydisp(1,jgp)+ydisp(2,jgp))/2.0
end
set dy_damp=rayl 0.05 20.0
his dytime
his pp i 1 j 2
his eff_stress
his settlement
his nstep 20
set ampl=0.005 freq=5.0
solve dyt=10.0
plot hold his 2 3 vs 1 skip 2
save martin.sav
; save byrne.sav

程序很好懂啊，主要是FISH语言的用法

这个问题问的好，发表程序的兄弟能把你的程序说明一下吗，有点看不懂，谢谢！

看看相关的东西能找到你要的不？在帮助文件里有

看看相关的东西能找到你要的不？在帮助文件里有