【求助】高手，斑竹请帮助，急急急！！！

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  
我认为先生成还是后生成，都应该是一样的，为什么会有差别？应该怎么样做呢

dynamax

肯定是不一样的
用 ini z add 的话
节点位置只是重合
实际上有两个节点
你可以查一下
gen merge 1.0e-5 可以合并重合节点
这样就一样了

chichilv

哦
这样时说我如果直接用gen zone成和的节点就覆盖了，而用移动的话就是两个独立的节点！

chichilv

但是要用到interface,这样的话必须要用ini z add 了，是吧，关于这个我是参考example 8 的例子做的，但是做出来的初始应力不收敛，肯定是我的问题，但我不知道问题在哪里，困扰我好多天了，请手指教，下面是我的命令
  
;;--------------------------------------------------------------------------
; gen geometry
;--------------------------------------------------------------------------
;gen radbrick
;zuo
gen zone brick p0 16 0 0 p1 80 0 0 p2 16 0 -18 p3 16 9 0 p4 80 0 -80 p5 16 9 -18 &
               p6 80 45 0 p7 80 45 -80 size 20 18 9 rat 1.1 1 1   
gen zone brick p0 16 9 0 p1 80 45 0 p2 16 9 -18 p3 16 10 0 p4 80 45 -80 p5 16 10 -18 &
               p6 80 50 0 p7 80 50 -80 size 20 18 2 rat 1.1 1 1
;you
gen zone brick p0 0 10 0 p1 15 10 0 p2 0 10 -18 p3 0 50 0 p4 15 10 -18 p5 0 50 -80 &
               p6 75 50 0 p7 75 50 -80 size 15 18 20 rat 1 1 1.1
gen zone brick p0 15 10 0 p1 16 10 0 p2 15 10 -18 p3 75 50 0 p4 16 10 -18 p5 75 50 -80 &
               p6 80 50 0 p7 80 50 -80 size 2 18 20 rat 1 1 1.1
;xia
gen zone brick p0 0 0 -18 p1 15 0 -18 p2 0 0 -80 p3 0 9 -18 p4 75 0 -80 p5 0 45 -80 &
               p6 15 9 -18 p7 75 45 -80 size 15 20 9 rat 1 1.1 1
gen zone brick p0 15 0 -18 p1 16 0 -18 p2 75 0 -80 p3 15 9 -18 p4 80 0 -80 p5 75 45 -80 &
               p6 16 9 -18 p7 80 45 -80 size 2 20 9 rat 1 1.1 1  
gen zone brick p0 0 9 -18 p1 15 9 -18 p2 0 45 -80 p3 0 10 -18 p4 75 45 -80 p5 0 50 -80 &
               p6 15 10 -18 p7 75 50 -80 size 15 20 2 rat 1 1.1 1  
gen zone brick p0 15 9 -18 p1 16 9 -18 p2 75 45 -80 p3 15 10 -18 p4 80 45 -80 p5 75 50 -80 &
               p6 16 10 -18 p7 80 50 -80 size 2 20 2 rat 1 1.1 1
;def interface
interface 1 face range x 15.9 16.1 y -0.1 10.1 z -18.1 0.1  
interface 2 face range x -0.1 16.1 y 9.9 10.1 z -18.1 0.1
;top  
gen zone brick p0 0 0 30 p1 15 0 30 p2 0 0 12 p3 0 9 30  size 15 18 9
gen zone brick p0 15 0 30 p1 16 0 30 p2 15 0 12 p3 15 9 30 size 2 18 9
gen zone brick p0 0 9 30 p1 15 9 30 p2 0 9 12 p3 0 10 30 size 15 18 2
gen zone brick p0 15 9 30 p1 16 9 30 p2 15 9 12 p3 15 10 30 size 2 18 2
group soil
  
ini z add -30 range z 12 30
;wall and cuttion
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                                       range group soil
pro bulk=6e7  shear=2e7  f=0 coh=30000    range group soil
mod elas                                      range group wall
pro bulk=6e7  shear=2e7    range group wall
;
interface 1 pro  kn 1e8  ks 1e8 fric 20 coh 30000
interface 2 pro  kn 1e8  ks 1e8 fric 20 coh 30000
;
ini dens 1600 range group soil
ini dens 1600 range group wall
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
save inisoil.sav
;---------ini stress-wall
mod elas                                     range group wall
pro bulk=13.9e9  shear=10.4e9                range group wall
ini dens 2400 range group wall
hist unbal
solve
save iniwall.sav
;--------------excavate--------
ini xdis 0 ydis 0 zdis 0
group excavate range x 0 15.5 y 0 9.5 z -9 0  
ini dens 1600 range group excavate
;
mod mohr                                range group excavate
pro bulk=6e6  shear=2e6  f=20  coh=1e4  range group excavate  
;excavate
mod null range  group excavate
set large
hist unbal
solve
这个会在2.*e-4左右，不在收敛，为什么？
烦劳多谢指教

chichilv

这是网格图

chichilv

这是interface

dynamax

gen merge 1.0e-5 ran z -18.1 17.9

chichilv

多谢多谢，dynamax  
我下午发现了这个问题了，现在可以算通了。谢谢，但是结果算出来还不太对
等我弄弄结果，请你指教

cybersky

请问chichilv 和dynamax ,下面的语句中
ini szz 0. grad 0 0 12300. range z -5.5 0.0
ini szz 17600 grad 0 0 15500 range z -8 -5.5
ini sxx 0.0 grad 0 0 5271.4 range z -5.5 0.0
ini sxx 7542.86 grad 0 0 6642.86 range z -8 -5.5
ini sxx add 31428.6 grad 0 0 5714.3 range z -8 -5.5
ini syy 0. grad 0 0 5271.4 range z -5.5 0.0
ini syy 7542.86 grad 0 0 6642.86 range z -8 -5.5
ini syy add 31428.6 grad 0 0 5714.3 range z -8 -5.5
除了12300和15500是水位分界面处上下土的重度外,其他数字是怎么得来的?
比如第三句ini sxx 0.0 grad 0 0 5271.4 range z -5.5 0.0中的5271.4是怎么得来的?

cybersky

怎么没人理我?(z)