格宾挡墙格宾单元的建模，求指点

jianndann

代码：

;====格宾挡墙====================================================================================================
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单元应力应变图时，如图提示
是不是geogrid单元不能定义竖向的
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

jianndann

有大神建过格宾单元的模型吗？可否提供代码以供学习

t413357530

我最近也在学习FLAC，从上面的命令流来看，首先，没有生成初始地应力，不知道土体有没有拉裂，处于塑性状态的，另外，格栅单元的建立只将范围确定到具体数值可不可行，是不是应该是个范围。有好多我也不懂，我们可以相互交流下。QQ745327686

mbhfz176cm

谢谢楼主的分享