一种常用的边坡支护形式--锚索框架的模拟命令流

li_2001 · 发表于 2009-11-13 11:39:27

问一下里面的节点是怎么找的

li_2001 · 发表于 2009-11-13 11:59:23

编写的是全锚索支护的情况一根锚索不需分段  现在我正在做一个有托盘的锚索  可是在运行过程中老是出现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 · 发表于 2009-11-13 11:59:40

编写的是全锚索支护的情况一根锚索不需分段  现在我正在做一个有托盘的锚索  可是在运行过程中老是出现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]

yshhui · 发表于 2010-3-10 13:28:08

还可以，命令流后有注释就好了

fanjunwei · 发表于 2010-3-10 22:00:15

期待完整的文件，不然没法学习！

在路上 · 发表于 2010-3-12 13:34:40

谢谢，学习学习

yshhui · 发表于 2010-3-17 20:42:46

很有借鉴意义谢谢

飞鱼游鹰 · 发表于 2010-3-17 22:07:36

楼主没发全啊!

hondu5 · 发表于 2010-3-20 21:29:58

期待，好好学习

haiyangzhou · 发表于 2010-3-26 15:03:43

先学习哈再说！！！

wangqihainl · 发表于 2010-8-3 12:23:32

谢谢。。。。我也是做这个的，，请以后多多指教

wangqihainl · 发表于 2010-8-4 10:36:26

期待楼主的完整命令流，万分感谢！

2007632001 · 发表于 2010-8-10 08:19:34

雪中送炭啊，最近也做锚索的模拟，可是锚索接触弄不明白，希望有用！

streamdream · 发表于 2011-3-5 01:42:44

wej101 · 发表于 2011-3-5 21:49:45

我下了，先看看

xinshe_cgs · 发表于 2011-4-5 20:35:19

支持分享！！！

wyj2005 · 发表于 2011-5-6 09:32:12

感谢楼主的奉献

ansys平安 · 发表于 2011-5-12 09:46:33

楼主你很有才锚索的直径133mm，犀利！！！

viking611 · 发表于 2011-5-16 22:09:12

不错，谢谢分享啦，呵呵

睡杭花梦 · 发表于 2011-5-23 00:43:43

很好，谢谢分享

账号		自动登录	找回密码
密码			注册

[工程实例] 一种常用的边坡支护形式--锚索框架的模拟命令流