FLAC2D手册中给出了速度波与应力波的换算公式。仅从关系式来看，似乎两者之间只存在量值上的倍数关系。但我将2个例题运行了一遍，然后查看了模型底部输入的应力波和输入波之间的关系，却发现其中一个符合这种关系，但另一个却并不是这么一回事。下面粘出例题3-2和3-3的命令流。
不知道有人注意到这种区别没？或者是我哪里没理解到，还请大家指教！
例题3-2（符合公式）
config dyn
grid 1,50
model elas
prop dens 1000 shear 1e7 bulk 2e7
def wave
if dytime > 1.0/freq
wave = 0.0
else
wave = 0.5 * (1.0 - cos(2.0*pi*freq * dytime))
endif
end
set freq=4.0
fix y
apply xquiet j=1
apply sxy -2e5 hist wave j=1
hist xvel i=1 j=1
hist xvel i=1 j=26
hist xvel i=1 j=51
hist sxy i=1 j=1
hist dytime
solve dytime 1.8
从上到下分别为应力波 速度波
例题3-3 （不符合公式）
; -- Free-field test --
config dyn
def wave
wave = 0.5 * (1.0 - cos(2*pi*dytime/period))
end
set period 0.015
grid 16 10
mod elas
gen line 6 10 8 6
gen line 8 6 10 10
mod null reg 7 10
prop bulk 66667 shear 40000 den 0.0025
set grav 10
fix x i=1
fix x i=17
fix y j=1
;
set dyn off
hist unbal
hist ydis i 5 j 5
hist ydis i 5 j 11
solve
save ff0.sav
;
set dyn on
;
apply ff
;
apply xquiet j=1
apply yquiet j=1
apply sxy -1.0 hist wave j 1
pause
;
set dytime 0
hist reset
hist dytime
hist xvel i 5 j 11
hist xvel i 18 j 11
hist xvel i 19 j 11
hist sxy i 8 j 1
hist xacc i 8 j 1
hist xvel i 8 j 1
;
solve dytime 0.02
从上到下分别为应力波 速度波
