本文介绍孔洞体积和孔隙率计算方法。
在拉伸、辐照或者一些力学的lammps模拟中，材料内部会出现孔洞等缺陷。
Ovito中有一个Construct Surface mesh命令，可以识别原子表面，使用这个功能，就可以直观的显示空洞的大小和形态。
ovito免费版可以直观地显示孔洞的表面积，但是不能计算体积，计算体积需要ovito pro版本。
但是，使用python调用ovito模块，则可以很方便地计算孔洞的体积。
下面是一个简单的案例代码：
from ovito.io import import_filefrom ovito.modifiers import ConstructSurfaceModifierpipeline = import_file('vol.xyz')pipeline.modifiers.append(ConstructSurfaceModifier(method = ConstructSurfaceModifier.Method.AlphaShape,radius = 1.8,identify_regions = True))data = pipeline.compute()print(f"空洞体积: {data.attributes['ConstructSurfaceMesh.void_volume']}")
运行这个python代码之后，结果如下，有了孔洞体积就可以计算孔隙率：
这个命令还可以输出其他的量，如下图所示，具体可到ovito官网查看手册说明。
关于python ovito的配置可参考：
也可考虑报名python ovito后处理课程。课程包含：Python ovito模块配置、CNA结构、位错密度、角度分布、键长分布、rdf等计算
扫码订阅课程(↑↑↑)
---------- End -----------
报名一对一培训课程可领完整版《lammps中文教程》，配套学习更轻松！
扫码领取粉丝专属优惠券↓↓↓
