请高手看看我的命令哪里有问题

kritty

我模拟的是一个桩网结构加固路基的简单模型。取了半结构，桩间距是2m。主要想观察桩的变形。
在填土时说有一个zone没定义模型。
盼高手给指正！谢谢！

kritty

这个是我写的命令，大家帮忙看看错在哪里了。
new
;确定各单元组
gen zone brick p0 0 0 0 p1 7.5 0 0 p2 0 6 0 p3 0 0 3 &
               p4 7.5 6 0 p5 0 6 3 p6 3 0 3 p7 3 6 3 size 25 20 5

group AB  ran z 0.6 3
group geogrid ran z 0 0.6
gen zone brick p0 0 0 -5 p1 7.5 0 -5 p2 0 6 -5 p3 0 0 0 size 25 20 10
group soil ran z -1 0
group clay ran z -5 -1
group pile ran x 0   0.6 y 0.4 1.6 z -0.5 0
group pile ran x 0   0.6 y 2.4 3.6 z -0.5 0
group pile ran x 0   0.6 y 4.4 5.6 z -0.5 0
group pile ran x 1.4 2.6 y 0.4 1.6 z -0.5 0
group pile ran x 1.4 2.6 y 2.4 3.6 z -0.5 0
group pile ran x 1.4 2.6 y 4.4 5.6 z -0.5 0
group pile ran x 0   0.3 y 0.7 1.3 z -4.5 -0.5  
group pile ran x 0   0.3 y 2.7 3.3 z -4.5 -0.5
group pile ran x 0   0.3 y 4.7 5.3 z -4.5 -0.5
group pile ran x 1.7 2.3 y 0.7 1.3 z -4.5 -0.5
group pile ran x 1.7 2.3 y 2.7 3.3 z -4.5 -0.5
group pile ran x 1.7 2.3 y 4.7 5.3 z -4.5 -0.5

save 各单元组.sav
plot blo gro
pause
;边界条件
fix x y z ran z -4.9 -5.1
fix x ran x -.1 .1
fix x ran x 7.4 7.6
fix y ran y -0.1 0.1
fix y ran y 5.9 6.1
;参数
model mohr ran z -5 0
model null ran z 0 3

prop bulk 2.8e7 shear 1.68e7 coh 10e10 tension 1e10 ran group soil
ini dens 2230 ran group soil
prop bulk 6.73e6 shear 2.24e6 coh 10e10 tension 1e10 ran group clay
ini dens 2040 ran group clay
set grav 0 0 -9.8
hist id=1 unbal

prop bulk 2.8e7 shear 1.68e7  coh 10e3 fric 20 ran group soil
prop bulk 6.73e6 shear 2.24e6   coh 32.07e3 fric 11.05 ran group clay
model mohr ran group pile
prop bulk 29.2e9 sh 13.46e9 coh 10e3 fric 15  ran group pile
ini dens 2000 ran group pile
solve

plot contour  szz
plot add hist 1


;建立各接触面（导来导去法）
;格栅接触面
del ran group geogrid not
interface 1 face  range x 0     y 0 6 z 0 0.6 &
                        x 0 7.5 y 0 6 z 0 &
                        x 0 6.6 y 0 6 z 0.6  
restore 各单元组.sav
del ran group geogrid
expgrid 各单元组.flac3d  
impgrid 各单元组.flac3d   
plo inter yel ske
pause                              
;桩接触面
del ran group pile not
interface 2 face ran x 0     y 0.4 1.6 z -0.5 0 &
                     x 0.6   y 0.4 1.6 z -0.5 0  &                  
                     x 0 0.6 y 0.4     z -0.5 0 &
                     x 0 0.6 y 1.6     z -0.5 0 &
                     x 0 0.6 y 0.4 1.6 z 0 &
                     x 0 0.6 y 0.4 1.6 z -0.5
interface 3 face ran x 0     y 2.4 3.6 z -0.5 0 &
                     x 0.6   y 2.4 3.6 z -0.5 0 &
                     x 0 0.6 y 2.4     z -0.5 0 &
                     x 0 0.6 y 3.6     z -0.5 0 &
                     x 0 0.6 y 2.4 3.6 z 0 &
                     x 0 0.6 y 2.4 3.6 z -0.5
interface 4 face ran x 0     y 4.4 5.6 z -0.5 0 &
                     x 0.6   y 4.4 5.6 z -0.5 0 &
                     x 0 0.6 y 4.4     z -0.5 0 &
                     x 0 0.6 y 5.6     z -0.5 0 &
                     x 0 0.6 y 4.4 5.6 z 0 &
                     x 0 0.6 y 4.4 5.6 z -0.5
interface 5 face ran x 1.4     y 0.4 1.6 z -0.5 0 &
                     x 2.6     y 0.4 1.6 z -0.5 0 &
                     x 1.4 2.6 y 0.4     z -0.5 0 &
                     x 1.4 2.6 y 1.6     z -0.5 0 &
                     x 1.4 2.6 y 0.4 1.6 z 0 &
                     x 1.4 2.6 y 0.4 1.6 z -0.5
interface 6 face ran x 1.4     y 2.4 3.6 z -0.5 0 &
                     x 2.6     y 2.4 3.6 z -0.5 0 &
                     x 1.4 2.6 y 2.4     z -0.5 0 &
                     x 1.4 2.6 y 3.6     z -0.5 0 &
                     x 1.4 2.6 y 2.4 3.6 z 0  &
                     x 1.4 2.6 y 2.4 3.6 z -0.5
interface 7 face ran x 1.4     y 4.4 5.6 z -0.5 0 &
                     x 2.6     y 4.4 5.6 z -0.5 0 &
                     x 1.4 2.6 y 4.4     z -0.5 0 &
                     x 1.4 2.6 y 5.6     z -0.5 0  &
                     x 1.4 2.6 y 4.4 5.6 z 0 &
                     x 1.4 2.6 y 4.4 5.6 z -0.5
            
interface 8 face ran x 0     y 0.7 1.3 z -4.5 -0.5 &
                     x 0.3   y 0.7 1.3 z -4.5 -0.5 &
                     x 0 0.3 y 0.7     z -4.5 -0.5 &
                     x 0 0.3 y 1.3     z -0.5 0 &
                     x 0 0.3 y 0.7 1.3 z -0.5 &
                     x 0 0.3 y 0.7 0.3 z -4.5
interface 9 face ran x 0    y 2.7 3.3 z -4.5 -0.5 &
                     x 0.3   y 2.7 3.3 z -4.5 -0.5  &
                     x 0 0.3 y 2.7     z -4.5 -0.5 &
                     x 0 0.3 y 3.3     z -0.5 0 &
                     x 0 0.3 y 2.7 3.3 z -0.5 &
                     x 0 0.3 y 2.7 3.3 z -4.5                                                            
interface 10 face ran x 0    y 4.7 5.3 z -4.5 -0.5 &
                     x 0.3   y 4.7 5.3 z -4.5 -0.5 &
                     x 0 0.3 y 4.7     z -4.5 -0.5 &
                     x 0 0.3 y 5.3     z -0.5 0 &
                     x 0 0.3 y 4.7 5.3 z -0.5 &
                     x 0 0.3 y 4.7 5.3 z -4.5
interface 11 face ran x 1.7    y 0.7 1.3 z -4.5 -0.5 &
                     x 2.3     y 0.7 1.3 z -4.5 -0.5 &
                     x 1.7 2.3 y 0.7     z -4.5 -0.5 &
                     x 1.7 2.3 y 1.3     z -0.5 0 &
                     x 1.7 2.3 y 0.7 1.3 z -0.5 &
                     x 1.7 2.3 y 0.7 1.3 z -4.5                     
interface 12 face ran x 1.7   y 2.7 3.3 z -4.5 -0.5 &
                     x 2.3     y 2.7 3.3 z -4.5 -0.5 &
                     x 1.7 2.3 y 2.7     z -4.5 -0.5 &
                     x 1.7 2.3 y 3.3     z -0.5 0 &
                     x 1.7 2.3 y 2.7 3.3 z -0.5 &
                     x 1.7 2.3 y 2.7 3.3 z -4.5                     
interface 13 face ran x 1.7    y 4.7 5.3 z -4.5 -0.5 &
                      x 2.3    y 4.7 5.3 z -4.5 -0.5 &
                     x 1.7 2.3 y 4.7     z -4.5 -0.5 &
                     x 1.7 2.3 y 5.3     z -0.5 0 &
                     x 1.7 2.3 y 4.7 5.3 z -0.5 &
                     x 1.7 2.3 y 4.7 5.3 z -4.5   
restore 各单元组.sav
del ran group pile
expgrid 各单元组.flac3d  
impgrid 各单元组.flac3d   
plo inter red ske
app nstress -200e5 ran x 0 6 y 0 6 z 3
solve
pause
;回填土
;第一次填筑(格栅)
model elastic  ran z 0 3
ini xdis 0 yd 0 zd 0
ini xvel 0 yvel 0  zvel 0
hist id=2 gp zdis 0 0 0
hist id=3 gp zdis 7.5 0 0
hist id=4 gp xdis 7.5 0 0

model elastic  ran z 0 0.6 group geogrid
prop bulk 10e6 sh 6e6 ran  z 0 0.6
ini dens 2000  ran   z 0 0.6
solve
save 001.sav
plot con zdis ou on
plot con xdis ou on
plot blo gro
;第二次填筑
model elastic  ran z 0.6 2
prop bulk 23.3e6 sh 14e6 ran z 0.6 2
ini dens 1500 ran z 0.6  2
solve
save 002.sav
;第三次填筑
model elastic  ran z 2 3
prop bulk 23.3e6 sh 14e6 ran z 2 3
ini dens 1500 ran z  2 3
solve
save 003.sav
plot con zdis
plot con xdis
plot add hist  3 & 4