FLAC3D模拟强夯冲击荷载时出现错误，高手们看看怎么回事

半葫酒

config dynamic
set dyn off
generate zone brick p0 0 0 0 p1 10 0 0 p2 0 10 0 p3 0 0 20
model mohr
prop bulk 5.8e6 shear 2.505e10 coh 16 fri 33.3  ran z 18.5 20
prop bulk 2.84e6 shear 3.905e10 coh 11.3 fric 11.3 ran z 13.65 18.5
prop bulk 3.21e6 shear 2.490e10 coh 32.7 fric 14.1 ran z 10.2 13.65
prop bulk 7.62e6 shear 6.503e10 coh 15.7 fric 19.0 ran z 0 10.2
fix x ran x 0
fix x ran x 10
fix y ran y 0
fix y ran y 10
fix z ran z 0
ini dens 1.8 ran z 0 18.5
ini dens 1.5 ran z 18.5 20
ini szz -22.5 grad 0 0 15 ran z 18.5 20
ini sxx -12.25 grad 0 0 7.5 ran z 18.5 20
ini syy -12.25 grad 0 0 7.5 ran z 18.5 20
ini szz -355.5 grad 0 0 18 ran z 0  18.5
ini sxx -177.75 grad 0 0 9 ran z 0  18.5
ini syy -177.75 grad 0 0 9 ran z 0  18.5
ini pp 185 grad 0 0 -10 ran z 0 18.5
set grav 0 0 -10
step 30000
hist unbal
solve
set dyn on
model finn
prop ff_latency=50 ff_switch=0 ff_c1=0.8 ff_c2=0.79 ff_c3=0.45 ff_c4=0.73
def setup
freq=20.0
omega=2.0*pi*freq
pulse=1.0/freq
end
set freq=50
setup
def wave
if dytime>pulse
wave=0.0
else
wave=5.0*sin(omega*dytime)
endif
end
apply zacc=5.0 hist wave range z 19.9 20.1
apply nquiet dquiet squiet range z 19.9 20.1
set dyn damp local 0.314
set dyn  time=0
set large
hist unbal
hist dytime
hist wave
hist gp pp 0 0 18
hist gp pp 0 0 16
hist gp pp 0 0 14
hist gp pp 0 0 12
hist gp pp 0 0 10.5
hist gp pp 0 0 8
sav dongli.sav
plot show
solve
运行后出现这个结果
zero stiffness in grid-point 1
有人能解释怎么回事不？

antonyzhu

建议可把该问题分享到ITASCA官方论坛。
里面的ITASCA工程师也许能给你一些建议。
http://bbs.itasca.cn（开放式注册）

半葫酒

谢谢 2# antonyzhu

lrf136

模型参数没有赋上去