求解: FLAC3D基础与工程实例中7.3节 应用实例——路基施工过程模拟
按照以下程序怎么运行不了???总是会弹出错误来,显示,Zero stiffness in grid-point 1 如附图:
new
gen zone brick p0 0 0 -50 p1 27.5 0 -50 p2 0 5 -50 p3 0 0 -10 size 8 1 10 & group caly
gen zone brick p0 27.5 0 -50 p1 100 0 -50 p2 27.5 5 -50 p3 27.5 0 -10 & ratio 1.1 1 1 size 12 1 10 gr clay
gen zone brick p0 0 0 -10 p1 27.5 0 -10 p2 0 5 -10 p3 0 0 0 ratio 1 1 0.8 & size 8 1 4 group soil
gen zone brick p0 27.5 0 -10 p1 100 0 -10 p2 27.5 5 -10 p3 27.5 0 0 &
ratio 1.1 1 0.8 size 12 1 4 group soil
gen zone brick p0 0 0 0 p1 27.5 0 0 p2 0 5 0 p3 0 0 5 p4 27.5 5 0 &
p5 0 5 5 p6 20 0 5 p7 20 5 5 size 8 1 5 group dam
fix x y z ran z -49.9 -50.1
fix x ran x -.1 .1
fix x ran x 99.9 100.1
fix y
mo mo ran z -50 0
mo null ran z 0 5
prop bulk 7.8e6 shear 3.0e6 coh 10e10 tens 1e10 ran gr soil
ini dens 1500 ran gr soil
prop bulk 3.91e6 shear 1.5e6 coh 10e10 tens 1e10 ran gr clay
ini dens 1800 ran gr clay
set grav 0 0 -9.8
hist id=1 unbal
solve
prop bulk 7.8e6 shear 3.0e6 coh 10e3 fric 15 ran gr soil
prop bulk 3.91e6 shear 1.5e6 coh 20e3 fric 20 ran gr clay
solve
save 初始应力.sav
pl con szz ou on
pl con sxx ou on
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
hist id=2 gp zdisp 0 0 0
hist id=3 gp zdisp 27.5 0 0
hist id=4 gp xdisp 27.5 0 0
mo elas ran z 0 1
prop bulk 7.8e6 shear 3.0e6 ran z 0 1
ini den 1500 ran z 0 1
solve
save fill-1.sav
pl con zdisp ou on
pl con xdisp ou on
mo elas ran z 1 2
prop bulk 7.8e6 shear 3.0e6 ran z 1 2
ini den 1500 ran z 1 2
solve
save fill-2.sav
mo elas ran z 2 3
prop bulk 7.8e6 shear 3.0e6 ran z 2 3
ini den 1500 ran z 2 3
solve
save fill-3.sav
mo elas ran z 3 4
prop bulk 7.8e6 shear 3.0e6 ran z 3 4
ini den 1500 ran z 3 4
solve
save fill-4.sav
mo elas ran z 4 5
prop bulk 7.8e6 shear 3.0e6 ran z 4 5
ini den 1500 ran z 4 5
solve
save fill-5.sav
pl con zdisp
pl con xdisp