关于波在土壤内部的传播

avazhou

set log on
set log nquiet.log
config dyn
grid 1 50
model elas
prop dens 1000 shear 1e7 bulk 2e7
def wave                                                                    ；定义wave
    if dytime > 1.0/freq
    wave = 0.0
    else
    wave = 0.5 * (1 - cos(2*pi*freq * dytime))
    endif
end
set freq=4
fix y
apply xquiet j=1
;apply sxy -2e5 hist wave j=1                                          ;把wave放在底部
apply sxy -2e5 hist wave j=11                                         ;把wave放在j=11的位置

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 sxy  i=1 j=26
hist sxy  i=1 j=51
hist dytime
solve dytime 1.8
plot hold hist 1 vs 7
plot hold hist 2 vs 7
plot hold hist 3 vs 7
plot hold hist 4 vs 7
plot hold hist 5 vs 7
plot hold hist 6 vs 7
set hisfile nquiet.his
hist write 1,2,3,4,5,6,7,7 nquiet.his
请教各位大侠，我现在做的是个很简单的波传递模型，模型高50m,颗粒波速是1m/s,剪切波速是100m/s.波从底部开始传递（j=1)时候我得到的结果是合理的。 但是当我把input wave放在j=11的时候，我发现在i=1 j=1;i=1 j=26;i=1 j=51三个点得到的速度值是一样的，形状是一样，但是没有时间的推迟。我假如把input wave放在j=11的地方，那么到达点i=1 j=1的时间应该是(10/100)=0.1s.但是程序里面得到的在i=1 j=1的速度还是从0开始，请教各位大侠，这个是什么原因？是因为程序有错吗？

还有当我把input wave放在j=26或者大于j=26的地方，得到的速度曲线已经方向改变，幅值为原来的两倍，而且在i=1 j=1的地方没有反射波，请教各位大侠，这个结果是合理的吗？如果不合理应该怎样改正？谢谢各位！感激不尽！