版权归作者所有，转载请注明出处
Flac3D的云图插值方法总结，contour
- 一、Flac3D的视图根据值的来源大致可以分为几种：
1：源自zone（其实是离散四面体，一个6面体离散5个四面体）
- 源自zone的值较多，比如zone-pp on的孔隙水压力，渗流场的渗透率，比如本构模型的力学参数（体积模量，剪切模量，内摩擦角、杨氏模量，泊松比等），加上密度参数density，基本上都是源自zone的，比如模型的strain，比如模型的屈服状态state等，都是源自zone，具体只要看说明书的部分。
2、源自节点Gridpoint的云图：
- 源自节点的值非常多，比如速度，饱和度、位移，force力等，节点在Flac3D内部是很重要的一个物件，承载了质量和运动的，间接影响到模型的变形，其中阻尼施加的地方就是节点，我们在做静力学的时候其实就是上图的一些迭代。
3、源自其他物件，比如结构单元，接触面等。
- 二、Flac3D云图的添加（以zone视图为例）
其中有部分的云图contour可以可选来源是节点还是zone的，比如渗流场的孔隙水压力，其余的基本上是确定好的。当然，视图里面没有的云图也可以通过extra额外增量进行自定义添加（31集附加），实现自定义云图插值。
- 三、Flac3D云图的插值方法：
大致有以上图片四种插值方法，一般情况我们做应力或者位移云图的时候，建议使用最后一种（外差法），这样子上下界面的值和我们算出来的是一致的，当然你的模型网格足够多的，几种插值方法的差别会较小。
- 另外一种云图方法是根据Tension的XYX和Magnitude来的，其实很容易理解，就是xyz方向或者开根号的xyz的平方和。Magnitude就是一个值，不存在方向性。四、友情提醒：6.0的云图里面的zone unbalance force云图是显示有错误的，但是用fish或者history调用zone或者节点的最大不平衡力是正确的，仅仅是显示问题，不影响计算和其他结果。
fish define cxj
f=zone.unbal()
end
@cxj
flac3d>[f]
2.212336e-01 (real)
其实真正的zone 不平衡力是0.22，而不是500的。显示有问题，完全可以用extra额外增量来自己做个图。下图是7.0运行的结果：
- 本内容为07集附加的部分内容节选，更多内容和源代码文件请看07集附加上下集。为了方便大家查阅做的文档，希望对大家有帮助。转账或者引用请注明出处，本文为Flac3D系列教程内容，谢谢大家。有错误希望大家指出。
