受邀谈论FLAC3D预应力锚杆模拟心得。

gnik77

1# tg0215 精品！顶

天外飞仙人

楼主牛人呀，

deyin007

我现在试一试看看好不好试

啥都懂的兔子

楼主 很佩服你！
请教个问题：做基坑支护的时候，上部是土钉，下部是桩加锚索，为什么计算稳定系数的时候出现双排锚稳定系数反而比单排锚稳定系数小呢？第一排锚的位置长度等都一样，只是在下部增加一排锚。第二排锚的参数与第一排锚相同。

tangjincai28

不顶都说不过去了，向tg0215学习

keepfit

不错，看看

zuoruo69

顶，顶，这东西 太好了

888438

那我想问，要是模拟基坑桩锚支护，该是怎么一个模拟过程那？如何让编写命令流那？

xunxunyuanyuan

强帖留名标记一下
很有用~~~

xunxunyuanyuan

tg0215，不知道你有没有模拟过压力型或者压力分散型锚索呢？
由于受力机理不同，cable单元赋予参数只能模拟一般的拉力型锚索，所以很想请教一下，你有什么方法可以较好的模拟压力型的么？

liuzhonghua

花了半天终于仔细看完了,真是获益匪浅啊.在次感谢那些牛人们.

li_2001

如果有建模部分的命令流就好了,没做过这方面的,俺想系统学一下

li_2001

高手 您好 您编写的是全锚索支护的情况 一根锚索不需分段  现在我正在做一个有托盘的锚索  可是在运行过程中老是出现soure node not near enough to an exsiting non-null zone 现在将程序传上来望您指点 谢谢

new
gen zone brick size 10 10 10 &
p0 0 0 0 p1 10 0 0 p3 0 0 10
model mohr
prop bulk 1e8 shear 0.3e8 fric 35 coh 1e10 ten 1e10
set grav 0 0 -9.81
ini dens 1000
fix x range x 0
fix x range x 10
fix y range y 0
fix y range y 10
fix z range z 0
hist unbal
set hist_rep 10
solve
save 初始平衡.sav

ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
model null range x 0 2 y 0 2 z 0 2
solve
save 1次开挖.sav

restore 初始平衡.sav
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
model null range x 0 2 y 0 2 z 0 2
def place_cables
    loop m(1,2)
    caseof m
    case 1
    x_d1=0.5
    x_d2=0.5
    x_d3=0.5
    y_d1=0.5
    y_d2=0.5
    y_d3=0.5
    z_d1=1.5
    z_d2=2.0
    z_d3=3.5
       a=1
       b=2
       c=3
       d=6
       id_11=1
       id_111=1
    case 2
    x_d1=0.5
    x_d2=0.5
    x_d3=0.5
    y_d1=1.5
    y_d2=1.5
    y_d3=1.5
    z_d1=1.5
    z_d2=2.0
    z_d3=3.5
       a=15
       b=16
       c=17
       d=20
       id_11=3
       id_111=15
   endcase
   loop n(1,2)
   x_d1=x_d1+1
   x_d2=x_d2+1
   x_d3=x_d3+1
   y_d1=y_d1
   y_d2=y_d2
   y_d3=y_d3
   z_d1=z_d1
   z_d2=z_d2
   z_d3=z_d3
       a=a+(n-1)*7
       b=b+(n-1)*7
       c=c+(n-1)*7
       d=d+(n-1)*7
        id_11=id_11+1
       id_111=id_111+(n-1)*7
  command
  sel cable id=id_11 beg x_d1 y_d1 z_d1 end x_d2 y_d2 z_d2 nseg 2
  sel cable id=id_11 beg x_d2 y_d2 z_d2 end x_d3 y_d3 z_d3 nseg 4
  sel cable id=id_11 prop emod 2e10 ytension 3.1e5 xcarea 0.0004906 &
   gr_coh 1 gr_k 1 gr_per 0.0785 range cid a,b
  sel cable id= id_11 prop emod 2e10 ytension 3.1e5 xcarea 0.0004906 &
   gr_coh 10e5 gr_k 2e7 gr_per 0.0785 range cid c,d
  sel delete link range id id_111
  sel link id=id_111 id_111 target zone
  sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id id_111
  sel cable id=id_11 pretension 60e3 range cid a,d
  endcommand
  end_loop
  end_loop
  end
  place_cables

li_2001

[高手 您好 您编写的是全锚索支护的情况 一根锚索不需分段  现在我正在做一个有托盘的锚索  可是在运行过程中老是出现soure node not near enough to an exsiting non-null zone 现在将程序传上来望您指点 谢谢

new
gen zone brick size 10 10 10 &
p0 0 0 0 p1 10 0 0 p3 0 0 10
model mohr
prop bulk 1e8 shear 0.3e8 fric 35 coh 1e10 ten 1e10
set grav 0 0 -9.81
ini dens 1000
fix x range x 0
fix x range x 10
fix y range y 0
fix y range y 10
fix z range z 0
hist unbal
set hist_rep 10
solve
save 初始平衡.sav

ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
model null range x 0 2 y 0 2 z 0 2
solve
save 1次开挖.sav

restore 初始平衡.sav
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
model null range x 0 2 y 0 2 z 0 2
def place_cables
    loop m(1,2)
    caseof m
    case 1
    x_d1=0.5
    x_d2=0.5
    x_d3=0.5
    y_d1=0.5
    y_d2=0.5
    y_d3=0.5
    z_d1=1.5
    z_d2=2.0
    z_d3=3.5
       a=1
       b=2
       c=3
       d=6
       id_11=1
       id_111=1
    case 2
    x_d1=0.5
    x_d2=0.5
    x_d3=0.5
    y_d1=1.5
    y_d2=1.5
    y_d3=1.5
    z_d1=1.5
    z_d2=2.0
    z_d3=3.5
       a=15
       b=16
       c=17
       d=20
       id_11=3
       id_111=15
   endcase
   loop n(1,2)
   x_d1=x_d1+1
   x_d2=x_d2+1
   x_d3=x_d3+1
   y_d1=y_d1
   y_d2=y_d2
   y_d3=y_d3
   z_d1=z_d1
   z_d2=z_d2
   z_d3=z_d3
       a=a+(n-1)*7
       b=b+(n-1)*7
       c=c+(n-1)*7
       d=d+(n-1)*7
        id_11=id_11+1
       id_111=id_111+(n-1)*7
  command
  sel cable id=id_11 beg x_d1 y_d1 z_d1 end x_d2 y_d2 z_d2 nseg 2
  sel cable id=id_11 beg x_d2 y_d2 z_d2 end x_d3 y_d3 z_d3 nseg 4
  sel cable id=id_11 prop emod 2e10 ytension 3.1e5 xcarea 0.0004906 &
   gr_coh 1 gr_k 1 gr_per 0.0785 range cid a,b
  sel cable id= id_11 prop emod 2e10 ytension 3.1e5 xcarea 0.0004906 &
   gr_coh 10e5 gr_k 2e7 gr_per 0.0785 range cid c,d
  sel delete link range id id_111
  sel link id=id_111 id_111 target zone
  sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id id_111
  sel cable id=id_11 pretension 60e3 range cid a,d
  endcommand
  end_loop
  end_loop
  end
  place_cables]

li_2001

请tg0215 答复 万分感谢

mengqb

好资料 谢谢分享，受益匪浅。

mengqb

好资料 谢谢分享，受益匪浅。

kjj0501

brick里面的dim指的是什么
手册上好像没有啊 180# hh199422

gaosongsong

精彩！希望楼主能举出一些锚索方面的例子。

liuhewu

4# tg0215

sel link id=100 1 target node tgt_num 18  
上面命令行中最后的"18“是如何定的