lammps模型晃动问题

大家好，我是小马老师。

在本文中，分享一个关于如何修正LAMMPS模拟中晃动问题的小技巧。

最近，在咨询过程中，多次遇到模型晃动的问题。

具体表现为在模拟过程中，模型无法保持稳定，可能会上下晃动或左右晃动。

这种情况会在后续的压痕或摩擦模拟中导致受力曲线出现剧烈波动。

造成这种现象的根本原因是模型未能弛豫到平衡状态，导致原子间的能量过大，结构无法稳定。

直接原因是初始模型未能接近平衡状态，并且边界层固定得过早。

一般来说，晶体模型的初始状态应该是平衡状态。

如果建模不合理导致原子间距离与平衡间距相差过大，那么在弛豫过程中会自动调整到合理的平衡位置。

然而，如果在弛豫之前就固定了底部的边界，导致底部原子无法调整原子间距，这将导致上层未固定原子无法完全调整到平衡状态，从而出现晃动现象。

解决方法相对简单，即在弛豫过程中先取消固定代码，使用"fix all npt aniso"方式对整个体系进行充分的弛豫。弛豫完成后，再重新分组并固定边界层。通过这种方法，可以完美地解决模型晃动问题。

下图展示了采用新方法后的弛豫结果，结构已经非常稳定。