格宾挡墙格宾单元的建模，求指点
代码：;====格宾挡墙====================================================================================================
new
set log on
;======地基土==============================================================================================
ge zo brickp0 0 0 -4p1 16.2 0 -4p2 0 2 -4p3 0 0 0size 81 2 4
ge zo brickp0 -10 0 -4p1 0 0 -4p2 -10 2 -4p3 -10 0 0size 50 2 4
ge zo brickp0 16.2 0 -4p1 21.2 0 -4 p2 16.2 2 -4p3 16.2 0 0size 25 2 4
ge zo brickp0 21.2 0 -4 p1 22.20 0 -4p2 21.2 2 -4p3 21.2 0 0size 5 2 4
ge zo brickp0 22.2 0 -4 p1 22.35 0 -4p2 22.2 2 -4p3 22.2 0 0size 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 -7p1 16.2 0 -7p2 0 2 -7p3 0 0 -4size 81 2 3
ge zo brick p0 -10 0 -7p1 0 0 -7p2 -10 2 -7p3 -10 0 -4size 50 2 3
ge zo brick p0 16.2 0 -7p1 21.2 0 -7p2 16.2 2 -7p3 16.2 0 -4size 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=0zdisplacement=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 1size 4 2 2
;=========结构填土表面格宾网==============================================================================
sel geogrid id=2 rangex 0 y 0 2z 0 1
sel node local xdir 0 0 1ydir 0 1 0 range x 0
sel node fix lsys range x 0
sel geogrid id=2 range x 0.8 y 0 2z 0 1
sel node local xdir 0 0 1ydir 0 1 0range x 0.8
sel node fix lsys range x 0.8
sel geogrid id=2 range x 0 0.8y 0 2z 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 rangez 0
sel node fix lsys rangez 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 1size 77 2 2
ge zo brick p0 16.2 0 0p1 21.2 0 0 p216.2 2 0p3 16.2 0 1size 25 2 2
ge zo brick p0 21.2 0 0p1 22.20 0 0p2 21.2 2 0 p3 21.2 0 1size 5 2 2
ge zo brick p0 22.2 0 0p1 22.35 0 0p2 22.2 2 0 p3 22.2 0 1size 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单元应力应变图时，如图提示
是不是geogrid单元不能定义竖向的
sel geogrid id=2 rangex 0 y 0 2z 0 1
sel node local xdir 0 0 1ydir 0 1 0 range x 0
sel node fix lsys range x 0
sel geogrid id=2 range x 0.8 y 0 2z 0 1
sel node local xdir 0 0 1ydir 0 1 0range x 0.8
sel node fix lsys range x 0.8 有大神建过格宾单元的模型吗？可否提供代码以供学习 我最近也在学习FLAC，从上面的命令流来看，首先，没有生成初始地应力，不知道土体有没有拉裂，处于塑性状态的，另外，格栅单元的建立只将范围确定到具体数值可不可行，是不是应该是个范围。有好多我也不懂，我们可以相互交流下。QQ745327686 谢谢楼主的分享
页:
[1]