大家好，我是小马老师。
本文介绍lammps六边形结构建模方法。
lammps自带的形状有长方体、圆柱体、圆锥体和球体，其他特殊的情况大多需要用组合建模方法。
如下图截面为六边形的纳米线，很难用一条命令建模。

主要原理和之前的锐角刀具建模原理一样，用region plane设置6个面，用这个6个面组成一个封闭区域，在封闭区域内生成原子。

建模代码：

atom_style atomic
units metal
boundary p p p
#外接圆半径
variable R equal 50
#建盒子
lattice fcc 3.61
region box block -20 20 -20 20 -40 40 
create_box 1 box
#斜面1
variable x1 equal cos(PI*30/180)
variable y1 equal sin(PI*30/180)  
region p1 plane $R 0 0 ${x1} ${y1} 0  side out units box
#斜面2
variable x1 equal cos(PI*30/180)
variable y1 equal -sin(PI*30/180)  
region p2 plane -$R 0 0 ${x1} ${y1} 0  side in units box
#斜面3
variable x1 equal -cos(PI*30/180)
variable y1 equal -sin(PI*30/180)  
region p3 plane -$R 0 0 ${x1} ${y1} 0  side out units box
#斜面4
variable x1 equal -cos(PI*30/180)
variable y1 equal sin(PI*30/180)  
region p4 plane $R 0 0 ${x1} ${y1} 0  side in units box
#平面1
variable y2 equal $R*cos(PI*30/180)
region p5 plane  0 ${y2} 0 0 -1 0  side in units box
#平面2
variable y3 equal -$R*cos(PI*30/180)
region p6 plane  0 ${y3} 0 0 1 0  side in units box
#合并6个面
region piece intersect 6 p1 p2 p3 p4 p5 p6
create_atoms 1 region piece
mass 1 64
write_data all.data

运行结果：