关于塑性区查找的问题

xingchen866

def zPlasticDepthInfo
  array bufff(1)
  p_z = zone_head
  loop while p_z # null
      z = 0
      n_zone=0
      mn = z_state(p_z,1)
    if mn = 1 then
        z1 = z_zcen(p_z)
      if z1 > z then
          z = z_zcent(p_z)
          n_zone = z_id(p_Z)  

      endif
   endif
   p_z=z_next(p_z)
endloop
bufff(1) = string(n_zone) +' '+string(z)
status = write(bufff,1)
end
def wrStress
  status = close
  status = open('zDepth.dat',1,1)
  if status = 0 then
     zPlasticDepthInfo
   status = close
  end_if
end
wrStress


这个段程序用于查找塑性单元的最大深度，写出这个单元的id和z坐标，但是为什么找不到呢？个人感觉问题出在Z_state（p_z，ind）中的ind取值或者mn的值上，不知道有没有高手指点一下啊，万分感激

faye100

1、ind=1好像只是代表shear_now吧。。。。
2、if z1 > z then....而你的z=0，问题是如果zcent（）都<0呢？

是不是应该写成if abs(z1) >= z

xingchen866

十分感谢提醒，z确实均小于0 1# xingchen866