代码：
;====格宾挡墙====================================================================================================
new
set log on
;======地基土==============================================================================================
ge zo brick p0 0 0 -4 p1 16.2 0 -4 p2 0 2 -4 p3 0 0 0 size 81 2 4
ge zo brick p0 -10 0 -4 p1 0 0 -4 p2 -10 2 -4 p3 -10 0 0 size 50 2 4
ge zo brick p0 16.2 0 -4 p1 21.2 0 -4 p2 16.2 2 -4 p3 16.2 0 0 size 25 2 4
ge zo brick p0 21.2 0 -4 p1 22.20 0 -4 p2 21.2 2 -4 p3 21.2 0 0 size 5 2 4
ge zo brick p0 22.2 0 -4 p1 22.35 0 -4 p2 22.2 2 -4 p3 22.2 0 0 size 1 2 4
group niantu_diji range x -10 22.35 y 0 2 z -4 0
model mohr range group niantu_diji
property bulk 1.39e6 shear 1.04e6 friction 12 cohesion 3e4 range group niantu_diji
ini density 1835 range group niantu_diji
ge zo brick p0 0 0 -7 p1 16.2 0 -7 p2 0 2 -7 p3 0 0 -4 size 81 2 3
ge zo brick p0 -10 0 -7 p1 0 0 -7 p2 -10 2 -7 p3 -10 0 -4 size 50 2 3
ge zo brick p0 16.2 0 -7 p1 21.2 0 -7 p2 16.2 2 -7 p3 16.2 0 -4 size 25 2 3
ge zo brick p0 21.2 0 -7 p1 22.20 0 -7 p2 21.2 2 -7 p3 21.2 0 -4 size 5 2 3
ge zo brick p0 22.2 0 -7 p1 22.35 0 -7 p2 22.2 2 -7 p3 22.2 0 -4 size 1 2 3
group shayan_diji range x -10 22.35 y 0 2 z -7 -4
model mohr range group shayan_diji
property bulk 2.26e10 shear 1.11e10 friction 28 cohesion 1.1e5 tension 1.58e6 range group shayan_diji
ini density 2599 range group shayan_diji
set grav 0 0 -9.81
fix x z range z -7
fix x range x -10
fix x range x 22.35
fix y range z -7 0
solve
pause
;continue
initial xdisplacement=0 ydisplacement=0 zdisplacement=0 range z -7 0
;========第一层施工========================================================================================
;======先底层钢丝网========================================================================================
sel geogrid id=2 range x 0 0.8 y 0 2 z 0
sel geogrid id=2 range x 0.8 3 y 0 2 z 0
;====结构填土=============================================================================================
ge zo brick p0 0 0 0 p1 0.8 0 0 p2 0 2 0 p3 0 0 1 size 4 2 2
;=========结构填土表面格宾网==============================================================================
sel geogrid id=2 range x 0 y 0 2 z 0 1
sel node local xdir 0 0 1 ydir 0 1 0 range x 0
sel node fix lsys range x 0
sel geogrid id=2 range x 0.8 y 0 2 z 0 1
sel node local xdir 0 0 1 ydir 0 1 0 range x 0.8
sel node fix lsys range x 0.8
sel geogrid id=2 range x 0 0.8 y 0 2 z 1
sel node local xdir 1 0 0 ydir 0 1 0 range z 1
sel node fix lsys range z 1
sel node local xdir 1 0 0 ydir 0 1 0 range z 0
sel node fix lsys range z 0
sel geogrid id=2 prop iso=(8.0e8,0.33) thick=2.7e-3 cs_sk=1.5e6 cs_scoh=1.03e4 cs_sfric=22
;=======面墙后填土======================================================================================
ge zo brick p0 0.8 0 0 p1 16.2 0 0 p2 0.8 2 0 p3 0.8 0 1 size 77 2 2
ge zo brick p0 16.2 0 0 p1 21.2 0 0 p2 16.2 2 0 p3 16.2 0 1 size 25 2 2
ge zo brick p0 21.2 0 0 p1 22.20 0 0 p2 21.2 2 0 p3 21.2 0 1 size 5 2 2
ge zo brick p0 22.2 0 0 p1 22.35 0 0 p2 22.2 2 0 p3 22.2 0 1 size 1 2 2
group T1_01 range x 0 22.35 y 0 2 z 0 1
model mohr range group T1_01
property bulk 5.14e7 shear 1.34e7 friction 38 cohesion 1.4e4 range group T1_01
ini density 2038.7 range group T1_01
set grav 0 0 -9.81
fix x range x 22.35 z 0 1
fix y range z 0 1
history nstep 5
history unbalance
history gp zdisplacement 1 0 0.5
plot history 1
plot history 2
solve
运行没问题，但在绘制geogrid单元应力应变图时，提示
，求大神给予指点
