重新链接结构单元

echowasd

; Flac3D 5.0
; 自动链接结构单元，不多解释，用过结构单元应该明白
; w.gang.89@foxmail.com
def relink(begin_link_id)
    dist_tol = 1e-2
    link_id = begin_link_id
    node_num = 0
    node_pnt0 = nd_head
    ; m_pnt = null
    loop while node_pnt0 # null
        ; node_num = node_num + 1
        ; m_pnt1 = get_mem(2)
        ; mem(m_pnt1) = m_pnt
        ; mem(m_pnt1+1) = node_pnt0
        ; m_pnt = m_pnt1
        ; node_pnt0 = nd_next(node_pnt0)
        node_num = node_num + 1
        node_pnt0 = nd_next(node_pnt0)
    endloop

    node_arr = get_array(node_num)
    node_num = 0
    node_pnt0 = nd_head
    loop while node_pnt0 # null
        node_num = node_num + 1
        node_arr(node_num) = node_pnt0
        node_pnt0 = nd_next(node_pnt0)
    endloop

    node_num_minus1 = node_num - 1
    loop ii(1,node_num_minus1)
        node_pnt1 = node_arr(ii)
        node_id1 = nd_id(node_pnt1)
        xxa = nd_pos(node_pnt1,2,1)
        yya = nd_pos(node_pnt1,2,2)
        zza = nd_pos(node_pnt1,2,3)
        ii_plus1 = ii + 1
        loop jj(ii_plus1,node_num)
            node_pnt2 = node_arr(jj)
            node_id2 = nd_id(node_pnt2)
            xxb = nd_pos(node_pnt2,2,1)
            yyb = nd_pos(node_pnt2,2,2)
            zzb = nd_pos(node_pnt2,2,3)
            node_dist = sqrt((xxa-xxb)^2+(yya-yyb)^2+(zza-zzb)^2)
            if node_dist <= dist_tol then
                link_pnt1 = nd_link(node_pnt1)
                link_pnt2 = nd_link(node_pnt2)
                if link_pnt1 = null then
                    link_id = link_id + 1
                    command
                        ; sel set link node_tol=dist_tol
                        sel link id=@link_id @node_id1 target=node tgt_num=@node_id2
                        sel link attach xdir=rigid ydir=rigid zdir=rigid &
                            xrdir=rigid yrdir=rigid zrdir=rigid range id=@link_id
                    endcommand
                else
                    if link_pnt2 = null then
                        link_id = link_id + 1
                        command
                            ; sel set link node_tol=dist_tol
                            sel link id=@link_id @node_id2 target=node tgt_num=@node_id1
                            sel link attach xdir=rigid ydir=rigid zdir=rigid &
                                xrdir=rigid yrdir=rigid zrdir=rigid range id=@link_id
                        endcommand
                    else
                        if lk_type(link_pnt1) = lk_type(link_pnt2) then
                        ; value {1,2} denotes node-to-node or node-to-zone link
                            temp = lk_delete(link_pnt2)
                            link_id = link_id + 1
                            command
                                sel link id=@link_id @node_id2 target = node tgt_num=@node_id1
                                sel link attach xdir=rigid ydir=rigid zdir=rigid &
                                    xrdir=rigid yrdir=rigid zrdir=rigid range id=@link_id
                            endcommand
                        endif
                    endif
                endif
            endif
        endloop
    endloop
    temp = lose_array(node_arr)
end
@relink(10000)

echowasd

涛声依旧2015 发表于 2015-8-6 13:14
有在flac3d3.0中自动创建结构单元链接的例子么，我最近在做基坑支护的过程中想在桩上建立锚杆，桩与锚杆之 ...

我这个是5.0版本的，基于前人版本补充修改的。
5.0下用法示例：

;@install_shell
sel node local xdir=(1,0,0) ydir=(0,1,0) range x 24 24.1
sel node fix lsys range x 24 24.1
sel node fix x yr zr range x 24 24.1
; 开挖西部第一层 z1=40.5， z2=42.5
@excav(40.5,42.5)
@install_cable_w(1)
@relink(30000)
solve
save quarter_excav_w1; end

不需要特殊的设置，安装好结构单元后，直接relink就可以了，程序自动查找同一位置节点并重新链接。
w.gang.89@foxmail.com这个邮箱可以直接加QQ，西南交通大学土木工程学院地下工程系。

leige135264

flac3d 3.0 能用不  试试看

涛声依旧2015

有在flac3d3.0中自动创建结构单元链接的例子么，我最近在做基坑支护的过程中想在桩上建立锚杆，桩与锚杆之间的链接都是通过查看链接号用笨方法建立的，程序写的很多，感觉自动建立程序挺简洁的，我想你肯定会，盼望你能在百忙之中提供个例子，谢谢。邮箱1650791281@qq.com

龙一乜

好好哦啊后好好好

herculis

这个我使用的时候总是说目标点和原点不够近，怎么回事啊？

落魄皇

thank you master

3Dgeosciences

感谢分享指导，受益多多

boliwutoubang

感谢分享指导

ys610710656

谢谢楼主，好好学习学习