跟大家请教并探讨输出某个单元应变的问题

hndfzhanghui

new
;建立模型
gen zone brick p0=(0,0,0) p1=(5,0,0) p2=(0,5,0) p3=(0,0,5) size 5 5 5
;模型参数
model mohr
def look_d
  p_z=zone_head
  loop while p_z # null
       z_prop(p_z,'bulk')=urand*1e8
       z_prop(p_z,'shear')=urand*0.3e8
       z_prop(p_z,'friction')=urand*35
       z_prop(p_z,'cohesion')=urand*1e10
       z_prop(p_z,'tension')=urand*1e10
       p_z=z_next(p_z)
  end_loop
end
look_d
;边界条件
fix z range z -0.1 0.1
;初始条件
set gravity 0,0,-9.81
initial density=1000
apply szz=-2e7 range z 4.9 5.1
;解
solve
;寻找应力应变最大块
def find_i
   array ar(6)
   p_z=zone_head
   loop while p_z # null
   dum=z_fsi(p_z,ar)
   max=ar(3)
   li=z_szz(p_z)
   smax=max*li
       if tt<smax then
       tt=smax
       zid=z_id(p_z）
       end if
    p_z=z_next(p_z)
    end loop
    ii=out('单元体zid号=’string(zid))
    kk=out('应力乘应变最大='tt)
end
find_i  
这事源代码，一下附文本文档，但当程序运行至
;寻找应力应变最大块
def find_i
   array ar(6)
   p_z=zone_head
   loop while p_z # null
   dum=z_fsi(p_z,ar)
   max=ar(3)
   li=z_szz(p_z)
   smax=max*li
处时会出现如图问题

当我输入pr fish后命令栏显示

我想请教各位牛人，这是什么情况
该如何解决
我程序的目的是赋给每个单元格随即的属性
然后寻找到z轴方向应力乘应变最大的那个单元格
很感谢