土工格栅加筋路堤命令流问题

gengmin

最近写了一个土工格栅加筋路堤的flac3d程序，但是发现命令中加了土工格栅层，可是运行起来好像没加土工格栅似的。不知道命令流哪写错了。望各位大侠赐教，谢谢。
命令流如下：
new

gen zone brick p0=(0,0,0) p1=(3,0,0) p2=(0,1,0) p3=(0,0,3) size 2 1 3
gen zone brick p0=(3,0,0) p1=(35,0,0) p2=(3,1,0) p3=(3,0,3) size 17 1 3
gen zone brick p0=(3,0,3) p1=(35,0,3) p2=(3,1,3) p3=(18,0,13) p4=(35,1,3) p5=(18,1,13) p6=(35,0,13) p7=(35,1,13) size 17 1 10

model mohr
property density 2000 bulk 2e8 shear 1e8 cohesion 1e10
property friction 30 dilation 30 tension 0

sel geogrid id=1 range x=(7.5,35) z=(12.99,13.01)
sel node init zpos add 6 range z=(12.99,13.01)
sel geogrid property iso=(26e9,0.33) thick=5e-3 cs_sk=2.3e6 cs_scoh=0.0
sel geogrid property cs_sfric=5.0 ; 0.0

fix x y z range z -0.1 0.1
fix x range x 34.9 35.1
fix x range x -0.1 0.1
fix y

model elastic
property density 2000 bulk 1e8 shear 3e7
set gravity 0,0,-10
solve
initial xdisplacement 0 ydisplacement 0 zdisplacement 0
initial xvelocity 0 yvelocity 0 zvelocity 0
solve
plot block state