flac3d出现illegal geometry和运行的时候zero devide
命令流如下，主要是要模拟一个桩基础下岩层纵向应力的分布情况，故先做了一个钻周土体，然后做了桩，并用interface连接，最后在做桩下岩层gen zone radcy size 3 15 5 20 p0 0 0 51 p1 40 0 51 p2 0 0 1 p3 0 40 51 p4 40 0 1 p5 0 40 1 p6 40 40 51 p7 40 40 1 p8 5 0 51 p9 0 5 51 p10 5 0 1 p11 0 5 1
gen zone reflect norm 1 0 0
gen zone reflect norm 0 1 0
interface 1 face range cylinder end1 0 0 0.9 end2 0 0 51.1 radius 5.1 cylinder end1 0 0 0.9 end2 0 0 51.1 radius 4.9 not
gen zone cyl p0 0 0 102 p1 5 0 102 p2 0 0 51 p3 0 5 102 p4 5 0 51 p5 0 5 51 size 1 20 3 group pile
gen zone reflect norm 1 0 0 range group pile
gen zone reflect norm 0 1 0 range group pile
ini z add -50 range group pile
group soil range z 20 51 group pile not
group rock1 range z 1 20 group pile not
gen zone brick size 10 10 4 p0 -40 -40 -6 p1 40 -40 -6 p2 -40 40 -6 p3 -40 -40 0 group rock2
ini z add 1.0 range group rock2
interface 2 face range cylinder end1 0 0 0.9 end2 0 0 1.1 radius 5.1
plot add grid yellow range group soil
plot add grid red range group rock1
plot add grid green range group rock2
plot add grid blue range group pile
plot add interface blue
plot show
model mohr range group pile not
pro bu 0.53e9 sh 0.33e9 coh 4.387e6 fric 38 tension 1.63e6 dens 2200 range group rock1
pro bu 0.54e9 sh 0.35e9 coh 4.154e6 fric 40 tension 2.43e6dens 2500 range group rock2
pro bu 0.43e9 sh 0.14e9 coh 0.55e6 fric 20 tension 1.15e6 dens 1700 range group soil
model elas range group pile
pro bu 0.77e9 sh 0.44e9 dens 3000 range group pile
fix x range x -39.9 -40.1 39.9 40.1
fix y range y -39.9 -40.1 39.9 40.1
fix x y z range z -5.1 -4.9
set gravity 0 0 -9.81
inter 1 prop kn 300e9 ks 300e9 tens 1e10 sbratio 1
inter 2 prop kn 300e9 ks 300e9 tens 1e10 sbratio 1
hist id 1 unba
hist id 2 gp zdisp 0 0 51
plot add hist 1
step 15000
save 7.sav
ini szz -3.5e6 range z 51.9 52.1
in xdis 0 ydis 0 zdis 0
in xvel 0 yvel 0 zvel 0
set large
step 15000
plot set plane ori 0 0 0 norm 0 1 0
plot contour szz plane
据说这个模型问题很多，我是初级使用者，不大了解问题出现在哪儿，怎么解决，希望得到各位的帮助
页:
[1]