为甚么在加筋土地基承载力计算中，加进土工格栅地基承载力没有提高
下面是我想到的几个原因，地基承载力监控的方式出错，土工格栅的参数设置错误或者是土体的参数设置出错，其他的原因暂时还没有想到，有哪位大神做过这方面的研究，可以运行以下下面的代码。
new
gen zone brick p0 0 0 0 p1 900 0 0 &
p2 0 900 0 p3 0 0 600 size 20 20 12 group shatu
model elas
prop bulk 5e9 shear 1e9
ini dens 1560
fix x range x -0.1 0.1
fix x range x 899.9 900.1
fix y range y -0.1 0.1
fix y range y 899.9 900.1
fix z range z -0.1 0.1
set gravity 0 0 -9.85
so
sel geogrid ele cst range x 100 800 y 100 800 z 599.9 600.1
sel node init zpos add -50
sel geogrid ele cst range x 100 800 y 100 800 z 599.9 600.1
sel node init zpos add -100 range z 599.9 600.1
sel geogrid ele cst range x 100 800 y 100 800 z 599.9 600.1
sel node init zpos add -150 range z 599.9 600.1
sel geogrid ele cst range x 100 800 y 100 800 z 599.9 600.1
sel node init zpos add -200 range z 599.9 600.1
sel geogrid prop iso = (1.1e9,0.3) cs_sk = &
3.65e7 cs_sfric = 34.0 cs_scoh = 3e4 thick = 5e-2
gen zone brick p0 375 375 600 p1 525 375 600 &
p2 375 525 600 p3 375 375 620 size 6 6 1 group jichu
ini z add -10 range group shatu
int 1 face range x 375 525 y 375 525 z 599.9 600.1
ini z add 10 range group shatu
int 1 prop ks 2e9 kn 2e9 fric 15
model mohr range group shatu
model elas range group jichu
prop bulk 2.75e11 shear 0.824e11 density 3000 range group jichu
prop bulk 17e6 shear 8.1e6 fric 34 dilation 34 coh 5e5 dens 1560 range group shatu
ini xvel 0 yvel 0 zvel 0
ini zdis 0 xdis 0 ydis 0
ini zvel -8e-4 range group jichu
fix zvel -8e-4 range group jichu
hist id = 1 gp zdis 450 450 600
hist id = 2 zone szz id = 4590
hist id = 3 zone szz id = 4591
hist id = 4 zone szz id = 4611
hist id = 5 zone szz id = 4828
hist id = 6 sel node xdis id = 269
hist id = 7 sel node xdis id = 267
hist id = 8 sel node xdis id = 263
hist id = 9 sel node xdis id = 257
hist id = 10 sel node xdis id = 272
step 10000
