求助zvel的方向对结果的影响？

linzheng200

有一个简单的例子如下：
new
conf dyn
gen zone brick size 3 3 3     
model elas
prop bulk 3e8 shear 1e8
ini dens 2000
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 dyn dt=2e-4
fix z ran z -0.1 0.1
fix x ran x -0.1 0.1
fix x ran x 2.9 3.1
fix y ran y -0.1 0.1
fix y ran y 2.9 3.1
apply zvel 1 hist wave  yvel 0 zvel 0  range z 3 x 1 2 y 1 2
set grav 0 0 -10
hist unbal
hist gp vel 0 0 3
hist gp vel 0 3 3
solve age 0.25
其中把wave改写成wave=0.5*(1.0 - cos(2.0*pi*freq*dytime))好像对结果没有影响，为什么？求教高手！

linzheng200

速度分量的方向问题，还有单位怎么定义的呢？

cloud328

应当不会出现这样的问题。