VTUinterface

VTUinterface 是一个用于方便访问由有限元软件（如OpenGeoSys）输出的 VTU/PVD 文件的 Python 包。它使用 VTK Python 包装器和时间步长与网格点之间的线性插值来访问模拟域内任何空间和时间点。

VTUinterface 与 ogs6py 一起可以在此处查看操作

0. 安装

注意：VTUinterface 需要Python的vtk包装器。或者，基于 MESHIO 的 版本 也正在开发中。克隆仓库并使用pip安装包

# pip install [--user] https://github.com/joergbuchwald/VTUinterface/archive/refs/heads/master.zip

2. VTUinterface 文档

您可以在 https://joergbuchwald.github.io/VTUinterface-doc/ 找到文档

3. 快速开始

注意：方法命名风格已更改（2021-05-20）

基本用法（python）

尽管 VTUinterface 是一个 Python 包，但它通过 PyCall 在 Julia 中也可以测试运行

基本用法（julia）

可以通过以下方式运行单元测试：

# python tests/test_vtuinterface.py

从项目根目录。

4. 常见问题/故障排除

安装

• 如果无法在 PyPI 上找到 vtk 轮子，您可以查看 https://github.com/pyvista/pyvista/discussions/2064 以获取非官方轮子。

由于输入数据使用 QHull 进行了三角剖分，线性插值可能会在边界处失败或如果给出错误输入维度。可能的解决方案

• 为了正确执行插值，提供问题正确的维度（通过 dim 关键字设置）至关重要。
• 由于dim关键字还指定了要使用的坐标，VTUinterface假定dim=1指的是x坐标，而默认情况下dim=2表示问题出在xy平面上。这可以通过指定一维的one_d_axis或二维的two_d_planenormal来改变。
• 对于某些网格，调整三角剖分考虑的点数可能会有所帮助，这可以通过使用nneighbors关键字来完成。默认值为20。
• 特别是在边界附近，使用QHULL方法的线性插值通常会失败，这可以通过使用最近邻插值来解决。
• 或者，您现在可以将interpolation_backend从scipy更改为vtk，并尝试不同的插值核。