求助，路基堆载预压模拟，预压一施加，沉降就变大，流固耦合算不过去。。。

小羊

[求助]FLAC3D模拟路基堆载预压固结，一施加预压就沉降过大，涉及到流固耦合（附带命令流）   Post By：2011-3-20 15:36:47

小弟要模拟路基分层堆载，并进行预压模拟，在分层堆载过程中沉降还算正常，但是一到预压施加，马上沉降就大起来了。。。堆载预压也是分层进行的，刚加第一层的时候沉降就变得非常的大，这个预压堆载的坡度跟路基的坡度是不一样的，建模的时候跟下面的是分开建的，小弟试过坡度不变的话施加堆载，并将两个建为一个块得话，沉降很正常，但是坡度变化了，分开建块的话，沉降立马就大了 小弟百思不得其解啊，请各位大牛帮帮小弟啊。。。。下面有附图。。。。       另外在计算流固耦合的时候，小弟设置的参数不知道对不对，就是算不过去，FLAC命令框里面的那些也不是很明白是什么意思，还望版主大侠给讲解下。。。。命令流是更新过的，里面有流固耦合的命令，王版主帮帮小弟。 下面是命令流： new config fluid gen zone brick p0 0 0 -25 p1 12.4 0 -25 p2 0 3 -25 p3 0 0 -20 size 8 3 4 group zsand plot add surface red plot add axes yellow plot set magnification 0.1 plot show   gen zone brick p0 12.4 0 -25 p1 40 0 -25 p2 12.4 3 -25 p3 12.4 0 -20 ratio 1.1 1 1 size 12 3 4 group zsand gen zone brick p0 0 0 -20 p1 12.4 0 -20 p2 0 3 -20 p3 0 0 -0.5 ratio 1 1 0.98 size 8 3 20 group clay gen zone brick p0 12.4 0 -20 p1 40 0 -20 p2 12.4 3 -20 p3 12.4 0 -0.5 ratio 1.1 1 0.98 size 12 3 20 group clay gen zone brick p0 0 0 -0.5 p1 12.4 0 -0.5 p2 0 3 -0.5 p3 0 0 0 size 8 3 1 group stone gen zone brick p0 12.4 0 -0.5 p1 40 0 -0.5 p2 12.4 3 -0.5 p3 12.4 0 0 ratio 1.1 1 0.95 size 12 3 1 group stone gen zone brick p0 0 0 0 p1 12.4 0 0 p2 0 3 0 p3 0 0 3.3 p4 12.4 3 0 &                p5 0 3 3.3 p6 7.45 0 3.3 p7 7.45 3 3.3 size 8 3 10 group huitian gen zone brick p0 0 0 3.3 p1 7.45 0 3.3 p2 0 3 3.3 p3 0 0 4.1 p4 7.45 3 3.3 &                p5 0 3 4.1 p6 6.25 0 4.1 p7 6.25 3 4.1 size 7 3 3 group huitian2 gen zone reflect n (-1 0 0) o = (0 0 0) plot block group fix x y z ran z -24.9 -25.1 fix x ran x -39.9 -40.1 fix x ran x 39.9 40.1 fix y ran y -0.1 0.1 fix y ran y 2.9 3.1 model mohr ran z -25 0 model null ran z 0 4.1 prop bulk 17.5e6 shear 13.125e6 coh 1e10 tension 1e10 ran group zsand ini dens 2100 ran group zsand prop bulk 7.7e6 shear 1.6e6 coh 1e10 tension 1e10 ran group clay ini dens 1950 ran group clay prop bulk 11e6 shear 8.25e6 coh 1e10 tension 1e10 ran group stone ini dens 2300 ran group stone model fl_iso ini pp 0 grad 0 0 -10e3 ran z -6 -25 ini ftens -1e-3 ini fdens 1000 ini fmod 2e8 prop perm 1e-11 poro 0.45 ran z -6 -25 fix pp 0 range z -5.9 -6.1 fix pp range z -25.1 -24.9 set fluid off set grav 0 0 -9.8 plot con zdis sk solve prop bulk 17.5e6 shear 13.125e6 coh 0 fric 36 ran group zsand prop bulk 7.7e6 shear 1.6e6 coh 10e3 fric 8 ran group clay prop bulk 11e6 shear 8.25e6 coh 0 fric 30 ran group stone solve ini xdis 0 ydis 0 zdis 0 ini xvel 0 yvel 0 zvel 0 ;土体中心创建桩 sel pile id 1 beg=(0.75 0.2 0.0) end (0.75 0.2 -18) nseg 18 sel pile id 1 beg=(2.25 0.2 0.0) end (2.25 0.2 -18) nseg 18 sel pile id 1 beg=(3.75 0.2 0.0) end (3.75 0.2 -18) nseg 18 sel pile id 1 beg=(5.25 0.2 0.0) end (5.25 0.2 -18) nseg 18 sel pile id 1 beg=(6.75 0.2 0.0) end (6.75 0.2 -18) nseg 18 sel pile id 1 beg=(8.25 0.2 0.0) end (8.25 0.2 -18) nseg 18 sel pile id 1 beg=(9.75 0.2 0.0) end (9.75 0.2 -18) nseg 18 sel pile id 1 beg=(11.25 0.2 0.0) end (11.25 0.2 -18) nseg 18 sel pile id 1 beg=(12.75 0.2 0.0) end (12.75 0.2 -18) nseg 18 sel pile id 1 beg=(0.0 1.5 0.0) end (0.0 1.5 -18) nseg 18 sel pile id 1 beg=(1.5 1.5 0.0) end (1.5 1.5 -18) nseg 18 sel pile id 1 beg=(3.0 1.5 0.0) end (3.0 1.5 -18) nseg 18 sel pile id 1 beg=(4.5 1.5 0.0) end (4.5 1.5 -18) nseg 18 sel pile id 1 beg=(6.0 1.5 0.0) end (6.0 1.5 -18) nseg 18 sel pile id 1 beg=(7.5 1.5 0.0) end (7.5 1.5 -18) nseg 18 sel pile id 1 beg=(9.0 1.5 0.0) end (9.0 1.5 -18) nseg 18 sel pile id 1 beg=(10.5 1.5 0.0) end (10.5 1.5 -18) nseg 18 sel pile id 1 beg=(12 1.5 0.0) end (12 1.5 -18) nseg 18 sel pile id 1 beg=(0.75 2.8 0.0) end (0.75 2.8 -18) nseg 18 sel pile id 1 beg=(2.25 2.8 0.0) end (2.25 2.8 -18) nseg 18 sel pile id 1 beg=(3.75 2.8 0.0) end (3.75 2.8 -18) nseg 18 sel pile id 1 beg=(5.25 2.8 0.0) end (5.25 2.8 -18) nseg 18 sel pile id 1 beg=(6.75 2.8 0.0) end (6.75 2.8 -18) nseg 18 sel pile id 1 beg=(8.25 2.8 0.0) end (8.25 2.8 -18) nseg 18 sel pile id 1 beg=(9.75 2.8 0.0) end (9.75 2.8 -18) nseg 18 sel pile id 1 beg=(11.25 2.8 0.0) end (11.25 2.8 -18) nseg 18 sel pile id 1 beg=(12.75 2.8 0.0) end (12.75 2.8 -18) nseg 18 sel pile id 1 beg=(-0.75 0.2 0.0) end (-0.75 0.2 -18) nseg 18 sel pile id 1 beg=(-2.25 0.2 0.0) end (-2.25 0.2 -18) nseg 18 sel pile id 1 beg=(-3.75 0.2 0.0) end (-3.75 0.2 -18) nseg 18 sel pile id 1 beg=(-5.25 0.2 0.0) end (-5.25 0.2 -18) nseg 18 sel pile id 1 beg=(-6.75 0.2 0.0) end (-6.75 0.2 -18) nseg 18 sel pile id 1 beg=(-8.25 0.2 0.0) end (-8.25 0.2 -18) nseg 18 sel pile id 1 beg=(-9.75 0.2 0.0) end (-9.75 0.2 -18) nseg 18 sel pile id 1 beg=(-11.25 0.2 0.0) end (-11.25 0.2 -18) nseg 18 sel pile id 1 beg=(-12.75 0.2 0.0) end (-12.75 0.2 -18) nseg 18 sel pile id 1 beg=(-1.5 1.5 0.0) end (-1.5 1.5 -18) nseg 18 sel pile id 1 beg=(-3.0 1.5 0.0) end (-3.0 1.5 -18) nseg 18 sel pile id 1 beg=(-4.5 1.5 0.0) end (-4.5 1.5 -18) nseg 18 sel pile id 1 beg=(-6.0 1.5 0.0) end (-6.0 1.5 -18) nseg 18 sel pile id 1 beg=(-7.5 1.5 0.0) end (-7.5 1.5 -18) nseg 18 sel pile id 1 beg=(-9.0 1.5 0.0) end (-9.0 1.5 -18) nseg 18 sel pile id 1 beg=(-10.5 1.5 0.0) end (-10.5 1.5 -18) nseg 18 sel pile id 1 beg=(-12 1.5 0.0) end (-12 1.5 -18) nseg 18 sel pile id 1 beg=(-0.75 2.8 0.0) end (-0.75 2.8 -18) nseg 18 sel pile id 1 beg=(-2.25 2.8 0.0) end (-2.25 2.8 -18) nseg 18 sel pile id 1 beg=(-3.75 2.8 0.0) end (-3.75 2.8 -18) nseg 18 sel pile id 1 beg=(-5.25 2.8 0.0) end (-5.25 2.8 -18) nseg 18 sel pile id 1 beg=(-6.75 2.8 0.0) end (-6.75 2.8 -18) nseg 18 sel pile id 1 beg=(-8.25 2.8 0.0) end (-8.25 2.8 -18) nseg 18 sel pile id 1 beg=(-9.75 2.8 0.0) end (-9.75 2.8 -18) nseg 18 sel pile id 1 beg=(-11.25 2.8 0.0) end (-11.25 2.8 -18) nseg 18 sel pile id 1 beg=(-12.75 2.8 0.0) end (-12.75 2.8 -18) nseg 18 sel pile id=1 prop emod=8e10 nu=0.3 xcarea=1.256e-3 & xcj=1.57e-4 xciy=7.85e-5 xciz=7.85e-5 & per=1.256 dens 2500 & cs_sk=1.3e11 cs_scoh=1e4 cs_sfric=12 & cs_nk=1.3e11 cs_ncoh=0 cs_nfric=0 & cs_ngap=off sel geogrid id 1  range x = (-11.9 11.9 ) z = (-0.01 0.01) sel node init zpos add -0.3 range z = (-0.01 0.01) sel geogrid id 2  range x = (-11.9 11.9 ) z = (-0.01 0.01) sel node init zpos add -0.15 range z = (-0.01 0.01) sel geogrid prop iso = (26e9 0.33) thick 5e-3 &     cs_sk = 2.3e6 cs_scoh = 0.0 sel geogrid prop cs_sfric 30 set fluid on set mech force 1.5e3 set mech subs 10000 auto set fluid subs 30 def day    day=fltime/24/3600 end day hist id=1 day model elastic ran z 0 0.8 ;激活0 m ~ 0.8 m的单元 prop bulk 5e6 shear 7.5e6 ran z 0 0.8 ini dens 1900 ran z 0 0.8 hist id=2 gp zdis 0 0 0         ;记录底部中心点的沉降 hist id=3 gp zdis 0 0 0.8         ;记录顶部中心点的沉降 hist id=4 gp zdis 12.4 0 0 ;记录坡脚处的沉降 hist id=5 gp xdis 12.4 0 0  ;记录坡脚处的水平位移 hist id=6 gp zdis 11.2 0 0.8 ;记录路肩处的沉降 hist id=7 gp xdis 11.2 0 0.8 ;记录路肩处的水平位移 solve age 7.776e6                 ;按软件默认精度求解 model elastic ran z 0.8 1.6 ;激活0.8m ~ 1.6m的单元 prop bulk 7e6 shear 7.5e6 ran z 0.8 1.6 ini dens 1900 ran z 0.8 1.6 hist id=8 gp zdis 0 0 0         ;记录底部中心点的沉降 hist id=9 gp zdis 0 0 1.6         ;记录顶部中心点的沉降 hist id=10 gp zdis 12.4 0 0 ;记录坡脚处的沉降 hist id=11 gp xdis 12.4 0 0  ;记录坡脚处的水平位移 hist id=12 gp zdis 10 0 1.6 ;记录路肩处的沉降 hist id=13 gp xdis 10 0 1.6  ;记录路肩处的水平位移 solve age 7.776e6                         ;按软件默认精度求解 model elastic ran z 1.6 2.4   ;激活1.6m ~ 2.4m的单元 prop bulk 7e6 shear 7.5e6 ran z 1.6 2.4 ini dens 1900 ran z 1.6 2.4 hist id=14 gp zdis 0 0 0        ;记录底部中心点的沉降 hist id=15 gp zdis 0 0 3        ;记录顶部中心点的沉降 hist id=16 gp zdis 12.4 0 0 ;记录坡脚处的沉降 hist id=17 gp xdis 12.4 0 0  ;记录坡脚处的水平位移 hist id=18 gp zdis 8.8 0 2.4 ;记录路肩处的沉降 hist id=19 gp zdis 8.8 0 2.4  ;记录路肩处的水平位移 solve age 7.776e6                         ;按软件默认精度求解 model elastic ran z 2.4 3.3 ;激活2.4m ~ 3.3m的单元 prop bulk 7e6 shear 7.5e6 ran z 2.4 3.3 ini dens 1900 ran z 2.4 3.3 hist id=20 gp zdis 0 0 0 ;记录底部中心点的沉降 hist id=21 gp zdis 0 0 3.3      ;记录顶部中心点的沉降 hist id=22 gp zdis 12.4 0 0 ;记录坡脚处的沉降 hist id=23 gp xdis 12.4 0 0  ;记录坡脚处的水平位移 hist id=24 gp zdis 7.45 0 3.3 ;记录路肩处的沉降 hist id=25 gp zdis 7.45 0 3.3  ;记录路肩处的水平位移 solve age 7.776e6                   ;按软件默认精度求解 model elastic ran z 3.3 3.95 ;激活3.3m ~ 3.95m的单元 prop bulk 7e6 shear 7.5e6 ran z 3.3 3.95 ini dens 1900 ran z 3.3 3.95 hist id=26 gp zdis 0 0 0 ;记录底部中心点的沉降 hist id=27 gp zdis 0 0 3.95     ;记录顶部中心点的沉降 hist id=28 gp zdis 12.4 0 0 ;记录坡脚处的沉降 hist id=29 gp xdis 12.4 0 0  ;记录坡脚处的水平位移 hist id=30 gp zdis 6.475 0 3.95 ;记录路肩处的沉降 hist id=31 gp zdis 6.475 0 3.95 ;记录路肩处的水平位移 solve age 7.776e6                           ;按软件默认精度求解

小羊

自己占个SF吧。。。

jun_118899

把命令用附件加在后面，要不很难检查~

jun_118899

不过看你那个图，也还算正常吧

小羊

4# jun_118899
恩 那个图是正常的图 一加上面的预压荷载之后就不正常了啊 沉降就变得很大啊。。。。然后是流固耦合的东西还不对。。。。请帮帮小弟吧。。。。

小羊

5# 小羊 命令流已经 传上去了 帮小弟看看吧 万分感谢啊

zhq19861225

你的初始地应力计算就有问题，刚开始没有把流体模量设置为零会引起超空隙水压力

hyyfj

你好，我也是做这方面的，请问你的问题解决了吗？能共享一下吗？谢谢！

颠倒

前辈啊高手！膜拜。。。
求指点。。。

binggehaoyang

bucuo     

jkidd23

学习中，研究一下

ling1227

做的挺好的了啊

anningjhj

我是来学习的

信徒

厉害厉害，我也想做