lammps后处理：多帧孔洞体积和孔隙率的计算

本文介绍lammps后处理技巧：多帧孔洞体积和孔隙率的计算方法。
在前面的专栏中，已经介绍了单帧孔洞体积的计算方法，有不少粉丝朋友咨询多帧孔洞体积的计算方法。

在上一次案例代码的基础上，稍加修改，添加一个for循环遍历所有的帧即可实现多帧孔洞体积的计算。
计算的结果保存到Vol.txt文本中，便于绘图。
代码如下：

from ovito.io import import_file
from ovito.modifiers import ConstructSurfaceModifier

pipeline = import_file('pka.xyz')
pipeline.modifiers.append