急求帮助！新手上路，一段关于双排桩支护的，结果显示不可以cycle，恳请各位点拨

liangqiumin

诚恳希望得到各位的点拨，以达对flac结构化的认识，为感！！双排桩，中间以连梁（beam单元模拟）相连
new
plot sk

;建立土体模型
gen zone brick size 30 11 30
group clay
mod mohr range group clay
pro bulk 5e6 shear 1.6e6 cohesion 15000 &
    friction 8 dilation 0. tension 0. range group clay


;设置边界条件和初始应力
fix x range x -.1 .1
fix x range x 29.9 30.1
fix y range y -.1 .1
fix y range y 10.9 11.1
fix z range z -0.1 0.1

ini szz -369000 grad 0 0 12300 range z 0 30
ini sxx -158130 grad 0 0 5271  range z 0 30
ini syy -158130 grad 0 0 5271  range z 0 30


;设置结构单元
sel        pile id 1 begin 10 0 30 end 10 0 10 nseg 17
sel        pile id 1 begin 11 0 30 end        11 0 10 nseg 17
sel        pile id 1 begin 10 2.5 30 end 10 2.5 10 nseg 17
sel        pile id 1 begin 11 2.5 30 end        11 2.5 10 nseg 17
sel        pile id 1 begin 10 5 30 end 10 5 10 nseg 17
sel        pile id 1 begin 11 5 30 end        11 5 10 nseg 17
sel        pile id 1 begin 10 7.5 30 end 10 7.5 10 nseg 17
sel        pile id 1 begin 11 7.5 30 end        11 7.5 10 nseg 17
sel        pile id 1 begin 10 10 30 end 10 10 10 nseg 17
sel        pile id 1 begin 11 10 30 end 11 10 10 nseg 17

sel nod id 501 10.5 0 30
sel nod id 502 10.5 2.5 30
sel nod id 503 10.5 5 30
sel nod id 504 10.5 7.5 30
sel nod id 505 10.5 10 30

sel nod id 601 10 1.25 30
sel nod id 602 10 3.75 30
sel nod id 603 10 6.25 30
sel nod id 604 10 8.75 30

sel nod id 701 11 1.25 30
sel nod id 702 11 3.75 30
sel nod id 703 11 6.25 30
sel nod id 704 11 8.75 30

def set_beam
   loop ii (1,5)
       iid1=1+36*(ii-1);节点号
       iid2=500+ii        
       iid3=19+36*(ii-1)
       command
       sel  beamsel  id 2  nodes iid1  iid2      ;刚性梁
       sel  beamsel  id 2  nodes iid2  iid3
       endcommand
   endloop
end
set_beam

def qianliang
   loop ii (1,4)
       iid1=1+36*(ii-1) ;节点号
       iid2=600+ii
       iid3=37+36*(ii-1);节点号
       command
       sel  beamsel  id 3  nodes iid1   iid2      ;前梁
       sel  beamsel  id 3  nodes iid2   iid3
       endcommand
   endloop
end
qianliang

def houliang
   loop ii (1,4)
       iid1=19+36*(ii-1);节点号
       iid2=700+ii
       iid3=55+36*(ii-1);节点号
       command
       sel  beamsel  id 3  nodes iid1   iid2      ;后梁
       sel  beamsel  id 3  nodes iid2   iid3
       endcommand
   endloop
end
houliang


sel pile id = 1 prop emod 3.0e10 nu 0.16 xcarea 0.785 &  ;半径设置不合理？？？
                xciy 0.049  xciz 0.049   xcj 0.092  &
                per 3.14 &
                cs_sk  0.11e9 cs_scoh 1e4 cs_sfric 13.46 &
                cs_nk  0.67e9 cs_ncoh 0.0 cs_nfric 0.0 cs_ngap off

sel beam id = 2 prop emod 3.0e10 nu 0.16 xcarea 0.785 &
               xciy 0.049  xciz 0.049   xcj 0.092
      

sel beam id = 3 prop emod 2.2e11 nu 0.13 xcarea 0.196 &
                 xciy 2e-4  xciz 2e-4   xcj 4e-4 &
                 ydir 1 0 0  

ini dens 2000 range group clay
set gravity 0 0 -10.0

hist unbal
solve ratio 1.e-6
solve

;开挖模拟
ini xdis 0 ydis 0 zdis 0

model null range x 0 10 y 0 11 z 20 30
step 100000     
save liang.sav

[ 本帖最后由 liangqiumin 于 2007-2-16 20:10 编辑 ]

liangqiumin

自己先顶起来！！！！！！！！！！！！！！！

friendpp

难题，需要高手解决呀。