如题，flac3d如何提取不规则形状模型任意平面或曲面上应力？比如彭文斌老师FLAC3D实用教程141页的那个用傅里叶条件做模型上表面的模型，我想做一个z=50的水平切面，并想把这个水平面上的力提取出来，大家有没有什么好的办法？
解决方法一：我知道在flac3d中应力是存储在单元体重心的，其余地方的应力在应力云图表现出来时都是插值表现出来的，但是，咱们有没有那种大神会这种插值方法？
解决方法二：根据flac3d建模时单元体编号是有规律分布的，根据这种分布规律来知道某一层所有单元体编号范围，进而通过单元体编号获得某一层（可以是平面也可以是曲面，哈哈）单元体的应力。这个比较容易实现。
解决方法三：将flac3d的数据用郑海棠老师的Flac3D2Tecplot_new.dat导入tecplot，编个python的代码，利用tecplot中的探针函数提取某一个切面上的数据（应力），但这个难在写python的代码。
