求FLAC3D模拟灌注桩后压浆的相关命令流

岩子228

各位前辈大家好，小弟先继续用FLAC3D模拟经过后压浆处理的灌注桩的静载试验过程，由于软件模拟经验较少，很难入手，如果各位大神谁有相关的命令流或者相关的指导，还望不吝赐教，定当不胜感激~先行谢过！如能提供有效帮助，赠给所有仿真币以表感激

岩子228

这么冷清，只能自己顶一下了。已慢慢开始上手建模了～自我鼓励一下！

岩子228

有兴趣的同学可以一起交流～现在已建立几何模型，正在土层分层，建立接触面

岩子228

坚持求教，也是对自己的一种鞭策！

岩子228

在遇到的问题是，多层土的地应力场形成。之考虑重力场的情况下，为什么第二层土的初始应力不是第一层土的自重应力呢？

岩子228

我现在编的命令流如下，请大家多指点！
;
gen zone radcyl p0 (0,0,0) p1 (10,0,0) p2 (0,0,-29) p3 (0,10,0) &
                p4 (10,0,-29) p5 (0,10,-29) p6 (10,10,0) p7 (10,10,-29) &
                p8 (0.5,0,0) p9 (0,0.5,0) p10 (0.5,0,-29) p11 (0,0.5,-29) &
                size 5 10 29 15 ratio 1 1 1 1.15
gen zone radcyl p0 (0,0,-29) p1 (10,0,-29) p2 (0,0,-58) p3 (0,10,-29) &
                p4 (10,0,-58) p5 (0,10,-58) p6 (10,10,-29) p7 (10,10,-58) &
                p8 (0.5,0,-29) p9 (0,0.5,-29) p10 (0.5,0,-58) p11 (0,0.5,-58) &
                size 5 8 29 15 ratio 1 1 1 1.15 fill
;镜像，建立1/2模型，待定
gen zone reflect dd 270 dip 90
;土体网格分层
group tu2 range z -4 0
group tu3 range z -10 -4
group tu4 range z -14 -10
group tu5 range z -21 -14
group tu6 range z -24 -21
group tu7 range z -27 -24
group tu8 range z -34 -27
group tu9 range z -58 -34
;建立接触面
interface 1 face range cylinder end1 (0,0,0.01) end2 (0,0,-4) radius 0.51 &
                       cylinder end1 (0,0,0.01) end2 (0,0,-4) radius 0.49 not
interface 1 face range cylinder end1 (0,0,-4) end2 (0,0,-10) radius 0.51 &
                       cylinder end1 (0,0,-4) end2 (0,0,-10) radius 0.49 not
interface 1 face range cylinder end1 (0,0,-10) end2 (0,0,-14) radius 0.51 &
                       cylinder end1 (0,0,-10) end2 (0,0,-14) radius 0.49 not
interface 1 face range cylinder end1 (0,0,-14) end2 (0,0,-21) radius 0.51 &
                       cylinder end1 (0,0,-14) end2 (0,0,-21) radius 0.49 not
interface 1 face range cylinder end1 (0,0,-21 end2 (0,0,-24) radius 0.51 &
                       cylinder end1 (0,0,-21) end2 (0,0,-24) radius 0.49 not
interface 1 face range cylinder end1 (0,0,-24) end2 (0,0,-27) radius 0.51 &
                       cylinder end1 (0,0,-24) end2 (0,0,-27) radius 0.49 not
interface 1 face range cylinder end1 (0,0,-27) end2 (0,0,-29) radius 0.51 &
                       cylinder end1 (0,0,-27) end2 (0,0,-29) radius 0.49 not
interface 2 face range cylinder end1 (0,0,-28.9) end2 (0,0,-29.1) radius 0.51
;建立桩的模型（悬空）
gen zone cyl p0 (0,0,30) p1 (0.5,0,30) p2 (0,0,1) p3 (0,0.5,30) &
             p4 (0.5,0,1) p5 (0,0.5,1) &
             size 5 58 6
gen zone cyl p0 (0,0,30.1) p1 (0.5,0,30.1) p2 (0,0,30) p3 (0,0.5,30.1) &
             p4 (0.5,0,30) p5 (0,0.5,30) &
             size 5 1 6
gen zone reflect dd 270 dip 90 range z 1 30.1
group pile rang z 1 30.1
ini z add -30.0 range group pile  
; 定义本构关系及材料特性
model mohr range group tu2
model mohr range group tu3
model mohr range group tu4
model mohr range group tu5
model mohr range group tu6
model mohr range group tu7
model mohr range group tu8
model mohr range group tu9
prop bulk 1.333e7 shear 6.15e6 coh 12000 fric 22 range group tu2
prop bulk 1.842e7 shear 8.22e6 coh 28000 fric 5 range group tu3
prop bulk 1.5e7 shear 6.92e6 coh 22000 fric 10 range group tu4
prop bulk 1.583e7 shear 7.31e6 coh 20000 fric 15 range group tu5
prop bulk 9.33e6 shear 5.6e6 coh 8000 fric 36 range group tu6
prop bulk 1.75e7 shear 8.08e6 coh 20000 fric 15 range group tu7
prop bulk 2.16e7 shear 1.426e7 coh 6000 fric 38 range group tu8
prop bulk 8.333e7 shear 3.846e7 coh 30000 fric 0 range group tu9
model elas range group pile
prop bulk 1.75e10 shear 1.3125e10 range group pile
interface 1 prop kn 2.6e8 ks 2.6e8 coh 9600 fric 17.6 range z -4 0.01
interface 1 prop kn 2.183e8 ks 2.183e8 coh 22400 fric 4 range z -10 -4
interface 1 prop kn 1.800e8 ks 1.800e8 coh 17600 fric 8 range z -14 -10
interface 1 prop kn 1.900e8 ks 1.900e8 coh 16000 fric 12 range z -21 -14
interface 1 prop kn 1.248e8 ks 1.248e8 coh 6400 fric 28.8 range z -24 -21
interface 1 prop kn 2.100e8 ks 2.100e8 coh 16000 fric 12 range z -27 -24
interface 1 prop kn 3.014e8 ks 3.014e8 coh 4800 fric 30.4 range z -29 -27
interface 2 prop kn 3.014e8 ks 3.014e8 coh 4800 fric 30.4
     ;初始条件
ini dens 2000 range group tu2
ini dens 1890 range group tu3
ini dens 1950 range group tu4
ini dens 2050 range group tu5
ini dens 2060 range group tu6
ini dens 2030 range group tu7
ini dens 2060 range group tu8
ini dens 2080 range group tu9
ini dens 2500 range group pile
model null range z -0.1 0.15
  ;边界条件
fix z range z -58.1 -57.9
fix x range x -10.1 -9.9
fix x range x 9.9 10.1
fix y range y 9.9 10.1
set grav 0 0 -9.81
ini szz 0. grad 0 0 19620 range group tu2
ini szz 78480. grad 0 0 18540.9 range group tu3
ini szz 152643.6 grad 0 0 19129.5 range group tu4
ini szz 229161.6 grad 0 0 20110.5 range group tu5
ini szz 369935.1 grad 0 0 20208.6 range group tu6
ini szz 430560.9 grad 0 0 19914.3 range group tu7
ini szz 490303.8 grad 0 0 20208.6 range group tu8
ini szz 510512.4 grad 0 0 20404.4 range group tu9  最后的这个初始应力的参数赋值肯定有问题，但是现在还是没查出来，下个土层的初始应力算法，求教！

岩子228

现在的问题是，计算时报部分接触面节点刚度为零，不知道问题出在哪了。这个有影响吗？在自问自答中，期待高手的指点！