请教：塑性单元的提取及重新定义性质？？？

53105189

请教高手一个问题：以下是我编写的一段fish程序，目的是想将塑性区中处于tensionnow状态下的塑性单元提取，并将它们定义为'tensionyieldnow'组，但试过之后并不能达到目的，望高手给与指点！！！不胜感激！！！

def showplasticzones
shearnow =  1
tensionnow = 2
shearpast = 4
tensionpast = 8

zp=zone_head
loop while zp # null
if  z_state(zp , 0) = 2
z_group(zp) = 'tensionyieldnow'
end_if
zp=z_next(zp)
end_loop
end
showplasticzones

haiqsim

来向楼主学习并讨论下，
我不大会用fish，不过偶有个想法想与楼主交流下。

您说的“将塑性区中处于tensionnow状态下的塑性单元提取，并将它们定义为'tensionyieldnow'组”，就是想将处于tensionnow状态下的塑性单元提取出来，并定义为一个名为tensionyieldnow的group吧？

而一个group中包含多个单元体，应该可以理解为一个group里包含多个单元体的ID号吧？

于是，您的想法是否可以这样实现：
遍历所有单元，找到处于tensionnow状态下的单元后，将这个单元赋予一个全新的ID号，这个ID号必须是模型中没有的，比如模型所有单元ID号为1~15000，而你的新赋予的ID好就可以从15001甚至20001开始。每找到一个符合要求的单元，就重新赋一个ID号，(可以设置这个新的ID号在循环后自动加1)，等到最后循环完毕后，采用命令“group tensionnow range id_start_new id_end_new”，这样应该就可以实现了吧？

我只是有这个想法，具体的函数不大熟，不知楼主以为是否可行？仅供参考！