你的问题是你建的模型很大,包含好几个采场,计算完毕后目前你只想提取其中一个采场的塑性区体积,可以这么理解吗?这个采场在一定空间里,在模型里有固定的位置,可以这么写:
p_z = zone_head
loop while p_z # null
*********
p_z = z_next(p_z)
假设这个采场的x坐标区间为(1,3),y坐标区间为(2,3),z坐标区间为(3,4),,上面的单元体地址 p_z 每循环一次做一个判断,判断这个单元体地址是不是在这个采场之内,若在这个采场内如果还满足塑性区条件则该单元体体积有效并进行累加,反之无效继续寻找,一直到null,累加起来的塑性区体积则为单个采场塑性区体积。fish程序编写原理就是这样的,详细程序我编写试试看。