请教钢拱架问题

xiaolongqdf001

各位同仁你们好，前几天看了galu的建立钢拱架的帖子，感觉非常好，但是由于水平有限，还处于初级水平阶段，有很多都是一知半解，自己也试着按照galu大哥的思路编了一段生成钢拱架的代码，可是就是不对，不知道错在什么地方了，期望能得到大家的帮助！在此谢谢各位了！
new
set log on
set logfile sdkw.log
gen zone radcyl p0 0 0 0 p1 20 0 0 p2 0 20 0 p3 0 0 20 &
p6 20 0 20 p8 10 0 0 p9 0 0 10 p10 10 20 0 p11 0 20 10 size 5 10 9 5 &
dimension 10 10 10 10
gen zone reflect dd 270 dip 90 origin 0 0 20
sel shell  range cyl end1 0 0 0 &
end2 0 20 0  rad 10 z 0  10
sel shell prop isotropic 10.5e9 0.25 thickness 0.02 density 2500
sel beam   beg 9.85 0 0 end 9.70035637 0 1.71043455
sel beam   beg  9.70035637 0 1.71043455  end 9.25597231 0 3.36889841
sel beam   beg  9.25597231 0 3.36889841 end 8.53035023 0 4.925
sel beam   beg  8.53035023 0 4.925 end 7.54553776 0 6.33145796
sel beam   beg  7.54553776 0 6.33145796 end 6.33145796 0 7.54553776
sel beam   beg  6.33145796 0 7.54553776 end 4.925 0 8.53035023
sel beam   beg 4.925 0 8.53035023 end 3.36889841 0 9.25597231
sel beam   beg  3.36889841 0 9.25597231 end 1.71043455 0 9.70035637
sel beam   beg  3.36889841 0 9.25597231 end 1.71043455 0 9.70035637
sel beam   beg  1.71043455 0 9.70035637 end 0 0 9.85
sel beam   beg   0 0 9.85 end  -1.71043455 0 9.70035637
sel beam   beg   -1.71043455 0 9.70035637 end -3.36889841 0 9.25597231
sel beam   beg   -3.36889841 0 9.25597231 end -4.925 0 8.53035023
sel beam   beg   -4.925 0 8.53035023 end  -6.33145796 0 7.54553776
sel beam   beg    -6.33145796 0 7.54553776 end -7.54553776 0 6.33145796
sel beam   beg  -7.54553776 0 6.33145796 end -8.53035023 0 4.925
sel beam   beg  -8.53035023 0 4.925 end -9.25597231 0 3.36889841
sel beam   beg  -9.25597231 0 3.36889841 end -9.70035637 0 1.71043455
sel beam   beg   -9.70035637 0 1.71043455 end -9.85 0 0
def parmater
   y_=0
   p=5
   Rbeam=9.85000
   Rshell=9.98  
    a=0
    det=10
end
parmater
def operate_link
       link_id=1000
        y=y_
       loop m(1,p)
          loop n(1,18)
             x_shellnode=Rshell*cos((a+n*det)*degrad)
             z_shellnode=Rshell*sin((a+n*det)*degrad)
             y_shellnode=y
             x_beamnode=Rbeam*cos((a+n*det)*degrad)
             z_beamnode=Rbeam*sin((a+n*det)*degrad)
             y_beamnode=y
             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
         y=y+0.8
      end_loop
end
operate_link

cloud328

你的问题我给你找出来了，首先你多定义了一段钢拱架，总共只有18段，你定义了19段，其次，你的竖直方向是y方向，而你的命令流里面是z方向，我把方向换了下，再次，我用fish批量定义了5柄钢拱架。建议新手，先认认真真看看书，把基本命令搞清楚。你的这个问题其实很简单。

xiaolongqdf001

非常希望各位版主和高手能够帮忙解决解决！

xiaolongqdf001

求求各位大侠帮帮我吧！时间紧迫啊…………

cloud328

说出你觉得问题所在啊。我看前面都没有问题，是不是后面的link fish有问题？

xiaolongqdf001

是啊！我运行的时候，一直出现target node and source node are the same，不知道什么原因？请指点指点。感谢！

xiaolongqdf001

版主说的没错，我真的应该看看书，一些基本问题都没搞明白。真是谢谢你了！

xxh345

我按这个加上参数和约束算了下，为什么加不加梁计算结果都一样呢？

yushibo1985429

8# xxh345

msh2010

版主的命令流还有和楼主的问题一样啊！！一直出现target node and source node are the same

benjackxu

奇了怪了，用cloud328 的命令运行没问题，一改参数就有问题。