大家好，我是小马老师。
对于一些金属体系的变形模拟，在lammps运算完成后，通常需要将轨迹文件导入到ovito进行分析。
ovito分析晶体变形的方法有很多，比如DXA、CNA等。
在一些论文中，会使用下图的方法分析晶内变形情况：
下面介绍一下这种图像的绘制方法。
（1）ovito打开轨迹文件
以Ni的拉伸模拟结果为例，拉伸轨迹文件导入到ovito中如下图所示：
（2）调用位错分析DXA
“Add modification”-“Dislocation analysis”
调用DXA后，主界面会以不同的颜色显示不同的相，并标出位错线位置。
下面通过删除部分相观察晶内变形情况，如删除FCC相。
（3）添加公式选择要删除的原子
“Add modification”-“Expression selection”
在下面的公式输入框中输入：
(StructureType ==1 || StructureType ==0) && (Position.X <= 105 && Position.X >= 0)
StructureType对应不同的相，例如1为绿色fcc相，0为白色other。
Position.X设置为只选择x坐标位于[0,105]之间的原子，这个根据个人需要设置，可用也可以不用。
（4）删除选定原子
“Add modification”-“Delete selected”
删除选定原子，为显示方便，可关闭“Simulation cell"、“Defect mesh”。
如果显示不清晰，可添加阴影渲染"Ambient occlusion"
这个案例给大家提供一种后处理思路，至于需要删除哪些原子以及保留哪部分数据，要根据模型具体情况具体分析，模型不同或者分析侧重点不同，处理方法可能会不同。
公众号定期分享lammps命令和模拟技巧，感谢分享和关注^_^。
END
公众号案例代码以及学习交流，请加QQ群：754749935
点击关注微信公众号：lammps加油站
