【讨论】ini z add 的作用，对模型的影响

chichilv

在2.0手册EXAMPLE APPLICATION 8  Axial and Lateral Loading of a Concrete Pile 里用到了"ini z add -6.0 range group pile"
但是我再用的过程中遇到了问题
在用到这个时在生成初始应力场时就不平衡，不知道试我用的不对还是？但我不用这个命令在想要的地方直接生成网格后计算有没了个这个问题，不明白是为什么；下面是我做的简单的命令
1）用到了 ini z add  
;gen radbrick
  gen zone radbrick  p0 0 0 0 p1 80 0 0 p2 0 0 -80  p3 0 50 0 p8 16 0 0 p9 0 0 -18 p10 0 10 0 &
     size 32 18 20 15 rat 1 1 1 1.08
  gen zone brick p0 0 0 30 p1 16 0 30 p2 0 0 12 p3 0 10 30  size 32 18 20
  ini z add -30 range x 0 16 y 0 10 z 12 30  
  group soil
  group wal1 range x 15.5 16 y 0 10 z -18 0
  group wal2 range x 0 15.5 y 9.5 10 z -18 0  
  group wall range  group wal1 any group wal2 any
mod mohr                                       
pro bulk=6e6  shear=2e6  f=20 coh=1e9  ten=1e9   
ini dens 1600  
ini szz 0.0 grad 0 0 1.6e4 range z -80 0
;--------------boun----------------
fix x     range x -0.1 0.1
fix x     range x 79.9 80.1  
fix y     range y -0.1 0.1
fix y     range y 49.9 50.1
fix x y z range z -80.1  -79.9
;
set g 0 0 -10
set small
hist unbal
solve
2)不用 ini z add  
gen zone radbrick  p0 0 0 0 p1 20 0 0 p2 0 0 -30  p3 0 20 0 p8 8 0 0 p9 0 0 -17 p10 0 8 0 &
     size 8 17 8 8 rat 1 1 1 1.15  
gen zone brick p0 0 0 0 p1 8 0 0 p2 0 0 -17 p3 0 8 0 size 8 17 8  
group soil
plot add surface yellow
group excavate range x 0 5 y 0 5 z -12 0  
group wall1  range x 5 6 y 0 6 z -14 0
group wall2  range x 0 5 y 5 6 z -14 0
group wall range  group wall1 any group wall2 any
model mohr
pro bulk=5e8   shear=6e8     fric=20
pro  coh=8e4  tens=8e4  
set gravity 0,0,-10
ini  density 1600 range group soil
ini density 2500 range group wall  
fix x range x -0.1  0.1
fix x range x 19.9 20.1
fix y range y -0.1 0.1
fix y range y 19.9 20.1
fix x y z range  z -29.9 -30.1
set g 0 0 -10
ini szz 0 grad 0 0 1.6e4 range z -30 0
hist unbal
set small
solve

chichilv

ruti