Python脚本提取Abaqus SPH结果粒子位置教程

Abaqus软件后处理中SPH粒子的渲染真的太差了，粒子大小、光照等都没有有效的调整方法，于是想着从ODB文件中读取出粒子位置信息，提取出来的信息还可以用于其他后处理及渲染等。

但是有很多问题：1. PC3D粒子的信息怎么输出？与节点关系如何对应？2. 可以在界面上选择节点编号查询节点，但是批量怎么读取？3. fieldOutput中物理量只有位移等，没有当前坐标。

想了一个笨办法：

1. 先在后处理模块中用“Display Group”功能显示所有流体粒子，然后利用查询功能获得所有当前粒子的label，存到一个列表中；
2. 在odb文件的assembly对象中获得上面粒子label编号列表对应的节点的坐标，就是初始坐标。
3. 在steps.frames对象中获得位移场变量，位移是一个容器，里面存储了所有节点的位移，找出流体对应节点的位移，位移与初始坐标的和就是当前坐标。
4. 将查询到的结果（包括节点当前坐标、节点位移值等）按照节点顺序写入vtk文件，利用paraview进行可视化。

如下为读取及写入vtk文件的Python代码

免责声明：本文系网络转载或改编，未找到原创作者，版权归原作者所有。如涉及版权，请联系删