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。