如何提取地层(group)单元体的应力值？

yjw_413

通过对单元遍历循环，可以提取对应的z_szz(p_z)，但是我的模型是分不同地层的，我想分别提取各不同地层对应单元的垂直应力，该如何做呢？
以下是我对整体单元体提取应力值的程序。因为分开提取是为了对不同地层的垂直应力乘以不同的数值！不知我描述的清楚不？请问有会的吗？帮助一下，自己没找到合适的资料！

1. 
   
2. 
   
3. def ini_stress ;存储初始地应力
   
4. 
   
5. p_z=zone_head
   
6. 
   
7. loop while p_z#null
   
8. n=z_id(p_z)
   
9. xtable(1,n)=z_szz(p_z)
   
10. p_z=z_next(p_z)
    
11. 
    
12. endloop
    
13. end

复制代码

yjw_413

我自己解决了通过z_group命令查询每个单元所属的几何状态，再赋不同的属性！

wanlang

define a_moment
p_z=zone_head
loop while p_z # null
if z_group(p_z)='pile'
n=z_id(p_z)
xtable(1,n)=z_szz(p_z)
x1=xtable(1,n)
ytable(1,n)=z_xcen(p_z)
y1=ytable(1,n)
mmz=x1*y1
xtable(2,n)=mmz
ytable(2,n)=n
endif
p_z=z_next(p_z)
endloop
end
@a_moment
set logfile 221.txt
set log on
print table 2

zyt2013

询问下如何获取指点坐标点的应力值，谢谢！