lammps教程：薄膜渗透过滤模拟--平衡弛豫过程

大家好，我是小马老师。
本文继续介绍lammps模拟薄膜过滤的案例代码。
在前面推文中，已经介绍了模型的建模过程。
薄膜建模方法

模型经过minimize能量最小化之后，进入弛豫阶段。
因弛豫阶段代码较长，不再一一介绍各部分代码的功能，仅对代码中的几个比较重要的功能进行详细介绍。
（1）薄膜固定及活塞墙的加压控制
在模拟过程中，过滤薄膜不能移动，使用fix setforce命令对其进行固定。
fix mysf1 membrane setforce 0 0 0
活塞墙在y、z方向需要固定，在x方向允许移动，因此，setforce命令中，x方向参数为NULL，表示不设置该方向受力，y、z方向为0，表示在y和z方向原子受力0，只能沿x方向移动。

fix		mysf2 piston_left setforce NULL 0 0
fix		mysf3 piston_right setforce NULL