斑竹一定帮帮我呀！土工格栅加筋挡墙模型计算出问题

刘德贵

gen zone brick size 36 2 12 p0 0 0 0 p1 18 0 0 p2 0 1 0 p3 0 0 6  
group base range z=-0.1 6.1
model mohr range group base
prop dens=2000 bulk=20.84e6 shear=9.62e6 fric=25 coh=0 dila=0
interface 1 face range x=6 18 z=5.9 6.1
interface 1 prop kn=1.12e9 ks=2e6 fric=30 coh=0 tens=0 dila=0
gen zone brick size 24 2 12 p0 6 0 6 p1 18 0 6 p2 6 1 6 p3 6 0 12
group re range z=6.0 12
model null range group re
fix z range z=-0.1 0.1
fix y range y=-0.1 0.1
fix y range y=0.9 1.1
fix x range x=-0.1 0.1
fix x range x=17.9 18.1
set grav 0 0 -9.81
ini szz=-1.2e5 grad 0 0 2.0e4 range z=0.0 6.0
ini sxx=-0.474e5 grad 0 0 0.79e4 range z=0.0 6.0
ini syy=-0.474e5 grad 0 0 0.79e4 range z=0.0 6.0
hist gp xdisp 6 5 6.0
hist gp xdisp 6 5 6.5
save reinstage1
restore reinstage1
model mohr range x=6.0 18.0 z=6.0 7.0
prop dens=1900 bulk=41.74e6 shear=19.23e6 fric=35 coh=0 dila=7
sel geogrid id=1 range x=6.0 11.0 z=6.99 7.01
sel node init zpos add -5.5 range z=6.99 7.01
sel geogrid prop iso(160e6 0.4) cs_sk=2.0e6 cs_scoh=0.0 cs_sfric=30 thick=5e3
hist gp xdisp 6 5 7.0
solve
save reinstage2
restore reinstage2
model mohr range x=6.0 18.0 z=7.0 8.0
prop dens=1900 bulk=41.74e6 shear=19.23e6 fric=35 coh=0 dila=7
sel geogrid id=2 range x=6.0 11.0 z=7.99 8.01
sel node init zpos add -1.0 range z=7.99 8.01
sel geogrid id=2 range x=6.0 11.0 z=7.99 8.01
sel node init zpos add -0.5 range z=7.99 8.01
sel geogrid prop iso(160e6 0.4) cs_sk=2.0e6 cs_scoh=0.0 cs_sfric=30 thick=5e3
hist gp xdisp 6 5 7.5
hist gp xdisp 6 5 8.0
solve
save reinstage3
restore reinstage3
model mohr range x=6.0 18.0 z=8.0 9.0
prop dens=1900 bulk=41.74e6 shear=19.23e6 fric=35 coh=0 dila=7
sel geogrid id=3 range x=6.0 11.0 z=8.99 9.01
sel node init zpos add -1.0 range z=8.99 9.01
sel geogrid id=3 range x=6.0 11.0 z=8.99 9.01
sel node init zpos add -0.5 range z=8.99 9.01
sel geogrid prop iso(160e6 0.4) cs_sk=2.0e6 cs_scoh=0.0 cs_sfric=30 thick=5e3
hist gp xdisp 6 5 8.5
hist gp xdisp 6 5 9.0
solve
save reinstage4
restore reinstage4
model mohr range x=6.0 18.0 z=9.0 10.0
prop dens=1900 bulk=41.74e6 shear=19.23e6 fric=35 coh=0 dila=7
sel geogrid id=4 range x=6.0 11.0 z=9.99 10.01
sel node init zpos add -1.0 range z=9.99 10.01
sel geogrid id=4 range x=6.0 11.0 z=9.99 10.01
sel node init zpos add -0.5 range z=9.99 10.01
sel geogrid prop iso(160e6 0.4) cs_sk=2.0e6 cs_scoh=0.0 cs_sfric=30 thick=5e3
hist gp xdisp 6 5 9.5
hist gp xdisp 6 5 10.0
solve
save reinstage5
restore reinstage5
model mohr range x=6.0 18.0 z=10.0 11.0
prop dens=1900 bulk=41.74e6 shear=19.23e6 fric=35 coh=0 dila=7
sel geogrid id=5 range x=6.0 11.0 z=10.99 11.01
sel node init zpos add -1.0 range z=10.99 11.01
sel geogrid id=5 range x=6.0 11.0 z=10.99 11.01
sel node init zpos add -0.5 range z=10.99 11.01
sel geogrid prop iso(160e6 0.4) cs_sk=2.0e6 cs_scoh=0.0 cs_sfric=30 thick=5e3
hist gp xdisp 6 5 10.5
hist gp xdisp 6 5 11.0
solve
save reinstage6
restore reinstage6
model mohr range x=6.0 18.0 z=11.0 12.0
prop dens=1900 bulk=41.74e6 shear=19.23e6 fric=35 coh=0 dila=7
sel geogrid id=6 range x=6.0 11.0 z=11.99 12.01
sel node init zpos add -1.0 range z=11.99 12.01
sel geogrid id=6 range x=6.0 11.0 z=11.99 12.01
sel node init zpos add -0.5 range z=11.99 12.01
sel geogrid prop iso(160e6 0.4) cs_sk=2.0e6 cs_scoh=0.0 cs_sfric=30 thick=5e3
hist gp xdisp 6 5 11.5
hist gp xdisp 6 5 12.0
solve
save reinstage7
以上是模型命令流，但是计算出来简直就乱完了，　不知道为什么／百思不得其解！请求斑竹及高手指点，帮我看看到底什么问题呢？
模型主要问题是出在分阶段填筑填土和土工格栅上！但是不知道是怎么回事！

huahua2005

有三点，
一是明确你的挡土墙形式，土工格栅加筋挡土墙主要有三种形式，一是面板和筋材组成，筋材与面板之间有连接，面板受力作用；二是只有筋材没有面板，但此时土工格栅需采用包裹式；三是面板和筋材组成，筋材包裹式，面板起防护作用，你的好像不在此列；
二是明确施工工序，加筋材前，其下土是否有一定压实，你的情况沉降已经太大了；
三是建议你在新加一层土前，将前一层土位移和土工格栅位移用FISH记录，然后清零，最后叠加，否则新加土与旧土之间的连接会有一些怪异的现象，不过这一点仁者见仁，智者见智，你可以不采纳 。

eagleming

　　你最好　能把命令流，写个文本，弄到附件上。
　我们好下下来，看看哪里有问题。
　你这样贴上　，你试试看怎么弄下来，运行

刘德贵

谢谢好的，我弄下来！

刘德贵

原帖由 刘德贵 于 2007-6-26 08:53 发表
谢谢好的，我弄下来！

多谢斑竹了　！

刘德贵

　　谢谢斑竹！你的意见我觉得很中肯，我现在是在试着建立模型看是不是符合实际情况！但是就您上面的意见，建议在新加一层土前，将前一层土位移和土工格栅位移用FISH记录，然后清零，最后叠加，否则新加土与旧土之间的连接会有一些怪异的现象，这一点我真还不知道怎么弄！可不可以教我看怎么弄！
　　再者，如果我如果在建立面板用是面板和筋材组成，筋材与面板之间有连接，面板受力作用，　那是不是就是建立一个面板模块，然后用interface和土与土工格栅连接就可以了呢！
　　我等会把我的不考虑施工工程的小程序给你看看，看是不是就那样就可以了！
　　非常感谢斑竹！

刘德贵

　　斑竹再帮我看看，我这样处理面板和土工格栅的接触和连接可不可以，接触连接就是用了interface和对面板的模型定义为弹性的材料！

刘德贵

原帖由 刘德贵 于 2007-6-26 15:17 发表
     斑竹再帮我看看，我这样处理面板和土工格栅的接触和连接可不可以，接触连接就是用了interface和对面板的模型定义为弹性的材料！

   

[ 本帖最后由 刘德贵 于 2007-6-26 15:27 编辑 ]

刘德贵

斑竹怎么往上面发图片呢，我的ＷＯＲＤ上编辑的图片怎么发不上去呢！

huahua2005

建议地基做初始平衡，其余的思路好像还可以，没仔细看，只运行了一遍，感觉沉降太夸张了~~

婴儿乐

看看，正在做。

zhangyuping

楼主能不能把你的全部命令流发给我看看啊，QQ：398446713