大家好，我是小马老师。

本文介绍使用lammps建立球壳模型的方法，如下图所示：

内核为Al球体，外面包裹一层Al球壳。

主要使用region sphere命令。

1.内核建模

使用region sphere命令生成一个球体区域

region ball_in sphere 0 0 0 60 units box

2.球壳建模

分别建立球壳内外两个球面，内球面添加side out关键词，指定球面的外侧区域为我们需要的区域。

region ball1 sphere 0 0 0 70 units box
region ball2 sphere 0 0 0 60 side out units box

使用intersect命令对两个球面求交集，得到shell区域，在shell区域内填充原子

region shell intersect 2 ball1 ball2

3.填充原子

在以上两个区域内使用create_atoms命令填充原子

create_atoms 2 region ball_in
create_atoms 1 region shell

保存结构文件为Al.data，导入到Ovito即可查看建模效果，至此，建模完成。

下面给出完整的建模in文件代码。

units   metal
atom_style atomic
boundary    p p p
lattice fcc 4.05
region box block -20 20 -20 20 -20 20 
create_box 2 box
#创建球壳区域
region ball1 sphere 0 0 0 70 units box
region ball2 sphere 0 0 0 60 side out units box
region shell intersect 2 ball1 ball2
#创建内核区域
region ball_in sphere 0 0 0 60 units box
#填充原子
create_atoms 2 region ball_in
create_atoms 1 region shell
#保存文件
mass    * 27
write_data Al.data

region命令可选功能比较丰富，灵活使用可创建比较复杂的模型。

in文件编写

公众号案例代码以及学习交流，请加QQ群：754749935

扫描关注微信公众号：lammps加油站。

lammps建模：球壳模型建模方法_lammps