请教高手，我建了一个四分之一的井壁与土体接触的模型，上部施加荷载，两个土层，为什么计算结果时，上面的土层单元压坏，而力却并没有向下传递？谁知道说一下，搞了一个星期了，还是没有找到原因，急呀！！先谢谢了！
下面是命令流：
;creat soil
gen zone radcyl p0 (0,0,0) p1 (30,0,0) p2 (0,10,0) p3 (0,0,30) &
size 4 10 8 8 dim 5.2 5.2 5.2 5.2 ratio 1 1 1 1.15 group soil1
gen zone radcyl p0 (0,10,0) p1 (30,10,0) p2 (0,16,0) p3 (0,10,30) &
size 4 6 8 8 dim 5.2 5.2 5.2 5.2 ratio 1 1 1 1.15 group soil2
gen zone cshell p0 (0,0,0) p1 (5.2,0,0) p2 (0,10,0) p3 (0,0,5.2) &
size 2 10 8 8 dim 4 4 4 4 ratio 1 1 1 1 group shaft1
gen zone cshell p0 (0,10,0) p1 (5.2,10,0) p2 (0,16,0) p3 (0,10,5.2) &
size 2 6 8 8 dim 4 4 4 4 ratio 1 1 1 1 group shaft2
save 1.sav
model mohr ;增大参数使自应力平衡
prop bulk 3e9 shear 2e9 cohesion 2.1e7 friction 45 dens 2500 ran y 0 10
prop bulk 81e6 shear 37e6 cohesion 2.1e7 friction 20 dens 2200 ran y 10 16
fix x y z range y -0.01 0.01
fix x range x -0.01 0.01 any x 29.99 30.01 any
fix z range x -0.01 0.01 any z 29.99 30.01 any
fix x range z -0.01 0.01 any x 29.99 30.01 any
fix z range z -0.01 0.01 any z 29.99 30.01 any
set gravity 0 -10 0
ini syy -1068500. grad 0 534250. 0 range y 0. 10
ini sxx -534250. grad 0 267125. 0 range y 0 10
ini szz -534250. grad 0 267125. 0 range y 0 10
ini syy -818500. grad 0 409250. 0 range y 10 16
ini sxx -409250. grad 0 204625. 0 range y 10 16
ini szz -409250. grad 0 204625. 0 range y 10 16
hist id=999 unbal
plot hist 999
solve
prop bulk 3e9 shear 2e9 cohesion 2.1e7 friction 45 dens 2500 ran group soil1
prop bulk 81e6 shear 37e6 cohesion 0.3e6 friction 20 dens 2200 ran group soil2
solve
save 地应力平衡.sav
ini xdis 0 ydis 0 zdis 0
ini xvel 0 yvel 0 zvel 0
ini state 0
rest 1.sav
del ran group shaft1
del ran group shaft2
interface 1 face range cylinder end1 (0,0,0) end2 (0,10,0) radius 5.3 &
cylinder end1 (0,0,0) end2 (0,16,0) radius 5.1 not
save 2.sav
rest 1.sav
del ran group soil1
del ran group soil2
expgrid 1.flac3d
rest 2.sav
impgrid 1.flac3d
plot inter yel ske
set log on
set logfile jierfenxi.log
set large
model mohr
interface 1 prop kn 30e6 ks 30e6 coh 0.01e6 friction 20
prop bulk 3e9 shear 2e9 cohesion 2.1e7 friction 45 dens 2500 ran group soil1
prop bulk 81e6 shear 37e6 cohesion 0.3e6 friction 20 dens 2200 ran group soil2
prop bulk 15e9 shear 12e9 cohesion 6e6 friction 32 dens 2500 ran group shaft1
prop bulk 15e9 shear 12e9 cohesion 6e6 friction 32 dens 2500 ran group shaft2
apply syy -2.5e6 range group soil2
apply syy -2.6e6 range group shaft2
set mech ratio 1e-5
solve
save jier.sav
