FLAC3d 中波动方程wave怎么定义

fg888134

我定义了一段wave函数，怎么plot的wave没有周期的变化：下附命令流，麻烦高手给看看。
new
config dyn
set dyn off
gen zone brick p0 6.5 0 14.9 p1 17.5 0 14.9 p2 6.5 40 14.9 p3 6.5 0 15.2 size 22 40 3 ;道床板层
gen zone brick p0 5.5 0 10 p1 18.5 0 10 p2 5.5 40 10 p3 5.5 0 14.9 size 26 40 49 ;梯形正方形层
gen zone wedge p0 18.5 0 10 p1 24 0 10 p2 18.5 40 10 p3 18.5 0 14.9 p4 24 40 10 p5 18.5 40 14.9 size 11 40 49 ;右楔形
gen zone wedge p2 5.5 0 10 p4 0 0 10 p0 5.5 40 10 p5 5.5 0 14.9 p1 0 40 10 p3 5.5 40 14.9 size 11 40 49 ;左楔形
gen zone brick p0 -3 0 0 p1 27 0 0 p2 -3 40 0 p3 -3 0 10 size 60 40 20 ;第 7 层

group 1 range z 14.9 15.2
group 2 range z 14.6 14.9
group 3 range z 13.9 14.6
group 4 range z 11.6 13.9
group 5 range z 10.6 11.6
group 6 range z 10 10.6
group 7 range z 0 10

;CFG 桩建立  nseg结构单元的数量


;弹性求解初始应力场
model elas
pro dens 2800 bulk 3.79e8 shear 3.21e8 range group 1 any
pro dens 2600 bulk 1.94e8 shear 1.46e8 range group 2 any
pro dens 2400 bulk 8.70e7 shear 4.72e7 range group 3 any
pro dens 2300 bulk 4.00e7 shear 2.40e7 range group 4 any
pro dens 2100 bulk 3.13e7 shear 1.79e7 range group 5 any
pro dens 2000 bulk 2.22e7 shear 1.67e7 range group 6 any
pro dens 1900 bulk 4.67e6 shear 1.92e6 range group 7 any

;施加边界约束条件
fix x range x -3.1 -2.9
fix x range x 26.9 27.1
fix y range y -0.1 0.1
fix y range y 39.9 40.1
fix z range z 15.19 15.21


ini dens 100
set gravity 0,0,-10
hist unbal
plot cont zdisp rev on
set mech ratio 3.0e-3
solve

;位移和速度归零
ini xdisp 0
ini ydisp 0
ini zdisp 0
ini xvel 0
ini yvel 0
ini zvel 0

save fosSEL1.sav

;定义材料参数
model elas
pro dens 2800 bulk 3.79e8 shear 3.21e8 range group 1 any
pro dens 2600 bulk 1.94e8 shear 1.46e8 range group 2 any
pro dens 2400 bulk 8.70e7 shear 4.72e7 range group 3 any
pro dens 2300 bulk 4.00e7 shear 2.40e7 range group 4 any
pro dens 2100 bulk 3.13e7 shear 1.79e7 range group 5 any
pro dens 2000 bulk 2.22e7 shear 1.67e7 range group 6 any
pro dens 1900 bulk 4.67e6 shear 1.92e6 range group 7 any

set dyn on
def iniwave
  per=0.01
end
iniwave
def wave
  wave=0.5*(1.0-cos(2*pi*dytime/per))
end
free x y z range z -.1 .1
apply nquiet squiet dquiet ran z -.1 .1
apply dstress 83.2e3 hist wave ran z 15.19 15.21
apply ff
group ff_corner
group ff_side range x -3 27
group ff_side range y 0 40
group main_grid ran x -3 27 y 0 40
set dyn time = 0
set dyn damp rayleigh 2 22.8 stiffness
hist reset
hist unbal
hist wave
history dytime
history gp zdis 12 20 15.2
history gp xvel 12 20 15.2
history gp zdis 12 20 14.9
history gp zdis 12 20 13.9
history gp zdis 12 20 11.6
history gp zdis 12 20 10
history gp zdis 12 20 9
plot hist 4 5 6 7 8 9 vs 3
plot add cont zdisp rev on
solve age 0.015
save DL.sav