lammps后处理：堆垛层错和孪晶的数量统计

大家好，我是小马老师。
本文介绍python统计堆垛层错和孪晶数量的方法。
lammps金属变形模拟中，可能会产生堆垛层错或者孪晶等结构，但是ovito免费版本无法统计这些结构的数量。
使用python和ovito软件包，则可以方便的统计出这些结构的数量。

代码如下：

from ovito.io import import_file
from ovito.modifiers import PolyhedralTemplateMatchingModifier, IdentifyFCCPlanarFaultsModifier#定义变量
steps =[]
counts_isf =