大家都来交流一下隧道初期支护中的“钢拱架”应该如何建立

galu

我先开个头吧，说说自己的体会：
1、钢拱架一般都采用多个beamsel小段来模拟（调整形状），同一榀钢拱架只要保证组成它的beamsel都属同一id就可以了。这一部分好解决。
2、设置每个node与网格的关系。我个人认为：实际的钢拱架都与围岩密贴接触（加垫块调整），这和shell模拟的喷混凝土与围岩的接触关系有异曲同工之处，区别仅在于网格不能带动node（钢拱架）旋转。因此node与zone之间的Link就应该只保留x,y,z，而将link中的xrot,yrot,zrot删除。

以上是我对钢拱架如何建立的一些浅显想法，有什么不足之处还请各位版主和大牛们提出中肯意见；以起到抛砖引玉的作用，大家共同交流经验，一起进步。

galu

补充一下：
3、如果钢拱架安装之前，已经有shell（喷混凝土）附在网格上。那么在生成beamsel时，两个端点node可以与shell的node设在同一位置（两个node，分属于shell和beam），然后删除node（beam）与zone的Link，最后建立beam  node与shell  node 的联结（仅是x,y,z,不包括xrot,yrot,zrot），这保证了beam和shell的实际联结情况——密贴但不会随之转动。

wylbest

楼主的帖子真好啊，应该建议加分阿！在下也要模拟钢拱架，无奈对楼主的““那么在生成beamsel时，两个端点node可以与shell的node设在同一位置（两个node，分属于shell和beam），然后删除node（beam）与zone的Link，最后建立beam  node与shell  node 的联结（仅是x,y,z,不包括xrot,yrot,zrot）””这句话不太了解，渴望得到楼主指点！

galu

你看看下面的程序吧，应该会对你有所启发。


按照这个思路我试着建立了模型，效果还是不错的.我把link的fish处理程序（一个处理Link的小思路）贴出来吧，大家分享一下，可能有的地方不是很合理，欢迎斧正。



restore 111.sav

def parmater
   ;起始位置
   z_=0
   ;榀数（纵向）
   p=5
   ;半径
   Rbeam=5.85000
   Rshell=6.0084(洞室半径)
   ;首个node的角度位置
    a=0
    det(a)=11.0589
end
parmater

;beam半径稍小于shell（洞室）半径

def operate_link
       link_id=1000
        z=z_
       loop m(1,p)
          loop n(0,13)
             x_shellnode=Rshell*cos((a+n*det(a))*degrad)
             y_shellnode=Rshell*sin((a+n*det(a))*degrad)
             z_shellnode=z

             x_beamnode=Rbeam*cos((a+n*det(a))*degrad)
             y_beamnode=Rbeam*sin((a+n*det(a))*degrad)
             z_beamnode=z

             sn_pointer=nd_near(x_shellnode,y_shellnode,z_shellnode)
             bn_pointer=nd_near(x_beamnode,y_beamnode,z_beamnode)
            
             shellnode_id=nd_id(sn_pointer)
             beamnode_id=nd_id(bn_pointer)

             command
                  sel set link node_tol 0.3
                  sel link    id   link_id    beamnode_id     target      node      tgt_num      shellnode_id
                  sel link attach    xdir=rigid    ydir=rigid    zdir=rigid   range     id    link_id
                  sel link attach    xrdir=free    yrdir=free    zrdir=free   range    id    link_id
             end_command

             link_id=link_id+1
         end_loop
         z=z+0.8
      end_loop
end
operate_link

wylbest

一看就知道有高手风范，不但自己会又分享给后来者，楼主，怎以句感谢了得阿，我下去一定好好体会一下你的程序！

wolf-4989

早看到就好，呵呵

lakewater

加分鼓励！

xhb

galu :
    你好！你提出的建立shell  与beem建立连接整体受力的方案很好。
但你帖出的“link的fish处理程序”中没有删除原来的shellnode或beemnode与zoom之间的连接（link ），而直接创建beamnode与shellnode之间的连接能行吗？一个节点已存在某种连接时不能再和其它节点创建连接吧。

cqsunhui

先下来学习下。

zqutyut

galu :
你好! 一看就是高手,可我们这些初学者很难看懂这些东西,虽然说明了意图,我也查了手册的相应的命令,但是资质有限,还是搞不懂,比如:x_shellnode=Rshell*cos((a+n*det(a))*degrad);  sel link    id   link_id    beamnode_id     target      node      tgt_num      shellnode_id;sel link attach    xdir=rigid    ydir=rigid    zdir=rigid   range     id    link_id 等是什么意思?

zqutyut

sel set link node_tol 0.3有什么作用?

galu

sel set link node_tol 0.3 表示只要两个node之间距离小于0.3，就产生link

x_shellnode=Rshell*cos((a+n*det(a))*degrad)  这是给出和beam的node产生link的node（shell）坐标值

sel link    id   link_id    beamnode_id     target      node      tgt_num      shellnode_id 是在将beam的node 和shell的node  link起来。手册中有这个命令行的解释。搜索Link章节即可查得。

sel link attach    xdir=rigid    ydir=rigid    zdir=rigid   range     id    link_id  是设置link的初始条件。手册中亦有这个命令行的解释。搜索Link章节即可查得。

galu

原帖由 xhb 于 2006-6-12 18:54 发表
galu :
    你好！你提出的建立shell  与beem建立连接整体受力的方案很好。
但你帖出的“link的fish处理程序”中没有删除原来的shellnode或beemnode与zoom之间的连接（link ），而直接创建beamnode与shellnode之 ...

        产生beamnode时，我控制了beam产生的半径略小于洞室半径，这样beam不会和zone产生Linke。另外，shellnode与zone之间的link不用删除。因为对于一个source node(shellnode)来说，它只能拥有一个属于它自己的Link（shellnode和zone之间的Link，对于这个Link来说，shellnode是source node），也就是说这个shellnode不能再去"link"其他的zone或者node。但是这个shellnode却可以被多个其他的zone或node去“link”，在上面的程序中，我设置beamnode去"link" shellnode。

teapotd

好人呀

lbt18872

先学习学习再发言吧

zqutyut

好人!

lianghlq

不知楼主对2D中的钢拱架的模拟有何高见，请教中

galu

2D中不是有很多beam单元用来模拟衬砌的实例么，把那些例子中的beam属性改成钢拱架的参数就可以了。

freddie_1999

是好帖子。
论坛就需要这样的步骤：问题的提出－讨论－最终的解决。
体现了提高的过程。

sch

好贴！考虑很细！请教：你设置beam半径稍小于shell（洞室）半径，他们厚度如何考虑？还请解释清楚！