[求助]flac3D怎么用fish实现多个结点力相加 [复制链接]

哪位兄弟对fish比较了解啊，讨教个问题
我想用loop循环，把x从0-15，y从0-15，z是50。16*16个结点力相加.

def get_ad
        i=1       
        loop while i<=15
        j=1
                loop while j<=15
                ad(i,j)=gp_near(i,j,50)  ;这句有问题,是不是gp_near()里不能包含变量?
                j=j+1
                end loop
        i=i+1
        end loop
end
get_ad
def load1
        load=0
                loop i(1,15)
                                loop j(1,15)                                 
                load=load+gp_zfunbal(ad(i,j)) ;这里感觉可能也有问题,但不知道咋改,请高手指点
                end loop 
        end loop
end
问题简单,还请高手不吝赐教.我现在实在没招了.翻手册找不到地方.

我把fish 改了下,现在编译通过了,但是运行几步发现load始终为0.快崩溃了,还以为问题已经解决了.下附完整命令.
gen zone brick size 5 5 5 & p0 0 0 0 p1 5 0 0 p2 0 5 0 p3 0 0 5
model mohr
prop bulk 2e+7 shea 1.200e+7  den 0 fri 10 coh 9.800e+4 ten 9.8e3
fix x y z range z  -.1 .1
apply zvel -0.001e-3 range z 4.9 5.1 x 0 5 y 0 5
;定义数组
def arry1
        array ad(6,6)
end
arry1
;把坐标存入数组
def get_ad
        i=1
        loop while i<=6
        j=1
                loop while j<=6
                m=i-1
                n=j-1
                ad(i,j)=gp_near(m,n,5)
                j=j+1
                end_loop
        i=i+1
        end_loop
end
get_ad
;计算各个结点力
def load1
        load=0
        loop i(1,6)
                loop j(1,6)
                load=load+gp_zfunbal(ad(i,j))
                end_loop
        end_loop
end
load1
his load
期待高手指点.

问题已解决 

你好，我想问一下，结点力与等效结点力有什么区别，能否告诉你是怎么实现的，谢谢Q272891858