lammps模拟——NiTi纳米压痕模拟PH曲线异常

问题描述

将经过升温模拟的NiTi模型作为初始模型进行纳米压痕模拟，模拟结束后发现PH曲线异常，其中有个别数据点异常大，如下图

正常的PH曲线应该如下图

原因分析

经过升温模拟后，原子排布不再是紧密规律的堆叠在一起，原子位置会有些移位。如下图红圈中所示。

所以划分牛顿层、恒温层、固定层时还按照初始基体模型的高来（原始模型高60埃米），就会导致有个别原子没有被包括进去，从而导致PH曲线表现异常。

# 模拟区域
region          fixed_layer block 0 120 0 120 0 9 units box     # 工件的边界区域
region          temp_layer block 0 120 0 120 9 18 units box  # 工件的温度层区域
region          newton_layer block 0 120 0 120 18 60 units box  # 工件的牛顿层区域
region          piece union 3 fixed_layer temp_layer newton_layer

修改后的代码

# 模拟区域
region          fixed_layer block 0 120 0 120 -1 9 units box     # 工件的边界区域
region          temp_layer block 0 120 0 120 9 18 units box  # 工件的温度层区域
region          newton_layer block 0 120 0 120 18 62 units box  # 工件的牛顿层区域
region          piece union 3 fixed_layer temp_layer newton_layer

 可以看到，我在原代码的基础上，将固定层、牛顿层的区域进行了扩大，尽可能的将原子包括进来，从而保证能够得到正确的PH曲线。