ParaView---FLAC3D和3DEC输出变量到VTK文件

1 引言

在《免费二维有限元分析软件ADONIS计算步骤》中曾经提及ADONIS的计算结果可以输出为VTK文件，事实上3DEC和FLAC3D的计算结果也可以输出为VTK文件。这个笔记简要讨论了输出步骤。

2 VTK命令

3DEC和FLAC3D使用block vtk和zone vtk命令输出vtk文件，可以输出的变量有：

(1) displacements: 网格点的位移矢量

(2) zone-extra: 单元索引i的额外变量

(3) gridpoint-extra: 网格点索引i的额外变量

(4) forces: 网格点不平衡力的信息，包括最后的局部应力比【FLAC/Slope的强度降低过程(Strength Reduction Procedure in FLAC/Slope)】

(5) groups: slot中单元设置的组

(6) model-fluid: 流体属性

(7) model-mechanical: 力学属性

(8) model-thermal: 热属性

(9) pore-pressure: 网格点的孔隙压力值

(10) range: 设置range, 用来控制哪些单元包括在vtk文件中

(11) saturation: 网格点饱和度的值

(11) stresses: 每个单元的平均应力。不输出内部四面体的应力。

(12) temperatures: 网格点的温度值

(13) velocities: 网格点的速度矢量

(14) filename: 设置vtk的文件名

默认情况下，只输出应力(stresses)，组(groups)和位移(displacements). vtk的文件格式是XML语言。

block vtk filename "slope.vtk"

3 VTK显示

下图所示的边坡共输出1491的网格点(points)和3184的单元(cells), 左边是3DEC计算的位移等值线图，右边是tvk产生的位移等值线图。

4 ParaView简介

VTK代表着Visualization ToolKit(VTK)，vtk文件由ParaView读取。ParaView是一个开源的用于可视化二维和三维数据集的程序。ParaView主要是为了处理具有空间代表性数据而设计的。这个程序与商业性软件Tecplot的功能类似([New] Tecplot Focus 2021 R1---数值模拟后处理软件)，但ParaView是开源的，完全使用Python语言编写。ParaView的开发得到了桑迪亚国家实验室(Sandia National Laboratories)，洛斯阿拉莫斯国家实验室(Los Alamos National Laboratories)和美国国家科学基金会(US National Science Foundation)以及其它一些工业项目的资助，目前的最新版本是5.10.0-RC1。