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

tg0215

pl sel geom link on id on range ....

wang_509

感谢斑竹劳动了,学习了!...................

baiqunli

楼主真是厉害，讲的也精彩！
楼主多发这么经典的文章啊！

haitang

大师。写的很详细，赞一个。不过三维模型中，锚杆打多了，计算速度狂慢。我的一个模型打了两万根，机子算了四五天。

tg0215

原帖由 haitang 于 2007-1-10 23:45 发表
大师。写的很详细，赞一个。不过三维模型中，锚杆打多了，计算速度狂慢。我的一个模型打了两万根，机子算了四五天。

哈，海棠兄来捧场了，荣幸啊。是啊，我的一个模型不知道打了多数根锚杆（没数），用hp服务器算了4天还没算完，最后为了求速度，把一些不重要区域的锚杆删除了。

chunhuo420

sel cable id=1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 &
gr_coh 1 gr_k 1 gr_per 0.0785 range cid 1,10
sel cable id=1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 &
                    gr_coh 10e5 gr_k 2e7 range cid 11,17


楼主，你好，第一个表示自由段，第二表示锚固段，为什么第二个里面没有 gr_per 0.0785 这一项，gr_per是孔径大小吗？

chunhuo420

搂住写得太好了，
能不能把所写的东西打包一下，因为复制有好多乱码

姚旭东天

根椐tg0215班主的第二种情况~做的预应力锚索:
却出现如下错误提示:
source node 2 already has a link
为什么node 1没有错误?~到底错在哪里??
;=============加入锚索结构单元===============
sel cable id=1 begin=(32.5642,3,38) end=(46.7724,3,33.1911) nseg=7
sel cable id=1 begin=(46.7724,3,33.1911) end=(56.2446,3,29.9851) nseg=6
sel cable id=2 begin=(32.5642,3,37) end=(43.5851,3,32.2525) nseg=7
sel cable id=2 begin=(43.5851,3,32.2525) end=(52.7692,3,28.2963) nseg=6
sel cable id=3 begin=(32.5642,9,38) end=(46.7724,9,33.1911) nseg=7
sel cable id=3 begin=(46.7724,9,33.1911) end=(56.2446,9,29.9851) nseg=6
sel cable id=4 begin=(32.5642,9,37) end=(43.5851,9,32.2525) nseg=7
sel cable id=4 begin=(43.5851,9,32.2525) end=(52.7692,9,28.2963) nseg=6
;==================================================================
def ad
h=10000
loop ad1(1,4)
h=h+1
command
sel cable id ad1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 &
                  gr_coh 1 gr_k 1 gr_per 0.0785 range cid 1,7
sel cable id ad1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 &
                gr_coh 10e5 gr_k 2e7 range cid 8,13
sel delete link range id ad1
sel link id h ad1 target zone
sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id ad1
sel cable id ad1 pretension 4e5 range cid 1,7
endcommand
endloop
end
ad

tg0215

原帖由 chunhuo420 于 2007-1-11 15:14 发表
sel cable id=1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 &
gr_coh 1 gr_k 1 gr_per 0.0785 range cid 1,10
sel cable id=1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 &
             ...

gr_per是指孔径周长，也即锚固圈周长，我没有给出是因为我认为锚固剂的剪切强度跟围压没有关系，gr_per是用来计算sbond的。

tg0215

to 姚旭东天
请检查node，把该删除的node删除掉，source node 2 already has a link 的意思是node 2已经存在一个链接，必须要把已有的link删除掉才能建立新的node。用plot sel geom node on id on 查看node。

chunhuo420

原帖由 tg0215 于 2007-1-11 16:11 发表


gr_per是指孔径周长，也即锚固圈周长，我没有给出是因为我认为锚固剂的剪切强度跟围压没有关系，gr_per是用来计算sbond的。

这么说是砂浆的厚度没有影响吗？？
不好理解啊 。。。

benjackxu

原帖由 tg0215 于 2007-1-11 16:14 发表
to 姚旭东天
请检查node，把该删除的node删除掉，source node 2 already has a link 的意思是node 2已经存在一个链接，必须要把已有的link删除掉才能建立新的node。用plot sel geom node on id on 查看node。

鉴于tg0215版主孜孜不倦地、勤勤恳恳地为他的精华帖提供售后服务，特此加分！

chunhuo420

QUOTE:
原帖由 tg0215 于 2007-1-11 16:11 发表


gr_per是指孔径周长，也即锚固圈周长，我没有给出是因为我认为锚固剂的剪切强度跟围压没有关系，gr_per是用来计算sbond的。
这么说是砂浆的厚度没有影响吗？？
不好理解啊 。。。


再顶。。。。。。。。。。

chunhuo420

急于想知道，。。。。。。。。我顶。。。。。。。。。。

tg0215

gr_per是用来求sbond的，如果给出了gr_per，又给出了fric，则sbond会随有效围压的增大而增大。但是现在没有文献证明锚固剂的最大抗剪强度会随有效围压增大而增大（我没有看到过，如果你有这方面资料，希望能共享一下）。实际情况中，砂浆厚度肯定影响锚固力，数值模拟中通过sbond和kbond来描述砂浆的锚固力，并不表明砂浆厚度对锚固力没影响。

chunhuo420

谢谢回答。。。。。。。。。。。。。
kbond是啥东东，手册中没找到啊。

chunhuo420

原帖由 tg0215 于 2007-1-14 16:12 发表
gr_per是用来求sbond的，如果给出了gr_per，又给出了fric，则sbond会随有效围压的增大而增大。但是现在没有文献证明锚固剂的最大抗剪强度会随有效围压增大而增大（我没有看到过，如果你有这方面资料，希望能共享一 ...

根据M-C准则，抗剪强度是随围压的增大而增大的啊。。。。

chunhuo420

看到一个论文中的参数

如下图，算出锚杆的半径为50mm

钻孔的直径为100mm

好像有点不对吧，还是我的理解有误。。。？？

chunhuo420

:(:(:(:(:(:(:(:(:(:(:(:hug:

hanfengleng

to：姚旭东天
sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id ad1
这个地方出错了，id 后面的ad1 要改成 h