关于一个动力作用地面的例子讨论
请各位高手指点一下小生，程序有什么错误吗？set log on
new
gen zon bri p0 0 0 -15 p1 15 0 -15 p2 0 15 -15 p3 0 0 0ratio 1.2 1.2 1 size 15 15 15
model mohr
prop bulk 3e7 shear 1e7 coh 1e10 fri 15 tension 1e10
fix z ran z -14.9 -15.1
fix x ran x -.1 .1
fix x ran x 14.9 15.1
fix y ran y -.1 .1
fix y ran y 14.9 15.1
ini dens 2000
set grav 0 0 -10
plot cont szz out on
solve
save 1.sav
config dyn
prop bulk 3e7 shear 1e7 coh 10e3 fri 15 ten 0
def setup
omega = 2.0 * pi * freq
pulse = 1.0 / freq
end
set freq=4.0
setup
def wave
if dytime > pulse
wave = 0.0
else
wave = 0.5 * (1.0 - cos(omega * dytime))
endif
end
apply dquiet squiet range x -0.1 0.1
apply dquiet squiet range x 14.9 15.1
apply dquiet squiet range y -0.1 0.1
apply dquiet squiet range y 14.9 15.1
apply dquiet squiet range z -14.9 -15.1
apply szz -40e5 hist wave syz 0.0 sxz 0.0 range z -0.1 0.1 x 6.5 7 y 6.5 7;-2e5的系数来源于 的值
apply zvel 0range x -0.1 0.1
apply zvel 0range x 14.9 15.1
apply zvel 0range y -0.1 0.1
apply zvel 0range y 14.9 15.1
apply zvel 0range z -14.9 -15.1
hist gp xvel 6.7,6.7,0
hist dytime
hist wave
plot create hhh
plot add hist 1
plot show
plot cont szz out on
solve age 2
sav 2.sav
页:
[1]