大家好，我是小马老师。
在前面的推文中，介绍了一维正弦粗糙界面的建模方法。
本文给出二维正弦粗糙界面建模代码。
二维曲面和一维曲面建模原理一样，只要能够构造出二维曲面的方程即可。
一维曲面沿y方向没有变化，而二维曲面y方向同样按照正弦曲线变化，因此，需要在方程中加入y的变化。
二维曲面方程形式为：
z=Asin(ωx)*sin(ωy)+k
对应的代码为：
variable v equal "v_zz < (15 * sin(v_xx * 2.0*PI/50)*sin(v_yy * 2.0*PI/50) +30 ) "
create_atoms 1 box var v set x xx set y yy set z zz
最终效果为：
公众号：lammps加油站
QQ群：227321877
