flac3d群桩模拟求助 [复制链接]

刚刚学习flac3d不久，模拟一个群桩，想计算桩顶沉降量、看看桩周围应力图。下面是我在一个算例的基础上改的，算着有问题，我也弄不太明白。恳请有经验得朋友，帮帮忙，看看哪里有问题，怎么修改才能实现我的模拟，十分感谢。ps：命令流中一些土体参数都是借鉴相关工程的，我之前学工民建，对于岩土没什么经验，也不知道设置的对不对。
以下是命令流：
n
;桩周土模型建立
gen zone radcyl p0 (2.,2.,0)   P1(3.,2.,0)    p2 (2.,2.,-15)    p3(2.,3.,0)    &    
                p4 (3.,2.,-15) p5 (2.,3.,-15) p6 (3.,3.,0)      p7 (3.,3.,-15) &
                p8 (2.25,2.,0) p9 (2.,2.25,0) p10 (2.25,2.,-15) p11 (2.,2.25,-15) &
                size 1 20 4 2 ratio 1 1 1 1   group Soil_Pc
plo blo gro                  
gen zone reflect dip 90 dd 90   origin = (2.,2.,0) range group Soil_Pc
gen zone reflect normal (0,1,0) origin = (2.,2.,0) range group Soil_Pc
;
gen zone radcyl p0 (2.,0,0)   P1(3.,0,0)     p2 (2.,0,-15)    p3 (2.,1.,0)      &    
                p4 (3.,0,-15) p5 (2.,1.,-15) p6 (3.,1.,0)     p7 (3.,1.,-15)    &
                p8 (2.25,0,0) p9 (2.,0.25,0) p10 (2.25,0,-15) p11 (2.,0.25,-15) &      
                size 1 20 4 2 ratio 1 1 1 1                                  
gen zone radcyl p0 (2.,0,0)    P1 (2.,1.,0)  p2 (2.,0,-15)     p3  (1.,0,0)     &    
                p4 (2.,1.,-15) p5 (1.,0,-15) p6 (1.,1.,0)      p7 (1.,1.,-15)   &  
                p8 (2.,0.25,0) p9 (1.75,0,0) p10 (2.,0.25,-15) p11 (1.75,0,-15) &  
                size 1 20 4 2 ratio 1 1 1 1      
gen zone radcyl p0 (0,2.,0)    p1 (1.,2.,0)  p2 (0,2.,-15)     p3 (0,3.,0) p4 (1.,2.,-15) &
                p5 (0,3.,-15)  p6 (1.,3.,0)  p7 (1.,3.,-15)  &
                p8 (0.25,2.,0) p9 (0,2.25,0) p10 (0.25,2.,-15) p11 (0,2.25,-15)&
                size 1 20 4 2 ratio 1 1 1 1  
gen zone radcyl p0 (0,2.,0)    p1 (0,1.,0)    p2 (0,2.,-15)    p3 (1.,2.,0) p4 (0,1.,-15) &
                p5 (1.,2.,-15) p6 (1.,1.,0)   p7 (1.,1.,-15) &
                p8 (0,1.75,0)  p9 (0.25,2.,0) p10 (0,1.75,-15) p11 (0.25,2.,-15) &
                size 1 20 4 2 ratio 1 1 1 1                
gen zone radcyl p0 (0,0,0)    p1 (1.,0,0)   p2 (0,0,-15)     p3 (0,1.,0) p4 (1.,0,-15) &
                p5 (0,1.,-15) p6 (1.,1.,0)  p7 (1.,1.,-15) &
                p8 (0.25,0,0) p9 (0,0.25,0) p10 (0.25,0,-15) p11 (0,0.25,-15) &
                size 1 20 4 2 ratio 1 1 1 1  
;桩下土模型建立
gen zone radcyl p0 (2.,2.,-15)   P1(3.,2.,-15)    p2 (2.,2.,-30)    p3(2.,3.,-15)    &          
                p4 (3.,2.,-30)   p5 (2.,3.,-30)   p6 (3.,3.,-15)    p7 (3.,3.,-30) &
                p8 (2.25,2.,-15) p9 (2.,2.25,-15) p10 (2.25,2.,-30) p11 (2.,2.25,-30) &
                size 1 5 4 2 ratio 1 1 1 1   group Base_Pc                    
gen zone cyl p0 (2.,2.,-15)   P1(2.25,2.,-15) p2 (2.,2.,-30) p3(2.,2.25,-15)  &
             p4 (2.25,2.,-30) p5 (2.,2.25,-30) &
             size 1 5 4  group Base_Pc                                  
gen zone reflect dip 90 dd 90   origin = (2.,2.,-15) range group Base_Pc
gen zone reflect normal (0,1,0) origin = (2.,2.,-15) range group Base_Pc
;
gen zone radcyl p0 (2.,0,-15)   P1(3.,0,-15)     p2 (2.,0,-30)    p3 (2.,1.,-15)      &          
                p4 (3.,0,-30)   p5 (2.,1.,-30)   p6 (3.,1.,-15)   p7 (3.,1.,-30)    &
                p8 (2.25,0,-15) p9 (2.,0.25,-15) p10 (2.25,0,-30) p11 (2.,0.25,-30) &      
                size 1 5 4 2 ratio 1 1 1 1           fill                      
gen zone radcyl p0 (2.,0,-15)    P1 (2.,1.,-15)  p2 (2.,0,-30)     p3  (1.,0,-15)     &                          
                p4 (2.,1.,-30)   p5 (1.,0,-30)   p6 (1.,1.,-15)    p7 (1.,1.,-30)   &            
                p8 (2.,0.25,-15) p9 (1.75,0,-15) p10 (2.,0.25,-30) p11 (1.75,0,-30) &  
                size 1 5 4 2 ratio 1 1 1 1     fill
gen zone radcyl p0 (0,2.,-15)    p1 (1.,2.,-15)  p2 (0,2.,-30)     p3 (0,3.,-15) p4 (1.,2.,-30) &
                p5 (0,3.,-30)    p6 (1.,3.,-15)  p7 (1.,3.,-30)  &
                p8 (0.25,2.,-15) p9 (0,2.25,-15) p10 (0.25,2.,-30) p11 (0,2.25,-30)&
                size 1 5 4 2 ratio 1 1 1 1     fill
gen zone radcyl p0 (0,2.,-15)    p1 (0,1.,-15)    p2 (0,2.,-30)    p3 (1.,2.,-15) p4 (0,1.,-30) &
                p5 (1.,2.,-30)   p6 (1.,1.,-15)   p7 (1.,1.,-30) &
                p8 (0,1.75,-15)  p9 (0.25,2.,-15) p10 (0,1.75,-30) p11 (0.25,2.,-30) &
                size 1 5 4 2 ratio 1 1 1 1      fill        
gen zone radcyl p0 (0,0,-15)    p1 (1.,0,-15)   p2 (0,0,-30)     p3 (0,1.,-15) p4 (1.,0,-30) &
                p5 (0,1.,-30)   p6 (1.,1.,-15)  p7 (1.,1.,-30) &
                p8 (0.25,0,-15) p9 (0,0.25,-15) p10 (0.25,0,-30) p11 (0,0.25,-30) &
                size 1 5 4 2 ratio 1 1 1 1      fill
;边界内其余土模型建立
gen zone radtunnel  p0 (0,0,0)     P1(15,0,0)     p2 (0,0,-15)   p3(0,15,0)     &
                    p4 (15,0,-15)  p5 (0,15,-15)  p6 (15,15,0)   p7 (15,15,-15) &
                    p8 (3.,0,0)    p9 (0,3.,0)    p10 (3.,0,-15) p11 (0,3.,-15) p12 (3.,3.,0) p13 (3.,3.,-15) &
                    size 6 20 6 15 ratio 1 1 1 1.1
gen zone radtunnel  p0 (0,0,-15)  P1(15,0,-15)  p2 (0,0,-30)   p3(0,15,-15)   &
                    p4 (15,0,-30) p5 (0,15,-30) p6 (15,15,-15) p7 (15,15,-30) &
                    p8 (3.,0,-15) p9 (0,3.,-15) p10 (3.,0,-30) p11 (0,3.,-30) p12 (3.,3.,-15) p13 (3.,3.,-30) &
                    size 6 5 6 15 ratio 1 1 1 1.1
;
attach face range z = -14.9 -15.1
;设置不同土层
group soft_soil ran z 0 -8.6
group sand ran z -8.6 -15
group clay ran z -16 -30
;接触面建立
interface 1 face range cylin end1 (2.,2.,0)     end2 (2.,2.,-15.1) radius .26 &
                       cylin end1 (2.,2.,0)     end2 (2.,2.,-15.1) radius .24 not        
interface 2 face range cylin end1 (2.,2.,-14.9) end2 (2.,2.,-15.1) radius .26
interface 3 face range cylin end1 (2.,0,0)      end2 (2.,0,-15.1) radius .26 &
                       cylin end1 (2.,0,0)      end2 (2.,0,-15.1) radius .24 not    
interface 4 face range cylin end1 (2.,0,-14.9)  end2 (2.,0,-15.1) radius .26
interface 5 face range cylin end1 (0,2.,0)      end2 (0,2.,-15.1) radius .26 &
                       cylin end1 (0,2.,0)      end2 (0,2.,-15.1) radius .24 not      
interface 6 face range cylin end1 (0,2.,-14.9)  end2 (0,2.,-15.1) radius .26
interface 7 face range cylin end1 (0,0,0)       end2 (0,0,-15.1) radius .26 &
                       cylin end1 (0,0,0)       end2 (0,0,-15.1) radius .24 not    
interface 8 face range cylin end1 (0,0,-14.9)   end2 (0,0,-15.1) radius .26
;桩模型建立
;角桩
gen zone cyl p0 (2.,2.,16)  P1(2.25,2.,16) p2 (2.,2.,1) p3(2.,2.25,16)  &
             p4 (2.25,2.,1) p5 (2.,2.25,1) &
             size 1 20 4  group Pc
gen zone reflect dip 90 dd 90   origin = (2.,2.,16) range group Pc
gen zone reflect normal (0,1,0) origin = (2.,2.,16) range group Pc
;边桩
gen zone cyl p0 (2.,0,16) P1(2.25,0,16) p2 (2.,0,1) p3(2.,0.25,16)  &
             p4 (2.25,0,1) p5 (2.,0.25,1) &
             size 1 20 4  
gen zone cyl p0 (2.,0,16) P1(2.,0.25,16) p2 (2.,0,1) p3(1.75,0,16)  &
             p4 (2.,0.25,1) p5 (1.75,0,1) &
             size 1 20 4  
gen zone cyl p0 (0,2.,16) P1(0.25,2.,16) p2 (0,2.,1) p3(0,2.25,16)  &
             p4 (0.25,2.,1) p5 (0,2.25,1) &
             size 1 20 4  
gen zone cyl p0 (0,2.,16) P1(0,1.75,16) p2 (0,2.,1) p3(0.25,2.,16)  &
             p4 (0,1.75,1) p5 (0.25,2.,1) &
             size 1 20 4                
;中心桩
gen zone cyl p0 (0,0,16) P1(.25,0,16) p2 (0,0,1) p3(0,.25,16)  &
             p4 (.25,0,1) p5 (0,.25,1) &
             size 1 20 4
group pile range z 1 16.1
ini z add -16.0 range group pile
;模型参数设置
model mohr range group soft_soil
prop bulk 5.556e6 shear 1.852e6 coh 15000 fric 13 range group soft_soil
model mohr range group sand
prop bulk 1.250e7 shear 5.769e6 coh 20000 fric 16 range group sand
model mohr range group clay
prop bulk 5.00e7 shear 2.308e7 coh 30000 fric 22 range group clay
model elas range group pile
prop bulk 5.556e6 shear 1.852e6 range group pile
;接触面参数
interface 1 prop kn 1e10 ks 1e10 fric 16.7 coh 0
interface 2 prop kn 1e10 ks 1e10 fric 16.7 coh 0
interface 3 prop kn 1e10 ks 1e10 fric 16.7 coh 0
interface 4 prop kn 1e10 ks 1e10 fric 16.7 coh 0
interface 5 prop kn 1e10 ks 1e10 fric 16.7 coh 0
interface 6 prop kn 1e10 ks 1e10 fric 16.7 coh 0
interface 7 prop kn 1e10 ks 1e10 fric 16.7 coh 0
interface 8 prop kn 1e10 ks 1e10 fric 16.7 coh 0
;密度设置
ini dens 2000 range group clay
ini dens 1800 range group sand
ini dens 1800 range group soft_soil
ini dens 1800 range group pile
;边界条件
fix z range z -30.1 -29.9
fix x range x -.1 .1
fix x range x 14.9 15.1
fix y range y -.1 .1
fix y range y 14.9 15.1
;施加重力作用
set grav 0 0 -10
;计算
plot con zd
solve
save PG_0.sav
;桩参数重新赋值
model elas range group pile
prop bulk 6.667e9 shear 5.0e9 range group pile
ini dens 2500 range group pile
plo cont zdisp
solve
save PG_1.sav
;位移清零
ini state 0
ini xdis 0 ydis 0 zdis 0
set logfile 1
set log on
;
rest PG_1.sav
;加载
apply szz -3.2e5 range z -0.1 0.1 group pile
hist gp zdis 0,0,0
plot con zd
;计算
solve
save Q.sav

给你参考一个例子

;一号桩外土体模型
gen zone radcyl p0 (0,0,0) p1 (.75,0,0) p2 (0,0,-22) p3 (0,.75,0) &
p4 (.75,0,-22) p5 (0,.75,-22) p6 (.75,.75,0) p7 (.75,.75,-22) &
p8 (.25,0,0) p9 (0,.25,0) p10 (.25,0,-22) p11 (0,.25,-22) &
size 2 22 6 3 ratio 1 1 1 1.25
gen zone radcyl p0 (0,0,-22) p1 (.75,0,-22) p2 (0,0,-45) p3 (0,.75,-22) &
p4 (.75,0,-45) p5 (0,.75,-45) p6 (.75,.75,-22) p7 (.75,.75,-45) &
p8 (.25,0,-22) p9 (0,.25,-22) p10 (.25,0,-45) p11 (0,.25,-45) &
size 2 11 6 3 ratio 1 1 1 1.25 fill
;二号桩外土体模型
gen zone copy 1.5 0 0
gen zone reflect normal 1 0 0 ori 1.5 0 0 ran x 1.5 2.25
;三号桩外土体模型
gen zone copy 0 1.5 0 ran x 0 0.75
gen zone reflect normal 0 1 0 ori 0 1.5 0 ran y 1.5 2.25
;四号桩外土体模型
gen zone copy 0 1.5 0 ran x 0.75 2.25
gen zone reflect normal 0 1 0 ori 0 1.5 0 ran x 0.75 2.25 y 1.5 2.25
;承台外土体模型
gen zone radtunnel p0 (0,0,0) p1 (6,0,0) p2 (0,0,-45) p3 (0,6,0) &
p4 (6,0,-45) p5 (0,6,-45) p6 (6,6,0) p7 (6,6,-45) &
p8 (2.25,0,0) p9 (0,2.25,0) p10 (2.25,0,-45) p11(0,2.25,-45) &
p12 (2.25,2.25,0) p13 (2.25,2.25,-45) &
size 9 45 9 5 ratio 1 1 1 1.3
group soil1 range z 0 -1.03
group soil2 range z -1.03 -2.63
group soil3 range z -2.63 -7.61
group soil4 range z -7.61 -13.29
group soil5 range z -13.29 -17.69
group soil6 range z -17.69 -33.55
group soil7 range z -33.55 -45
model mohr range group soil1
model mohr range group soil2
model mohr range group soil3
model mohr range group soil4
model mohr range group soil5
model mohr range group soil6
model mohr range group soil7
prop bulk 4.0e6  shear 1.25e6 fric 11.4 cohesion 13e3 ten 1e10 range group soil1
ini density 1870 range group soil1
prop bulk 6.67e6 shear 3.08e6 fric 10.8 cohesion 10e3 ten 1e10 range group soil2
ini density 1890 range group soil2
prop bulk 6.67e6 shear 3.33e6 fric 9.8  cohesion 16e3 ten 1e10 range group soil3
ini density 1920  range group soil3
prop bulk 13.33e6 shear 4.44e6 fric 12.5 cohesion 14e3 ten 1e10 range group soil4
ini density 1990 range group soil4
prop bulk 10e6 shear 6e6 fric 13.1 cohesion 12e3 ten 1e10 range group soil5
ini density 1970 range group soil5
prop bulk 16.67e6 shear 7.69e6 fric 16 cohesion 30e3 ten 1e10 range group soil6
ini density 2020  range group soil6
prop bulk 16.67e6 shear 7.69e6 fric 20.8 cohesion 43e3 ten 1e10 range group soil7
ini density 2040 range group soil7
;1号桩接触面
interface 1 face range cylinder end1(0,0,0.1) end2(0,0,-22.1) radius .26 &
cylinder end1(0,0,0.1) end2(0,0,-22.1) radius .24 not
interface 2 face range cylinder end1(0,0,-21.9) end2(0,0,-22.1) radius .26
;2号桩接触面
interface 1 face range cylinder end1(1.5,0,0.1) end2(1.5,0,-22.1) radius .26 &
cylinder end1(1.5,0,0.1) end2(1.5,0,-22.1) radius .24 not
interface 2 face range cylinder end1(1.5,0,-21.9) end2(1.5,0,-22.1) radius .26
;三号桩接触面
interface 1 face range cylinder end1(0,1.5,0.1) end2(0,1.5,-22.1) radius .26 &
cylinder end1(0,1.5,0.1) end2(0,1.5,-22.1) radius .24 not
interface 2 face range cylinder end1(0,1.5,-21.9) end2(0,1.5,-22.1) radius .26
;四号桩接触面
interface 1 face range cylinder end1(1.5,1.5,0.1) end2(1.5,1.5,-22.1) radius .26 &
cylinder end1(1.5,1.5,0.1) end2(1.5,1.5,-22.1) radius .24 not
interface 2 face range cylinder end1(1.5,1.5,-21.9) end2(1.5,1.5,-22.1) radius .26
interface 1 prop kn=2.7e7 ks=2.7e7 fric=8 coh=43e3
interface 2 prop kn=2.7e7 ks=2.7e7 fric=8 coh=43e3
;建立桩体模型
;一号桩
gen zone cyl p0 (0,0,23) p1 (.25,0,23) p2 (0,0,1) p3 (0,.25,23) &
p4 (0.25,0,1) p5 (0,0.25,1) &
size 1 22 6
;二号桩
gen zone copy 1.5 0 0 ran z 1 23
gen zone reflect normal 1 0 0 ori 1.5 0 0 ran x 1.5 2.25 z 1 23
;三号桩
gen zone copy 0 1.5 0 ran x 0 0.75 z 1 23
gen zone reflect normal 0 1 0 ori 0 1.5 0 ran y 1.5 2.25 z 1 23
;4号桩
gen zone copy 0 1.5 0 ran x 0.75 2.25 z 1 23
gen zone reflect normal 0 1 0 ori 0 1.5 0 ran x 0.75 2.25 y 1.5 2.25 z 1 23
group zhuang range z 1 23
ini z add -23.0 range group zhuang
model elas range group zhuang
prop bulk 15.56e6  shear 11.67e6
ini density 2500 range group zhuang
;建立承台模型
gen zone brick p0(0,0,0) p1(2.25,0,0) p2(0,0,-1) p3(0,2.25,0) size 1,1,1
group chengtai range x 0 2.25 y 0 2.25 z 0 -1
model elas range group chengtai
prop bulk 15.56e6  shear 11.67e6
ini density 2500 range group chengtai
gen merge 1e-4 ran cylinder end1 (0,0,-.01) end2 (0,0,.01) radius 0.16
gen merge 1e-4 ran cylinder end1 (1.5,0,-.01) end2 (1.5,0,.01) radius 0.16
gen merge 1e-4 ran cylinder end1 (0,1.5,-.01) end2 (0,1.5,.01) radius 0.16
gen merge 1e-4 ran cylinder end1 (1.5,1.5,-.01) end2 (1.5,1.5,.01) radius 0.16
fix z range z -45.1 -44.9
fix x range x -.1 .1
fix x range x 5.9 6.1
fix y range y -.1 .1
fix y range y 5.9 6.1
set grav 0 0 -10
ini szz -19261 grad 0 0 19261 range z 0 -1.03
ini szz -49801 grad 0 0 30240 range z -1.03 -2.63
ini szz -145417 grad 0 0 95616 range z -2.63 -7.61
ini szz -258449 grad 0 0 113032 range z -7.61 -13.29
ini szz -345129 grad 0 0 86680 range z -13.29 -17.69
ini szz -665501 grad 0 0 320372 range z -17.69 -33.55
ini szz -899081 grad 0 0 233580 range z -33.55 -45
apply szz -48e3 range z .05 .15
solve
plot gp disp ran id 0 0 0